A catalogue of splice junction and putative branch point sequences from plant introns.

Splice junction and possible branch point sequences have been collected from 177 plant introns. Consensus sequences for the 5' and 3' splice junctions and for possible branch points have been derived. The splice junction consensus sequences were virtually identical to those of animal introns except that the polypyrimidine stretch at the 3' splice junction was less pronounced in the plant introns. A search for possible branch points with sequences related to the yeast, vertebrate and fungal consensus sequences revealed a similar sequence in plant introns.     

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.     

Mutation of putative branchpoint consensus sequences in plant introns reduces splicing efficiency

Intron lariat formation between the 5' end of an intron and a branchpoint adenosine is a fundamental aspect of the first step in animal and yeast nuclear pre-mRNA splicing. Despite similarities in intron sequence requirements and the components of splicing, differences exist between the splicing of plant and vertebrate introns. The identification of AU-rich sequences as major functional elements in plant introns and the demonstration that a branchpoint consensus sequence was not required for splicing have led to the suggestion that the transition from AU-rich intron to GC-rich exon is a major potential signal by which plant pre-mRNA splice sites are recognized. The role of putative branchpoint sequences as an internal signal in plant intron recognition/definition has been re-examined. Single nucleotide mutations in putative branchpoint adenosines contained within CUNAN sequences in four different plant introns all significantly reduced splicing efficiency. These results provide the most direct evidence to date for preferred branchpoint sequences being required for the efficient splicing of at least some plant introns in addition to the important role played by AU sequences in dicot intron recognition. The observed patterns of 3' splice site selection in the introns studied are consistent with the scanning model described for animal intron 3' splice site selection. It is suggested that, despite the clear importance of AU sequences for plant intron splicing, the fundamental processes of splice site selection and splicing in plants are similar to those in animals.     

Conserved putative signals in 3' intron junctions in rodents

Several 3' splice signals are known todate. At the 3' splice site an AG doublet is frequently found. Just upstream of the splice site there is a string of 6-11 pyrimidines. More recently it has been found that one of the stages in the splicing process involves formation of a lariat, in which the 5' end of the intron forms a 2'-5' branch with an A residue located 18-37 nucleotides upstream of the 3' splice site. The branching-point consensus is weakly defined and consists of the sequence YNYTRAY, where Y is a pyrimidine, R a purine and N any base. The A in the sixth position is the one with which branching occurs. Here we present the results of extensive searches for additional putative signals around the branching-point consensus and the 3' splice site in rodent nuclear precursor mRNAs. The signals obtained for the over 370 rodent introns are compared with those found in a larger eukaryotic sample containing over 900 nuclear pre-mRNA introns. Of particular interest are GGGA and CCCA. In both analyses GGGA occurs about 60 nucleotides upstream and CCCA is found 3-40 nucleotides downstream from the 3' splice site. A model explaining some of the putative signals discussed here is also proposed. This model involves formation of alternate stem-loop structures around the branching point and 3' splice site. Such signals and structures can possibly aid in protein or nucleoprotein branching point and splice site recognition.     

Conserved Putative Signals in 3′ Intron Junctions in Rodents

Abstract

Several 3′ splice signals are known todate. At the 3′ splice site an AG doublet is frequently found. Just upstream of the splice site there is a string of 6–11 pyrimidines. More recently it has been found that one of the stages in the splicing process involves formation of a lariat, in which the 5′ end of the intron forms a 2′-5′ branch with an A residue located 18–37 nucleotides upstream of the 3′ splice site. The branching-point consensus is weakly defined and consists of the sequence YNYTRAY, where Y is a pyrimidine, R a purine and N any base. The A in the sixth position is the one with which branching occurs. Here we present the results of extensive searches for additional putative signals around the branching-point consensus and the 3′ splice site in rodent nuclear precursor mRNAs. The signals obtained for the over 370 rodent introns are compared with those found in a larger eukaryotic sample containing over 900 nuclear pre-mRNA introns. Of particular interest are GGGA and CCCA In both analyses GGGA occurs about 60 nucleotides upstream and CCCA is found 3–40 nucleotides downstream from the 3′ splice site. A model explaining some of the putative signals discussed here is also proposed. This model involves formation of alternate stem-loop structures around the branching point and 3′ splice site. Such signals and structures can possibly aid in protein or nucleoprotein branching point and splice site recognition.     

