Research Paper Volume 16, Issue 11 pp 9350—9368

Associations of childhood, adolescence, and midlife cognitive function with DNA methylation age acceleration in midlife

Junyu Chen1, , Leah Moubadder1, , Elizabeth S. Clausing1,2, , Katrina L. Kezios3, , Karen N. Conneely4, , Anke Hüls1,5, , Andrea Baccarelli6, , Pam Factor-Litvak3, , Piera Cirrillo7, , Rachel C. Shelton8, , Bruce G. Link9, , Shakira F. Suglia1, ,

  • 1 Department of Epidemiology, Rollins School of Public Health, Emory University, Atlanta, GA 30322, USA
  • 2 School of Global Integrative Studies, University of Nebraska, Lincoln, NE 68508, USA
  • 3 Department of Epidemiology, Mailman School of Public Health, Columbia University, New York, NY 10032, USA
  • 4 Department of Human Genetics, School of Medicine, Emory University, Atlanta, GA 30322, USA
  • 5 Gangarosa Department of Environmental Health, Rollins School of Public Health, Emory University, Atlanta, GA 30322, USA
  • 6 Department of Environmental Health Sciences, Mailman School of Public Health, Columbia University, New York, NY 10032, USA
  • 7 Child Health and Development Studies, Public Health Institute, Washington, DC 20024, USA
  • 8 Department of Sociomedical Sciences, Mailman School of Public Health, Columbia University, New York, NY 10032, USA
  • 9 Department of Sociology, University of California Riverside, Riverside, CA 92507, USA

Received: October 23, 2023       Accepted: April 25, 2024       Published: June 13, 2024      

https://doi.org/10.18632/aging.205943
How to Cite

Copyright: © 2024 Chen et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

Prior studies showed increased age acceleration (AgeAccel) is associated with worse cognitive function among old adults. We examine the associations of childhood, adolescence and midlife cognition with AgeAccel based on DNA methylation (DNAm) in midlife.

Data are from 359 participants who had cognition measured in childhood and adolescence in the Child Health and Development study, and had cognition, blood based DNAm measured during midlife in the Disparities study. Childhood cognition was measured by Raven’s Progressive Matrices and Peabody Picture Vocabulary Test (PPVT). Adolescent cognition was measured only by PPVT. Midlife cognition included Wechsler Test of Adult Reading (WTAR), Verbal Fluency (VF), Digit Symbol (DS). AgeAccel measures including Horvath, Hannum, PhenoAge, GrimAge and DunedinPACE were calculated from DNAm. Linear regressions adjusted for potential confounders were utilized to examine the association between each cognitive measure in relation to each AgeAccel.

There are no significant associations between childhood cognition and midlife AgeAccel. A 1-unit increase in adolescent PPVT, which measures crystalized intelligence, is associated with 0.048-year decrease of aging measured by GrimAge and this association is attenuated after adjustment for adult socioeconomic status. Midlife crystalized intelligence measure WTAR is negatively associated with PhenoAge and DunedinPACE, and midlife fluid intelligence measure (DS) is negatively associated with GrimAge, PhenoAge and DunedinPACE. AgeAccel is not associated with VF in midlife.

In conclusion, our study showed the potential role of cognitive functions at younger ages in the process of biological aging. We also showed a potential relationship of both crystalized and fluid intelligence with aging acceleration.

Introduction

Many neurodegenerative pathologies such as amyloid plaques and neurofibrillary tangles develop long before the clinical diagnosis of cognitive disorders [1]. This has led to many investigations of early-life cognitive outcomes, which have been demonstrated to be instructive in studying late-life cognitive impairment and dementia [2]. Recent research showed that both level of and change in language-based cognition during adolescence were associated with midlife cognition [3], suggesting that prevention of cognitive disorders may be able to start as early as adolescence.

DNA methylation (DNAm), which can regulate gene expression without altering the DNA segments by adding a methyl group at cytosine residues, has been found to be associated with cognitive abilities [4]. DNAm may contribute to learning and memory through a variety of mechanisms including synaptic plasticity and neurogenesis, stem cell function, immunosenescence, circadian rhythms, and the effects of cumulative stress through glucocorticoid signaling [5]. In recent years, epigenetic clocks were created based on DNAm to estimate the construct of ‘biological’ age [610]. DNAm age accelerations, defined as the difference between epigenetic clocks and chronological age or the residuals in epigenetic clocks independent of chronological age, were proposed to measure the discrepancy between chronological and biological age. These DNAm age accelerations have been repeatedly shown to be associated with early-life exposures and diverse age-related disorders and premature mortality [1115]. Investigating the relationship between cognitive functions and DNAm age acceleration could be useful for identification of mechanisms underlying the development of cognitive disorder independent of known neuro-pathologies.

To date, only a handful of studies have examined the relationship between cognitive function and DNAm age acceleration. While previous evidence from a small number of cross-sectional studies was inconsistent [5, 1619], three recent longitudinal studies on cognition at older age (>65 years old) have found that a faster age acceleration was associated with cognitive decline [5, 20, 21]. Two other studies showed that the negative associations between DNAm age acceleration and cognitive skills may exist in late adolescence [22, 23]. However, the relationship between cognitive function and DNAm age acceleration from adolescence to midlife hasn’t been explored. In the current study, we investigated whether cognitive function in childhood and adolescence is related to DNAm age acceleration assessed 35 years after the cognitive function assessment. Many studies have demonstrated the role of epigenetic modifications in cognitive aging at old age [5, 20, 21, 24]. Others have suggested DNAm age acceleration is a mediator for sex differences in verbal memory and processing speed [25, 26]. Thus, we hypothesized that the DNAm age acceleration influenced by cognitive skills at early age could in turn affect one’s cognition in later life and examined the cross-sectional associations of DNAm age acceleration and cognitive measures in midlife. By using this unique study design, we hope to make the first step to understand the relationship between cognition at younger age and aging, which might potentially help refine the strategy to prevent cognitive aging in the future.

Results

