Allelic imbalance and haploinsufficiency in MYBPC3-linked hypertrophic cardiomyopathy

The Genetic Factors Influencing Cardiomyopathies and Heart Failure across the Allele Frequency Spectrum

Heart failure (HF) remains a major cause of mortality and morbidity worldwide. Understanding the genetic basis of HF allows for the development of disease-modifying therapies, more appropriate risk stratification, and personalised management of patients. The advent of next-generation sequencing has enabled genome-wide association studies; moving beyond rare variants identified in a Mendelian fashion and detecting common DNA variants associated with disease. We summarise the latest GWAS and rare variant data on mixed and refined HF aetiologies, and cardiomyopathies. We describe the recent understanding of the functional impact of titin variants and highlight FHOD3 as a novel cardiomyopathy-associated gene. We describe future directions of research in this field and how genetic data can be leveraged to improve the care of patients with HF.

Bringing into focus the central domains C3-C6 of myosin binding protein C

Myosin binding protein C (MyBPC) is a multi-domain protein with each region having a distinct functional role in muscle contraction. The central domains of MyBPC have often been overlooked due to their unclear roles. However, recent research shows promise in understanding their potential structural and regulatory functions. Understanding the central region of MyBPC is important because it may have specialized function that can be used as drug targets or for disease-specific therapies. In this review, we provide a brief overview of the evolution of our understanding of the central domains of MyBPC in regard to its domain structures, arrangement and dynamics, interaction partners, hypothesized functions, disease-causing mutations, and post-translational modifications. We highlight key research studies that have helped advance our understanding of the central region. Lastly, we discuss gaps in our current understanding and potential avenues to further research and discovery.

Hypertrophic cardiomyopathy in MYBPC3 carriers in aging

Hypertrophic cardiomyopathy (HCM) is characterized by abnormal thickening of the myocardium, leading to arrhythmias, heart failure, and elevated risk of sudden cardiac death, particularly among the young. This inherited disease is predominantly caused by mutations in sarcomeric genes, among which those in the cardiac myosin binding protein-C3 (MYBPC3) gene are major contributors. HCM associated with MYBPC3 mutations usually presents in the elderly and ranges from asymptomatic to symptomatic forms, affecting numerous cardiac functions and presenting significant health risks with a spectrum of clinical manifestations. Regulation of MYBPC3 expression involves various transcriptional and translational mechanisms, yet the destiny of mutant MYBPC3 mRNA and protein in late-onset HCM remains unclear. Pathogenesis related to MYBPC3 mutations includes nonsense-mediated decay, alternative splicing, and ubiquitin-proteasome system events, leading to allelic imbalance and haploinsufficiency. Aging further exacerbates the severity of HCM in carriers of MYBPC3 mutations. Advancements in high-throughput omics techniques have identified crucial molecular events and regulatory disruptions in cardiomyocytes expressing MYBPC3 variants. This review assesses the pathogenic mechanisms that promote late-onset HCM through the lens of transcriptional, post-transcriptional, and post-translational modulation of MYBPC3, underscoring its significance in HCM across carriers. The review also evaluates the influence of aging on these processes and MYBPC3 levels during HCM pathogenesis in the elderly. While pinpointing targets for novel medical interventions to conserve cardiac function remains challenging, the emergence of personalized omics offers promising avenues for future HCM treatments, particularly for late-onset cases.

Identification of novel genetic risk factors of dilated cardiomyopathy: from canine to human

BACKGROUND

Dilated cardiomyopathy (DCM) is a life-threatening heart disease and a common cause of heart failure due to systolic dysfunction and subsequent left or biventricular dilatation. A significant number of cases have a genetic etiology; however, as a complex disease, the exact genetic risk factors are largely unknown, and many patients remain without a molecular diagnosis.

METHODS

We performed GWAS followed by whole-genome, transcriptome, and immunohistochemical analyses in a spontaneously occurring canine model of DCM. Canine gene discovery was followed up in three human DCM cohorts.

RESULTS

Our results revealed two independent additive loci associated with the typical DCM phenotype comprising left ventricular systolic dysfunction and dilatation. We highlight two novel candidate genes, RNF207 and PRKAA2, known for their involvement in cardiac action potentials, energy homeostasis, and morphology. We further illustrate the distinct genetic etiologies underlying the typical DCM phenotype and ventricular premature contractions. Finally, we followed up on the canine discoveries in human DCM patients and discovered candidate variants in our two novel genes.

CONCLUSIONS

Collectively, our study yields insight into the molecular pathophysiology of DCM and provides a large animal model for preclinical studies.

The penetrance of rare variants in cardiomyopathy-associated genes: A cross-sectional approach to estimating penetrance for secondary findings

Understanding the penetrance of pathogenic variants identified as secondary findings (SFs) is of paramount importance with the growing availability of genetic testing. We estimated penetrance through large-scale analyses of individuals referred for diagnostic sequencing for hypertrophic cardiomyopathy (HCM; 10,400 affected individuals, 1,332 variants) and dilated cardiomyopathy (DCM; 2,564 affected individuals, 663 variants), using a cross-sectional approach comparing allele frequencies against reference populations (293,226 participants from UK Biobank and gnomAD). We generated updated prevalence estimates for HCM (1:543) and DCM (1:220). In aggregate, the penetrance by late adulthood of rare, pathogenic variants (23% for HCM, 35% for DCM) and likely pathogenic variants (7% for HCM, 10% for DCM) was substantial for dominant cardiomyopathy (CM). Penetrance was significantly higher for variant subgroups annotated as loss of function or ultra-rare and for males compared to females for variants in HCM-associated genes. We estimated variant-specific penetrance for 316 recurrent variants most likely to be identified as SFs (found in 51% of HCM- and 17% of DCM-affected individuals). 49 variants were observed at least ten times (14% of affected individuals) in HCM-associated genes. Median penetrance was 14.6% (±14.4% SD). We explore estimates of penetrance by age, sex, and ancestry and simulate the impact of including future cohorts. This dataset reports penetrance of individual variants at scale and will inform the management of individuals undergoing genetic screening for SFs. While most variants had low penetrance and the costs and harms of screening are unclear, some individuals with highly penetrant variants may benefit from SFs.

