Clinical data and hearing phenotypes
In total, 18 male patients with IP-III malformations were recruited for this study (Fig. 1). The patients were from 14 unrelated Chinese families including thirteen Han and one Zang family (Family 08). The mean age at identification of hearing loss was 1.1 years (range: 0.3–3.6). All patients had bilateral, prelingual, sensorineural or mixed deafness. The results of click-ABR of 14 pediatric patients ranged from 50 to > 97 dB nHL. With respect to the tympanograms, 67% (12/18) of them presented a type “A” tympanogram, 17% (3/18) had type “C”, and 11% (2/18) presented type “B”. Among them, 4 had other problems in addition to hearing loss (Table 1). To note, patient 4 − 1 was diagnosed with mild autism and had undergone bilateral tympanotomy tube insertion when he was 7 months old to treat “secretory otitis media (SOM)”.
Table 1
Clinical and molecular genetic characteristics of patients with IP-III malformation
Family | Patient | Age at detection (y) | Age at visit (y) | Hearing loss level a | Tympanogram | Other problems | Hearing intervention (side and age at which CIb or HAc was fitted) | POU3F4 variants |
01 | 1–1 | 0.5 | 4.2 | L 85; R 60 dB nHL | A | / | HA (L + R, 1.7y) | c.962T > G (p.Val321Gly) |
02 | 2 − 1 | 0.4 | 2.4 | L 80; R 70 dB nHL | A | / | CI (L, 4.0y) + HA (R, 0.7y) | c.541C > T (p.Gln181*) |
03 | 3 − 1 | 0.4 | 2.3 | L 65; R 65 dB nHL | A | | CI (R,2.7y) + HA (L,0.5y) | c.699C > A (p.Cys233*) |
04 | 4 − 1 | 0.5 | 2.1 | L > 97; R 80 dB nHL | C | Mild autism | CI (L,4.5y) + HA (R,2.5y) | g.81548899_82006629delins ACCAATTGGTAGTACAAT |
05 | 5 − 1 | 0.6 | 3.8 | L 60; R 65 dB nHL | A | / | HA (L + R, 3.8y) | c.644T > G (p.Val215Gly) |
5 − 2 | 3.4 | 23.5 | L 65; R 67.5 dB HL | A | / | HA (L + R, 14.0y) | c.644T > G (p.Val215Gly) |
06 | 6 − 1 | 0.5 | 0.5 | L > 97; R > 97 dB nHL | A | / | CI (R, 1.0y) | g.81806051_82292259del |
6 − 2 | 3.6 | 67.2 | L 106; R 107.5 dB HL | A | / | None | g.81806051_82292259del |
07 | 7 − 1 | 0.5 | 2.0 | L > 97; R > 97 dB nHL | | / | CI (R, 2.5y) | c.648dupG (p.Leu217ValfsTer9) |
08 | 8 − 1 | 3.5 | 14.8 | L 115; R 116 dB HL | C | / | CI (R, 15.0y) | c.845G > A (p.Arg282Gln) |
8 − 2 | 3.0 | 12.0 | L 107.5; R 105 dB HL | C | Cataract | CI (R, 12.0y) | c.845G > A (p.Arg282Gln ) |
09 | 9 − 1 | 0.6 | 0.9 | L 70; R > 97 dB nHL | A | Developmental retardation | CI (R, 1.2y) + HA (L, 0.9y) | c.171G > A (p.Trp57*) d |
10 | 10 − 1 | 0.4 | 2.0 | L > 97; R > 97 dB nHL | A | / | CI (R, 3.0y) | g.81839469_82004841del |
10 − 2 | 0.4 | 2.0 | L > 97; R > 97 dB nHL | A | / | CI (R, 2.4y) | g.81839469_82004841del |
11 | 11 − 1 | 0.3 | 0.7 | L 65; R 60 dB nHL | A | / | HA (L + R, 0.7y) | c.946C > T (p.Gln316*) |
12 | 12 − 1 | 0.4 | 0.7 | L 50; R 55 dB nHL | B | / | HA (L + R, 1.0y) | c.903_912 delins TGCCA (p.Lys302AlafsTer25) |
13 | 13 − 1 | 0.3 | 4.5 | L 75; R 60 dB nHL | B | / | HA (L + R, 4.5y) | c.614_616delGAA (p.Arg205del) d |
14 | 14 − 1 | 0.5 | 0.6 | L 65; R 60 dB nHL | A | Atrial septal defect | HA (L + R, 0.6y) | g.81807331_81887213del d |
a dB nHL: click-ABR results, dB HL: PTA results. |
b CI: cochlear implantation. c HA: hearing aid. |
d De novo mutation. |
Among the mothers of all probands, subject 5 − 3 was 47 years of age and reported a history of left myringoplasty about two decades earlier. She presented mixed hearing loss in the left ear (PTA = 48.75 dB HL) and sensorineural hearing loss in the right ear (PTA = 40 dB HL). Meanwhile, HRCT scans showed unaffected inner ears.
Radiological findings
The HRCT images of all 18 patients showed typical IP-III anomalies: bilateral and symmetrical malformation, a relatively normal shape of outer coating with absence of cochlear modiolus and bony spiral lamina. There was a direct intercommunication between cochlea and IAC with a “cloudy like” characteristic (Fig. 2a I–III). There was a high jugular bulb in 83.3% (15/18) of patients (Fig. 2a IV). Stapes abnormalities were present in 11.1% (2/18) of patients and included a thickened stapes footplate and absence of fissula ante fenestram (Fig. 2a V). A vestibular aqueduct was occasionally visible, with a normal shape on axial CT images (Fig. 2a VI). The semicircular canals and vestibular aqueduct shapes were also normal. Postoperative HRCT was performed in six CI users. As examples, straight electrode arrays were successfully implanted on the right side in patients 7 − 1 and 10 − 2 (Fig. 2b I–II).
Mutations in POU3F4 gene by Sanger sequencing
From Sanger sequencing, 12 patients from 10 families had hemizygous variants in POU3F4 (Table 1) including 4 nonsense mutations, 3 missense variants, 2 frame-shift mutations and one indel mutation. Apart from c.648dupG identified in Family-07, the other 9 variants were novel. Besides, two variants p.Trp57* and p.Arg205del were de novo mutations for the patients from Family-09 and − 13 cause their mothers had no variants of POU3F4. All these mutations were not present in 200 normal hearing Chinese controls (150 Han and 50 Zang).
Four nonsense mutations, p.Gln181*, p.Cys233*, p.Trp57* and p.Gln316* were identified in Family-02, -03, -09 and − 11, respectively (Fig. 3a). Mutations p.Trp57* and p.Gln181* resulted in the loss of the entire POU-specific (POUs) domain and the POU homeodomain (POUHD). In particular, proband 9 − 1 had a de novo p.Trp57* hemizygous mutation. p.Cys233* was located in the POUs domain and p.Gln316* was present in the POUHD. Both caused premature termination in POU3F4. These residues are highly conserved in species including primates, mammals and other vertebrates (Fig.S1ab).
Three missense variants p.Val321Gly, p.Val215Gly and p.Arg282Gln were identified in Family-01, -05 and − 08, respectively (Fig. 3a). These residues are also highly conserved in different species (Fig.S1ab). These three variants were predicted to be disease causing and deleterious by the three prediction tools used in this study (Table.S1).
Two frame-shift mutations c.648dupG and c.903_912 delins TGCCA were identified in Family-07 and − 12, respectively. They caused amino acid mutations of p.Leu217ValfsTer9 and p.Lys302AlafsTer25, respectively (Fig. 3a), and resulted in premature termination of POU3F4.
The indel mutation c.614_616delGAA (p.Arg205del), identified in Family-13, was also a de novo mutation in this study. This residue is also highly conserved in different species (Fig.S1a). The predicted structures of wild-type and mutant proteins were observed and analyzed using RasMol (Fig.S1c).
Deletions identified in upstream of POU3F4 by Nanopore long-read single-molecule sequencing
Sanger sequencing of POU3F4 was negative in 4 families (Family-04, -06, -10 and − 14), so Nanopore long-read single-molecule sequencing was conducted. A total of 40–50 Gb of reads with a mean length of 11.6–22 kb was obtained in each sample with an average coverage of 95.3%, and an average depth of 13–16× for the whole genome. Four different types of DELs were identified (Fig. 3b). The breakpoint junctions of the DELs were verified by Sanger sequencing and the sequences of primers were listed in Table S2.
In proband 4 − 1, we detected a 458 kb DEL and an 18 bp INS located at ChrX (g.81548899_82006629delinsACCAATTGGTAGTACAAT), approximately 1502 kb upstream of POU3F4. A total of 7 reads supported the DEL. The breakpoint junction and insertion were successfully sequenced in the proband and his mother (Fig. 3c I).
In proband 6 − 1, we identified a 486 kb DEL (g.81806051_82292259del), 1217 kb upstream of POU3F4. A total of 7 reads supported the DEL. The junction was successfully sequenced in the proband, his mother and his grandfather (Fig. 3c II).
In patient 10 − 1, a 165 kb DEL (g.81839469_82004841del) was detected, about 1503 kb upstream of the POU3F4 gene. A total of 8 reads supported the DEL. The junction sequence was successfully amplified for the proband, his twin brother and their mother (Fig. 3c III).
In proband 14 − 1, we detected a 80 kb DEL (g.81807331_81887213del), 876 kb upstream of the POU3F4 gene, and this mutation was successfully sequenced (Fig. 3c IV). His mother’s genotype was normal. The proband had a de novo DEL located upstream of POU3F4.
Effect of hearing intervention and its correlation with genotypes
Appropriate hearing intervention was carried out in 94.4% (17/18) of patients in this study. Ten patients received unilateral CI with or without HA on the contralateral side. Seven patients wore bilateral HAs. Only one patient did not receive hearing intervention: Patient 6 − 2, 67 years of age, had congenital hearing loss, never had a HA and used sign language in daily life.
For 10 CI users, mean aided PTA was 40.3 ± 7.59 dB HL (range: 25.0–51.3 dB HL) 12 months after activation of the implant. Compared with the control group (age-matched CI recipients with normal cochlea, N = 20, PTA = 34.0 ± 5.74 dB HL), the result was not statistically different (p = 0.20). The genotypes of these IP-III patients consisted of 3 nonsense mutations (p.Gln181*, p.Cys233*, p.Trp57*), one frame-shift mutation (c.648dupG), one missense variant (p.Arg282Gln) and 3 upstream DELs (165–486 kb).
For 7 patients with bilateral HAs, mean aided PTA of the better ear was 41.1 ± 5.18 dB HL (range: 35.0–50.0 dB HL). The genotypes of them were p.Val215Gly (two siblings), p.Val321Gly, p.Arg205del, c.903_912 delins TGCCA, p.Gln316* and g.81807331_81887213del (80 kb DEL).
Based on the location of variants, 17 patients with HAs or CI were further divided into two groups: variants located in the exon region of POU3F4 (ER group, N = 12) or the upstream region (UR group, N = 5). Mean aided PTA of the ER group was 39.6 ± 6.31 dB HL, and that of the UR group was 43.0 ± 7.10 dB HL. The difference was not significant (p = 0.342).