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.

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.

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.

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.

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.

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.

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.

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.

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.

Transcriptome-wide association study of the plasma proteome reveals cis and trans regulatory mechanisms underlying complex traits.

Regulation of transcription and translation are mechanisms through which genetic variants affect complex traits. Expression quantitative trait locus (eQTL) studies have been more successful at identifying cis-eQTL (within 1 Mb of the transcription start site) than trans-eQTL. Here, we tested the cis component of gene expression for association with observed plasma protein levels to identify cis- and trans-acting genes that regulate protein levels. We used transcriptome prediction models from 49 Genotype-Tissue Expression (GTEx) Project tissues to predict the cis component of gene expression and tested the predicted expression of every gene in every tissue for association with the observed abundance of 3,622 plasma proteins measured in 3,301 individuals from the INTERVAL study. We tested significant results for replication in 971 individuals from the Trans-omics for Precision Medicine (TOPMed) Multi-Ethnic Study of Atherosclerosis (MESA). We found 1,168 and 1,210 cis- and trans-acting associations that replicated in TOPMed (FDR < 0.05) with a median expected true positive rate (π1) across tissues of 0.806 and 0.390, respectively. The target proteins of trans-acting genes were enriched for transcription factor binding sites and autoimmune diseases in the GWAS catalog. Furthermore, we found a higher correlation between predicted expression and protein levels of the same underlying gene (R = 0.17) than observed expression (R = 0.10, p = 7.50 × 10-11). This indicates the cis-acting genetically regulated (heritable) component of gene expression is more consistent across tissues than total observed expression (genetics + environment) and is useful in uncovering the function of SNPs associated with complex traits.

Genome-wide association analysis of composite sleep health scores in 413,904 individuals.

Recent genome-wide association studies (GWASs) of several individual sleep traits have identified hundreds of genetic loci, suggesting diverse mechanisms. Moreover, sleep traits are moderately correlated, and together may provide a more complete picture of sleep health, while also illuminating distinct domains. Here we construct novel sleep health scores (SHSs) incorporating five core self-report measures: sleep duration, insomnia symptoms, chronotype, snoring, and daytime sleepiness, using additive (SHS-ADD) and five principal components-based (SHS-PCs) approaches. GWASs of these six SHSs identify 28 significant novel loci adjusting for multiple testing on six traits (p<8.3e-9), along with 341 previously reported loci (p<5e-08). The heritability of the first three SHS-PCs equals or exceeds that of SHS-ADD (SNP-h2=0.094), while revealing sleep-domain-specific genetic discoveries. Significant loci enrich in multiple brain tissues and in metabolic and neuronal pathways. Post GWAS analyses uncover novel genetic mechanisms underlying sleep health and reveal connections to behavioral, psychological, and cardiometabolic traits.

Refinement of a published gene-physical activity interaction impacting HDL-cholesterol: role of sex and lipoprotein subfractions.

Large-scale gene-environment interaction (GxE) discovery efforts often involve compromises in the definition of outcomes and choice of covariates for the sake of data harmonization and statistical power. Consequently, refinement of exposures, covariates, outcomes, and population subsets may be helpful to establish often-elusive replication and evaluate potential clinical utility. Here, we used additional datasets, an expanded set of statistical models, and interrogation of lipoprotein metabolism via nuclear magnetic resonance (NMR)-based lipoprotein subfractions to refine a previously discovered GxE modifying the relationship between physical activity (PA) and HDL-cholesterol (HDL-C). This GxE was originally identified by Kilpeläinen et al., with the strongest cohort-specific signal coming from the Women's Genome Health Study (WGHS). We thus explored this GxE further in the WGHS (N = 23,294), with follow-up in the UK Biobank (UKB; N = 281,380), and the Multi-Ethnic Study of Atherosclerosis (MESA; N = 4,587). Self-reported PA (MET-hrs/wk), genotypes at rs295849 (nearest gene: LHX1), and NMR metabolomics data were available in all three cohorts. As originally reported, minor allele carriers of rs295849 in WGHS had a stronger positive association between PA and HDL-C (pint = 0.002). When testing a range of NMR metabolites (primarily lipoprotein and lipid subfractions) to refine the HDL-C outcome, we found a stronger interaction effect on medium-sized HDL particle concentrations (M-HDL-P; pint = 1.0×10-4) than HDL-C. Meta-regression revealed a systematically larger interaction effect in cohorts from the original meta-analysis with a greater fraction of women (p = 0.018). In the UKB, GxE effects were stronger both in women and using M-HDL-P as the outcome. In MESA, the primary interaction for HDL-C showed nominal significance (pint = 0.013), but without clear differences by sex and with a greater magnitude using large, rather than medium, HDL-P as an outcome. Towards reconciling these observations, further exploration leveraging NMR platform-specific HDL subfraction diameter annotations revealed modest agreement across all cohorts in the interaction affecting medium-to-large particles. Taken together, our work provides additional insights into a specific known gene-PA interaction while illustrating the importance of phenotype and model refinement towards understanding and replicating GxEs.

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.