An Unbiased Screen Identified the Hsp70-BAG3 Complex as a Regulator of Myosin-Binding Protein C3

Variants in the gene myosin-binding protein C3 (MYBPC3) account for approximately 50% of familial hypertrophic cardiomyopathy (HCM), leading to reduced levels of myosin-binding protein C3 (MyBP-C), the protein product made by gene MYBPC3. Elucidation of the pathways that regulate MyBP-C protein homeostasis could uncover new therapeutic strategies. Toward this goal, we screened a library of 2,426 bioactive compounds and identified JG98, an allosteric modulator of heat shock protein 70 that inhibits interaction with Bcl-2-associated athanogene (BAG) domain co-chaperones. JG98 reduces MyBP-C protein levels. Furthermore, genetic reduction of BAG3 phenocopies treatment with JG-98 by reducing MYBP-C protein levels.. Thus, an unbiased compound screen identified the heat shock protein 70-BAG3 complex as a regulator of MyBP-C stability.

Hypertrophic Cardiomyopathy: Genetic Foundations, Outcomes, Interconnections, and Their Modifiers

Hypertrophic cardiomyopathy (HCM) is the most prevalent heritable cardiomyopathy. HCM is considered to be caused by mutations in cardiac sarcomeric protein genes. Recent research suggests that the genetic foundation of HCM is much more complex than originally postulated. The clinical presentations of HCM are very variable. Some mutation carriers remain asymptomatic, while others develop severe HCM, terminal heart failure, or sudden cardiac death. Heterogeneity regarding both genetic mutations and the clinical course of HCM hinders the establishment of universal genotype-phenotype correlations. However, some trends have been identified. The presence of a mutation in some genes encoding sarcomeric proteins is associated with earlier HCM onset, more severe left ventricular hypertrophy, and worse clinical outcomes. There is a diversity in the mechanisms implicated in the pathogenesis of HCM. They may be classified into groups, but they are interrelated. The lack of known supplementary elements that control the progression of HCM indicates that molecular mechanisms that exist between genotype and clinical presentations may be crucial. Secondary molecular changes in pathways implicated in HCM pathogenesis, post-translational protein modifications, and epigenetic factors affect HCM phenotypes. Cardiac loading conditions, exercise, hypertension, diet, alcohol consumption, microbial infection, obstructive sleep apnea, obesity, and environmental factors are non-molecular aspects that change the HCM phenotype. Many mechanisms are implicated in the course of HCM. They are mostly interconnected and contribute to some extent to final outcomes.

Long-Term Prevalence of Systolic Dysfunction in MYBPC3 Versus MYH7-Related Hypertrophic Cardiomyopathy

BACKGROUND

The 2 sarcomere genes most commonly associated with hypertrophic cardiomyopathy (HCM), MYBPC3 (myosin-binding protein C3) and MYH7 (β-myosin heavy chain), are indistinguishable at presentation, and genotype-phenotype correlations have been elusive. Based on molecular and pathophysiological differences, however, it is plausible to hypothesize a different behavior in myocardial performance, impacting lifetime changes in left ventricular (LV) function.

METHODS

We reviewed the initial and final echocardiograms of 402 consecutive HCM patients with pathogenic or likely pathogenic MYBPC3 (n=251) or MYH7 (n=151) mutations, followed over 9±8 years.

RESULTS

At presentation, MYBPC3 patients were less frequently obstructive (15% versus 26%; P=0.005) and had lower LV ejection fraction compared with MYH7 (66±8% versus 68±8%, respectively; P=0.03). Both HCM patients harboring MYBPC3 and MYH7 mutations exhibited a small but significant decline in LV systolic function during follow-up; however, new onset of severe LV systolic dysfunction (LV ejection fraction, <50%) was greater among MYBPC3 patients (15% versus 5% among MYH7; P=0.013). Prevalence of grade II/III diastolic dysfunction at final evaluation was comparable between MYBPC3 and MYH7 patients (P=0.509). In a Cox multivariable analysis, MYBPC3-positive status (hazard ratio, 2.53 [95% CI, 1.09-5.82]; P=0.029), age (hazard ratio, 1.03 [95% CI, 1.00-1.06]; P=0.027), and atrial fibrillation (hazard ratio, 2.39 [95% CI, 1.14-5.05]; P=0.020) were independent predictors of severe systolic dysfunction. No statistically significant differences occurred with regard to incidence of atrial fibrillation, heart failure, appropriate implanted cardioverter defibrillator shock, or cardiovascular death.

CONCLUSIONS

MYBPC3-related HCM showed increased long-term prevalence of systolic dysfunction compared with MYH7, in spite of similar outcome. Such observations suggest different pathophysiology of clinical progression in the 2 subsets and may prove relevant for understanding of genotype-phenotype correlations in HCM.

