Widespread intron retention diversifies most cancer transcriptomes

1The discovery of high-frequency mutations affecting components of the RNA splicing machinery is one of the most unexpected results of cancer genome sequencing. ‘Spliceosomal mutations’ are enriched in diverse diseases, including myelodysplastic syndromes, lymphoid leukemias, and solid tumors of the lung, breast, pancreas, and eye, and most commonly cause specific missense changes to the SF3B1, SRSF2, and U2AF1 proteins. Mechanistic studies revealed that U2AF1 mutations alter the preferred 3′ splice site sequence both in vivo and in vitro, thereby influencing genome-wide recognition of alternative and constitutive 3′ splice sites. SRSF2 mutations similarly alter interactions between SRSF2 and pre-mRNA, resulting in altered exon recognition that promotes dysplastic hematopoiesis.

2In addition to the direct genetic link between abnormal RNA splicing and tumorigenesis provided by point mutations affecting the spliceosome, indirect evidence suggests that important differences distinguish RNA splicing in normal versus cancerous cells even in the absence of these mutations. Small molecules that inhibit splicing have antitumor activity; the SF3b component PHF5A is differentially required for constitutive splicing in glioblastoma versus normal neural stem cells; RNA splicing is reportedly noisier in cancers than normal cells; increased intron retention is associated with SETD2 mutations in kidney cancer and castration resistance in prostate cancer. These and other studies together suggest that common RNA processing differences may distinguish cancer and normal cells irrespective of tissue of origin. However, this hypothesis has not been systematically tested.

3Here, we took advantage of the comprehensive transcriptome data produced by The Cancer Genome Atlas (TCGA) to identify large-scale differences in RNA splicing between tumor and normal control samples across 16 distinct cancer types. While we observed no obvious biases in cassette exon recognition or 5′ or 3′ splice site recognition, almost all analyzed cancer types exhibited increased levels of intron retention relative to normal controls. The sole exception was breast cancer, for which intron retention characterized normal breast rather than cancer samples. Our results indicate that intron retention is a common correlate of tumorigenesis, and suggest that an abundance of intron-containing mRNAs in cancer cells may increase the diversity of many cancer transcriptomes.

4Unprocessed RNA-seq reads from TCGA were downloaded from CGHub, using all solid tumors with patient-matched samples from the adjacent normal tissue, as well as unmatched acute myeloid leukemia (AML) and breast cancer samples (the unmatched breast cancer samples were only used for the subgroup analysis involving all 1,080 cancer patients). Samples were identified using cgquery v2.1, with ‘state = live’, ‘library_strategy = RNA-Seq’, and ‘sample_type = 0*’ or ‘sample_type = 1*’ for cancer and normal samples, respectively, and the sequence data were downloaded with the GeneTorrent client software. For samples extracted from CGHub prior to November 2013, the raw reads were extracted in BAM format and converted to FASTQ format using sam2fastq v1.2 from UNC Bioinformatics Utilities. For samples extracted after November 2013, the reads were available directly in FASTQ format. All samples were sequenced using the Illumina Genome Analyzer or HiSeq, and reads were from unstranded paired-end libraries, with the exception of a subset of the uterine corpus endometrial carcinoma samples, which were single-end. Samples where the sequencing protocol was ‘TotalRNASeqV2’ were excluded, in order to include only poly(A)-selected RNA-seq libraries.

5RNA-seq reads from four healthy bone marrow donors were obtained from the NCBI Gene Expression Omnibus (GEO) under accession number GSE61410. The library characteristics of these samples match those of the AML RNA-seq samples (average read count: 75 M; paired-end libraries; with read length 2×49 or 50 nt). Subcellular fractionation RNA-seq reads from MCF-7 and K562 cells were obtained from GEO under accession number GSE30567, and restricted to poly(A)-selected libraries. RNA-seq data from breast cancer cell lines were obtained from GSE52643 and GSE48213.

6Alternative splicing events were classified as cassette exons, competing 5′ and 3′ splice sites, and retained introns, using annotations from MISO v2.0. Constitutive splice junctions and introns were defined as junctions that were not alternatively spliced in any isoform of the UCSC knownGene track. RNA annotation files were created for transcripts and splice junctions individually, to map the RNA-seq reads to each annotation set separately. The RNA transcript annotation is a combination of isoforms from MISO v2.0, UCSC knownGene, and the Ensembl 71 gene annotation. The RNA splice junction annotation was created using an enumerating of all possible combinations of annotated splice sites as previously described.