Splice site consensus sequences are preferentially accessible to nucleases in isolated adenovirus RNA.

The conformation of RNA sequences spanning five 3' splice sites and two 5' splice sites in adenovirus mRNA was probed by partial digestion with single-strand specific nucleases. Although cleavage of nucleotides near both 3' and 5' splice sites was observed, most striking was the preferential digestion of sequences near the 3' splice site. At each 3' splice site a region of very strong cleavage is observed at low concentrations of enzyme near the splice site consensus sequence or the upstream branch point consensus sequence. Additional sites of moderately strong cutting near the branch point consensus sequence were observed in those sequences where the splice site was the preferred target. Since recognition of the 3' splice site and branch site appear to be early events in mRNA splicing these observations may indicate that the local conformation of the splice site sequences may play a direct or indirect role in enhancing the accessibility of sequences important for splicing.     

Different effects of intron nucleotide composition and secondary structure on pre-mRNA splicing in monocot and dicot plants.

We have found previously that the sequences important for recognition of pre-mRNA introns in dicot plants differ from those in the introns of vertebrates and yeast. Neither a conserved branch point nor a polypyrimidine tract, found in yeast and vertebrate introns respectively, are required. Instead, AU-rich sequences, a characteristic feature of dicot plant introns, are essential. Here we show that splicing in protoplasts of maize, a monocot, differs significantly from splicing in a dicot, Nicotiana plumbaginifolia. As in the case of dicots, a conserved branch point and a polypyrimidine tract are not required for intron processing in maize. However, unlike in dicots, AU-rich sequences are not essential, although their presence facilitates splicing if the splice site sequences are not optimal. The lack of an absolute requirement for AU-rich stretches in monocot introns in reflected in the occurrence of GC-rich introns in monocots but not in dicots. We also show that maize protoplasts are able to process a mammalian intron and short introns containing stem--loops, neither of which are spliced in N.plumbaginifolia protoplasts. The ability of maize, but not of N.plumbaginifolia to process stem--loop-containing or GC-rich introns suggests that one of the functions of AU-rich sequences during splicing of dicot plant pre-mRNAs may be to minimize secondary structure within the intron.     

Distribution and consensus of branch point signals in eukaryotic genes: a computerized statistical analysis.

An intermediate stage in the process of eukaryotic RNA splicing is the formation of a lariat structure. It is anchored at an adenosine residue in intron between 10 and 50 nucleotides upstream of the 3' splice site. A short conserved sequence (the branch point sequence) functions as the recognition signal for the site of lariat formation. It has been generally assumed that the branch point is recognized mainly by the presence of its unique sequence where the lariat is formed. However, the known branch point consensus sequence is found to be distributed nearly randomly throughout the gene sequence with only a slightly higher frequency in the expected lariat region. Further, the known consensus sequence is found to be clearly inadequate to specify branch points. These observations have implications for understanding the mechanism of branch point recognition in the process of splicing, and the possible evolution of the branch point signal.     

Mutational analysis of the U12-dependent branch site consensus sequence

Highly conserved sequences at the 5′ splice site and branch site of U12-dependent introns are important determinants for splicing by U12-dependent spliceosomes. This study investigates the in vivo splicing phenotypes of mutations in the branch site consensus sequence of the U12-dependent intron F from a human NOL1 (P120) minigene. Intron F contains a fully consensus branch site sequence (UUCCUUAAC). Mutations at each position were analyzed for their effects on U12-dependent splicing in vivo. Mutations at most positions resulted in a significant reduction of correct U12-dependent splicing. Defects observed included increased unspliced RNA levels, the activation of cryptic U2-dependent 5′ and 3′ splice sites, and the activation of cryptic U12-dependent branch/3′ splice sites. A strong correlation was observed between the predicted thermodynamic stability of the branch site: U12 snRNA interaction and correct U12-dependent splicing. The lack of a polypyrimidine tract between the branch site and 3′ splice site of U12-dependent introns and the observed reliance on base-pairing interactions for correct U12-dependent splicing emphasize the importance of RNA/RNA interactions during U12-dependent intron recognition and proper splice site selection.     

