Skip to main content

Cardiovascular health and four epigenetic clocks

Abstract

Background

Cardiovascular health (CVH) was defined by the American Heart Association as an integrative idealness of seven clinical or lifestyle factors. Based on populations of European ancestry, recent studies have shown that ideal CVH is associated with a slower aging rate. The aging rate is measured by levels of epigenetic age acceleration (EAA), usually obtained from the residuals of regressing DNA methylation (DNAm) age on chronological age. However, little has been known about the association of CVH with biological aging in Asian populations.

Methods and results

We here analyzed blood DNAm data and clinical/lifestyle factors of 2474 Taiwan Biobank (TWB) participants, to investigate the association of CVH with EAA. CVH was assessed by seven components: smoking status, physical activity, dietary habits, body mass index, total cholesterol, fasting glucose, and blood pressure levels. Four measures of EAA were applied, among which two were based on the first-generation DNAm clocks (HannumEAA and IEAA) and two were based on the second-generation clocks (PhenoEAA and GrimEAA). After excluding 276 individuals with cardiovascular diseases, we regressed EAA on the CVH score (ranging from 0 to 7, integrating the abovementioned seven components) while adjusting for sex, drinking status, and educational attainment.

Our results showed that a decrease in one point in the CVH score was associated with a 0.350-year PhenoEAA (p = 4.5E−4) and a 0.499-year GrimEAA (p = 4.2E−15). By contrast, HannumEAA and IEAA were not significantly associated with the CVH score. We have obtained consistent results within each generation of epigenetic clocks.

Conclusions

This is one of the first studies to comprehensively investigate the associations of CVH with four epigenetic clocks. Our TWB data showed that ideal CVH is associated with lower levels of EAA calculated according to the second-generation epigenetic clocks (PhenoEAA and GrimEAA). Having an ideal CVH status can lower EAA and reduce the risk of aging-related disorders.

Introduction

Cardiovascular health (CVH) was defined by the American Heart Association (AHA) as an integrative idealness of four lifestyle factors (cigarette smoking, physical activity, body mass index (BMI), and dietary habits) and three clinical factors (fasting glucose, blood pressure, and total cholesterol level) [1]. Ideal CVH is associated with lower risks of cardiovascular diseases (CVDs) and chronic diseases and decreased all-cause mortality [1, 2]. Ideal CVH is also associated with more extended longevity and a better quality of life [3, 4].

Epigenetic change may reflect the idealness of CVH [5, 6]. One of the epigenetic mechanisms is DNA methylation (DNAm), which can dynamically express human aging-related physiological changes [7,8,9]. Therefore, DNAm has been used to estimate human biological age [9,10,11,12,13]. Epigenetic age (i.e., the so-called DNAm age/clock) measures biological age to predict people’s healthspan and lifespan. For example, Hannum et al. constructed a DNAm age estimator (i.e., the so-called “Hannum’s clock”) from whole-blood samples by selecting 71 cytosine-phosphate-guanine (CpG) sites that are highly predictive of chronological age [10]. Horvath further developed a multi-tissue DNAm clock with an average of 96% correlation between the DNAm age and chronological age [11]. Horvath used the elastic-net regularized regression [14] to select 353 CpGs by regressing a transformed version of chronological age on 21,369 CpG sites. These CpGs were therefore used to form a DNAm age (i.e., the so-called “Horvath’s clock”) [11]. Horvath showed that an average of 36 years of age acceleration was observed in 20 different cancer types [11].

Hannum’s clock and Horvath’s clock are regarded as the first-generation DNAm clocks because they are predictive of chronological age [15, 16]. The second-generation DNAm clocks, on the other hand, aim to estimate biological age that can better reflect disease morbidity and mortality [15,16,17]. For example, Levine et al. used a two-step approach to develop the “PhenoAge” [12]. They first selected 9 biomarkers and chronological age by regressing the hazard of aging-related mortality on 43 candidate markers (including 42 clinical biomarkers and chronological age) through a regularized Cox regression model. These 9 biomarkers and chronological age were then aggregated into a “phenotypic age,” reflecting one’s physiological conditions better than chronological age [12, 16]. In the second step, Levine et al. selected 513 CpGs from 20,169 CpGs as predictors of the phenotypic age. A weighted average of these 513 CpGs formed a DNAm clock (i.e., the so-called “PhenoAge”) [12].

Another second-generation DNAm clock is “GrimAge” [13]. A total of 7 DNAm-based surrogate markers of plasma proteins and smoking pack-years were found to be associated with time-to-death. These eight sets of surrogate markers were based on 1,030 unique CpGs. GrimAge was developed by aggregating the information of these 1,030 CpGs. Therefore, GrimAge is highly predictive of disease mortality and can reflect more aging-related physiological conditions [17].

Epigenetic age acceleration (EAA) is a novel measure to evaluate human aging rate, usually quantified by the residuals of regressing DNAm age on chronological age [18,19,20,21]. In this way, the EAA calculation will be robust to different measurement platforms and normalization methods [18,19,20,21]. A positive EAA implies that an individual is aging faster than the expected aging rate regarding his/her chronological age. Recent studies have shown that ideal CVH is associated with lower EAA in European ancestry populations [5, 6]. However, little has been known about the association of CVH with biological aging in Asian people.

