New maximum likelihood estimators for eukaryotic intron evolution

The evolution of spliceosomal introns remains poorly understood. Although many approaches have been used to infer intron evolution from the patterns of intron position conservation, the results to date have been contradictory. In this paper, we address the problem using a novel maximum likelihood method, which allows estimation of the frequency of intron insertion target sites, together with the rates of intron gain and loss. We analyzed the pattern of 10,044 introns (7,221 intron positions) in the conserved regions of 684 sets of orthologs from seven eukaryotes. We determined that there is an average of one target site per 11.86 base pairs (bp) (95% confidence interval, 9.27 to 14.39 bp). In addition, our results showed that: (i) overall intron gains are ~25% greater than intron losses, although specific patterns vary with time and lineage; (ii) parallel gains account for ~18.5% of shared intron positions; and (iii) reacquisition following loss accounts for ~0.5% of all intron positions. Our results should assist in resolving the long-standing problem of inferring the evolution of spliceosomal introns.     

Widespread occurrence of spliceosomal introns in the rDNA genes of ascomycetes

Spliceosomal (pre-mRNA) introns have previously been found in eukaryotic protein-coding genes, in the small nuclear RNAs of some fungi, and in the small- and large-subunit ribosomal DNA genes of a limited number of ascomycetes. How the majority of these introns originate remains an open question because few proven cases of recent and pervasive intron origin have been documented. We report here the widespread occurrence of spliceosomal introns (69 introns at 27 different sites) in the small- and large-subunit nuclear-encoded rDNA of lichen-forming and free-living members of the Ascomycota. Our analyses suggest that these spliceosomal introns are of relatively recent origin, i.e., within the Euascomycetes, and have arisen through aberrant reverse-splicing (in trans) of free pre-mRNA introns into rRNAs. The spliceosome itself, and not an external agent (e.g., transposable elements, group II introns), may have given rise to these introns. A nonrandom sequence pattern was found at sites flanking the rRNA spliceosomal introns. This pattern (AG-intron-G) closely resembles the proto-splice site (MAG-intron-R) postulated for intron insertions in pre-mRNA genes. The clustered positions of spliceosomal introns on secondary structures suggest that particular rRNA regions are preferred sites for insertion through reverse-splicing.     

Molecular evolution of eukaryotic genomes: hemiascomycetous yeast spliceosomal introns

As part of the exploratory sequencing program Génolevures, visual scrutinisation and bioinformatic tools were used to detect spliceosomal introns in seven hemiascomycetous yeast species. A total of 153 putative novel introns were identified. Introns are rare in yeast nuclear genes (<5% have an intron), mainly located at the 5′ end of ORFs, and not highly conserved in sequence. They all share a clear non-random vocabulary: conserved splice sites and conserved nucleotide contexts around splice sites. Homologues of metazoan snRNAs and putative homologues of SR splicing factors were identified, confirming that the spliceosomal machinery is highly conserved in eukaryotes. Several introns’ features were tested as possible markers for phylogenetic analysis. We found that intron sizes vary widely within each genome, and according to the phylogenetic position of the yeast species. The evolutionary origin of spliceosomal introns was examined by analysing the degree of conservation of intron positions in homologous yeast genes. Most introns appeared to exist in the last common ancestor of present day yeast species, and then to have been differentially lost during speciation. However, in some cases, it is difficult to exclude a possible sliding event affecting a pre-existing intron or a gain of a novel intron. Taken together, our results indicate that the origin of spliceosomal introns is complex within a given genome, and that present day introns may have resulted from a dynamic flux between intron conservation, intron loss and intron gain during the evolution of hemiascomycetous yeasts.     

USING RDNA GENES TO UNDERSTAND THE ORIGIN AND EVOLUTION OF SPLICEOSOMAL AND GROUP I INTRONS

Ribosomal DNA genes in lichen algae and lichen fungi are astonishingly rich in spliceosomal and group I introns. We use phylogenetic, secondary structure, and biochemical analyses to understand the evolution of these introns. Despite the widespread distribution of spliceosomal introns in nuclear pre-mRNA genes, their general mechanism of origin remains an open question because few proven cases of recent and pervasive intron origin have been documented. The lichen introns are valuable in this respect because they are undoubtedly of a "recent" origin and limited to the Euascomycetes. Our analyses suggest that rDNA spliceosomal introns have arisen through aberrant reverse-splicing (in trans) of free pre-mRNA introns into r RNAs. We propose that the spliceosome itself (and not an external agent; e.g. transposable elements, group II introns) has given rise to the introns. The rDNA introns are found most often between the flanking sequence G (78%) - intron-G (72%), and their clustered positions on secondary structures suggest that particular r RNA regions are preferred sites (i.e., proto-splice sites) for insertion. Mapping of intron positions on the newly available tertiary structures show that they are found most often in exposed regions of the ribosomes. This again is consistent with an intron origin through reverse-splicing. Remarkably, the distribution and phylogenetic relationships of most group I introns in nuclear rDNA genes are also consistent with a reverse-splicing origin. These data underline the value of lichens as a model system for understanding intron origin and stress the importance of RNA-level processes in the spread of these sequences in nuclear coding regions.     

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.     

