Meta-Analysis of Genome-Wide Association Studies Reveals Genetic Mechanisms of Supraventricular Arrhythmias.

BACKGROUND: Substantial data support a heritable basis for supraventricular tachycardias, but the genetic determinants and molecular mechanisms of these arrhythmias are poorly understood. We sought to identify genetic loci associated with atrioventricular nodal reentrant tachycardia (AVNRT) and atrioventricular accessory pathways or atrioventricular reciprocating tachycardia (AVAPs/AVRT). METHODS: We performed multiancestry meta-analyses of genome-wide association studies to identify genetic loci for AVNRT (4 studies) and AVAP/AVRT (7 studies). We assessed evidence supporting the potential causal effects of candidate genes by analyzing relations between associated variants and cardiac gene expression, performing transcriptome-wide analyses, and examining prior genome-wide association studies. RESULTS: Analyses comprised 2384 AVNRT cases and 106 489 referents, and 2811 AVAP/AVRT cases and 1,483 093 referents. We identified 2 significant loci for AVNRT, which implicates NKX2-5 and TTN as disease susceptibility genes. A transcriptome-wide association analysis supported an association between reduced predicted cardiac expression of NKX2-5 and AVNRT. We identified 3 significant loci for AVAP/AVRT, which implicates SCN5A, SCN10A, and TTN/CCDC141. Variant associations at several loci have been previously reported for cardiac phenotypes, including atrial fibrillation, stroke, Brugada syndrome, and electrocardiographic intervals. CONCLUSIONS: Our findings highlight gene regions associated with ion channel function (AVAP/AVRT), as well as cardiac development and the sarcomere (AVAP/AVRT and AVNRT) as important potential effectors of supraventricular tachycardia susceptibility.

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.

Titin-Truncating variants Predispose to Dilated Cardiomyopathy in Diverse Populations.

Importance: The effect of high percentage spliced in (hiPSI) TTN truncating variants (TTNtvs) on risk of dilated cardiomyopathy (DCM) has historically been studied among population subgroups defined by genetic similarity to European reference populations. This has raised questions about the effect of TTNtvs in diverse populations, especially among individuals genetically similar to African reference populations. Objective: To determine the effect of TTNtvs on risk of DCM in diverse population as measured by genetic distance (GD) in principal component (PC) space. Design: Cohort study. Setting: Penn Medicine Biobank (PMBB) is a large, diverse biobank. Participants: Participants were recruited from across the Penn Medicine healthcare system and volunteered to have their electronic health records linked to biospecimen data including DNA which has undergone whole exome sequencing. Main Outcomes and Measures: Risk of DCM among individuals carrying a hiPSI TTNtv. Results: Carrying a hiPSI TTNtv was associated with DCM among PMBB participants across a range of GD deciles from the 1000G European centroid; the effect estimates ranged from odds ratio (OR) = 3.29 (95% confidence interval [CI] 1.26 to 8.56) to OR = 9.39 (95% CI 3.82 to 23.13). When individuals were assigned to population subgroups based on genetic similarity to the 1000G reference populations, hiPSI TTNtvs conferred significant risk of DCM among those genetically similar to the 1000G European reference population (OR = 7.55, 95% CI 4.99 to 11.42, P<0.001) and individuals genetically similar to the 1000G African reference population (OR 3.50, 95% CI 1.48 to 8.24, P=0.004). Conclusions and Relevance: TTNtvs are associated with increased risk of DCM among a diverse cohort. There is no significant difference in effect of TTNtvs on DCM risk across deciles of GD from the 1000G European centroid, suggesting genetic background should not be considered when screening individuals for titin-related DCM.

Polygenic risk scores for cardiometabolic traits demonstrate importance of ancestry for predictive precision medicine.

