Long-range RNA pairings contribute to mutually exclusive splicing

Mutually exclusive splicing is an important means of increasing the protein repertoire, by which the Down''s syndrome cell adhesion molecule (Dscam) gene potentially generates 38,016 different isoforms in Drosophila melanogaster. However, the regulatory mechanisms remain obscure due to the complexity of the Dscam exon cluster. Here, we reveal a molecular model for the regulation of the mutually exclusive splicing of the serpent pre-mRNA based on competition between upstream and downstream RNA pairings. Such dual RNA pairings confer fine tuning of the inclusion of alternative exons. Moreover, we demonstrate that the splicing outcome of alternative exons is mediated in relative pairing strength-correlated mode. Combined comparative genomics analysis and experimental evidence revealed similar bidirectional structural architectures in exon clusters 4 and 9 of the Dscam gene. Our findings provide a novel mechanistic framework for the regulation of mutually exclusive splicing and may offer potentially applicable insights into long-range RNA–RNA interactions in gene regulatory networks.     

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 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.     

Autoregulation of Fox protein expression to produce dominant negative splicing factors

The Fox proteins are a family of regulators that control the alternative splicing of many exons in neurons, muscle, and other tissues. Each of the three mammalian paralogs, Fox-1 (A2BP1), Fox-2 (RBM9), and Fox-3 (HRNBP3), produces proteins with a single RNA-binding domain (RRM) flanked by N- and C-terminal domains that are highly diversified through the use of alternative promoters and alternative splicing patterns. These genes also express protein isoforms lacking the second half of the RRM (FoxΔRRM), due to the skipping of a highly conserved 93-nt exon. Fox binding elements overlap the splice sites of these exons in Fox-1 and Fox-2, and the Fox proteins themselves inhibit exon inclusion. Unlike other cases of splicing autoregulation by RNA-binding proteins, skipping the RRM exon creates an in-frame deletion in the mRNA to produce a stable protein. These FoxΔRRM isoforms expressed from cDNA exhibit highly reduced binding to RNA in vivo. However, we show that they can act as repressors of Fox-dependent splicing, presumably by competing with full-length Fox isoforms for interaction with other splicing factors. Interestingly, the Drosophila Fox homolog contains a nearly identical exon in its RRM domain that also has flanking Fox-binding sites. Thus, rather than autoregulation of splicing controlling the abundance of the regulator, the Fox proteins use a highly conserved mechanism of splicing autoregulation to control production of a dominant negative isoform.     

Role of six single nucleotide polymorphisms,risk factors in coronary disease,in OLR1 alternative splicing

The OLR1 gene encodes the oxidized low-density lipoprotein receptor (LOX-1), which is responsible for the cellular uptake of oxidized LDL (Ox-LDL), foam cell formation in atheroma plaques and atherosclerotic plaque rupture. Alternative splicing (AS) of OLR1 exon 5 generates two protein isoforms with antagonistic functions in Ox-LDL uptake. Previous work identified six single nucleotide polymorphisms (SNPs) in linkage disequilibrium that influence the inclusion levels of OLR1 exon 5 and correlate with the risk of cardiovascular disease. Here we use minigenes to recapitulate the effects of two allelic series (Low- and High-Risk) on OLR1 AS and identify one SNP in intron 4 (rs3736234) as the main contributor to the differences in exon 5 inclusion, while the other SNPs in the allelic series attenuate the drastic effects of this key SNP. Bioinformatic, proteomic, mutational and functional high-throughput analyses allowed us to define regulatory sequence motifs and identify SR protein family members (SRSF1, SRSF2) and HMGA1 as factors involved in the regulation of OLR1 AS. Our results suggest that antagonism between SRSF1 and SRSF2/HMGA1, and differential recognition of their regulatory motifs depending on the identity of the rs3736234 polymorphism, influence OLR1 exon 5 inclusion and the efficiency of Ox-LDL uptake, with potential implications for atherosclerosis and coronary disease.     

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.     

Multisite and bidirectional exonic splicing enhancer in CD44 alternative exon v3

The human CD44 gene encodes multiple isoforms of a transmembrane protein that differ in their extracellular domains as a result of alternative splicing of its variable exons. Expression of CD44 is tightly regulated according to the type and physiological status of a cell, with expression of high molecular weight isoforms by inclusion of variable exons and low molecular weight isoforms containing few or no variable exons. Human CD44 variable exon 3 (v3) can follow a specific alternative splicing route different from that affecting other variable exons. Here we map and functionally describe the splicing enhancer element within CD44 exon v3 which regulates its inclusion in the final mRNA. The v3 splicing enhancer is a multisite bipartite element consisting of a tandem nonamer, the XX motif, and an heptamer, the Y motif, located centrally in the exon. Each of the three sites of this multisite enhancer partially retains its splicing enhancing capacity independently from each other in CD44 and shows full enhancing function in gene contexts different from CD44. We further demonstrate that these motifs act cooperatively as at least two motifs are needed to maintain exon inclusion. Their action is differential with respect to the splice-site target abutting v3. The first X motif acts on the 3' splice site, the second X motif acts on both splice sites (as a bidirectional exonic splicing enhancer), and the Y motif acts on the 5' splice site. We also show that the multisite v3 splicing enhancer is functional irrespective of flanking intron length and spatial organization within v3.     

