Finding New Genes for Non-Syndromic Hearing Loss through an In Silico Prioritization Study

Abstract

At present, 51 genes are already known to be responsible for Non-Syndromic hereditary Hearing Loss (NSHL), but the knowledge of 121 NSHL-linked chromosomal regions brings to the hypothesis that a number of disease genes have still to be uncovered. To help scientists to find new NSHL genes, we built a gene-scoring system, integrating Gene Ontology, NCBI Gene and Map Viewer databases, which prioritizes the candidate genes according to their probability to cause NSHL. We defined a set of candidates and measured their functional similarity with respect to the disease gene set, computing a score () that relies on the assumption that functionally related genes might contribute to the same (disease) phenotype. A Kolmogorov-Smirnov test, comparing the pair-wise distribution on the disease gene set with the distribution on the remaining human genes, provided a statistical assessment of this assumption. We found at a p-value that the former pair-wise is greater than the latter, justifying a prioritization strategy based on the functional similarity of candidate genes respect to the disease gene set. A cross-validation test measured to what extent the ranking for NSHL is different from a random ordering: adding 15% of the disease genes to the candidate gene set, the ranking of the disease genes in the first eight positions resulted statistically different from a hypergeometric distribution with a p-value and a power. The twenty top-scored genes were finally examined to evaluate their possible involvement in NSHL. We found that half of them are known to be expressed in human inner ear or cochlea and are mainly involved in remodeling and organization of actin formation and maintenance of the cilia and the endocochlear potential. These findings strongly indicate that our metric was able to suggest excellent NSHL candidates to be screened in patients and controls for causative mutations.

About the Authors:

Matteo Accetturo

Contributed equally to this work with: Matteo Accetturo, Teresa M. Creanza, Claudia Santoro

Affiliation: GBS BAO Advanced Analytics Services and MBLab, IBM, Bari, Italy

Teresa M. Creanza

Contributed equally to this work with: Matteo Accetturo, Teresa M. Creanza, Claudia Santoro

Affiliation: GBS BAO Advanced Analytics Services and MBLab, IBM, Bari, Italy

Claudia Santoro

Contributed equally to this work with: Matteo Accetturo, Teresa M. Creanza, Claudia Santoro

Affiliation: GBS BAO Advanced Analytics Services and MBLab, IBM, Bari, Italy

Giancarlo Tria

Affiliation: GBS BAO Advanced Analytics Services and MBLab, IBM, Bari, Italy

Antonio Giordano

Affiliation: GBS BAO Advanced Analytics Services and MBLab, IBM, Bari, Italy

Simone Battagliero

Affiliation: GBS BAO Advanced Analytics Services and MBLab, IBM, Bari, Italy

Antonella Vaccina

Affiliation: GBS BAO Advanced Analytics Services and MBLab, IBM, Bari, Italy

Gaetano Scioscia

Affiliation: GBS BAO Advanced Analytics Services and MBLab, IBM, Bari, Italy

Pietro Leo

* E-mail: [email protected]

Affiliation: GBS BAO Advanced Analytics Services and MBLab, IBM, Bari, Italy

Introduction

Non Syndromic hereditary Hearing Loss (NSHL) is one of the most genetically heterogeneous disorders known. Indeed it can present an autosomal recessive, autosomal dominant, X-linked or mitochondrial pattern of inheritance; furthermore, mutations in the same gene may cause syndromic or non syndromic hearing loss, and recessive forms may be caused by a combination of two mutations in different genes from the same functional group [1].

Due to this tremendous genetic heterogeneity, the identification of genes and gene defects that affect the process of hearing is challenging [1]. At present 51 genes have been already identified to be responsible, if mutated, for this phenotype (see Table 1 for references); nevertheless not all these genes have been fully characterized. They usually are involved in the inner ear development or functionality, and their mutations generally cause hearing loss interfering in the process of the elaboration of sound.

[Figure omitted. See PDF.]

Table 1. NSHL disease genes.

https://doi.org/10.1371/journal.pone.0012742.t001

About 50% of cases of NSHL are due to mutations of GJB2, a gene coding for a gap-junction protein called connexin 26, involved in the cell-cell communication process. Another important gene responsible for NSHL is GJB6, belonging to the same family of GJB2 and adjacent to it. The identification of these two genes highlighted the role of connexins, and therefore of the cochlear gap-junction ion channels, in the auditory function [2], [3].

However the biology of hearing is extremely complex and many other different classes of genes are involved in NSHL. For instance, SLC26A4, associated with autosomal recessive NSHL [4] and Pendred syndrome, is a gene coding for pendrin, a chloride/iodide transporter; COCH, responsible for autosomal dominant non syndromic post-lingual with a progressive onset in adulthood [5], encodes for cochlin, a component of the extracellular matrix of the inner ear; POU3F4, responsible for an X-linked non syndromic progressive and profound sensorineural hearing loss [6], encodes for a transcription factor; while WFS1 associated with autosomal recessive Wolfram syndrome and autosomal dominant low frequency NSHL [7], [8], is a gene coding for the glycoprotein wolframin.

Moreover, several linkage studies over the years have shown that many chromosomal regions are involved in NSHL. At present 121 loci are known to be involved in this phenotype [9], and for many of them the genes causing NSHL have not been identified yet. Due to their often extremely large dimensions – they can even contain several hundreds of genes – it is not feasible to experimentally validate all the genes contained in each locus. In addition, some loci might contain more than one disease gene, as in the case of DFNA3 that harbors GJB2 and GJB6.

In this scenario, a bioinformatic approach to narrow down the list of possible candidate genes is an essential requirement in order to experimentally validate first those genes most likely associated with the disease.

Many strategies have been devised to address this issue, mostly sharing the common prioritization idea of ranking the candidate genes on the basis of their similarity with a set of training genes – genes already associated to the phenotype – relying on the main assumption that genes whose dysfunction contributes to a disease phenotype tend to be functionally related (see [10] and references within).

Quantifying the functional relatedness between two genes is not trivial; often existing information about gene function are exploited to infer functional relationships among genes. In this kind of approach an excellent means is provided by Gene Ontology (GO, The Gene Ontology Consortium, 2001) [11], which is the golden standard ontology in the field of genes and gene products.

Indeed one of the advantages of having genes annotated with GO terms is the possibility to compare them not only from a qualitative point of view (e.g. by searching for common terms with which they are annotated), but also by defining an explicit semantic similarity measure which reflects the closeness in meaning between the terms with which they are annotated [12], [13]. This semantic similarity measure gives in turn a measure of the functional similarity of the annotated gene products, as extensively discussed in Pesquita et al [12].

Briefly, when comparing two terms in an ontology, two main approaches are generally distinguished, the edge-based, which counts the edges in the graph path between two terms [14]–[18], and the node-based, which looks at the properties of the terms, their ancestors and descendants [19]–[25]. Most of the node-based similarity measures are functions of the information content (IC) of each term, and their most informative common ancestors [25]. IC is the amount of information a term contains, meaning that a term contains less information if it occurs very often; in this context the similarity between two terms is quantified looking at the amount of information they share. Very often gene products are annotated with multiple GO terms, in this case maximum [26], [27], average [13], [18], [28] or sum [29] of the GO term similarities may be taken as the gene similarity.

Here we define a new Semantic Similarity Measure (SSM) between gene products by directly extending to sets of concepts (the gene annotations) the Lin's idea [25] of quantifying the similarity between two concepts in an ontology. Our metric provides a measure of the functional similarity between two genes and its reliability is tested in this paper in the context of gene prioritization for NSHL. Indeed the overall aim of this paper is (i) to support researchers in search of new genes responsible for NSHL and (ii) provide indications about the main biological processes, molecular functions and cellular components to be explored to study NSHL, by defining a procedure to computationally prioritize candidate genes for their association with this phenotype. The availability of a good training gene set for NSHL – 51 genes already associated with this phenotype (disease genes) – allows to select new genes most likely responsible for this phenotype estimating their similarity with the disease gene set.

Finally we define a systematic and unbiased statistical assessment to validate the obtained results.

Results

The candidate genes prioritized for NSHL in this study were selected as described in the Methods section. They were prioritized against all the genes already known to be responsible for NSHL (disease genes, see Methods section for details on their selection), according to a score which is function of the Semantic Similarity Measure () estimated for each candidate-disease gene pair. All candidate genes were ranked by computing the for each candidate-disease gene pair; the final score used for prioritizing each candidate was obtained as the mean of the scores estimated for that candidate against all the disease genes and was defined Semantic Similarity Measure Average ().

Validation of the for NSHL gene prioritization

Before being able to assert that the ranking produced by is worthy of attention and therefore evaluating it from a biological point of view, we wanted to evaluate two main aspects concerning our prioritization methodology. We first wanted to test whether the main hypothesis upon which this and most of the prioritization studies are based – genes whose dysfunction contributes to a disease phenotype tend to be functionally related – is quantifiable in terms of semantic similarity, especially in the particular case of NSHL, where the complexity of the hearing process and the complexity of the genetics of the disease both play an important role. Second aspect is whether our metric is able to catch this functional relatedness. To test these two aspects is equivalent to answer the following question: are the disease genes more functionally related than two generic human genes according to ? A positive answer would yield a positive result for both aspects at the same time, implying that the more a candidate gene obtains a high score respect to the disease gene set, the higher is its probability to cause NSHL when mutated. To address this issue, we estimated the pair-wise SSM distribution on the disease gene set, and compared it with the pair-wise SSM distribution estimated on the entire human gene set. In Figure 1 a population pyramid shows the pair-wise SSM distribution across the disease genes and All-human-genes sets in two back-to-back histograms. It provides the graphical evidence that the majority of the disease gene pairs assume values in the range of 0.5–0.6, much greater than those assumed by the majority of all the remaining human genes (around 0.4). This clearly indicates that the NSHL genes are more functionally related in terms of similarity than two generic human genes. In order to statistically support this result, we formulated the following test: the null hypothesis is that the pair-wise SSM distribution in the disease genes set is equal to the pair-wise SSM distribution in the All-human-genes set, while the alternative hypothesis is that the former is greater than the latter, i.e. the cdf (cumulative distribution function) of the former population is smaller than the cdf of the latter population. The test was performed using the bootstrap version of the Kolmogorov-Smirnov test (ks.boot), which allows ties and is included in the R package Matching [30]. We found a p-value, confirming the hypothesis that the disease genes are indeed more similar according to than two generic human genes. This evidence shows the ability of our metric in capturing the functional relatedness of NSHL genes respect to the rest of all human genes, justifying therefore a gene prioritization strategy for association with NSHL based on the similarity of the candidate genes with respect to the disease gene set.

