Reduction of Filamin C Results in Altered Proteostasis, Cardiomyopathy, and Arrhythmias

BACKGROUND

Many cardiomyopathy-associated FLNC pathogenic variants are heterozygous truncations, and FLNC pathogenic variants are associated with arrhythmias. Arrhythmia triggers in filaminopathy are incompletely understood.

METHODS AND RESULTS

We describe an individual with biallelic FLNC pathogenic variants, p.Arg650X and c.970-4A>G, with peripartum cardiomyopathy and ventricular arrhythmias. We also describe clinical findings in probands with FLNC variants including Val2715fs87X, Glu2458Serfs71X, Phe106Leu, and c.970-4A>G with hypertrophic and dilated cardiomyopathy, atrial fibrillation, and ventricular tachycardia. Induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) were generated. The FLNC truncation, Arg650X/c.970-4A>G, showed a marked reduction in filamin C protein consistent with biallelic loss of function mutations. To assess loss of filamin C, gene editing of a healthy control iPSC line was used to generate a homozygous FLNC disruption in the actin binding domain. Because filamin C has been linked to protein quality control, we assessed the necessity of filamin C in iPSC-CMs for response to the proteasome inhibitor bortezomib. After exposure to low-dose bortezomib, FLNC-null iPSC-CMs showed an increase in the chaperone proteins BAG3, HSP70 (heat shock protein 70), and HSPB8 (small heat shock protein B8) and in the autophagy marker LC3I/II. FLNC null iPSC-CMs had prolonged electric field potential, which was further prolonged in the presence of low-dose bortezomib. FLNC null engineered heart tissues had impaired function after low-dose bortezomib.

CONCLUSIONS

FLNC pathogenic variants associate with a predisposition to arrhythmias, which can be modeled in iPSC-CMs. Reduction of filamin C prolonged field potential, a surrogate for action potential, and with bortezomib-induced proteasome inhibition, reduced filamin C led to greater arrhythmia potential and impaired function.

Managing differential performance of polygenic risk scores across groups: Real-world experience of the eMERGE Network

The differential performance of polygenic risk scores (PRSs) by group is one of the major ethical barriers to their clinical use. It is also one of the main practical challenges for any implementation effort. The social repercussions of how people are grouped in PRS research must be considered in communications with research participants, including return of results. Here, we outline the decisions faced and choices made by a large multi-site clinical implementation study returning PRSs to diverse participants in handling this issue of differential performance. Our approach to managing the complexities associated with the differential performance of PRSs serves as a case study that can help future implementers of PRSs to plot an anticipatory course in response to this issue.

Risk factors affecting polygenic score performance across diverse cohorts

Apart from ancestry, personal or environmental covariates may contribute to differences in polygenic score (PGS) performance. We analyzed effects of covariate stratification and interaction on body mass index (BMI) PGS (PGS BMI ) across four cohorts of European (N=491,111) and African (N=21,612) ancestry. Stratifying on binary covariates and quintiles for continuous covariates, 18/62 covariates had significant and replicable R 2 differences among strata. Covariates with the largest differences included age, sex, blood lipids, physical activity, and alcohol consumption, with R 2 being nearly double between best and worst performing quintiles for certain covariates. 28 covariates had significant PGS BMI -covariate interaction effects, modifying PGS BMI effects by nearly 20% per standard deviation change. We observed overlap between covariates that had significant R 2 differences among strata and interaction effects - across all covariates, their main effects on BMI were correlated with their maximum R 2 differences and interaction effects (0.56 and 0.58, respectively), suggesting high-PGS BMI individuals have highest R 2 and increase in PGS effect. Using quantile regression, we show the effect of PGS BMI increases as BMI itself increases, and that these differences in effects are directly related to differences in R 2 when stratifying by different covariates. Given significant and replicable evidence for context-specific PGS BMI performance and effects, we investigated ways to increase model performance taking into account non-linear effects. Machine learning models (neural networks) increased relative model R 2 (mean 23%) across datasets. Finally, creating PGS BMI directly from GxAge GWAS effects increased relative R 2 by 7.8%. These results demonstrate that certain covariates, especially those most associated with BMI, significantly affect both PGS BMI performance and effects across diverse cohorts and ancestries, and we provide avenues to improve model performance that consider these effects.

Selection, optimization and validation of ten chronic disease polygenic risk scores for clinical implementation in diverse US populations

Polygenic risk scores (PRSs) have improved in predictive performance, but several challenges remain to be addressed before PRSs can be implemented in the clinic, including reduced predictive performance of PRSs in diverse populations, and the interpretation and communication of genetic results to both providers and patients. To address these challenges, the National Human Genome Research Institute-funded Electronic Medical Records and Genomics (eMERGE) Network has developed a framework and pipeline for return of a PRS-based genome-informed risk assessment to 25,000 diverse adults and children as part of a clinical study. From an initial list of 23 conditions, ten were selected for implementation based on PRS performance, medical actionability and potential clinical utility, including cardiometabolic diseases and cancer. Standardized metrics were considered in the selection process, with additional consideration given to strength of evidence in African and Hispanic populations. We then developed a pipeline for clinical PRS implementation (score transfer to a clinical laboratory, validation and verification of score performance), and used genetic ancestry to calibrate PRS mean and variance, utilizing genetically diverse data from 13,475 participants of the All of Us Research Program cohort to train and test model parameters. Finally, we created a framework for regulatory compliance and developed a PRS clinical report for return to providers and for inclusion in an additional genome-informed risk assessment. The initial experience from eMERGE can inform the approach needed to implement PRS-based testing in diverse clinical settings.

The impact of damaging epilepsy and cardiac genetic variant burden in sudden death in the young

BACKGROUND

Sudden unexpected death in children is a tragic event. Understanding the genetics of sudden death in the young (SDY) enables family counseling and cascade screening. The objective of this study was to characterize genetic variation in an SDY cohort using whole genome sequencing.

METHODS

The SDY Case Registry is a National Institutes of Health/Centers for Disease Control and Prevention surveillance effort to discern the prevalence, causes, and risk factors for SDY. The SDY Case Registry prospectively collected clinical data and DNA biospecimens from SDY cases < 20 years of age. SDY cases were collected from medical examiner and coroner offices spanning 13 US jurisdictions from 2015 to 2019. The cohort included 211 children (median age 0.33 year; range 0-20 years), determined to have died suddenly and unexpectedly and from whom DNA biospecimens for DNA extractions and next-of-kin consent were ascertained. A control cohort consisted of 211 randomly sampled, sex- and ancestry-matched individuals from the 1000 Genomes Project. Genetic variation was evaluated in epilepsy, cardiomyopathy, and arrhythmia genes in the SDY and control cohorts. American College of Medical Genetics/Genomics guidelines were used to classify variants as pathogenic or likely pathogenic. Additionally, pathogenic and likely pathogenic genetic variation was identified using a Bayesian-based artificial intelligence (AI) tool.

RESULTS

The SDY cohort was 43% European, 29% African, 3% Asian, 16% Hispanic, and 9% with mixed ancestries and 39% female. Six percent of the cohort was found to harbor a pathogenic or likely pathogenic genetic variant in an epilepsy, cardiomyopathy, or arrhythmia gene. The genomes of SDY cases, but not controls, were enriched for rare, potentially damaging variants in epilepsy, cardiomyopathy, and arrhythmia-related genes. A greater number of rare epilepsy genetic variants correlated with younger age at death.

CONCLUSIONS

While damaging cardiomyopathy and arrhythmia genes are recognized contributors to SDY, we also observed an enrichment in epilepsy-related genes in the SDY cohort and a correlation between rare epilepsy variation and younger age at death. These findings emphasize the importance of considering epilepsy genes when evaluating SDY.

Identification of Clinical Drivers of Left Atrial Enlargement Through Genomics of Left Atrial Size

BACKGROUND

