Blood DNA methylation sites predict death risk in a longitudinal study of 12, 300 individuals

DNA methylation has fundamental roles in gene programming and aging that may help predict mortality. However, no large-scale study has investigated whether site-specific DNA methylation predicts all-cause mortality. We used the Illumina-HumanMethylation450-BeadChip to identify blood DNA methylation sites associated with all-cause mortality for 12, 300 participants in 12 Cohorts of the Heart and Aging Research in Genetic Epidemiology (CHARGE) Consortium. Over an average 10-year follow-up, there were 2,561 deaths across the cohorts. Nine sites mapping to three intergenic and six gene-specific regions were associated with mortality (P < 9.3x10^-7) independently of age and other mortality predictors. Six sites (cg14866069, cg23666362, cg20045320, cg07839457, cg07677157, cg09615688)—mapping respectively to BMPR1B, MIR1973, IFITM3, NLRC5, and two intergenic regions—were associated with reduced mortality risk. The remaining three sites (cg17086398, cg12619262, cg18424841)—mapping respectively to SERINC2, CHST12, and an intergenic region—were associated with increased mortality risk. DNA methylation at each site predicted 5%–15% of all deaths. We also assessed the causal association of those sites to age-related chronic diseases by using Mendelian randomization, identifying weak causal relationship between cg18424841 and cg09615688 with coronary heart disease. Of the nine sites, three (cg20045320, cg07839457, cg07677157) were associated with lower incidence of heart disease risk and two (cg20045320, cg07839457) with smoking and inflammation in prior CHARGE analyses. Methylation of cg20045320, cg07839457, and cg17086398 was associated with decreased expression of nearby genes (IFITM3, IRF, NLRC5, MT1, MT2, MARCKSL1) linked to immune responses and cardiometabolic diseases. These sites may serve as useful clinical tools for mortality risk assessment and preventative care.

Genome-wide meta-analysis of variant-by-diuretic interactions as modulators of lipid traits in persons of European and African ancestry

Hypertension (HTN) is a significant risk factor for cardiovascular morbidity and mortality. Metabolic abnormalities, including adverse cholesterol and triglycerides (TG) profiles, are frequent comorbid findings with HTN and contribute to cardiovascular disease. Diuretics, which are used to treat HTN and heart failure, have been associated with worsening of fasting lipid concentrations. Genome-wide meta-analyses with 39,710 European-ancestry (EA) individuals and 9925 African-ancestry (AA) individuals were performed to identify genetic variants that modify the effect of loop or thiazide diuretic use on blood lipid concentrations. Both longitudinal and cross sectional data were used to compute cohort-specific interaction results, which were then combined through meta-analysis in each ancestry. These ancestry-specific results were further combined through trans-ancestry meta-analysis. Analysis of EA data identified two genome-wide significant (p < 5 × 10-8) loci with single nucleotide variant (SNV)-loop diuretic interaction on TG concentrations (including COL11A1). Analysis of AA data identified one genome-wide significant locus adjacent to BMP2 with SNV-loop diuretic interaction on TG concentrations. Trans-ancestry analysis strengthened evidence of association for SNV-loop diuretic interaction at two loci (KIAA1217 and BAALC). There were few significant SNV-thiazide diuretic interaction associations on TG concentrations and for either diuretic on cholesterol concentrations. Several promising loci were identified that may implicate biologic pathways that contribute to adverse metabolic side effects from diuretic therapy.

Multi-ancestry GWAS of the electrocardiographic PR interval identifies 202 loci underlying cardiac conduction

The electrocardiographic PR interval reflects atrioventricular conduction, and is associated with conduction abnormalities, pacemaker implantation, atrial fibrillation (AF), and cardiovascular mortality. Here we report a multi-ancestry (N = 293,051) genome-wide association meta-analysis for the PR interval, discovering 202 loci of which 141 have not previously been reported. Variants at identified loci increase the percentage of heritability explained, from 33.5% to 62.6%. We observe enrichment for cardiac muscle developmental/contractile and cytoskeletal genes, highlighting key regulation processes for atrioventricular conduction. Additionally, 8 loci not previously reported harbor genes underlying inherited arrhythmic syndromes and/or cardiomyopathies suggesting a role for these genes in cardiovascular pathology in the general population. We show that polygenic predisposition to PR interval duration is an endophenotype for cardiovascular disease, including distal conduction disease, AF, and atrioventricular pre-excitation. These findings advance our understanding of the polygenic basis of cardiac conduction, and the genetic relationship between PR interval duration and cardiovascular disease.

Insulin resistance and reduced cardiac autonomic function in older adults: the Atherosclerosis Risk in Communities study

Background

