U5 snRNA interacts with exon sequences at 5' and 3' splice sites.

U5 snRNA is an essential pre-mRNA splicing factor whose function remains enigmatic. Specific mutations in a conserved single-stranded loop sequence in yeast U5 snRNA can activate cleavage of G1----A mutant pre-mRNAs at aberrant 5' splice sites and facilitate processing of dead-end lariat intermediates to mRNA. Activation of aberrant 5' cleavage sites involves base pairing between U5 snRNA and nucleotides upstream of the cleavage site. Processing of dead-end lariat intermediates to mRNA correlates with base pairing between U5 and the first two bases in exon 2. The loop sequence in U5 snRNA may therefore by intimately involved in the transesterification reactions at 5' and 3' splice sites. This pattern of interactions is strikingly reminiscent of exon recognition events in group II self-splicing introns and is consistent with the notion that U5 snRNA may be related to a specific functional domain from a group II-like self-splicing ancestral intron.     

ATP-dependent interaction of yeast U5 snRNA loop 1 with the 5' splice site.

Pre-messenger RNA splicing is a two-step process by which introns are removed and exons joined together. In yeast, the U5 snRNA loop 1 interacts with the 5' exon before the first step of splicing and with the 5' and 3' exons before the second step. In vitro studies revealed that yeast U5 loop 1 is not required for the first step of splicing but is essential for holding the 5' and 3' exons for ligation during the second step. It is critical, therefore, that loop 1 contacts the 5' exon before the first step of splicing to hold this exon following cleavage from the pre-mRNA. At present it is not known how U5 loop 1 is positioned on the 5' exon prior to the first step of splicing. To address this question, we have used site-specific photoactivated crosslinking in yeast spliceosomes to investigate the interaction of U5 loop 1 with the pre-mRNA prior to the first step of splicing. We have found that the highly conserved uridines in loop 1 make ATP-dependent contacts with an approximately 8-nt region at the 5' splice site that includes the invariant GU. These interactions are dependent on functional U2 and U6 snRNAs. Our results support a model where U5 snRNA loop 1 interacts with the 5' exon in two steps during its targeting to the 5' splice site.     

Proximity of conserved U6 and U2 snRNA elements to the 5' splice site region in activated spliceosomes

Major structural changes occur in the spliceosome during its catalytic activation, which immediately precedes the splicing of pre-mRNA. Whereas changes in snRNA conformation are well documented at the level of secondary RNA-RNA interactions, little is known about the tertiary structure of this RNA-RNA network, which comprises the spliceosome's catalytic core. Here, we have used the hydroxyl-radical probe Fe-BABE, tethered to the tenth nucleotide (U(+10)) of the 5' end of a pre-mRNA intron, to map RNA-RNA proximities in spliceosomes. These studies revealed that several conserved snRNA regions are close to U(+10) in activated spliceosomes, namely (i) the U6 snRNA ACAGAG-box region, (ii) portions of the U6 intramolecular stem-loop (U6-ISL) including a nucleotide implicated in the first catalytic step (U74), and (iii) the region of U2 that interacts with the branch point. These data constrain the relative orientation of these structural elements with respect to U(+10) in the activated spliceosome. Upon conversion of the activated spliceosome to complex C, the accessibility of U6-ISL to hydroxyl-radical cleavage is altered, suggesting rearrangements after the first catalytic step.     

Mutations at the 3' splice site can be suppressed by compensatory base changes in U1 snRNA in fission yeast.

U1 snRNA is an essential splicing factor known to base pair with 5' splice sites of premessenger RNAs. We demonstrate that pairing between the universally conserved CU just downstream from the 5' junction interaction region and the 3' splice site AG contributes to efficient splicing of Schizosaccharomyces pombe introns that typify the AG-dependent class described in mammals. Strains carrying mutations in the 3' AG of an artificial intron accumulate linear precursor, indicative of a first step block. Lariat formation is partially restored in these mutants by compensatory changes in nucleotides C7 and U8 of U1 snRNA. Consistent with a general role in fission yeast splicing, mutations at C7 are lethal, while U8 mutants are growth impaired and accumulate linear, unspliced precursor to U6 snRNA. U1 RNA-mediated recognition of the 3' splice site may have origins in analogous intramolecular interactions in an ancestral self-splicing RNA.     

Site-specific crosslinks of yeast U6 snRNA to the pre-mRNA near the 5' splice site.

