ExInt: an Exon/Intron database

The Exon/Intron (ExInt) database incorporates information on the exon/intron structure of eukaryotic genes. Features in the database include: intron nucleotide sequence, amino acid sequence of the corresponding protein, position of the introns at the amino acid level and intron phase. From ExInt, we have also generated four additional databases each with ExInt entries containing predicted introns, introns experimentally defined, organelle introns or nuclear introns. ExInt is accessible through a retrieval system with pointers to GenBank. The database can be searched by keywords, locus name, NID, accession number or length of the protein. ExInt is freely accessible at http://intron.bic.nus.edu.sg/exint/exint.html     

Generation of a database containing discordant intron positions in eukaryotic genes (MIDB)

MOTIVATION: Intron sliding is the relocation of intron-exon boundaries over short distances and is often also referred to as intron slippage or intron migration or intron drift. We have generated a database containing discordant intron positions in homologous genes (MIDB--Mismatched Intron DataBase). Discordant intron positions are those that are either closely located in homologous genes (within a window of 10 nucleotides) or an intron position that is present in one gene but not in any of its homologs. The MIDB database aims at systematically collecting information about mismatched introns in the genes from GenBank and organizing it into a form useful for understanding the genomics and dynamics of introns thereby helping understand the evolution of genes. RESULTS: Intron displacement or sliding is critically important for explaining the present distribution of introns among orthologous and paralogous genes. MIDB allows examining of intron movements and allows mapping of intron positions from homologous proteins onto a single sequence. The database is of potential use for molecular biologists in general and for researchers who are interested in gene evolution and eukaryotic gene structure. Partial analysis of this database allowed us to identify a few putative cases of intron sliding. AVAILABILITY: http://intron.bic.nus.edu.sg/midb/midb.html     

Extremely short 20-33 nucleotide introns are the standard length in Paramecium tetraurelia.

Paramecium tetraurelia has the shortest known introns as its standard intron length. Sequenced introns vary between 20 and 33 nucleotides in length. The intron sequences were discovered in genomic sequences coding for a variety of different proteins, including phosphatases, kinases, and low-molecular weight GTP-binding proteins. All intron sequences begin with the conserved dinucleotide GT and end with the conserved dinucleotide AG. The sequences are more AT rich than the Paramecium coding sequences. The identified sequences were confirmed as introns by sequencing several cDNA fragments. We report here analysis of the characteristics of 50 separate introns, including size, base composition, and a consensus sequence.     

Statistical analysis of vertebrate sequences reveals that long genes are scarce in GC-rich isochores

We compared the exon/intron organization of vertebrate genes belonging to different isochore classes, as predicted by their GC content at third codon position. Two main features have emerged from the analysis of sequences published in GenBank: (1) genes coding for long proteins (i.e., 500 aa) are almost two times more frequent in GC-poor than in GC-rich isochores; (2) intervening sequences (=sum of introns) are on average three times longer in GC-poor than in GC-rich isochores. These patterns are observed among human, mouse, rat, cow, and even chicken genes and are therefore likely to be common to all warm-blooded vertebrates. Analysis of Xenopus sequences suggests that the same patterns exist in cold-blooded vertebrates. It could be argued that such results do not reflect the reality because sequence databases are not representative of entire genomes. However, analysis of biases in GenBank revealed that the observed discrepancies between GC-rich and GC-poor isochores are not artifactual, and are probably largely underestimated. We investigated the distribution of microsatellites and interspersed repeats in introns of human and mouse genes from different isochores. This analysis confirmed previous studies showing that Ll repeats are almost absent from GC-rich isochores. Microsatellites and SINES (Alu, B1, B2) are found at roughly equal frequencies in introns from all isochore classes. Globally, the presence of repeated sequences does not account for the increased intron length in GC-poor isochores. The relationships between gene structure and global genome organization and evolution are discussed.     

Convergent intron loss of MRP1 in Drosophila and mosquito species

Previous study revealed that the MRP1 gene ortholog DMRP1/CG6214 of Drosophila melanogaster contains 12 exons in the coding region. In the current study, the genes of DMRP1/CG6214 from D. melanogaster and Drosophila virilis were compared, and the result indicated that D. virilis had an extra intron located in exon 2, implying that intron loss or gain might have occurred at this locus. To track the evolution of the extra intron (Intron Z), orthologous nucleotide sequences of 37 arthropod species were cloned or annotated. Based on phylogenetic analysis, we found that Intron Z should present in the common ancestor of arthropod species, more than 420 Ma. In addition, we found that Sophophora subgenus species and mosquito (Culex pipiens) lost Intron Z independently, showing evolutionary convergence.     

