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.     

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.     

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.     

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.     

Conserved intron positions in FGFR genes reflect the modular structure of FGFR and reveal stepwise addition of domains to an already complex ancestral FGFR

We have analyzed the evolution of fibroblast growth factor receptor (FGFR) tyrosine kinase genes throughout a wide range of animal phyla. No evidence for an FGFR gene was found in Porifera, but we tentatively identified an FGFR gene in the placozoan Trichoplax adhaerens. The gene encodes a protein with three immunoglobulin-like domains, a single-pass transmembrane, and a split tyrosine kinase domain. By superimposing intron positions of 20 FGFR genes from Placozoa, Cnidaria, Protostomia, and Deuterostomia over the respective protein domain structure, we identified ten ancestral introns and three conserved intron groups. Our analysis shows (1) that the position of ancestral introns correlates to the modular structure of FGFRs, (2) that the acidic domain very likely evolved in the last common ancestor of triploblasts, (3) that splicing of IgIII was enabled by a triploblast-specific insertion, and (4) that IgI is subject to substantial loss or duplication particularly in quickly evolving genomes. Moreover, intron positions in the catalytic domain of FGFRs map to the borders of protein subdomains highly conserved in other serine/threonine kinases. Nevertheless, these introns were introduced in metazoan receptor tyrosine kinases exclusively. Our data support the view that protein evolution dating back to the Cambrian explosion took place in such a short time window that only subtle changes in the domain structure are detectable in extant representatives of animal phyla. We propose that the first multidomain FGFR originated in the last common ancestor of Placozoa, Cnidaria, and Bilateria. Additional domains were introduced mainly in the ancestor of triploblasts and in the Ecdysozoa.     

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.     

Investigation of loss and gain of introns in the compact genomes of pufferfishes (Fugu and Tetraodon)

We have investigated intron evolution in the compact genomes of 2 closely related species of pufferfishes, Fugu rubripes and Tetraodon nigroviridis, that diverged about 32 million years ago (MYA). Analysis of 148,028 aligned intron positions in 13,547 gene pairs using human as an outgroup identified 57 and 24 intron losses in Tetraodon and fugu lineages, respectively, and no gain in either lineage. For comparison, we analyzed 144,545 intron positions in 12,866 orthologous pairs of genes in human and mouse that diverged about 61 MYA using dog as an outgroup and identified 51 intron losses in mouse and 3 losses in human and no gain. The rate of intron loss in Tetraodon is higher than that in fugu, mouse, and human but lower than the previous estimates for other eukaryotes. The introns lost in pufferfishes and mammals are significantly shorter than the mean size of introns in the genome. One intron deleted in fugu and another in Tetraodon have left behind 6 and 3 nucleotides, respectively, suggesting that they were lost due to genomic deletions. Such losses of introns are likely to be the result of a higher rate of DNA deletions experienced by the genomes of pufferfishes compared with mammals. The shorter generation time of Tetraodon compared with fugu, and the rich diversity and higher activity of transposable elements in pufferfishes compared with mammals, may be responsible for the higher rate of intron loss in Tetraodon. Our findings indicate that overall very little intron turnover has occurred in pufferfishes and mammals during recent evolution and that intron gain is an extremely rare event in vertebrate evolution.     

Intron loss and gain during evolution of the catalase gene family in angiosperms.

Angiosperms (flowering plants), including both monocots and dicots, contain small catalase gene families. In the dicot, Arabidopsis thaliana, two catalase (CAT) genes, CAT1 and CAT3, are tightly linked on chromosome 1 and a third, CAT2, which is more similar to CAT1 than to CAT3, is unlinked on chromosome 4. Comparison of positions and numbers of introns among 13 angiosperm catalase genomic sequences indicates that intron positions are conserved, and suggests that an ancestral catalase gene common to monocots and dicots contained seven introns. Arabidopsis CAT2 has seven introns; both CAT1 and CAT3 have six introns in positions conserved with CAT2, but each has lost a different intron. We suggest the following sequence of events during the evolution of the Arabidopsis catalase gene family. An initial duplication of an ancestral catalase gene gave rise to CAT3 and CAT1. CAT1 then served as the template for a second duplication, yielding CAT2. Intron losses from CAT1 and CAT3 followed these duplications. One subclade of monocot catalases has lost all but the 5''-most and 3''-most introns, which is consistent with a mechanism of intron loss by replacement of an ancestral intron-containing gene with a reverse-transcribed DNA copy of a fully spliced mRNA. Following this event of concerted intron loss, the Oryza sativa (rice, a monocot) CAT1 lineage acquired an intron in a novel position, consistent with a mechanism of intron gain at proto-splice sites.     

