Nucleic Acids Research, 2002, Vol. 30, No. 1 191-194
© 2002 Oxford University Press
ExInt: an Exon Intron Database
Ludwig Institute for Cancer Research, Sao Paulo Branch, Rua Prof. Antonio Prudente 109, 01509-010, Sao Paulo, Brazil, 1Bioinformatics Center, National University of Singapore, Singapore and 2Department of Ecology and Evolution, University of Chicago, IL, USA
Received September 19, 2001; Accepted September 26, 2001.
| ABSTRACT |
|---|
|
|
|---|
The Exon/Intron Database (ExInt) stores information of all GenBank eukaryotic entries containing an annotated intron sequence. Data are available through a retrieval system, as flat-files and as a MySQL dump file. In this report we discuss several implementations added to ExInt, which is accessible at http://intron.bic.nus.edu.sg/exint/newexint/exint.html.
| INTRODUCTION |
|---|
|
|
|---|
The exponential growth of sequence databases, especially due to genome and EST sequencing, has generated a parallel increase in the amount of sequences showing an intron/exon organization. We have recently developed a database containing all sequences in GenBank bearing in their annotation at least one exon/intron boundary (1). This, and other related databases (2,3), has been used in several studies approaching issues related to the exon/intron organization of eukaryotic genes (4,5).
In this report, we describe a series of implementations to the Exon/Intron Database (ExInt) as follows:
1. Relational database: data are now stored in a relational database (MySQL). The table structure is presented in Figure 1. Data from the database tables can be downloaded in a dump format, which allows direct incorporation in other MySQL relational databases.
|
2. Purged database: it is known that GenBank is extremely redundant. To avoid any potential bias, we have made available in this latest version of ExInt a non-redundant set of the data. Overall analysis of both redundant and non-redundant sets confirmed that most of the sequences (>80%) are redundant in current databases. Both datasets are available for download as Fasta libraries. They are also searchable using ExInt Blast engine.
3. Statistics link: several statistical features (for the whole database and models species) are available, such as number of genes, exons and introns before and after purging (Table 1); exon length distribution (Fig. 2); intron length distribution (Fig. 3) and intron phase distribution (Table 2).
|
|
|
|
4. Validation of predicted gene structure using EST data: we provided a validated subset for genes predicted in seven species: Homo sapiens, Mus musculus, Rattus sp., Caenorhabditis elegans, Drosophila melanogaster, Arabidopsis thaliana and Saccharomyces cerevisae (Table 3).
|
| METHODOLOGY |
|---|
|
|
|---|
We have used GenBank release 122 to construct a raw database containing all eukaryotic sequences with an exon/intron organization. The approach used to identify all intron-containing sequences in GenBank has been described previously (1). The same is true for the methodology used to construct the following derived databases: predicted introns, experimentally defined introns, organelle and nuclear genes (1). A purged database was constructed using a modification of the method of Long et al. (6), as follows. We performed an all-against-all protein sequence comparison using a PVM-version of Fasta in an eight-node cluster of PCs running Linux. When two protein sequences have an identity level
25% over at least 70% of the length of the shorter sequence, just one sequence is kept. These comparisons are exhaustive until a complete non-redundant database is obtained. As a representative of the gene cluster we have taken the sequence with the largest number of exons and introns. To validate the predicted gene structures, we take the predicted cDNA structure (keeping the positional information of all predicted introns) for all genes within seven model species and used Blast (7) to search them against the respective (same species) EST datasets. A script in PERL was written to parse the Blast output looking for cases where a predicted exon/exon boundary (by that we mean a region in the cDNA where a predicted intron is present at the genomic level) was confirmed by at least one EST.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
ExInt contains a wealth of relevant biological information. Here, we present some statistics that are important to the database construction and for a general evaluation of the data. Table 1 shows the number of genes, exons and introns for the redundant and non-redundant datasets and for seven model species. We note that there are, on average, 5.48 exons per gene with AL445795 having the higher number (96). Figures 2 and 3 show the exon and intron length distributions, respectively. We confirm an observation from Deutsch and Long (8) that invertebrate introns are on average smaller than human introns. As also seen by Deutsch and Long (8), we have observed a bimodal distribution of intron length for the whole dataset, which does not seem to be due to predicted introns, since the same pattern is also observed for the confirmed introns (Fig. 3). Positioning of introns along the coding region (Fig. 4) shows a bias distribution towards the C-terminal half of the protein molecule. This piece of information is important for interpretation of data related to gene structure. For example, it has recently been suggested that alternative splicing events are more frequent on the C-terminal half of proteins (9), a bias that can be due to the distribution shown in Figure 4.
|
The validation of predicted gene structures is probably the most important implementation to ExInt. It has been shown that gene prediction programs may generate a large amount of artefactual gene structures, and analysis using these datasets may draw incorrect conclusions (10). We have made use of the large amount of EST data available in dbEST to validate the predicted gene structure for sequences of seven different model species, H.sapiens, M.musculus, Rattus sp., C.elegans, D.melanogaster, A.thaliana and S.cerevisae. This validation step creates a sub-set of trusted predicted gene structure that may be important in a number of biological queries. The sub-set of validated intron/exon boundaries may also constitute a useful resource for developers of gene prediction programs. It is important to emphasize that the absence of validation does not imply that the predicted gene structure is wrong, since the coverage of the transcriptome by ESTs is not yet complete.
| AVAILABILITY |
|---|
|
|
|---|
ExInt is accessible via a World Wide Web interface at http://intron.bic.nus.edu.sg/exint/newexint/exint.html. Different features can be used as a query element such as: NID, locus name and keyword. The whole database, as well as derived databases, is available for download. Derived databases include: purged database, predicted intron, experimentally defined introns, organelle genes and nuclear genes. Users can also search all databases with a query sequence using Blast. ExInt will be updated twice a year.
| ACKNOWEDGEMENT |
|---|
|
|
|---|
F.P. is supported by Fapesp (00/02228-9).
|
| FOOTNOTES |
|---|
* To whom correspondence should be addressed. Tel: +55 11 32074922; Fax: +55 11 32077001; Email: sandro{at}compbio.ludwig.org.br The authors wish it to be known that, in their opinion, the first two authors should be regarded as joint First Authors
| REFERENCES |
|---|
|
|
|---|
-
1 Sakharkar,M., Long,M., Tan,T.W. and de Souza,S.J. (2000) ExInt: an Exon/Intron database. Nucleic Acids Res., 23, 191192.
2 Saxonov,S., Daizadeh,I., Fedorov,A. and Gilbert,W. (2000) EID: the ExonIntron Databasean exhaustive database of protein-coding intron-containing genes. Nucleic Acids Res., 28, 185190.
3 Schisler,N.J. and Palmer,J.D. (2000) The IDB and IEDB: intron sequence and evolution databases. Nucleic Acids Res., 28, 181184.
4 Sakharkar,M., Kangueane,P., Woon,T.W., Tan,T.W., Long,M., Kolatkar,P.R. and de Souza,S.J. (2000) IEKban exon intron knowledge base from databases. Bioinformatics, 16, 11511152.
5 Sakharkar,M., Tan,T.W. and de Souza,S.J. (2001). Generation of a database containing discordant intron positions in eukaryotic genes (MIDB). Bioinformatics, 17, 671675.
6 Long,M., Rosenberg,C. and Gilbert,W. (1995) Intron phase correlations and the evolution of the intron/exon structure of genes. Proc. Natl Acad. Sci. USA, 92, 1249512499.
7 Altschul,S.F., Madden,T.L., Schaffer,A.A., Zhang,J., Zhang,Z., Miller,W. and Lipman,D.J. (1997) Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res., 25, 33893402.
8 Deutsch,M and Long,M. (1999) Intronexon structures of eukaryotic model organisms. Nucleic Acids Res., 27, 32193228.
9 Modrek,B., Resch,A., Grasso,C. and Lee,C. (2001) Genome-wide detection of alternative splicing in expressed sequences of human genes. Nucleic Acids Res., 29, 28502859.
10 Rogic,S., Mackworth,A.K. and Ouellette,F.B. (2001) Evaluation of gene-finding programs on mammalian sequences. Genome Res., 409, 685690.
This article has been cited by other articles:
![]() |
V. Shepelev and A. Fedorov Advances in the Exon-Intron Database (EID) Brief Bioinform, June 1, 2006; 7(2): 178 - 185. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. Issac and G. P. S. Raghava EGPred: Prediction of Eukaryotic Genes Using Ab Initio Methods After Combining With Sequence Similarity Approaches Genome Res., September 1, 2004; 14(9): 1756 - 1766. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. Gopalan, T. W. Tan, B. T. K. Lee, and S. Ranganathan Xpro: database of eukaryotic protein-encoding genes Nucleic Acids Res., January 1, 2004; 32(90001): D59 - 63. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Lynch and A. Kewalramani Messenger RNA Surveillance and the Evolutionary Proliferation of Introns Mol. Biol. Evol., April 1, 2003; 20(4): 563 - 571. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Kaessmann, S. Zollner, A. Nekrutenko, and W.-H. Li Signatures of Domain Shuffling in the Human Genome Genome Res., November 1, 2002; 12(11): 1642 - 1650. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||