We have introduced a single photochemical crosslinking reagent into specific sites in the central domain of U6 to identify the sites that are in close proximity to the pre-mRNA substrate. Four distinct U6 snRNAs were synthesized with a single 4-thiouridine (4-thioU) at positions 46, 51, 54, and 57, respectively. Synthetic U6 RNA containing the 4-thioU modifications can functionally reconstitute splicing activity in cell-free yeast splicing extracts depleted of endogenous U6 snRNA. Upon photoactivation with UV (>300 nm), 4-thioU at position 46 forms crosslinks to pre-mRNA near the 5' splice site at nt +4, +5, +6, and +7 in the intron, whereas 4-thioU at position 51 crosslinks to the pre-mRNA at positions -2, -1, +1, +2, +3, and at the invariant G in the lariat intermediate. All crosslinks are dependent on the presence of ATP and the splicing substrate. The two crosslinks to the pre-mRNA from position 46 and 51 of U6 can also occur in prp2 heat-inactivated yeast splicing extracts blocked immediately prior to the first chemical step. Significantly, the crosslink from position 51 can undergo subsequent splicing when the mutant extract is complemented with functional Prp2 protein in a chase experiment, indicating that the crosslink reflects a functional interaction that is maintained during the first step. The crosslink to lariat intermediate appears when the mutant spliceosomes are complemented with functional Prp2 protein added exogenously. This experiment is a paradigm for future studies in which different mutant extracts are used to establish the stage in assembly at which particular RNA-RNA interactions defined by unique crosslinks occur.     

Activation of a cryptic 5' splice site by U1 snRNA

In the course of analyzing 5' splice site mutations in the second intron of Schizosaccharomyces pombe cdc2, we identified a cryptic 5' junction containing a nonconsensus nucleotide at position +2. An even more unusual feature of this cryptic 5' junction was its pattern of activation. By analyzing the profile of splicing products for an extensive series of cdc2 mutants in the presence and absence of compensatory U1 alleles, we have obtained evidence that the natural 5' splice site participates in activation of the cryptic 5' splice site, and that it does so via base pairing to U1 snRNA. Furthermore, the results of follow-up experiments strongly suggest that base pairing between U1 snRNA and the cryptic 5' junction itself plays a dominant role in its activation. Most remarkably, a mutant U1 can activate the cryptic 5' splice site even in the presence of a wild-type sequence at the natural 5' junction, providing unambiguous evidence that this snRNA redirects splicing via base pairing. Although previous work has demonstrated that U5 and U6 snRNAs can activate cryptic 5' splice sites through base pairing interactions, this is the first example in which U1 snRNA has been implicated in the final selection of a cryptic 5' junction.     

Alternative splicing regulation at tandem 3' splice sites

Alternative splicing (AS) constitutes a major mechanism creating protein diversity in humans. Previous bioinformatics studies based on expressed sequence tag and mRNA data have identified many AS events that are conserved between humans and mice. Of these events, ~25% are related to alternative choices of 3′ and 5′ splice sites. Surprisingly, half of all these events involve 3′ splice sites that are exactly 3 nt apart. These tandem 3′ splice sites result from the presence of the NAGNAG motif at the acceptor splice site, recently reported to be widely spread in the human genome. Although the NAGNAG motif is common in human genes, only a small subset of sites with this motif is confirmed to be involved in AS. We examined the NAGNAG motifs and observed specific features such as high sequence conservation of the motif, high conservation of ~30 bp at the intronic regions flanking the 3′ splice site and overabundance of cis-regulatory elements, which are characteristic of alternatively spliced tandem acceptor sites and can distinguish them from the constitutive sites in which the proximal NAG splice site is selected. Our findings imply that AS at tandem splice sites and constitutive splicing of the distal NAG are highly regulated.     

The 5' and 3' splice sites come together via a three dimensional diffusion mechanism.

We present evidence that the splice sites in mammalian pre-mRNAs are brought together via a three dimensional diffusion mechanism. We tested two mechanisms for splice site pairing: a lateral diffusion ('scanning') model and the currently favored three dimensional diffusion ('jumping') model. Two lines of evidence that distinguish between these two models are presented. The first utilized bipartite splicing substrates tethered by double-stranded RNA stems predicted to provide either a moderate or severe block to splice site pairing via a scanning mechanism. Splice site pairing via a jumping mechanism was expected to be unaffected or affected minimally. The second approach utilized a flexible poly(ethylene glycol) moiety within the intron. This insertion was predicted to reduce scanning efficiency but not the efficiency of a three dimensional diffusion mechanism. The best explanation for the data with the bipartite RNAs is that splice site pairing occurs through three dimensional diffusion. Kinetic analysis of the poly(ethylene glycol) containing substrate showed that neither the lag phase nor the initial rates of mRNA production and spliceosome assembly were affected by this insertion. Therefore, both experimental approaches supported the three dimensional diffusion model of splice site pairing.     