High rate of recent intron gain and loss in simultaneously duplicated Arabidopsis genes

We examined the gene structure of a set of 2563 Arabidopsis thaliana paralogous pairs that were duplicated simultaneously 20-60 MYA by tetraploidy. Out of a total of 23,164 introns in these genes, we found that 10,004 pairs have been conserved and 578 introns have been inserted or deleted in the time since the duplication event. This intron insertion/deletion rate of 2.7 x 10(-3) to 9.1 x 10(-4) per site per million years is high in comparison to previous studies. At least 56 introns were gained and 39 lost based on parsimony analysis of the phylogenetic distribution of these introns. We found weak evidence that genes undergoing intron gain and loss are biased with respect to gene ontology terms. Gene pairs that experienced at least 2 intron insertions or deletions show evidence of enrichment for membrane location and transport and transporter activity function. We do not find any relationship of intron flux to expression level or G + C content of the gene. Detection of a bias in the location of intron gains and losses within a gene depends on the method of measurement: an intragene method indicates that events (specifically intron losses) are biased toward the 3' end of the gene. Despite the relatively recent acquisition of these introns, we found only one case where we could identify the mechanism of intron origin--the TOUCH3 gene has experienced 2 tandem, partial, internal gene duplications that duplicated a preexisting intron and also created a novel, alternatively spliced intron that makes use of a duplicated pair of cryptic splice sites.     

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.     

Higher intron loss rate in Arabidopsis thaliana than A. lyrata is consistent with stronger selection for a smaller genome

The number of introns varies considerably among different organisms. This can be explained by the differences in the rates of intron gain and loss. Two factors that are likely to influence these rates are selection for or against introns and the mutation rate that generates the novel intron or the intronless copy. Although it has been speculated that stronger selection for a compact genome might result in a higher rate of intron loss and a lower rate of intron gain, clear evidence is lacking, and the role of selection in determining these rates has not been established. Here, we studied the gain and loss of introns in the two closely related species Arabidopsis thaliana and A. lyrata as it was recently shown that A. thaliana has been undergoing a faster genome reduction driven by selection. We found that A. thaliana has lost six times more introns than A. lyrata since the divergence of the two species but gained very few introns. We suggest that stronger selection for genome reduction probably resulted in the much higher intron loss rate in A. thaliana, although further analysis is required as we could not find evidence that the loss rate increased in A. thaliana as opposed to having decreased in A. lyrata compared with the rate in the common ancestor. We also examined the pattern of the intron gains and losses to better understand the mechanisms by which they occur. Microsimilarity was detected between the splice sites of several gained and lost introns, suggesting that nonhomologous end joining repair of double-strand breaks might be a common pathway not only for intron gain but also for intron loss.     

Tempo and Mode of Spliceosomal Intron Evolution in Actin of Foraminifera

Spliceosomal introns are present in almost all eukaryotic genes, yet little is known about their origin and turnover in the majority of eukaryotic phyla. There is no agreement whether most introns are ancestral and have been lost in some lineage or have been gained recently. We addressed this question by analyzing the spatial and temporal distribution of introns in actins of foraminifera, a group of testate protists whose exceptionally rich fossil record permits the calibration of molecular phylogenies to date intron origins. We identified 24 introns dispersed along the sequence of two foraminiferan actin paralogues and actin deviating proteins, an unconventional type of fast-evolving actin found in some foraminifera. Comparison of intron positions indicates that 20 of 24 introns are specific to foraminifera. Four introns shared between foraminifera and other eukaryotes were interpreted as parallel gains because they have been found only in single species belonging to phylogenetically distinctive lineages. Moreover, additional recent intron gain due to the transfer between the actin paralogues was observed in two cultured species. Based on a relaxed molecular clock timescale, we conclude that intron gains in actin took place throughout the evolution of foraminifera, with the oldest introns inserted between 550 and 500 million years ago and the youngest ones acquired less than 100 million years ago. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Debashish Bhattacharya]     

