Heart Failure Risk Among African-American Women With an ICAM1 Missense Variant

BACKGROUND

A common genetic variant of ICAM1 among African-American individuals (rs5491; p.K56M) is associated with heart failure (HF) hospitalization, but whether this risk is specific to heart failure with preserved ejection fraction (HFpEF) remains unclear. Older women are at high risk for HFpEF, and the relationship between rs5491 and HFpEF across the age spectrum is unknown.

OBJECTIVES

This study assessed risk of HF and its subtypes conferred by ICAM1 p.K56M (rs5491).

METHODS

Associations of rs5491 with risk of HF and its subtypes were estimated among African American individuals in WHI (Women's Health Initiative). The study evaluated whether the association between rs5491 and HF hospitalizations was modified by baseline age. Subsequently, African-American women in WHI and MESA (Multi-Ethnic Study of Atherosclerosis) were pooled and analyses were repeated.

RESULTS

Among 8,401 women in WHI, the minor allele frequency of rs5491 was 20.7%, and 731 HF hospitalizations occurred over 19.2 years. The rs5491 variant was not associated with HF or its subtypes across WHI. Interaction analyses suggested that age as a continuous variable modified the association of rs5491 with HFpEF hospitalization (interaction P = 0.04). Upon categorizing women into age decades, rs5491 conferred increased risk of HFpEF among women ≥70 years (HR per additional rs5491 allele: 1.82 [95% CI: 1.25-2.65]; P = 0.002) but was not associated with HFpEF risk among women <70 years. Pooling African-American women in WHI (n = 8,401) and MESA (n = 856) demonstrated that the effect modification by age on the association of rs5491 with HFpEF became more significant (interaction P = 0.009), with consistent HFpEF risk effect estimates among women ≥70 years.

CONCLUSIONS

ICAM1 p.K56M (rs5491) is associated with HFpEF among African-American women ≥70 years.

A genetic association study of circulating coagulation Factor VIII and von Willebrand Factor levels

Coagulation Factor VIII (FVIII) and its carrier protein von Willebrand factor (VWF) are critical to coagulation and platelet aggregation. We leveraged whole genome sequence data from the Trans-Omics for Precision Medicine (TOPMed) program along with TOPMed-based imputation of genotypes in additional samples to identify genetic associations with circulating FVIII and VWF levels in a single variant meta-analysis including up to 45,289 participants. Gene-based aggregate tests were implemented in TOPMed. We identified three candidate causal genes and tested their functional effect on FVIII release from human liver endothelial cells (HLECs) and VWF release from human umbilical vein endothelial cells (HUVECs). Mendelian randomization was also performed to provide evidence for causal associations of FVIII and VWF with thrombotic outcomes. We identified associations (P<5×10-9) at seven new loci for FVIII (ST3GAL4, CLEC4M, B3GNT2, ASGR1, F12, KNG1, and TREM1/NCR2) and one for VWF (B3GNT2). VWF, ABO, and STAB2 were associated with FVIII and VWF in gene-based analyses. Multi-phenotype analysis of FVIII and VWF identified another three new loci, including PDIA3. Silencing of B3GNT2 and the previously reported CD36 gene decreased release of FVIII by HLECs, while silencing of B3GNT2, CD36, and PDIA3 decreased release of VWF by HVECs. Mendelian randomization supports causal association of higher FVIII and VWF with increased risk of thrombotic outcomes. Seven new loci were identified for FVIII and one for VWF, with evidence supporting causal associations of FVIII and VWF with thrombotic outcomes. B3GNT2, CD36, and PDIA3 modulate the release of FVIII and/or VWF in vitro.

Clonal Hematopoiesis of Indeterminate Potential (CHIP) and Incident Type 2 Diabetes Risk

OBJECTIVE

Clonal hematopoiesis of indeterminate potential (CHIP) is an aging-related accumulation of somatic mutations in hematopoietic stem cells, leading to clonal expansion. CHIP presence has been implicated in atherosclerotic coronary heart disease (CHD) and all-cause mortality, but its association with incident type 2 diabetes (T2D) is unknown. We hypothesized that CHIP is associated with elevated risk of T2D.

RESEARCH DESIGN AND METHODS

CHIP was derived from whole-genome sequencing of blood DNA in the National Heart, Lung, and Blood Institute Trans-Omics for Precision Medicine (TOPMed) prospective cohorts. We performed analysis for 17,637 participants from six cohorts, without prior T2D, cardiovascular disease, or cancer. We evaluated baseline CHIP versus no CHIP prevalence with incident T2D, including associations with DNMT3A, TET2, ASXL1, JAK2, and TP53 variants. We estimated multivariable-adjusted hazard ratios (HRs) and 95% CIs with adjustment for age, sex, BMI, smoking, alcohol, education, self-reported race/ethnicity, and combined cohorts' estimates via fixed-effects meta-analysis.

RESULTS

Mean (SD) age was 63.4 (11.5) years, 76% were female, and CHIP prevalence was 6.0% (n = 1,055) at baseline. T2D was diagnosed in n = 2,467 over mean follow-up of 9.8 years. Participants with CHIP had 23% (CI 1.04, 1.45) higher risk of T2D than those with no CHIP. Specifically, higher risk was for TET2 (HR 1.48; CI 1.05, 2.08) and ASXL1 (HR 1.76; CI 1.03, 2.99) mutations; DNMT3A was nonsignificant (HR 1.15; CI 0.93, 1.43). Statistical power was limited for JAK2 and TP53 analyses.

CONCLUSIONS

CHIP was associated with higher incidence of T2D. CHIP mutations located on genes implicated in CHD and mortality were also related to T2D, suggesting shared aging-related pathology.

Differences in the Circulating Proteome in Individuals with versus without Sickle Cell Trait

BACKGROUND

Sickle cell trait affects approximately 8% of Black individuals in the United States, along with many other individuals with ancestry from malaria-endemic regions worldwide. While traditionally considered a benign condition, recent evidence suggests that sickle cell trait is associated with lower eGFR and higher risk of kidney diseases, including kidney failure. The mechanisms underlying these associations remain poorly understood. We used proteomic profiling to gain insight into the pathobiology of sickle cell trait.

METHODS

We measured proteomics ( N =1285 proteins assayed by Olink Explore) using baseline plasma samples from 592 Black participants with sickle cell trait and 1:1 age-matched Black participants without sickle cell trait from the prospective Women's Health Initiative cohort. Age-adjusted linear regression was used to assess the association between protein levels and sickle cell trait.

RESULTS

In age-adjusted models, 35 proteins were significantly associated with sickle cell trait after correction for multiple testing. Several of the sickle cell trait-protein associations were replicated in Black participants from two independent cohorts (Atherosclerosis Risk in Communities study and Jackson Heart Study) assayed using an orthogonal aptamer-based proteomic platform (SomaScan). Many of the validated sickle cell trait-associated proteins are known biomarkers of kidney function or injury ( e.g. , hepatitis A virus cellular receptor 1 [HAVCR1]/kidney injury molecule-1 [KIM-1], uromodulin [UMOD], ephrins), related to red cell physiology or hemolysis (erythropoietin [EPO], heme oxygenase 1 [HMOX1], and α -hemoglobin stabilizing protein) and/or inflammation (fractalkine, C-C motif chemokine ligand 2/monocyte chemoattractant protein-1 [MCP-1], and urokinase plasminogen activator surface receptor [PLAUR]). A protein risk score constructed from the top sickle cell trait-associated biomarkers was associated with incident kidney failure among those with sickle cell trait during Women's Health Initiative follow-up (odds ratio, 1.32; 95% confidence interval, 1.10 to 1.58).

CONCLUSIONS

We identified and replicated the association of sickle cell trait with a number of plasma proteins related to hemolysis, kidney injury, and inflammation.

Deleterious heteroplasmic mitochondrial mutations are associated with an increased risk of overall and cancer-specific mortality

Mitochondria carry their own circular genome and disruption of the mitochondrial genome is associated with various aging-related diseases. Unlike the nuclear genome, mitochondrial DNA (mtDNA) can be present at 1000 s to 10,000 s copies in somatic cells and variants may exist in a state of heteroplasmy, where only a fraction of the DNA molecules harbors a particular variant. We quantify mtDNA heteroplasmy in 194,871 participants in the UK Biobank and find that heteroplasmy is associated with a 1.5-fold increased risk of all-cause mortality. Additionally, we functionally characterize mtDNA single nucleotide variants (SNVs) using a constraint-based score, mitochondrial local constraint score sum (MSS) and find it associated with all-cause mortality, and with the prevalence and incidence of cancer and cancer-related mortality, particularly leukemia. These results indicate that mitochondria may have a functional role in certain cancers, and mitochondrial heteroplasmic SNVs may serve as a prognostic marker for cancer, especially for leukemia.

Missense Genetic Variation of ICAM1 and Incident Heart Failure

BACKGROUND

Intercellular adhesion molecule-1 (ICAM-1) is a cell surface protein that participates in endothelial activation and is hypothesized to play a central role in heart failure (HF). We evaluated associations of ICAM1 missense genetic variants with circulating ICAM-1 levels and with incident HF.

METHODS AND RESULTS

We identified 3 missense variants within ICAM1 (rs5491, rs5498 and rs1799969) and evaluated their associations with ICAM-1 levels in the Coronary Artery Risk Development in Young Adults Study and the Multi-Ethnic Study of Atherosclerosis (MESA). We determined the association among these 3 variants and incident HF in MESA. We separately evaluated significant associations in the Atherosclerosis Risk in Communities (ARIC) study. Of the 3 missense variants, rs5491 was common in Black participants (minor allele frequency [MAF] > 20%) and rare in other race/ethnic groups (MAF < 5%). In Black participants, the presence of rs5491 was associated with higher levels of circulating ICAM-1 at 2 timepoints separated by 8 years. Among Black participants in MESA (n = 1600), the presence of rs5491 was associated with an increased risk of incident HF with preserved ejection fraction (HFpEF; HR = 2.30; [95% CI 1.25-4.21; P = 0.007]). The other ICAM1 missense variants (rs5498 and rs1799969) were associated with ICAM-1 levels, but there were no associations with HF. In ARIC, rs5491 was significantly associated with incident HF (HR = 1.24 [95% CI 1.02 - 1.51]; P = 0.03), with a similar direction of effect for HFpEF that was not statistically significant.

CONCLUSIONS

A common ICAM1 missense variant among Black individuals may be associated with increased risk of HF, which may be HFpEF-specific.

Association of Clonal Hematopoiesis of Indeterminate Potential with Incident Heart Failure with Preserved Ejection Fraction

Background

Clonal hematopoiesis of indeterminate potential (CHIP) was recently identified as a risk factor for incident heart failure (HF). Whether CHIP is associated selectively with risk of heart failure with reduced ejection fraction (HFrEF) or heart failure with preserved ejection fraction (HFpEF) subtypes is unknown.

Objectives

To evaluate whether CHIP is associated with incident HF subtypes, HFrEF versus HFpEF.

Methods

We obtained CHIP status from whole genome sequencing of blood DNA in participants without prevalent HF from a multi-ethnic sample of post-menopausal women without prevalent HF (N=5,214) from the Women's Health Initiative (WHI). Cox proportional hazards models were performed, adjusting for demographic and clinical risk factors.

Results

CHIP was significantly associated with a 42% (95%CI 6%, 91%) increased risk of HFpEF (P=0.02). In contrast, there was no evidence of association between CHIP and risk of incident HFrEF. When the three most common CHIP subtypes were assessed individually, the risk of HFpEF was more strongly associated with TET2 (HR=2.5; 95%CI 1.54, 4.06; P<0.001), than DNMT3A or ASXL1.

Conclusion

CHIP, particularly mutations in TET2, represents a potential new risk factor for incident HFpEF.

Canonical correlation analysis for multi-omics: Application to cross-cohort analysis

Integrative approaches that simultaneously model multi-omics data have gained increasing popularity because they provide holistic system biology views of multiple or all components in a biological system of interest. Canonical correlation analysis (CCA) is a correlation-based integrative method designed to extract latent features shared between multiple assays by finding the linear combinations of features-referred to as canonical variables (CVs)-within each assay that achieve maximal across-assay correlation. Although widely acknowledged as a powerful approach for multi-omics data, CCA has not been systematically applied to multi-omics data in large cohort studies, which has only recently become available. Here, we adapted sparse multiple CCA (SMCCA), a widely-used derivative of CCA, to proteomics and methylomics data from the Multi-Ethnic Study of Atherosclerosis (MESA) and Jackson Heart Study (JHS). To tackle challenges encountered when applying SMCCA to MESA and JHS, our adaptations include the incorporation of the Gram-Schmidt (GS) algorithm with SMCCA to improve orthogonality among CVs, and the development of Sparse Supervised Multiple CCA (SSMCCA) to allow supervised integration analysis for more than two assays. Effective application of SMCCA to the two real datasets reveals important findings. Applying our SMCCA-GS to MESA and JHS, we identified strong associations between blood cell counts and protein abundance, suggesting that adjustment of blood cell composition should be considered in protein-based association studies. Importantly, CVs obtained from two independent cohorts also demonstrate transferability across the cohorts. For example, proteomic CVs learned from JHS, when transferred to MESA, explain similar amounts of blood cell count phenotypic variance in MESA, explaining 39.0% ~ 50.0% variation in JHS and 38.9% ~ 49.1% in MESA. Similar transferability was observed for other omics-CV-trait pairs. This suggests that biologically meaningful and cohort-agnostic variation is captured by CVs. We anticipate that applying our SMCCA-GS and SSMCCA on various cohorts would help identify cohort-agnostic biologically meaningful relationships between multi-omics data and phenotypic traits.

Centenarian clocks: epigenetic clocks for validating claims of exceptional longevity

Claims surrounding exceptional longevity are sometimes disputed or dismissed for lack of credible evidence. Here, we present three DNA methylation-based age estimators (epigenetic clocks) for verifying age claims of centenarians. The three centenarian clocks were developed based on n = 7039 blood and saliva samples from individuals older than 40, including n = 184 samples from centenarians, 122 samples from semi-supercentenarians (aged 105 +), and 25 samples from supercentenarians (aged 110 +). The oldest individual was 115 years old. Our most accurate centenarian clock resulted from applying a neural network model to a training set composed of individuals older than 40. An epigenome-wide association study of age in different age groups revealed that age effects in young individuals (age < 40) are correlated (r = 0.55) with age effects in old individuals (age > 90). We present a chromatin state analysis of age effects in centenarians. The centenarian clocks are expected to be useful for validating claims surrounding exceptional old age.

DNAmFitAge: biological age indicator incorporating physical fitness

Physical fitness is a well-known correlate of health and the aging process and DNA methylation (DNAm) data can capture aging via epigenetic clocks. However, current epigenetic clocks did not yet use measures of mobility, strength, lung, or endurance fitness in their construction. We develop blood-based DNAm biomarkers for fitness parameters gait speed (walking speed), maximum handgrip strength, forced expiratory volume in one second (FEV1), and maximal oxygen uptake (VO2max) which have modest correlation with fitness parameters in five large-scale validation datasets (average r between 0.16-0.48). We then use these DNAm fitness parameter biomarkers with DNAmGrimAge, a DNAm mortality risk estimate, to construct DNAmFitAge, a new biological age indicator that incorporates physical fitness. DNAmFitAge is associated with low-intermediate physical activity levels across validation datasets (p = 6.4E-13), and younger/fitter DNAmFitAge corresponds to stronger DNAm fitness parameters in both males and females. DNAmFitAge is lower (p = 0.046) and DNAmVO2max is higher (p = 0.023) in male body builders compared to controls. Physically fit people have a younger DNAmFitAge and experience better age-related outcomes: lower mortality risk (p = 7.2E-51), coronary heart disease risk (p = 2.6E-8), and increased disease-free status (p = 1.1E-7). These new DNAm biomarkers provide researchers a new method to incorporate physical fitness into epigenetic clocks.

DNA methylation GrimAge version 2

We previously described a DNA methylation (DNAm) based biomarker of human mortality risk DNAm GrimAge. Here we describe version 2 of GrimAge (trained on individuals aged between 40 and 92) which leverages two new DNAm based estimators of (log transformed) plasma proteins: high sensitivity C-reactive protein (logCRP) and hemoglobin A1C (logA1C). We evaluate GrimAge2 in 13,399 blood samples across nine study cohorts. After adjustment for age and sex, GrimAge2 outperforms GrimAge in predicting mortality across multiple racial/ethnic groups (meta P=3.6x10^-167 versus P=2.6x10^-144) and in terms of associations with age related conditions such as coronary heart disease, lung function measurement FEV1 (correlation= -0.31, P=1.1x10^-136), computed tomography based measurements of fatty liver disease. We present evidence that GrimAge version 2 also applies to younger individuals and to saliva samples where it tracks markers of metabolic syndrome. DNAm logCRP is positively correlated with morbidity count (P=1.3x10^-54). DNAm logA1C is highly associated with type 2 diabetes (P=5.8x10^-155). DNAm PAI-1 outperforms the other age-adjusted DNAm biomarkers including GrimAge2 in correlating with triglyceride (cor=0.34, P=9.6x10^-267) and visceral fat (cor=0.41, P=4.7x10^-41). Overall, we demonstrate that GrimAge version 2 is an attractive epigenetic biomarker of human mortality and morbidity risk.

Whole genome sequencing identifies structural variants contributing to hematologic traits in the NHLBI TOPMed program

Genome-wide association studies have identified thousands of single nucleotide variants and small indels that contribute to variation in hematologic traits. While structural variants are known to cause rare blood or hematopoietic disorders, the genome-wide contribution of structural variants to quantitative blood cell trait variation is unknown. Here we utilized whole genome sequencing data in ancestrally diverse participants of the NHLBI Trans Omics for Precision Medicine program (N = 50,675) to detect structural variants associated with hematologic traits. Using single variant tests, we assessed the association of common and rare structural variants with red cell-, white cell-, and platelet-related quantitative traits and observed 21 independent signals (12 common and 9 rare) reaching genome-wide significance. The majority of these associations (N = 18) replicated in independent datasets. In genome-editing experiments, we provide evidence that a deletion associated with lower monocyte counts leads to disruption of an S1PR3 monocyte enhancer and decreased S1PR3 expression.

Genetic diversity fuels gene discovery for tobacco and alcohol use

Tobacco and alcohol use are heritable behaviours associated with 15% and 5.3% of worldwide deaths, respectively, due largely to broad increased risk for disease and injury1-4. These substances are used across the globe, yet genome-wide association studies have focused largely on individuals of European ancestries5. Here we leveraged global genetic diversity across 3.4 million individuals from four major clines of global ancestry (approximately 21% non-European) to power the discovery and fine-mapping of genomic loci associated with tobacco and alcohol use, to inform function of these loci via ancestry-aware transcriptome-wide association studies, and to evaluate the genetic architecture and predictive power of polygenic risk within and across populations. We found that increases in sample size and genetic diversity improved locus identification and fine-mapping resolution, and that a large majority of the 3,823 associated variants (from 2,143 loci) showed consistent effect sizes across ancestry dimensions. However, polygenic risk scores developed in one ancestry performed poorly in others, highlighting the continued need to increase sample sizes of diverse ancestries to realize any potential benefit of polygenic prediction.

Circulating Soluble CD163, Associations With Cardiovascular Outcomes and Mortality, and Identification of Genetic Variants in Older Individuals: The Cardiovascular Health Study

Background Monocytes/macrophages participate in cardiovascular disease. CD163 (cluster of differentiation 163) is a monocyte/macrophage receptor, and the shed sCD163 (soluble CD163) reflects monocyte/macrophage activation. We examined the association of sCD163 with incident cardiovascular disease events and performed a genome-wide association study to identify sCD163-associated variants. Methods and Results We measured plasma sCD163 in 5214 adults (aged ≥65 years, 58.7% women, 16.2% Black) of the CHS (Cardiovascular Health Study). We used Cox regression models (associations of sCD163 with incident events and mortality); median follow-up was 26 years. Genome-wide association study analyses were stratified on race. Adjusted for age, sex, and race and ethnicity, sCD163 levels were associated with all-cause mortality (hazard ratio [HR], 1.08 [95% CI, 1.04-1.12] per SD increase), cardiovascular disease mortality (HR, 1.15 [95% CI, 1.09-1.21]), incident coronary heart disease (HR, 1.10 [95% CI, 1.04-1.16]), and incident heart failure (HR, 1.18 [95% CI, 1.12-1.25]). When further adjusted (eg, cardiovascular disease risk factors), only incident coronary heart disease lost significance. In European American individuals, genome-wide association studies identified 38 variants on chromosome 2 near MGAT5 (top result rs62165726, P=3.3×10^-18),19 variants near chromosome 17 gene ASGR1 (rs55714927, P=1.5×10^-14), and 18 variants near chromosome 11 gene ST3GAL4. These regions replicated in the European ancestry ADDITION-PRO cohort, a longitudinal cohort study nested in the Danish arm of the Anglo-Danish-Dutch study of Intensive Treatment Intensive Treatment In peOple with screeNdetcted Diabetes in Primary Care. In Black individuals, we identified 9 variants on chromosome 6 (rs3129781 P=7.1×10^-9) in the HLA region, and 3 variants (rs115391969 P=4.3×10^-8) near the chromosome 16 gene MYLK3_. Conclusions Monocyte function, as measured by sCD163, may be predictive of overall and cardiovascular-specific mortality and incident heart failure.

Whole-exome sequencing of 14 389 individuals from the ESP and CHARGE consortia identifies novel rare variation associated with hemostatic factors

Plasma levels of fibrinogen, coagulation factors VII and VIII and von Willebrand factor (vWF) are four intermediate phenotypes that are heritable and have been associated with the risk of clinical thrombotic events. To identify rare and low-frequency variants associated with these hemostatic factors, we conducted whole-exome sequencing in 10 860 individuals of European ancestry (EA) and 3529 African Americans (AAs) from the Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium and the National Heart, Lung and Blood Institute's Exome Sequencing Project. Gene-based tests demonstrated significant associations with rare variation (minor allele frequency < 5%) in fibrinogen gamma chain (FGG) (with fibrinogen, P = 9.1 × 10-13), coagulation factor VII (F7) (with factor VII, P = 1.3 × 10-72; seven novel variants) and VWF (with factor VIII and vWF; P = 3.2 × 10-14; one novel variant). These eight novel rare variant associations were independent of the known common variants at these loci and tended to have much larger effect sizes. In addition, one of the rare novel variants in F7 was significantly associated with an increased risk of venous thromboembolism in AAs (Ile200Ser; rs141219108; P = 4.2 × 10-5). After restricting gene-based analyses to only loss-of-function variants, a novel significant association was detected and replicated between factor VIII levels and a stop-gain mutation exclusive to AAs (rs3211938) in CD36 molecule (CD36). This variant has previously been linked to dyslipidemia but not with the levels of a hemostatic factor. These efforts represent the largest integration of whole-exome sequence data from two national projects to identify genetic variation associated with plasma hemostatic factors.

Genome-wide association analyses of physical activity and sedentary behavior provide insights into underlying mechanisms and roles in disease prevention

Although physical activity and sedentary behavior are moderately heritable, little is known about the mechanisms that influence these traits. Combining data for up to 703,901 individuals from 51 studies in a multi-ancestry meta-analysis of genome-wide association studies yields 99 loci that associate with self-reported moderate-to-vigorous intensity physical activity during leisure time (MVPA), leisure screen time (LST) and/or sedentary behavior at work. Loci associated with LST are enriched for genes whose expression in skeletal muscle is altered by resistance training. A missense variant in ACTN3 makes the alpha-actinin-3 filaments more flexible, resulting in lower maximal force in isolated type IIA muscle fibers, and possibly protection from exercise-induced muscle damage. Finally, Mendelian randomization analyses show that beneficial effects of lower LST and higher MVPA on several risk factors and diseases are mediated or confounded by body mass index (BMI). Our results provide insights into physical activity mechanisms and its role in disease prevention.

Insights From a Large-Scale Whole-Genome Sequencing Study of Systolic Blood Pressure, Diastolic Blood Pressure, and Hypertension

BACKGROUND

The availability of whole-genome sequencing data in large studies has enabled the assessment of coding and noncoding variants across the allele frequency spectrum for their associations with blood pressure.

METHODS

We conducted a multiancestry whole-genome sequencing analysis of blood pressure among 51 456 Trans-Omics for Precision Medicine and Centers for Common Disease Genomics program participants (stage-1). Stage-2 analyses leveraged array data from UK Biobank (N=383 145), Million Veteran Program (N=318 891), and Reasons for Geographic and Racial Differences in Stroke (N=10 643) participants, along with whole-exome sequencing data from UK Biobank (N=199 631) participants.

RESULTS

Two blood pressure signals achieved genome-wide significance in meta-analyses of stage-1 and stage-2 single variant findings (P<5×10^-8). Among them, a rare intergenic variant at novel locus, LOC100506274, was associated with lower systolic blood pressure in stage-1 (beta [SE]=-32.6 [6.0]; P=4.99×10^-8) but not stage-2 analysis (P=0.11). Furthermore, a novel common variant at the known INSR locus was suggestively associated with diastolic blood pressure in stage-1 (beta [SE]=-0.36 [0.07]; P=4.18×10^-7) and attained genome-wide significance in stage-2 (beta [SE]=-0.29 [0.03]; P=7.28×10^-23). Nineteen additional signals suggestively associated with blood pressure in meta-analysis of single and aggregate rare variant findings (P<1×10^-6 and P<1×10^-4, respectively).

DISCUSSION

We report one promising but unconfirmed rare variant for blood pressure and, more importantly, contribute insights for future blood pressure sequencing studies. Our findings suggest promise of aggregate analyses to complement single variant analysis strategies and the need for larger, diverse samples, and family studies to enable robust rare variant identification.

TOP-LD: A tool to explore linkage disequilibrium with TOPMed whole-genome sequence data

Current publicly available tools that allow rapid exploration of linkage disequilibrium (LD) between markers (e.g., HaploReg and LDlink) are based on whole-genome sequence (WGS) data from 2,504 individuals in the 1000 Genomes Project. Here, we present TOP-LD, an online tool to explore LD inferred with high-coverage (∼30×) WGS data from 15,578 individuals in the NHLBI Trans-Omics for Precision Medicine (TOPMed) program. TOP-LD provides a significant upgrade compared to current LD tools, as the TOPMed WGS data provide a more comprehensive representation of genetic variation than the 1000 Genomes data, particularly for rare variants and in the specific populations that we analyzed. For example, TOP-LD encompasses LD information for 150.3, 62.2, and 36.7 million variants for European, African, and East Asian ancestral samples, respectively, offering 2.6- to 9.1-fold increase in variant coverage compared to HaploReg 4.0 or LDlink. In addition, TOP-LD includes tens of thousands of structural variants (SVs). We demonstrate the value of TOP-LD in fine-mapping at the GGT1 locus associated with gamma glutamyltransferase in the African ancestry participants in UK Biobank. Beyond fine-mapping, TOP-LD can facilitate a wide range of applications that are based on summary statistics and estimates of LD. TOP-LD is freely available online.

Intake and Sources of Dietary Fiber, Inflammation, and Cardiovascular Disease in Older US Adults

IMPORTANCE

Higher intake of dietary fiber has been associated with lower inflammation, but whether there are differences in this association by source of dietary fiber (ie, cereal, vegetable, or fruit) has not been studied to date.

OBJECTIVES

To evaluate the associations of total fiber intake and source (ie, cereal, vegetable, and fruit fiber intake) with inflammation and to evaluate whether inflammation mediates the inverse association between dietary fiber intake and cardiovascular disease (CVD).

DESIGN, SETTING, AND PARTICIPANTS

At the baseline visit (1989-1990) of 4125 adults aged 65 years or older in an ongoing US cohort study, dietary intake was assessed by a food frequency questionnaire among study participants without prevalent CVD (stroke and myocardial infarction) at enrollment. Inflammation was assessed from blood samples collected at baseline with immunoassays for markers of inflammation. Multivariable linear regression models tested the association of dietary fiber intake with inflammation. Also assessed was whether each inflammatory marker and its composite derived from principal component analysis mediated the association of baseline cereal fiber intake with development of CVD (stroke, myocardial infarction, and atherosclerotic cardiovascular death) through June 2015. Data from June 1, 1989, through June 30, 2015, were analyzed.

EXPOSURES

Total fiber intake and sources of fiber (cereal, vegetable, and fruit).

MAIN OUTCOMES AND MEASURES

Systemic markers of inflammation. Cardiovascular disease was the outcome in the mediation analysis.

RESULTS

Of 4125 individuals, 0.1% (n = 3) were Asian or Pacific Islander, 4.4% (n = 183) were Black, 0.3% (n = 12) were Native American, 95.0% (n = 3918) were White, and 0.2% (n = 9) were classified as other. Among these 4125 individuals (2473 women [60%]; mean [SD] age, 72.6 [5.5] years; 183 Black individuals [4.4%]; and 3942 individuals of other races and ethnicitites [95.6%] [ie, race and ethnicity other than Black, self-classified by participant]), an increase in total fiber intake of 5 g/d was associated with significantly lower concentrations of C-reactive protein (adjusted mean difference, -0.05 SD; 95% CI, -0.08 to -0.01 SD; P = .007) and interleukin 1 receptor antagonist (adjusted mean difference, -0.04 SD; 95% CI, -0.07 to -0.01 SD; P < .02) but with higher concentrations of soluble CD163 (adjusted mean difference, 0.05 SD; 95% CI, 0.02-0.09 SD; P = .005). Among fiber sources, only cereal fiber was consistently associated with lower inflammation. Similarly, cereal fiber intake was associated with lower CVD incidence (adjusted hazard ratio, 0.90; 95% CI, 0.81-1.00; 1941 incident cases). The proportion of the observed association of cereal fiber with CVD mediated by inflammatory markers ranged from 1.5% for interleukin 18 to 14.2% for C-reactive protein, and 16.1% for their primary principal component.

CONCLUSIONS AND RELEVANCE

Results of this study suggest that cereal fiber intake was associated with lower levels of various inflammatory markers and lower risk of CVD and that inflammation mediated approximately one-sixth of the association between cereal fiber intake and CVD.

Whole-Genome Sequencing Association Analyses of Stroke and Its Subtypes in Ancestrally Diverse Populations From Trans-Omics for Precision Medicine Project

BACKGROUND AND PURPOSE

Stroke is the leading cause of death and long-term disability worldwide. Previous genome-wide association studies identified 51 loci associated with stroke (mostly ischemic) and its subtypes among predominantly European populations. Using whole-genome sequencing in ancestrally diverse populations from the Trans-Omics for Precision Medicine (TOPMed) Program, we aimed to identify novel variants, especially low-frequency or ancestry-specific variants, associated with all stroke, ischemic stroke and its subtypes (large artery, cardioembolic, and small vessel), and hemorrhagic stroke and its subtypes (intracerebral and subarachnoid).

METHODS

Whole-genome sequencing data were available for 6833 stroke cases and 27 116 controls, including 22 315 European, 7877 Black, 2616 Hispanic/Latino, 850 Asian, 54 Native American, and 237 other ancestry participants. In TOPMed, we performed single variant association analysis examining 40 million common variants and aggregated association analysis focusing on rare variants. We also combined TOPMed European populations with over 28 000 additional European participants from the UK BioBank genome-wide array data through meta-analysis.

RESULTS

In the single variant association analysis in TOPMed, we identified one novel locus 13q33 for large artery at whole-genome-wide significance (P<5.00×10-9) and 4 novel loci at genome-wide significance (P<5.00×10-8), all of which need confirmation in independent studies. Lead variants in all 5 loci are low-frequency but are more common in non-European populations. An aggregation of synonymous rare variants within the gene C6orf26 demonstrated suggestive evidence of association for hemorrhagic stroke (P<3.11×10-6). By meta-analyzing European ancestry samples in TOPMed and UK BioBank, we replicated several previously reported stroke loci including PITX2, HDAC9, ZFHX3, and LRCH1.

CONCLUSIONS

We represent the first association analysis for stroke and its subtypes using whole-genome sequencing data from ancestrally diverse populations. While our findings suggest the potential benefits of combining whole-genome sequencing data with populations of diverse genetic backgrounds to identify possible low-frequency or ancestry-specific variants, they also highlight the need to increase genome coverage and sample sizes.

Rare coding variants in RCN3 are associated with blood pressure

BACKGROUND

While large genome-wide association studies have identified nearly one thousand loci associated with variation in blood pressure, rare variant identification is still a challenge. In family-based cohorts, genome-wide linkage scans have been successful in identifying rare genetic variants for blood pressure. This study aims to identify low frequency and rare genetic variants within previously reported linkage regions on chromosomes 1 and 19 in African American families from the Trans-Omics for Precision Medicine (TOPMed) program. Genetic association analyses weighted by linkage evidence were completed with whole genome sequencing data within and across TOPMed ancestral groups consisting of 60,388 individuals of European, African, East Asian, Hispanic, and Samoan ancestries.

RESULTS

Associations of low frequency and rare variants in RCN3 and multiple other genes were observed for blood pressure traits in TOPMed samples. The association of low frequency and rare coding variants in RCN3 was further replicated in UK Biobank samples (N = 403,522), and reached genome-wide significance for diastolic blood pressure (p = 2.01 × 10- 7).

CONCLUSIONS

Low frequency and rare variants in RCN3 contributes blood pressure variation. This study demonstrates that focusing association analyses in linkage regions greatly reduces multiple-testing burden and improves power to identify novel rare variants associated with blood pressure traits.

Epigenome-wide association study of serum urate reveals insights into urate co-regulation and the SLC2A9 locus

Elevated serum urate levels, a complex trait and major risk factor for incident gout, are correlated with cardiometabolic traits via incompletely understood mechanisms. DNA methylation in whole blood captures genetic and environmental influences and is assessed in transethnic meta-analysis of epigenome-wide association studies (EWAS) of serum urate (discovery, n = 12,474, replication, n = 5522). The 100 replicated, epigenome-wide significant (p < 1.1E-7) CpGs explain 11.6% of the serum urate variance. At SLC2A9, the serum urate locus with the largest effect in genome-wide association studies (GWAS), five CpGs are associated with SLC2A9 gene expression. Four CpGs at SLC2A9 have significant causal effects on serum urate levels and/or gout, and two of these partly mediate the effects of urate-associated GWAS variants. In other genes, including SLC7A11 and PHGDH, 17 urate-associated CpGs are associated with conditions defining metabolic syndrome, suggesting that these CpGs may represent a blood DNA methylation signature of cardiometabolic risk factors. This study demonstrates that EWAS can provide new insights into GWAS loci and the correlation of serum urate with other complex traits.

Meta-analyses identify DNA methylation associated with kidney function and damage

Chronic kidney disease is a major public health burden. Elevated urinary albumin-to-creatinine ratio is a measure of kidney damage, and used to diagnose and stage chronic kidney disease. To extend the knowledge on regulatory mechanisms related to kidney function and disease, we conducted a blood-based epigenome-wide association study for estimated glomerular filtration rate (n = 33,605) and urinary albumin-to-creatinine ratio (n = 15,068) and detected 69 and seven CpG sites where DNA methylation was associated with the respective trait. The majority of these findings showed directionally consistent associations with the respective clinical outcomes chronic kidney disease and moderately increased albuminuria. Associations of DNA methylation with kidney function, such as CpGs at JAZF1, PELI1 and CHD2 were validated in kidney tissue. Methylation at PHRF1, LDB2, CSRNP1 and IRF5 indicated causal effects on kidney function. Enrichment analyses revealed pathways related to hemostasis and blood cell migration for estimated glomerular filtration rate, and immune cell activation and response for urinary albumin-to-creatinineratio-associated CpGs.

Soluble Urokinase Plasminogen Activator Receptor: Genetic Variation and Cardiovascular Disease Risk in Blacks

BACKGROUND

suPAR (Soluble urokinase plasminogen activator receptor) has emerged as an important biomarker of coagulation, inflammation, and cardiovascular disease (CVD) risk. The contribution of suPAR to CVD risk and its genetic influence in the Black population have not been evaluated.

METHODS

We measured suPAR in 3492 Blacks from the prospective, community-based JHS (Jackson Heart Study). Cross-sectional associations of suPAR with lifestyle and CVD risk factors were assessed, whole-genome sequence data were used to evaluate genetic associations of suPAR, and relationships of suPAR with incident CVD outcomes and overall mortality were estimated over follow-up.

RESULTS

In Cox models adjusted for traditional CVD risk factors, estimated glomerular filtration rate, and CRP (C-reactive protein), each 1-SD higher suPAR was associated with a 21% to 31% increased risk of incident coronary heart disease, heart failure, stroke, and mortality. In the genome-wide association study, 2 missense (rs399145 encoding p.Thr86Ala, rs4760 encoding p.Phe272Leu) and 2 noncoding regulatory variants (rs73935023 within an enhancer element and rs4251805 within the promoter) of PLAUR on chromosome 19 were each independently associated with suPAR and together explained 14% of suPAR phenotypic variation. The allele frequencies of each of the four suPAR-associated genetic variants differ considerably across African and European populations. We further show that PLAUR rs73935023 can alter transcriptional activity in vitro. We did not find any association between genetically determined suPAR and CVD in JHS or a larger electronic medical record-based analyses of Blacks or Whites.

CONCLUSIONS

Our results demonstrate the importance of ancestry-differentiated genetic variation on suPAR levels and indicate suPAR is a CVD biomarker in Black adults.

A System for Phenotype Harmonization in the National Heart, Lung, and Blood Institute Trans-Omics for Precision Medicine (TOPMed) Program

Genotype-phenotype association studies often combine phenotype data from multiple studies to increase statistical power. Harmonization of the data usually requires substantial effort due to heterogeneity in phenotype definitions, study design, data collection procedures, and data-set organization. Here we describe a centralized system for phenotype harmonization that includes input from phenotype domain and study experts, quality control, documentation, reproducible results, and data-sharing mechanisms. This system was developed for the National Heart, Lung, and Blood Institute's Trans-Omics for Precision Medicine (TOPMed) program, which is generating genomic and other -omics data for more than 80 studies with extensive phenotype data. To date, 63 phenotypes have been harmonized across thousands of participants (recruited in 1948-2012) from up to 17 studies per phenotype. Here we discuss challenges in this undertaking and how they were addressed. The harmonized phenotype data and associated documentation have been submitted to National Institutes of Health data repositories for controlled access by the scientific community. We also provide materials to facilitate future harmonization efforts by the community, which include 1) the software code used to generate the 63 harmonized phenotypes, enabling others to reproduce, modify, or extend these harmonizations to additional studies, and 2) the results of labeling thousands of phenotype variables with controlled vocabulary terms.

Supplemental Association of Clonal Hematopoiesis With Incident Heart Failure

BACKGROUND

Age-related clonal hematopoiesis of indeterminate potential (CHIP), defined as clonally expanded leukemogenic sequence variations (particularly in DNMT3A, TET2, ASXL1, and JAK2) in asymptomatic individuals, is associated with cardiovascular events, including recurrent heart failure (HF).

OBJECTIVES

This study sought to evaluate whether CHIP is associated with incident HF.

METHODS

CHIP status was obtained from whole exome or genome sequencing of blood DNA in participants without prevalent HF or hematological malignancy from 5 cohorts. Cox proportional hazards models were performed within each cohort, adjusting for demographic and clinical risk factors, followed by fixed-effect meta-analyses. Large CHIP clones (defined as variant allele frequency >10%), HF with or without baseline coronary heart disease, and left ventricular ejection fraction were evaluated in secondary analyses.

RESULTS

Of 56,597 individuals (59% women, mean age 58 years at baseline), 3,406 (6%) had CHIP, and 4,694 developed HF (8.3%) over up to 20 years of follow-up. CHIP was prospectively associated with a 25% increased risk of HF in meta-analysis (hazard ratio: 1.25; 95% confidence interval: 1.13-1.38) with consistent associations across cohorts. ASXL1, TET2, and JAK2 sequence variations were each associated with an increased risk of HF, whereas DNMT3A sequence variations were not associated with HF. Secondary analyses suggested large CHIP was associated with a greater risk of HF (hazard ratio: 1.29; 95% confidence interval: 1.15-1.44), and the associations for CHIP on HF with and without prior coronary heart disease were homogenous. ASXL1 sequence variations were associated with reduced left ventricular ejection fraction.

CONCLUSIONS

CHIP, particularly sequence variations in ASXL1, TET2, and JAK2, represents a new risk factor for HF.

Genome-wide association studies identify 137 genetic loci for DNA methylation biomarkers of aging

BACKGROUND

Biological aging estimators derived from DNA methylation data are heritable and correlate with morbidity and mortality. Consequently, identification of genetic and environmental contributors to the variation in these measures in populations has become a major goal in the field.

RESULTS

Leveraging DNA methylation and SNP data from more than 40,000 individuals, we identify 137 genome-wide significant loci, of which 113 are novel, from genome-wide association study (GWAS) meta-analyses of four epigenetic clocks and epigenetic surrogate markers for granulocyte proportions and plasminogen activator inhibitor 1 levels, respectively. We find evidence for shared genetic loci associated with the Horvath clock and expression of transcripts encoding genes linked to lipid metabolism and immune function. Notably, these loci are independent of those reported to regulate DNA methylation levels at constituent clock CpGs. A polygenic score for GrimAge acceleration showed strong associations with adiposity-related traits, educational attainment, parental longevity, and C-reactive protein levels.

CONCLUSION

This study illuminates the genetic architecture underlying epigenetic aging and its shared genetic contributions with lifestyle factors and longevity.

BinomiRare: A robust test for association of a rare genetic variant with a binary outcome for mixed models and any case-control proportion

Whole-genome sequencing (WGS) and whole-exome sequencing studies have become increasingly available and are being used to identify rare genetic variants associated with health and disease outcomes. Investigators routinely use mixed models to account for genetic relatedness or other clustering variables (e.g., family or household) when testing genetic associations. However, no existing tests of the association of a rare variant with a binary outcome in the presence of correlated data control the type 1 error where there are (1) few individuals harboring the rare allele, (2) a small proportion of cases relative to controls, and (3) covariates to adjust for. Here, we address all three issues in developing a framework for testing rare variant association with a binary trait in individuals harboring at least one risk allele. In this framework, we estimate outcome probabilities under the null hypothesis and then use them, within the individuals with at least one risk allele, to test variant associations. We extend the BinomiRare test, which was previously proposed for independent observations, and develop the Conway-Maxwell-Poisson (CMP) test and study their properties in simulations. We show that the BinomiRare test always controls the type 1 error, while the CMP test sometimes does not. We then use the BinomiRare test to test the association of rare genetic variants in target genes with small-vessel disease (SVD) stroke, short sleep, and venous thromboembolism (VTE), in whole-genome sequence data from the Trans-Omics for Precision Medicine (TOPMed) program.

Clonal hematopoiesis associated with epigenetic aging and clinical outcomes

Clonal hematopoiesis of indeterminate potential (CHIP) is a common precursor state for blood cancers that most frequently occurs due to mutations in the DNA‐methylation modifying enzymes DNMT3A or TET2. We used DNA‐methylation array and whole‐genome sequencing data from four cohorts together comprising 5522 persons to study the association between CHIP, epigenetic clocks, and health outcomes. CHIP was strongly associated with epigenetic age acceleration, defined as the residual after regressing epigenetic clock age on chronological age, in several clocks, ranging from 1.31 years (GrimAge, p < 8.6 × 10⁻⁷) to 3.08 years (EEAA, p < 3.7 × 10⁻¹⁸). Mutations in most CHIP genes except DNA‐damage response genes were associated with increases in several measures of age acceleration. CHIP carriers with mutations in multiple genes had the largest increases in age acceleration and decrease in estimated telomere length. Finally, we found that ~40% of CHIP carriers had acceleration >0 in both Hannum and GrimAge (referred to as AgeAccelHG+). This group was at high risk of all‐cause mortality (hazard ratio 2.90, p < 4.1 × 10⁻⁸) and coronary heart disease (CHD) (hazard ratio 3.24, p < 9.3 × 10⁻⁶) compared to those who were CHIP−/AgeAccelHG−. In contrast, the other ~60% of CHIP carriers who were AgeAccelHG− were not at increased risk of these outcomes. In summary, CHIP is strongly linked to age acceleration in multiple clocks, and the combination of CHIP and epigenetic aging may be used to identify a population at high risk for adverse outcomes and who may be a target for clinical interventions.

Soluble CD14 Levels in the Jackson Heart Study: Associations With Cardiovascular Disease Risk and Genetic Variants

[Figure: see text].

Multi-Ancestry Genome-wide Association Study Accounting for Gene-Psychosocial Factor Interactions Identifies Novel Loci for Blood Pressure Traits

Psychological and social factors are known to influence blood pressure (BP) and risk of hypertension and associated cardiovascular diseases. To identify novel BP loci, we carried out genome-wide association meta-analyses of systolic, diastolic, pulse, and mean arterial BP taking into account the interaction effects of genetic variants with three psychosocial factors: depressive symptoms, anxiety symptoms, and social support. Analyses were performed using a two-stage design in a sample of up to 128,894 adults from 5 ancestry groups. In the combined meta-analyses of Stages 1 and 2, we identified 59 loci (p value <5e-8), including nine novel BP loci. The novel associations were observed mostly with pulse pressure, with fewer observed with mean arterial pressure. Five novel loci were identified in African ancestry, and all but one showed patterns of interaction with at least one psychosocial factor. Functional annotation of the novel loci supports a major role for genes implicated in the immune response (PLCL2), synaptic function and neurotransmission (LIN7A, PFIA2), as well as genes previously implicated in neuropsychiatric or stress-related disorders (FSTL5, CHODL). These findings underscore the importance of considering psychological and social factors in gene discovery for BP, especially in non-European populations.

Multi-Ethnic Genome-Wide Association Study of Decomposed Cardioelectric Phenotypes Illustrates Strategies to Identify and Characterize Evidence of Shared Genetic Effects for Complex Traits

BACKGROUND

We examined how expanding electrocardiographic trait genome-wide association studies to include ancestrally diverse populations, prioritize more precise phenotypic measures, and evaluate evidence for shared genetic effects enabled the detection and characterization of loci.

METHODS

We decomposed 10 seconds, 12-lead electrocardiograms from 34 668 multi-ethnic participants (15% Black; 30% Hispanic/Latino) into 6 contiguous, physiologically distinct (P wave, PR segment, QRS interval, ST segment, T wave, and TP segment) and 2 composite, conventional (PR interval and QT interval) interval scale traits and conducted multivariable-adjusted, trait-specific univariate genome-wide association studies using 1000-G imputed single-nucleotide polymorphisms. Evidence of shared genetic effects was evaluated by aggregating meta-analyzed univariate results across the 6 continuous electrocardiographic traits using the combined phenotype adaptive sum of powered scores test.

RESULTS

We identified 6 novels (CD36, PITX2, EMB, ZNF592, YPEL2, and BC043580) and 87 known loci (adaptive sum of powered score test P<5×10-9). Lead single-nucleotide polymorphism rs3211938 at CD36 was common in Blacks (minor allele frequency=10%), near monomorphic in European Americans, and had effects on the QT interval and TP segment that ranked among the largest reported to date for common variants. The other 5 novel loci were observed when evaluating the contiguous but not the composite electrocardiographic traits. Combined phenotype testing did not identify novel electrocardiographic loci unapparent using traditional univariate approaches, although this approach did assist with the characterization of known loci.

CONCLUSIONS

Despite including one-third as many participants as published electrocardiographic trait genome-wide association studies, our study identified 6 novel loci, emphasizing the importance of ancestral diversity and phenotype resolution in this era of ever-growing genome-wide association studies.

Role of Rare and Low-Frequency Variants in Gene-Alcohol Interactions on Plasma Lipid Levels

BACKGROUND

Alcohol intake influences plasma lipid levels, and such effects may be moderated by genetic variants. We aimed to characterize the role of aggregated rare and low-frequency protein-coding variants in gene by alcohol consumption interactions associated with fasting plasma lipid levels.

METHODS

In the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium, fasting plasma triglycerides and high- and low-density lipoprotein cholesterol were measured in 34 153 individuals with European ancestry from 5 discovery studies and 32 277 individuals from 6 replication studies. Rare and low-frequency functional protein-coding variants (minor allele frequency, ≤5%) measured by an exome array were aggregated by genes and evaluated by a gene-environment interaction test and a joint test of genetic main and gene-environment interaction effects. Two dichotomous self-reported alcohol consumption variables, current drinker, defined as any recurrent drinking behavior, and regular drinker, defined as the subset of current drinkers who consume at least 2 drinks per week, were considered.

RESULTS

We discovered and replicated 21 gene-lipid associations at 13 known lipid loci through the joint test. Eight loci (PCSK9, LPA, LPL, LIPG, ANGPTL4, APOB, APOC3, and CD300LG) remained significant after conditioning on the common index single-nucleotide polymorphism identified by previous genome-wide association studies, suggesting an independent role for rare and low-frequency variants at these loci. One significant gene-alcohol interaction on triglycerides in a novel locus was significantly discovered (P=6.65×10^-6 for the interaction test) and replicated at nominal significance level (P=0.013) in SMC5.

CONCLUSIONS

In conclusion, this study applied new gene-based statistical approaches and suggested that rare and low-frequency genetic variants interacted with alcohol consumption on lipid levels.

Comparison of Proteomic Assessment Methods in Multiple Cohort Studies

Novel proteomics platforms, such as the aptamer-based SOMAscan platform, can quantify large numbers of proteins efficiently and cost-effectively and are rapidly growing in popularity. However, comparisons to conventional immunoassays remain underexplored, leaving investigators unsure when cross-assay comparisons are appropriate. The correlation of results from immunoassays with relative protein quantification is explored by SOMAscan. For 63 proteins assessed in two chronic obstructive pulmonary disease (COPD) cohorts, subpopulations and intermediate outcome measures in COPD Study (SPIROMICS), and COPDGene, using myriad rules based medicine multiplex immunoassays and SOMAscan, Spearman correlation coefficients range from -0.13 to 0.97, with a median correlation coefficient of ≈0.5 and consistent results across cohorts. A similar range is observed for immunoassays in the population-based Multi-Ethnic Study of Atherosclerosis and for other assays in COPDGene and SPIROMICS. Comparisons of relative quantification from the antibody-based Olink platform and SOMAscan in a small cohort of myocardial infarction patients also show a wide correlation range. Finally, cis pQTL data, mass spectrometry aptamer confirmation, and other publicly available data are integrated to assess relationships with observed correlations. Correlation between proteomics assays shows a wide range and should be carefully considered when comparing and meta-analyzing proteomics data across assays and studies.

Coagulation factor VIII: Relationship to cardiovascular disease risk and whole genome sequence and epigenome-wide analysis in African Americans

BACKGROUND

Prospective studies have suggested higher factor VIII (FVIII) levels are an independent risk factor for coronary heart disease (CHD) and stroke. However, limited information, including on genetic and epigenetic contributors to FVIII variation, is available specifically among African Americans (AAs), who have higher FVIII levels than Europeans.

OBJECTIVES

We measured FVIII levels in ~3400 AAs from the community-based Jackson Heart Study and assessed genetic, epigenetic, and epidemiological correlates of FVIII, as well as incident cardiovascular disease (CVD) associations.

METHODS

We assessed cross-sectional associations of FVIII with CVD risk factors as well as incident CHD, stroke, heart failure, and mortality associations. We additionally assessed associations with TOPMed whole genome sequencing data and an epigenome-wide methylation array.

RESULTS

Our results confirmed associations between FVIII and risk of incident CHD events and total mortality in AAs; mortality associations were largely independent of traditional risk factors. We also demonstrate an association of FVIII with incident heart failure, independent of B-type natriuretic peptide. Two genomic regions were strongly associated with FVIII (ABO and VWF). The index variant at VWF is specific to individuals of African descent and is distinct from the previously reported European VWF association signal. Epigenome-wide association analysis showed significant FVIII associations with several CpG sites in the ABO region. However, after adjusting for ABO genetic variants, ABO CpG sites were not significant.

CONCLUSIONS

Larger sample sizes of AAs will be required to discover additional genetic and epigenetic contributors to FVIII phenotypic variation, which may have consequences for CVD health disparities.

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.

Interferon gamma-induced protein 10 (IP-10) and cardiovascular disease in African Americans

Biomarkers of chronic inflammation (such as C-reactive protein) have long been associated with cardiovascular disease and mortality; however, biomarkers involved in antiviral cytokine induction and adaptive immune system activation remain largely unexamined. We hypothesized the cytokine interferon gamma inducible protein 10 (IP-10) would be associated with clinical and subclinical cardiovascular disease and all-cause mortality in African Americans. We assessed these associations in the Jackson Heart Study (JHS) cohort and the REasons for Geographic and Racial Differences in Stroke (REGARDS) study. There was a modest association of IP-10 with higher odds of left ventricular hypertrophy (OR = 1.20 (95% confidence interval (CI) 1.03, 1.41) per standard deviation (SD) higher natural log-transformed IP-10 in JHS). We did not observe associations with ankle brachial index, intima-media thickness, or arterial calcification. Each SD higher increment of ln-transformed IP-10 concentration was associated with incident heart failure (hazard ratio (HR) 1.26; 95% CI 1.11, 1.42, p = 4x10-4) in JHS, and with overall mortality in both JHS (HR 1.12 per SD, 95% CI 1.03, 1.21, p = 7.5x10-3) and REGARDS (HR 1.31 per SD, 95% CI 1.10, 1.55, p = 2.0 x 10-3), adjusting for cardiovascular risk factors and C-reactive protein. However, we found no association between IP-10 and stroke or coronary heart disease. These results suggest a role of IP-10 in heart failure and mortality risk independent of C-reactive protein. Further research is needed to investigate how the body's response to chronic viral infection may mediate heart failure and overall mortality risk in African Americans.

Ancestry-specific associations identified in genome-wide combined-phenotype study of red blood cell traits emphasize benefits of diversity in genomics

BACKGROUND

Quantitative red blood cell (RBC) traits are highly polygenic clinically relevant traits, with approximately 500 reported GWAS loci. The majority of RBC trait GWAS have been performed in European- or East Asian-ancestry populations, despite evidence that rare or ancestry-specific variation contributes substantially to RBC trait heritability. Recently developed combined-phenotype methods which leverage genetic trait correlation to improve statistical power have not yet been applied to these traits. Here we leveraged correlation of seven quantitative RBC traits in performing a combined-phenotype analysis in a multi-ethnic study population.

RESULTS

We used the adaptive sum of powered scores (aSPU) test to assess combined-phenotype associations between ~ 21 million SNPs and seven RBC traits in a multi-ethnic population (maximum n = 67,885 participants; 24% African American, 30% Hispanic/Latino, and 43% European American; 76% female). Thirty-nine loci in our multi-ethnic population contained at least one significant association signal (p < 5E-9), with lead SNPs at nine loci significantly associated with three or more RBC traits. A majority of the lead SNPs were common (MAF > 5%) across all ancestral populations. Nineteen additional independent association signals were identified at seven known loci (HFE, KIT, HBS1L/MYB, CITED2/FILNC1, ABO, HBA1/2, and PLIN4/5). For example, the HBA1/2 locus contained 14 conditionally independent association signals, 11 of which were previously unreported and are specific to African and Amerindian ancestries. One variant in this region was common in all ancestries, but exhibited a narrower LD block in African Americans than European Americans or Hispanics/Latinos. GTEx eQTL analysis of all independent lead SNPs yielded 31 significant associations in relevant tissues, over half of which were not at the gene immediately proximal to the lead SNP.

CONCLUSION

This work identified seven loci containing multiple independent association signals for RBC traits using a combined-phenotype approach, which may improve discovery in genetically correlated traits. Highly complex genetic architecture at the HBA1/2 locus was only revealed by the inclusion of African Americans and Hispanics/Latinos, underscoring the continued importance of expanding large GWAS to include ancestrally diverse populations.

Soluble CD14, Ischemic Stroke, and Coronary Heart Disease Risk in a Prospective Study: The REGARDS Cohort

Background

Soluble CD14 (sCD14), a circulating pattern recognition receptor, has been suggested as a cardiovascular disease risk factor. Prospective studies evaluating sCD14 with incident cardiovascular disease events are limited, particularly among racially diverse populations.

Methods and Results

Between 2003 and 2007, the REGARDS (Reasons for Geographic and Racial Differences in Stroke) study recruited 30 239 black and white participants across the United States. In a nested case–cohort study, sCD14 was measured in baseline serum from 548 cases of incident ischemic stroke, 612 cases of incident coronary heart disease (CHD), and a cohort random sample (n=1039). Cox models estimated hazards ratios (HR) of incident ischemic stroke or CHD per 1 SD higher sCD14, adjusting for cardiovascular disease risk factors. There was a differential association of sCD14 with ischemic stroke and CHD risk by race. Among blacks, the adjusted HR of stroke per SD increment of sCD14 was 1.42 (95% CI: 1.12, 1.80), with no association among whites (HR 1.02 [95% CI: 0.82, 1.27]). Higher sCD14 was associated with increased CHD risk in blacks but not whites, and relationships between sCD14 and CHD were stronger at younger ages. Adjusted for risk factors, the HR of CHD per SD higher sCD14 among blacks at age 45 years was 2.30 (95% CI: 1.45, 3.65) compared with 1.56 (95% CI: 0.94, 2.57) among whites. At age 65 years, the CHD HR was 1.51 (95% CI: 1.20, 1.91) among blacks and 1.02 (95% CI: 0.80, 1.31) among whites.

Conclusions

sCD14 may be a race‐specific stroke and CHD risk marker.

Soluble CD14 and Risk of Heart Failure and Its Subtypes in Older Adults

BACKGROUND

CD14 is a membrane glycoprotein primarily expressed by myeloid cells that plays a key role in inflammation. Soluble CD14 (sCD14) levels carry a poor prognosis in chronic heart failure (HF), but whether elevations in sCD14 precede HF is unknown. We tested the hypothesis that sCD14 is associated with HF incidence and its subtypes independent of major inflammatory biomarkers among older adults.

METHODS AND RESULTS

We included participants in the Cardiovascular Health Study without preexisting HF and available baseline sCD14. We evaluated the associations of sCD14, high-sensitivity C-reactive protein (hsCRP), interleukin (IL)-6, and white blood cell count (WBC) with incident HF and subtypes using Cox regression. Among 5217 participants, 1878 had incident HF over 13.6 years (609 classifiable as HF with preserved ejection fraction [HFpEF] and 419 as HF with reduced ejection fraction [HFrEF]). After adjusting for clinical and laboratory covariates, sCD14 was significantly associated with incident HF (hazard ratio [HR]: 1.56 per doubling, 95% confidence interval [CI]: 1.29-1.89), an association that was numerically stronger than for hsCRP (HR per doubling: 1.10, 95% CI: 1.06-1.15), IL-6 (HR: 1.18, 95% CI: 1.10-1.25), and WBC (HR: 1.24, 95% CI: 1.09-1.42), and that remained significant after adjustment for the other markers of inflammation. This association for sCD14 was observed with HFpEF (HR: 1.50, 95% CI: 1.07-2.10) but not HFrEF (HR: 0.99, 95% CI: 0.67-1.49).

CONCLUSIONS

Plasma sCD14 was associated with incident HF independently and numerically more strongly than other major inflammatory markers. This association was only observed with HFpEF in the subset with classifiable HF subtypes. Pending replication, these findings have potentially important therapeutic implications.

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 9,925 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⁻⁸) 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.

