Phenome-wide identification of therapeutic genetic targets, leveraging knowledge graphs, graph neural networks, and UK Biobank data.

The ongoing expansion of human genomic datasets propels therapeutic target identification; however, extracting gene-disease associations from gene annotations remains challenging. Here, we introduce Mantis-ML 2.0, a framework integrating AstraZeneca's Biological Insights Knowledge Graph and numerous tabular datasets, to assess gene-disease probabilities throughout the phenome. We use graph neural networks, capturing the graph's holistic structure, and train them on hundreds of balanced datasets via a robust semi-supervised learning framework to provide gene-disease probabilities across the human exome. Mantis-ML 2.0 incorporates natural language processing to automate disease-relevant feature selection for thousands of diseases. The enhanced models demonstrate a 6.9% average classification power boost, achieving a median receiver operating characteristic (ROC) area under curve (AUC) score of 0.90 across 5220 diseases from Human Phenotype Ontology, OpenTargets, and Genomics England. Notably, Mantis-ML 2.0 prioritizes associations from an independent UK Biobank phenome-wide association study (PheWAS), providing a stronger form of triaging and mitigating against underpowered PheWAS associations. Results are exposed through an interactive web resource.

Rare variant associations with plasma protein levels in the UK Biobank.

Integrating human genomics and proteomics can help elucidate disease mechanisms, identify clinical biomarkers and discover drug targets1-4. Because previous proteogenomic studies have focused on common variation via genome-wide association studies, the contribution of rare variants to the plasma proteome remains largely unknown. Here we identify associations between rare protein-coding variants and 2,923 plasma protein abundances measured in 49,736 UK Biobank individuals. Our variant-level exome-wide association study identified 5,433 rare genotype-protein associations, of which 81% were undetected in a previous genome-wide association study of the same cohort5. We then looked at aggregate signals using gene-level collapsing analysis, which revealed 1,962 gene-protein associations. Of the 691 gene-level signals from protein-truncating variants, 99.4% were associated with decreased protein levels. STAB1 and STAB2, encoding scavenger receptors involved in plasma protein clearance, emerged as pleiotropic loci, with 77 and 41 protein associations, respectively. We demonstrate the utility of our publicly accessible resource through several applications. These include detailing an allelic series in NLRC4, identifying potential biomarkers for a fatty liver disease-associated variant in HSD17B13 and bolstering phenome-wide association studies by integrating protein quantitative trait loci with protein-truncating variants in collapsing analyses. Finally, we uncover distinct proteomic consequences of clonal haematopoiesis (CH), including an association between TET2-CH and increased FLT3 levels. Our results highlight a considerable role for rare variation in plasma protein abundance and the value of proteogenomics in therapeutic discovery.

Genetic Determinants of the Interventricular Septum Are Linked to Ventricular Septal Defects and Hypertrophic Cardiomyopathy.

BACKGROUND: A large proportion of genetic risk remains unexplained for structural heart disease involving the interventricular septum (IVS) including hypertrophic cardiomyopathy and ventricular septal defects. This study sought to develop a reproducible proxy of IVS structure from standard medical imaging, discover novel genetic determinants of IVS structure, and relate these loci to diseases of the IVS, hypertrophic cardiomyopathy, and ventricular septal defect. METHODS: We estimated the cross-sectional area of the IVS from the 4-chamber view of cardiac magnetic resonance imaging in 32 219 individuals from the UK Biobank which was used as the basis of genome wide association studies and Mendelian randomization. RESULTS: Measures of IVS cross-sectional area at diastole were a strong proxy for the 3-dimensional volume of the IVS (Pearson r=0.814, P=0.004), and correlated with anthropometric measures, blood pressure, and diagnostic codes related to cardiovascular physiology. Seven loci with clear genomic consequence and relevance to cardiovascular biology were uncovered by genome wide association studies, most notably a single nucleotide polymorphism in an intron of CDKN1A (rs2376620; ß, 7.7 mm2 [95% CI, 5.8-11.0]; P=6.0×10-10), and a common inversion incorporating KANSL1 predicted to disrupt local chromatin structure (ß, 8.4 mm2 [95% CI, 6.3-10.9]; P=4.2×10-14). Mendelian randomization suggested that inheritance of larger IVS cross-sectional area at diastole was strongly associated with hypertrophic cardiomyopathy risk (pIVW=4.6×10-10) while inheritance of smaller IVS cross-sectional area at diastole was associated with risk for ventricular septal defect (pIVW=0.007). CONCLUSIONS: Automated estimates of cross-sectional area of the IVS supports discovery of novel loci related to cardiac development and Mendelian disease. Inheritance of genetic liability for either small or large IVS, appears to confer risk for ventricular septal defect or hypertrophic cardiomyopathy, respectively. These data suggest that a proportion of risk for structural and congenital heart disease can be localized to the common genetic determinants of size and shape of cardiovascular anatomy.

