Skip to main content

Associations of epigenetic age acceleration at birth and age 12 years with adolescent cardiometabolic risk: the HOME study

Abstract

Background

Cardiometabolic risk factors among youth are rising. Epigenetic age acceleration, a biomarker for aging and disease-risk, has been associated with adiposity in children, but its association with other cardiometabolic risk markers remains understudied. We employed data from the Health Outcomes and Measures of the Environment (HOME) study, a prospective pregnancy and birth cohort in the greater Cincinnati metropolitan area, to examine whether accelerated epigenetic age at birth as well as accelerated epigenetic age and faster pace of biological aging at age 12 years were associated with higher cardiometabolic risk in adolescents.

Results

After adjusting for potential confounders, including estimated cell type proportions, epigenetic gestational age acceleration at birth, derived from the Bohlin, Knight, and Haftorn clocks using cord blood DNA methylation data, was not associated with cardiometabolic risk z-scores or individual cardiometabolic risk score components (visceral fat, leptin to adiponectin ratio, HOMA-IR, triglycerides to HDL-C ratio, HbA1c, or systolic blood pressure) at age 12 years. We also did not observe any associations of epigenetic age acceleration, calculated with Horvath’s skin and blood, Hannum’s, and Wu’s epigenetic clocks using peripheral blood at age 12 years, with these same cardiometabolic risk markers. In contrast, faster pace of biological aging was associated with higher cardiometabolic risk [βs (95% CIs)] cardiometabolic risk score 0.25 (0.07, 0.42); visceral fat 0.21 (0.05, 0.38); and hemoglobin A1c 0.23 (0.05, 0.41) per standard deviation increase in pace of biological aging. Faster pace of biological aging was also positively associated with systolic blood pressure, triglycerides to HDL-C ratio, HOMA-IR, and leptin to adiponectin ratio, although these associations were not statistically significant.

Conclusions

Our findings provide evidence that faster pace of biological aging was associated with higher cardiometabolic risk score, visceral fat, and HbA1c at age 12 years. Further research is needed to determine whether these associations persist from adolescence through adulthood.

Background

The prevalence of factors associated with elevated cardiometabolic risk among youth is rising [1,2,3]. Among United States (U.S.) youth ages 12–19 years, it is estimated that nearly 25% of have elevated blood glucose, 18% elevated waist circumference, 20% low high-density lipoprotein cholesterol (HDL-C), 5% elevated triglycerides, and 4% metabolic syndrome [4]. These risk factors track from childhood to adulthood and contribute to adult cardiovascular disease, the leading cause of global deaths, thus underscoring the critical need for preventive measures and intervention strategies in early life [5, 6].

Several social, behavioral, and environmental factors have been associated with cardiometabolic risk across the lifespan [6]. A potential biological mechanism underlying these associations is epigenetic aging. Epigenetic aging is a biomarker of biological age increasingly used to predict age-related disease risk [7]. It is estimated from the DNA methylation at select CpG sites across the genome using epigenetic clocks [7]. Several epigenetic clocks, many which use different CpG sites, have been developed to measure epigenetic gestational age (EGA) in neonates and quantify epigenetic age (EA) in children and adults [8,9,10,11,12,13]. These epigenetic clocks were constructed using different DNA methylation profile microarrays and training populations, as well as various tissue types [8,9,10,11,12,13,14]. Moreover, different generations of the clocks have emerged. First-generation clocks were derived from CpG sites correlated with chronological age and are used to calculate epigenetic age acceleration (EAA; residual obtained from regressing chronological age on estimated epigenetic age) [8,9,10,11,12,13], while a recently developed third-generation clock, DunedinPACE, quantifies the pace of biological aging and was trained on 19 longitudinal biomarkers linked to the pulmonary, cardiovascular, metabolic, renal, hepatic, periodontal, and immune systems [15]. Multiple epigenetic clocks have been associated with several disease risk factors (e.g., smoking and low socioeconomic status) and adverse health outcomes (e.g., cardiovascular disease, cancer, and mortality), but these associations vary with each clock [16,17,18,19].

Among children and adolescents, results from prior studies examining the associations of EAA and cardiometabolic risk factors have been inconsistent and the associations of pace of biological aging with cardiometabolic risk have been scarce [20,21,22,23,24,25]. While two studies found that epigenetic gestational age acceleration (EGAA) estimated via the Bohlin clock was associated with lower likelihood of rapid weight gain during childhood [20, 21] and lower weight at 10 years of age [21], another study reported that EAA at birth, estimated via the Horvath clock, was associated with greater fat mass and faster gain in weight and BMI during childhood and adolescence [22]. However, this same study found that EAA at age 7 years was not associated with adiposity during childhood and adolescence [22]. In late adolescence, Hannum EAA at age 17 years was associated with greater lean mass, BMI, and waist circumference, while Horvath EAA was associated with higher leptin levels and HOMA-IR [23]. For blood pressure, both positive and inverse associations with EGAA at birth and EAA in childhood and adolescence have been reported, although the 95% confidence intervals included the null [23, 25]. Only one study to date has examined the associations of EAA with lipid biomarkers among children, with findings indicating that EAA at age 17 years was both inversely and positively associated, although not significantly, with total cholesterol, triglycerides, HDL-C, and low-density lipoprotein cholesterol (LDL-C) in late adolescence [23]. Further, the associations of pace of biological aging with cardiometabolic risk factors in adolescence are unknown.

Therefore, using six cardiometabolic risk biomarkers assessed at age 12 years, we examined whether accelerated EGA at birth as well as accelerated EA and faster pace of biological aging at age 12 years were associated with higher cardiometabolic risk in adolescence. Because of the variability in the design of the epigenetic clocks, we also assessed the deviations and correlations between the estimated epigenetic age and the observed age (gestational or chronological) across three widely used EGA clocks and three EA clocks.

Methods

Study design and participants

We used data from the Health Outcomes and Measures of the Environment (HOME) Study, a prospective pregnancy and birth cohort in the greater Cincinnati metropolitan area [26, 27]. Between March 2003 and January 2006, the HOME study enrolled 468 out of 1263 eligible pregnant women from nine prenatal practices affiliated with three hospitals in the Cincinnati, Ohio region [26]. Pregnant women were eligible for inclusion if they: 1) were at least age 18 years; 2) were 16 ± 3 weeks of gestation; 3) were fluent in English; 4) resided in the study region in a home built prior to 1978 (i.e., when the U.S federal government banned consumer used of lead-based paint) with goal of including children with increased lead exposure given that one of the aims of the HOME study was to test the efficacy of lead hazard controls on children’s lead exposure and health) [28]; 5) not residing in a mobile or trailer home; 6) expected to continue prenatal care and deliver at a collaborating hospital; 7) were HIV negative; 8) were not taking seizure, thyroid, or chemotherapy/radiation medications; 9) had no diagnosis of diabetes, bipolar disorder, schizophrenia, or cancer requiring radiation or chemotherapy [26].