Polygenic risk scores (PRS) have predominantly been derived from genome-wide association studies (GWAS) conducted in European ancestry (EUR) individuals. In this study, we present an in-depth evaluation of PRS based on multi-ancestry GWAS for five cardiometabolic phenotypes in the Penn Medicine BioBank (PMBB) followed by a phenome-wide association study (PheWAS). We examine the PRS performance across all individuals and separately in African ancestry (AFR) and EUR ancestry groups. For AFR individuals, PRS derived using the multi-ancestry LD panel showed a higher effect size for four out of five PRSs (DBP, SBP, T2D, and BMI) than those derived from the AFR LD panel. In contrast, for EUR individuals, the multi-ancestry LD panel PRS demonstrated a higher effect size for two out of five PRSs (SBP and T2D) compared to the EUR LD panel. These findings underscore the potential benefits of utilizing a multi-ancestry LD panel for PRS derivation in diverse genetic backgrounds and demonstrate overall robustness in all individuals. Our results also revealed significant associations between PRS and various phenotypic categories. For instance, CAD PRS was linked with 18 phenotypes in AFR and 82 in EUR, while T2D PRS correlated with 84 phenotypes in AFR and 78 in EUR. Notably, associations like hyperlipidemia, renal failure, atrial fibrillation, coronary atherosclerosis, obesity, and hypertension were observed across different PRSs in both AFR and EUR groups, with varying effect sizes and significance levels. However, in AFR individuals, the strength and number of PRS associations with other phenotypes were generally reduced compared to EUR individuals. Our study underscores the need for future research to prioritize 1) conducting GWAS in diverse ancestry groups and 2) creating a cosmopolitan PRS methodology that is universally applicable across all genetic backgrounds. Such advances will foster a more equitable and personalized approach to precision medicine.

Genome-wide association meta-analysis identifies risk loci for abdominal aortic aneurysm and highlights PCSK9 as a therapeutic target.

Abdominal aortic aneurysm (AAA) is a common disease with substantial heritability. In this study, we performed a genome-wide association meta-analysis from 14 discovery cohorts and uncovered 141 independent associations, including 97 previously unreported loci. A polygenic risk score derived from meta-analysis explained AAA risk beyond clinical risk factors. Genes at AAA risk loci indicate involvement of lipid metabolism, vascular development and remodeling, extracellular matrix dysregulation and inflammation as key mechanisms in AAA pathogenesis. These genes also indicate overlap between the development of AAA and other monogenic aortopathies, particularly via transforming growth factor ß signaling. Motivated by the strong evidence for the role of lipid metabolism in AAA, we used Mendelian randomization to establish the central role of nonhigh-density lipoprotein cholesterol in AAA and identified the opportunity for repurposing of proprotein convertase, subtilisin/kexin-type 9 (PCSK9) inhibitors. This was supported by a study demonstrating that PCSK9 loss of function prevented the development of AAA in a preclinical mouse model.

Relationships between body fat distribution and metabolic syndrome traits and outcomes: A mendelian randomization study.

BACKGROUND: Obesity is a complex, multifactorial disease associated with substantial morbidity and mortality worldwide. Although it is frequently assessed using BMI, many epidemiological studies have shown links between body fat distribution and obesity-related outcomes. This study examined the relationships between body fat distribution and metabolic syndrome traits using Mendelian Randomization (MR). METHODS/FINDINGS: Genetic variants associated with visceral adipose tissue (VAT), abdominal subcutaneous adipose tissue (ASAT), and gluteofemoral adipose tissue (GFAT), as well as their relative ratios, were identified from a genome wide association study (GWAS) performed with the United Kingdom BioBank. GWAS summary statistics for traits and outcomes related to metabolic syndrome were obtained from the IEU Open GWAS Project. Two-sample MR and BMI-controlled multivariable MR (MVMR) were performed to examine relationships between each body fat measure and ratio with the outcomes. Increases in absolute GFAT were associated with a protective cardiometabolic profile, including lower low density lipoprotein cholesterol (ß: -0.19, [95% CI: -0.28, -0.10], p < 0.001), higher high density lipoprotein cholesterol (ß: 0.23, [95% CI: 0.03, 0.43], p = 0.025), lower triglycerides (ß: -0.28, [95% CI: -0.45, -0.10], p = 0.0021), and decreased systolic (ß: -1.65, [95% CI: -2.69, -0.61], p = 0.0019) and diastolic blood pressures (ß: -0.95, [95% CI: -1.65, -0.25], p = 0.0075). These relationships were largely maintained in BMI-controlled MVMR analyses. Decreases in relative GFAT were linked with a worse cardiometabolic profile, with higher levels of detrimental lipids and increases in systolic and diastolic blood pressures. CONCLUSION: A MR analysis of ASAT, GFAT, and VAT depots and their relative ratios with metabolic syndrome related traits and outcomes revealed that increased absolute and relative GFAT were associated with a favorable cardiometabolic profile independently of BMI. These associations highlight the importance of body fat distribution in obesity and more precise means to categorize obesity beyond BMI.

