Metabolomic profiling identifies novel metabolites associated with cardiac dysfunction.

Metabolic comorbidities, such as obesity and diabetes, are associated with subclinical alterations in both cardiac structure/function and natriuretic peptides prior to the onset of heart failure (HF). Despite this, the exact metabolic pathways of cardiac dysfunction which precede HF are not well-defined. Among older individuals without HF in the Multi-Ethnic Study of Atherosclerosis (MESA), we evaluated the associations of 47 circulating metabolites measured by 1H-NMR with echocardiographic measures of cardiac structure and function. We then evaluated associations of significant metabolites with circulating N-terminal pro-B-type natriuretic peptide (NT-proBNP). In a separate cohort, we evaluated differences between top metabolites in patients with HF with preserved ejection fraction (HFpEF) and comorbidity-matched controls. Genetic variants associated with top metabolites (mQTLs) were then related to echocardiographic measures and NT-proBNP. Among 3440 individuals with metabolic and echocardiographic data in MESA (62 ± 10 years, 52% female, 38% White), 10 metabolites broadly reflective of glucose and amino acid metabolism were associated with at least 1 measure of cardiac structure or function. Of these 10 metabolites, 4 (myo-inositol, glucose, dimethylsulfone, carnitine) were associated with higher NT-proBNP and 2 (d-mannose, acetone) were associated with lower NT-proBNP. In a separate cohort, patients with HFpEF had higher circulating myo-inositol levels compared with comorbidity-matched controls. Genetic analyses revealed that 1 of 6 known myo-inositol mQTLs conferred risk of higher NT-proBNP. In conclusion, metabolomic profiling identifies several novel metabolites associated with cardiac dysfunction in a cohort at high risk for HF, revealing pathways potentially relevant to future HF risk.

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.

Identification of proteins associated with type 2 diabetes risk in diverse racial and ethnic populations.