[Figure omitted. See PDF.]

Figure 1. Similarity population pyramid.

Back-to-back histograms showing the asymmetry in frequencies of SSM values (in 0.1 bin interval between 0 and 1) among gene pairs, for disease genes (on the right) and the entire human gene set (on the left).

https://doi.org/10.1371/journal.pone.0012742.g001

In order to validate the reliability of in ranking the candidate genes with respect to their probability to play a causative role in NSHL manifestation, we designed a specific cross-validation procedure that quantifies how much the ranking obtained with our metric differs from a random ordering of the candidate genes. Indeed, due to the specific context we are dealing with, i.e. the gene prioritization, we could not use the classical cross-validation procedure, we in fact added 15% of the disease genes randomly drawn for 10000 times from the disease gene set to the candidates, and counted each time the number of the diseases that fell in the top four windows of 100, 75, 50 and 8 genes. Here the candidate gene set was used exclusively to produce noise, as the positions of the candidates in the ranking were never evaluated during the cross-validation procedure. We in fact tested if the number of disease genes ranked in the top windows were significantly greater than expected when a random extraction of 100, 75, 50 and 8 genes was performed from the total (candidates plus 15% of disease genes) gene set. In Figure 2 we report the distributions obtained from the cross-validation procedure (in blue) applied to the four top windows. In this figure we compare these distributions with the hypergeometric ones (in red), which mimic the random extraction of 100, 75, 50 and 8 genes from the 8748 genes (8740 candidate genes plus 8 disease genes). In all four cases the two distributions are clearly distinct (i.e. the overlapping regions are small). Moreover the means of the distributions for the cross-validation (blue triangles in the figure) result always greater than the means of the hypergeometric distributions (red triangles in the figure). This confirms that the ranking computed by our gene scoring system is significantly different from a random ordering. This is equivalent to assert that our scoring system is able to put at the top of the ranking those genes which are functionally more related to the NSHL genes and thus, more likely, potentially to cause the disease when mutated. This evidence is statistically supported as the p-value and the power of the test for each of the four windows (see Methods section) resulted always smaller than 0.01 and greater than 0.99, respectively (Table 2).

[Figure omitted. See PDF.]

Figure 2. Cross-validation and hypergeometric distributions in case of (a)100, (b)75, (c)50, (d)8 window widths.

In red the hypergeometric distributions with their expectation values (red filled triangles); in blue the distributions, estimated by cross-validation, of disease genes in the top-ranked genes with their mean values (blue filled triangles).

https://doi.org/10.1371/journal.pone.0012742.g002

[Figure omitted. See PDF.]

Table 2. Gene scoring system cross-validation.

https://doi.org/10.1371/journal.pone.0012742.t002

Analysis of the top-ranked candidate genes

The candidates ranked according to were then examined looking at their functions and expression sites. The twenty top-scored genes are reported in Table 3 together with a brief description of their functions. The number of 20 was arbitrarily chosen, mainly thinking about the intrinsic technical limitations of experimentally testing a great number of genes for disease association – this is actually the reason why such prioritization studies are becoming routine.

[Figure omitted. See PDF.]

Table 3. Top-ranked candidate genes.

https://doi.org/10.1371/journal.pone.0012742.t003

Half of them are reported in literature to be expressed in human inner ear or cochlea, despite the very limited availability of gene expression data for these tissues due to the technical difficulties of obtaining undamaged hair-cell samples for gene expression experiments. For the remaining genes, six are reported to be expressed in other organisms' inner ear or cochlea, mainly mouse or chicken, while four have no gene expression data for these tissues. Taken altogether, these are important indications supporting the goodness of the ranking we produced in respect to the NSHL, especially if we think that the initial candidate gene list was not a priori filtered by any criterium except that of being all annotated genes located in the susceptibility loci.

Moreover, looking at their functions, we found that most of the top-ranked genes play roles compatible with a possible involvement in NSHL phenotype. Among the most relevant, we identified a) processes of remodeling and organization of actin (WDR1, CLRN1, FLII), an essential component of the hair-cell bundle; b) formation and maintenance of cilia (ALMS1, USH1G, CC2D2A), the sensory organelles devoted to receive the mechanical stimulus; c) cycling and pH homeostasis in cochlear fluids (ATP6V0A4, KCNQ1, KCNE1L), essentials for the generation and maintenance of the endocochlear potential; d) signal transduction (PTPLA, PTPN11, TBL1X, TIMM8A). They are all important molecular mechanisms underlying the hearing process, which involve the hair cell capability to transduce the mechanical stimulus into electrical signal, as well as the endolymph production and maintenance.

Stronger evidences come from some of the top-ranked genes which are already linked to different syndromic forms of deafness: USH1G for instance is known to cause Usher syndrome type 1G [31], associated with sensorineural hearing impairment; for this gene a possible role in the development and maintenance of the stereocilia bundles is reported by Weil et al. [31]: it might in fact function as an anchoring/scaffolding protein in hair cells and could be involved in the functional network formed by USH1C, CDH23 and MYO7A that is required for cohesion of the growing hair bundle, making its role in the hearing impairment process quite easily explainable. Similarly, KCNE1L has been associated by Piccini et al. [32] to AMME syndrome (Alport syndrome - mental retardation - midface hypoplasia - elliptocytosis) whose symptoms include, among others, hearing loss, and analogous situations are reported also for TIMM8A, involved in Mohr-Tranebjaerg syndrome [33] and Jensen syndrome [34], and ALMS1, involved in Alström syndrome [35]. It is noteworthy that the association of some top-ranked genes to syndromic deafness forms does not exclude them from being good NSHL candidates, as clearly demonstrated by USH1C involved both in Usher syndrome type 1C [36], and NSHL [37], [38], depending on which mutations it undergoes.

Finally, we produced a graphical bidimensional representation of the 20 top-ranked genes together with the disease genes using Proxscal SPSS, which performs multidimensional scaling of similarity data to find a least squares representation of the objects in a low-dimensional space (Figure 3). The proximity of the two gene sets was in this way highlighted; this allowed identifying different groups of NSHL disease genes (red balls in the figure) – namely myosins, connexins, cadherins, ion channels and so forth – and mapping the best candidates within these groups. The inclusion of the top-scored candidate genes did not enlarge the area occupied by the disease genes and their membership to the relative subgroups was mantained in the graphical representation.

[Figure omitted. See PDF.]

Figure 3. Multidimensional scaling of similarity data to represent the disease and the 20 top-scored candidate genes in a bidimensional space.

https://doi.org/10.1371/journal.pone.0012742.g003

Overall, on the basis of these considerations, the majority of them seem to be excellent candidates for subsequent studies on NSHL patients and controls.

Functional characterization of candidate and disease genes using GO

In order to further investigate the obtained ranking and in order to have a more general picture of the molecular functions, biological processes and cellular components more associated to NSHL, as suggested by both the best candidates and disease genes, we designed and implemented two specific statistical tests that allowed to identify the GO terms more representative of NSHL, exploiting the score estimated by our gene scoring system. For the disease genes, we quantified and tested the enrichment of gene-sets defined by functional categories provided by Gene Ontology annotations in disease gene list. In this case the score was used to define the non-disease gene class (see Methods section). For the candidates, we analyzed all GO terms in their annotations, and evaluated the enrichment of the gene set annotated with each of them, by using the score obtained from our ranking to quantify their association with NSHL (see Materials and Methods section). In this case the score allowed us, starting from the GO terms associated to all the candidate genes, to identify those GO terms significantly associated with the best candidates, without making any a priori decision on which candidates should be considered as the “best” candidates.

This survey had the purpose to examine the ranking on a larger scale – extending the ranking examination to the whole candidate gene set – to possibly suggest non-obvious pathways to further look into when studying NSHL, hence it was devised as a way to look at the results from a different point of view (i.e. moving from a view of NSHL in terms of genes responsible of the disease to a view of NSHL in terms of biological processes, molecular functions and cellular components distinctive of the disease).

We considered as significantly descriptive of the best candidate and disease genes, those GO terms with a p-value and we ordered them according to their score, function of their p-value and specificity in the corpus of the GO annotations.

As for the candidate genes, the enriched terms, divided into biological processes, molecular functions and cellular components (Table 4), include expected concepts such as “auditory receptor cell stereocilium organization” (GO:0060088), “large conductance calcium-activated potassium channel activity” (GO:0060072), “sensory perception of sound” (GO:0007605), “auditory receptor cell stereocilium organization” (GO:0060088), consistent with hearing physiology, as well as less obvious functions or processes such as “regulation of circadian sleep/wake cycle, REM and non-REM sleep” (GO:0042320, GO:0045188), “response to cocaine” (GO:0042220), or “mu-type opioid receptor binding” (GO:0031852), that need further (experimental) investigations. This on the one hand supports again the goodness of the ranking, confirming that the top-scored genes are actually promising candidates for association with NSHL, on the other hand fulfils our initial requirement to suggest new prospective insights in NSHL.

[Figure omitted. See PDF.]

Table 4. Enriched biological processes, cellular components and molecular functions for candidate genes.

https://doi.org/10.1371/journal.pone.0012742.t004