Using a polygenic score to account for genomic risk factors in a model to detect individuals with dilated ascending thoracic aortas.

Background: Ascending thoracic aortic dilation is a complex trait that involves modifiable and non-modifiable risk factors and can lead to thoracic aortic aneurysm and dissection. Clinical risk factors have been shown to predict ascending thoracic aortic diameter. Polygenic scores (PGS) are increasingly used to assess clinical risk for multifactorial diseases. The degree to which a PGS can improve aortic diameter prediction is not known. In this study we tested the extent to which the addition of a PGS to clinical prediction algorithms improves the prediction of aortic diameter. Methods: The patient cohort comprised 6,790 Penn Medicine Biobank (PMBB) participants with available echocardiography and clinical data linked to genome-wide genotype data. Linear regression models were used to integrate PGS weights derived from a large genome wide association study of thoracic aortic diameter in the UK biobank and were compared to the performance of the standard and a reweighted variation of the recently published AORTA Score. Results: Cohort participants were 56% male, had a median age of 61 years (IQR 52-70) with a mean ascending aortic diameter of 3.4 cm (SD 0.5). Compared to the AORTA Score which explained 28.4% (95% CI 28.1% to 29.2%) of the variance in aortic diameter, AORTA Score + PGS explained 28.8%, (95% CI 28.1% to 29.6%), the reweighted AORTA score explained 30.4% (95% CI 29.6% to 31.2%), and the reweighted AORTA Score + PGS explained 31.0% (95% CI 30.2% to 31.8%). The addition of a PGS to either the AORTA Score or the reweighted AORTA Score improved model sensitivity for the identifying individuals with a thoracic aortic diameter ≥ 4 cm. The respective areas under the receiver operator characteristic curve for the AORTA Score + PGS (0.771, 95% CI 0.756 to 0.787) and reweighted AORTA Score + PGS (0.785, 95% CI 0.770 to 0.800) were greater than the standard AORTA Score (0.767, 95% CI 0.751 to 0.783) and reweighted AORTA Score (0.780 95% CI 0.765 to 0.795). Conclusions: We demonstrated that inclusion of a PGS to the AORTA Score results in a small but clinically meaningful performance enhancement. Further investigation is necessary to determine if combining genetic and clinical risk prediction improves outcomes for thoracic aortic disease.

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.

Carriers of rare damaging CCR2 genetic variants are at lower risk of atherosclerotic disease.

