Introduction

Physical activity, defined as any behaviour involving bodily movement that results in energy expenditure and increased physical fitness [1, 2], and its subcategory of physical exercise, which is movement performed in a structured and repeated manner in order to improve physical fitness [1], have emerged as promising modifiable lifestyle factors that support healthy development and ageing [3, 4]. In the context of neuroscience and brain health, physical activity and exercise have been identified as key lifestyle factors that influence brain plasticity [5, 6], with accumulating evidence suggesting potential for delaying cognitive decline and reducing the risk of dementia[7–9].

A large body of observational and experimental studies reports beneficial effects of physical activity on brain structure and function across the lifespan [10, 11]. Evidence from animal models has identified several mechanisms through which physical activity may support brain health and cognitive function, including the increased production of neurotrophic factors such as BDNF, enhanced cellular signalling promoting neural plasticity, and neurogenesis [12, 13]. Neuroimaging research has demonstrated the associations of physical activity with white matter integrity [14, 15] and grey matter volume, predominantly in the (left) hippocampus [16, 17], but also in cerebellum [18], frontal, and parietal regions [19]. Several studies using machine learning algorithms to derive markers of brain ageing have also linked physical activity to younger-appearing brains [20–22]. Cognitive function has similarly been found to benefit from physical activity, particularly in domains such as executive function, memory, and processing speed [23, 3, 24].

Nevertheless, a recent large-scale population-based study using the UK Biobank cohort reported negligible effects of physical activity on cognitive function, potentially due to sampling bias (volunteer-based sample), younger sample age, or reverse causality [25]. Similarly, no significant relationship was observed between physical activity and the neuroimaging-derived biomarker of brain ageing, “brain age gap”, within the same cohort [26]. In light of these contradictory findings, the present study aimed to estimate the cross-sectional and longitudinal effects of physical activity on brain structure and cognitive function in the population-based LIFE-Adult cohort of individuals from Leipzig, Germany [27, 28]. Based on a large volume of research, we hypothesised that there would be a positive effect of physical activity (both self-reported and assessed by accelerometer) on brain structure and cognitive function. BrainAGE (Brain Age Gap Estimate), a machine learning derived marker of structural brain aging, and hippocampal volume were used as brain structural outcomes, and processing speed and executive function (derived from Trail Making Test A & B, respectively) as cognitive outcomes. Additional analyses were conducted to address the potential age-bias effect and a quasi-cross-lagged panel model was estimated to investigate the possibility of reverse causality.

Methods

Participants

The participants in this study were part of the LIFE-Adult Study, a large population-based cohort study, whose baseline assessment was conducted from August 2011 to November 2014 by the Leipzig Research Centre for Civilization Diseases (LIFE) [28]. The LIFE-Adult Study was approved by the ethical board of the Medical Faculty of the University of Leipzig. All the participants provided their written informed consent and received a small remuneration for their participation.

Ten thousand randomly selected participants from Leipzig, Germany, underwent an extensive data collection at baseline, consisting of physical examinations, anthropometric data collection, interviews on medical and familiar history, self-report sociodemographic and lifestyle questionnaires, as well as cognitive testing [28]. Subsequently, a subset of participants wore an accelerometer to assess physical activity and sleep patterns. Structural and functional neuroimaging data were acquired for a large proportion of participants at a separate visit several weeks after the baseline visit of the study. The follow-up data of 1799 participants were collected from June 2018 to August 2021 with as much overlap between the measures from the first time point as possible [27]. Accelerometer data were not collected at the follow-up.

In the present study, the analyses were conducted on a base-line sample of 2576 participants (age range: 19-82, Mdn = 64.58 years, IQR = 22.41 years, 47% women). Two separate longitudinal path analyses were performed on two (overlapping) smaller subsamples for which data on physical activity was available at both time points (sample with cognitive scores: n = 342, Mdn = 63.57 years, IQR = 27.38 years, 40.1% women; and sample with structural brain age: n = 253, Mdn = 62.88 years, IQR = 30.42 years, 39.5% women). See Figure 1 for details of the selection process.

Figure 1.

Measures

Exposure: Physical Activity

