BACKGROUND
hearing loss, acquired or genetic, is a major worldwide public health concern.
numerous genes have been linked to hearing disorders  <cit> .
these disorders may be syndromic or nonsyndromic; conductive, sensorineural, or mixed; and prelingual or postlingual  <cit> .
the various genetic forms of hearing loss are distinguished based on otologic, audiologic and physical examination combined with linkage analysis.
some representative deafness genes that have been identified include the alport syndrome , branchio-oto-renal syndrome , mohr-tranebjaerg syndrome , pendred syndrome , jervell and lange-nielsen syndrome , usher syndrome with its several types, norrie disease , dfnb <dig> , dfn <dig> , dfnb <dig> , dfna6/ <dig> , and several others  <cit> .
the mutational analysis of genes such as gjb <dig>  and gjb <dig>   <cit>  has aided diagnosis and geneticcounselling.
syndromic hearing loss is associated with a variety of other clinical findings and is relatively less prevalent.
in contrast, nonsyndromic hearing loss accounts for more than 70% of deafness cases and involves autosomal as well as x or y -linked deafness phenotypes  <cit> .
the molecular causes of nearly all nonsyndromic hearing loss are associated with inner ear structural damage, and changes in both the inner and the middle ear  <cit> .
mutations in genes such as the actg <dig>  coch, col11a <dig>  dfna <dig>  eya <dig>  gjb <dig>  gjb <dig>  kcnq <dig>  myo <dig>  myo7a, tecta, tmc <dig>  and wfs <dig>  as well as altered expression of genes such as gjb <dig> and myo1a have been associated with the autosomal dominant types that are generally progressive and involve changes in inner ear  <cit> .
the autosomal recessive phenotypes are associated with mutations in genes such as the cdh <dig>  cldn <dig>  espn, gjb <dig>  gjb <dig>  myo15a, myo <dig>  myo7a, otof, pcdh <dig>  slc26a <dig>  strc, tecta, tmc <dig>  tmie, tmprss <dig>  and ush1c, as well as altered expression of gjb <dig>  <cit> .
the map locations of a large number of nonsyndromic autosomal recessive deafness phenotypes are known, but the specific genes responsible for all these phenotypes have not been identified  <cit> .
the cloning of genes involved in such phenotypes requires refinement of the suspected genomic interval to as short a region as possible by linkage analysis.
however, it is not always possible to map a gene within an interval that is amenable for mutation analysis.
the mutation analysis of all genes encoded by a large genomic interval is extremely labor-intensive.
we describe here a bioinformatic approach that can reduce the candidate genes to a manageable number for mutation analysis.
initially, all the genes from a particular locus are cross-referenced to the databases of expressed mouse inner ear genes and the expressed human cochlear genes.
the alternative procedure included a search for interacting proteins for the gene products mapping to the candidate region.
as presented here, this approach has led to a set of specific candidate genes.
RESULTS
the locations of  <dig> autosomal dominant and  <dig> autosomal recessive nonsyndromic deafness phenotypes mapped to several chromosomes downloaded from hereditary hearing loss homepage are shown in tables  <dig> and  <dig> <cit> .
additional loci for nonsyndromic conditions are mapped to chromosomes  <dig>   <dig>  x and y  <cit> .
the hereditary hearing loss homepage is updated on a regular basis.
the marker boundaries of these locations encompass between  <dig>  and  <dig>  million basepairs  for various loci.
to generate a set of candidate genes for the listed loci, a strategy schematically represented in figure  <dig> was followed.
the determination of coding sequences and/or genes in a genomic region was made by unigene  <cit> .
however, the genes encoded in a large genomic interval are too many to be characterized by mutational analysis in a gene-by-gene approach.
therefore, we used the human cochlea and mouse inner ear expression databases  <cit>  to eliminate from the candidate list certain genes that were not expressed in these organs.
such in silico expression analysis relies on the assumption that the expression databases are comprehensive.
however, the characterization of all transcripts expressed in the ear is far from complete.
we, therefore, introduced another step in our candidate gene strategy by taking advantage of the human protein reference database  and generated a list of interacting genes for every gene mapping to candidate deafness loci  <cit> .
the rationale for protein interaction is as follows.
if a gene encoded in the candidate region interacts with a gene that is either involved in inner ear development/function, or a protein shows interaction with more than one candidate genes mapping to different loci, then such a gene is likely to be involved in the phenotype in question.
the interaction pattern of the gene products from usher syndrome is a good example to illustrate this point.
the known gene products for several usher syndrome loci are known to form interactions in vivo  <cit> .
the mutation of each one of these genes affects protein interactions and influences usher type  <dig> phenotype  <cit> .
the five forms of usher syndrome have defects in myosin viia, harmonin, cadherin  <dig>  protocadherin  <dig>  and a putative scaffolding protein sans.
harmonin binds sans, and it also binds myosin viia and protocadherin  <dig>  <cit> .
the role of cadherins in mediating cell-cell interaction is well-characterized.
furthermore, interactions of harmonin  with ush2a, ush2c and ush2b are mediated by pdz domains  <cit> .
in retrospect, if the interacting protein strategy had been used to select candidate genes for usher syndrome subtypes, it is likely that several genes could have been eliminated from consideration.
therefore, it is reasonable to assume that physical interactions will exist between proteins that are involved in inner ear developmental pathway or inner ear signal transduction pathways, and mutations in any one protein of the pathway is likely to give the same altered phenotype.
if proteins of interacting networks can be identified or predicted, then such genes are natural candidates for a given phenotype.
the above hypothesis is the underlying rationale for incorporating interacting proteins as a criterion for selecting candidate genes presented in this paper.
briefly, the strategy is as follows.
first, assemble the genes encoded in all candidate intervals, list the proteins that interact with genes in the candidate region, and then search for candidates on different loci that interact with a common protein.
such a criterion will fulfil the rationale of putative involvement of proteins at two different loci involved in a common biological process, and by association the respective genes mapping to two different loci will be considered as candidates.
the application of candidate gene isolation is demonstrated for the autosomal dominant condition dfna <dig>  the gene is mapped to the genomic interval 4q <dig> spanning  <dig> mbp  <cit> .
this region codes for  <dig> known and  <dig> hypothetical proteins   <cit> .
the comparison of these genes to expression databases reduced the list to  <dig> genes from the human cochlear database and three found in the mouse inner ear   <cit> .
the possibility remained for the elimination of a stronger candidate just on the basis that it did not score a hit in expression databases.
to avoid such an error, we have assembled lists of interacting proteins by using the human protein reference database   <cit>  for every gene identified by generetriever® from the candidate region.
if an interactor of a gene in the candidate interval is expressed in inner ear then the gene is considered a candidate.
alternatively, the interacting genes from a specific locus list were compared against lists from other loci to identify if a hit was scored against proteins among two or more lists.
the original genes corresponding to such interactor were considered as candidates for the respective deafness loci.
the strong candidates, based on the above analyses, for various deafness loci are presented in table  <dig>
in principle, the interactions-based strategy can be targeted to identify candidates for deafness if a database for interacting proteins involved in inner ear development and function is available.
for example, oncomodulin and prestin are expressed in outer hair cells  <cit> .
the protein interaction approach could link the possible candidate genes to specific cochlear cells by identifying known interactants.
if the interactors happen to map to a region harboring a deafness gene, such interactors are obvious candidates for mutational analysis.
however, such an approach will require identification of interacting proteins.
the primary limitation of the in silico approach described here is inadequate description of interacting protein networks.
the strong candidate list includes genes such as various cadherins, collagens, some cytoskeletal components and a number of growth factors and inner ear specific transcripts.
for example, hat  from the dfna <dig> locus is known to enhance cell growth and il- <dig> production.
it has been implicated in induction of par- <dig> -mediated il- <dig> release in psoriasis vulgaris  <cit> .
because hat is expressed in the ear, and protease activated receptor  has the ability to activate g-proteins followed by an increase in calcium ion concentration, we consider hat as a candidate.
kdr, a vascular endothelial growth factor receptor-type  <dig>  from the same locus shows age-dependent expression in the inner ear  <cit> .
our analyses indicated that only a fraction  of genes mapping to various genomic intervals was expressed in the inner ear.
we attribute these observations to depth of inner ear libraries.
it is likely that the genes being scored in these libraries have multiplicity for certain transcripts and absence of other transcripts.
for example, out of  <dig> genes at the dfna <dig> locus, only  <dig> genes are present in the cochlear library.
we cannot reasonably rule out the expression of the remaining  <dig> genes in the inner ear.
therefore, the approach presented here will be more comprehensive if we do not include ear expression in this scheme.
consequently, in a second attempt to mine the protein-interaction data obtained from the hrpd, we analyzed all genes encoded in the candidate intervals for their interactors.
the interaction data were considerably exhaustive and resulted in many more possible candidates with their expression not reported in the ear expression library.
a summary of gene numbers at different loci before and after interacting proteins analysis using the ear-expression scheme is presented in table  <dig>  the mouse syntenic genes are also indicated in these results.
the number of unfiltered candidate genes for each locus obtained by interacting proteins analysis is shown in table  <dig>  to elucidate the relevance of genes not found in the ear-expression library as possible candidates, we performed a literature search cross-referencing the identified gene with any reported hearing-associated condition in humans or other model animals.
some of these genes were linked to ear-development or hearing impairment as a secondary or unrelated symptom of other conditions.
for example, neurod <dig> gene mapping to dfnb <dig> locus was not reported in any of the inner ear libraries.
however, it appears to participate in the development of the auditory system as neurod <dig> null mice exhibited severe reduction of sensory neurons in the cochlear-vestibular ganglion  <cit> .
e2f <dig>  a transcription factor of the e2f family mapping to the dfna <dig> locus, may be indirectly implicated by its ability to regulate cell proliferation possibly during the developmental stages  <cit> .
other candidate genes from the unfiltered candidate analysis for the various loci are listed in table  <dig>  thus the unfiltered strategy adds  <dig> candidates for  <dig> loci and expands the candidate list to  <dig> genes for further mutation analysis.
our approach indicated the presence of possible candidates within most of the mapped loci.
however, prediction of candidate genes was not easy for loci indicated by asterisks in table  <dig>  because the genes mapping to these loci did not fulfil the criteria we have employed.
we further examined these genes on the basis of their reported function.
the following description pertains to specific genes that are not indicated in the candidate lists.
within the dfna <dig> locus, scna <dig> and scna <dig>  both being voltage-gated sodium channels, can be considered candidates based on involvement of related sodium channels in hearing  <cit> .
similarly, atp2c <dig> in dfna <dig> locus is a likely candidate because mutations in a related atpase have been shown in mice that are profoundly deaf and have a balance defect  <cit> .
the ephb <dig> gene, within the dfna <dig> locus, plays a major role during the development of the inner ear in mice  <cit> .
the dfna <dig> locus has six <dig> gene that plays a pivotal role in the control of the mouse otic vesicle patterning  <cit> .
neugrin, mapping to dfna <dig> locus, appeared to be an appropriate candidate as it was shown to be up-regulated throughout neuronal differentiation  <cit> .
a possible candidate for the dfna <dig> locus is the transcription factor nfib, an essential player in the maturation of lungs and brain development  <cit> .
the splicing regulation carried out by pnn, mapping to the dfnb <dig> region, is a reasonable candidate  <cit> .
we believe the genes presented in this article may serve as starting candidates toward identifying molecular mechanisms for specific deafness phenotypes.
CONCLUSIONS
we have used an in silico strategy to assemble a list of candidate genes that are positionally linked to and could be causing specific nonsyndromic hereditary hearing loss conditions.
as presented here, a list of  <dig> genes mapping to various genomic intervals have been narrowed down to  <dig> genes as candidates.
these candidates may be analyzed for mutations in various deafness phenotypes in parallel with attempts to further narrow down the suspected region by genetic linkage analysis.
it warrants mention that the potential of the approach presented here will be better harnessed as more information becomes available about inner ear transcripts and protein interaction networks.

 hearing loss, acquired or genetic, is a major worldwide public health concern. numerous genes have been linked to hearing disorders. these disorders may be syndromic or nonsyndromic; conductive, sensorineural, or mixed; and prelingual or postlingual. the various genetic forms of hearing loss are distinguished based on otologic, audiologic and physical examination combined with linkage analysis. some representative deafness genes that have been identified include the alport syndrome, branchio-oto-renal syndrome, mohr-tranebjaerg syndrome, pendred syndrome, jervell and lange-nielsen syndrome, usher syndrome with its several types, norrie disease, dfnb, dfn, dfnb, dfna6/, and several others. the mutational analysis of genes such as gjb and gjb has aided diagnosis and geneticcounselling. syndromic hearing loss is associated with a variety of other clinical findings and is relatively less prevalent. in contrast, nonsyndromic hearing loss accounts for more than 70% of deafness cases and involves autosomal as well as x or y -linked deafness phenotypes. the molecular causes of nearly all nonsyndromic hearing loss are associated with inner ear structural damage, and changes in both the inner and the middle ear. mutations in genes such as the actg coch, col11a dfna eya gjb gjb kcnq myo myo7a, tecta, tmc and wfs as well as altered expression of genes such as gjb and myo1a have been associated with the autosomal dominant types that are generally progressive and involve changes in inner ear. the autosomal recessive phenotypes are associated with mutations in genes such as the cdh cldn espn, gjb gjb myo15a, myo myo7a, otof, pcdh slc26a strc, tecta, tmc tmie, tmprss and ush1c, as well as altered expression of gjb. the map locations of a large number of nonsyndromic autosomal recessive deafness phenotypes are known, but the specific genes responsible for all these phenotypes have not been identified. the cloning of genes involved in such phenotypes requires refinement of the suspected genomic interval to as short a region as possible by linkage analysis. however, it is not always possible to map a gene within an interval that is amenable for mutation analysis. the mutation analysis of all genes encoded by a large genomic interval is extremely labor-intensive. we describe here a bioinformatic approach that can reduce the candidate genes to a manageable number for mutation analysis. initially, all the genes from a particular locus are cross-referenced to the databases of expressed mouse inner ear genes and the expressed human cochlear genes. the alternative procedure included a search for interacting proteins for the gene products mapping to the candidate region. as presented here, this approach has led to a set of specific candidate genes.  the locations of autosomal dominant and autosomal recessive nonsyndromic deafness phenotypes mapped to several chromosomes downloaded from hereditary hearing loss homepage are shown in tables and. additional loci for nonsyndromic conditions are mapped to chromosomes x and y. the hereditary hearing loss homepage is updated on a regular basis. the marker boundaries of these locations encompass between and million basepairs for various loci. to generate a set of candidate genes for the listed loci, a strategy schematically represented in figure was followed. the determination of coding sequences and/or genes in a genomic region was made by unigene. however, the genes encoded in a large genomic interval are too many to be characterized by mutational analysis in a gene-by-gene approach. therefore, we used the human cochlea and mouse inner ear expression databases to eliminate from the candidate list certain genes that were not expressed in these organs. such in silico expression analysis relies on the assumption that the expression databases are comprehensive. however, the characterization of all transcripts expressed in the ear is far from complete. we, therefore, introduced another step in our candidate gene strategy by taking advantage of the human protein reference database and generated a list of interacting genes for every gene mapping to candidate deafness loci. the rationale for protein interaction is as follows. if a gene encoded in the candidate region interacts with a gene that is either involved in inner ear development/function, or a protein shows interaction with more than one candidate genes mapping to different loci, then such a gene is likely to be involved in the phenotype in question. the interaction pattern of the gene products from usher syndrome is a good example to illustrate this point. the known gene products for several usher syndrome loci are known to form interactions in vivo. the mutation of each one of these genes affects protein interactions and influences usher type phenotype. the five forms of usher syndrome have defects in myosin viia, harmonin, cadherin protocadherin and a putative scaffolding protein sans. harmonin binds sans, and it also binds myosin viia and protocadherin. the role of cadherins in mediating cell-cell interaction is well-characterized. furthermore, interactions of harmonin with ush2a, ush2c and ush2b are mediated by pdz domains. in retrospect, if the interacting protein strategy had been used to select candidate genes for usher syndrome subtypes, it is likely that several genes could have been eliminated from consideration. therefore, it is reasonable to assume that physical interactions will exist between proteins that are involved in inner ear developmental pathway or inner ear signal transduction pathways, and mutations in any one protein of the pathway is likely to give the same altered phenotype. if proteins of interacting networks can be identified or predicted, then such genes are natural candidates for a given phenotype. the above hypothesis is the underlying rationale for incorporating interacting proteins as a criterion for selecting candidate genes presented in this paper. briefly, the strategy is as follows. first, assemble the genes encoded in all candidate intervals, list the proteins that interact with genes in the candidate region, and then search for candidates on different loci that interact with a common protein. such a criterion will fulfil the rationale of putative involvement of proteins at two different loci involved in a common biological process, and by association the respective genes mapping to two different loci will be considered as candidates. the application of candidate gene isolation is demonstrated for the autosomal dominant condition dfna the gene is mapped to the genomic interval 4q spanning mbp. this region codes for known and hypothetical proteins. the comparison of these genes to expression databases reduced the list to genes from the human cochlear database and three found in the mouse inner ear. the possibility remained for the elimination of a stronger candidate just on the basis that it did not score a hit in expression databases. to avoid such an error, we have assembled lists of interacting proteins by using the human protein reference database for every gene identified by generetriever® from the candidate region. if an interactor of a gene in the candidate interval is expressed in inner ear then the gene is considered a candidate. alternatively, the interacting genes from a specific locus list were compared against lists from other loci to identify if a hit was scored against proteins among two or more lists. the original genes corresponding to such interactor were considered as candidates for the respective deafness loci. the strong candidates, based on the above analyses, for various deafness loci are presented in table in principle, the interactions-based strategy can be targeted to identify candidates for deafness if a database for interacting proteins involved in inner ear development and function is available. for example, oncomodulin and prestin are expressed in outer hair cells. the protein interaction approach could link the possible candidate genes to specific cochlear cells by identifying known interactants. if the interactors happen to map to a region harboring a deafness gene, such interactors are obvious candidates for mutational analysis. however, such an approach will require identification of interacting proteins. the primary limitation of the in silico approach described here is inadequate description of interacting protein networks. the strong candidate list includes genes such as various cadherins, collagens, some cytoskeletal components and a number of growth factors and inner ear specific transcripts. for example, hat from the dfna locus is known to enhance cell growth and il- production. it has been implicated in induction of par- -mediated il- release in psoriasis vulgaris. because hat is expressed in the ear, and protease activated receptor has the ability to activate g-proteins followed by an increase in calcium ion concentration, we consider hat as a candidate. kdr, a vascular endothelial growth factor receptor-type from the same locus shows age-dependent expression in the inner ear. our analyses indicated that only a fraction of genes mapping to various genomic intervals was expressed in the inner ear. we attribute these observations to depth of inner ear libraries. it is likely that the genes being scored in these libraries have multiplicity for certain transcripts and absence of other transcripts. for example, out of genes at the dfna locus, only genes are present in the cochlear library. we cannot reasonably rule out the expression of the remaining genes in the inner ear. therefore, the approach presented here will be more comprehensive if we do not include ear expression in this scheme. consequently, in a second attempt to mine the protein-interaction data obtained from the hrpd, we analyzed all genes encoded in the candidate intervals for their interactors. the interaction data were considerably exhaustive and resulted in many more possible candidates with their expression not reported in the ear expression library. a summary of gene numbers at different loci before and after interacting proteins analysis using the ear-expression scheme is presented in table the mouse syntenic genes are also indicated in these results. the number of unfiltered candidate genes for each locus obtained by interacting proteins analysis is shown in table to elucidate the relevance of genes not found in the ear-expression library as possible candidates, we performed a literature search cross-referencing the identified gene with any reported hearing-associated condition in humans or other model animals. some of these genes were linked to ear-development or hearing impairment as a secondary or unrelated symptom of other conditions. for example, neurod gene mapping to dfnb locus was not reported in any of the inner ear libraries. however, it appears to participate in the development of the auditory system as neurod null mice exhibited severe reduction of sensory neurons in the cochlear-vestibular ganglion. e2f a transcription factor of the e2f family mapping to the dfna locus, may be indirectly implicated by its ability to regulate cell proliferation possibly during the developmental stages. other candidate genes from the unfiltered candidate analysis for the various loci are listed in table thus the unfiltered strategy adds candidates for loci and expands the candidate list to genes for further mutation analysis. our approach indicated the presence of possible candidates within most of the mapped loci. however, prediction of candidate genes was not easy for loci indicated by asterisks in table because the genes mapping to these loci did not fulfil the criteria we have employed. we further examined these genes on the basis of their reported function. the following description pertains to specific genes that are not indicated in the candidate lists. within the dfna locus, scna and scna both being voltage-gated sodium channels, can be considered candidates based on involvement of related sodium channels in hearing. similarly, atp2c in dfna locus is a likely candidate because mutations in a related atpase have been shown in mice that are profoundly deaf and have a balance defect. the ephb gene, within the dfna locus, plays a major role during the development of the inner ear in mice. the dfna locus has six gene that plays a pivotal role in the control of the mouse otic vesicle patterning. neugrin, mapping to dfna locus, appeared to be an appropriate candidate as it was shown to be up-regulated throughout neuronal differentiation. a possible candidate for the dfna locus is the transcription factor nfib, an essential player in the maturation of lungs and brain development. the splicing regulation carried out by pnn, mapping to the dfnb region, is a reasonable candidate. we believe the genes presented in this article may serve as starting candidates toward identifying molecular mechanisms for specific deafness phenotypes.  we have used an in silico strategy to assemble a list of candidate genes that are positionally linked to and could be causing specific nonsyndromic hereditary hearing loss conditions. as presented here, a list of genes mapping to various genomic intervals have been narrowed down to genes as candidates. these candidates may be analyzed for mutations in various deafness phenotypes in parallel with attempts to further narrow down the suspected region by genetic linkage analysis. it warrants mention that the potential of the approach presented here will be better harnessed as more information becomes available about inner ear transcripts and protein interaction networks. 
