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.     

A genetic screen in C. elegans reveals roles for KIN17 and PRCC in maintaining 5’ splice site identity

Pre-mRNA splicing is an essential step of eukaryotic gene expression carried out by a series of dynamic macromolecular protein/RNA complexes, known collectively and individually as the spliceosome. This series of spliceosomal complexes define, assemble on, and catalyze the removal of introns. Molecular model snapshots of intermediates in the process have been created from cryo-EM data, however, many aspects of the dynamic changes that occur in the spliceosome are not fully understood. Caenorhabditis elegans follow the GU-AG rule of splicing, with almost all introns beginning with 5’ GU and ending with 3’ AG. These splice sites are identified early in the splicing cycle, but as the cycle progresses and “custody” of the pre-mRNA splice sites is passed from factor to factor as the catalytic site is built, the mechanism by which splice site identity is maintained or re-established through these dynamic changes is unclear. We performed a genetic screen in C. elegans for factors that are capable of changing 5’ splice site choice. We report that KIN17 and PRCC are involved in splice site choice, the first functional splicing role proposed for either of these proteins. Previously identified suppressors of cryptic 5’ splicing promote distal cryptic GU splice sites, however, mutations in KIN17 and PRCC instead promote usage of an unusual proximal 5’ splice site which defines an intron beginning with UU, separated by 1nt from a GU donor. We performed high-throughput mRNA sequencing analysis and found that mutations in PRCC, and to a lesser extent KIN17, changed alternative 5’ splice site usage at native sites genome-wide, often promoting usage of nearby non-consensus sites. Our work has uncovered both fine and coarse mechanisms by which the spliceosome maintains splice site identity during the complex assembly process.     

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.     

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.     

Evolutionarily emerged G tracts between the polypyrimidine tract and 3′ AG are splicing silencers enriched in genes involved in cancer

Background

The 3′ splice site (SS) at the end of pre-mRNA introns has a consensus sequence (Y)_nNYAG for constitutive splicing of mammalian genes. Deviation from this consensus could change or interrupt the usage of the splice site leading to alternative or aberrant splicing, which could affect normal cell function or even the development of diseases. We have shown that the position “N” can be replaced by a CA-rich RNA element called CaRRE1 to regulate the alternative splicing of a group of genes.

Results

Taking it a step further, we searched the human genome for purine-rich elements between the -3 and -10 positions of the 3′ splice sites of annotated introns. This identified several thousand such 3′SS; more than a thousand of them contain at least one copy of G tract. These sites deviate significantly from the consensus of constitutive splice sites and are highly associated with alterative splicing events, particularly alternative 3′ splice and intron retention. We show by mutagenesis analysis and RNA interference that the G tracts are splicing silencers and a group of the associated exons are controlled by the G tract binding proteins hnRNP H/F. Species comparison of a group of the 3′SS among vertebrates suggests that most (~87%) of the G tracts emerged in ancestors of mammals during evolution. Moreover, the host genes are most significantly associated with cancer.

Conclusion

We call these elements together with CaRRE1 regulatory RNA elements between the Py and 3′AG (REPA). The emergence of REPA in this highly constrained region indicates that this location has been remarkably permissive for the emergence of de novo regulatory RNA elements, even purine-rich motifs, in a large group of mammalian genes during evolution. This evolutionary change controls alternative splicing, likely to diversify proteomes for particular cellular functions.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1143) contains supplementary material, which is available to authorized users.     

Prp4 Kinase Grants the License to Splice: Control of Weak Splice Sites during Spliceosome Activation

The genome of the fission yeast Schizosaccharomyces pombe encodes 17 kinases that are essential for cell growth. These include the cell-cycle regulator Cdc2, as well as several kinases that coordinate cell growth, polarity, and morphogenesis during the cell cycle. In this study, we further characterized another of these essential kinases, Prp4, and showed that the splicing of many introns is dependent on Prp4 kinase activity. For detailed characterization, we chose the genes res1 and ppk8, each of which contains one intron of typical size and position. Splicing of the res1 intron was dependent on Prp4 kinase activity, whereas splicing of the ppk8 intron was not. Extensive mutational analyses of the 5’ splice site of both genes revealed that proper transient interaction with the 5’ end of snRNA U1 governs the dependence of splicing on Prp4 kinase activity. Proper transient interaction between the branch sequence and snRNA U2 was also important. Therefore, the Prp4 kinase is required for recognition and efficient splicing of introns displaying weak exon1/5’ splice sites and weak branch sequences.     

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.     

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

Despite a growing number of splicing mutations found in hereditary diseases, utilization of aberrant splice sites and their effects on gene expression remain challenging to predict. We compiled sequences of 346 aberrant 5′splice sites (5′ss) that were activated by mutations in 166 human disease genes. Mutations within the 5′ss consensus accounted for 254 cryptic 5′ss and mutations elsewhere activated 92 de novo 5′ss. Point mutations leading to cryptic 5′ss activation were most common in the first intron nucleotide, followed by the fifth nucleotide. Substitutions at position +5 were exclusively G>A transitions, which was largely attributable to high mutability rates of C/G>T/A. However, the frequency of point mutations at position +5 was significantly higher than that observed in the Human Gene Mutation Database, suggesting that alterations of this position are particularly prone to aberrant splicing, possibly due to a requirement for sequential interactions with U1 and U6 snRNAs. Cryptic 5′ss were best predicted by computational algorithms that accommodate nucleotide dependencies and not by weight-matrix models. Discrimination of intronic 5′ss from their authentic counterparts was less effective than for exonic sites, as the former were intrinsically stronger than the latter. Computational prediction of exonic de novo 5′ss was poor, suggesting that their activation critically depends on exonic splicing enhancers or silencers. The authentic counterparts of aberrant 5′ss were significantly weaker than the average human 5′ss. The development of an online database of aberrant 5′ss will be useful for studying basic mechanisms of splice-site selection, identifying splicing mutations and optimizing splice-site prediction algorithms.     

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.     

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.     

Cryptic splicing sites are differentially utilized in vivo

It has long been considered that cryptic splice sites are ignored by the splicing machinery in the context of intact genuine splice sites. In the present study, it is shown that cryptic splice sites are utilized in all circumstances, when the authentic site is intact, partially functional or completely abolished. Their use would therefore contribute to a background lack of fidelity in the context of the wild-type sequence. We also found that a mutation at the 5' splice site of beta-globin intron 1 accommodates multiple cryptic splicing pathways, including three previously reported pathways. Focusing on the two major cryptic 5' splice sites within beta-globin exon 1, we show that cryptic splice site selection ex vivo varies depending upon: (a) the cell stage of development during terminal erythroid differentiation; (b) the nature of the mutation at the authentic 5' splice site; and (c) the nature of the promoter. Finally, we found that the two major cryptic 5' splice sites are utilized with differential efficiencies in two siblings sharing the same beta-globin chromosome haplotype in the homozygous state. Collectively, these data suggest that intrinsic, sequence specific factors and cell genetic background factors both contribute to promote a subtle differential use of cryptic splice sites in vivo.     

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.