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.     

Functional selection and systematic analysis of intronic splicing elements identify active sequence motifs and associated splicing factors

Despite the critical role of pre-mRNA splicing in generating proteomic diversity and regulating gene expression, the sequence composition and function of intronic splicing regulatory elements (ISREs) have not been well elucidated. Here, we employed a high-throughput in vivo Screening PLatform for Intronic Control Elements (SPLICE) to identify 125 unique ISRE sequences from a random nucleotide library in human cells. Bioinformatic analyses reveal consensus motifs that resemble splicing regulatory elements and binding sites for characterized splicing factors and that are enriched in the introns of naturally occurring spliced genes, supporting their biological relevance. In vivo characterization, including an RNAi silencing study, demonstrate that ISRE sequences can exhibit combinatorial regulatory activity and that multiple trans-acting factors are involved in the regulatory effect of a single ISRE. Our work provides an initial examination into the sequence characteristics and function of ISREs, providing an important contribution to the splicing code.     

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.     

Does distance matter? Variations in alternative 3′ splicing regulation

Alternative splicing constitutes a major mechanism creating protein diversity in humans. This diversity can result from the alternative skipping of entire exons or by alternative selection of the 5′ or 3′ splice sites that define the exon boundaries. In this study, we analyze the sequence and evolutionary characteristics of alternative 3′ splice sites conserved between human and mouse genomes for distances ranging from 3 to 100 nucleotides. We show that alternative splicing events can be distinguished from constitutive splicing by a combination of properties which vary depending on the distance between the splice sites. Among the unique features of alternative 3′ splice sites, we observed an unexpectedly high occurrence of events in which a polypyrimidine tract was found to overlap the upstream splice site. By applying a machine-learning approach, we show that we can successfully discriminate true alternative 3′ splice sites from constitutive 3′ splice sites. Finally, we propose that the unique features of the intron flanking alternative splice sites are indicative of a regulatory mechanism that is involved in splice site selection. We postulate that the process of splice site selection is influenced by the distance between the competitive splice sites.     

A Novel Intra-U1 snRNP Cross-Regulation Mechanism: Alternative Splicing Switch Links U1C and U1-70K Expression

The U1 small nuclear ribonucleoprotein (snRNP)-specific U1C protein participates in 5′ splice site recognition and regulation of pre-mRNA splicing. Based on an RNA-Seq analysis in HeLa cells after U1C knockdown, we found a conserved, intra-U1 snRNP cross-regulation that links U1C and U1-70K expression through alternative splicing and U1 snRNP assembly. To investigate the underlying regulatory mechanism, we combined mutational minigene analysis, in vivo splice-site blocking by antisense morpholinos, and in vitro binding experiments. Alternative splicing of U1-70K pre-mRNA creates the normal (exons 7–8) and a non-productive mRNA isoform, whose balance is determined by U1C protein levels. The non-productive isoform is generated through a U1C-dependent alternative 3′ splice site, which requires an adjacent cluster of regulatory 5′ splice sites and binding of intact U1 snRNPs. As a result of nonsense-mediated decay (NMD) of the non-productive isoform, U1-70K mRNA and protein levels are down-regulated, and U1C incorporation into the U1 snRNP is impaired. U1-70K/U1C-deficient particles are assembled, shifting the alternative splicing balance back towards productive U1-70K splicing, and restoring assembly of intact U1 snRNPs. Taken together, we established a novel feedback regulation that controls U1-70K/U1C homeostasis and ensures correct U1 snRNP assembly and function.     

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.     

Evolution of Nova-dependent splicing regulation in the brain

