Evolution of Alternative Splicing Regulation: Changes in Predicted Exonic Splicing Regulators Are Not Associated with Changes in Alternative Splicing Levels in Primates

Alternative splicing is tightly regulated in a spatio-temporal and quantitative manner. This regulation is achieved by a complex interplay between spliceosomal (trans) factors that bind to different sequence (cis) elements. cis-elements reside in both introns and exons and may either enhance or silence splicing. Differential combinations of cis-elements allows for a huge diversity of overall splicing signals, together comprising a complex ‘splicing code’. Many cis-elements have been identified, and their effects on exon inclusion levels demonstrated in reporter systems. However, the impact of interspecific differences in these elements on the evolution of alternative splicing levels has not yet been investigated at genomic level. Here we study the effect of interspecific differences in predicted exonic splicing regulators (ESRs) on exon inclusion levels in human and chimpanzee. For this purpose, we compiled and studied comprehensive datasets of predicted ESRs, identified by several computational and experimental approaches, as well as microarray data for changes in alternative splicing levels between human and chimpanzee. Surprisingly, we found no association between changes in predicted ESRs and changes in alternative splicing levels. This observation holds across different ESR exon positions, exon lengths, and 5′ splice site strengths. We suggest that this lack of association is mainly due to the great importance of context for ESR functionality: many ESR-like motifs in primates may have little or no effect on splicing, and thus interspecific changes at short-time scales may primarily occur in these effectively neutral ESRs. These results underscore the difficulties of using current computational ESR prediction algorithms to identify truly functionally important motifs, and provide a cautionary tale for studies of the effect of SNPs on splicing in human disease.     

Cleavage and polyadenylation specificity factor 1 (CPSF1) regulates alternative splicing of interleukin 7 receptor (IL7R) exon 6

Interleukin 7 receptor, IL7R, is expressed exclusively on cells of the lymphoid lineage, and its expression is crucial for the development and maintenance of T cells. Alternative splicing of IL7R exon 6 results in membrane-bound (exon 6 included) and soluble (exon 6 skipped) IL7R isoforms. Interestingly, the inclusion of exon 6 is affected by a single-nucleotide polymorphism associated with the risk of developing multiple sclerosis. Given the potential association of exon 6 inclusion with multiple sclerosis, we investigated the cis-acting elements and trans-acting factors that regulate exon 6 splicing. We identified multiple exonic and intronic cis-acting elements that impact inclusion of exon 6. Moreover, we utilized RNA affinity chromatography followed by mass spectrometry to identify trans-acting protein factors that bind exon 6 and regulate its splicing. These experiments identified cleavage and polyadenylation specificity factor 1 (CPSF1) among protein-binding candidates. A consensus polyadenylation signal AAUAAA is present in intron 6 of IL7R directly downstream from the 5′ splice site. Mutations to this site and CPSF1 knockdown both resulted in an increase in exon 6 inclusion. We found no evidence that this site is used to produce cleaved and polyadenylated mRNAs, suggesting that CPSF1 interaction with intronic IL7R pre-mRNA interferes with spliceosome binding to the exon 6 5′ splice site. Our results suggest that competing mRNA splicing and polyadenylation regulate exon 6 inclusion and consequently determine the ratios of soluble to membrane-bound IL7R. This may be relevant for both T cell ontogeny and function and development of multiple sclerosis.     

HRP-2, the Caenorhabditis elegans Homolog of Mammalian Heterogeneous Nuclear Ribonucleoproteins Q and R, Is an Alternative Splicing Factor That Binds to UCUAUC Splicing Regulatory Elements