Identification,characterization and molecular phylogeny of U12-dependent introns in the Arabidopsis thaliana genome

U12-dependent introns are spliced by the minor U12-type spliceosome and occur in a variety of eukaryotic organisms, including Arabidopsis. In this study, a set of putative U12-dependent introns was compiled from a large collection of cDNA/EST- confirmed introns in the Arabidopsis thaliana genome by means of high-throughput bioinformatic analysis combined with manual scrutiny. A total of 165 U12-type introns were identified based upon stringent criteria. This number of sequences well exceeds the total number of U12-type introns previously reported for plants and allows a more thorough statistical analysis of U12-type signals. Of particular note is the discovery that the distance between the branch site adenosine and the acceptor site ranges from 10 to 39 nt, significantly longer than the previously postulated limit of 21 bp. Further analysis indicates that, in addition to the spacing constraint, the sequence context of the potential acceptor site may have an important role in 3′ splice site selection. Several alternative splicing events involving U12-type introns were also captured in this study, providing evidence that U12-dependent acceptor sites can also be recognized by the U2-type spliceosome. Furthermore, phylogenetic analysis suggests that both U12-type AT-AC and U12-type GT-AG introns occurred in Na⁺/H⁺ antiporters in a progenitor of animals and plants.     

The spliceosomal snRNAs of Caenorhabditis elegans.

Nematodes are the only group of organisms in which both cis- and trans-splicing of nuclear mRNAs are known to occur. Most Caenorhabditis elegans introns are exceptionally short, often only 50 bases long. The consensus donor and acceptor splice site sequences found in other animals are used for both cis- and trans-splicing. In order to identify the machinery required for these splicing events, we have characterized the C. elegans snRNAs. They are similar in sequence and structure to those characterized in other organisms, and several sequence variations discovered in the nematode snRNAs provide support for previously proposed structure models. The C. elegans snRNAs are encoded by gene families. We report here the sequences of many of these genes. We find a highly conserved sequence, the proximal sequence element (PSE), about 65 bp upstream of all 21 snRNA genes thus far sequenced, including the SL RNA genes, which specify the snRNAs that provide the 5' exons in trans-splicing. The sequence of the C. elegans PSE is distinct from PSE's from other organisms.     

Human GC-AG alternative intron isoforms with weak donor sites show enhanced consensus at acceptor exon positions

It has been previously observed that the intrinsically weak variant GC donor sites, in order to be recognized by the U2-type spliceosome, possess strong consensus sequences maximized for base pair formation with U1 and U5/U6 snRNAs. However, variability in signal strength is a fundamental mechanism for splice site selection in alternative splicing. Here we report human alternative GC-AG introns (for the first time from any species), and show that while constitutive GC-AG introns do possess strong signals at their donor sites, a large subset of alternative GC-AG introns possess weak consensus sequences at their donor sites. Surprisingly, this subset of alternative isoforms shows strong consensus at acceptor exon positions 1 and 2. The improved consensus at the acceptor exon can facilitate a strong interaction with U5 snRNA, which tethers the two exons for ligation during the second step of splicing. Further, these isoforms nearly always possess alternative acceptor sites and exhibit particularly weak polypyrimidine tracts characteristic of AG-dependent introns. The acceptor exon nucleotides are part of the consensus required for the U2AF³⁵-mediated recognition of AG in such introns. Such improved consensus at acceptor exons is not found in either normal or alternative GT-AG introns having weak donor sites or weak polypyrimidine tracts. The changes probably reflect mechanisms that allow GC-AG alternative intron isoforms to cope with two conflicting requirements, namely an apparent need for differential splice strength to direct the choice of alternative sites and a need for improved donor signals to compensate for the central mismatch base pair (C-A) in the RNA duplex of U1 snRNA and the pre-mRNA. The other important findings include (i) one in every twenty alternative introns is a GC-AG intron, and (ii) three of every five observed GC-AG introns are alternative isoforms.     