Of the 468 enrolled women, 412 women and their children were eligible for the 12-year visit after 67 women dropped out prior to delivery, but 11 of them re-enrolled at a later follow-up (Supplemental Fig. 1) [27]. Women commonly dropped out because of time commitment and family member’s unwillingness to participate [26]. Among 242 singleton children who completed the 12-years follow-up visit [27], 181 singleton children had epigenetic gestational age at birth, at least one cardiometabolic risk biomarker at age 12 years, and complete covariate data. A total of 183 singleton children had epigenetic age at age 12 years, at least one cardiometabolic risk biomarker, and complete data on relevant covariates.

HOME study protocols were approved by the Institutional Review Boards (IRBs) at the Cincinnati Children’s Hospital Medical Center (CCHMC) and collaborating delivery hospitals [26]. The Centers for Disease Control and Prevention and Brown University deferred to the CCHMC IRB as the IRB of record. Women provided written informed consent for themselves and their children [26]. At the 12-year visit, adolescents provided written informed assent [27].

DNA methylation

We collected venous blood samples at delivery (cord blood) and age 12 years, and stored them at − 80 °C prior to DNA methylation (DNAm) analysis as previously described [29]. Briefly, trained laboratory technicians quantified DNAm at > 850,000 CpG sites per sample using the Illumina EPIC Bead Chip (EPIC) v1.0 array [30]. Quality control and normalization of the DNAm data were performed using the minfi Bioconductor package in R (v1.48.0) [29, 31, 32]. We excluded samples with more than 5% high detection p-values (p-value > 1 × 10–6) and filtered out probes with at least one high detection p-value (p-value > 1 × 10–6), known to be cross-reactive, having single nucleotide polymorphisms, or mapping to sex chromosomes [29, 33]. A total of 669,622 probes were used for further analysis. We used normal-exponential convolution using out-of-band probes (Noob) for background and dye bias correction [34], Beta MIxture Quantile (BMIQ) normalization for probe design bias correction [35], and ComBat for batch effect adjustment [36].

Epigenetic age

We used the preprocessed DNAm to estimate: 1) epigenetic gestational age (EGA) at birth using three cord-blood-based EGA clocks—Bohlin, Knight, and Haftorn [8,9,10]; 2) epigenetic age (EA) at age 12 years using the Horvath skin and blood, Hannum, and Wu clocks [11,12,13]; and 3) pace of biological aging with the recently developed DunedinPACE clock [15]. We have selected these epigenetic clocks because they are frequently used in other studies, allowing us to compare our results. Moreover, these selected clocks were trained in tissue types (i.e., cord blood for EGA clocks and blood for EA clocks) consistent with the specimens collected in our cohort [8,9,10,11,12,13, 15]. Further, because we quantified DNA methylation using the Illumina EPIC array, we included at least one clock in each period (i.e., birth and age 12 years) and generation (i.e., first and third generation) trained using DNA methylation from the EPIC array. Blood-based pediatric-specific epigenetic clocks are limited [14], and as a result we examined epigenetic clocks at age 12 years even if the training population did not cover this age range. We also considered Levine’s PhenoAge, a second-generation epigenetic clock trained on chronological age and 9 clinical biomarkers, but excluded it due to a high median absolute error (medAE = 25.3 years) in our cohort [37].

Details of each epigenetic clock [8,9,10,11,12,13, 15] and comparison of their characteristics have been previously reported [14, 16]. Briefly, the first-generation epigenetic clocks—Bohlin, Knight, Haftorn, Horvath, Hannum, and Wu—are predictors of age and were developed by regressing gestational or chronological age on DNAm at CpG sites primarily via elastic net regression [8,9,10,11,12,13]. In contrast, DunedinPACE, a third-generation clock and predictor of the pace/rate of aging, was based on 19 biomarkers repeatedly measured from young to middle adulthood [15]. Across the epigenetic clocks, variability in the design, demographics of populations used for training the data, and selected CpG sites resulted in little to no overlap in CpGs between the clocks (Supplemental Fig. 2) [8,9,10,11,12,13, 15]. Therefore, we used multiple epigenetic clocks to assess consistency in the effect estimates from the regression analyses.

Epigenetic age (EGA and EA) and their respective acceleration measures (EGAA and EAA) were computed using the methylclock R package (v1.8.0) [38]. EGAA and EAA were estimated from the residuals obtained after regressing EGA on gestational age and EA on chronological age (methylclock v1.8.0) [38]. Positive values of EGAA and EAA indicate accelerated epigenetic aging, epigenetic age that is greater than chronological age. Due to the known associations of age acceleration with cell type proportions [39], we calculated intrinsic measures of EGAA and EAA by adjusting for estimated cell type proportions. We used the reference-based approach developed by Gervin et al. [40] for cord blood and by Salas et al. [41] for adult peripheral blood, as implemented by methylclock to estimate cell type proportions [38]. Pace of biological aging, defined as an average rate of change in one year of biological aging per year of chronological aging [15], was estimated using the DunedinPACE algorithm and the intrinsic measure was separately calculated in a linear regression model adjusted for estimated cell type proportions; neutrophils were excluded due to multicollinearity [41]. Missing CpG probes required for estimating each first-generation clock were recovered using the k-nearest neighbors (k = 10) imputation approach in methylclock consistent with imputation method used by the Knight, Wu, and Hannum epigenetic clocks [9, 12, 13, 38], while missing probes for DunedinPACE were rescued using the mean imputation procedure in the PACEProjector function of the DunedinPACE algorithm (v0.99.0) [15].

Adolescent cardiometabolic risk

We measured several cardiometabolic risk biomarkers at the 12-years visit [27]. Trained laboratory technicians at the Cincinnati Children’s Hospital Medical Center Clinical Translational Research Center Core Laboratory quantified serum insulin, glucose, HDL-C, triglycerides, adiponectin, leptin, and hemoglobin A1C (HbA1c) concentrations using immunoassays conducted on fasting blood samples. We obtained three sitting blood pressure measurements at 1-min intervals using a Dinamap Pro100 automated monitor and used the average of the second and third measures for analyses. Cross-sectional area of fat inside the abdominal cavity (cm2) was determined by dual-energy X-ray absorptiometry (DXA) using a Horizon densitometer (Hologic) [27]. We estimated the ratios of triglycerides to HDL-C and leptin to adiponectin. The homeostatic model assessment for insulin resistance (HOMA-IR) was calculated as: fasting insulin (mIU/L) × fasting glucose (mg/dL)/405 [42].

