Whole exome sequencing identifies a novel EMDmutation in a Chinese family with dilated cardiomyopathy
- Mingqiu Zhang†2,
- Jia Chen†1,
- Dayong Si3,
- Yu Zheng1,
- Haixu Jiao2,
- Zhaohui Feng2,
- Zhengmao Hu1Email author and
- Ranhui Duan1Email author
© Zhang et al.; licensee BioMed Central Ltd. 2014
Received: 21 April 2014
Accepted: 27 June 2014
Published: 5 July 2014
Variants in the emerin gene (EMD) were implicated in X-linked recessive Emery-Dreifuss muscular dystrophy (EDMD), characterized by early-onset contractures of tendons, progressive muscular weakness and cardiomyopathy. To date, 223 mutations have been reported in EMD gene and the majority of them caused a predominant skeletal muscular phenotype. In this study, we identified a novel deletion mutation in EMD exon 1, which results in almost a complete loss of emerin protein in a large Chinese family. However, the patients suffered severe dilated cardiomyopathy (DCM) but very mild skeletal muscle disorder.
Whole exome sequencing (WES) and linkage analysis were performed to identify the underlying mutation in a Chinese DCM family spanning five generations. A missense variation in the GPR50 gene was found co-segregated with the disease phenotype, whereas no functional alteration was detected in the variant GPR50 protein. When analyzing the failure sequences in the exome sequencing data, a novel deletion mutation (c.26_39delATACCGAGCTGACC) in EMD exon 1, was identified in this family. Different from the typical clinical features caused by most reported EMD mutations, patients in our study presented very mild skeletal muscle degeneration that had not been diagnosed until the mutation was found.
We described a family with rare clinical presentations caused by a novel EMD deletion mutation. Our findings broaden the heterogeneous spectrum of phenotypes attributed to EMD mutations and provide new insight to explain the genotype-phenotype correlations between EMD mutations and EDMD symptoms.
KeywordsDilated cardiomyopathy Emery-Dreifuss muscular dystrophy EMD gene Novel mutation Whole exome sequencing
The EMD gene, consisting of 6 exons and encoding a nuclear envelope protein emerin (254 amino acids), was involved in X-linked Emery-Dreifuss muscular dystrophy (EDMD) . Among the 223 reported EMD mutations in the Leiden Open Variation Database (LOVD), the overwhelming majority of mutations have been identified in EDMD patients who presented childhood-onset contractures of the elbows, Achilles tendons and post-cervical muscles, slowly progressive muscle weakness with humeroperoneal distribution in adolescence and cardiomyopathy with conduction block in early adulthood [1, 2].
Dilated cardiomyopathy (DCM) is a common cause of advanced heart failure and the primary indication for heart transplantation . Currently over 40 individual genes that encode for sarcomeric, cytoskeletal, nuclear membrane, dystrophin-associated glycoprotein complex and desmosomal proteins have been implicated in inherited DCM [3, 4]. One of the most common nonsyndromic DCM-causing genes is LMNA, which encodes the nuclear lamina proteins lamin A and C; mutations in this gene also cause autosomal dominant and recessive EDMD [1, 4]. Even though the functions of emerin and lamin A, C are very similar, the EMD gene is not listed in the routine screening of DCM because of lacking sufficient body evidences .
In this study, we conducted WES combined with linkage analysis to identify the underlying mutation in a large Chinese DCM family. A missense variation in the GPR50 gene was found co-segregated with the disease phenotype, whereas no functional alteration was detected in the variant GPR50 protein. When analyzing the failure sequences in the WES data, a novel deletion mutation (c.26_39delATACCGAGCTGACC) in EMD exon1, which results in almost a complete loss of emerin protein, was identified in this family. Different from the typical clinical features caused by most reported EMD mutations, the patients in our study presented severe DCM but very mild skeletal atrophy.
Materials and methods
Peripheral blood of 39 family members was obtained. Genomic DNA was extracted from the peripheral blood lymphocytes, using the standard phenol–chloroform method. A total of 22 subjects were included in the linkage study and followed fine mapping analysis. Genotyping was performed using 18 microsatellite markers (Panel 28, Version 2.5, ABI PRISM Linkage Mapping Set) in X chromosome. Five additional markers (DXS998, DXS1193, DXS8069, DXS8011, DXS8103) were chosen to refine mapping the gene locus. Two-point LOD scores were calculated by the program MLINK, and multipoint LOD scores with use of Linkmap. The disease was considered to be X-linked with a frequency of 0.0001 and 95% penetrance.
Five micrograms of DNA from each of two affected male individuals was applied for the construction of exome library using the Agilent SureSelect exome capture system and was sequenced on the Illumina Solexa GAIIx platform following the manufacturer’s instructions (Illumina, San Diego, CA). Raw image files were processed by the Illumina pipeline (version 1.3.4) for base calling and generating the reads set.
Reads were aligned to the human reference genome (University of California Santa Cruz, UCSC hg19) using SOAPaligner. Single nucleotide polymorphisms (SNPs) and indels (insertions and deletions) identification was performed by samtools and/or Genome Analysis Toolkit (GATK) and such SNPs with read depth > 4 and quality > 20 were reserved for subsequent analysis. SNPs and indels were annotated using SeattleSeq annotation. On the basis of SNP database (dbSNP) and 1000 genomes annotation and the supplied PolyPhen prediction, any non-synonymous variants not assigned a ‘benign’ prediction were considered to be damaging.
Analysis of candidate genes
DNA sequences were obtained from GenBank. Primers were designed with primer-3 software. Protein-encoding exons of candidate genes were amplified from genomic DNA. PCR products were purified with use of kits (golden polymerase, ABI, California) and sequenced with a dye-terminator cycle-sequencing system (Applied Biosystems, California).
All the individuals in the family were interviewed, and reported symptoms of cardiomyopathy disease such as chest pain, shortness of breath and fatigue; 8 males and 4 females were identified with clinical symptoms (Figure 1A). Sudden cardiac death was reported in five of eight male patients at the age of 25, 30, 33, 33 and 42 years old. The female patients presented milder heart failure and less life-threatening symptoms than the males. The proband A (III: 19) was initially diagnosed as sick sinus syndrome with chest tightness and early onset fatigue, especially after heavy activity, and his electrocardiogram (ECG) showed sinus bradycardia (sinus node dysfunction), with a minimum heart rate of 36 times/min at the age of 17 years. The symptoms progressively became more severe when the patient was 34 years old. The pacemaker implantation failed and the ultrasonic cardiogram (UCG) showed chamber enlargement (left ventricular end-diastolic diameter of 56.8 mm, the right atrium 65 × 83 mm, right room 26.4 mm, left atrial 54.2 mm) (Figure 1B), with a thinning of the chamber wall. The subject was hospitalized at the age of 41 years; the 24-h Holter detected paroxysmal ventricular tachycardia and sudden death occurred 1 month after discharge. A similar cardiac disorder has been identified in other patients in the family, but not among unaffected members.
Candidate genes analysis and Linkage study
Genetic and clinical evaluation of the pedigree provides the evidence of an X-linked recessive inheritance. DMD and TAZ are the two reported genes associated with X-linked nonsyndromic DCM . No deletion or insertion was observed in the DMD gene of proband A using the method of multiplex ligation-dependent probe amplification (MLPA) and Sanger sequencing did not detect any mutation in the exons and exon-intron junctions of both strands in the TAZ gene. Linkage analyses in the family identified a locus associated with DCM on chromosome X, band Xq28 (maximal lod score 1.55 at DXS1073). Five additional markers were analyzed to narrow the positive candidate region, defined a common haplotype in all the affected subjects, yielding a maximum 2-point LOD score of 2.43 with a marker located at DXS1193 (Figure 1A).
Whole exome sequencing
Whole exomes of two affected male individuals in the family (sample III: 11 and IV: 6) were sequenced, generating an average of 4.5 billion bases of sequence and a mean coverage of 65× for each affected individual. An average of 97% of the targeted bases was sufficiently covered to pass our thresholds for variant calling. Variants of 73 hereditary cardiomyopathies related genes were screened , the filtering conditions are as following: (1) same variants in both WES data; (2) missense, nonsense, insertion and deletion variants; (3) SNPs with minor allele frequency not more than 0.005 according to the SNP database of National Center for Biotechnology Information (NCBI). 7 variants passed the filtering conditions and none of them co-segregated with the disease phenotype. Then 140 genes in the linked locus (Xq28) were also screened according to the above conditions. 9 variants passed the conditions and only the GPR50 c.113C > T variant was found co-segregation with the disease phenotype.
Candidate variant analysis
Mutation detection and proband reexamination
In this study, we performed WES and linkage analysis to identify a novel mutation in the EMD gene in a large Chinese DCM family. The WES ruled out the mutations in the exons of all reported nonsyndromic DCM-causing genes and found a single-nucleotide variant in GPR50 gene which was completely co-segregated with the disease phenotype. But the following functional studies did not reveal any functional change between the wild-type and variant GPR50 protein. These results pushed us to screen the sequencing failure exons in the WES data and eventually the disease-causing mutation in EMD was found. This process highlighted the indispensability of functional study in identifying new disease-causing genes and the importance of analyzing the sequencing failure data of WES.
Most reported mutations in the EMD gene caused apparent skeletal muscle disorders with cardiomyopathy in the form of EDMD [1, 2]. A novel deletion mutation in EMD exon 1 which results in almost a complete loss of emerin protein was identified in the pedigree. However, different from the typical EDMD clinical features caused by most reported EMD mutations, the patients in our study did not present apparent skeletal muscle atrophy and contractures. To date, only three EMD mutations (an in-frame 3-bp deletion in exon 2, a 5-bp deletion and a nonsense mutation in exon 6) were reported in patients with predominant cardiac diseases and mild skeletal muscle disorder [7–9]. These cases together with our findings pointed out that environmental or genetic modifications, such as functional overlapping proteins complement, may contribute to the observed clinical variability caused by EMD mutations.
Emerin is a member of the LEM-domain protein family, which presents a group of inner nuclear membrane and intranuclear proteins such as lamina-associated polypeptide 2 (LAP2) and MAN antigen 1 (MAN1) . Key components of the Rb pathway, which plays a crucial role in E2F-mediated cell cycle regulation and MyoD-mediated induction of myogenesis, were specifically altered in human muscle biopsies of EDMD patients . Emerin-knockout mice also revealed that Rb1 and MyoD expression levels and pathways during muscle regeneration were disrupted even though no overt abnormality was observed . Consistently, the protein levels of MyoD decreased in the complete and muscle-specific knockout mice of LAP2α, which is a mammal-specific non-membrane-associated isoform of the LAP2 gene . A mutation in the C-terminal tail of LAP2α has also been linked to DCM . These findings provided evidences that emerin and LAP2α may have some overlapping functions in the process of muscle regeneration and differentiation. Considering common cardiac disorder but diverse skeletal muscle involvement observed among the patients carrying EMD mutations, the functional complement from LAP2α and environmental effect may partly contribute to skeletal muscle regeneration.
We applied WES and linkage analysis to identify a novel deletion mutation in EMD exon 1 (c.26_39delATACCGAGCTGACC) in a Chinese family spanning five generations. Different from the typical clinical features caused by the majority of reported EMD mutations, the patients in our study present severe DCM but very mild skeletal muscle disorder. Our study suggests the EMD gene, which was not listed in the routine screening of DCM, should be considered in patients with DCM of unknown etiology and further functional studies of emerin in skeletal and cardiac muscle cells are needed to illustrate the phenotypic heterogeneity caused by EMD mutations.
Whole exome sequencing
Emery-Dreifuss muscular dystrophy
Leiden open variation database
Multiplex ligation-dependent probe amplification.
We wish to thank all of the individuals in the family for kindly agreement to participate in this study and Professor Jiada Li1 for his generous help in predicting the three-dimensional structure of GPR50 protein.
This work was supported in part by National Natural Science Foundation of China (Grant number: 81071028, 81172513), Program for New Century Excellent Talents in University (Grant number: NCET-10-0832) and National Key Basic Research Program (973 program) (Grant number: 2012CB944601).
- Wagner N, Krohne G: LEM-domain proteins: new insights into lamin-interacting proteins. Int Rev Cytol. 2007, 261: 1-46.View ArticlePubMedGoogle Scholar
- Dunnen JD: Leiden Muscular Dystrophy pages of Leiden Open Variation Database. [http://www.dmd.nl/nmdb2/home.php?select_db=EMD]
- Hershberger RE, Siegfried JD: Update 2011: clinical and genetic issues in familial dilated cardiomyopathy. J Am Coll Cardiol. 2011, 2011 (57): 1641-1649.View ArticleGoogle Scholar
- Garcia-Pavia P, Cobo-Marcos M, Guzzo-Merello G, Gomez-Bueno1 M, Bornstein B, Lara-Pezzi E, Segovia1 J, Alonso-Pulpon L: Genetics in dilated cardiomyopathy. Biomarkers Med. 2013, 74: 517-533.View ArticleGoogle Scholar
- Kaplan JC, Hamroun D: Gene Table of Neuromuscular Disorders. [http://www.musclegenetable.fr/]
- Levoye A, Dam J, Ayoub MA, Guillaume JL, Couturier C, Delagrange P, Jockers R: The orphan GPR50 receptor specifically inhibits MT1 melatonin receptor function through heterodimerization. EMBO J. 2006, 25: 3012-3023. 10.1038/sj.emboj.7601193.View ArticlePubMedPubMed CentralGoogle Scholar
- Vohanka S, Vytopil M, Bednarik J, Lukas Z, Kadankaa Z, Schildberger J, Ricotti R, Bioned S, Toniolod D: A mutation in the X-linked Emery-Dreifuss muscular dystrophy gene in a patient affected with conduction cardiomyopathy. Neuromuscul Disord. 2001, 11: 411-413. 10.1016/S0960-8966(00)00206-6.View ArticlePubMedGoogle Scholar
- Karst ML, Herron KJ, Olson TM: X-linked nonsyndromic sinus node dysfunction and atrial fibrillation caused by emerin mutation. J Cardiovasc Electrophysiol. 2008, 195: 510-515.View ArticleGoogle Scholar
- Sakata K, Shimizu M, Ino H, Yamaguchi M, Terai H, Fujino N, Hayashi K, Kaneda T, Inoue M, Oda Y, Fujita T, Kaku B, Kanaya H, Mabuchi H: High incidence of sudden cardiac death with conduction disturbances and atrial cardiomyopathy caused by a nonsense mutation in the STA gene. Circulation. 2005, 111: 3352-3358. 10.1161/CIRCULATIONAHA.104.527184.View ArticlePubMedGoogle Scholar
- Bakay M, Wang Z, Melcon G, Schiltz L, Xuan J, Zhao P, Sartorelli V, Seo J, Pegoraro E, Angelini C, Shneiderman B, Escolar D, Chen YW, Winokur ST, Pachman LM, Fan C, Mandler R, Nevo Y, Gordon E, Zhu Y, Dong Y, Wang Y, Hoffman EP: Nuclear envelope dystrophies show a transcriptional fingerprint suggesting disruption of Rb-MyoD pathways in muscle regeneration. Brain. 2006, 129: 996-1013. 10.1093/brain/awl023.View ArticlePubMedGoogle Scholar
- Melcon G, Kozlov S, Cutler DA, Sullivan T, Hernandez L, Zhao P, Mitchell S, Nader G, Bakay M, Rottman JN, Hoffman EP, Stewart CL: Loss of emerin at the nuclear envelope disrupts the Rb1/E2F and MyoD pathways during muscle regeneration. Hum Mol Genet. 2006, 15: 637-651. 10.1093/hmg/ddi479.View ArticlePubMedGoogle Scholar
- Gotic I, Schmidt WM, Biadasiewicz K, Leschnik M, Spilka R, Braun J, Stewart CL, Foisner R: Loss of LAP2a delays satellite cell differentiation and affects postnatal fiber-type determination. Stem Cells. 2010, 28: 480-488. 10.1002/stem.292.View ArticlePubMedGoogle Scholar
- Taylor MRG, Slavov D, Gajewski A, Vlcek S, Ku L, Fain PR, Carniel E, Lenarda AD, Sinagra G, Boucek MM, Cavanaugh J, Graw SL, Ruegg P, Feiger J, Zhu X, Ferguson DA, Bristow MR, Gotzmann J, Foisner R, Mestroni L, the Familial Cardiomyopathy Registry Research Group: Thymopoietin (lamina-associated polypeptide 2) gene mutation associated with dilated cardiomyopathy. Hum Mutat. 2005, 26: 566-574. 10.1002/humu.20250.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/15/77/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.