Background: The CCL2/CCR2 axis governs monocyte trafficking and recruitment to atherosclerotic lesions. Human genetic analyses and population-based studies support an association between circulating CCL2 levels and atherosclerosis. Still, it remains unknown whether pharmacological targeting of CCR2, the main CCL2 receptor, would provide protection against human atherosclerotic disease. Methods: In whole-exome sequencing data from 454,775 UK Biobank participants (40-69 years), we identified predicted loss-of-function (LoF) or damaging missense (REVEL score >0.5) variants within the CCR2 gene. We prioritized variants associated with lower monocyte count (p<0.05) and tested associations with vascular risk factors and risk of atherosclerotic disease over a mean follow-up of 14 years. The results were replicated in a pooled cohort of three independent datasets (TOPMed, deCODE and Penn Medicine BioBank; total n=441,445) and the effect of the most frequent damaging variant was experimentally validated. Results: A total of 45 predicted LoF or damaging missense variants were identified in the CCR2 gene, 4 of which were also significantly associated with lower monocyte count, but not with other white blood cell counts. Heterozygous carriers of these variants were at a lower risk of a combined atherosclerosis outcome, showed a lower burden of atherosclerosis across four vascular beds, and were at a lower lifetime risk of coronary artery disease and myocardial infarction. There was no evidence of association with vascular risk factors including LDL-cholesterol, blood pressure, glycemic status, or C-reactive protein. Using a cAMP assay, we found that cells transfected with the most frequent CCR2 damaging variant (3:46358273:T:A, M249K, 547 carriers, frequency: 0.14%) show a decrease in signaling in response to CCL2. The associations of the M249K variant with myocardial infarction were consistent across cohorts (ORUKB: 0.62 95%CI: 0.39-0.96; ORexternal: 0.64 95%CI: 0.34-1.19; ORpooled: 0.64 95%CI: 0.450.90). In a phenome-wide association study, we found no evidence for higher risk of common infections or mortality among carriers of damaging CCR2 variants. Conclusions: Heterozygous carriers of damaging CCR2 variants have a lower burden of atherosclerosis and lower lifetime risk of myocardial infarction. In conjunction with previous evidence from experimental and epidemiological studies, our findings highlight the translational potential of CCR2-targeting as an atheroprotective approach.

Common- and rare-variant genetic architecture of heart failure across the allele frequency spectrum.

Heart failure (HF) is a complex trait, influenced by environmental and genetic factors, that affects over 30 million individuals worldwide. Historically, the genetics of HF have been studied in Mendelian forms of disease, where rare genetic variants have been linked to familial cardiomyopathies. More recently, genome-wide association studies (GWAS) have successfully identified common genetic variants associated with risk of HF. However, the relative importance of genetic variants across the allele-frequency spectrum remains incompletely characterized. Here, we report the results of common- and rare-variant association studies of all-cause heart failure, applying recently developed methods to quantify the heritability of HF attributable to different classes of genetic variation. We combine GWAS data across multiple populations including 207,346 individuals with HF and 2,151,210 without, identifying 176 risk loci at genome-wide significance (p < 5×10-8). Signals at newly identified common-variant loci include coding variants in Mendelian cardiomyopathy genes (MYBPC3, BAG3), as well as regulators of lipoprotein (LPL) and glucose metabolism (GIPR, GLP1R), and are enriched in cardiac, muscle, nerve, and vascular tissues, as well as myocyte and adipocyte cell types. Gene burden studies across three biobanks (PMBB, UKB, AOU) including 27,208 individuals with HF and 349,126 without uncover exome-wide significant (p < 3.15×10-6) associations for HF and rare predicted loss-of-function (pLoF) variants in TTN, MYBPC3, FLNC, and BAG3. Total burden heritability of rare coding variants (2.2%, 95% CI 0.99-3.5%) is highly concentrated in a small set of Mendelian cardiomyopathy genes, and is lower than heritability attributable to common variants (4.3%, 95% CI 3.9-4.7%) which is more diffusely spread throughout the genome. Finally, we demonstrate that common-variant background, in the form of a polygenic risk score (PRS), significantly modifies the risk of HF among carriers of pathogenic truncating variants in the Mendelian cardiomyopathy gene TTN. These findings suggest a significant polygenic component to HF exists that is not captured by current clinical genetic testing.

Genome-wide association study of thoracic aortic aneurysm and dissection in the Million Veteran Program.

The current understanding of the genetic determinants of thoracic aortic aneurysms and dissections (TAAD) has largely been informed through studies of rare, Mendelian forms of disease. Here, we conducted a genome-wide association study (GWAS) of TAAD, testing ~25 million DNA sequence variants in 8,626 participants with and 453,043 participants without TAAD in the Million Veteran Program, with replication in an independent sample of 4,459 individuals with and 512,463 without TAAD from six cohorts. We identified 21 TAAD risk loci, 17 of which have not been previously reported. We leverage multiple downstream analytic methods to identify causal TAAD risk genes and cell types and provide human genetic evidence that TAAD is a non-atherosclerotic aortic disorder distinct from other forms of vascular disease. Our results demonstrate that the genetic architecture of TAAD mirrors that of other complex traits and that it is not solely inherited through protein-altering variants of large effect size.