An Update on MYBPC3 Gene Mutation in Hypertrophic Cardiomyopathy

Hypertrophic cardiomyopathy (HCM) is the most prevalent genetically inherited cardiomyopathy that follows an autosomal dominant inheritance pattern. The majority of HCM cases can be attributed to mutation of the MYBPC3 gene, which encodes cMyBP-C, a crucial structural protein of the cardiac muscle. The manifestation of HCM's morphological, histological, and clinical symptoms is subject to the complex interplay of various determinants, including genetic mutation and environmental factors. Approximately half of MYBPC3 mutations give rise to truncated protein products, while the remaining mutations cause insertion/deletion, frameshift, or missense mutations of single amino acids. In addition, the onset of HCM may be attributed to disturbances in the protein and transcript quality control systems, namely, the ubiquitin-proteasome system and nonsense-mediated RNA dysfunctions. The aforementioned genetic modifications, which appear to be associated with unfavorable lifelong outcomes and are largely influenced by the type of mutation, exhibit a unique array of clinical manifestations ranging from asymptomatic to arrhythmic syncope and even sudden cardiac death. Although the current understanding of the MYBPC3 mutation does not comprehensively explain the varied phenotypic manifestations witnessed in patients with HCM, patients with pathogenic MYBPC3 mutations can exhibit an array of clinical manifestations ranging from asymptomatic to advanced heart failure and sudden cardiac death, leading to a higher rate of adverse clinical outcomes. This review focuses on MYBPC3 mutation and its characteristics as a prognostic determinant for disease onset and related clinical consequences in HCM.

cMyBP-C ablation in human engineered cardiac tissue causes progressive Ca2+-handling abnormalities

Truncation mutations in cardiac myosin binding protein C (cMyBP-C) are common causes of hypertrophic cardiomyopathy (HCM). Heterozygous carriers present with classical HCM, while homozygous carriers present with early onset HCM that rapidly progress to heart failure. We used CRISPR-Cas9 to introduce heterozygous (cMyBP-C+/-) and homozygous (cMyBP-C-/-) frame-shift mutations into MYBPC3 in human iPSCs. Cardiomyocytes derived from these isogenic lines were used to generate cardiac micropatterns and engineered cardiac tissue constructs (ECTs) that were characterized for contractile function, Ca2+-handling, and Ca2+-sensitivity. While heterozygous frame shifts did not alter cMyBP-C protein levels in 2-D cardiomyocytes, cMyBP-C+/- ECTs were haploinsufficient. cMyBP-C-/- cardiac micropatterns produced increased strain with normal Ca2+-handling. After 2 wk of culture in ECT, contractile function was similar between the three genotypes; however, Ca2+-release was slower in the setting of reduced or absent cMyBP-C. At 6 wk in ECT culture, the Ca2+-handling abnormalities became more pronounced in both cMyBP-C+/- and cMyBP-C-/- ECTs, and force production became severely depressed in cMyBP-C-/- ECTs. RNA-seq analysis revealed enrichment of differentially expressed hypertrophic, sarcomeric, Ca2+-handling, and metabolic genes in cMyBP-C+/- and cMyBP-C-/- ECTs. Our data suggest a progressive phenotype caused by cMyBP-C haploinsufficiency and ablation that initially is hypercontractile, but progresses to hypocontractility with impaired relaxation. The severity of the phenotype correlates with the amount of cMyBP-C present, with more severe earlier phenotypes observed in cMyBP-C-/- than cMyBP-C+/- ECTs. We propose that while the primary effect of cMyBP-C haploinsufficiency or ablation may relate to myosin crossbridge orientation, the observed contractile phenotype is Ca2+-mediated.

Is haploinsufficiency a sufficient mechanism for MYBPC3 truncating mutations?

Reduced expression of MYBPC3 causes early dysfunction in human cell culture models prior to reduced cMyBP-C levels.

Transcriptional bursts and heterogeneity among cardiomyocytes in hypertrophic cardiomyopathy

Transcriptional bursting is a common expression mode for most genes where independent transcription of alleles leads to different ratios of allelic mRNA from cell to cell. Here we investigated burst-like transcription and its consequences in cardiac tissue from Hypertrophic Cardiomyopathy (HCM) patients with heterozygous mutations in the sarcomeric proteins cardiac myosin binding protein C (cMyBP-C, MYBPC3) and cardiac troponin I (cTnI, TNNI3). Using fluorescence in situ hybridization (RNA-FISH) we found that both, MYBPC3 and TNNI3 are transcribed burst-like. Along with that, we show unequal allelic ratios of TNNI3-mRNA among single cardiomyocytes and unequally distributed wildtype cMyBP-C protein across tissue sections from heterozygous HCM-patients. The mutations led to opposing functional alterations, namely increasing (cMyBP-C_c.927−2A>G) or decreasing (cTnI_R145W) calcium sensitivity. Regardless, all patients revealed highly variable calcium-dependent force generation between individual cardiomyocytes, indicating contractile imbalance, which appears widespread in HCM-patients. Altogether, we provide strong evidence that burst-like transcription of sarcomeric genes can lead to an allelic mosaic among neighboring cardiomyocytes at mRNA and protein level. In HCM-patients, this presumably induces the observed contractile imbalance among individual cardiomyocytes and promotes HCM-development.

Population Prevalence of Premature Truncating Variants in Plakophilin-2 and Association With Arrhythmogenic Right Ventricular Cardiomyopathy: A UK Biobank Analysis

BACKGROUND

