- School of Information and Communication Engineering, University of Electronic Science and Technology of China, Chengdu, China
Introduction: Alzheimer's disease (AD) is a neurodegenerative disease that significantly impacts the quality of life of patients and their families. Neuroimaging-driven brain age prediction has been proposed as a potential biomarker to detect mental disorders, such as AD, aiding in studying its effects on functional brain networks. Previous studies have shown that individuals with AD display impaired resting-state functional connections. However, most studies on brain age prediction have used structural magnetic resonance imaging (MRI), with limited studies based on resting-state functional MRI (rs-fMRI).
Methods: In this study, we applied a graph neural network (GNN) model on controls to predict brain ages using rs-fMRI in patients with AD. We compared the performance of the GNN model with traditional machine learning models. Finally, the post hoc model was also used to identify the critical brain regions in AD.
Results: The experimental results demonstrate that our GNN model can predict brain ages of normal controls using rs-fMRI data from the ADNI database. Moreover the differences between brain ages and chronological ages were more significant in AD patients than in normal controls. Our results also suggest that AD is associated with accelerated brain aging and that the GNN model based on resting-state functional connectivity is an effective tool for predicting brain age.
Discussion: Our study provides evidence that rs-fMRI is a promising modality for brain age prediction in AD research, and the GNN model proves to be effective in predicting brain age. Furthermore, the effects of the hippocampus, parahippocampal gyrus, and amygdala on brain age prediction are verified.
1. Introduction
Alzheimer's disease (AD) is a neurodegenerative disease commonly occurring in older adults (American Psychiatric Association and Association, 2013). People with AD usually first develop symptoms, such as mild memory degradation, followed by a continuous decline in cognitive function (Reitz et al., 2011), and are eventually diagnosed with AD after experiencing a brief phase of mild cognitive impairment (MCI). As the disease worsens, patients with reduced self-care ability become more dependent on family members for care, which seriously affects the patients' and their families quality of life. Therefore, the early diagnosis method for AD is of utmost importance. However, the existing diagnostic methods mainly rely on psychological tests and clinical observation for middle and advanced patients, lacking objectivity and effective early diagnosis methods. Therefore, it is essential to develop more accurate identification methods and find more objective biomarkers for early diagnosis, which will help patients receive treatment earlier and improve cure rates (Tahami Monfared et al., 2022; Warren and Moustafa, 2023). One such potential biomarker is brain age estimation, which could help detect mental disorders such as AD. Moreover, the neuroimaging-driven brain age prediction will contribute to studying the effect of AD on brain structures and functional networks. Meanwhile, recent studies have also highlighted that the significant delaying of the progression of MCI to AD will reduce both the prevalence and cost of AD (Anderson, 2019). To this end, it is crucial to identify individuals with accelerated brain aging by brain age prediction for precise identification and intervention of AD. Therefore, identifying those individuals at risk of developing AD early through accurate and reliable brain age prediction will significantly pave the way for developing effective preventive strategies (Villemagne et al., 2013).
With the development of artificial intelligence and medical imaging technology, studies focusing on brain age estimation have increased in recent years (Frizzell et al., 2022). Most of these studies have used structural magnetic resonance imaging (MRI) data (Gaser et al., 2013; Sajedi and Pardakhti, 2019; Bashyam et al., 2020; Levakov et al., 2020; Lee et al., 2022). For example, Bashyam et al. trained DeepBrainNet with two-dimensional images obtained by T1-weighted MRI with minimal preprocessing steps from 11,729 health control (HC) subjects. They achieved accurate age prediction and revealed that a moderately fitted model of brain aging was more suitable for distinguishing between AD and HC. This proposed method ensured the broad applicability of the model in clinical settings through straightforward preprocessing steps. However, a significant limitation is that it does not identify the anatomical regions affecting the model's performance. Lee et al. (2022) suggested that occlusion sensitivity analysis enhanced the interpretability of the model. They further discovered that the sulci and white matter were positively correlated with the brain age gap. In contrast, the gyri and periventricular regions were negatively correlated with the brain age gap. By conducting these analyses, the authors shed light on the effects of AD on brain structures. However, it is essential to note that structural changes in hippocampal atrophy occur many years after the accumulation of beta-amyloid (Aβ) pathology in the brain (McKhann et al., 2011; Villemagne et al., 2013), and impairment of functional connections in AD can be detected almost synchronously with Aβ and tau measured using positron emission tomography (PET). Therefore, detecting brain functional changes from resting-state functional MRI (rs-fMRI) could be a more sensitive and earlier method for individuals at risk of AD compared to the brain structural changes (Gonneaud et al., 2021). Previous studies have also demonstrated the utility of rs-fMRI data in predicting brain age and proposed that AD leads to accelerated aging. One such study adopted a Gaussian process regression (GPR) and obtained a mean absolute error (MAE) of 8.195 and a root mean square error (RMSE) of 10.31 in the test set (Millar et al., 2022). However, most existing studies using fMRI data for brain age prediction rely on traditional machine learning methods, and only a few studies apply deep learning networks. To the best of our knowledge, no study employs deep learning models for age prediction trained on graph data derived from functional connectivity (FC) and further applies those for diagnosing AD.
To summarize, while accelerating aging has been extensively investigated as the biomarker in the diagnosis of AD, many unknown aspects still exist for further exploration. Most studies have focused on structural MRI, which reveals structural changes in the brain. However, studies have shown that the fMRI change occurs earlier than structural changes, and resting state FC may be a more sensitive way to detect brain changes in preclinical AD patients (Gonneaud et al., 2021). Furthermore, while most studies using rs-fMRI have employed machine learning models, deep learning models may be more suitable for learning non-linear relationships in brain imaging data (Abrol et al., 2021). Meanwhile, the machine-learned feature may be more informative than those extracted by experience-based artificial methods.
To address these concerns, we have developed an attention-based graph neural networks (GNN) framework to detect accelerated brain aging in patients with AD. First, we extracted a Pearson correlation matrix from rs-fMRI and constructed graph data. Then we trained the GNN model on the graph data to predict brain ages in the HC group. Meanwhile, we applied the model to predict brain ages in patients with MCI and AD. Eventually, the results of our study were compared with those obtained from traditional machine learning models. In this study, we hope to make early diagnosis possible by identifying people who are likely to develop AD as early as possible.
2. Materials and methods
2.1. The dataset
The rs-fMRI data used in our study were obtained from the Alzheimer's Disease Neuroimaging Initiative (adni.loni.usc.edu) (ADNI) (Weiner et al., 2010). The ADNI is a multi-site longitudinal data repository, allowing researchers to access publicly available data upon request and approval, making a significant contribution to promoting the research related to AD. The repository contains a large amount of MRI, PET, and other medical image data. In addition, clinical, genomic, and biomarker data are also provided. Our study used the first available rs-fMRI images obtained during the medical follow-up period of 1,006 subjects from 48 sites, including 535 patients and 471 normal controls (ages 51–97 years). The specific demographic information is shown in Table 1.
2.2. Data preprocessing
All the rs-fMRI images in the Digital Imaging and Communications in Medicine (DICOM) format were preprocessed using the Brainnetome fMRI Toolkit (Xu et al., 2018). The preprocessing steps include DICOM to Neuroimaging Informatics Technology Initiative (NIFTI) conversion, the deletion of the first 10 unstable time points, slice timing correction, registration, normalization to MNI standard space, and denoising. Then, based on the automatic anatomical labeling (AAL) atlas, which divides the human brain into 116 regions (Tzourio-Mazoyer et al., 2002), we extracted the time-series data from 116 brain regions and calculated Pearson correlation coefficients between two brain regions to obtain the Pearson correlation matrix. It is worth noting that the length of the time-series data was different because of different instruments of image acquisition. However, the calculation of the Pearson correlation matrix was not affected. Furthermore, we used Combat (Pomponio et al., 2020) on multi-site data to eliminate site differences. The ablation experiment was performed to verify the effect of Combat.
2.3. Graph data construction
In graph theory, a graph is generally denoted by a pair G = (V, E), where V refers to the vertex set and E represents the edge set. In this study, the FC matrix of each subject was constructed as graph-structured data. Specifically, 116 brain regions defined by the AAL atlas were assigned a corresponding vertex in the graph. The edges in the graph corresponded to the connections between brain regions. The FC values were treated as the features of vertices, while the FC values satisfying a threshold criterion were considered as the features of edges. A complete graph was obtained when the threshold was set to 0. Otherwise, edges that did not meet the threshold requirement were removed from the graph.
2.4. Brain prediction model based on the graph neural network
This study concerns the modeling of FC matrices as graph data, which characterizes the spatial correlations between diverse brain regions, as exemplified by Ktena et al. (2017). In this context, GNN is a prevailing method employed for feature learning of graph data by aggregating information from neighboring regions through convolution, thus displaying remarkable performance in graph representation learning (Zhang et al., 2019).
This study introduces a novel approach for brain age prediction using a GNN model based on transformer convolution (TransformerConv) (Shi et al., 2020). The proposed framework was trained and evaluated on the ADNI dataset, and the graphical representation of the model is depicted in Figure 1. First, the TransformerConv layer was utilized to aggregate and update the node features of the graph data. To be more specific, the input vectors from the graph data were fed into the self-attention layer and calculated the query, key, and value for each region. Then, the self-attention coefficient was calculated, which represented the similarity between the query and the keyword. Finally, the self-attention coefficient was used as the weight, and the output vector of the layer was the weighted sum of the value. Then, the LayerNorm layer (Ba et al., 2016) was applied to normalize the neurons in the middle layer, ensuring distribution stability. The following LeakyReLU (Maas et al., 2013) activation layer was used to enhance the learning ability of the model via non-linear mapping. Finally, the outputs of the two layers mentioned above were used as inputs to a multilayer perceptron (MLP) consisting of a fully connected layer, ReLU layer (Glorot et al., 2011), BatchNorm (Ioffe and Szegedy, 2015), and another fully connected layer, with 32 hidden nodes and one output node, to predict the brain age.
Figure 1. Overview of our proposed framework. First, the function connection matrices were extracted from the rs-fMRI data. Then, the graph data for each subject was constructed and fed into the network. The network contained TransformerConv layers, which updated the aggregated node features of the graph data. Finally, following one layer of MLP, the output was the predicted brain age.
We employed the deep learning approach for predicting age based on graph data from the HC group. Specifically, we utilized a batch size of 16, an initial learning rate of 0.001, and a weight attenuation of 0.001 with the Adam optimizer (Kingma and Ba, 2014) and the cosine annealing learning rate adjustment algorithm. Then, the model was trained based on the 10-fold cross-validation and the mean square error cost function as the loss function. The model's performance was evaluated by three criteria: MAE, RMSE, and the Pearson Correlation Coefficient (PCC) between predicted and chronological age. In addition, we compared our approach to the other six regression methods, namely, support vector regression (SVR) (Vapnik, 1998), GPR (Williams and Rasmussen, 1995), random forest regression (RFR) (Cutler et al., 2012), least absolute shrinkage and selection operator (LASSO) regression (LR) (Tibshirani, 1996), AlexNet (Krizhevsky et al., 2017), and autoencoder (AE) (Heinsfeld et al., 2018).
2.5. Estimated age difference in the patient group
We applied the model to analyze independent test samples of individuals categorized into distinct groups, including significant memory concern (SMC), early mild cognitive impairment (EMCI), MCI, late mild cognitive impairment (LMCI), and AD groups. Additionally, the brain age gap (BAG) as the difference between model-predicted and chronological age was calculated to investigate the potential impact of AD on the brain function.
2.6. Post hoc model based on perturbation
Deep learning is widely used in the diagnosis and treatment of AD. However, its lack of transparency hinders its clinical application as doctors need to understand the impact of changes in FC on patients with AD (Ahmedt-Aristizabal et al., 2021). GNN explainability methods have emerged as a promising solution to address this issue. In recent years, GNN explainability methods have developed rapidly. Existing methods can be divided into two categories: instance-level methods and model-level methods. The former includes the gradient/feature-based, perturbation-based, decomposition, and surrogate methods (Yuan et al., 2022), and the latter includes XGNN (Yuan et al., 2020). In this study, we used the post hoc model based on perturbation GNNExplainer (Ying et al., 2019) to explain networks and help doctors analyze pathology images related to AD. First, the features of each individual graph were perturbed, and the generated node feature mask was multiplied with the input graph element by element to obtain the masked graph. Then, the masked graph was used as the input of the trained model to evaluate and train the mask. The loss function was defined as the mean square error of predicted brain ages using the original graph and masked graph. When nodes with important information were retained, there was little change between the output before and after the perturbation, which would reveal the significant brain regions in the prediction task. Finally, we obtained the explainable weight matrix of each subject. The rows and columns represented the brain region corresponding to the node and the explainable weight of the brain region.
3. Results
3.1. Results of age prediction based on FC
This study evaluated a brain age prediction model based on GNN using a 10-fold cross-validation technique. A total of 471 data from the ADNI dataset were randomly divided into 10 groups for cross-validation. The test set comprised one group of data, while the remaining nine groups served as the training set, which was repeated 10 times. The performance of the model was assessed using the mean and standard deviation of the results obtained from 10-fold cross-validation. The results demonstrated that the model could effectively predict brain age, with an MAE of 5.92 ± 0.62, RMSE of 7.56 ± 0.78, and PCC of 0.44 ± 0.11.
Meanwhile, the results using the data before and after performing the Combat operation were compared. We also presented a comparative study of our GNN-based brain age prediction model against the other six methods (i.e., SVR, GPR, RFR, LR, Alexnet, and AE). It is important to note that traditional machine learning models utilize a one-dimensional vector as input. Thus, we flattened the upper triangular region of a symmetric 116*116 matrix by row, excluding the main diagonal. We then fed a one-dimensional vector with 6,670 features for each subject into the traditional models. Our results in Table 2 indicate that our GNN-based model outperformed other models in terms of MAE and RMSE, and the Combat operation enhanced the performance of the model. Furthermore, to ensure the consistency of different parameters, we also conducted experiments to investigate the influence of the number of heads on model performance. The results presented in Table 3 indicated that the model with a single head achieved the best performance. Additionally, we compared the performance of models with different numbers of self-attention layers, ranging from 1 to 4. The results, shown in Table 4, indicated that the model with two self-attention layers performed the best.
Table 3. Comparison of the performances of models with different numbers of heads in the self-attention layer.
3.2. Accelerated functional brain aging in patient groups
Our study aimed to investigate the relationship between BAG, defined as the difference between predicted and chronological age, and the diagnosis of AD using an independent test set comprising five groups of data, including SMC, EMCI, MCI, LMCI, and AD. We calculated criteria and BAG and analyzed the results using various statistical measures such as MAE, RMSE, and PCC. The findings of our study indicated a significant increase in the values of MAE, RMSE, and BAG in the AD group compared to the HC group, accompanied by a decrease in the PCC value, suggesting accelerated brain function aging in AD patients. Thus, BAG could serve as a valuable biomarker to evaluate the severity of AD and facilitate early diagnosis. The detailed results of each data group are summarized in Table 5.
3.3. The explainability of our model
To account for the neural substrates that impact the age prediction model and to elucidate functional changes that underlie deviations in brain age from typical aging trajectories in individuals with AD, we ranked the average explanatory weight coefficients in the HC and AD groups, respectively, and subsequently found the shared and distinct brain regions. We specifically focused on the top 30% of shared regions, which were primarily located in the medial and paracingulate gyrus, inferior and superior temporal gyrus, transverse temporal gyrus, anterior cingulate and paracingulate gyrus, anterior central gyrus, thalamus, and other brain regions, all of which made notable contributions to age prediction. Additionally, we observed that functional areas implicated in the diagnosis of AD, including the parahippocampal gyrus and amygdala, were also represented among these shared regions. These findings are reported in Table 6. Furthermore, our results revealed five distinct brain regions between the two groups shown in Figure 2 visualized using the BrainNet Viewer (Xia et al., 2013). Specifically, the HC group mainly focused on the inferior middle temporal gyrus, insula, globus pallidum, and two regions of the inferior cerebellum, which were not present in the AD group. In contrast, the AD group mainly focused on the inferior paracentral lobules, supplementary motor areas, fusiform gyrus, inferior cerebellum, and hippocampus, which were not observed in the HC group.
Figure 2. Visualization of the difference between HC and AD in the top 35 brain regions. Insula_L, Insula, left; Pallidum_R, Lenticular nucleus, pallidum, right; Temporal_Mid_L, Middle temporal gyrus, left; Cerebelum_Crus2_R, Cerebelum_Crus2, right; Cerebelum_7b_L, Cerebelum_7b, left; Supp_Motor_Area_L, Supplementary motor area, left; Hippocampus_L, Hippocampus, left; Fusiform_R, Fusiform gyrus, right; Paracentral_Lobule_L, Paracentral lobule, left; Cerebelum_10_R, Cerebelum_10, right.
4. Discussion
In this study, we proposed a prognostication model that leverages GNN augmented with an attention mechanism that can adeptly anticipate brain age. Compared with machine learning methods that employ one-dimensional features, our proposed model retains the topological data of brain regions in the rs-fMRI data, thereby fortifying its FC representation. Furthermore, the attention network manifests superior capability in extracting global features.
Recent studies have identified the accumulation of Aβ plaques and tau protein chain accumulation as pathological markers of AD that are strongly associated with accelerated brain aging (Gonneaud et al., 2021; Mecca and van Dyck, 2021). Although structural MRI is commonly used to detect brain atrophy in AD patients, rs-fMRI is more sensitive to preclinical brain changes that occur earlier in the disease progression (Gonneaud et al., 2021). In this study, rs-fMRI data were used to develop a model based on graph neural networks to predict brain age in HC and different groups of patients. The resulting model achieved relatively accurate predictions, with an MAE of 5.92 as reported in Tables 2–4. However, our study found that the accuracy of brain age prediction based on rs-fMRI data was generally lower than that of studies using structural MRI data. This finding is consistent with previous research indicating that structural MRI is usually more accurate in predicting brain age (Liem et al., 2017; Bashyam et al., 2020; Cole, 2020; Levakov et al., 2020; Dunas et al., 2021; Hwang et al., 2022; Lee et al., 2022). The reason is due to the dynamic nature of rs-fMRI, which captures functional aging changes that may occur earlier than structural changes in the normal population, leading to our predicted older brain age in the subjects.
In this study, we aimed to identify the most critical brain regions for predictive tasks and AD diagnosis by explaining the model. The regions of considerable interest in neuroscience include the cingulate cortex, hippocampal structure, precuneus, inferior temporal gyrus, angular gyrus, and pivotal brain areas involved in default mode networks. Additionally, the supraoccipital gyrus within visual networks, the anterior and posterior central gyrus within sensorimotor networks, as well as the thalamus within the limbic system and other brain regions have garnered considerable attention within the field (Yao et al., 2010; Lei et al., 2020; Ma et al., 2022; Han et al., 2023; Zhou et al., 2023). These regions are crucial for understanding various neural processes and their interconnectivity. In addition, the medial temporal lobe is a vital region implicated in brain senescence and cognitive deterioration in individuals with AD. Our investigation also validated the considerable impact of specific brain regions, including the hippocampus, parahippocampal cortex, and amygdala, on the predictive capacity of our model, as previously demonstrated by Hrybouski et al. (2023). Consistent with an earlier study by Libby et al. (2012) and Liu et al. (2016), our results emphasize the critical role of parahippocampal gyrus connectivity in memory function and its association with AD severity. Furthermore, a positive correlation between cognitive decline and reduced FC between the amygdala and specific brain regions was observed in the studies by Yao et al. (2013) and Yao et al. (2014). Previous studies have consistently shown that hippocampal structure and function changes are closely linked to core memory deficits in individuals with AD (Khazaee et al., 2017; Ibrahim et al., 2021). The middle temporal gyrus, fusiform gyrus, paraventral lobule, cerebellum, and auxiliary motor area have been identified as regions displaying functional alterations that may be associated with accelerated brain aging in individuals with AD (Frisoni et al., 2009; Yao et al., 2010; Brier et al., 2012; Eavani et al., 2018; Hojjati et al., 2019; Ibrahim et al., 2021). These findings validate the utility of our interpretation prediction model as a valuable tool for investigating the advancement of AD.
5. Conclusion
In this study, we introduced a GNN-based prediction model that utilized rs-fMRI data for brain age estimation. Our comparative experiments revealed that our proposed approach outperformed other machine learning and deep learning prediction methods. Additionally, we confirmed the potential of the BAG as a diagnostic marker for AD. Moreover, our post hoc interpretation of the model identified crucial functional brain regions that contributed to age prediction and investigated possible reasons for accelerated brain aging. We anticipate that our research will facilitate the early detection of AD.
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found here: https://adni.loni.usc.edu.
Ethics statement
The studies involving human participants were reviewed and approved by Institutional Review Board. The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author contributions
All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.
Funding
This study is supported by grants from the National Natural Science Foundation of China (Nos. 62276049, 61701078, and 62006038), the Sichuan Province Science and Technology Support Program (Nos. 2019YJ0193 and 2022YFS0180), and the Medico-Engineering Cooperation Funds from University of Electronic Science and Technology of China (No. ZYGX2021YGLH014).
Acknowledgments
Data collection and sharing for this project was funded by the Alzheimer's Disease Neuroimaging Initiative (ADNI) (National Institutes of Health Grant U01 AG024904) and DOD ADNI (Department of Defense award number W81XWH-12-2-0012). The ADNI is funded by the National Institute on Aging, the National Institute of Biomedical Imaging and Bioengineering, and through generous contributions from the following: AbbVie, Alzheimer's Association; Alzheimer's Drug Discovery Foundation; Araclon Biotech; BioClinica, Inc.; Biogen; Bristol-Myers Squibb Company; CereSpir, Inc.; Cogstate; Eisai Inc.; Elan Pharmaceuticals, Inc.; Eli Lilly and Company; EuroImmun; F. Hoffmann-La Roche Ltd. and its affiliated company Genentech, Inc.; Fujirebio; GE Healthcare; IXICO Ltd.; Janssen Alzheimer Immunotherapy Research & Development, LLC.; Johnson & Johnson Pharmaceutical Research & Development LLC.; Lumosity; Lundbeck; Merck & Co., Inc.; Meso Scale Diagnostics, LLC.; NeuroRx Research; Neurotrack Technologies; Novartis Pharmaceuticals Corporation; Pfizer Inc.; Piramal Imaging; Servier; Takeda Pharmaceutical Company; and Transition Therapeutics. The Canadian Institutes of Health Research is providing funds to support ADNI clinical sites in Canada. Private sector contributions are facilitated by the Foundation for the National Institutes of Health (www.fnih.org). The grantee organization is the Northern California Institute for Research and Education, and the study is coordinated by the Alzheimer's Therapeutic Research Institute at the University of Southern California. ADNI data are disseminated by the Laboratory for Neuro Imaging at the University of Southern California.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Abrol, A., Fu, Z., Salman, M., Silva, R., Du, Y., Plis, S., et al. (2021). Deep learning encodes robust discriminative neuroimaging representations to outperform standard machine learning. Nat. Commun. 12, 353. doi: 10.1038/s41467-020-20655-6
Ahmedt-Aristizabal, D., Armin, M. A., Denman, S., Fookes, C., and Petersson, L. (2021). Graph-based deep learning for medical diagnosis and analysis: past, present and future. Sensors (Basel) 21, 4758. doi: 10.3390/s21144758
American Psychiatric Association, D., and Association, A. P. (2013). Diagnostic and Statistical Manual of Mental Disorders: DSM-5. Washington, DC: American Psychiatric Association. doi: 10.1176/appi.books.9780890425596
Anderson, N. D. (2019). State of the science on mild cognitive impairment (MCI). CNS Spectr. 24, 78–87. doi: 10.1017/S1092852918001347
Ba, J. L., Kiros, J. R., and Hinton, G. E. (2016). Layer normalization. arXiv preprint arXiv: 06450.
Bashyam, V. M., Erus, G., Doshi, J., Habes, M., Nasrallah, I., Truelove-Hill, M., et al. (2020). MRI signatures of brain age and disease over the lifespan based on a deep brain network and 14 468 individuals worldwide. Brain 143, 2312–2324. doi: 10.1093/brain/awaa160
Brier, M. R., Thomas, J. B., Snyder, A. Z., Benzinger, T. L., Zhang, D., Raichle, M. E., et al. (2012). Loss of intranetwork and internetwork resting state functional connections with Alzheimer's disease progression. J. Neurosci. 32, 8890–8899. doi: 10.1523/JNEUROSCI.5698-11.2012
Cole, J. H. (2020). Multimodality neuroimaging brain-age in UK biobank: relationship to biomedical, lifestyle, and cognitive factors. Neurobiol. Aging 92, 34–42. doi: 10.1016/j.neurobiolaging.2020.03.014
Cutler, A., Cutler, D. R., and Stevens, J. R. (2012). “Random forests,” in Ensemble machine learning: Methods applications 157–175. doi: 10.1007/978-1-4419-9326-7_5
Dunas, T., Wahlin, A., Nyberg, L., and Boraxbekk, C. J. (2021). Multimodal image analysis of apparent brain age identifies physical fitness as predictor of brain maintenance. Cereb Cortex 31, 3393–3407. doi: 10.1093/cercor/bhab019
Eavani, H., Habes, M., Satterthwaite, T. D., An, Y., Hsieh, M. K., Honnorat, N., et al. (2018). Heterogeneity of structural and functional imaging patterns of advanced brain aging revealed via machine learning methods. Neurobiol. Aging 71, 41–50. doi: 10.1016/j.neurobiolaging.2018.06.013
Frisoni, G. B., Prestia, A., Rasser, P. E., Bonetti, M., and Thompson, P. M. (2009). In vivo mapping of incremental cortical atrophy from incipient to overt Alzheimer's disease. J. Neurol. 256, 916–924. doi: 10.1007/s00415-009-5040-7
Frizzell, T. O., Glashutter, M., Liu, C. C., Zeng, A., Pan, D., Hajra, S. G., et al. (2022). Artificial intelligence in brain MRI analysis of Alzheimer's disease over the past 12 years: A systematic review. Ageing Res. Rev. 77, 101614. doi: 10.1016/j.arr.2022.101614
Gaser, C., Franke, K., Kloppel, S., Koutsouleris, N., Sauer, H., and Alzheimer's Disease Neuroimaging, I. (2013). BrainAGE in mild cognitive impaired patients: predicting the conversion to Alzheimer's disease. PLoS ONE 8, e67346. doi: 10.1371/journal.pone.0067346
Glorot, X., Bordes, A., and Bengio, Y. (2011). “Deep sparse rectifier neural networks,” in Proceedings of the fourteenth International Conference on Artificial Intelligence and Statistics, JMLR Workshop and Conference Proceedings.
Gonneaud, J., Baria, A. T., Pichet Binette, A., Gordon, B. A., Chhatwal, J. P., Cruchaga, C., et al. (2021). Accelerated functional brain aging in pre-clinical familial Alzheimer's disease. Nat. Commun. 12, 5346. doi: 10.1038/s41467-021-25492-9
Han, H., Ge, S., and Wang, H. (2023). Prediction of brain age based on the community structure of functional networks. Biomed. Signal Proc. Control 79, 104151. doi: 10.1016/j.bspc.2022.104151
Heinsfeld, A. S., Franco, A. R., Craddock, R. C., Buchweitz, A., and Meneguzzi, F. (2018). Identification of autism spectrum disorder using deep learning and the ABIDE dataset. NeuroImage: Clin. 17, 16–23. doi: 10.1016/j.nicl.2017.08.017
Hojjati, S. H., Ebrahimzadeh, A., and Babajani-Feremi, A. (2019). Identification of the early stage of Alzheimer's disease using structural MRI and resting-state fMRI. Front. Neurol. 10, 904. doi: 10.3389/fneur.2019.00904
Hrybouski, S., Das, S. R., Xie, L., Wisse, L. E., Kelley, M., Lane, J., et al. (2023). Aging and Alzheimer's disease have dissociable effects on medial temporal lobe connectivity. medRxiv: 2023.2001. 2018, 23284749. doi: 10.1101/2023.01.18.23284749
Hwang, G., Abdulkadir, A., Erus, G., Habes, M., Pomponio, R., Shou, H., et al. (2022). Disentangling Alzheimer's disease neurodegeneration from typical brain ageing using machine learning. Brain Commun. 4, fcac117. doi: 10.1093/braincomms/fcac117
Ibrahim, B., Suppiah, S., Ibrahim, N., Mohamad, M., Hassan, H. A., Nasser, N. S., et al. (2021). Diagnostic power of resting-state fMRI for detection of network connectivity in Alzheimer's disease and mild cognitive impairment: A systematic review. Hum. Brain Mapp. 42, 2941–2968. doi: 10.1002/hbm.25369
Ioffe, S., and Szegedy, C. (2015). “Batch normalization: Accelerating deep network training by reducing internal covariate shift,” in International Conference on Machine Learning.
Khazaee, A., Ebrahimzadeh, A., Babajani-Feremi, A., and Initiative, A. s. D. N. (2017). Classification of patients with MCI and AD from healthy controls using directed graph measures of resting-state fMRI. Behav. Brain Res. 322, 339–350. doi: 10.1016/j.bbr.2016.06.043
Kingma, D. P., and Ba, J. (2014). Adam: A method for stochastic optimization. arXiv preprint arXiv: 6980.
Krizhevsky, A., Sutskever, I., and Hinton, G. E. (2017). Imagenet classification with deep convolutional neural networks. Commun. ACM 60, 84–90. doi: 10.1145/3065386
Ktena, S. I., Parisot, S., Ferrante, E., Rajchl, M., Lee, M., Glocker, B., et al. (2017). “Distance metric learning using graph convolutional networks: Application to functional brain networks,” in Medical Image Computing and Computer Assisted Intervention– MICCAI 2017, 20th International Conference (Quebec City, QC, Canada: Springer). doi: 10.1007/978-3-319-66182-7_54
Lee, J., Burkett, B. J., Min, H.-K., Senjem, M. L., Lundt, E. S., Botha, H., et al. (2022). Deep learning-based brain age prediction in normal aging and dementia. Nature Aging 2, 412–424. doi: 10.1038/s43587-022-00219-7
Lei, B., Cheng, N., Frangi, A. F., Tan, E. L., Cao, J., Yang, P., et al. (2020). Self-calibrated brain network estimation and joint non-convex multi-task learning for identification of early Alzheimer's disease. Med. Image Anal. 61, 101652. doi: 10.1016/j.media.2020.101652
Levakov, G., Rosenthal, G., Shelef, I., Raviv, T. R., and Avidan, G. (2020). From a deep learning model back to the brain-Identifying regional predictors and their relation to aging. Hum. Brain Mapp. 41, 3235–3252. doi: 10.1002/hbm.25011
Libby, L. A., Ekstrom, A. D., Ragland, J. D., and Ranganath, C. (2012). Differential connectivity of perirhinal and parahippocampal cortices within human hippocampal subregions revealed by high-resolution functional imaging. J. Neurosci. 32, 6550–6560. doi: 10.1523/JNEUROSCI.3711-11.2012
Liem, F., Varoquaux, G., Kynast, J., Beyer, F., Masouleh, S. K., Huntenburg, J. M., et al. (2017). Predicting brain-age from multimodal imaging data captures cognitive impairment. Neuroimage 148, 179–188. doi: 10.1016/j.neuroimage.2016.11.005
Liu, J., Zhang, X., Yu, C., Duan, Y., Zhuo, J., Cui, Y., et al. (2016). Impaired parahippocampus connectivity in mild cognitive impairment and Alzheimer's disease. J. Alzheimers Dis. 49, 1051–1064. doi: 10.3233/JAD-150727
Ma, J., Zhang, J., and Wang, Z. (2022). Multimodality Alzheimer's Disease Analysis in Deep Riemannian Manifold. Inf. Proc. Manage. 59, 102965. doi: 10.1016/j.ipm.2022.102965
Maas, A. L., Hannun, A. Y., and Ng, A. Y. (2013). Rectifier Nonlinearities Improve Neural Network Acoustic Models. Atlanta, Georgia, USA: Proc. ICML.
McKhann, G. M., Knopman, D. S., Chertkow, H., Hyman, B. T., Jack, C. R. Jr., Kawas, C. H., et al. (2011). The diagnosis of dementia due to Alzheimer's disease: recommendations from the National Institute on Aging-Alzheimer's Association workgroups on diagnostic guidelines for Alzheimer's disease. Alzheimers Dement 7, 263–269. doi: 10.1016/j.jalz.2011.03.005
Mecca, A. P., and van Dyck, C. H. (2021). Alzheimer's and dementia: the journal of the Alzheimer's Association. Alzheimer's Dement. 17, 327–406. doi: 10.1002/alz.12190
Millar, P. R., Luckett, P. H., Gordon, B. A., Benzinger, T. L. S., Schindler, S. E., Fagan, A. M., et al. (2022). Predicting brain age from functional connectivity in symptomatic and preclinical Alzheimer disease. Neuroimage 256, 119228. doi: 10.1016/j.neuroimage.2022.119228
Pomponio, R., Erus, G., Habes, M., Doshi, J., Srinivasan, D., Mamourian, E., et al. (2020). Harmonization of large MRI datasets for the analysis of brain imaging patterns throughout the lifespan. Neuroimage 208, 116450. doi: 10.1016/j.neuroimage.2019.116450
Reitz, C., Brayne, C., and Mayeux, R. (2011). Epidemiology of Alzheimer disease. Nat. Rev. Neurol. 7, 137–152. doi: 10.1038/nrneurol.2011.2
Sajedi, H., and Pardakhti, N. (2019). Age prediction based on brain MRI image: a survey. J. Med. Syst. 43, 1–30. doi: 10.1007/s10916-019-1401-7
Shi, Y., Huang, Z., Feng, S., Zhong, H., Wang, W., and Sun, Y. (2020). Masked label prediction: Unified message passing model for semi-supervised classification. arXiv preprint arXiv:2009.03509. doi: 10.24963/ijcai.2021/214
Tahami Monfared, A. A., Byrnes, M. J., White, L. A., and Zhang, Q. (2022). Alzheimer's disease: epidemiology and clinical progression. Neurol. Ther. 11, 553–569. doi: 10.1007/s40120-022-00338-8
Tibshirani, R. (1996). Regression shrinkage and selection via the lasso. J. Royal Stat. Soc. 58, 267–288. doi: 10.1111/j.2517-6161.1996.tb02080.x
Tzourio-Mazoyer, N., Landeau, B., Papathanassiou, D., Crivello, F., Etard, O., Delcroix, N., et al. (2002). Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain. Neuroimage 15, 273–289. doi: 10.1006/nimg.2001.0978
Vapnik, V. (1998). “The support vector method of function estimation,” in Nonlinear Modeling: Advanced Black-Box Techniques 55–85. doi: 10.1007/978-1-4615-5703-6_3
Villemagne, V. L., Burnham, S., Bourgeat, P., Brown, B., Ellis, K. A., Salvado, O., et al. (2013). Amyloid beta deposition, neurodegeneration, and cognitive decline in sporadic Alzheimer's disease: a prospective cohort study. Lancet Neurol. 12, 357–367. doi: 10.1016/S1474-4422(13)70044-9
Warren, S. L., and Moustafa, A. A. (2023). Functional magnetic resonance imaging, deep learning, and Alzheimer's disease: A systematic review. J. Neuroimag. 33, 5–18. doi: 10.1111/jon.13063
Weiner, M. W., Aisen, P. S., Jack, C. R. Jr., Jagust, W. J., Trojanowski, J. Q., Shaw, L., et al. (2010). The Alzheimer's disease neuroimaging initiative: progress report and future plans. Alzheimers Dement 6, 202–211.e207. doi: 10.1016/j.jalz.2010.03.007
Williams, C., and Rasmussen, C. (1995). “Gaussian processes for regression,” in Advances in Neural Information Processing Systems 8.
Xia, M., Wang, J., and He, Y. (2013). BrainNet Viewer: a network visualization tool for human brain connectomics. PLoS ONE 8, e68910. doi: 10.1371/journal.pone.0068910
Xu, K., Liu, Y., Zhan, Y., Ren, J., and Jiang, T. (2018). BRANT: A Versatile and Extendable Resting-State fMRI Toolkit. Front. Neuroinform. 12, 52. doi: 10.3389/fninf.2018.00052
Yao, H., Liu, Y., Zhou, B., Zhang, Z., An, N., Wang, P., et al. (2013). Decreased functional connectivity of the amygdala in Alzheimer's disease revealed by resting-state fMRI. Eur. J. Radiol. 82, 1531–1538. doi: 10.1016/j.ejrad.2013.03.019
Yao, H., Zhou, B., Zhang, Z., Wang, P., Guo, Y., Shang, Y., et al. (2014). Longitudinal alteration of amygdalar functional connectivity in mild cognitive impairment subjects revealed by resting-state FMRI. Brain Connect 4, 361–370. doi: 10.1089/brain.2014.0223
Yao, Z., Zhang, Y., Lin, L., Zhou, Y., Xu, C., Jiang, T., et al. (2010). Abnormal cortical networks in mild cognitive impairment and Alzheimer's disease. PLoS Comput. Biol. 6, e1001006. doi: 10.1371/journal.pcbi.1001006
Ying, Z., Bourgeois, D., You, J., Zitnik, M., and Leskovec, J. (2019). “Gnnexplainer: Generating explanations for graph neural networks,” in Advances in Neural Information Processing Systems 32.
Yuan, H., Tang, J., Hu, X., and Ji, S. (2020). “XGNN: towards model-level explanations of graph neural networks,” in Proceedings of the 26th ACM SIGKDD International Conference on Knowledge Discovery and Data Mining 430–438. doi: 10.1145/3394486.3403085
Yuan, H., Yu, H., Gui, S., and Ji, S. (2022). “Explainability in graph neural networks: A taxonomic survey,” in IEEE Transactions on Pattern Analysis Machine Intelligence. doi: 10.1109/TPAMI.2022.3204236
Zhang, S., Tong, H., Xu, J., and Maciejewski, R. (2019). Graph convolutional networks: a comprehensive review. Comput. Soc. Netw. 6, 1–23. doi: 10.1186/s40649-019-0069-y
Keywords: Alzheimer's disease, brain aging, resting-state functional magnetic resonance imaging, graph neural network, age prediction
Citation: Gao J, Liu J, Xu Y, Peng D and Wang Z (2023) Brain age prediction using the graph neural network based on resting-state functional MRI in Alzheimer's disease. Front. Neurosci. 17:1222751. doi: 10.3389/fnins.2023.1222751
Received: 15 May 2023; Accepted: 08 June 2023;
Published: 30 June 2023.
Edited by:
Ni Zhang, Sichuan University, ChinaCopyright © 2023 Gao, Liu, Xu, Peng and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zhengning Wang, emhlbmduaW5nLndhbmcmI3gwMDA0MDt1ZXN0Yy5lZHUuY24=
†Data used in preparation of this article were obtained from the Alzheimer's Disease Neuroimaging Initiative (ADNI) database (adni.loni.usc.edu). As such, the investigators within the ADNI contributed to the design and implementation of ADNI and/or provided data but did not participate in analysis or writing of this report. A complete listing of ADNI investigators can be found at: http://adni.loni.usc.edu/wp-content/uploads/how_to_apply/ADNI_Acknowledgement_List.pdf