A high density of ancient spliceosomal introns in oxymonad excavates

The Oxymonad Genome Displays Canonical Eukaryotic Complexity in the Absence of a Mitochondrion

The discovery that the protist Monocercomonoides exilis completely lacks mitochondria demonstrates that these organelles are not absolutely essential to eukaryotic cells. However, the degree to which the metabolism and cellular systems of this organism have adapted to the loss of mitochondria is unknown. Here, we report an extensive analysis of the M. exilis genome to address this question. Unexpectedly, we find that M. exilis genome structure and content is similar in complexity to other eukaryotes and less "reduced" than genomes of some other protists from the Metamonada group to which it belongs. Furthermore, the predicted cytoskeletal systems, the organization of endomembrane systems, and biosynthetic pathways also display canonical eukaryotic complexity. The only apparent preadaptation that permitted the loss of mitochondria was the acquisition of the SUF system for Fe-S cluster assembly and the loss of glycine cleavage system. Changes in other systems, including in amino acid metabolism and oxidative stress response, were coincident with the loss of mitochondria but are likely adaptations to the microaerophilic and endobiotic niche rather than the mitochondrial loss per se. Apart from the lack of mitochondria and peroxisomes, we show that M. exilis is a fully elaborated eukaryotic cell that is a promising model system in which eukaryotic cell biology can be investigated in the absence of mitochondria.

Reverse transcriptase and intron number evolution

BACKGROUND

Introns are universal in eukaryotic genomes and play important roles in transcriptional regulation, mRNA export to the cytoplasm, nonsense-mediated decay as both a regulatory and a splicing quality control mechanism, R-loop avoidance, alternative splicing, chromatin structure, and evolution by exon-shuffling.

METHODS

Sixteen complete fungal genomes were used 13 of which were sequenced and annotated by JGI. Ustilago maydis, Cryptococcus neoformans, and Coprinus cinereus (also named Coprinopsis cinerea) were from the Broad Institute. Gene models from JGI-annotated genomes were taken from the GeneCatalog track that contained the best representative gene models. Varying fractions of the GeneCatalog were manually curated by external users. For clarity, we used the JGI unique database identifier.

RESULTS

The last common ancestor of eukaryotes (LECA) has an estimated 6.4 coding exons per gene (EPG) and evolved into the diverse eukaryotic life forms, which is recapitulated by the development of a stem cell. We found a parallel between the simulated reverse transcriptase (RT)-mediated intron loss and the comparative analysis of 16 fungal genomes that spanned a wide range of intron density. Although footprints of RT (RTF) were dynamic, relative intron location (RIL) to the 5'-end of mRNA faithfully traced RT-mediated intron loss and revealed 7.7 EPG for LECA. The mode of exon length distribution was conserved in simulated intron loss, which was exemplified by the shared mode of 75 nt between fungal and Chlamydomonas genomes. The dominant ancient exon length was corroborated by the average exon length of the most intron-rich genes in fungal genomes and consistent with ancient protein modules being ~25 aa. Combined with the conservation of a protein length of 400 aa, the earliest ancestor of eukaryotes could have 16 EPG. During earlier evolution, Ascomycota's ancestor had significantly more 3'-biased RT-mediated intron loss that was followed by dramatic RTF loss. There was a down trend of EPG from more conserved to less conserved genes. Moreover, species-specific genes have higher exon-densities, shorter exons, and longer introns when compared to genes conserved at the phylum level. However, intron length in species-specific genes became shorter than that of genes conserved in all species after genomes experiencing drastic intron loss. The estimated EPG from the most frequent exon length is more than double that from the RIL method.

CONCLUSIONS

This implies significant intron loss during the very early period of eukaryotic evolution. De novo gene-birth contributes to shorter exons, longer introns, and higher exon-density in species-specific genes relative to conserved genes.

Characterization of newly gained introns in Daphnia populations