A large number of alternative exons are spliced with tissue-specific patterns, but little is known about how such patterns have evolved. Here, we study the conservation of the neuron-specific splicing factors Nova1 and Nova2 and of the alternatively spliced exons they regulate in mouse brain. Whereas Nova RNA binding domains are 94% identical across vertebrate species, Nova-dependent splicing silencer and enhancer elements (YCAY clusters) show much greater divergence, as less than 50% of mouse YCAY clusters are conserved at orthologous positions in the zebrafish genome. To study the relation between the evolution of tissue-specific splicing and YCAY clusters, we compared the brain-specific splicing of Nova-regulated exons in zebrafish, chicken, and mouse. The presence of YCAY clusters in lower vertebrates invariably predicted conservation of brain-specific splicing across species, whereas their absence in lower vertebrates correlated with a loss of alternative splicing. We hypothesize that evolution of Nova-regulated splicing in higher vertebrates proceeds mainly through changes in cis-acting elements, that tissue-specific splicing might in some cases evolve in a single step corresponding to evolution of a YCAY cluster, and that the conservation level of YCAY clusters relates to the functions encoded by the regulated RNAs.     

Alternative splicing regulation at tandem 3' splice sites

Alternative splicing (AS) constitutes a major mechanism creating protein diversity in humans. Previous bioinformatics studies based on expressed sequence tag and mRNA data have identified many AS events that are conserved between humans and mice. Of these events, ~25% are related to alternative choices of 3′ and 5′ splice sites. Surprisingly, half of all these events involve 3′ splice sites that are exactly 3 nt apart. These tandem 3′ splice sites result from the presence of the NAGNAG motif at the acceptor splice site, recently reported to be widely spread in the human genome. Although the NAGNAG motif is common in human genes, only a small subset of sites with this motif is confirmed to be involved in AS. We examined the NAGNAG motifs and observed specific features such as high sequence conservation of the motif, high conservation of ~30 bp at the intronic regions flanking the 3′ splice site and overabundance of cis-regulatory elements, which are characteristic of alternatively spliced tandem acceptor sites and can distinguish them from the constitutive sites in which the proximal NAG splice site is selected. Our findings imply that AS at tandem splice sites and constitutive splicing of the distal NAG are highly regulated.     

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 Alternative Polyadenylation by U1 snRNPs and SRp20

Although considerable information is currently available about the factors involved in constitutive vertebrate polyadenylation, the factors and mechanisms involved in facilitating communication between polyadenylation and splicing are largely unknown. Even less is known about the regulation of polyadenylation in genes in which 3′-terminal exons are alternatively recognized. Here we demonstrate that an SR protein, SRp20, affects recognition of an alternative 3′-terminal exon via an effect on the efficiency of binding of a polyadenylation factor to an alternative polyadenylation site. The gene under study codes for the peptides calcitonin and calcitonin gene-related peptide. Its pre-mRNA is alternatively processed by the tissue-specific inclusion or exclusion of an embedded 3′-terminal exon, exon 4, via factors binding to an intronic enhancer element that contains both 3′ and 5′ splice site consensus sequence elements. In cell types that preferentially exclude exon 4, addition of wild-type SRp20 enhances exon 4 inclusion via recognition of the intronic enhancer. In contrast, in cell types that preferentially include exon 4, addition of a mutant form of SRp20 containing the RNA-binding domain but missing the SR domain inhibits exon 4 inclusion. Inhibition is likely at the level of polyadenylation, because the mutant SRp20 inhibits binding of CstF to the exon 4 poly(A) site. This is the first demonstration that an SR protein can influence alternative polyadenylation and suggests that this family of proteins may play a role in recognition of 3′-terminal exons and perhaps in the communication between polyadenylation and splicing.     

Splicing Reporter Mice Revealed the Evolutionally Conserved Switching Mechanism of Tissue-Specific Alternative Exon Selection