Alternative splicing is regulated by cis sequences in the pre-mRNA that serve as binding sites for trans-acting alternative splicing factors. In a previous study, we used bioinformatics and molecular biology to identify and confirm that the intronic hexamer sequence UCUAUC is a nematode alternative splicing regulatory element. In this study, we used RNA affinity chromatography to identify trans factors that bind to this sequence. HRP-2, the Caenorhabditis elegans homolog of human heterogeneous nuclear ribonucleoproteins Q and R, binds to UCUAUC in the context of unc-52 intronic regulatory sequences as well as to RNAs containing tandem repeats of this sequence. The three Us in the hexamer are the most important determinants of this binding specificity. We demonstrate, using RNA interference, that HRP-2 regulates the alternative splicing of two genes, unc-52 and lin-10, both of which have cassette exons flanked by an intronic UCUAUC motif. We propose that HRP-2 is a protein responsible for regulating alternative splicing through binding interactions with the UCUAUC sequence.Alternative pre-mRNA splicing is a mechanism for generating multiple mRNA isoforms from a single gene. This process can allow a gene to encode for more than one protein isoform. For some genes, it is a mechanism for regulating message stability through production of alternative mRNA isoforms that are substrates for the nonsense-mediated mRNA decay pathway (1). The majority of human genes undergo alternative splicing (2), and the process can be regulated in tissue-specific and developmental stage-specific manners. Current models propose that cis elements on the pre-mRNA, in exons and introns, serve as recognition sites for trans-acting protein factors that bind to the pre-mRNA and regulate assembly of the splicing machinery, thus regulating splice site choice (3).In recent years, a number of groups have employed bioinformatics techniques to identify cis splicing regulatory elements (4). These techniques include using multiple interspecies sequence alignments to identify conserved intronic regions, identification of short sequences in exons that are bounded by weak consensus splice sites, and identification of common intronic sequences flanking similarly regulated alternative exons (5–9). These efforts have added many new sequences to the list of known and potential splicing regulators. The identification of the protein factor partners for these sequences will be important for understanding their function in alternative splicing regulation.Experimental approaches have identified alternative splicing factors that interact with specific cis elements (10), but the number of trans factors discovered still lags behind the number of newly identified cis element partners. Some examples of well-characterized cis element/trans-acting factor interactions include the NOVA K homology domain splicing factor binding to the sequence UCAY (11), the FOX splicing factors binding to the sequence UGCAUG (12–14), and hnRNP³ F/H proteins binding to the sequence GGGG (15, 16). By using cross-linking immunoprecipitation followed by large scale sequencing, entire catalogs of RNAs that the splicing factors NOVA, SF2/ASF, and FOX2 bind to in vivo have been determined (17–19). These approaches have led to models for how the proteins binding to their cis regulatory elements may alter splicing. These models include a role for the relative position of a cis element to an alternative cassette exon in determining alternative exon inclusion or skipping (18, 19).In a previous bioinformatics analysis of evolutionarily conserved intronic sequences flanking alternatively spliced exons, we identified the hexamer sequence UCUAUC as a novel splicing regulatory element (8). UCUAUC is found flanking both sides of alternative exon 16 of the unc-52 gene of Caenorhabditis elegans. Genetic analysis of a class of viable unc-52 mutants led to the discovery that exons 16–18 are alternative cassette exons and that every combination of skipping and inclusion of these three exons occurs (20). This splicing is regulated by the alternative splicing factor MEC-8 (21). Fig. 1A shows a schematic diagram of the alternatively spliced region of unc-52, with the MEC-8-enhanced alternative splicing events indicated. Using an unc-52 splicing reporter trans gene containing alternative exons 15–19, we previously reported that alternative splicing is regulated by the intronic motif UCUAUC in the intron downstream of exon 16 (8). In addition we showed that this element works cooperatively with a UGCAUG hexamer (the consensus FOX-1-binding site) in the upstream intron to regulate alternative splicing (8).Open in a separate windowFIGURE 1.RNA affinity chromatography identifies HRP-2 as binding to UCUAUC elements. A, schematic representation of the alternatively spliced region of unc-52 (adapted from Ref. 21). The alternative splicing events promoted by MEC-8 are indicated by bold lines. The lines next to introns 15 and 16 are the sites of the UCUAUC elements in those introns whose sequences were used in the RNA affinity chromatography. B, table showing sequences of RNAs immobilized to beads in the RNA affinity chromatography experiment. C, Coomassie-stained SDS-PAGE analysis of RNA affinity chromatography. C. elegans embryo extract was incubated with the different immobilized RNA substrates listed on top of the gel. Proteins identified by mass spectrometry are listed to the right of the gel, with arrows pointing to coincident protein bands. D, the left panel shows the silver stain result for the RNA affinity chromatography experiment. Each lane represents a different immobilized substrate, as indicated above. The band corresponding to HRP-2 is indicated by an arrow. The right panel is an immunoblot of the same gel using anti-HRP-2 polyclonal antibody. E, anti-HRP-2 immunoblot of an RNA affinity chromatography experiment for the indicated substrates.In this study, we report the results of a biochemical identification of a protein factor from C. elegans that binds to the UCUAUC intronic splicing regulatory element. We transcribed different short RNA sequences containing the UCUAUC element in its native intronic context, or as part of a repeating unit, and immobilized these onto agarose beads. After passing embryo extracts across these beads, we found that the protein HRP-2, the C. elegans homolog of the mammalian hnRNP Q/R proteins, binds to this sequence with high affinity. By using RNAi to reduce the level of HRP-2 in worms, we observed changes in alternative splicing of unc-52 and lin-10, two genes that contain UCUAUC elements in introns flanking alternative exons. We propose that HRP-2 is an alternative splicing factor that works through the UCUAUC intronic elements to regulate alternative splicing.     

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.     

An Intronic Sequence Element Mediates Both Activation and Repression of Rat Fibroblast Growth Factor Receptor 2 Pre-mRNA Splicing

