Novel MYO15A variants are associated with hearing loss in the two Iranian pedigrees
BMC Medical Genetics volume 21, Article number: 226 (2020)
Clinical genetic diagnosis of non-syndromic hearing loss (NSHL) is quite challenging. With regard to its high heterogeneity as well as large size of some genes, it is also really difficult to detect causative mutations using traditional approaches. One of the recent technologies called whole-exome sequencing (WES) has been thus developed in this domain to remove the limitations of conventional methods.
This study was a report on a research study of two unrelated pedigrees with multiple affected cases of hearing loss (HL). Accordingly, clinical evaluations and genetic analysis were performed in both families.
The results of WES data analysis to uncover autosomal recessive non-syndromic hearing loss (ARNSHL) disease-causing variants was reported in the present study. Initial analysis identified two novel variants of MYO15A i.e. c.T6442A:p.W2148R and c.10504dupT:p.C3502Lfs*15 correspondingly which were later confirmed by Sanger validations and segregation analyses. According to online prediction tools, both identified variants seemed to have damaging effects.
In this study, whole exome sequencing were used as a first approach strategy to identify the two novel variants in MYO15A in two Iranian families with ARNSHL.
Hearing impairment is considered as an etiologically heterogeneous sensory deficiency with incidence 1 in 1000 newborns around the world . In this regard, genetic hearing loss (HL) has been divided into syndromic and non-syndromic types. Considering the high rate of consanguineous marriages in the Middle East, autosomal recessive non-syndromic hearing loss (ARNSHL) is reportedly more prevalent western countries . However; due to the wide variety of pathogenic genes associated with non-syndromic hearing loss (NSHL), including both nuclear and mitochondrial ones, the disease includes diverse patterns of inheritance comprised of autosomal dominant, autosomal recessive, mitochondrial, and X-linked recessive. Disease heterogeneity has been admittedly recognized as the most important challenge in genetic diagnosis of NSHL. Diagnostic approaches which have been relied on conventional methods based on genetic testing of the most common genes, often fail to determine the exact genetic cause of the disease in many countries including Iran . In heterogeneous populations like Iran, the distribution of mutations in the gap junction beta-2 protein, also known as connexin 26, (i.e. GJB2) gene as a major cause of ARNSHL can be extremely diverse depending on patients’ ethnicities. The prevalence rate of GJB2-related hearing loss has been reported by 38.3% in northern Iran, but the percentage of such variations has been found very rarely in southern regions . Various frequencies of causative mutations, compound heterozygotes, as well as nuclear modifier genes can also render the molecular diagnosis of ARNSHL as a challenge [3, 4]. Moreover, mutation screening in large genes such as MYO15A (66 exons) can be impeded following the use of traditional approaches based on Sanger sequencing [5, 6]. Today, molecular genetic testing on the basis of multi-gene screening such as whole-exome sequencing (WES) are being used instead of traditional diagnostic procedures .
For first time, MYO15A (DFNB3 locus) mutation was reported from Indonesia population. Until now, there have been many reports of mutations of MYO15A causing ARNSHL in different countries of Asia such as Pakistan , Turkey  and Iran . In the present study, two novel MYO15A variants identified by WES from two Iranian families with ARNSHL is reported.
Subjects and clinical evaluations
The study was approved by the Ethics Committee of Shahrekord University of Medical Sciences (IR.SKUMS.REC.1397.008), Iran. Two Iranian families from Hamedan Province with hearing impairments, without any other additional symptoms were thus studied. Informed written consent was taken from both families. The proband from each family was further subjected to clinical evaluations of the inner ear accompanied by pure-tone audiometry (PTA).
WES was used to detect the deafness associated variants in the DNA sample in probands. Genomic deoxyribonucleic acid (DNA) was extracted from whole peripheral blood of each study subject utilizing DNA Extraction Kit DNP (Sinacolon, Iran) according to the manufacturer’s instructions. Purity and concentration of DNA samples were further measured via Thermo Scientific NanoDrop 2000c Spectrophotometer.
DNA samples from each pedigree’s proband (Fig. 1a, V-4 in family 1; Fig. 1c, II-3 in family 2) were then subjected to WES at Macrogen Online Sequencing Order System (Seoul, South Korea) on Genome Analyzer/HiSeq 2000 (Illumina, San Diego, CA, USA, 151-bp paired-end reads). It should be noted that the library had been prepared through SureSelect XT Library Prep Kit (Agilent Technologies, CA, USA). Data analysis was correspondingly performed using an in-house developed pipeline, adopted from Genome Analysis Tool Kit v3.6 and ANNOVAR software . Homozygous missense, start codon change, splice site, nonsense, stop loss, and indel variants with minor allele frequency < 1% were further filtered in dbSNP (version 138), 1000 Genomes Project, Exome Aggregation Consortium (ExAC), and NHLBI GO Exome Sequencing Project (ESP). Based on autosomal recessive inheritance, the homozygosity region of samples was determined using homozygosity mapping algorithms.
In order to prioritize the candidate functional variants, several online prediction software including MutationTaster2, FATHMM, PANTHER, SIFT, PROVEAN, MetaLR, PolyPhen-2, CADD, and ConSurf were also used to evaluate the pathogenic effects of the variants with a frequency less than 0.01. After observing the autosomal recessive inheritance pattern in both pedigrees, the homozygous variants were then prioritized. Next, the variants were investigated in the Human Gene Mutation Database, hereditaryhearingloss home page and the related literature to survey their association with a phenotype and novelty of the variants.
Besides, candidate variant segregation from exome data was evaluated through polymerase chain reaction (PCR)-based Sanger sequencing. Therefore, the following primers were synthesized: 5ʹ-GAACTACATCGTGCAGAAGG-3ʹ and 5ʹ-CCTATCCAGTCCCACTCACT-3ʹ for human MYO15A c.T6442A variant and 5ʹ-CCACCATTCGGCCTTCCA-3ʹ and 5ʹ-CTGCCTCCTCTTAGTGTCCTC-3ʹ for human MYO15A c.10504dupT variant.
Clinical and molecular findings
The first family pedigree is displayed in Fig. 1a. Accordingly, four members of the pedigree including two affected and two unaffected individuals who consented to be included in this study are indicated. In this family, the proband (V: 4) was a 21-year-old woman with congenital HL. No additional abnormal phenotypic features including visual impairments or any limb and facial malformations were observed in the proband. The parents were consanguineous and both showed normal hearing. According to the audiogram, the proband is suffering from congenital profound deafness (Fig. 1b).
A three-generation pedigree, depicted in Fig. 1c, was presented as the second family with ten members, six males and fourteen females, suffering from ARNSHL. The proband was a 25-year-old male individual, born as the second child of non-sanguineous healthy parents, who had been diagnosed for congenital HL when he was 1 year old. At the age of 5, he had gone through bilateral cochlear implant surgery based on physical examinations and audiometry testing. Audiogram analysis also confirmed HL to be profound in the proband (II: 3) (Fig. 1d).
Considering the all limitations in genetic diagnosis of ARNSHL, it was decided to perform WES as a first approach on proband’s genome DNA. The total number of bases, reads, GC (%), and Q30 (%) are calculated for the subjected sample. 70,760,838 reads were produced, and total read bases were 10.7G bp. The GC content (%) and Q30 were 51.87 and 95.69% respectively. Following the filtering step depicted in Fig. 2, two novel homozygous variants in MYO15A i.e. c.T6442A:p.W2148R and c.10504dupT:p.C3502Lfs*15 were prioritized in family I and II, respectively. Given that the proband of family I was the offspring of consanguineous parents, WES dataset revealed that MYO15A (c.T6442A) variants had resided in the large homozygous regions on chromosome 17 (Fig. 3). All in silico programs also predicted damaging effects of p.W2148R variant. Mutation taster further suggested that c.10504dupT variant had deleterious effects. Moreover, the first conservative amino acid alternation (p.C3502L) predicted to be pathogenic using online softwares (Table 1). Analysis of genotype-phenotype correlations revealed that patients with c.T6442A:p.W2148R and c.10504dupT:p.C3502Lfs*15 tended to have profound hearing loss.
Additionally, the results of Sanger sequencing confirmed the presence of p.W2148R variant in MYO15A gene in the proband and other affected members who were studied, but the unaffected sister (V.5) and her father were found heterozygous for the variant (Fig. 1e).
The findings of Sanger sequencing also revealed co-segregation of c.10504dupT variant in the second family (Fig. 1f). The affected proband was thus homozygous, whereas his unaffected sister (II:5) and his father (I:1) were heterozygous for this locus and his unaffected brother (II:2) lacked the variant.
It should be noted that these two variants were absent in 50 ethnically-matched control cases.
Based on WES data as well as segregation and genotype-phenotype correlational studies, mutations in MYO15A gene were identified as a main contributor of NSHL in the first and second surveyed families.
MYO15A, as a new branch of the myosin protein-coding gene superfamily, has a role in stereocilia formation . This gene is considered as the third leading cause of ARNSHL in many populations [12, 13], including Iran, with a prevalence rate ranged from 4.8–9.6% [14,15,16] . Moreover, mutations in this gene have been associated with severe-to-profound HL. Screening of 66 coding exons through Sanger sequencing is expensive and more time-consuming. On the contrary, high-throughput techniques can save time and money .
MYO15A is a different form of myosins protein with long N-terminal extensions following by the conserved motor domain, IQ motifs (calmodulin/ myosin light chain binding), MyTh4 domains (Myosin-Tail like Homology region 4), FERM motifs (4.1 protein, Ezrin, Radixin, and Moesin), SH3 domain (Src Homology 3), and the PDZ ligand domain (Post synaptic density protein (PSD95), Dlg1 (Drosophila disc large tumor suppressor), and zo-1 (Zonula occludens-1 protein) . The identified c.T6442A: p.W2148R variant in our study is located in the first (MyTH4) domain and is also conserved among different species (data not shown). The association of (MyTH4) domain mutations with hearing loss was firstly reported in 1998 . Due to documents, variants in this domain have related severe to profound hearing loss, which is consistent with our pedigree phenotype [5, 15, 16, 20,21,22]. To the best of our knowledge, p.P2073S, p.R2124Q  p.P2073L, p.V2114G , R2146Q , p.R2146W  and c.6273 + 1G > A  variants have been detected in Iranian population until now [15, 16, 22]. The p.R2124Q and p.P2073S were the first reported mutations in the MyTH4 domain of MYO15A protein which was located in conserved fourth helix of MyTH4. Variations in the MyTH4 domains interfere in forming of transmembrane actin microfilament assembly complex at the stereocilia tips .
Mehregan et al., reported p.Arg2146Gln in the fourth helix of the first MyTH4 core, which results in severe to profound hearing loss. Structural analysis of this variation has revealed that this substitution alters binding properties at the domain surface .
The substitution of a highly conserved amino acid (Table 1), hydrophobic non-polar tryptophan, with arginine can lead to a loss of hydrophobic pocket. The counterparts of p.W2148 in the more studied family members, i.e. MYO7A is W1192. Sans CEN2 (a scaffold protein) interacted with MyTH4 domain by extensive hydrogen bonding, hydrophobic contacts, and charge-charge interactions. The hydrophobic pocket in MYO7A comprised the conserved A1189, W1192, I1193, P1220, and Y1223 residues. Substitution of Ala1189 by Glu, leads to a.
~ 10-fold decrease of the binding affinity between MyTH4 domain and CEN. The complex of myosins and SANS linked cadherins to the actin cytoskeleton . Woo et al., explained that p.R2146Q in myosin 15A and R1190 of MYO7A had similar structure in MyTH4- FERM domains and interfere to binding CEN2 to this domain .
In addition, Myosin interacts with other scaffolding proteins (whirlin and Eps8 (Epidermal Growth Factor Receptor Pathway Substrate 8)) and can transport them to the tip of stereocilia to form a stereocilia tip complex, which can facilitate maturation of stereocilia. These scaffold proteins are essential for normal hearing in humans [27, 28]. MyTH4 domain contains the actin-binding sites. The overall surface of the microtubule (MT) is negatively charged. The positively charged motifs with surface-exposed hydrophobic side chains found in the myosin MyTH4 domain can serve as an MT binding site .
We speculate the, c.T6442A:p.W2148R variant interferes with the formation of the Myosin 15A-whirlin-Eps8-CEN2 complexs and microtubule-binding.
In the second studied family, WES could successfully detect a novel homozygous insertion variant i.e. c.10504dupT:p.C3502Lfs*15 in MYO15A gene, co-segregated with the disease within the pedigree. This variant between the second FERM and PDZ domains could also lead to a reading frame shift at position 10,504 and a stop codon (p.C3502Lfs*15) with truncation and translation of mRNA resulting in lack of its conserved amino acids at C-terminal (data not shown). Lezirovitz et al. reported a frameshift mutation c.10573delA to cause profound.
hearing loss. This mutation in the PDZ binding ligand of MYO15A altered interaction of this protein with whirlin. Thus, stereocilia elongation did not occurred .
Zhang et al. identified p.Leu3501Glu variant was associated with profound hearing loss .
Today, the topic of oligogenic inheritance traits has made a big wave in diagnostic medicine; since, in many monogenic diseases, it represents that there is not just one gene affected phenotype which causes new challenges in diagnosing these diseases and it can be more complicated for diseases with heterogenic pathology in which many genes are involved .
In summary, we identified two novel variants (p.W2148R and p.C3502Lfs*15), in the MYO15A gene in two Iranian families using whole-exome sequencing. Accordingly, we showed that these likely pathogenic variants were segregated with the profound hearing loss in both families. However, further functional analysis is required to confirm the results of the present study.
Availability of data and materials
The raw datasets generated and/or analyzed during the current study are not publicly available because it is possible that individual privacy could be compromised, but are available from the corresponding author on reasonable request.
Exome Aggregation Consortium
Exome Sequencing Project
Epidermal Growth Factor Receptor Pathway Substrate 8
4.1 protein, Ezrin, Radixin, and Moesin
Myosin-Tail like Homology region 4
Non-syndromic hearing loss
Src Homology 3
Whole exome sequencing
Zonula occludens-1 protein
Motavaf M, Soveizi M, Maleki M, Mahdieh N. MYO15A splicing mutations in hearing loss: a review literature and report of a novel mutation. Int J Pediatr Otorhinolaryngol. 2017;96:35–8.
Najmabadi H, Kahrizi K. Genetics of non-syndromic hearing loss in the Middle East. Int J Pediatr Otorhinolaryngol. 2014;78(12):2026–36.
Beheshtian M, Babanejad M, Azaiez H, Bazazzadegan N, Kolbe D, Sloan-Heggen C, et al. Heterogeneity of hereditary hearing loss in Iran: a comprehensive review. Arch Iran Med. 2016;19(10):720.
Khatami S, Rokni-Zadeh H, Mohsen-Pour N, Biglari A, Changi-Ashtiani M, Shahrooei M, et al. Whole exome sequencing identifies both nuclear and mitochondrial variations in an Iranian family with non-syndromic hearing loss. Mitochondrion. 2019;46:321–5.
Woo H-M, Park H-J, Baek J-I, Park M-H, Kim U-K, Sagong B, et al. Whole-exome sequencing identifies MYO15A mutations as a cause of autosomal recessive nonsyndromic hearing loss in Korean families. BMC Med Genet. 2013;14(1):72.
Asgharzade S, Chaleshtori MH, Tabatabaifar MA, Reisi S, Modaressi MH. Mutation in second exon of MYO15A gene cause of nonsyndromic hearing loss and its association in the Arab population in Iran. Genetika. 2016;48(2):587–96.
Seco CZ, Wesdorp M, Feenstra I, Pfundt R, Hehir-Kwa JY, Lelieveld SH, et al. The diagnostic yield of whole-exome sequencing targeting a gene panel for hearing impairment in the Netherlands. Eur J Hum Genet. 2017;25(3):308.
Liburd N, Ghosh M, Riazuddin S, Naz S, Khan S, Ahmed Z, et al. 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(5):535–41.
Kalay E, Uzumcu A, Krieger E, Çaylan R, Uyguner O, Ulubil-Emiroglu M, et al. MYO15A (DFNB3) mutations in Turkish hearing loss families and functional modeling of a novel motor domain mutation. Am J Med Genet A. 2007;143(20):2382–9.
Pourreza MR, Mohammadi H, Sadeghian L, Asgharzadeh S, Sehhati M, Tabatabaiefar MA. Applying two different Bioinformatic approaches to discover novel genes associated with hereditary hearing loss via whole-exome sequencing: ENDEAVOUR and HomozygosityMapper. Adv Biomed Res. 2018;7:141.
Rędowicz MJ. Myosins and deafness. J Muscle Res Cell Motil. 1999;20(3):241–8.
Rehman AU, Bird JE, Faridi R, Shahzad M, Shah S, Lee K, et al. Mutational spectrum of MYO15A and the molecular mechanisms of DFNB3 human deafness. Hum Mutat. 2016;37(10):991–1003.
Hilgert N, Smith RJ, Van Camp G. Forty-six genes causing nonsyndromic hearing impairment: which ones should be analyzed in DNA diagnostics? Mutat Res. 2009;681(2):189–96.
Fattahi Z, Shearer AE, Babanejad M, Bazazzadegan N, Almadani SN, Nikzat N, et al. Screening for MYO15A gene mutations in autosomal recessive nonsyndromic, GJB2 negative Iranian deaf population. Am J Med Genet A. 2012;158(8):1857–64.
Babanejad M, Fattahi Z, Bazazzadegan N, Nishimura C, Meyer N, Nikzat N, et al. A comprehensive study to determine heterogeneity of autosomal recessive nonsyndromic hearing loss in Iran. Am J Med Genet A. 2012;158(10):2485–92.
Sloan-Heggen CM, Babanejad M, Beheshtian M, Simpson AC, Booth KT, Ardalani F, et al. Characterising the spectrum of autosomal recessive hereditary hearing loss in Iran. J Med Genet. 2015;52(12):823–9.
Zhang J, Guan J, Wang H, Yin L, Wang D, Zhao L, et al. Genotype-phenotype correlation analysis of MYO15A variants in autosomal recessive non-syndromic hearing loss. BMC Med Genet. 2019;20(1):60.
Anderson DW, Probst FJ, Belyantseva IA, Fridell RA, Beyer L, Martin DM, et al. The motor and tail regions of myosin XV are critical for normal structure and function of auditory and vestibular hair cells. Hum Mol Genet. 2000;9(12):1729–38.
Wang A, Liang Y, Fridell RA, Probst FJ, Wilcox ER, Touchman JW, et al. Association of unconventional myosin MYO15 mutations with human nonsyndromic deafness DFNB3. Science. 1998;280(5368):1447–51.
Miyagawa M, Nishio S-Y, Hattori M, Moteki H, Kobayashi Y, Sato H, et al. Mutations in the MYO15A Gene Are a Significant Cause of Nonsyndromic Hearing Loss: Massively Parallel DNA Sequencing–Based Analysis. Ann Otol Rhinol Laryngol. 2015;124(1_suppl):158S–68S.
Yang T, Wei X, Chai Y, Li L, Wu H. Genetic etiology study of the non-syndromic deafness in Chinese Hans by targeted next-generation sequencing. Orphanet J Rare Dis. 2013;8(1):85.
Shearer AE, Hildebrand MS, Webster JA, Kahrizi K, Meyer NC, Jalalvand K, et al. Mutations in the first MyTH4 domain of MYO15A are a common cause of DFNB3 hearing loss. Laryngoscope. 2009;119(4):727–33.
Mehregan H, Mohseni M, Jalalvand K, Arzhangi S, Nikzat N, Banihashemi S, et al. Novel mutations in MYTH4-FERM domains of myosin 15 are associated with autosomal recessive nonsyndromic hearing loss. Int J Pediatr Otorhinolaryngol. 2019;117:115–26.
Mehregan H, Mohseni M, Jalalvand K, Arzhangi S, Nikzat N, Banihashemi S, et al. Novel mutations in MYTH4-FERM domains of myosin 15 are associated with autosomal recessive nonsyndromic hearing loss. Int J Pediatr Otorhinolaryngol. 2018;117:115–26.
Yan D, Tekin D, Bademci G, Foster J, Cengiz FB, Kannan-Sundhari A, et al. Spectrum of DNA variants for non-syndromic deafness in a large cohort from multiple continents. Hum Genet. 2016;135(8):953–61.
Wu L, Pan L, Wei Z, Zhang M. Structure of MyTH4-FERM domains in myosin VIIa tail bound to cargo. Science. 2011;331(6018):757–60.
Belyantseva IA, Boger ET, Naz S, Frolenkov GI, Sellers JR, Ahmed ZM, et al. Myosin-XVa is required for tip localization of whirlin and differential elongation of hair-cell stereocilia. Nat Cell Biol. 2005;7(2):148.
Behlouli A, Bonnet C, Abdi S, Bouaita A, Lelli A, Hardelin J-P, et al. EPS8, encoding an actin-binding protein of cochlear hair cell stereocilia, is a new causal gene for autosomal recessive profound deafness. Orphanet J Rare Dis. 2014;9(1):55.
Planelles-Herrero VJ, Blanc F, Sirigu S, Sirkia H, Clause J, Sourigues Y, et al. Myosin MyTH4-FERM structures highlight important principles of convergent evolution. Proc Natl Acad Sci. 2016;113(21):E2906–E15.
Lezirovitz K, Pardono E, de Mello Auricchio MT, FLdC ES, Lopes JJ, Abreu-Silva RS, et al. Unexpected genetic heterogeneity in a large consanguineous Brazilian pedigree presenting deafness. Eur J Hum Genet. 2008;16(1):89–96.
Badano JL, Katsanis N. Beyond Mendel: an evolving view of human genetic disease transmission. Nat Rev Genet. 2002;3(10):779–89.
We sincerely thank the patients and their families for participating in this study.
This work was financially supported by Sharekord University of Medical Sciences. Whole exome sequencing and its data collection and data analysis were supported by Sharekord University of Medical Sciences (grant No. 2685).
Ethics approval and consent to participate
The study was approved by the Review Board of Sharekord University of Medical Sciences (IR.SKUMS.REC.1397.008). Written informed consent was obtained from all of the participants in the study and a written consent to participate was obtained from the parents of the patient (younger than the age of 16).
Consent for publication
Written informed consent for publication of clinical details and clinical images was obtained from the all of the participants. The parents /legal guardian of the children (under the age of 18) has signed written informed consent for publication of clinical and genetic data.
All authors declare no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Khatami, S., Askari, M., Bahreini, F. et al. Novel MYO15A variants are associated with hearing loss in the two Iranian pedigrees. BMC Med Genet 21, 226 (2020). https://doi.org/10.1186/s12881-020-01168-x
- Whole exome sequencing
- First approach