Truncating variants in the desmosomal gene PKP2 (PKP2tv) cause arrhythmogenic right ventricular cardiomyopathy (ARVC) yet display varied penetrance and expressivity.

METHODS

We identified individuals with PKP2tv from the UK Biobank (UKB) and determined the prevalence of an ARVC phenotype and other cardiovascular traits based on clinical and procedural data. The PKP2tv minor allelic frequency in the UKB was compared with a second cohort of probands with a clinical diagnosis of ARVC (ARVC cohort), with a figure of 1:5000 assumed for disease prevalence. In silico predictors of variant pathogenicity (combined annotation-dependent depletion and Splice AI [Illumina, Inc.]) were assessed.

RESULTS

PKP2tv were identified in 193/200 643 (0.10%) UKB participants, with 47 unique PKP2tv. Features consistent with ARVC were present in 3 (1.6%), leaving 190 with PKP2tv without manifest disease (UKB cohort; minor allelic frequency 4.73×10-4). The ARVC cohort included 487 ARVC probands with 144 distinct PKP2tv, with 25 PKP2tv common to both cohorts. The odds ratio for ARVC for the 25 common PKP2tv was 0.047 (95% CI, 0.001-0.268; P=2.43×10-6), and only favored ARVC (odds ratio >1) for a single variant, p.Arg79*. In silico variant analysis did not differentiate PKP2tv between the 2 cohorts. Atrial fibrillation was over-represented in the UKB cohort in those with PKP2tv (7.9% versus 4.3%; odds ratio, 2.11; P=0.005).

CONCLUSIONS

PKP2tv are prevalent in the population and associated with ARVC in only a small minority, necessitating a more detailed understanding of how PKP2tv cause ARVC in combination with associated genetic and environmental risk factors.

Targeting the sarcomere in inherited cardiomyopathies

Variants in >12 genes encoding sarcomeric proteins can cause various cardiomyopathies. The two most common are hypertrophic cardiomyopathy (HCM) and dilated cardiomyopathy (DCM). Current therapeutics do not target the root causes of these diseases, but attempt to prevent disease progression and/or to manage symptoms. Accordingly, novel approaches are being developed to treat the cardiac muscle dysfunction directly. Challenges to developing therapeutics for these diseases include the diverse mechanisms of pathogenesis, some of which are still being debated and defined. Four small molecules that modulate the myosin motor protein in the cardiac sarcomere have shown great promise in the settings of HCM and DCM, regardless of the underlying genetic pathogenesis, and similar approaches are being developed to target other components of the sarcomere. In the setting of HCM, mavacamten and aficamten bind to the myosin motor and decrease the ATPase activity of myosin. In the setting of DCM, omecamtiv mecarbil and danicamtiv increase myosin activity in cardiac muscle (but omecamtiv mecarbil decreases myosin activity in vitro). In this Review, we discuss the therapeutic strategies to alter sarcomere contractile activity and summarize the data indicating that targeting one protein in the sarcomere can be effective in treating patients with genetic variants in other sarcomeric proteins, as well as in patients with non-sarcomere-based disease.

The mechanics of the heart: zooming in on hypertrophic cardiomyopathy and cMyBP-C

Hypertrophic cardiomyopathy (HCM), a disease characterized by cardiac muscle hypertrophy and hypercontractility, is the most frequently inherited disorder of the heart. HCM is mainly caused by variants in genes encoding proteins of the sarcomere, the basic contractile unit of cardiomyocytes. The most frequently mutated among them is MYBPC3, which encodes cardiac myosin-binding protein C (cMyBP-C), a key regulator of sarcomere contraction. In this review, we summarize clinical and genetic aspects of HCM and provide updated information on the function of the healthy and HCM sarcomere, as well as on emerging therapeutic options targeting sarcomere mechanical activity. Building on what is known about cMyBP-C activity, we examine different pathogenicity drivers by which MYBPC3 variants can cause disease, focussing on protein haploinsufficiency as a common pathomechanism also in nontruncating variants. Finally, we discuss recent evidence correlating altered cMyBP-C mechanical properties with HCM development.

Translation of New and Emerging Therapies for Genetic Cardiomyopathies

The primary etiology of a diverse range of cardiomyopathies is now understood to be genetic, creating a new paradigm for targeting treatments on the basis of the underlying molecular cause. This review provides a genetic and etiologic context for the traditional clinical classifications of cardiomyopathy, including molecular subtypes that may exhibit differential responses to existing or emerging treatments. The authors describe several emerging cardiomyopathy treatments, including gene therapy, direct targeting of myofilament function, protein quality control, metabolism, and others. The authors discuss advantages and disadvantages of these approaches and indicate areas of high potential for short- and longer term efficacy.

Incomplete Penetrance and Variable Expressivity: From Clinical Studies to Population Cohorts

The same genetic variant found in different individuals can cause a range of diverse phenotypes, from no discernible clinical phenotype to severe disease, even among related individuals. Such variants can be said to display incomplete penetrance, a binary phenomenon where the genotype either causes the expected clinical phenotype or it does not, or they can be said to display variable expressivity, in which the same genotype can cause a wide range of clinical symptoms across a spectrum. Both incomplete penetrance and variable expressivity are thought to be caused by a range of factors, including common variants, variants in regulatory regions, epigenetics, environmental factors, and lifestyle. Many thousands of genetic variants have been identified as the cause of monogenic disorders, mostly determined through small clinical studies, and thus, the penetrance and expressivity of these variants may be overestimated when compared to their effect on the general population. With the wealth of population cohort data currently available, the penetrance and expressivity of such genetic variants can be investigated across a much wider contingent, potentially helping to reclassify variants that were previously thought to be completely penetrant. Research into the penetrance and expressivity of such genetic variants is important for clinical classification, both for determining causative mechanisms of disease in the affected population and for providing accurate risk information through genetic counseling. A genotype-based definition of the causes of rare diseases incorporating information from population cohorts and clinical studies is critical for our understanding of incomplete penetrance and variable expressivity. This review examines our current knowledge of the penetrance and expressivity of genetic variants in rare disease and across populations, as well as looking into the potential causes of the variation seen, including genetic modifiers, mosaicism, and polygenic factors, among others. We also considered the challenges that come with investigating penetrance and expressivity.

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.