Genome-Wide Association Transethnic Meta-Analyses Identifies Novel Associations Regulating Coagulation Factor VIII and von Willebrand Factor Plasma Levels

BACKGROUND

Factor VIII (FVIII) and its carrier protein von Willebrand factor (VWF) are associated with risk of arterial and venous thrombosis and with hemorrhagic disorders. We aimed to identify and functionally test novel genetic associations regulating plasma FVIII and VWF.

METHODS

We meta-analyzed genome-wide association results from 46 354 individuals of European, African, East Asian, and Hispanic ancestry. All studies performed linear regression analysis using an additive genetic model and associated ≈35 million imputed variants with natural log-transformed phenotype levels. In vitro gene silencing in cultured endothelial cells was performed for candidate genes to provide additional evidence on association and function. Two-sample Mendelian randomization analyses were applied to test the causal role of FVIII and VWF plasma levels on the risk of arterial and venous thrombotic events.

RESULTS

We identified 13 novel genome-wide significant ( P≤2.5×10^-8) associations, 7 with FVIII levels ( FCHO2/TMEM171/TNPO1, HLA, SOX17/RP1, LINC00583/NFIB, RAB5C-KAT2A, RPL3/TAB1/SYNGR1, and ARSA) and 11 with VWF levels ( PDHB/PXK/KCTD6, SLC39A8, FCHO2/TMEM171/TNPO1, HLA, GIMAP7/GIMAP4, OR13C5/NIPSNAP, DAB2IP, C2CD4B, RAB5C-KAT2A, TAB1/SYNGR1, and ARSA), beyond 10 previously reported associations with these phenotypes. Functional validation provided further evidence of association for all loci on VWF except ARSA and DAB2IP. Mendelian randomization suggested causal effects of plasma FVIII activity levels on venous thrombosis and coronary artery disease risk and plasma VWF levels on ischemic stroke risk.

