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.     

Spliceosomal snRNAs in the unicellular eukaryote Trichomonas vaginalis are structurally conserved but lack a 5'-cap structure

Few genes in the divergent eukaryote Trichomonas vaginalis have introns, despite the unusually large gene repertoire of this human-infective parasite. These introns are characterized by extended conserved regulatory motifs at the 5' and 3' boundaries, a feature shared with another divergent eukaryote, Giardia lamblia, but not with metazoan introns. This unusual characteristic of T. vaginalis introns led us to examine spliceosomal small nuclear RNAs (snRNAs) predicted to mediate splicing reactions via interaction with intron motifs. Here we identify T. vaginalis U1, U2, U4, U5, and U6 snRNAs, present predictions of their secondary structures, and provide evidence for interaction between the U2/U6 snRNA complex and a T. vaginalis intron. Structural models predict that T. vaginalis snRNAs contain conserved sequences and motifs similar to those found in other examined eukaryotes. These data indicate that mechanisms of intron recognition as well as coordination of the two catalytic steps of splicing have been conserved throughout eukaryotic evolution. Unexpectedly, we found that T. vaginalis spliceosomal snRNAs lack the 5' trimethylguanosine cap typical of snRNAs and appear to possess unmodified 5' ends. Despite the lack of a cap structure, U1, U2, U4, and U5 genes are transcribed by RNA polymerase II, whereas the U6 gene is transcribed by RNA polymerase III.     

Structural features of the maize sus1 gene and protein.

Genomic clones, cDNA clones, and protein of the maize (Zea mays L.) Suc synthase1 (sus1) gene were isolated and sequenced. Termini (5' and 3') of the transcribed unit were identified. The SUS1 protein was purified from tissue culture cells as a phosphorylated protein. The overall structure of sus1 is virtually identical with that of the paralogous gene, shrunken1 (sh1); however, the last intron of sh1 is missing in sus1. This intron bears much sequence similarity with the adjacent exon, suggesting that the intron arose from an internal duplication. Although the placement of the other 14 introns is identical in both genes, the introns exhibit markedly greater differences in size and sequence relative to that shown by the exons. An explanation for the differential rate of divergence of exons and introns is selection pressure for gene function. Additionally, comparisons of coding regions of plant sucrose synthases show that sh1-like and sus1-like genes can be found in all monocots so far analyzed. These latter observations point to an important role played by both genes in this group of plants.     

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.     

Statistical characteristics of eukaryotic intron database

A database called eukaryotic intron database (EID) was developed based on the data from GenBank.Studies on the statistical characteristics of EID show that there were 103,848 genes,478,484 introns,and 582,332 exons,with an average of 4.61 introns and 5.61 exons per gene.Introns of 40-120 nt in length were abundant in the database.Results of the statistical analysis on the data from nine model species showed that in eukaryotes,higher species do not necessarily have more introns or exons in a gene than lower species.Furthermore,characteristics of EID,such as intron phase,distribution of different splice sites,and the relationship between genome size and intron proportion or intron density,have been studied.     

The structure of two fast-white myosin heavy chain promoters. A comparative study

Two complete myosin heavy chain genes were isolated from chicken genomic libraries, and shown to code for fast-white isoforms. Isoform specific probes were developed from the 5' nontranslated regions of the two genes and used to identify the developmental stages at which each of the genes are expressed. One of the genes is transcribed in the embryo and the other only in the adult. The 5' flanking regions of the two genes were sequenced along with the first three exons. The 5' untranslated sequences in both genes are not contiguous, one intron is present in the adult gene while the embryonic gene contains two. The promoters of both genes contain the conserved CAAT and TATA box elements observed in other eucaryotic genes. A computer assisted comparison was performed on the two genes at the nucleotide and amino acid levels. No homology could be detected in the 5' flanking regions of the genes except in and around the CAAT and TATA elements, however, structural sequences at the 5' ends were highly conserved as well as the position of the first three introns. The amino acids in and around the ATP binding site are completely conserved between the two isoforms.     

Can Codon Usage Bias Explain Intron Phase Distributions and Exon Symmetry?

More introns exist between codons (phase 0) than between the first and the second bases (phase 1) or between the second and the third base (phase 2) within the codon. Many explanations have been suggested for this excess of phase 0. It has, for example, been argued to reflect an ancient utility for introns in separating exons that code for separate protein modules. There may, however, be a simple, alternative explanation. Introns typically require, for correct splicing, particular nucleotides immediately 5 in exons (typically a G) and immediately 3 in the following exon (also often a G). Introns therefore tend to be found between particular nucleotide pairs (e.g., G|G pairs) in the coding sequence. If, owing to bias in usage of different codons, these pairs are especially common at phase 0, then intron phase biases may have a trivial explanation. Here we take codon usage frequencies for a variety of eukaryotes and use these to generate random sequences. We then ask about the phase of putative intron insertion sites. Importantly, in all simulated data sets intron phase distribution is biased in favor of phase 0. In many cases the bias is of the magnitude observed in real data and can be attributed to codon usage bias. It is also known that exons may carry either the same phase (symmetric) or different phases (asymmetric) at the opposite ends. We simulated a distribution of different types of exons using frequencies of introns observed in real genes assuming random combination of intron phases at the opposite sides of exons. Surprisingly the simulated pattern was quite similar to that observed. In the simulants we typically observe a prevalence of symmetric exons carrying phase 0 at both ends, which is common for eukaryotic genes. However, at least in some species, the extent of the bias in favor of symmetric (0,0) exons is not as great in simulants as in real genes. These results emphasize the need to construct a biologically relevant null model of successful intron insertion.Reviewing Editor: Dr. Manyuan Long     

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.     

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.     

Correlation between splicing sites within an intron and their sequence complementarity with U1 RNA

We have determined the nucleotide sequence of two short introns (respectively 215 and 90 nucleotides) in the chick alpha 2-collagen (type I) gene as well as parts of the adjacent exons. For one of these introns we find that the 5' end of U1 RNA is complementary not only to the two ends of the intron but also to one end of the intron and sequences inside this intron. These complementarities predict three potential internal splicing sites. By S1 mapping experiments we find three discrete RNA precursors in which different portions of this intron have been deleted. The sizes of the deleted segments are in good agreement with the location of the predicted splicing points inside the intron. The DNA sequence indicates that removal of one portion of the intron should still allow the subsequent elimination of the rest of the intron and the correct splicing of the coding segments located at each end of the intron. The new introns created by the first splicing events contain sequences at each end which are also complementary to U1 RNA. Our data indicate that in the intron which we have examined the sequences at the 3' end of the intron are removed before those at the 5' end.     

Patterns of intron gain and conservation in eukaryotic genes

Background:  

The presence of introns in protein-coding genes is a universal feature of eukaryotic genome organization, and the genes of multicellular eukaryotes, typically, contain multiple introns, a substantial fraction of which share position in distant taxa, such as plants and animals. Depending on the methods and data sets used, researchers have reached opposite conclusions on the causes of the high fraction of shared introns in orthologous genes from distant eukaryotes. Some studies conclude that shared intron positions reflect, almost entirely, a remarkable evolutionary conservation, whereas others attribute it to parallel gain of introns. To resolve these contradictions, it is crucial to analyze the evolution of introns by using a model that minimally relies on arbitrary assumptions.     

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.