We derived a summary cardiometabolic risk score as previously described [43]. In brief, we calculated age and sex standardized z-scores for HOMA-IR, HbA1c, triglycerides to HDL-C ratio, leptin to adiponectin ratio, and cross-sectional area of fat inside the abdominal cavity (visceral fat) from the studentized residuals of the linear regression for each component [43, 44]. Prior to standardization, we log2-transformed HOMA-IR, triglycerides to HDL-C ratio, leptin to adiponectin ratio, and visceral fat, because they were right-skewed. Systolic blood pressure was first nationally standardized according to sex, age, and height measures, and subsequently standardized internally to ensure comparability with the individual components of the cardiometabolic risk score [45, 46]. The summed standardized z-scores for HOMA-IR, triglycerides to HDL-C ratio, systolic blood pressure, leptin to adiponectin ratio, and visceral fat comprised the cardiometabolic risk score, with a higher score indicative of higher cardiometabolic risk [43].

Covariates

We used previous literature and directed acyclic graphs (DAG) to identify potential confounders associated with epigenetic age at birth or age at 12 years and early adolescent cardiometabolic risk (Supplemental Figs. 3 and 4). Trained research staff abstracted data from medical charts and used standardized computer-assisted interviews to obtain information on maternal age, maternal education, household income, gestational age at delivery, birthweight, as well as the child’s age, sex, and race/ethnicity [26, 27]. We calculated pre-pregnancy body mass index (BMI) from self-reported height and weight [47]. For mothers with missing self-reported weight, we imputed their weight with super learning [47, 48]. With the SuperLearner (v2.0–22) R package, we used a tenfold cross-validation approach to obtain weighted predictions based on a subset of candidate algorithms [47]. Prenatal tobacco smoke exposure was assessed using the average log10-transformed serum cotinine concentrations [49]. Serum cotinine concentrations were measured using high-performance liquid chromatography-tandem mass spectroscopy from blood samples collected from mothers via venipuncture at 16 or 26 weeks of gestation [49]. To obtain estimated cell type proportions, we used the reference-based algorithm developed by Gervin et al. [40] for cord blood and by Salas et al. [41] for adult peripheral blood. At the 12-year study visit, adolescents used images illustrating the Tanner stages of development to self-assess their pubertal status based on pubic hair growth [27, 50].

Statistical analyses

Using descriptive statistics, we summarized the distribution (median and interquartile range for continuous variables, frequency and proportion for categorical variables) of sociodemographic and perinatal characteristics of the mother and children in the study sample. We tested differences in characteristics for included and excluded participants using Pearson’s chi-squared or Fisher’s exact test for categorical variables, and Wilcoxon rank sum test for continuous measures. We reported the mean and standard deviation, overall and by covariates, of EGAA for the Bohlin, Knight, and Haftorn epigenetic clocks; EAA for the Wu, Hannum, and Horvath’s skin and blood epigenetic clocks; and pace of biological aging for DunedinPACE.

To evaluate the performance of the first-generation epigenetic clocks, we calculated Spearman correlations and median absolute error (medAE) between epigenetic age and age at birth or 12 years for each first-generation epigenetic clock.15 The median absolute error was calculated as the median absolute difference between epigenetic age and chronological age. We additionally used Spearman rank correlations to investigate the relationship between EGAA at birth and EAA at age 12 years.

To estimate the association of intrinsic EGAA at birth and EAA at age 12 years with cardiometabolic risk, we fit separate multivariable linear regression models with EGAA or EAA, respectively, as predictors and with cardiometabolic risk score, each of the individual cardiometabolic risk components (i.e., HOMA-IR, triglycerides to HDL-C ratio, systolic blood pressure, leptin to adiponectin ratio, and visceral fat), and HbA1c at age 12 years as outcomes. We standardized EGAA, EAA, and pace of biological aging using z-scores (i.e., mean and SD were set to 0 and 1, respectively) to ensure comparability of effect sizes across the different epigenetic clocks. We adjusted all models for pre-pregnancy BMI (continuous), gestational serum cotinine (< 0.015 ng/mL as unexposed; 0.015–3 ng/mL as secondhand exposure; ≥ 3 ng/mL as active smoking), maternal age at delivery (continuous, years), household income (continuous, USD/year), child age at outcome measurement (continuous, years), child sex (male; female), and child race/ethnicity (non-Hispanic, white; non-Hispanic, black; other). Estimated cell type proportions were accounted in the intrinsic accelerated epigenetic age measures of all first-generation clocks; and thus, were not adjusted for in the linear regression models [40, 41]. For pace of biological aging, we adjusted for estimated cell type proportions in the linear regression models. We excluded stage of pubertal development because it was a potential intermediate in the relationship between EAA and cardiometabolic risk in our DAG (Supplemental Fig. 4) [51, 52].

To test the robustness of our results, we conducted two sensitivity analyses: 1) we additionally adjusted models for birthweight in the analyses examining the associations of cardiometabolic risk with intrinsic EAA and pace of biological aging but excluded it in the associations of EGAA with cardiometabolic risk because it is a potential intermediate [22, 53]; and 2) we evaluated the associations of cardiometabolic risk with EGAA, EAA, and pace of biological aging without adjusting for cell type proportions.

Results were considered statistically significant when the two-sided p-value < 0.05. We did not adjust for multiplicity because we aimed to explore associations of epigenetic age measures with cardiometabolic risk biomarkers. We conducted statistical analyses using R (version 4.3.2) [54].

Results

Characteristics of study participants

The median gestational age at birth among the 181 singleton children included in our analysis was 39.4 weeks (IQR = 38.7–40.0), with 55% being female and 60% non-Hispanic white (Supplemental Table 1; Table 1). Nearly two-thirds of mothers had an annual household income > $40,000 during pregnancy (62%), and 60% had gestational serum cotinine concentrations indicative of secondhand tobacco smoke exposure (0.015 – 3 ng/mL, 60%). At the 12-year visit, the median age of children included in the analysis was 12.3 years (IQR = 11.9–12.8). While no significant differences in characteristics were detected between the included and excluded participants at birth and the 12-years visit, a higher proportion of children included in this analysis were non-Hispanic white with mothers who were > 25–35 years old at delivery (Supplemental Table 1).

Table 1 Epigenetic gestational age acceleration (EGAA) at birth across covariates among children from the HOME study (2003–2006), N = 181

Differences in EGAA and EAA across sociodemographic characteristics

The children in our sample were on average decelerated at birth (Table 1i.e., lower EGAA than expected), and had average EAA (i.e., EAA near zero) and slower pace of biological aging than expected at age 12 years (Table 2). Across all clocks, there was variability in mean EGAA and EAA across levels of some participants’ characteristics. For instance, male children had higher acceleration than females based on the Hannum clock, while female children had higher acceleration based on skinHorvath clock. In contrast, children who were non-Hispanic black and whose mothers had a pre-pregnancy BMI ≥ 30 had the lowest mean EGAA at birth based on all three EGA clocks (Table 1). Children with mothers > 35 years of age at delivery and with pre-pregnancy BMI ≥ 30 had the lowest mean EAA on the Wu, Hannum, and skinHorvath clocks (Table 2). Mean pace of biological aging was higher with lower annual household income and higher gestational serum cotinine concentrations (Table 2).