Evaluation of the Performance of the RECODe Equation with the Addition of Polygenic Risk Scores for Adverse Cardiovascular Outcomes in Individuals with Type II Diabetes.

Aims/Hypothesis: Individuals with T2D are at an increased risk of developing cardiovascular complications; early identification of individuals can lead to an alteration of the natural history of the disease. Current approaches to risk prediction tailored to individuals with T2D are exemplified by the RECODe algorithms which predict CVD outcomes among individuals with T2D. Recent efforts to improve CVD risk prediction among the general population have included the incorporation of polygenic risk scores (PRS). This paper aims to investigate the utility of the addition of a coronary artery disease (CAD), stroke and heart failure risk score to the current RECODe model for disease stratification. Methods: We derived PRS using summary statistics for ischemic stroke (IS) from the coronary artery disease (CAD) and heart failure (HF) and tested prediction accuracy in the Penn Medicine Biobank (PMBB). A Cox proportional hazards model was used for time-to-event analyses within our cohort, and we compared model discrimination for the RECODe model with and without a PRS using AUC. Results: The RECODe model alone demonstrated an AUC [95% CI] of 0.67 [0.62-0.72] for ASCVD; the addition of the three PRS to the model demonstrated an AUC [95% CI] of 0.66 [0.63-0.70]. A z-test to compare the AUCs of the two models did not demonstrate a detectable difference between the two models (p=0.97). Conclusions/Interpretation: In the present study, we demonstrate that although PRS associate with CVD outcomes independent of traditional risk factors among individuals with T2D, the addition of PRS to contemporary clinical risk models does not specifically improve the predictive performance as compared to the baseline model.

Session Introduction: SALUD: Scalable Applications of cLinical risk Utility and preDiction.

This PSB 2023 session discusses challenges in clinical implication and application of risk prediction models, which includes but is not limited to: implementation of risk models, responsible use of polygenic risk scores (PGS), and other risk prediction strategies. We focus on the development and use of new, scalable methods for harmonizing and refining risk prediction models by incorporating genetic and non-genetic risk factors, applying new phenotyping strategies, and integrating clinical factors and biomarkers. Lastly, we will discuss innovation in expanding the utility of these prediction models to underrepresented populations. This session focuses on the overarching theme of enabling early diagnosis, and treatment and preventive measures related to complex diseases and comorbidities.

Metabolite patterns associated with individual response to supervised exercise therapy in patients with intermittent claudication.

Objective: Supervised exercise therapy (SET) is the first line treatment for intermittent claudication owing to peripheral arterial disease. Despite multiple randomized controlled trials proving the efficacy of SET, there are large differences in individual patient's responses. We used plasma metabolomics to identify potential metabolic influences on the individual response to SET. Methods: Primary metabolites, complex lipids, and lipid mediators were measured on plasma samples taken at before and after Gardner graded treadmill walking tests that were administered before and after 12 weeks of SET. We used an ensemble modeling approach to identify metabolites or changes in metabolites at specific time points that associated with interindividual variability in the functional response to SET. Specific time points analyzed included baseline metabolite levels before SET, dynamic metabolomics changes before SET, the difference in pre- and post-SET baseline metabolomics, and the difference (pre- and post-SET) of the dynamic (pre- and post-treadmill). Results: High levels of baseline anandamide levels pre- and post-SET were associated with a worse response to SET. Increased arachidonic acid (AA) and decreased levels of the AA precursor dihomo-γ-linolenic acid across SET were associated with a worse response to SET. Participants who were able to tolerate large increases in AA during acute exercise had longer, or better, walking times both before and after SET. Conclusions: We identified two pathways of relevance to individual response to SET that warrant further study: anandamide synthesis may activate endocannabinoid receptors, resulting in worse treadmill test performance. SET may train patients to withstand higher levels of AA, and inflammatory signaling, resulting in longer walking times. Clinical Relevance: This manuscript describes the use of metabolomic techniques to measure the interindividual effects of SET in patients with peripheral artery disease (PAD). We identified high levels of AEA are linked to CB1 signaling and activation of inflammatory pathways. This alters energy expenditure in myoblasts by decreasing glucose uptake and may induce an acquired skeletal muscle myopathy. SET may also help participants tolerate increased levels of AA and inflammation produced during exercise, resulting in longer walking times. This data will enhance understanding of the pathophysiology of PAD and the mechanism by which SET improves walking intolerance.