Prior studies have shown insulin resistance is associated with reduced cardiac autonomic function measured at rest, but few studies have determined whether insulin resistance is associated with reduced cardiac autonomic function measured during daily activities.

Methods

We examined older adults without diabetes with 48-h ambulatory electrocardiography (n = 759) in an ancillary study of the Atherosclerosis Risk in Communities Study. Insulin resistance, the exposure, was defined by quartiles for three indexes: 1) the homeostatic model assessment of insulin resistance (HOMA-IR), 2) the triglyceride and glucose index (TyG), and 3) the triglyceride to high-density lipoprotein cholesterol ratio (TG/HDL-C). Low heart rate variability, the outcome, was defined by <25th percentile for four measures: 1) standard deviation of normal-to-normal R-R intervals (SDNN), a measure of total variability; 2) root mean square of successive differences in normal-to-normal R-R intervals (RMSSD), a measure of vagal activity; 3) low frequency spectral component (LF), a measure of sympathetic and vagal activity; and 4) high frequency spectral component (HF), a measure of vagal activity. Logistic regression was used to estimate odds ratios (OR) and 95% confidence intervals weighted for sampling/non-response, adjusted for age at ancillary visit, sex, and race/study-site. Insulin resistance quartiles 4, 3, and 2 were compared to quartile 1; high indexes refer to quartile 4 versus quartile 1.

Results

The average age was 78 years, 66% (n = 497) were women, and 58% (n = 438) were African American. Estimates of association were not robust at all levels of HOMA-IR, TyG, and TG/HDL-C, but suggest that high indexes were associated consistently with indicators of vagal activity. High HOMA-IR, high TyG, and high TG/HDL-C were consistently associated with low RMSSD (OR: 1.68 (1.00, 2.81), OR: 2.03 (1.21, 3.39), and OR: 1.73 (1.01, 2.91), respectively). High HOMA-IR, high TyG, and high TG/HDL-C were consistently associated with low HF (OR: 1.90 (1.14, 3.18), OR: 1.98 (1.21, 3.25), and OR: 1.76 (1.07, 2.90), respectively).

Conclusions

In older adults without diabetes, insulin resistance was associated with reduced cardiac autonomic function – specifically and consistently for indicators of vagal activity – measured during daily activities. Primary prevention of insulin resistance may reduce the related risk of cardiac autonomic dysfunction.

DNA methylation GrimAge strongly predicts lifespan and healthspan

It was unknown whether plasma protein levels can be estimated based on DNA methylation (DNAm) levels, and if so, how the resulting surrogates can be consolidated into a powerful predictor of lifespan. We present here, seven DNAm-based estimators of plasma proteins including those of plasminogen activator inhibitor 1 (PAI-1) and growth differentiation factor 15. The resulting predictor of lifespan, DNAm GrimAge (in units of years), is a composite biomarker based on the seven DNAm surrogates and a DNAm-based estimator of smoking pack-years. Adjusting DNAm GrimAge for chronological age generated novel measure of epigenetic age acceleration, AgeAccelGrim.Using large scale validation data from thousands of individuals, we demonstrate that DNAm GrimAge stands out among existing epigenetic clocks in terms of its predictive ability for time-to-death (Cox regression P=2.0E-75), time-to-coronary heart disease (Cox P=6.2E-24), time-to-cancer (P= 1.3E-12), its strong relationship with computed tomography data for fatty liver/excess visceral fat, and age-at-menopause (P=1.6E-12). AgeAccelGrim is strongly associated with a host of age-related conditions including comorbidity count (P=3.45E-17). Similarly, age-adjusted DNAm PAI-1 levels are associated with lifespan (P=5.4E-28), comorbidity count (P= 7.3E-56) and type 2 diabetes (P=2.0E-26). These DNAm-based biomarkers show the expected relationship with lifestyle factors including healthy diet and educational attainment.Overall, these epigenetic biomarkers are expected to find many applications including human anti-aging studies.

Particulate Matter and Albuminuria, Glomerular Filtration Rate, and Incident CKD

BACKGROUND AND OBJECTIVES

Exposure to particulate matter (PM) <2.5 μm in aerodynamic diameter (PM2.5) has been linked to detrimental health effects. This study aimed to describe the relationship between long-term PM2.5 exposure and kidney disease, including eGFR, level of albuminuria, and incident CKD.

DESIGN, SETTING, PARTICIPANTS, & MEASUREMENTS

