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.     

Biased exon/intron distribution of cryptic and de novo 3' splice sites

We compiled sequences of previously published aberrant 3′ splice sites (3′ss) that were generated by mutations in human disease genes. Cryptic 3′ss, defined here as those resulting from a mutation of the 3′YAG consensus, were more frequent in exons than in introns. They clustered in ~20 nt region adjacent to authentic 3′ss, suggesting that their under-representation in introns is due to a depletion of AG dinucleotides in the polypyrimidine tract (PPT). In contrast, most aberrant 3′ss that were induced by mutations outside the 3′YAG consensus (designated ‘de novo’) were in introns. The activation of intronic de novo 3′ss was largely due to AG-creating mutations in the PPT. In contrast, exonic de novo 3′ss were more often induced by mutations improving the PPT, branchpoint sequence (BPS) or distant auxiliary signals, rather than by direct AG creation. The Shapiro–Senapathy matrix scores had a good prognostic value for cryptic, but not de novo 3′ss. Finally, AG-creating mutations in the PPT that produced aberrant 3′ss upstream of the predicted BPS in vivo shared a similar ‘BPS-new AG’ distance. Reduction of this distance and/or the strength of the new AG PPT in splicing reporter pre-mRNAs improved utilization of authentic 3′ss, suggesting that AG-creating mutations that are located closer to the BPS and are preceded by weaker PPT may result in less severe splicing defects.     

Global control of aberrant splice-site activation by auxiliary splicing sequences: evidence for a gradient in exon and intron definition

Auxiliary splicing signals play a major role in the regulation of constitutive and alternative pre-mRNA splicing, but their relative importance in selection of mutation-induced cryptic or de novo splice sites is poorly understood. Here, we show that exonic sequences between authentic and aberrant splice sites that were activated by splice-site mutations in human disease genes have lower frequencies of splicing enhancers and higher frequencies of splicing silencers than average exons. Conversely, sequences between authentic and intronic aberrant splice sites have more enhancers and less silencers than average introns. Exons that were skipped as a result of splice-site mutations were smaller, had lower SF2/ASF motif scores, a decreased availability of decoy splice sites and a higher density of silencers than exons in which splice-site mutation activated cryptic splice sites. These four variables were the strongest predictors of the two aberrant splicing events in a logistic regression model. Elimination or weakening of predicted silencers in two reporters consistently promoted use of intron-proximal splice sites if these elements were maintained at their original positions, with their modular combinations producing expected modification of splicing. Together, these results show the existence of a gradient in exon and intron definition at the level of pre-mRNA splicing and provide a basis for the development of computational tools that predict aberrant splicing outcomes.     

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.     

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.     

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.     

Multiple features contribute to efficient constitutive splicing of an unusually large exon

Vertebrate internal exons are usually between 50 and 400 nt long; exons outside this size range may require additional exonic and/or intronic sequences to be spliced into the mature mRNA. The mouse polymeric immunoglobulin receptor gene has a 654 nt exon that is efficiently spliced into the mRNA. We have examined this exon to identify features that contribute to its efficient splicing despite its large size; a large constitutive exon has not been studied previously. We found that a strong 5′ splice site is necessary for this exon to be spliced intact, but the splice sites alone were not sufficient to efficiently splice a large exon. At least two exonic sequences and one evolutionarily conserved intronic sequence also contribute to recognition of this exon. However, these elements have redundant activities as they could only be detected in conjunction with other mutations that reduced splicing efficiency. Several mutations activated cryptic 5′ splice sites that created smaller exons. Thus, the balance between use of these potential sites and the authentic 5′ splice site must be modulated by sequences that repress or enhance use of these sites, respectively. Also, sequences that enhance cryptic splice site use must be absent from this large exon.     

Functional studies on the ATM intronic splicing processing element