Greater left atrial size is associated with a higher incidence of cardiovascular disease and mortality, but the full spectrum of diagnoses associated with left atrial enlargement in sex-stratified clinical populations is not well known. Our study sought to identify genetic risk mechanisms affecting left atrial diameter (LAD) in a clinical cohort.

METHODS

Using Vanderbilt deidentified electronic health record, we studied 6163 females and 5993 males of European ancestry who had at least 1 LAD measure and available genotyping. A sex-stratified polygenic score was constructed for LAD variation and tested for association against 1680 International Classification of Diseases code-based phenotypes. Two-sample univariable and multivariable Mendelian randomization approaches were used to assess etiologic relationships between candidate associations and LAD.

RESULTS

A phenome-wide association study identified 25 International Classification of Diseases code-based diagnoses in females and 11 in males associated with a polygenic score of LAD (false discovery rate q<0.01), 5 of which were further evaluated by Mendelian randomization (waist circumference [WC], atrial fibrillation, heart failure, systolic blood pressure, and coronary artery disease). Sex-stratified differences in the genetic associations between risk factors and a polygenic score for LAD were observed (WC for females; heart failure, systolic blood pressure, atrial fibrillation, and WC for males). By multivariable Mendelian randomization, higher WC remained significantly associated with larger LAD in females, whereas coronary artery disease, WC, and atrial fibrillation remained significantly associated with larger LAD in males.

CONCLUSIONS

In a clinical population, we identified, by genomic approaches, potential etiologic risk factors for larger LAD. Further studies are needed to confirm the extent to which these risk factors may be modified to prevent or reverse adverse left atrial remodeling and the extent to which sex modifies these risk factors.

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.

Selection, optimization, and validation of ten chronic disease polygenic risk scores for clinical implementation in diverse populations

Polygenic risk scores (PRS) have improved in predictive performance supporting their use in clinical practice. Reduced predictive performance of PRS in diverse populations can exacerbate existing health disparities. The NHGRI-funded eMERGE Network is returning a PRS-based genome-informed risk assessment to 25,000 diverse adults and children. We assessed PRS performance, medical actionability, and potential clinical utility for 23 conditions. Standardized metrics were considered in the selection process with additional consideration given to strength of evidence in African and Hispanic populations. Ten conditions were selected with a range of high-risk thresholds: atrial fibrillation, breast cancer, chronic kidney disease, coronary heart disease, hypercholesterolemia, prostate cancer, asthma, type 1 diabetes, obesity, and type 2 diabetes. We developed a pipeline for clinical PRS implementation, used genetic ancestry to calibrate PRS mean and variance, created a framework for regulatory compliance, and developed a PRS clinical report. eMERGE's experience informs the infrastructure needed to implement PRS-based implementation in diverse clinical settings.

Returning integrated genomic risk and clinical recommendations: The eMERGE study

PURPOSE

Assessing the risk of common, complex diseases requires consideration of clinical risk factors as well as monogenic and polygenic risks, which in turn may be reflected in family history. Returning risks to individuals and providers may influence preventive care or use of prophylactic therapies for those individuals at high genetic risk.

METHODS

To enable integrated genetic risk assessment, the eMERGE (electronic MEdical Records and GEnomics) network is enrolling 25,000 diverse individuals in a prospective cohort study across 10 sites. The network developed methods to return cross-ancestry polygenic risk scores, monogenic risks, family history, and clinical risk assessments via a genome-informed risk assessment (GIRA) report and will assess uptake of care recommendations after return of results.

RESULTS

GIRAs include summary care recommendations for 11 conditions, education pages, and clinical laboratory reports. The return of high-risk GIRA to individuals and providers includes guidelines for care and lifestyle recommendations. Assembling the GIRA required infrastructure and workflows for ingesting and presenting content from multiple sources. Recruitment began in February 2022.

CONCLUSION

Return of a novel report for communicating monogenic, polygenic, and family history-based risk factors will inform the benefits of integrated genetic risk assessment for routine health care.

The impact of damaging epilepsy and cardiac genetic variant burden in sudden death in the young

Background

Sudden unexpected death in children is a tragic event. Understanding the genetics of sudden death in the young (SDY) enables family counseling and cascade screening. The objective of this study was to characterize genetic variation in an SDY cohort using whole genome sequencing.

Methods

The SDY Case Registry is a National Institutes of Health/Centers for Disease Control surveillance effort to discern the prevalence, causes, and risk factors for SDY. The SDY Case Registry prospectively collected clinical data and DNA biospecimens from SDY cases <20 years of age. SDY cases were collected from medical examiner and coroner offices spanning 13 US jurisdictions from 2015-2019. The cohort included 211 children (mean age 1 year; range 0-20 years), determined to have died suddenly and unexpectedly and in whom DNA biospecimens and next-of-kin consent were ascertained. A control cohort consisted of 211 randomly sampled, sex-and ancestry-matched individuals from the 1000 Genomes Project. Genetic variation was evaluated in epilepsy, cardiomyopathy and arrhythmia genes in the SDY and control cohorts. American College of Medical Genetics/Genomics guidelines were used to classify variants as pathogenic or likely pathogenic. Additionally, genetic variation predicted to be damaging was identified using a Bayesian-based artificial intelligence (AI) tool.

Results

The SDY cohort was 42% European, 30% African, 17% Hispanic, and 11% with mixed ancestries, and 39% female. Six percent of the cohort was found to harbor a pathogenic or likely pathogenic genetic variant in an epilepsy, cardiomyopathy or arrhythmia gene. The genomes of SDY cases, but not controls, were enriched for rare, damaging variants in epilepsy, cardiomyopathy and arrhythmia-related genes. A greater number of rare epilepsy genetic variants correlated with younger age at death.

Conclusions

While damaging cardiomyopathy and arrhythmia genes are recognized contributors to SDY, we also observed an enrichment in epilepsy-related genes in the SDY cohort, and a correlation between rare epilepsy variation and younger age at death. These findings emphasize the importance of considering epilepsy genes when evaluating SDY.

Opposing effects of genetic variation in MTCH2 for obesity versus heart failure

Genetic variation in genes regulating metabolism may be advantageous in some settings but not others. The non-failing adult heart relies heavily on fatty acids as a fuel substrate and source of ATP. In contrast, the failing heart favors glucose as a fuel source. A bootstrap analysis for genes with deviant allele frequencies in cardiomyopathy cases versus controls identified the MTCH2 gene as having unusual variation. MTCH2 encodes an outer mitochondrial membrane protein, and prior genome-wide studies associated MTCH2 variants with body mass index, consistent with its role in metabolism. We identified the referent allele of rs1064608 (p.Pro290) as being overrepresented in cardiomyopathy cases compared to controls, and linkage disequilibrium analysis associated this variant with the MTCH2 cis eQTL rs10838738 and lower MTCH2 expression. To evaluate MTCH2, we knocked down Mtch in Drosophila heart tubes which produced a dilated and poorly functioning heart tube, reduced adiposity and shortened life span. Cardiac Mtch mutants generated more lactate at baseline, and they displayed impaired oxygen consumption in the presence of glucose but not palmitate. Treatment of cardiac Mtch mutants with dichloroacetate, a pyruvate dehydrogenase kinase inhibitor, reduced lactate and rescued lifespan. Deletion of MTCH2 in human cells similarly impaired oxygen consumption in the presence of glucose but not fatty acids. These data support a model in which MTCH2 reduction may be favorable when fatty acids are the major fuel source, favoring lean body mass. However, in settings like heart failure, where the heart shifts toward using more glucose, reduction of MTCH2 is maladaptive.

Quantifying factors that affect polygenic risk score performance across diverse ancestries and age groups for body mass index

