A Large-Scale Genome-Wide Gene-Sleep Interaction Study in 732,564 Participants Identifies Lipid Loci Explaining Sleep-Associated Lipid Disturbances.

We performed large-scale genome-wide gene-sleep interaction analyses of lipid levels to identify novel genetic variants underpinning the biomolecular pathways of sleep-associated lipid disturbances and to suggest possible druggable targets. We collected data from 55 cohorts with a combined sample size of 732,564 participants (87% European ancestry) with data on lipid traits (high-density lipoprotein [HDL-c] and low-density lipoprotein [LDL-c] cholesterol and triglycerides [TG]). Short (STST) and long (LTST) total sleep time were defined by the extreme 20% of the age- and sex-standardized values within each cohort. Based on cohort-level summary statistics data, we performed meta-analyses for the one-degree of freedom tests of interaction and two-degree of freedom joint tests of the main and interaction effect. In the cross-population meta-analyses, the one-degree of freedom variant-sleep interaction test identified 10 loci (P int <5.0e-9) not previously observed for lipids. Of interest, the ASPH locus (TG, LTST) is a target for aspartic and succinic acid metabolism previously shown to improve sleep and cardiovascular risk. The two-degree of freedom analyses identified an additional 7 loci that showed evidence for variant-sleep interaction (P joint <5.0e-9 in combination with P int <6.6e-6). Of these, the SLC8A1 locus (TG, STST) has been considered a potential treatment target for reduction of ischemic damage after acute myocardial infarction. Collectively, the 17 (9 with STST; 8 with LTST) loci identified in this large-scale initiative provides evidence into the biomolecular mechanisms underpinning sleep-duration-associated changes in lipid levels. The identified druggable targets may contribute to the development of novel therapies for dyslipidemia in people with sleep disturbances.

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.

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.

Participant-derived cell line transcriptomic analyses and mouse studies reveal a role for ZNF335 in plasma cholesterol statin response.

BACKGROUND: Statins lower circulating low-density lipoprotein cholesterol (LDLC) levels and reduce cardiovascular disease risk. Though highly efficacious in general, there is considerable inter-individual variation in statin efficacy that remains largely unexplained. METHODS: To identify novel genes that may modulate statin-induced LDLC lowering, we used RNA-sequencing data from 426 control- and 2 µM simvastatin-treated lymphoblastoid cell lines (LCLs) derived from European and African American ancestry participants of the Cholesterol and Pharmacogenetics (CAP) 40 mg/day 6-week simvastatin clinical trial (ClinicalTrials.gov Identifier: NCT00451828). We correlated statin-induced changes in LCL gene expression with plasma LDLC statin response in the corresponding CAP participants. For the most correlated gene identified (ZNF335), we followed up in vivo by comparing plasma cholesterol levels, lipoprotein profiles, and lipid statin response between wild-type mice and carriers of a hypomorphic (partial loss of function) missense mutation in Zfp335 (the mouse homolog of ZNF335). RESULTS: The statin-induced expression changes of 147 human LCL genes were significantly correlated to the plasma LDLC statin responses of the corresponding CAP participants in vivo (FDR = 5%). The two genes with the strongest correlations were zinc finger protein 335 (ZNF335 aka NIF-1, rho = 0.237, FDR-adj p = 0.0085) and CCR4-NOT transcription complex subunit 3 (CNOT3, rho = 0.233, FDR-adj p = 0.0085). Chow-fed mice carrying a hypomorphic missense (R1092W; aka bloto) mutation in Zfp335 had significantly lower non-HDL cholesterol levels than wild-type C57BL/6J mice in a sex combined model (p = 0.04). Furthermore, male (but not female) mice carrying the Zfp335R1092W allele had significantly lower total and HDL cholesterol levels than wild-type mice. In a separate experiment, wild-type mice fed a control diet for 4 weeks and a matched simvastatin diet for an additional 4 weeks had significant statin-induced reductions in non-HDLC (-43 ± 18% and -23 ± 19% for males and females, respectively). Wild-type male (but not female) mice experienced significant reductions in plasma LDL particle concentrations, while male mice carrying Zfp335R1092W allele(s) exhibited a significantly blunted LDL statin response. CONCLUSIONS: Our in vitro and in vivo studies identified ZNF335 as a novel modulator of plasma cholesterol levels and statin response, suggesting that variation in ZNF335 activity could contribute to inter-individual differences in statin clinical efficacy.

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.

