1. Introduction
An outbreak of a coronavirus in late 2019 (COVID-19) caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection quickly became pandemic. Based on the World Health Organization statistics, there have been 771,820,937 confirmed cases of COVID-19 globally by November 2023, including 6,978,175 deaths (WHO Coronavirus (COVID-19) Dashboard). Many studies indicate the persistence of COVID-19 symptoms following recovery of the acute infection despite clearance of the virus from the body [
1,
2,
3]. This phenomenon of the COVID-19 symptoms persistence more than 12 weeks following the infection has been characterized as post-COVID-19 syndrome or long-COVID-19 syndrome.
The long-term neurological and psychiatric complications of COVID-19 infection are being intensively studied. Complications of the COVID-19 include “brain fog”, fatigue, headache, sleep disorders, cognitive impairment, impaired sense of smell and taste, depression, anxiety, sleep disturbances, post-traumatic disorder and obsessive–compulsive symptoms [
4,
5,
6,
7,
8,
9,
10,
11,
12,
13,
14].
Among the complications related to COVID-19, depression represents one of the major public health concerns; it decreases quality-of-life outcomes and may cause disability [
15]. People who survived the acute phase of COVID-19 have been found to be at increased risk of mental health disorders, including being diagnosed with major depressive disorder (MDD) more often [
16]. There are still too few studies to confidently judge whether PCD has any specific features compared to MDD. Some publications described differences in symptoms of PCD and MDD. Simonetti et al. [
17] found high level of excitatory symptoms in post-COVID-19 syndrome, which is more likely corresponded to a mixed form of depression. In contrast, the comparative analysis of biomarkers suggested a similar etiopathogenesis and inflammatory hypothesis for post-COVID depression and MDD [
18].
The etiology of MDD is commonly considered to be multifactorial [
19]. The pathological mechanisms underlying post-COVID depressive symptoms are mainly related to the inflammation triggered by the peripheral immune-inflammatory response to the viral infection [
18,
20]. Several studies have also demonstrated neuroinflammation [
21,
22] in post-COVID syndrome. Neuroinflammation-associated reactivation of microglia and astrocytes is accompanied by the release of cytokines, pro-inflammatory and cytotoxic factors. This change in the extracellular environment negatively affects the functions of myelinating oligodendrocytes and oligodendrocyte precursors [
23,
24,
25]. This leads to the hypothesis that neuroinflammation and subsequent demyelination may be a major cause of post-COVID depression.
Several studies of brain demyelination in post-COVID patients support this hypothesis. White matter damage was seen in 57% of 4342 hospitalized COVID-19 patients [
26]. Reports indicated effect of COVID-19 infection in triggering demyelinating diseases, such as multiple sclerosis [
27], acute disseminated encephalomyelitis [
28], Guillain Barré syndrome [
29,
30]. The neurotoxic action of the virus is mediated by binding to angiotensin-converting enzyme-2 (ACE2) receptors or indirectly by inducing cytokine storm leading to disruption of the blood-brain barrier, immunological mediation, increasing blood coagulation and as a trigger for autoimmune-mediated demyelinating injuries in the CNS [
31,
32]. Recent MRI studies had shown a widespread alterations of white matter (WM) microstructure in the brain after COVID-19 and increase in radial diffusivity is indicative of demyelination and axonal damage [
33,
34,
35,
36,
37,
38,
39].
Considering the similarities between MDD and PCD symptoms [
40], another argument in favor of the significant role of demyelination in PCD is that MDD is often considered as a disorder of WM connectivity, which is caused by changes in brain myelination. Impaired WM integrity and demyelination in patients with MDD has been shown in several neuroimaging studies [
41,
42,
43,
44,
45,
46,
47].
Despite the large number of publications. related to brain demyelination in post-COVID-19 patients, we have not encountered a single study that systematically assessed demyelination in patients with diagnosed clinical depression [
36]. In this work, we applied the quantitative MRI method called macromolecular proton fraction (MPF) mapping to assess brain demyelination as a possible cause of post-COVID depression. MPF mapping is standing out among other MRI methods sensitive to myelin due its improved specificity to myelin, fewer physiological confounders based on diffusion, relaxation, and susceptibility [
48,
49,
50]. MPF has strongest correlations with myelin content assessed with histology [
51,
52,
53,
54,
55]. MPF values are independent of the magnetic field strength and provide solid method for myelin measurements across a variety of human and animal MRI platforms [
56,
57]. Finally, MPF maps can be obtained without modification of original manufacturers’ pulse sequences available on the standard clinical MRI scanners [
58,
59,
60] including for assessing demyelination in mental disorders [
58].
The current study aim was to quantify the degree of brain myelination in patients with clinically diagnosed post-COVID depression. We hypothesize that depression in post COVID-19 period is directly related to the decrease in brain myelination.
2. Materials and Methods
2.1. Study participants
Eighty seven participants were recruited by Mental Health Research Institute (Tomsk, Russia), Medica Diagnostic and Treatment Center (Tomsk, Russia), and Tomsk State University (Tomsk, Russia) between September 2022 and June 2023. The inclusion criteria were the following: age from 18 to 60 years, previous positive to COVID-19 PCR test and persistence of post-COVID complications (except for the control group), the absence of the history of traumatic brain injury, and the absence of any diagnosed neurologic or psychiatric condition prior to COVID-19. The exclusion criteria were: pregnancy, symptoms of acute infectious and somatic diseases, contraindications to MRI, inability to tolerate the MRI procedure, and self-withdrawal from the study. Written informed consent was obtained from all participants. The study design was approved by the local Ethical Committee of the Mental Health Research Institute (protocol №15/8.2022) and Bioethics Committee of Tomsk State University (№12/06.2022) following the guidelines of the Declaration of Helsinki
Patients with post-COVID complications were recruited by a neurologist at the Medica Diagnostic and Treatment Center and a psychiatrist at the Mental Health Research Institute. The control group was formed from employees and students of Tomsk State University, as well as their relatives and friends, who had not prior COVID-19 history. Patients and healthy volunteers were interviewed regarding eligibility for inclusion and exclusion criteria. Those who met these criteria signed an informed consent form and completed the questionnaire (Supplementary Materials). A few days later, participants underwent an MRI scanning. Data acquired from 5 participants who were newly diagnosed with brain pathologies based on MRI (cavernoma, glioma, vascular anomalies) were excluded from further analysis.
The severity of depressive symptoms was assessed using two scales, Hospital Anxiety and Depression Scale (HADS and Hamilton Depression Rating Scale (HDRS) [
61]. The screening testing with HADS was carried out by a clinical psychologist. the HDRS testing was carried out by a qualified clinical psychiatrist during a clinical interview. After that, the participants who scored higher (>8) on the HADS were assessed by a psychiatrist. A group of patients with PCD was formed by a psychiatrist based on a structured clinical interview for the International Classification of Diseases (ICD-10) and baseline assessment report, including socio-demographic characteristics, medical history, questionnaire regarding COVID-19, clinical and psychometric examination. The severity of the current depressive episode was assessed before the start of drug therapy using the Hamilton Rating Scale for Depression (HDRS) [
62,
63]. The total score is interpreted as follows: no depression (0-7); mild depression (8-16); moderate depression (17-23); and severe depression (≥24).
Twenty five individuals with diagnosed clinical depression (moderate depressive episode - F32.1, severe depressive episode without psychotic symptoms – F32.2, recurrent depressive disorder (first diagnosed at the time of the study), current episode moderate – F33.1, according to ICD-10) were included in the post-COVID depression (PCD) group. The participants (n=38) with neurological complications of COVID-19 and without clinical depression were included in the comparison group (noPCD group). The control group (n=19) included healthy volunteers who were not COVID-19 positive and did not experience symptoms of COVID-19 from the start of the pandemic until the time of examination. The demographic characteristics of participants are shown in
Table 1. The groups did not differ significantly in age, gender, and severity of COVID-19 (PCD and noPCD groups) according to Chi-square criteria.
2.3. Patient survey
All participants filled out a questionnaire regarding COVID-19 (Supplementary Materials). The questionnaire included questions about the number, severity, and date of illnesses, the PCR tests, vaccination, symptoms of the acute and post-COVID phases. As symptoms of the acute phase, patients were asked to note the presence or absence of anosmia, ageusia, fever, difficulty breathing, cough, muscle weakness, myalgia, headache, and dizziness. As symptoms of the post-COVID phase, patients were asked to note the presence or absence of headache, dizziness, brain fog, anosmia, ageusia, sensitivity, hypertensia/hypotensia, insomnia, fatigue, attention and memory deficit, myalgia, depression, panic attacks. Based on the results of the answers, the number of symptoms in the acute and post-COVID phases was calculated as the sum of symptoms (1 symptom – 1 point), for which positive answers were given for all diseases. This parameter has proven itself well as a predictor of post-COVID complications and for assessing the severity of post-COVID [
64,
65,
66,
67].
2.2. MRI data acquisition
All participants underwent MR imaging using 1.5T clinical scanner Magnetom Essenza (Siemens, Erlangen, Germany). The fast MPF mapping protocol [
58] included three 3D spoiled gradient-echo pulse sequences with following acquisition parameters:
− Magnetization-transfer-weighted: TR = 20 ms, echo time (TE) = 4.76 ms, flip angle (FA) = 8°, scan time 5 min 40 s;
− T1-weighted: TR =16 ms, TE = 4.76 ms, FA =18°, scan time 4 min 32 s;
− Proton-density-weighted: TR= 16 ms, TE = 4.76 ms, FA= 3°, scan time 4 min 32 s.
In addition, the following sequences were included in the protocol:
− 3D FLAIR: TR = 5000 ms, TE = 390 ms, TI = 1800 ms;
− 3D T1-weighted: TR = 16 ms, TE = 4.76 ms;
− 3D T2-weighted: TR=3000ms, TE=335ms.
All scans were acquired in the sagittal plane with a voxel size of 1.25 × 1.25 × 1.25 mm3 (matrix 192 × 192 × 160, field of view 240 × 240 × 200 mm3), single signal averaging.
The total scanning time was about 35 minutes.
2.3. Image processing
MPF maps were reconstructed using the previously developed software in the C++ language (available at
https://www.macromolecularmri.org/), which implements a single-point algorithm with a synthetic reference image [
68,
69].
Regional WM and GM segmentation was carried out using Advanced Normalization Tools (ANTs) [
70,
71] and Eve anatomical atlas [
72] as described in [
73]. T1 template image of Eve atlas was registered to individual MPF maps using antsRegistrationSyNQuick algorithm. Then the obtained deformation field was applied to Type-III Eve atlas segmentation [
72] to register the template atlas labels to individual MPF maps (
Figure 1).
The measurements on MPF maps were performed for 118 GM and WM structures (including measurements in the right and left hemispheres) using ITK-snap software. The list of structures included:
- 1)
juxtacortical (superficial) WM: superior parietal, superior, middle, and inferior frontal; precentral; postcentral; angular; pre-cuneus; cuneus; lingual; fusiform; superior, inferior, and middle occipital; superior, inferior, and middle temporal; lateral and middle fronto-orbital, supramarginal, rectus, cingulum (parts of cingulate gyrus and hippocampus);
- 2)
WM pathways and fasciculi: corticospinal tract (CST); medial lemniscus; anterior limb, posterior limb, and retrolenticular part of internal capsule (IC); inferior, superior, and middle cerebellar peduncles (CP); cerebral peduncles; posterior thalamic radiation; anterior, superior, and posterior corona radiata (CR); fornix (FX) (stria terminalis, column and body); superior longitudinal (SL) fasciculus; superior (SFO) and inferior fronto-occipital (IFO) fasciculi; uncinate fasciculus; sagittal stratum; external capsule; pontine crossing tract; genu, body, and splenium of corpus callosum (CC); tapetum;
- 3)
subcortical and allocortical GM structures: amygdala; hippocampus; entorhinal area; caudate nucleus; putamen; globus pallidus; thalamus;
- 4)
brainstem structures: midbrain; pons; medulla.
The measurements for left and right hemispheres were averaged for the midbrain, pons, and medulla. Other brain structures were analyzed, taking into account whether they belonged to the left or right hemisphere.
2.5. Statistical analysis
Statistical analysis was performed using Statistica 10.0 software. Differences in MPF between the PCD, noPCD, and control groups for each brain structure were analyzed using the repeated measures analysis of variance (ANOVA) following by post-hoc Fisher LSD tests. One-way ANOVA was used to assess between-group differences in MPF of brainstem structures, results of psychological tests, and parameters associated with COVID-19, such as severity, time after recovery, number of acute and post-COVID symptoms. Post-hoc Fisher LSD tests were performed to clarify differences between groups.
Multiple regression analysis was performed to identify the best predictors of clinical post-COVID depression. First, for all patients of the PCD and noPCD groups for all studied brain structures, the percentage of change in MPF relative to the control group was calculated using the formula: (MPFind - MPFcontrol) / MPFcontrol x 100, where MPFind is the individual MPF value for each brain structure, MPFcontrol is the average MPF for the same structure in the control group. 115 variables were obtained corresponding to the percentage changes in the MPF of 59 structures of the left and right hemispheres. To avoid errors associated with multicollinearity of data for multiple regression analysis, the values of percentage changes in MPF were examined using factor analysis. As a result, principal component analysis allowed to identify 19 independent factors with eigenvalues >1 that explained 87.3% of the variance in total. To interpret the obtained factors, the percentage changes in the MPF of structures with factor loadings >0.7 were considered. The individual factor scores were used for multiple regression analyses. The multiple regression analysis was performed for the PCD and noPCD groups separately and for the total sample of post-COVID patients. HDRS score was used as the dependent variable. Independent variables included 19 variables of individual factor scores, age, gender, time since the acute phase of the first and last disease, the number of symptoms of the acute and post-COVID phase. The quality of regression models was assessed using multiple R correlation coefficient and R2 determination coefficient. The contribution and significance of predictors were assessed using beta-coefficient and p-value.
The differences were considered statistically significant at p<0.05.
4. Discussion
In this work, brain demyelination was investigated as a possible cause of post-COVID depression. Patients with newly diagnosed post-COVID-19 clinical depression showed extensive brain demyelination. Changes in myelination were statistically significant in comparison of the post-COVID patients with depression (PCD) and controls without previous COVID-19 as well as in comparison with the patients with long-term complications after COVID-19 but without diagnosed depression. Patients in the PCD group showed extensive demyelination of the juxtacortical WM, most pronounced in the occipital lobe, but also including the frontal, parietal and temporal lobes. In addition, patients in this group showed demyelination of the WM tracts, the most prominent for the association pathways including the IFOF, sagittal stratum, and left external capsule. Projection (posterior thalamic radiation) and commissural (left tapetum) WM pathways also were affected. Besides, we found GM demyelination, including the hippocampus, putamen, left globus pallidus and amygdala. Significant demyelination was also observed in the noPCD group compared to the control group, but with lower magnitude and smaller affected area than for the PCD group.
Multiple regression analysis revealed only one factor as the main predictor of post-COVID depression: factor 7, which with the highest level of significance was included in both the regression equation for the total sample, which differentiates the PCD and noPCD groups (R2 = 0.41, p < 0.001), and in the regression equation predicting the severity of depression by Hamilton score in the PCD group (R2 = 0.68, p < 0.001). Factor 7 included two structures with the highest factor loadings: IFOF and the uncinate fasciculus of both hemispheres. IFOF also showed significant demyelination in the PCD group compared with both the non-PCD group and the control according to ANOVA results. The uncinate fasciculus showed only a trend (p = 0.08) toward decreased MPF for the PCD group compared with the noPCD group and no significant differences with the control group. Thus, according to the results of two types of analysis, IFOF demyelination can be considered as the best predictor of clinical post-COVID depression.
The regression equation for the total sample did not include any other MPF-related variables other than Factor 7, while the regression equation for the PCD group also included factor 12 and factor 15. Factor 12 was associated with left lingual WM, which showed no significant differences between groups according to ANOVA results. Factor 15 was not strictly associated with any of the studied structures (factor loadings < 0.7 for all structures) and demonstrated weak, although significant, correlations with changes in myelination in multiple regions of juxtacortical WM (superior and middle occipital, middle frontal, angular, inferior and superior temporal WM) and medulla. Most of these juxtacortical WM regions (superior and middle occipital, angular, superior temporal WM) showed significant decrease in MPF in the PCD group compared to controls.
The number of symptoms in the acute and post-COVID phases were also important for group classification. The equations for the total sample and noPCD group included the number of symptoms in the post-COVID phase, and this variable was the only significant predictor for the noPCD group. In contrast, the regression equation for the PCD group included the number of symptoms in the acute phase.
According to the literature, depression is often considered as a disorder of WM connectivity [
41,
42,
43,
44,
45,
46,
47] in which IFOF plays a significant role [
45,
46,
47]. The IFOF connects early visual processing in the cuneus and lingual gyrus as well as parts of the parietal lobe to frontal lobe regions and plays a critical role in semantic language processing, goal-oriented behavior, and visual switching tasks [
74,
75]. In addition, the tract includes the connections between the cingulo-opercular and frontoparietal networks related to executive function and goal-oriented behavior [
75,
76]. High angle diffusion spectrum imaging analysis identified five subcomponents of the IFOF, which primarily included connections from the frontal or fronto-orbital cortex to inferior, superior, and middle occipital lobes [
77]. The IFOF degeneration has been demonstrated in patients with Alzheimer's disease and neuropsychological behavioral disorders, including antisocial personality disorder and obsessive compulsive disorder [
78,
79,
80].
Demyelination of IFOF along with other brain structures has been described in MDD patients [
45,
47,
81,
82]. Lai et al. [
45] found lower fractional anisotropy (FA) in the bilateral IFOF, SLF, inferior longitudinal fasciculi, and CC in MDD patients compared to controls. Changes in FA were found in the left IFOF, uncinate fasciculus, anterior thalamic radiation, and bilateral CC compared to the patients with panic disorder. Liang et al. [
82] identified 3 subgroups of MDD patients based on the spatial localization of reduced FA: the first group with widespread WM disruption (decrease in 8 of 20 studied tracts, including IFOF), the second group with a predominant decreased FA in the CC and left cingulate, and the third group with no statistically significant tract disruption. Reduced FA in the genu of the corpus callosum, IFOF, and posterior thalamic radiation in MDD patients was found by Coloigner et al [
47].
Different studies reported different number of WM-demyelinated brain structures in MDD patients [
42,
43,
44,
45,
47,
81,
82,
83,
84]. Thus, Reppermund et al. [
85] found a significant decrease in FA in 45 brain regions, while the study by Hollocks et al. [
86] found no significant association between our WM parameters and depressive symptoms. Demyelination in MDD patients was observed most often in the CR [
42,
87], IFOF [
42,
45,
47,
81,
82,
87], uncinate fasciculus [
45,
87], posterior thalamic radiation [
47,
85,
88], cingulum [
42,
82,
84,
87], sagittal stratum [
42,
43,
87], IC [
42,
43,
47], and frontal lobe [
44,
85],which was confirmed by our results. Other demyelinated brain structured in MDD were CC [
42,
45,
47,
82,
83,
85], SLF [
85,
87,
88], and FX [
42,
43,
83]; in which we did not find significant changes compared to controls. Perhaps the reason for these discrepancies lies in the differences in the etiology of MDD and post-COVID depression.
The etiology of MDD is commonly considered as multifactorial, in other words, it might be caused by the interaction of biological, genetic, environmental, and psychosocial factors [
19]. Results of our current study point to COVID-19 as the main factor causing recent depressive episode. The multiple regression results identified that the number of symptoms of the acute and post-COVID phase as significant predictors of the presence and severity of clinical depression. Patients in the PCD group were significantly more likely to report ageusia, cough, and headache in the acute phase, as well as anosmia, ageusia, insomnia, fatigue, and attention deficit in the post-COVID phase compared to patients in the noPCD group. At the same time, factors such as age, gender and severity of COVID-19 were not among the significant predictors of post-COVID depression, although previous studies have shown these factors as predictors of post-COVID complications [
64,
67,
89,
90,
91,
92].
In addition to WM demyelination, we found MPF decrease in PCD patients for GM structures: hippocampus, left amygdala, putamen and left globus pallidus. Since all published studies measured only WM myelination, we cannot compare our results describing demyelination of GM brain structures in MDD patients, although some evidence suggests the involvement of the amygdala, hippocampus, and deep GM in depressive disorders [
42,
93]. Application of MPF mapping to myelin quantification was important advantage of our study in comparison with data published with use of DTI and other MRI methods. MPF mapping allows reliable quantification of weak GM myelination [
51,
52,
53,
55,
58] and is independent of iron accumulation in the basal ganglia [
94].
A few studies systematically examined myelination in post-COVID patients [
33,
34,
35,
36,
37,
38,
39]. The works by Huang et al. examined WM changes in the longitudinal MRI (DTI, DKI, NODDI) studies one [
34] and two [
33] years after COVID-19 recovery. They found abnormal diffusion metrics in the corona radiata, genu of the CC and left SLF in one year after recovery and in the CC, CR, CP, IC, posterior thalamic radiation, sagittal stratum, left external capsule, SLF, and CST in two years after recovery. Lower FA in the body of the corpus callosum was observed in the acute phase of COVID-19. Inflammation levels in the acute stage positively correlated with white matter abnormalities and negatively with cognitive function. Qin et al. [
37] revealed significant changes in the volumes of numerous WM structures and in FA in severe compared to mild patients, and in mild patients compared to controls patients. FA differences were found in the following WM tracts: anterior thalamic radiation, SLF, optic radiation, ILF, inferior longitudinal fasciculus, forceps minor, right IFOF, left FX, acoustic radiation, cingulum, and frontal aslant tract. MRI study by Bispo et al. [
35], in patients about 3 months after COVID-19 recovery, showed no changes in GM and lower fiber-specific apparent fiber density in the corona radiata, CST, CC, arcuate fasciculus, cingulate, fornix, IFOF, inferior longitudinal fasciculus, SLF, and uncinate fasciculus. Thus, there was a significant overlap of our results and the literature data regarding a number of brain structures affected by the disease, in particular, the IFOF, cingulum, corona radiata, IC, posterior thalamic radiation, sagittal stratum, external capsule, and uncinate fasciculus. The differences in results can be explained by variability in COVID-19 complications among patients.
We found only one work by Benedetti et al. [
36] that examined associations between the manifestations of post-COVID depression, brain myelination and functional connectivity. The study included voxel-based morphometry, DTI, and resting-state fMRI on 42 patients imaged 3 months after COVID-19. Self-rated depression inversely correlated with GM volumes in anterior cingulate cortex and insula, axial diffusivity, and was associated with functional connectivity. In this study, depressive psychopathology was self-rated on the Zung Self-Rating Depression Scale and high scores (>9) were observed in only 9 of 42 patients. This situation is critically different from our study of a homogeneous group with a clinical diagnosis of depression made by a psychiatrist. Unfortunately, in the study by Benedetti et al. the control group was missing. Our results partially overlap with the above in terms of decreased connectivity in singular WM, but differences in samples and the lack of a control group do not allow us to draw a clear conclusion.