Table 2 Epigenetic age acceleration (EAA) and pace of biological aging at age 12 years across covariates among children from the HOME study (2003–2006), N = 183

Differences between the epigenetic clocks in the correlations and deviations of epigenetic age and observed age

Gestational age at delivery (weeks) was more strongly correlated with the Bohlin EGA clock (ρ = 0.54) than the Haftorn EGA (ρ = 0.46, p < 0.001) and Knight clock (ρ = 0.36, p < 0.001) (Fig. 1). The median absolute error for the EGA clocks ranged from 0.80 to 1.22 weeks, with Bohlin and Haftorn having the lowest medAE and Knight having the highest medAE. Although the correlations of chronological age (years) with EA at age 12 years were weak (ρ range: 0.15–0.25), they were consistently positive across all EA clocks (Fig. 2). Among the EA clocks at age 12 years, Hannum EA had the lowest median absolute error (medAE = 2.48 years), whereas skinHorvath EA had the highest (medAE = 5.39 years).

Fig. 1
figure 1

Spearman rank correlations of gestational age and epigenetic gestational age at birth (weeks) among singleton children from the HOME study (N = 181)

Fig. 2
figure 2

Spearman rank correlations of chronological age and epigenetic age (years) at age 12 years among singleton children from the HOME study (N = 183)

Correlations of EGAA at birth with EAA and pace of biological aging at age 12 years

Among the 152 singleton children with DNA methylation data at birth and age 12 years, we observed consistent positive correlations between EGAA at birth and EAA at age 12 years (Fig. 3). EGAA at birth was inversely correlated with pace of biological aging (Bohlin ρ = − 0.09, p = 0.524; Knight ρ = − 0.17, p = 0.038; Haftorn ρ = − 0.05, p = 0.263), although the correlations of pace of biological aging with Bohlin EGAA and Haftorn EGAA were not statistically significant.

Fig. 3
figure 3

Spearman rank correlations of epigenetic gestational age acceleration (EGAA) at birth with epigenetic age acceleration (EAA) and pace of biological aging at age 12 years, among singleton children from the HOME study with DNA methylation data at birth and age 12 years (N = 152)

Adjusted associations of EGAA at birth and EAA/pace of biological aging at age 12 years with cardiometabolic risk in adolescence

In covariate-adjusted analysis, EGAA at birth had consistent positive associations with leptin to adiponectin ratio but inverse associations with visceral fat and systolic blood pressure across the 3 EGA clocks, although the 95% confidence intervals included the null (Fig. 4). Wu, Hannum, and Horvath EAA at age 12 years were consistently associated with lower cardiometabolic risk score, HOMA-IR, triglycerides to HDL-C ratio, and systolic blood pressure, but the 95% confidence intervals included the null. We did not observe any other consistent patterns of associations of EGAA at birth with cardiometabolic risk score, triglycerides to HDL-C ratio, and HOMA-IR; and of EAA at age 12 years with visceral fat and HbA1c.

Fig. 4
figure 4

Adjusted associations of intrinsic epigenetic gestational age acceleration (EGAA) at birth, and epigenetic age acceleration (EAA) and pace of biological aging at 12 years with adolescent cardiometabolic risk

Faster pace of biological aging was associated with higher cardiometabolic risk score, leptin to adiponectin ratio, visceral fat, HOMA-IR, HbA1c, triglycerides to HDL-C ratio, and systolic blood pressure at age 12 years (Fig. 4). After adjusting for covariates, as the pace of biological aging increased by one SD (0.08), the cardiometabolic risk score, visceral fat, and HbA1c increased, on average, by 0.25 SD (95% CI 0.07–0.42), 0.21 SD (95% CI 0.05–0.38), and 0.23 (95% CI 0.05–0.41), respectively. Further adjustment for birthweight in sensitivity analyses yielded similar results (Supplemental Fig. 5). When excluding the estimated cell type proportions in sensitivity analyses, we observed that the associations of pace of biological aging with cardiometabolic risk score, visceral fat, leptin to adiponectin ratio, and HOMA-IR were strengthened while the associations with triglycerides to HDL-C ratio, HbA1c, and systolic blood pressure were attenuated (Supplemental Fig. 6).

Each circle represents the effect estimate of each association and the error bars represent the corresponding 95% confidence intervals. Associations with p < 0.05 are displayed in purple.

Adjusted standardized difference in cardiometabolic risk score or individual cardiometabolic risk component at age 12 years per 1 SD increase in EGAA at birth or EAA/pace of biological aging at age 12 years. CI = Confidence Interval; Intrinsic = adjusted for estimated cell type proportions. Adjusted for pre-pregnancy BMI (continuous, kg/m2), gestational serum cotinine (< 0.015 ng/mL as unexposed; 0.015–3 ng/mL as secondhand exposure; ≥ 3 ng/mL as active smoking), maternal age at delivery (continuous, years), household income (continuous, USD/year), child age (continuous, years), child sex (male; female), child race/ethnicity (non-Hispanic White; non-Hispanic Black; other), and estimated cell type proportions (at birth: nucleated red blood cells, CD8+ T-cells, CD4+ T-cells, natural killer cells, B-cells, monocytes, and granulocytes; at 12-year visit: eosinophils, basophils, monocytes, naïve and memory B cells, naïve and memory CD4+ and CD8+ T cells, natural killer, and T regulatory cells; neutrophils were excluded due to multicollinearity).

Discussion

Our results indicate that faster pace of biological aging, as measured by DunedinPACE, was associated with elevated overall cardiometabolic risk score, as well as higher visceral fat and hemoglobin A1c. Adjusting for estimated cell type proportions, given that epigenetic clocks are known to vary with the immune cell composition [39], strengthened the associations of pace of biological aging with some cardiometabolic risk biomarkers including triglycerides to HDL-C ratio, HbA1c, and systolic blood pressure, but attenuated the associations with others, such as cardiometabolic risk score and visceral fat. Currently, only one study has investigated the relationship of pace of biological aging with biomarkers of cardiometabolic risk in adolescents [55]. In the Child Health Study, a cohort of 439 children ages 8–13 years with and without history of maltreatment based on the Pennsylvania’s Statewide Child Welfare Information System, higher BMI z-score was associated with faster Dunedin methylation pace of biological aging (DunedinPoA) [55]. However, the analysis assessed the influence of BMI on pace of biological aging, and not vice versa.

