Skip to main content
medRxiv
  • Home
  • About
  • Submit
  • ALERTS / RSS
Advanced Search

Association between allostatic load and accelerated white matter brain aging: findings from the UK Biobank

View ORCID ProfileLi Feng, Zhenyao Ye, Zewen Du, Yezhi Pan, Travis Canida, Hongjie Ke, Song Liu, Shuo Chen, L. Elliot Hong, Peter Kochunov, Jie Chen, David K.Y. Lei, Edmond Shenassa, View ORCID ProfileTianzhou Ma
doi: https://doi.org/10.1101/2024.01.26.24301793
Li Feng
1Department of Nutrition and Food Science, College of Agriculture & Natural Resources, University of Maryland, College Park, Maryland, United States of America
2Department of Epidemiology and Biostatistics, School of Public Health, University of Maryland, College Park, Maryland, United States of America
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Li Feng
Zhenyao Ye
3Maryland Psychiatric Research Center, Department of Psychiatry, School of Medicine, University of Maryland, Baltimore, Maryland, United States of America
4Division of Biostatistics and Bioinformatics, Department of Epidemiology and Public Health, School of Medicine, University of Maryland, Baltimore, Maryland, United States of America
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Zewen Du
5Department of Biostatistics, School of Global Public Health, New York University, New York, New York, United States of America
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Yezhi Pan
3Maryland Psychiatric Research Center, Department of Psychiatry, School of Medicine, University of Maryland, Baltimore, Maryland, United States of America
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Travis Canida
6Department of Mathematics, The college of Computer, Mathematical, and Natural Sciences, University of Maryland, College Park, Maryland, United States of America
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Hongjie Ke
2Department of Epidemiology and Biostatistics, School of Public Health, University of Maryland, College Park, Maryland, United States of America
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Song Liu
7School of Computer Science and Technology, Qilu University of Technology (Shandong Academy of Sciences), Jinan, Shandong, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Shuo Chen
3Maryland Psychiatric Research Center, Department of Psychiatry, School of Medicine, University of Maryland, Baltimore, Maryland, United States of America
4Division of Biostatistics and Bioinformatics, Department of Epidemiology and Public Health, School of Medicine, University of Maryland, Baltimore, Maryland, United States of America
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
L. Elliot Hong
8Louis A. Faillace Department of Psychiatry & Behavioral Sciences, McGovern Medical School, University of Texas Health Science Center at Houston, Houston, Texas, United States of America
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Peter Kochunov
8Louis A. Faillace Department of Psychiatry & Behavioral Sciences, McGovern Medical School, University of Texas Health Science Center at Houston, Houston, Texas, United States of America
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jie Chen
9Department of Health Policy and Management, School of Public Health, University of Maryland, College Park, Maryland, United States of America
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
David K.Y. Lei
1Department of Nutrition and Food Science, College of Agriculture & Natural Resources, University of Maryland, College Park, Maryland, United States of America
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Edmond Shenassa
2Department of Epidemiology and Biostatistics, School of Public Health, University of Maryland, College Park, Maryland, United States of America
10Maternal & Child Health Program, School of Public Health, University of Maryland, College Park, Maryland, United States of America
11Department of Epidemiology, School of Public Health, Brown University, Rhode Island, United States of America
12Department of Epidemiology & Public Health, School of Medicine, University of Maryland, Baltimore, Maryland, United States of America
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: tma0929{at}umd.edu shenassa{at}umd.edu
Tianzhou Ma
2Department of Epidemiology and Biostatistics, School of Public Health, University of Maryland, College Park, Maryland, United States of America
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Tianzhou Ma
  • For correspondence: tma0929{at}umd.edu shenassa{at}umd.edu
  • Abstract
  • Full Text
  • Info/History
  • Metrics
  • Supplementary material
  • Data/Code
  • Preview PDF
Loading

Abstract

White matter (WM) brain age, a neuroimaging-derived biomarker indicating WM microstructural changes, helps predict dementia and neurodegenerative disorder risks. The cumulative effect of chronic stress on WM brain aging remains unknown. In this study, we assessed cumulative stress using a multi-system composite allostatic load (AL) index based on inflammatory, anthropometric, respiratory, lipidemia, and glucose metabolism measures, and investigated its association with WM brain age gap (BAG), computed from diffusion tensor imaging data using a machine learning model, among 22 951 European ancestries aged 40 to 69 (51.40% women) from UK Biobank. Linear regression, Mendelian randomization, along with inverse probability weighting and doubly robust methods, were used to evaluate the impact of AL on WM BAG adjusting for age, sex, socioeconomic, and lifestyle behaviors. We found increasing one AL score unit significantly increased WM BAG by 0.29 years in association analysis and by 0.33 years in Mendelian analysis. The age- and sex-stratified analysis showed consistent results among participants 45-54 and 55-64 years old, with no significant sex difference. This study demonstrated that higher chronic stress was significantly associated with accelerated brain aging, highlighting the importance of stress management in reducing dementia and neurodegenerative disease risks.

Introduction

Human brain aging is characterized by continuous microstructure transformations, including a reduction in cerebral integrity of grey and white matter compartments and cognitive decline(1, 2). Specifically, white matter microstructural integrity can be an early biomarker for the prodromal stage of dementia and other neurodegenerative disorders(3–5). We propose that allostasis, a dynamic physiological adaptation to chronic exposure to stressors(6), provides a working model to understand how chronic physical and psychosocial stressors may affect brain aging. In our investigation, we employed a machine learning algorithm(7–9), leveraging multiple fractional anisotropy tract measurements obtained from diffusion tensor imaging data to predict white matter (WM) Brain Age Gap (BAG), representing the difference between an individual’s chronologic and biological brain age. The mechanisms by which chronic stress impacts brain aging are not fully understood. Hence, we examined the association between chronic stress, quantified by Allostatic Load (AL), and biological brain aging, assessed through WM BAG using multiple statistical methods, from linear regression to potential causality inferring Mendelian Randomization, inverse probability weighting, and doubly robust methods.