Assessment of the Contribution of a Thermodynamic and Mechanical Destabilization of Myosin-Binding Protein C Domain C2 to the Pathomechanism of Hypertrophic Cardiomyopathy-Causing Double Mutation MYBPC3Δ25bp/D389V

Mutations in the gene encoding cardiac myosin-binding protein-C (MyBPC), a thick filament assembly protein that stabilizes sarcomeric structure and regulates cardiac function, are a common cause for the development of hypertrophic cardiomyopathy. About 10% of carriers of the Δ25bp variant of MYBPC3, which is common in individuals from South Asia, are also carriers of the D389V variant on the same allele. Compared with noncarriers and those with MYBPC3^Δ25bp alone, indicators for the development of hypertrophic cardiomyopathy occur with increased frequency in MYBPC3^Δ25bp/D389V carriers. Residue D389 lies in the IgI-like C2 domain that is part of the N-terminal region of MyBPC. To probe the effects of mutation D389V on structure, thermostability, and protein–protein interactions, we produced and characterized wild-type and mutant constructs corresponding to the isolated 10 kDa C2 domain and a 52 kDa N-terminal fragment that includes subdomains C0 to C2. Our results show marked reductions in the melting temperatures of D389V mutant constructs. Interactions of construct C0–C2 D389V with the cardiac isoforms of myosin-2 and actin remain unchanged. Molecular dynamics simulations reveal changes in the stiffness and conformer dynamics of domain C2 caused by mutation D389V. Our results suggest a pathomechanism for the development of HCM based on the toxic buildup of misfolded protein in young MYBPC3^Δ25bp/D389V carriers that is supplanted and enhanced by C-zone haploinsufficiency at older ages.

Cardiomyocyte Dysfunction in Inherited Cardiomyopathies

Inherited cardiomyopathies form a heterogenous group of disorders that affect the structure and function of the heart. Defects in the genes encoding sarcomeric proteins are associated with various perturbations that induce contractile dysfunction and promote disease development. In this review we aimed to outline the functional consequences of the major inherited cardiomyopathies in terms of myocardial contraction and kinetics, and to highlight the structural and functional alterations in some sarcomeric variants that have been demonstrated to be involved in the pathogenesis of the inherited cardiomyopathies. A particular focus was made on mutation-induced alterations in cardiomyocyte mechanics. Since no disease-specific treatments for familial cardiomyopathies exist, several novel agents have been developed to modulate sarcomere contractility. Understanding the molecular basis of the disease opens new avenues for the development of new therapies. Furthermore, the earlier the awareness of the genetic defect, the better the clinical prognostication would be for patients and the better the prevention of development of the disease.

A genome-first approach to rare variants in hypertrophic cardiomyopathy genes MYBPC3 and MYH7 in a medical biobank

'Genome-first' approaches to analyzing rare variants can reveal new insights into human biology and disease. Because pathogenic variants are often rare, new discovery requires aggregating rare coding variants into 'gene burdens' for sufficient power. However, a major challenge is deciding which variants to include in gene burden tests. Pathogenic variants in MYBPC3 and MYH7 are well-known causes of hypertrophic cardiomyopathy (HCM), and focusing on these 'positive control' genes in a genome-first approach could help inform variant selection methods and gene burdening strategies for other genes and diseases. Integrating exome sequences with electronic health records among 41 759 participants in the Penn Medicine BioBank, we evaluated the performance of aggregating predicted loss-of-function (pLOF) and/or predicted deleterious missense (pDM) variants in MYBPC3 and MYH7 for gene burden phenome-wide association studies (PheWAS). The approach to grouping rare variants for these two genes produced very different results: pLOFs but not pDM variants in MYBPC3 were strongly associated with HCM, whereas the opposite was true for MYH7. Detailed review of clinical charts revealed that only 38.5% of patients with HCM diagnoses carrying an HCM-associated variant in MYBPC3 or MYH7 had a clinical genetic test result. Additionally, 26.7% of MYBPC3 pLOF carriers without HCM diagnoses had clear evidence of left atrial enlargement and/or septal/LV hypertrophy on echocardiography. Our study shows the importance of evaluating both pLOF and pDM variants for gene burden testing in future studies to uncover novel gene-disease relationships and identify new pathogenic loss-of-function variants across the human genome through genome-first analyses of healthcare-based populations.

Protein haploinsufficiency drivers identify MYBPC3 variants that cause hypertrophic cardiomyopathy