CONCLUSIONS

The meta-analysis identified 13 novel genetic loci regulating FVIII and VWF plasma levels, 10 of which we validated functionally. We provide some evidence for a causal role of these proteins in thrombotic events.

Multi-ancestry sleep-by-SNP interaction analysis in 126,926 individuals reveals lipid loci stratified by sleep duration

Both short and long sleep are associated with an adverse lipid profile, likely through different biological pathways. To elucidate the biology of sleep-associated adverse lipid profile, we conduct multi-ancestry genome-wide sleep-SNP interaction analyses on three lipid traits (HDL-c, LDL-c and triglycerides). In the total study sample (discovery + replication) of 126,926 individuals from 5 different ancestry groups, when considering either long or short total sleep time interactions in joint analyses, we identify 49 previously unreported lipid loci, and 10 additional previously unreported lipid loci in a restricted sample of European-ancestry cohorts. In addition, we identify new gene-sleep interactions for known lipid loci such as LPL and PCSK9. The previously unreported lipid loci have a modest explained variance in lipid levels: most notable, gene-short-sleep interactions explain 4.25% of the variance in triglyceride level. Collectively, these findings contribute to our understanding of the biological mechanisms involved in sleep-associated adverse lipid profiles.

Sleep duration is associated with an adverse lipid profile. Here, the authors perform genome-wide gene-by-sleep interaction analysis and find 49 previously unreported lipid loci when considering short or long total sleep time.

