The Genetic Determinants and Genomic Consequences of Non-Leukemogenic Somatic Point Mutations.

Clonal hematopoiesis (CH) is defined by the expansion of a lineage of genetically identical cells in blood. Genetic lesions that confer a fitness advantage, such as point mutations or mosaic chromosomal alterations (mCAs) in genes associated with hematologic malignancy, are frequent mediators of CH. However, recent analyses of both single cell-derived colonies of hematopoietic cells and population sequencing cohorts have revealed CH frequently occurs in the absence of known driver genetic lesions. To characterize CH without known driver genetic lesions, we used 51,399 deeply sequenced whole genomes from the NHLBI TOPMed sequencing initiative to perform simultaneous germline and somatic mutation analyses among individuals without leukemogenic point mutations (LPM), which we term CH-LPMneg. We quantified CH by estimating the total mutation burden. Because estimating somatic mutation burden without a paired-tissue sample is challenging, we developed a novel statistical method, the Genomic and Epigenomic informed Mutation (GEM) rate, that uses external genomic and epigenomic data sources to distinguish artifactual signals from true somatic mutations. We performed a genome-wide association study of GEM to discover the germline determinants of CH-LPMneg. After fine-mapping and variant-to-gene analyses, we identified seven genes associated with CH-LPMneg (TCL1A, TERT, SMC4, NRIP1, PRDM16, MSRA, SCARB1), and one locus associated with a sex-associated mutation pathway (SRGAP2C). We performed a secondary analysis excluding individuals with mCAs, finding that the genetic architecture was largely unaffected by their inclusion. Functional analyses of SMC4 and NRIP1 implicated altered HSC self-renewal and proliferation as the primary mediator of mutation burden in blood. We then performed comprehensive multi-tissue transcriptomic analyses, finding that the expression levels of 404 genes are associated with GEM. Finally, we performed phenotypic association meta-analyses across four cohorts, finding that GEM is associated with increased white blood cell count and increased risk for incident peripheral artery disease, but is not significantly associated with incident stroke or coronary disease events. Overall, we develop GEM for quantifying mutation burden from WGS without a paired-tissue sample and use GEM to discover the genetic, genomic, and phenotypic correlates of CH-LPMneg.

Proteogenomic analysis integrated with electronic health records data reveals disease-associated variants in Black Americans.

BACKGROUND: Most genome wide association studies (GWAS) of plasma proteomics have focused on White individuals of European ancestry, limiting biological insight from other ancestry enriched protein quantitative loci (pQTL). METHODS: We conducted a discovery GWAS of ~3,000 plasma proteins measured by the antibody based Olink platform in 1,054 Black adults from the Jackson Heart Study (JHS), and validated our findings in the Multi-Ethnic Study of Atherosclerosis (MESA). The genetic architecture of identified pQTLs were further explored through fine mapping and admixture association analysis. Finally, using our pQTL findings, we performed a phenome wide association study (PheWAS) across two large multi-ethnic electronic health record (EHR) systems in All of Us and BioMe. RESULTS: We identified 1002 pQTLs for 925 proteins. Fine mapping and admixture analyses suggested allelic heterogeneity of the plasma proteome across diverse populations. We identified associations for variants enriched in African ancestry, many in diseases that lack precise biomarkers, including cis-pQTLs for Cathepsin L (CTSL) and Siglec-9 that were linked with sarcoidosis and non-Hodgkin's lymphoma, respectively. We found concordant associations across clinical diagnoses and laboratory measurements, elucidating disease pathways, including a cis-pQTL associated with circulating CD58, white blood cell count, and multiple sclerosis. CONCLUSIONS: Our findings emphasize the value of leveraging diverse populations to enhance biological insights from proteomics GWAS, and we have made this resource readily available as an interactive web portal.

Determinants of mosaic chromosomal alteration fitness.

Clonal hematopoiesis (CH) is characterized by the acquisition of a somatic mutation in a hematopoietic stem cell that results in a clonal expansion. These driver mutations can be single nucleotide variants in cancer driver genes or larger structural rearrangements called mosaic chromosomal alterations (mCAs). The factors that influence the variations in mCA fitness and ultimately result in different clonal expansion rates are not well understood. We used the Passenger-Approximated Clonal Expansion Rate (PACER) method to estimate clonal expansion rate as PACER scores for 6,381 individuals in the NHLBI TOPMed cohort with gain, loss, and copy-neutral loss of heterozygosity mCAs. Our mCA fitness estimates, derived by aggregating per-individual PACER scores, were correlated (R2 = 0.49) with an alternative approach that estimated fitness of mCAs in the UK Biobank using population-level distributions of clonal fraction. Among individuals with JAK2 V617F clonal hematopoiesis of indeterminate potential or mCAs affecting the JAK2 gene on chromosome 9, PACER score was strongly correlated with erythrocyte count. In a cross-sectional analysis, genome-wide association study of estimates of mCA expansion rate identified a TCL1A locus variant associated with mCA clonal expansion rate, with suggestive variants in NRIP1 and TERT.

Methylation patterns associated with C-reactive protein in racially and ethnically diverse populations.

