ABSTRACT
Brain-age can be inferred from structural neuroimaging and compared to chronological age (brain-age delta) as a marker of biological brain aging. Accelerated aging has been found in neurodegenerative disorders like Alzheimer’s disease (AD), but its validation against markers of neurodegeneration and AD is lacking. Here, imaging-derived measures from the UK Biobank dataset (N=22,661) were used to predict brain-age in 2,314 cognitively unimpaired (CU) individuals at higher risk of AD and mild cognitive impaired (MCI) patients from four independent cohorts with available biomarker data: ALFA+, ADNI, EPAD and OASIS. Brain-age delta was associated with abnormal amyloid-β, more advanced stages (AT) of AD pathology and APOE-ε4 status. Brain-age delta was positively associated with plasma neurofilament light, a marker of neurodegeneration, and sex differences in the brain effects of this marker were found. These results validate brain-age delta as a non-invasive marker of biological brain aging related to markers of AD and neurodegeneration.
INTRODUCTION
Age is the main risk factor for Alzheimer’s Disease (AD) and most neurodegenerative diseases. However, the mechanisms underlying this association are still poorly understoodd (Fjell et al., 2014). Both normal aging and AD are associated with region-specific cerebral morphological changes characterized by the occurrence of atrophy (Bakkour et al., 2013; Fjell et al., 2014). Both aging and AD have differential and partially overlapping effects on specific regions of the cerebral cortex like, for instance, the dorsolateral prefrontal cortex (Bakkour et al., 2013; Fjell et al., 2014; Pichet Binette et al., 2020). Conversely, some regions are predominantly affected by age (e.g., calcarine cortex) and some others are predominantly affected by AD (e.g., medial temporal cortex) (Bakkour et al., 2013). A better understanding of the mechanistic links between the brain aging process and neurodegenerative diseases is an urgent priority to develop effective strategies to deal with their rising burden amid an ageing population (Franke & Gaser, 2019). Therefore, a growing amount of research is focusing on using neuroimaging techniques to develop a biomarker of biological brain aging. In this framework, the concept of brain-age has emerged as an appealing comprehensive marker that enables determining on an individual basis, the risk for age-associated brain diseases (James H. Cole et al., 2017; James H. Cole & Franke, 2017; Franke et al., 2010; Franke & Gaser, 2019). However, this is a challenging task because, even though the cerebral structural changes related to aging are well established, the older population is characterized by substantial variation in neurobiological aging trajectories (J. H. Cole et al., 2018; Fjell et al., 2014).
Recently, machine learning techniques have gained popularity as brain-age prediction models (James H. Cole et al., 2017; Dafflon et al., 2020; de Lange et al., 2019; Franke & Gaser, 2019), due to their ability in identifying relevant data-driven patterns within complex data (Zhavoronkov et al., 2019). These models learn the association between chronological age and cerebral morphological features derived from structural magnetic resonance imaging (MRI) in healthy individuals, yielding a predicted brain-age for each individual. Individuals with a predicted brain-age higher than their chronological age may have an “older” brain than expected, whereas an individual with an estimated brain-age lower than their chronological age has a “younger” brain. Subtracting chronological age from estimated brain-age hence provides an estimate of accelerated brain aging, namely the brain-age delta. Recent literature has shown the adequacy of using a brain-age predicted measurement in the assessment of the clinical severity of AD, by finding higher brain-age deltas in AD and individuals with mild cognitive impairment (MCI) with respect to cognitively unimpaired (CU) individuals (Beheshti et al., 2018; Kaufmann et al., 2019). A higher brain-age delta has also been reported in other diseases, such as multiple sclerosis, epilepsy and psychiatric disorders, with respect to healthy controls (Beheshti et al., 2018; Kaufmann et al., 2019). In addition, brain-age delta has also been associated with other biological measures such as: lifestyle factors (James H. Cole, 2020), cognition (Beheshti et al., 2018; James H. Cole, 2020) hypertension (de Lange, et al., 2020) and prediction of mortality (J. H. Cole et al., 2018).
Even though these studies support the association of brain-age delta as a biomarker of biological aging with relevance to various brain diseases, there are no comprehensive studies validating this measurement in association with specific biological markers of AD pathology (i.e. Amyloid-β [Aβ] and tau pathology), neurodegeneration and cerebrovascular disease. This is a very relevant aspect since the recent AD research framework criteria defines AD as a biological construct, namely the presence of both abnormal Aβ (A+) and tau (T+) biomarkers, regardless of clinical manifestations(Jack et al., 2018). The term “Alzheimer’s pathological change” is proposed whenever there is evidence of Aβ but not tau pathology (A+T-). The umbrella term “Alzheimer’s continuum” includes both “Alzheimer’s pathological change” (A+T-) and “Alzheimer’s Disease” (A+T+). Under this definition, A-T+ individuals would not fall into the AD continuum. Then, under this framework, neurodegeneration biomarkers (N) and cognitive status (i.e. CU, MCI and dementia syndromes) are used to stage disease progression.
A recent study used brain-age measurements to identify amnestic MCI (aMCI), the typical clinical presentation of prodromal AD, from other individuals with MCI, by studying the association with AD risk factors such as apolipoprotein E (APOE) and Aβ(Huang et al., 2021). Another study focusing on the impact of training the brain-age prediction model in individuals with Aβ pathology (Aβ+) showed that CU Aβ+ individuals had a higher brain-age delta than CU Aβ-individuals (Ly et al., 2020). Nonetheless, there remains a need to study the associations between brain-age prediction and AD as well as neurodegeneration biomarkers in preclinical stages in different and independent cohorts and in a larger sample size. Moreover, given that female individuals have a higher AD prevalence compared to males (Nebel et al., 2018) and display different lifetime trajectories in the brain morphological features (Gennatas et al., 2017), it is of interest to determine the effect of sex on brain age delta and its interaction with AD biomarkers. Literature describes sex differences in AD biomarkers, such as that females with abnormal Aβ who are APOE-ε4 carriers show greater subsequent increase in cerebrospinal fluid (CSF) tau than their male counterparts (Buckley et al., 2019), or that females with higher Aβ burden show higher entorhinal cortical tau than their male counterparts (Buckley et al., 2019). Conversely, levels of the neurodegeneration biomarker CSF neurofilament light (NfL) have been widely reported to be higher in males than in females (Mielke, 2020; Milà-Alomà et al., 2020). In line with this, AD risk factors are associated with greater brain aging in women than men (Subramaniapillai et al., 2021).
Therefore, in the present study, we aim to validate brain-age delta as a clinically relevant marker related to markers of AD and neurodegeneration. For this purpose, we determine the association between the predicted structural brain-age delta with biomarkers and risk factors for AD and neurodegeneration in non-demented individuals, as well as to study the effect of sex on these associations. We trained a model to predict the brain-age separately for females and males, using machine learning on imaging-derived measures of cortical thickness, cortical volume, and subcortical volume from the UK BioBank cohort (N=22,661). Using this model, we then estimated brain-age in four independent cohorts: ALFA+ (N=380), ADNI (N=719), EPAD (N=808) and OASIS (N=407). In each cohort, we studied the associations of brain-age delta with biomarkers of AD pathology (CSF Aβ and p-tau as continuous values, as well as categorized in AT stages), the APOE-ε4 genotype which is the main genetic risk factor for AD, neurodegeneration (CSF and plasma NfL), and small vessel disease (White Matter Hyperintensities [WMH]). Finally, we studied the sex differences in brain age prediction and the sex effects with these biomarkers on brain-age delta.
RESULTS
Participants’ Characteristics
Table 1 summarizes the demographic characteristics of the cohorts included in the study. ADNI and EPAD cohorts included both CU and MCI individuals, while the UK BioBank, ALFA+ and OASIS cohorts only included CU individuals. Table 2 summarizes the variables used to study the associations with brain-age delta, which included biomarkers for AD (Aβ positron emission tomography [PET] and CSF Aβ and p-tau), neurodegeneration (CSF and plasma NfL), and cerebrovascular pathology (WMH on MRI), as well as the aging signature composite(Bakkour et al., 2013), both cross-sectional and longitudinally. The aging signature composite is a map of specific brain regions that undergo cortical thinning in normal aging, which has been used as a proxy measurement for brain aging. These validation variables were correlated with chronological age for all cohorts (see Supplementary Table 1). Some of the participants for ALFA+ (N=25), ADNI (N=116) and EPAD (N=71) fell into the A-T+ group, corresponding to non-AD pathologic change. Since our aim was to specifically validate the brain-age delta measurements in the AD continuum, we excluded these participants from subsequent analyses; and they are reported within Table 1 and Table 2 solely for descriptive purposes. In addition, the number of MCI individuals with available data of CSF NfL and of aging signature change was relatively low and, therefore, these variables were excluded from the analysis in MCI individuals.
Brain-Age Prediction and Chronological Age
We trained the prediction model using the UK BioBank cohort and tested the model using four independent cohorts (ALFA+, ADNI, EPAD and OASIS), as shown in Supplementary Fig. 1. Table 3 shows the prediction accuracy before and after age-bias correction for the combined female and male predictions. The average prediction accuracy of the model run on UK BioBank using ten-fold cross-validation as measured by the mean absolute error (MAE) and by Pearson’s correlation were, originally, MAE=4.19 and R=0.71 and, after correction, MAE=2.95 and R=0.89 (Table 3 and Supplementary Fig. 2).
We then investigated the association of predicted brain-age with chronological age on each of the independent cohorts. All the cohorts showed a similar positive correlation and fitting performance metrics as measured by the mean absolute error (MAE), R and root mean squared error (RMSE) between chronological age and predicted brain-age. Correlation coefficients were not different between cohorts (P>0.05, for all comparisons, see Supplementary Table 2). As an example, after bias correction, the highest numerical difference was between ALFA+ and OASIS, with quite similar MAE=3.25 and R=0.729 and MAE=3.81 and R=0.910, respectively.
In order to study the effect of sex on brain age prediction, we also computed the performance metrics stratified by females and males (Supplementary Table 3 and 4). Correlations and fitting performance metrics were not significantly different between females and males (Pearson’s r (William’s test), P>0.05; RMSE (F-test) P>0.05), see Supplementary Table 5. Plots of the correlations between predicted brain-age and chronological age for females and males in each of the cohorts can be seen in the Supplementary Fig. 3.
Brain Regions Associated with Aging
We computed the SHapley Additive exPlanation (SHAP) values, which reflect the marginal contribution of each brain region to the brain-age prediction, using the UKBioBank dataset. SHAP values interpret the impact in the prediction of the values of volume or cortical thickness for a given brain region. In other words, they reflect the most important features that consistently influenced the prediction of brain-age and whether the decrease or increase of each region impacted into predicting a higher or lower brain-age. The SHAP values were computed separately for females and males. We compared the regions with higher SHAP values for females and males, and vice-versa, by averaging the SHAP values within each sex separately and then subtracting the mean SHAP of males to the mean SHAP of females.
There were regions whose SHAP values were high in both females and males, including the volumes of the amygdala, nucleus accumbens, cerebellar white matter, lateral ventricles and the insula, as well as the cortical thickness of the superior-temporal cortex. All the brain regions with consistent highest SHAP values for females and males are shown in Fig. 1a-b, as well as the effect of each region (larger or lower value) on predicting a higher brain-age. Conversely, the thickness of regions such as the transverse temporal cortex, the pars triangularis, the inferior parietal cortex and the left frontal pole thickness, as well as the volume of the left entorhinal cortex had higher SHAP values in females than in males, while the opposite occurred with the thickness of the left isthmus cingulate and the right cuneus and the cortical volume of the superior frontal and right rostral middle regions (Fig. 1c).
In Fig. 1d we can see the aging trajectories of three regions whose SHAP values were different for females and males. For example, the bi-lateral superior frontal volumes decreased more over the years within males than females. This result was seen as the interaction of sex with age (Pinteratcion<0.001). We also found an interaction effect of sex and age for the isthmus cingulate thickness (Pinteratcion<0.001), by which the thickness of males decreased more over the years than from the females. On the contrary, we also found regions, such as the middle temporal thickness, which followed the same trajectory over time for both sexes (Pinteraction=0.671), but which was lower for females than for males (P<0.001).
Associations with AD biomarkers and risk factors
We studied the association between brain-age delta and AD biomarker classifications (Aβ status, AT stages) and APOE-ε4 status in all the independent cohorts pooled together, with a linear model adjusting for the effect of age and sex (Figure 2 and Table 4). Aβ status was defined by CSF (ALFA+, ADNI and EPAD) or amyloid PET (OASIS) using pre-established cut-off values (Hansson et al., 2018; Milà-Alomà et al., 2020; Salvadó et al., 2019; Schindler et al., 2018). Brain-age delta was higher in MCI with respect to CU individuals (P<0.001). In both CU and MCI, a higher brain-age delta was significantly associated with abnormal Aβ status (CU: P<0.001 and MCI: P<0.001) and with progressive AT stages (CU: P<0.001 and MCI: P<0.001) (see Table 4 and Supplementary Table 6 for more details). The mean brain-age delta values for the different Aβ status and AT stages can be found in Supplementary Table 7. The brain-age effect on AT stages was progressive, as that of the A+T-group was larger than that of A-T-, while the brain-age delta of A+T+ was larger than those of the other two previous stages (Table 4 and Fig 2a). Brain-age delta was also significantly associated with APOE status (CU: P<0.001 and MCI: P=0.029). In particular, APOE-ε4 carriers had larger brain-age deltas (i.e., older-appearing brains than expected for their chronological age) compared to APOE-ε33 individuals for both CU (β=0.105, P=0.032) and MCI (β=0.266, P=0.005) (see Table 4 and Figure 2). The mean brain-age delta values for the different APOE status can be found in Supplementary Table 7. These results were consistent with the results from the within-cohort analyses (see Supplementary Table 8).
We next studied the association between brain-age delta and AD biomarkers and risk factors stratified by sex (Table 5). In general, the same associations found with the whole sample was seen for females and males separately. However, although a higher brain-age delta was significantly associated with progressive AT stages both for females (CU: P<0.001 and MCI: P<0.001) and males (CU: P=0.009 and MCI: P<0.001), brain-age delta of A+T+ was significantly larger than those of the other two previous stages (A-T-and A+T-) in CU females (β=0.431, P=0.001) but not in CU males (β=0.139, P=0.364). We conducted regression analyses to test the interaction effect of sex and AT stages on CU brain-age delta. Although we found a trend by which the proportion of A+T+ with larger brain-age delta was larger in females than in males, the interaction effect was not significant (Pinteraction=0.071).
Associations with neurodegeneration biomarkers
We next tested the associations between brain-age delta and neurodegeneration biomarkers (Fig. 2 and Table 4). CSF NfL, plasma NfL and longitudinal change of the aging signature were available in ALFA+ and ADNI. The positive associations between brain-age deltas and plasma NfL were significant within the CU (β=0.142, P=0.002) and MCI individuals (β=242, P<0.001). CSF NfL was not significantly associated with brain-age delta (β=0.077, P=0.122). The aging signature composite at both visits was negatively associated with brain-age delta (Visit 1: β=-0.366, P<0.001 and Visit 2: β=- 0.302, P<0.001). That is, larger brain-age delta was associated with reduced cortical thickness in aging-vulnerable regions. However, the association between the longitudinal change in the aging signature and brain-age delta was not statistically significant (β=0.014, P=0.799).
We next studied the association between brain-age delta and neurodegeneration biomarkers stratified by sex (Table 5). The associations between brain-age delta and CSF NfL were significant within the CU females (β=0.131, P=0.042), but not within the CU males (β=-0.004, P=0.959). However, the interaction effect of sex and CSF NfL on CU brain-age delta (Fig. 3a) did not reach significance (Pinteraction=0.170). In the same line, the associations between brain-age delta and plasma NfL were significant within the CU and MCI females (CU: β=0.193, P=0.001 and MCI: β=0.342, P=0.001), but not within the males (CU: β=0.079, P=0.254 and MCI: β=0.157, P=0.086). The interaction effect of sex and plasma NfL on brain-age delta (Fig. 3a) revealed a trend within CU individuals, by which, although not significant, plasma NfL was larger on females when the brain-age delta showed larger values (older-appearing brain) than in males (CU: Pinteraction=0.092 and MCI: Pinteraction=0.194).
In addition, we tested for the interaction effect of age and these biomarkers on brain-age delta. We found a significant interaction effect of age and CSF NfL on CU brain-age delta (Pinteraction< 0.001) within the CU individuals (Fig 4a), by which the measures of CSF NfL were higher with age and with larger brain-age deltas (older-appearing brain). When stratifying by sex, this interaction effect of age was seen in females (Pinteraction< 0.001), but not in males (Pinteraction=0.241). Regarding plasma NfL (Fig. 4b), although we found a similar trend by which the measures of plasma NfL were higher with age and with larger brain-age deltas for CU and MCI individuals, the interaction effects were not significant (CU: Pinteraction=0.136 and MCI: Pinteraction=0.145). When stratifying by sex, this interaction effect of age was seen in CU females (Pinteraction<0.016) and not in CU males (Pinteraction=0.656). On the contrary, this interaction effect of age and plasma NfL on brain-age delta was seen in MCI males (Pinteraction<0.017) and not in MCI females (Pinteraction=0.621).
Associations with markers of cerebrovascular disease
We lastly tested the associations between brain-age delta and markers of cerebrovascular disease WMH; WMH data were available in ALFA+, ADNI and EPAD. In both CU and MCI, brain-age delta was significantly associated with WMH (CU: β=0.171, P<0.001 and MCI: β=0.222, P<0.001) (see Table 4). These results were consistent with the results from the within-cohort analyses (see Supplementary Table 7).
When studying the association between brain-age delta and WMH stratified by sex (Table 5) we found that the brain-age delta was positively associated with WMH both in CU females (β=0.202, P<0.001) and CU males (β=0.132, P=0.007). The interaction effect of sex and WMH on CU brain-age delta (Fig. 3a) was not significant (Pinteraction=0.182). Conversely, we found that the brain-age delta was positively associated with WMH MCI males (β=0.154, P<0.001), but not in females (β=0.140, P=0.100). The interaction effect of sex and WMH on MCI brain-age delta (Fig. 3b) was also not significant (Pinteraction=0.112).
DISCUSSION
In this study, we show that, in non-demented individuals, the predicted brain-age delta is associated with specific AD biomarkers (amyloid-β PET, CSF Aβ42 and CSF pTau) and risk factors (APOE genotype), as well as with unspecific neurodegeneration biomarkers (plasma NfL), and markers of cerebrovascular disease (WMH volume). Our results also indicate that there are sex differences in the development of brain aging trajectories. Taken together, our findings validate the use of machine learning predicted brain-age deltas as biomarkers of brain aging and AD pathology.
We have studied, to our knowledge for the first time, the associations between brain-age delta and different biomarkers of AD pathology and neurodegeneration in CU individuals. We are aware of the complexity of disentangling the effects of aging and pathology in brain aging. In this work, we do not aim to disentangle to what extent the brain structural differences are caused by AD pathology (as measured by the biomarkers that we study) or aging. Here, we show that an unspecific estimation of biological brain aging, agnostic of the underlying mechanisms is associated with the specific biological process of AD.
Regarding the associations with AD biomarkers and risk factors, regression analyses revealed significant positive associations of brain-age delta with increased Aβ pathology and with AT stages for CU and MCI individuals. We also found significant associations with APOE status in the CU and MCI individuals, in which larger brain-age deltas were associated with the presence of APOE-ε4 allele. This result is in line with previous literature that has shown that APOE-ε4 carriership may accelerate AD-related brain atrophy (Evans et al., 2014; Filippini et al., 2011), as accelerated brain aging has also been found in MCI and AD patients (Beheshti et al., 2018; Kaufmann et al., 2019). The association of brain-age delta with APOE-ε4 was also previously studied, for which significant associations were found in MCI individuals (J. H. Cole et al., 2018; Löwe et al., 2016). Taken together, our results advocate for an effect of APOE-ε4 in physiological brain aging, albeit of a lesser magnitude than when AD pathology is present. These results with AD biomarkers and risk factors were highly reproducible in within-cohort analyses.
With the aim of studying the associations between brain-age delta and neurodegeneration, we computed the associations with NfL, a marker of neuro-axonal damage (Khalil et al., 2018) which can be measured both in CSF and in plasma, and with longitudinal changes in the aging signature composite as marker of age-related brain atrophy. The particular use of NfL in this context is supported by its correlation with age throughout the lifespan, as well as its strong association with all-cause mortality in the elderly (Kaeser et al., 2021). We found significant positive associations between brain-age delta and plasma NfL both in CU and MCI individuals, but we did not find significant associations between CSF NfL and CU brain-age delta. Even though we found NfL to be positively associated with chronological age in CU individuals, in line with previous studies (Beheshti et al., 2018; Kaufmann et al., 2019; Khalil et al., 2020; Milà-Alomà et al., 2020), the expected annual change of NfL in CU individuals whose mean age range was 65 years old is around 3.5% (Khalil et al., 2020). Therefore, we expected to find weak associations with brain-age delta. Still, we found a significant interaction effect of age and CSF NfL on CU brain-age delta, for which the individuals with larger brain-age delta had increased CSF NfL over the years, whereas the decelerated ones remained more stable. This trend was also seen in plasma NfL for CU and MCI individuals, although the interaction did not reach significance. Taken together, a strong association between brain-age delta and plasma NfL was observed across all individuals whereas the association on CU brain-age delta and CSF NfL was milder and could only be detected as an interaction with age. Overall, these mild associations between brain-age delta and NfL suggest that the morphological effects of aging in the brain are not fully driven by neurodegeneration. In this regard, it is worth noting that cortical thinning with age has also been linked to loss of volume of the neuropil and other non-neuronal processes which are not necessarily implicated in neurodegeneration (Vidal-Pineiro et al., 2020).
We also studied the associations between brain-age delta and cerebrovascular disease. Regression analyses revealed significant associations of brain-age delta with increased WMH for both CU and MCI individuals. These results were expected, as the increase in WMH with age has been previously studied (Maniega et al., 2015) and it has been shown that individuals with high WMH burden display spatial patterns of atrophy that partially overlap with those of brain aging (Brugulat-Serrat, Salvadó, et al., 2020; Habes et al., 2016). In addition, WMH have been linked to cognitive dysfunction and dementia (Brugulat-Serrat, Salvadó, et al., 2020; Brugulat-Serrat, Salvadó, et al., 2020; Maniega et al., 2015) and a potential pathway has been proposed, in which small vessel cerebrovascular disease affects cognition by promoting neurodegenerative changes (Rizvi et al., 2018). In summary, our results support an effect of cerebrovascular disease in physiological brain aging.
Brain structure aging-associated changes have been widely studied (Bakkour et al., 2013; Fjell et al., 2014). In our study, the brain regions that had highest impact on the machine learning prediction were similar to regions previously mentioned in literature (Arenaza-Urquijo et al., 2019; Bakkour et al., 2013). We found an overlap between some of our selected regions and regions included in the aging signature for both females and males, such as the precentral sulcus, insula, superior frontal and rostral middle frontal regions. In addition, the effect of sex on age-related changes in brain structure has also been studied in the recent years, with some studies reporting age–sex interactions in volumes of certain brain structures (Coffey et al., 1998; DeCarli et al., 2005), and others not finding such interactions (Greenberg et al., 2008). In our study, we found that, even though most of the regions with highest impact were the same for males and females, there were some regions that were sex specific. In particular, we found reduction in the superior-frontal, isthmus-cingulate and pars orbitalis regions within males and regions such as inferior-parietal, pars triangularis and paracentral within females. Most of these sex-specific regions were in concordance with previous studies (Armstrong et al., 2019; Podgórski et al., 2021). The mechanisms underlying these sex-specific brain aging differences are not well-known. Sexual hormones such as estrogen, progesterone and androgen could play a role in brain atrophy (Armstrong et al., 2019); in the WHIMS-MRI study (Resnick et al., 2009), women under menopausal hormone therapy were associated with greater brain atrophy. Others, however, have proposed that estrogen and progesterone may play a protective effect in women (Armstrong et al., 2019; Green & Simpkins, 2000). Other possible biological mechanisms influencing these results could be developmental (Baron-Cohen et al., 2005) or the influence of a greater presence of adverse lifestyle-related factors in men (DeCarli et al., 2005).
In line with the effect of sex on age-related changes in brain structure, we studied the effect of sex on the associations between brain-age delta and the above-mentioned variables. Regarding the AD biomarkers and risk factors, we found that the association between brain-age delta and a larger proportion of A+T+ was only seen in females for the CU individuals, but the interaction effect of sex and AT stages on brain-age delta was not significant. Regarding the neurodegeneration variables, we found a mild positive association between brain-age delta and CSF NfL in CU females, as well as a stronger positive association between brain-age delta and plasma NfL in CU and in MCI females. These associations were not seen in their male counterparts. However, the interaction effect of sex and CSF and plasma NfL on CU and MCI brain-age delta was not significant. In addition, we found a significant interaction effect of age and CSF and plasma NfL on brain-age delta, for which the CU females with larger brain-age delta had increased CSF and plasma NfL over the years, whereas the decelerated ones remained more stable. On the contrary, we found an effect of age interaction with plasma NfL on brain-age delta in MCI males, but not in MCI females. These results were expected, as females have higher chances of developing neurodegeneration and have showed to undergo faster cognitive decline than males (Ferretti et al., 2018). Although the role of sex hormones still needs to be clarified, it has been suggested that the menopausal drop of estrogen increases vulnerability to neurological events (Green & Simpkins, 2000; Maioli et al., 2021). On the contrary, results suggest that morphological effects of aging in the CU males’ brain are not fully driven by neurodegeneration, although these effects might increase with older age in MCI males. Lastly, regarding the cerebrovascular disease biomarkers, we found a positive association between brain-age delta and WMH for both CU females and males, while no interaction effect of sex was found. Conversely, in MCI individuals, we only found positive associations between brain-age delta and WMH on males, but no interaction effect of sex was found in MCI individuals. Overall, we found sex differences in the associations between brain-age delta and markers of neurodegeneration and cerebrovascular disease. NfL was only positively associated with brain-age delta in females and, although WMH were positively associated with both CU females and males, only MCI males showed this positive association. Positive associations between NfL and WMH have been previously demonstrated, for both CU and MCI (Osborn et al., 2018), and the different AD stages (Walsh et al., 2021). Moreover, it has been proposed that WMH may reflect two different pathological pathways, one including amyloid aggregation and another including axonal injury (Osborn et al., 2018). Our results may suggest that brain aging in males might be driven more strongly by the former pathway, while brain aging in females might be driven more by the latter one.
Our purpose was to study the clinical validity of using the brain-age delta as a proxy biomarker of brain aging associated to AD and neurodegeneration. Therefore, our main aim was studying the characteristics of the individuals whose brain age is more accelerated or decelerated. One of the strengths of this study was the robustness of the brain-age delta measurement using a widely used segmentation atlas such as the Desikan-Killiany. Notably, we demonstrated the robustness of our method by training our model with one cohort and testing independently on four independent cohorts. The similar results obtained in all cohorts allowed us to seek associations in a large sample of participants with biomarker data and to further stratify the data by sex. This aspect is critical for this type of analyses as the effects of biological aging are necessarily very small, particularly in CU individuals of limited age range. This may explain why we could not detect significant effects versus longitudinal brain atrophy, as the available sample size for these analyses were smaller since this variable was not available in all cohorts. Another strength of our work was that we were able to include a wide range of different biomarkers of AD pathology allowing us to perform an in-depth analysis of the effect of these measurements with the brain-age delta. Conversely, our model used a smaller number of features and a training set with a more limited age range than other models seen in literature recently (Liem et al., 2017; Peng et al., 2021), leading to a performance which cannot be compared against these state-of-the-art models. However, the performance of age prediction was similar to other publications that used similar methodologies (Beheshti et al., 2018; Dafflon et al., 2020; de Lange, et al., 2020) and, most importantly, was successful in studying the utility of brain-age delta as biomarker for AD and neurodegeneration. Future work should focus on developing a model with larger number of features or a 3D model and should study the effect of these validation measurements for AD and neurodegeneration with the brain-age delta more in depth.
In conclusion, we validated that machine-learning based brain age prediction obtained from a widely used segmentation atlas can be used as a biomarker of biological brain aging associated with AD pathology, risk factors and neurodegeneration. Moreover, our results confirm the presence of sex-related brain aging structural changes and suggest the prevalence of different neuropathological pathways involved in brain aging within females and males. Therefore, these results show the necessity to consider different approaches for assessing aging and neurodegeneration differently for each sex.
MATERIALS AND METHODS
Participants
We used a collection of T1-weighted brain MRI scans included in the UKBiobank (www.ukbiobank.ac.uk) dataset for training the proposed model and for calculating cross-validated brain age predictions. The dataset consisted of CU individuals (N = 22,661), after excluding subjects with ICD-9 and ICD-10 diagnosis, covering individuals of ages 44 to 81.
We also used four different cohorts to investigate the association between brain-age deltas with different sets of biomarker and AD risk factor measurements. Inclusion criteria for the independent cohorts consisted of: (i) availability of T1-weighted MRI brain scans; (ii) and availability of apolipoprotein E (APOE) categories and of CSF or PET measures for amyloid-β pathology acquired in less than a year from the MRI acquisition. These datasets included CU and MCI subjects from ADNI 1,2 and 3 (N = 751, CU = 253, MCI = 498), CU and MCI (as specified by a Clinical Dementia Rating = 0.5) subjects from the EPAD cohort (N = 808, CU = 653, MCI = 155), CU subjects from the ALFA+ cohort (N =380) and CU subjects from the OASIS cohort (N =407).
All the individuals had available data for the following clinical variables: chronological age, sex, MMSE and years of education, which will be referred as clinical variables from now on. A more detailed description of the clinical variables of these datasets is given in Table 1. Regarding AD-related variables, ALFA+, ADNI and EPAD cohorts included CSF Aβ42 measurements for categorizing Aβ pathology status, AT status determined by CSF Aβ42 and CSF p-tau, APOE categories and WMH. OASIS, meanwhile, only had data available for Aβ PET and APOE categories. In addition, ALFA+ and ADNI included biomarkers of neurodegeneration such as CSF NfL, plasma NfL and cortical atrophy measured by longitudinal changes in the so-called aging signature (Bakkour et al., 2013). The combination of available AD-related variables and neurodegeneration biomarkers will be referred as validation variables from so on. A more detailed description of the validation variables can be seen in Table 2.
Image Acquisition and Preprocessing
The UKBioBank, ADNI and OASIS datasets had available T1-weighted magnetic resonance (MR) images that had already been segmented with Freesurfer and had been parcellated using the FreeSurfer’s cortical Desikan-Killiany(Desikan et al., 2006) and subcortical aseg(Fischl et al., 2002) labeling pipelines, which had undergone a quality control procedure. Taking advantage of this available data, we decided to use the same segmentation pipeline with the ALFA+ and EPAD cohorts. All the image acquisition and preprocessing done is as follows.
The UKBiobank dataset consisted of T1-weighted magnetic resonance (MR) images, all collected using a 3T Siemens Skyra scanner and preprocessed as previously explained in more detail (https://biobank.ctsu.ox.ac.uk/crystal/crystal/docs/brain_mri.pdf). Images were previously segmented with Freesurfer 6.0 and underwent a quality control procedure.
For ADNI participants (Petersen et al., 2010), MRI acquisition methods are described in more detail elsewhere (http://adni.loni.usc.edu/methods/documents/). In brief, most of the T1-weighted MR were MP-RAGE, acquired with 1.5T or 3T scanners. Images were segmented with Freesurfer 5.1 and 6.0 and subjected to a quality control procedure. When possible, we also included a second T1-weighted MRI image sequence for the participants that underwent another MRI visit 3 years later. These scans were also segmented following the previously explained procedure.
For the OASIS subjects (Marcus et al., 2007), the MRI scans were acquired on 1.5 T or on 3.0 T scanners. T1-weighted magnetization-prepared rapid gradient echo (MP-RAGE) scans were obtained according to previously explained protocol (https://theunitedconsortium.com/wp-content/uploads/2021/07/OASIS-3_Imaging_Data_Dictionary_v1.8.pdf). All MRI sessions were segmented using FreeSurfer 5.1 or 5.3 and followed quality control measures. The PET images were acquired with [11C]PIB Pittsburgh’s compound 60-minute dynamic PET scan in 3D mode and the corresponding analysis analyses were performed using the PET unified pipeline (PUP, https://github.com/ysu001/PUP). Mean standardized uptake value (SUVR) values were converted to Centiloid scale as previously explained.
For the ALFA+ participants, a high-resolution 3D T1-weighted MRI sequence was acquired in a 3T Philips Ingenia CX scanner (TE/TR=4.6/9.9 ms, Flip Angle = 8º; voxel size= 0.75×0.75×0.75 mm). Images were segmented with Freesurfer 6.0 and subjected to a quality control procedure to identify and remove incidental findings (Brugulat-Serrat et al., 2017) and segmentation errors(Huguet et al., 2021). Some of these ALFA+ subjects (N=187) underwent a second MRI visit 3 years after the initial visit, where another T1-weighted MRI sequence was acquired and segmented following the same procedure as in the first visit.
For the EPAD cohort(Solomon et al., 2018), which is a multisite study, T1-weighted MRIs were inversion-recovery prepare 3D gradient-echo sequences, acquired with 3T scanners. Images were segmented with Freesurfer 6.0 and subjected to a quality control procedure (Lorenzini et al., 2021).
For all the cohorts, subsequent to the FreeSurfer segmentation, tissue regions were parcellated into 183 different anatomical regions of interest (ROI)s using the widely-used FreeSurfer’s cortical Desikan-Killiany(Desikan et al., 2006) and subcortical aseg(Fischl et al., 2002) labeling pipelines. As mentioned before, we used the available FreeSurfer segmentations from UKBioBank, ADNI and OASIS cohorts. All volumes were residualized with respect to total intracranial volume (TIV) and to scanning site, while all cortical thicknesses were residualized with respect to scanning site, using linear models. Lastly, we performed a standardization procedure by computing z-score measurements feature-wise within each cohort, as previously performed (Casamitjana et al., 2018; Subramaniapillai et al., 2021; Ten Kate et al., 2018). We then assessed that there were not statistical differences in mean cortical thickness and volumes between the cohorts (see Supplementary Fig. 4).
Biomarkers
CSF and plasma collection, processing and biomarkers measurements
CSF and blood collection, processing and storage in the ALFA+ study have been described previously (Milà-Alomà et al., 2020a; Suárez-Calvet et al., 2020). CSF p-tau181 was measured using the Elecsys® Phospho-Tau (181P) CSF electrochemiluminescence immunoassay on a fully automated cobas e 601 instrument (Roche Diagnostics International Ltd, Rotkreuz, Switzerland). CSF Aβ42 and NfL were measured with the NeuroToolKit on a cobas e 411 or cobas e 601 instrument (Roche Diagnostics International Ltd, Rotkreuz, Switzerland). Plasma NfL was measured using the commercial Quanterix® assay (Simoa® NF-light Kit cat. no. 103186) on a HD-X analyzer following the manufacturer’s instructions (Quanterix, Billerica, MA, USA). All these measurements were previously reported (Milà-Alomà et al., 2020; Suárez-Calvet et al., 2020). All measurements were performed at the Clinical Neurochemistry Laboratory, University of Gothenburg, Mölndal, Sweden, by laboratory technicians and scientists blinded to participants’ clinical information.
In the ADNI study, CSF samples were measured according to the kit manufacturer’s instructions and as described in previous studies (Bittner et al., 2016), using the Elecsys β-amyloid(1–42) CSF (Bittner et al., 2016), and the Elecsys Phospho-Tau (181P) and Elecsys Total-Tau CSF immunoassays on a cobas e 601 analyzer at the Biomarker Research Laboratory, University of Pennsylvania, USA. Plasma NfL was measured on an in-house immunoassay on the single-molecule array (Simoa) platform, using the same methodology as described previously, at the Clinical Neurochemistry Laboratory, University of Gothenburg, Mölndal, Sweden.
In the EPAD study, CSF was measured using the Elecsys β-amyloid (1–42) and the Elecsys Phospho-Tau (181P) CSF electrochemiluminescence immunoassay on a fully automated cobas e 601 instrument (Roche Diagnostics International Ltd.). All measurements were performed at the Clinical Neurochemistry Laboratory, University of Gothenburg, Mölndal, Sweden, by laboratory technicians and scientists blinded to participants’ clinical information. Concentrations of CSF Aβ42 and p-tau181 were determined according to the manufacturer’s instructions(Solomon et al., 2018).
Amyloid-β positivity cutoffs
For ALFA+, ADNI and EPAD, AT stages were defined by CSF Aβ42 and CSF p-tau, respectively. Previously used cut-offs were applied to each cohort, consisting of 1098 pg/mL for CSF Aβ42 for ALFA+ and EPAD (Schindler et al., 2018) and of 880 pg/mL for CSF Aβ42 for ADNI (Hansson et al., 2018) and of 24 pg/mL for p-tau for the three cohorts (Milà-Alomà et al., 2020). For OASIS, we used the cut-off value of 17 Centiloids from literature(Salvadó et al., 2019).
WMH Volumes
WMH volumes were generated for ALFA+ and EPAD cohorts using Bayesian Model Selection (BaMoS) procedure (Sudre et al., 2015), which has been provided previously. We also obtained the already available WMH volumes for ADNI cohort, in which the method of WMH volumetric quantification was performed using probabilistic models in a Markov Random Field framework, as previously provided (Schwarz et al., 2009). For each cohort, total WMH volumes were derived by summing and multiplying the number of labeled voxels by voxel dimensions. Total WMH volumes were natural log transformed and residualized with respect to TIV using linear models.
Aging signature measurements
For ALFA+ and ADNI, we computed the weighted Dickerson’s aging signature (Bakkour et al., 2013), which has been used as a proxy measurement for brain aging. The aging signature is a map of specific brain regions that undergo cortical atrophy in normal aging. This meta-ROI is composed of the surface-area weighted average of the mean cortical thickness in the following individual ROIs: calcarine, caudal fusiform, caudal insula, cuneus, inferior frontal gyrus, medial superior frontal and precentral cortices. A Z-score of this aging-specific measure was calculated based on the mean and standard deviation of the CU individuals, as done previously (Bakkour et al., 2013). This is referred as Aging Signature V1.
In addition, we also computed this measurement on the scanners from the second MRI visit, referred as Aging Signature V2. We then computed a longitudinal brain atrophy measurement by computing the aging signature change over the years between the MRI acquisitions. Therefore, longitudinal aging signature change was computed as: For another secondary analysis shown in Supplementary Table 9, we also computed the aging signature (Aging Signature V1) for the remaining independent cohorts: EPAD and OASIS.
Brain-Age Prediction
Regression Model
Model Workflow
For the current study, we used a gradient boosting framework: the XGBoost regressor model from the XGBoost python package (https://xgboost.readthedocs.io/en/latest/python) to run the brain age prediction. This regressor, which is based on a decision-tree based ensemble algorithm, was selected due to its speed and performance and its advanced regularization to reduce overfitting (Chen & Guestrin, 2016). In addition, large-scale brain age studies have demonstrated its adequacy (Bashyam et al., 2020; de Lange et al., 2019; de Lange, Barth, et al., 2020; Kaufmann et al., 2019). As it has been shown that there are sex-related trajectories in normal aging (Podgórski et al., 2021), we trained separate models for females and males. For each model, we first performed Bayesian parameter optimization based on a cross-validation scheme with ten folds and ten repeats per fold using the FreeSurfer volumes and thickness of the UK BioBank as input. For the optimization we used HYPEROPT (Bergstra et al., 2013), with which we scanned for maximum depth, number of estimators, learning rate, alpha regularization, lambda regularization, subsample, gamma and colsample by tree. The optimized parameters were maximum depth = 4, number of estimators = 800, learning rate = 0.03, alpha regularization = 4, lambda regularization =1, subsample= 0.36, gamma = 3 and colsample by tree = 0.89 for the males model; and maximum depth = 4, number of estimators = 850, learning rate = 0.03, alpha regularization = 8.5, lambda regularization =14.5, subsample= 0.449, gamma = 3.5 and colsample by tree = 0.72 for the females model. We trained these two models and performed the brain-age prediction on the independent cohorts. We decided to compute a ROI based model using these 183 FreeSurfer regions because they are widely used and available in most of the neuroimaging datasets. Therefore, our aim was not to compare our performance to the one achieved by a model trained with larger number of ROIs or with the full 3D images, but to study the generalizability and the relevance of our model in the AD field.
Contribution of Brain Regions in prediction
We computed SHAP (SHapley Additive exPlanation) values (https://github.com/slundberg/shap) to measure the contribution of each brain region in the prediction of age for each subject. SHAP assigns an importance value within the prediction to each feature (in this case, brain region), which is based on its unique consistent and locally accurate attribution (Lundberg et al., 2020). We calculated the average SHAP value for each region for all females and males of the UK BioBank cohort.
In addition, to assess that the regions with highest SHAP values were stable, we performed a permutation approach to study the significance of each region, separately for females and males. With this aim we compared the averaged SHAP value (region-specific) obtained when using the entire train set on the model to a null distribution calculated from 1,000 permutations performing subsample of the subjects, in which we trained and tested the model using 80% and 20% of the individuals, respectively.
Brain-Age Delta Estimation
We predicted brain-age on the independent cohorts separately: ALFA+, EPAD, ADNI and OASIS, using the previously trained model. To investigate the prediction performance, correlation analyses were run for predicted brain-age versus chronological age, R2, root mean square error (RMSE), and mean absolute error (MAE) were calculated for each independent cohort separately, as well as for females and males separate pooled from all independent cohorts. We also investigated the prediction performance on the UKBioBank cohort by computing the average latter metrics from a cross validation with ten splits and ten repetitions.
As recent research has shown that brain-age estimation involves a proportional bias (de Lange & Cole, 2020; Le et al., 2018; Liang et al., 2019; Smith et al., 2019), we applied a well-established age-bias correction procedure to our data (de Lange & Cole, 2020; Le et al., 2018). This correction, as originally proposed (de Lange & Cole, 2020; Le et al., 2018), consists of a linear regression between age (Ω) and brain-predicted age (Y) on each of the independent cohorts, Y = a × Ω + β. The derived values of slope (α) and intercept (β) from the training set were used to correct the predicted brain-age in each test set by applying: Corrected Predicted Brain Age = Predicted Brain Age + [Ω − (α × Ω + β)]. By subtracting the chronological age from the Corrected Predicted Brain Age, we obtained the brain-age delta which was used to test the associations with the validation measurements. The result of the correction is shown in Supplementary Fig.3.
Statistical Analyses
All statistical analyses were conducted using Python 3.7.0. We tested for normality of the distribution for each biomarker using the Kolmogorov-Smirnov test and visual inspection of histograms. CSF NfL and plasma NfL did not follow a normal distribution and were thus natural log transformed. In addition, to compare the measurements for CSF NfL and for plasma NfL coming from different cohorts (ALFA+ and ADNI), CSF and plasma NfL was converted to z-scores.
To study the performance and accuracy of the brain-age prediction for each cohort, correlation analyses were run for predicted brain-age versus chronological age. We also computed R2, RMSE and MAE for each cohort, as well as the age bias of the prediction after bias correction. We assessed statistically whether the accuracy of the predicted brain-age was different between cohorts by using Fisher’s z-transformation for correlation coefficients. In addition, we computed these performance metrics to assess the differences in the model for females and males in the pooled cohorts. Results from another secondary analysis are also shown in Supplementary Table 2, in which we assessed the performance and accuracy of the aging signature for all cohorts, by performing correlation analyses between the aging signature versus chronological age and computing the R2 and RMSE. In this secondary analyses, we also studied whether the performance obtained for the predicted brain-age was better than the aging signature by performing the William’s test (Williams, 1959) for the Pearson’s correlation coefficient and a F-test to assess which model was statistically better.
We used the brain-age delta as a measure of brain aging to study the associations between this measurement and the different AD and neurodegeneration biomarkers and risk factors. With this aim, we pooled all the subjects from all cohorts together and computed linear regression models for each validation variable, in which chronological age and sex were included as covariates. Local effect size of each of the continuous validation variables was calculated using Cohen’s f2 (Cohen, 2013). The mean brain age delta among Aβ pathology, AT stages and APOE status, were assessed by a one-way analysis of covariance (ANCOVA) adjusting for age and sex. Effect size of the different levels was calculated by dividing the estimated difference in the brain-age delta between the different categories by the estimated standard deviation. We also stratified the individuals by sex and studied the associations between brain-age delta and the different validation variables by computing linear regression models in which chronological age was included as covariate. We next tested for interactions between sex and the validation variables on brain-age delta using linear regression models and including chronological age as covariate. Lastly, we tested for interactions between age and the validation variables on brain-age delta for CU and MCI individuals.
We also studied the differences in volumes and cortical thickness between females and males in the UKBioBank for the brain regions that contributed the most to the prediction according to the SHAP values. With this aim we performed regression models for each ROI with sex as predictor variable, in which linear and quadratic expansions of age, site and TIV (only included for volume ROIs), were included as covariates.
As a secondary analysis we wanted to identify the individuals whose predicted brain-age deviate the most from chronological aging, i.e., individuals with the highest positive or lowest negative brain-age deltas, to study the above-mentioned associations. With this aim, we selected the individuals whose brain-age delta was included within the 10th and 90th percentile of the distribution for each independent cohort and studied the differences between these groups. The methodology and the results of this analysis can be found in Supplementary Appendix A.
Data Availability
UKBiobank data availability at www.ukbiobank.ac.uk ADNI data availability at www.adni.loni.usc.edu/data-samples/ EPAD data availability at www.ep-ad.org/ ALFA the datasets generated during and/or analysed during the current study may be available from the corresponding author on reasonable request.
CONFLICTS OF INTEREST
Ivonne Suridjan is a full-time employee and shareholder of Roche Diagnostics International Ltd; Gwendlyn Kollmorgen is a full-time employee of Roche Diagnostics GmbH; Anna Bayfield is a full-time employee and shareholder of Roche Diagnostics GmbH.
MSC has served as a consultant and at advisory boards for Roche Diagnostics International Ltd and has given lectures in symposia sponsored by Roche Diagnostics, S.L.U, Roche Farma, S.A and Roche Sistemas de Diagnósticos, Sociedade Unipessoal, Lda.
JL.M is currently a full-time employee of H. Lundbeck A/S and previously has served as a consultant or on advisory boards for the following for-profit companies or has given lectures in symposia sponsored by the following for-profit companies: Roche Diagnostics, Genentech, Novartis, Lundbeck, Oryzon, Biogen, Lilly, Janssen, Green Valley, MSD, Eisai, Alector, BioCross, GE Healthcare, and ProMIS Neurosciences.
ACKNOWLEDGEMENTS & FUNDING
This publication is part of the ALFA study (ALzheimer and FAmilies). The authors would like to express their most sincere gratitude to the ALFA project participants, without whom this research would have not been possible. Authors would like to thank Roche Diagnostics International Ltd. for kindly providing the kits for the CSF analysis of ALFA+ participants and GE Healthcare for kindly providing [18F]flutemetamol doses of ALFA+ participants. Collaborators of the ALFA Study are: Müge Akinci, Annabella Beteta, Alba Cañas, Irene Cumplido, Carme Deulofeu, Ruth Dominguez, Maria Emilio, Karine Fauria, Sherezade Fuentes, Oriol Grau-Rivera, Laura Hernandez, Gema Huesa, Jordi Huguet, Eider M Arenaza-Urquijo, Eva M Palacios, Paula Marne, Tania Menchón, Carolina Minguillon, Eleni Palpatzis, Cleofé Peña-Gómez, Albina Polo, Sandra Pradas, Blanca Rodríguez-Fernández, Aleix Sala-Vila, Gemma Salvadó, Mahnaz Shekari, Anna Soteras, Laura Stankeviciute, Marc Vilanova and Natalia Vilor-Tejedor.
The project leading to these results has received funding from “la Caixa” Foundation (ID 100010434), under agreement LCF/PR/GN17/50300004 and the Alzheimer’s Association and an international anonymous charity foundation through the TriBEKa Imaging Platform project (TriBEKa-17-519007). Additional support has been received from the Universities and Research Secretariat, Ministry of Business and Knowledge of the Catalan Government under the grant no. 2017-SGR-892 and the Spanish Research Agency (AEI) under project PID2020-116907RB-I00 of the call MCIN/ AEI /10.13039/501100011033. FB is supported by the NIHR biomedical research center at UCLH. MSC receives funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (Grant agreement No. 948677), the Instituto de Salud Carlos III (PI19/00155), and from a fellowship from “la Caixa” Foundation (ID 100010434) and from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No 847648 (LCF/BQ/PR21/11840004).
COBAS, COBAS E and ELECSYS are trademarks of Roche. All other product names and trademarks are the property of their respective owners. The Roche NeuroToolKit is a panel of exploratory prototype assays designed to robustly evaluate biomarkers associated with key pathologic events characteristic of AD and other neurological disorders, used for research purposes only and not approved for clinical use.