7All RNA-seq reads were processed using a standardized pipeline. Step 1: Map all reads to the UCSC hg19 (NCBI GRCh37) human genome assembly using Bowtie v1.0.0 and RSEM v.1.2.4 modified to call Bowtie with the -v 2 mapping strategy and invoked with the arguments --bowtie-m 100 --bowtie-chunkmbs 500 --calc-ci --output-genome-bam on the gene annotation file. Step 2: Filter the resulting BAM file to: (1) remove alignments with mapq scores of 0; and (2) require a minimum splice junction overhang of 6 bp. Step 3: Align all previously unaligned reads to the splice junction annotation file with TopHat v2.0.8b invoked with the arguments --bowtie1 --read-mismatches 3 --read-edit-dist 2 --no-mixed --no-discordant --min-anchor-length 6 --splice-mismatches 0 --min-intron-length 10 --max-intron-length 1,000,000 --min-isoform-fraction 0.0 --no-novel-juncs --no-novel-indels --raw-juncs. The parameters --mate-inner-dist and --mate-std-dev were determined by mapping to constitutive coding exons as determined with MISO’s exon_utils.py script. Step 4: Filter the resulting alignments as in step 2. Step 5: Merge the BAM files from TopHat and RSEM to generate a final set of all read alignments.

8RNA transcript levels from the 50 genes included in the PAM50 classifier were normalized using the trimmed mean of M values (TMM) method with the scaling factors calculated based on protein-coding transcripts only. The scaled centroids from pam50.robust in the ‘genefu’ R package were used to cluster all tumor samples into the five intrinsic molecular subtypes using intrinsic.cluster.predict with parameter do.prediction.strength = TRUE. The resulting subtype probabilities were assessed to confirm that less than 5 % of samples had a maximum probability below 0.5.

9MISO v2.0 was used to quantify isoform ratios for all annotated alternative splicing events (cassette exons, competing 5′ and 3′ splice sites, and retained introns). For constitutive junctions and retention of constitutive introns the isoform ratios were estimated using junction-spanning reads, as previously described. Individual tumor samples were analyzed by comparing them directly to the patient-matched normal sample (solid tumors) or to the median across normal bone marrow samples (AML). For each sample pair we restricted the analysis to splicing events that had at least 20 reads supporting either or both isoforms, and that were alternatively spliced in our data. Events were defined as differentially spliced within a sample pair if they satisfied the following criteria: (1) at least 20 relevant reads in both samples; (2) a change in isoform ratio of at least 10 %; and (3) a Bayes factor greater than or equal to 1 for differences in isoform ratios, calculated using Wagenmakers’s framework. Sets of tumor or normal samples were analyzed using a two-sided Wilcoxon rank-sum test with the total number of isoform reads from the test versus the reference set of samples (this method was only used for the mutation/subgroup analyses of AML and breast cancer). Events were differentially spliced if they had: (1) a total of at least 20 reads in both sample sets; (2) a change in isoform ratio of at least 10 %; and (3) a P value below 0.01.

10For each cancer type, we identified the parent genes of introns that were differentially retained in at least 20 % of samples. These genes were compared against all protein-coding genes using the R package goseq to test for enrichment of gene ontology (GO) Biological Process terms. The ‘Wallenius’ method was used, the results were corrected for gene length bias, and the resulting false discovery rates were corrected using the Benjamini-Hochberg approach. Only terms with at least two ancestors were tested, and only the most detailed terms were returned by the analysis, to eliminate parent terms associated with generic biological processes.

11All 5′ and 3′ splice site scores were calculated using the maximum entropy modeling method. Intron lengths and GC content were calculated using the hg19 human genome assembly, and the GC content was averaged across each individual intron.

12To calculate the proportion of the variation in intron retention that could be explained by mRNA levels, we used a generalized additive model (GAM) from the R-package mgcv v1.8.3. For each cancer type, we selected the most variable introns, defined as having a standard deviation of changes in intron retention for tumor-normal pairs across all patient samples exceeding 0.1 (AML: 852 introns, breast cancer: 621 introns, colon cancer: 740 introns). We extracted all genes associated with the gene ontology terms ‘mRNA catabolic process’ (GO:0006402), ‘mRNA splicing, via spliceosome’ (GO:0000398) and ‘mRNA transport’ (GO:0051028), and used the log2 fold-change of TMM normalized gene expression of the matched tumor-normal pairs. GAM models for all individual intron/mRNA combinations were calculated using a Gaussian distribution with ‘identity’ as the link function, and the method ‘GCV.Cp’ was used for estimating the smoothing parameters of the log2 fold-changes for the mRNAs.