Genome-wide association analysis of venous thromboembolism identifies new risk loci and genetic overlap with arterial vascular disease

Venous thromboembolism is a significant cause of mortality1, yet its genetic determinants are incompletely defined. We performed a discovery genome-wide association study in the Million Veteran Program and UK Biobank, with testing of approximately 13 million DNA sequence variants for association with venous thromboembolism (26,066 cases and 624,053 controls) and meta-analyzed both studies, followed by independent replication with up to 17,672 venous thromboembolism cases and 167,295 controls. We identified 22 previously unknown loci, bringing the total number of venous thromboembolism-associated loci to 33, and subsequently fine-mapped these associations. We developed a genome-wide polygenic risk score for venous thromboembolism that identifies 5% of the population at an equivalent incident venous thromboembolism risk to carriers of the established factor V Leiden p.R506Q and prothrombin G20210A mutations. Our data provide mechanistic insights into the genetic epidemiology of venous thromboembolism and suggest a greater overlap among venous and arterial cardiovascular disease than previously thought.

Phenome-wide association analysis of LDL-cholesterol lowering genetic variants in PCSK9

BACKGROUND

We characterised the phenotypic consequence of genetic variation at the PCSK9 locus and compared findings with recent trials of pharmacological inhibitors of PCSK9.

METHODS