AIMS/HYPOTHESIS: Several studies have reported associations between specific proteins and type 2 diabetes risk in European populations. To better understand the role played by proteins in type 2 diabetes aetiology across diverse populations, we conducted a large proteome-wide association study using genetic instruments across four racial and ethnic groups: African; Asian; Hispanic/Latino; and European. METHODS: Genome and plasma proteome data from the Multi-Ethnic Study of Atherosclerosis (MESA) study involving 182 African, 69 Asian, 284 Hispanic/Latino and 409 European individuals residing in the USA were used to establish protein prediction models by using potentially associated cis- and trans-SNPs. The models were applied to genome-wide association study summary statistics of 250,127 type 2 diabetes cases and 1,222,941 controls from different racial and ethnic populations. RESULTS: We identified three, 44 and one protein associated with type 2 diabetes risk in Asian, European and Hispanic/Latino populations, respectively. Meta-analysis identified 40 proteins associated with type 2 diabetes risk across the populations, including well-established as well as novel proteins not yet implicated in type 2 diabetes development. CONCLUSIONS/INTERPRETATION: Our study improves our understanding of the aetiology of type 2 diabetes in diverse populations. DATA AVAILABILITY: The summary statistics of multi-ethnic type 2 diabetes GWAS of MVP, DIAMANTE, Biobank Japan and other studies are available from The database of Genotypes and Phenotypes (dbGaP) under accession number phs001672.v3.p1. MESA genetic, proteome and covariate data can be accessed through dbGaP under phs000209.v13.p3. All code is available on GitHub ( https://github.com/Arthur1021/MESA-1K-PWAS ).

Obstructive sleep apnea mediates genetic risk of Diabetes Mellitus: The Hispanic Community Health Study/Study of Latinos.

Objective: We sought to evaluate whether obstructive sleep apnea (OSA), and other sleep disorders, increase genetic risk of developing diabetes mellitus (DM). Research Design and Methods: Using GWAS summary statistics from the DIAGRAM consortium and Million Veteran Program, we developed multi-ancestry Type 2 Diabetes (T2D) polygenic risk scores (T2D-PRSs) useful in admixed Hispanic/Latino individuals. We estimated the association of the T2D-PRS with cross-sectional and incident DM in the Hispanic Community Health Study/Study of Latinos (HCHS/SOL). We conducted a mediation analysis with T2D-PRSs as an exposure, incident DM as an outcome, and OSA as a mediator. Additionally, we performed Mendelian randomization (MR) analysis to assess the causal relationship between T2D and OSA. Results: Of 12,342 HCHS/SOL participants, at baseline, 48.4% were normoglycemic, 36.6% were hyperglycemic, and 15% had diabetes, and 50.9% identified as female. Mean age was 41.5, and mean BMI was 29.4. T2D-PRSs was strongly associated with baseline DM and with incident DM. At baseline, a 1 SD increase in the primary T2D-PRS had DM adjusted odds ratio (OR) = 2.67, 95% CI [2.40; 2.97] and a higher incident DM rate (incident rate ratio (IRR) = 2.02, 95% CI [1.75; 2.33]). In a stratified analysis based on OSA severity categories the associations were stronger in individuals with mild OSA compared to those with moderate to severe OSA. Mediation analysis suggested that OSA mediates the T2D-PRS association with DM. In two-sample MR analysis, T2D-PRS had a causal effect on OSA, OR = 1.03, 95% CI [1.01; 1.05], and OSA had a causal effect on T2D, with OR = 2.34, 95% CI [1.59; 3.44]. Conclusions: OSA likely mediates genetic effects on T2D.

Steroid Hormone Biosynthesis and Dietary Related Metabolites Associated with Excessive Daytime Sleepiness.

Background: Excessive daytime sleepiness (EDS) is a complex sleep problem that affects approximately 33% of the United States population. Although EDS usually occurs in conjunction with insufficient sleep, and other sleep and circadian disorders, recent studies have shown unique genetic markers and metabolic pathways underlying EDS. Here, we aimed to further elucidate the biological profile of EDS using large scale single- and pathway-level metabolomics analyses. Methods: Metabolomics data were available for 877 metabolites in 6,071 individuals from the Hispanic Community Health Study/Study of Latinos (HCHS/SOL) and EDS was assessed using the Epworth Sleepiness Scale (ESS) questionnaire. We performed linear regression for each metabolite on continuous ESS, adjusting for demographic, lifestyle, and physiological confounders, and in sex specific groups. Subsequently, gaussian graphical modelling was performed coupled with pathway and enrichment analyses to generate a holistic interactive network of the metabolomic profile of EDS associations. Findings: We identified seven metabolites belonging to steroids, sphingomyelin, and long chain fatty acids sub-pathways in the primary model associated with EDS, and an additional three metabolites in the male-specific analysis. The identified metabolites particularly played a role in steroid hormone biosynthesis. Interpretation: Our findings indicate that an EDS metabolomic profile is characterized by endogenous and dietary metabolites within the steroid hormone biosynthesis pathway, with some pathways that differ by sex. Our findings identify potential pathways to target for addressing the causes or consequences of EDS and related sleep disorders. Funding: Details regarding funding supporting this work and all studies involved are provided in the acknowledgments section.

Analysis of C-reactive protein omics-measures associates methylation risk score with sleep health and related health outcomes.

Introduction: DNA methylation (DNAm) predictors of high sensitivity C-reactive protein (CRP) offer a stable and accurate means of assessing chronic inflammation, bypassing the CRP protein fluctuations secondary to acute illness. Poor sleep health is associated with elevated inflammation (including elevated blood CRP levels) which may explain associations of sleep insufficiency with metabolic, cardiovascular and neurological diseases. Our study aims to characterize the relationships among sleep health phenotypes and CRP markers -blood, genetic, and epigenetic indicators-within the Hispanic Community Health Study/Study of Latinos (HCHS/SOL). Methods: In HCHS/SOL, methylation risk scores (MRS)-CRP and polygenetic risk score (PRS)-CRP were constructed separately as weighted sums of methylation beta values or allele counts, respectively, for each individual. Sleep health phenotypes were measured using self-reported questionnaires and objective measurements. Survey-weighted linear regression established the association between the multiple sleep phenotypes (obstructive sleep apnea (OSA), sleep duration, insomnia and excessive sleepiness symptom), cognitive assessments, diabetes and hypertension with CRP markers while adjusting for age, sex, BMI, study center, and the first five principal components of genetic ancestry in HCHS/SOL. Results: We included 2221 HCHS/SOL participants (age range 37-76 yrs, 65.7% female) in the analysis. Both the MRS-CRP (95% confidence interval (CI): 0.32-0.42, p = 3.3 × 10-38) and the PRS-CRP (95% CI: 0.15-0.25, p = 1 × 10-14) were associated with blood CRP level. Moreover, MRS-CRP was associated with sleep health phenotypes (OSA, long sleep duration) and related conditions (diabetes and hypertension), while PRS-CRP markers were not associated with these traits. Circulating CRP level was associated with sleep duration and diabetes. Associations between OSA traits and metabolic comorbidities weakened after adjusting for MRS-CRP, most strongly for diabetes, and least for hypertension. Conclusions: MRS-CRP is a promising estimate for systemic and chronic inflammation as reflected by circulating CRP levels, which either mediates or serves as a common cause of the association between sleep phenotypes and related comorbidities, especially in the presence of diabetes.

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.

Proteomics of left ventricular structure in the Multi-Ethnic Study of Atherosclerosis.

AIMS: Proteomic profiling offers an expansive approach to biomarker discovery and mechanistic hypothesis generation for LV remodelling, a critical component of heart failure (HF). We sought to identify plasma proteins cross-sectionally associated with left ventricular (LV) size and geometry in a diverse population-based cohort without known cardiovascular disease (CVD). METHODS AND RESULTS: Among participants of the Multi-Ethnic Study of Atherosclerosis (MESA), we quantified plasma abundances of 1305 proteins using an aptamer-based platform at exam 1 (2000-2002) and exam 5 (2010-2011) and assessed LV structure by cardiac magnetic resonance (CMR) at the same time points. We used multivariable linear regression with robust variance to assess cross-sectional associations between plasma protein abundances and LV structural characteristics at exam 1, reproduced findings in later-life at exam 5, and explored relationships of associated proteins using annotated enrichment analysis. We studied 763 participants (mean age 60 ± 10 years at exam 1; 53% female; 19% Black race; 31% Hispanic ethnicity). Following adjustment for renal function and traditional CVD risk factors, plasma levels of 3 proteins were associated with LV mass index at both time points with the same directionality (FDR < 0.05): leptin (LEP), renin (REN), and cathepsin-D (CTSD); 20 with LV end-diastolic volume index: LEP, NT-proBNP, histone-lysine N-methyltransferase (EHMT2), chordin-like protein 1 (CHRDL1), tumour necrosis factor-inducible gene 6 protein (TNFAIP6), NT-3 growth factor receptor (NTRK3), c5a anaphylatoxin (C5), neurogenic locus notch homologue protein 3 (NOTCH3), ephrin-B2 (EFNB2), osteomodulin (OMD), contactin-4 (CNTN4), gelsolin (GSN), stromal cell-derived factor 1 (CXCL12), calcineurin subunit B type 1 (PPP3R1), insulin-like growth factor 1 receptor (IGF1R), bone sialoprotein 2 (IBSP), interleukin-11 (IL-11), follistatin-related protein 1 (FSTL1), periostin (POSTN), and biglycan (BGN); and 4 with LV mass-to-volume ratio: RGM domain family member B (RGMB), transforming growth factor beta receptor type 3 (TGFBR3), ephrin-A2 (EFNA2), and cell adhesion molecule 3 (CADM3). Functional annotation implicated regulation of the PI3K-Akt pathway, bone morphogenic protein signalling, and cGMP-mediated signalling. CONCLUSIONS: We report proteomic profiling of LV size and geometry, which identified novel associations and reinforced previous findings on biomarker candidates for LV remodelling and HF. If validated, these proteins may help refine risk prediction and identify novel therapeutic targets for HF.

Metagenomic Study of the MESA: Detection of Gemella Morbillorum and Association With Coronary Heart Disease.

BACKGROUND: Inflammation is a feature of coronary heart disease (CHD), but the role of proinflammatory microbial infection in CHD remains understudied. METHODS AND RESULTS: CHD was defined in the MESA (Multi-Ethnic Study of Atherosclerosis) as myocardial infarction (251 participants), resuscitated arrest (2 participants), and CHD death (80 participants). We analyzed sequencing reads from 4421 MESA participants in the National Heart, Lung, and Blood Institute Trans-Omics for Precision Medicine program using the PathSeq workflow of the Genome Analysis Tool Kit and a 65-gigabase microbial reference. Paired reads aligning to 840 microbes were detected in >1% of participants. The association of the presence of microbe reads with incident CHD (follow-up, ~18 years) was examined. First, important variables were ascertained using a single regularized Cox proportional hazard model, examining change of risk as a function of presence of microbe with age, sex, education level, Life's Simple 7, and inflammation. For variables of importance, the hazard ratio (HR) was estimated in separate (unregularized) Cox proportional hazard models including the same covariates (significance threshold Bonferroni corrected P<6×10-5, 0.05/840). Reads from 2 microbes were significantly associated with CHD: Gemella morbillorum (HR, 3.14 [95% CI, 1.92-5.12]; P=4.86×10-6) and Pseudomonas species NFACC19-2 (HR, 3.22 [95% CI, 2.03-5.41]; P=1.58×10-6). CONCLUSIONS: Metagenomics of whole-genome sequence reads opens a possible frontier for detection of pathogens for chronic diseases. The association of G morbillorum and Pseudomonas species reads with CHD raises the possibilities that microbes may drive atherosclerotic inflammation and that treatments for specific pathogens may provide clinical utility for CHD reduction.

Sexually dimorphic DNA methylation and gene expression patterns in human first trimester placenta.

BACKGROUND: Fetal sex and placental development impact pregnancy outcomes and fetal-maternal health, but the critical timepoint of placenta establishment in first trimester is understudied in human pregnancies. METHODS: Pregnant subjects were recruited in late first trimester (weeks 10-14) at time of chorionic villus sampling, a prenatal diagnostic test. Leftover placenta tissue was collected and stored until birth outcomes were known, then DNA and RNA were isolated from singleton, normal karyotype pregnancies resulting in live births. DNA methylation was measured with the Illumina Infinium MethylationEPIC BeadChip array (n = 56). Differential methylation analysis compared 25 females versus 31 males using a generalized linear model on 743,461 autosomal probes. Gene expression sex differences were analyzed with RNA-sequencing (n = 74). An integrated analysis was performed using linear regression to correlate gene expression and DNA methylation in 51 overlapping placentas. RESULTS: Methylation analysis identified 151 differentially methylated probes (DMPs) significant at false discovery rate < 0.05, including 89 (59%) hypermethylated in females. Probe cg17612569 (GABPA, ATP5J) was the most significant CpG site, hypermethylated in males. There were 11 differentially methylated regions affected by fetal sex, with transcription factors ZNF300 and ZNF311 most significantly hypermethylated in males and females, respectively. RNA-sequencing identified 152 genes significantly sexually dimorphic at false discovery rate < 0.05. The 151 DMPs were associated with 18 genes with gene downregulation (P < 0.05) in the direction of hypermethylation, including 2 genes significant at false discovery rate < 0.05 (ZNF300 and CUB and Sushi multiple domains 1, CSMD1). Both genes, as well as Family With Sequence Similarity 228 Member A (FAM228A), showed significant correlation between DNA methylation and sexually dimorphic gene expression, though FAM228A DNA methylation was less sexually dimorphic. Comparison with other sex differences studies found that cg17612569 is male-hypermethylated across gestation in placenta and in human blood up to adulthood. CONCLUSIONS: Overall, sex dimorphic differential methylation with associated differential gene expression in the first trimester placenta is small, but there remain significant genes that may be regulated through methylation leading to differences in the first trimester placenta.

Fetal sex and placenta development affect pregnancy outcomes for both the fetus and mother throughout pregnancy, including risk of miscarriages, preterm birth, preeclampsia, and other outcomes. Epigenetics, the "overlay" of regulatory signals on DNA which affects how DNA is read, is not well understood in early pregnancy when critical placenta developments are happening that affect the rest of pregnancy. Here, we use leftover placenta biopsy samples (n = 56) donated by Cedars-Sinai patients with informed consent to learn about first trimester human placenta DNA methylation differences due to fetal sex. Out of the total 743,461 sites analyzed, we identified 151 sites significantly affected by fetal sex after correcting p-values to reduce false positives (false discovery rate < 0.05). We also performed an analysis to look at multiple sites and identified 11 regions across the genome with significant DNA methylation changes due to fetal sex. Furthermore, because DNA methylation is a regulatory mark on DNA which typically dampens gene expression, we also compared the DNA methylation sex differences to placental RNA-sequencing gene expression analysis using the same tissue from a mostly overlapping patient group (n = 74 total sequenced, n = 51 overlap). We identify 18 genes which show both significant DNA methylation differences and gene expression changes. The most significant gene was transcription factor ZNF300 with higher DNA methylation in males and reduced gene expression in males (and thus higher gene expression in females). This study identifies some sex differences that continue until later pregnancy and others that are unique to first trimester.

A methylation risk score for chronic kidney disease: a HyperGEN study.

Chronic kidney disease (CKD) impacts about 1 in 7 adults in the United States, but African Americans (AAs) carry a disproportionately higher burden of disease. Epigenetic modifications, such as DNA methylation at cytosine-phosphate-guanine (CpG) sites, have been linked to kidney function and may have clinical utility in predicting the risk of CKD. Given the dynamic relationship between the epigenome, environment, and disease, AAs may be especially sensitive to environment-driven methylation alterations. Moreover, risk models incorporating CpG methylation have been shown to predict disease across multiple racial groups. In this study, we developed a methylation risk score (MRS) for CKD in cohorts of AAs. We selected nine CpG sites that were previously reported to be associated with estimated glomerular filtration rate (eGFR) in epigenome-wide association studies to construct a MRS in the Hypertension Genetic Epidemiology Network (HyperGEN). In logistic mixed models, the MRS was significantly associated with prevalent CKD and was robust to multiple sensitivity analyses, including CKD risk factors. There was modest replication in validation cohorts. In summary, we demonstrated that an eGFR-based CpG score is an independent predictor of prevalent CKD, suggesting that MRS should be further investigated for clinical utility in evaluating CKD risk and progression.

Intracranial Volume Is Driven by Both Genetics and Early Life Exposures: The SOL-INCA-MRI Study.

Intracranial volume (ICV) reflects maximal brain development and is associated with later-life cognitive abilities. We quantified ICV among first- and second-generation Hispanic and Latino adults from the Study of Latinos-Investigation of Cognitive Aging - MRI (SOL-INCA-MRI), estimated ICV heritability, and tested its associations with previously reported genetic variants, both individually and as a genetic risk score (GRS). We also estimated the association of ICV with early life environmental measures: nativity or age of immigration and parental education. The estimated heritability of ICV was 19% (95% CI, 0.1%-56%) in n=1781 unrelated SOL-INCA-MRI individuals. Four of 10 tested genetic variants were associated with ICV and an increase of 1 SD of the ICV-GRS was associated with an increase of 10.37 cm3 in the ICV (95% CI, 5.29-15.45). Compared to being born in the continental United States, immigrating to the United States at age 11 years or older was associated with 24 cm3 smaller ICV (95% CI, -39.97 to -8.06). Compared to both parents having less than high-school education, at least 1 parent completing high-school education was associated with 15.4 cm3 greater ICV (95% CI, 4.46-26.39). These data confirm the importance of early life health on brain development.

Immune cell profiling of the ICAM1 p.K56M heart failure with preserved ejection fraction risk variant.

AIMS: Intercellular adhesion molecule-1 (ICAM-1) facilitates inflammation via leucocyte recruitment and has been implicated in heart failure with preserved ejection fraction (HFpEF). Approximately 35% of African American individuals carry a copy of an ICAM1 missense variant (rs5491; p.K56M), which is associated with an increased risk of HFpEF. The pathways by which rs5491 increases HFpEF risk are not well defined. We evaluated the circulating immune cell profile of rs5491. METHODS: Among African American individuals in the Multi-Ethnic Study of Atherosclerosis, we evaluated the associations of rs5491 with 29 circulating peripheral blood mononuclear cell subsets. The top immune cells were then related to echocardiographic measures of structure and function. RESULTS: Among 502 individuals with immune cell profiling (mean age 63 years, 51% female), 191 individuals (38%) had at least one copy of rs5491. Each additional rs5491 allele was significantly associated with higher proportions of Tc17 CD8+ cytotoxic T cells (ß = 1.34, SE = 0.45, P = 9.5 × 10-5) and Tc2 CD8+ cytotoxic T cells (ß = 1.19, SE = 0.44, P = 0.00012). There were no other associations noted between rs5491 and the remaining immune cells. A higher proportion of Tc17 cells was significantly associated with a higher left ventricular ejection fraction, E/e' average and right ventricular systolic pressure (RVSP), while a higher proportion of Tc2 cells was significantly associated with a higher RVSP. CONCLUSIONS: The ICAM1 p.K56M variant (rs5491) carries a distinct and inflammatory T-cell subset profile. These cytotoxic T cells are in turn associated with alterations in cardiac function and adverse haemodynamics later in life, thus providing insight into pathways by which rs5491 may increase the risk of HFpEF.

Plasma Proteomics of Exercise Blood Pressure and Incident Hypertension.

Importance: Blood pressure response during acute exercise (exercise blood pressure [EBP]) is associated with the future risk of hypertension and cardiovascular disease (CVD). Biochemical characterization of EBP could inform disease biology and identify novel biomarkers of future hypertension. Objective: To identify protein markers associated with EBP and test their association with incident hypertension. Design, Setting, and Participants: This study assayed 4977 plasma proteins in 681 healthy participants (from 763 assessed) of the Health, Risk Factors, Exercise Training and Genetics (HERITAGE; data collection from January 1993 to December 1997 and plasma proteomics from January 2019 to January 2020) Family Study at rest who underwent 2 cardiopulmonary exercise tests. Individuals were free of CVD at the time of recruitment. Individuals with resting SBP ≥160 mm Hg or DBP ≥100 mm Hg or taking antihypertensive drug therapy were excluded from the study. The association between resting plasma protein levels to both resting BP and EBP was evaluated. Proteins associated with EBP were analyzed for their association with incident hypertension in the Framingham Heart Study (FHS; n = 1177) and validated in the Jackson Heart Study (JHS; n = 772) and Multi-Ethnic Study of Atherosclerosis (MESA; n = 1367). Proteins associated with incident hypertension were tested for putative causal links in approximately 700 000 individuals using cis-protein quantitative loci mendelian randomization (cis-MR). Data were analyzed from January 2023 to January 2024. Exposures: Plasma proteins. Main Outcomes and Measures: EBP was defined as the BP response during a fixed workload (50 W) on a cycle ergometer. Hypertension was defined as BP ≥140/90 mm Hg or taking antihypertensive medication. Results: Among the 681 participants in the HERITAGE Family Study, the mean (SD) age was 34 (13) years; 366 participants (54%) were female; 238 (35%) were self-reported Black and 443 (65%) were self-reported White. Proteomic profiling of EBP revealed 34 proteins that would not have otherwise been identified through profiling of resting BP alone. Transforming growth factor ß receptor 3 (TGFBR3) and prostaglandin D2 synthase (PTGDS) had the strongest association with exercise systolic BP (SBP) and diastolic BP (DBP), respectively (TGFBR3: exercise SBP, ß estimate, -3.39; 95% CI, -4.79 to -2.00; P = 2.33 × 10-6; PTGDS: exercise DBP ß estimate, -2.50; 95% CI, -3.29 to -1.70; P = 1.18 × 10-9). In fully adjusted models, TGFBR3 was inversely associated with incident hypertension in FHS, JHS, and MESA (hazard ratio [HR]: FHS, 0.86; 95% CI, 0.75-0.97; P = .01; JHS, 0.87; 95% CI, 0.77-0.97; P = .02; MESA, 0.84; 95% CI, 0.71-0.98; P = .03; pooled cohort, 0.86; 95% CI, 0.79-0.92; P = 6 × 10-5). Using cis-MR, genetically predicted levels of TGFBR3 were associated with SBP, hypertension, and CVD events (SBP: ß, -0.38; 95% CI, -0.64 to -0.11; P = .006; hypertension: odds ratio [OR], 0.99; 95% CI, 0.98-0.99; P < .001; heart failure with hypertension: OR, 0.86; 95% CI, 0.77-0.97; P = .01; CVD: OR, 0.84; 95% CI, 0.77-0.92; P = 8 × 10-5; cerebrovascular events: OR, 0.77; 95% CI, 0.70-0.85; P = 5 × 10-7). Conclusions and Relevance: Plasma proteomic profiling of EBP identified a novel protein, TGFBR3, which may protect against elevated BP and long-term CVD outcomes.

Epigenetic and proteomic signatures associate with clonal hematopoiesis expansion rate.

Clonal hematopoiesis of indeterminate potential (CHIP), whereby somatic mutations in hematopoietic stem cells confer a selective advantage and drive clonal expansion, not only correlates with age but also confers increased risk of morbidity and mortality. Here, we leverage genetically predicted traits to identify factors that determine CHIP clonal expansion rate. We used the passenger-approximated clonal expansion rate method to quantify the clonal expansion rate for 4,370 individuals in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) cohort and calculated polygenic risk scores for DNA methylation aging, inflammation-related measures and circulating protein levels. Clonal expansion rate was significantly associated with both genetically predicted and measured epigenetic clocks. No associations were identified with inflammation-related lab values or diseases and CHIP expansion rate overall. A proteome-wide search identified predicted circulating levels of myeloid zinc finger 1 and anti-Müllerian hormone as associated with an increased CHIP clonal expansion rate and tissue inhibitor of metalloproteinase 1 and glycine N-methyltransferase as associated with decreased CHIP clonal expansion rate. Together, our findings identify epigenetic and proteomic patterns associated with the rate of hematopoietic clonal expansion.

