Bioinformatics analysis of plant orthologous introns: identification of an intronic tRNA-like sequence

Orthologous introns have identical positions relative to the coding sequence in orthologous genes of different species. By analyzing the complete genomes of five plants we generated a database of 40,512 orthologous intron groups of dicotyledonous plants, 28,519 orthologous intron groups of angiosperms, and 15,726 of land plants (moss and angiosperms). Multiple sequence alignments of each orthologous intron group were obtained using the Mafft algorithm. The number of conserved regions in plant introns appeared to be hundreds of times fewer than that in mammals or vertebrates. Approximately three quarters of conserved intronic regions among angiosperms and dicots, in particular, correspond to alternatively-spliced exonic sequences. We registered only a handful of conserved intronic ncRNAs of flowering plants. However, the most evolutionarily conserved intronic region, which is ubiquitous for all plants examined in this study, including moss, possessed multiple structural features of tRNAs, which caused us to classify it as a putative tRNA-like ncRNA. Intronic sequences encoding tRNA-like structures are not unique to plants. Bioinformatics examination of the presence of tRNA inside introns revealed an unusually long-term association of four glycine tRNAs inside the Vac14 gene of fish, amniotes, and mammals.     

Analysis of evolution of exon-intron structure of eukaryotic genes

The availability of multiple, complete eukaryotic genome sequences allows one to address many fundamental evolutionary questions on genome scale. One such important, long-standing problem is evolution of exon-intron structure of eukaryotic genes. Analysis of orthologous genes from completely sequenced genomes revealed numerous shared intron positions in orthologous genes from animals and plants and even between animals, plants and protists. The data on shared and lineage-specific intron positions were used as the starting point for evolutionary reconstruction with parsimony and maximum-likelihood approaches. Parsimony methods produce reconstructions with intron-rich ancestors but also infer lineage-specific, in many cases, high levels of intron loss and gain. Different probabilistic models gave opposite results, apparently depending on model parameters and assumptions, from domination of intron loss, with extremely intron-rich ancestors, to dramatic excess of gains, to the point of denying any true conservation of intron positions among deep eukaryotic lineages. Development of models with adequate, realistic parameters and assumptions seems to be crucial for obtaining more definitive estimates of intron gain and loss in different eukaryotic lineages. Many shared intron positions were detected in ancestral eukaryotic paralogues which evolved by duplication prior to the divergence of extant eukaryotic lineages. These findings indicate that numerous introns were present in eukaryotic genes already at the earliest stages of evolution of eukaryotes and are compatible with the hypothesis that the original, catastrophic intron invasion accompanied the emergence of the eukaryotic cells. Comparison of various features of old and younger introns starts shedding light on probable mechanisms of intron insertion, indicating that propagation of old introns is unlikely to be a major mechanism for origin of new ones. The existence and structure of ancestral protosplice sites were addressed by examining the context of introns inserted within codons that encode amino acids conserved in all eukaryotes and, accordingly, are not subject to selection for splicing efficiency. It was shown that introns indeed predominantly insert into or are fixed in specific protosplice sites which have the consensus sequence (A/C)AG|Gt.     

Origin and evolution of spliceosomal introns