As for the disease genes, as expected, the enriched terms are all consistent with the hearing physiology (Table 5). To give some examples, among the most relevant enriched biological processes we found “actin filament-based movement” (GO:0030048), “inner ear morphogenesis” (GO:0042472), “regulation of cell shape” (GO:0008360) and a group involving sensory perception (GO:0007605, GO:0007601, GO:0050957). Likewise, among the enriched cellular components, are “stereocilium” (GO:0032420), “myosin complex” (GO:0016459), “cell junction” (GO:0030054), and among the molecular functions, “actin binding” (GO:0003779), “actin filament binding” (GO:0051015), and so forth.

[Figure omitted. See PDF.]

Table 5. Enriched biological processes, cellular components and molecular functions for disease genes.

https://doi.org/10.1371/journal.pone.0012742.t005

Interestingly, among all the enriched terms – for both candidate and disease genes – there is a very small amount of overlapping. Only two biological processes are in fact shared between the two gene lists, “sensory perception of sound” (GO:0007605) and “equilibrioception” (GO:0050957), which are neverthless extremely specific terms – a very small number of gene products are annotated with these terms – both deeply linked to the inner ear function. Looking at the GO graph, however, many of the non-shared terms are interconnected with each other on a larger scale, sharing a common parent at different levels of specificity. This is due to the structure of our algorithm that favours the closeness in the graph of the terms in estimating the similarity between genes. It is noteworthy that with this approach we can think of NSHL from a different perspective, exploring portions of the graph that otherwise would have never been explored.

In Figure 4 we reported an elucidative example of this issue: by mapping some enriched disease and candidate biological processes to the GO graph, we observed that the addition of “inner ear receptor cell differentiation” (GO:0060113) to the list of NSHL possible biological processes clearly enlarges the NSHL subgraph covering a new branch of the “inner ear development” (GO:0048839) different from the “inner ear morphogenesis” (GO:0042472), while the addition of “auditory receptor cell stereocilium organization” (GO:0060088) narrows and specializes the concept “inner ear morphogenesis” to one of its components.

[Figure omitted. See PDF.]

Figure 4. GO subgraph of some disease and candidate gene enriched GO terms.

Red circles indicate terms enriched for the disease genes, green circles indicate terms enriched for the candidate genes. Dark blue arrows indicate is a relations, light blue arrows indicate part of relations between the terms.

https://doi.org/10.1371/journal.pone.0012742.g004

These findings, as a whole, on the one hand support again the goodness of the ranking, on the other hand they suggest that also some pathways apparently unrelated with NSHL, might deserve future attention by NSHL researchers.

Discussion

In the perspective of discovering new genes potentially involved in NSHL, we built a gene scoring system integrating Gene Ontology (GO), NCBI Gene and Map Viewer databases, which scores the candidate genes for NSHL by comparing them with the 51 NSHL disease genes already known, relying on the assumption that functionally related genes might contribute to the same (disease) phenotype.

We defined a set of candidate genes for NSHL as all the genes contained in the susceptibility loci known so far, and we prioritized them for the association with the disease, without making any a priori selection except that of being annotated with at least one GO term.

We first of all tested whether our metric, , was able to capture the above assumption, verifying that the disease genes are indeed more similar, according to the metric, than two generic human genes. We also demonstrated that our metric is able to pool the disease genes respect to the other human genes, implying that the former are indeed more closely functionally related than the latter: these results therefore justify a prioritization strategy based on the similarity of the candidate genes respect to the disease gene set.

Afterwards, we wanted to investigate to what extent our metric is reliable in ranking candidate genes for their potential role in NSHL manifestation. To this purpose we designed a cross-validation procedure and we obtained excellent results also considering the more disadvantageous condition of ranking eight disease genes in the first 8 positions of a list of more than 8700 genes.

Given these preliminary validations, we are extremely confident that the ranking we produced with respect to NSHL is worthy of attention for future NSHL research plan. Indeed, the top-scored candidate genes play all roles compatible with a possible involvement in NSHL phenotype, representing therefore excellent candidates for subsequent studies on NSHL patients and controls.

However two main limitations of this kind of approach should also be taken into account when looking at these data, both concerning the usage of Gene Ontology annotations to build the gene profiles on which the semantic similarity is measured. One is linked to the current knowledge about the human genome and its content in terms of genes. Indeed, the only prerequisite for a gene to be prioritized by our gene scoring system for a given disease is that of being annotated with at least one GO term, but, as clearly evidenced in this study, we are still far from the complete annotation of the entire human genome, as we were forced to exclude almost half of the possible candidates since they completely lacked GO annotations. This limitation obviously biases the results towards the best studied genes; however it will be progressively overcome in the future, due to the daily updates in this field. The other limitation regards the nature of the associations between GO terms and gene products. All the associations in Gene Ontology fall in five general categories indicating the evidences that support the annotation of a gene to a specific term. Four of these categories comprise exclusively manually-curated associations supported by experimental, computational analysis, author statements or curatorial statements. Unfortunately the great majority of GO associations does not fall in any of these manually-curated categories, being inferred from electronic annotation (IEA), which may open a debate on how reliable and precise they are. At present, given the high percentage of IEA associations in GO, it is not conceivable to discard them and consider only those manually-curated. Other solutions must therefore be devised to address this issue. Future developments of our gene scoring system could for instance take into account this problem by down-weighting the IEA associations respect to those manually-curated. However the quantification of the difference in weight between the manual and electronic associations is not trivial and requires an accurate study of the algorithms behind the electronic associations. We reserve in future to enhance our algorithm in this direction.

Final and essential step to confirm the results presented in this study is however the experimental validation. To this end two main aspects should be taken into account: (i) the accurate study and selection among the top-ranked genes of the most intringuing candidates for NSHL; we think for instance that the first one (WDR1) represents a good starting point, due to both what is known about its functions and structure – it is indeed involved in the organization of the actin, fundamental for the auditory process, and small enough to be quite easily sequenced in a large number of subjects; (ii) an equally accurate selection of the appropriate NSHL patients and controls to be screened for causative mutations; it is advisable for instance to screen these genes on a cohort of patients already excluded to carry mutations in GJB2, due to the high incidence of NSHL cases caused by mutations in this gene, and on a control set appropriately matched for their geographic origin, in order to take into account the geographic distribution of the human DNA sequence variation.

Methods

A total of 15727 genes (candidate genes) were prioritized for NSHL in this study. We chose as candidate genes all the genes contained in the NSHL susceptibility loci known so far (Tables S1, S2, S3 respectively for NSHL autosomal dominant, autosomal recessive and X-linked, Y-linked and modifier loci), so that all evidences coming from previous linkage analysis studies were taken into account.

We drew the complete disease gene list starting from the Hereditary Hearing Loss Homepage [9] and a team of experts (geneticists and molecular biologists) further analysed the literature to find additional advances in the field by performing multiple queries on PubMed. To the best of our knowledge, 51 genes belong to this category, as reported in Table 1.

For each disease and candidate gene, we extracted all their GO annotations using the file gene2go downloaded on 29th May 2009 from NCBI Entrez Gene ftp site [39]. One out of fifty-one disease genes – MIRN96 – had no GO annotations, therefore it was not included in this study, consequently narrowing the disease gene list to fifty genes. Likewise, 6987 out of 15727 candidate genes had no GO annotations, therefore the candidate gene list was consequently narrowed to 8740 genes.

Semantic similarity between two genes

As a node-based approach, our metric computes the similarity between two genes by comparing the GO terms describing them, their ancestors, and their descendants in the GO network. It is based on the Information Content (IC), which gives a measure of how specific and informative a term is. The IC of a term is quantified as the negative log likelihood(1)where is the probability of occurrence of in a specific corpus, which is normally estimated by the frequency of annotation of the term and its children in the GO structure [24], [40].

The concept of IC was used by Lin to quantify the semantic similarity between two terms in a tree-structured ontology, measuring the information they share normalized respect to the information contained in their total descriptions. According to Lin's metric [25] the similarity between two terms and is defined as:(2)where is the most informative common ancestor of the terms and , i.e. the common ancestor with the smallest .

However, two aspects of this metric limit its application:

* it applies only to trees, where a unique most informative common ancestor between two any given concepts exists;

* it measures the distance between single terms rather than set of terms.

For the first drawback, it is well known that in the case of a direct acyclic graph (DAG), such as GO, two terms can share parents by multiple paths, as multiple parents for each concept are allowed. Therefore, we chose, as , the minimum subsumer between and along all their independent paths to the graph root [13].

To address the second issue, we defined our Semantic Similarity Measure (SSM) by directly extending Lin's idea to quantify the similarity between two concepts to the comparison between two gene products, i.e. two sets of concepts, therefore measuring the IC of the common description of the two gene products, normalized respect to the IC of their global description.

Let betwo gene products, annotated respectively with and GO terms, that are to be compared. The idea is that each term in is an independent view of the gene and has to be compared with its counterpart in the gene annotation, namely the term in with maximum IC for the common description respect to it. In formulas, for the term , its counterpart in is defined as:The IC for each GO term c is estimated using its probability of occurrence in the corpus of all gene annotations provided by ENGINE database [41], [42]: in details, the probability is calculated for every term by counting the number of gene products associated with the term or any of its children, divided by the number of total associations between the GO terms and gene products.

Considering as independent the single views of a gene offered by each of its terms, the semantic similarity of respect to is estimated by the sum of the shared common description ICs between each term in and its counterpart in normalized with the IC of their global description:

The similarity of the gene respect to () is obtained by inverting the roles of and in the above formula. Finally, we defined our Semantic Similarity Measure between and , , as the mean between the similarity of respect to and the similarity of respect to :(3) generates normalized similarity values between 0 and 1: it's equal to 0 for genes annotated with terms that share only the root and equal to 1 for genes annotated with the same terms.

Validation of the for NSHL gene prioritization

A cross-validation procedure was used to check the reliability of the ranking of candidate genes for their involvement in NSHL. A random set of 8 disease genes was added to the set of candidate genes for 10000 times. Each time the values for this enlarged set of candidates were computed against the remaining disease genes and the number of disease genes in the first top-ranked positions was counted. The corresponding distributions of these countings were then compared with the probabilities of counting disease genes when a random drawn of 100, 75, 50, 8 genes, respectively, was performed from a set of 8748 genes (8740 candidate genes plus 8 disease genes): in the case of random drawns, the countings are described by a hypergeometric distribution with successes for d draws without replacement.