Evolutionary dynamics of spliceosomal intron revealed by <Emphasis Type="Italic">in silico</Emphasis> analyses of the P-Type ATPase superfamily genes

It has been long debated whether spliceosomal introns originated in the common ancestor of eukaryotes and prokaryotes. In this study, we tested the possibility that extant introns were inherited from the common ancestor of eukaryotes and prokaryotes using in silico simulation. We first identified 21 intron positions that are shared among different families of the P-Type ATPase superfamily, some of which are known to have diverged before the separation of prokaryotes and eukaryotes. Theoretical estimates of the expected number of intron positions shared by different genes suggest that the introns at those 21 positions were inserted independently. There seems to be no intron that arose from before the diversification of the P-Type ATPase superfamily. Namely, the present introns were inserted after the separation of eukaryotes and prokaryotes.     

Testing the "proto-splice sites" model of intron origin: evidence from analysis of intron phase correlations

A few nucleotide sites of nuclear exons that flank introns are often conserved. A hypothesis has suggested that these sites, called "proto-splice sites," are remnants of recognition signals for the insertion of introns in the early evolution of eukaryotic genes. This notion of proto-splice sites has been an important basis for the insertional theory of introns. This hypothesis predicts that the distribution of proto-splice sites would determine the distribution of intron phases, because the positions of introns are just a subset of the proto-splice sites. We previously tested this prediction by examining the proportions of the phases of proto-splice sites, revealing nothing in these proportion distributions similar to observed proportions of intron phases. Here, we provide a second independent test of the proto-splice site hypothesis, with regard to its prediction that the proto-splice sites would mimic intron phase correlations, using a CDS database we created from GenBank. We tested four hypothetical proto-splice sites G / G, AG / G, AG / GT, and C/AAG / R. Interestingly, while G / G and AG / GT site phase distributions are not consistent with actual introns, we observed that AG / G and C/AAG / R sites have a symmetric phase excess. However, the patterns of the excess are quite different from the actual intron phase distribution. In addition, particular amino acid repeats in proteins were found to partially contribute to the excess of symmetry at these two types of sites. The phase associations of all four sites are significantly different from those of intron phases. Furthermore, a general model of intron insertion into proto-splice sites was simulated by Monte Carlo simulation to investigate the probability that the random insertion of introns into AG / G and C/AAG / R sites could generate the observed intron phase distribution. The simulation showed that (1) no observed correlation of intron phases was statistically consistent with the phase distribution of proto-splice sites in the simulated virtual genes; (2) most conservatively, no simulation in 10,000 Monte Carlo experiments gave a pattern with an excess of symmetric (1, 1) exons larger than those of (0, 0) and (2, 2), a major statistical feature of intron phase distribution that is consistent with the directly observed cases of exon shuffling. Thus, these results reject the null hypothesis that introns are randomly inserted into preexisting proto-splice sites, as suggested by the insertional theory of introns.     

Evidence of splice signal migration from exon to intron during intron evolution

A comparison of the nucleotide sequences around the splice junctions that flank old (shared by two or more major lineages of eukaryotes) and new (lineage-specific) introns in eukaryotic genes reveals substantial differences in the distribution of information between introns and exons. Old introns have a lower information content in the exon regions adjacent to the splice sites than new introns but have a corresponding higher information content in the intron itself. This suggests that introns insert into nonrandom (proto-splice) sites but, during the evolution of an intron after insertion, the splice signal shifts from the flanking exon regions to the ends of the intron itself. Accumulation of information inside the intron during evolution suggests that new introns largely emerge de novo rather than through propagation and migration of old introns.     

A very high fraction of unique intron positions in the intron-rich diatom Thalassiosira pseudonana indicates widespread intron gain

Although spliceosomal introns are present in all characterized eukaryotes, intron numbers vary dramatically, from only a handful in the entire genomes of some species to nearly 10 introns per gene on average in vertebrates. For all previously studied intron-rich species, significant fractions of intron positions are shared with other widely diverged eukaryotes, indicating that 1) large numbers of the introns date to much earlier stages of eukaryotic evolution and 2) these lineages have not passed through a very intron-poor stage since early eukaryotic evolution. By the same token, among species that have lost nearly all of their ancestral introns, no species is known to harbor large numbers of more recently gained introns. These observations are consistent with the notion that intron-dense genomes have arisen only once over the course of eukaryotic evolution. Here, we report an exception to this pattern, in the intron-rich diatom Thalassiosira pseudonana. Only 8.1% of studied T. pseudonana intron positions are conserved with any of a variety of divergent eukaryotic species. This implies that T. pseudonana has both 1) lost nearly all of the numerous introns present in the diatom-apicomplexan ancestor and 2) gained a large number of new introns since that time. In addition, that so few apparently inserted T. pseudonana introns match the positions of introns in other species implies that insertion of multiple introns into homologous genic sites in eukaryotic evolution is less common than previously estimated. These results suggest the possibility that intron-rich genomes may have arisen multiple times in evolution. These results also provide evidence that multiple intron insertion into the same site is rare, further supporting the notion that early eukaryotic ancestors were very intron rich.     