Application of learning techniques to splicing site recognition

Most genes of eukaryotic genomes are disrupted by introns. The application of a learning technique which uses both statistic and syntactic analysis lead to the establishment of logical rules enabling the recognition of intron/exon junctions between uncoding and coding sequences. The rules were tested on rat actin gene sequences containing some or all of the introns and 50 exon nucleotides on either side of the intron. The results show good recognition of the excision site. This recognition is more ambiguous when the sequence is short; for the acceptor sequence it presents a good selection. The learning achieved with both the donor and acceptor sequence does not lead to recognition. This result indicates that it is not the relationship between donor and acceptor sites in the same intron which determines sequence selection or the splicing mechanism.     

The plastid aldolase gene fromChlamydomonas reinhardtii: Intron/exon organization,evolution, and promoter structure

Genomic clones encoding the plastidic fructose- 1,6-bisphosphate aldolase ofChlamydomonas reinhardtii were isolated and sequenced. The gene contains three introns which are located within the coding sequence for the mature protein. No introns are located within or near the sequence encoding the transit-peptide, in contrast to the genes for plastidic aldolases of higher plants. Neither the number nor the positions of the three introns of theC. reinhardtii aldolase gene are conserved in the plastidic or cytosolic aldolase genes of higher plants and animals. The 5′ border sequences of introns in the aldolase gene ofC. reinhardtii exhibit the conserved plant consensus sequence. The 3′ acceptor splice sites for introns 1 and 3 show much less similarity to the eukaryotic consensus sequences than do those of intron 2. The plastidic aldolase gene has two tandemly repeated CAAT box motifs in the promoter region. Genomic Southern blots indicate that the gene is encoded by a single locus in theC. reinhardtii genome.     

Nuclear pre-mRNA splicing in the fission yeast Schizosaccharomyces pombe strictly requires an intron-contained, conserved sequence element.

It has recently been argued that pre-mRNA splicing in the fission yeast Schizosaccharomyces pombe may be more similar to splicing in metazoan species than in the budding yeast Saccharomyces cerevisiae. In this report we show that, contrary to this assumption, the conserved sequence element 5'-CTPu APy-3' found in all S. pombe introns 6-18 nucleotides upstream of the 3' splice site is, like the TACTAAC box in S. cerevisiae, indispensable for efficient splicing. The conserved adenine residue of this sequence is used for branch formation and point mutations introduced into the CTPuAPy sequence abolish splicing and seem not to result in the recruitment of cryptic branch sites. We also show that an S. cerevisiae intron is correctly excised in S. pombe whereby the TACTAAC box is used in branch formation.     

An SF1 affinity model to identify branch point sequences in human introns

Splicing factor 1 (SF1) binds to the branch point sequence (BPS) of mammalian introns and is believed to be important for the splicing of some, but not all, introns. To help identify BPSs, particularly those that depend on SF1, we generated a BPS profile model in which SF1 binding affinity data, validated by branch point mapping, were iteratively incorporated into computational models. We searched a data set of 117,499 human introns for best matches to the SF1 Affinity Model above a threshold, and counted the number of matches at each intronic position. After subtracting a background value, we found that 87.9% of remaining high-scoring matches identified were located in a region upstream of 3'-splice sites where BPSs are typically found. Since U2AF65 recognizes the polypyrimidine tract (PPT) and forms a cooperative RNA complex with SF1, we combined the SF1 model with a PPT model computed from high affinity binding sequences for U2AF65. The combined model, together with binding site location constraints, accurately identified introns bound by SF1 that are candidates for SF1-dependent splicing.