Motivation: Alternative splicing (AS) is a pre-mRNA maturation process leading to the expression of multiple mRNA variants from the same primary transcript. More than 90% of human genes are expressed via AS. Therefore, quantifying the inclusion level of every exon is crucial for generating accurate transcriptomic maps and studying the regulation of AS.

Results: Here we introduce SpliceTrap, a method to quantify exon inclusion levels using paired-end RNA-seq data. Unlike other tools, which focus on full-length transcript isoforms, SpliceTrap approaches the expression-level estimation of each exon as an independent Bayesian inference problem. In addition, SpliceTrap can identify major classes of alternative splicing events under a single cellular condition, without requiring a background set of reads to estimate relative splicing changes. We tested SpliceTrap both by simulation and real data analysis, and compared it to state-of-the-art tools for transcript quantification. SpliceTrap demonstrated improved accuracy, robustness and reliability in quantifying exon-inclusion ratios.

Conclusions: SpliceTrap is a useful tool to study alternative splicing regulation, especially for accurate quantification of local exon-inclusion ratios from RNA-seq data.

Availability and Implementation: SpliceTrap can be implemented online through the CSH Galaxy server http://cancan.cshl.edu/splicetrap and is also available for download and installation at http://rulai.cshl.edu/splicetrap/.

Contact: michael.zhang@utdallas.edu

Supplementary Information: Supplementary data are available at Bioinformatics online.

The splicing-factor oncoprotein SRSF1 (also known as SF2/ASF) is upregulated in breast cancers. We investigated SRSF1’s ability to transform human and mouse mammary epithelial cells in vivo and in vitro. SRSF1-overexpressing COMMA-1D cells formed tumors, following orthotopic transplantation to reconstitute the mammary gland. In 3-D culture, SRSF1-overexpressing MCF-10A cells formed larger acini than control cells, reflecting increased proliferation and delayed apoptosis during acinar morphogenesis. These effects required the first RNA-recognition motif and nuclear functions of SRSF1. SRSF1 overexpression promoted alternative splicing of BIM and BIN1 isoforms that lack pro-apoptotic functions and contribute to the phenotype. Finally, SRSF1 cooperated specifically with MYC to transform mammary epithelial cells, in part by potentiating eIF4E activation, and these cooperating oncogenes are significantly co-expressed in human breast tumors. Thus, SRSF1 can promote breast cancer, and SRSF1 itself or its downstream effectors may be valuable targets for therapeutics development.

The SR protein splicing factor SRSF1 is a potent proto-oncogene that is frequently upregulated in cancer. Here we show that SRSF1 is a direct target of the transcription-factor oncoprotein MYC. These two oncogenes are significantly co-expressed in lung carcinomas, and MYC knockdown downregulates SRSF1 expression in lung-cancer cell lines. MYC directly activates transcription of SRSF1 through two non-canonical E-boxes in its promoter. The resulting increase in SRSF1 protein is sufficient to modulate alternative splicing of a subset of transcripts. In particular, MYC induction leads to SRSF1-mediated alternative splicing of the signaling kinase MKNK2 and the transcription factor TEAD1. SRSF1 knockdown reduces MYC’s oncogenic activity, decreasing proliferation and anchorage-independent growth. These results suggest a mechanism for SRSF1 upregulation in tumors with elevated MYC, and identify SRSF1 as a critical MYC target that contributes to its oncogenic potential by enabling MYC to regulate the expression of specific protein isoforms through alternative splicing.

Alternative splicing of the pyruvate kinase M gene (PK-M) can generate the M2 isoform and promote aerobic glycolysis and tumor growth. However, the cancer-specific alternative splicing regulation of PK-M is not completely understood. Here, we demonstrate that PK-M is regulated by reciprocal effects on the mutually exclusive exons 9 and 10, such that exon 9 is repressed and exon 10 is activated in cancer cells. Strikingly, exonic, rather than intronic, cis-elements are key determinants of PK-M splicing isoform ratios. Using a systematic sub-exonic duplication approach, we identify a potent exonic splicing enhancer in exon 10, which differs from its homologous counterpart in exon 9 by only two nucleotides. We identify SRSF3 as one of the cognate factors, and show that this serine/arginine-rich protein activates exon 10 and mediates changes in glucose metabolism. These findings provide mechanistic insights into the complex regulation of alternative splicing of a key regulator of the Warburg effect, and also have implications for other genes with a similar pattern of alternative splicing.