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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

The Pivotal Roles of TIA Proteins in 5′ Splice-Site Selection of Alu Exons and Across Evolution

More than 5% of alternatively spliced internal exons in the human genome are derived from Alu elements in a process termed exonization. Alus are comprised of two homologous arms separated by an internal polypyrimidine tract (PPT). In most exonizations, splice sites are selected from within the same arm. We hypothesized that the internal PPT may prevent selection of a splice site further downstream. Here, we demonstrate that this PPT enhanced the selection of an upstream 5′ splice site (5′ss), even in the presence of a stronger 5′ss downstream. Deletion of this PPT shifted selection to the stronger downstream 5′ss. This enhancing effect depended on the strength of the downstream 5′ss, on the efficiency of base-pairing to U1 snRNA, and on the length of the PPT. This effect of the PPT was mediated by the binding of TIA proteins and was dependent on the distance between the PPT and the upstream 5′ss. A wide-scale evolutionary analysis of introns across 22 eukaryotes revealed an enrichment in PPTs within ∼20 nt downstream of the 5′ss. For most metazoans, the strength of the 5′ss inversely correlated with the presence of a downstream PPT, indicative of the functional role of the PPT. Finally, we found that the proteins that mediate this effect, TIA and U1C, and in particular their functional domains, are highly conserved across evolution. Overall, these findings expand our understanding of the role of TIA1/TIAR proteins in enhancing recognition of exons, in general, and Alu exons, in particular.     

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.     

Vitamin A Metabolite,All-trans-retinoic Acid,Mediates Alternative Splicing of Protein Kinase C δVIII (PKCδVIII) Isoform via Splicing Factor SC35

Vitamin A metabolite, all-trans-retinoic acid (RA), induces cell growth, differentiation, and apoptosis and has an emerging role in gene regulation and alternative splicing events. Protein kinase Cδ (PKCδ), a serine/threonine kinase, has a role in cell proliferation, differentiation, and apoptosis. We reported an alternatively spliced variant of human PKCδ, PKCδVIII that functions as a pro-survival protein (1). RA regulates the splicing and expression of PKCδVIII via utilization of a downstream 5′ splice site of exon 10 on PKCδ pre-mRNA. Here, we further elucidate the molecular mechanisms involved in RA regulation of alternative splicing of PKCδVIII mRNA. Overexpression and knockdown of the splicing factor SC35 (i.e. SRp30b) indicated that it is involved in PKCδVIII alternative splicing. To identify the cis-elements involved in 5′ splice site selection we cloned a minigene, which included PKCδ exon 10 and its flanking introns in the pSPL3 splicing vector. Alternative 5′ splice site utilization in the minigene was promoted by RA. Further, co-transfection of SC35 with PKCδ minigene promoted selection of 5′ splice site II. Mutation of the SC35 binding site in the PKCδ minigene abolished RA-mediated utilization of 5′ splice splice II. RNA binding assays demonstrated that the enhancer element downstream of PKCδ exon 10 is a SC35 cis-element. We conclude that SC35 is pivotal in RA-mediated PKCδ pre-mRNA alternative splicing. This study demonstrates how a nutrient, vitamin A, via its metabolite RA, regulates alternative splicing and thereby gene expression of the pro-survival protein PKCδVIII.     

Cap-independent translation through the p27 5'-UTR

Several recent publications have explored cap-independent translation through an internal ribosome entry site (IRES) in the 5′-UTR of the mRNA encoding the cyclin-dependent kinase inhibitor p27. The major experimental tool used in these reports was the use of bicistronic reporter constructs in which the 5′-UTR was inserted between the upstream and downstream cistrons. None of these reports has completely ruled out the possibility that the 5′-UTR has either cryptic promoter activity or a cryptic splice acceptor site. Either of these possibilities could result in expression of a monocistronic mRNA encoding the downstream cistron and false identification of an IRES. Indeed, Liu et al. recently published data suggesting that the p27 5′-UTR harbors cryptic promoter activity which accounts for its putative IRES activity. In this report, we have explored this potential problem further using promoterless bicistronic constructs coupled with RNase protection assays, siRNA knockdown of individual cistrons, RT-PCR to detect mRNA encoded by the bicistronic reporter with high sensitivity, direct transfection of bicistronic mRNAs, and insertion of an iron response element into the bicistronic reporter. The results do not support the conclusion that the p27 5′-UTR has significant functional promoter activity or cryptic splice sites, but rather that it is able to support cap-independent initiation of translation.     

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.     

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.     

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.     

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.     