Here, we investigated the association of ideal CVH with four measures of EAA based on the DNAm data from 2474 Taiwan Biobank (TWB) participants. Two measures of EAA were calculated according to the first-generation DNAm clocks, among which HannumEAA was estimated from Hannum’s clock [11] and IEAA was obtained from Horvath’s clock [10]. The other two measures of EAA were calculated according to the second-generation DNAm clocks, among which PhenoEAA was gauged from PhenoAge [12] and GrimEAA was calculated from GrimAge [13]. We obtained four measures of EAA by uploading the DNAm data to the online DNAm Age Calculator created by Horvath’s laboratory (https://dnamage.genetics.ucla.edu/new). Similar to previous studies [5, 6], after excluding participants with extreme EAA levels and participants diagnosed with CVDs or leukemia, we regressed EAA (the response variable) on the CVH score (the explanatory variable) [1] while adjusting for sex, drinking status, and educational attainment. We aimed to test the associations of the CVH scores with EAA among Asian individuals without CVDs.

Both the 7-point and 14-point scales of the CVH score were applied according to the definition from the AHA [1]. However, the scoring criteria of BMI and ideal diet behaviors were adjusted to accommodate the standards of Asians (will be described in “Methods” section). Because only ~ 60% of the 2474 TWB participants provided their dietary habits, we also calculated 6-point and 12-point CVH scores by excluding the component of “ideal diet” (will be described in “Methods” section).

Results

Table 1 presents the basic characteristics of the 2,474 TWB participants stratified by tertiles of GrimEAA. The results stratified by tertiles of PhenoEAA, IEAA, and HannumEAA were shown in Additional file 1: Tables S1–S3, respectively. Participants were aged from 30 to 70 years. All seven components of the CVH score demonstrated significantly different distributions across the three GrimEAA tertile groups (p < 0.05, Table 1).

Table 1 Basic characteristics of the 2474 TWB participants stratified by tertiles of GrimEAA

The percentage of “ideal BMI” was larger in the lowest tertile (T1) than in the second (T2) or the highest tertile (T3) of GrimEAA (p < 0.001, Table 1). Similar results were also observed in PhenoEAA (p < 0.001, Additional file 1: Table S1), IEAA (p < 0.01 between T1 and T3, Additional file 1: Table S2), and HannumEAA (p < 0.01 between T1 and T2, p < 0.001 between T1 and T3, Additional file 1: Table S3). These results showed that an increased BMI value is associated with a faster biological aging rate, consistent with previous studies indicating that an increased BMI is associated with a shortened telomere length [22,23,24]. Telomere length is another well-known biological aging biomarker. A shortened telomere length can reflect an acceleration of the aging pace [25, 26].

Regarding the CVH component of “smoking status”, fewer never-smoking individuals were in the highest tertile than in the lowest tertile of GrimEAA (p < 0.001, Table 1), while more current or former smokers in the highest tertile than in the lowest tertile of GrimEAA (p < 0.001, Table 1). Similar results were also observed in the other three measures of EAA (Additional file 1: Tables S1–S3).

For the ideal fasting glucose level, more individuals were in the lowest tertile than in both the second (p < 0.01) and the highest tertile of GrimEAA (p < 0.001, Table 1). Significant differences between the lowest tertile and the second or highest tertile of GrimEAA were also observed in ideal blood pressure and ideal diet (p < 0.001, Table 1).

After excluding participants with CVDs (n = 276) or with extreme EAA values (n = 7 for HannumEAA; n = 1 for IEAA; n = 2 for PhenoEAA; n = 5 for GrimEAA; Methods), we performed multiple linear regression analysis to investigate the association between the CVH score and the four measures of EAA (HannumEAA, IEAA, PhenoEAA, and GrimEAA), respectively. Because our study sample contained no leukemia cases, no participants were excluded due to leukemia.

Table 2 shows the results of regressing the four measures of EAA on four different definitions of the CVH score while adjusting for sex, drinking status, and educational attainment (coefficients of covariates are presented in Additional file 1: Table S5). All models consistently demonstrated inverse associations between the CVH scores and the levels of EAA, as we can see from the negative regression coefficients of the CVH scores (Beta in Table 2). However, only EAA derived from the second-generation DNAm clocks (i.e., PhenoEAA and GrimEAA) led to significant inverse associations with all four definitions of the CVH score (p ≤ 4.5E−4).

Table 2 Results of regressing four measures of EAA on the CVH score

For the 7-point CVH score, a decrease in one point in the CVH score was associated with a 0.350-year PhenoEAA (p = 4.5E−4, 95% CI: 0.1550–0.5459; Table 2) and a 0.499-year GrimEAA (p = 4.2E−15, 95% CI: 0.3758–0.6222; Table 2). For the 14-point CVH score, a decrease in one point in the CVH score was associated with a 0.268-year PhenoEAA (p = 2.4E−5, 95% CI: 0.1439–0.3919; Table 2) and a 0.364-year GrimEAA (p = 1.5E−19, 95% CI: 0.2865–0.4419; Table 2).

By contrast, IEAA (p = 0.177 for 7-point; p = 0.120 for 14-point) and HannumEAA (p = 0.108 for 7-point; p = 0.087 for 14-point) were not significantly associated with the CVH scores (Table 2). We obtained similar results even when the CVH score was defined by six components (i.e., excluding diet information largely missing in the TWB data), as shown in the bottom two rows of Table 2. For the 6-point CVH score, a decrease in one point in the CVH score was associated with a 0.388-year PhenoEAA (p = 3.9E−6, 95% CI: 0.2238–0.5528) and a 0.526-year GrimEAA (p = 6.1E−23, 95% CI: 0.4222–0.6289). For the 12-point CVH score, a decrease in one point in the CVH score was associated with a 0.278-year PhenoEAA (p = 9.6E−8, 95% CI: 0.1761–0.3798) and a 0.377-year GrimEAA (p = 2.1E−30, 95% CI: 0.3136–0.4407). Again, IEAA and HannumEAA were not significantly associated with the six-component CVH scores (Table 2).

Discussion

In this study, we evaluated the associations of CVH with four measures of EAA through DNAm, clinical, and lifestyle data of 2,474 TWB participants. We found that, among the participants free of CVD, higher CVH scores were associated with lower EAA calculated based on PhenoEAA [12] and GrimEAA [13] (second-generation epigenetic age). In contrast, the associations were not significant between the CVH scores and EAA calculated based on Horvath’s [11] and Hannum’s [10] clocks (first-generation epigenetic age). This result was robust to various definitions of the CVH scores (7-point, 14-point, 6-point, and 12-point scales).

The lack of significance of the association when shifting the EAA measurement from the second-generation epigenetic clocks to the first-generation clocks may result from the different purposes of the two generations of clocks. The second-generation epigenetic clocks were trained to reflect aging-related physiological conditions, whereas the first-generation clocks aimed to predict human chronological age [16]. Increased chronological age is an important risk factor for many diseases and cancers. However, it does not account for the heterogeneity of physiological complexity among individuals of the same age. The EAA calculated based on the first-generation clocks would tend to “inherit” this “ignorance of physiological conditions.” By contrast, the second-generation clocks-based EAA is more sensitive to the impacts of CVH [13].

This is one of the first studies to investigate the association of CVH with EAA in Asian populations (specifically, the Taiwanese population). As most studies of DNAm age were focused on European ancestry populations [10,11,12,13, 15, 17, 27], our DNAm data from the 2,474 TWB participants were unique as far as the ethnicity is concerned.

Recently, Pottinger et al. [5] investigated the associations of the CVH score with the first-generation clocks-based EAA [10, 11] in 2,170 postmenopausal women (aged from 50 to 79 years) from the Women’s Health Initiative (WHI) cohort [28]. Unlike our results, Pottinger et al. [5] showed that ideal CVH was associated with the first-generation clocks-based EAA [10, 11]. The differences between the two data sets (WHI [28] and TWB) may result in these inconsistent results. First, among the 2,170 individuals analyzed by Pottinger et al. [5], 50% were non-Hispanic white, 25% were African Americans, and 15% were Hispanics. In contrast, most TWB individuals were of Han Chinese ancestry [29]. Second, the 2,170 individuals in Pottinger et al.’s study were all postmenopausal women aging from 50 to 79 years [5], whereas our data included 2,474 men (50.24%) or women (49.76%) aging from 30 to 70 years. The TWB data presented a broader spectrum of sex and age than the WHI cohort [28].

Joyce et al. also investigated the associations of the CVH score with GrimEAA at 3,224 participants from two cohort studies [6]. DNAm data were measured from 2,106 offspring participants of the Framingham Heart Study (FHS) at exam 8 (2005–2008) and 1,118 randomly-selected individuals of Coronary Artery Risk Development in Young Adults (CARDIA) study at Y15 (2000–2001) and Y20 (2005–2006), respectively. Their cross-sectional and longitudinal analyses both showed significant associations between the CVH score and GrimEAA (p < 0.01). Our results were in line with their study, indicating that the association of CVH with GrimEAA can be replicated in the Taiwanese population. Moreover, we here investigated the other three epigenetic clocks as well (Hannum et al.’s clock [10], Horvath’s clock [11], and Levine et al.’s PhenoAge [12]), and we have obtained consistent results within each generation of epigenetic clocks. Our study demonstrates that, for even CVD-free individuals, ideal CVH is important with its link to a slower aging rate in an Asian population.

Conclusions

This is one of the first studies to comprehensively investigate the associations of CVH with four epigenetic clocks. Ideal CVH is associated with lower levels of EAA calculated according to the second-generation epigenetic clocks in an Asian population. Having an ideal CVH status can lower EAA and reduce the risk of aging-related disorders.

Methods

Taiwan Biobank

The TWB has recruited ~ 153,543 community-based volunteers living in Taiwan since October 2012. The written informed consent of each participant was obtained before joining the study. Participants provided blood and urine samples and then underwent physical examinations. Clinical factors of participants were examined and recorded by the TWB. For example, serum glucose was measured with a Hitachi LST008 analyzer (Hitachi High-Technologies, Tokyo, Japan) after a fast for at least 6 h. Participants’ lifestyle factors such as dietary habits, cigarette smoking, physical exercise, educational attainment, and alcohol consumption were collected through a face-to-face interview with TWB researchers.

Because responding to all items in the original TWB questionnaire was time-consuming, a simplified version of the questionnaire was developed to speed up the interview process. Diet-related questions were not included in this simplified version of the questionnaire. As a result, the diet information of 42% of the 2474 individuals was missing because 1039 participants chose the simplified questionnaire.

From 2016 to 2021, the TWB randomly selected 2474 participants for DNAm quantification analysis. The Illumina Infinium MethylationEPIC BeadChip (Illumina, Inc., San Diego, CA) was used to analyze their blood DNAm levels.

The TWB was approved by the Ethics and Governance Council of Taiwan Biobank and the Institutional Review Board on Biomedical Science Research/IRB-BM at Academia Sinica. The TWB approved our access to the data on February 18, 2020, with an application number of TWBR10810-07. Our study was also approved by the Research Ethics Committee of National Taiwan University Hospital (NTUH-REC no. 201805050RINB). Written informed consent was obtained from each individual, following the principles of the Declaration of Helsinki and the institutional requirements.

Epigenetic age acceleration

The quality control and normalization of the DNAm data can be found in our previous works [20, 21]. Briefly speaking, the quality of all 2474 samples was satisfactory according to the quality control process of DNAm data [30]. Our intensity data were normalized with the normal-exponential out-of-band (noob) method [31] implemented by the preprocessNoob function in the R package minfi v1.36 [32].

To calculate epigenetic age, we uploaded the DNAm data quantified through the blood samples of the 2474 TWB participants to the online DNAm Age Calculator developed by Horvath’s laboratory (https://dnamage.genetics.ucla.edu/new). Four epigenetic clocks were applied, including Hannum et al.’s clock [10], Horvath’s clock [11], Levine et al.’s PhenoAge [12], and Lu et al.’s GrimAge [13]. Once the epigenetic age was calculated, four measures of EAA (HannumEAA, IEAA, PhenoEAA, and GrimEAA, respectively) were then obtained by the residuals of regressing the corresponding DNAm age on chronological age.

According to the definition of residuals, EAA is the difference between one’s DNAm age and the predicted DNAm age of individuals at the same chronological age. Positive EAA levels indicate that individuals are aging faster than the expected aging rate at the same chronological age. Contrarily, negative EAA levels imply that individuals are aging slower than the expected aging rate at the same chronological age.

Cardiovascular health score

CVH scores were calculated according to AHA’s “Life’s simple seven” CVH metrics [1]. The seven components include smoking status, physical activity, dietary habits, BMI, total cholesterol, blood pressure, and fasting glucose (Table 3). The CVH score integrates the seven components mentioned above, representing individuals’ health in cardiovascular conditions. Two different scales were used to evaluate the idealness of each CVH component. Each component was scored as 0–1 on a 2-level scale (ideal vs. not ideal) or 0–2 on a 3-level scale (poor, intermediate, or ideal) (Table 3). Therefore, a 7-point CVH score and a 14-point CVH score were obtained by summing up all scores from the seven components.

Table 3 Ideal CVH criterion (the 14-point CVH score)

A higher CVH score reflects a better overall CVH condition, associated with a lower incidence of CVD-related mortality and morbidity [33,34,35]. Generally, an individual’s CVH was categorized as favorable given at least 5 CVH components achieving the ideal criteria (i.e., CVH score ≥ 5 points on the 7-point scale, or CVH score ≥ 10 points on the 14-point scale); moderate given 3–4 components fulfilling the ideal criteria; unfavorable given ≤ 2 components achieving the ideal standards (i.e., CVH score ≤ 2 points on the 7-point scale, or CVH score ≤ 4 points on the 14-point scale) [36,37,38].

Some studies have suggested more stringent BMI cutoff levels for East Asians [39,40,41]. BMI > 24 has been found to be a risk factor for CVD, type 2 diabetes, and metabolic syndromes in Taiwanese populations [42, 43]. Therefore, according to the definition by the Ministry of Health and Welfare of Taiwan, BMI over 24 was defined as overweight, and BMI over 27 was regarded as obesity (Table 3).

To let “smoking status” be consistent with the TWB questionnaire, we adjusted the definition of a former smoker from quitting smoking for at least 12–6 months. Moreover, regarding “ideal diet,” because TWB did not survey the amount of various foods consumed every day, the scoring of ideal diet was modified according to 17 diet-related questions in the TWB questionnaire (Additional file 1: Table S4). Furthermore, 42% of the 2474 participants chose the simplified version of the TWB questionnaire, and therefore the diet-related information was only available for 58% of the individuals in our study. To not compromise the statistical power by only analyzing 58% of the data, we also calculated a 6-component CVH score while removing the “ideal diet” component. Therefore, we had a 6-point scale and a 12-point scale of the CVH score, following the 2-level and 3-level scoring rules, respectively.

Exclusion criteria

Similar to the previous study linking CVH with EAA (focusing on the WHI data) [5], we excluded participants diagnosed with leukemia or CVDs and participants with extreme EAA levels. Including leukemia cases into DNAm data analysis may introduce bias [5] and reduce the precision of DNAm clock models [44]. Therefore, participants with leukemia diagnoses were excluded from our research. However, because no individuals have been diagnosed with leukemia in our study sample, we did not exclude any participants due to this criterion.

Following the World Health Organization (WHO) definition, CVDs included the diagnosis of valvular heart disease, coronary artery disease, arrhythmia, cardiomyopathy, congenital heart disease, apoplexy or any other diseases involving blood vessels or the heart. We removed 276 participants with CVDs because we aimed to test the associations of the CVH scores with EAA among Asian individuals without CVDs.

Moreover, extreme outliers were defined if EAA \(> Q_{3} + 3 \times {\text{IQR}}\) or EAA \(< Q_{1} - 3 \times {\text{IQR}}\), where \(Q_{1}\) and \(Q_{3}\) were the 25th and 75th percentiles and \({\text{IQR}} = Q_{3} - Q_{1}\) (i.e., interquartile range). Among the 2474 TWB participants, 7, 1, 2, and 5 extreme outliers were excluded from our analyses for HannumEAA, IEAA, PhenoEAA, and GrimEAA, respectively.

Statistical analysis

To investigate the association of CVH with EAA, we regressed each of the four measures of EAA (HannumEAA [10], IEAA [11], PhenoEAA [12], and GrimEAA [13]) on the CVH score after excluding 276 participants with CVDs and participants with extreme EAA levels.

For the regression analysis, we considered the following model:

$${\text{EAA}} = \beta_{0} + \beta_{1} {\text{CVH}}_{{{\text{score}}}} + \beta_{2} {\text{SEX}} + \beta_{3} {\text{DRK}} + \beta_{4} {\text{EDU}} + \varepsilon ,$$
(1)

where EAA was served as the response variable, the CVH score was the primary explanatory variable, and \(\varepsilon\) was the random error term. Four different calculations of the CVH scores (7-point, 14-point, 6-point, and 12-point) were included in our analyses to evaluate whether our main conclusions may be sensitive to the CVH scoring scales or the exclusion of diet information. A total of 16 regression models were therefore fitted according to four measures of EAA (IEAA, HannumEAA, PhenoEAA, and GrimEAA) and four definitions of the CVH score (6-, 7-, 12-, and 14-point scales).

Similar to the two previous studies linking CVH with EAA [5, 6], all our models were adjusted for sex (SEX), alcohol drinking status (yes vs. no), and educational attainment (an integer ranging from 1 to 7). Drinking (DRK) was defined as an individual having more than 150 mL intake of alcohol per week for at least 6 months. Educational attainment (EDU) was classified into seven categories (1: Illiterate; 2: No formal education but literate; …; 7: Master’s or higher degree), as listed in Table 1.

We considered EDU a covariate because it was associated with EAA based on various DNAm clocks [45, 46]. Moreover, EDU was also adjusted in the models of the two previous studies linking CVH with EAA [5, 6]. In our analyses, EDU was inversely associated with GrimEAA in the 6-point and 12-point CVH score models (both p = 0.002, Additional file 1: Table S5) and PhenoEAA in all models (p < 0.05, Additional file 1: Table S5).

Variance inflation factor (VIF) values were calculated (Additional file 1: Table S6) to check multicollinearity. All VIF values were controlled under 1.2 (the largest VIF value was 1.1842). No multicollinearity among the explanatory variables was detected in all our models.

Regarding the model evaluation, we performed residual analyses (Additional file 1: Figs. S1–S4) to check the assumptions of normality and constant variance. No substantial violation of these two assumptions was observed in regression models based on IEAA (Additional file 1: Fig. S1), HannumEAA (Additional file 1: Fig. S2), and PhenoEAA (Additional file 1: Fig. S3). However, the quantile–quantile plots of the models based on GrimEAA (Additional file 1: Fig. S4) showed that the residuals followed distributions with heavier tails than the normal distribution.

Deviations from normality usually do not bias the regression coefficients (Beta in Table 2) [47, 48]. To evaluate whether our statistical significance (p-values) may be sensitive to the violation of the normality assumption, we further performed the rank-based inverse normal transformation (rank-based INT) on GrimEAA (Additional file 1: Fig. S8). Through this transformation, no substantial violation of the assumption of normality or constant variance was observed for any model (Additional file 1: Fig. S8). Results of using the rank-based INT on the four measures of EAA as the response variable in model (1) were shown in Additional file 1: Table S7. The statistical significance was similar to our original results (Table 2).

Availability of data and materials

The datasets used and analyzed during the current study are available from https://www.twbiobank.org.tw/.

Abbreviations

CVH:

Cardiovascular health

AHA:

American Heart Association

EAA:

Epigenetic age acceleration

DNAm:

DNA methylation

TWB:

Taiwan Biobank

IEAA:

Intrinsic EAA

EEAA:

Extrinsic EAA

CVH score:

Cardiovascular health score

BMI:

Body mass index

CVD:

Cardiovascular disease

CpG:

Cytosine-phosphate-guanine

WHI:

Women’s Health Initiative

CARDIA:

Coronary Artery Risk Development in Young Adults

FHS:

Framingham Heart Study

IRB-BM:

Institutional Review Board on Biomedical Science Research

NTU:

National Taiwan University

NTUH:

National Taiwan University Hospital

CI:

Confidence interval

SD:

Standard deviation

VIF:

Variance inflation factor

INT:

Inverse normal transformation

WHO:

World Health Organization

References

  1. Lloyd-Jones DM, Hong Y, Labarthe D, Mozaffarian D, Appel LJ, Van Horn L, Greenlund K, Daniels S, Nichol G, Tomaselli GF. Defining and setting national goals for cardiovascular health promotion and disease reduction: the American Heart Association’s strategic Impact Goal through 2020 and beyond. Circulation. 2010;121(4):586–613.

    Article  PubMed  Google Scholar 

  2. Gao B, Wang F, Zhu M, Wang J, Zhou M, Zhang L, Zhao M. Cardiovascular health metrics and all-cause mortality and mortality from major non-communicable chronic diseases among Chinese adult population. Int J Cardiol. 2020;313:123–8.

    Article  PubMed  Google Scholar 

  3. Agostinis-Sobrinho C, Werneck AO, Kievišienė J, Moreira C, Ramírez-Vélez R, Rosário R, Norkiene S, Lopes L, Mota J, Santos R. Ideal cardiovascular health status and health-related quality of life in adolescents: the LABMED physical activity study. Rev Paul Pediatr. 2021;39:e2019343.

    Article  PubMed  Google Scholar 

  4. Wu X, Liu X, Liao W, Kang N, Sang S, Abdulai T, Niu M, He Y, Zhai Z, Pan M, et al. The association between ideal cardiovascular health and health-related quality of life in adults: a population-based cross-sectional study. Int J Public Health. 2021;66:592043.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Pottinger TD, Khan SS, Zheng Y, Zhang W, Tindle HA, Allison M, Wells G, Shadyab AH, Nassir R, Martin LW, et al. Association of cardiovascular health and epigenetic age acceleration. Clin Epigenetics. 2021;13(1):42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Joyce BT, Gao T, Zheng Y, Ma J, Hwang SJ, Liu L, Nannini D, Horvath S, Lu AT, Bai Allen N, et al. Epigenetic age acceleration reflects long-term cardiovascular health. Circ Res. 2021;129(8):770–81.

    Article  CAS  PubMed  Google Scholar 

  7. Jones MJ, Goodman SJ, Kobor MS. DNA methylation and healthy human aging. Aging Cell. 2015;14(6):924–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Saul D, Kosinsky RL. Epigenetics of aging and aging-associated diseases. Int J Mol Sci. 2021;22(1):401.

    Article  CAS  PubMed Central  Google Scholar 

  9. He X, Liu J, Liu B, Shi J. The use of DNA methylation clock in aging research. Exp Biol Med (Maywood). 2021;246(4):436–46.

    Article  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

  11. Horvath S. DNA methylation age of human tissues and cell types. Genome Biol. 2013;14(10):R115.

    Article  PubMed  PubMed Central  Google Scholar 

  12. 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(4):573–91.

    Article  Google Scholar 

  13. 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(2):303–27.

    Article  CAS  Google Scholar 

  14. Zou H, Hastie T. Regularization and variable selection via the elastic net. J R Stat Soc Ser B Stat Methodol. 2005;67:301–20.

    Article  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

  16. Noroozi R, Ghafouri-Fard S, Pisarek A, Rudnicka J, Spólnicka M, Branicki W, Taheri M, Pośpiech E. DNA methylation-based age clocks: From age prediction to age reversion. Ageing Res Rev. 2021;68:101314.

    Article  CAS  PubMed  Google Scholar 

  17. McCrory C, Fiorito G, Hernandez B, Polidoro S, O’Halloran AM, Hever A, Ni Cheallaigh C, Lu AT, Horvath S, Vineis P, et al. GrimAge outperforms other epigenetic clocks in the prediction of age-related clinical phenotypes and all-cause mortality. J Gerontol A Biol Sci Med Sci. 2021;76(5):741–9.

    Article  PubMed  Google Scholar 

  18. White AJ, Kresovich JK, Xu Z, Sandler DP, Taylor JA. Shift work, DNA methylation and epigenetic age. Int J Epidemiol. 2019;48(5):1536–44.

    Article  PubMed  PubMed Central  Google Scholar 

  19. McEwen LM, Jones MJ, Lin DTS, Edgar RD, Husquin LT, MacIsaac JL, Ramadori KE, Morin AM, Rider CF, Carlsten C, et al. Systematic evaluation of DNA methylation age estimation with common preprocessing methods and the Infinium MethylationEPIC BeadChip array. Clin Epigenetics. 2018;10(1):123.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Lin WY, Wang YC, Teng IH, Liu C, Lou XY. Associations of five obesity metrics with epigenetic age acceleration: evidence from 2474 Taiwan Biobank participants. Obesity (Silver Spring). 2021;29(10):1731–8.

    Article  CAS  Google Scholar 

  21. Lin WY. Genome-wide association study for four measures of epigenetic age acceleration and two epigenetic surrogate markers using DNA methylation data from Taiwan biobank. Hum Mol Genet 2021.

  22. Batsis JA, Mackenzie TA, Vasquez E, Germain CM, Emeny RT, Rippberger P, Lopez-Jimenez F, Bartels SJ. Association of adiposity, telomere length and mortality: data from the NHANES 1999–2002. Int J Obes (Lond). 2018;42(2):198–204.

    Article  CAS  Google Scholar 

  23. Koriath M, Müller C, Pfeiffer N, Nickels S, Beutel M, Schmidtmann I, Rapp S, Münzel T, Westermann D, Karakas M, et al. Relative telomere length and cardiovascular risk factors. Biomolecules. 2019;9(5):192.

    Article  CAS  PubMed Central  Google Scholar 

  24. Wu Y, Liu Y, Ni N, Bao B, Zhang C, Lu L. High lead exposure is associated with telomere length shortening in Chinese battery manufacturing plant workers. Occup Environ Med. 2012;69(8):557–63.

    Article  CAS  PubMed  Google Scholar 

  25. Rizvi S, Raza ST, Mahdi F. Telomere length variations in aging and age-related diseases. Curr Aging Sci. 2014;7(3):161–7.

    Article  CAS  PubMed  Google Scholar 

  26. Shammas MA. Telomeres, lifestyle, cancer, and aging. Curr Opin Clin Nutr Metab Care. 2011;14(1):28–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. McCartney DL, Min JL, Richmond RC, Lu AT, Sobczyk MK, Davies G, Broer L, Guo X, Jeong A, Jung J, et al. Genome-wide association studies identify 137 genetic loci for DNA methylation biomarkers of aging. Genome Biol. 2021;22(1):194.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Design of the Women’s Health Initiative clinical trial and observational study. The Women’s Health Initiative Study Group. Control Clin Trials. 1998;19(1):61–109.

  29. Chen CH, Yang JH, Chiang CWK, Hsiung CN, Wu PE, Chang LC, Chu HW, Chang J, Song IW, Yang SL, et al. Population structure of Han Chinese in the modern Taiwanese population based on 10,000 participants in the Taiwan Biobank project. Hum Mol Genet. 2016;25(24):5321–31.

    CAS  PubMed  PubMed Central  Google Scholar 

  30. Maksimovic J, Phipson B, Oshlack A. A cross-package Bioconductor workflow for analysing methylation array data. F1000Res. 2016;5:1281.

    Article  PubMed  CAS  Google Scholar 

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Aryee MJ, Jaffe AE, Corrada-Bravo H, Ladd-Acosta C, Feinberg AP, Hansen KD, Irizarry RA. Minfi: a flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinformatics. 2014;30(10):1363–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Michos ED, Khan SS. Further understanding of ideal cardiovascular health score metrics and cardiovascular disease. Expert Rev Cardiovasc Ther. 2021;19(7):607–17.

    Article  CAS  PubMed  Google Scholar 

  34. Ramírez-Vélez R, Saavedra JM, Lobelo F, Celis-Morales CA, Pozo-Cruz BD, García-Hermoso A. Ideal cardiovascular health and incident cardiovascular disease among adults: a systematic review and meta-analysis. Mayo Clin Proc. 2018;93(11):1589–99.

    Article  PubMed  Google Scholar 

  35. Guo L, Zhang S. Association between ideal cardiovascular health metrics and risk of cardiovascular events or mortality: a meta-analysis of prospective studies. Clin Cardiol. 2017;40(12):1339–46.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Osondu CU, Aneni EC, Valero-Elizondo J, Salami JA, Rouseff M, Das S, Guzman H, Younus A, Ogunmoroti O, Feldman T, et al. Favorable cardiovascular health is associated with lower health care expenditures and resource utilization in a large US employee population: The Baptist Health South Florida Employee Study. Mayo Clin Proc. 2017.

  37. Roger VL, Sidney S, Fairchild AL, Howard VJ, Labarthe DR, Shay CM, Tiner AC, Whitsel LP, Rosamond WD. Recommendations for cardiovascular health and disease surveillance for 2030 and beyond: a policy statement from the American Heart Association. Circulation. 2020;141(9):e104–19.

    Article  PubMed  Google Scholar 

  38. Benjamin EJ, Muntner P, Alonso A, Bittencourt MS, Callaway CW, Carson AP, Chamberlain AM, Chang AR, Cheng S, Das SR, et al. Heart disease and stroke statistics-2019 update: a report from the American Heart Association. Circulation. 2019;139(10):e56–528.

    Article  PubMed  Google Scholar 

  39. Lin WY, Lee LT, Chen CY, Lo H, Hsia HH, Liu IL, Lin RS, Shau WY, Huang KC. Optimal cut-off values for obesity: using simple anthropometric indices to predict cardiovascular risk factors in Taiwan. Int J Obes Relat Metab Disord. 2002;26(9):1232–8.

    Article  PubMed  Google Scholar 

  40. Xu F, Wang YF, Lu L, Liang Y, Wang Z, Hong X, Li J. Comparison of anthropometric indices of obesity in predicting subsequent risk of hyperglycemia among Chinese men and women in Mainland China. Asia Pac J Clin Nutr. 2010;19(4):586–93.

    PubMed  Google Scholar 

  41. Zeng Q, He Y, Dong S, Zhao X, Chen Z, Song Z, Chang G, Yang F, Wang Y. Optimal cut-off values of BMI, waist circumference and waist:height ratio for defining obesity in Chinese adults. Br J Nutr. 2014;112(10):1735–44.

    Article  CAS  PubMed  Google Scholar 

  42. Hwang LC, Bai CH, Sun CA, Chen CJ. Prevalence of metabolically healthy obesity and its impacts on incidences of hypertension, diabetes and the metabolic syndrome in Taiwan. Asia Pac J Clin Nutr. 2012;21(2):227–33.

    CAS  PubMed  Google Scholar 

  43. Hwang LC, Bai CH, Chen CJ. Prevalence of obesity and metabolic syndrome in Taiwan. J Formos Med Assoc. 2006;105(8):626–35.

    Article  PubMed  Google Scholar 

  44. Finesso GE, McDevitt RA, Roy R, Brinster LR, Di Francesco A, Meade T, de Cabo R, Ferrucci L, Perdue KA. Impact of large granular lymphocyte leukemia on blood DNA methylation and epigenetic clock modeling in Fischer 344 rats. J Gerontol A Biol Sci Med Sci. 2021;77:956–63.

    Article  Google Scholar 

  45. 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(2):419–46.

    Article  CAS  Google Scholar 

  46. Zhao W, Ammous F, Ratliff S, Liu J, Yu M, Mosley TH, Kardia SLR, Smith JA. Education and lifestyle factors are associated with DNA methylation clocks in older African Americans. Int J Environ Res Public Health. 2019;16(17):3141.

    Article  CAS  PubMed Central  Google Scholar 

  47. Ramsey F, Schafer D. The statistical sleuth: a course in methods of data analysis. Brooks/Cole 2013.

  48. Williams MN, Grajales CAG, Kurkiewicz D. Assumptions of multiple regression: correcting two misconceptions. Pract Assess Res Eval. 2013;18:Article 11.

Download references

Acknowledgments

The authors would like to thank three anonymous reviewers for their insightful and constructive comments, and the Taiwan Biobank for approving our application to access the data.

Funding

This study was supported by the Ministry of Science and Technology of Taiwan (Grant Numbers MOST 107-2314-B-002-195-MY3 and 110-2314-B-002-071 to Wan-Yu Lin).

Author information

Authors and Affiliations

Authors

Contributions

This study was conceived, designed, and interpreted by YHL and WYL. YHL performed statistical analyses and drafted the manuscript. WYL supervised this study and substantively revised the manuscript. Both authors read and approved the final manuscript.

Corresponding author

Correspondence to Wan-Yu Lin.

Ethics declarations

Ethics approval and consent to participate

TWB was approved by the Institutional Review Board on Biomedical Science Research/IRB-BM, Academia Sinica, and also by the Ethics and Governance Council of Taiwan Biobank, Taiwan. Written informed consent was obtained from each participant in accordance with institutional requirements and the principles of the Declaration of Helsinki. Our use of the TWB research data was approved by TWB on February 18, 2020 (application number: TWBR10810-07). This study further received approval from the Research Ethics Committee of National Taiwan University Hospital (NTUH-REC no. 201805050RINB). Informed consent was obtained from all individual participants included in this study.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

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

Supplementary Information

Additional file 1. Table S1.

Basic characteristics of the 2474 TWB participants stratified by tertiles of PhenoEAA. Table S2. Basic characteristics of the 2474 TWB participants stratified by tertiles of IEAA. Table S3. Basic characteristics of the 2474 TWB participants stratified by tertiles of HannumEAA. Table S4. The 17 diet-related questions in the TWB questionnaire. Table S5. Regression coefficients of all the covariates included in statistical models. Table S6. Variance inflation factors (VIF) to check multicollinearity. Table S7. Regressing rank-based inverse normal transformation of the four measures of EAA on the CVH score. Figure S1. Residual and Normal Quantile-Quantile plots for regression models based on IEAA. Figure S2. Residual and Normal Quantile-Quantile plots for regression models based on HannumEAA. Figure S3. Residual and Normal Quantile-Quantile plots for regression models based on PhenoEAA. Figure S4. Residual and Normal Quantile-Quantile plots for regression models based on GrimEAA. Figure S5. Residual and Normal Quantile-Quantile plots for rank-based inverse normal transformation (rank-based INT) of IEAA. Figure S6. Residual and Normal Quantile-Quantile plots for rank-based inverse normal transformation (rank-based INT) of HannumEAA. Figure S7. Residual and Normal Quantile-Quantile plots for rank-based inverse normal transformation (rank-based INT) of PhenoEAA. Figure S8. Residual and Normal Quantile-Quantile plots for rank-based inverse normal transformation (rank-based INT) of GrimEAA.

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Lo, YH., Lin, WY. Cardiovascular health and four epigenetic clocks. Clin Epigenet 14, 73 (2022). https://doi.org/10.1186/s13148-022-01295-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13148-022-01295-7

Keywords