While we found no evidence of significant associations between accelerated epigenetic measures and early adolescent cardiometabolic risk markers, our data was compatible with EGAA at birth being positively associated with leptin to adiponectin ratio but inversely associated with systolic blood pressure. Among 1115 mother–offspring pairs of the Generation R Study, Bohlin EGAA at birth was also inversely associated with systolic blood pressure at age 10 years, although the results included the null [25]. The positive association of skinHorvath EAA with systolic blood pressure at age 10 years in the Generation R Study as well as Hannum EEAA and Horvath IEAA with systolic blood pressure at age 17 years in the Raine Study contrasted our data, which was most compatible with an inverse association of EAA with systolic blood pressure at age 12 years [23, 25]. Further, in two studies examining the relationship of age acceleration at birth with adiposity in children, Bohlin EGAA at birth was associated with elevated weight from birth to 9 months but lower weight at age 10 years, which was consistent with our findings of Bohlin EGAA with visceral fat at age 12 years [20, 21], but differed from findings of another study reporting that Horvath EAA at birth was associated with higher fat mass, and rapid gain in weight and BMI during childhood and adolescence [22].

Compared with the consistently positive associations of the cardiometabolic risk markers with pace of biological aging—calculated from third-generation clock DunedinPACE, the associations of EGAA and EAA, derived from first-generation clocks, with the cardiometabolic risk biomarkers reflected both positive and inverse associations. Differences in these associations may be due to how the first-generation clocks were trained compared to the third-generation clock DunedinPACE. DunedinPACE was trained on 19 clinical biomarkers with the goal of quantifying the rate of aging and capturing decline in organ system function [15]. Thus, it may be more sensitive at assessing disease risk. In contrast, the first-generation clocks were exclusively trained on gestational or chronological age, where accelerated epigenetic aging measures represent deviations between epigenetic and chronological age [8,9,10,11,12,13, 15].

We evaluated the deviations and correlations of estimated epigenetic age and observed age across six first-generation epigenetic clocks, three EGA clock at birth, and three EA clocks during early adolescence in our study due to variability in their design and prior studies showing that estimates may vary based on the clock [16, 20, 56,57,58,59,60,61]. In our study, the Bohlin clock had the strongest correlation and lowest median absolute error between gestational age and EGA at birth. These findings align with prior studies indicating that the Bohlin clock more accurately estimates EGA at birth compared with the Knight clock [20, 56,57,58] and Haftorn clock [59]. However, in some cohorts, the Bohlin and Haftorn clocks performed similarly [60, 61]. At age 12 years, we observed low correlations between chronological age and EA, with the Hannum clock displaying the strongest correlation. Although the Wu and skinHorvath clocks included children in their training populations, the broad age ranges spanning multiple stages of childhood—and adulthood for skinHorvath—may have affected the precision of their estimates and influenced their performance in our cohort [11, 14]. Taken together, differences in performance across the clocks could stem from the variability in how the epigenetic clocks were developed, including the demographics and sample size of the training population and selected CpGs, and results for them should be interpreted in the context of these differences [8,9,10,11,12,13].

Our findings suggest that pace of biological aging could be used as a biomarker to monitor cardiometabolic health as early as childhood and adolescence. Faster rate of biological aging has implications for increased risk of disease (e.g., cardiovascular disease) and mortality [15, 16]. Accelerated epigenetic aging in adolescence has been associated with higher risk of cardiovascular disease in middle adulthood, which may be consequential to premature mortality [23]. Therefore, identifying individuals aging faster could inform early intervention strategies and prevent the onset of adverse health outcomes in adulthood.

Our study has some strengths including the assessment of age acceleration during two critical windows of development in childhood using multiple epigenetic clocks at birth and at age 12 years as well as the inclusion of various cardiometabolic risk biomarkers in our analysis. However, this study has some limitations. First, the cross-sectional nature of the associations at the 12-years visit precludes us from determining their directionality and extrapolating causality. Prior cross-sectional analyses among children have found inverse associations of cardiometabolic risk factors with epigenetic age accelerations and pace of biological aging [55, 62]. Our findings should be validated in longitudinal analyses examining the pace of biological aging with changes in cardiometabolic risk through childhood and adolescence. Second, we tested the associations of multiple epigenetic-based measures with several cardiometabolic risk outcomes, which may have resulted in findings due to random chance. Therefore, we have focused on the patterns of associations. Third, estimates of epigenetic age using adult-specific clocks such as the Hannum clock may have reduced precision given that they were trained exclusively on adult populations. However, in our cohort, the Hannum epigenetic clock had stronger correlation and lower median absolute error than the pediatric-specific Wu epigenetic clock. Fourth, we assessed cardiometabolic risk using a composite score with equal weights allocated to each individual cardiometabolic risk component. Notably, findings were generally consistent across the individual components of the cardiometabolic risk score.

Further, there is potential for selection bias due to the exclusion of adolescents with missing data or from loss to follow-up. Sociodemographic characteristics of participants who returned for follow-up visits were comparable to those of the full cohort [27]. In addition, not all eligible pregnant women enrolled into the study, which may have induced self-selection bias. Reporting bias may have been introduced with self-reported data such as maternal weight and height. There may be residual confounding by factors such as race/ethnicity given the heterogeneity of individuals within racial categories. Finally, given our study eligibility criteria, our cohort may not be representative of births in the study region; thus, our findings may not be generalizable to other populations [27].

Conclusions

Our results suggest that faster pace of biological aging may influence cardiometabolic risk, visceral fat, and hemoglobin A1c during early adolescence. However, we did not find evidence that accelerated EGA at birth or EA at age 12 years were associated with early adolescent cardiometabolic risk. Future research should further explore this relationship using repeated measures of DunedinPACE and cardiometabolic risk markers at different stages during adolescence and adulthood.

Availability of data and materials

Data are available upon request. The HOME study Principal Investigators welcome new collaborations with other investigators and actively engage in collaborative data-sharing projects. Interested investigators should visit https://homestudy.research.cchmc.org/contact or contact JMB (joseph_braun_1@brown.edu) and KY (kimberly.yolton@cchmc.org) to obtain additional information about The HOME study, discuss collaborative opportunities, and request a project proposal form. The HOME study Protocol Review Committee reviews proposed research projects to ensure that they do not overlap with extant projects and are an efficient use of scarce resources (e.g., biospecimens).