Machine learning models for predicting blood pressure phenotypes by combining multiple polygenic risk scores.

We construct non-linear machine learning (ML) prediction models for systolic and diastolic blood pressure (SBP, DBP) using demographic and clinical variables and polygenic risk scores (PRSs). We developed a two-model ensemble, consisting of a baseline model, where prediction is based on demographic and clinical variables only, and a genetic model, where we also include PRSs. We evaluate the use of a linear versus a non-linear model at both the baseline and the genetic model levels and assess the improvement in performance when incorporating multiple PRSs. We report the ensemble model's performance as percentage variance explained (PVE) on a held-out test dataset. A non-linear baseline model improved the PVEs from 28.1 to 30.1% (SBP) and 14.3% to 17.4% (DBP) compared with a linear baseline model. Including seven PRSs in the genetic model computed based on the largest available GWAS of SBP/DBP improved the genetic model PVE from 4.8 to 5.1% (SBP) and 4.7 to 5% (DBP) compared to using a single PRS. Adding additional 14 PRSs computed based on two independent GWASs further increased the genetic model PVE to 6.3% (SBP) and 5.7% (DBP). PVE differed across self-reported race/ethnicity groups, with primarily all non-White groups benefitting from the inclusion of additional PRSs. In summary, non-linear ML models improves BP prediction in models incorporating diverse populations.

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.

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.

Multi-Ancestry Polygenic Risk Score for Coronary Heart Disease Based on an Ancestrally Diverse Genome-Wide Association Study and Population-Specific Optimization.