ABSTRACT: Evolution of exon-intron structure of eukaryotic genes has been a matter of long-standing, intensive debate. The introns-early concept, later rebranded 'introns first' held that protein-coding genes were interrupted by numerous introns even at the earliest stages of life's evolution and that introns played a major role in the origin of proteins by facilitating recombination of sequences coding for small protein/peptide modules. The introns-late concept held that introns emerged only in eukaryotes and new introns have been accumulating continuously throughout eukaryotic evolution. Analysis of orthologous genes from completely sequenced eukaryotic genomes revealed numerous shared intron positions in orthologous genes from animals and plants and even between animals, plants and protists, suggesting that many ancestral introns have persisted since the last eukaryotic common ancestor (LECA). Reconstructions of intron gain and loss using the growing collection of genomes of diverse eukaryotes and increasingly advanced probabilistic models convincingly show that the LECA and the ancestors of each eukaryotic supergroup had intron-rich genes, with intron densities comparable to those in the most intron-rich modern genomes such as those of vertebrates. The subsequent evolution in most lineages of eukaryotes involved primarily loss of introns, with only a few episodes of substantial intron gain that might have accompanied major evolutionary innovations such as the origin of metazoa. The original invasion of self-splicing Group II introns, presumably originating from the mitochondrial endosymbiont, into the genome of the emerging eukaryote might have been a key factor of eukaryogenesis that in particular triggered the origin of endomembranes and the nucleus. Conversely, splicing errors gave rise to alternative splicing, a major contribution to the biological complexity of multicellular eukaryotes. There is no indication that any prokaryote has ever possessed a spliceosome or introns in protein-coding genes, other than relatively rare mobile self-splicing introns. Thus, the introns-first scenario is not supported by any evidence but exon-intron structure of protein-coding genes appears to have evolved concomitantly with the eukaryotic cell, and introns were a major factor of evolution throughout the history of eukaryotes. This article was reviewed by I. King Jordan, Manuel Irimia (nominated by Anthony Poole), Tobias Mourier (nominated by Anthony Poole), and Fyodor Kondrashov. For the complete reports, see the Reviewers' Reports section.     

Molecular cloning and characterization of the complete acetylcholinesterase gene (Ace1) from the mosquito Aedes aegypti with implications for comparative genome analysis

Insensitive acetylcholinesterase (AChE) has been shown to be responsible for resistance to organophosphates and carbamates in a number of arthropod species. Some arthropod genomes contain a single Ace gene, while others including mosquitoes contain two genes, but only one confers insecticide resistance. Here we report the isolation of the full-length cDNA and characterization of the complete genomic DNA sequence for the Ace1 gene in the yellow fever mosquito, Aedes aegypti. The Ace1 homolog in other mosquito species has been associated with insecticide resistance. The full-length cDNA consists of 2721bp and contains a 2109bp open reading frame that encodes a 702 amino acid protein. The amino acid sequence is highly conserved with that of other mosquitoes, including greater than 90% identity with Culex spp. and about 80% identity with Anopheles gambiae. The genomic DNA sequence includes 138,970bp and consists of eight exons with seven introns ranging from 59 to 114,350bp. Exons 2 and 8 show reduced amino acid conservation across mosquito species, while exons 3-7 are highly conserved. The Ace1 introns in Ae. aegypti reflect a high frequency of repetitive sequences that comprise about 45% of the total intron sequence. The Ace1 locus maps to the p-arm of chromosome 3, which corresponds to the orthologous genome regions in Culex spp. and An. gambiae.     

Analysis of phylogenetic signal in protostomial intron patterns using Mutual Information

Many deep evolutionary divergences still remain unresolved, such as those among major taxa of the Lophotrochozoa. As alternative phylogenetic markers, the intron–exon structure of eukaryotic genomes and the patterns of absence and presence of spliceosomal introns appear to be promising. However, given the potential homoplasy of intron presence, the phylogenetic analysis of this data using standard evolutionary approaches has remained a challenge. Here, we used Mutual Information (MI) to estimate the phylogeny of Protostomia using gene structure data, and we compared these results with those obtained with Dollo Parsimony. Using full genome sequences from nine Metazoa, we identified 447 groups of orthologous sequences with 21,732 introns in 4,870 unique intron positions. We determined the shared absence and presence of introns in the corresponding sequence alignments and have made this data available in “IntronBase”, a web-accessible and downloadable SQLite database. Our results obtained using Dollo Parsimony are obviously misled through systematic errors that arise from multiple intron loss events, but extensive filtering of data improved the quality of the estimated phylogenies. Mutual Information, in contrast, performs better with larger datasets, but at the same time it requires a complete data set, which is difficult to obtain for orthologs from a large number of taxa. Nevertheless, Mutual Information-based distances proved to be useful in analyzing this kind of data, also because the estimation of MI-based distances is independent of evolutionary models and therefore no pre-definitions of ancestral and derived character states are necessary.     

