Accelerated DNA methylation age plays a role in the impact of cardiovascular risk factors on the human heart
Clinical Epigenetics volume 15, Article number: 164 (2023)
DNA methylation (DNAm) age acceleration (AgeAccel) and cardiac age by 12-lead advanced electrocardiography (A-ECG) are promising biomarkers of biological and cardiac aging, respectively. We aimed to explore the relationships between DNAm age and A-ECG heart age and to understand the extent to which DNAm AgeAccel relates to cardiovascular (CV) risk factors in a British birth cohort from 1946.
We studied four DNAm ages (AgeHannum, AgeHorvath, PhenoAge, and GrimAge) and their corresponding AgeAccel. Outcomes were the results from two publicly available ECG-based cardiac age scores: the Bayesian A-ECG-based heart age score of Lindow et al. 2022 and the deep neural network (DNN) ECG-based heart age score of Ribeiro et al. 2020. DNAm AgeAccel was also studied relative to results from two logistic regression-based A-ECG disease scores, one for left ventricular (LV) systolic dysfunction (LVSD), and one for LV electrical remodeling (LVER). Generalized linear models were used to explore the extent to which any associations between biological cardiometabolic risk factors (body mass index, hypertension, diabetes, high cholesterol, previous cardiovascular disease [CVD], and any CV risk factor) and the ECG-based outcomes are mediated by DNAm AgeAccel. We derived the total effects, average causal mediation effects (ACMEs), average direct effects (ADEs), and the proportion mediated [PM] with their 95% confidence intervals [CIs]. 498 participants (all 60–64 years) were included, with the youngest ECG heart age being 27 and the oldest 90. When exploring the associations between cardiometabolic risk factors and Bayesian A-ECG cardiac age, AgeAccelPheno appears to be a partial mediator, as ACME was 0.23 years [0.01, 0.52] p = 0.028 (i.e., PM≈18%) for diabetes, 0.34 [0.03, 0.74] p = 0.024 (i.e., PM≈15%) for high cholesterol, and 0.34 [0.03, 0.74] p = 0.024 (PM≈15%) for any CV risk factor. Similarly, AgeAccelGrim mediates ≈30% of the relationship between diabetes or high cholesterol and the DNN ECG-based heart age. When exploring the link between cardiometabolic risk factors and the A-ECG-based LVSD and LVER scores, it appears that AgeAccelPheno or AgeAccelGrim mediate 10–40% of these associations.
By the age of 60, participants with accelerated DNA methylation appear to have older, weaker, and more electrically impaired hearts. We show that the harmful effects of CV risk factors on cardiac age and health, appear to be partially mediated by DNAm AgeAccelPheno and AgeAccelGrim. This highlights the need to further investigate the potential cardioprotective effects of selective DNA methyltransferases modulators.
The biological aging process is complex and involves a range of cellular and molecular changes that accumulate over time, resulting in deterioration of physiological function and increased vulnerability to chronic diseases, and eventually increased mortality . Heterogeneity across human biological aging phenotypes cannot be solely explained by chronological age . Therefore, biological aging biomarkers capable of identifying individuals at risk of functional impairment beyond chronological age have been extensively researched.
The most promising biological aging biomarkers are those based on DNA methylation (DNAm) levels at specific cytosine-phosphate-guanine (CpG) sites . The 1st generation DNAm ages were created using ElasticNet penalized regression in which CpG sites were used to predict chronological age. DNAm AgeHannum  was based on blood cells, whilst DNAm AgeHorvath  incorporated methylation data from multiple tissues. Both were highly correlated with chronological age but poorly associated with CpG sites that captured lifespan and healthspan. As a result, 2nd generation DNAm biomarkers emerged: DNAm PhenoAge based on whole-blood CpG sites which associated with a composite of mortality-related clinical and physiological measures ; and DNAm GrimAge based on plasma protein estimates, smoking pack-years, chronological age and sex as a function of the time-to-death . After regressing these DNAm ages on chronological age, the residuals can be interpreted as a measure of age acceleration (AgeAccel), providing valuable insights into the speed of the epigenetic clock. A positive residual value is adverse as it suggests faster biological aging, while a negative value is beneficial as it suggests slower biological aging. The 2nd generation AgeAccel DNAm biomarkers have been shown to be more predictive of adverse health outcomes than their 1st generation DNAm counterparts [8, 9].
Cardiac ageing is naturally associated with a progressively increasing burden of cardiovascular diseases (CVD) across the human lifespan. Given the complexity of cardiac aging, chronological age is not an optimal proxy . Thus, cardiac ages based on imaging, electrocardiography (ECG) or cardiometabolic data have been proposed [11, 12]. Of these, advanced electrocardiography (A-ECG) which incorporates results from conventional ECG, vectorcardiography and waveform complexity within Bayesian statistical frameworks, or alternatively within deep neural network (DNN) ECG-based algorithms, has the highest translational potential given the low cost and the ubiquitous availability of ECG across healthcare systems [13,14,15].
Underpinning both DNAm and cardiac aging are a set of shared risk factors such as diet , smoking, exercise , lifetime psychological stress  and ambient air pollution , amongst others. Since many of these risk factors were shown to have different epigenetic methylation signatures at specific CpG sites [19, 20], the existence of a strong link between DNAm and cardiac aging was previously postulated . Indeed, although 1st generation DNAm ages appear be weakly linked to CVDs , the 2nd generation AgeAccel DNAm were shown to be more predictive of adverse cardiovascular outcomes [8, 9]. As distinct DNAm profiles have been associated with cardiometabolic risk factors (e.g., body mass index [BMI] , diabetes , high cholesterol , hypertension  and coronary artery disease ), the idea that DNAm might be a mediator between CV risk factors and cardiac age gained traction. Thus, the role of DNA methyltransferase inhibitors (DNMT-i) was studied in pre-clinical trials with encouraging effects in ameliorating cardiac hypertrophy , fibrosis  and atherosclerosis . However, to what extent AgeAccel DNAm might mediate the effect of biological CV risk factors on ECG-based cardiac age remains to be elucidated. Using prospectively collected life-course data from the 1946 Medical Research Council (MRC) British National Survey of Health and Development (NSHD) study, we sought to answer this question.
The MRC NSHD is a birth cohort study that includes 5,362 individuals (2,547 males and 2,815 females) who were born in one week in March 1946 in Britain. This cohort has been extensively followed up with periodic assessments of various aspects of their lives, including anthropometric measurements, socio-economic status, lifestyle factors, and health outcomes .
Between 2006 and 2010, when NSHD participants were 60–64 years old, those residing in the UK who had not been lost to follow-up or withdrawn were invited to attend a clinic-based assessment which included a standard 12-lead surface ECG. The 10-sec ECGs were stored in digital format to avoid the signal quality degradation that affects paper traces.
Based on the earlier method of Ball et al. , Lindow et al.  recently implemented a machine learning-based Bayesian-centric approach to predict cardiac age from multiple discrete features derivable from standard 12-lead ECGs, combining inputs from: (1) conventional ECG durations (e.g., P and QTc), amplitudes and axes (e.g., QRS and T); (2) the spatial QRS-T angles, spatial ventricular gradient, spatial QRS- and T-wave axes, azimuths, elevations, velocities, waveform amplitudes and areas from the derived, Frank X,Y, and Z lead vectorcardiogram; and (3) QRS- and T-waveform complexity obtained via singular value decomposition after signal averaging. Univariable linear regression models were used to select the ECG features, and multivariable linear regression models to estimate the cardiac age from the original Bayesian A-ECG heart age model of Ball et al. that had used higher-fidelity 5-min A-ECGs . The model of Lindow et al.  was utilized in the present study to derive the estimated Bayesian A-ECG cardiac age for each member of our own cohort from their respective, standard 12-lead ECGs.
The results from two logistic regression-based A-ECG scores for cardiac diseases were also evaluated: one for left ventricular (LV) systolic dysfunction (LVSD), and one for LV electrical remodeling (LVER) [33,34,35,36,37,38]. Both the LVSD and the LVER A-ECG scores prominently incorporate results from the spatial QRS-T angle, a measure also known to have important prognostic utility . While the results from these scores comprise continuous variables (as utilized in this study) rather than categorical variables, the presence of the given disease by A-ECG is usually also clinically defined as the score’s related probability exceeding 0.5 (50%). The presence of a positive LVSD score was originally designed to correspond to an imaging-proven left ventricular ejection fraction (LVEF) < 50% . However, ongoing research and subsequent clinical practice suggest that the LVSD score slightly better correlates with changes in global longitudinal strain than in LVEF [34, 35]. The presence of a positive LVER score was in turn designed to more accurately predict (versus strictly conventional ECG criteria) the presence of moderate or greater left ventricular hypertrophy by gold-standard imaging [33, 37, 38].
Finally, among others, Ribeiro et al. [13, 14] have also recently designed a deep neural network (DNN) to estimate the cardiac age based on raw, standard 12-lead ECG tracings. The DNN of Ribeiro et al. consists of 11 convolutional layers (with the last 10 organized into 5 blocks) whose weights were initialized via random sampling from a scaled normal variable. It was trained on a dataset of 1,558,415 patients from 811 counties in the state of Minas Gerais (Brazil) collected by the Telehealth Network of Minas Gerais (TNMG). Batch normalization was employed to rescale the output of each convolutional layer before being fed into a rectified linear activation unit, with an Adam optimizer being used to minimize the mean square error. Ribeiro et al. have made their DNN model publicly available [13, 14]. Thus DNN ECG-based cardiac ages were also derived for our own cohort by using Ribeiro et al.’s publicly accessible algorithm.
DNA methylation-based aging biomarkers
Blood samples from NSHD study members were collected in 1999 and again between 2006 and 2010 as previously described . Illumina Infinium Methylation EPIC BeadChips kits (Illumina, San Diego, California, US) were used to measure DNAm signals at > 850,000 CPG sites. The signals were processed for quality control (QC) in the R Enmix package , and beta-values were obtained using the noob normalization method in the R minfi package . Signals with a detection p-value > 10–6 and a number of beads < 3 were set to missing. We excluded: (1) samples with missing data in > 5% of the CpGs; (2) CpGs with missing data in > 5% of the samples; and (3) samples with outliers in bisulfite intensity, total intensity, or beta-values. Outliers were defined as values more than 3 standard deviations (SDs) from the mean or 3 interquartile ranges (IQRs) below the 1st or above the 3rd quartiles. Sample identity was verified by calculating the Pearson correlation coefficient between the 59 methylation bead chips SNPs and the imputed genotype data yielding coefficients > 0.90.
We calculated DNAm AgeHannum, DNAm AgeHorvath, DNAm PhenoAge and DNAm GrimAge using the methodology described by Horvath which is also available as online age calculator (https://dnamage.genetics.ucla.edu/home) . Age acceleration (in years) was defined as the residual produced by linearly regressing DNAm age on the chronological age , yielding the corresponding AgeAccelHannum, AgeAccelHorvath, AgeAccelPheno and AgeAccelGrim.
Cardiometabolic risk factors and covariates
Sex was assigned at birth as male or female. During the same clinic visit when ECG was recorded, participants’ weight and height were also measured and used to compute body mass index (BMI). Participants’ socioeconomic position (SEP) was evaluated at the time of echocardiography (60–64 years) or at 53 years where the former was not available, according to the UK Office of Population Censuses and Surveys Registrar General’s occupational-based social class dichotomized as manual or non-manual. Self-reported questionnaires at 60–64 years also provided information about smoking status (never smoked, ex-smoker and currently smoking), the average number of units of alcohol consumed per day and leisure time physical activity. The latter was dichotomized as inactive or active (exercises at least once per month). Blood samples at 60–64 were analyzed to provide the white cell counts: naïve and exhausted CD8+ T-lymphocytes, CD4+ T-lymphocytes, B-cells, natural killer cells, granulocytes, and monocytes. The presence of CVD, diabetes, high cholesterol or hypertension was recorded as 1 = present or 0 = absent as previously described . In addition, we defined the presence of ‘any CV risk factor’ as any participant having at least one out of diabetes, high cholesterol, hypertension, CVD, or a BMI > 30.
Statistical analysis was performed in R (version-4.2.1), and a two-tailed p value < 0.05 was considered statistically significant. Distributions of data were assessed on histograms and using the Shapiro–Wilk test. Continuous variables were expressed as mean ± 1 SD or median (IQR) as appropriate; and categorical variables, as counts and percentages.
In all analyses, the DNAm age or AgeAccel were the independent variables, whilst the A-ECG-based and DNN ECG-based cardiac ages or the A-ECG-based LVSD and LVER scores were the outcomes. To test for associations, we employed generalized linear models (glms) with Gaussian distributions and identity links. To mitigate the influence of confounders, the models were adjusted for chronological age, sex, SEP, smoking, alcohol consumption, and physical activity. To obtain better estimates of the actual epigenetic changes, the models were also adjusted for the white cell counts. We repeated the analyses for each DNAm age (i.e., AgeHannum, AgeHorvath, PhenoAge and GrimAge) as well as for the corresponding AgeAccel.
Firstly, we explored the associations between the DNAm ages and the A-ECG and DNN ECG-related outcomes. Secondly, we explored to what extent DNAm AgeAccel mediates the effect of CV risk factors (BMI, diabetes, high cholesterol, hypertension, CVD and ‘any CV risk factor’) on the A-ECG and DNN ECG-related outcomes using mediation analysis. The hypothetical mediation mechanism is presented in Fig. 1. We used the counterfactual framework methodology of causal inference developed by Imai, Tingley and Yamamoto which relies on the no-interaction (i.e., no exposure-mediator and mediator-outcome interactions) and sequential ignorability (SI; i.e., the absence of unmeasured confounding) assumptions [44,45,46]. To calculate the total effects (the direct model), we regressed the CV risk factors on the A-ECG and DNN ECG-related outcomes. To calculate the effect of the independent variable onto the mediator (the mediator model), we regressed the CV risk factors on the DNAm AgeAccel. To calculate the effect of the mediator on the dependent variable (the outcome model), we regressed the DNAm AgeAccel on the A-ECG and DNN ECG-related outcomes whilst adjusting for the CV risk factors. Then, we derived the total effects, average causal mediation effects (ACMEs) and average direct effects (ADEs) with their 95% confidence intervals (CIs) using nonparametric bootstrapping with 1000 Monte Carlo simulations . The proportion mediated (PM) was derived by dividing ACME by the total effect for each analysis. To investigate the robustness of our results, we used the Baron-Kenny procedure to calculate the correlation coefficients (ρ) between the residuals of the mediator and the regression outcomes from the linear structural equation models (LSEM) rather than glms as this was the only available implementation in R. We judged the robustness of the results by the magnitude of ρ required to reverse the sign of ACME.
Participant characteristics for the whole cohort, along with those stratified by sex, are showcased in Table 1. On average, male participants were more likely to have a non-manual-labor SEP (34.9% vs. 26.0%, p = 0.039), high cholesterol (29.0% vs. 13.4%, p < 0.0001), hypertension (57.5% vs. 43.9%, p = 0.003) and CVD (10.7% vs. 3.3%, p = 0.005). Conversely, females were more likely to have a lower DNAm AgeHannum (49.9 vs. 52.9, p < 0.0001), DNAm AgeHorvath (57.3 vs. 58.7, p < 0.0001), DNAm GrimAge (60.0 vs. 63.0, p < 0.0001) and hence a higher AgeAccel in these DNAm ages. Males were more likely to have a worse DNN ECG-related cardiac ages (62.7 vs. 59.6, p = 0.005), and females a worse Bayesian A-ECG cardiac age (70.2 vs. 68.0, p < 0.0001). Males had a worse LVSD score, whilst the LVER score was not statistically different based on sex.
Associations between DNAm ages and ECG-based cardiac ages and disease scores
The DNAm AgeHannum, AgeHorvath, PhenoAge and GrimAge were weakly correlated with both the Bayesian A-ECG and DNN ECG-based ages as well as with the LVSD and LVER disease scores (Table 2). A 1-year increase in the DNAm AgeHorvath or PhenoAge was associated with a ≈ 0.1 years (both p < 0.023) increase in the Bayesian A-ECG cardiac age, whilst a 1-year increase in DNAm GrimAge was associated with 0.21 years ([0.05, 0.37], p = 0.009) increase. Similarly, a 1-year increase in the DNAm AgeHorvath, PhenoAge or GrimAge was associated with a 0.3, 0.2, and 0.4 years respectively, increase in the DNN-ECG-based age. A higher DNAm AgeHorvath, PhenoAge and GrimAge were associated with worse LVSD and LVER scores.
DNAm as a mediator downstream of CV risk factors and upstream of ECG-based cardiac ages
Per Table 3, a 1-unit increase in BMI resulted in 0.32 years ([0.21, 0.44], p < 0.001) increase in the Bayesian A-ECG cardiac age. Moreover, having diabetes, high cholesterol, hypertension, or any CV risk factor, resulted in an increase of 1.30 ([0.09, 2.62] p = 0.034), 2.31 ([0.89, 3.57] p < 0.001), 2.67 ([1.70, 3.74] p < 0.001) and 2.70 years ([1.61, 3.74] p < 0.001) respectively, in the Bayesian A-ECG cardiac age. When adjusting for the AgeAccelPheno, the ACME for diabetes was 0.23 years ([0.01, 0.52] p = 0.028) and PM ≈18%, for high cholesterol 0.34 ([0.03, 0.74] p = 0.024) and PM ≈15%, and for any CV risk factor 0.17 ([0.01, 0.41] p = 0.048) and PM ≈ 6%. Similarly, AgeAccel Grim appears to mediate ≈19% of the relationship between diabetes and Bayesian A-ECG, and ≈7% of the hypertension-Bayesian A-ECG association. However, there was no significant ACMEs for AgeAccelHannum or AgeAccelHorvath.
A unit increase in BMI resulted in a 0.59 years ([0.37, 0.82] p < 0.001) increase in the DNN ECG-based cardiac age. Moreover, having hypertension, CVD or any CV risk factor resulted in an increase of 2.72 ([0.43, 5.01] p = 0.020), 4.16 ([0.18, 8.69] p = 0.040) or 3.20 years ([0.84, 5.74] p = 0.008) respectively, in the DNN ECG-based cardiac age. Although diabetes and high cholesterol were associated with higher DNN ECG-based cardiac ages, these were not statistically significant. When adjusting for AgeAccelGrim, the ACME for diabetes was 0.48 years ([0.02, 1.22] p = 0.028) and PM ≈37%, for high cholesterol 0.76 ([0.04, 1.71] p = 0.036) and PM ≈33%, and hypertension 0.46 ([0.01, 1.07] p = 0.046) and PM ≈ 17%. However, there were no significant ACMEs for AgeAccelHannum, AgeAccelHorvath, or AgeAccelPheno.
DNAm as a mediator downstream of CV risk factors and upstream of A-ECG disease scores
In general, when considering the associations between CV risk factors and the A-ECG-based LVSD and LVER scores, only AgeAccelPheno and AgeAccelGrim appear to be mediators (Table 4). AgeAccelPheno appears to mediate ≈15% of the relationships between BMI and the two A-ECG-based risk scores, ≈28% of the association between diabetes and the LVSD score, and ≈17% of the BMI-LVER score association. Similarly, AgeAccelGrim emerged as a significant mediator for the associations between BMI, diabetes, high cholesterol and hypertension, and the two A-ECG-based risk scores (all p < 0.05 for ACMEs). Importantly, AgeAccelGrim mediates ≈40% of the relationships between diabetes, high cholesterol and hypertension, and LVSD; and 30% of their assocaition with LVER disease score.
The ρ at which ACME = 0 is presented in Additional file 1: Supplementary Table S1. In general, even small deviations from the SI assumption can reverse the sign of ACME.
In this cross-sectional analysis, we show that the association between CV risk factors and ECG-based cardiac ages and disease scores could be partly mediated by the 2nd generation DNAm AgeAccel biomarkers. AgeAccelPheno and AgeAccelGrim appear to mediate the relationships between most CV risk factors and the Bayesian A-ECG cardiac age, and the LVSD and LVER disease scores AgeAccelGrim also appears to mediate the relationship between most CV risk factors and the DNN ECG-based cardiac age. However, this mediation appears to account for ≈10–40% of the total effects. AgeAccelHannum and AgeAccelHorvath appear to have a limited role.
Although epigenetics is a broad topic, our understanding revolves mainly on DNAm at CpG sites which are concentrated in the promoter regions of the genes albeit sparse in other parts of the genome. In general, promoter methylation at CpG sites can lead to gene silencing, whilst unmethylated promoters remain transcriptionally active. Previous epigenome-wide association studies have highlighted DNAm profiles associated with cardiometabolic risk factors (e.g., BMI , diabetes , high cholesterol , and hypertension ). In this study we show that DNAm AgeAccel derived from CpG methylation at specific sites in different tissues (especially blood) is a partial mediator downstream of the CV risk factors but upstream of the A-ECG phenotypes. This reinforces the theory that CV risk factors can lead to dynamic DNA changes with potentially adverse long-term cardiac phenotypic sequelae. Indeed, DNMT-i showed promising results in pre-clinical trials as they reduced pathological hypertension-related myocardial hypertrophy  and fibrosis , and ameliorated atherosclerosis . Interestingly, aspirin which has a proven benefit especially in the secondary prevention of CVDs, might also exert some of its effects by acting as a DNMT-i . In general, DNMT-i have a reversible effect suggesting that their safety, efficacy, and effectiveness in protecting against pathological age-related remodeling in human clinical trials might be possible to explore .
Although the 1st generation DNAm ages appeared to be only weak predictors of CVDs , the 2nd generation DNAm ages incorporating clinical and physiological prognostic methylation biomarkers were intended to act as better markers of healthspan (DNAm PhenoAge) and lifespan (DNAm GrimAge). Our findings suggest that only the 2nd generation DNAm AgeAccel metrics could act as mediators of the association between the CV risk factors and the ECG-based cardiac outcomes. DNAm GrimAge incorporates blood-based biomarkers related to extracellular matrix (ECM) remodeling (e.g., epidermal growth factor-containing fibulin-like ECM protein, plasminogen activator inhibitor 1, tissue inhibitor metalloprotease 1 etc.) . As LVSD and LVER are characterized by extensive remodeling, it is not surprising that DNAm AgeAccelGrim emerged as a significant mediator downstream of CV risk factors but upstream of these A-ECG disease scores. In contrast, DNAmPheno encompasses CpG sites related to immune system (e.g., white blood cell count), inflammation (e.g., C-reactive protein) and metabolism (e.g., glucose). Inflammaging  and immunosenescence  have both been associated with faster cardiac aging and increased susceptibility to CVDs as they augment endothelial damage, impair tissue repair, and promote insulin resistance and atherosclerosis . Given this strong relationship between immunity and metabolism, and cardiac ageing, DNAm AgeAccelPheno emerged as a consistent mediator in the association between CV risk factors and A-ECG ages and disease scores in our study (Fig. 2).
According to our results, only 10–40% of the total effects of the association between CV risk factors and A-ECG phenotypes are mediated by DNAm AgeAccel. However, the proportion mediated as a numerical quantity provides limited insights into the relative importance of the DNAm pathway. Whether DNAm has only a limited impact or whether it is a major stem which sets into motion a raft of cascading pathological cardiac ageing pathways, remains to be further elucidated. Moreover, the directionality of effect could not be firmly deduced using the analyses provided in this study. While it is theoretically possible that a more advanced cardiac age leads to a higher DNAm age, it is more biologically plausible that CV risk factors induce physiological stress driving DNAm. The advent of the International Human Epigenome Consortium promises to provide novel insights into the epigenetic changes most strongly associated with cardiac ageing . As CpG methylation is closely linked to gene transcription, the transcriptome and proteome, now measurable via high-throughput RNA sequencing and shotgun proteomics, respectively, will bolster our understanding of the mechanisms downstream of DNAm but upstream of the cardiac ageing phenotypes.
Strengths and limitations
NSHD was representative of a British-born population at the time of participant recruitment. The implicit age homogeneity of the birth cohort participants was another strength of the study, as it enables age-matching across analyses meaning that age-related confounding was minimized. In addition, participants were exposed to similar secular trends and risk factors, and similar access to diagnostic technologies and treatment facilities over time, which minimizes the bias of environmental factors.
An important limitation is that only participants who had DNAm and analyzable 12-lead ECGs which were collected as part of two different separate sub-studies were included in this study. Recently, the updated version of DNAm GrimAge was published, but this was not available at the time of the current analysis . As this study was retrospectively designed, selection bias may have influenced the observed associations. Moreover, the study has all the limitations inherent to cross-sectional studies (e.g., antecedent-consequent bias, susceptibility to transient effects etc.). In addition, the limited sample size meant that we were underpowered to significantly detect certain associations, with statistically significant ACMEs but not total effects found in some analyses. In that case, although we could not claim an association between the CV risk factor and the corresponding ECG-based outcome, if indeed there is one it is probably being mediated by the DNAm variable. Although sex-specific differences between DNAm and CV health have been reported, this study was underpowered for sex-stratified analyses, so they were not pursued. Moreover, we assumed that no exposure-mediator or mediator-outcome interactions exists. Whist mediation frameworks taking into account interactions exists (e.g., VanderWeel’s four-way decomposition ), they are more suitable for studies with larger sample sizes. Similarly, the results are sensitive to the violation of the SI assumption (Additional file 1: Supplementary Table S1) and the existence of unexplored confounding cannot be excluded although we adjusted for chronological age, sex, SEP, smoking, alcohol, and physical activity. Lastly, repeating the analyses using the sub-components of the DNAm ages would have provided a more comprehensive understanding of the association between CV risk factors and the ECG-outcomes, and this represents our plans for future work.
By the age of 60, individuals with accelerated DNA methylation appear to have older, weaker, and more electrically impaired hearts. The harmful effects of CV risk factors on cardiac age and health, appear to be partially mediated by the 2nd generation DNA methylation age biomarkers. This highlights the need for more research into the potentially cardioprotective roles of selective DNA methyltransferases modulators.
Availability of data and materials
NSHD data are available from: https://www.nshd.mrc.ac.uk/data. The mediation analysis code template used can be accessed here: https://cran.r-project.org/web/packages/mediation/mediation.pdf.
Harman D. Aging: overview. Ann N Y Acad Sci. 2001;928:1–21. https://doi.org/10.1111/j.1749-6632.2001.tb05631.x.
Lowsky DJ, Olshansky SJ, Bhattacharya J, Goldman DP. Heterogeneity in healthy aging. J Gerontol A Biol Sci Med Sci. 2014;69:640–9. https://doi.org/10.1093/gerona/glt162.
Jylhävä J, Pedersen NL, Hägg S. Biological age predictors. EBioMedicine. 2017;21:29–36. https://doi.org/10.1016/j.ebiom.2017.03.046.
Hannum G, Guinney J, Zhao L, Zhang L, Hughes G, Sadda S, Klotzle B, Bibikova M, Fan J-B, Gao Y, et al. Genome-wide methylation profiles reveal quantitative views of human aging rates. Mol Cell. 2013;49:359–67. https://doi.org/10.1016/j.molcel.2012.10.016.
Horvath S. DNA methylation age of human tissues and cell types. Genome Biol (Online Edn). 2013;14:R115–R115. https://doi.org/10.1186/gb-2013-14-10-r115.
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. https://doi.org/10.18632/aging.101414.
Lu AT, Quach A, Wilson JG, Reiner AP, Aviv A, Raj K, Hou L, Baccarelli AA, Li Y, Stewart JD, et al. DNA methylation GrimAge strongly predicts lifespan and healthspan. Aging (Albany, NY). 2019;11:303–27. https://doi.org/10.18632/aging.101684.
Maddock J, Castillo-Fernandez J, Wong A, Cooper R, Richards M, Ong KK, Ploubidis GB, Goodman A, Kuh D, Bell JT, et al. DNA methylation age and physical and cognitive ageing. 2019.
Maddock J, Castillo-Fernandez J, Wong A, Ploubidis GB, Kuh D, Bell JT, Hardy R. Childhood growth and development and DNA methylation age in mid-life. 2021.
Pavanello S, Campisi M, Fabozzo A, Cibin G, Tarzia V, Toscano G, Gerosa G. The biological age of the heart is consistently younger than chronological age. Sci Rep. 2020;10:10752–10752. https://doi.org/10.1038/s41598-020-67622-1.
Choi EJ, Lee SY. Index of Cardiac Age Index and cardiometabolic biomarkers in adults with metabolic syndrome. Metab Syndr Relat Disord. 2021;19:288–95. https://doi.org/10.1089/met.2020.0095.
Raisi-Estabragh Z, Salih A, Gkontra P, Atehortúa A, Radeva P, Boscolo Galazzo I, Menegaz G, Harvey NC, Lekadir K, Petersen SE. Estimation of biological heart age using cardiovascular magnetic resonance radiomics. Sci Rep. 2022;12:12805–12805. https://doi.org/10.1038/s41598-022-16639-9.
Ribeiro AH, Ribeiro MH, Paixão GMM, Oliveira DM, Gomes PR, Canazart JA, Ferreira MPS, Andersson CR, Macfarlane PW, Wagner M Jr, et al. Automatic diagnosis of the 12-lead ECG using a deep neural network. Nat Commun. 2020;11:1760–1760. https://doi.org/10.1038/s41467-020-15432-4.
Lima EM, Ribeiro AH, Paixão GMM, Ribeiro MH, Pinto-Filho MM, Gomes PR, Oliveira DM, Sabino EC, Duncan BB, Giatti L, et al. Deep neural network-estimated electrocardiographic age as a mortality predictor. Nat Commun. 2021. https://doi.org/10.1038/s41467-021-25351-7.
Lindow T, Palencia-Lamela I, Schlegel TT, Ugander M. Heart age estimated using explainable advanced electrocardiography. Sci Rep. 2022;12:9840–9840. https://doi.org/10.1038/s41598-022-13912-9.
Quach A, Levine ME, Tanaka T, Lu AT, Chen BH, Ferrucci L, Ritz B, Bandinelli S, Neuhouser ML, Beasley JM, et al. Epigenetic clock analysis of diet, exercise, education, and lifestyle factors. Aging (Albany, NY). 2017;9:419–46. https://doi.org/10.18632/aging.101168.
Zannas AS, Arloth J, Carrillo-Roa T, Iurato S, Röh S, Ressler KJ, Nemeroff CB, Smith AK, Bradley B, Heim C, et al. Lifetime stress accelerates epigenetic aging in an urban, African American cohort: Relevance of glucocorticoid signaling. Genome Biol. 2015;16:266–266. https://doi.org/10.1186/s13059-015-0828-5.
Ward-Caviness CK, Nwanaji-Enwerem JC, Wolf K, Wahl S, Colicino E, Trevisi L, Kloog I, Just AC, Vokonas P, Cyrys J, et al. Long-term exposure to air pollution is associated with biological aging. Oncotarget. 2016;7:74510–25. https://doi.org/10.18632/oncotarget.12903.
Joehanes R, Just AC, Marioni RE, Pilling LC, Reynolds LM, Mandaviya P, Guan WH, Xu T, Elks CE, Aslibekyan S, et al. Epigenetic signatures of cigarette smoking. Circ Cardiovasc Genet. 2016;9:436–47. https://doi.org/10.1161/CIRCGENETICS.116.001506.
Liu C, Marioni RE, Hedman Å, Pfeiffer L, Tsai PC, Reynolds LM, Just AC, Duan Q, Boer CG, Tanaka T, et al. A DNA methylation biomarker of alcohol consumption. 2018.
Niiranen TJ, Vasan RS. Epidemiology of cardiovascular disease: recent novel outlooks on risk factors and clinical approaches. Expert Rev Cardiovasc Ther. 2016;14:855–69. https://doi.org/10.1080/14779072.2016.1176528.
Roetker NSNS, Pankow JSJS, Bressler JJ, Morrison ACAC, Boerwinkle EE. A prospective study of epigenetic age acceleration and incidence of cardiovascular disease outcomes in the Atherosclerosis Risk in Communities (ARIC) Study. Circul Genom Precis Med. 2018;11:e001937–e001937. https://doi.org/10.1161/CIRCGEN.117.001937.
Dick KJP, Nelson CPP, Tsaprouni LP, Sandling JKP, Aïssi DM, Wahl SM, Meduri EP, Morange P-EP, Gagnon FP, Grallert HP, et al. DNA methylation and body-mass index: a genome-wide analysis. Lancet (British Edition). 2014;383:1990–8. https://doi.org/10.1016/S0140-6736(13)62674-4.
Chambers JCD, Loh MP, Lehne BP, Drong AP, Kriebel JP, Motta VP, Wahl SM, Elliott HRDB, Rota FP, Scott WRM, et al. Epigenome-wide association of DNA methylation markers in peripheral blood from Indian Asians and Europeans with incident type 2 diabetes: a nested case-control study. Lancet Diabetes Endocrinol. 2015;3:526–34. https://doi.org/10.1016/S2213-8587(15)00127-8.
Irvin MR, Zhi D, Joehanes R, Mendelson M, Aslibekyan S, Claas SA, Thibeault KS, Patel N, Day K, Jones LW, et al. Epigenome-wide association study of fasting blood lipids in the genetics of lipid-lowering drugs and diet network study. Circulation (New York, NY). 2014;130:565–72. https://doi.org/10.1161/CIRCULATIONAHA.114.009158.
Roberts ML, Kotchen TA, Pan X, Li Y, Yang C, Liu P, Wang T, Laud PW, Chelius TH, Munyura Y, et al. Unique associations of DNA methylation regions with 24-hour blood pressure phenotypes in black participants. Hypertension (Dallas, Tex 1979). 2022;79:761–72. https://doi.org/10.1161/HYPERTENSIONAHA.121.18584.
Xia Y, Brewer A, Bell JT. DNA methylation signatures of incident coronary heart disease: findings from epigenome-wide association studies. Clin Epigenet. 2021. https://doi.org/10.1186/s13148-021-01175-6.
Stenzig J, Schneeberger Y, Löser A, Peters BS, Schaefer A, Zhao R-R, Ng SL, Höppner G, Geertz B, Hirt MN, et al. Pharmacological inhibition of DNA methylation attenuates pressure overload-induced cardiac hypertrophy in rats. J Mol Cell Cardiol. 2018;120:53–63. https://doi.org/10.1016/j.yjmcc.2018.05.012.
Watson CJ, Horgan S, Neary R, Glezeva N, Tea I, Corrigan N, McDonald K, Ledwidge M, Baugh J. Epigenetic therapy for the treatment of hypertension-induced cardiac hypertrophy and fibrosis. J Cardiovasc Pharmacol Ther. 2016;21:127–37. https://doi.org/10.1177/1074248415591698.
Cao Q, Wang X, Jia L, Mondal AK, Diallo A, Hawkins GA, Das SK, Parks JS, Yu L, Shi H, et al. Inhibiting DNA methylation by 5-aza-2′-deoxycytidine ameliorates atherosclerosis through suppressing macrophage inflammation. Endocrinology (Philadelphia). 2014;155:4925–38. https://doi.org/10.1210/en.2014-1595.
Kuh D, Pierce M, Adams J, Deanfield J, Ekelund U, Friberg P, Ghosh AK, Harwood N, Hughes A, Macfarlane PW, et al. Cohort profile: updating the cohort profile for the MRC National Survey of Health and Development: a new clinic-based data collection for ageing research. Int J Epidemiol. 2011;40(1):e1–9.
Ball R, Feiveson A, Schlegel T, Starc V, Dabney A. Predicting, “Heart Age” using electrocardiography. J Person Med. 2014;4:65–78. https://doi.org/10.3390/jpm4010065.
Schlegel TT, Kulecz WB, Feiveson AH, Greco EC, Depalma JL, Starc V, Vrtovec B, Rahman MA, Bungo MW, Hayat MJ, et al. Accuracy of advanced versus strictly conventional 12-lead ECG for detection and screening of coronary artery disease, left ventricular hypertrophy and left ventricular systolic dysfunction. BMC Cardiovasc Disord. 2010;10:28. https://doi.org/10.1186/1471-2261-10-28.
Gleeson S, Liao Y-W, Dugo C, Cave A, Zhou L, Ayar Z, Christiansen J, Scott T, Dawson L, Gavin A, et al. ECG-derived spatial QRS-T angle is associated with ICD implantation, mortality and heart failure admissions in patients with LV systolic dysfunction. PLoS ONE. 2017;12:e0171069. https://doi.org/10.1371/journal.pone.0171069.
Gladding PA, Loader S, Smith K, Zarate E, Green S, Villas-Boas S, Shepherd P, Kakadiya P, Hewitt W, Thorstensen E, et al. Multiomics, virtual reality and artificial intelligence in heart failure. Future Cardiol. 2021;17:1335–47. https://doi.org/10.2217/fca-2020-0225.
Johnson K, Neilson S, To A, Amir N, Cave A, Scott T, Orr M, Parata M, Day V, Gladding P. Advanced electrocardiography identifies left ventricular systolic dysfunction in non-ischemic cardiomyopathy and tracks serial change over time. J Cardiovasc Dev Disease. 2015;2:93–107. https://doi.org/10.3390/jcdd2020093.
Maanja M, Schlegel TT, Kozor R, Lundin M, Wieslander B, Wong TC, Schelbert EB, Ugander M. The electrical determinants of increased wall thickness and mass in left ventricular hypertrophy. J Electrocardiol. 2020;58:80–6. https://doi.org/10.1016/j.jelectrocard.2019.09.024.
Maanja M, Schlegel TT, Kozor R, Bacharova L, Wong TC, Schelbert EB, Ugander M. Improved evaluation of left ventricular hypertrophy using the spatial QRS-T angle by electrocardiography. Sci Rep. 2022. https://doi.org/10.1038/s41598-022-16712-3.
Oehler A, Feldman T, Henrikson CA, Tereshchenko LG. QRS-T angle: a review. Ann Noninvasive Electrocardiol. 2014;19:534–42. https://doi.org/10.1111/anec.12206.
Xu Z, Niu L, Li L, Taylor JA. ENmix: a novel background correction method for Illumina HumanMethylation450 BeadChip. Nucleic Acids Res. 2016;44:e20–e20. https://doi.org/10.1093/nar/gkv907.
Fortin J-P, Triche TJ, Hansen KD. Preprocessing, normalization and integration of the Illumina HumanMethylationEPIC array with minfi. Bioinformatics. 2017;33:558–60. https://doi.org/10.1093/bioinformatics/btw691.
Horvath S, Gurven M, Levine ME, Trumble BC, Kaplan H, Allayee H, Ritz BR, Chen B, Lu AT, Rickabaugh TM, et al. An epigenetic clock analysis of race/ethnicity, sex, and coronary heart disease. Genome Biol. 2016;17:171–171. https://doi.org/10.1186/s13059-016-1030-0.
Topriceanu C-C, Moon JC, Hardy R, Chaturvedi N, Hughes AD, Captur G. Longitudinal birth cohort study finds that life-course frailty associates with later-life heart size and function. Sci Rep. 2021. https://doi.org/10.1038/s41598-021-85435-8.
Imai K, Tingley D, Yamamoto T. Experimental Designs for Identifying Causal Mechanisms. J R Stat Soc Ser A Stat Soc. 2013;176:5–51. https://doi.org/10.1111/j.1467-985x.2012.01032.x.
Imai K, Yamamoto T. Identification and sensitivity analysis for multiple causal mechanisms: revisiting evidence from framing experiments. Polit Anal. 2013;21:141–71. https://doi.org/10.1093/pan/mps040.
Imai K, Keele L, Tingley D. A general approach to causal mediation analysis. Psychol Methods. 2010;15:309–34. https://doi.org/10.1037/a0020761.
Tingley D, Yamamoto T, Hirose K, Keele L, Imai K. Mediation: R package for causal mediation analysis. J Stat Softw. 2014;59:1–38. https://doi.org/10.18637/jss.v059.i05.
Guay S-P, Légaré C, Houde A-A, Mathieu P, Bossé Y, Bouchard L. Acetylsalicylic acid, aging and coronary artery disease are associated with ABCA1 DNA methylation in men. Clin Epigenet. 2014;6:14–14. https://doi.org/10.1186/1868-7083-6-14.
Shi Y, Zhang H, Huang S, Yin L, Wang F, Luo P, Huang H. Epigenetic regulation in cardiovascular disease: mechanisms and advances in clinical trials. Signal Transduct Target Ther. 2022;7:200–200. https://doi.org/10.1038/s41392-022-01055-2.
Franceschi C, Garagnani P, Parini P, Giuliani C, Santoro A. Inflammaging: a new immune-metabolic viewpoint for age-related diseases. Nat Rev Endocrinol. 2018;14:576–90. https://doi.org/10.1038/s41574-018-0059-4.
Pawelec G. Age and immunity: What is “immunosenescence”? Exp Gerontol. 2018;105:4–9. https://doi.org/10.1016/j.exger.2017.10.024.
Ferrucci L, Fabbri E. Inflammageing: chronic inflammation in ageing, cardiovascular disease, and frailty. Nat Rev Cardiol. 2018;15:505–22. https://doi.org/10.1038/s41569-018-0064-2.
Stunnenberg HG, Hirst M, Abrignani S, Adams D, De Almeida M, Altucci L, Amin V, Amit I, Antonarakis SE, Aparicio S, et al. The International human epigenome consortium: a blueprint for scientific collaboration and discovery. Cell. 2016;167:1145–9. https://doi.org/10.1016/j.cell.2016.11.007.
Lu AT, Binder AM, Zhang J, Yan Q, Reiner AP, Cox SR, Corley J, Harris SE, Kuo P-L, Moore AZ, et al. DNA methylation GrimAge version 2. Aging (Albany, NY). 2022;14:9484–549. https://doi.org/10.18632/aging.204434.
Vanderweele TJ. A unification of mediation and interaction. Epidemiology. 2014;25:749–61. https://doi.org/10.1097/ede.0000000000000121.
This study was funded by the UK Medical Research Council (program codes MC_UU_12019/1; MC_UU_12019/4; MC_UU_12019/5). G.C. is supported by British Heart Foundation (MyoFit46 Special Programme Grant SP/20/2/34841), the Barts Charity HeartOME1000 project grant (MGU0427 / G-001411) and by the NIHR UCL Hospitals Biomedical Research Centre. J.C.M. is directly and indirectly supported by the UCL Hospitals NIHR BRC and Biomedical Research Unit at Barts Hospital respectively. AH receives support from the British Heart Foundation, the Economic and Social Research Council (ESRC), the Horizon 2020 Framework Programme of the European Union, the National Institute on Aging, the National Institute for Health Research University College London Hospitals Biomedical Research Centre, the UK Medical Research Council and works in a unit that receives support from the UK Medical Research Council. None of the funders was involved in the study design, the collection, the analysis, the interpretation of the data, and in the decision to submit the article for publication. For the purpose of open access, the authors have applied a creative commons attribution (CC BY) license to any author accepted manuscript version arising.
Ethics approval and consent to participate
The 2006–2010 NSHD data collection sweep included an in-depth cardiovascular assessment and was granted ethical approval from the Greater Manchester Local Research Ethics Committee and the Scotland Research Ethics Committee , with written informed consent obtained from all study participants. Our project was approved by the NSHD committee. All procedures performed were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards.
The views expressed in this article are those of the authors. TTS is owner and founder of Nicollier-Schlegel SARL, which performs ECG interpretation consultancy using software that can quantify the advanced ECG measures used in the current study. TTS and MU are owners and founders of Advanced ECG Systems, a company that is developing commercial applications of advanced ECG technology used in the current study. JCM is the chief executive officer of Myocardium AI and has served on advisory boards for Genzyme and Sanofi. The remaining authors declare that they have no potential conflicts of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Topriceanu, CC., Dev, E., Ahmad, M. et al. Accelerated DNA methylation age plays a role in the impact of cardiovascular risk factors on the human heart. Clin Epigenet 15, 164 (2023). https://doi.org/10.1186/s13148-023-01576-9