Somatic Copy Number Alterations in Colorectal Cancer Lead to a Differentially Expressed ceRNA Network (ceRNet)

Abstract

Colorectal cancer (CRC) represents the second deadliest malignancy worldwide. Around 75% of CRC patients exhibit high levels of chromosome instability that result in the accumulation of somatic copy number alterations. These alterations are associated with the amplification of oncogenes and deletion of tumor-ppressor genes and contribute to the tumoral phenotype in different malignancies. Even though this relationship is well known, much remains to be investigated regarding the effect of said alterations in long non-coding RNAs (lncRNAs) and, in turn, the impact these alterations have on the tumor phenotype. The present study aimed to evaluate the role of differentially expressed lncRNAs coded in regions with copy number alterations in colorectal cancer patient samples. We downloaded RNA-seq files of the Colorectal Adenocarcinoma Project from the The Cancer Genome Atlas (TCGA) repository (285 sequenced tumor tissues and 41 non-tumor tissues), evaluated differential expression, and mapped them over genome sequencing data with regions presenting copy number alterations. We obtained 78 differentially expressed (LFC > 1|< −1, padj < 0.05) lncRNAs, 410 miRNAs, and 5028 mRNAs and constructed a competing endogenous RNA (ceRNA) network, predicting significant lncRNA–miRNA–mRNA interactions. Said network consisted of 30 lncRNAs, 19 miRNAs, and 77 mRNAs. To understand the role that our ceRNA network played, we performed KEGG and GO analysis and found several oncogenic and anti-oncogenic processes enriched by the molecular players in our network. Finally, to evaluate the clinical relevance of the lncRNA expression, we performed survival analysis and found that C5orf64, HOTAIR, and RRN3P3 correlated with overall patient survival. Our results showed that lncRNAs coded in regions affected by SCNAs form a complex gene regulatory network in CCR.