Polygenic risk scores (PRS) have led to enthusiasm for precision medicine. However, it is well documented that PRS do not generalize across groups differing in ancestry or sample characteristics e.g., age. Quantifying performance of PRS across different groups of study participants, using genome-wide association study (GWAS) summary statistics from multiple ancestry groups and sample sizes, and using different linkage disequilibrium (LD) reference panels may clarify which factors are limiting PRS transferability. To evaluate these factors in the PRS generation process, we generated body mass index (BMI) PRS (PRSBMI) in the Electronic Medical Records and Genomics (eMERGE) network (N=75,661). Analyses were conducted in two ancestry groups (European and African) and three age ranges (adult, teenagers, and children). For PRSBMI calculations, we evaluated five LD reference panels and three sets of GWAS summary statistics of varying sample size and ancestry. PRSBMI performance increased for both African and European ancestry individuals using cross-ancestry GWAS summary statistics compared to European-only summary statistics (6.3% and 3.7% relative R2 increase, respectively, pAfrican=0.038, pEuropean=6.26x10-4). The effects of LD reference panels were more pronounced in African ancestry study datasets. PRSBMI performance degraded in children; R2 was less than half of teenagers or adults. The effect of GWAS summary statistics sample size was small when modeled with the other factors. Additionally, the potential of using a PRS generated for one trait to predict risk for comorbid diseases is not well understood especially in the context of cross-ancestry analyses - we explored clinical comorbidities from the electronic health record associated with PRSBMI and identified significant associations with type 2 diabetes and coronary atherosclerosis. In summary, this study quantifies the effects that ancestry, GWAS summary statistic sample size, and LD reference panel have on PRS performance, especially in cross-ancestry and age-specific analyses.

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.

Partial and complete loss of myosin binding protein H-like cause cardiac conduction defects

A premature truncation of MYBPHL in humans and a loss of Mybphl in mice is associated with dilated cardiomyopathy, atrial and ventricular arrhythmias, and atrial enlargement. MYBPHL encodes myosin binding protein H-like (MyBP-HL). Prior work in mice indirectly identified Mybphl expression in the atria and in small puncta throughout the ventricle. Because of its genetic association with human and mouse cardiac conduction system disease, we evaluated the anatomical localization of MyBP-HL and the consequences of loss of MyBP-HL on conduction system function. Immunofluorescence microscopy of normal adult mouse ventricles identified MyBP-HL-positive ventricular cardiomyocytes that co-localized with the ventricular conduction system marker contactin-2 near the atrioventricular node and in a subset of Purkinje fibers. Mybphl heterozygous ventricles had a marked reduction of MyBP-HL-positive cells compared to controls. Lightsheet microscopy of normal perinatal day 5 mouse hearts showed enrichment of MyBP-HL-positive cells within and immediately adjacent to the contactin-2-positive ventricular conduction system, but this association was not apparent in Mybphl heterozygous hearts. Surface telemetry of Mybphl-null mice revealed atrioventricular block and atrial bigeminy, while intracardiac pacing revealed a shorter atrial relative refractory period and atrial tachycardia. Calcium transient analysis of isolated Mybphl-null atrial cardiomyocytes demonstrated an increased heterogeneity of calcium release and faster rates of calcium release compared to wild type controls. Super-resolution microscopy of Mybphl heterozygous and homozygous null atrial cardiomyocytes showed ryanodine receptor disorganization compared to wild type controls. Abnormal calcium release, shorter atrial refractory period, and atrial dilation seen in Mybphl null, but not wild type control hearts, agree with the observed atrial arrhythmias, bigeminy, and atrial tachycardia, whereas the proximity of MyBP-HL-positive cells with the ventricular conduction system provides insight into how a predominantly atrial expressed gene contributes to ventricular arrhythmias and ventricular dysfunction.

Genetically Based Atrial Fibrillation: Current Considerations for Diagnosis and Management

Atrial fibrillation (AF) is the most common atrial arrhythmia and is subcategorized into numerous clinical phenotypes. Given its heterogeneity, investigations into the genetic mechanisms underlying AF have been pursued in recent decades, with predominant analyses focusing on early onset or lone AF. Linkage analyses, genome wide association studies (GWAS), and single gene analyses have led to the identification of rare and common genetic variants associated with AF risk. Significant overlap with genetic variants implicated in dilated cardiomyopathy syndromes, including truncating variants of the sarcomere protein titin, have been identified through these analyses, in addition to other genes associated with cardiac structure and function. Despite this, widespread utilization of genetic testing in AF remains hindered by the unclear impact of genetic risk identification on clinical outcomes and the high prevalence of variants of unknown significance (VUS). However, genetic testing is a reasonable option for patients with early onset AF and in those with significant family history of arrhythmia. While many knowledge gaps remain, emerging data support genotyping to inform selection of AF therapeutics. In this review we highlight the current understanding of the complex genetic basis of AF and explore the overlap of AF with inherited cardiomyopathy syndromes. We propose a set of criteria for clinical genetic testing in AF patients and outline future steps for the integration of genetics into AF care. This article is protected by copyright. All rights reserved.

South Asian-Specific MYBPC3 Δ25bp Deletion Carriers Display Hypercontraction and Impaired Diastolic Function Under Exercise Stress

Background: A 25-base pair (25bp) intronic deletion in the MYBPC3 gene enriched in South Asians (SAs) is a risk allele for late-onset left ventricular (LV) dysfunction, hypertrophy, and heart failure (HF) with several forms of cardiomyopathy. However, the effect of this variant on exercise parameters has not been evaluated. Methods: As a pilot study, 10 asymptomatic SA carriers of the MYBPC3 Δ25bp variant (52.9 ± 2.14 years) and 10 age- and gender-matched non-carriers (NCs) (50.1 ± 2.7 years) were evaluated at baseline and under exercise stress conditions using bicycle exercise echocardiography and continuous cardiac monitoring. Results: Baseline echocardiography parameters were not different between the two groups. However, in response to exercise stress, the carriers of Δ25bp had significantly higher LV ejection fraction (%) (CI: 4.57 ± 1.93; p < 0.0001), LV outflow tract peak velocity (m/s) (CI: 0.19 ± 0.07; p < 0.0001), and higher aortic valve (AV) peak velocity (m/s) (CI: 0.103 ± 0.08; p = 0.01) in comparison to NCs, and E/A ratio, a marker of diastolic compliance, was significantly lower in Δ25bp carriers (CI: 0.107 ± 0.102; p = 0.038). Interestingly, LV end-diastolic diameter (LVIDdia) was augmented in NCs in response to stress, while it did not increase in Δ25bp carriers (CI: 0.239 ± 0.125; p = 0.0002). Further, stress-induced right ventricular systolic excursion velocity s' (m/s), as a marker of right ventricle function, increased similarly in both groups, but tricuspid annular plane systolic excursion increased more in carriers (slope: 0.008; p = 0.0001), suggesting right ventricle functional differences between the two groups. Conclusions: These data support that MYBPC3 Δ25bp is associated with LV hypercontraction under stress conditions with evidence of diastolic impairment.

Genomic Autopsy of Sudden Deaths in Young Individuals

Importance

Postmortem genetic testing of young individuals with sudden death has previously identified pathogenic gene variants. However, prior studies primarily considered highly penetrant monogenic variants, often without detailed decedent and family clinical information.

Objective

To assess genotype and phenotype risk in a diverse cohort of young decedents with sudden death and their families.

Design, Setting, and Participants

Pathological and whole-genome sequence analysis was conducted in a cohort referred from a national network of medical examiners. Cases were accrued prospectively from May 2015 to March 2019 across 24 US states. Analysis began September 2016 and ended November 2020.

Exposures

Evaluation of autopsy and clinical data integrated with whole-genome sequence data and family member evaluation.

Results