Defining a 5' splice site by functional selection in the presence and absence of U1 snRNA 5' end

Pre-mRNA splicing in metazoans is mainly specified by sequences at the termini of introns. We have selected functional 5' splice sites from randomized intron sequences through repetitive rounds of in vitro splicing in HeLa cell nuclear extract. The consensus sequence obtained after one round of selection in normal extract closely resembled the consensus of natural occurring 5' splice sites, suggesting that the selection pressures in vitro and in vivo are similar. After three rounds of selection under competitive splicing conditions, the base pairing potential to the U1 snRNA increased, yielding a G100%U100%R94%A67%G89%U76%R83% intronic consensus sequence. Surprisingly, a nearly identical consensus sequence was obtained when the selection was performed in nuclear extract containing U1 snRNA with a deleted 5' end, suggesting that other factors than the U1 snRNA are involved in 5' splice site recognition. The importance of a consecutive complementarity between the 5' splice site and the U1 snRNA was analyzed systematically in the natural range for in vitro splicing efficiency and complex formation. Extended complementarity was inhibitory to splicing at a late step in spliceosome assembly when pre-mRNA substrates were incubated in normal extract, but favorable for splicing under competitive splicing conditions or in the presence of truncated U1 snRNA where transition from complex A to complex B occurred more rapidly. This suggests that stable U1 snRNA binding is advantageous for assembly of commitment complexes, but inhibitory for the entry of the U4/U6.U5 tri-snRNP, probably due to a delayed release of the U1 snRNP.     

Identification of alternative 5'/3' splice sites based on the mechanism of splice site competition

Alternative splicing plays an important role in regulating gene expression. Currently, most efficient methods use expressed sequence tags or microarray analysis for large-scale detection of alternative splicing. However, it is difficult to detect all alternative splice events with them because of their inherent limitations. Previous computational methods for alternative splicing prediction could only predict particular kinds of alternative splice events. Thus, it would be highly desirable to predict alternative 5'/3' splice sites with various splicing levels using genomic sequences alone. Here, we introduce the competition mechanism of splice sites selection into alternative splice site prediction. This approach allows us to predict not only rarely used but also frequently used alternative splice sites. On a dataset extracted from the AltSplice database, our method correctly classified approximately 70% of the splice sites into alternative and constitutive, as well as approximately 80% of the locations of real competitors for alternative splice sites. It outperforms a method which only considers features extracted from the splice sites themselves. Furthermore, this approach can also predict the changes in activation level arising from mutations in flanking cryptic splice sites of a given splice site. Our approach might be useful for studying alternative splicing in both computational and molecular biology.     

Natural base-pairing interactions between 5' splice site and branch site sequences affect mammalian 5' splice site selection.

In the murine gene encoding the neuronal cell adhesion molecule (NCAM), the integrity of the 5' splice site of exon 18 (E18) is essential for regulation of alternative splicing. To further study the contribution of 5' splice site sequences, we used a simple NCAM pre-mRNA containing a portion of E18 fused to E19 and separated by a shortened intron. This RNA is spliced in vitro to produce five sets of lariat intermediates and products, the most abundant set displaying aberrant migration in acrylamide/urea gels. Base pairing interactions between positions +5 and +8 of the intron and positions -3 and -6 from the branch point were responsible for the faster migration of this set of lariat molecules. To test whether the duplex structure forms earlier and contributes to 5' splice site selection, we used NCAM substrates carrying the 5' splice sites of E17 and E18 in competition for the 3' splice site of E19. Mutations upstream of the major branch site improve E18/E19 splicing in NIH3T3 extracts, whereas compensatory mutations at positions +7 and +8 neutralize the effect of branch site mutations and curtail E18/E19 splicing. Our data suggest that duplex formation occurs early and interferes with the assembly of complexes initiated on the 5' splice site of NCAM E18. This novel type of intron interaction may exist in the introns of other mammalian pre-mRNAs.     

Extensive interactions of PRP8 protein with the 5'' and 3'' splice sites during splicing suggest a role in stabilization of exon alignment by U5 snRNA.

