Abstract
Importance As SARS-CoV-2 infections have been shown to affect the central nervous system, it is crucial to investigate associated alterations of brain structure and neuropsychological sequelae to help address future health care needs.
Objective To determine whether a mild to moderate SARS-CoV-2 infection is associated with alteration of brain structure detected by magnetic resonance imaging (MRI) and neuropsychological deficits.
Design, Setting and Participants Following a case-control design, 223 non-vaccinated individuals with a positive polymerase chain reaction test (PCR) for SARS-CoV-2 obtained between 1 March and 31 December 2020 received MRI and neuropsychological assessments within the framework of the Hamburg City Health Study (median 9.7 months after testing). Two hundred twenty-three healthy controls, examined prior to the SARS-CoV-2 pandemic, were drawn from the main study and matched for age, sex, education and cardiovascular risk factors.
Exposure Infection with SARS-CoV-2 confirmed by a positive PCR.
Main Outcomes and Measures Primary study outcomes were advanced diffusion MRI measures of white matter microstructure, cortical thickness, white matter hyperintensity load and neuropsychological test scores.
Results The present analysis included 223 individuals recovered from mainly mild to moderate SARS-CoV-2 infections (100 female/123 male, age [years], mean ± SD, 55.54 ± 7.07) and 223 matched healthy controls (93 female/130 male, 55.74 ± 6.60). Among all 11 MR imaging markers tested, significant differences between groups were found in global measures of mean diffusivity and extracellular free-water which were both elevated in the white matter of post-SARS-CoV-2 individuals comparing to matched controls (free-water: 0.148 ± 0.018 vs. 0.142 ± 0.017, P<.001; mean diffusivity [10−3 mm2/s]: 0.747 ± 0.021 vs. 0.740 ± 0.020, P<.001). Classification accuracy for detecting post-SARS-CoV-2 individuals based on diffusion imaging markers was up to 80%. Neuropsychological test scores did not significantly differ between groups.
Conclusions and Relevance Our findings suggest that subtle changes in white matter extracellular water content may last beyond the acute infection with SARS-CoV-2. However, in our sample, a mild to moderate SARS-CoV-2 infection was not associated with neuropsychological deficits, significant changes in cortical structure or vascular lesions several months after recovery. External validation of our findings and longitudinal follow-up investigations are needed.
Question Do individuals who recovered from mild to moderate SARS-CoV-2 infection show differences in brain structure and neuropsychological assessments compared to matched healthy controls?
Findings In this case-control study, individuals recovered from SARS-CoV-2 showed significant alterations of the cerebral white matter identified by diffusion weighted imaging, such as global increases in extracellular free-water and mean diffusivity compared to healthy controls. No differences in performance in neuropsychological tests were detected.
Meaning Despite the observed brain white matter alterations in our sample, a mild to moderate SARS-CoV-2 infection was not associated with worse cognitive functions within the first year after recovery.
INTRODUCTION
As the number of patients recovering from an acute infection with the severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2) grows, the study of its long-term consequences on health outcomes has gained much attention.1–4
It is widely recognized that coronavirus disease 2019 (COVID-19) caused by SARS-CoV-2 not only leads to respiratory dysfunction, but also impacts various other organ systems during the acute phase and well beyond.1,5,6 Neurological symptoms, such as headache, fatigue, memory and attention deficits, may significantly impede well-being in individuals suffering from COVID-19 sequelae.4,7,8 Advancing our understanding of the underlying pathological mechanisms will be crucial for addressing future health care needs.
Different potential mechanisms have been suggested to be involved in the development and persistence of neurological symptoms in patients with COVID-19. Post-mortem histopathological and molecular studies have demonstrated viral neurotropism, signs of neuroinflammation,9,10 neurodegeneration,11 demyelination,12 axonal disruption,13 as well as micro- and macrovascular damage.14,15 However, most studies were conducted in patients with severe COVID-19, whereas histopathological findings from individuals with mild to moderate courses are lacking.
In vivo studies applying modern brain imaging joined by comprehensive clinical and neuropsychological assessment are scarce. Recent preliminary evidence from the UK Biobank suggests cortical thickness reductions in the olfactory and limbic network, as well as neurocognitive decline in former COVID-19 patients, although these findings still need to be replicated in independent datasets.11 The majority of remaining studies focused on visually apparent pathological findings such as intracranial hemorrhage, stroke or white matter hyperintensities in small case series or single case reports of more severely affected patients.16–19 Taken together, current evidence is of limited transferability to patients with a mild to moderate SARS-CoV-2 infection, therefore necessitating further investigations.
In order to address this research need, we studied 223 non-vaccinated individuals in median 289 days after recovery from mainly mild to moderate SARS-CoV-2 infections in a retrospective, cross-sectional case-control design. We leveraged advanced magnetic resonance imaging (MRI) techniques enabling the study of imaging phenotypes associated with neurodegeneration, atrophy, myelin/cellular disruption, inflammation, as well as vascular damage.20–23 Moreover, study participants received a comprehensive clinical and neuropsychological assessment. Building upon our previous multi-organ assessment in this cohort,1 here, we provide a detailed in vivo assessment of the cerebral white and gray matter, as well as neuropsychological outcomes in former COVID-19 patients.
METHODS
Study population
We examined participants of the Hamburg City Health Study (HCHS) COVID Program. A detailed description of the study design has been published previously.1 Post-SARS-CoV-2 participants (1) had a positive polymerase chain reaction (PCR) test for SARS-CoV-2 and (2) were aged between 45 and 74 at inclusion. Recruitment routes included both invitation upon laboratory-confirmed SARS-CoV-2 infection and self-referral of participants following newspaper announcement. Subsequent to recruitment, the participants underwent the study protocol of the HCHS24 – including cranial MR imaging, neuropsychological testing and a self-report questionnaire on COVID-19-associated symptoms. In addition, a healthy control group was sampled from the original HCHS cohort.24 The previously reported matching procedure1 was refined to account for confounds known to affect markers derived from structural and diffusion MR imaging: the groups were matched for age, sex, years of education as well as for the prevalence of arterial hypertension, diabetes mellitus, dyslipidemia and smoking behavior by propensity scores using the matchit (v4.3.3)25 R package.
Ethics approval
The local ethics committee of the Landesärztekammer Hamburg (State of Hamburg Chamber of Medical Practitioners, PV5131) approved the study and written informed consent was obtained from all participants.1,26
Clinical assessments
Cognitive testing was performed by a trained study nurse and included the Mini Mental State Examination,27 Trail Making Test A and B,28 Verbal Fluency and Word List Recall subtests of the Consortium to Establish a Registry for Alzheimer’s Disease Neuropsychological Assessment Battery (CERAD-Plus),29 as well as the clock drawing test.30 Psychosocial symptom burden was evaluated by the General Anxiety Disorder (GAD)-731 (anxiety) and Patient Health Questionnaire (PHQ)-9 (depression).32 Moreover, self-reported neurological symptoms (headache, dizziness, fatigue, sleep disturbances) were assessed by part of the PHQ-15.33
Brain imaging
Image acquisitions have been described in detail before.26 Put briefly, 3D T1-weighted rapid acquisition gradient-echo sequence (MPRAGE, 0.83×0.83×0.94mm), 3D T2-weighted FLAIR (0.75×0.75×0.9mm) and single-shell diffusion MRI (2×2×2mm, 64 noncollinear gradient directions, b=1000 s/mm2) were acquired on a single 3T Siemens Skyra MRI scanner (Siemens, Erlangen, Germany). Detailed parameters can be found in the Supplement.
An overview of the derived imaging markers for the grey and white matter can be found in Figure 1. For a detailed account on image preprocessing, derivation of morphometric and diffusion indices, as well as quality assurance (QA) please refer to the Supplement.34 All code is publicly available on GitHub (links can be found in the Supplement).
Following preprocessing, we derived conventional diffusion tensor imaging (DTI) markers of white matter microstructure, i.e., fractional anisotropy (FA) and mean diffusivity (MD), which have been extensively used in neuroscientific and neuropsychological research.35,36 Free-water imaging was employed to model an extracellular free-water compartment, sensitive to immune activation37 and atrophy38, as well as a cellular tissue compartment (FA of the tissue, FAT), more closely reflecting myelin and axonal alterations than their DTI equivalents.39 Fixel-based analysis, a novel multi-tissue model addressing more complex white matter compositions, was used to derive metrics of fiber density, fiber-bundle cross section (FC), fiber density and cross section (FDC), and complexity.40 For further statistical analysis, diffusion markers were averaged across a representative skeleton of the entire white matter derived by tract-based spatial statistics.41 Finally, peak width of skeletonized mean diffusivity (PSMD), a surrogate marker of cerebral small vessel disease, was calculated.23
After cortical surface reconstruction with the Computational Anatomy Toolbox for SPM (CAT12), mean cortical thickness was estimated as a proxy for neurodegenerative processes.20,42,43 Normalized volumes of white matter hyperintensities (WMH load) were obtained by FSL’s Brain Intensity AbNormality Classification Algorithm (BIANCA) with LOCally Adaptive Threshold Estimation (LOCATE).44,45
Statistical analysis
All statistical analyses were conducted in python 3.9.1,46,47 CAT12,42,43 as well as mrclusterstats48. Statistical tests were two-sided, with a P<0.05 as significance threshold. In the case of averaged imaging and clinical data, P-values were adjusted by Bonferroni-correction. Additionally, sensitivity analyses were performed by (1) excluding post-SARS-CoV-2 individuals who had been hospitalized, and (2) stratifying the post-SARS-CoV-2 group by recruitment strategy, following the same procedures as described below.
Phenotypical data
Sample characteristics were compared between healthy controls and post-SARS-CoV-2 participants using Χ2-tests (binary) and two-sample t-tests (continuous). Clinical variables were compared between groups in separate analyses of covariance (ANCOVA) adjusted for age, sex, and education.
Imaging
Statistical analysis of imaging parameters was conducted in two steps. First, global measures, i.e., mean skeletonized diffusion parameters, mean cortical thickness, WMH load and PSMD, were compared between post-SARS-CoV-2 individuals and healthy controls in separate ANCOVAs, adjusted for age, sex and education. In the case of FC and FDC, total intracranial volume served as an additional covariate. Next, in an effort to interrogate spatial patterns of brain structural alterations associated with a mild to moderate SARS-CoV-2 infection, we performed whole-brain voxel-wise permutation-based testing for skeletonized diffusion markers. Utilizing the same design matrices as in the ANCOVAs, we employed 5000 permutations, threshold-free cluster enhancement and family-wise error correction across multiple hypotheses. Vertex-wise cortical thickness was statistically analyzed in a general linear model as implemented in CAT12 with family-wise error correction and a cluster threshold of 10.
Machine learning prediction
To further evaluate their predictive capacities, all brain imaging markers calculated in the study were averaged within regions of interest where applicable (Desikan-Killiany cortical atlas parcels and TractSeg-derived anatomical white matter tracts) and propagated to a comparative supervised machine learning pipeline (scikit-learn v1.0.2).49–51 Per marker, multivariate logistic regression models were trained to predict past COVID-19. Models were scored with prediction accuracy and statistical significance was assessed via comparison to null model predictions. Further details are provided in the Supplement.
RESULTS
Sample characteristics
Imaging data was available for 230 post-SARS-CoV-2 individuals. Following quality assessment (QA), in total 7 post-SARS-CoV-2 individuals had to be excluded, leaving 223 cases for propensity-score matching with healthy controls who had passed QA. Results of the matching procedure are visualized in Figure S1. The final sample included 223 matched controls (93 female, age in years, mean ±standard deviation [SD], 55.74 ± 6.60) and 223 post-SARS-CoV-2 individuals (100 female, 55.54 ± 7.07, see table 1). Of the latter, the majority had a mild to moderate course of COVID-19 (without symptoms, n=7; mild symptoms, n=125; moderate symptoms, n=67), 18 were hospitalized and none required mechanical ventilation or intensive care unit treatment. There were no significant differences between post-SARS-CoV-2 individuals and matched controls regarding age, sex, years of education and cardiovascular risk factors.
Clinical data
Although post-SARS-CoV-2 individuals showed better test performances in VF, MMSTE, and CDT compared to matched controls, no significant group differences in neuropsychological test scores remained after Bonferroni-correction for multiple comparisons (Table 2).
Imaging
Post-SARS-CoV-2 individuals exhibited higher global free-water and MD relative to matched controls (mean ± SD, free-water: 0.148 ± 0.018 vs. 0.142 ± 0.017, F=18.47, Pbonf<.001; MD [10−3 mm2/s]: 0.747 ± 0.021 vs. 0.740 ± 0.020, F=17.28, Pbonf<0.001) (Figure 2 [left], Table S1). While PSMD (Puncorr=.005) and cortical thickness (Puncorr=.01) were nominally increased in post-SARS-CoV-2 individuals, both measures, as well as the remaining averaged imaging markers were not significantly different between groups after Bonferroni-correction (Figure S2, Table S1).
To detect spatial patterns of brain structural alterations, we additionally performed vertex- and voxel-wise analyses of gray and white matter imaging markers. Vertex-wise comparisons of cortical thickness did not reveal significant differences between matched controls and post-SARS-CoV-2 participants. Voxel-wise statistics on the entire white-matter skeleton revealed predominant free-water and MD increases in the white matter skeleton of post-SARS-CoV-2 subjects encompassing all brain lobes, comparing to very localized changes in other diffusion markers (Figure 2 [right], Figure S3, Table S2). More specifically, conventional diffusion tensor imaging markers FA and MD showed significant differences between groups, with FA increases in 0.8% and decreases in 1.2% of the skeleton in cases relative to healthy controls. MD was significantly increased in 41.3% and decreased in 0.1% of the skeleton of post-SARS-CoV-2 participants. Employing free-water imaging, post-SARS-CoV-2 individuals showed significant free-water elevations in 38.3% and reductions in 0.4% of the skeleton, as well as FAT elevations in 3.3% of the skeleton, but no FAT reductions. Alterations of the remaining diffusion markers were of even less spatial extent (<3%).
Prediction of a past SARS-CoV-2 infection based on imaging markers
We examined the predictive capacity of derived imaging markers employing a supervised machine learning approach (Figure 3). Free-water (80.21%) and MD (79.38%) achieved the strongest median prediction accuracies. The median cortical thickness score was 45.95%. All investigated metrics but cortical thickness scored significantly better than null models for which the group assignment was randomly permuted.
Sensitivity analyses
Analysis results remained stable if formerly hospitalized post-SARS-CoV-2 individuals were excluded and if post-SARS-CoV-2 individuals were stratified by recruitment route (see Supplement).
DISCUSSION
We investigated brain structural alterations and neuropsychological sequelae in a large sample of individuals who recovered from mainly mild to moderate COVID-19. At median 289 days after the acute infection, these individuals showed significantly higher average free-water and MD in the white matter compared to matched healthy controls. In contrast, cortical thickness and markers of cerebral small vessel disease were not significantly different between groups. In addition, white matter diffusion indices successfully predicted a past SARS-CoV-2 infection. We did not detect neuropsychological deficits post-SARS-CoV-2 individuals. Collectively, our study suggests that a mild to moderate SARS-CoV-2 infection is associated with subtle microstructural alterations in the cerebral white matter beyond the stage of acute infection.
A key aspect of COVID-19 neuropathology appears to be the dynamic response of the intrathecal immune system to the virus. Evidence of neuroinflammation was reported in histopathological and clinical studies of COVID-19 patients: virus invasion,52,53 activation of glial cells,9,54,55 and a cytokine response in the cerebrospinal fluid accompanying neurological and psychiatric COVID-19 symptoms.10 In our study, we observed widespread increases of extracellular free-water and MD in post-SARS-CoV-2 individuals encompassing all brain lobes. Both free-water and MD are sensitive to an activated immune response causing excessive extracellular free-water and thus increased diffusivity.56–58 More specifically, microglia and astrocytes emit cytokines upon activation, inducing osmosis of water from the blood into the extracellular space.59,60 Interestingly, endothelial dysfunction and subsequent vascular leakage due to persistent immune activation has been previously implicated in the pathophysiology of COVID-19.61,62 Taken together, it is conceivable that the observed increase in free-water and MD could be an indirect sign of a prolonged neuroinflammatory reaction to a SARS-CoV-2 infection. Nevertheless, other possible mechanisms for changes in the extracellular space need to be considered.
Volume increases in the extracellular compartment might be accompanied by structural damage like demyelination as well as axonal disruption secondary to neuroinflammation. Free-water corrected diffusion markers enable further guidance in microstructural interpretations. While analyses of overall mean value showed no significant group differences in free-water corrected FAT, voxel-wise statistics identified increased FAT in corresponding frontal areas of post-SARS-CoV-2 individuals. Yet, these changes included only ∼3% of the white matter skeleton, indicating either subtle or spatially limited effects localized to association tracts. Normal to increased FAT in the presence of elevated free-water suggests minor microstructural alterations like axonal compression or displacement rather than damage to myelin sheaths or axons which would rather lead to FAT decreases.21 Moreover, fixel markers, which also model properties of the tissue compartment,63,64 did not show group differences averaged across the entire white matter skeleton. Thus, in contrast to previous work demonstrating more widespread FA reductions in small samples of hospitalized COVID-19 cases,65–67 our findings suggest that white matter changes following a mild to moderate SARS-CoV-2 infection most likely reflect subtle increases in extracellular free-water as opposed to structural neural damage.
Based on previous histopathological reports of vasculopathy in COVID-19 and higher ACE-2 expression in cells of the blood-brain-barrier, we hypothesized that post-SARS-CoV-2 individuals would show alterations in imaging markers of small vessel disease burden.14 However, in our study, WMH load was not significantly different, indicating that a mild to moderate course of COVID-19 does not lead to visually accessible vascular lesions (WMH) as previously reported.68 PSMD, an established imaging marker of small vessel disease more sensitive to microstructural changes,23 showed nominally higher values in post-SARS-CoV-2 individuals, but differences did not survive Bonferroni-correction. Taken together, follow-up investigations are needed to understand whether subtle long-term vascular impairments will eventually increase the prevalence of cerebrovascular disease among COVID-19 convalescents.61
Alterations of cortical gray and white matter commonly co-occur in neurological diseases.69,70 Notably, this was not the case in our study. This is contrasted by a recent report on mildly-affected COVID-19 subjects in the UK Biobank demonstrating longitudinal volumetric reductions in the gray matter in olfactory networks.11 Besides general differences in recruitment strategies, this discrepancy might be due to our cross-sectional study design not being sufficiently sensitive to capture the subtle intraindividual changes.
By providing scores of prediction accuracy, our machine learning analysis evaluated brain imaging markers for their diagnostic relevance. Logistic regression models based on free-water and MD achieved a considerable accuracy of ∼80% in predicting a past SARS-CoV-2 infection, outperforming other imaging markers under study. Of note, cortical thickness achieved the lowest accuracy, not significantly differing from prediction by chance. The difference in accuracy between diffusion metrics and cortical thickness might highlight that, in mild to moderate COVID-19, pathophysiological aspects are better detected by diffusion imaging-based techniques. Finally, the higher accuracy of fiber cross-section compared to cortical thickness – both morphometric measures – might imply that COVID-19-associated alterations preferably occur in the white matter.
We want to emphasize, that our results represent average effects, i.e., not every mild to moderate affected COVID-19 patient may exhibit the reported changes. In addition, our results are based on a non-vaccinated cohort. As vaccination has been repetitively demonstrated to be a highly effective measure against COVID-19, vaccinated patients possibly exhibit less of the pathophysiological substrates identified in our study.71
It is important to put the observed brain white matter alterations into a clinical perspective. Reported persisting clinical sequelae of COVID-19 include executive dysfunction, anxiety, depression, fatigue, muscle weakness and sleep impairment.11,72–74 In contrast, we found no significant difference for any cognitive domain, depression, anxiety and neurological symptoms between groups. Besides the absence or mild expression of the respective symptoms, other reasons, such as the relatively long follow-up period in our investigation, as well as differing degrees of social deprivation as a result of country-specific pandemic control measures, may explain the discrepancy with other studies. Clearly, more research is needed to increase our understanding of factors underlying the persistence of neurological symptoms in a sub-group of COVID-19 patients.
Strengths of this work lie in its considerable sample size; high quality imaging and phenotypical data; a robust and reproducible image processing pipeline; the investigation happening at an early stage of the pandemic, potentially alleviating the problem of different COVID-19 strains and vaccinations as confounders; and a conservative statistical correction scheme to reduce the false-positive rate.
However, our study also exhibits limitations. Our investigation lacks information about SARS-CoV-2 strains as well as precise disease severity stratification beyond the self-reported information on hospitalization and subjective perception of disease intensity. In addition, we follow a cross-sectional study design unable to fully address premorbid group differences despite a rigorous matching procedure. Future longitudinal studies could not only elaborate on the trajectory of the identified microstructural white matter alterations, but also address the question whether these findings are markers of an increased susceptibility for the development of neurological sequelae.
CONCLUSION
We performed a comprehensive assessment of established neuroimaging markers for structural neural integrity to characterize neurobiological changes potentially underlying post-acute COVID-19 neuropsychological sequelae after a mainly mild to moderate disease course. Our findings support the notion of a prolonged neuroinflammatory response indicated by subtle, but widespread increases in extracellular free-water and mean diffusivity in the white matter of COVID-19 convalescents. In contrast, we did not observe signs of cortical atrophy or macrostructural vascular damage. Importantly, despite identifying this characteristic imaging footprint, the investigated sample exhibited no marked neuropsychological symptoms 10 months after SARS-CoV-2 infection. External validation and longitudinal investigations are needed to further clarify the clinical relevance of our findings.
Data Availability
All data produced in the present study are available upon reasonable request to the authors.
Competing interests
JG has received speaker fees from Lundbeck, Janssen-Cilag, Lilly, Otsuka and Boehringer outside the submitted work. JF reported receiving personal fees from Acandis, Cerenovus, Microvention, Medtronic, Phenox, and Penumbra; receiving grants from Stryker and Route 92; being managing director of eppdata; and owning shares in Tegus and Vastrax; all outside the submitted work. CG reports personal fees from Abbott, Amgen, Bayer, Boehringer Ingelheim, Daiichi Sankyo, GlaxoSmithKline, Novartis, Prediction Biosciences, outside the submitted work. GT has received fees as consultant or lecturer from Acandis, Alexion, Amarin, Bayer, Boehringer Ingelheim, BristolMyersSquibb/Pfizer, Daichi Sankyo, Portola, and Stryker outside the submitted work. The remaining authors declare no conflicts of interest.
Acknowledgments
We describe contributions to the paper using the CRediT contributor role taxonomy. Conceptualization: M.P, F.L.N., G.T., B.C.; Data Curation: M.P., F.L.N., C.M., M.S.; Formal analysis: M.P., F.L.N.; Investigation: M.P., F.L.N., C.M., M.S., E.P., S.K., J.G., J.F., O.P., J.M., M.G., C.G., R.T, G.T., B.C.; Methodology: M.P., F.L.N.; Software: M.P., F.L.N., O.P., C.M.; Supervision: G.T., B.C.; Funding Acquisition: G.T., B.C.; Visualization: M.P., F.L.N.; Writing—original draft: M.P., F.L.N.; Writing—review & editing: M.P., F.L.N., C.M., M.S., E.P., S.K., J.G., O.P, M.G., J.M., J.F., R.T., C.G., G.T., B.C.
Data access: M.P., F.L.N., G.T. and B.C. had full access to all the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis.
The authors wish to acknowledge all participants of the Hamburg City Health Study and cooperation partners, patrons and the Deanery from the University Medical Center Hamburg—Eppendorf for supporting the Hamburg City Health Study. Special thanks applies to the staff at the Epidemiological Study Center for conducting the study. The participating institutes and departments from the University Medical Center Hamburg-Eppendorf contribute all with individual and scaled budgets to the overall funding. The Hamburg City Health Study is also supported by Amgen, Astra Zeneca, Bayer, BASF, Deutsche Gesetzliche Unfallversicherung (DGUV), DIFE, the Innovative medicine initiative (IMI) under grant number No. 116074 and the Fondation Leducq under grant number 16 CVD 03., Novartis, Pfizer, Schiller, Siemens, Unilever and “Förderverein zur Förderung der HCHS e.V.”. This work was also funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation – SFB 936 – 178316478 – C2 (M.P., F.L.N., C.M., G.T., B.C.) & C7 (J.G., S.K.). This publication has been approved by the Steering Board of the Hamburg City Health Study.