A total of 103 decedents (mean [SD] age at death, 23.7 [11.9] years; age range, 1-44 years), their surviving family members, and 140 sex- and genetic ancestry-matched controls were analyzed. Among 103 decedents, autopsy and clinical data review categorized 36 decedents with postmortem diagnoses, 23 decedents with findings of uncertain significance, and 44 with sudden unexplained death. Pathogenic/likely pathogenic (P/LP) genetic variants in arrhythmia or cardiomyopathy genes were identified in 13 decedents (12.6%). A multivariable analysis including decedent phenotype, ancestry, and sex demonstrated that younger decedents had a higher burden of P/LP variants and select variants of uncertain significance (effect size, -1.64; P = .001). These select, curated variants of uncertain significance in cardiac genes were more common in decedents than controls (83 of 103 decedents [86%] vs 100 of 140 controls [71%]; P = .005), and decedents harbored more rare cardiac variants than controls (2.3 variants per individual vs 1.8 in controls; P = .006). Genetic testing of 31 parent-decedent trios and 14 parent-decedent dyads revealed 8 transmitted P/LP variants and 1 de novo P/LP variant. Incomplete penetrance was present in 6 of 8 parents who transmitted a P/LP variant.

Conclusions and Relevance

Whole-genome sequencing effectively identified P/LP variants in cases of sudden death in young individuals, implicating both arrhythmia and cardiomyopathy genes. Genomic analyses and familial phenotype association suggest potentially additive, oligogenic risk mechanisms for sudden death in this cohort.

Identification of Cardiac Fibrosis in Young Adults With a Homozygous Frameshift Variant in SERPINE1

Importance

Cardiac fibrosis is exceedingly rare in young adults. Identification of genetic variants that cause early-onset cardiomyopathy may inform novel biological pathways. Experimental models and a single case report have linked genetic deficiency of plasminogen activator inhibitor-1 (PAI-1), a downstream target of cardiac transforming growth factor β, with cardiac fibrosis.

Objective

To perform detailed cardiovascular phenotyping and genotyping in young adults from an Amish family with a frameshift variant (c.699_700dupTA) in SERPINE1, the gene that codes for PAI-1.

Design, Setting, and Participants

This observational study included participants from 3 related nuclear families from an Amish community in the primary analysis and participants from the extended family in the secondary analysis. Participants were recruited from May 2015 to December 2016, and analysis took place from June 2015 to June 2020.

Main Outcomes and Measures

(1) Multimodality cardiovascular imaging (transthoracic echocardiography and cardiac magnetic resonance imaging), (2) whole-exome sequencing, and (3) induced pluripotent stem cell-derived cardiomyocytes.

Results

Among 17 participants included in the primary analysis, the mean (interquartile range) age was 23.7 (20.9-29.9) years and 9 individuals (52.9%) were confirmed to be homozygous for the SERPINE1 c.699_700dupTA variant. Late gadolinium enhancement was present in 6 of 9 homozygous participants (67%) with absolute PAI-1 deficiency vs 0 of 8 in the control group (P = .001). Late gadolinium enhancement patterns tended to be dense and linear, usually subepicardial but also midmyocardial and transmural with noncoronary distributions. Targeted whole-exome sequencing analysis identified that homozygosity for c.699_700dupTA SERPINE1 was the only shared pathogenic variant or variant of uncertain significance after examination of cardiomyopathy genes among those with late gadolinium enhancement. Induced pluripotent stem cell-derived cardiomyocytes from participants homozygous for the SERPINE1 c.699_700dupTA variant exhibited susceptibility to cardiomyocyte injury in response to angiotensin II (increased transforming growth factor β1 secretion and release of lactate dehydrogenase) compared with control induced pluripotent stem cell-derived cardiomyocytes. In a secondary analysis based on echocardiography in 155 individuals across 3 generations in the extended family, no difference in global longitudinal strain was observed in carriers for the SERPINE1 c.699_700dupTA variant compared with wild-type participants, supporting an autosomal recessive inheritance pattern.

Conclusions and Relevance

In this study, a highly penetrant, autosomal recessive, cardiac fibrosis phenotype among young adults with homozygous frameshift variant for SERPINE1 was identified, suggesting an optimal range of PAI-1 levels are needed for cardiac homeostasis.

Mitochondrial cardiomyopathy and ventricular arrhythmias associated with biallelic variants in C1QBP

Patients with biallelic mutations in the nuclear-encoded mitochondrial gene C1QBP/p32 have been described with syndromic features and autosomal recessive cardiomyopathy. We describe the clinical course in two siblings who developed cardiomyopathy and ventricular fibrillation in infancy. We provide genomic analysis and clinical-pathologic correlation. Both siblings had profound cardiac failure with ventricular arrhythmia. One child died suddenly. The second sibling survived resuscitation but required extracorporeal cardiopulmonary support and died shortly afterward. On cardiac autopsy, the left ventricle was hypertrophied in both children. Histological examination revealed prominent cardiomyocyte cytoplasmic clearing, and electron microscopy confirmed abnormal mitochondrial structure within cardiomyocytes. DNA sequencing revealed compound heterozygous variants in C1QBP (p.Thr40Asnfs*45 and p.Phe204Leu) in both children. Family segregation analysis demonstrated each variant was inherited from an unaffected, heterozygous parent. Inherited loss of C1QBP/p32 is associated with recessive cardiomyopathy, ventricular fibrillation, and sudden death in early life. Ultrastructural mitochondrial evaluation in the second child was similar to findings in engineered C1qbp-deficient mice. Rapid trio analysis can define rare biallelic variants in genes that may be implicated in sudden death and facilitate medical management and family planning. (184/200).

Genomic Context Differs Between Human Dilated Cardiomyopathy and Hypertrophic Cardiomyopathy

Background Inherited cardiomyopathies display variable penetrance and expression, and a component of phenotypic variation is genetically determined. To evaluate the genetic contribution to this variable expression, we compared protein coding variation in the genomes of those with hypertrophic cardiomyopathy (HCM) and dilated cardiomyopathy (DCM). Methods and Results Nonsynonymous single-nucleotide variants (nsSNVs) were ascertained using whole genome sequencing from familial cases of HCM (n=56) or DCM (n=70) and correlated with echocardiographic information. Focusing on nsSNVs in 102 genes linked to inherited cardiomyopathies, we correlated the number of nsSNVs per person with left ventricular measurements. Principal component analysis and generalized linear models were applied to identify the probability of cardiomyopathy type as it related to the number of nsSNVs in cardiomyopathy genes. The probability of having DCM significantly increased as the number of cardiomyopathy gene nsSNVs per person increased. The increase in nsSNVs in cardiomyopathy genes significantly associated with reduced left ventricular ejection fraction and increased left ventricular diameter for individuals carrying a DCM diagnosis, but not for those with HCM. Resampling was used to identify genes with aberrant cumulative allele frequencies, identifying potential modifier genes for cardiomyopathy. Conclusions Participants with DCM had more nsSNVs per person in cardiomyopathy genes than participants with HCM. The nsSNV burden in cardiomyopathy genes did not correlate with the probability or manifestation of left ventricular measures in HCM. These findings support the concept that increased variation in cardiomyopathy genes creates a genetic background that predisposes to DCM and increased disease severity.

Genetic Variation in Enhancers Modifies Cardiomyopathy Gene Expression and Progression

BACKGROUND

Inherited cardiomyopathy associates with a range of phenotypes, mediated by genetic and nongenetic factors. Noninherited cardiomyopathy also displays varying progression and outcomes. Expression of cardiomyopathy genes is under the regulatory control of promoters and enhancers, and human genetic variation in promoters and enhancers may contribute to this variability.

METHODS

We superimposed epigenomic profiling from hearts and cardiomyocytes, including promoter-capture chromatin conformation information, to identify enhancers for 2 cardiomyopathy genes, MYH7 and LMNA. Enhancer function was validated in human cardiomyocytes derived from induced pluripotent stem cells. We also conducted a genome-wide search to ascertain genomic variation in enhancers positioned to alter cardiac expression and correlated one of these variants to cardiomyopathy progression using biobank data.

RESULTS

