Extended base pair complementarity between U1 snRNA and the 5' splice site does not inhibit splicing in higher eukaryotes, but rather increases 5' splice site recognition

Spliceosome formation is initiated by the recognition of the 5′ splice site through formation of an RNA duplex between the 5′ splice site and U1 snRNA. We have previously shown that RNA duplex formation between U1 snRNA and the 5′ splice site can protect pre-mRNAs from degradation prior to splicing. This initial RNA duplex must be disrupted to expose the 5′ splice site sequence for base pairing with U6 snRNA and to form the active spliceosome. Here, we investigated whether hyperstabilization of the U1 snRNA/5′ splice site duplex interferes with splicing efficiency in human cell lines or nuclear extracts. Unlike observations in Saccharomyces cerevisiae, we demonstrate that an extended U1 snRNA/5′ splice site interaction does not decrease splicing efficiency, but rather increases 5′ splice site recognition and exon inclusion. However, low complementarity of the 5′ splice site to U1 snRNA significantly increases exon skipping and RNA degradation. Although the splicing mechanisms are conserved between human and S.cerevisiae, these results demonstrate that distinct differences exist in the activation of the spliceosome.     

hnRNP H binding at the 5' splice site correlates with the pathological effect of two intronic mutations in the NF-1 and TSHbeta genes

We have recently reported a disease-causing substitution (+5G > C) at the donor site of NF-1 exon 3 that produces its skipping. We have now studied in detail the splicing mechanism involved in analyzing RNA–protein complexes at several 5′ splice sites. Characteristic protein patterns were observed by pulldown and band-shift/super-shift analysis. Here, we show that hnRNP H binds specifically to the wild-type GGGgu donor sequence of the NF-1 exon 3. Depletion analyses shows that this protein restricts the accessibility of U1 small nuclear ribonucleoprotein (U1snRNA) to the donor site. In this context, the +5G > C mutation abolishes both U1snRNP base pairing and the 5′ splice site (5′ss) function. However, exon recognition in the mutant can be rescued by disrupting the binding of hnRNP H, demonstrating that this protein enhances the effects of the +5G > C substitution. Significantly, a similar situation was found for a second disease-causing +5G > A substitution in the 5′ss of TSHβ exon 2, which harbors a GGgu donor sequence. Thus, the reason why similar nucleotide substitutions can be either neutral or very disruptive of splicing function can be explained by the presence of specific binding signatures depending on local contexts.     

Complex splicing control of the human Thrombopoietin gene by intronic G runs

The human thrombopoietin (THPO) gene displays a series of alternative splicing events that provide valuable models for studying splicing mechanisms. The THPO region spanning exon 1–4 presents both alternative splicing of exon 2 and partial intron 2 (IVS2) retention following the activation of a cryptic 3′ splice site 85 nt upstream of the authentic acceptor site. IVS2 is particularly rich in stretches of 3–5 guanosines (namely, G1–G10) and we have characterized the role of these elements in the processing of this intron. In vivo studies show that runs G7–G10 work in a combinatorial way to control the selection of the proper 3′ splice site. In particular, the G7 element behaves as the splicing hub of intron 2 and its interaction with hnRNP H1 is critical for the splicing process. Removal of hnRNP H1 by RNA interference promoted the usage of the cryptic 3′ splice site so providing functional evidence that this factor is involved in the selection of the authentic 3′ splice site of THPO IVS2.     

Evidence that U2/U6 helix I promotes both catalytic steps of pre-mRNA splicing and rearranges in between these steps