Regulation of BCL-X splicing reveals a role for the polypyrimidine tract binding protein (PTBP1/hnRNP I) in alternative 5′ splice site selection

Alternative splicing (AS) modulates many physiological and pathological processes. For instance, AS of the BCL-X gene balances cell survival and apoptosis in development and cancer. Herein, we identified the polypyrimidine tract binding protein (PTBP1) as a direct regulator of BCL-X AS. Overexpression of PTBP1 promotes selection of the distal 5′ splice site in BCL-X exon 2, generating the pro-apoptotic BCL-Xs splice variant. Conversely, depletion of PTBP1 enhanced splicing of the anti-apoptotic BCL-XL variant. In vivo cross-linking experiments and site-directed mutagenesis restricted the PTBP1 binding site to a polypyrimidine tract located between the two alternative 5′ splice sites. Binding of PTBP1 to this site was required for its effect on splicing. Notably, a similar function of PTBP1 in the selection of alternative 5′ splice sites was confirmed using the USP5 gene as additional model. Mechanistically, PTBP1 displaces SRSF1 binding from the proximal 5′ splice site, thus repressing its selection. Our study provides a novel mechanism of alternative 5′ splice site selection by PTBP1 and indicates that the presence of a PTBP1 binding site between two alternative 5′ splice sites promotes selection of the distal one, while repressing the proximal site by competing for binding of a positive regulator.     

A Genetic Screen for Suppressors of a Mutated 5′ Splice Site Identifies Factors Associated With Later Steps of Spliceosome Assembly