Physical activity at baseline was measured using a short version of the International Physical Activity Questionnaire (IPAQ-SH, [29]), a self-report measure of physical activity over the previous 7 days. Participants recorded how many days in the previous week they had spent walking or performing moderate or vigorous physical activities, as well as how many hours and minutes they had spent at each respective intensity level on each day. The evaluation was carried out based on the published evaluation guidelines (https://sites.google.com/view/ipaq/score), calculating the sum of the three reported activity areas: walking, moderate physical activity, and vigorous physical activity. The resulting measure is expressed in metabolic equivalent of task (MET) minutes per week [30]. In this study, we examined the effects of both total physical activity (the sum of all reported physical activities) and moderate to vigorous physical activity (the sum of moderate and vigorous physical activities), as previous studies have shown that moderate-to-vigorous, rather than total, PA drives changes in brain structure and/or function [31, 21, 32].

A subset of participants who were willing to participate in the sleep and activity part of the LIFE-Adult study also wore a 2D-axial accelerometer SenseWear Pro 3 (Bodymedia, Inc., Pittsburgh, Pennsylvania). Participants were asked to wear the accelerometer for seven consecutive days. Only the data from participants who had worn the accelerometer for at least four weekdays and one weekend day, and at most 14 days prior to the scanning session, were used in the analysis. The raw actigraphy data were processed using the SenseWear Professional software package version 6.1. The minute-by-minute data was subsequently exported to a Microsoft Excel template. The Excel template was used to customise the analysis windows according to the specific day-night-cycles of each participant (the standard SenseWear software uses a fixed time-window for activity analyses, e.g. 12 am to 11:59 pm). Based on the bedtime information provided in a sleep diary kept during the actigraphy recording, the respective night sleep intervals were identified manually. For the subsequent analyses, only daytime intervals with a minimum duration of 20 hours and an off-wrist duration of less than 15% of the daytime interval were used. Activity parameters were then recalculated for the daytime interval (defined as the period between two consecutive night sleep intervals). We used a measure of energy expenditure in MET-minutes that were individually averaged across the measuring days for the participants.

During the follow-up, physical exercise capacity rather than physical activity was measured using the Veteran Specific Activity Questionnaire (VSAQ, [33]). This questionnaire consists of a list of 13 activities arranged in order of increasing physical demand (e.g., from “Eating”, “Taking a shower” to “Any competitive activity”). Participants had to rate which activity would result in notable fatigue and heavy breathing. Their exercise capacity in MET was extrapolated based on the self-reported loading and age [33, 34].

Outcome: Brain structure

Brain structure was assessed from T1-weighted MRI scans, acquired on one 3T Siemens Verio scanner with a 32-channel coil. An MP-RAGE sequence with 1 mm isotropic voxels, 176 slices, TR=2300 ms, TE = 2.98 ms, and inversion time (TI) = 900 ms was used. The images were processed with the CAT12.9 [35] of SMP12 in Matlab 2024b. Two structural outcomes were devised: BrainAGE and hippocampal volume. The quality of the images was automatically tested using the CAT12 structural image quality rating (SIQR, [36]), with a cut-off threshold of 2.4, and by visual inspection, using the ‘Check Sample Homogeneity Tool’ in CAT12.

Brain age gap estimation (BrainAGE)

The Brain Age Gap Estimate (BrainAGE) represents a difference between a person’s chronological age and the age estimated by machine learning based on their brain scan. A positive or negative BrainAGE indicates that the brain appears older or younger than expected for the person’s age, respectively [37, 38]. The BrainAGE machine learning model (Gaussian process regression) previously described in Kalc et al. [39] was trained on a subset of scans from a large sample of healthy participants from various openly available neuroimaging datasets (n = 2446, M_Age = 57.69 ± 14.88 years, 58% women). It was then applied to baseline and follow-up brain images from the LIFE-Adult sample. As the brain age gap is always overestimated for younger subjects and underestimated for older ones, an age-bias correction was applied as described in Cole [26], resulting in a mean absolute error of 4.68 years and a correlation to age of 0.91. The age-bias correction is used to remove the dependency of the predicted brain age on age [40]. Therefore, we did not further control for the age in the regression analyses using the BrainAGE as the outcome. A detailed process of training and testing the algorithm is available in the Appendix.

Hippocampal volume

Hippocampal volume was extracted from the segmented grey matter regions using anatomical definitions provided by the Neuromorphometric atlas (Neuromorphometrics, Inc.). Specifically, the left and right hippocampal regions were identified and normalised by the total intracranial volume based on the analysis of covariance approach as described in Raz et al. [41]. The adjusted volumes were averaged into a single bilateral hippocampal volume measure.

Outcome: Cognitive function

Cognitive functioning of participants was assessed by a set of neuropsychological tests [28]. In this study, the time-to-complete scores from the Trail Making Test A and B (TMT, [42]) were available for both time points for the majority of the participants, and were therefore used as cognitive outcome measures. In TMT-A, participants had to connect the numbers on the paper in the correct order as quickly as possible. In part B, they had to connect numbers and letters of the alphabet, correctly alternating between both. The maximum time allowed for completion was 150 seconds and 300 seconds for part A and B, respectively. The time-to-complete scores (without accounting for errors) were used as a measure of processing speed and executive function of the participants (scores on TMT-A and TMT-B, respectively).

Covariates

Our covariate selection was based on the conceptual model of Campbell et al. [43], including sociodemographic factors and adult health behaviours. However, several possible covariates, such as factors that occurred prior to the birth or earlier in life, were not available in our dataset, therefore we opted for a smaller number. A minimal adjustment set of covariates was tested in the dagitty toolbox [44] (Figure S1 in the Appendix). Linear (and quadratic) polynomial terms of age, along with tobacco exposure (i.e., present and past smoker vs. non-smoker), alcohol consumption (g/day) and socioeconomic status (dummy coded) were included as covariates to account for confounding. Multicollinearity of the covariates was assessed using the formula for variance inflation factor: VIF = 1 / (1 - R²). The VIFs indicated no signs of multicollinearity, with all values being ≤ 1.757.

Statistical analysis

Statistical analyses were performed in R 4.4.3 [45], using RStudio version 2024.12.1+563. The cross-sectional multiple linear regression models were analysed with lavaan package version 0.6-19 [46], using a robust maximum likelihood estimation (MLR). Standard errors were calculated using a robust (Huber-White sandwich) estimator to better account for violations of normality. Missing data were handled by full information maximum likelihood (FIML) with the assumption that data was missing at random. We estimated the effect of self-reported total as well as moderate to vigorous physical activity as MET hours/day on two structural brain outcomes (i.e., BrainAGE and hippocampal volume) and two cognitive outcomes (i.e., attention and cognitive flexibility). As cognitive tests were performed before the accelerometer data collection, cognition was not included as an outcome in the analysis with accelerometer-derived physical activity measure. Two path analyses with a quasi-cross-lagged panel design were estimated for the subjects with available longitudinal physical capacity, BrainAGE, and cognitive scores to analyse the potential reverse causality.

Results

Sample characteristics for a cross-sectional baseline sample and its subsample with available actigraphy data, as well as for the two longitudinal subsamples with available cognitive scores (subsample 1) and structural MRI (subsample 2) are available in Table 1.

Table 1.

The following section presents the effects of physical activity on brain structural and cognitive outcomes in the cross-sectional analyses. We show the standardized estimates from unadjusted and adjusted models reflecting the effect of self-reported total PA on brain structure and cognitive outcome in Table 2. The effects of objectively measured PA on brain structure can be seen in Table 3. We present the effects of self-reported amount of MVPA on brain structure and cognition in Table 4.

Table 2.

Table 3.

Table 4.

The results indicate no statistically significant effects of self-reported total PA on brain structure. There is a statistically significant effect of total self-reported PA on cognitive function, indicating that higher levels of PA lead to higher time score on TMT-B. Similar results can be observed for MVPA, however, the results do not survive the correction for multiple comparisons.

The analysis of objectively measured PA indicated no statistically significant effects on brain structure.

The longitudinal effects of physical activity on cognition (executive function measured by TMT-B) and brain structural ageing (operationalised by BrainAGE) were estimated using path analysis. Self-reported total PA at baseline was used as a predictor of cognitive score or BrainAGE, and physical capacity at follow-up. The model testing the effects of PA on cognition showed no statistically significant effects of either PA at baseline on cognition at follow-up, nor of cognition at baseline on follow-up physical capacity. The model demonstrated acceptable fit to the data according to robust fit indices: χ² (1) = 2.21, p = .138, CFI = 0.995, TLI = 0.955, RMSEA [90% CI] = 0.077 [0.000, 0.22], SRMR = 0.026. In the model testing the effect of PA on brain structural ageing, we observed the effect of brain structural ageing at baseline on physical capacity at follow-up, but not the effect of PA at baseline on brain structural changes at follow-up. The robust fit indices of the model also suggested an acceptable fit (χ² (1) = 1.08, p = .300, CFI = 1.000, TLI = 0.998, RMSEA [90% CI] = 0.019 [0.000, 0.187], SRMR = 0.025). The path diagrams are shown in Figure 2.

Figure 2.

Sensitivity analyses

Given prior reports of lateralized effects of physical activity on left hippocampal volume [17], we conducted additional analyses focusing on this region only. However, no such effect was observed in our sample. The results are available in the supplementary table (Table S1).

Furthermore, due to the unexpected negative association between the hippocampal volume and self-reported PA, we examined the potential self-report bias. We compared the correlation coefficients between self-reported total and accelerometer-derived PA in the samples of subjects below and above the age of 60 using the package diffcor [47]. The age of 60 was chosen based on the characteristics of the sample, and to roughly divide between the working and retired population. There were significant differences in the correlation coefficients between the two groups (r_<60 = .52, r_>60 = .29, mean difference = 0.28, 95% CI [.002, .525]), showing systematic bias in self-reporting measures of PA in younger vs. older participants. We therefore conducted the analyses for the two age groups separately. The results are available in Tables S2 and S3. Total self-reported PA and MVPA had no effect on brain structure and cognition in either age group. However, when examining objectively measured PA, we observed a moderate effect of PA on the younger-appearing brain in the subset of adults aged between 30 and 60 years (n = 70), but we could not see the same effect in individuals over 60 years old (see Table S4 and S5).

Furthermore, we stratified the analyses based on biological sex, however, the effects were small and none survived the correction for multiple comparisons. The results are available in the Appendix in Tables S6–S11.

Discussion

A growing body of research shows the positive influence of physical activity on brain health and cognition [10, 11, 31]. In the present study, we investigated the effects of physical activity on brain structure and cognitive outcomes in a large, population-based cohort from the LIFE-Adult Study. We expected self-reported total and moderate to vigorous physical activity at baseline to have a positive effect on hippocampal volume and on making the brain appear younger (operationalised by smaller BrainAGE), as well as on better cognitive performance on the Trail Making Test. However, the results of multiple regression analyses did not support the hypothesised effect. Instead, the estimates showed the opposite direction (higher self-reported physical activity was negatively associated with hippocampal volume and positively with BrainAGE), although these effects were not statistically significant. There was a statistically significant negative effect of total self-reported physical activity on executive function, as measured by the TMT-B, showing an unexpected relationship between the exposure and outcome.

Self-reported measures of physical activity can be subject to reporting bias [48, 49]. A recent study examining sex-specific effects of self-reported physical activity on brain grey matter volume proposed that self-report measures of physical activity may not be sufficiently sensitive to investigate associations with brain health [50]. In line with this research, we found no effects of self-reported physical exercise on hippocampal volume or structural brain ageing. Furthermore, we did not find any effects in sex-stratified sensitivity analyses. However, we observed an age-related difference in the association between the self-reported physical activity measure and the objectively measured physical activity, suggesting a measurement error in the predictor, which probably affected our results.

To overcome the limitations of self-reported measures, we examined the effect of objectively measured PA by accelerometer on brain structure for a subset of participants. We observed no statistically significant effect of accelerometry-derived physical activity on either hippocampal volume (β = −0.04, p = 0.582) or BrainAGE (β = −0.076, p = 0.242). This is in contrast to some large-scale accelerometer studies, showing positive association of physical activity on brain structure [51, 52, 21]. However, the reported effects in those studies were small, so it is likely that our subsample lacked the statistical power to detect them. Furthermore, the subsample with accelerometer data in this study exhibited rather low activity levels [30], which may have fallen below the threshold of intensity and/or duration that the literature identifies as necessary to observe beneficial effects on brain health [53, 54]. Our sensitivity analyses have shown that in a subset of participants younger than 60 years, higher accelerometer-derived activity was associated with younger-appearing brains (lower BrainAGE), although these participants were on average more active than the older group and we could not see similar effects in a larger older group.

Contrary to a large body of research showing the positive influence of physical activity on brain structure and function, our results suggest no consistent effects in a large cohort of citizens from Leipzig, Germany. To date, only a small number of studies have reported no effects of physical activity on brain health or cognitive function (e.g., [50, 55–57]). Some of those studies were conducted on younger samples (at baseline), which may have reduced the impact of reverse causality. The majority of studies on older samples could have observed the effects of (prodromal) cognitive decline affecting physical activity rather than the other way around [25]. Similarly to Rodriguez-Ayllon et al. [58], who found a bidirectional association between physical activity and brain structure, with a more consistent pattern of brain structural measures affecting physical activity, the results of our path analysis partially supported the reverse causality explanation, indicating that baseline brain health influences follow-up physical capacity, rather than baseline physical activity affecting follow-up brain health. However, these results must be interpreted with caution, due to restrictions of the models and the possible attenuation of the effects by measurement error in the assessment of physical activity/capacity.

Strengths and limitations

The results of this cross-sectional observational study may be affected by various factors, including bias in self-reported physical activity [48, 49], accelerometer measurement error [59, 60] and reverse causality, among others. Moreover, our results may also be affected by the general medical status of the participants, since we did not control for other diseases within the sample. Although the longitudinal design is a major strength of this study, attrition of participants at follow-up may have affected our estimates. Furthermore, the use of cross-lagged panel model design in the longitudinal setting has frequently been criticised for not distinguishing between within-person changes and between-person differences [61, 62], and our adapted design suffers from these limitations, as well as others arising from the use of different instruments to measure the construct related to physical activity at each time point (IPAQ and VSAQ). Nevertheless, compared to large volunteer-based cohorts such as the UK Biobank, the registry-based recruitment strategy of the LIFE-Adult Study may be less susceptible to healthy volunteer bias. However, we cannot entirely eliminate the possibility of volunteer bias among the participants with accelerometry data in our case. Direct comparisons are limited in the absence of harmonised recruitment and assessment protocols. However, it is possible that the present study may potentially yield more generalisable estimates of the association between physical activity and brain structure and cognition for the white European population.

Conclusion

We found no consistent positive effects of physical activity on brain structure or cognition in the population-based LIFE-Adult cohort. Due to the aforementioned limitations, we cannot draw any definitive conclusions. However, the results are in line with those of observational longitudinal studies indicating minimal or no effect of physical activity on cognitive and/or brain health outcomes (e.g., [25, 56, 57]). Rather, they indicate a reverse pattern, whereby people with healthier brains tend to engage in more physical activities over time [63, 58]. Further longitudinal intervention studies with objectively measured physical activity in young and middle-aged adults are warranted to investigate the influence of physical activity on prospective brain and cognitive health.

Data Availability

We used data from the Leipzig Research Centre for Civilization Diseases (LIFE) under project agreement PV-801. All original data will be exclusively shared by LIFE (https://www.uniklinikum-leipzig.de/einrichtungen/life) based on a project proposal in accordance with the University of Leipzig’s data protection rules. The other publicly available neuroimaging datasets used in this study are accessible either openly or upon request: IXI (https://brain-development.org/ixi-dataset/), Cam-CAN (https://opendata.mrc-cbu.cam.ac.uk/projects/camcan/), Wayne State studies 10, 11, & EF datasets (https://fcon_1000.projects.nitrc.org/indi/retro/wayne_10.html, https://fcon_1000.projects.nitrc.org/indi/retro/wayne_11.html, https://fcon_1000.projects.nitrc.org/indi/retro/wayne_EF.html), the Enhanced NKI-Rockland sample (https://rocklandsample.org/for-researchers/step-3-accessing-the-data), Southwest University adult lifespan dataset (SALD; https://fcon_1000.projects.nitrc.org/indi/retro/sald.html), and Australian Imaging Biomarkers and Lifestyle study (AIBL; https://aibl.org.au/). The code used for this study is available on https://github.com/PolonaCa/PA_LIFE.

Acknowledgements

We would particularly like to thank the participants and researchers involved in the LIFE-Adult study, as well as other open neuroimaging databases that were used for this manuscript (IXI, Cam-CAN, Wayne studies, Enhanced NKI-Rockland sample, and Australian Imaging Biomarkers and Lifestyle flagship study of ageing - AIBL).

Additional information

Ethics Approval and consent to participate

All of the openly available datasets used in this manuscript originate from studies that have been approved by the relevant institutional ethics committees. Informed consent was obtained from all participants. The LIFE-Adult Study was approved by the ethical board of the Medical Faculty of the University of Leipzig. All the participants provided their written informed consent and received a small remuneration for their participation.

Funding

This work was supported by Carl Zeiss Stiftung as a part of the IMPULS project (IMPULS P2019-01-006), the Federal Ministry of Science and Education (BMBF) under the frame of ERA PerMed (Pattern-Cog ERAPERMED2021-127), and the Marie Skłodowska-Curie Innovative Training Network (SmartAge 859890 H2020-MSCA-ITN2019).

Author contributions

Polona Kalc: Conceptualization, Formal analysis, Writing - original draft, Writing - review and editing, Visualisation. Robert Dahnke: Software, Data curation, Visualisation, Writing - review and editing. Christian Sander: Resources, Data curation, Writing - review and editing. Frauke Beyer: Resources, Data curation, Writing - review and editing, Andrea Zülke: Writing - review and editing. Steffi G. Riedel-Heller: Resources, Data curation, Writing - review and editing. Veronica Witte: Resources, Data curation, Writing – review & editing. Christian Gaser: Funding acquisition, Supervision, Writing - review and editing.

Funding

Carl Zeiss Stiftung (IMPULS P2019-01-006)

• Christian Gaser

Bundesministerium für Forschung, Technologie und Raumfahrt (BMBF) (Pattern-Cog ERAPERMED2021-127)

• Christian Gaser

EC | Horizon Europe | Excellent Science | HORIZON EUROPE Marie Sklodowska-Curie Actions (MSCA)

https://doi.org/10.3030/859890

• Christian Gaser

Appendix 1

Directed acyclic graph (DAG)

We evaluated implied independencies in the DAG employing LOESS-based conditional independence testing with bootstrap resampling (500 iterations) in the dagitty R package [44]. The conditional independence tests confirmed DAG-implied independencies. No statistically significant associations were observed between age and alcohol consumption (estimate: 0.028, std. error: 0.020, 95% CI [−0.009, 0.067]), SES (estimate: −0.004, std. error: 0.024, 95% CI [−0.053, 0.043]), or smoking status (estimate: −0.022, std. error: 0.022, 95% CI [−0.070, 0.019]). A minimal adjustment set for the multiple linear regression was set to age, tobacco exposure, alcohol consumption, and SES.

Figure S1.

BrainAGE estimation

The supervised BrainAGE machine learning model uses a Gaussian Process Regression to learn the age-related pattern from affine-registered volumetric gray matter and white matter data and outputs the predicted brain age. The used model was trained on a subset of scans from a large collection of healthy participants from various openly available neuroimaging datasets. The scans from the following datasets were included: IXI (https://brain-development.org/ixi-dataset/), Cam-CAN [64], Wayne State studies 10, 11, & EF datasets (Wayne State University; https://fcon_1000.projects.nitrc.org/indi/retro/wayne_10.html, https://fcon_1000.projects.nitrc.org/indi/retro/wayne_11.html, https://fcon_1000.projects.nitrc.org/indi/retro/wayne_EF.html), the Enhanced NKI-Rockland sample [65], Southwest University adult lifespan dataset (SALD; https://fcon_1000.projects.nitrc.org/indi/retro/sald.html; [66]), and Australian Imaging Biomarkers and Lifestyle study (AIBL; [67]; https://aibl.org.au/). The subset for a training sample (n = 2446, M_Age = 57.69 ± 14.88 years, 58% women) was chosen to match the age and sex distribution of the participants in the LIFE-Adult study. The trained algorithm was then applied to a sample of 4297 brain images from the LIFE-Adult sample, resulting in a mean absolute error (MAE) of 5.97 years and a correlation to age of 0.91. The brain age gap is obtained by subtracting the chronological age from the predicted age. As the brain age gap is always overestimated for younger subjects and underestimated for older ones, an age-bias correction was applied as described in Cole [26]. Briefly, a random subset of the data was used to fit a linear regression between the chronological age as an independent and predicted age as a dependent variable. The predicted age in the LIFE-Adult dataset was then corrected by subtracting the intercept and dividing by the slope of the fitted model. This resulted in a MAE of 4.68 years. The age-bias correction is used to remove the dependency of the predicted age on age [40], we therefore did not further control for the age in the regression analyses using the BrainAGE as the outcome.

Table S1.

Table S2.

Table S3.

Table S4.

Table S5.

Table S6.

Table S7.

Table S8.

Table S9.

Table S10.

Table S11.