Discovery and analysis of evolutionarily conserved intronic splicing regulatory elements

Knowledge of the functional cis-regulatory elements that regulate constitutive and alternative pre-mRNA splicing is fundamental for biology and medicine. Here we undertook a genome-wide comparative genomics approach using available mammalian genomes to identify conserved intronic splicing regulatory elements (ISREs). Our approach yielded 314 ISREs, and insertions of ~70 ISREs between competing splice sites demonstrated that 84% of ISREs altered 5′ and 94% altered 3′ splice site choice in human cells. Consistent with our experiments, comparisons of ISREs to known splicing regulatory elements revealed that 40%–45% of ISREs might have dual roles as exonic splicing silencers. Supporting a role for ISREs in alternative splicing, we found that 30%–50% of ISREs were enriched near alternatively spliced (AS) exons, and included almost all known binding sites of tissue-specific alternative splicing factors. Further, we observed that genes harboring ISRE-proximal exons have biases for tissue expression and molecular functions that are ISRE-specific. Finally, we discovered that for Nova1, neuronal PTB, hnRNP C, and FOX1, the most frequently occurring ISRE proximal to an alternative conserved exon in the splicing factor strongly resembled its own known RNA binding site, suggesting a novel application of ISRE density and the propensity for splicing factors to auto-regulate to associate RNA binding sites to splicing factors. Our results demonstrate that ISREs are crucial building blocks in understanding general and tissue-specific AS regulation and the biological pathways and functions regulated by these AS events.     

Conserved sequence elements associated with exon skipping

One of the major forms of alternative splicing, which generates multiple mRNA isoforms differing in the precise combinations of their exon sequences, is exon skipping. While in constitutive splicing all exons are included, in the skipped pattern(s) one or more exons are skipped. The regulation of this process is still not well understood; so far, cis- regulatory elements (such as exonic splicing enhancers) were identified in individual cases. We therefore set to investigate the possibility that exon skipping is controlled by sequences in the adjacent introns. We employed a computer analysis on 54 sequences documented as undergoing exon skipping, and identified two motifs both in the upstream and downstream introns of the skipped exons. One motif is highly enriched in pyrimidines (mostly C residues), and the other motif is highly enriched in purines (mostly G residues). The two motifs differ from the known cis-elements present at the 5′ and 3′ splice site. Interestingly, the two motifs are complementary, and their relative positional order is conserved in the flanking introns. These suggest that base pairing interactions can underlie a mechanism that involves secondary structure to regulate exon skipping. Remarkably, the two motifs are conserved in mouse orthologous genes that undergo exon skipping.     

Wilms’ tumor 1-associating protein complex regulates alternative splicing and polyadenylation at potential G-quadruplex-forming splice site sequences

Wilms’ tumor 1-associating protein (WTAP) is a core component of the N6-methyladenosine (m6A)-methyltransferase complex, along with VIRMA, CBLL1, ZC3H13 (KIAA0853), RBM15/15B, and METTL3/14, which generate m6A, a key RNA modification that affects various processes of RNA metabolism. WTAP also interacts with splicing factors; however, despite strong evidence suggesting a role of Drosophila WTAP homolog fl(2)d in alternative splicing (AS), its role in splicing regulation in mammalian cells remains elusive. Here we demonstrate using RNAi coupled with RNA-seq that WTAP, VIRMA, CBLL1, and ZC3H13 modulate AS, promoting exon skipping and intron retention in AS events that involve short introns/exons with higher GC content and introns with weaker polypyrimidine-tract and branch points. Further analysis of GC-rich sequences involved in AS events regulated by WTAP, together with minigene assay analysis, revealed potential G-quadruplex formation at splice sites where WTAP has an inhibitory effect. We also found that several AS events occur in the last exon of one isoform of MSL1 and WTAP, leading to competition for polyadenylation. Proteomic analysis also suggested that WTAP/CBLL1 interaction promotes recruitment of the 3′-end processing complex. Taken together, our results indicate that the WTAP complex regulates AS and alternative polyadenylation via inhibitory mechanisms in GC-rich sequences.     

A Multi-Exon-Skipping Detection Assay Reveals Surprising Diversity of Splice Isoforms of Spinal Muscular Atrophy Genes

