Participants (see Table 1)
Patients were selected among all the patients from the Geneva University Hospitals (HUG) that showed evidence of a SARS-CoV-2 infection (between March 2020 and May 2021) either by positive polymerase chain reaction (PCR) from nasopharyngeal swab and/or by positive serology while being included according to the exclusion criteria (see below). Patients were divided into three groups and included to the study at 223.07 ± 41.69 days post-infection: 24 patients who had been admitted to ICU during the acute phase of the infection (severe); 42 patients who had been hospitalized but did not require mechanical ventilation (moderate); and 44 patients who had tested positive but had not been hospitalized (mild). Of these patients, 50 agreed to undergo MRI scans (severe: n = 9; moderate: n = 21; mild: n = 20).
The required number of participants in each group was determined by a power analysis involving the comparison of two means. This analysis was based on the literature evaluating the short-term neuropsychological effects of COVID-19 in mild patients 19. To achieve the desired statistical power (1 - β) of 90% and risk of Type I error (α) of 0.05, results indicated that for a one-sided hypothesis, 13 participants would be needed in each group and for a bilateral hypothesis 18. As we planned to perform nonparametric analyses, we had to increase the sample size by 15% 35, resulting in a minimum of 15 participants per group in the case of one-sided hypothesis and 21 participants per group in the case of bilateral hypothesis.
The mild and moderate groups were matched during the screening-inclusion process to the severe group for median age (mild = 57.50 years; moderate = 56.50 years; severe = 60 years), sociocultural level, and clinical variables (except for chronic renal failure) due to a limited number of available patients who were in ICU and met our exclusion criteria. Participants (n = 50) who underwent MRI were not matched during the screening-inclusion process, and all patients that agreed for the MRI study were included. Nevertheless, the groups were still comparable on sociodemographic characteristics (except gender) and severity. Participants were recruited via CoviCare program 36 following patients with post-COVID symptoms in Geneva, Switzerland (MN, OB and IG), as well as from registers from another study (LB). For each patient, we carried out a medical file review, followed by a telephone call inviting the patient to take part in the study, if all the eligibility criteria were met. Exclusion criteria were a history of neurological issues, psychiatric disorders (two of the included participants had had an episode of depression more than 10 years before their SARS-CoV-2 infection), cancer (to exclude possible chemotherapy- and radiotherapy-related cognitive impairment 37), neurodevelopmental pathologies, pregnancy, and age above 80 years (see Fig. 1).
Table 1. Sociodemographic data and medical history
|
Mild
n = 44
|
Moderate
n = 42
|
Severe
n = 24
|
p-value”
|
Mean age in years (± SD)
|
56.57 (± 7.23)
|
56.50 (± 9.58)
|
62.08 (± 12.03)
|
.078
|
Mean education level [1-3] (± SD)
|
2.72 (± 0.45)
|
2.64 (± 0.58)
|
2.50 (± 0.59)
|
.373
|
Gender (% women)
|
34.10
|
35.70
|
20.80
|
.420
|
Handedness (% right handed)
|
97.70
|
92.90
|
95.80
|
.553
|
Mean days of hospitalization (± SD)
|
-
|
12.00 (± 12.87)
|
40.13 (± 32.07)
|
-
|
Diabetes in %
|
2.30
|
9.50
|
20.80
|
.083
|
Smoking in %
|
11.40
|
2.40
|
4.20
|
.206
|
History of respiratory disorders in %
|
11.40
|
11.90
|
25.00
|
.259
|
History of cardiovascular disorders in %
|
13.60
|
14.30
|
25.00
|
.432
|
History of neurological disorders in %
|
0
|
0
|
0
|
1
|
History of psychiatric disorders in %
|
2.30+
|
2.40+
|
4.20+
|
.887
|
History of cancer in %
|
0
|
0
|
0
|
1
|
History of severe immunosuppression in %
|
0
|
0
|
0
|
1
|
History of developmental disorders in %
|
0
|
0
|
0
|
1
|
Chronic kidney disease in %
|
0
|
0
|
8.3
|
.026*
|
Sleep apnea syndrome in %
|
9.10
|
11.90
|
29.20
|
.067
|
Note. ns: not significant; SD: standard deviation; “Statistical analysis performed: Kruskal-Wallis or chi2; + treated depression more than 10 years prior to COVID-19.
General procedure and ethics
A flowchart displaying the successive stages of the study according to the eligibility criteria for each experimental group is provided in Fig. 1.
After being given a full description of the study, participants provided their writtenf informed consent. The study was conducted in accordance with the Declaration of Helsinki, and the study protocol was approved by the cantonal ethics committee of Geneva (CER-02186).
Neuropsychological assessment and other clinical outcomes
The experimental design and tests used are comparable to those used in a previous published study (Voruz et al., 2022). In addition, the detailed tests are available in SI 1.
Symptom validity and presence of noncredible symptoms
First, to validate our neuropsychological measurements, we checked the validity of patients’ symptoms. Both the measurement of symptom validity (i.e., congruence) and the measurement of noncredible symptoms with the BRIEF-A showed good-to-excellent results for all participants, validating the results of the neuropsychological tests and the psychiatric symptom questionnaires.
Neuroimaging processing
Image acquisition
A total of 50 participants (mild: n = 20; moderate: n = 21; severe: n = 9) underwent MRI scans at the CIBM Center for Biomedical Imaging in Geneva, on a Siemens Magnetom PrismaFit 3 tesla scanner. Analysis revealed no significant differences between the mild, moderate and severe groups on age (mild: 55.18 ± 8.58, moderate: 54.94 ± 12.93, severe: 57.80 ± 12.49, p = .885), sociocultural level (mild: 2.76 ± 0.44, moderate: 2.78 ± 0.43, severe: 2.80 ± 0.42, p = .978) or handedness (one left-handed in the mild group), whereas a significant difference was observed for gender (p = .049), with a higher proportion of men in the severe group as compared to mild and moderate. Intergroup analysis also failed to reveal any significant differences either on the interval between infection and MRI (mild: 254.18 ± 39.52 days; moderate: 287.17 ± 45.24 days; severe: 280.80 ± 54.06 days; p = .058) and the interval between neuropsychological testing and MRI (mild: 30.47 ± 20.66 days; moderate: 39.83 ± 26.23 days; severe: 51.39 ± 25.67 days; p = .112). Data from five patients were excluded due to high movement and/or poor registration. Structural images were obtained with a T1-weighted (T1w) magnetization-prepared rapid acquisition gradient echo sequence with an isotropic voxel size of 0.9375 x 0.9375 x 0.9 mm3 (SI 2). Resting-state functional images were acquired through a multiband accelerated echoplanar sequence with an isotropic voxel size of 2.5 x 2.5 x 2.5 mm3, 64 slices, and repetition time of 1 s for a total of 7 min 59 s of acquisition time (480 volumes; SI 3).
Preprocessing was performed using fMRIPrep 20.2.3 38, which is based on Nipype 1.6.1 39.
Anatomical preprocessing
Each T1w volume was corrected for intensity non-uniformity using N4BiasFieldCorrection v2.1.0 40, and skull-stripped using antsBrainExtraction.sh v2.1.0 (using the OASIS template). Spatial normalization to the ICBM 152 Nonlinear Asymmetrical template version 2009c 41 was performed through nonlinear registration with the antsRegistration tool of ANTs v2.1.0 42, using brain-extracted versions of both T1w volume and template. Brain tissue segmentation of cerebrospinal fluid (CSF), white matter (WM) and gray matter was performed on the brain-extracted T1w using fast 43 (FSL v5.0.9).
Functional preprocessing
Functional data were slice-time corrected using 3dTshift from AFNI v16.2.07 44, and motion corrected using mcflirt (FSL v5.0.9 45). This was followed by FLIRT (FSL) coregistration to the corresponding T1w images using boundary-based registration 46 with six degrees of freedom. Motion-correcting transformations, BOLD-to-T1w transformation and T1w-to-template (MNI) warp were concatenated and applied in a single step using antsApplyTransforms (ANTs v2.1.0), configured with Lanczos interpolation.
Physiological noise regressors were extracted with CompCor 47. Principal components were estimated for the temporal (tCompCor) and anatomical (aCompCor) variants. A mask to exclude signal with gray matter origin was obtained by eroding the brain mask, ensuring it only contained white matter and CSF structures. Six tCompCor components were then calculated including only the top 5% variable voxels within that subcortical mask. For aCompCor, six components were calculated within the intersection of the subcortical mask and the union of CSF and WM masks calculated in T1w space, after their projection to the native space of each functional run. Framewise displacement48 was calculated for each functional run using Nipype and volumes with a framewise displacement greater than 0.7 mm were excluded (SI 4).
Many internal operations of fMRIPrep use Nilearn 49, principally within the BOLD-processing workflow. For more details of the pipeline, see the section corresponding to workflows in the fMRIPrep documentation.
Behavioral statistical analyses
We compared the three groups (severe, moderate, mild) on the raw data for each neuropsychological, psychiatric, olfactory, fatigue, and dyspnea variable. Given the nonparametric distribution of the samples (as measured with Shapiro‑Wilks tests), we used nonparametric Kruskal‑Wallis tests. For significant (p < .050) measures, Mann‑Whitney U tests were performed for the 2 ´ 2 comparisons, with false discovery rate (FDR) corrections as function of each domain (cognition, psychiatry) and each Mann‑Whitney pairwise comparison (mild vs. severe; mild vs. moderate; moderate vs. severe).
Neuroimaging statistical analysis
Structural MRI inspection. First, the neuroimaging data were visually analyzed to look for noticeable brain lesions such as microbleeds and WM damages. Groups (SI 6) were compared on the total number of microbleeds and impact on WM, with the Wahlund scale 50.
fMRI statistical analysis. The processed functional time courses were averaged into 156 regions of interest (100 cortical regions 51 that can be associated with 17 resting-state networks 52, 34 cerebellar regions 53 and 22 regions from the basal ganglia (BG) 54 ), and the functional connectivity between pairs of regions was estimated by Pearson correlation. Measures of functional connectivity were converted into z scores with the Fisher z transformation and compared using two-sample t tests to investigate differences between groups. The normality of functional connectivity measures was confirmed with Shapiro-Wilk tests and p values were FDR corrected for multiple comparisons 55.
Relationship between neuropsychological scores and brain connectivity
A partial least squares correlation (PLSC) approach was used to evaluate multivariate associations between neuroimaging and behavioural data 56. This technique estimates latent components that consists out of linear combinations of brain functional connectivity and neuropsychological scores, respectively, to maximize their covariance across participants. The significance of the latent components was evaluated with permutation testing (1000 permutations), and the stability of the feature weights (called saliences) was assessed through bootstrapping (500 samples). Furthermore, we computed the imaging and behavioural loadings defined by the Pearson’s correlation between the original neuropsychological and functional connectivity values, and their corresponding PLSC weights. Only neuropsychological and functional connectivity scores surviving the FDR correction in the intergroup comparison were considered in this analysis. The PLSC analyses were performed on 26 brain and 7 behavioural scores (including 2 covariates to account for the effects of age and sociocultural level) in the three groups (group PLS), and using the myPLS toolbox (https://github.com/danizoeller/myPLS). In post-hoc analyses, the Pearson’s correlations between memory scores and the functional connectivity of specific brain-region pairs were compared across groups.
Data availability
At the end of the COVID-COG project, nonsensitive data will be made available in open access on a dedicated platform.