Hypertrophic cardiomyopathy (HCM) is the most common inherited cardiac disease. Variants in MYBPC3, the gene encoding cardiac myosin-binding protein C (cMyBP-C), are the leading cause of HCM. However, the pathogenicity status of hundreds of MYBPC3 variants found in patients remains unknown, as a consequence of our incomplete understanding of the pathomechanisms triggered by HCM-causing variants. Here, we examined 44 nontruncating MYBPC3 variants that we classified as HCM-linked or nonpathogenic according to cosegregation and population genetics criteria. We found that around half of the HCM-linked variants showed alterations in RNA splicing or protein stability, both of which can lead to cMyBP-C haploinsufficiency. These protein haploinsufficiency drivers associated with HCM pathogenicity with 100% and 94% specificity, respectively. Furthermore, we uncovered that 11% of nontruncating MYBPC3 variants currently classified as of uncertain significance in ClinVar induced one of these molecular phenotypes. Our strategy, which can be applied to other conditions induced by protein loss of function, supports the idea that cMyBP-C haploinsufficiency is a fundamental pathomechanism in HCM.

Computational prediction of protein subdomain stability in MYBPC3 enables clinical risk stratification in hypertrophic cardiomyopathy and enhances variant interpretation

Purpose

Variants in MYBPC3 causing loss of function are the most common cause of hypertrophic cardiomyopathy (HCM). However, a substantial number of patients carry missense variants of uncertain significance (VUS) in MYBPC3. We hypothesize that a structural-based algorithm, STRUM, which estimates the effect of missense variants on protein folding, will identify a subgroup of HCM patients with a MYBPC3 VUS associated with increased clinical risk.

Methods

Among 7,963 patients in the multicenter Sarcomeric Human Cardiomyopathy Registry (SHaRe), 120 unique missense VUS in MYBPC3 were identified. Variants were evaluated for their effect on subdomain folding and a stratified time-to-event analysis for an overall composite endpoint (first occurrence of ventricular arrhythmia, heart failure, all-cause mortality, atrial fibrillation, and stroke) was performed for patients with HCM and a MYBPC3 missense VUS.

Results

We demonstrated that patients carrying a MYBPC3 VUS predicted to cause subdomain misfolding (STRUM+, ΔΔG ≤ −1.2 kcal/mol) exhibited a higher rate of adverse events compared with those with a STRUM- VUS (hazard ratio = 2.29, P = 0.0282). In silico saturation mutagenesis of MYBPC3 identified 4,943/23,427 (21%) missense variants that were predicted to cause subdomain misfolding.

Conclusion

STRUM identifies patients with HCM and a MYBPC3 VUS who may be at higher clinical risk and provides supportive evidence for pathogenicity.

Transcriptome Sequencing of Patients With Hypertrophic Cardiomyopathy Reveals Novel Splice-Altering Variants in MYBPC3

BACKGROUND

Transcriptome sequencing can improve genetic diagnosis of Mendelian diseases but requires access to tissue expressing disease-relevant transcripts. We explored genetic testing of hypertrophic cardiomyopathy using transcriptome sequencing of patient-specific human induced pluripotent stem cell derived cardiomyocytes (hiPSC-CMs). We also explored whether antisense oligonucleotides (AOs) could inhibit aberrant mRNA splicing in hiPSC-CMs.

METHODS

We derived hiPSC-CMs from patients with hypertrophic cardiomyopathy due to MYBPC3 splice-gain variants, or an unresolved genetic cause. We used transcriptome sequencing of hiPSC-CM RNA to identify pathogenic splicing and used AOs to inhibit this splicing.

RESULTS

Transcriptome sequencing of hiPSC-CMs confirmed aberrant splicing in 2 people with previously identified MYBPC3 splice-gain variants (c.1090+453C>T and c.1224-52G>A). In a patient with an unresolved genetic cause of hypertrophic cardiomyopathy following genome sequencing, transcriptome sequencing of hiPSC-CMs revealed diverse cryptic exon splicing due to an MYBPC3 c.1928-569G>T variant, and this was confirmed in cardiac tissue from an affected sibling. Antisense oligonucleotide treatment demonstrated almost complete inhibition of cryptic exon splicing in one patient-specific hiPSC-CM line.

CONCLUSIONS

Transcriptome sequencing of patient specific hiPSC-CMs solved a previously undiagnosed genetic cause of hypertrophic cardiomyopathy and may be a useful adjunct approach to genetic testing. Antisense oligonucleotide inhibition of cryptic exon splicing is a potential future personalized therapeutic option.

Making waves: A proposed new role for myosin-binding protein C in regulating oscillatory contractions in vertebrate striated muscle

Myosin-binding protein C (MyBP-C) is a critical regulator of muscle performance that was first identified through its strong binding interactions with myosin, the force-generating protein of muscle. Almost simultaneously with its discovery, MyBP-C was soon found to bind to actin, the physiological catalyst for myosin’s activity. However, the two observations posed an apparent paradox, in part because interactions of MyBP-C with myosin were on the thick filament, whereas MyBP-C interactions with actin were on the thin filament. Despite the intervening decades since these initial discoveries, it is only recently that the dual binding modes of MyBP-C are becoming reconciled in models that place MyBP-C at a central position between actin and myosin, where MyBP-C alternately stabilizes a newly discovered super-relaxed state (SRX) of myosin on thick filaments in resting muscle and then prolongs the “on” state of actin on thin filaments in active muscle. Recognition of these dual, alternating functions of MyBP-C reveals how it is central to the regulation of both muscle contraction and relaxation. The purpose of this Viewpoint is to briefly summarize the roles of MyBP-C in binding to myosin and actin and then to highlight a possible new role for MyBP-C in inducing and damping oscillatory waves of contraction and relaxation. Because the contractile waves bear similarity to cycles of contraction and relaxation in insect flight muscles, which evolved for fast, energetically efficient contraction, the ability of MyBP-C to damp so-called spontaneous oscillatory contractions (SPOCs) has broad implications for previously unrecognized regulatory mechanisms in vertebrate striated muscle. While the molecular mechanisms by which MyBP-C can function as a wave maker or a wave breaker are just beginning to be explored, it is likely that MyBP-C dual interactions with both myosin and actin will continue to be important for understanding the new functions of this enigmatic protein.

