Whole exome sequencing in adult-onset hearing loss reveals a high load of predicted pathogenic variants in known deafness-associated genes and identifies new candidate genes
BMC Medical Genomics volume 11, Article number: 77 (2018)
Deafness is a highly heterogenous disorder with over 100 genes known to underlie human non-syndromic hearing impairment. However, many more remain undiscovered, particularly those involved in the most common form of deafness: adult-onset progressive hearing loss. Despite several genome-wide association studies of adult hearing status, it remains unclear whether the genetic architecture of this common sensory loss consists of multiple rare variants each with large effect size or many common susceptibility variants each with small to medium effects. As next generation sequencing is now being utilised in clinical diagnosis, our aim was to explore the viability of diagnosing the genetic cause of hearing loss using whole exome sequencing in individual subjects as in a clinical setting.
We performed exome sequencing of thirty patients selected for distinct phenotypic sub-types from well-characterised cohorts of 1479 people with adult-onset hearing loss.
Every individual carried predicted pathogenic variants in at least ten deafness-associated genes; similar findings were obtained from an analysis of the 1000 Genomes Project data unselected for hearing status. We have identified putative causal variants in known deafness genes and several novel candidate genes, including NEDD4 and NEFH that were mutated in multiple individuals.
The high frequency of predicted-pathogenic variants detected in known deafness-associated genes was unexpected and has significant implications for current diagnostic sequencing in deafness. Our findings suggest that in a clinic setting, efforts should be made to a) confirm key sequence results by Sanger sequencing, b) assess segregations of variants and phenotypes within the family if at all possible, and c) use caution in applying current pathogenicity prediction algorithms for diagnostic purposes. We conclude that there may be a high number of pathogenic variants affecting hearing in the ageing population, including many in known deafness-associated genes. Our findings of frequent predicted-pathogenic variants in both our hearing-impaired sample and in the larger 1000 Genomes Project sample unselected for auditory function suggests that the reference population for interpreting variants for this very common disorder should be a population of people with good hearing for their age rather than an unselected population.
Hearing loss is one of the most common sensory deficits in the human population, and it has a strong genetic component . However, although more than 140 human non-syndromic hearing impairment loci have been mapped and over 100 genes identified, most underlie childhood deafness or early-onset hearing loss. The vast majority of genes involved in hearing remain unknown, including those associated with adult-onset, age-related progressive hearing loss. Age-related hearing loss (ARHL) affects 1 in 3 people over the age of sixty, often leading to social isolation and depression, is associated with subsequent cognitive decline [2,3,4] and a predictor of dementia . The heritability of ARHL has been estimated to be between 30 and 50%, similar to other common complex disorders [1, 6, 7].
Although several ARHL genome-wide association studies (GWAS) have been carried out [8,9,10,11] and promising candidate genes identified, such as SIK3 and ESRRG [10, 11], only five loci have been associated with hearing status at the genome-wide significance level: GRM7 , PCDH20 and SLC28A3 , and ISG20 or ACAN and TRIOBP . Furthermore, single genes can underlie progressive hearing loss with post-lingual onset, including in middle-age, particularly genes underlying dominantly-inherited deafness [15, 16]. Thus, adult-onset hearing loss may result from either rare Mendelian gene variants with large effect size or multiple variants each making a small contribution to hearing loss. It is also unclear whether these variants are in novel genes involved in maintenance of auditory pathways or whether they are milder variants of the same genes that are mutated in congenital deafness.
Here, we have taken a more in-depth approach than GWAS, using whole-exome sequencing (WES) to study thirty patients carefully selected from a total sample of 1479 patients with a variety of adult-onset hearing loss phenotypes to represent the mixed phenotypes and varied genetic aetiology that might be present in a clinical scenario, targeting specific sub-phenotypes to maximise power to detect shared variants. Our aim was to establish to what extent exome sequencing is an effective and appropriate tool for genetic diagnosis of hearing loss in a clinic setting, where there is usually only a single adult patient involved and family members are not available for segregation analysis. Whole exome and genome sequencing are beginning to be used in this scenario for diagnosis of adult-onset hearing loss with the clinician faced with challenges in evaluating the candidate variants identified. Our results demonstrate the value of targeting well-characterised phenotypic subtypes and cross-species data comparison in exome sequencing analysis, and highlight issues which need to be considered in interpreting genetic variants of unknown pathogenicity in current genetic diagnosis and gene discovery studies, in particular the finding that many individuals have multiple predicted-pathogenic variants in different genes known to underlie deafness.
Recruitment of patients
Twenty patients (seven males, thirteen females) with non-syndromic sensorineural adult-onset hearing loss (self-reported age of onset between 20 and 50) were selected from a larger group of 700 patients recruited from the adult hearing aid clinic at the Royal National Throat Nose and Ear Hospital, London, U.K. (described in ). The twenty were chosen based on a family history of hearing loss and an age of onset in middle age. Air and bone conduction thresholds at 0.25, 0.5, 1, 2, 4 and 8 kHz and 0.5, 1, 2 and 4 kHz, respectively were measured with masking as indicated according to BSA Recommended Procedures  (Additional file 1: Figs. S1, and Additional file 2: Figs. S2). Ethical approval for this project was granted from the Royal Free Local Research Ethics Committee (reference 6202).
A second group of 10 older individuals were selected from the 779 people in the database of the longitudinal study of ARHL being conducted at the Medical University of South Carolina (MUSC) since 1987 (described in ). These individuals were selected for the current analyses on the basis of age (> 60 and < 79 years), negative or limited self-reported occupational and recreational noise history, available DNA samples, and audiometric phenotype (five metabolic: 4 females, 1 male; mean age 69.8 years and five sensory: 2 females, 3 males; mean age 68.3 years) (Additional file 3: Fig. S3). The protocols for this study were approved by the Institutional Review Board at MUSC.
All patients were recruited by written informed consent.
DNA was submitted for WES using either the Agilent SureSelect Human All Exon V3 kit or a custom library designed by Agilent for human whole exome sequence capture (which predates the SureSelect kit). Sequencing was carried out on either the Illumina Genome Analyzer IIx or Illumina HiSeq 2000 platform as paired-end 54 bp, 75 bp or 76 bp reads according to the manufacturer’s protocol.
The data were filtered on quality, Minor Allele Frequency (MAF) in the ExAC non-Finnish European population , most severe consequence and predicted pathogenicity using the pipeline shown in Fig. 1 (Fig. 1, Additional file 4: Table S1). Where no data were recorded in ExAC but a MAF was available from the 1000 Genomes European population , e.g. for rs1813100, this was used. The impact of filtering on variant number is shown in Additional file 5: Fig. S4. Following filtering, variants were examined both individually and by gene (using only genes and known miRNAs present in Ensembl, accessed May 2016), and the gene lists analysed to find those mutated genes shared between individuals.
1000 Genomes project data
Genotype and annotation data from the 1000 Genomes project [21, 24] were used to create files of gene variants from each of the 2504 individuals sequenced in the project. Variants were processed using the same pipeline (Fig. 1), with the exception of the quality filter (quality of variant call ≥30, read depth ≥ 20, mapping quality ≥45), which was omitted because these variant calls have already been validated and filtered .
Results and discussion
In order to determine the possibility of using WES to reveal the genetic basis of adult onset hearing loss in a typical clinical scenario, we identified 30 individuals with different phenotypic sub-types based on family history, age of onset and audiogram shape from large well characterised cohorts. Twenty patients with family histories were subdivided into ten with probable dominant hearing loss and ten with presumed recessive hearing loss. Ten older adults without a family history of hearing loss were subdivided into five people with a metabolic phenotype of ARHL and five with a sensory phenotype of ARHL based on audiogram shape . These sub-groups are henceforth referred to as Dominant, Recessive, Metabolic and Sensory respectively.
Because hearing loss is common, we selected a low-stringency allele frequency filter of 10%, to detect common risk variants in our phenotypic sub-groups. We compiled a list of known deafness genes comprising all genes listed in the Hereditary Hearing Loss Homepage  plus the human orthologues of all the mouse deafness genes from the Hereditary Hearing Impairment in Mice website (described in ); 357 genes in total (Additional file 7: Table S3). Mutant mice continue to be valuable tools for discovering genes required for hearing , many of which have subsequently been shown to also underlie deafness in humans, such as WBP2  and MIR96 [27, 28].
To identify potential false positives due to platform errors we utilised a list of 507 such genes described for the Illumina Genome Analyser IIx by Fuentes Fajardo et al.  (Additional file 8: Table S4). We have not excluded these genes from our analysis but have marked them as “candidates for exclusion” where they occur.
First analysis: Common variants
Here, common variants that were predicted to be pathogenic (≤10% MAF in the non-Finnish European population) were analysed to find genes common within each sub-group.
Known deafness genes in all four sub-groups combined
We first examined known deafness genes in all thirty people together and found that every person had at least ten known deafness genes with one or more predicted pathogenic variants (Fig. 2A, pale blue bars). The same analysis on data from the 1000 Genomes project, which includes exome sequences from 2504 people with unknown auditory function, produced a very similar distribution of predicted pathogenic variants to the 30 people with adult-onset hearing loss (Fig. 2B). As this high frequency of predicted pathogenic variants was an unexpected finding, we asked if the same pattern of distribution was present in genes known to be involved in retinal disease, another sensory deficit with a large number of single genes known to be involved. We repeated the analysis using 265 retinal disease genes (from RetNet [30, 31], Additional file 9: Table S5) instead of our list of deafness genes. Again, we found a very similar distribution of predicted pathogenic variants in both our thirty patients with hearing loss and in the 2504 samples from the 1000 Genomes project (Fig. 2C, D). We then looked at human dominant deafness genes only (n = 33, Additional file 7: Table S3), where pathogenic variants would be expected to have an effect even when only present in one allele. We found once again a similar distribution of variants in our thirty patients and the 1000 Genomes samples (Fig. 2E, F) and 27 of our 30 patients had at least one predicted pathogenic variant in a dominant deafness gene (Fig. 2C).
To test whether this finding was due to having too many common variants, we repeated the analysis using a 1% MAF cut-off (instead of 10%), and found the same pattern (Fig. 2, pink bars). We then tested a more stringent filter for pathogenicity, selecting only those variants predicted to be pathogenic by all the predictors used instead of the majority of predictors. The distributions were again similar for all groups of people and sets of disease-associated genes with both MAF filters (Fig. 2, dark blue and orange bars for MAF < 10% and MAF < 1% respectively). Finally, we looked for only loss of function variants (loss of transcript, essential splice site, frameshift indel, stop gain or start loss, see Fig. 1), and again observed a similar distribution (Fig. 2, purple and red bars for MAF < 10% and MAF < 1% respectively). Two of our patients and 152 people from the 1000 Genomes project had a loss of function variant in a known human dominant deafness gene (Fig. 2E, F). The two patients were both from the dominant sub-group; one had a single base pair insertion (causing a frameshift) in MYO6 and the other had a nonsense mutation in HOMER2.
Genes common to individual sub-groups
We then looked at individual genes with predicted pathogenic variants (MAF < 10%, predicted by a majority of predictors) present in all members of each sub-group. No known deafness genes were common to any of the sub-groups. We found four other such genes: TTN (a candidate for exclusion) bore variants in all ten individuals with recessive hearing loss, MON1B was mutated in all five individuals with metabolic deafness (but see later), and the five people with sensory hearing loss all had variants in NEDD4 and ZAN (Table 1).
Implications for WES analyses in common disease
Our common variant analysis has highlighted several potential problems with detecting common variants of small effect size in a common disorder. Although each individual had many deafness genes with predicted pathogenic variants, the overall spread didn’t look any different to that observed in the 1000 Genomes dataset  (Fig. 2). However, the 1000 Genomes data do not exclude people with hearing loss, particularly adult-onset hearing loss which may not be evident at the time of sampling an individual. This highlights the need for good controls in this type of analysis, in both clinics and research; suitable controls in this case might be older adults with good hearing typical of a younger adult. Other exome sequencing projects have reported similar results for age-related macular degeneration  and in the ExAC data an abundance of rare, functional variants were reported in many disease genes . It has been hypothesised that both false pathogenicity reports and incomplete penetrance contribute to this over-reporting, but whatever the reason, putative causative variants must be treated with caution until proof of pathogenicity has been obtained, preferably by functional studies and linkage analysis.
Another factor to consider is whether variants which are rare in the most relevant ethnic populations and passed the MAF filter but are more prevalent in other populations should be retained; 1234 of the 10,482 unique variants in the thirty patients which passed all our filtering steps were present in other ExAC populations at a MAF of more than 10%. To give one extreme example, it is worth considering whether a variant which is present in 60% of the African ExAC population should be included even if it is very rare (0.3%) in the non-Finnish European population, as for a T > C missense variant in LPP (chr3:g.188327555 T > C).
Second analysis: Very rare variants
Since the pattern of common variants in our thirty patients with hearing loss did not differ from that observed in the data from the 1000 Genomes project, we then filtered for very rare variants with a minor allele frequency of < 0.5%, predicted to be deleterious to protein function by a majority of predictors to identify any likely causal variants. This is the recommended MAF filter for recessive hearing loss (the recommended MAF for dominant hearing loss is 0.05%) . No individuals from any of the groups, including those with a strong family history consistent with recessive inheritance, were homozygous for a variant in a known deafness gene. Apparently-homozygous variants were found in 29 other genes in only one of the individuals, and two genes bore homozygous variants in multiple people: SIRPA (a candidate for exclusion) and ZAN (Additional file 10: Table S6).
Many known deafness genes were mutated in multiple samples in the heterozygous state (Table 2); for example, 7 individuals bore variants in GPR98. One individual from the dominant group had five predicted variants in WFS1 (Fig. 3), and many people bore more than one variant in the same gene, including three individuals from the recessive group who bore two heterozygous variants in a deafness gene (Additional file 11: Table S7). These variants might explain the hearing loss seen in these three patients, but with one exception, without segregation analysis it isn’t possible to tell whether the variants are on the same chromosome or were inherited one from each parent. The exception is the two variants in GPR98 in sample 11,813, which are close enough to fall within the same sequencing read, and can be confirmed to originate from the same chromosome. Furthermore, several very rare variants were also found to be present in more than one individual (Additional file 12: Tables S8 and Additional file 13: Table S9).
When we examined heterozygous variants in all genes, we found several genes bearing rare variants in six or more people (Additional file 13: Table S9). One of particular interest was NEFH, which encodes the neurofilament heavy chain and is strongly expressed in the spiral ganglion neurons innervating the inner hair cells .
We also examined the eight strongest candidate genes (PCDH20, GRM7, ESRRG, SIK3, SLC28A3, TRIOBP, ISG20 and ACAN) linked to adult hearing ability or ARHL in GWAS [10,11,12,13,14]. Five rare variants were found: a novel deletion predicted to cause a frameshift in GRM7, a missense variant in SIK3, and three missense variants in ACAN, each present in one individual except for one of the variants in ACAN, which is present in two people (Additional file 14: Table S10).
We chose 29 variants for confirmation by Sanger sequencing, covering fourteen genes, for our final quality control step. We focussed on variants in the genes common to phenotype-specific subgroups (NEDD4, MON1B, ZAN, Table 1), genes with rare variants in multiple people (SIRPA, ZAN, MON1B, NEFH, ADC, GPR98, LRBA, PAX2 and DNAH2, Table 2 and Additional file 10: Table S6), genes with identical variants in multiple people (PAX2, RBPJ, MON1B, NEFH, Additional file 12: Tables S8 and Additional file 13: Table S9), variants in GWAS genes (GRM7, SIK3, Additional file 14: Table S10) and WFS1, which had 5 variants in one individual (Additional file 11: Table S7). Where genes had multiple predicted variants, we focussed on those present in more than one individual.
Approximately two thirds of our selected variants were confirmed (Table 3 and Additional file 15: Table S11), although the predicted zygosity was not always correct. Most failures were single base pair indels called by Dindel, but three were SNVs called by SAMtools (Table 3). Four additional indels called by Dindel were present but found to be different to the Dindel prediction. For example, eight patients were predicted to bear a single base pair insertion in NEFH, which would result in a frameshift, but in fact Sanger sequencing detected an 18 bp insertion, which is not a frameshift. This miscalling also leads to incorrect minor allele frequencies being associated with each variant (compare the predicted MAF with the MAF of the confirmed variants in Table 3). None of the confirmed indels were detected by SAMtools, even though it is capable of calling small indels. Our findings suggest that while it is valuable to include Dindel, its output should be used with care in variant calling pipelines. In summary, we found a surprisingly high level of false calls from the exome sequencing, confirming that Sanger sequencing should always be used to verify important variants.
Eight genes had confirmed variants, including four candidate novel deafness genes; NEDD4, ZAN, DNAH2, and NEFH. The four known deafness genes with confirmed variants are GPR98, WFS1, GRM7 and SIK3.
Details of the 9 predicted variants in GPR98 are described in Additional file 16: Table S12, including those found in multiple individuals (Table 2), as well as 2 additional variants found in individuals who carry multiple variants in this gene (Additional file 11: Table S7). GPR98 is a large gene (90 exons encoding 6307 amino acids) encoding a G protein coupled receptor, and frameshift mutations in this gene cause one form of Usher’s syndrome, USH2C, an autosomal recessive disorder  causing congenital hearing impairment and retinitis pigmentosa (OMIM #605472). The variants described here in patients without retinitis pigmentosa are heterozygous missense variants, plus one in a splice site (Additional file 14: Table S10), and these are spread throughout the protein from exons 9 to 89. The two missense SNVs we sequenced were confirmed in the three patients predicted to carry them.
All five predicted variants in WFS1 in individual 2590 were confirmed by resequencing. Mutations in WFS1 can cause either a dominantly-inherited non-syndromic hearing loss typically affecting the low frequencies (below 2 kHz) or Wolfram Syndrome (OMIM #222300), a recessive neurodegenerative disease which can include mental retardation, childhood diabetes, optic atrophy and deafness which typically is progressive and affects the high frequencies . The family history, medical history and audiogram shape of individual 2590 are consistent with a dominant non-syndromic hearing loss rather than Wolfram Syndrome (Fig. 3). All five mutations are predicted to be pathogenic but only one of the mutations (p.His313Tyr) has been reported previously, in the heterozygous state in 4 families with Wolfram Syndrome [37, 38]. Individuals previously reported with this variant are all deaf, some also had mental retardation but not all, this mutation is reported as a probably pathogenic mutation in the Wolfram Syndrome Mutation Database (https://lovd.euro-wabb.org/home.php?select_db=WFS1).
Both variants in two of the genes linked to hearing loss by previous GWAS reports were confirmed, in two patients from the metabolic subgroup, one with a deletion in GRM7 and one with a SIK3 SNV (Table 3). GRM7 codes for a metabotropic glutamate receptor, and is expressed in the spiral ganglion and the cochlear and vestibular hair cells in the mouse . The variant was predicted to be a deletion of a T, but we found instead a deletion of TG 34 bp 5′ of the prediction and 22 bp 5′ of the acceptor splice site for exon 6 (transcript ENST00000357716), which makes it an intronic variant and unlikely to affect protein function. SIK3 is a salt-inducible kinase, and is also expressed in the mouse spiral ganglion and cochlear hair cells, although only in young mice . The contribution of these heterozygous variants to the hearing loss of the patients bearing them is unclear.
NEDD4 is a particularly interesting novel candidate; it is known to be widely expressed in the cochlear duct , and the encoded protein is a ubiquitin ligase that binds to and ubiquitylates the products of several deafness genes, including WBP2  and KCNQ1 , both of which are implicated in sensorineural deafness in humans and mice [26, 42,43,44] (Fig. 4). It is involved in AMPA receptor ubiquitination, playing a critical role in AMPA receptor trafficking in rat neurons . The two variants are both missense, one resulting in a substitution of arginine for serine in exon 1, and the other is a substitution of valine for methionine towards the start of exon 1. Patients bearing these variants are all heterozygotes, and fall into the dominant and sensory subgroups.
NEFH, the neurofilament heavy chain, is expressed in the rat cochlear nucleus  and the spiral ganglion neurons . It is known to bind to OTOF , which is involved in deafness in humans and mice [48,49,50], and its expression is affected by several other deafness genes (Fig. 4). We confirmed an 18 bp duplication (rs147489453) in eight people (homozygotes and heterozygotes, from the recessive, sensory and dominant subgroups), and a missense variant (rs59371099) in four people (heterozygotes and one homozygote, from the dominant and recessive subgroups) (Table 3). Four people carried both variants, three from the dominant subgroup and the one individual from the recessive subgroup who was heterozygous for both. The other two patients from the recessive subgroup with a variant in NEFH were homozygous for the 18 bp duplication (Table 3). Both variants are in the last exon of NEFH (exon 4), which has only one isoform. The 18 bp duplication is predicted to duplicate six amino acids in a low complexity region. The missense variant results in a substitution of lysine, which is positively charged, for glutamate, which is negatively charged, which may affect protein function, but without functional studies it is hard to predict what difference either variant would make to the function of the protein.
ZAN is a multiple-domain transmembrane protein found on the apical region of the sperm head , which functions to bind the sperm to the zona pellucida . The MAM domains, D domains, EGF-like domains and mucin-like domains it contains all play a role in adhesion in other proteins, typically in cell-cell or cell-extracellular matrix binding . Although ZAN is thought to be testis-specific, it may also be involved in cell-cell adhesion in the organ of Corti. Variants in ZAN were found in patients from each subgroup, all in the heterozygote state except for three dominant patients who were homozygous for a 35 bp deletion (g.7:100385563_100385597del).
DNAH2 is an axonemal dynein heavy chain, about which very little is known. Six missense variants, one intronic splice site SNV and one single base pair insertion were predicted in seven patients (MAF < 0.5%), and six further missense variants with a MAF < 10% were predicted in 5 patients, of which we have confirmed one missense variant present in the heterozygous state in two patients from the dominant and recessive subgroups.
Expression and mouse mutations of novel candidate deafness genes
We examined the expression of our novel candidate deafness genes using the gEAR portal , which displays data from the mouse organ of Corti at postnatal day (P)0  – P7 , and found that both Nedd4 and Zan were detected during this period. Nedd4 is expressed at reasonably high levels in sensory cells, and only slightly lower levels in non-sensory cells, in accordance with previous publications . Zan has very little expression in the hair cells but is strongly expressed in the supporting cells at P0 . Nefh is expressed in the neurons under the inner hair cells , and is commonly used as a neuronal marker . Dnah2 is expressed in the cochlear duct at P0, but not as strongly as Zan or Nedd4.
There are several mouse lines bearing mutations in Nedd4, Zan, and Nefh, but the only mouse line with a variant in Dnah2 has a large inversion on chromosome 11 covering 1898 genes. No hearing phenotype has been reported for any of the mutant alleles of these genes (as recorded in the Mouse Genome Informatics database ), but these mice may not have been tested for auditory function.
Our analysis has identified several candidate variants and genes for involvement in adult-onset progressive hearing loss, in particular variants in NEFH and NEDD4. It is perhaps surprising that our relatively small sample size of thirty individuals was able to identify good candidate variants and genes. By targeting our approach to specific phenotypic subtypes within large, well-characterised patient cohorts based on audiogram shape, family history and similar age of onset, we have sought to increase power to detect causal genetic variants. In addition, by cross-referencing our data in the filtering pipeline with data from the mouse we were able to prioritise the strongest candidate variants.
The two strongest candidates from our analysis, the variants in NEDD4 and NEFH, should be followed up in larger cohorts and by functional studies to confirm whether they are causal mutations.
Our study also highlights the potential pitfalls of using targeted sequencing to diagnose the cause of adult-onset hearing loss in a typical clinical scenario, where relatives are not available for segregation analysis especially on a gene by gene basis. Of the 30 individuals WES was only able to identify the likely causal mutation in one individual with five WFS1 variants. For the remaining patients the variants identified are of uncertain pathogenicity without further validation. Interpretation of these variants in single individuals is extremely challenging given that even when we limited our analysis to dominantly-inherited human deafness genes, 2334 of the 2504 individuals in the 1000 Genomes Project data carried at least one predicted pathogenic variant in a known dominant human deafness gene. These findings reveal the need for allele frequency databases from carefully-selected controls with good hearing for their age rather than the existing unselected general population controls in studies of highly-prevalent disorders such as hearing loss. Here, we have explored the use of WES in undiagnosed individuals, based on our results it would be interesting to pursue a similar study in individuals who have received a genetic diagnosis to ascertain the number of other predicted pathogenic variants that are present in deafness genes. Furthermore, our study demonstrates that confirmation of candidate variants by Sanger sequencing is always a necessary step. Our findings further suggest that for diagnostics in a clinical situation, novel candidate variants identified by sequencing should be investigated by family analysis if possible, to look for segregation. Without these additional steps, our data suggest that it is not possible to determine with confidence the causative mutation responsible for a patient’s hearing impairment from exome sequence alone. Even with these secondary steps great caution should be exercised in interpreting predicted disease-causing variants, given our findings of the incidence of such variants in every genome.
Age-related hearing loss
European Genome-Phenome Archive
Genome-wide association studies
Minor allele frequency
Medical University of South Carolina
Gates GA, Couropmitree NN, Myers RH. Genetic associations in age-related hearing thresholds. Arch otolaryngol Head Neck Surg. 1999;125:654–9.
Fellinger J, Holzinger D, Pollard R. Mental health of deaf people. Lancet. 2012;379:1037–44.
Karpa MJ, Gopinath B, Beath K, Rochtchina E, Cumming RG, Wang JJ, et al. Associations between hearing impairment and mortality risk in older persons: the Blue Mountains hearing study. Ann Epidemiol. 2010;20:452–9.
Mick P, Kawachi I, Lin FR. The association between hearing loss and social isolation in older adults. Otolaryngol Head Neck Surg. 2014;150:378–84.
Livingston G, Sommerlad A, Orgeta V, Costafreda SG, Huntley J, Ames D, et al. Dementia prevention, intervention, and care. Lancet. 2017;390:2673–734.
Karlsson KK, Harris JR, Svartengren M. Description and primary results from an audiometric study of male twins. Ear Hear. 1997;18:114–20.
Wolber LE, Steves CJ, Spector TD, Williams FM. Hearing ability with age in northern European women: a new web-based approach to genetic studies. PLoS One. 2012;7:e35500.
Fransen E, Bonneux S, Corneveaux JJ, Schrauwen I, Di Berardino F, White CH, et al. Genome-wide association analysis demonstrates the highly polygenic character of age-related hearing impairment. Eur J Hum Genet. 2015;23:110–5.
Huyghe JR, Van Laer L, Hendrickx JJ, Fransen E, Demeester K, Topsakal V, et al. Genome-wide SNP-based linkage scan identifies a locus on 8q24 for an age-related hearing impairment trait. Am J Hum Genet. 2008;83:401–7.
Nolan LS, Maier H, Hermans-Borgmeyer I, Girotto G, Ecob R, Pirastu N, et al. Estrogen-related receptor gamma and hearing function: evidence of a role in humans and mice. Neurobiol Aging. 2013;34:2077 e2071–9.
Wolber LE, Girotto G, Buniello A, Vuckovic D, Pirastu N, Lorente-Canovas B, et al. Salt-inducible kinase 3, SIK3, is a new gene associated with hearing. Hum Mol Genet. 2014;23:6407–18.
Friedman RA, Van Laer L, Huentelman MJ, Sheth SS, Van Eyken E, Corneveaux JJ, et al. GRM7 variants confer susceptibility to age-related hearing impairment. Hum Mol Genet. 2009;18:785–96.
Vuckovic D, Dawson S, Scheffer DI, Rantanen T, Morgan A, Di Stazio M, et al. Genome-wide association analysis on normal hearing function identifies PCDH20 and SLC28A3 as candidates for hearing function and loss. Hum Mol Genet. 2015;24:5655–64.
Hoffmann TJ, Keats BJ, Yoshikawa N, Schaefer C, Risch N, Lustig LR. A large genome-wide association study of age-related hearing impairment using electronic health records. PLoS Genet. 2016;12:e1006371.
Van Camp G, Smith RJ: Hereditary Hearing Loss Homepage. http://hereditaryhearingloss.org/. Accessed Jan 2016.
Van Camp G, Willems PJ, Smith RJ. Nonsyndromic hearing impairment: unparalleled heterogeneity. Am J Hum Genet. 1997;60:758–64.
Nolan LS, Cadge BA, Gomez-Dorado M, Dawson SJ. A functional and genetic analysis of SOD2 promoter variants and their contribution to age-related hearing loss. Mech Ageing Dev. 2013;134:298–306.
British Society of Audiology: BSA Recommended Procedure: Pure tone air-conduction and bone-conduction threshold audiometry with and without masking. http://www.thebsa.org.uk/resources/. Accessed July 2016.
Dubno JR, Eckert MA, Lee FS, Matthews LJ, Schmiedt RA. Classifying human audiometric phenotypes of age-related hearing loss from animal models. J Assoc Res Otolaryngol. 2013;14:687–701.
Lek M, Karczewski K, Minikel E, Samocha K, Banks E, Fennell T, et al. Analysis of protein-coding genetic variation in 60,706 humans. bioRxiv. 2015;
Auton A, Brooks LD, Durbin RM, Garrison EP, Kang HM, Korbel JO, et al. A global reference for human genetic variation. Nature. 2015;526:68–74.
Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M, et al. Primer3--new capabilities and interfaces. Nucleic Acids Res. 2012;40:e115.
Bonfield JK, Smith K, Staden R. A new DNA sequence assembly program. Nucleic Acids Res. 1995;23:4992–9.
The International Genome Sample Resource: 1000 Genomes data. http://ftp.1000genomes.ebi.ac.uk/vol1/ftp/release/20130502/. Accessed June 2016.
Ohlemiller KK, Jones SM, Johnson KR. Application of mouse models to research in hearing and balance. J Assoc Res Otolaryngol. 2016;17:493–523.
Buniello A, Ingham NJ, Lewis MA, Huma AC, Martinez-Vega R, Varela-Nieto I, et al. Wbp2 is required for normal glutamatergic synapses in the cochlea and is crucial for hearing. EMBO Mol Med. 2016;8:191–207.
Lewis MA, Quint E, Glazier AM, Fuchs H, De Angelis MH, Langford C, et al. An ENU-induced mutation of miR-96 associated with progressive hearing loss in mice. Nat Genet. 2009;41:614–8.
Mencia A, Modamio-Hoybjor S, Redshaw N, Morin M, Mayo-Merino F, Olavarrieta L, et al. Mutations in the seed region of human miR-96 are responsible for nonsyndromic progressive hearing loss. Nat Genet. 2009;41:609–13.
Fuentes Fajardo KV, Adams D, Mason CE, Sincan M, Tifft C, Toro C, et al. Detecting false-positive signals in exome sequencing. Hum Mutat. 2012;33:609–13.
Daiger SP, Rossiter BJF, Greenberg J, Christoffels A, Hide W. Data services and software for identifying genes and mutations causing retinal degeneration. Invest Opthalmol Vis Sci. 1998;39:S295.
Daiger SP, Sullivan LS, Bowne SJ: RetNet Retinal Information Network https://sph.uth.edu/retnet/home.htm. Accessed June 2016.
Tabor HK, Auer PL, Jamal SM, Chong JX, Yu JH, Gordon AS, et al. Pathogenic variants for Mendelian and complex traits in exomes of 6,517 European and African Americans: implications for the return of incidental results. Am J Hum Genet. 2014;95:183–93.
Shearer AE, Eppsteiner RW, Booth KT, Ephraim SS, Gurrola J 2nd, Simpson A, et al. Utilizing ethnic-specific differences in minor allele frequency to recategorize reported pathogenic deafness variants. Am J Hum Genet. 2014;95:445–53.
Hafidi A, Despres G, Romand R. Cochlear innervation in the developing rat: an immunocytochemical study of neurofilament and spectrin proteins. J Comp Neurol. 1990;300:153–61.
Jacobson SG, Cideciyan AV, Aleman TS, Sumaroka A, Roman AJ, Gardner LM, et al. Usher syndromes due to MYO7A, PCDH15, USH2A or GPR98 mutations share retinal disease mechanism. Hum Mol Genet. 2008;17:2405–15.
Cryns K, Sivakumaran TA, Van den Ouweland JM, Pennings RJ, Cremers CW, Flothmann K, et al. Mutational spectrum of the WFS1 gene in Wolfram syndrome, nonsyndromic hearing impairment, diabetes mellitus, and psychiatric disease. Hum Mutat. 2003;22:275–87.
Hansen L, Eiberg H, Barrett T, Bek T, Kjaersgaard P, Tranebjaerg L, et al. Mutation analysis of the WFS1 gene in seven Danish Wolfram syndrome families; four new mutations identified. Eur J Hum Genet. 2005;13:1275–84.
Rohayem J, Ehlers C, Wiedemann B, Holl R, Oexle K, Kordonouri O, et al. Diabetes and neurodegeneration in Wolfram syndrome: a multicenter study of phenotype and genotype. Diabetes Care. 2011;34:1503–10.
Zhong SX, Liu ZH. Expression patterns of Nedd4 isoforms and SGK1 in the rat cochlea. Acta Otolaryngol. 2009;129:935–9.
Persaud A, Alberts P, Amsen EM, Xiong X, Wasmuth J, Saadon Z, et al. Comparison of substrate specificity of the ubiquitin ligases Nedd4 and Nedd4-2 using proteome arrays. Mol Syst Biol. 2009;5:333.
Jespersen T, Membrez M, Nicolas CS, Pitard B, Staub O, Olesen SP, et al. The KCNQ1 potassium channel is down-regulated by ubiquitylating enzymes of the Nedd4/Nedd4-like family. Cardiovasc Res. 2007;74:64–74.
Casimiro MC, Knollmann BC, Ebert SN, Vary JC Jr, Greene AE, Franz MR, et al. Targeted disruption of the Kcnq1 gene produces a mouse model of Jervell and Lange-Nielsen syndrome. Proc Natl Acad Sci U S A. 2001;98:2526–31.
Lee MP, Ravenel JD, Hu RJ, Lustig LR, Tomaselli G, Berger RD, et al. Targeted disruption of the Kvlqt1 gene causes deafness and gastric hyperplasia in mice. J Clin Invest. 2000;106:1447–55.
Neyroud N, Tesson F, Denjoy I, Leibovici M, Donger C, Barhanin J, et al. A novel mutation in the potassium channel gene KVLQT1 causes the Jervell and Lange-Nielsen cardioauditory syndrome. Nat Genet. 1997;15:186–9.
Lin A, Hou Q, Jarzylo L, Amato S, Gilbert J, Shang F, et al. Nedd4-mediated AMPA receptor ubiquitination regulates receptor turnover and trafficking. J Neurochem. 2011;119:27–39.
Friedland DR, Popper P, Eernisse R, Ringger B, Cioffi JA. Differential expression of cytoskeletal genes in the cochlear nucleus. Anat Rec A Discov Mol Cell Evol Biol. 2006;288:447–65.
Duncker SV, Franz C, Kuhn S, Schulte U, Campanelli D, Brandt N, et al. Otoferlin couples to clathrin-mediated endocytosis in mature cochlear inner hair cells. J Neurosci. 2013;33:9508–19.
Longo-Guess C, Gagnon LH, Bergstrom DE, Johnson KR. A missense mutation in the conserved C2B domain of otoferlin causes deafness in a new mouse model of DFNB9. Hear Res. 2007;234:21–8.
Roux I, Safieddine S, Nouvian R, Grati M, Simmler MC, Bahloul A, et al. Otoferlin, defective in a human deafness form, is essential for exocytosis at the auditory ribbon synapse. Cell. 2006;127:277–89.
Yasunaga S, Grati M, Cohen-Salmon M, El-Amraoui A, Mustapha M, Salem N, et al. A mutation in OTOF, encoding otoferlin, a FER-1-like protein, causes DFNB9, a nonsyndromic form of deafness. Nat Genet. 1999;21:363–9.
Gao Z, Garbers DL. Species diversity in the structure of zonadhesin, a sperm-specific membrane protein containing multiple cell adhesion molecule-like domains. J Biol Chem. 1998;273:3415–21.
Lea IA, Sivashanmugam P, O'Rand MG. Zonadhesin: characterization, localization, and zona pellucida binding. Biol Reprod. 2001;65:1691–700.
Hertzano RP, Orvis J. gEAR Portal. https://umgear.org/. Accessed June 2016.
Cai T, Jen HI, Kang H, Klisch TJ, Zoghbi HY, Groves AK. Characterization of the transcriptome of nascent hair cells and identification of direct targets of the Atoh1 transcription factor. J Neurosci. 2015;35:5870–83.
Scheffer DI, Shen J, Corey DP, Chen ZY. Gene expression by mouse inner ear hair cells during development. J Neurosci. 2015;35:6366–80.
The Jackson Laboratory: Mouse Genome Informatics. www.informatics.jax.org. Accessed Sept 2017.
EGA Consortium: European Genome-Phenome Archive. http://ega-archive.org. Accessed June 2016.
Camarero G, Avendano C, Fernandez-Moreno C, Villar A, Contreras J, de Pablo F, et al. Delayed inner ear maturation and neuronal loss in postnatal Igf-1-deficient mice. J Neurosci. 2001;21:7630–41.
Brady ST, Witt AS, Kirkpatrick LL, de Waegh SM, Readhead C, Tu PH, et al. Formation of compact myelin is required for maturation of the axonal cytoskeleton. J Neurosci. 1999;19:7278–88.
Zer C, Sachs G, Shin JM. Identification of genomic targets downstream of p38 mitogen-activated protein kinase pathway mediating tumor necrosis factor-alpha signaling. Physiol Genomics. 2007;31:343–51.
Fujiyoshi T, Hood L, Yoo TJ. Restoration of brain stem auditory-evoked potentials by gene transfer in shiverer mice. Ann Otolo Rhinol Laryngol. 1994;103:449–56.
Oishi N, Chen J, Zheng HW, Hill K, Schacht J, Sha SH. Tumor necrosis factor-alpha-mutant mice exhibit high frequency hearing loss. J Assoc Res Otolaryngol. 2013;14:801–11.
Fukushima T, Yoshihara H, Furuta H, Kamei H, Hakuno F, Luan J, et al. Nedd4-induced monoubiquitination of IRS-2 enhances IGF signalling and mitogenic activity. Nat Commun. 2015;6:6780.
Edwin F, Anderson K, Patel TB. HECT domain-containing E3 ubiquitin ligase Nedd4 interacts with and ubiquitinates Sprouty2. J Biol Chem. 2010;285:255–64.
Wang X, Trotman LC, Koppie T, Alimonti A, Chen Z, Gao Z, et al. NEDD4-1 is a proto-oncogenic ubiquitin ligase for PTEN. Cell. 2007;128:129–39.
Trotman LC, Wang X, Alimonti A, Chen Z, Teruya-Feldstein J, Yang H, et al. Ubiquitination regulates PTEN nuclear import and tumor suppression. Cell. 2007;128:141–56.
Murillo-Cuesta S, Camarero G, Gonzalez-Rodriguez A, De La Rosa LR, Burks DJ, Avendano C, et al. Insulin receptor substrate 2 (IRS2)-deficient mice show sensorineural hearing loss that is delayed by concomitant protein tyrosine phosphatase 1B (PTP1B) loss of function. Mol Med. 2012;18:260–9.
Dong Y, Sui L, Yamaguchi F, Kamitori K, Hirata Y, Hossain MA, et al. Phosphatase and tensin homolog deleted on chromosome 10 regulates sensory cell proliferation and differentiation of hair bundles in the mammalian cochlea. Neuroscience. 2010;170:1304–13.
Shim K, Minowada G, Coling DE, Martin GR. Sprouty2, a mouse deafness gene, regulates cell fate decisions in the auditory sensory epithelium by antagonizing FGF signaling. Dev Cell. 2005;8:553–64.
(ESP) NGESP: Exome Variant Server. http://evs.gs.washington.edu/EVS/. Accessed July 2018.
We thank Carol Scott, Wellcome Trust Sanger Institute, for help and support with the exome sequencing study and variant calling pipeline.
This work was supported by the following: the Medical Research Council, the Wellcome Trust [grant numbers 098051; 100669 to KPS]; National Institutes of Health/National Institute on Deafness and Other Communication Disorders [grant number P50 DC000422]; the South Carolina Clinical & Translational Research Institute, the Medical University of South Carolina [grant numbers NIH/NCATS UL1 TR000062 and UL1 TR001450]; the Haigh Fellowship in Age Related Deafness and Deafness Research UK, now Action on Hearing Loss [grant number 444:UEI:SD to SJD]. The London ARHL Cohort collection was initiated by funding from Research into Ageing [grant number 223 to SJD]; Teresa Rosenbaum Golden Charitable Trust, Telethon Foundation [grant number GGP09037]. The funders played no part in the design of the study or the collection, analysis and interpretation of the data, or the writing of the manuscript.
Availability of data and materials
The datasets generated and/or analyzed during the current study are available in European Genome-Phenome Archive . The accession number for the twenty UK adults is EGA:EGAS00001000295, and the accession number for the ten older individuals from the US is EGA:EGAS00001000157. Access to the data may be requested through the EGA. The 1000 Genomes data is publically available from the 1000 Genomes repository [http://ftp.1000genomes.ebi.ac.uk/vol1/ftp/release/20130502/] .
Ethics approval and consent to participate
Ethical approval for the patients recruited in the U.K. was granted from the Royal Free Local Research Ethics Committee (reference 6202). For patients recruited in the U.S. the protocols for this study were approved by the Institutional Review Board at MUSC. All participants were recruited by written informed consent.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. showing the audiograms of each participant in the recessive patient group. (PDF 486 kb)
Figure S2. showing the audiograms of each participant in the dominant patient group. (PDF 485 kb)
Figure S3. showing the audiograms of each participant in the metabolic and sensory patient groups. (PDF 428 kb)
Table S1. which lists the software prediction tools used to call and annotate variants, with relevant references. (XLS 30 kb)
Figure S4. and legend detailing counts at each filtering step for common variant analysis. (PDF 259 kb)
Table S2. detailing primers used for confirming selected variants. (DOCX 20 kb)
Table S3. listing genes reported to underlie deafness in humans and/or mice. (XLS 45 kb)
Table S4. and lists candidates for exclusion - genes with predicted pathogenic variants in many exome projects or across multiple families exhibiting a variety of phenotypes. This list consists of the genes in Supplementary Tables S1 and S2 from . (XLS 46 kb)
Table S6. listing individuals homozygous for very rare variants in any gene. (DOCX 15 kb)
Table S7. listing individuals with more than one very rare mutation in the same gene. (DOCX 20 kb)
Table S8. giving details of the single very rare variants in known deafness genes found in multiple samples. (DOCX 16 kb)
Table S9. giving details of genes with identical variants found in more than one person. (DOCX 16 kb)
Table S10. detailing very rare mutations in GWAS candidates. (DOCX 20 kb)
Table S11. which contains further details of the validated variants listed in Table 3. (XLS 32 kb)
Table S12. giving details of very rare mutations identified in GPR98. (DOCX 19 kb)
About this article
Cite this article
Lewis, M.A., Nolan, L.S., Cadge, B.A. et al. Whole exome sequencing in adult-onset hearing loss reveals a high load of predicted pathogenic variants in known deafness-associated genes and identifies new candidate genes. BMC Med Genomics 11, 77 (2018). https://doi.org/10.1186/s12920-018-0395-1