The study included 10,997 participants from the Atherosclerosis Risk in Communities cohort who were followed from 1996-1998 through 2016. Monthly mean PM2.5 concentrations (μg/m3) were estimated at geocoded participant addresses using geographic information system-based, spatiotemporal generalized additive mixed models-including geospatial covariates such as land use-and then averaged over the 12-month period preceding participant examination. Covariate-adjusted, cross-sectional associations of PM2.5, baseline eGFR, and urinary albumin-creatinine ratio (UACR) were estimated using linear regression. PM2.5 and incident CKD (defined as follow-up eGFR <60 ml/min per 1.73 m2 with ≥25% eGFR decline relative to baseline, CKD-related hospitalization or death based on International Classification of Diseases 9/10 codes, or development of ESKD) associations were estimated using Cox proportional hazards regression. Modeling was stratified by study site, and stratum-specific estimates were combined using random-effects meta-analyses.

RESULTS

Baseline mean participant age was 63 (±6) years and eGFR was 86 (±16) ml/min per 1.73 m2. There was no significant PM2.5-eGFR association at baseline. Each 1-μg/m3 higher annual average PM2.5 was associated with higher UACR after adjusting for demographics, socioeconomic status, and clinical covariates (percentage difference, 6.6%; 95% confidence interval [95% CI], 2.6% to 10.7%). Each 1-μg/m3 higher annual average PM2.5 was associated with a significantly higher risk of incident CKD (hazard ratio, 1.05; 95% CI, 1.01 to 1.10).

CONCLUSIONS

Exposure to higher annual average PM2.5 concentrations was associated with a higher level of albuminuria and higher risk for incident CKD in a community-based cohort.

Race/Ethnicity, Gender, and Trajectories of Depressive Symptoms Across Early- and Mid-Life Among the Add Health Cohort

While disparities in depressive symptoms by race/ethnicity and gender have been documented, left unclear is how such status characteristics intersect to influence mental health, particularly across early life and among a diverse set of population subgroups. This study investigates how intra- and inter-individual trends in depressive symptoms unfold across a 30-year span (ages 12-42) and are structured by the intersection of race/ethnicity and gender among White, Black, Hispanic, and Asian American young adults (N = 18,566). Analyses use data from the National Longitudinal Study of Adolescent to Adult Health, a nationally representative study of adolescents who have been followed through their fourth decade of life. We draw on Waves I-IV and a representative subsample of the brand new Wave V data. Growth curve models indicated depressive symptoms decreased across adolescence and young adulthood before increasing in the early 30s. Racial/ethnic minorities reported more depressive symptoms than Whites. Women reported more depressive symptoms than men and experienced especially steep increases in their late 30s. Racial/ethnic-gender disparities remained stable with age, except for Hispanic-White disparities among women and Asian American-White disparities among men, which narrowed with age. Overall, findings demonstrate dynamic inequalities across a longer period of the life span than was previously known, as well as heterogeneity in trajectories of poor mental health within and between racial/ethnic-gender groups. Results also suggest that Black and Asian American women experience the highest mental health risks and that interventions for reducing disparities in depressive symptoms should focus on adults in their late 20s/early 30s, particularly women of color.

Leukocyte Traits and Exposure to Ambient Particulate Matter Air Pollution in the Women's Health Initiative and Atherosclerosis Risk in Communities Study

BACKGROUND

Inflammatory effects of ambient particulate matter (PM) air pollution exposures may underlie PM-related increases in cardiovascular disease risk and mortality, although evidence of PM-associated leukocytosis is inconsistent and largely based on small, cross-sectional, and/or unrepresentative study populations.

OBJECTIVES

Our objective was to estimate PM-leukocyte associations among U.S. women and men in the Women's Health Initiative and Atherosclerosis Risk in Communities study (n = 165,675 ).

METHODS

We based the PM-leukocyte estimations on up to four study visits per participant, at which peripheral blood leukocytes and geocoded address-specific concentrations of PM ≤ 10 , ≤ 2.5 , and 2.5 - 10 μ m in diameter (PM 10 , PM 2.5 , and PM 2.5 - 10 , respectively) were available. We multiply imputed missing data using chained equations and estimated PM-leukocyte count associations over daily to yearly PM exposure averaging periods using center-specific, linear, mixed, longitudinal models weighted for attrition and adjusted for sociodemographic, behavioral, meteorological, and geographic covariates. In a subset of participants with available data (n = 8,457 ), we also estimated PM-leukocyte proportion associations in compositional data analyses.

RESULTS

We found a 12   cells / μ L (95% confidence interval: - 9 , 33) higher leukocyte count, a 1.2% (0.6%, 1.8%) higher granulocyte proportion, and a - 1.1 % (- 1.9 % , - 0.3 % ) lower CD 8 + T-cell proportion per 10 - μ g / m 3 increase in 1-month mean PM 2.5 . However, shorter-duration PM 10 exposures were inversely and only modestly associated with leukocyte count.

DISCUSSION

The PM 2.5 -leukocyte estimates, albeit imprecise, suggest that among racially, ethnically, and environmentally diverse U.S. populations, sustained, ambient exposure to fine PM may induce subclinical, but epidemiologically important, inflammatory effects. https://doi.org/10.1289/EHP5360.