More in details, we computed the p-value and the power of a statistical test on the hypothesis of equal distributions against the hypothesis of a greater number of disease genes in the first positions for the ranked ordering respect to the random ordering. The p-value measures the probability to obtain, by random extraction, a number of disease genes greater than the mean value of the number of disease genes found in the top-ranked positions on the 10000 cross-validations:(4)The p-value is our estimate of the probability of rejecting when is true: whenever the p-value was less than the significance level , we maintained that the number of disease genes found in the top-ranked positions was statistically significantly greater than that found in random orderings.

The knowledge of the empirical distribution of estimated through the cross-validation procedure, allowed us to estimate the power of the test with level : indicating with the -quantile of the hypergeometric distribution, is computed as follows:(5)where is the number of disease genes found in the first positions for the randomization. The larger is the percentage of values obtained in cross-validation that are greater than , the more effective is the gene prioritization system.

Functional characterization of candidate and disease genes for NSHL

Candidate genes.

The statistical test used to identify the most representative GO terms associated with the candidate genes was designed as follows: the null hypothesis is that candidate genes annotated with a particular GO category have an average score equal to the average score expected for a random list of candidate genes with the same size, whereas the alternative hypothesis is that the GO category list has a higher average score and, therefore, is supposed to be more associated with the disease than a random candidate gene list. After selecting all the GO terms associated with all the candidate genes, we computed a p-value which scores each GO term according to the following strategy: the higher is the number of its associated candidates which obtained in our ranking a high value, the more the GO term is considered enriched in the candidate gene list. This implies that, choosing a significance threshold of 0.01, the GO terms with p-value can be considered significantly descriptive of the best candidate genes and consequently significantly associated with the disease. This provides directions for the NSHL researchers about the functions to be more deeply investigated in future laboratory experiments.

The p-value for the i-th GO term is computed as follows:(6)where is the average of the scores resulting for candidate genes annotated with , is the number of candidate genes annotated with the category, is the empirical cumulative distribution for the scores, averaged on lists of candidate genes of size . was computed by drawing 10000 random lists of candidate genes of size and averaging the respective gene scores.

Disease genes.

After selecting all the GO terms used to annotate the disease genes, we computed for each GO term a Fisher's exact test p-value which scores the GO category (GO Term C) highly for enrichment if many more disease genes than expected belong to the category. The contingency table (Table 6) is built by counting the disease and non-disease genes associated and not associated with the GO category.

[Figure omitted. See PDF.]

Table 6. Contingency table.

https://doi.org/10.1371/journal.pone.0012742.t006

The definition of the non-disease class is not trivial, as it is not possible to know in advance which candidate genes will be discovered as responsible for NSHL in the future – i.e. it is not possible to discriminate disease and non-disease genes among the candidates. To address this issue we decided to use the distribution of SSM scores in the class of candidate genes to define the non-disease class. We considered as non-disease genes the candidate genes with a score less than the 95th percentile of the distribution of candidate gene scores.

The GO terms with a Fisher's test p-value smaller than 0.01 are considered significantly over-represented in the list of the disease genes. These provide indications about the main functions and biological processes involved in the hearing mechanisms, taking into account the scores computed for our candidate gene list against the NSHL genes at present known.

score.

For both candidate and disease gene lists their over-represented GO terms are weighted taking into account their specificity in the corpus of the GO annotations as follows:where IC is estimated using the probability of occurrence of the GO terms in the corpus of all gene annotations provided by ENGINE database [41], [42].

Supporting Information

[Figure omitted. See PDF.]

Table S1.

NSHL autosomal dominant loci. Locus names and chromosomal locations have been inferred from literature. References are relative to the articles where the locus association to NSHL was identified.

https://doi.org/10.1371/journal.pone.0012742.s001

(0.05 MB PDF)

Table S2.

NSHL autosomal recessive loci. Locus names and chromosomal locations have been inferred from literature. References are relative to the articles where the locus association to NSHL was identified.

https://doi.org/10.1371/journal.pone.0012742.s002

(0.05 MB PDF)

Table S3.

NSHL X-linked, Y-linked and modifier loci. Locus names and chromosomal locations have been inferred from literature. References are relative to the articles where the locus association to NSHL was identified.

https://doi.org/10.1371/journal.pone.0012742.s003

(0.02 MB PDF)

Acknowledgments

We thank Andreas Gisel and Angelica Tulipano from CNR – ITB Bari for useful discussions and for providing probabilities for IC estimation. We thank Luigi Doronzo from IBM for useful insights and revision of the manuscript. Finally, we thank the members of the research group at Medical Genetics Service, Hospital CSS-IRCCS, San Giovanni Rotondo, Italy, Massimo Carella and Maria Stella Alemanno for useful discussions around the pathology and the gene validation protocol.

Author Contributions

Conceived and designed the experiments: MA TMC CS GS PL. Performed the experiments: MA TMC CS. Analyzed the data: MA TMC CS AV GS. Wrote the paper: MA TMC CS. Implemented the mathematical model used for gene prioritization: SB.

Citation: Accetturo M, Creanza TM, Santoro C, Tria G, Giordano A, Battagliero S, et al. (2010) Finding New Genes for Non-Syndromic Hearing Loss through an In Silico Prioritization Study. PLoS ONE5(9): e12742. https://doi.org/10.1371/journal.pone.0012742

References

1. Schrijver I (2004) Hereditary non-syndromic sensorineural hearing loss: transforming silence to sound. J Mol Diagn 6: 275–284.I. Schrijver2004Hereditary non-syndromic sensorineural hearing loss: transforming silence to sound.J Mol Diagn6275284

2. Kelsell DP, Dunlop J, Stevens HP, Lench NJ, Liang JN, et al. (2009) Connexin 26 mutations in hereditary non-syndromic sensorineural deafness. Nature 387: 80–83.DP KelsellJ. DunlopHP StevensNJ LenchJN Liang2009Connexin 26 mutations in hereditary non-syndromic sensorineural deafness.Nature3878083

3. Grifa A, Wagner CA, D'Ambrosio L, Melchionda S, Bernardi F, et al. (1999) Mutations in GJB6 cause nonsyndromic autosomal dominant deafness at DFNA3 locus. Nat Genet 23: 16–18.A. GrifaCA WagnerL. D'AmbrosioS. MelchiondaF. Bernardi1999Mutations in GJB6 cause nonsyndromic autosomal dominant deafness at DFNA3 locus.Nat Genet231618

4. Li XC, Everett LA, Lalwani AK, Desmukh D, Friedman TB, et al. (1998) A mutation in PDS causes non-syndromic recessive deafness. Nat Genet 18: 215–217.XC LiLA EverettAK LalwaniD. DesmukhTB Friedman1998A mutation in PDS causes non-syndromic recessive deafness.Nat Genet18215217

5. Robertson NG, Lu L, Heller S, Merchant SN, Eavey RD, et al. (1998) Mutations in a novel cochlear gene cause DFNA9, a human nonsyndromic deafness with vestibular dysfunction. Nat Genet 20: 299–303.NG RobertsonL. LuS. HellerSN MerchantRD Eavey1998Mutations in a novel cochlear gene cause DFNA9, a human nonsyndromic deafness with vestibular dysfunction.Nat Genet20299303

6. De Kok YJ, Van der Maarel SM, Bitner-Glindzicz M, Huber I, Monaco AP, et al. (1995) Association between X-linked mixed deafness and mutations in the POU domain gene POU3F4. Science 267: 685–688.YJ De KokSM Van der MaarelM. Bitner-GlindziczI. HuberAP Monaco1995Association between X-linked mixed deafness and mutations in the POU domain gene POU3F4.Science267685688

7. Bespalova IN, Van Camp G, Bom SJ, Brown DJ, Cryns K, et al. (2001) Mutations in the Wolfram syndrome 1 gene (WFS1) are a common cause of low frequency sensorineural hearing loss. Hum Mol Genet 10: 2501–2508.IN BespalovaG. Van CampSJ BomDJ BrownK. Cryns2001Mutations in the Wolfram syndrome 1 gene (WFS1) are a common cause of low frequency sensorineural hearing loss.Hum Mol Genet1025012508

8. Young TL, Ives E, Lynch E, Person R, Snook S, et al. (2001) Non-syndromic progressive hearing loss DFNA38 is caused by heterozygous missense mutation in the Wolfram syndrome gene WFS1. Hum Mol Genet 10: 2509–2514.TL YoungE. IvesE. LynchR. PersonS. Snook2001Non-syndromic progressive hearing loss DFNA38 is caused by heterozygous missense mutation in the Wolfram syndrome gene WFS1.Hum Mol Genet1025092514

9. Hereditary Hearing Loss home page. Hereditary Hearing Loss home page.Available: http://herditaryhearingloss.org/. Accessed April 2010. Available: http://herditaryhearingloss.org/. Accessed April 2010.

10. Linghu B, Snitkin ES, Hu Z, Xia Y, DeLisi C (2009) Genome-wide prioritization of disease genes and identification of disease-disease associations from an integrated human functional linkage network. Genome Biology 10: R91.B. LinghuES SnitkinZ. HuY. XiaC. DeLisi2009Genome-wide prioritization of disease genes and identification of disease-disease associations from an integrated human functional linkage network.Genome Biology10R91

11. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, et al. (2000) Gene Ontology: tool for the unification of biology. the Gene Ontology Consortium. Nat Genet 25: 25–29.M. AshburnerCA BallJA BlakeD. BotsteinH. Butler2000Gene Ontology: tool for the unification of biology. the Gene Ontology Consortium.Nat Genet252529

