Genomic HEXploring allows landscaping of novel potential splicing regulatory elements

Effective splice site selection is critically controlled by flanking splicing regulatory elements (SREs) that can enhance or repress splice site use. Although several computational algorithms currently identify a multitude of potential SRE motifs, their predictive power with respect to mutation effects is limited. Following a RESCUE-type approach, we defined a hexamer-based ‘HEXplorer score’ as average Z-score of all six hexamers overlapping with a given nucleotide in an arbitrary genomic sequence. Plotted along genomic regions, HEXplorer score profiles varied slowly in the vicinity of splice sites. They reflected the respective splice enhancing and silencing properties of splice site neighborhoods beyond the identification of single dedicated SRE motifs. In particular, HEXplorer score differences between mutant and reference sequences faithfully represented exonic mutation effects on splice site usage. Using the HIV-1 pre-mRNA as a model system highly dependent on SREs, we found an excellent correlation in 29 mutations between splicing activity and HEXplorer score. We successfully predicted and confirmed five novel SREs and optimized mutations inactivating a known silencer. The HEXplorer score allowed landscaping of splicing regulatory regions, provided a quantitative measure of mutation effects on splice enhancing and silencing properties and permitted calculation of the mutationally most effective nucleotide.     

Identification of a prognostic alternative splicing signature in oral squamous cell carcinoma

Alternative splicing (AS) is critically associated with tumorigenesis and patient's prognosis. Here, we systematically analyzed survival-associated AS signatures in oral squamous cell carcinoma (OSCC) and evaluated their prognostic predictive values. Survival-related AS events were identified by univariate and multivariate Cox regression analyses using OSCC data from the TCGA head neck squamous cell carcinoma data set. The Percent Spliced In calculated by SpliceSeq from 0 to 1 was used to quantify seven types of AS events. A predictive model based on AS events was constructed by least absolute shrinkage and selection operator Cox regression assay and further validated using a training-testing cohort design. Patient survival was estimated using the Kaplan–Meier method and compared with Log-rank test. The receiver operating characteristics curve area under the curves was used to evaluate the predictive abilities of these predictive models. Furthermore, gene–gene interaction networks and the splicing factors (SFs)-AS regulatory network was generated by Cytoscape. A total of 825 survival-related AS events within 719 genes were identified in OSCC samples. The integrative predictive model was better at predicting outcomes of patients as compared to those models built with the individual AS event. The predictive model based on three AS-related genes also effectively predicted patients’ survival. Moreover, seven survival-related SFs were detected in OSCC including RBM4, HNRNPD, and HNRNPC, which have been linked to tumorigenesis. The SF-AS network revealed a significant correlation between survival-related AS genes and these SFs. Our findings revealed a systemic portrait of survival-associated AS events and the splicing network in OSCC, suggesting that AS events might serve as novel prognostic biomarkers and therapeutic targets for OSCC.     

Computational identification of tissue-specific alternative splicing elements in mouse genes from RNA-Seq

Tissue-specific alternative splicing is a key mechanism for generating tissue-specific proteomic diversity in eukaryotes. Splicing regulatory elements (SREs) in pre-mature messenger RNA play a very important role in regulating alternative splicing. In this article, we use mouse RNA-Seq data to determine a positive data set where SREs are over-represented and a reliable negative data set where the same SREs are most likely under-represented for a specific tissue and then employ a powerful discriminative approach to identify SREs. We identified 456 putative splicing enhancers or silencers, of which 221 were predicted to be tissue-specific. Most of our tissue-specific SREs are likely different from constitutive SREs, since only 18% of our exonic splicing enhancers (ESEs) are contained in constitutive RESCUE-ESEs. A relatively small portion (20%) of our SREs is included in tissue-specific SREs in human identified in two recent studies. In the analysis of position distribution of SREs, we found that a dozen of SREs were biased to a specific region. We also identified two very interesting SREs that can function as an enhancer in one tissue but a silencer in another tissue from the same intronic region. These findings provide insight into the mechanism of tissue-specific alternative splicing and give a set of valuable putative SREs for further experimental investigations.     

VERSE: a varying effect regression for splicing elements discovery

Identification of splicing regulatory elements (SREs) deserves special attention because these cis-acting short sequences are vital parts of splicing code. The fact that a variety of other biological signals cooperatively govern the splicing pattern indicates the necessity of developing novel tools to incorporate information from multiple sources to improve splicing factor binding sites prediction. Under this context, we proposed a Varying Effect Regression for Splicing Elements (VERSE) to discover intronic SREs in the proximity of exon junctions by integrating other biological features. As a result, 1562 intronic SREs were identified in 16 human tissues, many of which overlapped with experimentally verified binding motifs for several well-known splicing factors, including FOX-1, PTB, hnRNP A/B, hnRNP F/H, and so on. The discovered tissue, region, and conservation preferences of the putative motifs demonstrate that splice site selection is a complicated process that needs subtle and delicate regulation. VERSE may serve as a powerful tool to not only discover SREs by incorporating additional informative signals but also precisely quantify their varying contribution under different biological contexts.     

Systematic characterization of short intronic splicing-regulatory elements in SMN2 pre-mRNA

