ABSTRACT
Purpose Functional magnetic resonance imaging (fMRI) and functional connectivity (FC) have been used to follow aging in both children and older adults. Robust changes have been observed in children, where high connectivity among all brain regions changes to a more modular structure with maturation. In this work, we examine changes in FC in older adults after two years of aging in the UK Biobank longitudinal cohort.
Approach We process data using the Power264 atlas, then test whether FC changes in the 2,722-subject longitudinal cohort are statistically significant using a Bonferroni-corrected t-test. We also compare the ability of Power264 and UKB-provided, ICA-based FC to determine which of a longitudinal scan pair is older.
Results We find a 6.8% average increase in SMT-VIS connectivity from younger to older scan (from ρ = 0.39 to ρ = 0.42) that occurs in male, female, older subject (> 65 years old), and younger subject (< 55 years old) groups. Among all inter-network connections, this average SMT-VIS connectivity is the best predictor of relative scan age, accurately predicting which scan is older 57% of the time. Using the full FC and a training set of 2,000 subjects, one is able to predict which scan is older 82.5% of the time using either the full Power264 FC or the UKB-provided ICA-based FC.
Conclusions We conclude that SMT-VIS connectivity increases in the longitudinal cohort, while resting state FC increases generally with age in the cross-sectional cohort. However, we consider the possibility of a change in resting state scanner task between UKB longitudinal data acquisitions.
1. INTRODUCTION
Functional magnetic resonance imaging (fMRI) is a non-invasive technique that has proven indispensable for investigating human neural processes in vivo.1 For example, it has been used to localize the areas associated with vision,2 attention,3,4 emotion,5,6,7 and language8 to specific regions in the cortex, or at least find the regions that are most significantly involved in a specific task. Functional connectivity (FC) is a quantity derived from fMRI that measures the time correlation of blood oxygen level-dependent (BOLD) signal between different regions in the brain.9 FC has recently been used to predict age,10,11 sex,12,13 race,14 psychiatric disease status,15,16 and preclinical Alzheimer’s disease.17 Efforts to predict general fluid intelligence, although common,18,13 are thought by some to be confounded by differential achievement score distribution among ethnicities and the strong presence of race signal in FC.14 FC has proven effective in predictive studies because of its simplicity and its robust representation of complex BOLD signal activity, as evidenced by high subject identifiability across different scanner tasks and across time.19,20,21
Besides being used as a predictive tool, FC has been observed to undergo changes throughout the lifespan. For example, connectivity in young children is generally very high between all brain regions and decreases while also becoming more modularized during and after puberty.22 The FC of males and females is also quantitatively different, with females having higher intra-DMN connectivity and males having relatively greater connectivity between the DMN and other networks, although there is a wide degree of individual variation.23,24 Meanwhile, studies have shown that changes occur in the DMN during late middle and old age,25 although the exact direction of change in FC does not always appear constant.26 In addition, various studies have examined age-related changes in the cingulum27 and medial temporal lobe.28 Given the recent interest in using fMRI to predict pre-clinical Alzheimer’s disease,26,17 we believe a knowledge of ordinary changes in FC during old age is essential. This is especially true because it has been shown that a confounder can easily be mistaken for a true signal indicative of, e.g., general fluid intelligence or achievement score.14
This study uses the longitudinal cohort of the UKB29 to examine changes in the FC of individuals after an average of two years, the time between longitudinal scans. The UKB population of subjects with fMRI scans is predominantly (98%) Caucasian, ruling out race as a possible confounding effect. Additionally, we investigate changes in FC in longitudinal sub-populations based on subject age and sex. We find that average FC between SMT-VIS networks increases on average from the first scan to the second, and that SMT and VIS-related connectivities are more predictive of scan age than those of other networks. The complete FC, or a large subset, is still required to attain the best accuracy.
2. METHODS
We first describe the UKB dataset and the longitudinal subset used for our analysis. We then describe pre-processing of the fMRI data and conversion into FC. Finally, we discuss prediction of older vs younger scan in the longitudinal cohort and detail our methods for analysis of FC changes.
2.1 UK Biobank Longitudinal Cohort
The UKB contains various data of more than 500,000 subjects in the UK, of who more than 40,000 have fMRI scans.29 We processed two longitudinal resting state scans for 2,722 subjects, taken approximately two years apart. These subjects are approximately equally split between male and female, and have significant numbers of younger and older adults. The longitudinal cohort is composed of 1,289 genetic males and 1,369 genetic females, with the rest not having genetic sex information. The ethnicity of the subset of the UKB with fMRI scans is 98% Caucasian. Besides the 2,722 subjects we processed, an additional 154 subjects have the second longitudinal scan but not the first, resulting either from missing original source data or a failure in our SPM12-based preprocessing pipeline.
2.2 fMRI Preprocessing
The original scan acquisition parameters are described elsewhere,30,31 but consist of both resting state and task fMRI scans with a repetition time of TR = 0.735 sec. For this study, we examined the resting state scans only. All resting state 4D fMRI volumes were processed with SPM12, including co-registration and warping to MNI space (http://www.fil.ion.ucl.ac.uk/spm/software/spm12/). BOLD signal was extracted using the Power264 atlas,32 which consists of 264 ROIs grouped into 14 functional networks and represented by 5mm radius spheres. The resulting timeseries were bandpass filtered between 0.01 and 0.15 Hz to remove scanner drift, noise, heartbeat, and some breathing signal. Pearson correlation of the filtered timeseries was used to create subject-specific FC matrices, which were reduced to the unique entries in the upper right triangle and vectorized. The entire procedure is summarized in Figure 1.
In contrast to the Power264 atlas-derived FC constructed by us, the original UKB data provided the unique part of 21-region and 55-region FC and partial correlation-based connectivity (PC) matrices based on ICA in vectorized format.30 These matrices were calculated through the use of PCA on whole cohort fMRI data followed by ICA,30 meaning that regions overlap in an unpredictable way and are not associated with specific functional networks. Although prediction using 55-component ICA-based FC and PC is often as good as and sometimes better than prediction using Power264 atlas-derived FC, the resulting connectivities are uninterpretable with regards to BOLD signal within specific regions. Additionally, in predicting which scan is older, Power264 asymptotes to a higher predictive accuracy than either of the ICA-derived measures (see Figure 4).
2.3 Prediction of Scan Order and Analysis of FC
Prediction of scan age in the UKB longitudinal cohort was carried out by logistic regression (https://scikit-learn. org/stable/) models with 20 bootstrapping repetitions, using the scikit-learn implementation.34 The regularization parameter was fixed to C = 1, which was found to be near the optimal value for all training set sizes using grid search. It was found that a simple difference of scan FCs gave the best prediction results compared to concatenation or difference and concatenation, using either logistic regression or MLP. The training set was created with randomization of whether older scan was subtracted from younger scan or younger scan was subtracted from older scan. Our code for computing prediction accuracy can be found online (https://github.com/aorliche/ukb-longitudinal-smt-vis). However, UKB data sharing policy precludes us from posting the longitudinal data itself; interested researchers may contact us with any questions.
Analysis of FC was performed by finding the mean (Figure 2) and standard deviation (Figure 9) of older scan FC minus younger scan FC for the longitudinal cohort. Additionally, prediction of scan order was carried out using the average connectivity between each of the Power264 networks, each network consisting of many individual ROIs. As before, logistic regression with 20 bootstrap repetitions and C = 1 was used for this purpose. A Bonferroni-corrected two-sided t-test was applied to the 105 average inter-network connectivity differences (from the complete graph of 14 functional networks) of the 2,722 longitudinal subjects to determine if they were significantly different from zero (Figure 5 Bottom).
3. RESULTS
We first describe trends in FC changes during the average of 2 years between longitudinal scans, summarize the ability of simple machine learning models to identify older vs younger scan, and investigate the ability of specific inter-network connectivities to predict scan order. We then summarize the statistical significance of inter-network FC changes with aging, in both the longitudinal and cross-sectional cohorts of the UKB. Finally, we consider the possibility that the observed longitudinal changes are due to a change in scanner task by presenting inter-task FC differences in the Philadelphia Neurodevelopmental Cohort (PNC) dataset.35
3.1 Inter-Network FC Changes
In Figure 2, we show that, on average, SMT-VIS connectivity increases from younger scan to older scan. The right hand side of Figure 2 displays divisions of the 14 functional networks included in the Power264 atlas. Network labels and abbreviations are listed in Table 1. The increase in connectivity is large and distinct over the majority of SMT-VIS FCs compared to other non-SMT and non-VIS FCs. Many FCs involving the VIS network appear to increase in connectivity from the first scan to the second. The average change in FC in the SMT-VIS connection is 6.8%, corresponding to a mean change µΔρ = +0.03, compared to a standard deviation of σΔρ = 0.26. Figure 5 shows that, although small compared to the standard deviation, this difference is very significant.
Figure 3 displays the same analysis, i.e., the average change from first scan to second, for four subsets of the cohort. These subsets are male subjects, female subjects, young (< 55 years old) subjects, and old (> 65 years old) subjects. All four subsets observed the same effect as the whole cohort, thus we rule out very old age or gender as confounding factors.
3.2 Predicting Older Scan of Pair
In Figure 4, one can see that is possible to predict which scan of a longitudinal pair is older with the Power264 atlas at an accuracy of 82.5%, having 2,000 subjects in the training set and the rest in the test set. This measurement was repeated with 20 bootstrap iterations and averaged. The entire 34,716-feature upper right triangle of the FC matrix was used ot make the prediction. One can also see that the ICA FC/PC matrices provided pre-processed along with the UKB data are also able to predict scan order, although at a slightly reduced accuracy. Prediction is possible at an accuracy of 60-70% using only 100-200 training set subjects.
3.3 Prediction of Older Scan Using Specific Inter-Network Connections
In Figure 5, we rank average inter-network FCs in their ability to predict scan order. As expected from the mean change in FC (Figure 2), the SMT-VIS connection is the most predictive of longitudinal scan age. Furthermore, SMT and VIS networks are included among the next several most predictive inter-network connections. In Figure 5 bottom, we plot the predictive ability of all 105 inter-network connections, along with a p-value for the inter-network FC change being significantly different from zero. The raw p-value has been multiplied by 105 to account for multiple comparisons. It is highly significant for the first 10 or so most predictive inter-network connections.
Table 3.3 lists the number of subjects whose FC increased or decreased for the SMT-VIS connection and over the entire brain. The table is divided among the four subsets of the longitudinal cohort mentioned previously. Additionally, we correlated several dozen subject phenotypes and longitudinally-tracked variables with changes in FC and report the most significant in Appendix B. In that section, we find an interesting but small correlation with hand grip strength, body mass index (BMI), and basal metabolic rate. In Section 3.4, we find that average resting state FC increases with age across most inter-network connections in the much larger UKB cross-sectional cohort.
3.4 FC Changes with Age in the UKB Cross-Sectional Cohort
We find that average resting state FC has a significant increase in almost all inter-network connections in the UKB cross-sectional cohort. Average maps of FC change are shown in Figure 7. We fail to find a higher SMT-VIS change compared to other connections; however, almost all inter-network regions have a large positive increase in FC with aging. We give precise numbers for four inter-network connections as well as total FC in Table 3.
In total, there are 9,387 older males (> 65 years old), 2,425 younger males (< 55 years old), 8,728 older females (> 65 years old), and 3,132 younger females (< 55 years old) in the UKB cross-sectional cohort.
3.5 Comparison with FC Differences Between Scanner Tasks in the PNC Dataset
We consider the possibility that the difference in SMT-VIS connectivity between the two scans of the longitudinal cohort is due to a change in scanner task. In Figure 8, we show the average FC differences between 3 different tasks in the PNC dataset.35 This dataset contains 1,345 children and young adults having all of three different scanner tasks: resting state, working memory, and emotion identification. The preprocessing and FC creation steps for this dataset have been described elsewhere.10 We note that the VIS-VIS is most different for change in task, but that the SMT-VIS is not qualitatively more different that the rest of FC. Also, the magnitude of change in FC in the PNC dataset between tasks is much larger than in the UKB longitudinal cohort.
4. DISCUSSION
Farràs-Permanyer et al. (2019)36 find that mean resting state FC may increase throughout the entire brain for the oldest subject (> 80 years old) group. As shown in Appendix A, we confirm a small, statistically insignificant increase in total longitudinal FC in the healthy controls of the ADNI dataset,37 another elderly population with multiple longitudinal fMRI scans. In Section 3.4, we show that there is a large, statistically significant increase in average resting state FC across almost all inter-network connections in the UKB cross-sectional cohort with increased age. This cross-sectional cohort is much larger than the longitudinal cohort we describe in the main part of this paper. The fact that SMT-VIS FC also increases in the cross-sectional cohort, but not disproportionately compared to the rest of FC, raises the possibility of a change in resting state scanner task during the second longitudinal scan. We believe this longitudinal change is not an artifact of our pre-processing methods. Credence should increase in our pre-processing methods since the UKB-provided ICA-based FC and PC is also able to predict longitudinal scan ordering at almost the same level as our Power264-based approach, although the ICA FC and PC matrices are not interpretable.
Many studies have focused on examining connectivity in the DMN associated with aging.38,39 These studies find areas of increased connectivity as well as areas of decreased connectivity. There are two problems with such studies. First, they are for the most part cross-sectional and do not follow a single subject across a multi-year period. Second, they mostly use small numbers of subjects, the majority of studies enrolling fewer than 50, making it impossible to identify small effects. On the other hand, one study performed on a cohort of more than 2,000 older subjects in Rotterdam found age-related changes in connectivity to be complicated, drawing no firm conclusions.40 We note that the Rotterdam study was not longitudinal but cross-sectional.
We conjecture the fact that most studies only focus on DMN and report decreased connectivity38 in aging populations may be related to the large number of ROIs in the DMN and an implicit bias inherent in the word “connectivity.” Naturally, as we reach very old age we expect physical connections to degenerate, not become stronger. In fact, FC is really the synchronization of BOLD signal between regions, and does not imply a direct physical link between regions. Young children are known to have higher average FC than young adults;22,41 thus older subjects may been as reverting to a less optimal state as they age.
On the other hand, as we describe in Appendix B, physical observables such as hand grip strength in the UKB longitudinal cohort are weakly correlated with an increase in FC in SMT-CB and VIS-CB connectivity. Additionally, we find BMI and basal metabolic rate are weakly correlated with the longitudinal increase in SMT-VIS connectivity (see Appendix B). This may suggest a small health related effect that is found throughout the study cohort and includes male, female, younger, and older subjects. Finally, as discussed in Section 3.5, we cannot rule out the possibility that changes in longitudinal FC are caused by a change in scanner task between acquisition of the two longitudinal timepoints.
We show in this work that the average connectivity increase in the SMT-VIS connection is small but highly statistically significant. The average change in FC in this connection is only 6.8%, corresponding to a mean change µΔρ = +0.03, compared to a standard deviation of change from subject to subject of σΔρ = 0.26 (see Figure 9). However, using our longitudinal sample of 2,722 subjects, we find the average SMT-VIS connectivity change from younger scan to older scan is significant as level of p < 10−15 after Bonferroni correction for multiple comparisons (Figure 5). Finding such small effects is helped by the use of large number of subjects and longitudinal data.
5. CONCLUSION
In this work, we pre-process a 2,722 subject longitudinal subset of the UK Biobank dataset and examine FC using the Power264 atlas. We find that in scans taken an average of two years apart, the average functional connectivity between SMT and VIS network regions tends to increase. This occurs in male, female, younger (< 55 years old), and older (> 65 years old) subjects. We verify the ability of this average FC increase to predict scan ordering using simple machine learning models. The identification of an increase in connectivity with non-pathological aging, in longitudinal as well as cross-sectional cohorts, and specifically in the SMT-VIS synchronization of BOLD signal, may lead to novel insights about brain function in old age. Additionally, we identify an effect that could possibly show up as a confounder in studies of dementia or neurodegenerative diseases. Nonetheless, we remain open to the idea of a change in resting state scanner task during acquisition of the longitudinal data in the UKB being partly responsible for this effect.
DISCLOSURES
The authors have no conflicts of interest to report.
CODE AND DATA AVAILABILITY
fMRI and phenotype data came from the UK Biobank (application ID 61915), available via application to qualified researchers. Additional data used in Appendix A came from the Alzheimer’s Disease Network Initiative (ADNI), available via application from https://adni.loni.usc.edu/. Additional neuroimaging data came from the Neurodevelopmental Genomics: Trajectories of Complex Phenotypes database of genotypes and phenotypes repository, dbGaP Study Accession ID phs000607.v3.p2.
All code used in this study is available from GitHub at https://github.com/aorliche/ukb-longitudinal-smt-vis. We do not have permission to post original subject data; however, it may be obtained via application from the sources listed above.
Biography
Anton Orlichenko received his B.S. in Electrical and Computer Engineering in 2010 from the Illinois Institute of Technology. He is currently a Ph.D. candidate at Tulane University, working on mathematical models of brain function and genetics.
Kuan-Jui Su received his Ph.D. from Tulane University School of Medicine in 2023. He currently works at the Center for Biomedical Informatics and Genomics at Tulane University.
Qing Tian holds an M.S. degree and works with Dr. Deng and Dr. Shen as the database administrator for the Center for Biomedical Informatics and Genomics at Tulane University.
Hui Shen received his B.S. from the University of Science and Technology in China and his Ph.D. from Creighton University. He is currently a geneticist at the Center for Biomedical Informatics and Genomics at Tulane University.
Hong-Wen Deng received his B.S. and M.S. from Peking University. He received a further M.S. and Ph.D. degrees from the University of Oregon. He currently runs the Center for Biomedical Informatics and Genomics at the Tulane University School of Medicine.
Yu-Ping Wang received his B.S. from Tianjin University in 1990, and his M.S. and PhD degrees from Xian Jiaotong University in 1993 and 1996, respectively. He currently runs the Multiscale Bioimaging and Bioinformatics Lab at Tulane University.
ACKNOWLEDGMENTS
The authors would like acknowledge the NIH (grants R01 GM109068, R01 MH104680, R01 MH107354, P20 GM103472, R01 EB020407, R01 EB006841, R56 MH124925) and NSF (#1539067) for partial funding support.
This research was supported in part using high performance computing (HPC) resources and services provided by Information Technology at Tulane University, New Orleans, LA.
APPENDIX A. LONGITUDINAL FC CHANGES IN THE ADNI DATASET
We examined the longitudinal change in FC of healthy controls in the Alzheimer’s Disease Neuroimaging Initiative (ADNI) dataset37 (age matched subjects who do not develop AD pathology). We used scans taken an average of one year apart. We confirm a small, statistically insignificant increase in total FC but fail to find the same SMT-VIS increase relative to the rest of FC as in the UKB. Statistics are given in Table 4 and the average FC change is shown in Figure 10.
APPENDIX B. correlation OF CHANGE IN FC WITH LONGITUDINAL OUTCOMES IN THE UKB
We identified several correlations between longitudinal change in FC and changes in clinical outcomes associated with the two scan timepoints in the UKB dataset. These are presented below, along with the UKB field identifiers of the outcomes. All p-values are Bonferroni-corrected with n = 105 multiple comparisons (one for each average inter-network connectivity).
B.1 SMT Hand, VIS, and CB Connectivity and Grip Strength (f.46.2.0, f.46.3.0, f.47.2.0, f.47.3.0)
We find a marginally significant association between change in hand grip strength and VIS-CB and SMT-CB connectivity change (Figure 11).
B.2 Body Mass Index and Basal Metabolic Rate (f.23104.2.0, f.23104.3.0, f.23105.2.0, f.23105.3.0)
We find a not statistically significant but suggestive association between BMI and basal metabolic rate change and SMT-VIS connectivity change (Figure 12).
Footnotes
Manuscript now includes analysis of PNC dataset demonstrating change in functional connectivity (FC) with change in scanner task (Section 3.5). The manuscript also moves data about overall increase in FC with aging into the main body of the paper (Section 3.4). The abstract is slightly revised.