The Depths of Despair Among US Adults Entering Midlife

OBJECTIVES

To test whether indicators of despair are rising among US adults as they age toward midlife and whether this rise is concentrated among low-educated Whites and in rural areas.

METHODS

We used data from the National Longitudinal Study of Adolescent to Adult Health, a nationally representative study of US adolescents in 1994. Our sample was restricted to individuals who participated in 1 or more of 5 waves (1994-2017) and self-identified as non-Hispanic White, non-Hispanic Black, or Hispanic (n = 18 446). We examined change in indicators of despair from adolescence to adulthood using multilevel regression analysis, testing for differences by race/ethnicity, education, and rurality.

RESULTS

We found evidence of rising despair among this cohort over the past decade. This increase was not restricted to low-educated Whites or to rural areas.

CONCLUSIONS

Results suggest that generally rising despair among the young adult cohort now reaching midlife that cuts across racial/ethnic, educational, and geographic groups may presage rising midlife mortality for these subgroups in the next decade.

Genome-Wide Association Study of Apparent Treatment-Resistant Hypertension in the CHARGE Consortium: The CHARGE Pharmacogenetics Working Group

BACKGROUND

Only a handful of genetic discovery efforts in apparent treatment-resistant hypertension (aTRH) have been described.

METHODS

We conducted a case-control genome-wide association study of aTRH among persons treated for hypertension, using data from 10 cohorts of European ancestry (EA) and 5 cohorts of African ancestry (AA). Cases were treated with 3 different antihypertensive medication classes and had blood pressure (BP) above goal (systolic BP ≥ 140 mm Hg and/or diastolic BP ≥ 90 mm Hg) or 4 or more medication classes regardless of BP control (nEA = 931, nAA = 228). Both a normotensive control group and a treatment-responsive control group were considered in separate analyses. Normotensive controls were untreated (nEA = 14,210, nAA = 2,480) and had systolic BP/diastolic BP < 140/90 mm Hg. Treatment-responsive controls (nEA = 5,266, nAA = 1,817) had BP at goal (<140/90 mm Hg), while treated with one antihypertensive medication class. Individual cohorts used logistic regression with adjustment for age, sex, study site, and principal components for ancestry to examine the association of single-nucleotide polymorphisms with case-control status. Inverse variance-weighted fixed-effects meta-analyses were carried out using METAL.

RESULTS

The known hypertension locus, CASZ1, was a top finding among EAs (P = 1.1 × 10-8) and in the race-combined analysis (P = 1.5 × 10-9) using the normotensive control group (rs12046278, odds ratio = 0.71 (95% confidence interval: 0.6-0.8)). Single-nucleotide polymorphisms in this locus were robustly replicated in the Million Veterans Program (MVP) study in consideration of a treatment-responsive control group. There were no statistically significant findings for the discovery analyses including treatment-responsive controls.

CONCLUSION

This genomic discovery effort for aTRH identified CASZ1 as an aTRH risk locus.

Methylome-wide association study provides evidence of particulate matter air pollution-associated DNA methylation

BACKGROUND

DNA methylation (DNAm) may contribute to processes that underlie associations between air pollution and poor health. Therefore, our objective was to evaluate associations between DNAm and ambient concentrations of particulate matter (PM) ≤2.5, ≤10, and 2.5-10 μm in diameter (PM2.5; PM10; PM2.5-10).

METHODS

We conducted a methylome-wide association study among twelve cohort- and race/ethnicity-stratified subpopulations from the Women's Health Initiative and the Atherosclerosis Risk in Communities study (n = 8397; mean age: 61.5 years; 83% female; 45% African American; 9% Hispanic/Latino American). We averaged geocoded address-specific estimates of daily and monthly mean PM concentrations over 2, 7, 28, and 365 days and 1 and 12 months before exams at which we measured leukocyte DNAm in whole blood. We estimated subpopulation-specific, DNAm-PM associations at approximately 485,000 Cytosine-phosphate-Guanine (CpG) sites in multi-level, linear, mixed-effects models. We combined subpopulation- and site-specific estimates in fixed-effects, inverse variance-weighted meta-analyses, then for associations that exceeded methylome-wide significance and were not heterogeneous across subpopulations (P < 1.0 × 10-7; PCochran's Q > 0.10), we characterized associations using publicly accessible genomic databases and attempted replication in the Cooperative Health Research in the Region of Augsburg (KORA) study.

RESULTS