Validation of human telomere length multi-ancestry meta-analysis association signals identifies POP5 and KBTBD6 as human telomere length regulation genes.

Genome-wide association studies (GWAS) have become well-powered to detect loci associated with telomere length. However, no prior work has validated genes nominated by GWAS to examine their role in telomere length regulation. We conducted a multi-ancestry meta-analysis of 211,369 individuals and identified five novel association signals. Enrichment analyses of chromatin state and cell-type heritability suggested that blood/immune cells are the most relevant cell type to examine telomere length association signals. We validated specific GWAS associations by overexpressing KBTBD6 or POP5 and demonstrated that both lengthened telomeres. CRISPR/Cas9 deletion of the predicted causal regions in K562 blood cells reduced expression of these genes, demonstrating that these loci are related to transcriptional regulation of KBTBD6 and POP5. Our results demonstrate the utility of telomere length GWAS in the identification of telomere length regulation mechanisms and validate KBTBD6 and POP5 as genes affecting telomere length regulation.

Genetic Subtypes of Prediabetes, Healthy Lifestyle, and Risk of Type 2 Diabetes.

Prediabetes is a heterogenous metabolic state with various risks for development of type 2 diabetes (T2D). In this study, we used genetic data on 7,227 US Hispanic/Latino participants without diabetes from the Hispanic Community Health Study/Study of Latinos (HCHS/SOL) and 400,149 non-Hispanic White participants without diabetes from the UK Biobank (UKBB) to calculate five partitioned polygenetic risk scores (pPRSs) representing various pathways related to T2D. Consensus clustering was performed in participants with prediabetes in HCHS/SOL (n = 3,677) and UKBB (n = 16,284) separately based on these pPRSs. Six clusters of individuals with prediabetes with distinctive patterns of pPRSs and corresponding metabolic traits were identified in the HCHS/SOL, five of which were confirmed in the UKBB. Although baseline glycemic traits were similar across clusters, individuals in cluster 5 and cluster 6 showed an elevated risk of T2D during follow-up compared with cluster 1 (risk ratios [RRs] 1.29 [95% CI 1.08, 1.53] and 1.34 [1.13, 1.60], respectively). Inverse associations between a healthy lifestyle score and risk of T2D were observed across different clusters, with a suggestively stronger association observed in cluster 5 compared with cluster 1. Among individuals with a healthy lifestyle, those in cluster 5 had a similar risk of T2D compared with those in cluster 1 (RR 1.03 [0.91, 1.18]). This study identified genetic subtypes of prediabetes that differed in risk of progression to T2D and in benefits from a healthy lifestyle.