Humans have two near identical copies of Survival Motor Neuron gene: SMN1 and SMN2. Loss of SMN1 coupled with the predominant skipping of SMN2 exon 7 causes spinal muscular atrophy (SMA), a neurodegenerative disease. SMA patient cells devoid of SMN1 provide a powerful system to examine splicing pattern of various SMN2 exons. Until now, similar system to examine splicing of SMN1 exons was unavailable. We have recently screened several patient cell lines derived from various diseases, including SMA, Alzheimer’s disease, Parkinson’s disease and Batten disease. Here we report a Batten disease cell line that lacks functional SMN2, as an ideal system to examine pre-mRNA splicing of SMN1. We employ a multiple-exon-skipping detection assay (MESDA) to capture simultaneously skipping of multiple exons. Our results show surprising diversity of splice isoforms and reveal novel splicing events that include skipping of exon 4 and co-skipping of three adjacent exons of SMN. Contrary to the general belief, MESDA captured oxidative-stress induced skipping of SMN1 exon 5 in several cell types, including non-neuronal cells. We further demonstrate that the predominant SMN2 exon 7 skipping induced by oxidative stress is modulated by a combinatorial control that includes promoter sequence, endogenous context, and the weak splice sites. We also show that an 8-mer antisense oligonucleotide blocking a recently described GC-rich sequence prevents SMN2 exon 7 skipping under the conditions of oxidative stress. Our findings bring new insight into splicing regulation of an essential housekeeping gene linked to neurodegeneration and infant mortality.     

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.     

Pyrvinium pamoate changes alternative splicing of the serotonin receptor 2C by influencing its RNA structure

The serotonin receptor 2C plays a central role in mood and appetite control. It undergoes pre-mRNA editing as well as alternative splicing. The RNA editing suggests that the pre-mRNA forms a stable secondary structure in vivo. To identify substances that promote alternative exons inclusion, we set up a high-throughput screen and identified pyrvinium pamoate as a drug-promoting exon inclusion without editing. Circular dichroism spectroscopy indicates that pyrvinium pamoate binds directly to the pre-mRNA and changes its structure. SHAPE (selective 2′-hydroxyl acylation analysed by primer extension) assays show that part of the regulated 5′-splice site forms intramolecular base pairs that are removed by this structural change, which likely allows splice site recognition and exon inclusion. Genome-wide analyses show that pyrvinium pamoate regulates >300 alternative exons that form secondary structures enriched in A–U base pairs. Our data demonstrate that alternative splicing of structured pre-mRNAs can be regulated by small molecules that directly bind to the RNA, which is reminiscent to an RNA riboswitch.     

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.     

Stabilization of the μ-Opioid Receptor by Truncated Single Transmembrane Splice Variants through a Chaperone-like Action

The μ-opioid receptor gene, OPRM1, undergoes extensive alternative pre-mRNA splicing, as illustrated by the identification of an array of splice variants generated by both 5′ and 3′ alternative splicing. The current study reports the identification of another set of splice variants conserved across species that are generated through exon skipping or insertion that encodes proteins containing only a single transmembrane (TM) domain. Using a Tet-Off system, we demonstrated that the truncated single TM variants can dimerize with the full-length 7-TM μ-opioid receptor (MOR-1) in the endoplasmic reticulum, leading to increased expression of MOR-1 at the protein level by a chaperone-like function that minimizes endoplasmic reticulum-associated degradation. In vivo antisense studies suggested that the single TM variants play an important role in morphine analgesia, presumably through modulation of receptor expression levels. Our studies suggest the functional roles of truncated receptors in other G protein-coupled receptor families.     

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.     

Rbm20 regulates titin alternative splicing as a splicing repressor

Titin, a sarcomeric protein expressed primarily in striated muscles, is responsible for maintaining the structure and biomechanical properties of muscle cells. Cardiac titin undergoes developmental size reduction from 3.7 megadaltons in neonates to primarily 2.97 megadaltons in the adult. This size reduction results from gradually increased exon skipping between exons 50 and 219 of titin mRNA. Our previous study reported that Rbm20 is the splicing factor responsible for this process. In this work, we investigated its molecular mechanism. We demonstrate that Rbm20 mediates exon skipping by binding to titin pre-mRNA to repress the splicing of some regions; the exons/introns in these Rbm20-repressed regions are ultimately skipped. Rbm20 was also found to mediate intron retention and exon shuffling. The two Rbm20 speckles found in nuclei from muscle tissues were identified as aggregates of Rbm20 protein on the partially processed titin pre-mRNAs. Cooperative repression and alternative 3′ splice site selection were found to be used by Rbm20 to skip different subsets of titin exons, and the splicing pathway selected depended on the ratio of Rbm20 to other splicing factors that vary with tissue type and developmental age.     