Diversity in isochore structure among cold-blooded vertebrates based on GC content of coding and non-coding sequences

To review the general consideration about the different compositional structure of warm and cold-blooded vertebrates genomes, we used of the increasing number of genetic sequences, including coding (exons) and non-coding (introns) regions, that have been deposited on the databases throughout last years. The nucleotide distributions of the third codon positions (GC3) have been analyzed in 1510 coding sequences (CDS) of fish, 1414 CDS of amphibians and 320 CDS of reptiles. Also, the relationship between GC content of 74, 56 and 25 CDS of fish, amphibians and reptiles, respectively and that of their corresponding introns (GCI) have been considerated. In accordance with recent data, sequence analysis showed the presence of very GC3-rich CDS in these poikilotherm vertebrates. However, very high diversity in compositional patterns among different orders of fish, amphibians and reptiles was found. Significant positive correlations between GC3 and GCI was also confirmed for the genes analyzed. Nevertheless, introns resulted to be poorer in GC than their corresponding CDS, this difference being larger than in human genome. Because the limited number of available sequences including exons and introns we must be cautious about the results derived from them. However, the indicious of higher GC richness of coding sequences than of their corresponding introns could aid to understand the discrepancy of sequence analysis with the ultracentrifugation studies in cold-blooded vertebrates that did not predict the existence of GC-rich isochores.     

Nucleotide sequence of the gene for the b subunit of human factor XIII

Factor XIII (Mr 320,000) is a blood coagulation factor that stabilizes and strengthens the fibrin clot. It circulates in blood as a tetramer composed of two a subunits (Mr 75,000 each) and two b subunits (Mr 80,000 each). The b subunit consists of 641 amino acids and includes 10 tandem repeats of 60 amino acids known as GP-I structures, short consensus repeats (SCR), or sushi domains. In the present study, the human gene for the b subunit has been isolated from three different genomic libraries prepared in lambda phage. Fifteen independent phage with inserts coding for the entire gene were isolated and characterized by restriction mapping, Southern blotting, and DNA sequencing. The gene was found to be 28 kilobases in length and consisted of 12 exons (I-XII) separated by 11 intervening sequences. The leader sequence was encoded by exon I, while the carbonyl-terminal region of the protein was encoded by exon XII. Exons II-XI each coded for a single sushi domain, suggesting that the gene evolved through exon shuffling and duplication. The 12 exons in the gene ranged in size from 64 to 222 base pairs, while the introns ranged in size from 87 to 9970 nucleotides and made up 92% of the gene. The introns contained four Alu repetitive sequences, one each in introns A, E, I, and J. A fifth Alu repeat was present in the flanking 3' end of the gene. Two partial KpnI repeats were also found in the introns, including one in intron I and one in intron J. The KpnI repeat in intron J was 89% homologous to a sequence of approximately 2200 nucleotides flanking the gene coding for human beta globin and approximately 3800 nucleotides from the L1 insertion present in the gene for human factor VIII. Intron H also contained an "O" family repeat, while two potential regions for Z-DNA were identified within introns G and J. One nucleotide change was found in the coding region of the gene when its sequence was compared to that of the cDNA. This difference, however, did not result in a change in the amino acid sequence of the protein.     

Transterm: a database of messenger RNA components and signals

Transterm facilitates studies of messenger RNAs and translational control signals. Each messenger RNA (mRNA) from GenBank is extracted and broken into its functional components, its coding sequence, initiation context, termination context, flanking sequence representing its 5' UTR (untranslated region), 3' UTR and translational signals. In addition, numerical parameters characterising each coding region in Transterm, including codon and GC bias, are available. For each species in Transterm, the initiation and termination regions are aligned by their start or stop codons and presented as base frequency matrices and tables of the information content of the bases in the alignments. Users can obtain summaries of characteristics of the mRNAs for species of their choice and search for translational signals both in the Transterm database and in their own sequence. The current release contains data from over 10 000 species, including the complete genomes of 20 prokaryotes and three eukaryotes. Both flat-file and relational database forms of Transterm are accessible via the WWW at http://biochem.otago.ac.nz/Transterm/     

Exon-primed, intron-crossing (EPIC) loci for five nuclear genes in deep-sea protobranch bivalves: primer design, PCR protocols and locus utility

