Skip Navigation

This Article
Right arrow Abstract Freely available
Right arrow Print PDF (178K) Freely available
Right arrow Supplementary Data
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (29)
Right arrowRequest Permissions
Right arrow Commercial Re-use Guidelines
for Open Access NAR Content
Google Scholar
Right arrow Articles by Hishiki, T.
Right arrow Articles by Okubo, K.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Hishiki, T.
Right arrow Articles by Okubo, K.
Social Bookmarking
Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

Nucleic Acids Research, 2000, Vol. 28, No. 1 136-138
© 2000 Oxford University Press

BodyMap: a human and mouse gene expression database

Teruyoshi Hishiki1,2, Shoko Kawamoto1, Shinichi Morishita2 and Kousaku Okubo1,*

1Institute for Molecular and Cellular Biology, Osaka University, 1-3 Yamadaoka, Suita, Osaka 565-0871, Japan and 2Department of Genome Knowledge Discovery System, Institute of Medical Science, University of Tokyo, 4-6-1 Shiroganedai, Minato, Tokyo 108-8639, Japan

Received September 3, 1999; Revised and Accepted October 6, 1999.


   ABSTRACT
TOP
 ABSTRACT
INTRODUCTION
 CONSTRUCTING THE DATABASE
 DATABASE USAGE
 FUTURE DEVELOPMENT
 SUPPLEMENTARY MATERIAL
 REFERENCES
 
BodyMap is a human and mouse gene expression database that has been maintained since 1993. It is based on site-directed 3'-ESTs collected from non-biased cDNA libraries constructed at Osaka University and contains >270 000 sequences from 60 human and 38 mouse tissues. The site-directed nature of the sequence tags allows unequivocal grouping of tags representing the same transcript and provides abundance information for each transcript in different parts of the body. Our collection of ESTs was compared periodically with other public databases for cross referencing. The histological resolution of source tissues and unique cloning strategy that minimized cloning bias enabled BodyMap to support three unique mRNA based experiments in silico. First, the recurrence information for clones in each library provides a rough estimate of the mRNA composition of each source tissue. Second, a user can search the entire data set with nucleotide sequences or keywords to assess expression patterns of particular genes. Third, and most important, BodyMap allows a user to select genes that have a desired expression pattern in humans and mice. BodyMap is accessible through the WWW at http://bodymap.ims.u-tokyo.ac.jp


   INTRODUCTION
TOP
 ABSTRACT
 INTRODUCTION
CONSTRUCTING THE DATABASE
 DATABASE USAGE
 FUTURE DEVELOPMENT
 SUPPLEMENTARY MATERIAL
 REFERENCES
 
The genome contains not only a series of blueprints for the building of individual proteins but also a coordinated program of protein synthesis and the means for controlling its execution. To deduce the functions of genes from their nucleotide sequences, decoding of these two independent pieces of genomic information is necessary because the former, the structural information, defines the action of each gene product and the latter, the regulatory information, determines the consequence of the action in the context of life. Knowledge of the codon rules, secondary structure prediction methods, and libraries of functional peptide motifs have made the extraction of the structural information from genomic sequences possible. In contrast, the paucity of information about cis-elements that define the execution of genetic programs have made it difficult to understand regulation of expression based on genomic sequence alone, especially in higher eukaryotes. Accordingly, the importance of systematic collection of gene expression data in parallel with the genome sequencing effort has been emphasized (1).

BodyMap is the first systematic effort to identify genes and collect gene expression information for the human and mouse genomes (2–4). During collection of expressed sequence tags (EST) for construction of the BodyMap database, non-structural information contained in the mRNA, transcript abundance and anatomical distribution were carefully preserved. The libraries were constructed from well characterized sources so as to minimize the differences in cloning efficiencies among transcripts, and libraries were never amplified prior to sequencing. Accordingly, BodyMap has characteristics distinct from those of other public EST data sets, which favor gene variety at the cost of quantitative expression information by normalizing libraries and using highly complex sources.


   CONSTRUCTING THE DATABASE
TOP
 ABSTRACT
 INTRODUCTION
 CONSTRUCTING THE DATABASE
DATABASE USAGE
 FUTURE DEVELOPMENT
 SUPPLEMENTARY MATERIAL
 REFERENCES
 
