key: cord-0930928-jr2ukudx authors: Benedetti, Francesco; Palladini, Mariagrazia; Paolini, Marco; Melloni, Elisa; Vai, Benedetta; De Lorenzo, Rebecca; Furlan, Roberto; Rovere-Querini, Patrizia; Falini, Andrea; Mazza, Mario Gennaro title: Brain correlates of depression, post-traumatic distress, and inflammatory biomarkers in COVID-19 survivors: A multimodal magnetic resonance imaging study date: 2021-11-02 journal: Brain Behav Immun Health DOI: 10.1016/j.bbih.2021.100387 sha: 3924c9df1b635a019d0dce7b9c03ab106d227cf2 doc_id: 930928 cord_uid: jr2ukudx Psychiatric sequelae substantially contribute to the post-acute burden of disease associated with COVID-19, persisting months after clearance of the virus. Brain imaging shows white matter (WM) hypodensities/hyperintensities, and the involvement of grey matter (GM) in prefrontal, anterior cingulate (ACC) and insular cortex after COVID, but little is known about brain correlates of persistent psychopathology. With a multimodal approach, we studied whole brain voxel-based morphometry, diffusion-tensor imaging, and resting-state connectivity, to correlate MRI measures with depression and post-traumatic distress (PTSD) in 42 COVID-19 survivors without brain lesions, at 90.59 ± 54.66 days after COVID. Systemic immune-inflammation index (SII) measured in the emergency department, which reflects the immune response and systemic inflammation based on peripheral lymphocyte, neutrophil, and platelet counts, predicted worse self-rated depression and PTSD, widespread lower diffusivity along the main axis of WM tracts, and abnormal functional connectivity (FC) among resting state networks. Self-rated depression and PTSD inversely correlated with GM volumes in ACC and insula, axial diffusivity, and associated with FC. We observed overlapping associations between severity of inflammation during acute COVID-19, brain structure and function, and severity of depression and post-traumatic distress in survivors, thus warranting interest for further study of brain correlates of the post-acute COVID-19 syndrome. Beyond COVID-19, these findings support the hypothesis that regional GM, WM microstructure, and FC could mediate the relationship between a medical illness and its psychopathological sequelae, and are in agreement with current perspectives on the brain structural and functional underpinnings of depressive psychopathology. Consistent prospective clinical studies affirm that neuropsychiatric sequelae substantially contribute to the post-acute burden of disease associated with COVID-19 (1, 2) . Over 30% of COVID-19 survivors show clinical depression and post-traumatic distress that persist for months after virus clearance and discharge from the hospital, and leads to the need for psychiatric intervention and psychotropic medication (3) (4) (5) . In the first months after recovery from the acute phase, psychopathology associates with neurocognitive impairment in multiple domains, both predicted by the systemic inflammation during the acute phase and by its pattern of change over time (2, 4, 6, 7) . Beyond COVID, influenza and viral pneumonia are recognized among common physical diseases preceding the onset of a major depressive episode (8, 9) , and persistent lowgrade inflammation in the absence of active medical illness is increasingly recognized as a core biological underpinning of depressive psychopathology in 'endogenous' mood disorders, such as major depressive disorder (MDD) and bipolar disorder (BD), affecting both pathogenesis and therapy outcomes (10) (11) (12) . The brain magnetic resonance imaging (MRI) correlates of depressive psychopathology have been extensively studied (13) . During a major depressive episode, patients show abnormalities in neural response to emotional stimuli (14) (15) (16) (17) (18) , while recovery from depression associates with the normalization of brain cortico-limbic neural responses and functional connectivity in brain areas implicated in the control of emotions and in the cognitive generation of affect (19) (20) (21) (22) (23) (24) . Abnormal neural responses are paralleled by reduced regional grey matter (GM) volumes at MRI voxel-based morphometry (VBM), affecting in particular the anterior cingulate cortex (ACC) (25) (26) (27) (28) , and by abnormal resting state functional connectivity suggesting circuit dysfunction mainly involving the default mode network (DMN), dorsolateral prefrontal cortex, orbitofrontal cortex, ACC, insula, amygdala, hippocampus, basal ganglia, thalamus, and cerebellum (13, 29, 30) . Consistent findings in patients with post-traumatic stress disorder (PTSD) also affirm aberrant GM reduction (31, 32) and resting-state connectivity (33, 34) in corticolimbic structures, largely overlapping with those implicated in mood disorders. Moreover, GM volume reduction and altered connectivity associate with neuropsychiatric symptoms in inflammatory medical conditions (35, 36) , suggesting that regional GM microstructure and function could mediate the relationship between a medical illness and its psychopathological sequelae. A consistent literature also associates mood disorders and PTSD with white matter microstructure as studied with Diffusion Tensor Imaging (DTI), associating both conditions with a decrease in fractional anisotropy of water diffusion (FA) in both conditions, and a decrease in diffusivity along the main axis of the fibers (axial diffusivity, AD) with an increase in mean and radial diffusivity (MD, RD) in mood disorders (37) (38) (39) (40) , spread in large brain regions and including cortico-limbic circuitries. These measures reflect the myelination, orientational coherence, and microtubular axonal structure of fibers, and associate with peripheral markers of inflammation in mood disorders, suggesting disruption of myelin sheaths and axonal pathology (11, 41) . Notwithstanding the consistent literature about brain correlates of neurological symptoms in patients who suffered brain damage during the acute phase of COVID-19 illness, little is known about the brain correlates of depression and post-traumatic distress in COVID-19 survivors (1). Beyond the abnormalities in the olfactory system, systematic reviews showed diffuse cerebral WM hypodensities/hyperintensities, and the involvement of prefrontal, ACC and insular cortex (42, 43): these regions have been consistently implicated in mood and anxiety disorders, thus raising the hypothesis that psychopathological symptoms after COVID-19 could correlate with their structure and function. A very recent study associated PTSD after COVID with reduced volumes of the left hippocampus and amygdala (44) . With a multimodal brain imaging approach, here we studied whole brain VBM, DTI, and resting-state connectivity, to correlate MRI measures with psychopathological ratings in COVID-19 survivors. We studied 50 participants surviving COVID-19, from January 7 to July 13 2021, during an ongoing prospective cohort study on COVID-19 sequelae at IRCCS San Raffaele Hospital in Milan (45) . Diagnosis of COVID-19 was made after clinical and radiological findings at the admission to the Emergency Department, as confirmed by positive real-time reverse-transcriptase polymerase chain reaction from a nasopharyngeal and/or throat swab. To keep a naturalistic study design, exclusion criteria were limited to exceeding age range 18-70, having needed Intensive Care Unit treatment during the acute phase; and the presence of coarse brain disorder prior to the onset of COVID-19, or recent brain lesions typical of COVID such as multifocal cortical infarcts or microhaemorrhagic lesions (46) . Written informed consent was obtained from all participants, and the institutional review board approved the study in accordance with the principles in the Declaration of Helsinki. Depressive psychopathology was self-rated on the Zung Self-Rating Depression Scale (ZSDS) (47) , which thanks to its behavioral and somatic orientation in assessing depressive symptomatology proved a valid instrument to differentiate depressed from non-depressed groups in the general population (48) , including COVID-19 survivors (2, 4, 6, 7, 49) , and on the 13-item Beck Depression Inventory (BDI). Post-traumatic symptomatology was self-rated on the Impact of Event Scale-Revised (IES-R) (50) , which is reliable and valid to measure subjective distress in the dimensions of intrusion, avoidance, and hyperarousal after traumatic experiences (51) , and to detect acute and persistent distress in COVID-19 survivors (2, 4) . MRI was obtained within one week after the clinical assessment. Inflammatory markers at hospital admission and discharge during acute COVID-19 were extracted from charts levels: C-reactive protein (CRP), and systemic immune-inflammation index (SII) (SII = platelets X neutrophils / lymphocytes), which previous prospective research associated with severity of psychopathology up to three months after hospital discharge (4, 6). All imaging was performed on a 3.0 Tesla scanner (Ingenia CX, Philips, The Netherlands) using a 32-channel sensitivity encoding (SENSE) head coil. T1 and T2-weighted FLAIR images were visually inspected as part of the quality check procedure. WM hyperintensities were first observed upon visual inspection and later confirmed via the SPM12 toolbox LST (v. 3.0.0), using the Lesion growth algorithm LGA (52) with the default threshold 0.3. Brain atrophy was observed and quantified via the FSL 6 package Sienax v. 2.6 (53). For VBM, structural T1-weighted MRI scans were acquired at baseline using a T1-weighted magnetizationprepared rapid gradient-echo (MP-RAGE) sequence (TR 2500 ms, TE 4.6 ms, yielding 220 transversal slices with a thickness of 0.8 mm), and processed with the CAT12 toolbox (http://dbm.neuro.uni-jena.de/cat), integrated in the SPM12 software package (Statistical Parametric Mapping, Institute of Neurology, London, UK). Voxel-based morphometry (VBM) statistics were carried out within the general linear model (GLM) framework, to associate markers of inflammation and severity of psychopathology with regional GM volumes. Age, sex, and total intracranial volume were added as nuisance covariates. We searched the whole brain for differences surviving a threshold, family-wise error (FWE) corrected for multiple comparisons, of pFWE<0.05 at the cluster level. For DTI, we used SE Echo-planar imaging (EPI) and the following parameters: TR/TE=5900/78 ms, FoV (mm) 240 (ap), 129 (fh), 232 (rl); acquisition matrix 2.14×2.73×2.30; 56 contiguous, 2.3mm thick axial slices reconstructed with in-plane pixel size 1.88×1.88x2.30 mm; SENSE acceleration factor=2; 1 b0 and 40 noncollinear directions of the diffusion gradients; b value=1000 s/mm 2 . Fat saturation was performed to avoid chemical shift artefacts. Image analyses and tensor calculations were done using the "Oxford Center J o u r n a l P r e -p r o o f for Functional Magnetic Resonance Imaging of the Brain Software Library" (FSL 6.0; www.fmrib.ox.ac.uk/fsl/index.html) (54, 55) . Anisotropy of water diffusion was estimated through the application of diffusion-sensitizing gradients and calculation of elements of the diffusion tensor matrix, i.e. the three eigenvalues 1, 2 and 3 (56-58), representing the tendency to diffuse along the principal direction of the fiber (axial diffusivity, AD, 1); perpendicular to axonal walls and myelin sheaths (radial diffusivity, RD, the average of 2 and 3); independent of tissue directionality (average molecular motion: mean diffusivity, MD, average of 1, 2 and 3). Fractional anisotropy (FA) was calculated as the square root of the sum of squares (SRSS) of the diffusivity differences, divided by the SRSS of the three diffusivities. Estimated FA, eigenvector, and eigenvalue maps (least-square fits) were skeletonized and transformed into a common space as used in Tract-Based Spatial Statistics (TBSS). All volumes were nonlinearly warped to the FMRIB58_FA template, and normalized to the Montreal Neurological Institute (MNI) space, to create a mean FA skeleton and warping individual FA values onto it. The resulting tract invariant skeletons were fed into voxelwise permutation-based cross-subject statistics. Multiple regressions on DTI measures were performed with voxelwise nonparametric permutation-based testing, with markers of inflammation and severity of psychopathology as factors, and with age, sex, and time elapsed from COVID diagnosis to the MRI scan as covariates. Threshold-free cluster enhancement (TFCE) (59) was used to avoid defining arbitrary cluster-forming thresholds and smoothing levels. Data were tested against an empirical null distribution generated by 5000 permutations for each contrast, thus providing statistical maps fully corrected for multiple comparisons across space. Corrected p<0.05 in a minimum cluster size of k=100 was considered significant. Scanning sessions for resting-state fMRI images comprised 200 sequential T2*-weighted volumes (interleaved ascending transverse slices covering whole brain, tilted 30° downward with respect to bicommissural line to reduce susceptibility artifacts in orbitofrontal region), acquired using an EPI pulse sequence (TR = 2000 ms; TE = 30 ms; flip angle = 85°; field of view = 192 mm; number of slices = 38; slice thickness = 3.7 mm; matrix size = 64 x 62 reconstructed up to 96 x 96 pixels). Six dummy scans before fMRI acquisition allowed obtaining longitudinal magnetization equilibrium. Total time acquisition was 6 min and 56 s. We studied whole-brain functional connectivity (FC) patterns with multivariate pattern analysis (MVPA), which can detect neuroimaging-based biomarkers of mood disorders (60) and the effects of biomarkers of brain damage (61) , as implemented in the CONN Toolbox (http://web.mit.edu.sanraffaele.idm.oclc.org/swg/software.htm), to create multivariate connectivity maps for each subject by deriving seeds based on the data, and computing functional connectivity patterns for each voxel via a principal component analysis (PCA) on the connectivity patterns between each single voxel to all other voxels. We performed separate analyses exploring effects of SII, CPR, ZSDS, BDI, IES-R. All models included age, sex, (and days after blood draw for inflammatory markers) as nuisance covariates. Analyses were thresholdet at peak level: p<0.001, uncorrected; cluster level: p<0.05 FWE-corrected). This estimates between-subject variance in connectivity patterns in relation to severity of psychopathology and inflammatory markers, performed across all voxels' connectivity patterns, to identify clusters of voxels that displayed a similar between-subject variance of their spatial connectivity. Post-hoc analysis using significant MVPA clusters as seeds was performed to test which areas of the brain change their connectivity in relation to the independent variables. Considering the sample size, 4 PCA components were kept for each voxel. Finally, to test their effects on severity of psychopathology, and considering the a priori expected significant interaction of independent factors (biomarkers of systemic inflammation, age, sex) and the non-normal distribution of inflammatory biomarkers, independent variables were entered into a Generalized Linear Model (GLZM) analysis of homogeneity of variances with an identity link function (62) . Parameter estimates were obtained with iterative re-weighted least squares maximum likelihood procedures (63, 64) . The quality of the statistical model was checked using the entropy maximization principle of the Akaike information criterion (AIC) (65) . Clinical and demographic characteristics of the studied sample are resumed in Table 1 . Eight patients were excluded from the study after MRI acquisitions: 2 because of multifocal cortical lesions, suggestive of micro-thromboses during acute COVID-19; 3 because of signs of cerebrovascular disease with a diffuse pattern of WM hyperintensities (WMH volumes of the excluded subjects: 8.7 ml, 10.13 ml, 10.34 ml; average WMH volume of the included sample: 0.7±0.9 ml); 2 because of brain atrophy; 1 because of movement artifacts. *************** Table 1 here *************** Severity of current psychopathology associated with regional GM volumes at VBM (Figure 1 ). Both ZSDS, BDI, and IES-R scores negatively associated with GM in large clusters in bilateral ACC, encompassing Brodmann areas (BA) 24 and 32 ( Table 2 ). BDI and IES-R scores also associated with GM in bilateral insular cortex. IES-R also associated with precuneus, the effects being driven by the hyperarousal and intrusion symptom dimensions, which correlated with GM in the same areas, while avoidance showed no effect. *************** Figure 1 here *************** *************** Table 2 here *************** Inflammatory markers in the ED predicted post-acute psychopathology. A GLZM analysis of variance on ZSDS scores showed that in the best fitting model (forward entry/backward removal, AIC criterion) SII at hospital admission positively (b=0.411, Wald W 2 =9.02, p=0.0027), and its decrease during treatment negatively (b=-0.469, Wald W 2 =7.85, p=0.0051) predicted them; accordingly, a highly significant effect was detected for SII at hospital discharge (b=0.337, Wald W 2 =13.07, p=0.0003). Similar opposite effects were detected for IES-R scores: a positive effect of SII at hospital discharge (b=0.479, Wald W 2 =7.967, p=0.0048), and a negative effect of delta SII (b=-0.053, Wald W 2 =3.997, p=0.0456). In all these models, female sex robustly associated with worse severity (ZSDS: W 2 =10.842, p=0.0010; IES-R: W 2 =9.332, p=0.0023). Effects in the same direction were observed for CRP at hospital admission on ZSDS (Wald W 2 =7.954, p=0.0048), but not on IES-R. Inflammatory markers showed no significant effect on VBM, but predicted DTI measures of WM microstructure (Figure 2 ). SII showed significant negative brain widespread associations with measures of water diffusion (AD, MD), with maximal effects in proportionally reducing AD in large portions of the WM skeleton. CRP also negatively associated with AD in spread brain WM tracts, with effects overlapping to those of SII (see detailed effects in Supplementary Table S1 ). *************** Figure 2 here *************** Severity of psychopathology also associated with WM microstructure in the same direction of inflammatory biomarkers. IES-R scores negatively associated with AD in several WM tracts in both hemispheres, with signal peaks in superior and posterior corona radiata and superior longitudinal fasciculus, but also involving Inferior longitudinal fasciculus, external capsule, anterior thalamic radiation. BDI negatively associated with J o u r n a l P r e -p r o o f AD in a small cluster including left superior corona radiata, superior longitudinal fasciculus, and posterior corona radiata. (see detailed effects in Supplementary Table S2) *************** Figure 3 here *************** Inflammatory markers also influenced resting state FC (Table 3, Figure 4 ). MVPA identified 5 seeds where SII significantly associated with rs-FC, located in left dorsolateral prefrontal cortex (salience network), inferior parietal sulcus (dorsal attentional network), cerebellum, right parietal lobe, and precentral gyrus (dorsal attentional network). Severity of inflammation negatively associated with rs-FC between these seeds and target clusters in the default mode network, language and frontoparietal networks; and positively associated with FC in dorsal attentional, salience, and sensorimotor network; with an opposite connectivity pattern for the seed in the cerebellum. CRP was negatively associated with rs-FC between right caudate and two clusters in right temporal cortex in the language network (see detailed effects in Supplementary Table S3 and S4) . Moreover, MVPA identified a seed in the dorsal cingulate cortex where severity of post-traumatic distress (IES-R score) negatively associated with connectivity to medial prefrontal cortex, and positively influenced connectivity to bilateral posterior cingulate cortex/precuneus in DMN, superior temporal cortex in sensorimotor network, and right middle frontal gyrus (see detailed effects in Supplementary Table S5 ). *************** Figure 4 here *************** This is the first study to detect robust associations between anxiety and depression in the post-acute COVID-19 syndrome, inflammatory markers during acute COVID, brain regional GM volumes, DTI measures of WM microstructure, and resting-state functional connectivity. The main findings are that severity of depressive psychopathology associates with decreasing GM volumes in the ACC, while post-traumatic symptoms associate with decreasing GM volumes both, in the ACC and in bilateral insular cortex; and both associate with WM microstructure. Inflammation during acute COVID-19 predicted both, WM microstructure and current psychopathology several weeks after clearance of the virus. Moreover, restingstate FC associated both, with inflammation and psychopathology, confirming the brain functional implication of the strucural effect. These findings support the hypothesis that regional GM and WM microstructure and functional and structural connectivity could mediate the relationship between a medical illness and its psychopathological sequelae, and are in agreement with current perspectives on the brain structural and functional underpinnings of depressive psychopathology. At the core of cortico-limbic structures involved in the control of emotions, in the cognitive generation of affect, and in the attribution of salience, the ACC and the insular cortex have been widely implicated in mood and post-traumatic stress disorders. In healthy subjects the ACC associates with the detection of unfavourable outcomes, response errors, response conflict, and decision uncertainty (66) . Smaller ACC volumes robustly associate with major depressive disorder (25) (26) (27) (28) , but are also proportional to the severity of mild depressive symptoms in the general population, as self-rated on the BDI (67), or observed (68). Smaller ACC and insula GM volumes consistently associate with PTSD (69, 70) , are negatively impacted by the cumulative exposure to stress (71, 72) , and can predict the occurrence of PTSD after mild brain trauma (73) , with insular volumes decreasing over time together with PTSD worsening (74) . Our data now suggest that GM volumes in these structures also associate with post-acute COVID psychopathology, leading to the development of both, depressive and post-traumatic symptoms. Previous prospective studies by our group showed that markers of inflammation during hospitalization predicted anxiety and depression up to three months after virus clearance (4, 6) . Here we confirmed the same effects, with higher SII and CRP at hospitalization predicting the severity of both, persistent depressive symptoms and post-traumatic distress. These same markers predicted widespread lower AD, with overlapping effects in several WM tracts. Effects were spread in large portions of the WM skeleton, as observed in patients with mood disorders (37) (38) (39) (40) , in particular for the effects of inflammatory biomarkers on WM microstructure (41) . Axial diffusivity represents the water diffusivity parallel to the axonal fibers, reflecting the greater freedom of water to diffuse along the principal fiber axis rather than to travel across the surrounding myelin sheaths. Both myelin and axonal microstructure, including microtubules and neurofilaments (75) , contribute to this diffusion anisotropy, with studies on neurodevelopment associating AD with fiber diameter or organization (76) , and animal models associating a reduction of AD with axonal injury (77-79). We previously showed that neurons, and their axons, suffer during COVID-19, as revealed by higher blood levels of neurofilament light chain (NfL), a cytoskeletal intermediate filament protein of central and peripheral neurons which is released during disease progression in a variety of neurological diseases, and which can mark fatal brain damage during COVID-19 (80) . It is then tempting to surmise that our observation of inverse relationship between AD post-COVID, and the severity of previous acute inflammation in the ED, could be the consequence of a slowly recovering axonal damage due to acute COVID-19. Further research and modelling, also using MRI techniques capable of discerning intra-from extra-neurite water diffusion (81) , are needed to clarify this issue. These measures of WM microstructure also associated with severity of psychopathology. Both depressive symptoms and post-traumatic distress, as self-rated on BDI and IES-R, again associated with reduced AD. This observation supports the hypothesis that reduction in axonal integrity or directionality of WM tracts could influence post-COVID psychopathology, by mediating at the level of the brain the effects of inflammation on behavior. In patients with mood disorders the spread disruption of WM microstructure can change in parallel to the severity of depression and ameliorate after successful antidepressant treatment (82) . Further research will clarify if recovery from post-COVID symptomatology will be associated with changes in AD. Effects of inflammation on DTI measures of WM microstructure were paralleled by effects on functional connectivity. Inflammation increased rs-FC within the dorsal attentional and salience networks, and between these networks and the sensorimotor network, which have all been implicated in depression and in PTSD. Interestingly, increased connectivity within Salience and Sensorimotor Networks has been previously associated to persistent physical pain and medical illness (83, 84) , possibly reflecting an increased reliance of body 'threat assessment' (83) and underpinning effects of persistent medical and psychopathological symptoms after COVID. Higher SII also associated with a reduction of rs-FC between (i) DMN, language and frontoparietal networks, and (ii) the salience and dorsal attentional networks. Salience and dorsal attentional network exert an inhibitory influence on the DMN, inducing anticorrelations (85), whereas frontoparietal system has been shown to mediate their interactions (86, 87) . Our results suggest that peripheral inflammatory status may be associated to an inhibition of DMN activity. Reduced rs-FC in this network was consistently associated to psychopathological conditions such as depression and PTSD (33, 88) . In line with these findings, in our sample post-traumatic symptomatology was associated with higher rs-FC between PCC/precuneus in the DMN and the mid portion of cingulate cortex, including dorsal ACC, a crucial brain region of the salience network (89) . Abnormal rs-FC in these regions was proposed as core neural underpinning of with PTSD (33), where increased FC in PCC/precuneus, rated soon after the traumatic experience, also predicted the onset of PTSD and its symptom severity several weeks after trauma (90, 91) . Finally, we detected that high level of CPR was associated to a reduced rs-FC between caudate and temporal cortex. Previous findings detected an inverse association between this inflammatory index and striatum rs-FC among patients with major depressive disorder (92) , this effect may parallel the reduction of striatal dopamine release induced by inflammatory cytokines, as demonstrated in a monkey model of cytokine-induced depression (93) . Strengths of the present study include a focused research question, state-of-the-art imaging methods, a whole-brain study approach, and straightforward effects, but our results must be viewed in light of some methodological limitations. The relative small sample size and the monocentric recruitment in a single ethnic group, could affect the possibility of population stratifications limiting the generalizability of the findings. The widely heterogeneous drug treatments administered during the course of the illness could have influenced MRI measures. We obtained an excellent power to study the association of single psychopathological measures with MRI, but not to investigate their interaction with clinical and pharmacological variables. Finally, a word of caution is needed because of the correlational nature of neuroimaging studies, hampering interpretation of the neurobiological basis of GM and WM MRI measures (94) . In conclusion, these limitations do not bias the main finding of overlapping associations of severity of inflammation during acute COVID-19, brain structure and function, and severity of depression and posttraumatic distress in survivors, thus warranting interest for further study of brain correlates of the postacute COVID-19 syndrome. Beyond COVID-19, our results support current perspectives on the association between inflammation and depression, and on the usefulness of multimodal brain imaging to investigate its underpinnings. J o u r n a l P r e -p r o o f . Effects of inflammatory biomarkers (SII, CRP) and severity of post-traumatic distress (IES-R score) in reducing (blue) or enhancing (red) resting-state functional connectivity from MVPA seeds to target clusters. Table 1 . Clinical and demographic characteristics of the sample; severity of depressive symptoms and posttraumatic distress; and medical comorbidities and cardiovascular risk factors. Table 2 . Results of VBM analysis. In all the listed brain areas, severity of psychopathology, as rated on ZSDS, BDI, and IES-R total and subscale scores, negatively associated with GM volumes. L=left; R=right; L/R= bilateral. Brodmann Post-acute COVID-19 syndrome Anxiety and depression in COVID-19 survivors: Role of inflammatory and clinical predictors. Brain, behavior, and immunity Bidirectional associations between COVID-19 and psychiatric disorder: retrospective cohort studies of 62 354 COVID-19 cases in the USA. The Lancet Psychiatry Persistent psychopathology and neurocognitive impairment in COVID-19 survivors: Effect of inflammatory biomarkers at three-month follow-up 6-month neurological and psychiatric outcomes in 236 379 survivors of COVID-19: a retrospective cohort study using electronic health records. The Lancet Psychiatry Can Cytokine Blocking Prevent Depression in COVID-19 Survivors? Correlation between immune response and selfreported depression during convalescence from COVID-19 Comprehensive textbook of psychiatry Association of seropositivity for influenza and coronaviruses with history of mood disorders and suicide attempts Brain-immune crosstalk in the treatment of major depressive disorder Neuroinflammation in bipolar depression The role of inflammation in depression: from evolutionary imperative to modern treatment target Toward circuit mechanisms of pathophysiology in depression Mood-congruent recall of affectively toned stimuli: A metaanalytic review The neuropsychology of mood disorders The neural basis of mood-congruent processing biases in depression Emotional information processing in mood disorders: a review of behavioral and neuroimaging findings Selective attention to emotional stimuli in a verbal go/no-go task: an fMRI study The neural substrates of affective processing in depressed patients treated with venlafaxine Depression: perspectives from affective neuroscience Neural and genetic correlates of antidepressant response to sleep deprivation -A functional magnetic resonance imaging study of moral valence decision, in bipolar depression Successful antidepressant chronotherapeutics enhance fronto-limbic neural responses and connectivity in bipolar depression Changes in medial prefrontal cortex neural responses parallel successful antidepressant combination of venlafaxine and light therapy Fronto-limbic effective connectivity as possible predictor of antidepressant response to SSRI administration Distinctive neuroanatomical substrates for depression in bipolar disorder versus major depressive disorder Predicting differential diagnosis between bipolar and unipolar depression with multiple kernel learning on multimodal structural neuroimaging Common and distinct patterns of greymatter volume alteration in major depression and bipolar disorder: evidence from voxel-based metaanalysis Gray matter abnormalities in major depressive disorder: a meta-analysis of voxel based morphometry studies All roads lead to the defaultmode network-global source of DMN abnormalities in major depressive disorder Cortico-limbic connectivity as a possible biomarker for bipolar disorder: where are we now? Grey matter reduction associated with posttraumatic stress disorder and traumatic stress Trauma-specific grey matter alterations in PTSD Aberrant resting-state brain activity in posttraumatic stress disorder: A meta-analysis and systematic review Deep Learning Model of fMRI Connectivity Predicts PTSD Symptom Trajectories in Recent Trauma Survivors A multi-modal MRI study of the central response to inflammation in rheumatoid arthritis Investigation into the neural correlates of emotional augmentation of clinical pain Structural white matter changes in adults and children with posttraumatic stress disorder: A systematic review and meta-analysis Recent findings on the role of white matter pathology in bipolar disorder White matter microstructural alterations in posttraumatic stress disorder: An ROI and whole-brain based meta-analysis Widespread white matter microstructural abnormalities in bipolar disorder: evidence from mega-and meta-analyses across 3033 individuals Inflammatory cytokines influence measures of white matter integrity in Bipolar Disorder Brain abnormalities in COVID-19 acute/subacute phase: a rapid systematic review. Brain, behavior, and immunity Brain Imaging in Patients with COVID-19: A Systematic Review. Brain, behavior, & immunity-health Post-traumatic stress symptoms in COVID-19 survivors: a self-report and brain imaging follow-up study Residual clinical damage after COVID-19: A retrospective and prospective observational cohort study Lifting the mask on neurological manifestations of COVID-19 A self-rating depression scale Self-Rating Depression Scale (SDS) Anxiety and depression and its correlates in patients with coronavirus disease 2019 in Wuhan Psychometric properties of the Impact of Event Scale -Revised Impact of Event Scale An automated tool for detection of FLAIR-hyperintense white-matter lesions in multiple sclerosis Accurate, robust, and automated longitudinal and cross-sectional brain change analysis Advances in functional and structural MR image analysis and implementation as FSL Bayesian analysis of neuroimaging data in FSL Looking into the functional architecture of the brain with diffusion MRI Diffusion tensor imaging: background, potential, and utility in psychiatric research MR diffusion tensor spectroscopy and imaging Threshold-free cluster enhancement: addressing problems of smoothing, threshold dependence and localisation in cluster inference Identifying major depression using whole-brain functional connectivity: a multivariate pattern analysis Functional resting-state fMRI connectivity correlates with serum levels of the S100B protein in the acute phase of traumatic brain injury Generalized linear models An Introduction to Categorical Data Analysis An introduction to generalized linear models A new look at the statistical model identification The role of the medial frontal cortex in cognitive control Reduced gray matter volume in the anterior cingulate, orbitofrontal cortex and thalamus as a function of mild depressive symptoms: a voxel-based morphometric analysis Depressive symptoms and neuroanatomical structures in community-dwelling women: A combined voxel-based morphometry and diffusion tensor imaging study with tract-based spatial statistics A systematic review and metaanalysis of magnetic resonance imaging measurement of structural volumes in posttraumatic stress disorder Meta-analysis of 89 structural MRI studies in posttraumatic stress disorder and comparison with major depressive disorder Preliminary evidence for the impact of combat experiences on gray matter volume of the posterior insula Military deployment correlates with smaller prefrontal gray matter volume and psychological symptoms in a subclinical population Smaller regional brain volumes predict posttraumatic stress disorder at 3 months after mild traumatic brain injury. Biological psychiatry: cognitive neuroscience and neuroimaging Trauma-focused psychotherapy response in youth with posttraumatic stress disorder is associated with changes in insula volume Apparent diffusion coefficient on rat brain and nerves intoxicated with methylmercury Diffusion tensor imaging detects and differentiates axon and myelin degeneration in mouse optic nerve after retinal ischemia Noninvasive detection of cuprizone induced axonal damage and demyelination in the mouse corpus callosum Assessment of lesion pathology in a new animal model of MS by multiparametric MRI and DTI Blood neurofilament light chain and total tau levels at admission predict death in COVID-19 patients Along-tract statistics of neurite orientation dispersion and density imaging diffusion metrics to enhance MR tractography quantitative analysis in healthy controls and in patients with brain tumors Changes of white matter microstructure after successful treatment of bipolar depression Salience network functional connectivity is spatially heterogeneous across sensorimotor cortex in healthy humans Brain signature and functional impact of centralized pain: a multidisciplinary approach to the study of chronic pelvic pain (MAPP) network study The hierarchical organization of the default, dorsal attention and salience networks in adolescents and young adults Evidence for a frontoparietal control system revealed by intrinsic functional connectivity Default network activity, coupled with the frontoparietal control network, supports goal-directed cognition Reduced default mode network functional connectivity in patients with recurrent major depressive disorder Dissociable intrinsic connectivity networks for salience processing and executive control Default mode network connectivity as a predictor of post-traumatic stress disorder symptom severity in acutely traumatized subjects A preliminary study of alterations in default network connectivity in post-traumatic stress disorder patients following recent trauma Inflammation is associated with decreased functional connectivity within corticostriatal reward circuitry in depression Chronic interferon-α decreases dopamine 2 receptor binding and striatal dopamine release in association with anhedonia-like behavior in nonhuman primates Natural killer cells protect white matter integrity in bipolar disorder