BACKGROUND: Predictive performance of polygenic risk scores (PRS) varies across populations. To facilitate equitable clinical use, we developed PRS for coronary heart disease (CHD; PRSCHD) for 5 genetic ancestry groups. METHODS: We derived ancestry-specific and multi-ancestry PRSCHD based on pruning and thresholding (PRSPT) and ancestry-based continuous shrinkage priors (PRSCSx) applied to summary statistics from the largest multi-ancestry genome-wide association study meta-analysis for CHD to date, including 1.1 million participants from 5 major genetic ancestry groups. Following training and optimization in the Million Veteran Program, we evaluated the best-performing PRSCHD in 176,988 individuals across 9 diverse cohorts. RESULTS: Multi-ancestry PRSPT and PRSCSx outperformed ancestry-specific PRSPT and PRSCSx across a range of tuning values. Two best-performing multi-ancestry PRSCHD (ie, PRSPTmult and PRSCSxmult) and 1 ancestry-specific (PRSCSxEUR) were taken forward for validation. PRSPTmult demonstrated the strongest association with CHD in individuals of South Asian ancestry and European ancestry (odds ratio per 1 SD [95% CI, 2.75 [2.41-3.14], 1.65 [1.59-1.72]), followed by East Asian ancestry (1.56 [1.50-1.61]), Hispanic/Latino ancestry (1.38 [1.24-1.54]), and African ancestry (1.16 [1.11-1.21]). PRSCSxmult showed the strongest associations in South Asian ancestry (2.67 [2.38-3.00]) and European ancestry (1.65 [1.59-1.71]), lower in East Asian ancestry (1.59 [1.54-1.64]), Hispanic/Latino ancestry (1.51 [1.35-1.69]), and the lowest in African ancestry (1.20 [1.15-1.26]). CONCLUSIONS: The use of summary statistics from a large multi-ancestry genome-wide meta-analysis improved the performance of PRSCHD in most ancestry groups compared with single-ancestry methods. Despite the use of one of the largest and most diverse sets of training and validation cohorts to date, improvement of predictive performance was limited in African ancestry. This highlights the need for larger genome-wide association study datasets of underrepresented populations to enhance the performance of PRSCHD.

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 association identifies novel ROP risk loci in a multiethnic cohort.

We conducted a genome-wide association study (GWAS) in a multiethnic cohort of 920 at-risk infants for retinopathy of prematurity (ROP), a major cause of childhood blindness, identifying 1 locus at genome-wide significance level (p < 5×10-8) and 9 with significance of p < 5×10-6 for ROP ≥ stage 3. The most significant locus, rs2058019, reached genome-wide significance within the full multiethnic cohort (p = 4.96×10-9); Hispanic and European Ancestry infants driving the association. The lead single nucleotide polymorphism (SNP) falls in an intronic region within the Glioma-associated oncogene family zinc finger 3 (GLI3) gene. Relevance for GLI3 and other top-associated genes to human ocular disease was substantiated through in-silico extension analyses, genetic risk score analysis and expression profiling in human donor eye tissues. Thus, we identify a novel locus at GLI3 with relevance to retinal biology, supporting genetic susceptibilities for ROP risk with possible variability by race and ethnicity.

Increased Plasma Branched Short-Chain Fatty Acids and Improved Glucose Homeostasis: The Microbiome and Insulin Longitudinal Evaluation Study (MILES).

Short-chain fatty acids (SCFAs) have been extensively studied for potential beneficial roles in glucose homeostasis and risk of diabetes; however, most of this research has focused on butyrate, acetate, and propionate. The effect on metabolism of branched SCFAs (BSCFAs; isobutyrate, isovalerate, and methylbutyrate) is largely unknown. In a cohort of 219 non-Hispanic White participants and 126 African American participants, we examined the association of BSCFA with dysglycemia (prediabetes and diabetes) and oral glucose tolerance test-based measures of glucose and insulin homeostasis, as well as with demographic, anthropometric, lifestyle, and lipid traits, and other SCFAs. We observed a bimodal distribution of BSCFAs, with 25 individuals having high levels (H-BSCFA group) and 320 individuals having lower levels (L-BSCFA group). The prevalence of dysglycemia was lower in the H-BSCFA group compared with the L-BSCFA group (16% vs. 49%; P = 0.0014). This association remained significant after adjustment for age, sex, race, BMI, and levels of other SCFAs. Consistent with the lower rate of dysglycemia, fasting and postprandial glucose levels were lower and the disposition index was higher in the H-BSCFA group. Additional findings in H-BSCFA versus L-BSCFA included lower fasting and postprandial C-peptide levels and lower insulin clearance without differences in insulin levels, insulin sensitivity, insulin secretion, or other variables examined, including diet and physical activity. As one of the first human studies associating higher BSCFA levels with lower odds of dysglycemia and improved glucose homeostasis, this study sets the stage for further investigation of BSCFA as a novel target for prevention or treatment of diabetes.

A statistical framework for powerful multi-trait rare variant analysis in large-scale whole-genome sequencing studies.

Large-scale whole-genome sequencing (WGS) studies have improved our understanding of the contributions of coding and noncoding rare variants to complex human traits. Leveraging association effect sizes across multiple traits in WGS rare variant association analysis can improve statistical power over single-trait analysis, and also detect pleiotropic genes and regions. Existing multi-trait methods have limited ability to perform rare variant analysis of large-scale WGS data. We propose MultiSTAAR, a statistical framework and computationally-scalable analytical pipeline for functionally-informed multi-trait rare variant analysis in large-scale WGS studies. MultiSTAAR accounts for relatedness, population structure and correlation among phenotypes by jointly analyzing multiple traits, and further empowers rare variant association analysis by incorporating multiple functional annotations. We applied MultiSTAAR to jointly analyze three lipid traits (low-density lipoprotein cholesterol, high-density lipoprotein cholesterol and triglycerides) in 61,861 multi-ethnic samples from the Trans-Omics for Precision Medicine (TOPMed) Program. We discovered new associations with lipid traits missed by single-trait analysis, including rare variants within an enhancer of NIPSNAP3A and an intergenic region on chromosome 1.

Lac-Phe mediates the anti-obesity effect of metformin.

Metformin is a widely prescribed anti-diabetic medicine that also reduces body weight. The mechanisms that mediate metformin's effects on energy balance remain incompletely defined. Here we show that metformin is a powerful pharmacological inducer of the anorexigenic metabolite Lac-Phe in mice as well as in two independent human cohorts. In cell culture, metformin drives Lac-Phe biosynthesis via inhibition of complex I, increased glycolytic flux, and intracellular lactate mass action. Other biguanides and structurally distinct inhibitors of oxidative phosphorylation also increase Lac-Phe levels in vitro. Genetic ablation of CNDP2, the principal biosynthetic enzyme for Lac-Phe, in mice renders animals resistant to metformin's anorexigenic and anti-obesity effects. Mediation analyses also support a role for Lac-Phe in metformin's effect on body mass index in humans. These data establish the CNDP2/Lac-Phe pathway as a critical mediator of the effects of metformin on energy balance.

Rare variants in long non-coding RNAs are associated with blood lipid levels in the TOPMed whole-genome sequencing study.

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions in lipid metabolism. Large-scale whole-genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess more associations between rare variants in lncRNA genes and complex traits across the genome. In this study, we used high-coverage WGS from 66,329 participants of diverse ancestries with measurement of blood lipids and lipoproteins (LDL-C, HDL-C, TC, and TG) in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) program to investigate the role of lncRNAs in lipid variability. We aggregated rare variants for 165,375 lncRNA genes based on their genomic locations and conducted rare-variant aggregate association tests using the STAAR (variant-set test for association using annotation information) framework. We performed STAAR conditional analysis adjusting for common variants in known lipid GWAS loci and rare-coding variants in nearby protein-coding genes. Our analyses revealed 83 rare lncRNA variant sets significantly associated with blood lipid levels, all of which were located in known lipid GWAS loci (in a ±500-kb window of a Global Lipids Genetics Consortium index variant). Notably, 61 out of 83 signals (73%) were conditionally independent of common regulatory variation and rare protein-coding variation at the same loci. We replicated 34 out of 61 (56%) conditionally independent associations using the independent UK Biobank WGS data. Our results expand the genetic architecture of blood lipids to rare variants in lncRNAs.