Primary data
The 3'-EST data accumulated in BodyMap were all collected from locally-produced 3'-directed libraries that contain only the most 3'-terminal MboI fragments of cDNAs. They all started with GATC and were read toward the polyA tail. As of July 1999, BodyMap contained 159 429 human and 123 468 mouse ESTs. The steps in construction of the database that are described below are illustrated in Figure 1.



View larger version (37K):
[in this window]
[in a new window]
  		Figure 1. Steps in construction of the database

 
Elimination/masking of uninformative sequences
Before storing the sequences in the database, uninformative sequences were identified and eliminated or separated to maintain the quality of the data. Those sequences with >5% Ns, not starting with GATC, or having more than one GATC were eliminated. We then eliminated those sequences having >90% similarity in a overlap longer than 50 bp or 70% of the EST length with libraries of vector and ribosomal sequences. Sequences for mitochondrial transcripts were separated from the analysis flow. The lengths of ESTs in BodyMap were determined primarily by the location of MboI sites in the cDNAs, and those sequences with a GATC site located within 20 bp of the polyA tail were separated out because they were considered too short for gene identification. Lastly, sequences were compared with a library of repetitive sequences based on REPBASE (5) using BLAST (6), and repetitive regions were masked. The repeat masking task was carried out as follows. Matching regions of no less than 20 bp overlap were identified. To simulate our judgement and the tendency of repeats to have many variations, we defined the threshold for similarity as a function of the length of the overlap. It started at 95% for 20 bp, and the threshold declined until it reached 70% for overlaps of >100 bp. Finally, the matched regions were masked.

The inspection of most sequences was carried out automatically, but ~1% of the sequences were reviewed manually. The inspected sequences are stored in the database.

Transcript counting/EST clustering
Newly-submitted tag sequences, ~1000 from each new library, were compared using FASTA (7). When the similarity exceeded 95% in a overlap longer than 50 bp or 70% of tag length and the overlap started at a GATC, they were considered the same tag and clustered (primary cluster). From each cluster one representative tag sequence was selected and compared with the representative sequences from previously generated clusters.

Using the same criteria, clusters of the same tags were grouped, and a new representative tag was selected from the new cluster (secondary cluster). A five figure cluster ID referred to as the Gene Signature (GS) number was assigned to each independent cluster.

Cross references/searching for equivalents in other data sets
Representative sequences for the GS clusters were compared periodically with the primate division of GenBank and ESTs used in the latest version of UniGene. The criteria for identity were the same as those for clustering. Based on the identities, links to one GenBank entry and one Hs_sequence in UniGene were made.

Expression matrix
As a result of the similarity searches, all GS clusters are indicated with their recurrence in all libraries in BodyMap, their identity in GenBank, and their corresponding UniGene sequence. All data are stored in a relational database (Sybase). A spreadsheet showing the occurrence of each tag in each library (GS versus sources), referred to as the expression matrix, was generated periodically. For efficiency, every query from the user interface concerning expression patterns is processed using the latest expression matrix.

GS card
Detailed information for each GS cluster is presented in the form of a ‘GS card’ that shows the representative sequence, presence or absence of repetitive sequences, the multiple alignment of the tag sequences represented by the GS, and the distribution across all analyzed libraries.


   DATABASE USAGE
TOP
 ABSTRACT
 INTRODUCTION
 CONSTRUCTING THE DATABASE
 DATABASE USAGE
FUTURE DEVELOPMENT
 SUPPLEMENTARY MATERIAL
 REFERENCES
 
BodyMap serves as computerized multi-tissue mRNA. Three types of sheets are prepared for each query. For every query, a portion of the compressed expression matrix is returned. From each compressed matrix, detailed information for each GS can be selected (Fig. 2).



View larger version (34K):
[in this window]
[in a new window]
  		Figure 2. The navigation map in the database

 
Composition of mRNA
In the query sheet, by selecting one source tissue from a menu, a GS list is created in descending order of recurrence in the selected library. This list shows the most abundant transcripts in the selected tissue and their rough concentrations in total mRNA. Such data are useful as a reference for microarray hybridizations or for promoter selection in gene transfer experiments.

Expression patterns of genes
GS can be searched by nucleotide sequences, keywords, GenBank accession numbers or UniGene IDs. Accession numbers and keywords are converted to nucleotide sequences using GenBank annotations prior to searching. Matches are returned in the form of an expression matrix from which a user can go into the GS cards for details.