Many alleles of human disease genes have mutations within splicing consensus sequences that activate cryptic splice sites. In Caenorhabditis elegans, the unc-73(e936) allele has a G-to-U mutation at the first base of the intron downstream of exon 15, which results in an uncoordinated phenotype. This mutation triggers cryptic splicing at the −1 and +23 positions and retains some residual splicing at the mutated wild-type (wt) position. We previously demonstrated that a mutation in sup-39, a U1 snRNA gene, suppresses e936 by increasing splicing at the wt splice site. We report here the results of a suppressor screen in which we identify three proteins that function in cryptic splice site choice. Loss-of-function mutations in the nonessential splicing factor smu-2 suppress e936 uncoordination through changes in splicing. SMU-2 binds SMU-1, and smu-1(RNAi) also leads to suppression of e936. A dominant mutation in the conserved C-terminal domain of the C. elegans homolog of the human tri-snRNP 27K protein, which we have named SNRP-27, suppresses e936 uncoordination through changes in splicing. We propose that SMU-2, SMU-1, and SNRP-27 contribute to the fidelity of splice site choice after the initial identification of 5′ splice sites by U1 snRNP.PRE-mRNA splicing takes place in a large ribonucleoprotein complex called the spliceosome (Burge et al. 1999). Components of this splicing machinery assemble at conserved signal sequences within the pre-mRNA. The 5′ splice site consensus sequence M₋₃A₋₂G₋₁ | G₊₁U₊₂R₊₃A₊₄G₊₅U₊₆ and the 3′ splice site consensus sequence Y₋₃A₋₂G₋₁ | R₊₁ (M is either A or C; R is a purine, and Y is a pyrimidine) define the limits of the intron. Base-pairing interactions between the 5′ end of the U1 snRNA and the 5′ splice site consensus sequence occur early in spliceosome assembly. It is the nearly invariable GU dinucleotide at the first two positions of the 5′ end of the intron that defines the beginning of the intron. The 5′ consensus sequence is essential but insufficient for splice site selection, as 5′ splice sites with weaker consensus matches may require additional determinants for proper activation (Sanford et al. 2005).Mutations that disrupt the 5′ consensus splice signal can lead to genetic disease in humans (Nelson and Green 1990; Cohen et al. 1994). Approximately 15% of point mutations that cause genetic diseases affect pre-mRNA splicing consensus sequences (Krawczak et al. 1992). For some specific disease genes, as many as 50% of the known heritable alleles alter splicing (Teraoka et al. 1999; Ars et al. 2000; Roca et al. 2003; Pagenstecher et al. 2006). Among all the positions of the 5′ splice site consensus sequence, the highest proportion of human disease mutations occur at the +1G position (Buratti et al. 2007). The fidelity of pre-mRNA splice site choice is largely disrupted by this defect, since this mutation causes splicing at this site to be either abolished or outcompeted by the activation of nearby cryptic 5′ splice sites (Nelson and Green 1990; Cohen et al. 1994). Cryptic splice sites are used only when the wild-type splice donor is disrupted by mutation, as they tend to have very weak splice donor consensus sequences outside of a 5′-GU dinucleotide that defines the beginning of the intron (Roca et al. 2003). Suppression of mutations to the 5′ splice site consensus sequence in vivo has been achieved through the expression of U1 snRNAs containing compensatory base substitutions (Zhuang and Weiner 1986); however, suppression of mutations to the +1 position of the intron using reverse genetic approaches has not been successful (Newman et al. 1985; Nelson and Green 1990; Cohen et al. 1994).We have used a specific allele of the Caenorhabditis elegans unc-73 gene, e936, which contains a G-to-U mutation at the first nucleotide of intron 16 (Steven et al. 1998), as a model for studying cryptic splice site choice (Roller et al. 2000; Zahler et al. 2004). unc-73 encodes a RAC guanine nucleotide exchange factor that is expressed in neurons and is important for axon guidance (Steven et al. 1998). The e936 allele induces the use of three different cryptic 5′ splice sites (Figure 1A). Two of these 5′ splice sites, located at the −1 and +23 positions, define introns beginning with GU. The third 5′ splice site used is at the mutated wild-type (wt) position and is referred to as “wt” since splicing at this site still produces wild-type unc-73 mRNA and protein, even though the intron begins with UU (Roller et al. 2000). Use of either the −1 or the +23 cryptic site causes a shift in the reading frame and loss of gene function. In e936 animals, 90% of the stable messages of unc-73 are out-of-frame, yet the phenotype is not as severe as for other alleles in this gene. This indicates that the 10% of steady-state messages that are in frame have some functional role.Open in a separate windowFigure 1.—(A) Diagram of the unc-73 gene between exons 15 and 16. The positions of the −1 and +23 cryptic 5′ splice sites are indicated by arrows. The intronic e936 (+1G → U) point mutation is highlighted. (B) γ-³²P-labeled RT–PCR results across the cryptic splicing region of unc-73(e936) for different strains. Lanes 1, 2, and 3 are loaded with RT–PCR reactions from wild type (N2), unc-73(e936);sup-39(je5), and unc-73(e936) RNA, respectively. The lines carrying the suppressor alleles and e936 follow in lanes 4–10 as indicated. (C) The unc-73 genomic sequence from exon 15 (uppercase letters) and intron 15 (lowercase letters). The locations of the az23 and e936 mutational substitutions are indicated below. The position of the −9 cryptic splice donor activated in e936az23 is indicated by an arrow above.In a previous genetic screen for extragenic suppressors of e936 movement defects, Way and colleagues identified sup-39 (Run et al. 1996). It was subsequently shown that mutations in sup-39 alter cryptic splice site choice of e936 (Roller et al. 2000). sup-39 encodes a U1 snRNA gene with a compensatory mutation at the position that normally base pairs with the +1G. This allows sup-39 to base pair with an intron with a +1U (Zahler et al. 2004). This dominant suppressor increases usage of the mutated splice site and improves the fraction of in-frame messages from e936 from 10 to 33%, with a dramatic improvement in coordination. A similar mutant U1 snRNA suppressor with a different compensatory substitution, sup-6(st19), was found to suppress the intronic +1G to A transition of unc-13(e309) to allow for splicing at the mutated wild-type site, even though the intron begins with AU instead of GU (Zahler et al. 2004).We are interested in identifying additional factors that play a role in cryptic 5′ splice site choice. To do this, we took advantage of unc-73(e936), in which modest increases in the use of the wt splice site lead to dramatic increases in coordination, as a sensitive screen for changes in cryptic splice site choice. In this article we report that the proteins SMU-1 and SMU-2, which are nonessential factors previously shown to have a role in alternative splicing (Spartz et al. 2004), have a role in selection of cryptic 5′ splice sites. We also report the identification of a new dominant suppressor of cryptic splicing, snrp-27, which encodes a C. elegans homolog of the human tri-snRNP 27K protein.