Multiple enhancers were identified and validated for LMNA and MYH7, including a key enhancer that regulates the switch from MYH6 expression to MYH7 expression. Deletion of this enhancer resulted in a dose-dependent increase in MYH6 and faster contractile rate in engineered heart tissues. We searched for genomic variation in enhancer sequences across the genome, with a focus on nucleotide changes that create or interrupt transcription factor binding sites. The sequence variant, rs875908, disrupts a T-Box Transcription Factor 5 binding motif and maps to an enhancer region 2 kilobases from the transcriptional start site of MYH7. Gene editing to remove the enhancer that harbors this variant markedly reduced MYH7 expression in human cardiomyocytes. Using biobank-derived data, rs875908 associated with longitudinal echocardiographic features of cardiomyopathy.

CONCLUSIONS

Enhancers regulate cardiomyopathy gene expression, and genomic variation within these enhancer regions associates with cardiomyopathic progression over time. This integrated approach identified noncoding modifiers of cardiomyopathy and is applicable to other cardiac genes.

Altered Enhancer and Promoter Usage Leads to Differential Gene Expression in the Normal and Failed Human Heart

BACKGROUND

The failing heart is characterized by changes in gene expression. However, the regulatory regions of the genome that drive these gene expression changes have not been well defined in human hearts.

METHODS

To define genome-wide enhancer and promoter use in heart failure, cap analysis of gene expression sequencing was applied to 3 healthy and 4 failed human hearts to identify promoter and enhancer regions used in left ventricles. Healthy hearts were derived from donors unused for transplantation and failed hearts were obtained as discarded tissue after transplantation.

RESULTS

Cap analysis of gene expression sequencing identified a combined potential for ≈23 000 promoters and ≈5000 enhancers active in human left ventricles. Of these, 17 000 promoters and 1800 enhancers had additional support for their regulatory function. Comparing promoter usage between healthy and failed hearts highlighted promoter shifts which altered aminoterminal protein sequences. Enhancer usage between healthy and failed hearts identified a majority of differentially used heart failure enhancers were intronic and primarily localized within the first intron, revealing this position as a common feature associated with tissue-specific gene expression changes in the heart.

CONCLUSIONS

This data set defines the dynamic genomic regulatory landscape underlying heart failure and serves as an important resource for understanding genetic contributions to cardiac dysfunction. Additionally, regulatory changes contributing to heart failure are attractive therapeutic targets for controlling ventricular remodeling and clinical progression.

Prevalence of Abnormal Heart Weight After Sudden Death in People Younger than 40 Years of Age

Background After sudden cardiac death in people aged <40 years, heart weight is a surrogate for cardiomegaly and a marker for cardiomyopathy. However, thresholds for cardiomegaly based on heart weight have not been validated in a cohort of cases of sudden cardiac death in young people. Methods and Results We surveyed medical examiner offices to determine which tools were used to assess heart weight norms. The survey determined that there was no gold standard for cardiomegaly (52 centers reported 22 different methods). We used a collection of heart weight data from sudden deaths in the Northwestern Sudden Death Collaboration (NSDC) to test the 22 methods. We found that the methods reported in our survey had little consistency: they classified between 18% and 81% of NSDC hearts with cardiomegaly. Therefore, we obtained biometric and postmortem data from a reference population of 3398 decedents aged <40 years. The reference population was ethnically diverse and had no known cardiac pathology on autopsy or histology. We derived and validated a multivariable regression model to predict normal heart weights and a threshold for cardiomegaly (upper 95% CI limit) in the young reference population (the Chicago model). Using the new model, the prevalence of cardiomegaly in hearts from the NSDC was 19%. Conclusions Medical examiner offices use a variety of tools to classify cardiomegaly. These approaches produce inconsistent results, and many overinterpret cardiomegaly. We recommend the model proposed to classify postmortem cardiomegaly in cases of sudden cardiac death in young people.

Abstract 439: Trajectory Analysis of Left Ventricular Dimensions From Biobank Data Uncovers Novel Genetic Associations

Background: Approximately 6 million adults in the United States have heart failure. The progression of heart failure is variable arising from differences in sex, age, genetic background including ancestry, and medication response. Many population-based genetic studies of heart failure have been cross-sectional in nature, failing to gain additional power from longitudinal analyses. As heart failure is known to change over time, using longitudinal data trajectories as a quantitative trait will increase power in genome wide association studies (GWAS).

Methods: We used the electronic health record in a racially and ethnically diverse medical biobank from a single, metropolitan US center. We used whole genome data from 896 unrelated participants analyzed, including 494 who had at least 1 electrocardiogram and 324 who had more than 1 echocardiogram (average of 3 observations per person). A mixture model based semiparametric latent growth curve model was used to cluster outcome measures used for genome-wide analyses.

Results: GWAS on the trajectory probability of QTc interval identified significant associations with variants in regulatory regions proximal to the WLS gene, which encodes Wntless, a Wnt ligand secretion mediator. WLS was previously associated with QTc and myocardial infarction, thus confirming the power of the method. GWAS on the trajectory probability of left ventricular diameter (LVIDd) identified significant associations with variants in regulatory regions near MYO10 , which encodes unconventional Myosin-10. MYO10 was previously associated with obesity and metabolic syndrome.

Conclusions: This is the first study to show an association with variants in or near MYO10 and left ventricular dimension changes over time. Further, we found that using trajectory probabilities can provide increased power to find novel associations with longitudinal data. This reduces the need for larger cohorts, and increases yield from smaller, well-phenotyped cohorts, such as those found in biobanks. This approach should be useful in the study of rare diseases and underrepresented populations.

Pathogenic and Uncertain Genetic Variants Have Clinical Cardiac Correlates in Diverse Biobank Participants

Background

Genome sequencing coupled with electronic heath record data can uncover medically important genetic variation. Interpretation of rare genetic variation and its role in mediating cardiovascular phenotypes is confounded by variants of uncertain significance.

Methods and Results

We analyzed the whole genome sequence of 900 racially and ethnically diverse biobank participants selected from a single US center. Participants were equally divided among European, African, Hispanic, and mixed races/ethnicities. We evaluated the American College of Medical Genetics and Genomics medically actionable list of 59 genes, focusing on the cardiac genes. Variation was interpreted using the most recent reports in ClinVar, a database of medically relevant human variation. We identified 19 individuals with pathogenic or likely pathogenic variants in cardiac actionable genes (2%) and found evidence of related clinical correlates in the electronic health record. Participants of African ancestry, compared with those of European ancestry, had more variants of uncertain significance in the medically actionable genes including the 30 cardiac actionable genes, even when normalized to total variant count per person. Longitudinal measures of left ventricle size from ≈400 biobank participants (1723 patient‐years) were correlated with genetic findings. The presence of ≥1 uncertain variant in the actionable cardiac genes and a cardiomyopathy diagnosis correlated with increased left ventricular internal diameter in diastole and in systole. In particular, MYBPC3 was identified as a gene with excess variants of uncertain significance.

Conclusions

These data indicate that a subset of uncertain genetic variants may confer risk and should not be considered benign.

Association of Cardiomyopathy With MYBPC3 D389V and MYBPC3Δ25bpIntronic Deletion in South Asian Descendants

Importance

The genetic variant MYBPC3Δ25bp occurs in 4% of South Asian descendants, with an estimated 100 million carriers worldwide. MYBPC3 Δ25bp has been linked to cardiomyopathy and heart failure. However, the high prevalence of MYBPC3Δ25bp suggests that other stressors act in concert with MYBPC3Δ25bp.

Objective

To determine whether there are additional genetic factors that contribute to the cardiomyopathic expression of MYBPC3Δ25bp.

Design, Setting, andParticipants

South Asian individuals living in the United States were screened for MYBPC3Δ25bp, and a subgroup was clinically evaluated using electrocardiograms and echocardiograms at Loyola University, Chicago, Illinois, between January 2015 and July 2016.

Main Outcomes and Measures

Next-generation sequencing of 174 cardiovascular disease genes was applied to identify additional modifying gene mutations and correlate genotype-phenotype parameters. Cardiomyocytes derived from human-induced pluripotent stem cells were established and examined to assess the role of MYBPC3Δ25bp.