References

  1. Noubiap JJ, Nansseu JR, Lontchi-Yimagou E, Nkeck JR, Nyaga UF, Ngouo AT, Tounouga DN, Tianyi FL, Foka AJ, Ndoadoumgue AL, Bigna JJ. Global, regional, and country estimates of metabolic syndrome burden in children and adolescents in 2020: a systematic review and modelling analysis. Lancet Child Adolesc Health. 2022;6:158–70.

    Article  PubMed  Google Scholar 

  2. Bendor CD, Bardugo A, Pinhas-Hamiel O, Afek A, Twig G. Cardiovascular morbidity, diabetes and cancer risk among children and adolescents with severe obesity. Cardiovasc Diabetol. 2020;19:79.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Skinner AC, Ravanbakht SN, Skelton JA, Perrin EM, Armstrong SC. Prevalence of obesity and severe obesity in US children, 1999–2016. Pediatrics 2018;141.

  4. Liu J, Ma J, Orekoya O, Vangeepuram N, Liu J. Trends in Metabolic Syndrome Among US Youth, From 1999 to 2018. JAMA Pediatr. 2022;176:1043–5.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Di Cesare M, Perel P, Taylor S, Kabudula C, Bixby H, Gaziano TA, McGhie DV, Mwangi J, Pervan B, Narula J, et al. The heart of the world. Global Heart 2024;19.

  6. Drozdz D, Alvarez-Pitti J, Wojcik M, Borghi C, Gabbianelli R, Mazur A, Herceg-Cavrak V, Lopez-Valcarcel BG, Brzezinski M, Lurbe E, Wuhl E. Obesity and cardiometabolic risk factors: from childhood to adulthood. Nutrients. 2021;13:4176.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Horvath S, Raj K. DNA methylation-based biomarkers and the epigenetic clock theory of ageing. Nat Rev Genet. 2018;19:371–84.

    Article  CAS  PubMed  Google Scholar 

  8. Bohlin J, Haberg SE, Magnus P, Reese SE, Gjessing HK, Magnus MC, Parr CL, Page CM, London SJ, Nystad W. Prediction of gestational age based on genome-wide differentially methylated regions. Genome Biol. 2016;17:207.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Knight AK, Craig JM, Theda C, Baekvad-Hansen M, Bybjerg-Grauholm J, Hansen CS, Hollegaard MV, Hougaard DM, Mortensen PB, Weinsheimer SM, et al. An epigenetic clock for gestational age at birth based on blood methylation data. Genome Biol. 2016;17:206.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Haftorn KL, Lee Y, Denault WRP, Page CM, Nustad HE, Lyle R, Gjessing HK, Malmberg A, Magnus MC, Naess O, et al. An EPIC predictor of gestational age and its application to newborns conceived by assisted reproductive technologies. Clin Epigenetics. 2021;13:82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Horvath S, Oshima J, Martin GM, Lu AT, Quach A, Cohen H, Felton S, Matsuyama M, Lowe D, Kabacik S, et al. Epigenetic clock for skin and blood cells applied to Hutchinson Gilford Progeria Syndrome and ex vivo studies. Aging (Albany NY). 2018;10:1758–75.

    Article  CAS  PubMed  Google Scholar 

  12. Hannum G, Guinney J, Zhao L, Zhang L, Hughes G, Sadda S, Klotzle B, Bibikova M, Fan JB, Gao Y, et al. Genome-wide methylation profiles reveal quantitative views of human aging rates. Mol Cell. 2013;49:359–67.

    Article  CAS  PubMed  Google Scholar 

  13. Wu X, Chen W, Lin F, Huang Q, Zhong J, Gao H, Song Y, Liang H. DNA methylation profile is a quantitative measure of biological aging in children. Aging (Albany NY). 2019;11:10031–51.

    Article  CAS  PubMed  Google Scholar 

  14. Fang F, Zhou L, Perng W, Marsit CJ, Knight AK, Cardenas A, Aung MT, Hivert MF, Aris IM, Goodrich JM, et al. Evaluation of pediatric epigenetic clocks across multiple tissues. Clin Epigenet. 2023;15:142.

    Article  Google Scholar 

  15. Belsky DW, Caspi A, Corcoran DL, Sugden K, Poulton R, Arseneault L, Baccarelli A, Chamarti K, Gao X, Hannon E, et al. DunedinPACE, a DNA methylation biomarker of the pace of aging. Elife 2022;11.

  16. Oblak L, van der Zaag J, Higgins-Chen AT, Levine ME, Boks MP. A systematic review of biological, social and environmental factors associated with epigenetic clock acceleration. Ageing Res Rev. 2021;69:101348.

    Article  CAS  PubMed  Google Scholar 

  17. Harris KM, Levitt B, Gaydosh L, Martin C, Meyer JM, Mishra AA, Kelly AL, Aiello AE. Sociodemographic and lifestyle factors and epigenetic aging in US young adults: NIMHD social epigenomics program. JAMA Netw Open. 2024;7:e2427889.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Irvin MR, Aslibekyan S, Do A, Zhi D, Hidalgo B, Claas SA, Srinivasasainagendra V, Horvath S, Tiwari HK, Absher DM, Arnett DK. Metabolic and inflammatory biomarkers are associated with epigenetic aging acceleration estimates in the GOLDN study. Clin Epigenet. 2018;10:56.

    Article  Google Scholar 

  19. Lin L, Kiryakos J, Ammous F, Ratliff SM, Ware EB, Faul JD, Kardia SLR, Zhao W, Birditt KS, Smith JA. Epigenetic age acceleration is associated with blood lipid levels in a multi-ancestry sample of older U.S. adults. BMC Med Genom. 2024;17:146.

    Article  CAS  Google Scholar 

  20. Alfano R, Zugna D, Barros H, Bustamante M, Chatzi L, Ghantous A, Herceg Z, Keski-Rahkonen P, de Kok TM, Nawrot TS, et al. Cord blood epigenome-wide meta-analysis in six European-based child cohorts identifies signatures linked to rapid weight growth. BMC Med. 2023;21:17.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Bright HD, Howe LD, Khouja JN, Simpkin AJ, Suderman M, O’Keeffe LM. Epigenetic gestational age and trajectories of weight and height during childhood: a prospective cohort study. Clin Epigenet. 2019;11:194.

    Article  CAS  Google Scholar 

  22. Simpkin AJ, Howe LD, Tilling K, Gaunt TR, Lyttleton O, McArdle WL, Ring SM, Horvath S, Smith GD, Relton CL. The epigenetic clock and physical development during childhood and adolescence: longitudinal analysis from a UK birth cohort. Int J Epidemiol. 2017;46:549–58.

    PubMed  PubMed Central  Google Scholar 

  23. Huang RC, Lillycrop KA, Beilin LJ, Godfrey KM, Anderson D, Mori TA, Rauschert S, Craig JM, Oddy WH, Ayonrinde OT, et al. Epigenetic age acceleration in adolescence associates with bmi, inflammation, and risk score for middle age cardiovascular disease. J Clin Endocrinol Metab. 2019;104:3012–24.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Phang M, Ross J, Raythatha JH, Dissanayake HU, McMullan RL, Kong Y, Hyett J, Gordon A, Molloy P, Skilton MR. Epigenetic aging in newborns: role of maternal diet. Am J Clin Nutr. 2020;111:555–61.

    Article  PubMed  Google Scholar 

  25. Monasso GS, Jaddoe VWV, Kupers LK, Felix JF. Epigenetic age acceleration and cardiovascular outcomes in school-age children: the generation R study. Clin Epigenet. 2021;13:205.

    Article  CAS  Google Scholar 

  26. Braun JM, Kalloo G, Chen A, Dietrich KN, Liddy-Hicks S, Morgan S, Xu Y, Yolton K, Lanphear BP. Cohort profile: the health outcomes and measures of the environment (HOME) study. Int J Epidemiol. 2017;46:24.

    PubMed  Google Scholar 

  27. Braun JM, Buckley JP, Cecil KM, Chen A, Kalkwarf HJ, Lanphear BP, Xu Y, Woeste A, Yolton K. Adolescent follow-up in the health outcomes and measures of the environment (HOME) study: cohort profile. BMJ Open. 2020;10:e034838.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Braun JM, Hornung R, Chen A, Dietrich KN, Jacobs DE, Jones R, Khoury JC, Liddy-Hicks S, Morgan S, Vanderbeek SB, et al. Effect of residential lead-hazard interventions on childhood blood lead concentrations and neurobehavioral outcomes. JAMA Pediatr. 2018;172:934.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Liu Y, Eliot MN, Papandonatos GD, Kelsey KT, Fore R, Langevin S, Buckley J, Chen A, Lanphear BP, Cecil KM, et al. Gestational perfluoroalkyl substance exposure and DNA methylation at birth and 12 years of age: A longitudinal epigenome-wide association study. Environ Health Perspect. 2022;130:37005.

    Article  CAS  PubMed  Google Scholar 

  30. Solomon O, MacIsaac J, Quach H, Tindula G, Kobor MS, Huen K, Meaney MJ, Eskenazi B, Barcellos LF, Holland N. Comparison of DNA methylation measured by Illumina 450K and EPIC BeadChips in blood of newborns and 14-year-old children. Epigenetics. 2018;13:655–64.

    Article  PubMed  PubMed Central  Google Scholar 

  31. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Wilhelm-Benartzi CS, Koestler DC, Karagas MR, Flanagan JM, Christensen BC, Kelsey KT, Marsit CJ, Houseman EA, Brown R. Review of processing and analysis methods for DNA methylation array data. Br J Cancer. 2013;109:1394–402.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Zhou W, Laird PW, Shen H. Comprehensive characterization, annotation and innovative use of Infinium DNA methylation BeadChip probes. Nucleic Acids Res. 2017;45:e22.

    PubMed  Google Scholar 

  34. Triche TJ Jr, Weisenberger DJ, Van Den Berg D, Laird PW, Siegmund KD. Low-level processing of Illumina Infinium DNA Methylation BeadArrays. Nucleic Acids Res. 2013;41:e90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Teschendorff AE, Marabita F, Lechner M, Bartlett T, Tegner J, Gomez-Cabrero D, Beck S. A beta-mixture quantile normalization method for correcting probe design bias in Illumina Infinium 450 k DNA methylation data. Bioinformatics. 2013;29:189–96.

    Article  CAS  PubMed  Google Scholar 

  36. Johnson WE, Li C, Rabinovic A. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics. 2007;8:118–27.

    Article  PubMed  Google Scholar 

  37. Levine ME, Lu AT, Quach A, Chen BH, Assimes TL, Bandinelli S, Hou L, Baccarelli AA, Stewart JD, Li Y, et al. An epigenetic biomarker of aging for lifespan and healthspan. Aging (Albany NY). 2018;10:573–91.

    Article  PubMed  Google Scholar 

  38. Pelegi-Siso D, de Prado P, Ronkainen J, Bustamante M, Gonzalez JR. methylclock: a bioconductor package to estimate DNA methylation age. Bioinformatics. 2021;37:1759–60.

    Article  CAS  PubMed  Google Scholar 

  39. Zhang Z, Reynolds SR, Stolrow HG, Chen JQ, Christensen BC, Salas LA. Deciphering the role of immune cell composition in epigenetic age acceleration: insights from cell-type deconvolution applied to human blood epigenetic clocks. Aging Cell. 2024;23:e14071.

    Article  CAS  PubMed  Google Scholar 

  40. Gervin K, Salas LA, Bakulski KM, van Zelm MC, Koestler DC, Wiencke JK, Duijts L, Moll HA, Kelsey KT, Kobor MS, et al. Systematic evaluation and validation of reference and library selection methods for deconvolution of cord blood DNA methylation data. Clin Epigenet. 2019;11:125.

    Article  Google Scholar 

  41. Salas LA, Zhang Z, Koestler DC, Butler RA, Hansen HM, Molinaro AM, Wiencke JK, Kelsey KT, Christensen BC. Enhanced cell deconvolution of peripheral blood using DNA methylation for high-resolution immune profiling. Nat Commun. 2022;13:761.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Matthews DR, Hosker JP, Rudenski AS, Naylor BA, Treacher DF, Turner RC. Homeostasis model assessment: insulin resistance and beta-cell function from fasting plasma glucose and insulin concentrations in man. Diabetologia. 1985;28:412–9.

    Article  CAS  PubMed  Google Scholar 

  43. Li N, Liu Y, Papandonatos GD, Calafat AM, Eaton CB, Kelsey KT, Cecil KM, Kalkwarf HJ, Yolton K, Lanphear BP, et al. Gestational and childhood exposure to per- and polyfluoroalkyl substances and cardiometabolic risk at age 12 years. Environ Int. 2021;147:106344.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Eisenmann JC. On the use of a continuous metabolic syndrome score in pediatric research. Cardiovasc Diabetol. 2008;7:17.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Expert Panel on Integrated Guidelines for Cardiovascular H, Risk Reduction in C, Adolescents, National Heart L Blood I. Expert panel on integrated guidelines for cardiovascular health and risk reduction in children and adolescents: summary report. Pediatrics. 2011;128(Suppl 5):S213–56.

    Google Scholar 

  46. National High Blood Pressure Education Program Working Group on High Blood Pressure in C Adolescents. The fourth report on the diagnosis, evaluation, and treatment of high blood pressure in children and adolescents. Pediatrics. 2004;114:555–76.

    Article  Google Scholar 

  47. Braun JM, Eliot M, Papandonatos GD, Buckley JP, Cecil KM, Kalkwarf HJ, Chen A, Eaton CB, Kelsey K, Lanphear BP, Yolton K. Gestational perfluoroalkyl substance exposure and body mass index trajectories over the first 12 years of life. Int J Obes (Lond). 2021;45:25–35.

    Article  CAS  PubMed  Google Scholar 

  48. Van Der Laan MJ, Polley EC, Hubbard AE. Super learner. Stat Appl Genet Mol Biol 2007;6.

  49. Braun JM, Daniels JL, Poole C, Olshan AF, Hornung R, Bernert JT, Xia Y, Bearer C, Barr DB, Lanphear BP. A prospective cohort study of biomarkers of prenatal tobacco smoke exposure: the correlation between serum and meconium and their association with infant birth weight. Environ Health. 2010;9:53.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Yayah Jones NH, Khoury JC, Xu Y, Newman N, Kalkwarf HJ, Braun JM, Lanphear B, Chen A, Cecil KM, Rose SR, Yolton K. Comparing adolescent self staging of pubertal development with hormone biomarkers. J Pediatr Endocrinol Metab. 2021;34:1531–41.

    Article  CAS  PubMed  Google Scholar 

  51. Suarez A, Lahti J, Czamara D, Lahti-Pulkkinen M, Girchenko P, Andersson S, Strandberg TE, Reynolds RM, Kajantie E, Binder EB, Raikkonen K. The epigenetic clock and pubertal, neuroendocrine, psychiatric, and cognitive outcomes in adolescents. Clin Epigenetics. 2018;10:96.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Binder AM, Corvalan C, Mericq V, Pereira A, Santos JL, Horvath S, Shepherd J, Michels KB. Faster ticking rate of the epigenetic clock is associated with faster pubertal development in girls. Epigenetics. 2018;13:85–94.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Khouja JN, Simpkin AJ, O’Keeffe LM, Wade KH, Houtepen LC, Relton CL, Suderman M, Howe LD. Epigenetic gestational age acceleration: a prospective cohort study investigating associations with familial, sociodemographic and birth characteristics. Clin Epigenet. 2018;10:86.

    Article  Google Scholar 

  54. Team RC. R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing; 2023.

    Google Scholar 

  55. Etzel L, Hastings WJ, Hall MA, Heim CM, Meaney MJ, Noll JG, O’Donnell KJ, Pokhvisneva I, Rose EJ, Schreier HMC, et al. Obesity and accelerated epigenetic aging in a high-risk cohort of children. Sci Rep. 2022;12:8328.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Daredia S, Huen K, Van Der Laan L, Collender PA, Nwanaji-Enwerem JC, Harley K, Deardorff J, Eskenazi B, Holland N, Cardenas A. Prenatal and birth associations of epigenetic gestational age acceleration in the Center for the Health Assessment of Mothers and Children of Salinas (CHAMACOS) cohort. Epigenetics. 2022;17:2006–21.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Bozack AK, Rifas-Shiman SL, Gold DR, Laubach ZM, Perng W, Hivert MF, Cardenas A. DNA methylation age at birth and childhood: performance of epigenetic clocks and characteristics associated with epigenetic age acceleration in the Project Viva cohort. Clin Epigenet. 2023;15:62.

    Article  CAS  Google Scholar 

  58. Song AY, Feinberg JI, Bakulski KM, Croen LA, Fallin MD, Newschaffer CJ, Hertz-Picciotto I, Schmidt RJ, Ladd-Acosta C, Volk HE. Prenatal exposure to ambient air pollution and epigenetic aging at birth in newborns. Front Genet. 2022;13:929416.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Deng WQ, Pigeyre M, Azab SM, Wilson SL, Campbell N, Cawte N, Morrison KM, Atkinson SA, Subbarao P, Turvey SE, et al. Consistent cord blood DNA methylation signatures of gestational age between South Asian and white European cohorts. Clin Epigenet. 2024;16:74.

    Article  CAS  Google Scholar 

  60. Xu K, Li S, Muskens IS, Elliott N, Myint SS, Pandey P, Hansen HM, Morimoto LM, Kang AY, Ma X, et al. Accelerated epigenetic aging in newborns with Down syndrome. Aging Cell. 2022;21:e13652.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Salontaji K, Haftorn KL, Sanders F, Page CM, Walton E, Felix JF, Bekkhus M, Bohlin J, Tiemeier H, Cecil CAM. Gestational epigenetic age and ADHD symptoms in childhood: a prospective, multi-cohort study. Mol Psychiatry. 2024;29:2911–8.

    Article  PubMed  Google Scholar 

  62. Raffington L, Belsky DW, Kothari M, Malanchini M, Tucker-Drob EM, Harden KP. Socioeconomic disadvantage and the pace of biological aging in children. Pediatrics. 2021;147(6):e2020024406. https://doiorg.publicaciones.saludcastillayleon.es/10.1542/peds.2020-024406.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

