Rare coding variant analysis for human diseases across biobanks and ancestries.

Large-scale sequencing has enabled unparalleled opportunities to investigate the role of rare coding variation in human phenotypic variability. Here, we present a pan-ancestry analysis of sequencing data from three large biobanks, including the All of Us research program. Using mixed-effects models, we performed gene-based rare variant testing for 601 diseases across 748,879 individuals, including 155,236 with ancestry dissimilar to European. We identified 363 significant associations, which highlighted core genes for the human disease phenome and identified potential novel associations, including UBR3 for cardiometabolic disease and YLPM1 for psychiatric disease. Pan-ancestry burden testing represented an inclusive and useful approach for discovery in diverse datasets, although we also highlight the importance of ancestry-specific sensitivity analyses in this setting. Finally, we found that effect sizes for rare protein-disrupting variants were concordant between samples similar to European ancestry and other genetic ancestries (ßDeming = 0.7-1.0). Our results have implications for multi-ancestry and cross-biobank approaches in sequencing association studies for human disease.

Convergence of coronary artery disease genes onto endothelial cell programs.

Linking variants from genome-wide association studies (GWAS) to underlying mechanisms of disease remains a challenge1-3. For some diseases, a successful strategy has been to look for cases in which multiple GWAS loci contain genes that act in the same biological pathway1-6. However, our knowledge of which genes act in which pathways is incomplete, particularly for cell-type-specific pathways or understudied genes. Here we introduce a method to connect GWAS variants to functions. This method links variants to genes using epigenomics data, links genes to pathways de novo using Perturb-seq and integrates these data to identify convergence of GWAS loci onto pathways. We apply this approach to study the role of endothelial cells in genetic risk for coronary artery disease (CAD), and discover 43 CAD GWAS signals that converge on the cerebral cavernous malformation (CCM) signalling pathway. Two regulators of this pathway, CCM2 and TLNRD1, are each linked to a CAD risk variant, regulate other CAD risk genes and affect atheroprotective processes in endothelial cells. These results suggest a model whereby CAD risk is driven in part by the convergence of causal genes onto a particular transcriptional pathway in endothelial cells. They highlight shared genes between common and rare vascular diseases (CAD and CCM), and identify TLNRD1 as a new, previously uncharacterized member of the CCM signalling pathway. This approach will be widely useful for linking variants to functions for other common polygenic diseases.

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.

A multi-ancestry polygenic risk score improves risk prediction for coronary artery disease.

Identification of individuals at highest risk of coronary artery disease (CAD)-ideally before onset-remains an important public health need. Prior studies have developed genome-wide polygenic scores to enable risk stratification, reflecting the substantial inherited component to CAD risk. Here we develop a new and significantly improved polygenic score for CAD, termed GPSMult, that incorporates genome-wide association data across five ancestries for CAD (>269,000 cases and >1,178,000 controls) and ten CAD risk factors. GPSMult strongly associated with prevalent CAD (odds ratio per standard deviation 2.14, 95% confidence interval 2.10-2.19, P < 0.001) in UK Biobank participants of European ancestry, identifying 20.0% of the population with 3-fold increased risk and conversely 13.9% with 3-fold decreased risk as compared with those in the middle quintile. GPSMult was also associated with incident CAD events (hazard ratio per standard deviation 1.73, 95% confidence interval 1.70-1.76, P < 0.001), identifying 3% of healthy individuals with risk of future CAD events equivalent to those with existing disease and significantly improving risk discrimination and reclassification. Across multiethnic, external validation datasets inclusive of 33,096, 124,467, 16,433 and 16,874 participants of African, European, Hispanic and South Asian ancestry, respectively, GPSMult demonstrated increased strength of associations across all ancestries and outperformed all available previously published CAD polygenic scores. These data contribute a new GPSMult for CAD to the field and provide a generalizable framework for how large-scale integration of genetic association data for CAD and related traits from diverse populations can meaningfully improve polygenic risk prediction.

Genome-wide association analysis and Mendelian randomization proteomics identify drug targets for heart failure.

We conduct a large-scale meta-analysis of heart failure genome-wide association studies (GWAS) consisting of over 90,000 heart failure cases and more than 1 million control individuals of European ancestry to uncover novel genetic determinants for heart failure. Using the GWAS results and blood protein quantitative loci, we perform Mendelian randomization and colocalization analyses on human proteins to provide putative causal evidence for the role of druggable proteins in the genesis of heart failure. We identify 39 genome-wide significant heart failure risk variants, of which 18 are previously unreported. Using a combination of Mendelian randomization proteomics and genetic cis-only colocalization analyses, we identify 10 additional putatively causal genes for heart failure. Findings from GWAS and Mendelian randomization-proteomics identify seven (CAMK2D, PRKD1, PRKD3, MAPK3, TNFSF12, APOC3 and NAE1) proteins as potential targets for interventions to be used in primary prevention of heart failure.

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.

