Epigenetic age acceleration and cardiovascular outcomes in school-age children: The Generation R Study
Clinical Epigenetics volume 13, Article number: 205 (2021)
Hypertension and atherosclerosis may partly originate in early life. Altered epigenetic aging may be a mechanism underlying associations of early-life exposures and the development of cardiovascular risk factors in childhood. A discrepancy between chronological age and age predicted from neonatal DNA methylation data is referred to as age acceleration. It may either be positive, if DNA methylation age is older than clinical age, or negative, if DNA methylation age is younger than chronological age. We examined associations of age acceleration at birth (‘gestational age acceleration’), and of age acceleration at school-age, with blood pressure and with intima-media thickness and distensibility of the common carotid artery, as markers of vascular structure and function, respectively, measured at age 10 years.
This study was embedded in the Generation R Study, a population-based prospective cohort study. We included 1115 children with information on cord blood DNA methylation and blood pressure, carotid intima-media thickness or carotid distensibility. Gestational age acceleration was calculated using the Bohlin epigenetic clock, which was developed specifically for cord blood DNA methylation data. It predicts gestational age based on methylation levels of 96 CpGs from HumanMethylation450 BeadChip. We observed no associations of gestational age acceleration with blood pressure, carotid intima-media thickness or carotid distensibility at age 10 years. In analyses among children with peripheral blood DNA methylation measured at age 6 (n = 470) and 10 (n = 449) years, we also observed no associations of age acceleration at these ages with the same cardiovascular outcomes, using the ‘skin and blood clock,’ which predicts age based on methylation levels at 391 CpGs from HumanMethylation450 BeadChip.
Our findings do not provide support for the hypothesis that altered epigenetic aging during the earliest phase of life is involved in the development of cardiovascular risk factors in childhood.
Cardiovascular risk factors such as hypertension and atherosclerosis partly originate in the earliest phases of life and track into adulthood [1,2,3,4]. Prenatal exposures, such as maternal smoking, have been associated with higher blood pressure and carotid intima-media thickness, and lower carotid distensibility in childhood [5, 6]. Arterial intima-media thickness and distensibility are related, noninvasive markers of arterial structure and function that predict cardiovascular events in adults [7, 8]. In addition to prenatal exposures, postnatal early-life exposures, such as higher weight, have been associated with suboptimal blood pressure, intima-media thickness and distensibility at school-age [9,10,11].
An underlying mechanism for the associations of early-life exposures with cardiovascular health later in childhood may be altered epigenetic aging [12, 13]. In recent years, multiple ‘epigenetic clocks’ have been developed [12,13,14,15,16]. These biomarkers estimate age from DNA methylation levels at selected cytosine-phosphate-guanine sites (CpGs). Any discrepancy between chronological age and age predicted from DNA methylation data is referred to as epigenetic age acceleration . Positive age acceleration refers to an older DNA methylation age than chronological age. Negative age acceleration refers to a younger DNA methylation age than chronological age . In adults, positive age acceleration has been associated with adverse outcomes, including higher risk of cancer, cardiovascular disease and all-cause mortality [14, 17, 18]. Two specific cord blood-based clocks have been developed to assess ‘gestational age acceleration’ at birth [12, 13]. Contrary to adults, there is no clear pattern in the directionality of associations of beneficial or detrimental prenatal exposures and gestational age acceleration . For child ages, two clocks have been developed [20, 21]. However, in blood DNA methylation samples, these clocks were outperformed by the ‘skin and blood clock,’ which is developed among participants with a broad age range [15, 22].
Previous cross-sectional studies in adolescents and adults did not report consistently on associations of epigenetic age acceleration with blood pressure, with some studies reporting positive, and others no associations [23,24,25,26,27,28]. For aortic and carotid intima-media thickness, positive associations with epigenetic age acceleration were reported in two cross-sectional studies, one in preterm born neonates and one in adults, respectively [24, 29]. Associations with arterial distensibility have not been studied. It is not known whether age acceleration during the earliest phase of life is associated with blood pressure and carotid intima-media thickness and carotid distensibility at school-age.
We hypothesized that positive age acceleration at birth and at school-age would be associated with higher blood pressure and carotid intima-media thickness, and lower carotid distensibility at age 10 years. In a population-based study among 1115 children, we examined associations of gestational age acceleration, estimated from cord blood DNA methylation data, with blood pressure, carotid intima-media thickness and carotid distensibility at age 10 years. We also assessed associations of age acceleration at the ages of 6 and 10 years, estimated from blood DNA methylation data, with the same outcomes, among 470 and 449 children, respectively.
In the current study, which was embedded in the Generation R Study, a prospective cohort study from pregnancy onwards, 1115 children with information on cord blood DNA methylation and at least one of the relevant outcomes at age 10 years were included for the analyses at birth . We included 470 children at age 6 years and 449 children at age 10 years with information on DNA methylation measured in peripheral blood and at least one relevant outcome at age 10 years. Of these, 458 and 435 children, respectively, were also included in the analyses at birth. Additional file 1: Figure S1 shows a flowchart of the study population. Table 1 and Additional file 2: Table S1 shows subject characteristics before and after imputation for covariates, respectively. DNA methylation gestational age was estimated from cord blood DNA methylation data, using Bohlin’s epigenetic clock . This clock estimates gestational age based on DNA methylation levels of 96 CpGs from HumanMethylation450 BeadChip selected trough Lasso-regression . At birth, children had older clinically determined gestational age (median 40.2 weeks (95% range 37.0, 42.4)) than DNA methylation gestational age (median 39.4 weeks (95% range 37.2, 40.8)). Consequently, raw gestational age acceleration (DNA methylation age minus clinical age) had a negative median value. Pearson’s correlation coefficient between clinical gestational age and DNA methylation gestational age was 0.77 (Additional file 1: Figure S2). Among 297 children born to mothers with optimal clinical pregnancy dating, based on a regular menstrual cycle and a known first date of last menstrual period, this correlation increased slightly (r = 0.80). A non-response analysis suggested that as compared to children included in the analyses at birth, the 281 non-included children were more often boys and had younger, lower educated mothers, who less frequently used folic acid supplements and smoked more often during pregnancy (Additional file 2: Table S2).
At the age of 6 years, chronological age (median 6.0; 95% range: 5.7, 7.0 years) was older than DNA methylation age (median 5.6; 95% range: 4.3, 7.9 years) estimated by the skin and blood clock (Table 1) . This clock is developed on fibroblasts, keratinocytes, buccal cells, endothelial cells, blood, and saliva from newborns, children and adults and predicts age based on the methylation levels at 391 CpGs from the HumanMethylation450 BeadChip selected through elastic net regression . Negative age acceleration was also observed at the age of 10 years (median chronological age 9.8; 95% range 9.2, 10.3 years; median DNA methylation age 8.5; 95% range: 6.7, 11.5; Table 1). In childhood, Pearson’s correlations between chronological age and DNA methylation age were low: r = 0.35 (6 years) and r = 0.20 (10 years) (Additional file 1: Figure S2).
Gestational age acceleration and cardiovascular outcomes at school-age
We used linear regression models to examine the associations of gestational age acceleration with blood pressure and carotid intima-media and carotid distensibility measured at age 10 years. To examine whether the timing of outcome measurement influenced our results, we also examined associations of gestational age acceleration with blood pressure measured at age 6 years, in the subgroup of included children with these data available. In the full population, neither raw nor residual (residuals from regressing DNA methylation age on clinical age) gestational age acceleration by Bohlin’s method were associated with blood pressure, carotid intima-media thickness or carotid distensibility in the main model, which was adjusted for child sex, batch, child age at outcome measurement, cell types and maternal confounders (age, education, pre-pregnancy body mass index and both folic acid supplementation and smoking during pregnancy) (all P values ≥ 0.05; Tables 2 and 3). This was not different in the subgroup of 295 children born to mothers with optimal pregnancy dating. The results from basic models (not adjusted for maternal confounders) and reduced main models (not corrected for cell type proportions) are shown in Additional file 2: Table S3. Linear mixed effect models for repeated outcome measures showed that gestational age acceleration was not associated with systolic or diastolic blood pressure at school-age (Additional file 2: Table S4).
A sensitivity analysis showed that results obtained for the right and left common carotid artery separately were similar, indicating that there were no differences based on the side of carotid artery measurement (Additional file 2: Table S5). We also estimated gestational age acceleration using another cord blood-based epigenetic clock, developed by Knight. This clock estimates gestational age based on methylation levels of 148 CpGs from the HumanMethylation27 BeadChip selected through elastic net regression . Pearson’s correlation coefficients between clinical gestational age and DNA methylation gestational age by Knight in the full population and subgroup with optimal pregnancy dating are shown in Additional file 1: Figure S2. We also observed no significant associations of gestational age acceleration by Knight with blood pressure, carotid intima-media thickness or carotid distensibility at age 10 years (Additional file 2: Tables S6 and S7).
Childhood age acceleration and cardiovascular outcomes at school-age
Age acceleration estimated by the ‘skin and blood’ clock at ages 6 and 10 years was not associated with blood pressure at either age. Similarly, it was not associated with carotid intima-media thickness and carotid distensibility at age 10 years (main model, Tables 4, 5) . Conditional regression analyses that accounted for the correlations between repeated age acceleration measurements showed that age acceleration at age 6 years was not associated with systolic or diastolic blood pressure, carotid intima-media thickness or carotid distensibility at age 10 years. Also, independent of epigenetic age acceleration at age 6 years, age acceleration at age 10 years was not associated with these cardiovascular outcomes (Additional file 2: Table S8). In neither the analyses at birth nor the analyses at school-age did the results change substantially upon additional adjustment for maternal history of hypertension and family history of cardiovascular disease (data not shown).
In this population-based prospective cohort study, we observed that age acceleration at birth or in childhood is not associated with blood pressure, carotid intima-media thickness and carotid distensibility at school-age. These correlated cardiovascular risk factors partly originate in the earliest phase of life [1,2,3,4,5,6, 9,10,11]. Altered epigenetic aging may be an underlying mechanism for associations of early-life exposures with cardiovascular health at school-age [12,13,14,15,16, 19]. We hypothesized that positive epigenetic age acceleration at birth and at school-age would be associated with higher blood pressure and carotid intima-media thickness, and lower carotid distensibility at the age of 10 years.
We observed no associations of gestational age acceleration with blood pressure and carotid intima-media thickness and carotid distensibility in childhood. These null findings were not explained by inaccurate clinical pregnancy dating, as we obtained similar results among children born to mothers with a regular menstrual cycle and pregnancy dating based on last menstrual period. Previous cross-sectional studies reported inconsistently on associations of epigenetic age acceleration with blood pressure [23,24,25,26,27,28, 31, 32]. Both an American and a Scottish study among 969 and 5100 adults, respectively, reported that positive age acceleration as estimated by the clock of Hannum, but not that of Horvath, was associated with hypertension [26, 28]. A study among 1100 mostly hypertensive African-Americans reported that positive age acceleration was associated with blood pressure, using Horvath’s and Hannum’s clocks . In addition, this study, as well as two studies among 4178 women from the Women’s Health Initiative, only reported associations with systolic but not diastolic blood pressure when using the PhenoAge and GrimAge clocks [27, 31, 32]. Three further studies in adults and one in adolescents reported no associations of age acceleration with blood pressure [23,24,25, 33].
For intima-media thickness, an Australian study among 169 neonates reported that gestational age acceleration estimated in saliva using the clock of Horvath was positively associated with aortic intima-media thickness measured within 24 h postpartum in preterm, but not term, born neonates . Another cross-sectional study, among 2543 middle-aged African-Americans, reported positive associations of age acceleration estimated by both the clocks of Hannum and Horvath with carotid intima-media thickness . Associations of age acceleration with arterial distensibility have not been studied, to the best of our knowledge. Overall, previous cross-sectional studies show some evidence for positive associations of age acceleration with blood pressure and carotid intima-media thickness. Results from our population-based cohort study, however, suggest no prospective or cross-sectional associations of age acceleration with these outcomes, nor with carotid distensibility, at school-age.
Our null findings may be interpreted in two ways. First, it could be that DNA methylation is simply not a mechanism underlying the associations of early-life exposures with childhood health. Some exposures, such as maternal smoking and folate concentrations during pregnancy, have been related to both differential DNA methylation and childhood cardiovascular health [5, 34,35,36]. This does not necessarily imply that DNA methylation is on the biological pathway between these exposures and childhood cardiovascular health. DNA methylation changes could be epiphenomena rather than consequences of early-life exposures . Second, DNA methylation may in fact be on the causal pathway from early-life exposures to childhood cardiovascular health, but at specific CpG sites that are not captured by the three clocks we used in this study. Only one CpG (cg07141002) and three CpGs (cg01243072, cg06690548, cg18933331) included in the clock of Knight and the skin and blood clock, respectively, have been individually associated with blood pressure in adults [38, 39]. The only CpG (cg05575921) that has been associated with carotid intima-media thickness in adults is not included in any of the three clocks, whereas arterial distensibility has not been studied in relation to DNA methylation [38, 40]. In adults, different epigenetic clocks seem to capture different age-related biological processes . In addition, previous studies in newborns and adults observed different associations of age acceleration with cardiometabolic outcomes, depending on which clock was used [26,27,28, 41, 42]. The clocks of Bohlin and Knight and the skin and blood clock may thus capture aspects of age acceleration that are not a reflection of cardiovascular adaptations in utero or in childhood.
The null findings may also be explained by some limitations of our study. First, although we had a relatively large study population, we may have been underpowered to find associations of smaller magnitude, especially in the subgroup with optimal pregnancy dating. Although the correlation between clinical gestational age and DNA methylation gestational age increased in this subgroup, the decrease in power may not have been outweighed by increased precision in the estimation of clinical gestational age. Second, our study population was relatively healthy. Limited variation in either age acceleration or in the outcomes may have prevented detection of meaningful associations. It needs confirmation whether associations are observed among populations showing more extreme variation in exposure or outcomes. Despite these limitations, our study also had strengths. These include its implementation in a large observational prospective birth cohort, the availability of repeated DNA methylation measurements and detailed cardiovascular outcomes, and the replication of the null findings of the main analysis at birth using another epigenetic clock.
Epigenetic age acceleration at birth or at school-age was not associated with blood pressure at ages 6 and 10 years or with carotid intima-media thickness and carotid distensibility at age 10 years. Our study does not support altered epigenetic aging as an underlying mechanism for associations of exposures during the earliest phase of life with childhood cardiovascular health. Further research is needed to confirm this.
This study was embedded in the Generation R Study, a population-based prospective cohort study from fetal life onwards in Rotterdam, the Netherlands . The Medical Ethical Committee of Erasmus MC, University Medical Center Rotterdam, approved the study (MEC 198.782/2001/31). Pregnant women with an expected delivery date between April 2002 and January 2006 living in Rotterdam were eligible to participate. We obtained written informed consent from all participants. In a relatively homogeneous European-ancestry subgroup taken from the total of 9901 live-born participating children, we measured genome-wide DNA methylation in cord blood taken at birth (n = 1396) and in peripheral blood at 6 years (n = 493) and 10 years (n = 464), when the children visited the research center for follow-up visits. Additionally, at age 10 years, we measured systolic and diastolic blood pressure and both common carotid artery intima-media thickness and carotid distensibility. In the current study, we included children with DNA methylation measured at any time point, and at least one of these cardiovascular outcomes measured (Additional file 1: Figure S1). There were no twins among our population for analysis. For each sibling pair, we included only one child by selecting based on completeness of outcome data, and if equal, randomly.
DNA methylation data
We used the salting-out method to extract DNA from cord blood and peripheral blood samples in childhood. Five-hundred nanograms of DNA were bisulfite converted using the EZ-96 DNA Methylation kit (Shallow) (Zymo Research Corporation, Irvine, USA). DNA methylation samples were plated onto 96-well plates in no specific order. Samples were processed with the Illumina Infinium HumanMethylation450 BeadChip (Illumina Inc., San Diego, USA). Quality control and normalization were performed using the CPACOR workflow . Probes with a detection P value ≥ 1E-16 were set to missing. Intensity values were quantile normalized. We removed arrays with technical problems, a call rate ≤ 95%, or a mismatch between the expected sex of participant and sex determined by chromosome X and Y probe intensities. Probes on the sex chromosomes were removed before the analyses. We used untransformed beta-values as measures of DNA methylation. The final DNA methylation dataset contained information on 458,563 CpGs.
Gestational age estimation
Mothers of participating children visited the research center in the first trimester of the index pregnancy for fetal ultrasound. During this visit, we established the ‘clinical’ estimate of gestational age . For those mothers with a known and reliable first day of the last menstrual period and a regular (28 ± 4 days) menstrual cycle, gestational age was based on last menstrual period . If mothers did not know their last menstrual period exactly, or had an irregular menstrual cycle, we established gestational age by ultrasound examination. Gestational age at birth was assessed from midwife or obstetric records. We calculated DNA methylation gestational age at birth (in weeks) using the epigenetic clock of Bohlin as primary method, as it was developed among newborns comparable to those in our study population . This clock estimates gestational age from DNA methylation levels at 96 CpGs from HumanMethylation450 BeadChip selected trough Lasso-regression . Bohlin’s DNA methylation gestational age was calculated using the GAprediction package version 1.16.0 in R 3.6.1 [13, 45]. We excluded children with missing values for one or more of the 96 required CpGs (total study population: n = 11; subgroup with optimal pregnancy dating: n = 2). As a sensitivity analysis, to examine the robustness of our findings, we estimated DNA methylation gestational age by another cord blood based epigenetic clock, developed by Knight . This clock estimates gestational age from DNA methylation levels at 148 CpGs that are available on the HumanMethylation450 BeadChip, selected through elastic net regression [12, 46]. Knight’s DNA methylation gestational age was calculated using the methylclock package 0.5.0 in R 3.6.1 (R Core Team, Vienna, Austria) . All children had information on the 148 required CpGs.
Childhood age estimation
Of the 470 and 449 children included in the analyses at age 6 and age 10 years, respectively, twelve (6 years) and fourteen (10 years) children were not included in the analyses at birth. Childhood DNA methylation age (in years) was estimated by the ‘skin and blood clock’, using the methylclock package 0.5.0 in R 3.6.1 (R Core Team, Vienna, Austria) [15, 47]. A previous study reported that this well-established clock outperformed two child-specific clocks when estimating DNA methylation age in blood [20,21,22]. The skin and blood clock estimates age from DNA methylation levels at 391 CpGs measured in fibroblasts, keratinocytes, buccal cells, endothelial cells, blood, and saliva from participants aged 0–94 years . These CpGs were automatically selected trough elastic net regression that were available on the HumanMethylation450 BeadChip . We obtained chronological (‘clinical’) age (in years) at the time of the child’s follow-up visit to the research center.
Age acceleration at birth and in childhood
Based on the epigenetic clocks of Bohlin and Knight at birth and the skin and blood clock at school-age, we calculated both raw and residual age acceleration, in line with previous research [46,47,48]. Raw age acceleration was obtained by subtracting clinical age from DNA methylation age. It does not take into account the potential confounding effect of clinical age on DNA methylation age, as variance is shared between both measures . We obtained residual age acceleration from the residuals from a linear regression of DNA methylation age on clinical age. By definition, residual acceleration is uncorrelated with clinical age . Positive age acceleration indicates older DNA methylation age than clinical age. Negative age acceleration indicates younger DNA methylation age than clinical age.
Cardiovascular outcomes at school-age
When children visited the research facility at age 6 and 10 years, we assessed blood pressure at the right brachial artery four times with one minute intervals using the validated automatic sphygmomanometer Datascope Accutorr Plus (Paramus, New Jersey, United States) . For children with complete data on blood pressure measurements available, the mean systolic and diastolic blood pressure was calculated from the last three measurements. Further, during the visit at age 10 years, we measured carotid intima-media thickness and carotid distensibility three times at both common carotid arteries (n = 5,746) using the Logiq E9 (GE Medical Systems, Wauwatosa, WI, USA) device. Children were in the supine position, with the head tilted slightly away from the transducer. The common carotid artery was identified in a longitudinal plane, ~ 10 mm proximal from the carotid bifurcation. We obtained six recordings that each ideally included multiple heart cycles. The analyses were performed offline and semi-automatically, using the application Carotid Studio (Cardiovascular Suite (Quipu srl, Pisa, Italy)). For each recording, at all R-waves of the simultaneous ECG, carotid intima-media thickness was computed at the far wall, as the average distance between lumen-intima and media-adventitia borders. The average carotid intima-media thickness of all frames of the acquired image sequence was computed. Carotid distensibility was defined as the relative change in lumen area during systole for a given peripheral brachial artery pressure change. The lumen diameter of the carotid artery was computed as the average distance between the far and near media‐adventitia interfaces, for each frame of the acquired image sequence. Distension was calculated as the difference between maximal (diastolic) and minimal (systolic) lumen diameter. Per recording, the average distension and diameter values were used to compute the average carotid distensibility. During these offline analyses, we excluded 516 and 704 children, without any valid carotid intima-media thickness or carotid distensibility measurement, respectively. Reasons were no appropriate recording, insufficient quality of the recording, recording of the heart only, or no blood pressure measurement available to calculate carotid distensibility. Further data processing for the remaining 5,230 and 5,042 children with carotid intima-media thickness and carotid distensibility data, respectively, was performed using R (R Core Team, Vienna, Austria). We excluded 9 children with unreliable low or high carotid distensibility values. We calculated the overall mean carotid intima-media thickness (in mm) and carotid distensibility (in kPa−1*10–3), both based on up to six measurements. In a reproducibility study performed among 47 subjects, the interobserver and intraobserver intraclass correlation coefficient were > 0.85 for carotid distensibility and > 0.94 for carotid intima-media thickness.
We selected potential covariates based on the previous literature. From maternal questionnaires sent out at intake, we obtained information on the following maternal covariates: age (in years), educational level (primary or secondary education, versus college or higher), pre-pregnancy body mass index (kg/m2) and both folic acid supplementation (no supplementation, started before ten weeks gestational age and started preconception) and smoking (no smoking or quit before second trimester, versus sustained smoking) during pregnancy. We also obtained information on maternal hypertension and family history of cardiovascular disease, which we defined as hypertension, myocardial infarction before the age of 65, or cerebrovascular attack in a first-degree relative. From midwife and medical records, we obtained information on child sex. We adjusted for batch effects to correct for systematic technical variation in the DNA methylation data by including sample plate number. For the analyses of gestational age acceleration, we used the “Salas” reference set to estimate cell type proportion in the ‘FlowSorted.CordBlood.Combined.450 K’ Bioconductor package. It includes CD8+ T cells, CD4+ T cells, Natural Killer cells, B cells, Monocytes, Granulocytes and Nucleated red blood cells . For the analyses in childhood, we used the Houseman method based on the Reinius reference set for cell type estimation (CD8+ T-cells, CD4+ T-cells, Natural Killer cells, B cells, Monocytes, Granulocytes) in the ‘FlowSorted.Blood.450 K’ package [51, 52].
First, we performed a non-response analysis using Student’s t-tests, Mann–Whitney tests and Chi-square tests. Among those children with cord blood DNA methylation measured, we compared characteristics of children included in the analyses at birth, to those of non-included children because they had no information on any cardiovascular outcome of interest, or had a sibling included in the analysis. We used multiple imputations for covariates with missing values, using the Markov Chain Monte Carlo method. We created five datasets and pooled analysis estimates from all datasets . Second, we examined the associations of raw and residual gestational age acceleration estimated by both the method of Bohlin, with all cardiovascular outcomes, using multivariable linear regression models. To compare effect estimates, the outcomes were analyzed in SDS, after natural log transformation of carotid distensibility, which had a skewed distribution. Basic models were adjusted for sample plate number (batch), estimated cell types, child sex and age at outcome measurement. Main models were additionally adjusted for maternal age, education, pre-pregnancy body mass index, folic acid supplementation and smoking. In line with previous research, the main model was also run without cell type adjustment (reduced main model), to examine the specific influence of variation in cell type proportions . We performed three sensitivity analyses. First, in an attempt to minimize the chance of measurement error in the calculation of gestational age acceleration, we restricted the analyses at birth to children of mothers with optimal pregnancy dating, based on last menstrual period. We reasoned that such error may result from inaccuracy in the prediction of clinical gestational age by fetal ultrasound, which does not take into account variation in early fetal growth . Second, as recommended in previous literature, we assessed the associations of gestational age acceleration with carotid intima-media thickness and carotid distensibility measured at the right versus left common carotid artery . Third, both in the full population and in the subgroup with optimal pregnancy dating, we reran the main model after calculating gestational age acceleration based on the clock of Knight. In childhood, we assessed the associations of raw and residual age acceleration at age 6 and 10 years with all outcomes (main models only). Further, we explored whether the timing of outcome measurement is critical for associations of age acceleration with blood pressure, on which we had information at age 6 years. To do this, we analyzed associations of age acceleration at birth and at 6 years with systolic and diastolic blood pressure. Finally, as a secondary analysis, in the full population and for each time point and all outcomes, we additionally adjusted the main model for both maternal hypertension and family history of cardiovascular disease. We considered P values < (0.05/4 assessed outcomes), so < 0.00125 statistically significant. We did not stratify the analyses on sex because the interaction terms between age acceleration measures and sex were not significant. We examined potential nonlinear associations between gestational age acceleration and all outcomes by adding a quadratic term for raw or residual gestational age acceleration to the main models. None of the quadratic terms were significant (all P > 0.05). All statistical analyses were performed using the Statistical Package of Social Sciences version 25.0 for Windows (SPSS IBM, Chicago, Illinois, United States).
Making use of the repeated measurements in our data, we performed linear mixed effects models using the R package nlme to examine the associations of gestational age acceleration with repeatedly measured systolic and diastolic blood pressure at ages 6 and 10 years. To decrease the number of covariates, we first obtained standardized residuals for both raw and residual gestational age acceleration from the regression of gestational age acceleration on batch and all cell types and used these residuals as the exposure. All models included a random intercept. We additionally included an interaction term between age acceleration and child age to allow the exposure-outcome association to change across childhood. We also examined the independent associations of age acceleration at ages 6 and 10 years with all cardiovascular outcomes at age 10 years using conditional multivariable linear regression to account for the correlations between epigenetic age measurements using SPSS . For these analyses, only the 329 children with information on age acceleration at both age 6 and age 10 years were included. First, for both ages, we obtained standardized residuals for both raw and residual age acceleration from the regression of age acceleration on batch and all cell types. Second, the obtained residuals at 10 years were made independent of those at 6 years by regressing the 10-year residuals on those at 6 years. Third, both measures were simultaneously included in the models.
Availability of data and materials
The datasets generated and/or analyzed during the current study are not publicly available due to privacy restrictions, but are available from the corresponding author on reasonable request, subject to the Generation R Study data access procedures.
Epigenome-wide association study
Standard deviation score
Raitakari OT, Juonala M, Kähönen M, Taittonen L, Laitinen T, Mäki-Torkko N, et al. Cardiovascular risk factors in childhood and carotid artery intima-media thickness in adulthood: the Cardiovascular Risk in Young Finns Study. JAMA. 2003;290(17):2277–83.
Allen NB, Krefman AE, Labarthe D, Greenland P, Juonala M, Kähönen M, et al. Cardiovascular health trajectories from childhood through middle age and their association with subclinical atherosclerosis. JAMA Cardiol. 2020;5(5):557–66.
Li S, Chen W, Srinivasan SR, Bond MG, Tang R, Urbina EM, et al. Childhood cardiovascular risk factors and carotid vascular changes in adulthood: the Bogalusa Heart Study. JAMA. 2003;290(17):2271–6.
Jonna J, Mervi O, Costan GM, Vera M, Niina S, Eero J, et al. Childhood physical, environmental, and genetic predictors of adult hypertension. Circulation. 2012;126(4):402–9.
Lawlor DA, Najman JM, Sterne J, Williams GM, Ebrahim S, Davey SG. Associations of parental, birth, and early life characteristics with systolic blood pressure at 5 years of age: findings from the Mater-University study of pregnancy and its outcomes. Circulation. 2004;110(16):2417–23.
Geerts CC, Bots ML, van der Ent CK, Grobbee DE, Uiterwaal CS. Parental smoking and vascular damage in their 5-year-old children. Pediatrics. 2012;129(1):45–54.
Yuan C, Wang J, Ying M. Predictive value of carotid distensibility coefficient for cardiovascular diseases and all-cause mortality: a meta-analysis. PLoS ONE. 2016;11(4):e0152799.
Lorenz MW, Markus HS, Bots ML, Rosvall M, Sitzer M. Prediction of clinical cardiovascular events with carotid intima-media thickness. Circulation. 2007;115(4):459–67.
Alexander J, Marietta C, Emanuela F, Aroon DH, Nicholas F, Stefano M, et al. Adipose and height growth through childhood and blood pressure status in a large prospective cohort study. Hypertension. 2012;59(5):919–25.
Geerts CC, Evelein AMV, Bots ML, van der Ent CK, Grobbee DE, Uiterwaal CSPM. Body fat distribution and early arterial changes in healthy 5-year-old children. Ann Med. 2012;44(4):350–9.
Hanna M, Katja P, Harri N, Tapani R, Jorma SAV, Antti J, et al. Cardiometabolic determinants of carotid and aortic distensibility from childhood to early adulthood. Hypertension. 2017;70(2):452–60.
Knight AK, Craig JM, Theda C, Bækvad-Hansen M, Bybjerg-Grauholm J, Hansen CS, et al. An epigenetic clock for gestational age at birth based on blood methylation data. Genome Biol. 2016;17(1):206.
Bohlin J, Håberg SE, Magnus P, Reese SE, Gjessing HK, Magnus MC, et al. Prediction of gestational age based on genome-wide differentially methylated regions. Genome Biol. 2016;17(1):207.
Horvath S. DNA methylation age of human tissues and cell types. Genome Biol. 2013;14(10):3156.
Horvath S, Oshima J, Martin GM, Lu AT, Quach A, Cohen H, et al. Epigenetic clock for skin and blood cells applied to Hutchinson Gilford Progeria Syndrome and ex vivo studies. Aging (Albany NY). 2018;10(7):1758–75.
Liu Z, Leung D, Thrush K, Zhao W, Ratliff S, Tanaka T, et al. Underlying features of epigenetic aging clocks in vivo and in vitro. Aging Cell. 2020;19(10):e13229.
Chen BH, Marioni RE, Colicino E, Peters MJ, Ward-Caviness CK, Tsai PC, et al. DNA methylation-based measures of biological age: meta-analysis predicting time to death. Aging. 2016;8(9):1844–65.
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):1–7.
Dieckmann L, Lahti-Pulkkinen M, Kvist T, Lahti J, DeWitt PE, Cruceanu C, et al. Characteristics of epigenetic aging across gestational and perinatal tissues. Clin Epigenetics. 2021;13(1):97.
McEwen LM, O’Donnell KJ, McGill MG, Edgar RD, Jones MJ, MacIsaac JL, et al. The PedBE clock accurately estimates DNA methylation age in pediatric buccal cells. Proc Natl Acad Sci. 2020;117(38):23329–35.
Wu X, Chen W, Lin F, Huang Q, Zhong J, Gao H, et al. DNA methylation profile is a quantitative measure of biological aging in children. Aging (Albany NY). 2019;11(22):10031–51.
Kling T, Wenger A, Carén H. DNA methylation-based age estimation in pediatric healthy tissues and brain tumors. Aging. 2020;12(21):21037–56.
Quach A, Levine ME, Tanaka T, Lu AT, Chen BH, Ferrucci L, et al. Epigenetic clock analysis of diet, exercise, education, and lifestyle factors. Aging. 2017;9(2):419–46.
Roetker NS, Pankow JS, Bressler J, Morrison AC, Boerwinkle E. Prospective study of epigenetic age acceleration and incidence of cardiovascular disease outcomes in the ARIC study (Atherosclerosis risk in communities). Circ Genomic Precis Med. 2018;11(3):e001937.
Huang R, Lillycrop KA, Beilin LJ, Godfrey KM, Anderson D, Mori AT, 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(7):3012–24.
McCartney DL, Stevenson A, Walker RM, Gibson J, Morris SW, Campbell A, et al. Investigating the relationship between DNA methylation age acceleration and risk factors for Alzheimer’s disease. Alzheimers Dement (Amst). 2018;10:429–427.
Ammous F, Zhao W, Ratliff SM, Mosley TH, Bielak LF, Zhou X, et al. Epigenetic age acceleration is associated with cardiometabolic risk factors and clinical cardiovascular disease risk scores in African Americans. Clin Epigenetics. 2021;13(1):55.
Horvath S, Gurven M, Levine ME, Trumble BC, Kaplan H, Allayee H, et al. An epigenetic clock analysis of race/ethnicity, sex, and coronary heart disease. Genome Biol. 2016;17(1):171.
Phang M, Ross J, Raythatha JH, Dissanayake HU, McMullan RL, Kong Y, et al. Epigenetic aging in newborns: role of maternal diet. Am J Clin Nutr. 2020;111(3):555–61.
Kooijman MN, Kruithof CJ, Van Duijn CM, Duijts L, Franco OH, Van Ijzendoorn MH, et al. The Generation R Study: design and cohort update 2017. Eur J Epidemiol. 2016;31(12):1243–64.
Levine ME, Lu AT, Quach A, Chen BH, Assimes TL, Bandinelli S, et al. An epigenetic biomarker of aging for lifespan and healthspan. Aging. 2018;10(4):573–91.
Lu AT, Quach A, Wilson JG, Reiner AP, Aviv A, Raj K, et al. DNA methylation GrimAge strongly predicts lifespan and healthspan. Aging. 2019;11(2):303–27.
Smith JA, Raisky J, Ratliff SM, Liu J, Kardia SLR, Turner ST, et al. Intrinsic and extrinsic epigenetic age acceleration are associated with hypertensive target organ damage in older African Americans. BMC Med Genomics. 2019;12(1):141.
Joubert BR, Felix JF, Yousefi P, Bakulski KM, Just AC, Breton C, et al. DNA methylation in newborns and maternal smoking in pregnancy: genome-wide consortium meta-analysis. Am J Hum Genet. 2016;98(4):680–96.
Joubert BR, Den Dekker HT, Felix JF, Bohlin J, Ligthart S, Beckett E, et al. Maternal plasma folate impacts differential DNA methylation in an epigenome-wide meta-analysis of newborns. Nat Commun. 2016;10(7):10577.
Monasso GS, Santos S, Geurtsen ML, Heil SG, Felix JF, Jaddoe VWV. Associations of early pregnancy and neonatal circulating folate, vitamin B-12, and homocysteine concentrations with cardiometabolic risk factors in children at 10 y of age. J Nutr. 2021;151(6):1628–36.
Felix JF, Cecil C. Population DNA methylation studies in the Developmental Origins of Health and Disease (DOHaD) framework. J Dev Orig Health Dis. 2019;10(3):306–13.
Battram T, Yousefi P, Crawford G, Prince C, Babei MS, Sharp G, et al. The EWAS Catalog: a database of epigenome-wide association studies. OSF Preprints; 2021. Available from: https://osf.io/837wn.
Richard MA, Huan T, Ligthart S, Gondalia R, Jhun MA, Brody JA, et al. DNA methylation analysis identifies loci for blood pressure regulation. Am J Hum Genet. 2017;101(6):888–902.
Portilla-Fernández E, Hwang S-J, Wilson R, Maddock J, Hill WD, Teumer A, et al. Meta-analysis of epigenome-wide association studies of carotid intima-media thickness (Online ahead of print). Eur J Epidemiol. 2021. https://doi.org/10.1007/s10654-021-00759-z.
Simpkin AJ, Howe LD, Tilling K, Gaunt TR, Lyttleton O, McArdle W, et al. The epigenetic clock and physical development during childhood and adolescence: longitudinal analysis from a UK birth cohort. Int J Epidemiol. 2017;2(46):549–58.
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 Epigenetics. 2019;11(1):194.
Lehne B, Drong AW, Loh M, Zhang W, Scott WR, Tan ST, et al. A coherent approach for analysis of the Illumina HumanMethylation450 BeadChip improves data quality and performance in epigenome-wide association studies. Genome Biol. 2015;16(1):37.
Gaillard R, Steegers EA, de Jongste JC, Hofman A, Jaddoe VWV. Tracking of fetal growth characteristics during different trimesters and the risks of adverse birth outcomes. Int J Epidemiol. 2014;43(4):1140–53.
Bohlin J. GAprediction: prediction of gestational age with Illumina HumanMethylation450 data. R package version 1.16.0 ed2020.
Girchenko P, Lahti J, Czamara D, Knight AK, Jones MJ, Suarez A, et al. Associations between maternal risk factors of adverse pregnancy and birth outcomes and the offspring epigenetic clock of gestational age at birth. Clin Epigenetics. 2017;8(9):49.
Pelegí-Sisó D, de Prado P, Ronkainen J, Bustamante M, González JR. Methylclock: a bioconductor package to estimate DNA methylation age. Bioinformatics. 2020;37:1759–60.
Chen L, Wagner CL, Dong Y, Wang X, Shary JR, Huang Y, et al. Effects of maternal vitamin D3 supplementation on offspring epigenetic clock of gestational age at birth: a post-hoc analysis of a randomized controlled trial. Epigenetics. 2020;15(8):830–40.
Wong SN, Tz Sung RY, Leung LC. Validation of three oscillometric blood pressure devices against auscultatory mercury sphygmomanometer in children. Blood Press Monit. 2006;11:281–91.
Gervin K, Salas LA, Bakulski KM, van Zelm MC, Koestler DC, Wiencke JK, et al. Systematic evaluation and validation of reference and library selection methods for deconvolution of cord blood DNA methylation data. Clin Epigenetics. 2019;11(1):125.
Houseman EA, Molitor J, Marsit CJ. Reference-free cell mixture adjustments in analysis of DNA methylation data. Bioinformatics. 2014;30(10):1431–9.
Reinius LE, Acevedo N, Joerink M, Pershagen G, Dahlén S-E, Greco D, et al. Differential DNA methylation in purified human blood cells: implications for cell lineage and studies on disease susceptibility. PLoS ONE. 2012;7(7):e41361.
Sterne JA, White IR, Carlin JB, Spratt M, Royston P, Kenward MG, et al. Multiple imputation for missing data in epidemiological and clinical research: potential and pitfalls. BMJ. 2009;338:b2393.
Khouja JN, Simpkin AJ, O’Keeffe LM, Wade KH, Houtepen LC, Relton CL, et al. Epigenetic gestational age acceleration: a prospective cohort study investigating associations with familial, sociodemographic and birth characteristics. Clin Epigenetics. 2018;10(86):1–9.
Dalla Pozza R, Ehringer-Schetitska D, Fritsch P, Jokinen E, Petropoulos A, Oberhoffer R. Intima media thickness measurement in children: a statement from the Association for European Paediatric Cardiology (AEPC) Working Group on Cardiovascular Prevention endorsed by the Association for European Paediatric Cardiology. Atherosclerosis. 2015;238(2):380–7.
Gishti O, Gaillard R, Manniesing R, Abrahamse-Berkeveld M, van der Beek EM, Heppe DH, et al. Fetal and infant growth patterns associated with total and abdominal fat distribution in school-age children. J Clin Endocrinol Metab. 2014;99(7):2557–66.
The Generation R Study is conducted by the Erasmus Medical Center in close collaboration with the School of Law and Faculty of Social Sciences of the Erasmus University Rotterdam, the Municipal Health Service Rotterdam area, Rotterdam, the Rotterdam Homecare Foundation, Rotterdam and the Stichting Trombosedienst & Artsenlaboratorium Rijnmond (STAR-MDC), Rotterdam. We gratefully acknowledge the contribution of children and parents, general practitioners, hospitals, midwives and pharmacies in Rotterdam. The generation and management of the Illumina 450K methylation array data (EWAS data) for the Generation R Study was executed by the Human Genotyping Facility of the Genetic Laboratory of the Department of Internal Medicine, Erasmus MC, and the Netherlands. We thank Mr. Michael Verbiest, Ms. Mila Jhamai, Ms. Sarah Higgins, Mr. Marijn Verkerk and Dr. Lisette Stolk for their help in creating the EWAS database. We thank Dr. Alexander Teumer for his work on the quality control and normalization scripts.
The general design of the Generation R Study is made possible by financial support from the Erasmus MC, University Medical Centre Rotterdam, Erasmus University Rotterdam, the Netherlands Organization for Health Research and Development (ZonMw), the Netherlands Organization for Scientific Research (NWO), the Ministry of Health, Welfare and Sport and the Ministry of Youth and Families. The EWAS data were funded by a grant to VWVJ from the Netherlands Genomics Initiative (NGI)/Netherlands Organization for Scientific Research (NWO) Netherlands Consortium for Healthy Aging (NCHA; project number 050-060-810), by funds from the Genetic Laboratory of the Department of Internal Medicine, Erasmus MC, and by a grant from the National Institute of Child and Human Development (R01HD068437). VWVJ received funding from the European Research Council (ERC-2014-CoG-648916). The project was supported by funding from the European Union's Horizon 2020 research and innovation program under grant agreements No 733206 (LifeCycle) and 874739 (LongITools), and from the European Joint Programming Initiative “A Healthy Diet for a Healthy Life” (JPI HDHL, NutriPROGRAM project, ZonMw the Netherlands no.529051022 and PREcisE project ZonMw the Netherlands no.529051023).
Ethics approval and consent to participate
The Medical Ethical Committee of Erasmus MC, University Medical Center Rotterdam, approved the study protocol of the Generation R study (MEC 198.782/2001/31). Written informed consent was obtained for all participants.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. Flowchart of the study population. Figure S2. Pearson’s correlation between clinical age and DNA methylation age.
. Maternal and child characteristics based on imputed data (n = 1115). Table S2. Non-response analysis. Table S3. Associations of gestational age acceleration by the epigenetic clock of Bohlin with cardiovascular outcomes in children aged ten years (basic model and reduced main model). Table S4. Associations of gestational age acceleration by the epigenetic clock of Bohlin with repeated blood pressure measurements (linear mixed effect models). Table S5. Associations of gestational age acceleration by the epigenetic clock of Bohlin with carotid intima-media thickness and carotid distensibility in children aged ten years, for the right and left common carotid artery separately (main model). Table S6. Associations of gestational age acceleration by the epigenetic clock of Knight with blood pressure in children aged six years (main model). Table S7. Associations of gestational age acceleration by the epigenetic clock of Knight with cardiovascular outcomes in children aged ten years (main model). Table S8. Associations of repeated epigenetic age acceleration measurements at age 6 and 10 years with cardiovascular outcomes in children aged 10 years (conditional regression analyses)
About this article
Cite this article
Monasso, G.S., Jaddoe, V.W.V., Küpers, L.K. et al. Epigenetic age acceleration and cardiovascular outcomes in school-age children: The Generation R Study. Clin Epigenet 13, 205 (2021). https://doi.org/10.1186/s13148-021-01193-4
- DNA methylation
- Epigenetic clock
- Gestational age
- Cohort study
- Cardiovascular disease
- Carotid intima-media thickness