Abstract
Background Our behavioral traits, and subsequent actions, could affect the risk of exposure to the coronavirus disease of 2019 (COVID-19). The current study aimed to determine whether unique brain endophenotypes predict the COVID-19 infection risk.
Methods This research was conducted using the UK Biobank Resource. Functional magnetic resonance imaging scans acquired before the COVID-19 pandemic in a cohort of general population older adults (n=3,662) were used to compute the whole-brain functional connectomes. A network-informed machine learning approach was used to identify connectome and nodal fingerprints that predicted positive COVID-19 status during the pandemic up to February 4th, 2021.
Results Brain scans, acquired an average of 3 years before COVID-19 testing, significantly predicted the infection results. The predictive models successfully identified 6 fingerprints that were associated with COVID-19 positive, compared to negative status (all p values < 0.005). Overall, lower integration across the brain modules and increased segregation, as reflected by internal within module connectivity, were associated with higher infection rates. More specifically, COVID-19 infections were predicted by 1) reduced connectivity between the central executive and ventral salience, as well as between the dorsal salience and default mode networks; 2) increased internal connectivity within the default mode, ventral salience, subcortical and sensorimotor networks; and 3) increased connectivity between the ventral salience, subcortical and sensorimotor networks.
Conclusion Individuals are at increased risk of COVID-19 infections if their brain connectome is consistent with reduced connectivity in the top-down attention and executive networks, along with increased internal connectivity in the introspective and instinctive networks.
Introduction
“What’s natural is the microbe. All the rest—health, integrity, purity (if you like)—is a product of the human will, of a vigilance that must never falter. The good man, the man who infects hardly anyone, is the man who has the fewest lapses of attention.” – Albert Camus, The Plague
Chronic stress pathology has been increasingly recognized as a major factor in the pathophysiology of neuropsychiatric disorders 1. For certain mental illnesses—for example posttraumatic stress and major depression—trauma and stress could be the triggering and/or the perpetuating factors. While for others, such as anxiety and paranoia, chronic stress is a detrimental outcome that may exacerbate the underlying pathology 2,3. Furthermore, biological correlates of stress-related disorders may reflect predisposing components and/or an outcome of chronic stress pathology. For example, reduced hippocampal volume is believed to be both a predisposing factor as well as an outcome of posttraumatic stress disorder (PTSD) 4-6. Hence, although it is important to determine the brain correlates of chronic stress, it is critical to disentangle the predisposing markers from the consequences of stress. Capitalizing on the fortuitous large neuroimaging dataset from the UK Biobank 7, the aim of the current report is to identify the brain signatures that predisposed the general population to a major worldwide stressor, the coronavirus disease of 2019 (COVID-19).
In early 2020, COVID-19 cases spreading globally instigated a devastating pandemic 8. To date, more than 115,000,000 cases of COIVD-19 were identified, and more than 2,500,000 deaths were related to COVID-19. From life-threatening hospitalizations and the loss of loved ones, to lockdowns, isolation, and increased unemployment and domestic conflict, the impact of the pandemic has been overwhelming 8,9. Moreover, COVID-19 can cause direct damage to the brain through encephalopathy 10. Thus, in the coming years, it will be essential for the field to assess the long-term impact of the pandemic on mental health and brain function. Equally important is the need to determine the brain signatures that predate the pandemic but correlate with higher COVID-19 infection rate. Identifying these biomarkers will help us disentangle the sequelae of COVID-19 from its predisposing brain markers. It will also provide a greater understanding of the brain functions role in the spread of the COVID-19 infections, which may assist in developing future preventive strategies.
This report will focus on the role of the brain intrinsic connectivity networks, using functional connectome fingerprinting. This machine-learning approach allows full assessment of the brain connectome, while providing network informed results 11,12. It is a combination of the network-restricted strength (NRS) 13 and connectome-based predictive modeling approaches 14. Functional connectome fingerprints (CFPs) were reported to predict behavior in the general population 15,16 and treatment response in depressed patients 12,17,18. This NRS predictive model (PM) approach has several major strengths. First, predictive features can be back-translated to the original space, which is not often the case in machine-learning algorithms. Thus, instead of establishing a “black box” computational algorithm that is predictive of the outcome but is undiscernible, the NRS-PM works to identify the brain biomarker that is associated, significantly and consistently, with the outcome of interest (e.g., infection status) regardless of the intensity of prediction. Second, the NRS-PM approach could enhance reproducibility by providing protection against overfitting, which is an issue with traditional interpretive statistics. Third, the multivariate pattern analysis permits the full assessment of the connectome, without the inherent increase of Type I error due to univariate multiple comparisons or the need to restrict the analysis to a limited selection of seeds and targets. Finally, the NRS-PM results are network-based by design, which both informs the neurobiological models and facilitates the integration of findings 17,19.
Based on the Akiki-Abdallah (AA) hierarchical connectivity atlas 11,12, the brain connectome is divided into 7 canonical networks: 1) central executive (CE); 2) default mode (DM); 3) ventral salience (VS); 4) dorsal salience (DS); 5) subcortical (SC); 6) sensorimotor (SM); and 7) visual (VI) 11, 12. In an environment with multiple priorities and stimuli competing for our attention, two brain systems, the DS and VS, dictate which stimuli is deserving of our attention 20. The DS, sometimes called the dorsal attention network, is involved in top-down voluntary attention to salient stimuli. In contrast, the VS network is primarily responsible for reorienting brain resources in response to involuntary salient (i.e., important or conspicuous) external and internal stimuli 21. These brain systems interact with the DM and CE networks, which are responsible for internally and externally directed cognitions, respectively 22. While the function of the SC network remains unknown, it is hierarchically derived from the salience system and was previously found to complement connectivity changes in the CE 12,17. The current study conducted a data-driven approach assessing all whole brain networks. However, considering the hypothesized role of the brain networks 23, it is conceivable to anticipate increased COVID-19 infections in individuals with reduced connectivity in the top-down, attention and executive, control networks (i.e., DS and CE).
Methods
Data used in this study were extracted from the UK Biobank data repository under application number 42826. All study procedures were approved by Institutional Review Boards and all participants completed an informed consent process.
Participants
The UK Biobank is a prospective epidemiological study of approximately 500,000 participants. Details of the UK Biobank resource and procedures can be found online (https://www.ukbiobank.ac.uk) and in previous reports 24. Briefly, between 2006 and 2010, community-dwelling general population individuals across the United Kingdom (n = 502,536; 40 to 69 years of age at the time of recruitment) provided extensive genetic, physical, and health data 24. In 2016, a followup imaging study was funded to scan 100,000 participants from the existing cohort (including brain, abdomen, heart, and whole-body scans) 7.
Imaging Data
This study used three UK Biobank brain imaging modalities acquired on Siemens Skyra 3T magnet (see 7,25 for more details), including structural high-resolution MRI (T1; 1×1×1 mm), resting-state fMRI (2.4×2.4×2.4 mm; 490 frames in 6min.), and task fMRI (2.4×2.4×2.4 mm; 332 frames in 4 min.; Hariri faces/shapes “emotion” task 26). We used the UK Biobank preprocessed NIfTI files (see 7 for more details). Briefly, only “usable” data (i.e., following manual review and auto quality checks 25) were used. For all modalities: quality check scores were generated, based on alignments and signal-to-noise ratios; gradient distortion correction was applied; and nonlinear transformations between native and standard spaces were generated. B0 fieldmaps were used to correct EPI distortion for fMRI. In addition, structural MRI preprocessing included tissue-type segmentation using FAST (FMRIB’s Automated Segmentation Tool 27) and subcortical structure modeling using FIRST (FMRIB’s Integrated Registration and Segmentation Tool 28). Preprocessing of fMRI scans included: motion correction, grand-mean intensity normalization, high-pass temporal filtering (sigma=50s), and structured artefact removal by ICA+FIX processing (Independent Component Analysis followed by FMRIB’s ICA-based X-noiseifier 29).
Connectome and Nodal Predive Models
Full details of the network restricted strength predictive model (NRS-PM) methods were previously reported 12,13,17. Briefly, individual specific FAST and FIRST segmentations were used to extract the average time series of 424 nodes that cover the whole-brain gray matter based on the A424 atlas 12,30-32. The Akiki-Abdallah hierarchical connectivity at 50 modules (AA-50; Fig. S1), 24 modules (AA-24; Fig. S2), and 7 modules (AA-7) were used to determine the network affiliation of the A424 nodes (https://github.com/emergelab). The full connectome is the Fisher-Z transformation of the pairwise correlation coefficients. NRS connectome is the pairwise average connectivity of all modules at AA-50, AA-24 and AA-7 12. Nodal strength (nS) is the average connectivity of a node to all other nodes. Nodal internal NRS (niNRS) is the average connectivity between each node and all other nodes within the same canonical connectivity network (i.e., AA-7). Nodal external NRS (neNRS) is the average connectivity between each node and all other nodes outside its canonical connectivity network 12. The predictive models used were adapted from the connectome-based predictive model approach 14, as previously detailed 12. All NRS-PM functions used in the current study are publicly available at https://github.com/emergelab. The modeling includes feature selection in training subsamples, followed by fitting a linear predictive model, then applying the model to the test subsample 14. Finally, 200 iterations of ten-fold cross-validation (CV) were conducted to ensure the stability of the models and to determine the statistical significance; that is by comparing true and random predictions 12. The predictive model included both resting and task fMRI connectome data to improve the study predictions 33.
Statistical analyses
Descriptive statistics were calculated prior to statistical analysis. Data distributions were checked using normal probability plots. The statistical significance threshold was set at 0.05 (2-tailed tests). MATLAB (2018a; Mathworks Inc.) and the Statistical Package for the Social Sciences (version 24; IBM) software were used for the analyses. False Discovery Rate (FDR; q < 0.05) was used to correct for multiple comparisons. The connectivity fingerprints (CFPs) were examined at AA-50, AA-24, and AA-7. The nodal fingerprints (NFPs) were determined for nS, niNRS, and neNRS. FDR was applied on all 6 outcome measures to determine statistical significance.
As in previous reports 12,17, connectivity per fingerprint was computed by multiplying the connectivity features (e.g., NRS at AA-50) by the corresponding weighted fingerprint masks (e.g., the COVID-19 CFP at AA-50). Thus, the CFP total connectivity is the sum of weighted estimates per subject per CFP. To facilitate the comparison across measures, the CFP connectivity values were standardized (z-scored). Follow-up analyses covarying for age and sex used general linear models with the 6 fingerprints’ total connectivity as dependent variables and COVID-19 status as fixed factor. The study atlases, code, and predictive models will be made publicly available at https://github.com/emergelab.
Results
The brain imaging data were based on a package downloaded on July 8th, 2020. At the time the data were downloaded, preprocessed brain imaging data from 40,681 participants were available for this report. The COVID-19 results were downloaded on February 4th, 2021. The COVID-19 data included 60,446 UK Biobank participants, of which 3,662 had successful structural MRI, and resting and task fMRI scans. These 3,662 individuals were investigated in the current report. They were 52% females (n=1896). Their average age was 63 years (SEM=0.13) at the time of the brain scan and 66 years (SEM=0.13) at the time of the COVID-19 testing. A total of 921 (25%) tested positive for COVID-19. All imaging data used in the current study were acquired prior to the COVID-19 pandemic.
The predictive models successfully identified 6 fingerprints that were significantly associated with COVID-19 positive, compared to negative status: (1) AA-50 CFP (r = 0.13, CV = 10, iterations = 200, p < 0.005, q < 0.05; Fig. 1A); (2) AA-24 CFP (r = 0.13, CV = 10, iterations = 200, p < 0.005, q < 0.05; Fig. 1B); (3) AA-7 CFP (r = 0.11, CV = 10, iterations = 200, p < 0.005, q < 0.05; Fig. 1C); (4) nS NFP (r = 0.09, CV = 10, iterations = 200, p < 0.005, q < 0.05; Fig. 2A-B); (5) neNRS NFP (r = 0.10, CV = 10, iterations = 200, p < 0.005, q < 0.05, Fig. 2C); and (6) niNRS NFP (r = 0.12, CV = 10, iterations = 200, p < 0.005, q < 0.05; Fig. 2D).
As shown in Fig. 1, positive COVID-19 tests were predicted by increased internal connectivity within modules but reduced external connectivity between the brain networks. In particular, positive COVID-19 tests were associated with reduced connections between the central executive (CE) and ventral salience (VS), as well as between the dorsal salience (DS) and default mode (DM) modules. In contrast, increased interference from the VS to the sensorimotor (SM) and subcortical (SC) networks predicted positive COVID-19 results (Fig. 1C).
The external to internal connectivity shifts observed in the CFPs were translated into overall reduced neNRS but increased niNRS as shown in Fig. 2. Independent of network constraints, positive COVID-19 results were associated with increased overall functional connectivity strength (i.e., nS) in the insula and surrounding regions, as well as in the medial frontal area (Fig. 2B).
To account for the effect of age, a general linear model examined the effects of COVID-19 status on the 6 fingerprints’ total connectivity covarying for age. This multivariate test showed statistically significant effects of COVID-19 status on the fingerprints’ connectivity (F = 2.7, p = 0.01). Post-hoc univariate tests of between-subject effects were significant for nS (F = 14.3, p < 0.001), niNRS (F = 10.8, p = 0.001), and neNRS (F = 9.6, p = 0.002), but not for CFPs at AA-50 (F = 3.5, p = 0.060), AA-24 (F = 3.6, p = 0.059), and AA-7 (F = 1.0, p = 0.31). Covarying for sex did not affect the main study results with all 6 fingerprints retaining significance (all p values < 0.001).
Discussion
This report successfully identified pre-pandemic brain functional connectivity markers that significantly predict COVID-19 status in a relatively large general population sample of older adults. The results showed that individuals with positive COVID-19 status tended to have lower integration across the brain connectivity networks but increased internal connectivity within modules. Together, these findings indicate a shift toward increased segregation between the brain networks. In particular, COVID-19 infections were predicted by reduced connectivity between the top-down, attention (DS) and executive (CE) control networks and the introspective (DM) and instinctive (VS) networks, respectively. Importantly, the DS-DM and CE-VS reductions in top-down connectivity were accompanied by increased internal connectivity in both the DM and VS networks. This shift from external to internal connectivity was evidenced in the frontoparietal reduction in neNRS and increased niNRS in the insular and medial frontal brain regions. At the level of global connectivity as measured by nS, COVID-19 infections were predicted by increased connectivity in areas within the VS and DM networks.
Considering that the study scans were acquired an average of 3 years prior to the COVID-19 testing, the results suggest that these connectivity fingerprints may reflect personality traits of related phenotypes. Investigating the behavioral correlates of the COVID-19 CFPs and NFPs in future studies will be essential to better understand the role of the brain in mitigating infection risk and perhaps develop new approaches to limit the spread of infections in future epidemics.
Among the limitations of the current report is that the causal relationship between the network alterations and infection status cannot be established in this associative study. However, the longitudinal design rules out the possibility that the identified fingerprints are the consequence of COVID-19 infections. These brain biomarkers could be an underlying factor to increasing infections. They may also be a confound or compensatory response to an unknown underlying causal factor. While the COVID-19 NFPs retained significance when controlling for age and sex, the study failed to rule out the possibility that the variance in the COVID-19 CFPs may be at least partially affected by age as a confound. Another limitation is that the study cannot determine whether the results will generalize to younger adults. Future studies should further investigate the role of age in the COVID-19 brain connectivity signatures.
The study has several strengths, including: 1) large cohort; 2) longitudinal data; 3) high quality imaging acquisition and preprocessing using standardized methods; 4) data-driven approach that is not limited to a biased selection of seeds; 5) network informed design to facilitate the interpretation and guide future follow-up studies; and 6) urgently needed data to better understand a devastating pandemic and ultimately devise additional preventive measures to improve global health and reduce suffering. It is rare to have such a large neuroimaging sample collected prior to unpredictable traumatic stressors like disasters or epidemics. Not only does this allow future studies of post-pandemic neural alterations disentangle consequence of disease and chronic stress from predisposing characteristics, it also underscores the tremendous value of large, prospective, longitudinal, multi-modal epidemiological efforts like the UK Biobank, and the critical need to continue funding similar projects that cover the life-span and other geographic regions.
Data Availability
Data are publicly available through the UK Biobank.
Funding
This work was supported by the Department of Veterans Affairs, the VA National Center for PTSD and the Beth K and Stuart Yudofsky Chair in the Neuropsychiatry of Military Post Traumatic Stress Syndrome.
The funding sources have had no involvement in the study design, the analysis and interpretation of data, the writing of this report, or the decision to submit this article for publication.
Declaration of Conflicting Interests
Dr. Abdallah has served as a consultant, speaker and/or on advisory boards for Genentech, Janssen, Psilocybin Labs, Lundbeck, Guidepoint, and FSV7, and as editor of Chronic Stress for Sage Publications, Inc. He also filed a patent for using mTORC1 inhibitors to augment the effects of antidepressants (Aug 20, 2018).
Acknowledgement
This research has been conducted using the UK Biobank Resource under Application Number 42826. I would like to thank all the UK Biobank investigators, staff and participants for contributing to this tremendous effort. I also would like to thank the Emerge Research Program (emerge.care) staff and investigators, who made possible the line of research conducted in this study.