We provide an overview of the study population and workflow in Figure 1. Briefly, data from 359 participants who had cognition measured in childhood and adolescence in the Child Health and Development study (CHDS), and had cognition, blood based DNAm measured during midlife in the Disparities study (DISPAR) were used. Childhood (9–11 years) fluid intelligence was measured by Raven Colored Progressive Matrices (RCPM-9) and crystallized intelligence was measured by Peabody Picture Vocabulary Test (PPVT-9). Adolescent (15–17) cognition was measured only by Peabody Picture Vocabulary Test (PPVT-15). Midlife crystallized intelligence was measured by Wechsler Test of Adult Reading (WTAR) and fluid intelligence was measured by Verbal Fluency (VF), Digit Symbol (DS). Midlife DNAm age acceleration measures including Horvath, Hannum, PhenoAge, GrimAge and DunedinPACE were calculated from DNAm. Linear regressions adjusted for potential confounders were utilized to examine the associations between each cognitive measure in relation to each DNAm age acceleration. To align with our hypotheses mentioned above, childhood/adolescent cognitive measures were used as independent variables in the childhood/adolescent analyses while midlife cognitive measures were used as dependent variables in the midlife analyses.

Overview of study population and workflow.

Figure 1. Overview of study population and workflow.

Characteristics of analytic sample

Sociodemographic and clinical characteristics of the participants are summarized in Table 1. Approximately half of the study population are female (46.8%). Over half are non-Hispanic White (56.5%) and 61 (34.3%) are non-Hispanic Black based on self-identified race of the participants in the DISPAR study. Participants who did not self-identify as non-Hispanic White or Black are included as “Other race/ethnicity” (n = 33).

Table 1. Demographic and clinical characteristics of the child health and development studies (CHDS) and the racial disparities (DISPAR) Aging study participants in childhood, adolescence, and midlife (n = 359).

ChildhoodAdolescenceMidlife
Female, n (%)168 (46.8%)
Race/Ethnicity, n (%)
 Non-Hispanic White203 (56.5%)
 Non-Hispanic Black123 (34.3%)
 Other Race/Ethnicity33 (9.2%)
Age, mean (SD)10.0 (0.9)16.6 (0.6)49.3 (1.3)
RCPM, mean (SD)32.3 (10.2)
PPVT, mean (SD)81.3 (11.2)114 (16.5)
WTAR, mean (SD)33.8 (10.2)
VF, mean (SD)23.4 (6.1)
DS, mean (SD)58.8 (18.1)
Socioeconomic status1, n (%)
 Low119 (33.1%)119 (33.0%)
 Medium178 (49.6%)155 (43.2%)
 High62 (17.3%)85 (23.7%)
BMI, mean (SD)17.7 (2.7)22.0 (4.10)29.6 (6.4)
Lived with someone that smoked, n (%)215 (59.9%)215 (59.9%)
Cigarette use2, n (%)
 Never307 (85.5%)192 (53.5%)
 Past46 (12.8%)90 (25.1%)
 Current6 (1.7%)77 (21.4%)
Alcohol use2, n (%)
 Never86 (24.0%)83 (23.1%)
 Past77 (21.4%)
 Current273 (76.0%)199 (55.4%)
Hypertension, n (%)98 (27.3%)
Horvath3, mean (SD)0.0 (4.34)
Hannum3, mean (SD)0.0 (4.11)
GrimAge3, mean (SD)0.0 (4.39)
PhenoAge3, mean (SD)0.0 (5.94)
DunedinPACE, mean (SD)1.0 (0.136)
Abbreviation: RCPM: Raven Colored Progressive Matrices; PPVT: Peabody Picture Vocabulary Test; WTAR: Wechsler Test of Adult Reading; VF: Verbal Fluency; DS: Digit Symbol; SES: Socioeconomic Status; BMI: Body Mass Index. 1One SES measure was created based on maternal education at birth, paternal occupation, and family income at birth, ages 9–11 and 15–17 to represent SES at both childhood and adolescence. Midlife SES index was developed by combining information on college education, self-reported family income and occupational status. 2Information on cigarette and alcohol use in adolescence were measured and defined in difference matrices from those in midlife. 3Horvath, Hannum, GrimAge and PhenoAge are the residuals in the corresponding epigenetic clocks after regressing out chronological age.

All 5 DNAm age acceleration measures follow a normal distribution (Supplementary Figure 1) but DunedinPACE has a much smaller standard deviation (SD) (Table 1). We show that DunedinPACE has a moderate correlation with GrimAge (R = 0.6) but neither measure correlates with Horvath or Hannum (R ≤ 0.15) (Supplementary Figure 2C). PhenoAge has a weak positive correlation with all other age acceleration measures (R ≈ 0.4).

All cognitive scores were comparable to average scores in the general population. All cognitive measures in childhood and adolescence are positively correlated with each other before and after multiple imputation (Supplementary Figure 2A, 2B). Distributions of midlife cognitive function are similar before and after multiple imputation (Supplementary Figure 3). VF and DS follow a normal distribution and WTAR is moderately left-skewed. WTAR shows moderate correlations with all childhood and adolescent cognitive measures, especially PPVT-15 (R = 0.77) as they all measure crystallized intelligence that’s developed at an early age (Supplementary Figure 2B). The three midlife cognitive measures have low correlations with each other (R ranges 0.21–0.4), given they assess different components of adult cognition [2730].

Childhood/adolescent cognitive function and midlife DNAm age acceleration

There are no statistically significant associations between childhood cognitive function and midlife DNAm age acceleration measures (Figure 2). Adolescent cognitive measure (PPVT-15) is not significantly associated with Horvath, Hannam and PhenoAge age acceleration but GrimAge and DunedinPACE. A 1-unit increase in PPVT-15 is associated with 0.048-year (1.1% SD) decrease of aging measured by GrimAge and 0.001 decrease (0.71% SD) in the pace of aging measured by DunedinPACE (Figure 2). However, the latter is not statistically significant after adjusting for multiple testing. Associations from complete case analyses are similar to the ones from analyses after multiple imputation (Supplementary Figure 4). After adjusting for adulthood SES, the absolute value of Beta for the associations of PPVT-15 with GrimAge and DunedinPACE has shrunk to half (Table 2).

