Participants were recruited from a monitoring program for essential workers who work in construction and law enforcement agencies 24 with an existing neuroimaging program 25, 26. In this program, participants receive an annual personalized monitoring examination that assesses physical, cognitive, and mental health problems. In addition, during the COVID-19 pandemic, the program began to determine COVID-19 status using a barrage of text messaging, phone calls, and letters, and by taking both virtual and in-clinic infection and vaccination histories whenever possible 27. In 2020, it became obvious that some individuals were experiencing persistent symptoms for months after initial infections like brain fog and difficulty concentrating so we developed a clinical neuroimaging study to determine whether these neuropsychiatric symptoms being reported alongside PASC were associated with indicators of white matter health. Participants were asked to participate in the neuroimaging program if they fit the eligibility criteria.
For all participants, we collected history of clinical severity (asymptomatic, mild, moderate, and severe), types of acute symptoms during the index COVID-19 case linked with the neuro-PASC diagnosis, the length of time since COVID-19 infection when an infection was present, and the duration of neurological PASC. Participants self-reported their COVID-19 status during their initial phone screening, and we verified information to the best of our ability; we completed a blood draw for COVID-19 antibody testing with all participants. We acquired copies of other positive antibody, PCR, antigen, and molecular test results acquired outside of our clinic, most frequently from local urgent care facilities. We used the date of infections to calculate months between COVID-19 infection and neuroimaging among those who reported a COVID-19 diagnosis.
Neuroimaging analyses included participants who were recruited from three separate groups of participants: neurological PASC: a history of positive COVID-19 testing (antigen, molecular, PCR, or antibody) from 2020–2021 paired with neuropsychiatric symptoms of brain fog and difficulty concentrating lasting ≥ 6 months. Controls included both “uninfected controls” who had a history of negative COVID-19 test results and no COVID-19 symptoms before imaging and “acute COVID-19” controls who reported positive COVID-19 test results but who had an asymptomatic, mild, or moderate presentation of COVID-19 from 2020–2021 but had no evidence of PASC.
We linked COVID-19-related information to demographic and clinical data from our clinical database. All participants were male so demographics included age, education in years, and body mass expressed in kg/m2. We measured premorbid crystallized cognition using the Wide Range Achievement Test. As physical functional markers, we included measures of maximal handgrip strength, in 100s of pounds, comfortable walking speed over 12 feet averaged across two trials, and chair rise speed averaged over five trials expressed in rises per second. As measures of fluid cognition, episodic memory was measured using the total score on the Hopkin’s Verbal Language Test, processing speed was measured using the Trails B test and expressed in connections per second, language function was measured using the Boston naming test, working memory was measured using the Symbol Digit Modalities Test, and Verbal Fluency was measured using the Controlled Word Association Test.
Eligibility Criteria
Due to occupational inclusion, approximately 94% of program participants are male. Since COVID-19 was more severe in men, and because the measures used here are sensitive to participant sex, only males were eligible for this study. Participants were excluded if, during a pre-imaging screening visit, they met any additional exclusionary criteria as discussed below. We were able to verify infection and clinical details for 84.5% of potential COVID-19 cases. Participants whose case status was verified needed to have a body mass index ≤ 40 kg/m2 to fit comfortably into the MRI scanner. Participants also needed the capacity to provide informed consent, and a willingness to undergo a blood draw totaling about 100 ml of blood. Participants were fluent in English for neuropsychological testing, and there was no upper age restriction for inclusion, though participants need to be ≥ 18 years of age.
Participants undergoing magnetic resonance imaging (MRI) were excluded if they presented with a history of psychosis, a history of stroke, a history of serious head trauma as defined by loss of consciousness accompanied by confusion, slurred speech, or amnesia, or other neurological disorders such as epilepsy. We also excluded participants with a history of brain cancer, chronic autoimmune diseases, or heart failure, and participants currently in renal failure or receiving dialysis, those who had a myocardial infarction in the past year, and/or an indication of unmanaged diabetes, and those who had evidence of severe liver disease or hepatitis. Participants receiving cognitively active medications, or anti-inflammatory, or immunomodulatory drugs were also excluded. Participants with claustrophobia, embedded ferromagnetic metal implants, pacemakers, shrapnel, wires, and/or other MRI-unsafe surgical implants were excluded.
MRI Acquisition
All images were acquired on a 3T Siemens Biograph mMR scanner (V.VE11P) using the vendor-provided 20-channel head/neck coil. Diffusion MRI images were acquired using a state-of-the-art multi-band diffusion-weighted imaging sequence 28, 29 (obtained via C2P, Center for Magnetic Resonance Research, University of Minnesota) with TE/TR = 121.4/6300 ms, voxel size = 2x2x2 mm3, multi-band factor = 3. A multi-shell diffusion scheme was used with multiple b-values = 1000, 2000, 3000, and 4000 s/mm² with diffusion sampling directions 64, 32, 32, and 32, respectively. T1-weighted magnetization-prepared rapid gradient echo (MPRAGE) images were also acquired (TE/TR/TI = 2.49/1900/900 ms, flip angle = 9°, isotropic voxel size = 0.9x0.9x0.9 mm3, grappa factor = 2, Turbo factor = 192).
MRI image processing
Brain parcellation was performed with FreeSurfer (7.3.0) using T1-MPRAGE images. Two subjects were removed from the analysis because, upon visual inspection, there was evidence of chronic cerebral infarctions. The diffusion MRI images were processed using DSI Studio (V.06142023). FSL’s eddy was used to correct for eddy current distortion. We calculated fractional anisotropy (FA), mean diffusivity (MD), axial diffusivity (AD), and radial diffusivity (RD) using diffusion tensor imaging with diffusion images acquired with b-value = 0; 1,000 s/mm². To generate the spin distribution function 30, Q-space diffeomorphic reconstruction 31 was performed in the MNI (Montreal Neurological Institute) space with a diffusion sampling length ratio of 1.25. Output resolution in diffeomorphic reconstruction was an isotropic 2 mm. The model-free technique separates isotropic from anisotropic diffusion, thereby reducing partial volume effects and crossing fibers. This approach has been shown to have improved reliability against edema32 and crossing fibers33. It allows the quantification of isotropic diffusion (ISO, isotropic value of the spin distribution function). Quantitative anisotropy (QA) was extracted as a local connectome fingerprint 34, and used in connectometry analyses. Whole-brain average white matter diffusion parameters were calculated using masks generated by combining left and right cerebral white matter in the FreeSurfer DKT atlas. Cortical thickness and cerebral volume were extracted from the T1MPRAGE images using the standard pipeline in Freesurfer. All image analyses were performed on a Mac mini with Apple Silicon M2 Pro (32GB memory) running Ventura 13.5.
Statistical Analyses
For gross comparisons between group characteristics, the Fisher exact test was utilized for categorical variables, while Welch’s t-test was utilized for continuous variables. To report effect sizes after adjustment for covariables we used generalized linear modeling and estimated standardized beta coefficients (b). All modeling assumptions were tested, and models were adjusted for potential unmatched confounders. When examining symptoms-based correlations within COVID-19-infected individuals, models adjusted for age, and intracranial volume. When examining cognitive outcomes, we additionally adjusted for schooling and estimated premorbid cognitive ability. Heat maps were generated, and estimated effect sizes were overlaid to aid in interpretation.
Correlational tractography was used to examine white matter tractography in COVID-19-infected versus uninfected participants and to compare those with acute COVID-19 to those with PASC. We used a nonparametric Spearman partial correlation to derive the correlation because of the small sample size and the potential for non-Gaussian distributions. As an effort to show differences in effect size, we reported images with a length threshold of 15 voxels (where T = 1), that were tracked using a deterministic fiber tracking algorithm with whole brain seeding. Generated tracts were then filtered by topology-informed pruning with 16 iterations. We use the area under the receiver-operating curve (AUC) to indicate the extent to which neuroimaging measures differed between diagnostic categories.
Since cerebral atrophy can indicate increased severity and poorer prognoses in neurological conditions, we examined correlations between gray and white matter parameters to determine if any white matter changes were related to neurodegenerative differences. Volumetric estimates were adjusted for total intracranial volume in all models.
Statistical significance was determined using a two-tailed p-value (a = 0.05) and adjusted for the false discovery rate (FDR) 35. FDR in tractography analyses was estimated using a total of 4,000 randomized permutations applied to the case label to obtain the null distribution of the tract length; tracts significant after adjusting for multiple comparisons (FDR = 0.05) were reported. Group-wise analyses were performed using Stata MP/17 (StataCorp), and Tractography was analyzed using DSI-Studio (dsi-studio.labsolver.org).
Sensitivity Analyses
Some factors might influence vulnerability to other cerebrovascular or neurodegenerative diseases, so we examined the potential for diffusion and volumetric parameters to relate to vaccination status or APOE4 status.
Ethics
The [Institution] internal review board reviewed all study procedures. Participants provided informed written consent to participate in this study.