During pre-mRNA splicing, the spliceosome must configure the substrate, catalyze 5′ splice site cleavage, reposition the substrate, and catalyze exon ligation. The highly conserved U2/U6 helix I, which adjoins sequences that define the reactive sites, has been proposed to configure the substrate for 5′ splice site cleavage and promote catalysis. However, a role for this helix at either catalytic step has not been tested rigorously and previous observations question its role at the catalytic steps. Through a comprehensive molecular genetic study of U2/U6 helix I, we found that weakening U2/U6 helix I, but not mutually exclusive structures, compromised splicing of a substrate limited at the catalytic step of 5′ splice site cleavage, providing the first compelling evidence that this helix indeed configures the substrate during 5′ splice site cleavage. Further, mutations that we proved weaken only U2/U6 helix I suppressed a mutation in PRP16, a DEAH-box ATPase required after 5′ splice site cleavage, providing persuasive evidence that helix I is destabilized by Prp16p and suggesting that this structure is unwound between the catalytic steps. Lastly, weakening U2/U6 helix I also compromised splicing of a substrate limited at the catalytic step of exon ligation, providing evidence that U2/U6 helix I reforms and functions during exon ligation. Thus, our data provide evidence for a fundamental and apparently dynamic role for U2/U6 helix I during the catalytic stages of splicing.     

Exonic Sequences in the 5′ Untranslated Region of α-Tubulin mRNA Modulate trans Splicing in Trypanosoma brucei

Previous studies have identified a conserved AG dinucleotide at the 3′ splice site (3′SS) and a polypyrimidine (pPy) tract that are required for trans splicing of polycistronic pre-mRNAs in trypanosomatids. Furthermore, the pPy tract of the Trypanosoma brucei α-tubulin 3′SS region is required to specify accurate 3′-end formation of the upstream β-tubulin gene and trans splicing of the downstream α-tubulin gene. Here, we employed an in vivo cis competition assay to determine whether sequences other than those of the AG dinucleotide and the pPy tract were required for 3′SS identification. Our results indicate that a minimal α-tubulin 3′SS, from the putative branch site region to the AG dinucleotide, is not sufficient for recognition by the trans-splicing machinery and that polyadenylation is strictly dependent on downstream trans splicing. We show that efficient use of the α-tubulin 3′SS is dependent upon the presence of exon sequences. Furthermore, β-tubulin, but not actin exon sequences or unrelated plasmid sequences, can replace α-tubulin exon sequences for accurate trans-splice-site selection. Taken together, these results support a model in which the informational content required for efficient trans splicing of the α-tubulin pre-mRNA includes exon sequences which are involved in modulation of trans-splicing efficiency. Sequences that positively regulate trans splicing might be similar to cis-splicing enhancers described in other systems.     

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.     

The Exon Splicing Silencer in Human Immunodeficiency Virus Type 1 Tat Exon 3 Is Bipartite and Acts Early in Spliceosome Assembly

Inefficient splicing of human immunodeficiency virus type 1 (HIV-1) RNA is necessary to preserve unspliced and singly spliced viral RNAs for transport to the cytoplasm by the Rev-dependent pathway. Signals within the HIV-1 genome that control the rate of splicing include weak 3′ splice sites, exon splicing enhancers (ESE), and exon splicing silencers (ESS). We have previously shown that an ESS present within tat exon 2 (ESS2) and a suboptimal 3′ splice site together act to inhibit splicing at the 3′ splice site flanking tat exon 2. This occurs at an early step in spliceosome assembly. Splicing at the 3′ splice site flanking tat exon 3 is regulated by a bipartite element composed of an ESE and an ESS (ESS3). Here we show that ESS3 is composed of two smaller elements (AGAUCC and UUAG) that can inhibit splicing independently. We also show that ESS3 is more active in the context of a heterologous suboptimal splice site than of an optimal 3′ splice site. ESS3 inhibits splicing by blocking the formation of a functional spliceosome at an early step, since A complexes are not detected in the presence of ESS3. Competitor RNAs containing either ESS2 or ESS3 relieve inhibition of splicing of substrates containing ESS3 or ESS2. This suggests that a common cellular factor(s) may be required for the inhibition of tat mRNA splicing mediated by ESS2 and ESS3.     

Four exons of the serotonin receptor 4 gene are associated with multiple distant branch points