As one of the few known species in an active phase of intron proliferation, the microcrustacean Daphnia pulex is an especially attractive system for interrogating the gain and loss of introns in natural populations. In this study, we used a comparative population-genomic approach to identify and characterize 90 recently gained introns in this species. Molecular clock analyses indicate that these introns arose between 3.9 × 10(5) and 1.45 × 10(4) years ago, with a spike in intron proliferation approximately 5.2 × 10(4) to 1.22 × 10(5) years ago. Parallel gains at homologous positions contribute to 47.8% (43/90) of discovered new introns. A disproportionally large number of new introns were found in historically isolated populations in Oregon. Nonetheless, derived, intron-bearing alleles were also identified in a wide range of geographic locations, suggesting intron gain and, to a lesser degree, intron loss are important sources of genetic variation in natural populations of Daphnia. A majority (55/90 or 61.1%) of the identified neointrons have associated internal direct repeats with lengths and compositions that are unlikely to occur by chance, suggesting repeated bouts of staggered double-strand breaks (DSBs) during their evolution. Accordingly, internal, staggered DSBs may contribute to a passive trend toward increased length and sequence diversity in nascent introns.

Single-Cell DNA barcoding using sequences from the small subunit rRNA and internal transcribed spacer region identifies new species of Trichonympha and Trichomitopsis from the hindgut of the termite Zootermopsis angusticollis

To aid in their digestion of wood, lower termites are known to harbour a diverse community of prokaryotes as well as parabasalid and oxymonad protist symbionts. One of the best-studied lower termite gut communities is that of Zootermopsis angusticollis which has been known for almost 100 years to possess 3 species of Trichonympha (T. campanula, T. collaris, and T. sphaerica), 1 species of Trichomitopsis (T. termopsidis), as well as smaller flagellates. We have re-assessed this community by sequencing the small subunit (SSU) rRNA gene and the internal transcribed spacer (ITS) region from a large number of single Trichonympha and Trichomitopsis cells for which morphology was also documented. Based on phylogenetic clustering and sequence divergence, we identify 3 new species: Trichonympha postcylindrica, Trichomitopsis minor, and Trichomitopsis parvus spp. nov. Once identified by sequencing, the morphology of the isolated cells for all 3 new species was re-examined and found to be distinct from the previously described species: Trichonympha postcylindrica can be morphologically distinguished from the other Trichonympha species by an extension on its posterior end, whereas Trichomitopsis minor and T. parvus are smaller than T. termopsidis but similar in size to each other and cannot be distinguished based on morphology using light microscopy. Given that Z. angusticollis has one of the best characterized hindgut communities, the near doubling of the number of the largest and most easily identifiable symbiont species suggests that the diversity of hindgut symbionts is substantially underestimated in other termites as well. Accurate descriptions of the diversity of these microbial communities are essential for understanding hindgut ecology and disentangling the interactions among the symbionts, and molecular barcoding should be a priority for these systems.

Phylogenomic evidence for the presence of a flagellum and cbb(3) oxidase in the free-living mitochondrial ancestor

The initiation of the intracellular symbiosis that would give rise to mitochondria and eukaryotes was a major event in the history of life on earth. Hypotheses to explain eukaryogenesis fall into two broad and competing categories: those proposing that the host was a phagocytotic proto-eukaryote that preyed upon the free-living mitochondrial ancestor (hereafter FMA), and those proposing that the host was an archaebacterium that engaged in syntrophy with the FMA. Of key importance to these hypotheses are whether the FMA was motile or nonmotile, and the atmospheric conditions under which the FMA thrived. Reconstructions of the FMA based on genome content of Rickettsiales representatives-generally considered to be the closest living relatives of mitochondria-indicate that it was nonmotile and aerobic. We have sequenced the genome of Candidatus Midichloria mitochondrii, a novel and phylogenetically divergent member of the Rickettsiales. We found that it possesses unique gene sets found in no other Rickettsiales, including 26 genes associated with flagellar assembly, and a cbb(3)-type cytochrome oxidase. Phylogenomic analyses show that these genes were inherited in a vertical fashion from an ancestral α-proteobacterium, and indicate that the FMA possessed a flagellum, and could undergo oxidative phosphorylation under both aerobic and microoxic conditions. These results indicate that the FMA played a more active and potentially parasitic role in eukaryogenesis than currently appreciated and provide an explanation for how the symbiosis could have evolved under low levels of oxygen.