Longitudinal change in mitochondrial heteroplasmy exhibits positive selection for deleterious variants.

A common feature of human aging is the acquisition of somatic mutations, and mitochondria are particularly prone to mutation due to their inefficient DNA repair and close proximity to reactive oxygen species, leading to a state of mitochondrial DNA heteroplasmy1,2. Cross-sectional studies have demonstrated that detection of heteroplasmy increases with participant age3, a phenomenon that has been attributed to genetic drift4-7. In this first large-scale longitudinal study, we measured heteroplasmy in two prospective cohorts (combined n=1405) at two timepoints (mean time between visits, 8.6 years), demonstrating that deleterious heteroplasmies were more likely to increase in variant allele fraction (VAF). We further demonstrated that increase in VAF was associated with increased risk of overall mortality. These results challenge the claim that somatic mtDNA mutations arise mainly due to genetic drift, instead demonstrating positive selection for predicted deleterious mutations at the cellular level, despite an negative impact on overall mortality.

Time-to-Event Genome-Wide Association Study for Incident Cardiovascular Disease in People With Type 2 Diabetes.

OBJECTIVE: To identify genetic risk factors for incident cardiovascular disease (CVD) among people with type 2 diabetes (T2D). RESEARCH DESIGN AND METHODS: We conducted a multiancestry time-to-event genome-wide association study for incident CVD among people with T2D. We also tested 204 known coronary artery disease (CAD) variants for association with incident CVD. RESULTS: Among 49,230 participants with T2D, 8,956 had incident CVD events (event rate 18.2%). We identified three novel genetic loci for incident CVD: rs147138607 (near CACNA1E/ZNF648, hazard ratio [HR] 1.23, P = 3.6 × 10-9), rs77142250 (near HS3ST1, HR 1.89, P = 9.9 × 10-9), and rs335407 (near TFB1M/NOX3, HR 1.25, P = 1.5 × 10-8). Among 204 known CAD loci, 5 were associated with incident CVD in T2D (multiple comparison-adjusted P < 0.00024, 0.05/204). A standardized polygenic score of these 204 variants was associated with incident CVD with HR 1.14 (P = 1.0 × 10-16). CONCLUSIONS: The data point to novel and known genomic regions associated with incident CVD among individuals with T2D.

Lac-Phe mediates the effects of metformin on food intake and body weight.

Metformin is a widely prescribed anti-diabetic medicine that also reduces body weight. There is ongoing debate about the mechanisms that mediate metformin's effects on energy balance. Here, we show that metformin is a powerful pharmacological inducer of the anorexigenic metabolite N-lactoyl-phenylalanine (Lac-Phe) in cells, in mice and two independent human cohorts. Metformin drives Lac-Phe biosynthesis through the inhibition of complex I, increased glycolytic flux and intracellular lactate mass action. Intestinal epithelial CNDP2+ cells, not macrophages, are the principal in vivo source of basal and metformin-inducible Lac-Phe. Genetic ablation of Lac-Phe biosynthesis in male mice renders animals resistant to the effects of metformin on food intake and body weight. Lastly, mediation analyses support a role for Lac-Phe as a downstream effector of metformin's effects on body mass index in participants of a large population-based observational cohort, the Multi-Ethnic Study of Atherosclerosis. Together, these data establish Lac-Phe as a critical mediator of the body weight-lowering effects of metformin.