Genetic architecture of spatial electrical biomarkers for cardiac arrhythmia and relationship with cardiovascular disease.

The 3-dimensional spatial and 2-dimensional frontal QRS-T angles are measures derived from the vectorcardiogram. They are independent risk predictors for arrhythmia, but the underlying biology is unknown. Using multi-ancestry genome-wide association studies we identify 61 (58 previously unreported) loci for the spatial QRS-T angle (N = 118,780) and 11 for the frontal QRS-T angle (N = 159,715). Seven out of the 61 spatial QRS-T angle loci have not been reported for other electrocardiographic measures. Enrichments are observed in pathways related to cardiac and vascular development, muscle contraction, and hypertrophy. Pairwise genome-wide association studies with classical ECG traits identify shared genetic influences with PR interval and QRS duration. Phenome-wide scanning indicate associations with atrial fibrillation, atrioventricular block and arterial embolism and genetically determined QRS-T angle measures are associated with fascicular and bundle branch block (and also atrioventricular block for the frontal QRS-T angle). We identify potential biology involved in the QRS-T angle and their genetic relationships with cardiovascular traits and diseases, may inform future research and risk prediction.

Large scale genome-wide association analyses identify novel genetic loci and mechanisms in hypertrophic cardiomyopathy.

Hypertrophic cardiomyopathy (HCM) is an important cause of morbidity and mortality with both monogenic and polygenic components. We here report results from the largest HCM genome-wide association study (GWAS) and multi-trait analysis (MTAG) including 5,900 HCM cases, 68,359 controls, and 36,083 UK Biobank (UKB) participants with cardiac magnetic resonance (CMR) imaging. We identified a total of 70 loci (50 novel) associated with HCM, and 62 loci (32 novel) associated with relevant left ventricular (LV) structural or functional traits. Amongst the common variant HCM loci, we identify a novel HCM disease gene, SVIL, which encodes the actin-binding protein supervillin, showing that rare truncating SVIL variants cause HCM. Mendelian randomization analyses support a causal role of increased LV contractility in both obstructive and non-obstructive forms of HCM, suggesting common disease mechanisms and anticipating shared response to therapy. Taken together, the findings significantly increase our understanding of the genetic basis and molecular mechanisms of HCM, with potential implications for disease management.

Effects of protein-coding variants on blood metabolite measurements and clinical biomarkers in the UK Biobank.

