Figures
Abstract
Background
Mutations in SLC26A4, which encodes pendrin, are a common cause of deafness. SLC26A4 mutations are responsible for Pendred syndrome and non-syndromic enlarged vestibular aqueduct (EVA). The mutation spectrum of SLC26A4 varies widely among ethnic groups. To investigate the incidence of EVA in Chinese population and to provide appropriate genetic testing and counseling to patients with SLC26A4 variants, we conducted a large-scale molecular epidemiological survey of SLC26A4.
Methods
A total of 2352 unrelated non-syndromic hearing loss patients from 27 different regions of China were included. Hot spot regions of SLC26A4, exons 8, 10 and 19 were sequenced. For patients with one allelic variant in the hot spot regions, the other exons were sequenced one by one until two mutant alleles had been identified. Patients with SLC26A4 variants were then examined by temporal bone computed tomography scan for radiological diagnosis of EVA. Ten SLC26A4 variants were cloned for functional study. Confocal microscopy and radioisotope techniques were used to examine the membrane expression of pendrin and transporter function.
Results
Of the 86 types of variants found, 47 have never been reported. The ratio of EVA in the Chinese deaf population was at least 11%, and that in patients of Han ethnicity reached at least 13%. The mutational spectrum and mutation detection rate of SLC26A4 are distinct among both ethnicities and regions of Mainland China. Most of the variants caused retention of pendrin in the intracellular region. All the mutant pendrins showed significantly reduced transport capability.
Conclusion
An overall description of the molecular epidemiological findings of SLC26A4 in China is provided. The functional assessment procedure can be applied to identification of pathogenicity of variants. These findings are valuable for genetic diagnosis, genetic counseling, prenatal testing and pre-implantation diagnosis in EVA families.
Citation: Yuan Y, Guo W, Tang J, Zhang G, Wang G, Han M, et al. (2012) Molecular Epidemiology and Functional Assessment of Novel Allelic Variants of SLC26A4 in Non-Syndromic Hearing Loss Patients with Enlarged Vestibular Aqueduct in China. PLoS ONE 7(11): e49984. https://doi.org/10.1371/journal.pone.0049984
Editor: Iris Schrijver, Stanford University School of Medicine, United States of America
Received: August 3, 2012; Accepted: October 19, 2012; Published: November 21, 2012
Copyright: © 2012 Yuan et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This work was supported by Chinese National Nature Science Foundation Research Grant (30728030,31071109,81230020), State 863 High Technology Research & Development Key Project of China(SQ2010AA0221634001) and grant from the Minister of Science and Technology (2012BAI09B02) to PD, the National Basic Research Program of China (973 Program) (#2012CB967900, #2012CB967901, 2011CBA01000) to SY. Chinese National Nature Science Foundation Research Grant (30801285) and Beijing Nova programme (2009B34) to YY, Chinese National Nature Science Foundation Research Grant (81000415) and Beijing Nova programme (2010B081) to GW, National Institutes of Health grant R01 DC 004696 from the National Institute on Deafness and Other Communication Disorders to JT and DZZH. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
SLC26A4 encodes an anion transporter transmembrane protein, pendrin, which is expressed in the thyroid, kidney, and cochlea [1]. Pendrin is a member of the anion transporter family SLC26, which mediates the exchange of anions including Cl-, HCO3-, OH-, I-, or formate [2]. In thyrocytes, iodide and sodium are brought into the cells via the basolaterally located sodium–iodide symporter. Apically located pendrin seems to be responsible for the efflux of iodide into the follicular lumen [3]. In the kidney, pendrin is suspected to mediate Cl−/HCO3- exchange in the acid–base regulating B- and non-A–non-B-intercalated cells [4]. Similarly, in the inner ear, pendrin is thought to mediate Cl−/HCO3- exchange, and is therefore involved in the conditioning of endolymphatic fluid, presumably due to HCO3- secretion [5]. Malfunction of pendrin leads to Pendred syndrome (PS) and non-syndromic DFNB4 deafness with EVA [6].
Pendred syndrome (PS) is the most common form of syndromic deafness and accounts for about 10% of hereditary hearing impairment [7]. It is an autosomal recessive disorder caused by biallelic mutations in SLC26A4 resulting in hearing loss, enlargement of the vestibular aqueduct (EVA) and iodine organification defect in the thyroid gland [8]. EVA, diagnosed based on the criteria of a greater than 1.5-mm diameter at the midpoint between the common crus and the external aperture, is detected in the ears of patients with PS by computed tomography (CT) and magnetic resonance imaging (MRI) [9]. EVA is the most common form of inner ear malformations associated with prelingual or postlingual sensorineural hearing loss, and is an important feature of PS [10]. EVA may occur alone or in combination with an incomplete partition of the apical turn of the cochlea as part of a Mondini deformity. PS is differentiated from non-syndromic hearing loss with EVA by the presence of goiter, which usually develops later at around the time of puberty. Since environmental and other genetic factors may modulate the effects of SLC26A4 mutations on the development of goiter, the expression of goiter in PS patients is variable and may have incomplete penetrance [11]. Hearing loss in PS and DFNB4 occurs either congenitally or after some mild head injury, showing fluctuating and progressive hearing loss.
To-date, 174 SLC26A4 mutations have been reported [12]. The mutation spectrum varies widely among ethnic groups [8], [11], [13]–[21]. Park and Pryor observed that patients with PS were always associated with two mutant alleles in SLC26A4, consistent with autosomal recessive disorder, whereas patients with non-syndromic hearing loss and EVA might have one or no mutant allele [8], [11]. In a Caucasian non-syndromic EVA cohort, about one-third of the patients had two mutant alleles, one-third had one mutant allele and the remaining one-third had zero [11]. In Japanese and Korean, the proportions of EVA patients with two identified mutant alleles in SLC26A4 are 57% and 81%, respectively [14], [20]. In China, 97.9% of the EVA patients in simplex families were detected with either biallelic or monoallelic mutations, of which 88.4% carried biallelic variants and 9.5% had a monoallelic mutation. Only 2.1% of the Chinese patients with EVA had no mutant SLC26A4 allele detected [17].
SLC26A4 mutations are the second-most common cause of deafness in China. The incidence of non-syndromic EVA accounts for 13.73% of hereditary hearing loss [22]. However, Pendred syndrome is relatively rare in the Chinese population. Molecular epidemiology surveys showed that there are more than 100 SLC26A4 variants in the Chinese hearing loss population, of which at least half are novel. With genetic testing of deafness becoming more popular in China, verification of the pathogenicity of novel variants is essential for both genetic and prenatal diagnosis. However, the two parameters used thus far–(i) low incidence of the mutation in the control population and (ii) substitution of evolutionarily conserved amino acids by the mutation–are not sufficient for clinical genetic diagnosis. To investigate the incidence of EVA in Chinese hearing loss populations and provide appropriate genetic testing and counseling to patients with SLC26A4 variants, we conducted a large-scale molecular epidemiological survey of SLC26A4 in China. In addition, in order to assess the pathogenicity and understand the molecular mechanism(s) underlying EVA and hearing loss, we examined the membrane expression of pendrin and the transporter function of 10 variants that comprised six novel variants and four reported only in Chinese patients.
Materials and Methods
Subjects
A total of 2352 unrelated non-syndromic hearing loss patients from 27 different regions of China were included in this study. To sample as many regions of China as possible, we included the remote northwestern provinces, including Xinjiang (translates as ‘New Territory’), Tibet, and Qinghai, where minorities make up a significant proportion of the local population. We also included samples from the southwestern provinces of Yunnan and Guizhou, where the populations comprise a number of minorities originating from various native tribes. The patients consisted of 1314 males and 1038 females ranging in age from 1.8 to 24 years with an average age of 14.17±3.41. The majority of patients were Han Chinese (1903), followed by minorities in the Southwest region (125), Tibetan (119), Hui (89), minorities in Xinjiang (69), Mongolian (33), Maan (12) and Korean (2). Ethnic subgroup designations were based on permanent residency documentation.
The subjects in this study were from special educational schools within each region. All patients showed moderate to profound bilateral sensorineural hearing impairment on audiograms. Careful medical examinations revealed no clinical features other than hearing impairment. In addition to the 2352 patients, 200 Han individuals with normal hearing (87 females and 113 males, aged 15 to 30 years) were recruited as control subjects from Beijing (Northern) and Jiangsu Province (Eastern), two densely populated regions consisting of 98 to 99% Han Chinese. Auditory function was evaluated by pure-tone air-and bone-conduction threshold audiometry.
Written informed consent was obtained from all subjects or guardians prior to blood sampling and genetic testing. The study protocol including the consent procedure was performed with the approval of the Ethics Committee of the Chinese PLA General Hospital.
Mutational Analysis of SLC26A4
DNA was extracted from peripheral blood leukocytes using a commercially available DNA extraction kit (Watson Biotechnologies Inc., Shanghai, China). The open reading frame (ORF) of SLC26A4 includes 2–20 exons. DNA sequence analyses of SLC26A4 were performed by polymerase chain reaction (PCR) amplification of the coding exons plus approximately 50–100 bp of the flanking intron regions followed by Big Dye sequencing and analysis using an ABI 3100 DNA sequencing machine (ABI, Foster City, USA.) and ABI 3100 Analysis Software v.3.7 NT, according to the manufacturer’s instructions. Because about 90% of Chinese EVA patients carried at least one allelic mutation in exons 8, 19 or 10 of SLC26A4 [17], [18], [23], we designated the three exons as hot spot regions and sequenced these first in all 2352 patients. For the patients with one allelic mutation or variant in exons 8, 19 or 10, the other exons of SLC26A4 were sequenced one by one until two mutant alleles had been identified.
Normal-hearing controls were screened in exons 2–21 of SLC26A4 for mutations or variants by DHPLC followed by sequencing analysis.
CT Scan and MRI
Patients with mutations or variants in SLC26A4 were examined by temporal bone CT scan for diagnosis of inner ear malformation or EVA based on a diameter of >1.5 mm at the midpoint between the common crus and the external aperture. Some patients with EVA were further examined using magnetic resonance imaging (MRI).
Pendrin cDNA and Site-directed Mutagenesis
The coding region of pendrin cDNA was amplified by PCR from the cDNA library of human thyroid tissue and subcloned into the pEGFP-N1 vector (Invitrogen, Carlsbad, CA, USA) using the XhoI and KpnI restriction sites. The pendrin cDNA sequence was verified by nucleotide sequencing, and was identical to a registered sequence (GenBank NM_000441). Oligonucleotide directed mutagenesis was performed to create constructs of 10 mutant plasmids using the QuikChange® Site-Directed Mutagenesis Kit (Stratagene, La Jolla, California, USA), according to the manufacturer’s protocol. The cloning and mutagenic primers used in this study are listed in Table S1.
Cell Culture and Transient Transfection
Human embryonic kidney (HEK) cells were cultured in DMEM solution (Invitrogen, Carlsbad, CA, USA) supplemented with 10% fetal bovine serum. Constructs of pendrin and its mutated plasmids were introduced to the cells using lipofectamine 2000 (Invitrogen, Carlsbad, CA, USA). The amount of DNA used per 35-mm dish was 4 µg, mixed with 10-µl lipofectamine. For radioisotope uptake experiments, the cells were passaged into 24-well plates 24-h before transfection, at 2×105 cells per well. Cells were quantified using a hemocytometer (Fisher Scientific Inc., Pittsburgh, PA). The quantity of DNA used per well was 0.8 µg, added to 1.6 µl lipofectamine. Fluorescence microscopy was used to evaluate membrane targeting during cell-selection experiments.
Transport Function Assessment
Conventional radioisotope techniques were used to measure transport functions from HEK cells transfected with plasmids encoding for pendrin and its mutated products, as described in previous studies [24]–[26]. Pendrin cDNA- and pEGFP-transfected cells were used as the positive and negative controls, respectively. To improve sensitivity, we measured transport function only from those cells that were positively transfected. Fluorescence-based flow cytometry was used for cell sorting. Cells with a fluorescence intensity 4–200 times greater than non-fluorescent cells were selected. In total, 600,000 events (cells) were collected for each sample. To measure [14C] formate uptake, sorted cells in the 24-well culture cluster were first incubated for 30 min in a solution containing 130 mM NaCl, 20 mM HEPES, 5 mM KCl, 5 mM glucose, 2 mM CaCl2 and 1 mM MgCl2 (pH 7.3 and 305 Osm/L). Cells were then incubated at RT for 12 min in a solution containing 140 mM K-gluconate, 20 mM HEPES and 5 mM glucose (pH 7.3 and 305 Osm/L). [14C] formate (Moravek Biochemicals, Inc., Brea, CA) was added to this solution at a concentration of 20 µM. Cells were then washed three times with cold K-gluconate solution without [14C] formate, lysed with 200 µl 0.5 M NaOH, and neutralized with 0.5 M HCl. The lysate was used for liquid scintillation counting to determine the [14C] formate uptake. In each run, three wells were assayed for each plasmid, and the experiments were repeated three times independently. Therefore, the data in each group represented a sample size (n) of 3×3 = 9 trials for each plasmid and control. Student’s t-test was used for statistical analysis.
Results
By sequencing exons 8, 10, and 19 of SLC26A4 in 2352 patients, we identified 342 patients with mutations or variants, of which 116 were homozygotes and 226 were heterozygotes. For the 226 patients with monoallelic mutations or variants in exons 8, 10, or 19 of SLC26A4, 132 were verified as compound heterozygotes by the subsequent whole ORF analysis, and the other 94 were verified as heterozygotes. The total variant detection rate was 14.54% (342/2352), including a homozygous mutation detection rate of 4.93% (116/2352), a compound heterozygous mutation detection rate of 5.44% (128/2352), and a monoallelic mutation detection rate of 4.17% (98/2352). For patients of Han ethnicity, accounting for 80.52% of our sample, the total variant detection rate, homozygous mutation detection rate, compound heterozygous mutation detection rate and monoallelic mutation detection rates were 16.55% (315/1903), 5.68% (108/1903), 6.46% (123/1903) and 4.41% (84/1903), respectively. Eighty-six types of variants in SLC26A4 were found, including 7 nonsense variants, 13 frameshift variants, 48 missense variants, 5 splicing site variants, 6 silent variants, 1 untranslated region (UTR) variant and 6 intron variants (Table 1). The missense variants accounted for 60.76% (48/79, excluding the UTR variant and intron variants) of all the SLC26A4 variants. Of the 86 variants, 47 were novel, which comprised 4 nonsense variants, 6 frameshift variants, 22 missense variants, 2 splicing site variants, 6 silent variants, 1 UTR variant and 6 intron variants (Table 2).
The mutational spectrum of SLC26A4 differed among the various ethnicities. IVS7-2A>G was the most frequent hot-spot mutation in Han, Hui, Mongolian and minorities in the Southwest region, with an allele frequency of 8.01% (377/4704) and a mutant allele frequency of 61.6% (377/612). p.H723R was the second-most frequent hot-spot mutation in Han, Hui, and Mongolian ethnicities with an allele frequency of 1.51% (71/4704) and a mutant allele frequency of 11.6% (71/612). In contrast, neither IVS7-2A>G nor p.H723R were found in the Tibetan and Uygur patients.
The mutation detection rates were distinct in the hearing loss populations of different ethnicities, with the highest detection rate in the Han population (16.55%), followed by the Hui population (13.48%; homozygous mutation detection rate: 4.49%, compound heterozygous mutation detection rate: 3.37% and monoallelic mutation detection rate: 5.62%) and the Mongolian population (12.12%; homozygous mutation detection rate: 6.06%, compound heterozygous mutation detection rate: 0 and monoallelic mutation detection rate: 6.06%). In contrast, the mutation detection rate in the Tibetan population, and minorities in the Southwest region and Uygur, were low at 4.2, 2.4 and 1.45%, respectively.
The SLC26A4 mutation detection rates differed geographically in China. Henan province ranked first with a detection rate of 31.58%. Detection rates in Gansu, Heibei, Hubei, Shaanxi and Heilongjiang province were 26.32, 25, 18.75, 18.37 and 16.67%, respectively. The SLC26A4 mutation detection rate in Guangxi province, the Tibet Autonomous Region and Guizhou province were the lowest at 3.3, 5.08 and 5.4%, respectively. In Guangxi province, no biallelic SLC26A4 mutations were found in the 90 patients enrolled. In the Tibet Autonomous Region, biallelic SLC26A4 mutations were detected only in patients of Han ethnicity. The geographic distribution and the proportion of patients carrying two SLC26A4 mutant alleles in each region are shown in Figure S1. Our results revealed that the incidence of EVA varied geographically in Mainland China. Surprisingly, we found that no Tibetan patients with EVA living in the Tibet plateau were diagnosed either by SLC26A4 genetic testing or by high resolute temporal bone CT scan.
Four IVS7-2A>G heterozygotes, one missense c.2176A>G (p.I726V) heterozygote, one silent variant c.2217A>G (p.Q739Q), one silent variant c.225C>G (p.L75L) and one intron variant (intron4-12T>A) were found in the control group. Although this control population may be too small to reach a final conclusion, the carrier rate of SLC26A4 mutations in China is estimated to be about 2.5%. Polymorphisms in SLC26A4 appear to be rare in the general population in comparison to those in GJB2 [22].
Phenotype
Of the 342 hearing loss patients with mutations or variants in SLC26A4, 320 were examined by temporal CT scan or MRI. The follow-up of the remaining 22 patients (4 with biallelic mutations and 18 with monoallelic mutations) failed. All patients (240/240) with biallelic mutations or biallelic variants (not including silent variants) and 23.75% (19/80) of those with monoallelic mutations or variants in the ORF of SLC26A4 were confirmed to have EVA by temporal bone CT scan or MRI (Figure S2). All patients diagnosed with EVA in our study showed moderate to profound bilateral sensorineural hearing impairment on audiograms (Figure S3). Because none of our patients showed goiter upon physical examination, the perchlorate discharge test was not performed in our study.
Functional Study on the 10 SLC26A4 Variants
We chose 10 SLC26A4 variants for functional study using the heterologous expression system. Of these, c.279T>A (p.S93R), c.665G>T (p.G222V), c.941C>T (p.S314L), c.1517T>G (p.L506R), c.1985G>A (p.C662Y), and c.2326C>G (p.R776G) have not been reported before. c.259G>T (p.D87Y), c.1079C>T (p.A360V), c.1225C>T (p.R409C) and c.1991C>T (p.A664V) have only been detected in ethnic Chinese [27]–[29]. Patients with the above 10 variants also carried an identified mutation, such as IVS7-2 A>G or c.2168 A>G (p.H723R) or c.1174 A>T (p.N392Y) on the other allele, and all were diagnosed with EVA by CT scan. The molecular epidemiological information and sift prediction of the pathogenicity of the ten variants are summarized in Table 3. The sift prediction of seven variants (p.D87Y, p.S93R, p.G222V, p.S314L, p.R409C, p.L506R, p.C662Y) was damaging, while the other three (p.A360V, p.A664V and p.R776G) was tolerated. All of the missense SLC26A4 variants, except p.L506R, p.C662Y and p.A664V, affect amino acids located in highly conserved regions of the SLC26A4 orthologs. However, L506, C662 and A664 are located in conserved regions of most of the SLC26A4 orthologs (Figure S4). None of the 10 SLC26A4 variants was found in the normal-hearing control group.
Membrane targeting is essential for the functional expression of a membrane-bound protein. We questioned whether these mutations could alter the folding of the 10 variants, which could result in reduced or loss of membrane expression. It has been demonstrated that pendrin shows robust membrane-targeted expression in HEK 293 cells. To determine whether these 10 variants were expressed in the membrane, we used confocal microscopy to examine the expression of EGFP and a membrane-based dye, di-8-amino-naphthylethenyl-pyridinium (di-8-ANEPPS) in HEK cells after transfecting the cells with plasmids containing one of the 10 variants tagged by EGFP. Di-8-ANEPPS is a member of the ANEP (aminoaphthyl ethenyl pyridinium) class of membrane potential dyes. Co-localization of EGFP expression and the di-8-ANEPPS dye in the membrane would suggest proper membrane targeting of the proteins [25]. Figure S5 shows representative images of membrane expression patterns of the 10 mutant pendrins. For comparison, membrane expression of the normal pendrin protein is also presented. Although it appeared that p.D87Y, p.S93R, p.G222V, p.R409C, p.L506R, p.C662Y, p.A664V and p.R776G all reached the plasma membrane, their expression levels were reduced when compared to normal pendrin. Other mutant pendrins (p.S314L, p.A360V, p.R409C, p.L506R, p.C662Y and p.R776G) were also localized to the cytosol, indicating retention in the intracellular region.
It has been suggested that the defect in Cl−/HCO3- exchange activities at the apical membrane of inner ear epithelial cells is the key factor that causes EVA and deafness [30]. We evaluated the transport activity of each mutant pendrin and compared it with that of WT-pendrin. A conventional radioisotope technique was used to measure [14C] formate uptake after sorting of EGFP-positive cells by fluorescence-based flow cytometry. Cells transfected with EGFP-vector only were used as a negative control. As shown in Figure S6, all mutant pendrins showed significantly decreased transport function compared to that of normal pendrin (P<0.01). Pfarr et al. reported that the p.R776C mutation does not compromise transporter activity [31]. We showed here that the transporter activity of p.R776G was significantly reduced, although the reduction was less dramatic than in other variants.
Discussion
Mutations in SLC26A4, which encodes pendrin, are a common cause of deafness, and are responsible for both syndromic and non-syndromic hearing loss. The mutation spectrum of SLC26A4 varies widely among ethnic groups [12]. In addition, the prevalence of mutations varies among ethnic groups. Campbell et al. reported p.T416P and IVS8+1G>A as the two most frequent mutations in northern European populations [6], whereas Blons et al. reported a completely different mutation spectrum that was extremely heterogeneous [13]. In a Japanese population, p.H723R accounted for 53% of the mutant alleles, and in a Korean population, p.H723R and the IVS7-2A>G mutation were the most prevalent, accounting for 45.5% of patients with PS or EVA [11], [20]. In China, the IVS7-2A>G mutation was the most common form, accounting for 57.63% of the mutant alleles [17].
All of the above studies focused on patients with EVA or Pendred syndrome. Few large-scale molecular epidemiological surveys on the SLC26A4 gene in hearing loss populations have been reported. Previous studies verified that about 90% of Chinese patients with EVA carried at least one allelic mutation in exons 8, 19 or 10 of SLC26A4 [17], [18], [23]. Because SLC26A4 is a large gene with 20 exons, hot-spot region screening in a molecular epidemiological survey is an economical method. SLC26A4 mutations were detected in nearly 15% of 2352 non-syndromic Chinese hearing impairment patients, with IVS7-2A>G being the most prevalent. By screening the above three exons of SLC26A4, we verified that the ratio of EVA in the Chinese deaf population was at least 11%, and that in the Han ethnicity reached at least 13%. Of our patients with EVA, 23.75% carried only one missense mutation. Whether the missense mutation causes a dominant negative effect and/or specifies a different phenotype is not clear. It is possible that the second mutant allele has not yet been identified due to the location of mutations deep in introns or promoter regions that were not sequenced, intragenic exon deletions, or the involvement of mutations in genes other than SLC26A4, such as FOXI and KCNJ10 (i.e., digenic synergistic mutations). However, none of our patients with one SLC26A4 missense mutation carried any of the reported FOXI or KCNJ10 mutations.
We found that different ethnicities that originated from various regions of China showed distinct SLC26A4 mutational spectra and detection rates, which may be explained by the genetic heterogeneity in various ethnicities. For example, the mutation detection rates in the Tibetan population and minorities in the Southwest region and Uygur populations were significantly lower than that in the Han, Hui and Mongolian populations. Although IVS7-2A>G and p.H723R were hot mutations in Han, Hui and Mongolian patients, they were not found in the Tibetan and Uygur patients. Our results also revealed that the SLC26A4 detection rate and incidence of EVA varied geographically. The incidence of inner ear malformation in the Tibetans living in the Tibet plateau was almost as high as that in Chinese Han patients, but the types of malformation in the Tibetan population differed greatly [32]. The most-common inner ear deformity, EVA, is rare in the Chinese Tibetan hearing-loss population. Hypoxia may be one of the causes of development of inner ear deformity, but further studies are required to determine the genetic etiology of hearing loss in Tibetan patients. We also found that the detection rate of another common gene, GJB2, was significantly lower in Guangxi province (6.05%) than the average rate in China (21.01%) [33]. These results may be explained by the influence of regional or environmental on the etiology of hearing loss in Guangxi.
A striking achievement of this study is the establishment of a new strategy for detecting SLC26A4 mutations prior to temporal bone CT scan for identifying patients with EVA. In China, the cost of a temporal CT scan is 200 to 300 RMB (approximately 3 days’ salary of an average Chinese). Because of the relatively high cost, it is not possible to perform a CT scan in every hearing loss patient for diagnosis of EVA. Because 97.9% of Chinese patients with EVA carry an SLC26A4 mutation [17], SLC26A4 mutations in hearing loss patients indicate a high possibility of EVA. This model has the unique advantage of identifying patients with EVA in an epidemiologic study in large-scale deaf populations.
The discovery of 47 novel variants enriches the mutational spectrum of SLC26A4. Functional tests of mutated pendrin allelic variants found in patients with Pendred syndrome or non-syndromic EVA revealed that the pathology is linked to a reduction or a loss of function in the ion transport activity of pendrin [34], [35]. Subcellular localization studies revealed that the mutated proteins are often retained in subcellular compartments and are unable to reach the plasma membrane [36], [37]. These mutations most likely cause protein misfolding, followed by impaired trafficking and subsequent degradation. Furthermore, the transport function of mutated proteins that do reach the membrane is impaired [34], [38], [39]. Besides clinical and radiological assessments, the functional evaluation of pendrin mutations is essential for correct diagnosis of Pendred syndrome and non-syndromic EVA [40]. Indeed, the clinical condition of pseudo-Pendred syndromes [41]–[43] and the high incidence of functional pendrin polymorphisms in some populations [30] could lead to an incorrect assignment of the pathological conditions in cases where proper functional characterizations are not possible. Whereas a broad spectrum of mutations is found in Caucasian populations, and Israeli and Palestinian populations have been functionally characterized [44], [45], little about the variants that occur in Chinese is known.
We used confocal microscopy and radioisotope techniques to examine the membrane expression and transporter function of 10 mutant proteins in human embryonic kidney cells. We showed that all 10 mutations resulted in a significant reduction in membrane expression of the proteins. Measurement using a radioisotope technique showed that their transport activity was significantly reduced. Our results suggest that the pathogenesis underlying the hearing loss associated with these mutations is the result of reduced membrane expression and decreased transport capability of the proteins.
Pera et al. [30] suggested that certain biochemical parameters explained some functional impairment. For example, transport function was impaired if an amino acid bearing a fixed charge was missing or introduced (i.e., aspartic acid (D) with fixed negative charge (K) or glutamic acid (EK), lysine (K) with positive fixed charge (C) or arginine (RC), and histidine (HC) with the positive charge of histidine, depending on the pH). Furthermore, the loss or inclusion of a proline (P; proline acts as a structural disruptor of regular secondary structures such as α-helices or β-sheets) in the SLC26A4 sequence is detrimental to its transport function. However, Pfarr et al. [31] described p.R776C as an exception, with an arginine located on the extreme C-terminus of SLC26A4 (a 780 amino acid protein). It is thought that the N- and C-termini of a protein are not structural and the mutations occurring in these areas probably have little or no functional impact [30]. However, with regard to the EVA phenotype, the p.R776G mutation in Chinese was found at very low allele frequency in the control group.
In conclusion, we provide an overview of the molecular epidemiological status of SLC26A4 in China. We showed that the prevalence of EVA in the Chinese deaf population was at least 11%, which reached at least 13% in the Han majority population. The mutational spectrum and mutation detection rate of SLC26A4 are distinct among both ethnicities and regions of mainland China. Our results suggest that the pathogenesis underlying hearing loss in patients with EVA is the result of reduced membrane expression and the decreased transport activity of the mutant pendrins. The functional assessment procedure can be applied to genetic diagnosis of EVA patients and prenatal diagnosis for their families. These findings are valuable for genetic counseling and risk assessment for families with EVA in China.
Supporting Information
Figure S1.
Geographic distribution and the proportion of patients carrying two SLC26A4 mutant alleles in each region studied. NE: not examined
https://doi.org/10.1371/journal.pone.0049984.s001
(TIF)
Figure S2.
Examples of CT and MRI images of normal children and children with SLC26A4 mutations. A: CT of temporal bone of a normal child. Arrows show the normal vestibular aquaduct. B: CT of a patient with EVA. Arrows mark the enlarged vestibular aquaduct. C: MRI of child with normal inner ear. Arrows indicate the normal endolymphatic sac. D: MRI of a patient with EVA. Arrows mark the hydrop of endolymphatic sac.
https://doi.org/10.1371/journal.pone.0049984.s002
(TIF)
Figure S3.
Representative pure-tone audiograms of children with SLC26A4 mutations. A: c.259G>T/c.2168A>G. B: c.941C>T/c.1174A>G. C: c.1517T>G/IVS7-2A>G. D: c.1991C>T/IVS7-2A>G
https://doi.org/10.1371/journal.pone.0049984.s003
(TIF)
Figure S4.
Protein sequences of 13 SLC26A4 orthologs, including: (1) Homo sapiens (NP_000432.1), (2) Nomascus leucogenys (XP_003268184.1), (3) Pan troglodytes (XP_519308.2), (4) Macaca mulatta (XP_001094049.1), (5) Callithrix jacchus (XP_002751785.1), (6) Sus scrofa (XP_003357559.1), (7) Canis lupus familiaris (XP_540382.3), (8) Rattus norvegicus (NP_062087.1), (9) Mus musculus (NP_035997.1), (10) Oryctolagus cuniculus (XP_002712085.1), (11) Loxodonta africana (XP_003407255.1), (12) Xenopus (Silurana) tropicalis (NP_001107135.1), and (13) Monodelphis domestica (XP_001363598.1).
https://doi.org/10.1371/journal.pone.0049984.s004
(TIF)
Figure S5.
Heterologous expression of wild-type and mutant pendrins. HEK cells were transfected with each plasmid of the 10 variants tagged by EGFP. A membrane-based dye (orange color), di-8-ANEPPS, was also added in the solution before confocal microscopy. Co-localization of EGFP expression and the di-8-ANEPPS dye in the membrane would suggest proper membrane targeting of the proteins.
https://doi.org/10.1371/journal.pone.0049984.s005
(TIF)
Figure S6.
Transport activity of wild-type pendrin and 10 mutant pendrins.
https://doi.org/10.1371/journal.pone.0049984.s006
(TIF)
Methods S1.
SLC26A4 Site-Directed Mutagenesis Primers. The mRNA from normal human adenoid specimens was extracted, and cDNA was synthesized. PCR was done to catch the target gene. Primers were as follows: XhoI -F:5′- CCGCTCGAGATGGCAGCGCCAGGCGGCAG -3′ KpnI -R:5′- GGGGTACCGTGGATGCAAGTGTACGCATAG -3′ XhoI -F contains protective base and XhoI restriction sites, its 5 'end partial sequence for PCR amplify target gene. KpnI -R contains protective base and KpnI restriction sites,its 3 'end partial sequence for PCR amplify target gene. The entire ORF of the human SLC26A4 gene was cloned into the expression vector pEGFP-N1 (Invitrogen, Carlsbad, CA, USA). The point mutation-specific primers were synthesized as shown in Table S1.
https://doi.org/10.1371/journal.pone.0049984.s007
(DOCX)
Table S1.
SLC26A4 Site-Directed Mutation Primers.
https://doi.org/10.1371/journal.pone.0049984.s008
(DOCX)
Acknowledgments
We thank all the normal-hearing subjects who collaborated in this study. We thank Dr.Dongyi Han and Dr.Shasha Huang who gave us detailed counseling on the genotype and phenotype correlation of SLC26A4.
Author Contributions
Conceived and designed the experiments: PD. Performed the experiments: YY WWG JT GZZ GW MYH. Analyzed the data: YY JT DZZH SY XZ. Contributed reagents/materials/analysis tools: PD YY DZZH. Wrote the paper: YY DZZH.
References
- 1. Everett LA, Morsli H, Wu DK, Green ED (1999) Expression pattern of the mouse ortholog of the Pendred's syndrome gene (Pds) suggests a key role for pendrin in the inner ear. Proc Natl Acad Sci U S A 96: 9727–9732.
- 2. Mount DB, Romero MF (2004) The SLC26 gene family of multifunctional anion exchangers. Pflugers Arch 447: 710–721.
- 3. Yoshida A, Hisatome I, Taniguchi S, Sasaki N, Yamamoto Y, et al. (2004) Mechanism of iodide/chloride exchange by pendrin. Endocrinology 145: 4301–4308.
- 4. Kim YH, Kwon TH, Frische S, Kim J, Tisher CC, et al. (2002) Immunocytochemical localization of pendrin in intercalated cell subtypes in rat and mouse kidney. Am J Physiol Renal Physiol 283: F744–754.
- 5. Wangemann P, Nakaya K, Wu T, Maganti RJ, Itza EM, et al. (2007) Loss of cochlear HCO3- secretion causes deafness via endolymphatic acidification and inhibition of Ca2+ reabsorption in a Pendred syndrome mouse model. Am J Physiol Renal Physiol 292: F1345–1353.
- 6. Campbell C, Cucci RA, Prasad S, Green GE, Edeal JB, et al. (2001) Pendred syndrome, DFNB4, and PDS/SLC26A4 identification of eight novel mutations and possible genotype-phenotype correlations. Hum Mutat 17: 403–411.
- 7. Reardon W, Coffey R, Phelps PD, Luxon LM, Stephens D, et al. (1997) Pendred syndrome–100 years of underascertainment? QJM 90: 443–447.
- 8. Pryor SP, Madeo AC, Reynolds JC, Sarlis NJ, Arnos KS, et al. (2005) SLC26A4/PDS genotype-phenotype correlation in hearing loss with enlargement of the vestibular aqueduct (EVA): evidence that Pendred syndrome and non-syndromic EVA are distinct clinical and genetic entities. J Med Genet 42: 159–165.
- 9. Phelps PD, Coffey RA, Trembath RC, Luxon LM, Grossman AB, et al. (1998) Radiological malformations of the ear in Pendred syndrome. Clin Radiol 53: 268–273.
- 10. Valvassori GE, Clemis JD (1978) The large vestibular aqueduct syndrome. Laryngoscope 88: 723–728.
- 11. Park HJ, Shaukat S, Liu XZ, Hahn SH, Naz S, et al. (2003) Origins and frequencies of SLC26A4 (PDS) mutations in east and south Asians: global implications for the epidemiology of deafness. J Med Genet 40: 242–248.
- 12.
Chang E, Kölln KA, Nishimura C, Fischer S, Smith RJH Mutation location within SLC26A4. [http://www.healthcare.uiowa.edu/labs/pendredandbor/slcMutations.htm]
- 13. Blons H, Feldmann D, Duval V, Messaz O, Denoyelle F, et al. (2004) Screening of SLC26A4 (PDS) gene in Pendred's syndrome: a large spectrum of mutations in France and phenotypic heterogeneity. Clin Genet 66: 333–340.
- 14. Park HJ, Lee SJ, Jin HS, Lee JO, Go SH, et al. (2005) Genetic basis of hearing loss associated with enlarged vestibular aqueducts in Koreans. Clin Genet 67: 160–165.
- 15. Prasad S, Kolln KA, Cucci RA, Trembath RC, Van Camp G, et al. (2004) Pendred syndrome and DFNB4-mutation screening of SLC26A4 by denaturing high-performance liquid chromatography and the identification of eleven novel mutations. Am J Med Genet A 124A: 1–9.
- 16. Albert S, Blons H, Jonard L, Feldmann D, Chauvin P, et al. (2006) SLC26A4 gene is frequently involved in nonsyndromic hearing impairment with enlarged vestibular aqueduct in Caucasian populations. Eur J Hum Genet 14: 773–779.
- 17. Wang QJ, Zhao YL, Rao SQ, Guo YF, Yuan H, et al. (2007) A distinct spectrum of SLC26A4 mutations in patients with enlarged vestibular aqueduct in China. Clin Genet 72: 245–254.
- 18. Dai P, Yuan Y, Huang D, Zhu X, Yu F, et al. (2008) Molecular etiology of hearing impairment in Inner Mongolia: mutations in SLC26A4 gene and relevant phenotype analysis. J Transl Med 6: 74.
- 19. Lopez-Bigas N, Melchionda S, de Cid R, Grifa A, Zelante L, et al. (2001) Identification of five new mutations of PDS/SLC26A4 in Mediterranean families with hearing impairment. Hum Mutat 18: 548.
- 20. Tsukamoto K, Suzuki H, Harada D, Namba A, Abe S, et al. (2003) Distribution and frequencies of PDS (SLC26A4) mutations in Pendred syndrome and nonsyndromic hearing loss associated with enlarged vestibular aqueduct: a unique spectrum of mutations in Japanese. Eur J Hum Genet 11: 916–922.
- 21. Wu CC, Yeh TH, Chen PJ, Hsu CJ (2005) Prevalent SLC26A4 mutations in patients with enlarged vestibular aqueduct and/or Mondini dysplasia: a unique spectrum of mutations in Taiwan, including a frequent founder mutation. Laryngoscope 115: 1060–1064.
- 22. Yuan Y, You Y, Huang D, Cui J, Wang Y, et al. (2009) Comprehensive molecular etiology analysis of nonsyndromic hearing impairment from typical areas in China. J Transl Med 7: 79.
- 23. Yuan YY WG, Huang DL, Kang DY, Dai P (2010) Analysis of the SLC26A4 gene hot spot mutation regions in Chinese Enlarged Vestibular Aqueduct Patients. Chinese Journal of Otology 8: 292–295.
- 24. Bai JP, Surguchev A, Montoya S, Aronson PS, Santos-Sacchi J, et al. (2009) Prestin's anion transport and voltage-sensing capabilities are independent. Biophys J 96: 3179–3186.
- 25. Tang J, Pecka JL, Tan X, Beisel KW, He DZ (2011) Engineered pendrin protein, an anion transporter and molecular motor. J Biol Chem 286: 31014–31021.
- 26. Tan X, Pecka JL, Tang J, Okoruwa OE, Zhang Q, et al. (2011) From zebrafish to mammal: functional evolution of prestin, the motor protein of cochlear outer hair cells. J Neurophysiol 105: 36–44.
- 27. Lai CC, Chiu CY, Shiao AS, Tso YC, Wu YC, et al. (2007) Analysis of the SLC26A4 gene in patients with Pendred syndrome in Taiwan. Metabolism 56: 1279–1284.
- 28. Zhao J, Wu LQ, Feng Y, Pan Q, Zhao K, et al. (2009) [Mutational screening of the SLC26A4 gene in patients with nonsyndromic hearing loss by denaturing high-performance liquid chromatography]. Zhonghua Yi Xue Yi Chuan Xue Za Zhi 26: 21–25.
- 29. Huang S, Han D, Yuan Y, Wang G, Kang D, et al. (2011) Extremely discrepant mutation spectrum of SLC26A4 between Chinese patients with isolated Mondini deformity and enlarged vestibular aqueduct. J Transl Med 9: 167.
- 30. Pera A, Dossena S, Rodighiero S, Gandia M, Botta G, et al. (2008) Functional assessment of allelic variants in the SLC26A4 gene involved in Pendred syndrome and nonsyndromic EVA. Proc Natl Acad Sci U S A 105: 18608–18613.
- 31. Pfarr N, Borck G, Turk A, Napiontek U, Keilmann A, et al. (2006) Goitrous congenital hypothyroidism and hearing impairment associated with mutations in the TPO and SLC26A4/PDS genes. J Clin Endocrinol Metab 91: 2678–2681.
- 32. Yuan Y, Zhang X, Huang S, Zuo L, Zhang G, et al. (2012) Common molecular etiologies are rare in nonsyndromic Tibetan Chinese patients with hearing impairment. PLoS One 7: e30720.
- 33. Dai P, Yu F, Han B, Liu X, Wang G, et al. (2009) GJB2 mutation spectrum in 2,063 Chinese patients with nonsyndromic hearing impairment. J Transl Med 7: 26.
- 34. Taylor JP, Metcalfe RA, Watson PF, Weetman AP, Trembath RC (2002) Mutations of the PDS gene, encoding pendrin, are associated with protein mislocalization and loss of iodide efflux: implications for thyroid dysfunction in Pendred syndrome. J Clin Endocrinol Metab 87: 1778–1784.
- 35. Yoon JS, Park HJ, Yoo SY, Namkung W, Jo MJ, et al. (2008) Heterogeneity in the processing defect of SLC26A4 mutants. J Med Genet 45: 411–419.
- 36. Rotman-Pikielny P, Hirschberg K, Maruvada P, Suzuki K, Royaux IE, et al. (2002) Retention of pendrin in the endoplasmic reticulum is a major mechanism for Pendred syndrome. Hum Mol Genet 11: 2625–2633.
- 37. Dai P, Stewart AK, Chebib F, Hsu A, Rozenfeld J, et al. (2009) Distinct and novel SLC26A4/Pendrin mutations in Chinese and U.S. patients with nonsyndromic hearing loss. Physiol Genomics 38: 281–290.
- 38. Choi BY, Stewart AK, Madeo AC, Pryor SP, Lenhard S, et al. (2009) Hypo-functional SLC26A4 variants associated with nonsyndromic hearing loss and enlargement of the vestibular aqueduct: genotype-phenotype correlation or coincidental polymorphisms? Hum Mutat 30: 599–608.
- 39. Dossena S, Vezzoli V, Cerutti N, Bazzini C, Tosco M, et al. (2006) Functional characterization of wild-type and a mutated form of SLC26A4 identified in a patient with Pendred syndrome. Cell Physiol Biochem 17: 245–256.
- 40. Fugazzola L, Mannavola D, Cerutti N, Maghnie M, Pagella F, et al. (2000) Molecular analysis of the Pendred's syndrome gene and magnetic resonance imaging studies of the inner ear are essential for the diagnosis of true Pendred's syndrome. J Clin Endocrinol Metab 85: 2469–2475.
- 41. Fugazzola L, Cerutti N, Mannavola D, Crino A, Cassio A, et al. (2002) Differential diagnosis between Pendred and pseudo-Pendred syndromes: clinical, radiologic, and molecular studies. Pediatr Res 51: 479–484.
- 42. Kara C, Kilic M, Ucakturk A, Aydin M (2010) Congenital goitrous hypothyroidism, deafness and iodide organification defect in four siblings: Pendred or pseudo-Pendred syndrome? J Clin Res Pediatr Endocrinol 2: 81–84.
- 43. Davis N, Lunardi C, Shield JP (2006) Sensori-neural deafness and hypothyroidism: autoimmunity causing 'pseudo-Pendred syndrome'. Horm Res 65: 267–268.
- 44. Dossena S, Rodighiero S, Vezzoli V, Nofziger C, Salvioni E, et al. (2009) Functional characterization of wild-type and mutated pendrin (SLC26A4), the anion transporter involved in Pendred syndrome. J Mol Endocrinol 43: 93–103.
- 45. Dossena S, Nofziger C, Brownstein Z, Kanaan M, Avraham KB, et al. (2011) Functional characterization of pendrin mutations found in the Israeli and Palestinian populations. Cell Physiol Biochem 28: 477–484.