Published and individual participant level data (300,000+ participants) were combined to construct a weighted PCSK9 gene-centric score (GS). Seventeen randomized placebo controlled PCSK9 inhibitor trials were included, providing data on 79,578 participants. Results were scaled to a one mmol/L lower LDL-C concentration.

RESULTS

The PCSK9 GS (comprising 4 SNPs) associations with plasma lipid and apolipoprotein levels were consistent in direction with treatment effects. The GS odds ratio (OR) for myocardial infarction (MI) was 0.53 (95% CI 0.42; 0.68), compared to a PCSK9 inhibitor effect of 0.90 (95% CI 0.86; 0.93). For ischemic stroke ORs were 0.84 (95% CI 0.57; 1.22) for the GS, compared to 0.85 (95% CI 0.78; 0.93) in the drug trials. ORs with type 2 diabetes mellitus (T2DM) were 1.29 (95% CI 1.11; 1.50) for the GS, as compared to 1.00 (95% CI 0.96; 1.04) for incident T2DM in PCSK9 inhibitor trials. No genetic associations were observed for cancer, heart failure, atrial fibrillation, chronic obstructive pulmonary disease, or Alzheimer's disease - outcomes for which large-scale trial data were unavailable.

CONCLUSIONS

Genetic variation at the PCSK9 locus recapitulates the effects of therapeutic inhibition of PCSK9 on major blood lipid fractions and MI. While indicating an increased risk of T2DM, no other possible safety concerns were shown; although precision was moderate.