Genome-wide association studies (GWASs) have established the contribution of common and low-frequency variants to metabolic blood measurements in the UK Biobank (UKB). To complement existing GWAS findings, we assessed the contribution of rare protein-coding variants in relation to 355 metabolic blood measurements-including 325 predominantly lipid-related nuclear magnetic resonance (NMR)-derived blood metabolite measurements (Nightingale Health Plc) and 30 clinical blood biomarkers-using 412,393 exome sequences from four genetically diverse ancestries in the UKB. Gene-level collapsing analyses were conducted to evaluate a diverse range of rare-variant architectures for the metabolic blood measurements. Altogether, we identified significant associations (p < 1 × 10-8) for 205 distinct genes that involved 1,968 significant relationships for the Nightingale blood metabolite measurements and 331 for the clinical blood biomarkers. These include associations for rare non-synonymous variants in PLIN1 and CREB3L3 with lipid metabolite measurements and SYT7 with creatinine, among others, which may not only provide insights into novel biology but also deepen our understanding of established disease mechanisms. Of the study-wide significant clinical biomarker associations, 40% were not previously detected on analyzing coding variants in a GWAS in the same cohort, reinforcing the importance of studying rare variation to fully understand the genetic architecture of metabolic blood measurements.

Genome-Wide Analysis of Left Ventricular Maximum Wall Thickness in the UK Biobank Cohort Reveals a Shared Genetic Background With Hypertrophic Cardiomyopathy.

BACKGROUND: Left ventricular maximum wall thickness (LVMWT) is an important biomarker of left ventricular hypertrophy and provides diagnostic and prognostic information in hypertrophic cardiomyopathy (HCM). Limited information is available on the genetic determinants of LVMWT. METHODS: We performed a genome-wide association study of LVMWT measured from the cardiovascular magnetic resonance examinations of 42 176 European individuals. We evaluated the genetic relationship between LVMWT and HCM by performing pairwise analysis using the data from the Hypertrophic Cardiomyopathy Registry in which the controls were randomly selected from UK Biobank individuals not included in the cardiovascular magnetic resonance sub-study. RESULTS: Twenty-one genetic loci were discovered at P<5×10-8. Several novel candidate genes were identified including PROX1, PXN, and PTK2, with known functional roles in myocardial growth and sarcomere organization. The LVMWT genetic risk score is predictive of HCM in the Hypertrophic Cardiomyopathy Registry (odds ratio per SD: 1.18 [95% CI, 1.13-1.23]) with pairwise analyses demonstrating a moderate genetic correlation (rg=0.53) and substantial loci overlap (19/21). CONCLUSIONS: Our findings provide novel insights into the genetic underpinning of LVMWT and highlight its shared genetic background with HCM, supporting future endeavours to elucidate the genetic etiology of HCM.

DrugnomeAI is an ensemble machine-learning framework for predicting druggability of candidate drug targets.