Modulation of exon skipping and inclusion by heterogeneous nuclear ribonucleoprotein A1 and pre-mRNA splicing factor SF2/ASF.

The essential splicing factor SF2/ASF and the heterogeneous nuclear ribonucleoprotein A1 (hnRNP A1) modulate alternative splicing in vitro of pre-mRNAs that contain 5' splice sites of comparable strengths competing for a common 3' splice site. Using natural and model pre-mRNAs, we have examined whether the ratio of SF2/ASF to hnRNP A1 also regulates other modes of alternative splicing in vitro. We found that an excess of SF2/ASF effectively prevents inappropriate exon skipping and also influences the selection of mutually exclusive tissue-specific exons in natural beta-tropomyosin pre-mRNA. In contrast, an excess of hnRNP A1 does not cause inappropriate exon skipping in natural constitutively or alternatively spliced pre-mRNAs. Although hnRNP A1 can promote alternative exon skipping, this effect is not universal and is dependent, e.g., on the size of the internal alternative exon and on the strength of the polypyrimidine tract in the preceding intron. With appropriate alternative exons, an excess of SF2/ASF promotes exon inclusion, whereas an excess of hnRNP A1 causes exon skipping. We propose that in some cases the ratio of SF2/ASF to hnRNP A1 may play a role in regulating alternative splicing by exon inclusion or skipping through the antagonistic effects of these proteins on alternative 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.     

Characterization of exon skipping mutants of the COP1 gene from Arabidopsis

The removal of introns from pre-mRNA requires accurate recognition and selection of the intron splice sites. Mutations which alter splice site selection and which lead to skipping of specific exons are indicative of intron/exon recognition mechanisms involving an exon definition process. In this paper, three independent mutants to the COP1 gene in Arabidopsis which show exon skipping were identified and the mutations which alter the normal splicing pattern were characterized. The mutation in cop1–1 was a G→A change 4 nt upstream from the 3′ splice site of intron 5, while the mutation in cop1–2 was a G→A at the first nucleotide of intron 6, abolishing the conserved G within the 5′ splice site consensus. The effect of these mutations was skipping of exon 6. The mutation in cop1–8 was G→A in the final nucleotide of intron 10 abolishing the conserved G within the 3′ splice site consensus and leading to skipping of exon 11. The splicing patterns surrounding exons 6 and 11 of COP1 in these three mutant lines of Arabidopsis provide evidence for exon definition mechanisms operating in plant splicing.     

Characterization of RNase R-digested cellular RNA source that consists of lariat and circular RNAs from pre-mRNA splicing

Besides linear RNAs, pre-mRNA splicing generates three forms of RNAs: lariat introns, Y-structure introns from trans-splicing, and circular exons through exon skipping. To study the persistence of excised introns in total cellular RNA, we used three Escherichia coli 3′ to 5′ exoribonucleases. Ribonuclease R (RNase R) thoroughly degrades the abundant linear RNAs and the Y-structure RNA, while preserving the loop portion of a lariat RNA. Ribonuclease II (RNase II) and polynucleotide phosphorylase (PNPase) also preserve the lariat loop, but are less efficient in degrading linear RNAs. RNase R digestion of the total RNA from human skeletal muscle generates an RNA pool consisting of lariat and circular RNAs. RT–PCR across the branch sites confirmed lariat RNAs and circular RNAs in the pool generated by constitutive and alternative splicing of the dystrophin pre-mRNA. Our results indicate that RNase R treatment can be used to construct an intronic cDNA library, in which majority of the intron lariats are represented. The highly specific activity of RNase R implies its ability to screen for rare intragenic trans-splicing in any target gene with a large background of cis-splicing. Further analysis of the intronic RNA pool from a specific tissue or cell will provide insights into the global profile of alternative splicing.