The authors would like to thank the HOME study participants for the time they have given to the study. This work was supported by grants from the National Institutes of Health/National Institute of Environmental Health Sciences (R01ES032836 03S1, R01 ES032386, R01 ES025214, P01 ES011261, R01 ES014575, R01 ES020349, R01 ES027224, and R01 ES033252). The content is solely the responsibility of the authors and do not necessarily represent the official views of the National Institutes of Health.

Funding

This work was supported by Grants from the National Institutes of Health/National Institute of Environmental Health Sciences (R01ES032836 03S1, R01 ES032386, R01 ES025214, P01 ES011261, R01 ES014575, R01 ES020349, R01 ES027224, and R01 ES033252). The content is solely the responsibility of the authors and do not necessarily represent the official views of the National Institutes of Health.

Author information

Authors and Affiliations

Authors

Contributions

JA and JMB contributed to the conceptualization of the study. JA led the formal analysis, created the visualizations, and wrote the original draft of the manuscript. JMB supervised this project and contributed to the methodology. GP and KK contributed to methodology and data visualizations. KC, AC, SL, BL, and KY contributed to data and funding acquisition. All authors edited, reviewed, and approved the final manuscript.

Corresponding author

Correspondence to Jennifer L. Arzu.

Ethics declarations

Ethics approval and consent to participate

The institutional review boards of Cincinnati Children’s Hospital Medical Center, the participating delivery hospitals, and the CDC approved this study. During a face-to-face visit, research assistants explained study protocols to each prospective participant and completed a checklist to ensure women were fully informed. At the 12-year visit, adolescents provided written informed assent after study protocols were explained and their questions answered.

Consent for publication

Not applicable.

Competing interests

Dr. Joseph Braun has served as an expert witness for plaintiffs involved in litigation related to PFAS-contaminated drinking water. The other authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Arzu, J.L., Kelsey, K.T., Papandonatos, G.D. et al. Associations of epigenetic age acceleration at birth and age 12 years with adolescent cardiometabolic risk: the HOME study. Clin Epigenet 16, 163 (2024). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13148-024-01779-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13148-024-01779-8

Keywords