ABSTRACT
A number of patients with colon cancer with local or local advanced disease suffer from recurrence and there is an urgent need for better prognostic biomarkers in this setting. Here, the transcriptomic landscape of mRNAs, long noncoding RNAs, snRNAs, small nucleolar RNAs (snoRNAs), small Cajal body-specific RNAs, pseudogenes, and circular RNAs, as well as RNAs denoted as miscellaneous RNAs, was profiled by total RNA sequencing. In addition to well-known coding and noncoding RNAs, differential expression analysis also uncovered transcripts that have not been implicated previously in colon cancer, such as RNA5SP149, RNU4-2, and SNORD3A. Moreover, there was a profound global up-regulation of snRNA pseudogenes, snoRNAs, and rRNA pseudogenes in more advanced tumors. A global down-regulation of circular RNAs in tumors relative to normal tissues was observed, although only a few were expressed differentially between tumor stages. Many previously undescribed transcripts, including RNU6-620P, RNU2-20P, VTRNA1-3, and RNA5SP60, indicated strong prognostic biomarker potential in receiver operating characteristics analyses. In summary, this study unveiled numerous differentially expressed RNAs across various classes between recurrent and nonrecurrent colon cancer. Notably, there was a significant global up-regulation of snRNA pseudogenes, snoRNAs, and rRNA pseudogenes in advanced tumors. Many of these newly discovered candidates demonstrate a strong prognostic potential for stage II colon cancer.
Subject(s)
Colonic Neoplasms , Gene Expression Regulation, Neoplastic , Neoplasm Recurrence, Local , Humans , Colonic Neoplasms/genetics , Colonic Neoplasms/pathology , Colonic Neoplasms/metabolism , Neoplasm Recurrence, Local/genetics , Neoplasm Recurrence, Local/pathology , Biomarkers, Tumor/genetics , Biomarkers, Tumor/metabolism , RNA, Untranslated/genetics , RNA, Small Nucleolar/genetics , RNA, Small Nucleolar/metabolism , Transcriptome/genetics , Male , Gene Expression Profiling/methods , FemaleABSTRACT
Circular RNAs (circRNA) are covalently closed molecules that can play important roles in cancer development and progression. Hundreds of differentially expressed circRNAs between tumors and adjacent normal tissues have been identified in studies using RNA sequencing or microarrays, emphasizing a strong translational potential. Most previous studies have been performed using RNA from bulk tissues and lack information on the spatial expression patterns of circRNAs. Here, we showed that the majority of differentially expressed circRNAs from bulk tissue analyses of colon tumors relative to adjacent normal tissues were surprisingly not differentially expressed when comparing cancer cells directly with normal epithelial cells. Manipulating the proliferation rates of cells grown in culture revealed that these discrepancies were explained by circRNAs accumulating to high levels in quiescent muscle cells due to their high stability; on the contrary, circRNAs were diluted to low levels in the fast-proliferating cancer cells due to their slow biogenesis rates. Thus, different subcompartments of colon tumors and adjacent normal tissues exhibited striking differences in circRNA expression patterns. Likewise, the high circRNA content in muscle cells was also a strong confounding factor in bulk analyses of circRNAs in bladder and prostate cancers. Together, these findings emphasize the limitations of using bulk tissues for studying differential circRNA expression in cancer and highlight a particular need for spatial analysis in this field of research. SIGNIFICANCE: The abundance of circRNAs varies systematically between subcompartments of solid tumors and adjacent tissues, implying that differentially expressed circRNAs discovered in bulk tissue analyses may reflect differences in cell type composition between samples.