Results

In this genotype-phenotype study, individuals of South Asian descent living in the United States from both sexes (36.23% female) with a mean population age of 48.92 years (range, 18-84 years) were recruited. Genetic screening of 2401 US South Asian individuals found an MYBPC3Δ25bpcarrier frequency of 6%. A higher frequency of missense TTN variation was found in MYBPC3Δ25bp carriers compared with noncarriers, identifying distinct genetic backgrounds within the MYBPC3Δ25bp carrier group. Strikingly, 9.6% of MYBPC3Δ25bp carriers also had a novel MYBPC3 variant, D389V. Family studies documented D389V was in tandem on the same allele as MYBPC3Δ25bp, and D389V was only seen in the presence of MYBPC3Δ25bp. In contrast to MYBPC3Δ25bp, MYBPC3Δ25bp/D389V was associated with hyperdynamic left ventricular performance (mean [SEM] left ventricular ejection fraction, 66.7 [0.7%]; left ventricular fractional shortening, 36.6 [0.6%]; P < .03) and stem cell-derived cardiomyocytes exhibited cellular hypertrophy with abnormal Ca2+ transients.

Conclusions and Relevance

MYBPC3Δ25bp/D389V is associated with hyperdynamic features, which are an early finding in hypertrophic cardiomyopathy and thought to reflect an unfavorable energetic state. These findings support that a subset of MYBPC3Δ25bp carriers, those with D389V, account for the increased risk attributed to MYBPC3Δ25bp.

Harmonizing Clinical Sequencing and Interpretation for the eMERGE III Network

The advancement of precision medicine requires new methods to coordinate and deliver genetic data from heterogeneous sources to physicians and patients. The eMERGE III Network enrolled >25,000 participants from biobank and prospective cohorts of predominantly healthy individuals for clinical genetic testing to determine clinically actionable findings. The network developed protocols linking together the 11 participant collection sites and 2 clinical genetic testing laboratories. DNA capture panels targeting 109 genes were used for testing of DNA and sample collection, data generation, interpretation, reporting, delivery, and storage were each harmonized. A compliant and secure network enabled ongoing review and reconciliation of clinical interpretations, while maintaining communication and data sharing between clinicians and investigators. A total of 202 individuals had positive diagnostic findings relevant to the indication for testing and 1,294 had additional/secondary findings of medical significance deemed to be returnable, establishing data return rates for other testing endeavors. This study accomplished integration of structured genomic results into multiple electronic health record (EHR) systems, setting the stage for clinical decision support to enable genomic medicine. Further, the established processes enable different sequencing sites to harmonize technical and interpretive aspects of sequencing tests, a critical achievement toward global standardization of genomic testing. The eMERGE protocols and tools are available for widespread dissemination.

Abstract 469: Genomic Context Predicts Dilated but Not Hypertrophic Cardiomyopathy

Rationale: Cardiomyopathies are a major cause of heart failure and have a strong heritable component.

Objective: Determine the role of both common and rare nonsynonymous genetic variation in hypertrophic (HCM) and dilated familial cardiomyopathy (DCM).

Methods and Results: Whole genome sequencing was used to determine common and rare nonsynonymous genetic variation in familial cases of HCM (n=56) or DCM (n=71). Variation was evaluated in 102 cardiomyopathy genes routinely assayed in clinical and commercial gene testing panels. We used cardiac gene expression data from GTEx (Genotype-Tissue Expression database) to define additional genes expressed in the heart. The number of nonsynonymous single nucleotide variants (nsSNVs), the majority of which were missense variants, was correlated with echocardiographic measurements. Principal component analysis (PCA) of left ventricular measures separated HCM and DCM. Regression of the first principal component using all nonsynonymous single nucleotide variants (nsSNVs) in cardiomyopathy genes showed the number of nsSNVs predicted DCM but not HCM. DCM probability in the cohort significantly increased as the number of cardiomyopathy gene nsSNVs increased (p<0.02). The increase in nsSNVs in cardiomyopathy genes significantly associated with reduced left ventricular ejection fraction and increased left ventricular diameter in DCM. Resampling methods identified genes with deviant cumulative allele frequencies, identifying potential modifier genes for cardiomyopathy.

Conclusions: DCM subjects carry a greater burden of nsSNVs in cardiomyopathy genes. This genomic burden translates to impaired systolic cardiac function in DCM. In contrast, nsSNV burden in cardiomyopathy genes did not correlate with the probability or manifestation of left ventricular measures in HCM. These findings support a complex inheritance for DCM where increased variation in cardiomyopathy genes creates a genetic background that predisposes to DCM and increased disease severity. The distinct genetic landscapes of HCM and DCM suggest that greater genetic variation in cardiac genes provokes unfavorable ventricular remodeling with reduced systolic function.

Prelamin A causes aberrant myonuclear arrangement and results in muscle fiber weakness

Physiological and premature aging are frequently associated with an accumulation of prelamin A, a precursor of lamin A, in the nuclear envelope of various cell types. Here, we aimed to underpin the hitherto unknown mechanisms by which prelamin A alters myonuclear organization and muscle fiber function. By experimentally studying membrane-permeabilized myofibers from various transgenic mouse lines, our results indicate that, in the presence of prelamin A, the abundance of nuclei and myosin content is markedly reduced within muscle fibers. This leads to a concept by which the remaining myonuclei are very distant from each other and are pushed to function beyond their maximum cytoplasmic capacity, ultimately inducing muscle fiber weakness.

Abstract 116: Mybphl is a Novel Myofilament Component Implicated in Arrhythmia and Dilated Cardiomyopathy

Dilated cardiomyopathy is under significant genetic influence. Using whole-genome sequencing, we found a premature stop variant (R255X) in the MYBPHL gene in a family with DCM and conduction system disease. A Mybphl null mouse model revealed systolic dysfunction with atrial and ventricular conduction system (VCS) abnormalities. We found MyBP-HL protein highly expressed in both human and mouse atria and in a subset of ventricular cardiomyocytes. We hypothesize that MyBP-HL regulates VCS function and loss of MyBP-HL causes ventricular arrhythmia and cardiac dysfunction. We now examined hearts and isolated VCMs from Mybphl WT, het, and null mice using immunofluorescence microscopy (IFM). Clusters of MyBP-HL+ ventricular cardiomyocytes were evident in WT hearts (~1/200,000), with only 10% as many MyBP-HL+ cardiomyocytes in het samples. In WT hearts, MyBP-HL+ cells localized in the right ventricular free wall, with significantly fewer MyBP-HL+ cells in the right ventricular free wall in het hearts. The VCS marker contactin-2 (Cntn2) was used to assess MyBP-HL co-localization in the VCS. MyBP-HL co-localized with Cntn2+ cardiomyocytes in the transition zone between the atria and ventricle near the AV node, as well as in a subset of Cntn2+ Purkinje fibers. MyBP-HL+ ventricular cardiomyocytes were also found outside the VCS, suggesting other roles for those cells in regulating ventricular function. To discern MyBP-HL myofilament localization, we used super-resolution structured illumination microscopy on atrial cardiomyocytes. MyBP-HL co-localized with cMyBP-C in the C-zone of the A-band but also occupied space in the inner portion of the A-band near the M-line. This localization was not altered in het mice; however, the ratio of MyBP-HL in the inner portion of the A-band was diminished in het mice compared to WT. This suggests that MyBP-HL may bind the light meromyosin coiled-coil tail that originates from the M-line. Following the discovery that loss of MyBP-HL is associated with DCM and arrhythmia, we identified that MyBP-HL associates with the VCS and reduced MyBP-HL levels reduce the amount of MyBP-HL positive VCMs. These data suggest that loss of MyBP-HL may impair the development and function of the VCS, resulting in arrhythmias.