Allostatic load was developed to quantify the “wear and tear” on the body, stemming from the continuous activation of stress responses via the hypothalamic–pituitary–adrenal axis and the sympathetic-adrenal-medullary system(10). A cross-sectional study of Lothian Birth Cohort (b. 1936) showed evidence of a positive association between AL and brain age by using T1-weighted MRI scans of gray matter and white matter among cognitively healthy older adults (mean age 73 years) residing in Edinburgh Scotland(11). Other investigations showed a similar positive association between elevated AL and white matter volume and integrity with limited sample size from 349 to 731 participants in mainly local cross-sectional cohorts(12–15) (See Table S1 for details). The selection of a brain aging indicator is essential. Recent studies highlight the efficacy of advanced machine learning techniques applied to neuroimaging data, offering the age-adjusted BAG estimates with minimal predictive biases(16, 17) proven to be linked to physical and brain health (9, 18–22). The single scalar metric BAG is a robust and easily interpretable outcome, simplifying the challenges of handling multiple variables and addressing dependencies typically encountered in complex multivariate brain imaging models(7, 8, 16).

To meet the challenge of establishing a link between AL and WM BAG while accounting for the multifactorial nature of WM BAG influenced by genetic, environmental, biological, and psychosocial factors(23–25), besides a linear regression model, we employed a Mendelian randomization method. Mendelian randomization method has been increasingly applied in observational studies by leveraging genetic variants as instrumental variables to establish a potential causal relationship between exposure and outcome(26). Compared to the previous methods, genetic variants are considered useful instrumental variables because they are less prone to be influenced by reverse causation or confounding factors, making them particularly strong for causal inference(27). Additionally, we considered a dichotomized AL risk and applied inverse probability weighting(28–30) and doubly robust methods to estimate the effect of AL on BAG for more robust inference(31).

In this study, we evaluated the association between AL and WM BAG among the general population as well as stratified by sex and age in the UK Biobank cohort. We hypothesized that an increase in chronic stress burden would accelerate white matter brain aging. The estimated effects derived from our study shed light on the long-term physiological toll of chronic stress on brain health. This research is a critical step for developing early interventions that can reduce the global burden of dementia by addressing chronic stress risk factors before they lead to irreversible brain damage.

Methods

Study population and variables

The data used in the present study were from the UK Biobank, a large-scale population-based longitudinal cohort study recruiting over 500 000 individuals aged 40 to 69 years in 22 recruitment centers across the UK at baseline from 2006-2010 and collecting comprehensive physical, environmental, genomic, and health data. Beginning in 2014, the brain imaging phenotypic data (N=∼40k) were collected(32). Imaging visit (instance 2, 2014+) where our outcome of interest WM BAG was measured happened six years after the baseline where our exposure the AL biomarkers were measured, providing the basis for evaluating the potentially causal relationship.

In this study, we focus on the non-pregnant, predominantly European white (self-reported) population in UK Biobank. We study white participants only to avoid bias of imbalanced data and reduce cross-population heterogeneity in training the brain age prediction model and we have shown in our previous studies that the white-only model estimates brain age well with high prediction accuracy (9, 21) The participants who had all five AL biomarkers at baseline (N = 348 361) were included to calculate the AL index. The participants who had fractional anisotropy data available but did not have extreme WM hyperintensities (N=30 375) were included to estimate the age-adjusted WM BAG in the testing dataset. Extreme WM hyperintensities can significantly alter brain structure and skew the brain age prediction(33). After merging, the final analytic sample includes N=22 951 participants for the analysis using linear regression, inverse probability weighting and doubly robust methods. For one-sample Mendelian Randomization (MR) analysis, we further merged the final analytical sample with genetic data and excluded participants with kinship, excessive heterozygosity, and sex discordance based on genetic data, resulting in an analytical sample of N=15 678 participants. These exclusion criteria are standard in Genome-Wide Association Studies (GWAS) to improve data reliability, maintain the integrity and accuracy of the genetic analysis(34, 35), and ensure more robust MR analyses. The details of the inclusion and exclusion criteria, data processing and analytical steps can be found in Figure 1.

Figure 1.
  • Download figure
  • Open in new tab
Figure 1. Flow chart of participants selection and data processing procedures.

Flowchart of our main analysis procedures and the number of subjects included at each step of the analysis. The main analyses are in three parts: (1) Estimation of age bias corrected white matter (WM) brain age gap (BAG) using machine learning model (Gray); (2) Linear association analysis between allostatic load (AL) index and BAG (Blue); (3) mendelian randomization (MR) method (green) to evaluate the potencial causal effect of AL on BAG; (4) inverse probability weighting (IPW) method (orange); (5) doubly robust (DR) method (pink).

Allostatic load (AL)

Our AL index is comprised of five biomarkers measured at baseline representing distinct physical systems: C-reactive protein (mg/L) representing immunological/inflammatory, waist-to-height-ratio representing anthropometric, peak expiratory flow (litres/min) representing respiratory, high-density lipoprotein cholesterol (mmol/L) representing lipidemia, glycated hemoglobin (mmol/mol) representing glucose metabolism. These five biomarkers were found to be available in most studies and predict mortality and longer, more elaborate AL batteries in recent meta-analysis(36).

Individuals at the highest quartile for C-reactive protein, waist-to-height-ratio, and glycated hemoglobin and the lowest quartile for high-density lipoprotein cholesterol and peak expiratory flow were marked as 1 and indicated as health risk among each 5-year and sex-specific distribution and otherwise as 0. Utilizing a combination of five distinct biomarkers representing various physiological systems, our AL index provides a comprehensive assessment of the cumulative effects of chronic stress on an individual’s health, offering a precise understanding of AL’s impact on brain aging(37). We further validated our proposed AL index by 1) testing its correlation with the subclinical health outcomes/measures (depressive score, self-reported health, left/right-hand grip strength, usual walking pace, and duration of walking for pleasure) as positive control and 2) testing its correlation with non-health-related variables (water intake each day, milk type used, usual side of head for mobile phone use) as negative control. Results of the additional validation tests can be found in Supplementary Table S2 and S3.

White matter brain age gap (WM BAG)