MYBPC3 Haplotype Linked to Hypertrophic Cardiomyopathy in Rhesus Macaques (Macaca mulatta)

In humans, abnormal thickening of the left ventricle of the heart clinically defines hypertrophic cardiomyopathy (HCM), a common inherited cardiovascular disorder that can precede a sudden cardiac death event. The wide range of clinical presentations in HCM obscures genetic variants that may influence an individual's susceptibility to sudden cardiac death. Although exon sequencing of major sarcomere genes can be used to detect high-impact causal mutations, this strategy is successful in only half of patient cases. The incidence of left ventricular hypertrophy (LVH) in a managed research colony of rhesus macaques provides an excellent comparative model in which to explore the genomic etiology of severe HCM and sudden cardiac death. Because no rhesus HCM-associated mutations have been reported, we used a next-generation genotyping assay that targets 7 sarcomeric rhesus genes within 63 genomic sites that are orthologous to human genomic regions known to harbor HCM disease variants. Amplicon sequencing was performed on 52 macaques with confirmed LVH and 42 unrelated, unaffected animals representing both the Indian and Chinese rhesus macaque subspecies. Bias-reduced logistic regression uncovered a risk haplotype in the rhesus MYBPC3 gene, which is frequently disrupted in both human and feline HCM; this haplotype implicates an intronic variant strongly associated with disease in either homozygous or carrier form. Our results highlight that leveraging evolutionary genomic data provides a unique, practical strategy for minimizing population bias in complex disease studies.

Under construction: The dynamic assembly, maintenance, and degradation of the cardiac sarcomere

The sarcomere is the basic contractile unit of striated muscle and is a highly ordered protein complex with the actin and myosin filaments at its core. Assembling the sarcomere constituents into this organized structure in development, and with muscle growth as new sarcomeres are built, is a complex process coordinated by numerous factors. Once assembled, the sarcomere requires constant maintenance as its continuous contraction is accompanied by elevated mechanical, thermal, and oxidative stress, which predispose proteins to misfolding and toxic aggregation. To prevent protein misfolding and maintain sarcomere integrity, the sarcomere is monitored by an assortment of protein quality control (PQC) mechanisms. The need for effective PQC is heightened in cardiomyocytes which are terminally differentiated and must survive for many years while preserving optimal mechanical output. To prevent toxic protein aggregation, molecular chaperones stabilize denatured sarcomere proteins and promote their refolding. However, when old and misfolded proteins cannot be salvaged by chaperones, they must be recycled via degradation pathways: the calpain and ubiquitin-proteasome systems, which operate under basal conditions, and the stress-responsive autophagy-lysosome pathway. Mutations to and deficiency of the molecular chaperones and associated factors charged with sarcomere maintenance commonly lead to sarcomere structural disarray and the progression of heart disease, highlighting the necessity of effective sarcomere PQC for maintaining cardiac function. This review focuses on the dynamic regulation of assembly and turnover at the sarcomere with an emphasis on the chaperones involved in these processes and describes the alterations to chaperones - through mutations and deficient expression - implicated in disease progression to heart failure.

A Novel "Cut and Paste" Method for In Situ Replacement of cMyBP-C Reveals a New Role for cMyBP-C in the Regulation of Contractile Oscillations

RATIONALE

cMyBP-C (cardiac myosin-binding protein-C) is a critical regulator of heart contraction, but the mechanisms by which cMyBP-C affects actin and myosin are only partly understood. A primary obstacle is that cMyBP-C localization on thick filaments may be a key factor defining its interactions, but most in vitro studies cannot duplicate the unique spatial arrangement of cMyBP-C within the sarcomere.

OBJECTIVE

The goal of this study was to validate a novel hybrid genetic/protein engineering approach for rapid manipulation of cMyBP-C in sarcomeres in situ.

METHODS AND RESULTS

