Abstract:Alternative splicing (AS), in higher eukaryotes, is one of the mechanisms of post-transcriptional regulation that generate multiple transcripts from the same gene. One particular mode of AS is the skipping event where an exon may be alternatively excluded or constitutively included in the resulting mature mRNA. Both transcript isoforms from this skipping event site, i.e. in which the exon is either included (inclusion isoform) or excluded (skipping isoform), are typically present in one cell, and maintain a subtle balance that is vital to cellular function and dynamics. However, how the prevailing conditions dictate which isoform is expressed and what biological factors might influence the regulation of this process remain areas requiring further exploration. In this study, we have developed a novel computational method, graph-based exon-skipping scanner (GESS), for de novo detection of skipping event sites from raw RNA-seq reads without prior knowledge of gene annotations, as well as for determining the dominant isoform generated from such sites. We have applied our method to publicly available RNA-seq data in GM12878 and K562 cells from the ENCODE consortium and experimentally validated several skipping site predictions by RT-PCR. Furthermore, we integrated other sequencing-based genomic data to investigate the impact of splicing activities, transcription factors (TFs) and epigenetic histone modifications on splicing outcomes. Our computational analysis found that splice sites within the skipping-isoform-dominated group (SIDG) tended to exhibit weaker MaxEntScan-calculated splice site strength around middle, ‘skipping’, exons compared to those in the inclusion-isoform-dominated group (IIDG). We further showed the positional preference pattern of splicing factors, characterized by enrichment in the intronic splice sites immediately bordering middle exons. Finally, our analysis suggested that different epigenetic factors may introduce a variable obstacle in the process of exon–intron boundary establishment leading to skipping events.

Transcriptome sequencing suggests that pre-mRNA splicing counteracts widespread intronic cleavage and polyadenylation

Interplay between splicing and transcriptional pausing exerts genome-wide control over alternative polyadenylation

Alternative polyadenylation by sequential activation of distal and proximal PolyA sites

Computational Analysis Reveals a Correlation of Exon-Skipping Events with Splicing, Transcription and Epigenetic Factors.

Genetic regulation of nascent RNA maturation revealed by direct RNA nanopore sequencing

Post-transcriptional Splicing of Nascent RNA Contributes to Widespread Intron Retention in Plants

by PAS-Seq Complex and dynamic landscape of RNA polyadenylation revealed Material Supplemental

Genome-wide Identification and Predictive Modeling of Tissue-Specific Alternative Polyadenylation.

Epigenetic Regulation of Alternative Splicing.

Pre-mRNA splicing inhibits m6A deposition, allowing longer mRNA half-life and flexible protein coding

<scp>RNA</scp> structures in alternative splicing and back‐splicing

Intronic Alus Influence Alternative Splicing

Distinct Polyadenylation Landscapes of Diverse Human Tissues Revealed by A Modified Pa-Seq Strategy

Dichotomous intronic polyadenylation profiles reveal multifaceted gene functions in the pan-cancer transcriptome

Complex and dynamic landscape of RNA polyadenylation revealed by PAS-Seq.

Study of the RNA splicing kinetics via in vivo 5-EU labeling

Rewriting endogenous human transcripts with -splicing

Interactions between RNA m6A modification, alternative splicing, and poly(A) tail revealed by MePAIso-seq2

Transcriptome Dynamics Through Alternative Polyadenylation in Developmental and Environmental Responses in Plants Revealed by Deep Sequencing

The RNA-binding protein hnRNPLL induces a T cell alternative splicing program delineated by differential intron retention in polyadenylated RNA

Alternative Polyadenylation A Mechanism Maximizing Transcriptome Diversity In Higher Eukaryotes