- Open Access
Novel compound heterozygous mutations in the MYO15A gene in autosomal recessive hearing loss identified by whole-exome sequencing
Journal of Translational Medicine volume 11, Article number: 284 (2013)
Inherited genetic defects play an important role in congenital hearing loss, contributing to about 60% of deafness occurring in infants. Hereditary nonsyndromic hearing loss is highly heterogeneous, and most patients with a presumed genetic etiology lack a specific molecular diagnosis.
By whole exome sequencing, we identified responsible gene of family 4794 with autosomal recessively nonsyndromic hearing loss (ARNSHL). We also used DNA from 56 Chinese familial patients with ARNSHL (autosomal recessive nonsyndromic hearing loss) and 108 ethnicity-matched negative samples to perform extended variants analysis.
We identified MYO15A c.IVS25 + 3G > A and c.8375 T > C (p.V2792A) as the disease-causing mutations. Both mutations co-segregated with hearing loss in family 4794, but were absent in the 56 index patients and 108 ethnicity-matched controls.
Our results demonstrated that the hearing loss of family 4794 was caused by novel compound heterozygous mutations in MYO15A.
Hearing loss is a common sensory defect that can significantly impact quality of life. The majority of congenital cases are attributable to genetic factors . Nonsyndromic hereditary forms, in which the hearing loss is the only clinical sign, are genetically heterogeneous. Autosomal recessive nonsyndromic hearing loss (ARNSHL) is the most common type and accounts for ~80% of cases of inherited hearing loss. To date, 61 genes and more than 100 genetic loci have been implicated in ARNSHL (see URL1). For many decades, linkage analysis has been the most powerful and widely used strategy to identify the gene defects responsible for inherited disorders. However, this approach is time consuming and requires the availability of cohorts of homogeneous and informative, possibly large families. As the molecular basis of deafness in most of the Chinese pedigrees is unsolved, we predict that many new deafness genes and mutations remain to be identified.
Whole-exome sequencing (WES) has become a highly efficient strategy for identifying novel causative genes and mutations involved in heritable disease . Both simple nonsyndromic and complex syndromic forms of hearing loss can be resolved efficiently using WES, especially in small families with distinct and interesting phenotypes that were once too small to map . To date, nine syndromic or nonsyndromic deafness genes have been identified using targeted genomic enrichment and next-generation sequencing: TPRN, GPSM2, CEACAM16, SMPX, HSD17B4, HARS2, MASP1, DNMT1, and TSPEAR[4–12]. As WES typically identifies many thousands of exomic variants, a selection strategy is important to facilitate the identification of the mutation that causes the disease [13, 14]. With ARSNHL, WES of several affected and healthy members from one family may be powerful for quickly identifying new susceptibility genes or mutations.
MYO15A has 66 exons and its coding protein, myosin XVa, is composed of 3,530 amino acids, and is critical for the formation of stereocilia in hair cells of the cochlea. In the organ of Corti, myosin XVa is localized exclusively at the tips of stereocilia and is a motor protein that uses energy from ATP hydrolysis to move along actin filaments . The tip of a stereocilium is the site of stereocilia growth and one of the proposed sites of mechano-electrical transduction . Myosin XVa-deficient mice were reported to lack any links between stereocilia, suggesting a complete disruption of the mechanotransduction machinery . Therefore, it was assumed that myosin XVa is required for proper formation and function of the mechanotransduction apparatus. Mutations in this gene cause DFNB3 hearing loss in individuals from different populations worldwide .
Here, we report a family with two siblings affected by severe-to-profound sensorineural hearing loss. Mutations in GJB2 and SLC26A4 were previously excluded. As the family is small and non-consanguineous, neither linkage analysis nor homozygosity mapping would have been informative for identifying the causative gene. Therefore, we used WES to identify the gene responsible for deafness in this family. WES was carried out in four members of family 4794 (two affected siblings and their parents), followed by validation. The results identified two novel compound heterozygous disease-segregating mutations, c.IVS25 + 3G > A and c.8375 T > C (p.V2792A), in the MYO15A gene. To exclude the possibility that these mutations were polymorphisms, DNA samples of 56 affected and 108 unaffected individuals were also analyzed. This is the first study to identify MYO15A mutations as the ARNSHL-associated gene in the Chinese deaf population.
Methods and materials
Family 4794 is a two-generation Chinese family with autosomal recessive prelingual non-syndromic sensorineural hearing loss. To screen for candidate mutations, we used 108 ethnicity-matched controls and 56 affected DNA samples from the Department of Otolaryngology, Head and Neck Surgery, Chinese PLA General Hospital. The 56 index patients were from families presenting with ARNSHL and in whom mutations of GJB2 and SLC26A4 had been previously excluded. Fully informed written consent was obtained from each subject. The study was approved by the Chinese PLA General Hospital Research Ethics Committee. Medical histories of the members of family 4794 were obtained using a questionnaire regarding the following aspects of the condition Otoscopy, physical examination, and pure tone audiometric examination (at frequencies from 250 to 8000 Hz) were performed to identify the phenotype. Immittance testing was applied to evaluate middle-ear pressure, ear canal volumes, and tympanic membrane mobility. Physical examination of all members revealed no signs of systemic illness or dysmorphic features. Computed tomography (CT) scan of the temporal bone was performed. The diagnosis of profound sensorineural hearing impairment was made according to the ICD-10 criteria based on audiometric examination. Tandem gait and Romberg tests were performed to evaluate balance.
Preparation of DNA
All genomic DNA was extracted from peripheral blood using a blood DNA extraction kit according to the protocol provided by the manufacturer (TianGen, Beijing, China).
Whole-exome capture and sequencing
Whole-exome capture and sequencing have been described in detail previously . In brief, genomic DNA was captured using the NimblegenSeqCap EZ Library (64 Mb for I1, I2, II:1, II:2). DNA was sheared, ligated to adaptors, extracted, amplified by ligation-mediated PCR (LM-PCR), and then hybridized to the Nimblegen SeqCap EZ Library for enrichment. The each captured library was loaded onto the Illumina Hiseq2000 platform. Raw image files were processed by Illumina base calling software 1.7.
Read mapping, variant detection, filtering, and annotation have been described in detail previously . Under the assumption of an autosomal-recessive pattern of inheritance, only variants that were homozygous or compound heterozygous in the two affected siblings and heterozygous in their parents were selected as candidates.
After filtering against multiple databases, Sanger sequencing was used to determine if any of the potential novel mutations in known causative genes of ARNSHL co-segregated with the disease phenotype in this family. Direct polymerase chain reaction (PCR) products were sequenced using Bigdye terminator v3.1 cycle sequencing kits (Applied Biosystems, Foster City, CA) and analyzed using an ABI 3700XL Genetic Analyzer.
Segregation of the mutations was studied in all family members. In addition, 108 negative samples and 56 ARNSHL families were also screened for the mutations by direct sequencing. Genotyping for c.IVS25 + 3G > A and c.8375 T > C was performed by PCR (Additional file 1: Table S1) and detected by bidirectional sequencing of the amplified fragments using an automated DNA sequencer (ABI 3100; Applied Biosystems). Nucleotide alteration(s) were identified by sequence alignment with the MYO15A GenBank sequence (NM_016239) using Genetool software.
Multiple sequence alignment
Multiple sequence alignment was performed according to a Homologene program with default settings and the sequences NP_057323.3 (H. sapiens), XP_0002888352.1 (M. mulatta), XP_003315632.1 (P. troglodytes), XP_002707808.1 (B. taurus), XP_536660.3 (C. lupus), NP_034992.2 (M. musculus), XP_577100.2 (R. norvegicus), XP_414818.3 (G. gallus), and XP_001924051 (D. rerio).
Splice site prediction
We used Fruitfly (see URL2) analysis to evaluate the effect of c.IVS25 + 3 G > A on splice site. A 2.0 kb genomic DNA fragment encompassing the region harboring the mutation was involved.
Model building and structural-based analysis
Three-dimensional (3D) modeling of the human wild type and p.V2792A mutation were performed using SWISS-MODEL, an automated homology modeling program (see URL3). This study used the automatic modeling approach to apply the complete protein sequence of human myosin XVa, including its 3519 amino acids and its mutation, which are available in the NCBI GenBank (NP_057323.3) in FASTA format. Data obtained by the homology models were visualized using Swiss-PdbViewer 4.1.
Clinical presentation of family 4794
Family 4794 includes two affected siblings (31 and 25 years old) and two unaffected parents (Figure 1A). Audiograms of the affected siblings showed that the hearing loss was bilateral and severe to profound (Figure 1B). The affected individuals did not have delays in gross motor development; neither did they have balance problems. Tandem walking was normal, and the Romberg test was negative. CT scan of the temporal bone in the proband excluded inner-ear malformations. The parents reported no history of stillbirth or miscarriage. Physical examination of all family members revealed no signs of systemic illness or dysmorphic features. The remaining examination results were normal. This phenotype is consistent with that reported previously for DFNB3 . No mutations were found in SLC26A4 and GJB2, and we proceeded to sequence the whole exomes of the two patients.
We applied whole-exome sequencing to the two affected siblings and two unaffected parents: I1, I2, II:1, and II:2. The reads were aligned with the human genome reference sequence. On average, we generated 8.71 billion bases of sequence with mean target region depths of 74.5. Approximately 99.2% (4805 mb in length) of the targeted bases were covered sufficiently to pass our thresholds for calling SNPs and short insertions or deletions (indels). After identification of variants, we focused only on non-synonymous variants, splice acceptor and donor site mutations, and frameshift coding indels, which were more likely to be pathogenic than others, especially those in homozygous or compound heterozygous mode (Additional file 2: Table S2).
On average, we identified 150,306 SNPs in the coding regions and 3,259 variants in introns that may affect splicing (within 10 bp of the intron/exon junction) for each sample (Additional file 3: Table S3). Then we compared these variants with the dbSNP135 (see URL4), HapMap project (see URL5), 1000 Genome Project (see URL6), and YH database (see URL7). We identified potential candidate variants by selecting the homozygous or compound heterozygous variants found in both affected siblings while unaffected parents contributed equally. Under the autosomal recessive mode of inheritance, we identified 67 genes with compound heterozygous variants and 16 genes with homozygous variants. Then we compared these variants to reported nonsyndromic hereditary hearing loss genes; 2 mutations (c.IVS25 + 3G > A, c.8375 T > C) were found in a previously reported deafness-related gene, MYO15A.
Mutation detection and analysis
According to the Berkeley Drosophila Genome Project (Fruitfly), a change in the splice donor sequence from GG to AG in intron 25 (c.IVS25 + 3 G > A) of MYO15A, detected in the present family, is predicted to result in an increase in splice site recognition by the splicing factor from 0.58 to 0.96, probably influencing the formation of the splice donor site.
The c.8375 T > C substitution occurs in exon 47 of MYO15A and results in a single amino acid substitution, valine to alanine (p.V2792A), within the FERMa domain. The V2792 amino acid is conserved across mammalian species, including human, chimpanzee, dog, mouse, and rat, as depicted in Figure 2. c.8375 T > C (p.V2792A) was predicted to be damaging by SIFT and Polyphen2 (see URL8,9).
Both mutations affect the highly preserved residues (Figure 2). Sanger sequencing validated the two affected siblings’ mutations and demonstrated that the two parents were heterozygous carriers of c.IVS25 + 3G > A (mother) and c.8375 T > C (father), showing complete co-segregation of the mutation with the phenotype (Figure 1C). The two mutations were absent both in 108 negative samples and in 56 ARNSHL families. These observations suggest that the mutation V2792A is likely to have a detrimental effect on the protein.
Structure modeling of p.V2792
A molecular model of myosin XVa was constructed based on the crystal structure of the talin head FERM domain (PDB ID:3ivf), a protein domain that binds F-actin in vitro and in vivo at a specific site within the actin filament. The constructed model covered the target sequence of myosin XVa (residues 2685–2896). Although the sequence identity between the target and the template was 14.16%, lower than the average 25%, the target area’s sequence similarity is high. Using the Swiss-Pdb Viewer 4.1, we show that the mutation is predicted to perturb the amino acid side chain because of the substitution of valine to alanine. This region of the protein is predicted to be highly hydrophobic, as previously shown in a hydrophobicity plot (Figure 3).
These data, together with the clinical presentation of the two affected siblings, clearly indicate that MYO15A mutations are the cause of the ARNSHL in this family.
Despite previous intensive linkage analysis and candidate gene screening, a large proportion of ARNSHL remain genetically unexplained. WES provides unprecedented opportunities to identify causative DNA mutations in rare heritable disorders. This report describes a Chinese family whose hearing impairment is caused by novel compound heterozygous mutations in MYO15A gene. Mutations in this gene have previously been described to cause autosomal recessive hearing loss.
Human myosin XVa consists of an N-terminal domain (amino acids 1–1,223), motor domain, neck region, and tail domain. The motor domain (amino acids 1,224–1,899) contains two binding sites for ATP and actin, the switches I and II helix, a relay helix, an SH3 helix, an SH1/SH2 helix, and a converter domain. The neck region contains two light-chain binding motifs (1,909–1,942). The tail region contains two MyTH4 (myosin tail homology 4) domains (2,066–2,174 and 3,051–3,161), two FERM domains (2,687–2,867 and 3,217–3,497) that have a role in binding of cytoskeletal proteins to cytoplasmic domains of transmembrane proteins, a putative SH3 domain (2,865–2,959), and a C-terminal class I PDZ-binding motif [20, 21]. The myosin XVa protein has been shown to be integral for development and elongation of the stereocilia through delivery of whirlin to the tips of the stereocilia. Whirlin binds to the SH3-MYTH4-FERM-domain-containing region of the myosin XVa protein and regulates actin filament elongation. There have been several reports of mutations in MYO15A causing hearing loss [17, 21–27]. Forty-seven mutations have previously been reported in MYO15A, and primarily occur in the motor domain.
The c.8375 T > C (p.V2792A) mutation is located within exon 47 and predicted to lie in the first FERM domain of myosin XVa. Three mutations in this domain (Q2716H, D2720H, D2823N) have been reported. The FERM domain is a widespread protein module involved in localizing protein to the plasma membrane, and its amino acid sequence is highly conserved .
In this family, a novel splice site mutation c.IVS25 + 3G > A was found in intron 25 of MYO15A and is predicted to change the highly conserved donor splice site of intron 25 from GG to AG. Mutations at splice sites are commonly associated with diminished amounts or abnormal maturation of mRNA, resulting in either exon skipping or cryptic splice site activation . These possibilities have not been investigated, and RNA studies are not possible in this instance because we cannot obtain fresh peripheral blood samples.
In summary, we have reported the clinical and genetic characteristics of a non-consanguineous Chinese family with ARNSHL by WES. To identify pathogenic variants, we consecutively filtered these variants by subjecting them to an analytical pipeline for high-confidence variant calling and annotation and identified novel compound heterozygous mutations in MYO15A. The identification of additional novel mutations in MYO15A further confirms the crucial role of MYO15A in auditory function. These results support sequence analysis of MYO15A in the clinical diagnostic testing of individuals with ARNSHL. Pre-implantation genetic diagnosis might be available for these families.
The URLs presented are as follows:
Hereditary Hearing Loss, http://hereditaryhearingloss.org/
HapMap project, ftp://ftp.ncbi.nlm.nih.gov/hapmap
1000 Genome Project, ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp
YH database, http://yh.genomics.org.cn/
XueGao, Qing-yan Zhu and Yue-Shuai Song are listed as co-first authors.
Autosomal recessive nonsyndromic hearing loss
Polymerase chain reaction.
Morton CC, Nance WE: Newborn hearing screening–a silent revolution. N Engl J Med. 2006, 354: 2151-2164. 10.1056/NEJMra050700.
Teer JK, Mullikin JC: Exome sequencing: the sweet spot before whole genomes. Hum Mol Genet. 2010, 19: R145-151. 10.1093/hmg/ddq333.
Diaz-Horta O, Duman D, Foster J, Sirmaci A, Gonzalez M, Mahdieh N, Fotouhi N, Bonyadi M, Cengiz FB, Menendez I: Whole-exome sequencing efficiently detects rare mutations in autosomal recessive nonsyndromic hearing loss. PLoS One. 2012, 7: e50628-10.1371/journal.pone.0050628.
Delmaghani S, Aghaie A, Michalski N, Bonnet C, Weil D, Petit C: Defect in the gene encoding the EAR/EPTP domain-containing protein TSPEAR causes DFNB98 profound deafness. Hum Mol Genet. 2012, 21: 3835-3844. 10.1093/hmg/dds212.
Li Y, Pohl E, Boulouiz R, Schraders M, Nurnberg G, Charif M, Admiraal RJ, von Ameln S, Baessmann I, Kandil M: Mutations in TPRN cause a progressive form of autosomal-recessive nonsyndromic hearing loss. Am J Hum Genet. 2010, 86: 479-484. 10.1016/j.ajhg.2010.02.003.
Walsh T, Shahin H, Elkan-Miller T, Lee MK, Thornton AM, Roeb W, Abu Rayyan A, Loulus S, Avraham KB, King MC, Kanaan M: Whole exome sequencing and homozygosity mapping identify mutation in the cell polarity protein GPSM2 as the cause of nonsyndromic hearing loss DFNB82. Am J Hum Genet. 2010, 87: 90-94. 10.1016/j.ajhg.2010.05.010.
Zheng J, Miller KK, Yang T, Hildebrand MS, Shearer AE, DeLuca AP, Scheetz TE, Drummond J, Scherer SE, Legan PK: Carcinoembryonic antigen-related cell adhesion molecule 16 interacts with alpha-tectorin and is mutated in autosomal dominant hearing loss (DFNA4). Proc Natl Acad Sci USA. 2011, 108: 4218-4223. 10.1073/pnas.1005842108.
Schraders M, Haas SA, Weegerink NJ, Oostrik J, Hu H, Hoefsloot LH, Kannan S, Huygen PL, Pennings RJ, Admiraal RJ: Next-generation sequencing identifies mutations of SMPX, which encodes the small muscle protein, X-linked, as a cause of progressive hearing impairment. Am J Hum Genet. 2011, 88: 628-634. 10.1016/j.ajhg.2011.04.012.
Pierce SB, Walsh T, Chisholm KM, Lee MK, Thornton AM, Fiumara A, Opitz JM, Levy-Lahad E, Klevit RE, King MC: Mutations in the DBP-deficiency protein HSD17B4 cause ovarian dysgenesis, hearing loss, and ataxia of Perrault Syndrome. Am J Hum Genet. 2010, 87: 282-288. 10.1016/j.ajhg.2010.07.007.
Pierce SB, Chisholm KM, Lynch ED, Lee MK, Walsh T, Opitz JM, Li W, Klevit RE, King MC: Mutations in mitochondrial histidyl tRNA synthetase HARS2 cause ovarian dysgenesis and sensorineural hearing loss of Perrault syndrome. Proc Natl Acad Sci USA. 2011, 108: 6543-6548. 10.1073/pnas.1103471108.
Sirmaci A, Walsh T, Akay H, Spiliopoulos M, Sakalar YB, Hasanefendioglu-Bayrak A, Duman D, Farooq A, King MC, Tekin M: MASP1 mutations in patients with facial, umbilical, coccygeal, and auditory findings of Carnevale, Malpuech, OSA, and Michels syndromes. Am J Hum Genet. 2010, 87: 679-686. 10.1016/j.ajhg.2010.09.018.
Klein CJ, Botuyan MV, Wu Y, Ward CJ, Nicholson GA, Hammans S, Hojo K, Yamanishi H, Karpf AR, Wallace DC: Mutations in DNMT1 cause hereditary sensory neuropathy with dementia and hearing loss. Nat Genet. 2011, 43: 595-600. 10.1038/ng.830.
van den Berg H, van der Gaag M, Hendriks H: Influence of lifestyle on vitamin bioavailability. Int J Vitam Nutr Res. 2002, 72: 53-59. 10.1024/0300-98184.108.40.206.
Sierksma A, van der Gaag MS, Kluft C, Hendriks HF: Effect of moderate alcohol consumption on fibrinogen levels in healthy volunteers is discordant with effects on C-reactive protein. Ann N Y Acad Sci. 2001, 936: 630-633.
Belyantseva IA, Boger ET, Friedman TB: Myosin XVa localizes to the tips of inner ear sensory cell stereocilia and is essential for staircase formation of the hair bundle. Proc Natl Acad Sci USA. 2003, 100: 13958-13963. 10.1073/pnas.2334417100.
Beyer LA, Odeh H, Probst FJ, Lambert EH, Dolan DF, Camper SA, Kohrman DC, Raphael Y: Hair cells in the inner ear of the pirouette and shaker 2 mutant mice. J Neurocytol. 2000, 29: 227-240. 10.1023/A:1026515619443.
Wang A, Liang Y, Fridell RA, Probst FJ, Wilcox ER, Touchman JW, Morton CC, Morell RJ, Noben-Trauth K, Camper SA, Friedman TB: Association of unconventional myosin MYO15 mutations with human nonsyndromic deafness DFNB3. Science. 1998, 280: 1447-1451. 10.1126/science.280.5368.1447.
Gao X, Su Y, Guan LP, Yuan YY, Huang SS, Lu Y, Wang GJ, Han MY, Yu F, Song YS: Novel Compound Heterozygous TMC1 Mutations Associated with Autosomal Recessive Hearing Loss in a Chinese Family. PLoS One. 2013, 8: e63026-10.1371/journal.pone.0063026.
Sirmaci A, Duman D, Ozturkmen-Akay H, Erbek S, Incesulu A, Ozturk-Hismi B, Arici ZS, Yuksel-Konuk EB, Tasir-Yilmaz S, Tokgoz-Yilmaz S: Mutations in TMC1 contribute significantly to nonsyndromic autosomal recessive sensorineural hearing loss: a report of five novel mutations. Int J Pediatr Otorhinolaryngol. 2009, 73: 699-705. 10.1016/j.ijporl.2009.01.005.
Garcia-Alvarez B, de Pereda JM, Calderwood DA, Ulmer TS, Critchley D, Campbell ID, Ginsberg MH, Liddington RC: Structural determinants of integrin recognition by talin. Mol Cell. 2003, 11: 49-58. 10.1016/S1097-2765(02)00823-7.
Kalay E, Uzumcu A, Krieger E, Caylan R, Uyguner O, Ulubil-Emiroglu M, Erdol H, Kayserili H, Hafiz G, Baserer N: MYO15A (DFNB3) mutations in Turkish hearing loss families and functional modeling of a novel motor domain mutation. Am J Med Genet A. 2007, 143A: 2382-2389. 10.1002/ajmg.a.31937.
Friedman TB, Sellers JR, Avraham KB: Unconventional myosins and the genetics of hearing loss. Am J Med Genet. 1999, 89: 147-157. 10.1002/(SICI)1096-8628(19990924)89:3<147::AID-AJMG5>3.0.CO;2-6.
Liburd N, Ghosh M, Riazuddin S, Naz S, Khan S, Ahmed Z, Liang Y, Menon PS, Smith T, Smith AC: Novel mutations of MYO15A associated with profound deafness in consanguineous families and moderately severe hearing loss in a patient with Smith-Magenis syndrome. Hum Genet. 2001, 109: 535-541. 10.1007/s004390100604.
Lezirovitz K, Pardono E, de Mello Auricchio MT, de Carvalho ESFL, Lopes JJ, Abreu-Silva RS, Romanos J, Batissoco AC, Mingroni-Netto RC: Unexpected genetic heterogeneity in a large consanguineous Brazilian pedigree presenting deafness. Eur J Hum Genet. 2008, 16: 89-96. 10.1038/sj.ejhg.5201917.
Belguith H, Aifa-Hmani M, Dhouib H, Said MB, Mosrati MA, Lahmar I, Moalla J, Charfeddine I, Driss N, Arab SB: Screening of the DFNB3 locus: identification of three novel mutations of MYO15A associated with hearing loss and further suggestion for two distinctive genes on this locus. Genet Test Mol Biomarkers. 2009, 13: 147-151. 10.1089/gtmb.2008.0077.
Shearer AE, Hildebrand MS, Webster JA, Kahrizi K, Meyer NC, Jalalvand K, Arzhanginy S, Kimberling WJ, Stephan D, Bahlo M: Mutations in the first MyTH4 domain of MYO15A are a common cause of DFNB3 hearing loss. Laryngoscope. 2009, 119: 727-733. 10.1002/lary.20116.
Cengiz FB, Duman D, Sirmaci A, Tokgoz-Yilmaz S, Erbek S, Ozturkmen-Akay H, Incesulu A, Edwards YJ, Ozdag H, Liu XZ, Tekin M: Recurrent and private MYO15A mutations are associated with deafness in the Turkish population. Genet Test Mol Biomarkers. 2010, 14: 543-550. 10.1089/gtmb.2010.0039.
Chishti AH, Kim AC, Marfatia SM, Lutchman M, Hanspal M, Jindal H, Liu SC, Low PS, Rouleau GA, Mohandas N: The FERM domain: a unique module involved in the linkage of cytoplasmic proteins to the membrane. Trends Biochem Sci. 1998, 23: 281-282. 10.1016/S0968-0004(98)01237-7.
Krawczak M, Reiss J, Cooper DN: The mutational spectrum of single base-pair substitutions in mRNA splice junctions of human genes: causes and consequences. Hum Genet. 1992, 90: 41-54.
This work was supported by grants from the Project of the National Natural Science Foundation of China (Nos. 30801285, 31071099, 81230020, 81200751, 81070792, 81000414, 81000415), State 863 High Technology R&D Key Project of China (No.2011AA02A112), China Postdoctoral Science Foundation (No.2012 M521878) a grant from Minister of Science and Technology of China (No.2012BAI09B02), and a grant from Minister of Health of China (No.201202005). We sincerely thank all the family members for their participation and cooperation in this study. We are also grateful to Dr. Guangchun Bai of Albany Medical College, New York State Department of Health for his valuable comments and suggestions on protein structure analysis.
The authors declare that they have no competing interests.
Conceived and designed the experiments: PD XG. Performed the experiments: XG QYZ YSS GJW FX DYK. Analyzed the data: XG YYY SSH MYH LPG. Contributed reagents/materials/analysis tools: JGZ. Wrote the paper: XG PD. All authors read and approved the final manuscript.
Xue Gao, Qing-yan Zhu, Yue-Shuai Song contributed equally to this work.