We designed a novel cut and paste approach for removal and replacement of cMyBP-C N'-terminal domains (C0-C7) in detergent-permeabilized cardiomyocytes from gene-edited Spy-C mice. Spy-C mice express a TEVp (tobacco etch virus protease) cleavage site and a SpyTag (st) between cMyBP-C domains C7 and C8. A cut is achieved using TEVp which cleaves cMyBP-C to create a soluble N'-terminal γC0C7 (endogenous [genetically encoded] N'-terminal domains C0 to C7 of cardiac myosin binding protein-C) fragment and an insoluble C'-terminal SpyTag-C8-C10 fragment that remains associated with thick filaments. Paste of new recombinant (r)C0C7 domains is achieved by a covalent bond formed between SpyCatcher (-sc; encoded at the C'-termini of recombinant proteins) and SpyTag. Results show that loss of γC0C7 reduced myofilament Ca²⁺ sensitivity and increased cross-bridge cycling (k_tr) at submaximal [Ca²⁺]. Acute loss of γC0C7 also induced auto-oscillatory contractions at submaximal [Ca²⁺]. Ligation of rC0C7 (exogenous [recombinant] N'-terminal domains C0 to C7 of cardiac myosin binding protein-C)-sc returned pCa₅₀ and k_tr to control values and abolished oscillations, but phosphorylated (p)-rC0C7-sc did not completely rescue these effects.

CONCLUSIONS

We describe a robust new approach for acute removal and replacement of cMyBP-C in situ. The method revealed a novel role for cMyBP-C N'-terminal domains to damp sarcomere-driven contractile waves (so-called spontaneous oscillatory contractions). Because phosphorylated (p)-rC0C7-sc was less effective at damping contractile oscillations, results suggest that spontaneous oscillatory contractions may contribute to enhanced contractility in response to inotropic stimuli.

Human Induced Pluripotent Stem-Cell-Derived Cardiomyocytes as Models for Genetic Cardiomyopathies

In the last few decades, many pathogenic or likely pathogenic genetic mutations in over hundred different genes have been described for non-ischemic, genetic cardiomyopathies. However, the functional knowledge about most of these mutations is still limited because the generation of adequate animal models is time-consuming and challenging. Therefore, human induced pluripotent stem cells (iPSCs) carrying specific cardiomyopathy-associated mutations are a promising alternative. Since the original discovery that pluripotency can be artificially induced by the expression of different transcription factors, various patient-specific-induced pluripotent stem cell lines have been generated to model non-ischemic, genetic cardiomyopathies in vitro. In this review, we describe the genetic landscape of non-ischemic, genetic cardiomyopathies and give an overview about different human iPSC lines, which have been developed for the disease modeling of inherited cardiomyopathies. We summarize different methods and protocols for the general differentiation of human iPSCs into cardiomyocytes. In addition, we describe methods and technologies to investigate functionally human iPSC-derived cardiomyocytes. Furthermore, we summarize novel genome editing approaches for the genetic manipulation of human iPSCs. This review provides an overview about the genetic landscape of inherited cardiomyopathies with a focus on iPSC technology, which might be of interest for clinicians and basic scientists interested in genetic cardiomyopathies.

SarcTrack

RATIONALE

Human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) in combination with CRISPR/Cas9 genome editing provide unparalleled opportunities to study cardiac biology and disease. However, sarcomeres, the fundamental units of myocyte contraction, are immature and nonlinear in hiPSC-CMs, which technically challenge accurate functional interrogation of contractile parameters in beating cells. Furthermore, existing analysis methods are relatively low-throughput, indirectly assess contractility, or only assess well-aligned sarcomeres found in mature cardiac tissues.

OBJECTIVE

We aimed to develop an analysis platform that directly, rapidly, and automatically tracks sarcomeres in beating cardiomyocytes. The platform should assess sarcomere content, contraction and relaxation parameters, and beat rate.

METHODS AND RESULTS

We developed SarcTrack, a MatLab software that monitors fluorescently tagged sarcomeres in hiPSC-CMs. The algorithm determines sarcomere content, sarcomere length, and returns rates of sarcomere contraction and relaxation. By rapid measurement of hundreds of sarcomeres in each hiPSC-CM, SarcTrack provides large data sets for robust statistical analyses of multiple contractile parameters. We validated SarcTrack by analyzing drug-treated hiPSC-CMs, confirming the contractility effects of compounds that directly activate (CK-1827452) or inhibit (MYK-461) myosin molecules or indirectly alter contractility (verapamil and propranolol). SarcTrack analysis of hiPSC-CMs carrying a heterozygous truncation variant in the myosin-binding protein C ( MYBPC3) gene, which causes hypertrophic cardiomyopathy, recapitulated seminal disease phenotypes including cardiac hypercontractility and diminished relaxation, abnormalities that normalized with MYK-461 treatment.

CONCLUSIONS

SarcTrack provides a direct and efficient method to quantitatively assess sarcomere function. By improving existing contractility analysis methods and overcoming technical challenges associated with functional evaluation of hiPSC-CMs, SarcTrack enhances translational prospects for sarcomere-regulating therapeutics and accelerates interrogation of human cardiac genetic variants.

Altered force generation and cell-to-cell contractile imbalance in hypertrophic cardiomyopathy

Hypertrophic cardiomyopathy (HCM) is mainly caused by mutations in sarcomeric proteins. Thirty to forty percent of identified mutations are found in the ventricular myosin heavy chain (β-MyHC). A common mechanism explaining how numerous mutations in several different proteins induce a similar HCM-phenotype is unclear. It was proposed that HCM-mutations cause hypercontractility, which for some mutations is thought to result from mutation-induced unlocking of myosin heads from a so-called super-relaxed state (SRX). The SRX was suggested to be related to the "interacting head motif," i.e., pairs of myosin heads folded back onto their S2-region. Here, we address these structural states of myosin in context of earlier work on weak binding cross-bridges. However, not all HCM-mutations cause hypercontractility and/or are involved in the interacting head motif. But most likely, all mutations alter the force generating mechanism, yet in different ways, possibly including inhibition of SRX. Such functional-hyper- and hypocontractile-changes are the basis of our previously proposed concept stating that contractile imbalance due to unequal fractions of mutated and wildtype protein among individual cardiomyocytes over time will induce cardiomyocyte disarray and fibrosis, hallmarks of HCM. Studying β-MyHC-mutations, we found substantial contractile variability from cardiomyocyte to cardiomyocyte within a patient's myocardium, much higher than in controls. This was paralleled by a similarly variable fraction of mutant MYH7-mRNA (cell-to-cell allelic imbalance), due to random, burst-like transcription, independent for mutant and wildtype MYH7-alleles. Evidence suggests that HCM-mutations in other sarcomeric proteins follow the same disease mechanism.