当前位置: 首页 > 期刊 > 《核酸研究》 > 2005年第We期 > 正文
编号:11369741
GeneSeeker: extraction and integration of human disease-related inform
http://www.100md.com 《核酸研究医学期刊》
     Centre for Molecular and Biomolecular Informatics, Radboud University Nijmegen PO Box 9010, 6500GL Nijmegen, The Netherlands 1Dalicon BV PO Box 354, 6700AJ Wageningen, The Netherlands 2Genomics Laboratory, University Medical Centre Utrecht PO Box 85060, 3508AB Utrecht, The Netherlands 3Wageningen University and Research Centre Wageningen, The Netherlands 4Department of Human Genetics, University Medical Centre Nijmegen PO Box 9101, 6500HB Nijmegen, The Netherlands

    *To whom correspondence should be addressed. Tel: +31 24 36 53391; Fax: +31 24 36 52977; Email: G.Vriend@cmbi.ru.al

    ABSTRACT

    The identification of genes underlying human genetic disorders requires the combination of data related to cytogenetic localization, phenotypes and expression patterns, to generate a list of candidate genes. In the field of human genetics, it is normal to perform this combination analysis by hand. We report on GeneSeeker (http://www.cmbi.ru.nl/GeneSeeker/), a web server that gathers and combines data from a series of databases. All database searches are performed via the web interfaces provided with the original databases, guaranteeing that the most recent data are queried, and obviating data warehousing. GeneSeeker makes the same selection of candidate genes as the human geneticists would have performed, and thus reducing the time-consuming process to a few minutes. GeneSeeker is particularly well suited for syndromes in which the disease gene displays altered expression patterns in the affected tissue(s).

    INTRODUCTION

    The identification of causative genes in human genetic disorders will be accelerated by the wealth of ‘omics’ information being generated. Geneticists consult a number databases to search for these genes. Each database concentrates on a different (molecular) aspect. In addition, databases have their own user interface, different formats to present the data and sometimes even their own ontologies. Data, such as gene localization and expression patterns, may be distributed over multiple databases.

    Geneticists normally collect phenotypic and/or expression data and the genes in the chromosomal region(s) of interest, and combine these to get a list of candidate genes. The rationale for this is that the gene that causes a disease is most probably expressed in the tissues affected by that disease (1–3). Using model organisms, such as the mouse, it is often possible to obtain information on genes, proteins, protein interactions and other functional attributes that can be transferred to Homo sapiens by means of synteny and protein homology relationships. The use of data from other species (such as mouse) often proves helpful in identifying the location or function of the equivalent human gene (4). GeneSeeker mimics this multi-species identification strategy (5).

    MATERIALS AND METHODS

    Databases used

    Table 1 lists the databases that GeneSeeker queries. These are divided over database groups (DB-groups). All databases are accessed through their standard WWW interfaces except MIMMAP and OXFORD. MIMMAP is a reformatted version of the OMIM (6) gene mapping information. OXFORD is used to translate human to mouse chromosomal locations, and is described in more detail in the pre-processing section. We use SRS (Sequence Retrieval System, Lion Biosciences, Cambridge, UK) to access these two databases (7). The SRS parser was modified to allow searches for chromosomal ranges.

    Table 1 Databases accessed by the GeneSeeker

    Data processing

    The layout of the GeneSeeker web server is shown in Figure 1. The user query consists of a chromosomal band range using standard nomenclature (e.g. 7p15–p21). This cytogenetic localization is passed through DB-group 1. Syntenic regions in the mouse are sought in DB-group 2 using an Oxford-grid. Tissues of interest or phenotypic features of a syndrome can be specified by the user as a Boolean expression that is split up and processed by DB-group 3. This modular set-up makes it easy to add extra DB-groups in the future. For every database, a plug-in was designed to perform all tasks from user-query pre-processing to query-result post-processing. These plug-ins deal with a series of technical topics, such as query reformatting, generating the correct URL, filling in the form on that database's web interface, requesting all hits rather than in chunks, parsing the database HTML output and so on.

    Figure 1 Overview of GeneSeeker. The query, which consists of a cytogenetic localization, a phenotypic description and expression data, is divided over the three DB-groups that use the database-specific plug-ins to deal with all topics ranging from user-query pre-processing to post-processing of the query output. Results from each DB-group are merged with a Boolean OR. The results of the three DB-groups are combined as specified in the user query.

    The name of a gene can vary from database to database. The gene for the multi-drug resistance-associated protein 1, for example, is stored as ABCC1, MRP or MRP1, depending on the database used. These gene nomenclature problems have to be solved because GeneSeeker depends on the gene names in the combination steps. For each DB-group the results are integrated with a Boolean OR. The resulting gene lists of the three DB-groups are combined according to the Boolean logic specified in the user query.

    Implementation issues

    Parallelization

    The database plug-ins run in parallel to minimize the waiting time. A queuing system prevents excessive loads on remote servers. The plug-ins return the results of the queries to GeneSeeker as a list containing the gene names and corresponding database hyperlinks.

    Mouse–human synteny

    An Oxford grid (8) is used to find the homologous genes and gene regions in the mouse genome for all human chromosome locations entered by the user. A human chromosomal band range is translated into the corresponding mouse chromosome locations. Two mouse locations are combined if the genetic distance is shorter than a user-specified value (defaults to 10 cM). We regenerate this Oxford grid weekly to ensure that the latest synteny information is used in each query.

    Gene nomenclature

    Inconsistent gene nomenclature is resolved using gene synonym information from UniProt database (9). We use the MGD human homologues information to interconvert mouse and human gene names. We maintain local copies of these conversion tables because nearly all queries require that gene nomenclature problems be solved.

    User interface

    The GeneSeeker interface consists of the query form shown in Figure 2 and an options form that usually requires no user input. A genetic localization and the phenotypic/expression terms should be entered for a meaningful search. Databases that generate more noise than signal can be removed from the query. The user can also suppress the display of housekeeping genes or a specified list of genes. The options form contains a thesaurus (10) that can help the user to select the correct expression terms: for example, when the user is interested in a genetic trait that results in abnormalities in the brain, selection of the ‘brain’ category returns the hints ‘brain or hindbrain or forebrain...’. Hints for the genetic localization data can be found in a table containing frequently aberrant chromosomal bands in specific disorders taken from literature (11,12). The user can be notified on request about the completion of GeneSeeker searches by email. All parameters are linked to help screens. The results are presented in four tables (Figure 3).

    Figure 2 An example of a GeneSeeker query. Analyses of Trismus-Pseudocamptodactyly syndrome (TPC; MIM 158 300) has been linked to 17p12–p13.1 (13). TPC is characterized by defects in muscle tissue mainly in limb and/or mouth. The options form is data not shown.

    Figure 3 The output of GeneSeeker for the Trismus-Pseudocamptodactyly syndrome query (see Figure 2). It has been shown that mutations in the MYH8 gene can cause TPC (13). Top left table: genes that agree perfectly with the user query. Top right table: genes found in mouse syntenic regions that cannot be mapped automatically on the human genome, but match the expression pattern. Bottom left table: genes found in mouse syntenic regions that match the expression pattern, but map on the human genome outside the candidate cytogenetic region. Bottom right table: human genes in the candidate cytogenetic region that do not match the phenotype/expression pattern. All genes are hyperlinked to the underlying database, and, when possible, to GeneCards (14).

    RESULTS AND DISCUSSION

    The GeneSeeker offers a user-friendly quick scan of several databases that are commonly used by geneticists to identify candidate genes for specific Mendelian diseases. As such, GeneSeeker uses those databases that are most appropriate for the questions asked. Several aspects are likely to change in the near future as genomics and genetics develop. For example, our usage of an Oxford grid can be improved or replaced as soon as consensus is reached about the localization of genes on the mouse and human genomes among the various databases. Expression pattern information (e.g. microarray data) is growing rapidly, and is expected to become useful for GeneSeeker in the near future. At the moment, publicly available expression information is still sparse, scattered and not yet standardized.

    In its present form, GeneSeeker is best suited for syndromes in which one can assume aberrant or absent gene expression in the affected tissues. GeneSeeker allows the user to query heterogeneous databases and obtain good candidate genes for the disease of interest based on positional, expression and model data (5). With the present hardware set-up GeneSeeker can perform 1000 searches per day.

    ACKNOWLEDGEMENTS

    We thank David Thomas for helpful corrections to the manuscript. This work was supported by NWO/Unilever, the Irene Kinderziekenhuis Foundation and the EU FP6 Programme (LHSG-CT-2003-503265). Funding to pay the Open Access publication charges for this article was provided by Radboud University Nijmegen.

    REFERENCES

    Blackshaw, S., Fraioli, R.E., Furukawa, T., Cepko, C.L. (2001) Comprehensive analysis of photoreceptor gene expression and the identification of candidate retinal disease genes Cell, 107, 579–589 .

    den Hollander, A.I., van Driel, M.A., de Kok, Y.J., van de Pol, D.J., Hoyng, C.B., Brunner, H.G., Deutman, A.F., Cremers, F.P. (1999) Isolation and mapping of novel candidate genes for retinal disorders using suppression subtractive hybridization Genomics, 58, 240–249 .

    Dryja, T.P. (1997) Gene-based approach to human gene–phenotype correlations Proc. Natl Acad. Sci. USA, 94, 12117–12121 .

    Chiang, A.P., Nishimura, D., Searby, C., Elbedour, K., Carmi, R., Ferguson, A.L., Secrist, J., Braun, T., Casavant, T., Stone, E.M., et al. (2004) Comparative genomic analysis identifies an ADP-ribosylation factor-like gene as the cause of Bardet–Biedl syndrome (BBS3) Am. J. Hum. Genet., 75, 475–484 .

    van Driel, M.A., Cuelenaere, K., Kemmeren, P.P., Leunissen, J.A., Brunner, H.G. (2003) A new web-based data mining tool for the identification of candidate genes for human genetic disorders Eur. J. Hum. Genet., 11, 57–63 .

    Hamosh, A., Scott, A.F., Amberger, J., Bocchini, C., Valle, D., McKusick, V.A. (2002) Online Mendelian Inheritance in Man (OMIM), a knowledgebase of human genes and genetic disorders Nucleic Acids Res., 30, 52–55 .

    Etzold, T., Ulyanov, A., Argos, P. (1996) SRS: information retrieval system for molecular biology data banks Methods Enzymol., 266, 114–128 .

    Edwards, J.H. (1991) The Oxford Grid Ann. Hum. Genet., 55, 17–31 .

    Apweiler, R., Bairoch, A., Wu, C.H., Barker, W.C., Boeckmann, B., Ferro, S., Gasteiger, E., Huang, H., Lopez, R., Magrane, M., et al. (2004) UniProt: the Universal Protein knowledgebase Nucleic Acids Res., 32, D115–D119 .

    van Steensel, M.A., Celli, J., van Bokhoven, J.H., Brunner, H.G. (1999) Probing the gene expression database for candidate genes Eur. J. Hum. Genet., 7, 910–919 .

    Brewer, C., Holloway, S., Zawalnyski, P., Schinzel, A., FitzPatrick, D. (1998) A chromosomal deletion map of human malformations Am. J. Hum. Genet., 63, 1153–1159 .

    Brewer, C., Holloway, S., Zawalnyski, P., Schinzel, A., FitzPatrick, D. (1999) A chromosomal duplication map of malformations: regions of suspected haplo- and triplolethality—and tolerance of segmental aneuploidy—in humans Am. J. Hum. Genet., 64, 1702–1708 .

    Veugelers, M., Bressan, M., McDermott, D.A., Weremowicz, S., Morton, C.C., Mabry, C.C., Lefaivre, J.F., Zunamon, A., Destree, A., Chaudron, J.M., et al. (2004) Mutation of perinatal myosin heavy chain associated with a Carney complex variant N. Engl. J. Med., 351, 460–469 .

    Safran, M., Chalifa-Caspi, V., Shmueli, O., Olender, T., Lapidot, M., Rosen, N., Shmoish, M., Peter, Y., Glusman, G., Feldmesser, E., et al. (2003) Human Gene-Centric Databases at the Weizmann Institute of Science: GeneCards, UDB, CroW 21 and HORDE Nucleic Acids Res., 31, 142–146 .

    Blake, J.A., Richardson, J.E., Bult, C.J., Kadin, J.A., Eppig, J.T. (2003) MGD: the Mouse Genome Database Nucleic Acids Res., 31, 193–195 .

    Letovsky, S.I., Cottingham, R.W., Porter, C.J., Li, P.W. (1998) GDB: the Human Genome Database Nucleic Acids Res., 26, 94–99 .

    Ringwald, M., Eppig, J.T., Begley, D.A., Corradi, J.P., McCright, I.J., Hayamizu, T.F., Hill, D.P., Kadin, J.A., Richardson, J.E. (2001) The Mouse Gene Expression Database (GXD) Nucleic Acids Res., 29, 98–101 .

    Woychik, R.P., Wassom, J.S., Kingsbury, D., Jacobson, D.A. (1993) TBASE: a computerized database for transgenic animals and targeted mutations Nature, 363, 375–376 .(M. A. van Driel, K. Cuelenaere1, P. P. C)