Forest plot showing associations of childhood and adolescent cognitive function with midlife DNA methylation age accelerations. Betas multiplied by 10 were shown for associations with DunedinPACE for better visualization. Column “Change” is showing that 1-unit change in childhood or adolescent cognitive function is associated with X.XX% standard deviation (SD) change in DNA methylation age acceleration. SD for each DNA methylation age can be found in Table 1.

Figure 2. Forest plot showing associations of childhood and adolescent cognitive function with midlife DNA methylation age accelerations. Betas multiplied by 10 were shown for associations with DunedinPACE for better visualization. Column “Change” is showing that 1-unit change in childhood or adolescent cognitive function is associated with X.XX% standard deviation (SD) change in DNA methylation age acceleration. SD for each DNA methylation age can be found in Table 1.

Table 2. Comparison of associations of adolescent cognitive function (PPVT-15) with midlife DNA methylation age accelerations before and after adjustment for adult SES.

Before adjustment for adult SESAfter adjustment for adult SES
Age accelerationBetaSEAge accelerationBetaSEAge accelerationBetaSE
Horvath−0.0180.0160.283−0.406%−0.0140.0170.413−0.323%
Hannum0.0120.0130.3730.285%0.0110.0140.4210.268%
GrimAge−0.0480.0151.693E-03−1.101%−0.0250.0150.105−0.566%
PhenoAge−0.0090.0210.665−0.154%−0.0050.0220.804−0.092%
DunedinPACE−9.653E-034.453E-030.031−0.710%−4.521E-034.549E-030.321−0.332%
Column “% SD change” is showing that 1-unit change in PPVT-15 is associated with X.XXX% standard deviation (SD) change in DNA methylation age acceleration. SD for each DNA methylation age can be found in Table 1.

Midlife cognitive function and midlife DNAm age acceleration

WTAR shows negative associations with all 5 age acceleration measures but only two are statistically significant with FDR-q value less than 0.05 (Figure 3). Specifically, 1-SD increase in PhenoAge is associated with 0.89-unit decrease in WTAR score. A 1-SD increase in DunedinPACE is associated with 1.12-unit decrease in DunedinPACE. Associations between WTAR on PhenoAge and DunedinPACE are no longer in complete case analysis than analysis after multiple imputation (Supplementary Figure 5).

Forest plot showing associations of midlife DNA methylation age acceleration with midlife cognitive function. Betas divided by 100 were shown for associations with DunedinPACE for better visualization. Column “Change” is showing that 1-standard deviation (SD) change in DNA methylation age acceleration is associated with X.XX-unit change in midlife cognitive function. SD for each DNA methylation age can be found in Table 1.

Figure 3. Forest plot showing associations of midlife DNA methylation age acceleration with midlife cognitive function. Betas divided by 100 were shown for associations with DunedinPACE for better visualization. Column “Change” is showing that 1-standard deviation (SD) change in DNA methylation age acceleration is associated with X.XX-unit change in midlife cognitive function. SD for each DNA methylation age can be found in Table 1.

VF test score decreases by 0.94 per SD increase of pace of aging measured by DunedinPACE and remains significant after adjusting for multiple testing (Figure 3). Although the estimated beta coefficient is similar to the one in the complete case analysis, it is no longer significant in the complete case analysis (Supplementary Figure 5). VF is not significantly associated with other age acceleration measures.

DS is not significantly associated with Horvath and Hannum age acceleration (Figure 3). A 1-SD year increase in aging measured by GrimAge, PhenoAge and DunedinPACE is associated with 3.57, 2.91 and 4.1 decrease in DS respectively. These three associations remain significant even after multiple testing adjustment, however, the precision of the estimates is low (wide 95% confidence interval). In complete case analysis, associations are larger than analysis after multiple imputation, with 1-SD year increase in GrimAge, PhenoAge or DunedinPACE associated with more than 4-unit decrease in DS.

Discussion

We investigated the associations between childhood, adolescent and midlife cognitive function with DNAm age acceleration in midlife. Our findings show that those with poorer crystalized intelligence (PPVT-15) in adolescence are more likely to have accelerated GrimAge in their midlife. We also show that accelerated PhenoAge and DunedinPACE in midlife is associated with both reduced crystalized intelligence (WTAR) and fluid intelligence (DS), while accelerated GrimAge in midlife is only associated with crystalized intelligence (WTAR). In the complete case analyses, the magnitude of the associations is similar, however, only the associations of DNAm age acceleration with fluid intelligence remain statistically significant, which could be due to the reduced sample size.

In previous studies among adults older than 65 years old [5, 20, 21], it’s hypothesized that DNAm age acceleration influences cognitive functions through numerous mechanisms such as synaptic plasticity and neurogenesis and immune dysfunction [31, 32]. Our study showed that crystallized intelligence measured by PPVT-15 is associated with age acceleration measured 35 years later, demonstrating the possibility of a reverse causation or even a bidirectional relationship, that is, the possibility that cognitive function at a younger age could be a cause rather than just the consequences of biological aging. Further research with more time points of cognitive function and DNAm age acceleration is needed to establish causality.

The effect of early crystallized intelligence on later biological aging could be through its influence in daily stress processes, behavior, lifestyle or health [33], which might be starting to accumulate during adolescence. For example, a higher crystallized cognition may increase the likelihood of exposure to daily stressors [33], which could influence DNAm age acceleration [34]. The attenuated association between adolescent crystallized intelligence and GrimAge acceleration after adjusting for adult SES suggests that adult SES could be a potential mediator (Supplementary Figure 6A) and improving SES could be helpful for slowing down some aspects of biological aging associated with cognition at younger ages. However, the attenuated association could also be due to unmeasured confounders in our study that induce confounding or collider bias through adult SES (Supplementary Figure 6B, 6C).