Genome-wide association meta-analysis identifies 17 loci associated with nonalcoholic fatty liver disease.

Nonalcoholic fatty liver disease (NAFLD) is common and partially heritable and has no effective treatments. We carried out a genome-wide association study (GWAS) meta-analysis of imaging (n = 66,814) and diagnostic code (3,584 cases versus 621,081 controls) measured NAFLD across diverse ancestries. We identified NAFLD-associated variants at torsin family 1 member B (TOR1B), fat mass and obesity associated (FTO), cordon-bleu WH2 repeat protein like 1 (COBLL1)/growth factor receptor-bound protein 14 (GRB14), insulin receptor (INSR), sterol regulatory element-binding transcription factor 1 (SREBF1) and patatin-like phospholipase domain-containing protein 2 (PNPLA2), as well as validated NAFLD-associated variants at patatin-like phospholipase domain-containing protein 3 (PNPLA3), transmembrane 6 superfamily 2 (TM6SF2), apolipoprotein E (APOE), glucokinase regulator (GCKR), tribbles homolog 1 (TRIB1), glycerol-3-phosphate acyltransferase (GPAM), mitochondrial amidoxime-reducing component 1 (MARC1), microsomal triglyceride transfer protein large subunit (MTTP), alcohol dehydrogenase 1B (ADH1B), transmembrane channel like 4 (TMC4)/membrane-bound O-acyltransferase domain containing 7 (MBOAT7) and receptor-type tyrosine-protein phosphatase δ (PTPRD). Implicated genes highlight mitochondrial, cholesterol and de novo lipogenesis as causally contributing to NAFLD predisposition. Phenome-wide association study (PheWAS) analyses suggest at least seven subtypes of NAFLD. Individuals in the top 10% and 1% of genetic risk have a 2.5-fold to 6-fold increased risk of NAFLD, cirrhosis and hepatocellular carcinoma. These genetic variants identify subtypes of NAFLD, improve estimates of disease risk and can guide the development of targeted therapeutics.

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