12. Pesquita C, Faria D, Falcao AO, Lord P, Couto F (2009) Semantic similarity in biomedical ontologies. PLoS Comput Biol 5: e1000443.C. PesquitaD. FariaAO FalcaoP. LordF. Couto2009Semantic similarity in biomedical ontologies.PLoS Comput Biol5e1000443

13. Lord PW, Stevens RD, Brass A, Goble CA (2003) Semantic similarity measures as tools for exploring the Gene Ontology. Pacific Symposium on Biocomputing 8: 601–612.PW LordRD StevensA. BrassCA Goble2003Semantic similarity measures as tools for exploring the Gene Ontology.Pacific Symposium on Biocomputing8601612

14. Wu H, Su Z, Mao F, Olman V, Xu Y (2005) Prediction of functional modules based on comparative genome analysis and gene ontology application. Nucleic Acids Res 33: 2822–2837.H. WuZ. SuF. MaoV. OlmanY. Xu2005Prediction of functional modules based on comparative genome analysis and gene ontology application.Nucleic Acids Res3328222837

15. Wu X, Zhu L, Guo J, Zhang DY, Lin K (2006) Prediction of yeast protein-protein interaction network: insights from the gene ontology and annotations. Nucleic Acids Res 34: 2137–2150.X. WuL. ZhuJ. GuoDY ZhangK. Lin2006Prediction of yeast protein-protein interaction network: insights from the gene ontology and annotations.Nucleic Acids Res3421372150

16. Yu H, Gao L, Tu K, Guo Z (2005) Broadly predicting specific gene functions with expression similarity and taxonomy similarity. Gene 352: 75–81.H. YuL. GaoK. TuZ. Guo2005Broadly predicting specific gene functions with expression similarity and taxonomy similarity.Gene3527581

17. Cheng J, Cline M, Martin J, Finkelstein D, Awad T, et al. (2004) A knowledge-based clustering algorithm driven by gene ontology. Journal of Biopharmaceutical Statistics 14: 687–700.J. ChengM. ClineJ. MartinD. FinkelsteinT. Awad2004A knowledge-based clustering algorithm driven by gene ontology.Journal of Biopharmaceutical Statistics14687700

18. del Pozo A, Pazos F, Valencia A (2008) Defining functional distances over gene ontology. BMC Bioinformatics 9: 50.A. del PozoF. PazosA. Valencia2008Defining functional distances over gene ontology.BMC Bioinformatics950

19. Resnik P (1995) Using information content to evaluate semantic similarity in a taxonomy. Proc of the 14th International Joint Conference on Artificial Intelligence 448–453.P. Resnik1995Using information content to evaluate semantic similarity in a taxonomy.Proc of the 14th International Joint Conference on Artificial Intelligence448453

20. Couto FM, Silva MJ, Coutinho PM (2005) Semantic similarity over the gene ontology: Family correlation and selecting disjunctive ancestors. pp. 343–344.FM CoutoMJ SilvaPM Coutinho2005Semantic similarity over the gene ontology: Family correlation and selecting disjunctive ancestors.343344In: Proc. of the ACM Conference in Information and Knowledge Management as a short paper. Bremen, Germany.

21. Othman R, Deris S, Illias R (2007) A genetic similarity algorithm for searching the gene ontology terms and annotating anonymous protein sequences. J Biomed Inform 23: 529–538.R. OthmanS. DerisR. Illias2007A genetic similarity algorithm for searching the gene ontology terms and annotating anonymous protein sequences.J Biomed Inform23529538

22. Jiang J, Conrath D (1997) Semantic similarity based on corpus statistics and lexical taxonomy. pp. 19–33.J. JiangD. Conrath1997Semantic similarity based on corpus statistics and lexical taxonomy.1933In: Proc. of the 10th International Conference on Research on Computational Linguistics. Taiwan.

23. Schlicker A, Domingues FS, Rahnenführer J, Lengauer T (2006) A new measure for functional similarity of gene products based on gene ontology. BMC Bioinformatics 7: 302.A. SchlickerFS DominguesJ. RahnenführerT. Lengauer2006A new measure for functional similarity of gene products based on gene ontology.BMC Bioinformatics7302

24. Bodenreider O, Aubry M, Burgun A (2005) Non-lexical approaches to identifying associative relations in the gene ontology. Pac Symp Biocomput 7: 91–102.O. BodenreiderM. AubryA. Burgun2005Non-lexical approaches to identifying associative relations in the gene ontology.Pac Symp Biocomput791102

25. Lin D (1998) An information-theoretic definition of similarity. Proc. of the 15th International Conference on Machine Learning. San Francisco, CA: Morgan Kaufmann. pp. 296–304.D. Lin1998An information-theoretic definition of similarity.Proc. of the 15th International Conference on Machine LearningSan Francisco, CAMorgan Kaufmann296304

26. Sevilla JL, Segura V, Podhorski A, Guruceaga E, Mato JM, et al. (2005) Correlation between gene expression and go semantic similarity. IEEE/ACM Transactions on Computational Biology and Bioinformatics 2: 330–338.JL SevillaV. SeguraA. PodhorskiE. GuruceagaJM Mato2005Correlation between gene expression and go semantic similarity.IEEE/ACM Transactions on Computational Biology and Bioinformatics2330338

27. Riensche RM, Baddeley BL, Sanfilippo AP, Posse C, Gopalan B (2007) Xoa: Web-enabled cross-ontological analytics. Services, IEEE Congress 99–105.RM RienscheBL BaddeleyAP SanfilippoC. PosseB. Gopalan2007Xoa: Web-enabled cross-ontological analytics.Services, IEEE Congress99105

28. Wang JZZ, Du Z, Payattakool R, Yu PSS, Chen CFF (2007) A new method to measure the semantic similarity of go terms. Bioinformatics 23: 1274–1281.JZZ WangZ. DuR. PayattakoolPSS YuCFF Chen2007A new method to measure the semantic similarity of go terms.Bioinformatics2312741281

29. Lei Z, Dai Y (2006) Assessing protein similarity with gene ontology and its use in subnuclear localization prediction. BMC Bioinformatics 7: 491.Z. LeiY. Dai2006Assessing protein similarity with gene ontology and its use in subnuclear localization prediction.BMC Bioinformatics7491

30. R Package MatchingR Package MatchingAvailable: http://sekhon.berkeley.edu/matching. Available: http://sekhon.berkeley.edu/matching.

31. Weil D, El-Amraoui A, Masmoudi S, Mustapha M, Kikkawa Y, et al. (2003) Usher syndrome type IG (USH1G) is caused by mutations in the gene encoding SANS, a protein that associates with the USH1C protein, harmonin. Hum Molec Genet 12: 463–471.D. WeilA. El-AmraouiS. MasmoudiM. MustaphaY. Kikkawa2003Usher syndrome type IG (USH1G) is caused by mutations in the gene encoding SANS, a protein that associates with the USH1C protein, harmonin.Hum Molec Genet12463471

32. Piccini M, Vitelli F, Seri M, Galietta LJV, Moran O, et al. (1999) KCNE1-like gene is deleted in AMME contiguous gene syndrome: Identification and characterization of the human and mouse homologs. Genomics 60: 251–257.M. PicciniF. VitelliM. SeriLJV GaliettaO. Moran1999KCNE1-like gene is deleted in AMME contiguous gene syndrome: Identification and characterization of the human and mouse homologs.Genomics60251257

33. Jin H, May M, Tranebjaerg L, Kendall E, Fontan G, et al. (1996) A novel X-linked gene, DDP, shows mutations in families with deafness (DFN-1), dystonia, mental deficiency and blindness. Nature Genet 14: 177–180.H. JinM. MayL. TranebjaergE. KendallG. Fontan1996A novel X-linked gene, DDP, shows mutations in families with deafness (DFN-1), dystonia, mental deficiency and blindness.Nature Genet14177180

34. Tranebjaerg L, van Ghelue M, Nilssen O, Hodes ME, Dlouhy SR, et al. (1997) Jensen syndrome is allelic to Mohr-Tranebjaerg syndrome and both are caused by stop mutations in the DDP gene. (Abstract) Am J Hum Genet 61: suppl.A349 only.L. TranebjaergM. van GhelueO. NilssenME HodesSR Dlouhy1997Jensen syndrome is allelic to Mohr-Tranebjaerg syndrome and both are caused by stop mutations in the DDP gene.(Abstract) Am J Hum Genet61suppl.A349 only

35. Collin GB, Marshall JD, Ikeda A, So WV, Russell-Eggitt I, et al. (2002) Mutations in ALMS1 cause obesity, type 2 diabetes and neurosensory degeneration in Alström syndrome. Nature Genet 31: 74–78.GB CollinJD MarshallA. IkedaWV SoI. Russell-Eggitt2002Mutations in ALMS1 cause obesity, type 2 diabetes and neurosensory degeneration in Alström syndrome.Nature Genet317478

36. Verpy E, Leibovici M, Zwaenepoel I, Liu XZ, Gal A, et al. (2000) A defect in harmonin, a PDZ domain-containing protein expressed in the inner ear sensory hair cells, underlies Usher syndrome type 1C. Nature Genet 26: 51–55.E. VerpyM. LeiboviciI. ZwaenepoelXZ LiuA. Gal2000A defect in harmonin, a PDZ domain-containing protein expressed in the inner ear sensory hair cells, underlies Usher syndrome type 1C.Nature Genet265155

37. Ouyang XM, Xia XJ, Verpy E, Du LL, Pandya A, et al. (2002) Mutations in the alternatively spliced exons of USH1C cause non-syndromic recessive deafness. Hum Genet 111: 26–30.XM OuyangXJ XiaE. VerpyLL DuA. Pandya2002Mutations in the alternatively spliced exons of USH1C cause non-syndromic recessive deafness.Hum Genet1112630

38. Ahmed ZM, Smith TN, Riazuddin S, Makishima T, Ghosh M, et al. (2002) Nonsyndromic recessive deafness DFNB18 and Usher syndrome type 1C are allelic mutations of USH1C. Hum Genet 110: 527–531.ZM AhmedTN SmithS. RiazuddinT. MakishimaM. Ghosh2002Nonsyndromic recessive deafness DFNB18 and Usher syndrome type 1C are allelic mutations of USH1C.Hum Genet110527531