On the other hand, accelerated aging is associated with both crystallized and fluid intelligence in midlife. Crystallized intelligence measured by WTAR in midlife is statistically associated with PhenoAge and DunedinPACE but not GrimAge. This observed heterogeneity may be explained by the inherent difference in the methodology used to construct the epigenetic clocks. Early cognition may be more related to biological aging in response to physiological risk factors and stress factors [9], while later cognition may be more related to biological aging in terms of phenotypic changes and organ failures [8, 35]. GrimAge, PhenoAge and DunedinPACE are significantly associated with fluid intelligence measured by DS, and all associations are robust to sensitivity analyses. This suggests that it is possible that fluid intelligence, which refers to the capacity to process complex information involved in reasoning and problem-solving tasks, may be more affected by accelerated aging than language processing skills. Midlife DNAm age acceleration may influence midlife crystallized and fluid intelligence through inflammatory or oxidative stress pathways [36].

Previous studies in older population found that faster decline of memory and executive function are associated with Horvath and Hannum DNAm age acceleration [5, 20], but neither of these measures are predictive of any cognitive function measure in our study. This could be due to a smaller sample size or usage of different cognition measures in our study. Conversely, our results suggest that a more progressive DNAm age acceleration measure that indicates mortality or morbidity such as GrimAge, PhenoAge and DunedinPACE may be more sensitive to changes in cognitive function compared to DNAm age acceleration that only considers chronological age.

To the best of our knowledge, this is the first study to show the association of cognition at younger age with midlife age acceleration, and associations between midlife age acceleration measures and cognitive function that are independent of childhood and adolescent cognition. This aligns with our hypothesis that cognitive function at early age may affect aging in midlife and in turn affect cognitive function in midlife. However, we only have one time point of DNAm in midlife and no information on cognition, DNAm or other covariates during the 35-year gap. This limits the inference on the causal relationship between cognitive function and DNAm age acceleration. Selection of participants into the DISPAR study at midlife was limited to those who continued to live in California at midlife, which may have introduced selection bias [3]. Another limitation that warrants mention is that the fluid intelligence was only administered in childhood, so we could not assess its association with DNAm age acceleration in midlife. Additionally, the measures of crystalized and fluid intelligence were different between childhood/adolescence and midlife, so we were not able to efficiently describe cognitive function trajectories over time. The differences in measures may also account for the observed difference in their associations with DNAm age acceleration.

In conclusion, our study brings attention to the potential influence of adolescent crystalized intelligence on age-related DNAm at older age. We showed a potential relationship of both crystalized and fluid intelligence with aging acceleration in midlife. Moreover, we demonstrated that GrimAge, PhenoAge, DunedinPACE age acceleration measures could prove valuable for studying the relationship between aging and cognitive function. Future studies with longitudinal cognitive functions and DNAm age accelerations at more time points across the lifespan are needed to decipher their relationships.

Materials and Methods

Study population

Our study utilized data from adult offspring of mothers that were recruited in the Child Health and Development Studies (CHDS), which was a prospective study of pregnant women seeking prenatal care from the Kaiser Foundation Health Plan in Oakland, California from 1959 to 1966 [37]. This original cohort was comprised of 19,044 live births from participating mothers with diverse race/ethnic background. Extensive sociodemographic, behavioral and clinical information on the offspring and their mothers were collected through in-person interviews and medical records [37]. Follow-up assessments of cognitive function in specific subsets of offspring were conducted subsequently at age 5, 9–11, and 15–17 years.

After approximately 35 years, a subset of the adult offspring who had assessments during adolescence were selected for another follow-up examination and interview for the CHDS Disparities (DISPAR) study [38]. Among 985 eligible pool that could be located and recruited by phone, 605 were successfully enrolled and participated in a telephone interview, 510 completed a home visit and 497 completed a self-administered questionnaire. At the home visit, participants completed a short cognitive battery and a blood draw that resulted in 400 viable serology samples after excluding participants who refused or had a technical issue during the draw [39].

The current study, DISPAR Aging, builds from the CHDS DISPAR study by using the blood samples it collected to profile for DNA methylation. From the sample of 400, we further excluded adult offspring who refused DNA analysis or use of their blood samples for future studies, resulting in a final analytic sample of 359.

Cognitive function measures

Childhood cognitive function was measured by two tests, the Raven Colored Progressive Matrices test at the age of 9–11 years (RCPM-9) and the Peabody Picture Vocabulary test at the age of 9–11 years (PPVT-9) [3]. RCPM-9 consisted of series of items each containing a matrix display of 12 graphical patterns where one of them was missing. The children were provided with 6 alternative patterns and asked to choose the best that fits the missing element. It is a test for logical and perceptual reasoning and is typically used to assess fluid intelligence [40]. In PPVT-9 test, participants were shown a series of sheets with 4 pictures and were asked to indicate which best represents a target word said by the examiner. This is a test for vocabulary and language processing, which is typically used to assess crystallized intelligence [41]. Adolescent cognitive function was only measured by the Peabody Picture Vocabulary Test at the age of 15–17 years (PPVT-15) [3]. The RCPM-9 raw score can range from 0–60. The average range for PPVT raw scores is 72.28–92.18 at the age of 11 years and 90.98–119.24 at the age of 17 years.

Three cognitive function measures were assessed in midlife, including the Wechsler Test of Adult Reading (WTAR) [27], Verbal Fluency (VF) test [28], the Wechsler Adult Intelligence Scale-Revised Digit Symbol (DS) test [29]. WTAR intended to measure participants’ intellectual functioning, reading and recognition abilities by asking the participants to pronounce a list of 50 phonetically incorrect words [27]. VF was typically used to evaluate participants’ executive functions, semantic fluency, and speed of processing by asking participants to name as many animals as possible within 1 minute [28]. For DS test, each participant was presented with a sheet with a number-to-symbol key at the top, along with a grid of numbers and blank boxes, and was asked to write down as many correct corresponding symbols underneath each number. This test measures participants’ executive functions and speed of processing [29]. The WTAR and DS test can range from 0–50 and 0–133 respectively. There is no limit for the maximum VF test score, however, a VF test score less than 15 could indicate potential cognitive impairment [42]. Details of the cognitive measures have been previously published [3]. Spearman correlations were calculated between cognitive function measures.

DNAm age acceleration in midlife