Remarkable interkingdom conservation of intron positions and massive,lineage-specific intron loss and gain in eukaryotic evolution

Sequencing of eukaryotic genomes allows one to address major evolutionary problems, such as the evolution of gene structure. We compared the intron positions in 684 orthologous gene sets from 8 complete genomes of animals, plants, fungi, and protists and constructed parsimonious scenarios of evolution of the exon-intron structure for the respective genes. Approximately one-third of the introns in the malaria parasite Plasmodium falciparum are shared with at least one crown group eukaryote; this number indicates that these introns have been conserved through >1.5 billion years of evolution that separate Plasmodium from the crown group. Paradoxically, humans share many more introns with the plant Arabidopsis thaliana than with the fly or nematode. The inferred evolutionary scenario holds that the common ancestor of Plasmodium and the crown group and, especially, the common ancestor of animals, plants, and fungi had numerous introns. Most of these ancestral introns, which are retained in the genomes of vertebrates and plants, have been lost in fungi, nematodes, arthropods, and probably Plasmodium. In addition, numerous introns have been inserted into vertebrate and plant genes, whereas, in other lineages, intron gain was much less prominent.     

The complete sequence of the mouse skeletal alpha-actin gene reveals several conserved and inverted repeat sequences outside of the protein-coding region.

The complete nucleotide sequence of a genomic clone encoding the mouse skeletal alpha-actin gene has been determined. This single-copy gene codes for a protein identical in primary sequence to the rabbit skeletal alpha-actin. It has a large intron in the 5'-untranslated region 12 nucleotides upstream from the initiator ATG and five small introns in the coding region at codons specifying amino acids 41/42, 150, 204, 267, and 327/328. These intron positions are identical to those for the corresponding genes of chickens and rats. Similar to other skeletal alpha-actin genes, the nucleotide sequence codes for two amino acids, Met-Cys, preceding the known N-terminal Asp of the mature protein. Comparison of the nucleotide sequences of rat, mouse, chicken, and human skeletal muscle alpha-actin genes reveals conserved sequences (some not previously noted) outside of the protein-coding region. Furthermore, several inverted repeat sequences, partially within these conserved regions, have been identified. These sequences are not present in the vertebrate cytoskeletal beta-actin genes. The strong conservation of the inverted repeat sequences suggests that they may have a role in the tissue-specific expression of skeletal alpha-actin genes.     

Near intron positions are reliable phylogenetic markers: an application to holometabolous insects

Today, the reconstruction of the organismal evolutionary tree is based mainly on molecular sequence data. However, sequence data are sometimes insufficient to reliably resolve in particular deep branches. Thus, it is highly desirable to find novel, more reliable types of phylogenetic markers that can be derived from the wealth of genomic data. Here, we consider the gain of introns close to older preexisting ones. Because correct splicing is impeded by very small exons, nearby pairs of introns very rarely coexist, that is, the gain of the new intron is nearly always associated with the loss of the old intron. Both events may even be directly connected as in cases of intron migration. Therefore, it should be possible to identify one of the introns as ancient (plesiomorphic) and the other as novel (derived or apomorphic). To test the suitability of such near intron pairs (NIPs) as a marker class for phylogenetic analysis, we undertook an analysis of the evolutionary positions of bees and wasps (Hymenoptera) and beetles (Coleoptera) in relation to moths (Lepidoptera) and dipterans (Diptera) using recently completed genome project data. By scanning 758 putatively orthologous gene structures of Apis mellifera (Hymenoptera) and Tribolium castaneum (Coleoptera), we identified 189 pairs of introns, one from each species, which are located less than 50 nt from each other. A comparison with genes from 5 other holometabolan and 9 metazoan outgroup genomes resulted in 22 shared derived intron positions found in beetle as well as in butterflies and/or dipterans. This strongly supports a basal position of hymenopterans in the holometabolous insect tree. In addition, we found 31 and 12 intron positions apomorphic for A. mellifera and T. castaneum, respectively, which seem to represent changes inside these branches. Another 12 intron pairs indicate parallel intron gains or extraordinarily small exons. In conclusion, we show here that the analysis of phylogenetically nested, nearby intron pairs is suitable to identify evolutionarily younger intron positions and to determine their relative age, which should be of equal importance for the understanding of intron evolution and the reconstruction of the eukaryotic tree.     