13Somatic mutations from each AML tumor sample were extracted from cBioPortal using the CGDS-R package. For all genes mutated in more than five samples, intron retention was analyzed in the mutated sample set compared to all wild-type tumor samples.

14To identify systematic differences in splicing between cancer and normal cells, we analyzed the transcriptomes of 16 solid and liquid tumor types that were sequenced as part of TCGA (Table 1). For the solid tumors, we restricted to samples for which patient-matched adjacent normal tissue was available in order to control for potential genetic differences in splicing, as well as reduce artifacts arising from different handling of unmatched samples. For the acute myeloid leukemia (AML) samples, we used four bone marrow samples from healthy donors as normal controls, as patient-matched controls were not available due to the circulating nature of leukemic cells. To help ensure that we studied the mature mRNA pool rather than pre-mRNA or RNA degradation products, we restricted to samples for which both the cancer and adjacent normal samples were sequenced as poly(A)-selected libraries. A total of 805 matched samples were available for solid cancers originating from the bladder, breast, colon, head and neck, kidney, liver, lung, prostate, rectum, stomach, thyroid, and endometrium, as well as for acute myeloid leukemia (Table 1). We quantified global patterns of splicing across these cancer and normal samples, which were sequenced to an average coverage of 74 million paired-end reads (+/- 22 million), using a database of approximately 125,000 annotated alternative splicing events and approximately 160,000 constitutive splice junctions (see Methods).

15We first identified annotated alternative splicing events that were differentially spliced in cancer versus normal cells. We defined differentially spliced events as those exhibiting a difference in isoform ratio of ≥10 % between a solid tumor sample and its matched normal, or alternately between an acute myeloid leukemia sample and the median across four normal bone marrow samples. The relative frequency of differential splicing between cancer and normal tissues differed substantially between cancer types, with AML and liver cancers exhibiting the most and fewest differentially spliced events on average. Cassette exons comprised approximately 50 % to 60 % of all differentially spliced events, while competing 5′ and 3′ splice sites and retained introns comprised less than 20 % each across all cancer types (Fig. 1a-d, Table 1), consistent with patterns observed in studies of tissue-specific alternative splicing. Differential splicing of cassette exons and competing 5′ and 3′ splice sites was balanced, in the sense that cassette exons were not preferentially included/excluded and intron-proximal/distal 5′ or 3′ splice sites were not preferentially used.

16In contrast to cassette exon recognition or competing splice site usage, intron retention was markedly imbalanced. All cancer types exhibited a strong enrichment (median of 1- to 20-fold) for increased retention of alternative introns in cancers, with the notable exception of breast cancer, for which adjacent normal tissue exhibited a four-fold enrichment for increased intron retention relative to cancer samples (Fig. 1d). We therefore tested whether constitutive as well as alternative introns were retained. Again requiring a difference in isoform ratio of ≥10 %, we found that constitutively spliced introns were frequently unspliced in cancers relative to normal controls (with the exception of breast cancer), although the enrichment was more modest than for alternative introns (Fig. 1e).

17Increased intron retention could potentially be a side effect of the reportedly high levels of noisy splicing in cancers rather than a specific failure to remove intronic sequences. To test this hypothesis, we estimated levels of noisy splicing by quantifying mis-splicing of constitutive junctions, such as splicing from the 5′ splice site of a constitutive junction to a different 3′ splice site of the gene. In contrast to intron retention, where all cancer types except for breast cancer exhibited increased levels of unspliced mRNA, we did not observe uniform increases or decreases in mis-splicing of constitutive junctions in cancer relative to normal transcriptomes (Fig. 1f). Specific introns of a gene were typically retained while others were removed with high efficiency (Fig. 1g).

18While intron retention was a readily detectable feature of all cancer types except for breast, the degree of intron retention in any particular sample was highly variable, ranging from increases of 2- to 40-fold in AML and colon cancer (Fig. 2a). The average number of differentially retained introns within a cancer type correlated with the absolute magnitude of the change in intron retention levels (Fig. 2b). The enrichment for retention of alternative versus constitutive introns was well-correlated within each sample for all cancer types (Fig. 2c), suggesting that the propensity to retain introns is a patient-specific characteristic of the transcriptional environment. Here and subsequently, we focused on AML and colon cancer as exemplars of cancers exhibiting strong intron retention, and on breast cancer as the sole example of a cancer exhibiting decreased intron retention relative to normal controls.

