Intron size, abundance, and distribution within untranslated regions of genes

Most research concerning the evolution of introns has largely considered introns within coding sequences (CDSs), without regard for introns located within untranslated regions (UTRs) of genes. Here, we directly determined intron size, abundance, and distribution in UTRs of genes using full-length cDNA libraries and complete genome sequences for four species, Arabidopsis thaliana, Drosophila melanogaster, human, and mouse. Overall intron occupancy (introns/exon kbp) is lower in 5' UTRs than CDSs, but intron density (intron occupancy in regions containing introns) tends to be higher in 5' UTRs than in CDSs. Introns in 5' UTRs are roughly twice as large as introns in CDSs, and there is a sharp drop in intron size at the 5' UTR-CDS boundary. We propose a mechanistic explanation for the existence of selection for larger intron size in 5' UTRs, and outline several implications of this hypothesis. We found introns to be randomly distributed within 5' UTRs, so long as a minimum required exon size was assumed. Introns in 3' UTRs were much less abundant than in 5' UTRs. Though this was expected for human and mouse that have intron-dependent nonsense-mediated decay (NMD) pathways that discourage the presence of introns within the 3' UTR, it was also true for A. thaliana and D. melanogaster, which may lack intron-dependent NMD. Our findings have several implications for theories of intron evolution and genome evolution in general.     

Intron loss and gain in Drosophila

Although introns were first discovered almost 30 years ago, their evolutionary origin remains elusive. In this work, we used multispecies whole-genome alignments to map Drosophila melanogaster introns onto 10 other fully sequenced Drosophila genomes. We were able to find 1,944 sites where an intron was missing in one or more species. We show that for most (>80%) of these cases, there is no leftover intronic sequence or any missing exonic sequence, indicating exact intron loss or gain events. We used parsimony to classify these differences as 1,754 intron loss events and 213 gain events. We show that lost and gained introns are significantly shorter than average and flanked by longer than average exons. They also display quite distinct phase distributions and show greater than average similarity between the 5' splice site and its 3' partner splice site. Introns that have been lost in one or more species evolve faster than other introns, occur in slowly evolving genes, and are found adjacent to each other more often than would be expected for independent single losses. Our results support the cDNA recombination mechanism of intron loss, suggest that selective pressures affect site-specific loss rates, and show conclusively that intron gain has occurred within the Drosophila lineage, solidifying the "introns-middle" hypothesis and providing some hints about the gain mechanism.     

The intronome of budding yeasts

Whatever their abundance in genomes, spliceosomal introns are the signature of eukaryotic genes. The sequence of Saccharomyces cerevisiae, achieved fifteen years ago, revealed that this yeast has very few introns, but conserved intron boundaries typical for an intron definition mechanism. With the improvement and the development of new sequencing technologies, yeast genomes have been extensively sequenced during the last decade. We took advantage of this plethora of data to compile and assess the intron content of the protein-coding genes of 13 genomes representative of the evolution of hemiascomycetous yeasts. We first observed that intron paucity is a general rule and that the fastest evolving genomes tend to lose their introns more rapidly (e.g. S. cerevisiae versus Yarrowia lipolytica). Noticeable differences were also confirmed for 5' splice sites and branch point sites (BP) as well as for the relative position of the BP. These changes seemed to be correlated with the lineage specific evolution of splicing factors.     

Nuclear pre-mRNA processing in plants: distinct modes of 3''-splice-site selection in plants and animals.

The report that human growth hormone pre-mRNA is not processed in transgenic plant tissues (A. Barta, K. Sommergruber, D. Thompson, K. Hartmuth, M.A. Matzke, and A.J.M. Matzke, Plant Mol. Biol. 6:347-357, 1986) has suggested that differences in mRNA splicing processes exist between plants and animals. To gain more information about the specificity of plant pre-mRNA processing, we have compared the splicing of the soybean leghemoglobin pre-mRNA with that of the human beta-globin pre-mRNA in transfected plant (Orychophragmus violaceus and Nicotiana tabacum) protoplasts and mammalian (HeLa) cells. Of the three introns of leghemoglobin pre-mRNA, only intron 2 was correctly and efficiently processed in HeLa cells. The 5' splice sites of the remaining two introns were faithfully recognized, but correct processing of the 3' sites took place only rarely (intron 1) or not at all (intron 3); cryptic 3' splice sites were used instead. While the first intron in human beta-globin pre-mRNA was not spliced in transfected plant protoplasts, intron 2 processing occurred at a low level, indicating that some mammalian introns can be recognized by the plant intron-splicing machinery. However, excision of intron 2 proved to be incorrect, involving the authentic 5' splice site and a cryptic 3' splice site. Our results indicate that the mechanism of 3'-splice-site selection during intron excision differs between plants and animals. This conclusion is supported by analysis of the 3'-splice-site consensus sequences in animal and plant introns which revealed that polypyrimidine tracts, characteristic of animal introns, are not present in plant pre-mRNAs. It is proposed that an elevated AU content of plant introns is important for their processing.     

Intron splicing: a conserved internal signal in introns of Drosophila pre-mRNAs.

The introns of Drosophila pre-mRNAs have been analysed for conserved internal sequence elements near the 3' intron boundary similar to the T-A-C-T-A-A-C in yeast introns and the C/T-T-A/G-A-C/T in introns of other organisms. Such conserved internal elements are the 3' splice signals recognized in intron splicing. In the lariat splicing mechanism, the G at the 5' end of an intron joins covalently to the last A of a 3' splice signal to form a branch point in a splicing intermediate. Analysis of 39 published sequences of Drosophila introns reveals that potential 3' splice signals with the consensus C/T-T-A/G-A-C/T are present in 18 cases. In 17 of the remaining cases signals are present which vary from this consensus just in the middle or last position. In Drosophila introns the 3' splice signal is usually located in a discrete region between 18 and 35 nucleotides upstream from the 3' splice point. We note that the Drosophila small nuclear U2-RNA has sequences complementary to C-T-G-A-T, one variant of the signal, and to C-A-G, one variant of the 3' terminus of an intron. We also note that the absence of any A-G between -3 and -19 from the 3' splice point may be an essential feature of a strong 3' boundary.     

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.     

Group II introns deleted for multiple substructures retain self-splicing activity.

Group II introns can be folded into highly conserved secondary structures with six major substructures or domains. Domains 1 and 5 are known to play key roles in self-splicing, while the roles of domains 2, 3, 4, and 6 are less clear. A trans assay for domain 5 function has been developed which indicates that domain 5 has a binding site on the precursor RNA that is not predicted from any secondary structure element. In this study, the self-splicing group II intron 5 gamma of the coxI gene of yeast mitochondrial DNA was deleted for various intron domains, singly and in combinations. Those mutant introns were characterized for self-splicing reactions in vitro as a means of locating the domain 5 binding site. A single deletion of domain 2, 3, 4, or 6 does not block in vitro reactions at either splice junction, though the deletion of domain 6 reduces the fidelity of 3' splice site selection somewhat. Even the triple deletion lacking domains 2, 4, and 6 retains some self-splicing activity. The deletion of domains 2, 3, 4, and 6 blocks the reaction at the 3' splice junction but not at the 5' junction. From these results, we conclude that the binding site for domain 5 is within domain 1 and that the complex of 5' exon, domain 1, and domain 5 (plus short connecting sequences) constitutes the essential catalytic core of this intron.