Bioinformatic and Molecular Analysis
Identified variants in this study include missense variants (27), insertion/deletion (23) , nonsense (6), synonymous (1), and non-coding variants (6) (
Figure 1. E). Ten known NSHL genes were identified in 27 of the 43 families (
Table 1). These genes included
GJB2,
SLC26A4,
MYO7A,
DIAPH3,
PTPRQ,
LOXHD1,
CDH23,
MPZL2,
EYA4, and
OTOA. Five known SHL genes were identified in 6 of the 43 families (
Table 2). These genes included
TCOF1,
EYA1,
FDXR,
AIFM1, and
SOX10.
ATP7B was segregate in SHL family as a candidate SHL gene (
Table 3). Variants in
GJB2 (14 families, 32.6%) and
SLC26A4 (5 families, 11.6%) accounted for the majority (44.2%).
EYA1 variants were found in two families, and variants in 13 other genes were only observed in one family. Twenty of the variants have been reported to be associated with the occurrence of hearing loss, while thirteen of the novel variants were found to segregate in 11/43 HL families, inlcude
DIAPH3 c.2256_2257insT,
PTPRQ c.6293T>C,
CDH23 c.4859T>A,
MPZL2 c.393_436+21del,
TCOF1 c.3997_4007del,
LOXHD1 c.2438T>A and c.1759C>T,
EYA1 c.1350_1353delTAATinsCAGACA,
FDXR c.1069G>T and c.364C>T, ATP7B c.4014T>A,
SOX10 c.133del,
AIFM1 c.1771-14T>A.
Non-syndromic hearing loss gene variants. Fourteen families segregated variants in
GJB2 (OMIM: 121011), while five families segregated variants in
SLC26A4 (OMIM: 605646), which are the most common HL genes in East Asia [
17,
18]. All the variants (N=12) identified in these families were classified as pathogenic and likely pathogenic based on the ACMG-AMP classification guidelines for HL.
MYO7A (OMIM: 276903) variants c.1183C>T (p.Arg395Cys) and c.3696_3706del (p.Arg1232SerTer72) were segregated from NSHL phenotype of proband in Fam 20. It has been Fam 21 was found to have a
DIAPH3 (OMIM: 603550) c.2256_2257insT (p.Ser752SerfsTer12) variant segregating with the HL phenotype. In Fam 21, the mother of the proband was confirmed by Sanger sequencing to carry
DIAPH3 c.2256_2257insT (
Figure 2.A). The mother's three brothers also affected HL, although they did not undergo genetic testing.
PTPRQ (OMIM: 603317) is associated with ADHL (OMIM: 617663) and ARNSHL (OMIM: 613391). Novel heterozygous variants
PTPRQ c.6293T>C (p.Leu2098Ser) were segregating in Fam 22 (
Figure 2.B). Proband showed normal vestibular function until now.
LOXHD1 (OMIM: 613072) is associated with ARNSHL (OMIM: 613079). Novel compound heterozygous variants
LOXHD1 c.2438T>A (p.Leu813Ter) and c.1759C>T (p.Arg587Trp) were segregating in Fam 23 (
Figure 2.C).
CDH23 (OMIM: 605516) is involved in Usher syndrome 1D (OMIM: 601067) and ARNSHL (OMIM: 601386). Fam 24 was found to have Homozygous
CDH23 c.4859T>A (p.Val1620Glu) segregating with the profound HL phenotype (
Figure 2.D). The proband was 13 years old at the time of testing and had not yet exhibited vestibular dysfunction or pigmentary retinopathy. Fam 25 segregates compound heterozygous variants
MPZL2 (OMIM: 604873) c.393_436+21del and c.220C>T (p.Gln74Ter) (
Figure 2.E).
MPZL2 has been associated with ARNSHL(OMIM: 618145). The feature of Fam26 are an HL father and two HL boys, which both carring
EYA4 (OMIM: 603550) c.1759C>T ( p.Arg587Ter) variant (
Figure 3.F).
EYA4 has been associated with ARNSHL(OMIM: 601316). No clinical signs were observed after the genetics tests that would indicate that affected family members have disorder affecting other organs relate to syndromic HL. We also identified a novel compound heterozygous variants c.2359G>T(p.Glu787Ter) and c.2353A>C(p.Thr785Pro) in
OTOA (
Figure 2.G).
OTOA (OMIM: 607038) has been associated with ARNSHL(607039).
Syndromic hearing loss gene variants. Novel heterozygous variants
TCOF1 (OMIM: 606847) c.3997_4007del (p.Ser1333GlnfsTer16) were identified to be likely cause of SHL in Fam 28 (
Figure 3.A). The proband in Fam28 has been clinically diagnosed with Treacher Collins syndrome (OMIM:154500). It was reported that clinical features of include downward-slanting palpebral fissures, malar hypoplasia, conductive deafness, mandibular hypoplasia, atresia of external ear canal, microtia, coloboma of the lower eyelid, asymmetry, projection of scalp hair onto the lateral cheek, cleft palate, choanal stenosis/atresia, cardiac malformation. Two families (Fam29 and Fam30) segregate variants
EYA1 (OMIM: 601653) c.1225C>T (p.Asn451ArgfsTer18) and c.1350_1353delTAATinsCAGACA (p.Pro458Ala) (
Figure 3.B,C).
EYA1 is associated with Branchiootic syndrome 1 (OMIM: 602588). Both probands were observed to have clinical features of branchial cervical fistulae, preauricular pits and cup shaped Pinnae. Novel compound heterozygous variants
FDXR c.1069G>T (p.Val357Leu) and c.364C>T (p.Arg122Cys) were identified to be likely cause of SHL in Fam 31 (
Figure 3.D).
FDXR (OMIM: 103270) variants
FDXR c.364C>T (p.Arg122Cys) and c.1069G>T (p.Val357Leu) could cause Auditory neuropathy and optic atrophy (OMIM: 617717) or Multiple mitochondrial dysfunctions syndrome 9B (OMIM: 620887). Proband with severe HL in Fam 31 have bilateral optic nerve atrophy and retinitis pigmentosa, while the brother of the probands in Fam 31 have Down syndrome(OMIM: 190685). Sanger sequencing confirms that the mother of the proband carries a heterozygous variant c.364C>T (p.Arg122Cys) in
FDXR, and the father of the proband does not carry
FDXR variant. In Fam 33, novel heterozygous variants
SOX10 (OMIM: 602229) C.133del (p.Gly38AlafsTer71) and heterozygous variants
SLC26A4 (OMIM: 605646) c.1975G>C were both identified in proband (
Figure 3.E).
SOX10 variants has been associated with Waardenburg syndrome 2E (OMIM:611584) and Waardenburg syndrome 4C (OMIM:613266). The proband have blue eyes and profound HL. The mother of the proband in Fam 33, also have blue eyes and HL, and was identfied carring heterozygous variant C.133del in
SOX10 by WES sequencing. The father of the proband in Fam 33, with the HL phenotype, and was identfied carring compound heterozygous variants c.1975G>C and c.1919G>A (p.Trp640Ter) in
SLC26A4. In addition, the younger brother of the proband's father in Fam 33 also suffers from HL and was identified to carry a compound heterozygous variant c.1975G>C and c.1919G>A in
SLC26A4. Additionally, the proband in Fam 34 was identified to carry novel variants
AIFM1 c.1771-14T>A (
Figure 3.F).
AIFM1 (OMIM: 300169) have been associated with X-linked deafness-5 with peripheral neuropathy (DFNX5, OMIM: 300614), Cowchock syndrome (OMIM: 300169), and Combined oxidative phosphorylation deficiency 6 (OMIM: 300816). The proband in Fam 34 showed delayed development at 6 months after birth. Then the proband underwent a series of clinical examinations and was diagnosed with profound sensorineural hearing loss. MRI shows bilateral widening ventricular. The proband has multiple epileptic seizures. By age 5 years, the proband could not stand or walk alone, has poor eye contact, and is almost unresponsive to external sound stimuli. At present, the proband underwent cochlear implantation at the age of 6 and received continuous speech and motor rehabilitation training. However, the parents of the proband do not believe that cochlear implantation and speech rehabilitation training increase the proband's response to external sound stimuli. The parents of the proband do not carry the variant in
AIFM1. The
AIFM1 variant in the proband may be spontaneous mutation.
Hearing loss candidate genes variants. ATP7B (OMIM: 606882) is with autosomal recessive disorder Wilson disease (OMIM: 277900). In Fam 32, the proband was identified to carry compound heterozygous variant of ATP7B c.4014T>A (p.Ile1338Ile) and c.3446G>A (p.Gly1149Glu) (
Figure 4. A). Variant
ATP7B c.3446G>A was predicted to be likely pathogenic [
19], while the clinical features of
ATP7B c.4014T>A carrier could be conflicting interpretat according to gnomAD. It has been reported that synonymous mutation c.4014T>A(p.Ile1338Ile) could cause exon skipping in the
ATP7B mRNA transcript [
20].
Expression of HL Candidate Genes ATP7B in the Mouse Inner Ear
We used various publicly available RNA sequencing and microarray data sets to investigated expression of HL candidate gene
ATP7B in the developing and adult mouse cochlear.
ATP7B shows a low mRNA expression in the developing cochlear epithelium include cochlear hair cells(HC) and cochlear surrounding cells(SC) throughout E14, P0, P4 and P7 (
Figure 4. B). The expressions of
ATP7B mRNA were detected in pillar cells (PC) at adult CBA/J mice cochlear, while less
ATP7B mRNA were detected in inner hair cells (IHC), outter hair cells (OHC), and Deiters’ cells (DC) (
Figure 4. C).
The onset of mice hearing arround P12 and the hearing threshold of mice remains stable around P18. To detect the expression of ATP7B protein in the adult mouse cochlea, we performed immunofluorescence staining on cochlear tissues of C57BL/6 mice at P25. ATP7B were detect to mainly expressed in the spiral ganglions (SGN) and organ of Corti (OC), include IHC, OHC, PC, DC, and Claudius' cells (CC) (
Figure 5. A-L). Immunofluorescence results also demonstrates that ATP7B were expressed in the cytoplasmic region of the cochlear outer hair cells (
Figure 5. M-O) and spiral ganglions cell (
Figure 6. A-H).
3.2. Figures, Tables and Schemes
Figure 1.
Principal component analysis and overview of gene variants in this cohort. A: Frequency of the known HI and candidate genes in chinese families included in this cohort. B: Proportions of non-syndromic hearing loss and syndromic hearing loss families in this cohort. C: Proportions of hereditary mode of variants in known and novel candidate HL genes were identified. D: Proportions of pre-lingual and post-lingual hearing loss of the probands in this cohort. E: Proportions of variant types in known and novel candidate HL genes were identified.
Figure 1.
Principal component analysis and overview of gene variants in this cohort. A: Frequency of the known HI and candidate genes in chinese families included in this cohort. B: Proportions of non-syndromic hearing loss and syndromic hearing loss families in this cohort. C: Proportions of hereditary mode of variants in known and novel candidate HL genes were identified. D: Proportions of pre-lingual and post-lingual hearing loss of the probands in this cohort. E: Proportions of variant types in known and novel candidate HL genes were identified.
Figure 2.
Pedigrees of families segregation novel non-sydromic hearing loss genes variants. The segregation of non-sydromic hearing loss genes was shown in the respective families. A: Fam21, B: Fam22, C: Fam23, D: Fam24, E: Fam25, F: Fam26, G: Fam27. * Individuals whose audiological examination results can not obtained.
Figure 2.
Pedigrees of families segregation novel non-sydromic hearing loss genes variants. The segregation of non-sydromic hearing loss genes was shown in the respective families. A: Fam21, B: Fam22, C: Fam23, D: Fam24, E: Fam25, F: Fam26, G: Fam27. * Individuals whose audiological examination results can not obtained.
Figure 3.
Pedigrees of families segregation sydromic hearing loss genes variants. The segregation of sydromic hearing loss genes was shown in the respective families. A: Fam28, B: Fam29, C: Fam30, D: Fam31, E: Fam33, F: Fam34.
Figure 3.
Pedigrees of families segregation sydromic hearing loss genes variants. The segregation of sydromic hearing loss genes was shown in the respective families. A: Fam28, B: Fam29, C: Fam30, D: Fam31, E: Fam33, F: Fam34.
Figure 4.
Pedigrees of Fam32 segregation ATP7B variants and expression of novel candidate gene ATP7B in the development and adult cochlear. A: Pedigrees of Fam32 segregation ATP7B variants. B: ATP7B mRNA expression in the hair cells(HC) and surrounding cells(SC) of mouse cochlea and utricle during four different developmental stages: E16, P0, P4, and P7. RNA-sequencing data of hair cells (GFP+) and surrounding cells (GFP−) from the cochlear and utricles of mouse expressing EGFP under the Pou4f3 promoter during developmental stages. C: ATP7B mRNA expression in the inner ear cells of adult CBA/J mice. ATP7B mRNA expression in Deiters’ cells (DC), Pillar cells (PC), Inner Hair Cells (IHC), and Outer Hair Cells (OHC) were evaluated. The y-axis represents the gene expression normalized to transcripts per million (TPM). Data for (A,B) were obtained from Gene Expression Omnibus (GEO) database and SHIELD (Shared Harvard Inner-Ear La-boratory Database).
Figure 4.
Pedigrees of Fam32 segregation ATP7B variants and expression of novel candidate gene ATP7B in the development and adult cochlear. A: Pedigrees of Fam32 segregation ATP7B variants. B: ATP7B mRNA expression in the hair cells(HC) and surrounding cells(SC) of mouse cochlea and utricle during four different developmental stages: E16, P0, P4, and P7. RNA-sequencing data of hair cells (GFP+) and surrounding cells (GFP−) from the cochlear and utricles of mouse expressing EGFP under the Pou4f3 promoter during developmental stages. C: ATP7B mRNA expression in the inner ear cells of adult CBA/J mice. ATP7B mRNA expression in Deiters’ cells (DC), Pillar cells (PC), Inner Hair Cells (IHC), and Outer Hair Cells (OHC) were evaluated. The y-axis represents the gene expression normalized to transcripts per million (TPM). Data for (A,B) were obtained from Gene Expression Omnibus (GEO) database and SHIELD (Shared Harvard Inner-Ear La-boratory Database).
Figure 5.
Immunostaining of ATP7B in the cochlear of C57BL/6 mice at P18. A-C: Whole mount immunostaining of the cochlear of P18 C57BL/6 mice. Immunostaining of ATP7B (Green) is observed in both spiral ganglions cell and organ of Corti. Scale bar: 100μm. D-L: Immunostaining of ATP7B (Green) in regions of the OC at the apical, medial and basal turns of the cochlea specimens of C57BL/6 mice at P18. Immunostaining of ATP7B (Green) is observed in both IHCs, OHCs, DCs, PCs, and claudius cell. Scale bar: 100μm. M-O: Immunostaining of ATP7B (Green) in cytoplasmic region of the cochlear outer hair cells of C57BL/6 mice at P18. Nucleus was stained with DAPI (Blue). F-actin was stained with phalloidin (red). Scale bar:10μm.
Figure 5.
Immunostaining of ATP7B in the cochlear of C57BL/6 mice at P18. A-C: Whole mount immunostaining of the cochlear of P18 C57BL/6 mice. Immunostaining of ATP7B (Green) is observed in both spiral ganglions cell and organ of Corti. Scale bar: 100μm. D-L: Immunostaining of ATP7B (Green) in regions of the OC at the apical, medial and basal turns of the cochlea specimens of C57BL/6 mice at P18. Immunostaining of ATP7B (Green) is observed in both IHCs, OHCs, DCs, PCs, and claudius cell. Scale bar: 100μm. M-O: Immunostaining of ATP7B (Green) in cytoplasmic region of the cochlear outer hair cells of C57BL/6 mice at P18. Nucleus was stained with DAPI (Blue). F-actin was stained with phalloidin (red). Scale bar:10μm.
Figure 6.
Immunostaining of ATP7B in the spiral ganglion cell of C57BL/6 mice at P18. A-H: Whole mount immunostaining of ATP7B in the spiral ganglion cell of P18 C57BL/6 mice. Nucleus was stained with DAPI (Blue). F-actin was stained with phalloidin (red). Immunostaining of ATP7B (Green) is observed in spiral ganglions cell at apical, middle and basic turn. Scale bar: 50μm.
Figure 6.
Immunostaining of ATP7B in the spiral ganglion cell of C57BL/6 mice at P18. A-H: Whole mount immunostaining of ATP7B in the spiral ganglion cell of P18 C57BL/6 mice. Nucleus was stained with DAPI (Blue). F-actin was stained with phalloidin (red). Immunostaining of ATP7B (Green) is observed in spiral ganglions cell at apical, middle and basic turn. Scale bar: 50μm.