19As we restricted our analysis to poly(A)-selected cDNA libraries, increased RNA-seq coverage of introns was likely due to the presence of incompletely spliced mRNA rather than inefficiently degraded spliced introns. To confirm that the signal did indeed arise from incompletely spliced mRNA, we quantified intron retention using only reads crossing the exon-intron boundary, which cannot arise from spliced introns or lariats. We observed similar enrichment for intron retention with this conservative measure, indicating that the intronic signal arises from incompletely spliced mRNA (Figs. 1d and 2d).

20Of the 16 cancer types that we studied, only breast cancer was associated with decreased, rather than increased, intron retention relative to normal controls. Infrequent intron retention in breast cancer versus normal breast tissue could be due to efficient intron removal in breast cancer, or instead inefficient intron removal in normal breast. To distinguish between these possibilities, we measured intron retention in breast cancer versus each of the 14 other solid cancer types, as well as in normal breast versus each of the 14 other normal control tissues. We did not observe consistently increased or decreased intron retention in breast cancer versus other cancers. In contrast, there was a strong bias towards increased intron retention in normal breast versus other control tissues, wherein normal breast exhibited more frequent intron retention than any other control tissue except for endometrium (Fig. 2e). Median gene expression of MKI67, encoding the proliferation marker Ki-67, was >10-fold higher in breast cancer samples than in adjacent controls, indicating that enhanced intron retention in normal breast is unlikely to be explained by contamination of adjacent control tissue by cancerous cells.

21Breast cancer’s status as an outlier in our analysis could potentially be due to the different cell types represented in the cancer versus normal samples. Breast cancer predominantly arises from the mammary ductal or (occasionally) lobular cells, whereas the surrounding stroma consists of fibroblasts, adipocytes, and cells from the immune system. We therefore compared intron retention levels in 50 tumorigenic and six non-tumorigenic breast cell line models (Fig. 2f). While the efficiency of intron removal varied across the cell lines, we did not observe a consistent association between tumorigenic/non-tumorigenic status and degree of intron retention. This cell line analysis is consistent with the hypothesis that cell type differences between tumor and normal samples may contribute to breast cancer’s outlier status, although the known propensity of cell lines to adopt ‘cancer-like’ RNA processing characteristics such as preferential usage of short 3′ UTRs suggests that the same could occur for intron retention.

22Select introns, such as two adjacent introns in FUS, were recurrently retained in cancers from many different tissues of origin (Fig. 3a). Quantifying this genome-wide, we found that while most introns were retained in just a few samples of a particular cancer, 1,205 and 171 introns, respectively, exhibited increased or decreased retention in >10 % of cancer samples relative to normal controls across all cancer types (Fig. 3b). An unsupervised cluster analysis confirmed that while some retained introns are shared across most cancer types, the majority are either specific to the cancer of origin, or present at low frequencies in multiple cancers. For example, cancers arising from similar tissues of origin (for example, kidney) clustered together, indicating that they exhibit similar genome-wide patterns of intron retention (Fig. 3c). The preponderance of differentially retained introns in any given cancer sample is typically specific to that cancer, although a subset of introns are frequently retained across diverse cancers.

23We next tested whether genes containing the small subset of introns that are frequently retained across diverse cancers are preferentially involved in particular biological processes (Fig. 3d). We identified Gene Ontology (GO) terms that were enriched for genes exhibiting changes in intron retention affecting at least 20 % of samples within each cancer type. Strikingly, the most enriched terms were all involved in RNA processing, with mRNA export as the most enriched term (Fig. 3d). This association between frequently retained introns and mRNA export is primarily mediated by nine genes that exhibit differential intron retention in more than one-third of all samples on average (Fig. 3e). These genes encode SR proteins, DEAD box proteins associated with the nuclear export factor NXF1, and components of the TREX mRNA export complex.