Analyses identified significant DNAm-PM associations at three CpG sites. Twenty-eight-day mean PM10 was positively associated with DNAm at cg19004594 (chromosome 20; MATN4; P = 3.33 × 10-8). One-month mean PM10 and PM2.5-10 were positively associated with DNAm at cg24102420 (chromosome 10; ARPP21; P = 5.84 × 10-8) and inversely associated with DNAm at cg12124767 (chromosome 7; CFTR; P = 9.86 × 10-8). The PM-sensitive CpG sites mapped to neurological, pulmonary, endocrine, and cardiovascular disease-related genes, but DNAm at those sites was not associated with gene expression in blood cells and did not replicate in KORA.

CONCLUSIONS

Ambient PM concentrations were associated with DNAm at genomic regions potentially related to poor health among racially, ethnically and environmentally diverse populations of U.S. women and men. Further investigation is warranted to uncover mechanisms through which PM-induced epigenomic changes may cause disease.

Short-term exposure to air pollution and incidence of stroke in the Women's Health Initiative

BACKGROUND

Evidence of the association between daily variation in air pollution and risk of stroke is inconsistent, potentially due to the heterogeneity in stroke etiology.

OBJECTIVES

To estimate the associations between daily variation in ambient air pollution and risk of stroke and its subtypes among participants of the Women's Health Initiative, a large prospective cohort study in the United States.

METHODS

We used national-scale, log-normal ordinary kriging models to estimate daily concentrations of fine particulate matter (PM_2.5), respirable particulate matter (PM₁₀), nitrogen dioxide (NO₂), nitrogen oxides (NOx), sulphur dioxide, and ozone at participant addresses. Stroke was adjudicated by trained neurologists and classified as ischemic or hemorrhagic. Ischemic strokes were further classified according to the Trial of ORG 10172 in Acute Stroke Treatment (TOAST) classification. We used a time-stratified case-crossover approach to estimate the odds ratio (OR) of the risk of stroke associated with an interquartile range (IQR) increase in concentrations of each air pollutant. We performed stratified analysis to examine whether associations varied across subgroups defined by age at stroke onset, US census region, smoking status, body mass index, and prior history of diabetes mellitus, hypertension, heart or circulation problems, or arterial fibrillation at enrollment.

RESULTS

Among 5417 confirmed strokes between 1993 and 2012, 4300 (79.4%) were classified as ischemic and 924 (17.1%) as hemorrhagic. No association was observed between day-to-day variation in any pollutant and risk of total stroke, ischemic stroke, or specific etiologies of ischemic stroke. We observed a positive association between risk of hemorrhagic stroke and NO₂ and NO_x in the 3 days prior to stroke with OR of 1.24 (95% CI: 1.01, 1.52) and 1.18 (95% CI: 1.03, 1.34) per IQR increase, respectively. The observed associations with hemorrhagic stroke were more pronounced among non-obese participants.

CONCLUSIONS

In this large cohort of post-menopausal US women, daily NO₂ and NOx were associated with higher risk of hemorrhagic stroke, but ambient levels of four other air pollutants were not associated with higher risk of total stroke, ischemic stroke, or ischemic stroke subtypes.

The prevalence of atrial fibrillation on 48-hour ambulatory electrocardiography in African Americans compared to Whites: The Atherosclerosis Risk in Communities (ARIC) study

BACKGROUND

A lower prevalence of atrial fibrillation (AF), but paradoxically higher burden of cardiovascular disease risk factors, has been observed among African Americans compared to Whites in studies of AF identified by mostly 12-lead electrocardiograms (ECGs) and clinically.

METHODS

We performed 48-hour ambulatory electrocardiography (aECG) in a biracial sample of 1,193 participants in the Atherosclerosis Risk in Communities (ARIC) (mean age = 78 years, 62% African Americans, 64% female). Atrial fibrillation was identified from aECG, study visit ECGs, and discharge codes from cohort hospitalizations. We used covariate-adjusted logistic regression to estimate prevalence odds ratios (ORs) for AF in African Americans versus Whites, with adjustment for sampling and nonresponse.

RESULTS

African Americans were more likely than Whites to have hypertension and diabetes but less likely to have coronary heart disease. The prevalence of AF detected by aECG or ARIC study ECG (adjusted for age and coronary heart disease) was lower in African Americans than Whites (2.7% vs 5.0%). White men had a higher (although not significant) AF prevalence of 7.8% compared to the other race and gender groups at 2.3%-2.8%. The adjusted OR for AF was 0.49 (0.24-0.99) comparing African Americans to Whites. Findings were similar when AF was defined to include prior AF hospitalizations (OR = 0.42, 0.25-0.72). There were no significant differences by race for asymptomatic or paroxysmal AF.

CONCLUSIONS

Atrial fibrillation was less prevalent in African American than white older adults, regardless of detection method. Although overall detection of new AF cases with aECG was low, future studies should consider longer-term monitoring to characterize AF by race.