Our main outcome, the WM BAG, was computed based on the fractional anisotropy measures collected from diffusion tensor imaging data from the UK Biobank imaging assessment (instance 2, 2014+)(9). This study focuses on 39 WM fractional anisotropy tracks that cover multiple brain regions (see Table S4). Following the ENIGMA protocol, the average value for each white matter tract in the brain was computed using Track-Based Spatial Statistics analysis applied to the diffusion tensor imaging fractional anisotropy images(38). The fractional anisotropy measurement, which falls within the range of 0 to 1, reflects the level of anisotropy in a diffusion process and the integrity of the cortical white matter(39).

We applied a machine learning model (random forest regression) to predict brain age and estimate BAG from fractional anisotropy data. We followed a similar data splitting scheme in previous studies to split the data into a training set (N=8570) and a testing (N=30 375) (9, 40). The training set includes healthy participants only (excluding who were previous or current smokers, had hypertension, CVD, diabetes, or brain diseases, etc.). We employed random forest regression with fractional anisotropy features as predictors and chronological age as the outcome. The model’s parameters were optimized using a 5-fold cross-validation, focusing on maximizing the correlation between chronological age and estimated brain age, and minimizing mean absolute error (MAE). This approach also identified key fractional anisotropy features influencing brain aging. The optimized model was then applied to the testing set for brain age prediction. We calculated the WM BAG by subtracting the chronological age from the predicted brain age. To address age-dependent bias noted in previous studies(22, 41), we used linear regression for bias correction, assessing the correction method’s performance via MAE(42). The adjusted WM BAG from the testing set was the primary outcome in subsequent analyses.

Potential confounders and instrumental variables (IVs)

We included potential confounding factors from baseline in the UK Biobank, including sociodemographic factors: age, sex (females and males), education (A levels/AS levels or equivalent, College or University degree, CSEs or equivalent, NVQ or HND or HNC or equivalent, O levels/GCSEs or equivalent, and other professional qualifications), Townsend Deprivation Index (an area-based score of social deprivation (accounting for unemployment, overcrowding, non-car ownership, and non-home ownership)), household income (less than 18k, 18k-31k, 31k-52k, 52k-100k, more than 100k); lifestyle/behavior factors (e.g. diet, smoking, physical activity; all were calculated based on American Heart Association Life’s Essential 8’s definition, ranged 0-100 for each score (Higher score means healthier behavior))(43). We assumed missing at random for these covariates and performed multiple imputations by chained equations method to impute the covariates data(44).

For MR analysis, we further incorporated the genotype data to select the genetic variants as candidate instrumental variables (IVs). The genotype data were assayed with UK BiLEVE Axiom Array and with UK Biobank Axiom Array(45, 46). We performed quality control of genotype data and only kept the genetic variants with: minor allele frequency ≥ 0.01, Hardy-Weinberg equilibrium exact test P value ≥ 0.001, missing genotype rate ≤ 0.05 and missingness per individual ≤ 0.2. In the Mendelian randomization analysis, we further filtered out subjects who had missing genotypes for more than two valid SNPs and used random imputation to impute the missing genotypes(47).

Statistical analysis

Descriptive statistics

Variables of baseline characteristics are shown as % (N) if categorical, mean (SD) if continuous. The mean (SD) distribution of the AL index and WM BAG by baseline categorical variables was compared using an ANOVA test. Pearson correlation tests were applied to test the association between the continuous variables and the AL index and WM BAG. P and r values were shown for the significance of both tests and the Pearson correlation coefficient, respectively (Table 1). The baseline characteristics for the age- and sex-stratified groups are detailed in Table S5. Table S6 compares the baseline characteristics of participants included in the AL/BAG final analytic sample vs. those excluded due to pregnancy or extreme WMH; Table S7 compares the baseline characteristics of those included in the MR analytic sample vs. excluded due to kinship, excessive heterozygosity, and sex discordance.

View this table:
  • View inline
  • View popup
Table 1. Characteristics of Participants for linear regression, inverse probability weighting, and doubly robust estimation.

Association analysis between AL index and BAG

We used four statistical models to investigate the association between AL and WM BAG (Figure 2) among the general population as well as stratified by sex and age groups (early middle-aged (40–44), late middle-aged I (45–54), late middle-aged II (55–64), and early older adults (65–69)) in the UK Biobank cohort. We first applied a linear regression model to test the association between WM BAG and AL, controlling for the aforementioned confounders. We also tested for age and sex interactions with AL in both the linear regression and one-sample MR to determine whether these factors significantly influenced the association between AL and WM BAG.

Figure 2.
  • Download figure
  • Open in new tab
Figure 2.

Directed Acyclic Graphs to investigate the potential causal relationship between allostatic load (AL) and white matter (WM) brain age gap (BAG). Multiple models: model 1 (Multiple linear regression), model 2 (one-sample Mendelian Randomization), model 3 (inverse probability weighting method), model 4 (doubly robust method). Confounders: age, sex, education, household income, Townsend deprivation score, smoke score, diet score, and physical activity score.