The Naegleria genome: a free-living microbial eukaryote lends unique insights into core eukaryotic cell biology

Naegleria gruberi, a free-living protist, has long been treasured as a model for basal body and flagellar assembly due to its ability to differentiate from crawling amoebae into swimming flagellates. The full genome sequence of Naegleria gruberi has recently been used to estimate gene families ancestral to all eukaryotes and to identify novel aspects of Naegleria biology, including likely facultative anaerobic metabolism, extensive signaling cascades, and evidence for sexuality. Distinctive features of the Naegleria genome and nuclear biology provide unique perspectives for comparative cell biology, including cell division, RNA processing and nucleolar assembly. We highlight here exciting new and novel aspects of Naegleria biology identified through genomic analysis.

Broadly sampled multigene analyses yield a well-resolved eukaryotic tree of life

An accurate reconstruction of the eukaryotic tree of life is essential to identify the innovations underlying the diversity of microbial and macroscopic (e.g., plants and animals) eukaryotes. Previous work has divided eukaryotic diversity into a small number of high-level "supergroups," many of which receive strong support in phylogenomic analyses. However, the abundance of data in phylogenomic analyses can lead to highly supported but incorrect relationships due to systematic phylogenetic error. Furthermore, the paucity of major eukaryotic lineages (19 or fewer) included in these genomic studies may exaggerate systematic error and reduce power to evaluate hypotheses. Here, we use a taxon-rich strategy to assess eukaryotic relationships. We show that analyses emphasizing broad taxonomic sampling (up to 451 taxa representing 72 major lineages) combined with a moderate number of genes yield a well-resolved eukaryotic tree of life. The consistency across analyses with varying numbers of taxa (88-451) and levels of missing data (17-69%) supports the accuracy of the resulting topologies. The resulting stable topology emerges without the removal of rapidly evolving genes or taxa, a practice common to phylogenomic analyses. Several major groups are stable and strongly supported in these analyses (e.g., SAR, Rhizaria, Excavata), whereas the proposed supergroup "Chromalveolata" is rejected. Furthermore, extensive instability among photosynthetic lineages suggests the presence of systematic biases including endosymbiotic gene transfer from symbiont (nucleus or plastid) to host. Our analyses demonstrate that stable topologies of ancient evolutionary relationships can be achieved with broad taxonomic sampling and a moderate number of genes. Finally, taxon-rich analyses such as presented here provide a method for testing the accuracy of relationships that receive high bootstrap support (BS) in phylogenomic analyses and enable placement of the multitude of lineages that lack genome scale data.

The genome of Naegleria gruberi illuminates early eukaryotic versatility

Genome sequences of diverse free-living protists are essential for understanding eukaryotic evolution and molecular and cell biology. The free-living amoeboflagellate Naegleria gruberi belongs to a varied and ubiquitous protist clade (Heterolobosea) that diverged from other eukaryotic lineages over a billion years ago. Analysis of the 15,727 protein-coding genes encoded by Naegleria's 41 Mb nuclear genome indicates a capacity for both aerobic respiration and anaerobic metabolism with concomitant hydrogen production, with fundamental implications for the evolution of organelle metabolism. The Naegleria genome facilitates substantially broader phylogenomic comparisons of free-living eukaryotes than previously possible, allowing us to identify thousands of genes likely present in the pan-eukaryotic ancestor, with 40% likely eukaryotic inventions. Moreover, we construct a comprehensive catalog of amoeboid-motility genes. The Naegleria genome, analyzed in the context of other protists, reveals a remarkably complex ancestral eukaryote with a rich repertoire of cytoskeletal, sexual, signaling, and metabolic modules.

Mystery of intron gain: new data and new models

Despite their ubiquity, the mechanisms and evolutionary forces responsible for the origins of spliceosomal introns remain mysterious. Recent molecular evidence supports the idea that intronic RNAs can reverse splice into RNA transcripts, a crucial step for an influential model of intron gain. However, a paradox attends this model because the rate of intron gain is expected to be orders of magnitude lower than the rate of intron loss in general, in contrast to findings from several lineages. We suggest two possible resolutions to this paradox, based on steric considerations and on the possibility of co-option by specific introns of retroelement transposition pathways, respectively. In addition, we introduce two potential mechanisms for intron creation, based on hybrid RNA-DNA reverse splicing and on template switching errors by reverse transcriptase.