Analysis of ribosomal protein gene structures: implications for intron evolution

Many spliceosomal introns exist in the eukaryotic nuclear genome. Despite much research, the evolution of spliceosomal introns remains poorly understood. In this paper, we tried to gain insights into intron evolution from a novel perspective by comparing the gene structures of cytoplasmic ribosomal proteins (CRPs) and mitochondrial ribosomal proteins (MRPs), which are held to be of archaeal and bacterial origin, respectively. We analyzed 25 homologous pairs of CRP and MRP genes that together had a total of 527 intron positions. We found that all 12 of the intron positions shared by CRP and MRP genes resulted from parallel intron gains and none could be considered to be “conserved,” i.e., descendants of the same ancestor. This was supported further by the high frequency of proto-splice sites at these shared positions; proto-splice sites are proposed to be sites for intron insertion. Although we could not definitively disprove that spliceosomal introns were already present in the last universal common ancestor, our results lend more support to the idea that introns were gained late. At least, our results show that MRP genes were intronless at the time of endosymbiosis. The parallel intron gains between CRP and MRP genes accounted for 2.3% of total intron positions, which should provide a reliable estimate for future inferences of intron evolution.     

Smoke without fire: most reported cases of intron gain in nematodes instead reflect intron losses

Identification of recently gained spliceosomal introns would provide crucial evidence in the continuing debate concerning the age and evolutionary significance of introns. A previously published genomic analysis reported to have identified 122 introns that had been gained since the divergence of the nematodes Caenorhabidits elegans and Caenorhabditis briggsae approximately 100 MYA. However, using newly available genomic sequence from additional Caenorhabditis species, we show that 74% (60/81) of the reported gains in C. elegans are present in a C. briggsae relative. This pattern indicates that these introns represent losses in C. briggsae, not gains in C. elegans. In addition, 61% (25/41) of the reported gains in C. briggsae are present in the more distant C. briggsae relative, in a pattern suggesting that additional reported gains in C. elegans and/or C. briggsae may in fact represent unrecognized losses. These results underscore the dominance of intron loss over intron gain in recent eukaryotic evolution, the pitfalls associated with parsimony in inferring intron gains, and the importance of genomic sequencing of clusters of closely related species for drawing accurate inferences about genome evolution.     

Evidence for extensive recent intron transposition in closely related fungi

Though spliceosomal introns are a major structural component of most eukaryotic genes and intron density varies by more than three orders of magnitude among eukaryotes [1-3], the origins of introns are poorly understood, and only a few cases of unambiguous intron gain are known [4-8]. We utilized population genomic comparisons of three closely related fungi to identify crucial transitory phases of intron gain and loss. We found 74 intron positions showing intraspecific presence-absence polymorphisms (PAPs) for the entire intron. Population genetic analyses identified intron PAPs at different stages of fixation and showed that intron gain or loss was very recent. We found direct support for extensive intron transposition among unrelated genes. A substantial proportion of highly similar introns in the genome either were recently gained or showed a transient phase of intron PAP. We also identified an intron transfer among paralogous genes that created a new intron. Intron loss was due mainly to homologous recombination involving reverse-transcribed mRNA. The large number of intron positions in transient phases of either intron gain or loss shows that intron evolution is much faster than previously thought and provides an excellent model to study molecular mechanisms of intron gain.     

Piecing together the puzzle of parasitic plant plastome evolution

The importance of photosynthesis as a mode of energy production has put plastid genomes of plants under a constant purifying selection. This has shaped the characteristic features of plastid genomes across the entire spectrum of photosynthetic plants and has led to a highly uniform and conserved plastid genome with respect to structure, size, gene order, intron and editing site positions and coding capacity. Parasitic species that have dropped photosynthesis as the main energy provider share striking deviations from the plastid genome norm: multiple rearrangements within the circular chromosome, pseudogenization and gene deletions, promoter losses, intron losses as well as the extensive loss of mRNA editing competence have been reported. The collective loss of larger sets of functionally related genes like those for the plastid NADH–dehydrogenase complex and concomitant losses of RNA polymerase genes together with their target promoters point to “domino effects” where an initial loss might have triggered others. An example, which will be discussed in more detail, is the concomitant loss of the intron maturase gene matK and all introns that are supposedly subject to MatK-dependent splicing in two Cuscuta species.