Association analysis using linear regression models does not directly infer causality. We then performed a one-sample MR analysis by treating genetic variants as IVs to evaluate the potential causal relationship between AL and WM BAG following the guidelines by Burgess et al. (2023)(48). First, from a larger non-overlapping sample (N=217 606) with valid AL index but no fractional anisotropy data in the UK Biobank, we performed a GWAS (35) (P < 5X 10-8) adjusting for population stratification followed by linkage disequilibrium clumping (r2=0.5 within window 250 kb) using Plink (version 1.9, www.cog-genomics.org/plink/1.9/) to select the variants associated with AL index as candidate IVs (IV assumption (i))(49). In one-sample mendelian randomization analysis, we further removed variants significantly associated with the covariates (false discovery rate adjusted P >0.05; IV assumption (ii)) and then removed the horizontal pleiotropic variants by regressing BAG on AL index with the genetic variables (false discovery rate adjusted P >0.1; IV assumption (iii)). We annotated the final selected IVs using ANNOVAR (https://annovar.openbioinformatics.org/en/latest/) and used R packages “gwascat” to find relevant traits associated with these IVs. We also validated the IV selection by testing the IV effect on AL using the first stage F-statistics and the IV x age or IV x sex interaction effects. We mainly used the inverse-variance weighted (IVW) MR method and other popular MR methods (simple median, weighted median, penalized weighted median, inverse variance weighted, Egger, etc.) were also applied for sensitivity analysis.

To further verify our findings, we split the participants into AL high-risk (AL index >2; N=6208) and low-risk group (AL index <2; N=16 743) and estimated the effect of AL risk on BAG using the inverse probability weighting (IPW) and doubly robust (DR) estimation methods. We also conducted sensitivity analyses by removing the outlier participants (N=1223) with extreme propensity scores, or redefining the high-risk AL group (AL index > 1, 3, 4, as representing the top 59.28%, 10.68%, 3.14% of the population, respectively) and rerunning IPW method for robustness.

All statistical analyses were conducted using R (version 4.0.5)(50). R packages, including “Hmisc” (version 5.1.0), “MICE” (version 3.14.0), “ipw” (version 1.2), “MendelianRandomization” (version 0.5.1)(51) were used to perform multiple imputation by chained equations, inverse probability weighting, doubly robust, and Mendelian randomization analyses. Except for IV selection, P<0.05 was regarded as statistically significant for all other analyses unless otherwise specified. For IV selection in MR analyses, Benjamini-Hochberg method was used to adjust for multiple comparisons(52).

Results

The optimal random forest regression model selected 16 fractional anisotropy measures to predict brain age (Table S2). The age-adjusted predicted brain age achieved a Pearson correlation coefficient (R)= 0.89, mean absolute error (MAE) = 2.74 years in the testing dataset (Figure S1), implying an excellent prediction performance of the proposed machine learning model. In the final analytical sample of 22 951 participants (mean age 55.37, 51.40% women), the mean AL index (ranging from 0-5) is 1.01 (SD, 1.09), and the mean of WM BAG is 0.13 (SD, 3.45).

Participants with a lower AL index were older, had higher education, more household income, lower Townsend deprivation index, and higher diet, smoke, and physical activity scores, but no sex difference (Table 1). We also observed that older groups, especially older females, had lower WM BAG, AL index, Townsend deprivation index, and higher diet, smoke, and physical activity scores (Table 1 and Table S5). The baseline characteristics are largely alike between the included and the excluded participants (Table S6-7) and no clear selection bias has been observed in our study.

Overall, all the four models showed significant positive relationship between AL and WM BAG (Table 2). In the linear regression model, when adjusting for age, sex, education, household income, Townsend Deprivation Index, smoking, diet, and physical activity, WM BAG is estimated to increase by 0.29 years (106 days, 95% CI: 0.25 to 0.33 years, P<0.001) per one AL score increment.

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 2. Results of 1) association analysis by linear regression model between continuous allostatic load (AL) index and white matter (WM) brain age gap (BAG), 2) one-sample Mendelian randomization (MR) analysis between continuous AL index and WM BAG, 3) inverse probability weighting (IPW) between High/low-risk AL groups and WM BAG, 4) doubly robust (DR) between High/low-risk AL groups and WM BAG among overall participants. All models were adjusted for age, sex, education, Townsend deprivation index, household income, smoke score, diet score, and physical score.

In one-sample MR model, the GWAS analysis identified 2636 genetic variants significantly (P< 5 X10-8) associated with AL index (see Manhattan plot in Figure S2). A final set of 238 variants were identified as valid IVs after checking the IV assumptions (Table S8). All these IVs had first stage F-statistics above the acceptable threshold and showed no significant interactions with age or sex (Table S9). Many of these IVs were found to be associated with biomarkers including C-reactive protein, hemoglobin and high-density lipoprotein cholesterol collected in GWAS Catalog(53), aligning well with the five biomarkers used to define our multi-system AL index (Table S9). Using these IVs, we found an overall significant positive effect of AL index (0.33 years/120 days, 95% CI=0.01 to 0.655 years per one AL score increment, P=0.042) on WM BAG (Table 2).

For dichotomized AL risk group comparison, both the IPW and DR methods showed that participants with high-risk AL were older in WM BAG than those with low-risk AL (IPW: 0.60 years/219 days, 95% CI: 0.51 to 0.69 years, P <0.001; DR: 0.60 years/219 days, 95% CI: 0.50 to 0.70 years, P <0.001), confirming a consistent finding of the positive relationship between AL and WM BAG (Table 2).

The results of the age- and sex-stratified analyses can be found in Figure 3, Table S10. No significant interactions between AL and age or sex were observed (Table S11). In the age-stratified analysis, the positive association between AL index and WM BAG was consistent in all four age strata (early middle age, late middle age I, late middle age II, late adulthood) from linear regression and inverse probability weighting method but not one-sample Mendelian randomization method, while we observed that the late middle-aged participants (45-54 years, 55-64 years) had the most significant effects. Sex showed no significant difference overall, but some sex disparities were found in different age strata.

Figure 3.
  • Download figure
  • Open in new tab
Figure 3.

(Left) Associations between adjusted white matter (WM) BAG (brain age gap) and allostatic load (AL) index (continuous) by sex and age group (Gray area, 95% CI). (Right) inverse probability weighting (IPW)-handled association analyses between WM BAG and High-risk AL (binary) by sex and age group. The estimate represents the absolute WM BAG difference between High-risk AL and Low-risk AL.

The sensitivity analysis results from other MR methods matched well with that by the main IVW method (Table S12). The results were not materially changed when removing extreme propensity scores (Figure S3, Table S13), or redefining high-risk vs. low-risk AL groups (Table S14) for IPW method.

Discussion

This study reveals a consistently positive association and potentially causal relationship between AL and WM BAG among the participants of White ethnicity aged 40-69 in the UK Biobank across a range of statistical models, with the strongest effect observed among participants aged 45-64 years old. To the best of our knowledge, this study marks the first exploration of the potential causal relationship between AL and brain aging in a large community-dwelling population, encompassing both the overall population and age- and sex-specific groups. Our proposed AL index represents a universal and consensus set of biomarkers from multiple systems that capture the physiological ‘wear and tear’ and were further validated by using positive and negative controls. It is noteworthy that our machine learning-based BAG estimation exhibits remarkable predictive accuracy. These results illuminate the positive relationship between chronic stress and brain aging, offering insights into potential strategies for promoting healthy brain aging and mitigating dementia and neurodegenerative diseases burden.