39. Maglott D, Ostell J, Pruitt KD, Tatusova T (2007) Entrez Gene: gene-centered information at NCBI. Nucleic Acids Res 35: D26–31.D. MaglottJ. OstellKD PruittT. Tatusova2007Entrez Gene: gene-centered information at NCBI.Nucleic Acids Res35D2631

40. Tlili A, Männikkö M, Charfedine I, Lahmar I, Benzina Z, et al. (2005) A novel autosomal recessive non-syndromic deafness locus, DFNB66, maps to chromosome 6p21.2–22.3 in a large Tunisian consanguineous family. Hum Hered 60: 123–128.A. TliliM. MännikköI. CharfedineI. LahmarZ. Benzina2005A novel autosomal recessive non-syndromic deafness locus, DFNB66, maps to chromosome 6p21.2–22.3 in a large Tunisian consanguineous family.Hum Hered60123128

41. Tulipano A, Donvito G, Licciulli F, Maggi G, Gisel A (2007) Gene analogue finder: a GRID solution for finding functionally analogous gene products. BMC Bioinformatics 8: 329.A. TulipanoG. DonvitoF. LicciulliG. MaggiA. Gisel2007Gene analogue finder: a GRID solution for finding functionally analogous gene products.BMC Bioinformatics8329

42. GENe AnaloGue FINdEr ENGINE database: A repository for functional analogous gene products. GENe AnaloGue FINdEr ENGINE database: A repository for functional analogous gene products.Available: http://spank.ba.itb.cnr.it/engine/. Accessed June 2009. Available: http://spank.ba.itb.cnr.it/engine/. Accessed June 2009.

43. Lynch ED, Lee MK, Morrow JE, Welcsh PL, León PE, et al. (1997) Nonsyndromic deafness DFNA1 associated with mutation of a human homolog of the Drosophila gene diaphanous. Science 278: 1315–1318.ED LynchMK LeeJE MorrowPL WelcshPE León1997Nonsyndromic deafness DFNA1 associated with mutation of a human homolog of the Drosophila gene diaphanous.Science27813151318

44. Xia JH, Liu CY, Tang BS, Pan Q, Huang L, et al. (1998) Mutations in the gene encoding gap junction protein beta-3 associated with autosomal dominant hearing impairment. Nature Genet 20: 370–373.JH XiaCY LiuBS TangQ. PanL. Huang1998Mutations in the gene encoding gap junction protein beta-3 associated with autosomal dominant hearing impairment.Nature Genet20370373

45. Kubisch C, Schroeder BC, Friedrich T, Lütjohann B, El-Amraoui A, et al. (1999) KCNQ4, a novel potassium channel expressed in sensory outer hair cells, is mutated in dominant deafness. Cell 96: 437–446.C. KubischBC SchroederT. FriedrichB. LütjohannA. El-Amraoui1999KCNQ4, a novel potassium channel expressed in sensory outer hair cells, is mutated in dominant deafness.Cell96437446

46. del Castillo I, Villamar M, Moreno-Pelayo MA, del Castillo FJ, Alvarez A, et al. (2002) A deletion involving the connexin 30 gene in nonsyndromic hearing impairment. N Engl J Med 346: 243–249.I. del CastilloM. VillamarMA Moreno-PelayoFJ del CastilloA. Alvarez2002A deletion involving the connexin 30 gene in nonsyndromic hearing impairment.N Engl J Med346243249

47. Donaudy F, Snoeckx R, Pfister M, Zenner HP, Blin N, et al. (2004) Nonmuscle myosin heavy-chain gene MYH14 is expressed in cochlea and mutated in patients affected by autosomal dominant hearing impairment (DFNA4). Am J of Hum Genet 74: 770–776.F. DonaudyR. SnoeckxM. PfisterHP ZennerN. Blin2004Nonmuscle myosin heavy-chain gene MYH14 is expressed in cochlea and mutated in patients affected by autosomal dominant hearing impairment (DFNA4).Am J of Hum Genet74770776

48. Van Laer L, Huizing EH, Verstreken M, van Zuijlen D, Wauters JG, et al. (1998) Nonsyndromic hearing impairment is associated with a mutation in DFNA5. Nature Genet 20: 194–197.L. Van LaerEH HuizingM. VerstrekenD. van ZuijlenJG Wauters1998Nonsyndromic hearing impairment is associated with a mutation in DFNA5.Nature Genet20194197

49. Verhoeven K, Van Laer L, Kirschhofer K, Legan PK, Hughes DC, et al. (1998) Mutations in the human alpha-tectorin gene cause autosomal dominant non-syndromic hearing impairment. Nature Genet 19: 60–62.K. VerhoevenL. Van LaerK. KirschhoferPK LeganDC Hughes1998Mutations in the human alpha-tectorin gene cause autosomal dominant non-syndromic hearing impairment.Nature Genet196062

50. Wayne S, Robertson NG, DeClau F, Chen N, Verhoeven K (2001) Mutations in the transcriptional activator EYA4 cause late-onset deafness at the DFNA10 locus. Hum Mol Genet 10: 195–200.S. WayneNG RobertsonF. DeClauN. ChenK. Verhoeven2001Mutations in the transcriptional activator EYA4 cause late-onset deafness at the DFNA10 locus.Hum Mol Genet10195200

51. Liu XZ, Walsh J, Mburu P, Kendrick-Jones J, Cope MJ, et al. (1997) Mutations in the myosin VIIA gene cause non-syndromic recessive deafness. Nat Genet 16: 188–190.XZ LiuJ. WalshP. MburuJ. Kendrick-JonesMJ Cope1997Mutations in the myosin VIIA gene cause non-syndromic recessive deafness.Nat Genet16188190

52. McGuirt WT, Prasad SD, Griffith AJ, Kunst HP, Green GE, et al. (1999) Mutations in COL11A2 cause non-syndromic hearing loss (DFNA13). Hum Mol Genet 23: 413–419.WT McGuirtSD PrasadAJ GriffithHP KunstGE Green1999Mutations in COL11A2 cause non-syndromic hearing loss (DFNA13).Hum Mol Genet23413419

53. Vahava O, Morell R, Lynch ED, Weiss S, Kagan ME, et al. (1998) Mutation in transcription factor POU4F3 associated with inherited progressive hearing loss in humans. Science 279: 1950–1954.O. VahavaR. MorellED LynchS. WeissME Kagan1998Mutation in transcription factor POU4F3 associated with inherited progressive hearing loss in humans.Science27919501954

54. Lalwani AK, Goldstein JA, Kelley MJ, Luxford W, Castelein CM, et al. (2000) Human nonsyndromic hereditary deafness DFNA17 is due to a mutation in nonmuscle myosin MYH9. Am J Hum Genet 67: 1121–1128.AK LalwaniJA GoldsteinMJ KelleyW. LuxfordCM Castelein2000Human nonsyndromic hereditary deafness DFNA17 is due to a mutation in nonmuscle myosin MYH9.Am J Hum Genet6711211128

55. Zhu M, Yang T, Wei S, DeWan AT, Morell RJ, et al. (2003) Mutations in the gamma-actin gene (ACTG1) are associated with dominant progressive deafness (DFNA20/26). Am J Hum Genet 73: 1082–1091.M. ZhuT. YangS. WeiAT DeWanRJ Morell2003Mutations in the gamma-actin gene (ACTG1) are associated with dominant progressive deafness (DFNA20/26).Am J Hum Genet7310821091

56. van Wijk E, Krieger E, Kemperman MH, De Leenheer EM, Huygen PL, et al. (2003) A mutation in the gamma actin 1 (ACTG1) gene causes autosomal dominant hearing loss (DFNA20/26). J Med Genet 40: 879–884.E. van WijkE. KriegerMH KempermanEM De LeenheerPL Huygen2003A mutation in the gamma actin 1 (ACTG1) gene causes autosomal dominant hearing loss (DFNA20/26).J Med Genet40879884

57. Melchionda S, Ahituv N, Bisceglia L, Sobe T, Glaser F, et al. (2001) MYO6, the human homologue of the gene responsible for deafness in Snell's waltzer mice, is mutated in autosomal dominant nonsyndromic hearing loss. Am J Hum Genet 69: 635–640.S. MelchiondaN. AhituvL. BiscegliaT. SobeF. Glaser2001MYO6, the human homologue of the gene responsible for deafness in Snell's waltzer mice, is mutated in autosomal dominant nonsyndromic hearing loss.Am J Hum Genet69635640

58. Peters LM, Anderson DW, Griffith AJ, Grundfast KM, San Agustin TB, et al. (2002) Mutation of a transcription factor, TFCP2L3, causes progressive autosomal dominant hearing loss, DFNA28. Hum Mol Genet 11: 2877–2885.LM PetersDW AndersonAJ GriffithKM GrundfastTB San Agustin2002Mutation of a transcription factor, TFCP2L3, causes progressive autosomal dominant hearing loss, DFNA28.Hum Mol Genet1128772885

59. Kurima K, Peters LM, Yang Y, Riazuddin S, Ahmed ZM, et al. (2002) Dominant and recessive deafness caused by mutations of a novel gene, TMC1, required for cochlear hair-cell function. Nat Genet 30: 277–284.K. KurimaLM PetersY. YangS. RiazuddinZM Ahmed2002Dominant and recessive deafness caused by mutations of a novel gene, TMC1, required for cochlear hair-cell function.Nat Genet30277284

60. Abe S, Katagiri T, Saito-Hisaminato A, Usami S, Inoue Y (2003) Identification of CRYM as a candidate responsible for nonsyndromic deafness, through cDNA microarray analysis of human cochlear and vestibular tissues. Am J Hum Genet 72: 73–82.S. AbeT. KatagiriA. Saito-HisaminatoS. UsamiY. Inoue2003Identification of CRYM as a candidate responsible for nonsyndromic deafness, through cDNA microarray analysis of human cochlear and vestibular tissues.Am J Hum Genet727382