Intronic splicing enhancers and silencers (ISEs and ISSs) are two groups of splicing-regulatory elements (SREs) that play critical roles in determining splice-site selection, particularly for alternatively spliced introns or exons. SREs are often short motifs; their mutation or dysregulation of their cognate proteins frequently causes aberrant splicing and results in disease. To date, however, knowledge about SRE sequences and how they regulate splicing remains limited. Here, using an SMN2 minigene, we generated a complete pentamer-sequence library that comprises all possible combinations of 5 nucleotides in intron 7, at a fixed site downstream of the 5′ splice site. We systematically analyzed the effects of all 1023 mutant pentamers on exon 7 splicing, in comparison to the wild-type minigene, in HEK293 cells. Our data show that the majority of pentamers significantly affect exon 7 splicing: 584 of them are stimulatory and 230 are inhibitory. To identify actual SREs, we utilized a motif set enrichment analysis (MSEA), from which we identified groups of stimulatory and inhibitory SRE motifs. We experimentally validated several strong SREs in SMN1/2 and other minigene settings. Our results provide a valuable resource for understanding how short RNA sequences regulate splicing. Many novel SREs can be explored further to elucidate their mechanism of action.     

Polypyrimidine tract‐binding protein homologues from Arabidopsis underlie regulatory circuits based on alternative splicing and downstream control

Alternative splicing (AS) of precursor mRNAs is a widespread phenomenon in plants; however, many questions, especially regarding its regulation and functional implications, remain to be elucidated. In vertebrates, polypyrimidine tract‐binding proteins (PTBs) have been identified as key splicing factors influencing splice site selection and orchestrating coordinated splicing programmes during developmental processes. Here, we analysed three PTB homologues from Arabidopsis thaliana and provide evidence for their gene regulatory potential based on AS and a splicing‐independent mechanism. Our data reveal that Arabidopsis PTB homologues are subject to extensive auto‐ and cross‐regulation via AS‐coupled nonsense‐mediated decay, thereby establishing a basis for interlinking their expression. Furthermore, the multiple modes of action of Arabidopsis PTB homologues are reflected in their subcellular localization in the nucleus, cytosol and processing bodies. This work provides insight into the regulation of AS in plants and highlights the regulatory potential of the multifunctional plant PTB homologues, which might have important implications in diverse biological processes.     

Identification of a cis element for tissue-specific alternative splicing of chloroplast ascorbate peroxidase pre-mRNA in higher plants

Alternative splicing events in the 3'-terminal region of chloroplast ascorbate peroxidase (chlAPX) pre-mRNA in spinach and tobacco, which produced four types of mRNA variants, one form (tAPX-I) encoding thylakoid-bound APX (tAPX) and three forms (sAPX-I, -II, and -III) encoding stromal APX (sAPX), were regulated in a tissue-specific manner. The ratio of the level of sAPX mRNAs (sAPX-I, -II, and -III) to tAPX-I mRNA was close to 1 in leaf, whereas the ratio in root was greatly elevated due to an increase in sAPX-III and a decrease in tAPX-I resulting from the alternative excision of intron 11 and intron 12, respectively. A putative splicing regulatory cis element (SRE), which is highly conserved in the sequences of chlAPX genes of higher plants, was identified upstream of the acceptor site in intron 12. The deletion of the SRE sequence diminished the splicing efficiency of intron 12 in tobacco leaf in vivo. Gel-shift analysis showed that SRE interacts strongly with a nuclear protein from leaves but not those from the roots of spinach and tobacco. These results indicate that the tissue-specific alternative splicing of chlAPX pre-mRNA is regulated by the splicing enhancer SRE.     

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.     

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.     

The secondary structure of the human immunodeficiency virus type 1 transcript modulates viral splicing and infectivity

Splicing of human immunodeficiency virus type 1 (HIV-1) exon 6D is regulated by the presence of a complex splicing regulatory element (SRE) sequence that interacts with the splicing factors hnRNP H and SC35. In this work, we show that, in the context of the wild-type viral sequence, hnRNP H acts as a repressor of exon 6D inclusion independent of its binding to the SRE. However, hnRNP H binding to the SRE acts as an enhancer of exon 6D inclusion in the presence of a critical T-to-C mutation. These seemingly contrasting functional properties of hnRNP H appear to be caused by a change in the RNA secondary structure induced by the T-to-C mutation that affects the spatial location of bound hnRNP H with respect to the exon 6D splicing determinants. We propose a new regulatory mechanism mediated by RNA folding that may also explain the dual properties of hnRNP H in splicing regulation.     

Integrative genome-wide analysis reveals cooperative regulation of alternative splicing by hnRNP proteins

Understanding how RNA binding proteins control the splicing code is fundamental to human biology and disease. Here, we present a comprehensive study to elucidate how heterogeneous nuclear ribonucleoparticle (hnRNP) proteins, among the most abundant RNA binding proteins, coordinate to regulate alternative pre-mRNA splicing (AS) in human cells. Using splicing-sensitive microarrays, crosslinking and immunoprecipitation coupled with high-throughput sequencing (CLIP-seq), and cDNA sequencing, we find that more than half of all AS events are regulated by multiple hnRNP proteins and that some combinations of hnRNP proteins exhibit significant synergy, whereas others act antagonistically. Our analyses reveal position-dependent RNA splicing maps, in vivo consensus binding sites, a surprising level of cross- and autoregulation among hnRNP proteins, and the coordinated regulation by hnRNP proteins of dozens of other RNA binding proteins and genes associated with cancer. Our findings define an unprecedented degree of complexity and compensatory relationships among hnRNP proteins and their splicing targets that likely confer robustness to cells.