A multi-cohort genome-wide association study in African ancestry individuals reveals risk loci for primary open-angle glaucoma.

Primary open-angle glaucoma (POAG), the leading cause of irreversible blindness worldwide, disproportionately affects individuals of African ancestry. We conducted a genome-wide association study (GWAS) for POAG in 11,275 individuals of African ancestry (6,003 cases; 5,272 controls). We detected 46 risk loci associated with POAG at genome-wide significance. Replication and post-GWAS analyses, including functionally informed fine-mapping, multiple trait co-localization, and in silico validation, implicated two previously undescribed variants (rs1666698 mapping to DBF4P2; rs34957764 mapping to ROCK1P1) and one previously associated variant (rs11824032 mapping to ARHGEF12) as likely causal. For individuals of African ancestry, a polygenic risk score (PRS) for POAG from our mega-analysis (African ancestry individuals) outperformed a PRS from summary statistics of a much larger GWAS derived from European ancestry individuals. This study quantifies the genetic architecture similarities and differences between African and non-African ancestry populations for this blinding disease.

Genome-wide study investigating effector genes and polygenic prediction for kidney function in persons with ancestry from Africa and the Americas.

Chronic kidney disease is a leading cause of death and disability globally and impacts individuals of African ancestry (AFR) or with ancestry in the Americas (AMS) who are under-represented in genome-wide association studies (GWASs) of kidney function. To address this bias, we conducted a large meta-analysis of GWASs of estimated glomerular filtration rate (eGFR) in 145,732 AFR and AMS individuals. We identified 41 loci at genome-wide significance (p < 5 × 10-8), of which two have not been previously reported in any ancestry group. We integrated fine-mapped loci with epigenomic and transcriptomic resources to highlight potential effector genes relevant to kidney physiology and disease, and reveal key regulatory elements and pathways involved in renal function and development. We demonstrate the varying but increased predictive power offered by a multi-ancestry polygenic score for eGFR and highlight the importance of population diversity in GWASs and multi-omics resources to enhance opportunities for clinical translation for all.

Interaction molecular QTL mapping discovers cellular and environmental modifiers of genetic regulatory effects.

Bulk-tissue molecular quantitative trait loci (QTLs) have been the starting point for interpreting disease-associated variants, and context-specific QTLs show particular relevance for disease. Here, we present the results of mapping interaction QTLs (iQTLs) for cell type, age, and other phenotypic variables in multi-omic, longitudinal data from the blood of individuals of diverse ancestries. By modeling the interaction between genotype and estimated cell-type proportions, we demonstrate that cell-type iQTLs could be considered as proxies for cell-type-specific QTL effects, particularly for the most abundant cell type in the tissue. The interpretation of age iQTLs, however, warrants caution because the moderation effect of age on the genotype and molecular phenotype association could be mediated by changes in cell-type composition. Finally, we show that cell-type iQTLs contribute to cell-type-specific enrichment of diseases that, in combination with additional functional data, could guide future functional studies. Overall, this study highlights the use of iQTLs to gain insights into the context specificity of regulatory effects.