61. Modamio-Høybjør S, Mencia A, Goodyear R, del Castillo I, Richardson G, et al. (2007) A mutation in CCDC50, a gene encoding an effector of epidermal growth factor-mediated cell signaling, causes progressive hearing loss. Am J Hum Genet 80: 1076–1089.S. Modamio-HøybjørA. MenciaR. GoodyearI. del CastilloG. Richardson2007A mutation in CCDC50, a gene encoding an effector of epidermal growth factor-mediated cell signaling, causes progressive hearing loss.Am J Hum Genet8010761089

62. Donaudy F, Ferrara A, Esposito L, Hertzano R, Ben-David O, et al. (2003) Multiple mutations of MYO1A, a cochlear-expressed gene, in sensorineural hearing loss. Hum Mol Genet 72: 1571–1577.F. DonaudyA. FerraraL. EspositoR. HertzanoO. Ben-David2003Multiple mutations of MYO1A, a cochlear-expressed gene, in sensorineural hearing loss.Hum Mol Genet7215711577

63. Yang T, Gurrola JGn, Wu H, Chiu SM, Wangemann P, et al. (2009) Mutations of KCNJ10 together with mutations of SLC26A4 cause digenic nonsyndromic hearing loss associated with enlarged vestibular aqueduct syndrome. Am J Hum Genet 84: 651–657.T. YangJGn GurrolaH. WuSM ChiuP. Wangemann2009Mutations of KCNJ10 together with mutations of SLC26A4 cause digenic nonsyndromic hearing loss associated with enlarged vestibular aqueduct syndrome.Am J Hum Genet84651657

64. Mencía A, Modamio-Høybjør S, Redshaw N, Morín M, Mayo-Merino F (2009) Mutations in the seed region of human miR-96 are responsible for nonsyndromic progressive hearing loss. Nat Genet 41: 609–613.A. MencíaS. Modamio-HøybjørN. RedshawM. MorínF. Mayo-Merino2009Mutations in the seed region of human miR-96 are responsible for nonsyndromic progressive hearing loss.Nat Genet41609613

65. Wang A, Liang Y, Fridell RA, Probst FJ, Wilcox ER, et al. (1998) Association of unconventional myosin MYO15 mutations with human nonsyndromic deafness DFNB3. Science 280: 1447–1451.A. WangY. LiangRA FridellFJ ProbstER Wilcox1998Association of unconventional myosin MYO15 mutations with human nonsyndromic deafness DFNB3.Science28014471451

66. Naz S, Giguere CM, Kohrman DC, Mitchem KL, Riazuddin S, et al. (2002) Mutations in a novel gene, TMIE, are associated with hearing loss linked to the DFNB6 locus. Am J Hum Genet 71: 632–636.S. NazCM GiguereDC KohrmanKL MitchemS. Riazuddin2002Mutations in a novel gene, TMIE, are associated with hearing loss linked to the DFNB6 locus.Am J Hum Genet71632636

67. Scott HS, Kudoh J, Wattenhofer M, Shibuya K, Berry A, et al. (2002) Insertion of beta-satellite repeats identifies a transmembrane protease causing both congenital and childhood onset autosomal recessive deafness. Nat Genet 27: 59–63.HS ScottJ. KudohM. WattenhoferK. ShibuyaA. Berry2002Insertion of beta-satellite repeats identifies a transmembrane protease causing both congenital and childhood onset autosomal recessive deafness.Nat Genet275963

68. Yasunaga S, Grati M, Cohen-Salmon M, El-Amraoui A, Mustapha M, et al. (1999) A mutation in OTOF, encoding otoferlin, a FER-1-like protein, causes DFNB9, a nonsyndromic form of deafness. Nat Genet 21: 363–369.S. YasunagaM. GratiM. Cohen-SalmonA. El-AmraouiM. Mustapha1999A mutation in OTOF, encoding otoferlin, a FER-1-like protein, causes DFNB9, a nonsyndromic form of deafness.Nat Genet21363369

69. Bork JM, Peters LM, Riazuddin S, Bernstein SL, Ahmed ZM, et al. (2001) Usher syndrome 1D and nonsyndromic autosomal recessive deafness DFNB12 are caused by allelic mutations of the novel cadherin-like gene CDH23. Am J Hum Genet 68: 26–37.JM BorkLM PetersS. RiazuddinSL BernsteinZM Ahmed2001Usher syndrome 1D and nonsyndromic autosomal recessive deafness DFNB12 are caused by allelic mutations of the novel cadherin-like gene CDH23.Am J Hum Genet682637

70. Verpy E, Masmoudi S, Zwaenepoel I, Leibovici M, Hutchin TP, et al. (2001) Mutations in a new gene encoding a protein of the hair bundle cause non-syndromic deafness at the DFNB16 locus. Nat Genet 29: 345–349.E. VerpyS. MasmoudiI. ZwaenepoelM. LeiboviciTP Hutchin2001Mutations in a new gene encoding a protein of the hair bundle cause non-syndromic deafness at the DFNB16 locus.Nat Genet29345349

71. Zwaenepoel I, Mustapha M, Leibovici M, Verpy E, Goodyear R, et al. (2002) Otoancorin, an inner ear protein restricted to the interface between the apical surface of sensory epithelia and their overlying acellular gels, is defective in autosomal recessive deafness DFNB22. Proc Natl Acad Sci U S A 99: 6240–6245.I. ZwaenepoelM. MustaphaM. LeiboviciE. VerpyR. Goodyear2002Otoancorin, an inner ear protein restricted to the interface between the apical surface of sensory epithelia and their overlying acellular gels, is defective in autosomal recessive deafness DFNB22.Proc Natl Acad Sci U S A9962406245

72. Ahmed ZM, Riazuddin S, Ahmad J, Bernstein SL, Guo Y, et al. (2003) PCDH15 is expressed in the neurosensory epithelium of the eye and ear and mutant alleles are responsible for both USH1F and DFNB23. Hum Mol Genet 12: 3215–3223.ZM AhmedS. RiazuddinJ. AhmadSL BernsteinY. Guo2003PCDH15 is expressed in the neurosensory epithelium of the eye and ear and mutant alleles are responsible for both USH1F and DFNB23.Hum Mol Genet1232153223

73. Khan SY, Ahmed ZM, Shabbir MI, Kitajiri S, Kalsoom S, et al. (2007) Mutations of the RDX gene cause nonsyndromic hearing loss at the DFNB24 locus. Hum Mutat 28: 417–423.SY KhanZM AhmedMI ShabbirS. KitajiriS. Kalsoom2007Mutations of the RDX gene cause nonsyndromic hearing loss at the DFNB24 locus.Hum Mutat28417423

74. Shahin H, Walsh T, Sobe T, Abu Sa'ed J, Abu Rayan A, et al. (2006) Mutations in a novel isoform of TRIOBP that encodes a filamentous-actin binding protein are responsible for DFNB28 recessive nonsyndromic hearing loss. Am J Hum Genet 78: 144–152.H. ShahinT. WalshT. SobeJ. Abu Sa'edA. Abu Rayan2006Mutations in a novel isoform of TRIOBP that encodes a filamentous-actin binding protein are responsible for DFNB28 recessive nonsyndromic hearing loss.Am J Hum Genet78144152

75. Riazuddin S, Khan SN, Ahmed ZM, Ghosh M, Caution K, et al. (2006) Mutations in TRIOBP, which encodes a putative cytoskeletal-organizing protein, are associated with nonsyndromic recessive deafness. Am J Hum Genet 78: 137–143.S. RiazuddinSN KhanZM AhmedM. GhoshK. Caution2006Mutations in TRIOBP, which encodes a putative cytoskeletal-organizing protein, are associated with nonsyndromic recessive deafness.Am J Hum Genet78137143

76. Wilcox ER, Burton QL, Naz S, Riazuddin S, Smith TN, et al. (2001) Mutations in the gene encoding tight junction claudin-14 cause autosomal recessive deafness DFNB29. Cell 104: 165–172.ER WilcoxQL BurtonS. NazS. RiazuddinTN Smith2001Mutations in the gene encoding tight junction claudin-14 cause autosomal recessive deafness DFNB29.Cell104165172

77. Walsh T, Walsh V, Vreugde S, Hertzano R, Shahin H, et al. (2002) From flies' eyes to our ears: mutations in a human class iii myosin cause progressive nonsyndromic hearing loss DFNB30. Proc Natl Acad Sci U S A 99: 7518–7523.T. WalshV. WalshS. VreugdeR. HertzanoH. Shahin2002From flies' eyes to our ears: mutations in a human class iii myosin cause progressive nonsyndromic hearing loss DFNB30.Proc Natl Acad Sci U S A9975187523

78. Mburu P, Mustapha M, Varela A, Weil D, El-Amraoui A, et al. (2003) Defects in whirlin, a PDZ domain molecule involved in stereocilia elongation, cause deafness in the whirler mouse and families with DFNB31. Nat Genet 34: 421–428.P. MburuM. MustaphaA. VarelaD. WeilA. El-Amraoui2003Defects in whirlin, a PDZ domain molecule involved in stereocilia elongation, cause deafness in the whirler mouse and families with DFNB31.Nat Genet34421428

79. Ansar M, Din MA, Arshad M, Sohail M, Faiyaz-Ul-Haque M, et al. (2003) A novel autosomal recessive non-syndromic deafness locus (DFNB35) maps to 14q24.1–14q24.3 in large consanguineous kindred from Pakistan. Eur J Hum Genet 11: 77–80.M. AnsarMA DinM. ArshadM. SohailM. Faiyaz-Ul-Haque2003A novel autosomal recessive non-syndromic deafness locus (DFNB35) maps to 14q24.1–14q24.3 in large consanguineous kindred from Pakistan.Eur J Hum Genet117780