Computer identification of snoRNA genes using a Mammalian Orthologous Intron Database

Based on comparative genomics, we created a bioinformatic package for computer prediction of small nucleolar RNA (snoRNA) genes in mammalian introns. The core of our approach was the use of the Mammalian Orthologous Intron Database (MOID), which contains all known introns within the human, mouse and rat genomes. Introns from orthologous genes from these three species, that have the same position relative to the reading frame, are grouped in a special orthologous intron table. Our program SNO.pl searches for conserved snoRNA motifs within MOID and reports all cases when characteristic snoRNA-like structures are present in all three orthologous introns of human, mouse and rat sequences. Here we report an example of the SNO.pl usage for searching a particular pattern of conserved C/D-box snoRNA motifs (canonical C- and D-boxes and the 6 nt long terminal stem). In this computer analysis, we detected 57 triplets of snoRNA-like structures in three mammals. Among them were 15 triplets that represented known C/D-box snoRNA genes. Six triplets represented snoRNA genes that had only been partially characterized in the mouse genome. One case represented a novel snoRNA gene, and another three cases, putative snoRNAs. Our programs are publicly available and can be easily adapted and/or modified for searching any conserved motifs within mammalian introns.     

Intron position as an evolutionary marker of thioredoxins and thioredoxin domains

In contrast to prokaryotes, which typically possess one thioredoxin gene per genome, three different thioredoxin types have been described in higher plants. All are encoded by nuclear genes, but thioredoxins m and f are chloroplastic while thioredoxins h have no transit peptide and are probably cytoplasmic. We have cloned and sequencedArabidopsis thaliana genomic fragments encoding the five previously described thioredoxins h, as well as a sixth gene encoding a new thioredoxin h. In spite of the high divergence of the sequences, five of them possess two introns at positions identical to the previously sequenced tobacco thioredoxin h gene, while a single one has only the first intron. The recently published sequence ofChlamydomonas thioredoxin h shows three introns, two at the same positions as in higher plants. This strongly suggests a common origin for all cytoplasmic thioredoxins of plants and green algae. In addition, we have cloned and sequenced pea DNA genomic fragments encoding thioredoxins m and f. The thioredoxin m sequence shows only one intron between the regions encoding the transit peptide and the mature protein, supporting the prokaryotic origin of this sequence and suggesting that its association with the transit peptide has been facilitated by exon shuffling. In contrast, the thioredoxin f sequence shows two introns, one at the same position as an intron in various plant and animal thioredoxins and the second at the same position as an intron in thioredoxin domains of disulfide isomerases. This strongly supports the hypothesis of a eukaryotic origin for chloroplastic thioredoxin f.     

Conservation versus parallel gains in intron evolution

Orthologous genes from distant eukaryotic species, e.g. animals and plants, share up to 25–30% intron positions. However, the relative contributions of evolutionary conservation and parallel gain of new introns into this pattern remain unknown. Here, the extent of independent insertion of introns in the same sites (parallel gain) in orthologous genes from phylogenetically distant eukaryotes is assessed within the framework of the protosplice site model. It is shown that protosplice sites are no more conserved during evolution of eukaryotic gene sequences than random sites. Simulation of intron insertion into protosplice sites with the observed protosplice site frequencies and intron densities shows that parallel gain can account but for a small fraction (5–10%) of shared intron positions in distantly related species. Thus, the presence of numerous introns in the same positions in orthologous genes from distant eukaryotes, such as animals, fungi and plants, appears to reflect mostly bona fide evolutionary conservation.     

Nucleotide sequence of the triosephosphate isomerase gene from Aspergillus nidulans: implications for a differential loss of introns