Genomic DNA was extracted from whole blood samples using with the PureLink® Genomic DNA Mini Kit by Invitrogen (Carlsbad, CA, USA). DNA methylation was measured using 850K BeadChips following the standard protocol. Intensity files in IDAT format were imported in R using minfi package [43]. Intensity values with a detection P-value ≥ 0.001 were set to missing for each DNAm site and DNAm sites missing in more than 5% of the samples were removed. All samples passed a call rate of 95%. Normal-exponential out-of-band background correction was used to account for technical variation in the intensity values, which were then used to generate a β value for each DNAm site. We calculated the proportions of 6 cell types including CD4+ T cells, CD8+ T cells, natural killer T cells, B cells, monocytes, and granulocytes in blood by applying algorithm developed by Houseman et al to the DNAm data [44].

We considered 5 different measures of DNAm age acceleration, including Horvath [10], Hannum [7], GrimAge [9], PhenoAge [8] and DunedinPACE [35]. For Horvath [10] and Hannum [7] epigenetic clocks, the difference in DNAm between younger and older people represents biological processes of aging while GrimAge and PhenoAge epigenetic clocks are using the difference in DNAm patterns that predict mortality and health, respectively, to represent biological processes of aging [8, 9]. These four epigenetic clocks were generated by uploading the DNAm β values of 30,084 prespecified DNAm sites and information on age and sex to the online DNA Methylation Age Calculator (https://dnamage.genetics.ucla.edu/) with the “normalization” and “advanced analysis in blood” options. DNAm age accelerations were calculated as the residuals in epigenetic clocks after regressing out chronological age, thus a positive DNAm age acceleration reflects more rapid biological aging.

Recently, Belsky et al. applied elastic-net regression to develop a DNA methylation model that predicts the pace of aging, which was quantified from 19 organ-system integrity data of four time points spanning two decades [35]. This resulted in a new epigenetic biomarker named DunedinPACE, with values larger than 1 indicating faster aging and values less than 1 indicating slower aging. The idea behind DunedinPACE is that the DNAm responses to the variation in function decline of multiple organ systems represents the biological processes of aging. DunedinPACE measures someone’s pace of aging in 20 years while the other clock-based age acceleration measures the accumulated years of aging at the time of assessment in a cross-sectional setting. We calculated DunedinPACE using the DunedinPACE R package. Spearman correlations were calculated between DNAm age acceleration measures.

Additional covariates

Childhood and adolescence

Participants’ exact age, school grade level and BMI at the time of cognition assessment were recorded. A composite score for childhood/adolescence socioeconomic status (SES) was created based on maternal education at birth, paternal occupation, and family income at birth, ages 9–11 and 15–17 [38]. At the age of 15–17 years, participants were asked how often they had consumed alcoholic beverage in the past 6 months and those who answered at least one drink were defined as current drinkers. They were also asked if they had ever ‘regularly’ smoked at least 1 cigarette daily and if they currently smoked cigarettes. Alcohol and cigarette use at age 9–11 was not recorded. Additionally, during interviews in midlife, participants were asked to recall if they had lived with a regular cigarette, cigar, or pipe smoker who smoked in their home when they were young. We note that information on BMI, alcohol use and personal cigarette use in childhood and adolescence was only used in multiple imputation but not in analyses examining cognition and age acceleration associations (see multiple imputation and statistical analyses method section for more details).

Midlife

Participants’ exact age and BMI were collected at the time of midlife cognition assessment. Midlife SES index was developed by combining information on college education, self-reported family income and occupational status. Details on how SES indices were developed have been described in detail in a prior publication [39]. Participants were asked about alcohol and cigarette use during the midlife interviews in DISPAR. If they had ever consumed any alcoholic beverages at least once a month for six months over the past 12 months, they were defined as current drinkers. If they had done so but not over the past 12 months, they were defined as past drinkers. Otherwise, they were categorized as “never drink”. Participants were also asked if they had ever smoked at least 1 cigarette per day for one month or more and if they smoked at any time in the past 12 months. If they answered yes to both questions, they were defined as current smokers. If they answered yes to the first question but no to the second question, they were defined as past smokers. Otherwise, they were categorized as “never smokers”. Blood pressure was measured during the home visit and hypertension was defined as having systolic blood pressure over or equal to 120 mm Hg or diastolic blood pressure over or equal to 80 mm Hg, or currently taking hypertension medication.

Multiple imputation

Of the 359 participants who had DNAm data, 26 were missing RCPM-9, 34 were missing PPVT-9, 61 were missing PPVT-15 and 19 were missing WTAR. We found that RCPM-9, PPVT-9, PPVT-15 and WTAR were positively correlated (Supplementary Figure 2A), thus, we used multiple imputation by chained equations to impute missing values on cognitive function and other covariates. Multiple imputation was performed among participants that were successfully enrolled in DISPAR, with or without DNAm data, and had at least one measure of RCPM-9, PPVT-9, PPVT-15 or WTAR (n = 599, Figure 1). We assumed missing at random and performed multiple imputation using Predictive mean matching by MICE package in R [45]. We assumed variables at the same time point can predict each other and variables from different time points cannot predict each other except for the variables that measure related information, such as adolescent and midlife smoking status. The predictor matrix is presented in Supplementary Table 1. Twenty complete datasets were imputed, and estimates obtained by statistical analyses were combined using “pool” function in the MICE package in R [45]. The convergence of cognitive function measures is shown in Supplementary Figure 7.

Statistical analyses

All analyses were conducted in R (4.3.0). Potential confounders were identified through literature review and directed acyclic graphs (DAG) (Supplementary Figure 8A8C).

Childhood and adolescence

Linear regressions were used to assess associations of each childhood or adolescent cognitive function measure (independent variable) with each DNAm age acceleration measure in midlife (dependent variable). All models were adjusted for age and school grade level at assessment, sex, race, household smoking status and cell type proportions. According to the DAG for adolescent cognitive function, BMI, alcohol use and personal cigarette use in adolescence could be potential mediators [4649] (Supplementary Figure 8B). Thus, these covariates were not included in the linear regressions but were used in multiple imputation. In a separate model, we additionally adjusted for adult SES to assess its impact on the associations between adolescent cognition and midlife aging. In order to adjust for multiple testing, Benjamini-Hochberg false discovery rate (FDR) adjustment was applied and FDR-q-values < 0.05 were used to define statistical significance. Sensitivity analyses by complete case analyses (n = 286) were performed to test the robustness of the associations.

Midlife

Another set of linear regressions were used to assess associations of each DNAm age acceleration measure in midlife with each midlife cognitive function measure. Similar to the majority of the current literature, we hypothesize that age acceleration in midlife reflects general deterioration and has an impact on cognitive function in midlife. Thus, in these analyses, age acceleration was treated as the independent variable and cognitive function was treated as the dependent variable. All models were adjusted for age at blood draw, sex, race, SES, BMI, alcohol and cigarette use, RCPM-9, PPVT-9, PPVT-15 and cell type proportions. We included RCPM-9, PPVT-9 and PPVT-15 in the models because we wanted to evaluate the associations in midlife independent of the influence from cognition at younger age. Unlike in childhood and adolescence, BMI, alcohol, and cigarette use in midlife are likely to be confounders (Supplementary Figure 8C), and thus were adjusted in the models. Sensitivity analyses by complete case analyses (n = 211) were performed to test the robustness of the associations.

Author Contributions

JC and SS designed the study; JC performed the analyses; JC and SS wrote the manuscript; LM cleaned the phenotypic data and assisted with designing the study; ES cleaned the methylation data; KC, AH, AB provided feedback on analyses; KLK, PFL, PC, RCS and BGL oversaw data collection and provided feedback on study design; SS supervised the study and provided funding for analysis; All authors provided critical comments on the final version of the manuscript.

Acknowledgments

The authors thank the staff and participants of the CHDS and DISPAR study. The authors also thank Barbara Cohn for assistance on an earlier version of the manuscript.

Conflicts of Interest

The authors declare no conflicts of interest related to this study.

Ethical Statement and Consent

This study was approved by the Institutional Review Board at Columbia University and the Public Health Institute, Oakland California. The Emory IRB approval number for DISPAR is IRB00100708. All CHDS and DISPAR participants provided informed consent.

Funding

This research was funded in part by the National Institute for Child Health and Development grant R01HD058515 and the National Institute of Aging Grant R01AG058704.

References

  • 1. Launer LJ. The epidemiologic study of dementia: a life-long quest? Neurobiol Aging. 2005; 26:335–40. https://doi.org/10.1016/j.neurobiolaging.2004.03.016 [PubMed]
  • 2. Novotný JS, Gonzalez-Rivas JP, Medina-Inojosa JR, Lopez-Jimenez F, Geda YE, Stokin GB. Investigating cognition in midlife. Alzheimers Dement (N Y). 2021; 7:e12234. https://doi.org/10.1002/trc2.12234 [PubMed]
  • 3. Agha G, Kezios K, Baccarelli AA, DuBois Bowman F, Rauh V, Susser ES, Cohn B, Cirillo P, Link BG, Factor-Litvak P, Staudinger UM. Cognition level and change in cognition during adolescence are associated with cognition in midlife. Ann Epidemiol. 2019; 35:48–52.e2. https://doi.org/10.1016/j.annepidem.2019.03.004 [PubMed]
  • 4. Marioni RE, McRae AF, Bressler J, Colicino E, Hannon E, Li S, Prada D, Smith JA, Trevisi L, Tsai PC, Vojinovic D, Simino J, Levy D, et al. Meta-analysis of epigenome-wide association studies of cognitive abilities. Mol Psychiatry. 2018; 23:2133–44. https://doi.org/10.1038/s41380-017-0008-y [PubMed]
  • 5. Vaccarino V, Huang M, Wang Z, Hui Q, Shah AJ, Goldberg J, Smith N, Kaseer B, Murrah N, Levantsevych OM, Shallenberger L, Driggers E, Bremner JD, Sun YV. Epigenetic Age Acceleration and Cognitive Decline: A Twin Study. J Gerontol A Biol Sci Med Sci. 2021; 76:1854–63. https://doi.org/10.1093/gerona/glab047 [PubMed]
  • 6. Horvath S, Raj K. DNA methylation-based biomarkers and the epigenetic clock theory of ageing. Nat Rev Genet. 2018; 19:371–84. https://doi.org/10.1038/s41576-018-0004-3 [PubMed]
  • 7. Hannum G, Guinney J, Zhao L, Zhang L, Hughes G, Sadda S, Klotzle B, Bibikova M, Fan JB, Gao Y, Deconde R, Chen M, Rajapakse I, et al. Genome-wide methylation profiles reveal quantitative views of human aging rates. Mol Cell. 2013; 49:359–67. https://doi.org/10.1016/j.molcel.2012.10.016 [PubMed]
  • 8. Levine ME, Lu AT, Quach A, Chen BH, Assimes TL, Bandinelli S, Hou L, Baccarelli AA, Stewart JD, Li Y, Whitsel EA, Wilson JG, Reiner AP, et al. An epigenetic biomarker of aging for lifespan and healthspan. Aging (Albany NY). 2018; 10:573–91. https://doi.org/10.18632/aging.101414 [PubMed]
  • 9. Lu AT, Quach A, Wilson JG, Reiner AP, Aviv A, Raj K, Hou L, Baccarelli AA, Li Y, Stewart JD, Whitsel EA, Assimes TL, Ferrucci L, Horvath S. DNA methylation GrimAge strongly predicts lifespan and healthspan. Aging (Albany NY). 2019; 11:303–27. https://doi.org/10.18632/aging.101684 [PubMed]
  • 10. Horvath S. DNA methylation age of human tissues and cell types. Genome Biol. 2013; 14:R115. https://doi.org/10.1186/gb-2013-14-10-r115 [PubMed]
  • 11. Marioni RE, Shah S, McRae AF, Chen BH, Colicino E, Harris SE, Gibson J, Henders AK, Redmond P, Cox SR, Pattie A, Corley J, Murphy L, et al. DNA methylation age of blood predicts all-cause mortality in later life. Genome Biol. 2015; 16:25. https://doi.org/10.1186/s13059-015-0584-6 [PubMed]
  • 12. Perna L, Zhang Y, Mons U, Holleczek B, Saum KU, Brenner H. Epigenetic age acceleration predicts cancer, cardiovascular, and all-cause mortality in a German case cohort. Clin Epigenetics. 2016; 8:64. https://doi.org/10.1186/s13148-016-0228-z [PubMed]
  • 13. Salameh Y, Bejaoui Y, El Hajj N. DNA Methylation Biomarkers in Aging and Age-Related Diseases. Front Genet. 2020; 11:171. https://doi.org/10.3389/fgene.2020.00171 [PubMed]
  • 14. Chen BH, Marioni RE, Colicino E, Peters MJ, Ward-Caviness CK, Tsai PC, Roetker NS, Just AC, Demerath EW, Guan W, Bressler J, Fornage M, Studenski S, et al. DNA methylation-based measures of biological age: meta-analysis predicting time to death. Aging (Albany NY). 2016; 8:1844–65. https://doi.org/10.18632/aging.101020 [PubMed]
  • 15. Dhingra R, Nwanaji-Enwerem JC, Samet M, Ward-Caviness CK. DNA Methylation Age-Environmental Influences, Health Impacts, and Its Role in Environmental Epidemiology. Curr Environ Health Rep. 2018; 5:317–27. https://doi.org/10.1007/s40572-018-0203-2 [PubMed]
  • 16. Maddock J, Castillo-Fernandez J, Wong A, Cooper R, Richards M, Ong KK, Ploubidis GB, Goodman A, Kuh D, Bell JT, Hardy R. DNA Methylation Age and Physical and Cognitive Aging. J Gerontol A Biol Sci Med Sci. 2020; 75:504–11. https://doi.org/10.1093/gerona/glz246 [PubMed]
  • 17. Starnawska A, Tan Q, Lenart A, McGue M, Mors O, Børglum AD, Christensen K, Nyegaard M, Christiansen L. Blood DNA methylation age is not associated with cognitive functioning in middle-aged monozygotic twins. Neurobiol Aging. 2017; 50:60–3. https://doi.org/10.1016/j.neurobiolaging.2016.10.025 [PubMed]
  • 18. Marioni RE, Shah S, McRae AF, Ritchie SJ, Muniz-Terrera G, Harris SE, Gibson J, Redmond P, Cox SR, Pattie A, Corley J, Taylor A, Murphy L, et al. The epigenetic clock is correlated with physical and cognitive fitness in the Lothian Birth Cohort 1936. Int J Epidemiol. 2015; 44:1388–96. https://doi.org/10.1093/ije/dyu277 [PubMed]
  • 19. Sugden K, Caspi A, Elliott ML, Bourassa KJ, Chamarti K, Corcoran DL, Hariri AR, Houts RM, Kothari M, Kritchevsky S, Kuchel GA, Mill JS, Williams BS, et al, and Alzheimer's Disease Neuroimaging Initiative*. Association of Pace of Aging Measured by Blood-Based DNA Methylation With Age-Related Cognitive Impairment and Dementia. Neurology. 2022; 99:e1402–13. https://doi.org/10.1212/WNL.0000000000200898 [PubMed]
  • 20. Beydoun MA, Shaked D, Tajuddin SM, Weiss J, Evans MK, Zonderman AB. Accelerated epigenetic age and cognitive decline among urban-dwelling adults. Neurology. 2020; 94:e613–25. https://doi.org/10.1212/WNL.0000000000008756 [PubMed]
  • 21. Reed RG, Carroll JE, Marsland AL, Manuck SB. DNA methylation-based measures of biological aging and cognitive decline over 16-years: preliminary longitudinal findings in midlife. Aging (Albany NY). 2022; 14:9423–44. https://doi.org/10.18632/aging.204376 [PubMed]
  • 22. Niccodemi G, Menta G, Turner J, D'Ambrosio C. Pace of aging, family environment and cognitive skills in children and adolescents. SSM Popul Health. 2022; 20:101280. https://doi.org/10.1016/j.ssmph.2022.101280 [PubMed]
  • 23. Raffington L, Tanksley PT, Sabhlok A, Vinnik L, Mallard T, King LS, Goosby B, Harden KP, Tucker-Drob EM. Socially Stratified Epigenetic Profiles Are Associated With Cognitive Functioning in Children and Adolescents. Psychol Sci. 2023; 34:170–85. https://doi.org/10.1177/09567976221122760 [PubMed]
  • 24. Day JJ, Sweatt JD. DNA methylation and memory formation. Nat Neurosci. 2010; 13:1319–23. https://doi.org/10.1038/nn.2666 [PubMed]
  • 25. O'Shea DM, Alaimo H, Davis JD, Galvin JE, Tremont G. A comparison of cognitive performances based on differing rates of DNA methylation GrimAge acceleration among older men and women. Neurobiol Aging. 2023; 123:83–91. https://doi.org/10.1016/j.neurobiolaging.2022.12.011 [PubMed]
  • 26. O'Shea DM, Maynard T, Tremont G. DNA Methylation "GrimAge" Acceleration Mediates Sex/Gender Differences in Verbal Memory and Processing Speed: Findings From the Health and Retirement Study. J Gerontol A Biol Sci Med Sci. 2022; 77:2402–12. https://doi.org/10.1093/gerona/glac133 [PubMed]
  • 27. Venegas J, Clark E. Wechsler Test of Adult Reading. In: Kreutzer JS, DeLuca J, Caplan B, eds. Encyclopedia of Clinical Neuropsychology. Springer: New York; 2011; 2693–4. https://doi.org/10.1007/978-0-387-79948-3_1500
  • 28. Brickman AM, Paul RH, Cohen RA, Williams LM, MacGregor KL, Jefferson AL, Tate DF, Gunstad J, Gordon E. Category and letter verbal fluency across the adult lifespan: relationship to EEG theta power. Arch Clin Neuropsychol. 2005; 20:561–73. https://doi.org/10.1016/j.acn.2004.12.006 [PubMed]
  • 29. Wechsler D. Wechsler Adult Intelligence Scale-Revised. New York, NY: Psychological.
  • 30. Cadar D, Pikhart H, Mishra G, Stephen A, Kuh D, Richards M. The role of lifestyle behaviors on 20-year cognitive decline. J Aging Res. 2012; 2012:304014. https://doi.org/10.1155/2012/304014 [PubMed]
  • 31. Grolleau-Julius A, Ray D, Yung RL. The role of epigenetics in aging and autoimmunity. Clin Rev Allergy Immunol. 2010; 39:42–50. https://doi.org/10.1007/s12016-009-8169-3 [PubMed]
  • 32. Jiang Y, Langley B, Lubin FD, Renthal W, Wood MA, Yasui DH, Kumar A, Nestler EJ, Akbarian S, Beckel-Mitchener AC. Epigenetics in the nervous system. J Neurosci. 2008; 28:11753–9. https://doi.org/10.1523/JNEUROSCI.3797-08.2008 [PubMed]
  • 33. Stawski RS, Mogle JA, Sliwinski MJ. Associations among fluid and crystallized cognition and daily stress processes in older adults. Psychol Aging. 2013; 28:57–63. https://doi.org/10.1037/a0029813 [PubMed]
  • 34. Suglia SF, Clausing ES, Shelton RC, Conneely K, Prada-Ortega D, DeVivo I, Factor-Litvak P, Cirillo P, Baccarelli AA, Cohn B, Link BG. Cumulative Stress Across the Life Course and Biological Aging in Adulthood. Psychosom Med. 2024; 86:137–45. https://doi.org/10.1097/PSY.0000000000001284 [PubMed]
  • 35. Belsky DW, Caspi A, Corcoran DL, Sugden K, Poulton R, Arseneault L, Baccarelli A, Chamarti K, Gao X, Hannon E, Harrington HL, Houts R, Kothari M, et al. DunedinPACE, a DNA methylation biomarker of the pace of aging. Elife. 2022; 11:e73420. https://doi.org/10.7554/eLife.73420 [PubMed]
  • 36. Cencioni C, Spallotta F, Martelli F, Valente S, Mai A, Zeiher AM, Gaetano C. Oxidative stress and epigenetic regulation in ageing and age-related diseases. Int J Mol Sci. 2013; 14:17643–63. https://doi.org/10.3390/ijms140917643 [PubMed]
  • 37. van den Berg BJ, Christianson RE, Oechsli FW. The California Child Health and Development Studies of the School of Public Health, University of California at Berkeley. Paediatr Perinat Epidemiol. 1988; 2:265–82. https://doi.org/10.1111/j.1365-3016.1988.tb00218.x [PubMed]
  • 38. Link BG, Susser ES, Factor-Litvak P, March D, Kezios KL, Lovasi GS, Rundle AG, Suglia SF, Fader KM, Andrews HF, Johnson E, Cirillo PM, Cohn BA. Disparities in self-rated health across generations and through the life course. Soc Sci Med. 2017; 174:17–25. https://doi.org/10.1016/j.socscimed.2016.11.035 [PubMed]
  • 39. Kezios KL, Suglia SF, Doyle DM, Susser E, Bradwin G, Cirillo P, Cohn B, Link B, Factor-Litvak P. Comparing different operationalizations of allostatic load measured in mid-life and their patterning by race and cumulative life course socioeconomic status. Psychoneuroendocrinology. 2022; 139:105689. https://doi.org/10.1016/j.psyneuen.2022.105689 [PubMed]
  • 40. Raven JC. Guide to using the Coloured Progressive Matrices. H. K. Lewis & Co.; 1958; 40.
  • 41. Dunn L. Expanded manual for the peabody picture. Circle Pines, Minn: American Guidance Service, Inc.; 1965.
  • 42. Canning SJ, Leach L, Stuss D, Ngo L, Black SE. Diagnostic utility of abbreviated fluency measures in Alzheimer disease and vascular dementia. Neurology. 2004; 62:556–62. https://doi.org/10.1212/wnl.62.4.556 [PubMed]
  • 43. Aryee MJ, Jaffe AE, Corrada-Bravo H, Ladd-Acosta C, Feinberg AP, Hansen KD, Irizarry RA. Minfi: a flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinformatics. 2014; 30:1363–9. https://doi.org/10.1093/bioinformatics/btu049 [PubMed]
  • 44. Houseman EA, Kelsey KT, Wiencke JK, Marsit CJ. Cell-composition effects in the analysis of DNA methylation array data: a mathematical perspective. BMC Bioinformatics. 2015; 16:95. https://doi.org/10.1186/s12859-015-0527-y [PubMed]
  • 45. Zhang Z. Multiple imputation with multivariate imputation by chained equation (MICE) package. Ann Transl Med. 2016; 4:30. https://doi.org/10.3978/j.issn.2305-5839.2015.12.63 [PubMed]
  • 46. Li C, Wang Z, Hardy T, Huang Y, Hui Q, Crusto CA, Wright ML, Taylor JY, Sun YV. Association of Obesity with DNA Methylation Age Acceleration in African American Mothers from the InterGEN Study. Int J Mol Sci. 2019; 20:4273. https://doi.org/10.3390/ijms20174273 [PubMed]
  • 47. Beach SR, Dogan MV, Lei MK, Cutrona CE, Gerrard M, Gibbons FX, Simons RL, Brody GH, Philibert RA. Methylomic Aging as a Window onto the Influence of Lifestyle: Tobacco and Alcohol Use Alter the Rate of Biological Aging. J Am Geriatr Soc. 2015; 63:2519–25. https://doi.org/10.1111/jgs.13830 [PubMed]
  • 48. Karlsson IK, Gatz M, Arpawong TE, Dahl Aslan AK, Reynolds CA. The dynamic association between body mass index and cognition from midlife through late-life, and the effect of sex and genetic influences. Sci Rep. 2021; 11:7206. https://doi.org/10.1038/s41598-021-86667-4 [PubMed]
  • 49. Valentine G, Sofuoglu M. Cognitive Effects of Nicotine: Recent Progress. Curr Neuropharmacol. 2018; 16:403–14. https://doi.org/10.2174/1570159X15666171103152136 [PubMed]