24As intron-containing mRNAs are frequently retained in the nucleus and/or subject to degradation, preferential intron retention in cancer cells could potentially suppress gene expression in the absence of compensatory transcriptional upregulation. Consistent with this hypothesis, a recent study reported that intron retention was associated with lower expression of the parent genes, although other studies have reported the opposite trend. We therefore tested whether differential intron retention was associated with altered gene expression by simply computing the Pearson correlation between differences in intron retention and parent gene expression for each cancer sample and matched normal control. For all cancers with the exception of AML, we observed a weak negative correlation between differential intron retention and expression of the corresponding parent genes (Fig. 3f). The origin of AML’s outlier status is unclear, although it may be due to the lack of patient-matched controls. While some tumor-normal pairs exhibited consistent anti-correlation between intron retention and gene expression for many genes, in many cases the anti-correlation was driven by a few outliers with unusually prominent changes in both intron retention and gene expression (Fig. 3g). Our data are consistent with the hypothesis that intron retention contributes to alterations in gene expression, although the changes that we observed were generally modest.

25Previous genome-wide studies have reported that cis-acting sequence features distinguish introns susceptible to retention, including weak splice sites, shorter lengths, and higher GC content. Consistent with these studies, introns that were preferentially retained in either normal or cancer samples had weak 5′ and 3′ splice sites, were very short, and had high GC content relative to constitutively spliced introns (Fig. 4a-d). These trends persisted even after controlling for intron length. Trends for introns that were preferentially retained or removed in cancer cells were qualitatively similar, suggesting that similar sequence features predispose introns to inefficient splicing in both normal and cancer cells.

26In addition to cis-acting sequence features, we investigated the extent to which trans-acting factors may affect differential intron retention. We used a generalized additive model to identify relationships between changes in expression of genes encoding trans-acting factors relevant to the RNA life cycle and intron retention. While the majority of RNA processing factors were not systematically predictive of changes in intron retention, expression of a subset of genes was strongly associated with intron retention (Fig. 4e). Many of these genes, which were involved in RNA splicing, transport, and degradation, were associated with intron retention across AML, colon cancer, and breast cancer (Fig. 4f). These genes encoded proteins involved in 3′ splice site selection (for example, SF3A1, SF3B1, SF3B2), the nuclear pore complex (NUP60, NUPL1, NPIPA1, RANBP2) and RNA degradation and nonsense-mediated decay (UPF1, SMG1, PAN2, XRN1).

27We next sought to identify clinical variables that potentially contributed to the high levels of inter-sample variability in intron retention that we observed (Fig. 2a, Table 1). We first tested whether specific somatic mutations were associated with enhanced or diminished intron retention relative to normal controls. We focused on AML due to the mounting evidence that presence/absence of specific driver mutations is clinically and therapeutically relevant. For each of 22 genes that were somatically mutated in at least five of the 169 AML samples, we computed the number of differentially retained introns in wild-type versus mutant samples. Mutations in most genes were not associated with increased or decreased intron retention relative to other AML samples (although all AML samples exhibited strong intron retention relative to normal bone marrow; Fig. 2a). Mutations in RUNX1, IDH1, and IDH2, in contrast, were associated with a strongly increased intron retention relative to wild-type AML samples (Fig. 4g). For IDH1 and IDH2, the augmented intron retention that we observed is likely an underestimate of the true effect. We treated IDH1 and IDH2 separately so that they could serve as a rough validation of each other, and so many of the IDH1 wild-type samples contained IDH2 mutations and vice versa.

28We next tested whether specific molecular subtypes of tumors exhibited stronger or weaker signals of intron retention. We focused on breast cancer, which is frequently divided into subtypes that correlate strongly with prognosis and response to therapy. Using all 1,080 poly(A)-selected breast cancer samples available in TCGA, we computed the number of differentially retained introns in each molecular subtype relative to all available normal controls (Fig. 4h). (As we performed a subgroup rather than per-sample analysis, we did not restrict to samples with matched tumor-normal pairs.) Normal-like samples exhibited the fewest differentially retained introns relative to all normal controls, consistent with their known similarity to normal breast tissue, and also exhibited relatively balanced increases and decreases in intron retention, as did luminal A and B samples. In contrast, basal-like and Her2 amplified or over-expressing samples exhibited the most differential intron retention relative to normal tissue, as well as the strongest bias towards decreased intron retention. As basal-like and Her2 amplified or over-expressing samples are associated with the worst prognosis of the subtypes analyzed, it is interesting to note that a recent study observed high levels of unspliced mRNA in bone metastases from castration-resistant prostate cancer relative to primary prostate cancer.

29Alternative splicing of specific genes has long been known to contribute to cancer initiation, progression, and metastasis. The recent discovery of spliceosomal mutations in diverse cancers, coupled with the identification of spliceosomal components that are differentially required in cancer versus normal cells and the ongoing development of RNA splicing inhibitors as potential antitumor drugs, suggests that the RNA splicing process may frequently differ substantially between cancer and normal cells. Our study bolsters this conjecture by demonstrating that cancer transcriptomes from diverse tissues of origin frequently exhibit marked increases in intron retention relative to adjacent normal tissue even in the absence of mutational insults to the splicing machinery.