Alternative splicing of fibroblast growth factor receptor 2 (FGF-R2) is an example of highly regulated alternative splicing in which exons IIIb and IIIc are utilized in a mutually exclusive manner in different cell types. The importance of this splicing choice is highlighted by studies which indicate that deregulation of the FGF-R2 splicing is associated with progression of prostate cancer. Loss of expression of a IIIb exon-containing isoform of FGF-R2 [FGF-R2 (IIIb)] accompanies the transition of a well-differentiated, androgen-dependent rat prostate cancer cell line, DT3, to the more aggressive, androgen-independent AT3 cell line. We have used transfection of rat FGF-R2 minigenes into DT3 and AT3 cancer cell lines to study the mechanisms that control alternative splicing of rat FGF-R2. Our results support a model in which an important cis-acting element located in the intron between these alternative exons mediates activation of splicing using the upstream IIIb exon and repression of the downstream IIIc exon in DT3 cells. This element consists of 57 nucleotides (nt) beginning 917 nt downstream of the IIIb exon. Analysis of mutants further demonstrates that an 18-nt “core sequence” within this element is most crucial for its function. Based on our observations, we have termed this sequence element ISAR (for intronic splicing activator and repressor), and we suggest that factors which bind this sequence are required for maintenance of expression of the FGF-R2 (IIIb) isoform.     

Depletion of TDP 43 overrides the need for exonic and intronic splicing enhancers in the human apoA-II gene

Exon 3 of the human apolipoprotein A-II (apoA-II) gene is efficiently included in the mRNA although its acceptor site is significantly weak because of a peculiar (GU)₁₆ tract instead of a canonical polypyrimidine tract within the intron 2/exon 3 junction. Our previous studies demonstrated that the SR proteins ASF/SF2 and SC35 bind specifically an exonic splicing enhancer (ESE) within exon 3 and promote exon 3 splicing. In the present study, we show that the ESE is necessary only in the proper context. In addition, we have characterized two novel sequences in the flanking introns that modulate apoA-II exon 3 splicing. There is a G-rich element in intron 2 that interacts with hnRNPH1 and inhibits exon 3 splicing. The second is a purine rich region in intron 3 that binds SRp40 and SRp55 and promotes exon 3 inclusion in mRNA. We have also found that the (GU) repeats in the apoA-II context bind the splicing factor TDP-43 and interfere with exon 3 definition. Significantly, blocking of TDP-43 expression by small interfering RNA overrides the need for all the other cis-acting elements making exon 3 inclusion constitutive even in the presence of disrupted exonic and intronic enhancers. Altogether, our results suggest that exonic and intronic enhancers have evolved to balance the negative effects of the two silencers located in intron 2 and hence rescue the constitutive exon 3 inclusion in apoA-II mRNA.     

A nonsense exon in the Tpm1 gene is silenced by hnRNP H and F

As well as generating protein isoform diversity, in some cases alternative splicing generates RNAs that harbor premature termination codons and that are subject to nonsense-mediated decay (NMD). We previously identified an apparent pseudo-exon in the rat α-tropomyosin (Tpm1) gene as a probable genuine alternatively spliced exon that causes NMD when spliced into Tpm1 RNA. Here, we report the analysis of cis-acting splicing regulatory elements within this “nonsense exon.” Guided by the data set of predicted splicing enhancer and silencer elements compiled by Zhang and Chasin, we made a series of mutations through the nonsense exon and found that like authentic exons it is densely packed with enhancer and silencer elements. Strikingly, 11 of 13 tested mutations behaved as predicted computationally. In particular, we found that a G-rich silencer at the 5′ end, which is crucial for skipping of the nonsense exon, functions by binding hnRNP-H and F.     

Dual Role of G-runs and hnRNP F in the Regulation of a Mutation-Activated Pseudoexon in the Fibrinogen Gamma-Chain Transcript

Most pathological pseudoexon inclusion events originate from single activating mutations, suggesting that many intronic sequences are on the verge of becoming exons. However, the precise mechanisms controlling pseudoexon definition are still largely unexplored. Here, we investigated the cis-acting elements and trans-acting regulatory factors contributing to the regulation of a previously described fibrinogen gamma-chain (FGG) pseudoexon, which is activated by a deep-intronic mutation (IVS6-320A>T). This pseudoexon contains several G-run elements, which may be bound by heterogeneous nuclear ribonucleoproteins (hnRNPs) F and H. To explore the effect of these proteins on FGG pseudoexon inclusion, both silencing and overexpression experiments were performed in eukaryotic cells. While hnRNP H did not significantly affect pseudoexon splicing, hnRNP F promoted pseudoexon inclusion, indicating that these two proteins have only partially redundant functions. To verify the binding of hnRNP F and the possible involvement of other trans-acting splicing modulators, pulldown experiments were performed on the region of the pseudoexon characterized by both a G-run and enrichment for exonic splicing enhancers. This 25-bp-long region strongly binds hnRNP F/H and weakly interacts with Serine/Arginine-rich protein 40, which however was demonstrated to be dispensable for FGG pseudoexon inclusion in overexpression experiments. Deletion analysis, besides confirming the splicing-promoting role of the G-run within this 25-bp region, demonstrated that two additional hnRNP F binding sites might instead function as silencer elements. Taken together, our results indicate a major role of hnRNP F in regulating FGG pseudoexon inclusion, and strengthen the notion that G-runs may function either as splicing enhancers or silencers of the same exon.     