We describe PCR primers and amplification protocols developed to obtain introns from conserved nuclear genes in deep-sea protobranch bivalves. Because almost no sequence data for protobranchs are publically available, mollusk and other protostome sequences from GenBank were used to design degenerate primers, making these loci potentially useful in other invertebrate taxa. Amplification and sequencing success varied across the test group of 30 species, and we present five loci spanning this range of outcomes. Intron presence in the targeted regions also varied across genes and species, often within single genera; for instance, the calmodulin and β-tubulin loci contained introns with high frequency, whereas the triose phosphate isomerase locus never contained an intron. In introns for which we were able to obtain preliminary estimates of polymorphism levels in single species, polymorphism was greater than traditional mitochondrial loci. These markers will greatly increase the ability to assess population structure in the ecologically important protobranchs, and may prove useful in other taxa as well.     

The translational signal database, TransTerm: more organisms, complete genomes.

TransTerm is a database of initiation and termination sequence contexts from more than 250 organisms listed in GenBank, including the four complete genomes:Haemophilus influenzae, Methanococcus jannaschii, Mycoplasma genitalium,and Saccharomyces cerevisiae. For the current release, more than 60 000 coding sequences were analysed. The tabulated data include initiation and termination contexts organised by species along with quantitative parameters about individual coding sequences (length, %GC, GC3, Nc and CAI). There are also tables of initiation- and termination-region nucleotide-frequencies, codon usage tables and summaries of stop signal usage. TransTerm is available on the World Wide Web at: http://biochem.otago.ac.nz:800/Transterm/homepage.h tml     

Selective constraints on intron evolution in Drosophila

Intron sizes show an asymmetrical distribution in a number of organisms, with a large number of "short" introns clustered around a minimal intron length and a much broader distribution of longer introns. In Drosophila melanogaster, the short intron class is centered around 61 bp. The narrow length distribution suggests that natural selection may play a role in maintaining intron size. A comparison of 15 orthologous introns among species of the D. melanogaster subgroup indicates that, in general, short introns are not under greater DNA sequence or length constraints than long introns. There is a bias toward deletions in all introns (deletion/insertion ratio is 1.66), and the vast majority of indels are of short length (<10 bp). Indels occurring on the internal branches of the phylogenetic tree are significantly longer than those occurring on the terminal branches. These results are consistent with a compensatory model of intron length evolution in which slightly deleterious short deletions are frequently fixed within species by genetic drift, and relatively rare larger insertions that restore intron length are fixed by positive selection. A comparison of paralogous introns shared among duplicated genes suggests that length constraints differ between introns within the same gene. The janusA, janusB, and ocnus genes share two short introns derived from a common ancestor. The first of these introns shows significantly fewer indels than the second intron, although the two introns show a comparable number of substitutions. This indicates that intron-specific selective constraints have been maintained following gene duplication, which preceded the divergence of the D. melanogaster species subgroup.     

Sequences and Evolution of Human and Squirrel Monkey Blue Opsin Genes

The sequences of the entire blue opsin gene in the squirrel monkey (Saimiri boliviensis) and the five introns of the human blue opsin gene were obtained. Intron 3 of these genes contains an Alu sequence and intron 4 contains a partial mer13 sequence. A comparison of the squirrel monkey opsin sequence with published mammalian opsin sequences shows that features believed to be functionally critical are all conserved. However, the blue opsin has evolved twice as fast as rhodopsin and is only as conservative as the β globin, which has evolved at the average rate of mammalian proteins. Interestingly, the interhelical loops are, on average, actually more conservative than the transmembrane α helical regions. The introns of the blue opsin gene have evolved at the average rate of introns in primate genes. Received: 5 August 1996 / Accepted: 2 October 1996     

Remarkable sequence conservation of the last intron in the PKD1 gene

The last intron of the PKD1 gene (intron 45) was found to have exceptionally high sequence conservation across four mammalian species: human, mouse, rat, and dog. This conservation did not extend to the comparable intron in pufferfish. Pairwise comparisons for intron 45 showed 91% identity (human vs. dog) to 100% identity (mouse vs. rat) for an average for all four species of 94% identity. In contrast, introns 43 and 44 of the PKD1 gene had average pairwise identities of 57% and 54%, and exons 43, 44, and 45 and the coding region of exon 46 had average pairwise identities of 80%, 84%, 82%, and 80%. Intron 45 is 90 to 95 bp in length, with the major region of sequence divergence being in a central 4-bp to 9-bp variable region. RNA secondary structure analysis of intron 45 predicts a branching stem-loop structure in which the central variable region lies in one loop and the putative branch point sequence lies in another loop, suggesting that the intron adopts a specific stem-loop structure that may be important for its removal. Although intron 45 appears to conform to the class of small, G-triplet-containing introns that are spliced by a mechanism utilizing intron definition, its high sequence conservation may be a reflection of constraints imposed by a unique mechanism that coordinates splicing of this last PKD1 intron with polyadenylation.