Epigenetic clock for skin and blood cells applied to Hutchinson Gilford Progeria Syndrome and ex vivo studies

DNA methylation (DNAm)-based biomarkers of aging have been developed for many tissues and organs. However, these biomarkers have sub-optimal accuracy in fibroblasts and other cell types used in ex vivo studies. To address this challenge, we developed a novel and highly robust DNAm age estimator (based on 391 CpGs) for human fibroblasts, keratinocytes, buccal cells, endothelial cells, lymphoblastoid cells, skin, blood, and saliva samples. High age correlations can also be observed in sorted neurons, glia, brain, liver, and even bone samples. Gestational age correlates with DNAm age in cord blood. When used on fibroblasts from Hutchinson Gilford Progeria Syndrome patients, this age estimator (referred to as the skin & blood clock) uncovered an epigenetic age acceleration with a magnitude that is below the sensitivity levels of other DNAm-based biomarkers. Furthermore, this highly sensitive age estimator accurately tracked the dynamic aging of cells cultured ex vivo and revealed that their proliferation is accompanied by a steady increase in epigenetic age. The skin & blood clock predicts lifespan and it relates to many age-related conditions. Overall, this biomarker is expected to become useful for forensic applications (e.g. blood or buccal swabs) and for a quantitative ex vivo human cell aging assay.

Age-related DNA hydroxymethylation is enriched for gene expression and immune system processes in human peripheral blood

DNA methylation (DNAm) has a well-established association with age in many tissues, including peripheral blood mononuclear cells (PBMCs). Compared to DNAm, the closely related epigenetic modification known as DNA hydroxymethylation (DNAhm) was much more recently discovered in mammals. Preliminary investigations have observed a positive correlation between gene body DNAhm and cis-gene expression. While some of these studies have observed an association between age and global DNAhm, none have investigated region-specific age-related DNAhm in human blood samples. In this study, we investigated DNAhm and gene expression in PBMCs of 10 young and 10 old, healthy female volunteers. Thousands of regions were differentially hydroxymethylated in the old vs. young individuals in gene bodies, exonic regions, enhancers, and promoters. Consistent with previous work, we observed directional consistency between age-related differences in DNAhm and gene expression. Further, age-related DNAhm and genes with high levels of DNAhm were enriched for immune system processes which may support a role of age-related DNAhm in immunosenescence.

Blood Leukocyte DNA Methylation Predicts Risk of Future Myocardial Infarction and Coronary Heart Disease

BACKGROUND

DNA methylation is implicated in coronary heart disease (CHD), but current evidence is based on small, cross-sectional studies. We examined blood DNA methylation in relation to incident CHD across multiple prospective cohorts.

METHODS

Nine population-based cohorts from the United States and Europe profiled epigenome-wide blood leukocyte DNA methylation using the Illumina Infinium 450k microarray, and prospectively ascertained CHD events including coronary insufficiency/unstable angina, recognized myocardial infarction, coronary revascularization, and coronary death. Cohorts conducted race-specific analyses adjusted for age, sex, smoking, education, body mass index, blood cell type proportions, and technical variables. We conducted fixed-effect meta-analyses across cohorts.

RESULTS

Among 11 461 individuals (mean age 64 years, 67% women, 35% African American) free of CHD at baseline, 1895 developed CHD during a mean follow-up of 11.2 years. Methylation levels at 52 CpG (cytosine-phosphate-guanine) sites were associated with incident CHD or myocardial infarction (false discovery rate<0.05). These CpGs map to genes with key roles in calcium regulation (ATP2B2, CASR, GUCA1B, HPCAL1), and genes identified in genome- and epigenome-wide studies of serum calcium (CASR), serum calcium-related risk of CHD (CASR), coronary artery calcified plaque (PTPRN2), and kidney function (CDH23, HPCAL1), among others. Mendelian randomization analyses supported a causal effect of DNA methylation on incident CHD; these CpGs map to active regulatory regions proximal to long non-coding RNA transcripts.

CONCLUSION

Methylation of blood-derived DNA is associated with risk of future CHD across diverse populations and may serve as an informative tool for gaining further insight on the development of CHD.

DNA methylation-based estimator of telomere length