Genome-wide association meta-analysis of spontaneous coronary artery dissection identifies risk variants and genes related to artery integrity and tissue-mediated coagulation.

Spontaneous coronary artery dissection (SCAD) is an understudied cause of myocardial infarction primarily affecting women. It is not known to what extent SCAD is genetically distinct from other cardiovascular diseases, including atherosclerotic coronary artery disease (CAD). Here we present a genome-wide association meta-analysis (1,917 cases and 9,292 controls) identifying 16 risk loci for SCAD. Integrative functional annotations prioritized genes that are likely to be regulated in vascular smooth muscle cells and artery fibroblasts and implicated in extracellular matrix biology. One locus containing the tissue factor gene F3, which is involved in blood coagulation cascade initiation, appears to be specific for SCAD risk. Several associated variants have diametrically opposite associations with CAD, suggesting that shared biological processes contribute to both diseases, but through different mechanisms. We also infer a causal role for high blood pressure in SCAD. Our findings provide novel pathophysiological insights involving arterial integrity and tissue-mediated coagulation in SCAD and set the stage for future specific therapeutics and preventions.

Clinical and genetic associations of deep learning-derived cardiac magnetic resonance-based left ventricular mass.

Left ventricular mass is a risk marker for cardiovascular events, and may indicate an underlying cardiomyopathy. Cardiac magnetic resonance is the gold-standard for left ventricular mass estimation, but is challenging to obtain at scale. Here, we use deep learning to enable genome-wide association study of cardiac magnetic resonance-derived left ventricular mass indexed to body surface area within 43,230 UK Biobank participants. We identify 12 genome-wide associations (1 known at TTN and 11 novel for left ventricular mass), implicating genes previously associated with cardiac contractility and cardiomyopathy. Cardiac magnetic resonance-derived indexed left ventricular mass is associated with incident dilated and hypertrophic cardiomyopathies, and implantable cardioverter-defibrillator implant. An indexed left ventricular mass polygenic risk score ≥90th percentile is also associated with incident implantable cardioverter-defibrillator implant in separate UK Biobank (hazard ratio 1.22, 95% CI 1.05-1.44) and Mass General Brigham (hazard ratio 1.75, 95% CI 1.12-2.74) samples. Here, we perform a genome-wide association study of cardiac magnetic resonance-derived indexed left ventricular mass to identify 11 novel variants and demonstrate that cardiac magnetic resonance-derived and genetically predicted indexed left ventricular mass are associated with incident cardiomyopathy.

Predictive Utility of a Coronary Artery Disease Polygenic Risk Score in Primary Prevention.

Importance: The clinical utility of polygenic risk scores (PRS) for coronary artery disease (CAD) has not yet been established. Objective: To investigate the ability of a CAD PRS to potentially guide statin initiation in primary prevention after accounting for age and clinical risk. Design, Setting, and Participants: This was a longitudinal cohort study with enrollment starting on January 1, 2006, and ending on December 31, 2010, with data updated to mid-2021, using data from the UK Biobank, a long-term population study of UK citizens. A replication analysis was performed in Biobank Japan. The analysis included all patients without a history of CAD and who were not taking lipid-lowering therapy. Data were analyzed from January 1 to June 30, 2022. Exposures: Polygenic risk for CAD was defined as low (bottom 20%), intermediate, and high (top 20%) using a CAD PRS including 241 genome-wide significant single-nucleotide variations (SNVs). The pooled cohort equations were used to estimate 10-year atherosclerotic cardiovascular disease (ASCVD) risk and classify individuals as low (<5%), borderline (5-<7.5%), intermediate (7.5-<20%), or high risk (≥20%). Main Outcomes and Measures: Myocardial infarction (MI) and ASCVD events (defined as incident clinical CAD [including MI], stroke, or CV death). Results: A total of 330 201 patients (median [IQR] age, 57 [40-74] years; 189 107 female individuals [57%]) were included from the UK Biobank. Over the 10-year follow-up, 4454 individuals had an MI. The CAD PRS was significantly associated with the risk of MI in all age groups but had significantly stronger risk prediction at younger ages (age <50 years: hazard ratio [HR] per 1 SD of PRS, 1.72; 95% CI, 1.56-1.89; age 50-60 years: HR, 1.46; 95% CI, 1.38-1.53; age >60 years: HR, 1.42; 95% CI, 1.37-1.48; P for interaction <.001). In patients younger than 50 years, those with high PRS had a 3- to 4-fold increased associated risk of MI compared with those in the low PRS category. A significant interaction between CAD PRS and age was replicated in Biobank Japan. When CAD PRS testing was added to the clinical ASCVD risk score in individuals younger than 50 years, 591 of 4373 patients (20%) with borderline risk were risk stratified into intermediate risk, warranting initiation of statin therapy and 3198 of 7477 patients (20%) with both borderline or intermediate risk were stratified as low risk, thus not warranting therapy. Conclusions and Relevance: Results of this cohort study suggest that the predictive ability of a CAD PRS was greater in younger individuals and can be used to better identify patients with borderline and intermediate clinical risk who should initiate statin therapy.