Evolution of GHF5 endoglucanase gene structure in plant-parasitic nematodes: no evidence for an early domain shuffling event

Background

Endo-1,4-beta-glucanases or cellulases from the glycosyl hydrolase family 5 (GHF5) have been found in numerous bacteria and fungi, and recently also in higher eukaryotes, particularly in plant-parasitic nematodes (PPN). The origin of these genes has been attributed to horizontal gene transfer from bacteria, although there still is a lot of uncertainty about the origin and structure of the ancestral GHF5 PPN endoglucanase. It is not clear whether this ancestral endoglucanase consisted of the whole gene cassette, containing a catalytic domain and a carbohydrate-binding module (CBM, type 2 in PPN and bacteria) or only of the catalytic domain while the CBM2 was retrieved by domain shuffling later in evolution. Previous studies on the evolution of these genes have focused primarily on data of sedentary nematodes, while in this study, extra data from migratory nematodes were included.

Results

Two new endoglucanases from the migratory nematodes Pratylenchus coffeae and Ditylenchus africanus were included in this study. The latter one is the first gene isolated from a PPN of a different superfamily (Sphaerularioidea); all previously known nematode endoglucanases belong to the superfamily Tylenchoidea (order Rhabditida). Phylogenetic analyses were conducted with the PPN GHF5 endoglucanases and homologous endoglucanases from bacterial and other eukaryotic lineages such as beetles, fungi and plants. No statistical incongruence between the phylogenetic trees deduced from the catalytic domain and the CBM2 was found, which could suggest that both domains have evolved together. Furthermore, based on gene structure data, we inferred a model for the evolution of the GHF5 endoglucanase gene structure in plant-parasitic nematodes. Our data confirm a close relationship between Pratylenchus spp. and the root knot nematodes, while some Radopholus similis endoglucanases are more similar to cyst nematode genes.

Conclusion

We conclude that the ancestral PPN GHF5 endoglucanase gene most probably consisted of the whole gene cassette, i.e. the GHF5 catalytic domain and the CBM2, rather than that it evolved by domain shuffling. Our evolutionary model for the gene structure in PPN GHF5 endoglucanases implies the occurrence of an early duplication event, and more recent gene duplications at genus or species level.

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.

Environmental PCR survey to determine the distribution of a non-canonical genetic code in uncultivable oxymonads

The universal genetic code is conserved throughout most living systems, but a non-canonical code where TAA and TAG encode glutamine has evolved in several eukaryotes, including oxymonad protists. Most oxymonads are uncultivable, so environmental RT-PCR and PCR was used to examine the distribution of this rare character. A total of 253 unique isolates of four protein-coding genes were sampled from the hindgut community of the cockroach, Cryptocercus punctulatus, an environment rich in diversity from two of the five subgroups of oxymonad, saccinobaculids and polymastigids. Four alpha-tubulins were found with non-canonical glutamine codons. Environmental RACE confirmed that these and related genes used only TGA as stop codons, as expected for the non-canonical code, whereas other genes used TAA or TAG as stop codons, as expected for the universal code. We characterized alpha-tubulin from manually isolated Saccinobaculus ambloaxostylus, confirming it uses the universal code and suggesting, by elimination, that the non-canonical code is used by a polymastigid. HSP90 and EF-1alpha phylogenies also showed environmental sequences falling into two distinct groups, and are generally consistent with previous hypotheses that polymastigids and Streblomastix are closely related. Overall, we propose that the non-canonical genetic code arose once in a common ancestor of Streblomastix and a subgroup of polymastigids.

Intron mis-splicing: no alternative?

Eukaryotes compensate for inefficient splicing by mechanisms that prevent the translation of mis-spliced or unspliced mRNAs.