In disease-associated genes, the understanding of the functional significance of deep intronic nucleotide variants may represent a difficult challenge. We have previously reported a new disease-causing mechanism that involves an intronic splicing processing element (ISPE) in ATM, composed of adjacent consensus 5′ and 3′ splice sites. A GTAA deletion within ISPE maintains potential adjacent splice sites, disrupts a non-canonical U1 snRNP interaction and activates an aberrant exon. In this paper, we demonstrate that binding of U1 snRNA through complementarity within a ~40 nt window downstream of the ISPE prevents aberrant splicing. By selective mutagenesis at the adjacent consensus ISPE splice sites, we show that this effect is not due to a resplicing process occurring at the ISPE. Functional comparison of the ATM mouse counterpart and evaluation of the pre-mRNA splicing intermediates derived from affected cell lines and hybrid minigene assays indicate that U1 snRNP binding at the ISPE interferes with the cryptic acceptor site. Activation of this site results in a stringent 5′–3′ order of intron sequence removal around the cryptic exon. Artificial U1 snRNA loading by complementarity to heterologous exonic sequences represents a potential therapeutic method to prevent the usage of an aberrant CFTR cryptic exon. Our results suggest that ISPE-like intronic elements binding U1 snRNPs may regulate correct intron processing.     

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.     

Delineation of the mechanisms of aberrant splicing caused by two unusual intronic mutations in the RSK2 gene involved in Coffin-Lowry syndrome

Coffin–Lowry syndrome (CLS) is caused by mutations in the RSK2 gene encoding a protein kinase of the Ras signalling pathway. We have studied two point mutations which cause aberrant splicing but do not concern the invariant GT or AG nucleotides of splice sites. The first, an A→G transition at position +3 of the 5′ splice site of exon 6, results in vivo and in vitro in exon skipping and premature translation termination. The natural 5′ splice site, although intrinsically weak, is not transactivated under normal conditions. Consequently, replacement of an A/U by a G/U base pairing with U1 snRNA reduces its strength below a critical threshold. The second mutation, an A→G transition 11 nt upstream of exon 5, creates a new AG near the natural 3′ splice site. In vitro this synthetic 3′ AG is used exclusively by the splicing machinery. In vivo this splicing event is also observed, but is underestimated because the resulting RSK2 mRNA contains premature stop codons which trigger the nonsense-mediated decay process. We show that a particular mechanism is involved in the aberrant splicing of exon 5, implying involvement of the natural 3′ AG during the first catalytic step and the new 3′ AG during the second step. Thus, our results explain how these mutations cause severe forms of CLS.     

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.     

Insights into the selective activation of alternatively used splice acceptors by the human immunodeficiency virus type-1 bidirectional splicing enhancer

The guanosine-adenosine-rich exonic splicing enhancer (GAR ESE) identified in exon 5 of the human immunodeficiency virus type-1 (HIV-1) pre-mRNA activates either an enhancer-dependent 5′ splice site (ss) or 3′ ss in 1-intron reporter constructs in the presence of the SR proteins SF2/ASF2 and SRp40. Characterizing the mode of action of the GAR ESE inside the internal HIV-1 exon 5 we found that this enhancer fulfils a dual splicing regulatory function (i) by synergistically mediating exon recognition through its individual SR protein-binding sites and (ii) by conferring 3′ ss selectivity within the 3′ ss cluster preceding exon 5. Both functions depend upon the GAR ESE, U1 snRNP binding at the downstream 5′ ss D4 and the E42 sequence located between these elements. Therefore, a network of cross-exon interactions appears to regulate splicing of the alternative exons 4a and 5. As the GAR ESE-mediated activation of the upstream 3′ ss cluster also is essential for the processing of intron-containing vpu/env-mRNAs during intermediate viral gene expression, the GAR enhancer substantially contributes to the regulation of viral replication.     

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.     

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.     

The Emergence of Alternative 3′ and 5′ Splice Site Exons from Constitutive Exons