A functional cDNA from Aspergillus nidulans encoding triosephosphate isomerase (TPI) was isolated by its ability to complement a tpi1 mutation in Saccharomyces cerevisiae. This cDNA was used to obtain the corresponding gene, tpiA. Alignment of the cDNA and genomic DNA nucleotide sequences indicated that tpiA contains five introns. The intron positions in the tpiA gene were compared with those in the TPI genes of human, chicken, and maize. One intron is present at an identical position in all four organisms, two other introns are located in similar positions in A. nidulans and maize, and the remaining two introns are unique to A. nidulans. These Aspergillus-specific introns are located in regions of the protein that were predicted to be interrupted by introns based on analysis of a Go plot of chicken TPI. These comparisons are discussed in relation to the evolution of introns within TPI genes.     

Structure of the gene encoding potato cytosolic pyruvate kinase.

The polymerase chain reaction (PCR) has been used to generate a series of overlapping genomic clones representing 43 bp of 5' untranslated sequence, 63 bp of 3' untranslated sequence and the entire coding sequence of the gene encoding potato cytosolic pyruvate kinase (PKc). This portion of the gene is approximately 4.5 kb in length and is interrupted by three introns, one of which is present in the 5' untranslated region. Southern blot analysis indicates that PKc is encoded by a small gene family, and sequence data from a number of PCR-derived genomic clones indicate that there are as many as six PKc genes. Sequence differences between the PCR-generated genomic clones and a PKc cDNA clone are discussed with respect to the fidelity of Taq polymerase. An alignment of intron placement in the potato PKc gene with intron placement in PK genes from other sources indicates that two of the potato introns correspond to intron positions in other species.     

<Emphasis Type="Italic">Oscheius tipulae</Emphasis> as an Example of eEF1A Gene Diversity in Nematodes

We characterized four eEF1A genes in the alternative rhabditid nematode model organism Oscheius tipulae. This is twice the copy number of eEF1A genes in C. elegans, C. briggsae, and, probably, many other free-living and parasitic nematodes. The introns show features remarkably different from those of other metazoan eEF1A genes. Most of the introns in the eEF1A genes are specific to O. tipulae and are not shared with any of the other genes described in metazoans. Most of the introns are phase 0 (inserted between two codons), and few are inserted in protosplice sites (introns inserted between the nucleotide sequence A/CAG and G/A). Two of these phase 0 introns are conserved in sequence in two or more of the four eEF1A gene copies, and are inserted in the same position in the genes. Neither of these characteristics has been detected in any of the nematode eEF1A genes characterized to date. The coding sequences were also compared with other eEF1A cDNAs from 11 different nematodes to determine the variability of these genes within the phylum Nematoda. Parsimony and distance trees yielded similar topologies, which were similar to those created using other molecular markers. The presence of more than one copy of the eEF1A gene with nearly identical coding regions makes it difficult to define the orthologous cDNAs. As shown by our data on O. tipulae, careful and extensive examination of intron positions in the eEF1A gene across the phylum is necessary to define their potential for use as valid phylogenetic markers.     

The intron-containing L3 ribosomal protein gene (RPL3): sequence analysis and identification of U43 and of two novel intronic small nucleolar RNAs

Isolation and sequencing of bovine and human intron-containing L3 ribosomal protein genes are here reported. They exhibit very similar organisation, both comprising 10 exons and nine introns. A polymorphic locus, involving a 19-bp deletion, was found in intron 6 of the human gene. The frequency of the two alleles has been estimated in 200 haploid genomes. In bovine and human genes intron sequences are rather different, except for limited regions, located in corresponding positions, which show a surprisingly high degree of identity. All these regions contain conserved features defining the box C/D class of small nucleolar RNAs. Demonstration is given that U43 small nucleolar RNA is encoded within the first intron of both bovine and human genes. Single nucleotide sequences, encoding two novel species of small nucleolar RNAs (U82, U83a and U83b), are located in introns 3, 5 and 7. Their expression has been investigated and a possible role of these molecules in 2'-O-ribose methylation of rRNAs is discussed.