In vitro and in silico analysis reveals an efficient algorithm to predict the splicing consequences of mutations at the 5' splice sites

We have found that two previously reported exonic mutations in the PINK1 and PARK7 genes affect pre-mRNA splicing. To develop an algorithm to predict underestimated splicing consequences of exonic mutations at the 5′ splice site, we constructed and analyzed 31 minigenes carrying exonic splicing mutations and their derivatives. We also examined 189 249 U2-dependent 5′ splice sites of the entire human genome and found that a new variable, the SD-Score, which represents a common logarithm of the frequency of a specific 5′ splice site, efficiently predicts the splicing consequences of these minigenes. We also employed the information contents (R_i) to improve the prediction accuracy. We validated our algorithm by analyzing 32 additional minigenes as well as 179 previously reported splicing mutations. The SD-Score algorithm predicted aberrant splicings in 198 of 204 sites (sensitivity = 97.1%) and normal splicings in 36 of 38 sites (specificity = 94.7%). Simulation of all possible exonic mutations at positions −3, −2 and −1 of the 189 249 sites predicts that 37.8, 88.8 and 96.8% of these mutations would affect pre-mRNA splicing, respectively. We propose that the SD-Score algorithm is a practical tool to predict splicing consequences of mutations affecting the 5′ splice site.     

Optimization of a weak 3' splice site counteracts the function of a bovine papillomavirus type 1 exonic splicing suppressor in vitro and in vivo

Alternative splicing is a critical component of the early to late switch in papillomavirus gene expression. In bovine papillomavirus type 1 (BPV-1), a switch in 3' splice site utilization from an early 3' splice site at nucleotide (nt) 3225 to a late-specific 3' splice site at nt 3605 is essential for expression of the major capsid (L1) mRNA. Three viral splicing elements have recently been identified between the two alternative 3' splice sites and have been shown to play an important role in this regulation. A bipartite element lies approximately 30 nt downstream of the nt 3225 3' splice site and consists of an exonic splicing enhancer (ESE), SE1, followed immediately by a pyrimidine-rich exonic splicing suppressor (ESS). A second ESE (SE2) is located approximately 125 nt downstream of the ESS. We have previously demonstrated that the ESS inhibits use of the suboptimal nt 3225 3' splice site in vitro through binding of cellular splicing factors. However, these in vitro studies did not address the role of the ESS in the regulation of alternative splicing. In the present study, we have analyzed the role of the ESS in the alternative splicing of a BPV-1 late pre-mRNA in vivo. Mutation or deletion of just the ESS did not significantly change the normal splicing pattern where the nt 3225 3' splice site is already used predominantly. However, a pre-mRNA containing mutations in SE2 is spliced predominantly using the nt 3605 3' splice site. In this context, mutation of the ESS restored preferential use of the nt 3225 3' splice site, indicating that the ESS also functions as a splicing suppressor in vivo. Moreover, optimization of the suboptimal nt 3225 3' splice site counteracted the in vivo function of the ESS and led to preferential selection of the nt 3225 3' splice site even in pre-mRNAs with SE2 mutations. In vitro splicing assays also showed that the ESS is unable to suppress splicing of a pre-mRNA with an optimized nt 3225 3' splice site. These data confirm that the function of the ESS requires a suboptimal upstream 3' splice site. A surprising finding of our study is the observation that SE1 can stimulate both the first and the second steps of splicing.     

Exon mutations uncouple 5' splice site selection from U1 snRNA pairing

It has previously been shown that a mutation of yeast 5' splice junctions at position 5 (GUAUGU) causes aberrant pre-mRNA cleavages near the correct 5' splice site. We show here that the addition of exon mutations to an aberrant cleavage site region transforms it into a functional 5' splice site both in vivo and in vitro. The aberrant mRNAs are translated in vivo. The results suggest that the highly conserved G at the 5' end of introns is necessary for the second step of splicing. Further analyses indicate that the location of the U1 snRNA-pre-mRNA pairing is not affected by the exon mutations and that the precise 5' splice site is selected independent of this pairing.     

A U1 snRNA:pre-mRNA base pairing interaction is required early in yeast spliceosome assembly but does not uniquely define the 5'' cleavage site.

We analyzed the effects of suppressor mutations in the U1 snRNA (SNR19) gene from Saccharomyces cerevisiae on the splicing of mutant pre-mRNA substrates. The results indicate that pairing between U1 snRNA and the highly conserved position 5 (GTATGT) of the intron occurs early in spliceosome assembly in vitro. This pairing is important for efficient splicing both in vitro and in vivo. However, pairing at position 5 does not appear to influence 5' splice site selection in vivo, indicating that the previously described U1 snRNA:5' splice junction base pairing interaction is not sufficient to define the 5' cleavage site.     

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.