Systemic low-grade inflammation is a feature of chronic disease. C-reactive protein (CRP) is a common biomarker of inflammation and used as an indicator of disease risk; however, the role of inflammation in disease is not completely understood. Methylation is an epigenetic modification in the DNA which plays a pivotal role in gene expression. In this study we evaluated differential DNA methylation patterns associated with blood CRP level to elucidate biological pathways and genetic regulatory mechanisms to improve the understanding of chronic inflammation. The racially and ethnically diverse participants in this study were included as 50% White, 41% Black or African American, 7% Hispanic or Latino/a, and 2% Native Hawaiian, Asian American, American Indian, or Alaska Native (total n = 13,433) individuals. We replicated 113 CpG sites from 87 unique loci, of which five were novel (CADM3, NALCN, NLRC5, ZNF792, and cg03282312), across a discovery set of 1,150 CpG sites associated with CRP level (p < 1.2E-7). The downstream pathways affected by DNA methylation included the identification of IFI16 and IRF7 CpG-gene transcript pairs which contributed to the innate immune response gene enrichment pathway along with NLRC5, NOD2, and AIM2. Gene enrichment analysis also identified the nuclear factor-kappaB transcription pathway. Using two-sample Mendelian randomization (MR) we inferred methylation at three CpG sites as causal for CRP levels using both White and Black or African American MR instrument variables. Overall, we identified novel CpG sites and gene transcripts that could be valuable in understanding the specific cellular processes and pathogenic mechanisms involved in inflammation.

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.

ABSTRACT: 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 3 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. 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 7 new loci for FVIII (ST3GAL4, CLEC4M, B3GNT2, ASGR1, F12, KNG1, and TREM1/NCR2) and 1 for VWF (B3GNT2). VWF, ABO, and STAB2 were associated with FVIII and VWF in gene-based analyses. Multiphenotype analysis of FVIII and VWF identified another 3 new loci, including PDIA3. Silencing of B3GNT2 and the previously reported CD36 gene decreased release of FVIII by HLECs, whereas 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 1 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.

Individual and Neighborhood-level Socioeconomic Status and Somatic Mutations Associated With Increased Risk of Cardiovascular Disease and Mortality: A Cross-Sectional Analysis in the Women's Health Initiative.

BACKGROUND: Clonal hematopoiesis of indeterminate potential (CHIP), the expansion of leukemogenic mutations in white blood cells, has been associated with increased risk of atherosclerotic cardiovascular diseases, cancer, and mortality. OBJECTIVE: We examined the relationship between individual- and neighborhood-level socioeconomic status (SES) and CHIP and evaluated effect modification by interpersonal and intrapersonal resources. METHODS: The study population included 10,799 postmenopausal women from the Women's Health Initiative without hematologic malignancy or antineoplastic medication use. Individual- and neighborhood (Census tract)-level SES were assessed across several domains including education, income, and occupation, and a neighborhood-level SES summary z-score, which captures multiple dimensions of SES, was generated. Interpersonal and intrapersonal resources were self-reports. CHIP was ascertained based on a prespecified list of leukemogenic driver mutations. Weighted logistic regression models adjusted for covariates were used to estimate risk of CHIP as an odds ratio (OR) and 95% confidence interval (95% CI). RESULTS: The interval-scale neighborhood-level SES summary z-score was associated with a 3% increased risk of CHIP: OR (95% CI) = 1.03 (1.00-1.05), p = .038. Optimism significantly modified that estimate, such that among women with low/medium and high levels of optimism, the corresponding ORs (95% CIs) were 1.03 (1.02-1.04) and 0.95 (0.94-0.96), pInteraction < .001. CONCLUSIONS: Our findings suggest that reduced risk of somatic mutation may represent a biological pathway by which optimism protects contextually advantaged but at-risk women against age-related chronic disease and highlight potential benefits of long-term, positive psychological interventions.

Determinants of mosaic chromosomal alteration fitness.

Clonal hematopoiesis (CH) is characterized by the acquisition of a somatic mutation in a hematopoietic stem cell that results in a clonal expansion. These driver mutations can be single nucleotide variants in cancer driver genes or larger structural rearrangements called mosaic chromosomal alterations (mCAs). The factors that influence the variations in mCA fitness and ultimately result in different clonal expansion rates are not well-understood. We used the Passenger-Approximated Clonal Expansion Rate (PACER) method to estimate clonal expansion rate for 6,381 individuals in the NHLBI TOPMed cohort with gain, loss, and copy-neutral loss of heterozygosity mCAs. Our estimates of mCA fitness were correlated (R 2 = 0.49) with an alternative approach that estimated fitness of mCAs in the UK Biobank using a theoretical probability distribution. Individuals with lymphoid-associated mCAs had a significantly higher white blood cell count and faster clonal expansion rate. In a cross-sectional analysis, genome-wide association study of estimates of mCA expansion rate identified TCL1A , NRIP1 , and TERT locus variants as modulators of mCA clonal expansion rate.

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.

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.

The Relationship of Duffy Gene Polymorphism, High Sensitivity C-Reactive Protein, and Long-term Outcomes.

Background: Black adults have higher incidence of all-cause death and worse cardiovascular outcomes when compared to other populations. The Duffy chemokine receptor is not expressed in a large majority of Black adults and the clinical implications of this are unclear. Methods: Here, we investigated the relationship of Duffy receptor status, high-sensitivity C-reactive protein (hs-CRP), and long-term cardiovascular outcomes in Black members of two contemporary, longitudinal cohort studies (the Jackson Heart Study and Multi-Ethnic Study of Atherosclerosis). Data on 4,307 Black participants (2,942 Duffy null and 1,365 Duffy receptor positive, as defined using Single Nucleotide Polymorphism (SNP) rs2814778) were included in this analysis. Results: Duffy null was not independently associated with elevated levels of serum hs-CRP levels once conditioning for known CRP locus alleles in linkage disequilibrium with the Duffy gene. Duffy null status was not found to be independently associated with higher incidence of all-cause mortality or secondary outcomes after adjusting for possible confounders in Black participants. Conclusions: These findings suggest that increased levels of hs-CRP found in Duffy null individuals is due to co-inheritance of CRP alleles known to influence circulating levels hs-CRP and that Duffy null status was not associated with worse adverse outcomes over the follow-up period in this cohort of well-balanced Black participants.

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.

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.

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.

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.

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.