Efficient exon skipping of SGCG mutations mediated by phosphorodiamidate morpholino oligomers

Exon skipping uses chemically modified antisense oligonucleotides to modulate RNA splicing. Therapeutically, exon skipping can bypass mutations and restore reading frame disruption by generating internally truncated, functional proteins to rescue the loss of native gene expression. Limb-girdle muscular dystrophy type 2C is caused by autosomal recessive mutations in the SGCG gene, which encodes the dystrophin-associated protein γ-sarcoglycan. The most common SGCG mutations disrupt the transcript reading frame abrogating γ-sarcoglycan protein expression. In order to treat most SGCG gene mutations, it is necessary to skip 4 exons in order to restore the SGCG transcript reading frame, creating an internally truncated protein referred to as Mini-Gamma. Using direct reprogramming of human cells with MyoD, myogenic cells were tested with 2 antisense oligonucleotide chemistries, 2'-O-methyl phosphorothioate oligonucleotides and vivo-phosphorodiamidate morpholino oligomers, to induce exon skipping. Treatment with vivo-phosphorodiamidate morpholino oligomers demonstrated efficient skipping of the targeted exons and corrected the mutant reading frame, resulting in the expression of a functional Mini-Gamma protein. Antisense-induced exon skipping of SGCG occurred in normal cells and those with multiple distinct SGCG mutations, including the most common 521ΔT mutation. These findings demonstrate a multiexon-skipping strategy applicable to the majority of limb-girdle muscular dystrophy 2C patients.

Novel nesprin-1 mutations associated with dilated cardiomyopathy cause nuclear envelope disruption and defects in myogenesis

Nesprins-1 and -2 are highly expressed in skeletal and cardiac muscle and together with SUN (Sad1p/UNC84)-domain containing proteins and lamin A/C form the LInker of Nucleoskeleton-and-Cytoskeleton (LINC) bridging complex at the nuclear envelope (NE). Mutations in nesprin-1/2 have previously been found in patients with autosomal dominant Emery–Dreifuss muscular dystrophy (EDMD) as well as dilated cardiomyopathy (DCM). In this study, three novel rare variants (R8272Q, S8381C and N8406K) in the C-terminus of the SYNE1 gene (nesprin-1) were identified in seven DCM patients by mutation screening. Expression of these mutants caused nuclear morphology defects and reduced lamin A/C and SUN2 staining at the NE. GST pull-down indicated that nesprin-1/lamin/SUN interactions were disrupted. Nesprin-1 mutations were also associated with augmented activation of the ERK pathway in vitro and in hearts in vivo. During C2C12 muscle cell differentiation, nesprin-1 levels are increased concomitantly with kinesin light chain (KLC-1/2) and immunoprecipitation and GST pull-down showed that these proteins interacted via a recently identified LEWD domain in the C-terminus of nesprin-1. Expression of nesprin-1 mutants in C2C12 cells caused defects in myoblast differentiation and fusion associated with dysregulation of myogenic transcription factors and disruption of the nesprin-1 and KLC-1/2 interaction at the outer nuclear membrane. Expression of nesprin-1α₂ WT and mutants in zebrafish embryos caused heart developmental defects that varied in severity. These findings support a role for nesprin-1 in myogenesis and muscle disease, and uncover a novel mechanism whereby disruption of the LINC complex may contribute to the pathogenesis of DCM.

Experimental Modeling Supports a Role for MyBP-HL as a Novel Myofilament Component in Arrhythmia and Dilated Cardiomyopathy

BACKGROUND

Cardiomyopathy and arrhythmias are under significant genetic influence. Here, we studied a family with dilated cardiomyopathy and associated conduction system disease in whom prior clinical cardiac gene panel testing was unrevealing.

METHODS

Whole-genome sequencing and induced pluripotent stem cells were used to examine a family with dilated cardiomyopathy and atrial and ventricular arrhythmias. We also characterized a mouse model with heterozygous and homozygous deletion of Mybphl.

RESULTS

Whole-genome sequencing identified a premature stop codon, R255X, in the MYBPHL gene encoding MyBP-HL (myosin-binding protein-H like), a novel member of the myosin-binding protein family. MYBPHL was found to have high atrial expression with low ventricular expression. We determined that MyBP-HL protein was myofilament associated in the atria, and truncated MyBP-HL protein failed to incorporate into the myofilament. Human cell modeling demonstrated reduced expression from the mutant MYBPHL allele. Echocardiography of Mybphl heterozygous and null mouse hearts exhibited a 36% reduction in fractional shortening and an increased diastolic ventricular chamber size. Atria weight normalized to total heart weight was significantly increased in Mybphl heterozygous and null mice. Using a reporter system, we detected robust expression of Mybphl in the atria, and in discrete puncta throughout the right ventricular wall and septum, as well. Telemetric electrocardiogram recordings in Mybphl mice revealed cardiac conduction system abnormalities with aberrant atrioventricular conduction and an increased rate of arrhythmia in heterozygous and null mice.

CONCLUSIONS

The findings of reduced ventricular function and conduction system defects in Mybphl mice support that MYBPHL truncations may increase risk for human arrhythmias and cardiomyopathy.

Abstract 74: Whole Genome Sequencing as a Diagnostic Tool for Cardiomyopathy

Background: Congestive heart failure (CHF) is an increasing medical problem that is disabling and costly, and CHF currently affects more than 5 million Americans. A leading cause of CHF is cardiomyopathy, a disorder with a high heritable component. Genetic studies of familial cardiomyopathy have identified more than 50 responsible genes. Knowledge of the genetic underpinnings of disease provides prognostic information and can guide clinical decision- making. Clinical genetic testing employs gene panels for diagnosis, and typical panels include 20-100 genes. We assessed the viability of whole genome sequencing (WGS) to identify potentially pathogenic cardiomyopathy-associated genetic variants in established and novel genes.

Methods: WGS was conducted on 91 families (99 individuals) with cardiomyopathy. In 46/91 (50%) families, there was a prior history of unrevealing panel-based clinical genetic testing. WGS was performed using Illumina next generation sequencing. Variants were called and annotated using MegaSeq and MegaSeq2. MegaSeq2 refines the number of candidate variants from ~4 million to <30 using variant effect prediction, conservation at variant site, frequency, and expression data.

Results: MegaSeq2 was applied to 99 cardiomyopathy genomes and in 45/91 families (49%), a potentially pathogenic or likely pathogenic variant was identified. Here we describe identification and segregation analysis in 4 families. Two families have unique RBM20 variants that segregate with cardiomyopathy. One family was found to carry a premature stop in SCN5A and exhibits variable phenotypes including sudden death, arrhythmia and cardiomyopathy. The fourth family has cardiomyopathy that segregates with a BAG3 early termination variant.

Conclusions: WGS was able to identify potentially pathogenic variants in nearly half (45/91). While some of the variants identified were in novel genes, the majority of pathogenic variants were in established cardiomyopathy genes despite previous panel-based testing. More comprehensive analysis, like WGS, allows for re-evaluation of variation as additional genetic information becomes available.

Abstract 391: The Pathophysiological Role of MYBPHL in Dilated Cardiomyopathy

Background: Cardiomyopathy is a leading cause of heart failure and is highly heritable. One common form of cardiomyopathy is dilated cardiomyopathy (DCM), which currently has over 70 identified genes that have been described as causative for the disease. Genetic testing for DCM employs gene panels and has a sensitivity of mutation detection less than 50%, indicating that additional genes contribute to DCM. Here, we employed whole genome sequencing (WGS) in a family with DCM and heart block who had previously undergone unrevealing genetic testing. We identified a premature stop codon in the MYBPHL gene, a gene that has not previously been linked to DCM as a likely cause of DCM in this family. Myosin binding protein H Like (MyBP-HL) is a muscle-expressed protein bearing structural similarity to myosin binding protein C (MyBP-C), which is commonly mutated gene in cardiomyopathies.

Objective: Determine the physiological and pathophysiological role of Mybphl .