The Penn Medicine BioBank: Towards a Genomics-Enabled Learning Healthcare System to Accelerate Precision Medicine in a Diverse Population.

The Penn Medicine BioBank (PMBB) is an electronic health record (EHR)-linked biobank at the University of Pennsylvania (Penn Medicine). A large variety of health-related information, ranging from diagnosis codes to laboratory measurements, imaging data and lifestyle information, is integrated with genomic and biomarker data in the PMBB to facilitate discoveries and translational science. To date, 174,712 participants have been enrolled into the PMBB, including approximately 30% of participants of non-European ancestry, making it one of the most diverse medical biobanks. There is a median of seven years of longitudinal data in the EHR available on participants, who also consent to permission to recontact. Herein, we describe the operations and infrastructure of the PMBB, summarize the phenotypic architecture of the enrolled participants, and use body mass index (BMI) as a proof-of-concept quantitative phenotype for PheWAS, LabWAS, and GWAS. The major representation of African-American participants in the PMBB addresses the essential need to expand the diversity in genetic and translational research. There is a critical need for a "medical biobank consortium" to facilitate replication, increase power for rare phenotypes and variants, and promote harmonized collaboration to optimize the potential for biological discovery and precision medicine.

Exome-wide association analysis of CT imaging-derived hepatic fat in a medical biobank.

Nonalcoholic fatty liver disease is common and highly heritable. Genetic studies of hepatic fat have not sufficiently addressed non-European and rare variants. In a medical biobank, we quantitate hepatic fat from clinical computed tomography (CT) scans via deep learning in 10,283 participants with whole-exome sequences available. We conduct exome-wide associations of single variants and rare predicted loss-of-function (pLOF) variants with CT-based hepatic fat and perform cross-modality replication in the UK Biobank (UKB) by linking whole-exome sequences to MRI-based hepatic fat. We confirm single variants previously associated with hepatic fat and identify several additional variants, including two (FGD5 H600Y and CITED2 S198_G199del) that replicated in UKB. A burden of rare pLOF variants in LMF2 is associated with increased hepatic fat and replicates in UKB. Quantitative phenotypes generated from clinical imaging studies and intersected with genomic data in medical biobanks have the potential to identify molecular pathways associated with human traits and disease.

Genome-wide association and multi-trait analyses characterize the common genetic architecture of heart failure.

Heart failure is a leading cause of cardiovascular morbidity and mortality. However, the contribution of common genetic variation to heart failure risk has not been fully elucidated, particularly in comparison to other common cardiometabolic traits. We report a multi-ancestry genome-wide association study meta-analysis of all-cause heart failure including up to 115,150 cases and 1,550,331 controls of diverse genetic ancestry, identifying 47 risk loci. We also perform multivariate genome-wide association studies that integrate heart failure with related cardiac magnetic resonance imaging endophenotypes, identifying 61 risk loci. Gene-prioritization analyses including colocalization and transcriptome-wide association studies identify known and previously unreported candidate cardiomyopathy genes and cellular processes, which we validate in gene-expression profiling of failing and healthy human hearts. Colocalization, gene expression profiling, and Mendelian randomization provide convergent evidence for the roles of BCKDHA and circulating branch-chain amino acids in heart failure and cardiac structure. Finally, proteome-wide Mendelian randomization identifies 9 circulating proteins associated with heart failure or quantitative imaging traits. These analyses highlight similarities and differences among heart failure and associated cardiovascular imaging endophenotypes, implicate common genetic variation in the pathogenesis of heart failure, and identify circulating proteins that may represent cardiomyopathy treatment targets.