Telomere length (TL) is associated with several aging-related diseases. Here, we present a DNA methylation estimator of TL (DNAmTL) based on 140 CpGs. Leukocyte DNAmTL is applicable across the entire age spectrum and is more strongly associated with age than measured leukocyte TL (LTL) (r ~-0.75 for DNAmTL versus r ~ -0.35 for LTL). Leukocyte DNAmTL outperforms LTL in predicting: i) time-to-death (p=2.5E-20), ii) time-to-coronary heart disease (p=6.6E-5), iii) time-to-congestive heart failure (p=3.5E-6), and iv) association with smoking history (p=1.21E-17). These associations are further validated in large scale methylation data (n=10k samples) from the Framingham Heart Study, Women's Health Initiative, Jackson Heart Study, InChianti, Lothian Birth Cohorts, Twins UK, and Bogalusa Heart Study. Leukocyte DNAmTL is also associated with measures of physical fitness/functioning (p=0.029), age-at-menopause (p=0.039), dietary variables (omega 3, fish, vegetable intake), educational attainment (p=3.3E-8) and income (p=3.1E-5). Experiments in cultured somatic cells show that DNAmTL dynamics reflect in part cell replication rather than TL per se. DNAmTL is not only an epigenetic biomarker of replicative history of cells, but a useful marker of age-related pathologies that are associated with it.

GWAS of QRS duration identifies new loci specific to Hispanic/Latino populations

BACKGROUND

The electrocardiographically quantified QRS duration measures ventricular depolarization and conduction. QRS prolongation has been associated with poor heart failure prognosis and cardiovascular mortality, including sudden death. While previous genome-wide association studies (GWAS) have identified 32 QRS SNPs across 26 loci among European, African, and Asian-descent populations, the genetics of QRS among Hispanics/Latinos has not been previously explored.

METHODS

We performed a GWAS of QRS duration among Hispanic/Latino ancestry populations (n = 15,124) from four studies using 1000 Genomes imputed genotype data (adjusted for age, sex, global ancestry, clinical and study-specific covariates). Study-specific results were combined using fixed-effects, inverse variance-weighted meta-analysis.

RESULTS

We identified six loci associated with QRS (P<5x10-8), including two novel loci: MYOCD, a nuclear protein expressed in the heart, and SYT1, an integral membrane protein. The top SNP in the MYOCD locus, intronic SNP rs16946539, was found in Hispanics/Latinos with a minor allele frequency (MAF) of 0.04, but is monomorphic in European and African descent populations. The most significant QRS duration association was with intronic SNP rs3922344 (P = 1.19x10-24) in SCN5A/SCN10A. Three other previously identified loci, CDKN1A, VTI1A, and HAND1, also exceeded the GWAS significance threshold among Hispanics/Latinos. A total of 27 of 32 previously identified QRS duration SNPs were shown to generalize in Hispanics/Latinos.

CONCLUSIONS

Our QRS duration GWAS, the first in Hispanic/Latino populations, identified two new loci, underscoring the utility of extending large scale genomic studies to currently under-examined populations.

An integrative cross-omics analysis of DNA methylation sites of glucose and insulin homeostasis

Despite existing reports on differential DNA methylation in type 2 diabetes (T2D) and obesity, our understanding of its functional relevance remains limited. Here we show the effect of differential methylation in the early phases of T2D pathology by a blood-based epigenome-wide association study of 4808 non-diabetic Europeans in the discovery phase and 11,750 individuals in the replication. We identify CpGs in LETM1, RBM20, IRS2, MAN2A2 and the 1q25.3 region associated with fasting insulin, and in FCRL6, SLAMF1, APOBEC3H and the 15q26.1 region with fasting glucose. In silico cross-omics analyses highlight the role of differential methylation in the crosstalk between the adaptive immune system and glucose homeostasis. The differential methylation explains at least 16.9% of the association between obesity and insulin. Our study sheds light on the biological interactions between genetic variants driving differential methylation and gene expression in the early pathogenesis of T2D.

Using Animations of Risk Functions to Visualize Trends in US All-Cause and Cause-Specific Mortality, 1968-2016

OBJECTIVES

To use dynamic visualizations of mortality risk functions over both calendar year and age as a way to estimate and visualize patterns in US life spans.

METHODS

We built 49 synthetic cohorts, 1 per year 1968 to 2016, using National Center for Health Statistics (NCHS) mortality and population data. Within each cohort, we estimated age-specific probabilities of dying from any cause (all-cause analysis) or from a particular cause (cause-specific analysis). We then used Kaplan-Meier (all-cause) or Aalen-Johansen (cause-specific) estimators to obtain risk functions. We illustrated risk functions using time-lapse animations.

RESULTS

Median age at death increased from 75 years in 1970 to 83 years in 2015. Risk by age 100 years of cardiovascular mortality decreased (from a risk of 55% in 1970 to 32% in 2015), whereas risk attributable to other (i.e., nonrespiratory and noncardiovascular) causes increased in compensation.

CONCLUSIONS

Our findings were consistent with the trends published in the NCHS 2015 mortality report, and our dynamic animations added an efficient, interpretable tool for visualizing US mortality trends over age and calendar time.

Association between DNA methylation in obesity-related genes and body mass index percentile in adolescents