Association of APOL1 With Heart Failure With Preserved Ejection Fraction in Postmenopausal African American Women

Importance

APOL1 genotypes are associated with kidney diseases in African American individuals and may influence cardiovascular disease and mortality risk, but findings have been inconsistent.

Objective

To discern whether high-risk APOL1 genotypes are associated with cardiovascular disease and stroke in postmenopausal African American women, who are at high risk for these outcomes.

Design, Setting, and Participants

The Women's Health Initiative is a prospective cohort that enrolled 161 838 postmenopausal women into clinical trials and an observational study between 1993 and 1998. This study includes 11 137 African American women participants who had a clinical event from enrollment to June 2014. Data analyses were completed from January 2017 to August 2017.

Exposures

The variants of APOL1 were genotyped or imputed from whole-exome sequencing.

Main Outcomes and Measures

Incident coronary heart disease, stroke and heart failure subtypes, and overall and cause-specific mortality were adjudicated from hospital records and death certificates. Estimated incidence rates were determined for each outcome and hazard ratios (HR) and 95% CIs for the associations of APOL1 groups with outcomes.

Results

The mean (SD) age of participants was 61.7 (7.1) years. Carriers of high-risk APOL1 variants (n = 1370; 12.3%) had higher prevalence of hypertension, use of cholesterol-lowering medications, and reduced estimated glomerular filtration rate (eGFR). After a mean (SD) of 11.0 (3.6) years, carriers of high-risk APOL1 variants had a higher incidence rate of hospitalized heart failure with preserved ejection fraction (HFpEF) than low-risk carriers did but showed no differences for other outcomes. In adjusted models, there was a significant 58% increased hazard of hospitalized HFpEF (HR, 1.58 [95% CI, 1.03-2.41]) among carriers of high-risk APOL1 variants compared with carriers of low-risk APOL1 variants. The association with HFpEF was attenuated (HR = 1.50 [95% CI, 0.98-2.30]) and no longer significant when adjusting for baseline eGFR.