Distribution of rRNA introns in the three-dimensional structure of the ribosome

More than 1200 introns have been documented at over 150 unique sites in the small and large subunit ribosomal RNA genes (as of February 2002). Nearly all of these introns are assigned to one of four main types: group I, group II, archaeal and spliceosomal. This sequence information has been organized into a relational database that is accessible through the Comparative RNA Web Site (http://www.rna.icmb.utexas.edu/) While the rRNA introns are distributed across the entire tree of life, the majority of introns occur within a few phylogenetic groups. We analyzed the distributions of rRNA introns within the three-dimensional structures of the 30S and 50S ribosomes. Most sites in rRNA genes that contain introns contain only one type of intron. While the intron insertion sites occur at many different coordinates, the majority are clustered near conserved residues that form tRNA binding sites and the subunit interface. Contrary to our expectations, many of these positions are not accessible to solvent in the mature ribosome. The correlation between the frequency of intron insertions and proximity of the insertion site to functionally important residues suggests an association between intron evolution and rRNA function.     

Phylogenetic distribution of intron positions in alpha-amylase genes of bilateria suggests numerous gains and losses

Most eukaryotes have at least some genes interrupted by introns. While it is well accepted that introns were already present at moderate density in the last eukaryote common ancestor, the conspicuous diversity of intron density among genomes suggests a complex evolutionary history, with marked differences between phyla. The question of the rates of intron gains and loss in the course of evolution and factors influencing them remains controversial. We have investigated a single gene family, alpha-amylase, in 55 species covering a variety of animal phyla. Comparison of intron positions across phyla suggests a complex history, with a likely ancestral intronless gene undergoing frequent intron loss and gain, leading to extant intron/exon structures that are highly variable, even among species from the same phylum. Because introns are known to play no regulatory role in this gene and there is no alternative splicing, the structural differences may be interpreted more easily: intron positions, sizes, losses or gains may be more likely related to factors linked to splicing mechanisms and requirements, and to recognition of introns and exons, or to more extrinsic factors, such as life cycle and population size. We have shown that intron losses outnumbered gains in recent periods, but that "resets" of intron positions occurred at the origin of several phyla, including vertebrates. Rates of gain and loss appear to be positively correlated. No phase preference was found. We also found evidence for parallel gains and for intron sliding. Presence of introns at given positions was correlated to a strong protosplice consensus sequence AG/G, which was much weaker in the absence of intron. In contrast, recent intron insertions were not associated with a specific sequence. In animal Amy genes, population size and generation time seem to have played only minor roles in shaping gene structures.     

The evolution of spliceosomal introns in alveolates

Many issues concerning the evolution of spliceosomal introns remain poorly understood. In this respect, the reconstruction of the evolution of introns in deep branching species such as alveolates is of special significance. In this study, we inferred the intron evolution in alveolates using 3,368 intron positions in 162 orthologs from 10 species (9 alveolates and 1 outgroup, Homo sapiens). We found that although very few intron gains and losses have occurred in Theileria and Plasmodium recently, many intron gains and losses have occurred in the evolution of alveolates. Thus, the rates of intron gain and loss in alveolates have varied greatly across time and lineage. Our results seem to support the notion that massive intron gains and losses have occurred during short episodes, perhaps coinciding with major evolutionary events.     

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.     

The recent origins of spliceosomal introns revisited

Does the intron/exon structure of eukaryotic genes belie their ancient assembly by exon-shuffling or have introns been inserted into preformed genes during eukaryotic evolution? These are the central questions in the ongoing ‘introns-early’ versus ‘introns-late’ controversy. The phylogenetic distribution of spliceosomal introns continues to strongly favor the intronslate theory. The introns-early theory, however, has claimed support from intron phase and protein structure correlations.     

Sm/Lsm Genes Provide a Glimpse into the Early Evolution of the Spliceosome

The spliceosome, a sophisticated molecular machine involved in the removal of intervening sequences from the coding sections of eukaryotic genes, appeared and subsequently evolved rapidly during the early stages of eukaryotic evolution. The last eukaryotic common ancestor (LECA) had both complex spliceosomal machinery and some spliceosomal introns, yet little is known about the early stages of evolution of the spliceosomal apparatus. The Sm/Lsm family of proteins has been suggested as one of the earliest components of the emerging spliceosome and hence provides a first in-depth glimpse into the evolving spliceosomal apparatus. An analysis of 335 Sm and Sm-like genes from 80 species across all three kingdoms of life reveals two significant observations. First, the eukaryotic Sm/Lsm family underwent two rapid waves of duplication with subsequent divergence resulting in 14 distinct genes. Each wave resulted in a more sophisticated spliceosome, reflecting a possible jump in the complexity of the evolving eukaryotic cell. Second, an unusually high degree of conservation in intron positions is observed within individual orthologous Sm/Lsm genes and between some of the Sm/Lsm paralogs. This suggests that functional spliceosomal introns existed before the emergence of the complete Sm/Lsm family of proteins; hence, spliceosomal machinery with considerably fewer components than today's spliceosome was already functional.     

Very little intron gain in Entamoeba histolytica genes laterally transferred from prokaryotes

The evolution of spliceosomal introns remains intensely debated. We studied 96 Entamoeba histolytica genes previously identified as having been laterally transferred from prokaryotes, which were presumably intronless at the time of transfer. Ninety out of the 96 are also present in the reptile parasite Entamoeba invadens, indicating lateral transfer before the species' divergence approximately 50 MYA. We find only 2 introns, both shared with E. invadens. Thus, no intron gains have occurred in approximately 50 Myr, implying a very low rate of intron gain of less than one gain per gene per approximately 4.5 billion years. Nine other predicted introns are due to annotation errors reflecting apparent mistakes in the E. histolytica genome assembly. These results underscore the massive differences in intron gain rates through evolution.     

Proto-splice site model of intron origin

It is proposed that nuclear pre-mRNA introns (classical introns) were first generated as by-products during the evolution of alternative splicing. They were formed whenever two splice sites within the coding sequence of ancestral genes were used at a frequency that removed the coding constraint from the intervening sequence. Once introns had evolved, it is suggested that they were spread by the splicing machinery which inserted them into proto or cryptic-splice sites of other genes by reverse splicing, so giving rise to genes that have introns yet are not alternatively spliced. It is argued that 5' and 3' splice sites evolved from common ancestral splice sites, referred to as proto-splice sites, that were bidirectional and had a core consensus sequence of C or A, A, G, R, which remains today as the immediate flanking sequence of most introns. The ancestral splicing machinery, although inefficient, would have been capable of generating vast mRNA diversity by splicing between proto-splice sites. Natural selection would be expected to have preserved mutations that increased the amounts of advantageously spliced mRNA. It is argued that this process drove the evolution of present 5' and 3' splice sites from a subset of proto-splice sites and also drove the evolution of a more efficient splicing machinery. The positions of most introns that evolved directly from the coding sequence would be expected to correlate with protein structure.     

Insertion of spliceosomal introns in proto-splice sites: the case of secretory signal peptides

Analysis of the exon-intron structures of 2208 human genes has revealed that there is a statistically highly significant excess of phase 1 introns in the vicinity of the signal peptide cleavage sites. It is suggested that amino acid sequences surrounding signal peptide cleavage sites are significantly enriched in phase 1 proto-splice sites and this has favored insertion of spliceosomal introns in these sites.     

Widespread intron loss suggests retrotransposon activity in ancient apicomplexans

Several facets of spliceosomal intron in apicomplexans remain mysterious. First, intron numbers vary across species by 2 orders of magnitude, indicating massive intron loss and/or gain. Second, previous studies have shown very different evolutionary patterns over different timescales, with 100-fold higher rates of intron loss/gain between genera than within genera. Third, the timing and dynamics of nearly complete intron loss in Cryptosporidium species, as well as reasons for retention of the few remaining introns, remain unknown. We compared intron positions in 785 orthologous genes between 3 moderate to intron-rich apicomplexan species. We estimate that the Theileria-Plasmodium ancestor had 4.5 times as many introns as modern Plasmodium species and 38% more than modern Theileria species, and that subsequent intron losses have outnumbered intron gains by 5.8 to 1 in Theileria and by some 56 to 1 in Plasmodium. Several patterns suggest that these intron losses occurred by recombination with reverse-transcribed mRNAs. Intriguingly, this finding suggests significant retrotransposon activity in the lineages leading to both Theileria and Plasmodium, in contrast to the dearth of known retrotransposons and intron loss within modern species from both genera. We also compared genomes from Cryptosporidium parvum and C. hominis and found no evidence of ongoing intron loss, nor of intron gain. By contrast, Cryptosporidium introns are less evolutionary conserved with Toxoplasma than are introns from other apicomplexans; thus the few remaining introns are not simply indispensable ancestral introns.