The druggability of targets is a crucial consideration in drug target selection. Here, we adopt a stochastic semi-supervised ML framework to develop DrugnomeAI, which estimates the druggability likelihood for every protein-coding gene in the human exome. DrugnomeAI integrates gene-level properties from 15 sources resulting in 324 features. The tool generates exome-wide predictions based on labelled sets of known drug targets (median AUC: 0.97), highlighting features from protein-protein interaction networks as top predictors. DrugnomeAI provides generic as well as specialised models stratified by disease type or drug therapeutic modality. The top-ranking DrugnomeAI genes were significantly enriched for genes previously selected for clinical development programs (p value < 1 × 10-308) and for genes achieving genome-wide significance in phenome-wide association studies of 450 K UK Biobank exomes for binary (p value = 1.7 × 10-5) and quantitative traits (p value = 1.6 × 10-7). We accompany our method with a web application ( http://drugnomeai.public.cgr.astrazeneca.com ) to visualise the druggability predictions and the key features that define gene druggability, per disease type and modality.

Human genetics uncovers MAP3K15 as an obesity-independent therapeutic target for diabetes.

We performed collapsing analyses on 454,796 UK Biobank (UKB) exomes to detect gene-level associations with diabetes. Recessive carriers of nonsynonymous variants in MAP3K15 were 30% less likely to develop diabetes (P = 5.7 × 10-10) and had lower glycosylated hemoglobin (ß = -0.14 SD units, P = 1.1 × 10-24). These associations were independent of body mass index, suggesting protection against insulin resistance even in the setting of obesity. We replicated these findings in 96,811 Admixed Americans in the Mexico City Prospective Study (P < 0.05)Moreover, the protective effect of MAP3K15 variants was stronger in individuals who did not carry the Latino-enriched SLC16A11 risk haplotype (P = 6.0 × 10-4). Separately, we identified a Finnish-enriched MAP3K15 protein-truncating variant associated with decreased odds of both type 1 and type 2 diabetes (P < 0.05) in FinnGen. No adverse phenotypes were associated with protein-truncating MAP3K15 variants in the UKB, supporting this gene as a therapeutic target for diabetes.

Cancer-driving mutations are enriched in genic regions intolerant to germline variation.

Large reference datasets of protein-coding variation in human populations have allowed us to determine which genes and genic subregions are intolerant to germline genetic variation. There is also a growing number of genes implicated in severe Mendelian diseases that overlap with genes implicated in cancer. We hypothesized that cancer-driving mutations might be enriched in genic subregions that are depleted of germline variation relative to somatic variation. We introduce a new metric, OncMTR (oncology missense tolerance ratio), which uses 125,748 exomes in the Genome Aggregation Database (gnomAD) to identify these genic subregions. We demonstrate that OncMTR can significantly predict driver mutations implicated in hematologic malignancies. Divergent OncMTR regions were enriched for cancer-relevant protein domains, and overlaying OncMTR scores on protein structures identified functionally important protein residues. Last, we performed a rare variant, gene-based collapsing analysis on an independent set of 394,694 exomes from the UK Biobank and find that OncMTR markedly improves genetic signals for hematologic malignancies.

Gene-SCOUT: identifying genes with similar continuous trait fingerprints from phenome-wide association analyses.

Large-scale phenome-wide association studies performed using densely-phenotyped cohorts such as the UK Biobank (UKB), reveal many statistically robust gene-phenotype relationships for both clinical and continuous traits. Here, we present Gene-SCOUT, a tool used to identify genes with similar continuous trait fingerprints to a gene of interest. A fingerprint reflects the continuous traits identified to be statistically associated with a gene of interest based on multiple underlying rare variant genetic architectures. Similarities between genes are evaluated by the cosine similarity measure, to capture concordant effect directionality, elucidating clusters of genes in a high dimensional space. The underlying gene-biomarker population-scale association statistics were obtained from a gene-level rare variant collapsing analysis performed on over 1500 continuous traits using 394 692 UKB participant exomes, with additional metabolomic trait associations provided through Nightingale Health's recent study of 121 394 of these participants. We demonstrate that gene similarity estimates from Gene-SCOUT provide stronger enrichments for clinical traits compared to existing methods. Furthermore, we provide a fully interactive web-resource (http://genescout.public.cgr.astrazeneca.com) to explore the pre-calculated exome-wide similarities. This resource enables a user to examine the biological relevance of the most similar genes for Gene Ontology (GO) enrichment and UKB clinical trait enrichment statistics, as well as a detailed breakdown of the traits underpinning a given fingerprint.

Incremental value of left atrial booster and reservoir strain in predicting atrial fibrillation in patients with hypertrophic cardiomyopathy: a cardiovascular magnetic resonance study.

BACKGROUND: Left atrial (LA) size and function are known predictors of new onset atrial fibrillation (AF) in hypertrophic cardiomyopathy (HCM) patients. Components of LA deformation including reservoir, conduit, and booster function provide additional information on atrial mechanics. Whether or not LA deformation can augment our ability to predict the risk of new onset AF in HCM patients beyond standard measurements is unknown. METHODS: We assessed LA size, function, and deformation on cardiovascular magnetic resonance (CMR) in 238 genotyped HCM patients and compared this with twenty age, sex, blood pressure and body mass index matched control subjects. We further evaluated the determinants of new onset AF in HCM patients. RESULTS: Compared to control subjects, HCM patients had higher LA antero-posterior diameter, lower LA ejection fraction and lower LA reservoir (19.9 [17.1, 22.2], 21.6 [19.9, 22.9], P = 0.047) and conduit strain (10.6 ± 4.4, 13.7 ± 3.3, P = 0.002). LA booster strain did not differ between healthy controls and HCM patients, but HCM patients who developed new onset AF (n = 33) had lower booster strain (7.6 ± 3.3, 9.5 ± 3.0, P = 0.001) than those that did not (n = 205). In separate multivariate models, age, LA ejection fraction, and LA booster and reservoir strain were each independent determinants of AF. Age ≥ 55 years was the strongest determinant (HR 6.62, 95% CI 2.79-15.70), followed by LA booster strain ≤ 8% (HR 3.69, 95% CI 1.81-7.52) and LA reservoir strain ≤ 18% (HR 2.56, 95% CI 1.24-5.27). Conventional markers of HCM phenotypic severity, age and sudden death risk factors were associated with LA strain components. CONCLUSIONS: LA strain components are impaired in HCM and, together with age, independently predicted the risk of new onset AF. Increasing age and phenotypic severity were associated with LA strain abnormalities. Our findings suggest that the routine assessment of LA strain components and consideration of age could augment LA size in predicting risk of AF, and potentially guide prophylactic anticoagulation use in HCM.

Rare variant contribution to human disease in 281,104 UK Biobank exomes.

Genome-wide association studies have uncovered thousands of common variants associated with human disease, but the contribution of rare variants to common disease remains relatively unexplored. The UK Biobank contains detailed phenotypic data linked to medical records for approximately 500,000 participants, offering an unprecedented opportunity to evaluate the effect of rare variation on a broad collection of traits1,2. Here we study the relationships between rare protein-coding variants and 17,361 binary and 1,419 quantitative phenotypes using exome sequencing data from 269,171 UK Biobank participants of European ancestry. Gene-based collapsing analyses revealed 1,703 statistically significant gene-phenotype associations for binary traits, with a median odds ratio of 12.4. Furthermore, 83% of these associations were undetectable via single-variant association tests, emphasizing the power of gene-based collapsing analysis in the setting of high allelic heterogeneity. Gene-phenotype associations were also significantly enriched for loss-of-function-mediated traits and approved drug targets. Finally, we performed ancestry-specific and pan-ancestry collapsing analyses using exome sequencing data from 11,933 UK Biobank participants of African, East Asian or South Asian ancestry. Our results highlight a significant contribution of rare variants to common disease. Summary statistics are publicly available through an interactive portal ( http://azphewas.com/ ).

Shared genetic pathways contribute to risk of hypertrophic and dilated cardiomyopathies with opposite directions of effect.

The heart muscle diseases hypertrophic (HCM) and dilated (DCM) cardiomyopathies are leading causes of sudden death and heart failure in young, otherwise healthy, individuals. We conducted genome-wide association studies and multi-trait analyses in HCM (1,733 cases), DCM (5,521 cases) and nine left ventricular (LV) traits (19,260 UK Biobank participants with structurally normal hearts). We identified 16 loci associated with HCM, 13 with DCM and 23 with LV traits. We show strong genetic correlations between LV traits and cardiomyopathies, with opposing effects in HCM and DCM. Two-sample Mendelian randomization supports a causal association linking increased LV contractility with HCM risk. A polygenic risk score explains a significant portion of phenotypic variability in carriers of HCM-causing rare variants. Our findings thus provide evidence that polygenic risk score may account for variability in Mendelian diseases. More broadly, we provide insights into how genetic pathways may lead to distinct disorders through opposing genetic effects.

Common genetic variants and modifiable risk factors underpin hypertrophic cardiomyopathy susceptibility and expressivity.

Hypertrophic cardiomyopathy (HCM) is a common, serious, genetic heart disorder. Rare pathogenic variants in sarcomere genes cause HCM, but with unexplained phenotypic heterogeneity. Moreover, most patients do not carry such variants. We report a genome-wide association study of 2,780 cases and 47,486 controls that identified 12 genome-wide-significant susceptibility loci for HCM. Single-nucleotide polymorphism heritability indicated a strong polygenic influence, especially for sarcomere-negative HCM (64% of cases; h2g = 0.34 ± 0.02). A genetic risk score showed substantial influence on the odds of HCM in a validation study, halving the odds in the lowest quintile and doubling them in the highest quintile, and also influenced phenotypic severity in sarcomere variant carriers. Mendelian randomization identified diastolic blood pressure (DBP) as a key modifiable risk factor for sarcomere-negative HCM, with a one standard deviation increase in DBP increasing the HCM risk fourfold. Common variants and modifiable risk factors have important roles in HCM that we suggest will be clinically actionable.

Secondary findings in inherited heart conditions: a genotype-first feasibility study to assess phenotype, behavioural and psychosocial outcomes.

Disclosing secondary findings (SF) from genome sequencing (GS) can alert carriers to disease risk. However, evidence around variant-disease association and consequences of disclosure for individuals and healthcare services is limited. We report on the feasibility of an approach to identification of SF in inherited cardiac conditions (ICC) genes in participants in a rare disease GS study, followed by targeted clinical evaluation. Qualitative methods were used to explore behavioural and psychosocial consequences of disclosure. ICC genes were analysed in genome sequence data from 7203 research participants; a two-stage approach was used to recruit genotype-blind variant carriers and matched controls. Cardiac-focused medical and family history collection and genetic counselling were followed by standard clinical tests, blinded to genotype. Pathogenic ICC variants were identified in 0.61% of individuals; 20 were eligible for the present study. Four variant carriers and seven non-carrier controls participated. One variant carrier had a family history of ICC and was clinically affected; a second was clinically unaffected and had no relevant family history. One variant, in two unrelated participants, was subsequently reclassified as being of uncertain significance. Analysis of qualitative data highlights participant satisfaction with approach, willingness to follow clinical recommendations, but variable outcomes of relatives' engagement with healthcare services. In conclusion, when offered access to SF, many people choose not to pursue them. For others, disclosure of ICC SF in a specialist setting is valued and of likely clinical utility, and can be expected to identify individuals with, and without a phenotype.

Reevaluation of the South Asian MYBPC3Δ25bp Intronic Deletion in Hypertrophic Cardiomyopathy.

BACKGROUND: The common intronic deletion, MYBPC3Δ25, detected in 4% to 8% of South Asian populations, is reported to be associated with cardiomyopathy, with ≈7-fold increased risk of disease in variant carriers. Here, we examine the contribution of MYBPC3Δ25 to hypertrophic cardiomyopathy (HCM) in a large patient cohort. METHODS: Sequence data from 2 HCM cohorts (n=5393) was analyzed to determine MYBPC3Δ25 frequency and co-occurrence of pathogenic variants in HCM genes. Case-control and haplotype analyses were performed to compare variant frequencies and assess disease association. Analyses were also undertaken to investigate the pathogenicity of a candidate variant MYBPC3 c.1224-52G>A. RESULTS: Our data suggest that the risk of HCM, previously attributed to MYBPC3Δ25, can be explained by enrichment of a derived haplotype, MYBPC3Δ25/-52, whereby a small subset of individuals bear both MYBPC3Δ25 and a rare pathogenic variant, MYBPC3 c.1224-52G>A. The intronic MYBPC3 c.1224-52G>A variant, which is not routinely evaluated by gene panel or exome sequencing, was detected in ≈1% of our HCM cohort. CONCLUSIONS: The MYBPC3 c.1224-52G>A variant explains the disease risk previously associated with MYBPC3Δ25 in the South Asian population and is one of the most frequent pathogenic variants in HCM in all populations; genotyping services should ensure coverage of this deep intronic mutation. Individuals carrying MYBPC3Δ25 alone are not at increased risk of HCM, and this variant should not be tested in isolation; this is important for the large majority of the 100 million individuals of South Asian ancestry who carry MYBPC3Δ25 and would previously have been declared at increased risk of HCM.

Distinct Subgroups in Hypertrophic Cardiomyopathy in the NHLBI HCM Registry.

BACKGROUND: The HCMR (Hypertrophic Cardiomyopathy Registry) is a National Heart, Lung, and Blood Institute-funded, prospective registry of 2,755 patients with hypertrophic cardiomyopathy (HCM) recruited from 44 sites in 6 countries. OBJECTIVES: The authors sought to improve risk prediction in HCM by incorporating cardiac magnetic resonance (CMR), genetic, and biomarker data. METHODS: Demographic and echocardiographic data were collected. Patients underwent CMR including cine imaging, late gadolinium enhancement imaging (LGE) (replacement fibrosis), and T1 mapping for measurement of extracellular volume as a measure of interstitial fibrosis. Blood was drawn for the biomarkers N-terminal pro-B-type natriuretic peptide (NT-proBNP) and high-sensitivity cardiac troponin T (cTnT), and genetic analysis. RESULTS: A total of 2,755 patients were studied. Mean age was 49 ± 11 years, 71% were male, and 17% non-white. Mean ESC (European Society of Cardiology) risk score was 2.48 ± 0.56. Eighteen percent had a resting left ventricular outflow tract (LVOT) gradient ≥30 mm Hg. Thirty-six percent had a sarcomere mutation identified, and 50% had any LGE. Sarcomere mutation-positive patients were more likely to have reverse septal curvature morphology, LGE, and no significant resting LVOT obstruction. Those that were sarcomere mutation negative were more likely to have isolated basal septal hypertrophy, less LGE, and more LVOT obstruction. Interstitial fibrosis was present in segments both with and without LGE. Serum NT-proBNP and cTnT levels correlated with increasing LGE and extracellular volume in a graded fashion. CONCLUSIONS: The HCMR population has characteristics of low-risk HCM. Ninety-three percent had no or only mild functional limitation. Baseline data separated patients broadly into 2 categories. One group was sarcomere mutation positive and more likely had reverse septal curvature morphology, more fibrosis, but less resting obstruction, whereas the other was sarcomere mutation negative and more likely had isolated basal septal hypertrophy with obstruction, but less fibrosis. Further follow-up will allow better understanding of these subgroups and development of an improved risk prediction model incorporating all these markers.

Analysis of 51 proposed hypertrophic cardiomyopathy genes from genome sequencing data in sarcomere negative cases has negligible diagnostic yield.

PURPOSE: Increasing numbers of genes are being implicated in Mendelian disorders and incorporated into clinical test panels. However, lack of evidence supporting the gene-disease relationship can hinder interpretation. We explored the utility of testing 51 additional genes for hypertrophic cardiomyopathy (HCM), one of the most commonly tested Mendelian disorders. METHODS: Using genome sequencing data from 240 sarcomere gene negative HCM cases and 6229 controls, we undertook case-control and individual variant analyses to assess 51 genes that have been proposed for HCM testing. RESULTS: We found no evidence to suggest that rare variants in these genes are prevalent causes of HCM. One variant, in a single case, was categorized as likely to be pathogenic. Over 99% of variants were classified as a variant of uncertain significance (VUS) and 54% of cases had one or more VUS. CONCLUSION: For almost all genes, the gene-disease relationship could not be validated and lack of evidence precluded variant interpretation. Thus, the incremental diagnostic yield of extending testing was negligible, and would, we propose, be outweighed by problems that arise with a high rate of uninterpretable findings. These findings highlight the need for rigorous, evidence-based selection of genes for clinical test panels.