Splicing of vertebrate introns involves recognition of three consensus elements at the 3′ end. The branch point (BP) and polypyrimidine tract (PPT) are usually located within 40 nucleotides (nt) of the 3′ splice site (3′ ss), AG, but can be much more distant. A characteristic of the region between distant BPs (dBPs) and the 3′ ss is the absence of intervening AG dinucleotides, leading to its designation as the “AG exclusion zone” (AGEZ). The human HTR4 gene, which encodes serotonin receptor 4 and has been associated with schizophrenia, bipolar disease, and gastrointestinal disorders, has four exons with extensive AGEZs. We have mapped the BPs for HTR4 exons 3, 4, 5, and g generated by in vitro splicing, and validated them by mutagenesis in exon-trapping vectors. All exons used dBPs up to 273 nt upstream of the exon. Strikingly, exons 4 and 5 used combinations of both distant and conventionally located BPs, suggesting that successful splicing of these exons can occur by distinct pathways. Our results emphasize the importance for single nucleotide polymorphism resequencing projects to take account of potential dBPs, as the extended AGEZs are vulnerable to mutations that could affect splicing itself or regulation of alternative splicing.     

The U1, U2 and U5 snRNAs crosslink to the 5' exon during yeast pre-mRNA splicing

Activation of pre-messenger RNA (pre-mRNA) splicing requires 5′ splice site recognition by U1 small nuclear RNA (snRNA), which is replaced by U5 and U6 snRNA. Here we use crosslinking to investigate snRNA interactions with the 5′ exon adjacent to the 5′ splice site, prior to the first step of splicing. U1 snRNA was found to interact with four different 5′ exon positions using one specific sequence adjacent to U1 snRNA helix 1. This novel interaction of U1 we propose occurs before U1-5′ splice site base pairing. In contrast, U5 snRNA interactions with the 5′ exon of the pre-mRNA progressively shift towards the 5′ end of U5 loop 1 as the crosslinking group is placed further from the 5′ splice site, with only interactions closest to the 5′ splice site persisting to the 5′ exon intermediate and the second step of splicing. A novel yeast U2 snRNA interaction with the 5′ exon was also identified, which is ATP dependent and requires U2-branchpoint interaction. This study provides insight into the nature and timing of snRNA interactions required for 5′ splice site recognition prior to the first step of pre-mRNA splicing.     

Tra2-Mediated Recognition of HIV-1 5′ Splice Site D3 as a Key Factor in the Processing of vpr mRNA

Small noncoding HIV-1 leader exon 3 is defined by its splice sites A2 and D3. While 3′ splice site (3′ss) A2 needs to be activated for vpr mRNA formation, the location of the vpr start codon within downstream intron 3 requires silencing of splicing at 5′ss D3. Here we show that the inclusion of both HIV-1 exon 3 and vpr mRNA processing is promoted by an exonic splicing enhancer (ESE_vpr) localized between exonic splicing silencer ESSV and 5′ss D3. The ESE_vpr sequence was found to be bound by members of the Transformer 2 (Tra2) protein family. Coexpression of these proteins in provirus-transfected cells led to an increase in the levels of exon 3 inclusion, confirming that they act through ESE_vpr. Further analyses revealed that ESE_vpr supports the binding of U1 snRNA at 5′ss D3, allowing bridging interactions across the upstream exon with 3′ss A2. In line with this, an increase or decrease in the complementarity of 5′ss D3 to the 5′ end of U1 snRNA was accompanied by a higher or lower vpr expression level. Activation of 3′ss A2 through the proposed bridging interactions, however, was not dependent on the splicing competence of 5′ss D3 because rendering it splicing defective but still competent for efficient U1 snRNA binding maintained the enhancing function of D3. Therefore, we propose that splicing at 3′ss A2 occurs temporally between the binding of U1 snRNA and splicing at D3.     

A novel function for the U2AF 65 splicing factor in promoting pre-mRNA 3'-end processing

Splicing and 3′-end processing (including cleavage and polyadenylation) of vertebrate pre-mRNAs are tightly coupled events that contribute to the extensive molecular network that coordinates gene expression. Sequences within the terminal intron of genes are essential to stimulate pre-mRNA 3′-end processing, although the factors mediating this effect are unknown. Here, we show that the pyrimidine tract of the last splice acceptor site of the human β-globin gene is necessary to stimulate mRNA 3′-end formation in vivo and binds the U2AF 65 splicing factor. Naturally occurring β-thalassaemia-causing mutations within the pyrimidine tract reduces both U2AF 65 binding and 3′-end cleavage efficiency. Significantly, a fusion protein containing U2AF 65, when tethered upstream of a cleavage/polyadenylation site, increases 3′-end cleavage efficiency in vitro and in vivo. Therefore, we propose that U2AF 65 promotes 3′-end processing, which contributes to 3′-terminal exon definition.     

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.     