80. Collin RW, Kalay E, Tariq M, Peters T, van der Zwaag B, et al. (2008) Mutations of ESRRB encoding estrogen-related receptor beta cause autosomal-recessive nonsyndromic hearing impairment DFNB35. Am J Hum Genet 82: 125–138.RW CollinE. KalayM. TariqT. PetersB. van der Zwaag2008Mutations of ESRRB encoding estrogen-related receptor beta cause autosomal-recessive nonsyndromic hearing impairment DFNB35.Am J Hum Genet82125138

81. Naz S, Griffith AJ, Riazuddin S, Hampton LL, Battey JFJ, et al. (2004) Mutations of ESPN cause autosomal recessive deafness and vestibular dysfunction. J Med Genet 41: 591–595.S. NazAJ GriffithS. RiazuddinLL HamptonJFJ Battey2004Mutations of ESPN cause autosomal recessive deafness and vestibular dysfunction.J Med Genet41591595

82. Schultz JM, Khan SN, Ahmed ZM, Riazuddin S, Waryah AM, et al. (2009) Noncoding mutations of HGF are associated with nonsyndromic hearing loss, DFNB39. Am J Hum Genet 85: 25–39.JM SchultzSN KhanZM AhmedS. RiazuddinAM Waryah2009Noncoding mutations of HGF are associated with nonsyndromic hearing loss, DFNB39.Am J Hum Genet852539

83. Abe S, Usami S, Nakamura Y (2003) Mutations in the gene encoding KIAA1199 protein, an inner-ear protein expressed in deiters' cells and the fibrocytes, as the cause of nonsyndromic hearing loss. J Hum Genet 48: 564–570.S. AbeS. UsamiY. Nakamura2003Mutations in the gene encoding KIAA1199 protein, an inner-ear protein expressed in deiters' cells and the fibrocytes, as the cause of nonsyndromic hearing loss.J Hum Genet48564570

84. Riazuddin S, Ahmed ZM, Fanning AS, Lagziel A, Kitajiri S (2006) Tricellulin is a tight-junction protein necessary for hearing. Am J Hum Genet 79: 1040–1051.S. RiazuddinZM AhmedAS FanningA. LagzielS. Kitajiri2006Tricellulin is a tight-junction protein necessary for hearing.Am J Hum Genet7910401051

85. Delmaghani S, del Castillo FJ, Michel V, Leibovici M, Aghaie A (2006) Mutations in the gene encoding pejvakin, a newly identified protein of the afferent auditory pathway, cause DFNB59 auditory neuropathy. Nat Genet 38: 770–778.S. DelmaghaniFJ del CastilloV. MichelM. LeiboviciA. Aghaie2006Mutations in the gene encoding pejvakin, a newly identified protein of the afferent auditory pathway, cause DFNB59 auditory neuropathy.Nat Genet38770778

86. Liu XZ, Ouyang XM, Xia XJ, Zheng J, Pandya A, et al. (2003) Prestin, a cochlear motor protein, is defective in non-syndromic hearing loss. Hum Mol Genet 12: 1155–1162.XZ LiuXM OuyangXJ XiaJ. ZhengA. Pandya2003Prestin, a cochlear motor protein, is defective in non-syndromic hearing loss.Hum Mol Genet1211551162

87. Ahmed ZM, Masmoudi S, Kalay E, Belyantseva IA, Mosrati MA, et al. (2008) Mutations of LRTOMT, a fusion gene with alternative reading frames, cause nonsyndromic deafness in humans. Nat Genet 40: 1335–40.ZM AhmedS. MasmoudiE. KalayIA BelyantsevaMA Mosrati2008Mutations of LRTOMT, a fusion gene with alternative reading frames, cause nonsyndromic deafness in humans.Nat Genet40133540

88. Shabbir MI, Ahmed ZM, Khan SY, Riazuddin S, Waryah AM (2006) Mutations of human TMHS cause recessively inherited non-syndromic hearing loss. J Med Genet 43: 634–640.MI ShabbirZM AhmedSY KhanS. RiazuddinAM Waryah2006Mutations of human TMHS cause recessively inherited non-syndromic hearing loss.J Med Genet43634640

89. Kalay E, Li Y, Uzumcu A, Uyguner O, Collin RW (2006) Mutations in the lipoma HMGIC fusion partner-like 5 (LHFPL5) gene cause autosomal recessive nonsyndromic hearing loss. Hum Mutat 27: 633–639.E. KalayY. LiA. UzumcuO. UygunerRW Collin2006Mutations in the lipoma HMGIC fusion partner-like 5 (LHFPL5) gene cause autosomal recessive nonsyndromic hearing loss.Hum Mutat27633639

90. Liu X, Han D, Li J, Han B, Ouyang X (2010) Loss-of-function mutations in the PRPS1 gene cause a type of nonsyndromic X-linked sensorineural deafness, DFN2. Am J Hum Genet 86: 65–71.X. LiuD. HanJ. LiB. HanX. Ouyang2010Loss-of-function mutations in the PRPS1 gene cause a type of nonsyndromic X-linked sensorineural deafness, DFN2.Am J Hum Genet866571

91. Schultz JM, Yang Y, Caride AJ, Filoteo AG, Penheiter AR (2005) Modification of human hearing loss by plasma-membrane calcium pump PMCA2. N Engl J Med 352: 1557–1564.JM SchultzY. YangAJ CarideAG FiloteoAR Penheiter2005Modification of human hearing loss by plasma-membrane calcium pump PMCA2.N Engl J Med35215571564

92. Ficarella R, Di Leva F, Bortolozzi M, Ortolano S, Donaudy F (2007) A functional study of plasma-membrane calcium-pump isoform 2 mutants causing digenic deafness. Proc Natl Acad Sci U S A 104: 1516–1521.R. FicarellaF. Di LevaM. BortolozziS. OrtolanoF. Donaudy2007A functional study of plasma-membrane calcium-pump isoform 2 mutants causing digenic deafness.Proc Natl Acad Sci U S A10415161521

93. Li G, Vega R, Nelms K, Gekakis N, Goodnow C, et al. (2007) A role for Alström syndrome protein, alms1, in kidney ciliogenesis and cellular quiescence. PLoS Genet 3: e8.G. LiR. VegaK. NelmsN. GekakisC. Goodnow2007A role for Alström syndrome protein, alms1, in kidney ciliogenesis and cellular quiescence.PLoS Genet3e8

94. Karamatic Crew V, Burton N, Kagan A, Green CA, Levene C, et al. (2004) CD151, the first member of the tetraspanin (TM4) superfamily detected on erythrocytes, is essential for the correct assembly of human basement membranes in kidney and skin. Blood 104: 2217–2223.V. Karamatic CrewN. BurtonA. KaganCA GreenC. Levene2004CD151, the first member of the tetraspanin (TM4) superfamily detected on erythrocytes, is essential for the correct assembly of human basement membranes in kidney and skin.Blood10422172223

95. Stover EH, Borthwick KJ, Bavalia C, Eady N, Fritz DM, et al. (2002) Novel ATP6V1B1 and ATP6V0A4 mutations in autosomal recessive distal renal tubular acidosis with new evidence for hearing loss. J Med Genet 39: 796–803.EH StoverKJ BorthwickC. BavaliaN. EadyDM Fritz2002Novel ATP6V1B1 and ATP6V0A4 mutations in autosomal recessive distal renal tubular acidosis with new evidence for hearing loss.J Med Genet39796803

96. Wangemann P (2006) Supporting sensory transduction: cochlear fluid homeostasis and the endocochlear potential. J Physiol 576: 11–21.P. Wangemann2006Supporting sensory transduction: cochlear fluid homeostasis and the endocochlear potential.J Physiol5761121

97. Bulfone A, Menguzzato E, Broccoli V, Marchitiello A, Gattuso C, et al. (2000) Barhl1, a gene belonging to a new subfamily of mammalian homeobox genes, is expressed in migrating neurons of the CNS. Hum Mol Genet 9: 1443–1452.A. BulfoneE. MenguzzatoV. BroccoliA. MarchitielloC. Gattuso2000Barhl1, a gene belonging to a new subfamily of mammalian homeobox genes, is expressed in migrating neurons of the CNS.Hum Mol Genet914431452

98. NCBI UniGene: Organized view of the transcriptomeNCBI UniGene: Organized view of the transcriptomeAvailable: http://www.ncbi.nlm.nih.gov/unigene. Accessed April 2010. Available: http://www.ncbi.nlm.nih.gov/unigene. Accessed April 2010.

99. UniProtKB (Universal Protein resource KnowledgeBase)UniProtKB (Universal Protein resource KnowledgeBase)Available: http://www.uniprot.org/. Accessed April 2010. Available: http://www.uniprot.org/. Accessed April 2010.

100. HPRD (Human Protein Reference Database)HPRD (Human Protein Reference Database)Available: http://www.hprd.org/. Accessed April 2010. Available: http://www.hprd.org/. Accessed April 2010.

101. The Morton Lab: Cochlear EST databaseThe Morton Lab: Cochlear EST databaseAvailable: http://www.brighamandwomens.org/bwh_hearing/default.aspx. Accessed January 2009. Available: http://www.brighamandwomens.org/bwh_hearing/default.aspx. Accessed January 2009.

102. NCBI GEO: Gene Expression Omnibus DatabaseNCBI GEO: Gene Expression Omnibus DatabaseAvailable: http://www.ncbi.nlm.nih.gov/gds. Accessed April 2010. Available: http://www.ncbi.nlm.nih.gov/gds. Accessed April 2010.

103. IHR (Institute of Hearing Research)Table of gene expression in the developing ear. IHR (Institute of Hearing Research)Table of gene expression in the developing ear.Available: http://www.ihr.mrc.ac.uk/. Accessed April 2010. Available: http://www.ihr.mrc.ac.uk/. Accessed April 2010.

104. Bgee (A dataBase for Gene Expression Evolution)Bgee (A dataBase for Gene Expression Evolution)Available: http://bgee.unil.ch. Accessed April 2010. Available: http://bgee.unil.ch. Accessed April 2010.