Previous studies(12–15, 54, 55) have identified a negative association between cortisol levels (i.e., serum, hair) and brain volume, changes in the microstructure of gray or white matter, and impaired cognitive functions. For chronic stress, our results are consistent with evidence indicating a positive link between chronic stress (AL index) and WM volume(12, 13, 54), WM tracks(15, 55), and WM microstructure(14). In contrast to earlier studies that used multiple singular measurements, WM BAG offers a holistic estimation of brain aging that is derived from 39 brain-wide white matter fractional anisotropy tracks using a machine learning model, providing a single measure of brain aging that is predictive of cognitive decline(23) and standardized for use across diverse studies(56).

The complex mechanism underlying the impact of chronic stress on brain aging and the microstructure of white matter remains an active area of research. At the core of this relationship is the dysregulation of the hypothalamic-pituitary-adrenal axis, a primary characteristic of an impaired adaptive stress response commonly observed in brain aging(57). Chronic stress precipitates a cascade of endocrine disruptions, causing sustained elevation of cortisol—a glucocorticoid linked to deleterious effects on cerebral architecture and functionality, manifesting as compromised white matter integrity(58, 59). This may be partly due to disrupted myelination process of oligodendrogenesis that maintained repair white matter. Chronic stress targets the prefrontal cortex, amygdala, and hippocampus(60), and can alter oligodendrogenesis and thus the balance between myelin damage and repair(61, 62). Structural changes in white matter can compromise its fundamental role in facilitating signal transmission between brain regions(63), hence indicating the changes in white matter may be a mediator between chronic stress and cognitive function and emotion regulation. The relation between chronic stress and persistent systemic inflammation, alongside an increase in oxidative stress, is potentially an additional pathway contributing to the brain’s vascular alterations, neurodegeneration(64), and accelerated aging of the white matter.

Collectively, these findings suggest that chronic stress acts through both hormonal dysregulation and inflammatory pathways may elicit changes that are detrimental to the brain’s white matter, potentially hastening the natural aging process. This study has provided empirical evidence that AL contributes to WM BAG. However, the precise mechanisms by which AL influences the biological brain age remains to be further elucidated, which is crucial for developing targeted interventions.

This study exhibits several strengths in design and analysis compared to other assessments of the impact of AL on brain aging. Firstly, we leveraged a large-scale UK Biobank cohort and conducted various statistical analyses, enabling us to unravel a potential causal relationship between chronic stress and brain aging across overall and different age and sex groups. Secondly, informed from the meta-analysis(36), we incorporated five robust indicators into one AL index, representing vital physiological systems validated from positive and negative controls (Supplementary Tables S3 and S4). Thirdly, we employed an optimal machine learning model to calculate the WM BAG and adjusted age bias, resulting in an accurate estimation of biological brain age. Fourthly, we utilized multiple robust causality inference methods, including Mendelian randomization, inverse probability weighting and doubly robust approaches, which mitigated the impact of confounding variables and reverse causation, thus enhancing the validity of the study. Lastly, we conducted multiple sensitivity analyses on AL measurements and High-risk AL definitions, as well as alternative Mendelian randomization analyses, ensuring the reliability and validity of our findings.

The study has certain limitations that offer opportunities for future research. First, the UK Biobank comprises a relatively healthy population, which may skew the sample towards healthier individuals and potentially introduce selection bias(65). The general population may include a larger proportion of subjects more vulnerable to the damaging effects of AL (e.g. those with mental illness(37)). This could compromise the generalizability of our results and potentially lead to an underestimation of the association between allostatic load and brain aging. Second, the construct of AL index in this study was limited to biomarkers available through the UK Biobank. To date, over 50 clinical and non-clinical biomarkers were found to define AL index before(36), yet no consensus has been reached regarding the optimal set of biomarkers. Future research could leverage high-throughput omics technologies to explore the potentially common molecular underpinnings and mechanism of chronic stress in the body, mitigating the controversy surrounding AL definition. Lastly, by design of UKB cohort, we cannot further track each participant for their change in chronic stress and brain aging over time, ongoing longitudinal cohorts such as Adolescent Brain Cognitive Development study and Baltimore Longitudinal Study of Aging are to be investigated in the future study to justify the causal relationship(66, 67).

Conclusion

In conclusion, in this pioneering examination of the association and potential causal relationship between chronic stress and brain aging using a robust set of statistical methods, we observed a positive relationship between allostatic load and WM BAG among the participants of White ethnicity aged 40-69 in the UK Biobank, which highlights the critical role of chronic stress management in preventing early-onset cognitive decline and advancing public health strategies against dementia and neurodegenerative diseases.

Author Contributions

L.F. took the lead in performing the analysis and wrote the manuscript. Z.Y. and Z.D. performed the one-sample Mendelian randomization analysis. T.M. and E.S. developed the idea, supervised the project and took the lead in editing the manuscript. Y.P., T.C., H.K.,S.C., E.H., P.K., J.C., D.L. contributed to manuscript writing and polishing. All authors provided critical feedback and helped to shape the research, analysis, and manuscript.

Source of Funding

This work was supported by the University of Maryland Grand Challenge grant, University of Maryland MPower Brain Health and Human Performance seed grant, University of Maryland Department of Epidemiology and Biostatistics pilot award and National Institutes of Health under award numbers 1DP1DA048968, R01MH123163, R01EB015611, and S10OD023696.

Disclosures

None.

Data Availability Statement