A recent report reveals widespread mis-splicing of RNA transcripts in eukaryotes, with mis-spliced RNA destroyed by nonsense-mediated mRNA decay. This striking inefficiency deepens the mystery of the proliferation and persistence of introns.

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.

Functional and evolutionary analysis of alternatively spliced genes is consistent with an early eukaryotic origin of alternative splicing

Background

Alternative splicing has been reported in various eukaryotic groups including plants, apicomplexans, diatoms, amoebae, animals and fungi. However, whether widespread alternative splicing has evolved independently in the different eukaryotic groups or was inherited from their last common ancestor, and may therefore predate multicellularity, is still unknown. To better understand the origin and evolution of alternative splicing and its usage in diverse organisms, we studied alternative splicing in 12 eukaryotic species, comparing rates of alternative splicing across genes of different functional classes, cellular locations, intron/exon structures and evolutionary origins.

Results

For each species, we find that genes from most functional categories are alternatively spliced. Ancient genes (shared between animals, fungi and plants) show high levels of alternative splicing. Genes with products expressed in the nucleus or plasma membrane are generally more alternatively spliced while those expressed in extracellular location show less alternative splicing. We find a clear correspondence between incidence of alternative splicing and intron number per gene both within and between genomes. In general, we find several similarities in patterns of alternative splicing across these diverse eukaryotes.

Conclusion

Along with previous studies indicating intron-rich genes with weak intron boundary consensus and complex spliceosomes in ancestral organisms, our results suggest that at least a simple form of alternative splicing may already have been present in the unicellular ancestor of plants, fungi and animals. A role for alternative splicing in the evolution of multicellularity then would largely have arisen by co-opting the preexisting process.

Intron-rich gene structure in the intracellular plant parasite Plasmodiophora brassicae

Plasmodiophora brassicae, a pathogen of Brassicaceae plants, is grouped within the eukaryotic supergroup, the Rhizaria. Although a large diversity of protists is found in the Rhizaria, genomes of organisms within the group have barely been examined. In this study, we identified DNA sequences spanning or flanking 24 P. brassicae genes, eventually sequencing close to 44 kb of genomic DNA. Evidence from this preliminary genome survey suggested that splicing is an important feature of P. brassicae gene expression; the P. brassicae genes were rich in spliceosomal introns and two mini-exons of less than 20 bp were identified. Consensus splice sites and branch-point sequences in P. brassicae introns were similar to those found in other eukaryotes. Examination of the promoter and transcription start sites of genes indicated that P. brassicae transcription is likely to begin from initiator elements rather than TATA-box containing promoters. Where neighbouring genes were confirmed, intergenic distances were short, ranging from 44 to 470 bp, but a number of larger DNA fragments containing no obvious genes were also sequenced.

Protecting exons from deleterious R-loops: a potential advantage of having introns

Background

Accumulating evidence indicates that the nascent RNA can invade and pair with one strand of DNA, forming an R-loop structure that threatens the stability of the genome. In addition, the cost and benefit of introns are still in debate.

Results

At least three factors are likely required for the R-loop formation: 1) sequence complementarity between the nascent RNA and the target DNA, 2) spatial juxtaposition between the nascent RNA and the template DNA, and 3) accessibility of the template DNA and the nascent RNA. The removal of introns from pre-mRNA reduces the complementarity between RNA and the template DNA and avoids the spatial juxtaposition between the nascent RNA and the template DNA. In addition, the secondary structures of group I and group II introns may act as spatial obstacles for the formation of R-loops between nearby exons and the genomic DNA.

Conclusion

Organisms may benefit from introns by avoiding deleterious R-loops. The potential contribution of this benefit in driving intron evolution is discussed. I propose that additional RNA polymerases may inhibit R-loop formation between preceding nascent RNA and the template DNA. This idea leads to a testable prediction: intermittently transcribed genes and genes with frequently prolonged transcription should have higher intron density.

Reviewers

This article was reviewed by Dr. Eugene V. Koonin, Dr. Alexei Fedorov (nominated by Dr. Laura F Landweber), and Dr. Scott W. Roy (nominated by Dr. Arcady Mushegian).

Coevolution of genomic intron number and splice sites

