Evolutionary convergence on highly-conserved 3' intron structures in intron-poor eukaryotes and insights into the ancestral eukaryotic genome

The presence of spliceosomal introns in eukaryotes raises a range of questions about genomic evolution. Along with the fundamental mysteries of introns' initial proliferation and persistence, the evolutionary forces acting on intron sequences remain largely mysterious. Intron number varies across species from a few introns per genome to several introns per gene, and the elements of intron sequences directly implicated in splicing vary from degenerate to strict consensus motifs. We report a 50-species comparative genomic study of intron sequences across most eukaryotic groups. We find two broad and striking patterns. First, we find that some highly intron-poor lineages have undergone evolutionary convergence to strong 3' consensus intron structures. This finding holds for both branch point sequence and distance between the branch point and the 3' splice site. Interestingly, this difference appears to exist within the genomes of green alga of the genus Ostreococcus, which exhibit highly constrained intron sequences through most of the intron-poor genome, but not in one much more intron-dense genomic region. Second, we find evidence that ancestral genomes contained highly variable branch point sequences, similar to more complex modern intron-rich eukaryotic lineages. In addition, ancestral structures are likely to have included polyT tails similar to those in metazoans and plants, which we found in a variety of protist lineages. Intriguingly, intron structure evolution appears to be quite different across lineages experiencing different types of genome reduction: whereas lineages with very few introns tend towards highly regular intronic sequences, lineages with very short introns tend towards highly degenerate sequences. Together, these results attest to the complex nature of ancestral eukaryotic splicing, the qualitatively different evolutionary forces acting on intron structures across modern lineages, and the impressive evolutionary malleability of eukaryotic gene structures.     

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.     

Hard tick calreticulin (CRT) gene coding regions have only one intron with conserved positions and variable sizes

Calreticulin (CRT) is a unique eukaryotic gene. The CRT gene product, calreticulin, was first identified as a calcium binding protein in 1974, but further investigations have indicated that CRT protein performs many functions in cells, including involvement in evading the host's immune system by parasites. Many studies of CRT have been published since the molecule was first discovered; however, the CRT gene exon-intron structure is only known for a limited number of ectoparasite species. In this study, we compared tick CRT genomic sequences to the corresponding cDNA from 28 species and found that 2 exons and 1 intron are present in the tick CRT gene. The intron position is conserved in 28 hard ticks, but intron size and nucleotide sequences vary. Three tick introns possess duplicated fragments and are twice as long as other introns. All tick CRT introns obey the GT-AG rule in the splice-site junctions and are phase 1 introns. By comparing tick CRT introns to those of fruit fly, mouse, and human, we conclude that tick CRT introns belong to the intron-late type. The number and size of CRT introns have increased through the evolution of eukaryotes.     

Spliceosomal introns as tools for genomic and evolutionary analysis

Over the past 5 years, the availability of dozens of whole genomic sequences from a wide variety of eukaryotic lineages has revealed a very large amount of information about the dynamics of intron loss and gain through eukaryotic history, as well as the evolution of intron sequences. Implicit in these advances is a great deal of information about the structure and evolution of surrounding sequences. Here, we review the wealth of ways in which structures of spliceosomal introns as well as their conservation and change through evolution may be harnessed for evolutionary and genomic analysis. First, we discuss uses of intron length distributions and positions in sequence assembly and annotation, and for improving alignment of homologous regions. Second, we review uses of introns in evolutionary studies, including the utility of introns as indicators of rates of sequence evolution, for inferences about molecular evolution, as signatures of orthology and paralogy, and for estimating rates of nucleotide substitution. We conclude with a discussion of phylogenetic methods utilizing intron sequences and positions.     

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     

Aberrant pre-mRNA maturation is caused by LINE insertions into introns of the white gene of Drosophila melanogaster.

Insertional mutagenesis screens have provided thousands of mutant alleles for analysing genes of varied functions in Drosophila melanogaster. We here document mechanisms of insertional mutagenesis by a LINE element, the I factor, by determining the molecular structure of RNAs produced from two alleles of the white gene of D.melanogaster, wIR1 and wIR6. These alleles result from insertion of the I factor into introns of the gene. We show that sequences present within the element direct aberrant splicing and termination events. When the I factor is inserted within the white first intron it may lead to the use of a cryptic 3' splice site which does not contain the dinucleotide AG. This splicing gives rise to a chimeric messenger RNA whose synthesis is controlled differently in tissues where the mutated gene is expressed. When the I factor is inserted within the white last intron it induces synthesis of truncated mRNAs. These results provide, for the first time, mechanisms for I factor insertional mutagenesis. They are discussed in the more general context of RNA processing in Drosophila and the evolution of eukaryotic gene 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.     

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.     

Patterns of intron loss and gain in plants: intron loss-dominated evolution and genome-wide comparison of O. sativa and A. thaliana