Alternative 3′ and 5′ splice site (ss) events constitute a significant part of all alternative splicing events. These events were also found to be related to several aberrant splicing diseases. However, only few of the characteristics that distinguish these events from alternative cassette exons are known currently. In this study, we compared the characteristics of constitutive exons, alternative cassette exons, and alternative 3′ss and 5′ss exons. The results revealed that alternative 3′ss and 5′ss exons are an intermediate state between constitutive and alternative cassette exons, where the constitutive side resembles constitutive exons, and the alternative side resembles alternative cassette exons. The results also show that alternative 3′ss and 5′ss exons exhibit low levels of symmetry (frame-preserving), similar to constitutive exons, whereas the sequence between the two alternative splice sites shows high symmetry levels, similar to alternative cassette exons. In addition, flanking intronic conservation analysis revealed that exons whose alternative splice sites are at least nine nucleotides apart show a high conservation level, indicating intronic participation in the regulation of their splicing, whereas exons whose alternative splice sites are fewer than nine nucleotides apart show a low conservation level. Further examination of these exons, spanning seven vertebrate species, suggests an evolutionary model in which the alternative state is a derivative of an ancestral constitutive exon, where a mutation inside the exon or along the flanking intron resulted in the creation of a new splice site that competes with the original one, leading to alternative splice site selection. This model was validated experimentally on four exons, showing that they indeed originated from constitutive exons that acquired a new competing splice site during evolution.     

Function of a Bovine Papillomavirus Type 1 Exonic Splicing Suppressor Requires a Suboptimal Upstream 3′ Splice Site

Alternative splicing is an important mechanism for the regulation of bovine papillomavirus type 1 (BPV-1) gene expression during the virus life cycle. Previous studies in our laboratory have identified two purine-rich exonic splicing enhancers (ESEs), SE1 and SE2, located between two alternative 3′ splice sites at nucleotide (nt) 3225 and nt 3605. Further analysis of BPV-1 late-pre-mRNA splicing in vitro revealed a 48-nt pyrimidine-rich region immediately downstream of SE1 that inhibits utilization of the nt 3225 3′ splice site. This inhibitory element, which we named an exonic splicing suppressor (ESS), has a U-rich 5′ end, a C-rich central part, and an AG-rich 3′ end (Z. M. Zheng, P. He, and C. C. Baker, J. Virol. 70:4691–4699, 1996). The present study utilized in vitro splicing of both homologous and heterologous pre-mRNAs to further characterize the ESS. The BPV-1 ESS was inserted downstream of the 3′ splice site in the BPV-1 late pre-mRNA, Rous sarcoma virus src pre-mRNA, human immunodeficiency virus tat-rev pre-mRNA, and Drosophila dsx pre-mRNA, all containing a suboptimal 3′ splice site, and in the human β-globin pre-mRNA, which contains a constitutive 3′ splice site. These studies demonstrated that suppression of splicing by the BPV-1 ESS requires an upstream suboptimal 3′ splice site but not an upstream ESE. Furthermore, the ESS functions when located either upstream or downstream of BPV-1 SE1. Mutational analyses demonstrated that the function of the ESS is sequence dependent and that only the C-rich region of the ESS is essential for suppression of splicing in all the pre-mRNAs tested.     

Intron Definition and a Branch Site Adenosine at nt 385 Control RNA Splicing of HPV16 E6*I and E7 Expression

HPV16 E6 and E7, two viral oncogenes, are expressed from a single bicistronic pre-mRNA. In this report, we provide the evidence that the bicistronic pre-mRNA intron 1 contains three 5′ splice sites (5′ ss) and three 3′ splice sites (3′ ss) normally used in HPV16⁺ cervical cancer and its derived cell lines. The choice of two novel alternative 5′ ss (nt 221 5′ ss and nt 191 5′ ss) produces two novel isoforms of E6E7 mRNAs (E6*V and E6*VI). The nt 226 5′ ss and nt 409 3′ ss is preferentially selected over the other splice sites crossing over the intron to excise a minimal length of the intron in RNA splicing. We identified AACAAAC as the preferred branch point sequence (BPS) and an adenosine at nt 385 (underlined) in the BPS as a branch site to dictate the selection of the nt 409 3′ ss for E6*I splicing and E7 expression. Introduction of point mutations into the mapped BPS led to reduced U2 binding to the BPS and thereby inhibition of the second step of E6E7 splicing at the nt 409 3′ ss. Importantly, the E6E7 bicistronic RNA with a mutant BPS and inefficient splicing makes little or no E7 and the resulted E6 with mutations of ⁹¹QYNK⁹⁴ to ⁹¹PSFW⁹⁴ displays attenuate activity on p53 degradation. Together, our data provide structural basis of the E6E7 intron 1 for better understanding of how viral E6 and E7 expression is regulated by alternative RNA splicing. This study elucidates for the first time a mapped branch point in HPV16 genome involved in viral oncogene expression.