Select genes by expression patterns
The query sheet allows users to define expression patterns by complex logical expressions regarding presence or absence of the tags in each library.

Although the total numbers of analyzed clones are ~1000 for each library, the well-defined nature and the relative homogeneity of cell types allows joining of multiple libraries into a hypothetical tissue or tissues that makes the selection more sensitive. For example, to select genes active preferentially in nervous tissues, the corresponding query would be ‘find those genes present in sources from nervous system and absent in the rest’. In the same manner, users can select the transcripts unique to some systems as well as those common across different systems. To facilitate writing such complex expressions, the web page provides a form for pattern definition that can be filled quickly by checks of radio buttons. On the form, each library is provided with a set of boxes representing one of three conditions: absent in any, present in some and present in all.


   FUTURE DEVELOPMENT
TOP
 ABSTRACT
 INTRODUCTION
 CONSTRUCTING THE DATABASE
 DATABASE USAGE
 FUTURE DEVELOPMENT
SUPPLEMENTARY MATERIAL
 REFERENCES
 
Analysis of regulated transcription based on statistical models (8) is currently in progress. In the coming update, prior to gene selection users will be able to set a confidence level for the inference about the differential expression between two libraries. For example, presence of a tag in brain will not affect the selection of the gastrin gene as preferentially expressed in gastric mucosa.

Incorporation of the emerging expression information into the coordinates of BodyMap expression matrix is also scheduled. PCR-based gene-by-gene analysis of expression pattern data that is being collected at Osaka University will be the first data to be included.


   SUPPLEMENTARY MATERIAL
TOP
 ABSTRACT
 INTRODUCTION
 CONSTRUCTING THE DATABASE
 DATABASE USAGE
 FUTURE DEVELOPMENT
 SUPPLEMENTARY MATERIAL
REFERENCES
 
Supplementary Material available at NAR Online consists of four figures representing (i) GS cluster; (ii) expression matrix; (iii) GS card and (iv) a summary of statistics.


   FOOTNOTES
 
* To whom correspondence should be addressed. Tel: +81 6 6879 7992; Fax: +81 6 6877 1922; Email: kousaku@imcb.osaka-u.ac.jp Back


   REFERENCES
TOP
 ABSTRACT
 INTRODUCTION
 CONSTRUCTING THE DATABASE
 DATABASE USAGE
 FUTURE DEVELOPMENT
 SUPPLEMENTARY MATERIAL
 REFERENCES
 

    1 Okubo,K. and Matsubara,K. (1997) FEBS Lett., 403, 225–229.[Web of Science][Medline]

    2 Bortoluzzi,S. and Danieli,G.A. (1999) Trends Genet., 15, 118–119.[Web of Science][Medline]

    3 Okubo,K., Hori,N., Matoba,R., Niiyama,T., Fukushima,A., Kojima,Y. and Matsubara,K. (1992) Nature Genet., 2, 173–179.[Web of Science][Medline]

    4 Kawamoto,S., Matsumoto,Y., Mizuno,K., Okubo,K. and Matsubara,K. (1996) Gene, 174, 151–158.[Web of Science][Medline]

    5 Jurka,J. (1998) Curr. Opin. Struct. Biol., 8, 333–337.[Web of Science][Medline]

    6 Altschul,S.F., Madden,T.L., Schaffer,A.A., Zhang,J., Zhang,Z., Miller,W. and Lipman,D.J. (1997) Nucleic Acids Res., 25, 3389–3402.[Abstract/Free Full Text]

    7 Pearson,W.R. and Lipman,D.J. (1988) Proc. Natl Acad. Sci. USA, 85, 2444–2448.[Abstract/Free Full Text]

    8 Audic,S. and Claverie,J.M. (1997) Genome Res., 7, 986–995.[Abstract/Free Full Text]


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?


This article has been cited by other articles:


Home page
 IOVSHome page
K. Yamasaki, S. Kawasaki, R. D. Young, H. Fukuoka, H. Tanioka, M. Nakatsukasa, A. J. Quantock, and S. Kinoshita
Genomic Aberrations and Cellular Heterogeneity in SV40-Immortalized Human Corneal Epithelial Cells
Invest. Ophthalmol. Vis. Sci., February 1, 2009; 50(2): 604 - 613.
[Abstract] [Full Text] [PDF]