Childhood obesity remains an epidemic in the U.S. and worldwide. However, little is understood regarding the epigenetic basis of obesity in adolescents. To investigate the cross-sectional association between DNA methylation level in obesity-related genes and body mass index (BMI) percentile, data from 263 adolescents in the population-based Penn State Child Cohort follow-up exam was analysed. Using DNA extracted from peripheral leukocytes, epigenome-wide single nucleotide resolution of DNA methylation in cytosine-phosphate-guanine (CpG) sites and surrounding regions was obtained. We used multivariable-adjusted linear regression models to assess the association between site-specific methylation level and age- and sex-specific BMI percentile. Hypergeometric and permutation tests were used to determine if obesity-related genes were significantly enriched among all intragenic sites that achieved a p < 0.05 throughout the epigenome. Among the 5,669 sites related to BMI percentile with p < 0.05, 28 were identified within obesity-related genes. Obesity-related genes were significantly enriched among 103,466 intragenic sites (P_{hypergeometric} = 0.006; P_permutation = 0.006). Moreover, increased methylation on one site within SIM1 was significantly related to higher BMI percentile (P = 4.2E-05). If externally validated, our data would suggest that DNA methylation in obesity-related genes may relate to obesity risk in adolescents.

Genome-wide meta-analysis of SNP-by9-ACEI/ARB and SNP-by-thiazide diuretic and effect on serum potassium in cohorts of European and African ancestry

We evaluated interactions of SNP-by-ACE-I/ARB and SNP-by-TD on serum potassium (K+) among users of antihypertensive treatments (anti-HTN). Our study included seven European-ancestry (EA) (N = 4835) and four African-ancestry (AA) cohorts (N = 2016). We performed race-stratified, fixed-effect, inverse-variance-weighted meta-analyses of 2.5 million SNP-by-drug interaction estimates; race-combined meta-analysis; and trans-ethnic fine-mapping. Among EAs, we identified 11 significant SNPs (P < 5 × 10-8) for SNP-ACE-I/ARB interactions on serum K+ that were located between NR2F1-AS1 and ARRDC3-AS1 on chromosome 5 (top SNP rs6878413 P = 1.7 × 10-8; ratio of serum K+ in ACE-I/ARB exposed compared to unexposed is 1.0476, 1.0280, 1.0088 for the TT, AT, and AA genotypes, respectively). Trans-ethnic fine mapping identified the same group of SNPs on chromosome 5 as genome-wide significant for the ACE-I/ARB analysis. In conclusion, SNP-by-ACE-I /ARB interaction analyses uncovered loci that, if replicated, could have future implications for the prevention of arrhythmias due to anti-HTN treatment-related hyperkalemia. Before these loci can be identified as clinically relevant, future validation studies of equal or greater size in comparison to our discovery effort are needed.

Rural-Urban Residence and Stage at Breast Cancer Diagnosis Among Postmenopausal Women: The Women's Health Initiative

BACKGROUND

Although social exposures have complex and dynamic relationships and interactions, the existing literature on the impact of rural-urban residence on stage at breast cancer diagnosis does not examine heterogeneity of effect. We examined the joint effect of social support, social relationship strain, and rural-urban residence on stage at breast cancer diagnosis.

METHODS

Using data from the Women's Health Initiative (WHI) (n = 161,808), we describe the distribution of social, behavioral, and clinical factors by rural-urban residence among postmenopausal women with incident breast cancer (n = 7,120). We used rural-urban commuting area (RUCA) codes to categorize baseline residential addresses as urban, large rural city/town, or small rural town, and the surveillance, epidemiology, and end results staging system to categorize breast cancer stage at diagnosis (dichotomized as early or late). We then used univariable and multivariable logistic regression to estimate odds ratios (ORs) and associated 95% confidence intervals (95% CI) for the relationship between rural-urban residence and stage at breast cancer diagnosis. We included separate interaction terms between rural-urban residence and social strain and social support to test for statistical interaction.

RESULTS

Of the social, behavioral, and clinical factors we examined, only younger age at WHI enrollment screening was significantly associated with late stage at breast cancer diagnosis (p = 0.003). Contrary to our hypothesis, rural-urban residence was not significantly associated with stage at breast cancer diagnosis among postmenopausal women ([adjusted OR, 95% CI] for urban compared with small town: 1.08 [0.76-1.53]; large town compared with small town: 1.16 [0.74-1.84]; and urban compared with large town: 0.93 [0.68-1.26]).The associations did not vary by social support or social strain (p for interaction between RUCA and social strain and social support, respectively: 0.99 and 0.17).

CONCLUSIONS

Future studies should examine other potential effect modifiers to identify novel factors predictive or protective for late stage at breast cancer diagnosis associated with rural-urban residence.