# Hearing Screening Combined with Target Gene Panel Testing Increased Etiological Diagnostic Yield in Deaf Children

Authors: Le Xie; Yue Qiu; Yuan Jin; Kai Xu; Xue Bai; Xiao-Zhou Liu; Xiao-Hui Wang; Sen Chen; Yu Sun
Journal: Neural Plasticity
 (2021)

---

## Abstract

Genetic testing is the gold standard for exploring the etiology of congenital hearing loss. Here, we enrolled 137 Chinese patients with congenital hearing loss to describe the molecular epidemiology by using 127 gene panel testing or 159 variant testing. Sixty-three deaf children received 127 gene panel testing, while seventy-four patients received 159 variant testing. By use of 127 gene panel testing, more mutant genes and variants were identified. The most frequent mutant genes wereGJB2, SLC26A4, MYO15A, CDH23, and OTOF. By analyzing the patients who received 127 gene panel testing, we found that 51 deaf children carried variants which were not included in 159 variant testing. Therefore, a large number of patients would be misdiagnosed if only 159 variant testing is used. This study highlights the advantage of 127 gene panel testing, and it suggests that broader genetic testing should be done to identify the genetic etiology of congenital hearing loss.


---

## Body

## 1. Introduction

Congenital hearing loss (HL) is a common disease, and about 1 to 2 per 1000 live births suffer from congenital hearing loss in the world. The prevalence of hearing loss increases with age. Among the kids in the primary school, the prevalence is 2.83/1000. Furthermore, in adolescents, it will rise to 3.5/1000 [1–3]. According to the latest date, approximately 50%-60% of congenital HL is caused by genetic factors [4, 5]. However, genetic disorder can also cause late-onset deafness in children and adolescents [2]. Both congenital or late-onset deafness caused by genetic disorder can be attributed to hereditary hearing loss. To date, more than 150 deafness-related genes and 6000 variants have been identified [6] (https://hereditaryhearingloss.org). The hearing loss genes and hot variants showed distinctly in different ethnic groups. Among these hearing loss genes, GJB2 and SLC26A4 gene variants were most common in Chinese deaf population [4]. For GJB2 gene variants, the most frequent etiological factor of nonsyndromic hereditary hearing loss, the hotspots were described as c.235delC and c.109G>A, common in Asian, c.35delG, found in European, and c.71G > A, predominant in the Indian [7, 8]. Although GJB2 and SLC26A4 mutations account for a large part of causes of hereditary deafness, there are still a large number of known or unknown gene mutations that cause hearing loss. How to quickly and accurately detect causes of patients with suspected hereditary deafness is an important issue in the diagnosis and treatment of deafness.Genetic testing is the gold standard for exploring the etiology of congenital and late-onset deafness. In addition, medical history and other auditory physiology tests are also very important. For hereditary deafness, genetic testing is a crucial step. Except for finding causes, it also could help diagnosis and intervention of syndromic hearing loss because some children only exhibited hearing loss without other symptoms at a young age. Besides, genetic testing was beneficial to protecting children with some variants such as m.1555A > G or variants inKCNQ1 gene from avoidable risk factors [9–12]. DNA sequencing first described 44 years ago has evolved from single-mutation sequencing (Sanger sequencing) to high-throughput sequencing (Next generation sequencing, NGS) [5, 13]. Compared to Sanger sequencing, NGS was able to simultaneously sequence millions of small fragments at a reasonable cost and reduced runtime [13]. Targeted gene panel, a common NGS technique, has been widely applied in deaf population. It was targeted to the detection of variants in the massively known hearing loss genes. Considering the cost and time consuming, different panels were designed in clinical applications [1, 4, 12, 14, 15].As the spectrum of gene mutations is distinct in different ethnic group, studying mutated character of Chinese deaf children and choosing the appropriate gene panel are critical for clinical diagnosis and treatment. In this study, we enrolled 137 Chinese children with sensorineural hearing loss to describe the molecular epidemiology by using 127 gene panel testing and 159 variants in 22 gene testing.


## 2. Materials and Methods

### 2.1. Patients and Samples

A total of 137 patients under the age of twelve with sensorineural hearing loss were recruited from the Wuhan Union Hospital from 2018 to 2020. All patients failed pass the neonatal hearing screening and underwent 127 gene panel testing (targeting the exon regions and exon-intron boundaries of 127 known deafness-causing nuclear genes as well as deafness-causing mitochondrial regions) or 159 variants in 22 gene testing (detecting a total of 159 hotspot variants in 22 known deafness-causing genes). Informed consent was obtained from all patients or their legal guardians. This study was approved by the Institutional Review Board of the BGI in accordance with the Declaration of Helsinki (1964). Written informed consent, clinical evaluations, and blood samples were obtained from all the participants or their legal guardians. The study of the protocols was approved by the review boards of the ethics committees of the Tongji Medical College of Huazhong University of Science and Technology.Peripheral venous blood samples were collected from all recruited patients. Genomic DNA was obtained and purified using a QIAamp DSP DNA Blood Mini Kit (61104, Qiagen Inc., Valencia, CA, USA).


### 2.2. Library Preparation, Sequencing, and Bioinformatics

Fragmented DNA (a size of 350–400 base pairs) were prepared using Covaris LE220 ultrasonicator (Covaris Inc., Woburn, Massachusetts, USA). Then end-repair and adaptor ligation were performed for library construction. After array hybridization, elution, and postcapture amplification, targeted DNA fragments were sequenced on the BGISEQ-500 platform. After sequencing and quality control, the clean reads derived from targeted high-throughput sequencing were aligned to the GRCh37/hg19 using the BWA (Burrows Wheeler Aligner) Multi-Vision software package. After alignment, the variants of single-nucleotide variants (SNVs) and inserts and deletions (InDels) were detected by GATK software. Then, filtered SNVs and InDels were compared with NCBI GenBank database, 1000 Genomes, ESP6500, dbSNP, HGMD, and ExAC. Candidate variants were classified into pathogenic variants, likely pathogenic variants, variants of uncertain significance (VUS), likely benign variants, and benign variants according to the American College of Medical Genetics and Genomics–Association for Molecular Pathology (ACMG–AMP) guideline. The method described in this study partly reproduces the wording in our previous article [16–18].


## 2.1. Patients and Samples

A total of 137 patients under the age of twelve with sensorineural hearing loss were recruited from the Wuhan Union Hospital from 2018 to 2020. All patients failed pass the neonatal hearing screening and underwent 127 gene panel testing (targeting the exon regions and exon-intron boundaries of 127 known deafness-causing nuclear genes as well as deafness-causing mitochondrial regions) or 159 variants in 22 gene testing (detecting a total of 159 hotspot variants in 22 known deafness-causing genes). Informed consent was obtained from all patients or their legal guardians. This study was approved by the Institutional Review Board of the BGI in accordance with the Declaration of Helsinki (1964). Written informed consent, clinical evaluations, and blood samples were obtained from all the participants or their legal guardians. The study of the protocols was approved by the review boards of the ethics committees of the Tongji Medical College of Huazhong University of Science and Technology.Peripheral venous blood samples were collected from all recruited patients. Genomic DNA was obtained and purified using a QIAamp DSP DNA Blood Mini Kit (61104, Qiagen Inc., Valencia, CA, USA).


## 2.2. Library Preparation, Sequencing, and Bioinformatics

Fragmented DNA (a size of 350–400 base pairs) were prepared using Covaris LE220 ultrasonicator (Covaris Inc., Woburn, Massachusetts, USA). Then end-repair and adaptor ligation were performed for library construction. After array hybridization, elution, and postcapture amplification, targeted DNA fragments were sequenced on the BGISEQ-500 platform. After sequencing and quality control, the clean reads derived from targeted high-throughput sequencing were aligned to the GRCh37/hg19 using the BWA (Burrows Wheeler Aligner) Multi-Vision software package. After alignment, the variants of single-nucleotide variants (SNVs) and inserts and deletions (InDels) were detected by GATK software. Then, filtered SNVs and InDels were compared with NCBI GenBank database, 1000 Genomes, ESP6500, dbSNP, HGMD, and ExAC. Candidate variants were classified into pathogenic variants, likely pathogenic variants, variants of uncertain significance (VUS), likely benign variants, and benign variants according to the American College of Medical Genetics and Genomics–Association for Molecular Pathology (ACMG–AMP) guideline. The method described in this study partly reproduces the wording in our previous article [16–18].


## 3. Result

### 3.1. Study Samples

A total of 137 patients were included in this research, seventy-four patients were male and sixty-three were female. Sixty-three patients (46%) received genetic testing by 127 gene panel, while seventy-four patients (54%) received 159 variant testing (Table1). Tested genes are provided in supplement (Table S1 and S2.)Table 1
The basic information of patients recruited in this research.
127 gene panel159 variantsTotalMale353974Female283563Total6374137


### 3.2. Genetic Spectrum

Among the 63 received 127 gene panel test, ninety-two variants in 32 genes were detected, the top five genes wereGJB2 (33 patients), SLC26A4 (16 patients), MYO15A (5 patients), CDH23 (5 patients), and OTOF (5 patients) (Table 2).Table 2
The genetic spectrum of hereditary hearing loss patients detected by 127 gene panel testing.
GenePatients carried variantPathogenic variantLikely pathogenic variantVUS variantGJB2333201SLC26A4161501OTOF5131MYO15A5113CDH235023TECTA4013MYO7A4004WFS13102PTPRQ2101COCH2020GJB32011OTOA2011USH2A2011ESPN2002TRIOBP2002GRXCR11100ALMS11100COL4A31010ILDR11010ADGR V11001CLDN141001CRYM1001DIAPH31001EDNRB1001LOXHD11001MITF1001MYH141001PCDH151001PDSS11001SLC17A81001TJP21001TMC11001Among the 74 who received 159 variant testing, twenty-one variants in 4 genes were detected; the patients with mutation ofGJB2, SLC26A4, GJB3, and MT-RNR1 were 45, 27, 2, and 1, respectively (Table 3).Table 3
The genetic spectrum of hereditary hearing loss patients detected by 159 variant testing.
GenePatients with mutant genePercentageGJB2450.608108SLC26A4270.364865GJB320.027027MT-RNR110.013514


### 3.3. Different Diagnosis between 127 Gene Panel Testing and 159 Variant Testing

Among the 63 patients who received 127 gene testing, fifty-one patients were found to carry variants which were not included in 159 variant testing, thirty-six of them had been found to carry pathogenic or likely pathogenic variants. If these patients received 159 variant testing only, in 36 of them, the genetic etiology of the hearing loss would be missed, and 26 of the patients were found to carry the pathogenic or likely pathogenic variants (Figure1). Twenty-two deafness gene mutation which was not included in 159 variant testing was found in 25 patients; nine of them carried pathogenic or likely pathogenic variants of these 22 deafness genes. In these 22 deafness genes which were not included in 159 variant testing, eighty-four variants were found, ten of them were pathogenic or likely pathogenic (Table 4).Figure 1
Analyze the genetic etiology of patients who received 127 gene panel testing compared with 159 variants. (a) 36 patients carried pathogenic or likely pathogenic variants which were not included in 159 variants, fifteen patients carried VUS variants which were not included in 159 variants, and 12 patients do not carry the variants which were not included in 159 variants. (b) 26 patients who carried pathogenic or likely pathogenic variants would misdiagnose if they were only detected by 159 variant testing, ten patients who carried VUS might miss diagnosis, and twenty-seven patients would not miss diagnosis.
(a)(b)Table 4
The variants would be missed if tested by 159 variant testing.
GeneVariantEffect on proteinCharacterization of variantADGR V1c.11704A>Gp.Met3902ValVUSc.2899-10T>AVUSALMS1c.2035C>Tp.Arg679TerPathogenicCDH23c.5067+1G>ALikely pathogenicc.6604G>Ap.Asp2202AsnLikely pathogenicc.1765G>Ap.Asp589AsnVUSc.2368A>Gp.Met790ValVUSc.3262G>Ap.Val1088MetVUSc.4859T>Ap.Val1620GluVUSc.5051G>Ap.Arg1684HisVUSc.805C>Tp.Arg269TrpVUSCLDN14c.449C>Tp.Pro150LeuVUSc.694G>Ap.Gly232ArgVUSCOCHc.812_813insGp.Val271Valfs X5Likely pathogenicCOL4A3c.4755+1G>ALikely pathogenicCRYMc.849C>Tp.His283HisVUSDIAPH3c.3543_3544i nsCp.Pro1181 Profs∗4VUSEDNRBc.758G>Ap.Arg253GlnVUSESPNc.1828G>Cp.Ala610ProVUSc.2069C>Tp.Ser690LeuVUSGJB2c.109G>Ap.Val37IlePathogenicc.139G>Tp.Glu47TerPathogenicc.88A>Gp.Ile30ValVUSGJB3c.580G>Ap.Ala194ThrLikely pathogenicc.595A>Gp.Ile199ValVUSGRXCR1c.784C>Tp.Arg262TerPathogenicILDR1c.206C>Ap.Pro69HisLikely pathogenicLOXHD1EX31 DUPVUSMITFc.950G>Cp.Arg317ThrVUSMYH14c.4641C>Tp.Asp1547AspVUSMYO15Ac.10251_10253 delCTTPathogenicc.4642G>Ap.Ala1548ThrPathogenicc.8702_8703i nsTp.Pro2901Pro fsX25Likely pathogenicc.3239G>Ap.Arg1080HisVUSc.3743G>Ap.Arg1248GlnVUSc.5681T>Cp.Leu1894ProVUSc.6340G>Ap.Val2114MetVUSMYO7Ac.1901G>Ap.Arg634GlnVUSc.1945C>Tp.Arg649TrpVUSc.6235C>Tp.Arg2079TrpVUSc.6470T>Cp.Ile2157ThrVUSOTOAc.943delTp.Ser315Profs∗5Likely pathogenicc.1172C>Tp.Ser391LeuVUSc.2359G>Tp.Glu787TerVUSOTOFc.5816G>Ap.Arg1939GlnPathogenicc.4023+1G>ALikely pathogenicc.5098G>Cp.Glu1700GlnLikely pathogenicc.1194T>Ap.Asp398GluVUSc.3429G>Ap.Arg1143ArgVUSPCDH15c.2563C>Ap.Arg855ArgVUSPDSS1EX1 DUPVUSPTPRQc.6475C>Tp.Arg2159TerPathogenicc.6025-4_6025-2delACAVUSc.6293T>Cp.Leu2098SerVUSSLC17A8c.1117T>Cp.Leu373LeuVUSc.243G>Ap.Met81IleVUSSLC26A4c.1079C>Tp.Ala360ValPathogenicc.1339delAp.Lys447Serfs∗8Pathogenicc.1519delTp.Leu507∗Pathogenicc.164+1G>CPathogenicc.2000T>Cp.Phe667SerPathogenicc.1001+5G>TVUSc.1087A>Cp.Ile363LeuVUSc.1340A>Tp.Lys447MetVUSc.208C>Ap.Pro70ThrVUSc.678T>Cp.Ala226AlaVUSc.718C>Tp.Leu240PheVUSc.765+5G>AVUSTECTAc.640G>Tp.Gly214TerLikely pathogenicc.235G>Cp.Val79LeuVUSc.2936T>Cp.Phe979SerVUSc.576G>Ap.Thr192ThrVUSTJP2c.474G>Tp.Arg158SerVUSTMC1c.1810C>Gp.Arg604GlyVUSc.589G>Ap.Gly197ArgVUSTRIOBPc.2133T>Cp.Pro711ProVUSc.2314T>Cp.Cys772ArgVUSc.-60-1G>CVUSUSH2Ac.6937G>Tp.Gly2313CysLikely pathogenicc.10740+7G>AVUSc.7230A>Tp.Val2410ValVUSWFS1c.2051C>Tp.Ala684ValPathogenicc.125G>Tp.Arg42LeuVUSc.2458G>Ap.Gly820SerVUS


### 3.4.GJB2 Variants

In 127 gene panel testing, six variants ofGJB2 were identified; c.109G>A, c.235delC, c.299_300del AT, c.176_191del GCTGCAAGAACGTGTG, c.139G>T, and c.88A>G were detected 18, 15, 5, 2, 1, and 1 times, respectively (Table 5). The patients who carried heterozygous variant of c.109G>A were the most, and secondary was heterozygous c.235delC mutation in 11 patients, and then was the homozygous c.109G>A mutation in 5 patients and homozygous c.235delC mutation in 4 patients. In the 11 patients who carried heterozygous c.235delC variants, four of them carried c.299_300del AT, 2 carried c.109G>A, and 2 carried c.176_191del GCTGCAAG AACGTGTG compound heterozygotes variants, respectively. In our study, a patient who carried a novel heterozygous variant was identified, the nucleotide change is c.88A>G, the protein change is p.Ile30Val, and the characterization of variant is VUS; this patient carried c.919-2A>G and c.1229C>T compound heterozygotes variants of SLC26A4 gene as well.Table 5
Variants frequency inGJB2 gene.
VariantsHomozygousHeterozygousTotalTesting by 127 gene panelc.109G>A51318c.235delC41115c.299_300del AT145c.176_191del GCTGCAAG AACGTGTG022c.139G>T101c.88A>G011Testing by 159 variant testingc.235delC171936c.299_300delAT099c.176_191delGCTGCAAGAACGTGTG055c.139G>T202c.187G>T011c.427C>T101c.428G>A011c.257C>G011In 159 variant testing, eight variants ofGJB2 were identified; the frequency of these variants from high to low were c.235delC, c.299_300delAT, c.176_191del GCTGCAAGAACGTGTG, c.139G>T, c.187G>T, c.257C>G, c.427C>T, and c.428G>A; the check-out times of these variants were 36, 9, 5, 2, 1, 1, 1, and 1, respectively. Nineteen patients carried heterozygous c.235delC variants, five of them carried c.299_300delAT compound heterozygotes variants, three patients carried c.176_191del GCTGCAAGAACGTGTG compound heterozygotes variants, and one carried c.257C>G compound heterozygotes variants. The homozygous c.235delC patients were seventeen. c.299_300delAT heterozygous variant was identified in 9 patients, two of them carried c.176_191delGCTGCAAGAACGTGTG compound heterozygotes variants.


### 3.5.SLC26A4 Variants

In 127 gene panel testing, seventeen variants were identified inSLC26A4 gene, the most frequent variant detected was c.919-2A>G, eleven patients carried heterozygous c.919-2A>G variant, nine of them were compound heterozygotes variants with other SLC26A4 variants, and one homozygous c.919-2A>G was identified. Eight novel variants were identified in our research; c.1339delA, c.1519delT, c.164+1G>C, and c.2000T>C were likely pathogenic variants, while c.1340A>T, c.718C>T, c.765+5G>A, and c.208C>A were VUS variants (Table 6).Table 6
Variants frequency inSLC26A4 gene.
VariantsHomozygousHeterozygousTotalTesting by 127 gene panelc.919-2A>G11112c.2000T>C022c.1087A>C011c.1174A>T011c.1339delA011c.164+1G>C011c.718C>T011c.1226G>A011c.678T>C011c.1519delT011c.1229C>T011c.1340A>T011Testing by 159 variant testingc.919-2A>G51419c.1229C>T055c.2168A>G044c.1174A>T011c.1226G>A011c.1336C>T011c.1343C>A011c.1343C>T011c.1594A>C011c.1692dupA011c.1707+5G >A011c.589G>A011Twelve variants ofSLC26A4 were identified in 159 variant testing; the most common variant was c.919-2A>G same as the 127 gene panel testing, then was c.1229C>T and c.2168A>G. Nineteen patients carried c.919-2A>G variant, fourteen of them were heterozygous while 5 were homozygous.


## 3.1. Study Samples

A total of 137 patients were included in this research, seventy-four patients were male and sixty-three were female. Sixty-three patients (46%) received genetic testing by 127 gene panel, while seventy-four patients (54%) received 159 variant testing (Table1). Tested genes are provided in supplement (Table S1 and S2.)Table 1
The basic information of patients recruited in this research.
127 gene panel159 variantsTotalMale353974Female283563Total6374137


## 3.2. Genetic Spectrum

Among the 63 received 127 gene panel test, ninety-two variants in 32 genes were detected, the top five genes wereGJB2 (33 patients), SLC26A4 (16 patients), MYO15A (5 patients), CDH23 (5 patients), and OTOF (5 patients) (Table 2).Table 2
The genetic spectrum of hereditary hearing loss patients detected by 127 gene panel testing.
GenePatients carried variantPathogenic variantLikely pathogenic variantVUS variantGJB2333201SLC26A4161501OTOF5131MYO15A5113CDH235023TECTA4013MYO7A4004WFS13102PTPRQ2101COCH2020GJB32011OTOA2011USH2A2011ESPN2002TRIOBP2002GRXCR11100ALMS11100COL4A31010ILDR11010ADGR V11001CLDN141001CRYM1001DIAPH31001EDNRB1001LOXHD11001MITF1001MYH141001PCDH151001PDSS11001SLC17A81001TJP21001TMC11001Among the 74 who received 159 variant testing, twenty-one variants in 4 genes were detected; the patients with mutation ofGJB2, SLC26A4, GJB3, and MT-RNR1 were 45, 27, 2, and 1, respectively (Table 3).Table 3
The genetic spectrum of hereditary hearing loss patients detected by 159 variant testing.
GenePatients with mutant genePercentageGJB2450.608108SLC26A4270.364865GJB320.027027MT-RNR110.013514


## 3.3. Different Diagnosis between 127 Gene Panel Testing and 159 Variant Testing

Among the 63 patients who received 127 gene testing, fifty-one patients were found to carry variants which were not included in 159 variant testing, thirty-six of them had been found to carry pathogenic or likely pathogenic variants. If these patients received 159 variant testing only, in 36 of them, the genetic etiology of the hearing loss would be missed, and 26 of the patients were found to carry the pathogenic or likely pathogenic variants (Figure1). Twenty-two deafness gene mutation which was not included in 159 variant testing was found in 25 patients; nine of them carried pathogenic or likely pathogenic variants of these 22 deafness genes. In these 22 deafness genes which were not included in 159 variant testing, eighty-four variants were found, ten of them were pathogenic or likely pathogenic (Table 4).Figure 1
Analyze the genetic etiology of patients who received 127 gene panel testing compared with 159 variants. (a) 36 patients carried pathogenic or likely pathogenic variants which were not included in 159 variants, fifteen patients carried VUS variants which were not included in 159 variants, and 12 patients do not carry the variants which were not included in 159 variants. (b) 26 patients who carried pathogenic or likely pathogenic variants would misdiagnose if they were only detected by 159 variant testing, ten patients who carried VUS might miss diagnosis, and twenty-seven patients would not miss diagnosis.
(a)(b)Table 4
The variants would be missed if tested by 159 variant testing.
GeneVariantEffect on proteinCharacterization of variantADGR V1c.11704A>Gp.Met3902ValVUSc.2899-10T>AVUSALMS1c.2035C>Tp.Arg679TerPathogenicCDH23c.5067+1G>ALikely pathogenicc.6604G>Ap.Asp2202AsnLikely pathogenicc.1765G>Ap.Asp589AsnVUSc.2368A>Gp.Met790ValVUSc.3262G>Ap.Val1088MetVUSc.4859T>Ap.Val1620GluVUSc.5051G>Ap.Arg1684HisVUSc.805C>Tp.Arg269TrpVUSCLDN14c.449C>Tp.Pro150LeuVUSc.694G>Ap.Gly232ArgVUSCOCHc.812_813insGp.Val271Valfs X5Likely pathogenicCOL4A3c.4755+1G>ALikely pathogenicCRYMc.849C>Tp.His283HisVUSDIAPH3c.3543_3544i nsCp.Pro1181 Profs∗4VUSEDNRBc.758G>Ap.Arg253GlnVUSESPNc.1828G>Cp.Ala610ProVUSc.2069C>Tp.Ser690LeuVUSGJB2c.109G>Ap.Val37IlePathogenicc.139G>Tp.Glu47TerPathogenicc.88A>Gp.Ile30ValVUSGJB3c.580G>Ap.Ala194ThrLikely pathogenicc.595A>Gp.Ile199ValVUSGRXCR1c.784C>Tp.Arg262TerPathogenicILDR1c.206C>Ap.Pro69HisLikely pathogenicLOXHD1EX31 DUPVUSMITFc.950G>Cp.Arg317ThrVUSMYH14c.4641C>Tp.Asp1547AspVUSMYO15Ac.10251_10253 delCTTPathogenicc.4642G>Ap.Ala1548ThrPathogenicc.8702_8703i nsTp.Pro2901Pro fsX25Likely pathogenicc.3239G>Ap.Arg1080HisVUSc.3743G>Ap.Arg1248GlnVUSc.5681T>Cp.Leu1894ProVUSc.6340G>Ap.Val2114MetVUSMYO7Ac.1901G>Ap.Arg634GlnVUSc.1945C>Tp.Arg649TrpVUSc.6235C>Tp.Arg2079TrpVUSc.6470T>Cp.Ile2157ThrVUSOTOAc.943delTp.Ser315Profs∗5Likely pathogenicc.1172C>Tp.Ser391LeuVUSc.2359G>Tp.Glu787TerVUSOTOFc.5816G>Ap.Arg1939GlnPathogenicc.4023+1G>ALikely pathogenicc.5098G>Cp.Glu1700GlnLikely pathogenicc.1194T>Ap.Asp398GluVUSc.3429G>Ap.Arg1143ArgVUSPCDH15c.2563C>Ap.Arg855ArgVUSPDSS1EX1 DUPVUSPTPRQc.6475C>Tp.Arg2159TerPathogenicc.6025-4_6025-2delACAVUSc.6293T>Cp.Leu2098SerVUSSLC17A8c.1117T>Cp.Leu373LeuVUSc.243G>Ap.Met81IleVUSSLC26A4c.1079C>Tp.Ala360ValPathogenicc.1339delAp.Lys447Serfs∗8Pathogenicc.1519delTp.Leu507∗Pathogenicc.164+1G>CPathogenicc.2000T>Cp.Phe667SerPathogenicc.1001+5G>TVUSc.1087A>Cp.Ile363LeuVUSc.1340A>Tp.Lys447MetVUSc.208C>Ap.Pro70ThrVUSc.678T>Cp.Ala226AlaVUSc.718C>Tp.Leu240PheVUSc.765+5G>AVUSTECTAc.640G>Tp.Gly214TerLikely pathogenicc.235G>Cp.Val79LeuVUSc.2936T>Cp.Phe979SerVUSc.576G>Ap.Thr192ThrVUSTJP2c.474G>Tp.Arg158SerVUSTMC1c.1810C>Gp.Arg604GlyVUSc.589G>Ap.Gly197ArgVUSTRIOBPc.2133T>Cp.Pro711ProVUSc.2314T>Cp.Cys772ArgVUSc.-60-1G>CVUSUSH2Ac.6937G>Tp.Gly2313CysLikely pathogenicc.10740+7G>AVUSc.7230A>Tp.Val2410ValVUSWFS1c.2051C>Tp.Ala684ValPathogenicc.125G>Tp.Arg42LeuVUSc.2458G>Ap.Gly820SerVUS


## 3.4.GJB2 Variants

In 127 gene panel testing, six variants ofGJB2 were identified; c.109G>A, c.235delC, c.299_300del AT, c.176_191del GCTGCAAGAACGTGTG, c.139G>T, and c.88A>G were detected 18, 15, 5, 2, 1, and 1 times, respectively (Table 5). The patients who carried heterozygous variant of c.109G>A were the most, and secondary was heterozygous c.235delC mutation in 11 patients, and then was the homozygous c.109G>A mutation in 5 patients and homozygous c.235delC mutation in 4 patients. In the 11 patients who carried heterozygous c.235delC variants, four of them carried c.299_300del AT, 2 carried c.109G>A, and 2 carried c.176_191del GCTGCAAG AACGTGTG compound heterozygotes variants, respectively. In our study, a patient who carried a novel heterozygous variant was identified, the nucleotide change is c.88A>G, the protein change is p.Ile30Val, and the characterization of variant is VUS; this patient carried c.919-2A>G and c.1229C>T compound heterozygotes variants of SLC26A4 gene as well.Table 5
Variants frequency inGJB2 gene.
VariantsHomozygousHeterozygousTotalTesting by 127 gene panelc.109G>A51318c.235delC41115c.299_300del AT145c.176_191del GCTGCAAG AACGTGTG022c.139G>T101c.88A>G011Testing by 159 variant testingc.235delC171936c.299_300delAT099c.176_191delGCTGCAAGAACGTGTG055c.139G>T202c.187G>T011c.427C>T101c.428G>A011c.257C>G011In 159 variant testing, eight variants ofGJB2 were identified; the frequency of these variants from high to low were c.235delC, c.299_300delAT, c.176_191del GCTGCAAGAACGTGTG, c.139G>T, c.187G>T, c.257C>G, c.427C>T, and c.428G>A; the check-out times of these variants were 36, 9, 5, 2, 1, 1, 1, and 1, respectively. Nineteen patients carried heterozygous c.235delC variants, five of them carried c.299_300delAT compound heterozygotes variants, three patients carried c.176_191del GCTGCAAGAACGTGTG compound heterozygotes variants, and one carried c.257C>G compound heterozygotes variants. The homozygous c.235delC patients were seventeen. c.299_300delAT heterozygous variant was identified in 9 patients, two of them carried c.176_191delGCTGCAAGAACGTGTG compound heterozygotes variants.


## 3.5.SLC26A4 Variants

In 127 gene panel testing, seventeen variants were identified inSLC26A4 gene, the most frequent variant detected was c.919-2A>G, eleven patients carried heterozygous c.919-2A>G variant, nine of them were compound heterozygotes variants with other SLC26A4 variants, and one homozygous c.919-2A>G was identified. Eight novel variants were identified in our research; c.1339delA, c.1519delT, c.164+1G>C, and c.2000T>C were likely pathogenic variants, while c.1340A>T, c.718C>T, c.765+5G>A, and c.208C>A were VUS variants (Table 6).Table 6
Variants frequency inSLC26A4 gene.
VariantsHomozygousHeterozygousTotalTesting by 127 gene panelc.919-2A>G11112c.2000T>C022c.1087A>C011c.1174A>T011c.1339delA011c.164+1G>C011c.718C>T011c.1226G>A011c.678T>C011c.1519delT011c.1229C>T011c.1340A>T011Testing by 159 variant testingc.919-2A>G51419c.1229C>T055c.2168A>G044c.1174A>T011c.1226G>A011c.1336C>T011c.1343C>A011c.1343C>T011c.1594A>C011c.1692dupA011c.1707+5G >A011c.589G>A011Twelve variants ofSLC26A4 were identified in 159 variant testing; the most common variant was c.919-2A>G same as the 127 gene panel testing, then was c.1229C>T and c.2168A>G. Nineteen patients carried c.919-2A>G variant, fourteen of them were heterozygous while 5 were homozygous.


## 4. Discussion

Hearing scanning combined with molecular genetic testing is an effective approach to diagnose hereditary hearing loss, while because of the high cost of genetic testing, such limited the usage of more comprehensive genetic testing. Through 127 gene panel testing, more genetic etiology of deaf children was identified. In 127 gene panel testing, thirty-two deafness genes were identified in patients. The five most frequent genes wereGJB2 (33 patients), SLC26A4 (16 patients), MYO15A (5 patients), OTOF (5 patients), and CDH23 (5 patients), while only four deafness genes were identified by 159 variant testing, and the frequency of these four genes was 45 (GJB2), 27 (SLC26A4), 2 (GJB3), and 1 (MT-RNR1), respectively. The genetic spectrum of 127 gene testing is consistent with the previous research in China [19]. The diagnostic rate of GJB2 and SLC26A4 is higher in 159 variant testing which may cause more mutant genes identified in 127 gene panel testing. Furthermore, eighty-four variants in 32 genes which were not included in 159 variant testing were identified, and ten variants of these were pathogenic or likely pathogenic. This result demonstrates that even though 159 variant testing contain hotspot deafness variants, it would still misdetect abundant of gene mutations.Broader genetic testing can improve the diagnostic sensitivity and accuracy. By analyzing the patients who received 127 gene panel testing, we found that 25 patients had mutant genes which were not included in 159 variant testing, and nine of them carried pathogenic or likely pathogenic variants. Fifty-one patients carried variants that were not included in 159 variant testing, and thirty-six of them carried pathogenic or likely pathogenic variants. If these patients only received 159 variant testing, the genetic etiology would not be identified. This result revealed that even if the patients pass the traditional gene testing, they still could have genetic causes. Therefore, expanded gene testing should be done to identify the genetic etiology.Moreover, accurate genetic testing can benefit patients in several aspects. Identifying genetic etiology can guide clinicians to determine treatment strategies. In our research, we found several patients carried variants inCDH23, WFS1, and PCDH15 by 127 gene panel testing. As reported in previous research, some patients with above gene mutation may have a poor prognosis of cochlear implantation (CI) [20]. Besides, early genetic diagnostic of hearing loss would provide more information on nature, mode of inheritance, and implications of genetic disorders, which would help individuals and families make informed medical and personal decisions [14, 21]. Furthermore, the relationship between genotype and phenotype of hereditary hearing loss is complicated; finding novel variants would broaden the understanding of hereditary hearing loss.Genetic testing should be optimized according to multiple dimensions. Firstly, the common detection of deafness genes and variants should vary by ethnicity and region. For example, in 127 gene panel testing, c.109G>A was the most frequent variant, then was the c.235delC. In our research, we did not select the ethnicity of the patients, which would influence the frequency of these variants. A previous research in Qinghai (China) has shown that c.109G>A was more common in minority patients, while c.235delC was more common in Han nationality [22]. Secondly, some gene mutation which might be associated with poor prognosis of CI such as POU3F4, TMPRSS3, PJVK, CDH23, and PCDH15 should be added in genetic testing. The full sequence of these genes or at least the exon regions should be detected. If pathological mutations in these genes are missed, it may bring catastrophic consequences to patients and their families.


## 5. Conclusion

Congenital deafness is the number one congenital disease that endangers human health. The identification of the cause through hearing screening combined with genetic testing is an important basis for mechanism research, clinical intervention, and genetic counseling. With the development of sequencing technology, deeper and broader genetic testing by next-generation sequencing are efficient and affordable. More accurate and convenient genetic testing should be designed and recommended. It can increase the detection rate of patients with hereditary deafness, guide clinical treatment strategies, expand new deafness genes or variants, and provide help for future clinical work.



---

*Source: 6151973-2021-07-23.xml*