Cross-Ancestry Investigation of Venous Thromboembolism Genomic Predictors.

BACKGROUND: Venous thromboembolism (VTE) is a life-threatening vascular event with environmental and genetic determinants. Recent VTE genome-wide association studies (GWAS) meta-analyses involved nearly 30 000 VTE cases and identified up to 40 genetic loci associated with VTE risk, including loci not previously suspected to play a role in hemostasis. The aim of our research was to expand discovery of new genetic loci associated with VTE by using cross-ancestry genomic resources. METHODS: We present new cross-ancestry meta-analyzed GWAS results involving up to 81 669 VTE cases from 30 studies, with replication of novel loci in independent populations and loci characterization through in silico genomic interrogations. RESULTS: In our genetic discovery effort that included 55 330 participants with VTE (47 822 European, 6320 African, and 1188 Hispanic ancestry), we identified 48 novel associations, of which 34 were replicated after correction for multiple testing. In our combined discovery-replication analysis (81 669 VTE participants) and ancestry-stratified meta-analyses (European, African, and Hispanic), we identified another 44 novel associations, which are new candidate VTE-associated loci requiring replication. In total, across all GWAS meta-analyses, we identified 135 independent genomic loci significantly associated with VTE risk. A genetic risk score of the significantly associated loci in Europeans identified a 6-fold increase in risk for those in the top 1% of scores compared with those with average scores. We also identified 31 novel transcript associations in transcriptome-wide association studies and 8 novel candidate genes with protein quantitative-trait locus Mendelian randomization analyses. In silico interrogations of hemostasis and hematology traits and a large phenome-wide association analysis of the 135 GWAS loci provided insights to biological pathways contributing to VTE, with some loci contributing to VTE through well-characterized coagulation pathways and others providing new data on the role of hematology traits, particularly platelet function. Many of the replicated loci are outside of known or currently hypothesized pathways to thrombosis. CONCLUSIONS: Our cross-ancestry GWAS meta-analyses identified new loci associated with VTE. These findings highlight new pathways to thrombosis and provide novel molecules that may be useful in the development of improved antithrombosis treatments.

Stroke genetics informs drug discovery and risk prediction across ancestries.

Previous genome-wide association studies (GWASs) of stroke - the second leading cause of death worldwide - were conducted predominantly in populations of European ancestry1,2. Here, in cross-ancestry GWAS meta-analyses of 110,182 patients who have had a stroke (five ancestries, 33% non-European) and 1,503,898 control individuals, we identify association signals for stroke and its subtypes at 89 (61 new) independent loci: 60 in primary inverse-variance-weighted analyses and 29 in secondary meta-regression and multitrait analyses. On the basis of internal cross-ancestry validation and an independent follow-up in 89,084 additional cases of stroke (30% non-European) and 1,013,843 control individuals, 87% of the primary stroke risk loci and 60% of the secondary stroke risk loci were replicated (P < 0.05). Effect sizes were highly correlated across ancestries. Cross-ancestry fine-mapping, in silico mutagenesis analysis3, and transcriptome-wide and proteome-wide association analyses revealed putative causal genes (such as SH3PXD2A and FURIN) and variants (such as at GRK5 and NOS3). Using a three-pronged approach4, we provide genetic evidence for putative drug effects, highlighting F11, KLKB1, PROC, GP1BA, LAMC2 and VCAM1 as possible targets, with drugs already under investigation for stroke for F11 and PROC. A polygenic score integrating cross-ancestry and ancestry-specific stroke GWASs with vascular-risk factor GWASs (integrative polygenic scores) strongly predicted ischaemic stroke in populations of European, East Asian and African ancestry5. Stroke genetic risk scores were predictive of ischaemic stroke independent of clinical risk factors in 52,600 clinical-trial participants with cardiometabolic disease. Our results provide insights to inform biology, reveal potential drug targets and derive genetic risk prediction tools across ancestries.