Alternative splicing of fibronectin mRNAs in chondrosarcoma cells: role of far upstream intron sequences

The fibronectin (FN) gene encodes multiple mRNAs through the process of alternative splicing, and production of certain isoforms is characteristic of a given cell type. Chondrocytes produce FNs that completely lack alternative exon EIIIA, and loss of inclusion of the exon is tightly linked to chondrogenic condensation of mesenchymal cells. The inclusion of a second exon, EIIIB, is high in embryonic cartilage, but declines with age. Multiple exons are omitted to produce the (V + C)-form that is highly specific for cartilage and chondrocytes. A rat chondrosarcoma cell line, RCS, was identified that preserves key features of the cartilage-specific splicing phenotype. RCS cells, which exclude exon EIIIA, and HeLa cells, which include exon EIIIA similar to mesenchymal cells, were used to assess the contribution of intron sequences flanking exon EIIIA to splicing regulation. Deletion of most of the intron downstream of the exon had little effect on splicing in either cell type. However, deletions within upstream intron 32-A reduced inclusion of the alternative exon in both cell types. The sequences involved lie more than 200 nucleotides away from the exon, but could not be localized to a single region by deletion mapping. These intronic sequences contribute to the efficiency of exon EIIIA recognition, but not to cell-type specific regulation. The normally inhibitory factor polypyrimidine tract binding protein promotes exon EIIIA inclusion in a manner that is partially dependent on the regulatory sequences within intron 32-A.     

Pressure-Overload Cardiac Hypertrophy Is Associated with Distinct Alternative Splicing Due to Altered Expression of Splicing Factors

Chronic pressure-overload cardiac hypertrophy is associated with an increased risk of morbidity/mortality, largely due to maladaptive remodeling and dilatation that progresses to dilated cardiomyopathy. Alternative splicing is an important biological mechanism that generates proteomic complexity and diversity. The recent development of next-generation RNA sequencing has improved our understanding of the qualitative signatures associated with alternative splicing in various biological conditions. However, the role of alternative splicing in cardiac hypertrophy is yet unknown. The present study employed RNA-Seq and a bioinformatic approach to detect the RNA splicing regulatory elements involved in alternative splicing during pressure-overload cardiac hypertrophy. We found GC-rich exonic motifs that regulate intron retention in 5′ UTRs and AT-rich exonic motifs that are involved in exclusion of the AT-rich elements that cause mRNA instability in 3′ UTRs. We also identified motifs in the intronic regions involved in exon exclusion and inclusion, which predicted splicing factors that bind to these motifs. We found, through Western blotting, that the expression levels of three splicing factors, ESRP1, PTB and SF2/ASF, were significantly altered during cardiac hypertrophy. Collectively, the present results suggest that chronic pressure-overload hypertrophy is closely associated with distinct alternative splicing due to altered expression of splicing factors.     

RBM4 promotes neuronal differentiation and neurite outgrowth by modulating Numb isoform expression

RBM4 participates in cell differentiation by regulating tissue-specific alternative pre-mRNA splicing. RBM4 also has been implicated in neurogenesis in the mouse embryonic brain. Using mouse embryonal carcinoma P19 cells as a neural differentiation model, we observed a temporal correlation between RBM4 expression and a change in splicing isoforms of Numb, a cell-fate determination gene. Knockdown of RBM4 affected the inclusion/exclusion of exons 3 and 9 of Numb in P19 cells. RBM4-deficient embryonic mouse brain also exhibited aberrant splicing of Numb pre-mRNA. Using a splicing reporter minigene assay, we demonstrated that RBM4 promoted exon 3 inclusion and exon 9 exclusion. Moreover, we found that RBM4 depletion reduced the expression of the proneural gene Mash1, and such reduction was reversed by an RBM4-induced Numb isoform containing exon 3 but lacking exon 9. Accordingly, induction of ectopic RBM4 expression in neuronal progenitor cells increased Mash1 expression and promoted cell differentiation. Finally, we found that RBM4 was also essential for neurite outgrowth from cortical neurons in vitro. Neurite outgrowth defects of RBM4-depleted neurons were rescued by RBM4-induced exon 9–lacking Numb isoforms. Therefore our findings indicate that RBM4 modulates exon selection of Numb to generate isoforms that promote neuronal cell differentiation and neurite outgrowth.     