the human genome sequence and expression profiles of transcripts in the inner ear have aided positional cloning approaches. we have employed a bioinformatic approach to assemble the genes encoded by genomic regions that harbor various deafness loci. we have established a list of strong candidate genes encoded by the regions linked to various nonsyndromic hereditary hearing loss phenotypes by using a novel bioinformatic approach. the advantages and shortcomings of this bioinformatic approach are discussed. the candidates presented here provide a starting point for mutational analysis in well-characterized families along with genetic linkage to refine the loci. such analyses have narrowed a list of genes from suspected regions of the genome to a manageable number of about for further analysis. the knowledge of protein interactions offers additional advantages in selecting candidate genes within a mapped region.
 there are more than genes for autosomal dominant and autosomal recessive nonsyndromic hereditary deafness that are yet to be cloned. the human genome sequence and expression profiles of transcripts in the inner ear have aided positional cloning approaches. the knowledge of protein interactions offers additional advantages in selecting candidate genes within a mapped region.  we have employed a bioinformatic approach to assemble the genes encoded by genomic regions that harbor various deafness loci. the genes were then in silico analyzed for their candidacy by expression pattern and ability to interact with other proteins. such analyses have narrowed a list of genes from suspected regions of the genome to a manageable number of about for further analysis.  we have established a list of strong candidate genes encoded by the regions linked to various nonsyndromic hereditary hearing loss phenotypes by using a novel bioinformatic approach. the candidates presented here provide a starting point for mutational analysis in well-characterized families along with genetic linkage to refine the loci. the advantages and shortcomings of this bioinformatic approach are discussed. 