Genetic architecture of heart failure with preserved versus reduced ejection fraction.

Pharmacologic clinical trials for heart failure with preserved ejection fraction have been largely unsuccessful as compared to those for heart failure with reduced ejection fraction. Whether differences in the genetic underpinnings of these major heart failure subtypes may provide insights into the disparate outcomes of clinical trials remains unknown. We utilize a large, uniformly phenotyped, single cohort of heart failure sub-classified into heart failure with reduced and with preserved ejection fractions based on current clinical definitions, to conduct detailed genetic analyses of the two heart failure sub-types. We find different genetic architectures and distinct genetic association profiles between heart failure with reduced and with preserved ejection fraction suggesting differences in underlying pathobiology. The modest genetic discovery for heart failure with preserved ejection fraction (one locus) compared to heart failure with reduced ejection fraction (13 loci) despite comparable sample sizes indicates that clinically defined heart failure with preserved ejection fraction likely represents the amalgamation of several, distinct pathobiological entities. Development of consensus sub-phenotyping of heart failure with preserved ejection fraction is paramount to better dissect the underlying genetic signals and contributors to this highly prevalent condition.

Discovery and systematic characterization of risk variants and genes for coronary artery disease in over a million participants.

The discovery of genetic loci associated with complex diseases has outpaced the elucidation of mechanisms of disease pathogenesis. Here we conducted a genome-wide association study (GWAS) for coronary artery disease (CAD) comprising 181,522 cases among 1,165,690 participants of predominantly European ancestry. We detected 241 associations, including 30 new loci. Cross-ancestry meta-analysis with a Japanese GWAS yielded 38 additional new loci. We prioritized likely causal variants using functionally informed fine-mapping, yielding 42 associations with less than five variants in the 95% credible set. Similarity-based clustering suggested roles for early developmental processes, cell cycle signaling and vascular cell migration and proliferation in the pathogenesis of CAD. We prioritized 220 candidate causal genes, combining eight complementary approaches, including 123 supported by three or more approaches. Using CRISPR-Cas9, we experimentally validated the effect of an enhancer in MYO9B, which appears to mediate CAD risk by regulating vascular cell motility. Our analysis identifies and systematically characterizes >250 risk loci for CAD to inform experimental interrogation of putative causal mechanisms for CAD.

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.

Clinical Implementation of Combined Monogenic and Polygenic Risk Disclosure for Coronary Artery Disease.

BACKGROUND: State-of-the-art genetic risk interpretation for a common complex disease such as coronary artery disease (CAD) requires assessment for both monogenic variants-such as those related to familial hypercholesterolemia-as well as the cumulative impact of many common variants, as quantified by a polygenic score. OBJECTIVES: The objective of the study was to describe a combined monogenic and polygenic CAD risk assessment program and examine its impact on patient understanding and changes to clinical management. METHODS: Study participants attended an initial visit in a preventive genomics clinic and a disclosure visit to discuss results and recommendations, primarily via telemedicine. Digital postdisclosure surveys and chart review evaluated the impact of disclosure. RESULTS: There were 60 participants (mean age 51 years, 37% women, 72% with no known CAD), including 30 (50%) referred by their cardiologists and 30 (50%) self-referred. Two (3%) participants had a monogenic variant pathogenic for familial hypercholesterolemia, and 19 (32%) had a high polygenic score in the top quintile of the population distribution. In a postdisclosure survey, both the genetic test report (in 80% of participants) and the discussion with the clinician (in 89% of participants) were ranked as very or extremely helpful in understanding the result. Of the 42 participants without CAD, 17 or 40% had a change in management, including statin initiation, statin intensification, or coronary imaging. CONCLUSIONS: Combined monogenic and polygenic assessments for CAD risk provided by preventive genomics clinics are beneficial for patients and result in changes in management in a significant portion of patients.