Numerous previous studies have elucidated 2 surprising patterns of spliceosomal intron evolution in diverse eukaryotes over the past roughly 100 Myr. First, rates of recent intron gain in a wide variety of eukaryotic lineages have been surprisingly low, far too low to explain modern intron densities. Second, intron losses have outnumbered intron gains over a variety of lineages. For several reasons, land plants might be expected to have comparatively high rates of intron gain and thus to represent a possible exception to this pattern. However, we report several studies that indicate low rates of intron gain and an excess of intron losses over intron gains in a variety of plant lineages. We estimate that intron losses have outnumbered intron gains in recent evolution in Arabidopsis thaliana (roughly 12.6 times more losses than gains), Oryza sativa (9.8 times), the green alga Chlamydomonas reinhardtii (5.1 times), and the Bigelowiella natans nucleomorph, an enslaved green algal nucleus (2.8 times). We estimate that during recent evolution, A. thaliana and O. sativa have experienced very low rates of intron gain of around one gain per gene per 2.6-8.0 billion years. In addition, we compared 8,258 pairs of putatively orthologous A. thaliana-O. sativa genes. We found that 5.3% of introns in conserved coding regions are species-specific. Observed species-specific A. thaliana and O. sativa introns tend to be exact and to lie adjacent to each other along the gene, in a pattern suggesting mRNA-mediated intron loss. Our results underscore that low intron gain rates and intron number reduction are common features of recent eukaryotic evolution. This pattern implies that rates of intron creation were higher during earlier periods of evolution and further focuses attention on the causes of initial intron proliferation.     

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.     

Ovomucoid intervening sequences specify functional domains and generate protein polymorphism

Two thirds of the natural chicken ovomucoid gene has been sequenced, including all exons and the intron sequences surrounding all fourteen intron/ exon junctions. The junction sequences surrounding four of the introns are redundant; however, the sequences surrounding the other three introns contain no redundancies and thus the splicing sites at either end of these three introns are unambiguous. The splicing in all cases conforms to the GT-AG rule. The ovomucoid gene sequence around intron F can be used to predict the cause of an internal deletion polymorphism in the ovomucoid protein, which is an apparent error in the processing of the ovomucoid pre-mRNA. We also compare the structural organization of the ovomucoid gene with the ovomucoid protein sequence to examine theories of the evolution of ovomucoids as well as the origin of intervening sequences. This analysis suggests that the present ovomucoid gene evolved from a primordial ovomucoid gene by two separate intragenic duplications. Furthermore, sequence analyses suggest that introns were present in the primordial ovomucoid gene before birds and mammals diverged, about 300 million years ago. Finally, the positions of the introns within the ovomucoid gene support the theory that introns separate gene segments that code for functional domains of proteins and provide insight on the manner by which eucaryotic genes were constructed during the process of evolution.     

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.     

Prevalence of intron gain over intron loss in the evolution of paralogous gene families

The mechanisms and evolutionary dynamics of intron insertion and loss in eukaryotic genes remain poorly understood. Reconstruction of parsimonious scenarios of gene structure evolution in paralogous gene families in animals and plants revealed numerous gains and losses of introns. In all analyzed lineages, the number of acquired new introns was substantially greater than the number of lost ancestral introns. This trend held even for lineages in which vertical evolution of genes involved more intron losses than gains, suggesting that gene duplication boosts intron insertion. However, dating gene duplications and the associated intron gains and losses based on the molecular clock assumption showed that very few, if any, introns were gained during the last ~100 million years of animal and plant evolution, in agreement with previous conclusions reached through analysis of orthologous gene sets. These results are generally compatible with the emerging notion of intensive insertion and loss of introns during transitional epochs in contrast to the relative quiet of the intervening evolutionary spans.     

A 25-bp ancient spliceosomal intron in the TvRab1a gene of Trichomonas vaginalis

Spliceosomal introns play a key role in eukaryotic genome evolution and protein diversity. A large Rab GTPase family has been identified in a unicellular eukaryote Trichomonas vaginalis. However, the characteristics of introns in Rab genes of T. vaginalis have not been investigated previously. In this study, we identified a 25-bp spliceosomal intron in the T. vaginalis Rab1a (TvRab1a) gene, the smallest intron in T. vaginalis to be characterized to date. This intron contains a canonical splice site at both 5' (GT) and 3' (AG) ends, and a putative branch-point sequence (TCTAAC) that matches the Trichomonad consensus sequence of ACTAAC except for the first nucleotide. The position and phase of the TvRab1a intron are evolutionarily conserved in Rab1 homologous genes across at least five eukaryotic supergroups, including Opisthokonta, Amoebozoa, Excavata, Chromalveolata, and Plantae. These results strongly suggest that the TvRab1a intron is likely to be an ancient spliceosomal intron, and it can therefore be used as a phylogenetic marker to evaluate particular eukaryotic groupings. Identification and characterization of the TvRabla intron may provide an insight into the evolution of the large Rab repertoire in T. vaginalis.     

Intron mobility in phage T4 is dependent upon a distinctive class of endonucleases and independent of DNA sequences encoding the intron core: mechanistic and evolutionary implications

Although mobility of the phylogenetically widespread group I introns appears to be mechanistically similar, the phage T4 intron-encoded endonucleases that promote mobility of the td and sunY introns are different from their eukaryotic counterparts. Most notably, they cleave at a distance from the intron insertion sites. The td enzyme was shown to cleave 23-26 nt 5' and the sunY endonuclease 13-15 nt 3' to the intron insertion site to generate 3-nt or 2-nt 3'-OH extensions, respectively. The absolute coconversion of exon markers between the distant cleavage and insertion sites is consistent with the double-strand-break repair model for intron mobility. As a further critical test of the model we have demonstrated that the mobility event is independent of DNA sequences that encode the catalytic intron core structure. Thus, in derivatives in which the lacZ or kanR coding sequences replace the intron, these marker genes are efficiently inserted into intron-minus alleles when the cognate endonuclease is provided in trans. The process is therefore endonuclease-dependent, rather than dependent on the intron per se. These findings, which imply that the endonucleases rather than the introns themselves were the primordial mobile elements, are incorporated into a model for the evolution of mobile introns.