The raw genetic and phenotypic data used for this study can be found in the UK Biobank (http://www.ukbiobank.ac.uk/)

Footnotes

  • We have added some sensitivity analyses and some explanations in detail. We have also minimized the "causal" word in the manuscript.

References

  1. 1.↵
    Coelho A, Fernandes HM, Magalhaes R, et al. Signatures of white-matter microstructure degradation during aging and its association with cognitive status. Sci Rep 2021;11(1):4517.
    OpenUrlCrossRefPubMed
  2. 2.↵
    Power MC, Su D, Wu A, et al. Association of white matter microstructural integrity with cognition and dementia. Neurobiol Aging 2019;83:63–72.
    OpenUrlCrossRef
  3. 3.↵
    Kantarci K, Schwarz CG, Reid RI, et al. White matter integrity determined with diffusion tensor imaging in older adults without dementia: influence of amyloid load and neurodegeneration. JAMA Neurol 2014;71(12):1547–54.
    OpenUrl
  4. 4.↵
    Gold BT, Johnson NF, Powell DK, et al. White matter integrity and vulnerability to Alzheimer’s disease: preliminary findings and future directions. Biochim Biophys Acta 2012;1822(3):416–22.
    OpenUrlCrossRefPubMedWeb of Science
  5. 5.↵
    Xiao D, Wang K, Theriault L, et al. White matter integrity and key structures affected in Alzheimer’s disease characterized by diffusion tensor imaging. Eur J Neurosci 2022;56(8):5319–31.
    OpenUrl
  6. 6.↵
    McEwen BS. Allostasis and allostatic load: implications for neuropsychopharmacology. Neuropsychopharmacology 2000;22(2):108–24.
    OpenUrlCrossRefPubMedWeb of Science
  7. 7.↵
    Cole JH, Franke K. Predicting Age Using Neuroimaging: Innovative Brain Ageing Biomarkers. Trends Neurosci 2017;40(12):681–90.
    OpenUrlCrossRefPubMed
  8. 8.↵
    Cole JH, Marioni RE, Harris SE, et al. Brain age and other bodily ‘ages’: implications for neuropsychiatry. Mol Psychiatry 2019;24(2):266–81.
    OpenUrl
  9. 9.↵
    Feng L, Ye Z, Mo C, et al. Elevated blood pressure accelerates white matter brain aging among late middle-aged women: a Mendelian Randomization study in the UK Biobank. J Hypertens 2023.
  10. 10.↵
    McEwen BS. Neurobiological and Systemic Effects of Chronic Stress. Chronic Stress (Thousand Oaks) 2017;1.
  11. 11.↵
    Cole JH, Ritchie SJ, Bastin ME, et al. Brain age predicts mortality. Mol Psychiatry 2018;23(5):1385–92.
    OpenUrlPubMed
  12. 12.↵
    Booth T, Royle NA, Corley J, et al. Association of allostatic load with brain structure and cognitive ability in later life. Neurobiol Aging 2015;36(3):1390–9.
    OpenUrlCrossRefPubMed
  13. 13.↵
    Ritchie SJ, Tucker-Drob EM, Cox SR, et al. Risk and protective factors for structural brain ageing in the eighth decade of life. Brain Struct Funct 2017;222(8):3477–90.
    OpenUrlCrossRef
  14. 14.↵
    Zsoldos E, Filippini N, Mahmood A, et al. Allostatic load as a predictor of grey matter volume and white matter integrity in old age: The Whitehall II MRI study. Sci Rep 2018;8(1):6411.
    OpenUrl
  15. 15.↵
    Ottino-Gonzalez J, Jurado MA, Garcia-Garcia I, et al. Allostatic load and disordered white matter microstructure in overweight adults. Sci Rep 2018;8(1):15898.
    OpenUrl
  16. 16.↵
    Cole JH, Poudel RPK, Tsagkrasoulis D, et al. Predicting brain age with deep learning from raw imaging data results in a reliable and heritable biomarker. NeuroImage 2017;163:115–24.
    OpenUrlCrossRefPubMed
  17. 17.↵
    Smith SM, Vidaurre D, Alfaro-Almagro F, et al. Estimation of brain age delta from brain imaging. Neuroimage 2019;200:528–39.
    OpenUrl
  18. 18.↵
    Hogan AM, Vargha-Khadem F, Saunders DE, et al. Impact of frontal white matter lesions on performance monitoring: ERP evidence for cortical disconnection. Brain 2006;129(8):2177–88.
    OpenUrlCrossRefPubMedWeb of Science
  19. 19.
    Gaser C, Franke K, Kloppel S, et al. BrainAGE in Mild Cognitive Impaired Patients: Predicting the Conversion to Alzheimer’s Disease. PLoS One 2013;8(6):e67346.
    OpenUrlCrossRefPubMed
  20. 20.
    Kolbeinsson A, Filippi S, Panagakis Y, et al. Accelerated MRI-predicted brain ageing and its associations with cardiometabolic and brain disorders. Scientific Reports 2020;10(1):1–9.
    OpenUrl
  21. 21.↵
    Ye Z, Mo C, Liu S, et al. Deciphering the causal relationship between blood pressure and regional white matter integrity: A two-sample Mendelian randomization study. J Neurosci Res 2023.
  22. 22.↵
    Mo C, Wang J, Ye Z, et al. Evaluating the causal effect of tobacco smoking on white matter brain aging: a two-sample Mendelian randomization analysis in UK Biobank. Addiction 2023;118(4):739–49.
    OpenUrl
  23. 23.↵
    Franke K, Gaser C. Ten Years of BrainAGE as a Neuroimaging Biomarker of Brain Aging: What Insights Have We Gained? Front Neurol 2019;10:789.
  24. 24.
    Bocklandt S, Lin W, Sehl ME, et al. Epigenetic predictor of age. PLoS One 2011;6(6):e14821.
    OpenUrlCrossRefPubMed
  25. 25.↵
    Rando TA, Chang HY. Aging, rejuvenation, and epigenetic reprogramming: resetting the aging clock. Cell 2012;148(1-2):46–57.
    OpenUrlCrossRefPubMedWeb of Science
  26. 26.↵
    Gupta V, Walia GK, Sachdeva MP. ’Mendelian randomization’: an approach for exploring causal relations in epidemiology. Public Health 2017;145:113–9.
    OpenUrlPubMed
  27. 27.↵
    Davies NM, Holmes MV, Davey Smith G. Reading Mendelian randomisation studies: a guide, glossary, and checklist for clinicians. BMJ 2018;362:k601.
  28. 28.↵
    Shiba K, Kawahara T. Using Propensity Scores for Causal Inference: Pitfalls and Tips. J Epidemiol 2021;31(8):457–63.
    OpenUrlCrossRefPubMed
  29. 29.
    Robins JM, Hernan MA, Brumback B. Marginal structural models and causal inference in epidemiology. Epidemiology 2000;11(5):550–60.
    OpenUrlCrossRefPubMedWeb of Science
  30. 30.↵
    Chesnaye NC, Stel VS, Tripepi G, et al. An introduction to inverse probability of treatment weighting in observational research. Clin Kidney J 2022;15(1):14–20.
    OpenUrlCrossRefPubMed
  31. 31.↵
    Funk MJ, Westreich D, Wiesen C, et al. Doubly robust estimation of causal effects. Am J Epidemiol 2011;173(7):761–7.
    OpenUrlCrossRefPubMedWeb of Science
  32. 32.↵
    Sudlow C, Gallacher J, Allen N, et al. UK biobank: an open access resource for identifying the causes of a wide range of complex diseases of middle and old age. PLoS Med 2015;12(3):e1001779.
    OpenUrlCrossRefPubMed
  33. 33.↵
    Maniega SM, Valdes Hernandez MC, Clayden JD, et al. White matter hyperintensities and normal-appearing white matter integrity in the aging brain. Neurobiol Aging 2015;36(2):909–18.
    OpenUrlCrossRefPubMed
  34. 34.↵
    Barsh GS, Copenhaver GP, Gibson G, et al. Guidelines for genome-wide association studies. PLoS Genet 2012;8(7):e1002812.
    OpenUrlCrossRefPubMed
  35. 35.↵
    Uffelmann E, Huang, Q.Q., Munung, N.S., et al. Genome-wide association studies. Nature Reviews Methods Primers 2021.
  36. 36.↵
    McCrory C, McLoughlin S, Layte R, et al. Towards a consensus definition of allostatic load: a multi-cohort, multi-system, multi-biomarker individual participant data (IPD) meta-analysis. Psychoneuroendocrinology 2023;153:106117.
  37. 37.↵
    Lenart-Bugla M, Szczesniak D, Bugla B, et al. The association between allostatic load and brain: A systematic review. Psychoneuroendocrinology 2022;145:105917.
  38. 38.↵
    Smith SM, Jenkinson M, Johansen-Berg H, et al. Tract-based spatial statistics: voxelwise analysis of multi-subject diffusion data. Neuroimage 2006;31(4):1487–505.
    OpenUrlCrossRefPubMedWeb of Science
  39. 39.↵
    Timpe JC, Rowe KC, Matsui J, et al. White matter integrity, as measured by diffusion tensor imaging, distinguishes between impaired and unimpaired older adult decision-makers: A preliminary investigation. J Cogn Psychol (Hove) 2011;23(6):760–7.
    OpenUrl
  40. 40.↵
    Mo C, Wang J, Ye Z, et al. Evaluating the causal effect of tobacco smoking on white matter brain aging: a two-sample Mendelian randomization analysis in UK Biobank. Addiction 2022.
  41. 41.↵
    Le TT, Kuplicki RT, McKinney BA, et al. A Nonlinear Simulation Framework Supports Adjusting for Age When Analyzing BrainAGE. Front Aging Neurosci 2018;10:317.
  42. 42.↵
    Beheshti I, Nugent S, Potvin O, et al. Bias-adjustment in neuroimaging-based brain age frameworks: A robust scheme. Neuroimage Clin 2019;24:102063.
  43. 43.↵
    Lloyd-Jones DM, Allen NB, Anderson CAM, et al. Life’s Essential 8: Updating and Enhancing the American Heart Association’s Construct of Cardiovascular Health: A Presidential Advisory From the American Heart Association. Circulation 2022;146(5):e18–e43.
    OpenUrlCrossRefPubMed
  44. 44.↵
    Azur MJ, Stuart EA, Frangakis C, et al. Multiple imputation by chained equations: what is it and how does it work? Int J Methods Psychiatr Res 2011;20(1):40–9.
    OpenUrlCrossRefPubMed
  45. 45.↵
    Wain LV, Shrine N, Miller S, et al. Novel insights into the genetics of smoking behaviour, lung function, and chronic obstructive pulmonary disease (UK BiLEVE): a genetic association study in UK Biobank. The Lancet Respiratory Medicine 2015;3(10):769–81.
    OpenUrl
  46. 46.↵
    Bycroft C, Freeman C, Petkova D, et al. The UK Biobank resource with deep phenotyping and genomic data. Nature 2018;562(7726):203–9.
    OpenUrlCrossRefPubMed
  47. 47.↵
    Frank E, Harrell CD. Jr. 2023-09-11; 2023. (https://cran.r-project.org/web/packages/Hmisc/index.html). (Accessed).
  48. 48.↵
    Burgess S, Davey Smith G, Davies NM, et al. Guidelines for performing Mendelian randomization investigations: update for summer 2023. Wellcome Open Res 2019;4:186.
  49. 49.↵
    Chang CC, Chow CC, Tellier LC, et al. Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience 2015;4:7.
  50. 50.↵
    R.C. T. R: A language and environment for statistical computing. 2013.
  51. 51.↵
    Yavorska OO, Burgess S. MendelianRandomization: an R package for performing Mendelian randomization analyses using summarized data. Int J Epidemiol 2017;46(6):1734–9.
    OpenUrlCrossRefPubMed
  52. 52.↵
    Benjamini Y. HY. Controlling the false discovry rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B Methodol 1995;57(1):289–300.
    OpenUrlCrossRef
  53. 53.↵
    Sollis E, Mosaku A, Abid A, et al. The NHGRI-EBI GWAS Catalog: knowledgebase and deposition resource. Nucleic Acids Res 2023;51(D1):D977–D85.
    OpenUrlCrossRef
  54. 54.↵
    Cole JH, Ritchie SJ, Bastin ME, et al. Brain age predicts mortality. Molecular psychiatry 2018;23(5):1385–92.
    OpenUrlPubMed
  55. 55.↵
    Savransky A, Chiappelli J, Fisseha F, et al. Elevated allostatic load early in the course of schizophrenia. Transl Psychiatry 2018;8(1):246.
    OpenUrl
  56. 56.↵
    Soumya Kumari LK, Sundarrajan R. A review on brain age prediction models. Brain Res 2023;1823:148668.
    OpenUrl
  57. 57.↵
    Mattson MP, Arumugam TV. Hallmarks of Brain Aging: Adaptive and Pathological Modification by Metabolic States. Cell Metab 2018;27(6):1176–99.
    OpenUrlCrossRefPubMed
  58. 58.↵
    Herman JP, McKlveen JM, Ghosal S, et al. Regulation of the Hypothalamic-Pituitary-Adrenocortical Stress Response. Compr Physiol 2016;6(2):603–21.
    OpenUrl
  59. 59.↵
    Echouffo-Tcheugui JB, Conner SC, Himali JJ, et al. Circulating cortisol and cognitive and structural brain measures: The Framingham Heart Study. Neurology 2018;91(21):e1961–e70.
    OpenUrlCrossRef
  60. 60.↵
    McEwen BS, Bowles NP, Gray JD, et al. Mechanisms of stress in the brain. Nat Neurosci 2015;18(10):1353–63.
    OpenUrlCrossRefPubMed
  61. 61.↵
    Chetty S, Friedman AR, Taravosh-Lahn K, et al. Stress and glucocorticoids promote oligodendrogenesis in the adult hippocampus. Mol Psychiatry 2014;19(12):1275–83.
    OpenUrlCrossRefPubMed
  62. 62.↵
    El Waly B, Macchi M, Cayre M, et al. Oligodendrogenesis in the normal and pathological central nervous system. Front Neurosci 2014;8:145.
  63. 63.↵
    Long KLP, Chao LL, Kazama Y, et al. Regional gray matter oligodendrocyte- and myelin-related measures are associated with differential susceptibility to stress-induced behavior in rats and humans. Transl Psychiatry 2021;11(1):631.
    OpenUrl
  64. 64.↵
    Fischer R, Maier O. Interrelation of oxidative stress and inflammation in neurodegenerative disease: role of TNF. Oxid Med Cell Longev 2015;2015:610813.
    OpenUrlCrossRefPubMed
  65. 65.↵
    Fry A, Littlejohns TJ, Sudlow C, et al. Comparison of Sociodemographic and Health-Related Characteristics of UK Biobank Participants With Those of the General Population. Am J Epidemiol 2017;186(9):1026–34.
    OpenUrlCrossRefPubMed
  66. 66.↵
    Casey BJ, Cannonier T, Conley MI, et al. The Adolescent Brain Cognitive Development (ABCD) study: Imaging acquisition across 21 sites. Dev Cogn Neurosci 2018;32:43–54.
    OpenUrlCrossRefPubMed
  67. 67.↵
    Ferrucci L. The Baltimore Longitudinal Study of Aging (BLSA): a 50-year-long journey and plans for the future. J Gerontol A Biol Sci Med Sci 2008;63(12):1416–9.
    OpenUrlCrossRefPubMed
Back to top
PreviousNext
Posted May 09, 2024.
Download PDF

Supplementary Material

Data/Code
Email

Thank you for your interest in spreading the word about medRxiv.

NOTE: Your email address is requested solely to identify you as the sender of this article.

Enter multiple addresses on separate lines or separate them with commas.
Association between allostatic load and accelerated white matter brain aging: findings from the UK Biobank
(Your Name) has forwarded a page to you from medRxiv
(Your Name) thought you would like to see this page from the medRxiv website.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Association between allostatic load and accelerated white matter brain aging: findings from the UK Biobank
Li Feng, Zhenyao Ye, Zewen Du, Yezhi Pan, Travis Canida, Hongjie Ke, Song Liu, Shuo Chen, L. Elliot Hong, Peter Kochunov, Jie Chen, David K.Y. Lei, Edmond Shenassa, Tianzhou Ma
medRxiv 2024.01.26.24301793; doi: https://doi.org/10.1101/2024.01.26.24301793
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools
Association between allostatic load and accelerated white matter brain aging: findings from the UK Biobank
Li Feng, Zhenyao Ye, Zewen Du, Yezhi Pan, Travis Canida, Hongjie Ke, Song Liu, Shuo Chen, L. Elliot Hong, Peter Kochunov, Jie Chen, David K.Y. Lei, Edmond Shenassa, Tianzhou Ma
medRxiv 2024.01.26.24301793; doi: https://doi.org/10.1101/2024.01.26.24301793

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Subject Area

  • Epidemiology
Subject Areas
All Articles
  • Addiction Medicine (349)
  • Allergy and Immunology (668)
  • Allergy and Immunology (668)
  • Anesthesia (181)
  • Cardiovascular Medicine (2648)
  • Dentistry and Oral Medicine (316)
  • Dermatology (223)
  • Emergency Medicine (399)
  • Endocrinology (including Diabetes Mellitus and Metabolic Disease) (942)
  • Epidemiology (12228)
  • Forensic Medicine (10)
  • Gastroenterology (759)
  • Genetic and Genomic Medicine (4103)
  • Geriatric Medicine (387)
  • Health Economics (680)
  • Health Informatics (2657)
  • Health Policy (1005)
  • Health Systems and Quality Improvement (985)
  • Hematology (363)
  • HIV/AIDS (851)
  • Infectious Diseases (except HIV/AIDS) (13695)
  • Intensive Care and Critical Care Medicine (797)
  • Medical Education (399)
  • Medical Ethics (109)
  • Nephrology (436)
  • Neurology (3882)
  • Nursing (209)
  • Nutrition (577)
  • Obstetrics and Gynecology (739)
  • Occupational and Environmental Health (695)
  • Oncology (2030)
  • Ophthalmology (585)
  • Orthopedics (240)
  • Otolaryngology (306)
  • Pain Medicine (250)
  • Palliative Medicine (75)
  • Pathology (473)
  • Pediatrics (1115)
  • Pharmacology and Therapeutics (466)
  • Primary Care Research (452)
  • Psychiatry and Clinical Psychology (3432)
  • Public and Global Health (6527)
  • Radiology and Imaging (1403)
  • Rehabilitation Medicine and Physical Therapy (814)
  • Respiratory Medicine (871)
  • Rheumatology (409)
  • Sexual and Reproductive Health (410)
  • Sports Medicine (342)
  • Surgery (448)
  • Toxicology (53)
  • Transplantation (185)
  • Urology (165)