This article has Open Peer Review reports available.
A substitution mutation in cardiac ubiquitin ligase, FBXO32, is associated with an autosomal recessive form of dilated cardiomyopathy
© Al-Hassnan et al. 2016
Received: 27 February 2015
Accepted: 6 January 2016
Published: 14 January 2016
Familial dilated cardiomyopathy (DCM) is genetically heterogeneous. Mutations in more than 40 genes have been identified in familial cases, mostly inherited in an autosomal dominant pattern. DCM due to recessive mutations is rarely observed. In consanguineous families, homozygosity mapping and whole exome sequencing (WES) can be utilized to identify the genetic defects in recessively inherited DCM.
In a consanguineous family with four affected siblings with severe DCM, we combined homozygosity mapping, linkage analysis and WES, to uncover the genetic defect.
A region of homozygosity (ROH) on chromosome 8q24.13–24.23 was found to be shared by all of the four affected siblings. WES detected ~47,000 variants that were filtered to a homozygous mutation (p.Gly243Arg) in the FBXO32 gene, located within the identified ROH. The mutation segregated with the phenotype, replaced a highly-conserved amino acid, and was not detected in 1986 ethnically-matched chromosomes. FBXO32, which encodes a muscle-specific ubiquitin ligase, has been implicated in the pathogenesis of cardiomyopathy through the ubiquitin proteasome system (UPS). In addition, FBXO32-knockout mice manifest with cardiomyopathy. Screening the index patient for all of the WES variants in 48 genes known to be implicated in hypertrophic and dilated cardiomyopathy was negative.
Our data suggest that FBXO32 is a candidate gene for recessive DCM. Acting as a cardiac ubiquitin ligase, mutated FBXO32 could perturb the degradation of target proteins in the UPS, the impairment of which has been observed in cardiomyopathy. Our work proposes that genes encoding other ubiquitin ligases could also be implicated in familial cardiomyopathy.
Dilated cardiomyopathy (DCM; OMIM 115200), characterized by impaired systolic function and left ventricular enlargement, is considered to be the most common indication for cardiac transplantation [1, 2]. A familial pattern, mostly inherited as autosomal dominant, has been observed in about 20–35 % of DCM cases with remarkable genetic heterogeneity; disease-associated mutations have been identified in more than 40 genes [3, 4]. A recessively inherited form of familial DCM is rarely observed and is often associated with extra-cardiac manifestations such as skeletal myopathy, hypotonia, hepatic encephalopathy, and impaired fatty acid oxidation defects.
The ubiquitin proteasome system (UPS) is crucial for the regulation of cellular protein degradation. The target proteins for degradation undergo an initial step of polyubiquitination before their destruction by the 26S proteasome . The UPS cascade consists of three components: a ubiquitin-activating enzyme (E1), a ubiquitin-conjugating enzyme (E2), and a ubiquitin-ligase enzyme (E3), which is the substrate-specific component of the UPS. Impaired UPS has recently been implicated in the pathogenesis of heart failure, cardiac ischemia, and cardiomyopathy (CMP) [6, 7]. The protein encoded by FBXO32, also known as atrogin-1 and MAFbx, is an E3 ligase that is expressed selectively in skeletal muscles and cardiomyocytes and plays a critical role in muscle atrophy [8, 9], as well as in the development of cardiac hypertrophy and atrophy [10, 11]. In support of an association between deficient FBXO32 protein and CMP, a recent study has shown that atrogin-1 knockout mice develop CMP with intracellular protein accumulation and cardiomyocyte apoptosis .
Subjects with cardiomyopathy
Echocardiographic findings of the four affected individuals. EF: ejection fraction; FS: fractional shortening; LV: left ventricle; LVDd: left ventricular diastolic diameter; MR: mitral regurgitation; RV: right ventricle; TR: tricuspid regurgitation; yrs: years
severely dilated LV. severely reduced LV systolic function. moderate MR.
moderately dilated LV. severely reduced LV systolic function. trace MR.
moderately dilated LV. severely reduced LV systolic function. moderate MR.
severely dilated LV and RV. severely reduced LV and RV systolic function. severe MR and TR
Homozygosity mapping and linkage analysis
A recessive pattern of inheritance was assumed, because the family was consanguineous, the parents had unremarkable cardiac evaluations, and the disease affected both sexes. Genomic DNA was extracted from whole-blood samples obtained from all affected patients, their parents, and unaffected siblings. SNP-based genotyping was performed using AxiomTM CEU Human Array from AffymetrixR (Santa Clara, CA, USA). To detect regions of homozygosity (ROHs), SNPs were analyzed using AutoSNPa software . Conventionally, ROH are defined as fragments where SNPs are homozygous for a stretch of consecutive alleles in affected individuals, and heterozygous, or homozygous for the other alleles, in unaffected members of the same family. Multipoint parametric linkage analysis was performed using Gene Hunter Easy Linkage Analysis Software module 4.0. A fully penetrant recessive model of inheritance was used, with a population disease allele frequency of 0.0001 and using Asian SNP allele frequencies.
Whole exome sequencing
Candidate gene search and mutation screening
Using the UCSC and Ensembl genome browsers [14, 15], the OMIM database , and the GeneDistiller2 software , candidate genes were analyzed based on their function, expression, biological pathways, or animal models. Genomic DNA of the affected patients, their parents, and siblings was amplified by PCR using intronic primers that were designed to flank (50–100 bp) the coding exons of the identified gene, as defined by the Ensembl Genome Browser . PCR was performed in a final volume of 20 μl containing approximately 10 ng of genomic DNA, using standard conditions (primer sequences and conditions are available on request). Purified PCR amplicons covering the entire coding region of the selected genes were directly sequenced with the dideoxy chain-termination method using an ABI Prism Big Dye Terminator v3.1 Cycle Sequencing Kit, following the manufacturer’s instructions, and processed on an ABI 3730XL capillary sequencer (Applied Biosystems, CA, USA). Sequence analysis was performed using the SeqMan 6.1 module of the Lasergene (DNA Star Inc. WI, USA) software package, then compared with the reference GenBank sequence. Numbering commenced with the A of the ATG initiation codon as +1.
Histopathology and immunohistochemistry
Histopathological and immunohistochemical studies were performed on ventricular myocardial tissue from the explanted heart of the proband, and on a control myocardial specimen from a non-cardiomyopathy case. To assess FBXO32 expression, a mouse monoclonal anti-human FBXO32 antibody raised against amino acids 1–300 (Santa Cruz Biotechnology, Dallas, Texas) was used for immunohistochemical staining. Antibody dilution was performed according to the manufacturer’s protocol.
Homozygosity mapping and linkage analysis
Whole exome sequencing
We have also screened all the WES variants that overlap 48 OMIM genes implicated in dilated and hypertrophic familial CMP (Additional file 1: Table S1). All of the identified 62 variants have been previously reported in the dbSNP and/or 1000 Genomes databases as polymorphisms, except for two intronic variants in the TTN and SGCD genes.
Bioinformatics-based prediction analysis for functional consequences of p.Gly243Arg
Several bioinformatics tools (Polyphen2 , PANTHER , SNPs&GO , and I-Mutant ) were utilized to predict the functional consequences of p.Gly243Arg (c.727G > C). The Polyphen2 program predicted the variant to be probably damaging. According to the PANTHER cSNP estimation, the change has a deleterious effect on protein function with a score of −3.624. SNPs&GO analysis predicted the p.Gly243Arg change as “disease” with a reliability score of 1, while I-Mutant predicted that the change will reduce the stability of the protein with a score of −0.67 (the “minus” score indicates reduced stability). Protein sequence alignment of FBXO32 orthologs demonstrated that the glycine residue is a conserved amino acid down to zebrafish , suggesting that this amino acid has an essential function (Fig. 5b). We also modeled p.Gly243Arg of FBXO32 on a template structure adapted from 2ovq (Fig. 5c). The analysis suggested that the replacement of the small and more hydrophobic glycine with the large, polar and less hydrophobic arginine may not only hinder normal protein folding, but also destabilize the protein, as was also predicted by the I-Mutant algorithm.
Histopathology and immunohistochemistry
In this work, we combined homozygosity mapping, linkage analysis, and WES, to identify the genetic alteration in a consanguineous family with a recessive form of DCM. The identified homozygous mutation in the FBXO32 gene was uncovered after homozygosity mapping, an approach that we adopted in this case in view of the high probability of identifying a homozygous mutation in a consanguineous family. The parents were asymptomatic with normal echocardiograms, and the disease affected both sexes.
The phenotype in the proband was remarkable, with severe CMP progressing to heart failure that necessitated heart transplantation. Screening family members by echocardiogram revealed variable presentation: two siblings were asymptomatic while the youngest affected child had suggestive cardiac symptoms. This observation underscores the importance of screening first-degree relatives in cases of DCM.
The FXO32 gene is one of more than 38 members of the FBXO (F-Box Only) family of proteins that have an F-box domain characterized by approximately 50 amino acids, which functions as a site for protein-protein interaction . FBXO32 is one of the components of the E3 ubiquitin ligase SCF (Skp1, Cullin-1, F-box), which binds target proteins for degradation by the UPS. The F-box domain links the F-box protein to the other SCF components by binding Skp1, leading to the destruction of the target proteins. In animals, through enhancing protein degradation, Fbxo32 was initially discovered to play a critical role in inducing muscle atrophy [8, 9].
Serving as a cardiac ubiquitin ligase, FBXO32 has been implicated in the pathogenesis of CMP.6,7 FBXO32 has been shown to repress calcineurin through assembling with the UPS complex, which promotes cardiac hypertrophy in response to pathologic stimuli . Li et al. have also found that over-expressed Fbxo32 in neonatal rat cardiomyocytes disrupts the Akt-dependent pathway responsible for physiological cardiac hypertrophy.11 Interestingly, down-regulation of Fbxo32 in knockout mice produces the opposite effect, with the inhibition of cardiac hypertrophy in response to pressure overload , suggesting that there are different mechanisms by which FBXO32 induces atrophy and hypertrophy in cardiomyocytes . In a recent study on atrogin-1 knockout mice, charged multivesicular body protein 2B (CHMP2B), which is part of an endosomal sorting complex required for autophagy, has been identified as a target of atrogin-1-mediated degradation.12 Mice lacking atrogin-1 fail to degrade CHMP2B, resulting in autophagy impairment, intracellular protein aggregate accumulation, unfolded protein response activation and subsequent cardiomyocyte apoptosis, leading, ultimately, to CMP and premature death.
Furthermore, in support of a potential link between FBXO32 and DCM in this family, other cardiac ubiquitin ligases (MuRF1, MuRF2, MuRF3, CHIP, MDM2) have also been implicated in the pathogenesis of cardiac hypertrophy, atrophy, and ischemia reperfusion injury.5 Recently, mutations in the cardiac ubiquitin ligase TRIM63, which encodes MuRF1, have been identified in patients with familial hypertrophic CMP . The expression of mutant TRIM63 was associated with impaired UPS-mediated protein degradation in cardiomyocytes. This finding, together with our work, suggests that the dysfunction of other proteins in the UPS complexes may also be implicated in the pathogenesis of Mendelian forms of CMP.
Although causality is yet to be established, several factors pertinent to the identified mutation (p.Gly243Arg) support the suggestion that FBXO32 is a likely candidate gene for DCM. There was a clear co-segregation of the mutation with the clinical phenotype. All of the affected patients, but none of their parents and unaffected siblings, were homozygous for the mutant allele. In addition, the mutation, which replaces a highly-conserved aliphatic non-polar amino acid with a polar positively charged one, was not present in 1986 chromosomes, and was predicted to destabilize the protein.
The identified homoallelic variant in the F-box domain is likely to be a loss-of-function mutation. Being within the F-box domain, this mutation may cause loss of substrate specificity, leading to the premature degradation of functional proteins. It may also cause loss of efficient recruitment of the components of the UPS. However, a gain-of-function mechanism is still a possibility, where the mutant allele could lead to the accumulation of proteins that are damaging. Both mechanisms would theoretically perturb the well-controlled system of protein homeostasis within cardiomyocytes, and lead to CMP.
It is noteworthy that none of the affected patients had skeletal muscle weakness. The mutation identified in this family may have damaged a cardiac-specific protein, and hence the skeletal muscles are spared. However, further studies in cell or animal models harboring this particular mutation may provide a plausible explanation of the mechanism of cardiac disease and of the apparent lack of skeletal involvement.
The methodology that we have adopted in this study has recognizable limitations. The family pedigree was highly suggestive of a recessive pattern of inheritance, and therefore we implemented the approach of homozygosity mapping and WES, followed by excluding non-homozygous variants. Although the recessive inheritance conforms to the segregation of the disease in this family, a dominant pattern with reduced penetrance remains a possibility. In addition, the unaffected individuals may still be in a presymptomatic phase of a late-onset CMP. A long-term follow-up with periodic evaluation may clarify if the identified variant is not fully penetrant. Last, as our study was conducted in a single family, screening FBXO32 in other families with CMP of unknown etiology may provide more insight and establish the causal relationship between FBXO32 and CMP.
Our work suggests that FBXO32 is a candidate gene for recessive DCM. Further studies are still needed to establish the causal relationship. Genetic analysis for FBXO32 may be considered in families with recessively inherited DCM of unknown genetic cause.
The authors would like to thank the patients and their family members for their participation in the study. We would also like to thank members of the sequencing and genotyping core facilities in the Department of Genetics at KFSH&RC. This work was supported by a grant from King Abdulaziz City for Science and Technology (KACST) and the National Comprehensive Plan for Science and Technology (NCPST) (grant numbers: 08-MED489-20 and 11-MED1439-20).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Taylor DO, Edwards LB, Boucek MM, Trulock EP, Aurora P, Christie J, et al. Registry of the International Society for Heart and Lung Transplantation: twenty-fourth official adult heart transplant report—2007. J Heart Lung Transplant. 2007;26:769–81.View ArticlePubMedGoogle Scholar
- Boucek MM, Aurora P, Edwards LB, Taylor DO, Trulock EP, Christie J, et al. Registry of the International Society for Heart and Lung Transplantation: tenth official pediatric heart transplantation report--2007. J Heart LungTransplant. 2007;26:796–807.View ArticleGoogle Scholar
- Hershberger RE, Morales A, Siegfried JD. Clinical and genetic issues in dilated cardiomyopathy: a review for genetics professionals. Genet Med. 2010;12:655–67.View ArticlePubMedPubMed CentralGoogle Scholar
- Online Mendelian Inheritance in Man. http://omim.org. (accessed 22 June 2014).
- Willis MS, Townley-Tilson WH, Kang EY, Homeister JW, Patterson C. Sent to destroy: the ubiquitin proteasome system regulates cell signaling and protein quality control in cardiovascular development and disease. Circ Res. 2010;106(3):463–78.View ArticlePubMedPubMed CentralGoogle Scholar
- Schlossarek S, Carrier L. The ubiquitin-proteasome system in cardiomyopathies. Curr Opin Cardiol. 2011;26(3):190–95.View ArticlePubMedGoogle Scholar
- Mearini G, Gedicke C, Schlossarek S, Witt CC, Kramer E, Cao P, et al. Atrogin-1 andMuRF1 regulate cardiac MyBP-C levels via different mechanisms. Cardiovasc Res. 2010;85(2):357–66.View ArticlePubMedGoogle Scholar
- Bodine SC, Latres E, Baumhueter S, Lai VK, Nunez L, Clarke BA, et al. Identification of ubiquitin ligases required for skeletal muscle atrophy. Science. 2001;294:1704–08.View ArticlePubMedGoogle Scholar
- Gomes MD, Lecker SH, Jagoe RT, Navon A, Goldberg AL. Atrogin-1, a muscle-specific Fbox protein highly expressed during muscle atrophy. Proc Natl Acad Sci USA. 2001;98:14440–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Li HH, Kedar V, Zhang C, McDonough H, Arya R, Wang DZ, et al. F-box inhibits calcineurin-dependent cardiac hypertrophy by participating in an SCF ubiquitin ligase complex. J Clin Invest. 2004;114(8):1058–71.View ArticlePubMedPubMed CentralGoogle Scholar
- Li HH, Willis MS, Lockyer P, Miller N, McDonough H, Glass DJ, et al. Atrogin-1inhibits Akt dependent cardiac hypertrophy in mice via ubiquitin-dependent coactivation of Forkhead proteins. J Clin Invest. 2007;117:3211–23.View ArticlePubMedPubMed CentralGoogle Scholar
- Zaglia T, Milan G, Ruhs A, Franzoso M, Bertaggia E, Pianca N, et al. Atrogin-1 deficiency promotes cardiomyopathy and premature death via impaired autophagy. J Clin Invest. 2014;124(6):2410–24.View ArticlePubMedPubMed CentralGoogle Scholar
- Carr IM, Flintoff KJ, Taylor GR, Markham AF, Bonthron DT. Interactive visual analysis of SNP data for rapid autozygosity mapping in consanguineous families. Hum Mut. 2006;27:1041–6.View ArticlePubMedGoogle Scholar
- UCSC Genome Browser. http://genome.ucsc.edu. (accessed 22 June 2014).
- Ensembl Genome Browser. http://www.ensembl.org/index.html. (accessed 22 June 2014).
- GeneDistiller2 software; http://www.genedistiller.org. (accessed 22 June 2014).
- Wang Y, Kowalski PE, Thalmann I, Ornitz DM, Mager DL, Thalmann R. Otoconin-90, the mammalian otoconial matrix protein, contains two domains of homology to secretory phospholipase A2. Proc Natl Acad Sci USA. 1998;95(26):15345–50.View ArticlePubMedPubMed CentralGoogle Scholar
- Mouse Genome Informatics URL: http://www.informatics.jax.org. (accessed 22 June 2014).
- Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, et al. A method and server for predicting damaging missense mutations. Nat Methods. 2010;7(4):248–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Mi H, Guo N, Kejariwal A, Thomas PD. PANTHER version 6: protein sequence and function evolution data with expanded representation of biological pathways. Nucleic Acids Res. 2007;35(Database issue):D247–52.View ArticlePubMedGoogle Scholar
- Calabrese R, Capriotti E, Fariselli P, Martelli PL, Casadio R. Functional annotations improve the predictive score of human disease-related mutations in proteins. Hum Mutat. 2009;30(8):1237–44.View ArticlePubMedGoogle Scholar
- Capriotti E, Fariselli P, Casadio R. I-Mutant2.0: predicting stability changes upon mutation from the protein sequence or structure. Nucleic Acids Res. 2005;33(WebServer issue):W306–10.View ArticlePubMedPubMed CentralGoogle Scholar
- Kipreos ET, Pagano M. The F-box protein family. Genome Biol. 2000;1(5):REVIEWS3002.View ArticlePubMedPubMed CentralGoogle Scholar
- Usui S, Maejima Y, Pain J, Hong C, Cho J, Park JY, et al. Endogenous muscle atrophy F-box mediates pressure overload-induced cardiac hypertrophy through regulation of nuclear factor-kappaB. Circ Res. 2011;109(2):161–71.View ArticlePubMedPubMed CentralGoogle Scholar
- Lee D, Goldberg A. Atrogin1/MAFbx: what atrophy, hypertrophy, and cardiac failure have in common. Circ Res. 2011;109(2):123–6.View ArticlePubMedGoogle Scholar
- Chen SN, Czernuszewicz G, Tan Y, Lombardi R, Jin J, Willerson JT, et al. Human Molecular Genetic and Functional Studies Identify TRIM63, Encoding Muscle RING Finger Protein 1, as a Novel Gene for Human Hypertrophic Cardiomyopathy. Circ Res. 2012;111(7):907–19.View ArticlePubMedPubMed CentralGoogle Scholar