Conclusions and Relevance

Status as a carrier of a high-risk APOL1 genotype was associated with HFpEF hospitalization among postmenopausal women, which is partly accounted for by baseline kidney function. These findings do not support an association of high-risk APOL1 genotypes with coronary heart disease, stroke, or mortality in postmenopausal African American women.

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.

A Cross-Sectional Analysis of Telomere Length and Sleep in the Women's Health Initiative

Telomere length is a heritable marker of cellular age that is associated with morbidity and mortality. Poor sleep behaviors, which are also associated with adverse health events, may be related to leukocyte telomere length (LTL). We studied a subpopulation of 3,145 postmenopausal women (1,796 European-American (EA) and 1,349 African-American (AA)) enrolled in the Women's Health Initiative in 1993-1998 with data on Southern blot-measured LTL and self-reported usual sleep duration and sleep disturbance. LTL-sleep associations were analyzed separately for duration and disturbance using weighted and confounder-adjusted linear regression models in the entire sample (AAs + EAs; adjusted for race/ethnicity) and in racial/ethnic strata, since LTL differs by ancestry. After adjustment for covariates, each additional daily hour of sleep beyond 5 hours, approximately, was associated with a 27-base-pair (95% confidence interval (CI): 6, 48) longer LTL in the entire sample. Associations between sleep duration and LTL were strongest among AAs (adjusted β = 37, 95% CI: 4, 70); a similar, nonsignificant association was observed for EAs (adjusted β = 20, 95% CI: -7, 48). Sleep disturbance was not associated with LTL in our study. Our models did not show departure from linearity (quadratic sleep terms: P ≥ 0.55). Our results suggest that longer sleep duration is associated with longer LTL in postmenopausal women.

Epigenome-wide association study of leukocyte telomere length

Telomere length is associated with age-related diseases and is highly heritable. It is unclear, however, to what extent epigenetic modifications are associated with leukocyte telomere length (LTL). In this study, we conducted a large-scale epigenome-wide association study (EWAS) of LTL using seven large cohorts (n=5,713) - the Framingham Heart Study, the Jackson Heart Study, the Women's Health Initiative, the Bogalusa Heart Study, the Lothian Birth Cohorts of 1921 and 1936, and the Longitudinal Study of Aging Danish Twins. Our stratified analysis suggests that EWAS findings for women of African ancestry may be distinct from those of three other groups: males of African ancestry, and males and females of European ancestry. Using a meta-analysis framework, we identified DNA methylation (DNAm) levels at 823 CpG sites to be significantly associated (P<1E-7) with LTL after adjusting for age, sex, ethnicity, and imputed white blood cell counts. Functional enrichment analyses revealed that these CpG sites are near genes that play a role in circadian rhythm, blood coagulation, and wound healing. Weighted correlation network analysis identified four co-methylation modules associated with LTL, age, and blood cell counts. Overall, this study reveals highly significant relationships between two hallmarks of aging: telomere biology and epigenetic changes.