Precursor RNAs containing 4-thiouridine at specific sites were used with UV-crosslinking to map the binding sites of the yeast protein splicing factor PRP8. PRP8 protein interacts with a region of at least eight exon nucleotides at the 5' splice site and a minimum of 13 exon nucleotides and part of the polypyrimidine tract in the 3' splice site region. Crosslinking of PRP8 to mutant and duplicated 3' splice sites indicated that the interaction is not sequence specific, nor does it depend on the splice site being functional. Binding of PRP8 to the 5' exon was established before step 1 and to the 3' splice site region after step 1 of splicing. These interactions place PRP8 close to the proposed catalytic core of the spliceosome during both transesterification reactions. To date, this represents the most extensive mapping of the binding site(s) of a splicing factor on the substrate RNA. We propose that the large binding sites of PRP8 stabilize the intrinsically weaker interactions of U5 snRNA with both exons at the splice sites for exon alignment by the U5 snRNP.     

Conserved signals around the 5' splice sites in eukaryotic nuclear precursor mRNAs: G-runs are frequent in the introns and C in the exons near both 5' and 3' splice sites

It is known that the GT doublet is well conserved at the 5' exon/intron splice junction and is frequently embedded in the AGGT quartet. Although only the underlined G is invariable, splicing and ligation are accurately executed. In this work we search for additional conserved potential signals which may aid in 5' splice site recognition. Extensive searches which are not limited to a preconceived consensus sequence are carried out. We investigate the distributions of the 256 quartets in a 1000 nucleotide span around the 5' splice sites in approximately 1700 eukaryotic nuclear precursor mRNAs. Several potential signals are noted. Of particular interest are quartets containing runs of G, e.g., G4, G3T, G3C, G3A and AG3 in the intron immediately downstream and some C-containing quartets in the exon upstream of the 5' splice site. In an analogous calculation, (A)GGG(A) has also been found to be frequent in the intron, 60 nucleotides upstream and (A)CCC(A) in the exon downstream of the 3' splice site. These results are consistent with the recent indications that exon sequences may play a role in efficient splicing. Some models are proposed.     

Intrinsic differences between authentic and cryptic 5' splice sites

Cryptic splice sites are used only when use of a natural splice site is disrupted by mutation. To determine the features that distinguish authentic from cryptic 5′ splice sites (5′ss), we systematically analyzed a set of 76 cryptic 5′ss derived from 46 human genes. These cryptic 5′ss have a similar frequency distribution in exons and introns, and are usually located close to the authentic 5′ss. Statistical analysis of the strengths of the 5′ss using the Shapiro and Senapathy matrix revealed that authentic 5′ss have significantly higher score values than cryptic 5′ss, which in turn have higher values than the mutant ones. β-Globin provides an interesting exception to this rule, so we chose it for detailed experimental analysis in vitro. We found that the sequences of the β-globin authentic and cryptic 5′ss, but not their surrounding context, determine the correct 5′ss choice, although their respective scores do not reflect this functional difference. Our analysis provides a statistical basis to explain the competitive advantage of authentic over cryptic 5′ss in most cases, and should facilitate the development of tools to reliably predict the effect of disease-associated 5′ss-disrupting mutations at the mRNA level.     

Mutations in yeast U5 snRNA alter the specificity of 5' splice-site cleavage

Recognition of 5' splice sites in pre-mRNA splicing is achieved in part by base pairing with U1 snRNA. We have used interactive suppression in the yeast Saccharomyces cerevisiae to look for other factors involved in 5' splice-site recognition. This approach identified an extragenic suppressor that activates a cryptic 5' splice site. The suppressor is a gene for U5 snRNA (snR7) with a single base mutation in a strictly conserved 9 base sequence. This suggests that U5 snRNA can play a part in determining the position of 5' splice-site cleavage. Consistent with this, we have been able to isolate other mutations in the 9 base element in U5 snRNA that specifically activate a second cryptic 5' splice site nearby.     

Efficient internal exon recognition depends on near equal contributions from the 3' and 5' splice sites

Pre-mRNA splicing is carried out by the spliceosome, which identifies exons and removes intervening introns. In vertebrates, most splice sites are initially recognized by the spliceosome across the exon, because most exons are small and surrounded by large introns. This gene architecture predicts that efficient exon recognition depends largely on the strength of the flanking 3' and 5' splice sites. However, it is unknown if the 3' or the 5' splice site dominates the exon recognition process. Here, we test the 3' and 5' splice site contributions towards efficient exon recognition by systematically replacing the splice sites of an internal exon with sequences of different splice site strengths. We show that the presence of an optimal splice site does not guarantee exon inclusion and that the best predictor for exon recognition is the sum of both splice site scores. Using a genome-wide approach, we demonstrate that the combined 3' and 5' splice site strengths of internal exons provide a much more significant separator between constitutive and alternative exons than either the 3' or the 5' splice site strength alone.