Rna-Seq Differential Transcription Analysis Of Tcga Colorectal Cancer (Crc) Transcriptomes Reveals Subtype-Specific Isoform Usage

Background: Differential transcription analysis may reveal transcriptome aberrations associated with cancer progression and help identify important drivers in molecularly distinct subtypes of a disease. Methods: We developed a data-driven bioinformatics pipeline (DiffSplice) that performs differential transcription analysis at the alternative splicing (AS) level and allows accurate transcriptome reconstruction through joint analysis of all RNA-seq samples. In this study, we applied DiffSplice to RNA-seq data from TCGA CRC (n = 686) and evaluated the expression patterns of the detected AS events across the four consensus molecular subtypes (CMS) of the disease (n = 511): CMS1 - MSI Immune (hypermutated, n = 76); CMS2 - Canonical (epithelial, WNT/MYC activation, n = 220); CMS3 - Metabolic (epithelial, enriched for KRAS mutations, n = 72); and CMS4 - Mesenchymal (stromal invasion, TGFB activation, n = 143). Results: A total of 15,613 AS events were detected, including 5,128 events that involved novel splice junctions not cataloged in UCSC GAF2.0 transcriptome annotation. Exon skipping was the most common pattern with 6,482 occurrences, followed by 4,962 alternative splice sites, 2,689 alternative 5′/3′ transcription sites and 1,262 retained introns. In the subset of samples with CMS classification, a cross-subtype comparative analysis revealed 1,835 differentially transcribed loci from 1,634 genes, including 759 novel junctions (FDR Discussion: Here we present an in-depth transcriptome analysis of TCGA CRC data across the molecular subtypes of the disease. We were able to identify differential transcription patterns in an unbiased way by using a scalable and efficient computational approach that models both annotated and novel splice variants simultaneously. Our pipeline revealed a number of promising candidates that may be acting as key modulators in specific disease subsets and deserve further investigation. Citation Format: Yin Hu, Rodrigo Dienstmann, Justin Guinney. RNA-seq differential transcription analysis of TCGA colorectal cancer (CRC) transcriptomes reveals subtype-specific isoform usage. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 1932. doi:10.1158/1538-7445.AM2015-1932