30Intriguingly, our simple correlation of clinical variables with intron retention suggested that mutations affecting factors not canonically involved in RNA processing, such as IDH1 and IDH2, may contribute to splicing dysregulation in cancer. While our analysis cannot distinguish between direct and indirect effects, previous reports suggest that aberrant methylation induced by IDH1 and IDH2 mutations could indeed alter intron recognition. Exons exhibit preferential DNA methylation relative to introns and changes in DNA methylation can influence alternative splicing through differential CTCF binding. We were unable to test whether differentially retained introns in IDH1 or IDH2-mutant samples exhibited differential DNA methylation due to the low probe coverage of introns on the Illumina HumanMethylation450K BeadChip used by TCGA. However, taken together with a recent report that mutations affecting the epigenetic factor SETD2 are associated with loss of histone H3 lysine 36 tri-methylation (H3K36me3) and increased intron retention in clear cell renal cell carcinoma, our results suggest that crosstalk between epigenetics and RNA splicing may link a diverse spectrum of mutations to RNA processing defects in cancer.

31While our study and others make increasingly clear that intron retention is much more common in human cells than previously believed, the functional consequences of this intron retention remain unknown. Impaired nuclear export and cytoplasmic RNA degradation both provide barriers to stable expression of intron-containing mRNAs, and likely reduce many such mRNAs to dead-end products of incomplete splicing. As TCGA data are from whole-cell transcriptomes, we were unable to test whether intron-containing mRNAs are present in the cytoplasmic mRNA pools of primary cancers. We therefore instead used subcellular fractionation data from MCF-7 (breast cancer) and K562 (erythroleukemic) cell lines as rough approximations of the cell types represented in the TCGA primary breast cancer and AML samples. We restricted to introns that were differentially retained in breast cancer or AML relative to normal controls, and quantified the relative abundance of mRNAs containing these introns in the nuclear and cytoplasmic fractions of MCF-7 or K562 cells (Fig. 5). Most intron-containing mRNAs were present at higher levels in the nuclear fraction, as expected, but many retained introns were nonetheless present at high levels in both fractions. Approximately 75 % or 52 % of introns that were differentially retained in breast cancer or AML relative to controls were retained at rates ≥10 % in the cytoplasmic fraction (meaning that the intron-containing mRNA constituted ≥10 % of the mRNAs from the parent gene). With the caveat that data from MCF-7 and K562 cell lines may not closely correspond to in vivo processes, we conclude that biased intron retention may generate a diversity of intron-containing mRNAs that are exported to the cytoplasm and sufficiently stable to be readily detectable.

32If translated, intron-containing mRNAs could produce novel peptides with unknown consequences, ranging from a potentially adaptive expansion of the proteomic repertoire to a deleterious triggering of tumor immunogenicity. While genes involved in many biological processes are affected by intron retention, it is interesting to note that the small subset of introns that are preferentially retained across many cancers are frequently found in genes encoding RNA processing factors involved in RNA export and splicing (Fig. 3d). Finally, it is tempting to speculate that the apparent abundance of intron-containing mRNAs in tumor cells may contribute to the observed antitumor activity of compounds that inhibit RNA splicing catalysis. Future work will reveal whether intron-containing mRNAs are important contributors to tumor cell biology, host-tumor interaction, and sensitivity to antitumor agents, or instead simply by-products of the widespread molecular dysregulation that accompanies oncogenesis.

33Our data show that intron retention characterizes the transcriptomes of many primary cancers, indicating that globally abnormal RNA splicing is a common correlate of tumorigenesis even in the absence of direct mutational insults to the RNA splicing machinery. Different introns were preferentially retained in specific samples and cancer types, although a small minority of introns were frequently subject to retention in diverse cancers. The degree of intron retention appears to be influenced by a complex interplay of cis-acting sequence features and RNA processing factors acting in trans, rather than by a single mechanism. Many introns that were preferentially retained in primary cancers were detectable in the cytoplasmic fractions of cancer cell lines, suggesting that at least some intron-containing mRNAs are sufficiently stable to contribute to the transcriptional diversity of cancer cells. Further investigation will determine whether these intron-containing mRNAs are translated into potentially functional proteins, or instead dead-end products of inefficient splicing.