Spliceosomal intron numbers and boundary sequences vary dramatically in eukaryotes. We found a striking correspondence between low intron number and strong sequence conservation of 5' splice sites (5'ss) across eukaryotic genomes. The phylogenetic pattern suggests that ancestral 5'ss were relatively weakly conserved, but that some lineages independently underwent both major intron loss and 5'ss strengthening. It seems that eukaryotic ancestors had relatively large intron numbers and 'weak' 5'ss, a pattern associated with frequent alternative splicing in modern organisms.

Evolutionary conservation of UTR intron boundaries in Cryptococcus

Despite significant progress, the general functional and evolutionary significance of the untranslated regions (UTRs) of eukaryotic transcripts remain mysterious. Particularly mysterious is the common occurrence of spliceosomal introns in transcript UTRs because UTR splicing is not necessary for restoration of transcript coding sequence. In general, it is not known to what extent such splicing performs an important function or merely represents spliceosomal "noise." We conducted the first analysis of evolutionary conservation of UTR splicing. Among 4 species from Cryptococcus neoformans species complex, we find high levels of conservation of UTR intron boundary sequences, strongly suggesting that UTR intron splicing is conserved by purifying selection. We estimate that 50-90% of splice boundaries are maintained by selection. Donor site sequences are more highly conserved than acceptor sequences, and splicing boundaries are more conserved in 5' UTRs than in 3' UTRs. In addition, we report a variety of differences between patterns of UTR splicing in Cryptococcus and corresponding patterns in animals and plants. These results focus attention on the functional roles of eukaryotic UTRs and deepen the mystery of UTR intron splicing.

Conserved intron positions in ancient protein modules

Background

The timing of the origin of introns is of crucial importance for an understanding of early genome architecture. The Exon theory of genes proposed a role for introns in the formation of multi-exon proteins by exon shuffling and predicts the presence of conserved splice sites in ancient genes. In this study, large-scale analysis of potential conserved splice sites was performed using an intron-exon database (ExInt) derived from GenBank.

Results

A set of conserved intron positions was found by matching identical splice sites sequences from distantly-related eukaryotic kingdoms. Most amino acid sequences with conserved introns were homologous to consensus sequences of functional domains from conserved proteins including kinases, phosphatases, small GTPases, transporters and matrix proteins. These included ancient proteins that originated before the eukaryote-prokaryote split, for instance the catalytic domain of protein phosphatase 2A where a total of eleven conserved introns were found. Using an experimental setup in which the relation between a splice site and the ancientness of its surrounding sequence could be studied, it was found that the presence of an intron was positively correlated to the ancientness of its surrounding sequence. Intron phase conservation was linked to the conservation of the gene sequence and not to the splice site sequence itself. However, no apparent differences in phase distribution were found between introns in conserved versus non-conserved sequences.

Conclusion

The data confirm an origin of introns deep in the eukaryotic branch and is in concordance with the presence of introns in the first functional protein modules in an 'Exon theory of genes' scenario. A model is proposed in which shuffling of primordial short exonic sequences led to the formation of the first functional protein modules, in line with hypotheses that see the formation of introns integral to the origins of genome evolution.

Reviewers

This article was reviewed by Scott Roy (nominated by Anthony Poole), Sandro de Souza (nominated by Manyuan Long), and Gáspár Jékely.

Origins and Evolution of Spliceosomal Introns

Research into the origins of introns is at a critical juncture in the resolution of theories on the evolution of early life (which came first, RNA or DNA?), the identity of LUCA (the last universal common ancestor, was it prokaryotic- or eukaryotic-like?), and the significance of noncoding nucleotide variation. One early notion was that introns would have evolved as a component of an efficient mechanism for the origin of genes. But alternative theories emerged as well. From the debate between the "introns-early" and "introns-late" theories came the proposal that introns arose before the origin of genetically encoded proteins and DNA, and the more recent "introns-first" theory, which postulates the presence of introns at that early evolutionary stage from a reconstruction of the "RNA world." Here we review seminal and recent ideas about intron origins. Recent discoveries about the patterns and causes of intron evolution make this one of the most hotly debated and exciting topics in molecular evolutionary biology today.