Aberrant 3' splice sites in human disease genes: mutation pattern, nucleotide structure and comparison of computational tools that predict their utilization

The frequency distribution of mutation-induced aberrant 3′ splice sites (3′ss) in exons and introns is more complex than for 5′ splice sites, largely owing to sequence constraints upstream of intron/exon boundaries. As a result, prediction of their localization remains a challenging task. Here, nucleotide sequences of previously reported 218 aberrant 3′ss activated by disease-causing mutations in 131 human genes were compared with their authentic counterparts using currently available splice site prediction tools. Each tested algorithm distinguished authentic 3′ss from cryptic sites more effectively than from de novo sites. The best discrimination between aberrant and authentic 3′ss was achieved by the maximum entropy model. Almost one half of aberrant 3′ss was activated by AG-creating mutations and ~95% of the newly created AGs were selected in vivo. The overall nucleotide structure upstream of aberrant 3′ss was characterized by higher purine content than for authentic sites, particularly in position −3, that may be compensated by more stringent requirements for positive and negative nucleotide signatures centred around position −11. A newly developed online database of aberrant 3′ss will facilitate identification of splicing mutations in a gene or phenotype of interest and future optimization of splice site prediction tools.     

Oriented Scanning Is the Leading Mechanism Underlying 5′ Splice Site Selection in Mammals

Splice site selection is a key element of pre-mRNA splicing. Although it is known to involve specific recognition of short consensus sequences by the splicing machinery, the mechanisms by which 5′ splice sites are accurately identified remain controversial and incompletely resolved. The human F7 gene contains in its seventh intron (IVS7) a 37-bp VNTR minisatellite whose first element spans the exon7–IVS7 boundary. As a consequence, the IVS7 authentic donor splice site is followed by several cryptic splice sites identical in sequence, referred to as 5′ pseudo-sites, which normally remain silent. This region, therefore, provides a remarkable model to decipher the mechanism underlying 5′ splice site selection in mammals. We previously suggested a model for splice site selection that, in the presence of consecutive splice consensus sequences, would stimulate exclusively the selection of the most upstream 5′ splice site, rather than repressing the 3′ following pseudo-sites. In the present study, we provide experimental support to this hypothesis by using a mutational approach involving a panel of 50 mutant and wild-type F7 constructs expressed in various cell types. We demonstrate that the F7 IVS7 5′ pseudo-sites are functional, but do not compete with the authentic donor splice site. Moreover, we show that the selection of the 5′ splice site follows a scanning-type mechanism, precluding competition with other functional 5′ pseudo-sites available on immediate sequence context downstream of the activated one. In addition, 5′ pseudo-sites with an increased complementarity to U1snRNA up to 91% do not compete with the identified scanning mechanism. Altogether, these findings, which unveil a cell type–independent 5′−3′-oriented scanning process for accurate recognition of the authentic 5′ splice site, reconciliate apparently contradictory observations by establishing a hierarchy of competitiveness among the determinants involved in 5′ splice site selection.     

Temperature-dependent splicing of beta-globin pre-mRNA

A T→G mutation at nucleotide 705 of human β-globin intron 2 creates an aberrant 5′ splice site and activates a cryptic 3′ splice site upstream. In consequence, the pre-mRNA is spliced via aberrant splice sites, despite the presence of the still functional correct sites. Surprisingly, when IVS2-705 HeLa or K562 cells were cultured at temperatures below 30°C, aberrant splicing was inhibited and correct splicing was restored. Similar temperature effects were seen for another β-globin pre-mRNA, IVS2-745, and in a construct in which a β-globin intron was inserted into a coding sequence of EGFP. Temperature-induced alternative splicing was affected by the nature of the internal aberrant splice sites flanking the correct sites and by exonic sequences. The results indicate that in the context of thalassemic splicing mutations and possibly in other alternatively spliced pre-mRNAs, temperature is one of the parameters that affect splice site selection.