Since alternative splicing of pre-mRNAs is essential for generating tissue-specific diversity in proteome, elucidating its regulatory mechanism is indispensable to understand developmental process or tissue-specific functions. We have been focusing on tissue-specific regulation of mutually exclusive selection of alternative exons because this implies the typical molecular mechanism of alternative splicing regulation and also can be good examples to elicit general rule of “splice code”. So far, mutually exclusive splicing regulation has been explained by the outcome from the balance of multiple regulators that enhance or repress either of alternative exons discretely. However, this “balance” model is open to questions of how to ensure the selection of only one appropriate exon out of several candidates and how to switch them. To answer these questions, we generated an original bichromatic fluorescent splicing reporter system for mammals using fibroblast growth factor-receptor 2 (FGFR2) gene as model. By using this splicing reporter, we demonstrated that FGFR2 gene is regulated by the “switch-like” mechanism, in which key regulators modify the ordered splice-site recognition of two mutually exclusive exons, eventually ensure single exon selection and their distinct switching. Also this finding elucidated the evolutionally conserved “splice code,” in which combination of tissue-specific and broadly expressed RNA binding proteins regulate alternative splicing of specific gene in a tissue-specific manner. These findings provide the significant cue to understand how a number of spliced genes are regulated in various tissue-specific manners by a limited number of regulators, eventually to understand developmental process or tissue-specific functions.     

Unusual intron conservation near tissue-regulated exons found by splicing microarrays

Alternative splicing contributes to both gene regulation and protein diversity. To discover broad relationships between regulation of alternative splicing and sequence conservation, we applied a systems approach, using oligonucleotide microarrays designed to capture splicing information across the mouse genome. In a set of 22 adult tissues, we observe differential expression of RNA containing at least two alternative splice junctions for about 40% of the 6,216 alternative events we could detect. Statistical comparisons identify 171 cassette exons whose inclusion or skipping is different in brain relative to other tissues and another 28 exons whose splicing is different in muscle. A subset of these exons is associated with unusual blocks of intron sequence whose conservation in vertebrates rivals that of protein-coding exons. By focusing on sets of exons with similar regulatory patterns, we have identified new sequence motifs implicated in brain and muscle splicing regulation. Of note is a motif that is strikingly similar to the branchpoint consensus but is located downstream of the 5′ splice site of exons included in muscle. Analysis of three paralogous membrane-associated guanylate kinase genes reveals that each contains a paralogous tissue-regulated exon with a similar tissue inclusion pattern. While the intron sequences flanking these exons remain highly conserved among mammalian orthologs, the paralogous flanking intron sequences have diverged considerably, suggesting unusually complex evolution of the regulation of alternative splicing in multigene families.     

A Short Sequence within Two Purine-Rich Enhancers Determines 5′ Splice Site Specificity

Purine-rich enhancers are exon sequences that promote inclusion of alternative exons, usually via activation of weak upstream 3′ splice sites. A recently described purine-rich enhancer from the caldesmon gene has an additional activity by which it directs selection of competing 5′ splice sites within an alternative exon. In this study, we have compared the caldesmon enhancer with another purine-rich enhancer from the chicken cardiac troponin T (cTNT) gene for the ability to regulate flanking splice sites. Although similar in sequence and length, the two enhancers demonstrated strikingly different specificities towards 5′ splice site choice when placed between competing 5′ splice sites in an internal exon. The 32-nucleotide caldesmon enhancer caused effective usage of the exon-internal 5′ splice site, whereas the 30-nucleotide cTNT enhancer caused effective usage of the exon-terminal 5′ splice site. Both enhancer-mediated splicing pathways represented modulation of the default pathway in which both 5′ splice sites were utilized. Each enhancer is multipartite, consisting of two purine-rich sequences of a simple (GAR)_n repeat interdigitated with two enhancer-specific sequences. The entire enhancer was necessary for maximal splice site selectivity; however, a 5- to 7-nucleotide region from the 3′ end of each enhancer dictated splice site selectivity. Mutations that interchanged this short region of the two enhancers switched specificity. The portion of the cTNT enhancer determinative for 5′ splice site selectivity was different than that shown to be maximally important for activation of a 3′ splice site, suggesting that enhancer environment can have a major impact on activity. These results are the first indication that individual purine-rich enhancers can differentiate between flanking splice sites. Furthermore, localization of the specificity of splice site choice to a short region within both enhancers indicates that subtle differences in enhancer sequence can have profound effects on the splicing pathway.