Rare and Common Genetic Variation Underlying the Risk of Hypertrophic Cardiomyopathy in a National Biobank.

Importance: Hypertrophic cardiomyopathy (HCM) is a leading cause of sudden cardiac death in young people. Although rare genetic variants are well-established contributors to HCM risk, common genetic variants have recently been implicated in disease pathogenesis. Objective: To assess the contributions of rare and common genetic variation to risk of HCM in the general population. Design, Setting, and Participants: This cohort study of the UK Biobank (data from 2006-2010) and the Mass General Brigham Biobank (2010-2019) assessed the relative and joint contributions of rare genetic variants and a common variant (polygenic) score to risk of HCM. Both rare and common variant predictors were then evaluated in the context of relevant clinical risk factors. Data analysis was conducted from May 2021 to February 2022. Exposures: Pathogenic rare variants, common-variant (polygenic) score, and clinical risk factors. Main Outcomes and Measures: Risk of HCM. Results: The primary study population comprised 184 511 individuals from the UK Biobank. Mean (SD) age was 56 (8) years, 83 690 (45%) of participants were men, and 204 (0.1%) participants had HCM. Of 51 genes included in clinical genetic testing panels for HCM, pathogenic or likely pathogenic variants in 14 core genes (designated by the American College of Medical Genetics and Genomics [ACMG]) were associated with 55-fold higher odds (95% CI, 35-83) of HCM, while those in the remaining 37 non-ACMG genes were not significantly associated with HCM (OR, 1.8; 95% CI, 0.6-4.0). ClinVar pathogenic or likely pathogenic mutations in MYBPC3 (OR, 72; 95% CI, 39-124) and MYH7 (OR, 61; 95% CI, 26-121) were strongly associated with HCM, as were loss-of-function variants in ALPK3 (OR, 13; 95% CI, 4.4-28). A polygenic score was strongly associated with HCM (OR per SD increase in score, 1.6; 95% CI, 1.4-1.8), with concordant results in the Mass General Brigham Biobank. Genetic factors enhanced clinical risk prediction for HCM: addition of rare variant carrier status and the polygenic score to clinical risk factors (obesity, hypertension, atrial fibrillation, and coronary artery disease) improved the area under the receiver operator characteristic curve from 0.71 (95% CI, 0.65-0.77) to 0.82 (95% CI, 0.77-0.87). Conclusions and Relevance: Both rare and common genetic variants contribute substantially to HCM susceptibility in the general population and improve HCM risk prediction beyond that achieved with clinical factors.

Genome-wide and phenome-wide analysis of ideal cardiovascular health in the VA Million Veteran Program.

BACKGROUND: Genetic studies may help identify causal pathways; therefore, we sought to identify genetic determinants of ideal CVH and their association with CVD outcomes in the multi-population Veteran Administration Million Veteran Program. METHODS: An ideal health score (IHS) was calculated from 3 clinical factors (blood pressure, total cholesterol, and blood glucose levels) and 3 behavioral factors (smoking status, physical activity, and BMI), ascertained at baseline. Multi-population genome-wide association study (GWAS) was performed on IHS and binary ideal health using linear and logistic regression, respectively. Using the genome-wide significant SNPs from the IHS GWAS, we created a weighted IHS polygenic risk score (PRSIHS) which was used (i) to conduct a phenome-wide association study (PheWAS) of associations between PRSIHS and ICD-9 phenotypes and (ii) to further test for associations with mortality and selected CVD outcomes using logistic and Cox regression and, as an instrumental variable, in Mendelian Randomization. RESULTS: The discovery and replication cohorts consisted of 142,404 (119,129 European American (EUR); 16,495 African American (AFR)), and 45,766 (37,646 EUR; 5,366 AFR) participants, respectively. The mean age was 65.8 years (SD = 11.2) and 92.7% were male. Overall, 4.2% exhibited ideal CVH based on the clinical and behavioral factors. In the multi-population meta-analysis, variants at 17 loci were associated with IHS and each had known GWAS associations with multiple components of the IHS. PheWAS analysis in 456,026 participants showed that increased PRSIHS was associated with a lower odds ratio for many CVD outcomes and risk factors. Both IHS and PRSIHS measures of ideal CVH were associated with significantly less CVD outcomes and CVD mortality. CONCLUSION: A set of high interest genetic variants contribute to the presence of ideal CVH in a multi-ethnic cohort of US Veterans. Genetically influenced ideal CVH is associated with lower odds of CVD outcomes and mortality.