Identification of mRNA isoform switching in breast cancer

Alternative splicing provides a major mechanism to generate protein diversity. Increasing evidence suggests a link of dysregulation of splicing associated with cancer. While previous genomic-based studies demonstrated the expression of a handful of tumor-specific isoforms, genome-wide alterations in the balance between isoforms and cancer subtypes is understudied.

We systematically analyzed the isoform-level expression patterns and isoform switching events of 819 breast tumor and normal samples assayed by mRNA-seq from TCGA project. On average, 2.2 isoforms per gene were detected and 67.5 % of detected genes (i.e. expressed) showed 1-2 isoforms only. While the majority of isoforms for a given gene were positively correlated with each other and the overall gene level, 470 pairs of isoforms displayed an inverse correlation suggesting a switching event. Most of the isoform switching events were associated with molecular subtypes, including a Basal-like-associated switching in CTNND1. 88 genes showed switching independent of subtypes, among which the isoform pattern of PRICKLE1 was associated with a large genomic signature of biological significance.

Our results reveal that the majority of genes do not undergo complex mRNA splicing within breast cancers, and that there is a general concordance in isoform and gene expression levels in breast tumors. We identified hundreds of isoform switching events across breast tumors, most of which were associated with differences in tumor subtypes. As exemplified by the detailed analysis of CTNND1 and PRICKLE1, these isoform switching events potentially provide new insights into the post-transcriptional regulatory mechanisms of tumor subtypes and cancer biology.