BACKGROUND: Type 2 diabetes mellitus (T2D) confers a two- to three-fold increased risk of cardiovascular disease (CVD). However, the mechanisms underlying increased CVD risk among people with T2D are only partially understood. We hypothesized that a genetic association study among people with T2D at risk for developing incident cardiovascular complications could provide insights into molecular genetic aspects underlying CVD. METHODS: From 16 studies of the Cohorts for Heart & Aging Research in Genomic Epidemiology (CHARGE) Consortium, we conducted a multi-ancestry time-to-event genome-wide association study (GWAS) for incident CVD among people with T2D using Cox proportional hazards models. Incident CVD was defined based on a composite of coronary artery disease (CAD), stroke, and cardiovascular death that occurred at least one year after the diagnosis of T2D. Cohort-level estimated effect sizes were combined using inverse variance weighted fixed effects meta-analysis. We also tested 204 known CAD variants for association with incident CVD among patients with T2D. RESULTS: A total of 49,230 participants with T2D were included in the analyses (31,118 European ancestries and 18,112 non-European ancestries) which consisted of 8,956 incident CVD cases over a range of mean follow-up duration between 3.2 and 33.7 years (event rate 18.2%). We identified three novel, distinct genetic loci for incident CVD among individuals with T2D that reached the threshold for genome-wide significance (P<5.0×10-8): rs147138607 (intergenic variant between CACNA1E and ZNF648) with a hazard ratio (HR) 1.23, 95% confidence interval (CI) 1.15 - 1.32, P=3.6×10-9, rs11444867 (intergenic variant near HS3ST1) with HR 1.89, 95% CI 1.52 - 2.35, P=9.9×10-9, and rs335407 (intergenic variant between TFB1M and NOX3) HR 1.25, 95% CI 1.16 - 1.35, P=1.5×10-8. Among 204 known CAD loci, 32 were associated with incident CVD in people with T2D with P<0.05, and 5 were significant after Bonferroni correction (P<0.00024, 0.05/204). A polygenic score of these 204 variants was significantly associated with incident CVD with HR 1.14 (95% CI 1.12 - 1.16) per 1 standard deviation increase (P=1.0×10-16). CONCLUSIONS: The data point to novel and known genomic regions associated with incident CVD among individuals with T2D.

Protein-metabolite association studies identify novel proteomic determinants of metabolite levels in human plasma.

Although many novel gene-metabolite and gene-protein associations have been identified using high-throughput biochemical profiling, systematic studies that leverage human genetics to illuminate causal relationships between circulating proteins and metabolites are lacking. Here, we performed protein-metabolite association studies in 3,626 plasma samples from three human cohorts. We detected 171,800 significant protein-metabolite pairwise correlations between 1,265 proteins and 365 metabolites, including established relationships in metabolic and signaling pathways such as the protein thyroxine-binding globulin and the metabolite thyroxine, as well as thousands of new findings. In Mendelian randomization (MR) analyses, we identified putative causal protein-to-metabolite associations. We experimentally validated top MR associations in proof-of-concept plasma metabolomics studies in three murine knockout strains of key protein regulators. These analyses identified previously unrecognized associations between bioactive proteins and metabolites in human plasma. We provide publicly available data to be leveraged for studies in human metabolism and disease.