Distinct Types of Disorder in the Human Proteome: Functional Implications for Alternative Splicing

Intrinsically disordered regions have been associated with various cellular processes and are implicated in several human diseases, but their exact roles remain unclear. We previously defined two classes of conserved disordered regions in budding yeast, referred to as “flexible” and “constrained” conserved disorder. In flexible disorder, the property of disorder has been positionally conserved during evolution, whereas in constrained disorder, both the amino acid sequence and the property of disorder have been conserved. Here, we show that flexible and constrained disorder are widespread in the human proteome, and are particularly common in proteins with regulatory functions. Both classes of disordered sequences are highly enriched in regions of proteins that undergo tissue-specific (TS) alternative splicing (AS), but not in regions of proteins that undergo general (i.e., not tissue-regulated) AS. Flexible disorder is more highly enriched in TS alternative exons, whereas constrained disorder is more highly enriched in exons that flank TS alternative exons. These latter regions are also significantly more enriched in potential phosphosites and other short linear motifs associated with cell signaling. We further show that cancer driver mutations are significantly enriched in regions of proteins associated with TS and general AS. Collectively, our results point to distinct roles for TS alternative exons and flanking exons in the dynamic regulation of protein interaction networks in response to signaling activity, and they further suggest that alternatively spliced regions of proteins are often functionally altered by mutations responsible for cancer.     

MBNL and CELF proteins regulate alternative splicing of the skeletal muscle chloride channel CLCN1

The expression and function of the skeletal muscle chloride channel CLCN1/ClC-1 is regulated by alternative splicing. Inclusion of the CLCN1 exon 7A is aberrantly elevated in myotonic dystrophy (DM), a genetic disorder caused by the expansion of a CTG or CCTG repeat. Increased exon 7A inclusion leads to a reduction in CLCN1 function, which can be causative of myotonia. Two RNA-binding protein families—muscleblind-like (MBNL) and CUG-BP and ETR-3-like factor (CELF) proteins—are thought to mediate the splicing misregulation in DM. Here, we have identified multiple factors that regulate the alternative splicing of a mouse Clcn1 minigene. The inclusion of exon 7A was repressed by MBNL proteins while promoted by an expanded CUG repeat or CELF4, but not by CUG-BP. Mutation analyses suggested that exon 7A and its flanking region mediate the effect of MBNL1, whereas another distinct region in intron 6 mediates that of CELF4. An exonic splicing enhancer essential for the inclusion of exon 7A was identified at the 5′ end of this exon, which might be inhibited by MBNL1. Collectively, these results provide a mechanistic model for the regulation of Clcn1 splicing, and reveal novel regulatory properties of MBNL and CELF proteins.     

Intronic alternative splicing regulators identified by comparative genomics in nematodes

Many alternative splicing events are regulated by pentameric and hexameric intronic sequences that serve as binding sites for splicing regulatory factors. We hypothesized that intronic elements that regulate alternative splicing are under selective pressure for evolutionary conservation. Using a Wobble Aware Bulk Aligner genomic alignment of Caenorhabditis elegans and Caenorhabditis briggsae, we identified 147 alternatively spliced cassette exons that exhibit short regions of high nucleotide conservation in the introns flanking the alternative exon. In vivo experiments on the alternatively spliced let-2 gene confirm that these conserved regions can be important for alternative splicing regulation. Conserved intronic element sequences were collected into a dataset and the occurrence of each pentamer and hexamer motif was counted. We compared the frequency of pentamers and hexamers in the conserved intronic elements to a dataset of all C. elegans intron sequences in order to identify short intronic motifs that are more likely to be associated with alternative splicing. High-scoring motifs were examined for upstream or downstream preferences in introns surrounding alternative exons. Many of the high- scoring nematode pentamer and hexamer motifs correspond to known mammalian splicing regulatory sequences, such as (T)GCATG, indicating that the mechanism of alternative splicing regulation is well conserved in metazoans. A comparison of the analysis of the conserved intronic elements, and analysis of the entire introns flanking these same exons, reveals that focusing on intronic conservation can increase the sensitivity of detecting putative splicing regulatory motifs. This approach also identified novel sequences whose role in splicing is under investigation and has allowed us to take a step forward in defining a catalog of splicing regulatory elements for an organism. In vivo experiments confirm that one novel high-scoring sequence from our analysis, (T)CTATC, is important for alternative splicing regulation of the unc-52 gene.     

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.