Home page
 Nucleic Acids ResHome page
Y. Zhang, S.-M. Luoh, L. S. Hon, R. Baertsch, W. I. Wood, and Z. Zhang
GeneHub-GEPIS: digital expression profiling for normal and cancer tissues based on an integrated gene database
Nucleic Acids Res., July 13, 2007; 35(suppl_2): W152 - W158.
[Abstract] [Full Text] [PDF]

Home page
 Nucleic Acids ResHome page
O. Ogasawara, M. Otsuji, K. Watanabe, T. Iizuka, T. Tamura, T. Hishiki, S. Kawamoto, and K. Okubo
BodyMap-Xs: anatomical breakdown of 17 million animal ESTs for cross-species comparison of gene expression
Nucleic Acids Res., January 1, 2006; 34(suppl_1): D628 - D631.
[Abstract] [Full Text] [PDF]

Home page
 Nucleic Acids ResHome page
M. Tanino, M.-A. Debily, T. Tamura, T. Hishiki, O. Ogasawara, K. Murakawa, S. Kawamoto, K. Itoh, S. Watanabe, S. J. de Souza, et al.
The Human Anatomic Gene Expression Library (H-ANGEL), the H-Inv integrative display of human gene expression across disparate technologies and platforms
Nucleic Acids Res., January 1, 2005; 33(suppl_1): D567 - D572.
[Abstract] [Full Text] [PDF]

Home page
 Proc. Natl. Acad. Sci. USAHome page
H. R. Ueda, S. Hayashi, S. Matsuyama, T. Yomo, S. Hashimoto, S. A. Kay, J. B. Hogenesch, and M. Iino
Universality and flexibility in gene expression from bacteria to human
PNAS, March 16, 2004; 101(11): 3765 - 3769.
[Abstract] [Full Text] [PDF]

Home page
 Genome ResHome page
S. Nelander, P. Mostad, and P. Lindahl
Prediction of Cell Type-Specific Gene Modules: Identification and Initial Characterization of a Core Set of Smooth Muscle-Specific Genes
Genome Res., August 1, 2003; 13(8): 1838 - 1854.
[Abstract] [Full Text] [PDF]

Home page
 Physiol. GenomicsHome page
M. Ho, E. Yang, G. Matcuk, D. Deng, N. Sampas, A. Tsalenko, R. Tabibiazar, Y. Zhang, M. Chen, S. Talbi, et al.
Identification of endothelial cell genes by combined database mining and microarray analysis
Physiol Genomics, May 13, 2003; 13(3): 249 - 262.
[Abstract] [Full Text] [PDF]

Home page
 Nucleic Acids ResHome page
J. Sese, H. Nikaidou, S. Kawamoto, Y. Minesaki, S. Morishita, and K. Okubo
BodyMap incorporated PCR-based expression profiling data and a gene ranking system
Nucleic Acids Res., January 1, 2001; 29(1): 156 - 158.
[Abstract] [Full Text] [PDF]

Home page
 Genome ResHome page
S. Kawamoto, J. Yoshii, K. Mizuno, K. Ito, Y. Miyamoto, T. Ohnishi, R. Matoba, N. Hori, Y. Matsumoto, T. Okumura, et al.
BodyMap: A Collection of 3' ESTs for Analysis of Human Gene Expression Information
Genome Res., November 1, 2000; 10(11): 1817 - 1827.
[Abstract] [Full Text]

Home page
 IOVSHome page
S. Bortoluzzi, F. d’Alessi, and G. A. Danieli
A Novel Resource for the Study of Genes Expressed in the Adult Human Retina
Invest. Ophthalmol. Vis. Sci., October 1, 2000; 41(11): 3305 - 3308.
[Abstract] [Full Text]

This Article
Right arrow Abstract Freely available
Right arrow Print PDF (178K) Freely available
Right arrow Supplementary Data
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (29)
Right arrowRequest Permissions
Right arrow Commercial Re-use Guidelines
for Open Access NAR Content
Google Scholar
Right arrow Articles by Hishiki, T.
Right arrow Articles by Okubo, K.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Hishiki, T.
Right arrow Articles by Okubo, K.
Social Bookmarking
Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?