Co-regulation of alternative splicing by diverse splicing factors in Caenorhabditis elegans

Regulation of alternative splicing is controlled by pre-mRNA sequences (cis-elements) and trans-acting protein factors that bind them. The combinatorial interactions of multiple protein factors with the cis-elements surrounding a given alternative splicing event lead to an integrated splicing decision. The mechanism of multifactorial splicing regulation is poorly understood. Using a splicing-sensitive DNA microarray, we assayed 352 Caenorhabditis elegans alternative cassette exons for changes in embryonic splicing patterns between wild-type and 12 different strains carrying mutations in a splicing factor. We identified many alternative splicing events that are regulated by multiple splicing factors. Many splicing factors have the ability to behave as splicing repressors for some alternative cassette exons and as splicing activators for others. Unexpectedly, we found that the ability of a given alternative splicing factor to behave as an enhancer or repressor of a specific splicing event can change during development. Our observations that splicing factors can change their effects on a substrate during development support a model in which combinatorial effects of multiple factors, both constitutive and developmentally regulated ones, contribute to the overall splicing decision.     

Biological Characterization of Rev Variation in Equine Infectious Anemia Virus

Sequence analysis identified significant variation in the second exon of equine infectious anemia virus (EIAV) rev. Functional analysis indicated that limited amino acid variation in Rev significantly altered the export activity of the protein but did not affect Rev-dependent alternative splicing. EIAV Rev can mediate export through two independent cis-acting Rev-responsive elements (RREs), and differences among Rev variants were more pronounced when both RREs were present. Variation in Rev may be an important mechanism for regulation of virus replication in vivo and may contribute to changes in clinical disease.     

Characterization of sequences and mechanisms through which ISE/ISS-3 regulates FGFR2 splicing

Alternative splicing of fibroblast growth factor receptor-2 (FGFR2) mutually exclusive exons IIIb and IIIc results in highly cell-type-specific expression of functionally distinct receptors, FGFR2-IIIb and FGFR2-IIIc. We previously identified an RNA cis-element, ISE/ISS-3, that enhanced exon IIIb splicing and silenced exon IIIc splicing. Here, we have performed comprehensive mutational analysis to define critical sequence motifs within this element that independently either enhance splicing of upstream exons or repress splicing of downstream exons. Such analysis included use of a novel fluorescence-based splicing reporter assay that allowed quantitative determination of relative functional activity of ISE/ISS-3 mutants using flow cytometric analysis of live cells. We determined that specific sequences within this element that mediate splicing enhancement also mediate splicing repression, depending on their position relative to a regulated exon. Thus, factors that bind the element are likely to be coordinately involved in mediating both aspects of splicing regulation. Exon IIIc silencing is dependent upon a suboptimal branchpoint sequence containing a guanine branchpoint nucleotide. Previous studies of exon IIIc splicing in HeLa nuclear extracts demonstrated that this guanine branchsite primarily impaired the second step of splicing suggesting that ISE/ISS-3 may block exon IIIc inclusion at this step. However, results presented here that include use of newly developed in vitro splicing assays of FGFR2 using extracts from a cell line expressing FGFR2-IIIb strongly suggest that cell-type-specific silencing of exon IIIc occurs at or prior to the first step of splicing.     

Identification of an exonic splicing silencer in exon 6A of the human VEGF gene

Background  

The different isoforms of vascular endothelial growth factor (VEGF) play diverse roles in vascular growth, structure and function. Alternative splicing of the VEGF gene results in the expression of three abundant isoforms: VEGF121, VEGF165 and VEGF189. The mRNA for VEGF189 contains the alternatively spliced exon 6A whereas the mRNA for VEGF165 lacks this exon. The objective of this study was to identify the cis elements that control utilization of exon 6A. A reporter minigene was constructed (pGFP-E6A) containing the coding sequence for GFP whose translation was dependent on faithful splicing for removal of the VEGF exon 6A. To identify cis-acting splicing elements, sequential deletions were made across exon 6A in the pGFP-E6A plasmid.