Results: RNA-seq and qPCR from mouse hearts revealed that Mybphl is highly expressed in the right and left atria with lower expression in the ventricle and virtually no expression in skeletal muscle. As MyBP-HL shares a high homology with the myofilament proteins cardiac myosin binding protein-C and H, we investigated if MyBP-HL is also myofilament-associated. We determined that MyBP-HL protein is myofilament-associated in the atria although not clearly so in ventricle. To assess the requirement of MyBP-HL in cardiac function, we used a mouse model with an insertional disruption of the Mybphl gene. These mice have deficits in in vivo cardiac function, with reduced fractional shortening. In addition, ECG recordings from the Mybphl null mice show conduction system abnormalities affecting atrioventricular conduction.

Conclusions: WGS identified a premature stop codon in MYBPHL in human DCM. A mouse model with a disrupted Mybphl gene showed similar pathophysiological features as the humans with reduced ventricular function and cardiac conduction system abnormalities. MyBP-HL is an important protein for normal cardiac function.

Genetic Variation in Cardiomyopathy and Cardiovascular Disorders

With the wider deployment of massively-parallel, next-generation sequencing, it is now possible to survey human genome data for research and clinical purposes. The reduced cost of producing short-read sequencing has now shifted the burden to data analysis. Analysis of genome sequencing remains challenged by the complexity of the human genome, including redundancy and the repetitive nature of genome elements and the large amount of variation in individual genomes. Public databases of human genome sequences greatly facilitate interpretation of common and rare genetic variation, although linking database sequence information to detailed clinical information is limited by privacy and practical issues. Genetic variation is a rich source of knowledge for cardiovascular disease because many, if not all, cardiovascular disorders are highly heritable. The role of rare genetic variation in predicting risk and complications of cardiovascular diseases has been well established for hypertrophic and dilated cardiomyopathy, where the number of genes that are linked to these disorders is growing. Bolstered by family data, where genetic variants segregate with disease, rare variation can be linked to specific genetic variation that offers profound diagnostic information. Understanding genetic variation in cardiomyopathy is likely to help stratify forms of heart failure and guide therapy. Ultimately, genetic variation may be amenable to gene correction and gene editing strategies.

Disruption of the lamin A and matrin-3 interaction by myopathic LMNA mutations

The nuclear face of the nuclear membrane is enriched with the intermediate filament protein lamin A. Mutations in LMNA, the gene encoding lamin A, lead to a diverse set of inherited conditions including myopathies that affect both the heart and skeletal muscle. To gain insight about lamin A protein interactions, binding proteins associated with the tail of lamin A were characterized. Of 130 nuclear proteins found associated with the lamin A tail, 17 (13%) were previously described lamin A binding partners. One protein not previously linked to lamin A, matrin-3, was selected for further study, because like LMNA mutations, matrin-3 has also been implicated in inherited myopathy. Matrin-3 binds RNA and DNA and is a nucleoplasmic protein originally identified from the insoluble nuclear fraction, referred to as the nuclear matrix. Anti-matrin-3 antibodies were found to co-immunoprecipitate lamin A, and the lamin-A binding domain was mapped to the carboxy-terminal half of matrin-3. Three-dimensional mapping of the lamin A-matrin-3 interface showed that the LMNA truncating mutation Δ303, which lacks the matrin-3 binding domain, was associated with an increased distance between lamin A and matrin-3. LMNA mutant cells are known to have altered biophysical properties and the matrin-3-lamin A interface is positioned to contribute to these defects.

The genetic landscape of cardiomyopathy and its role in heart failure

Heart failure is highly influenced by heritability, and nearly 100 genes link to familial cardiomyopathy. Despite the marked genetic diversity that underlies these complex cardiovascular phenotypes, several key genes and pathways have emerged. Hypertrophic cardiomyopathy is characterized by increased contractility and a greater energetic cost of cardiac output. Dilated cardiomyopathy is often triggered by mutations that disrupt the giant protein titin. The energetic consequences of these mutations offer molecular targets and opportunities for new drug development and gene correction therapies.

Targeted analysis of whole genome sequence data to diagnose genetic cardiomyopathy

BACKGROUND

Cardiomyopathy is highly heritable but genetically diverse. At present, genetic testing for cardiomyopathy uses targeted sequencing to simultaneously assess the coding regions of >50 genes. New genes are routinely added to panels to improve the diagnostic yield. With the anticipated $1000 genome, it is expected that genetic testing will shift toward comprehensive genome sequencing accompanied by targeted gene analysis. Therefore, we assessed the reliability of whole genome sequencing and targeted analysis to identify cardiomyopathy variants in 11 subjects with cardiomyopathy.

METHODS AND RESULTS

Whole genome sequencing with an average of 37× coverage was combined with targeted analysis focused on 204 genes linked to cardiomyopathy. Genetic variants were scored using multiple prediction algorithms combined with frequency data from public databases. This pipeline yielded 1 to 14 potentially pathogenic variants per individual. Variants were further analyzed using clinical criteria and segregation analysis, where available. Three of 3 previously identified primary mutations were detected by this analysis. In 6 subjects for whom the primary mutation was previously unknown, we identified mutations that segregated with disease, had clinical correlates, and had additional pathological correlation to provide evidence for causality. For 2 subjects with previously known primary mutations, we identified additional variants that may act as modifiers of disease severity. In total, we identified the likely pathological mutation in 9 of 11 (82%) subjects.

CONCLUSIONS

These pilot data demonstrate that ≈30 to 40× coverage whole genome sequencing combined with targeted analysis is feasible and sensitive to identify rare variants in cardiomyopathy-associated genes.

Genetic profiling for risk reduction in human cardiovascular disease

Cardiovascular disease is a major health concern affecting over 80,000,000 people in the U.S. alone. Heart failure, cardiomyopathy, heart rhythm disorders, atherosclerosis and aneurysm formation have significant heritable contribution. Supported by familial aggregation and twin studies, these cardiovascular diseases are influenced by genetic variation. Family-based linkage studies and population-based genome-wide association studies (GWAS) have each identified genes and variants important for the pathogenesis of cardiovascular disease. The advent of next generation sequencing has ushered in a new era in the genetic diagnosis of cardiovascular disease, and this is especially evident when considering cardiomyopathy, a leading cause of heart failure. Cardiomyopathy is a genetically heterogeneous disorder characterized by morphologically abnormal heart with abnormal function. Genetic testing for cardiomyopathy employs gene panels, and these panels assess more than 50 genes simultaneously. Despite the large size of these panels, the sensitivity for detecting the primary genetic defect is still only approximately 50%. Recently, there has been a shift towards applying broader exome and/or genome sequencing to interrogate more of the genome to provide a genetic diagnosis for cardiomyopathy. Genetic mutations in cardiomyopathy offer the capacity to predict clinical outcome, including arrhythmia risk, and genetic diagnosis often provides an early window in which to institute therapy. This discussion is an overview as to how genomic data is shaping the current understanding and treatment of cardiovascular disease.

Supercomputing for the parallelization of whole genome analysis

MOTIVATION

The declining cost of generating DNA sequence is promoting an increase in whole genome sequencing, especially as applied to the human genome. Whole genome analysis requires the alignment and comparison of raw sequence data, and results in a computational bottleneck because of limited ability to analyze multiple genomes simultaneously.

RESULTS

We now adapted a Cray XE6 supercomputer to achieve the parallelization required for concurrent multiple genome analysis. This approach not only markedly speeds computational time but also results in increased usable sequence per genome. Relying on publically available software, the Cray XE6 has the capacity to align and call variants on 240 whole genomes in ∼50 h. Multisample variant calling is also accelerated.

AVAILABILITY AND IMPLEMENTATION

The MegaSeq workflow is designed to harness the size and memory of the Cray XE6, housed at Argonne National Laboratory, for whole genome analysis in a platform designed to better match current and emerging sequencing volume.