Comprehensive analysis of microarray expression profiles of circRNAs and lncRNAs with associated co-expression networks in human colorectal cancer

Abstract

Increasing data demonstrate that circular RNAs (circRNAs) and long non-coding RNAs (lncRNAs) play important roles in tumorigenesis. However, the mechanisms in colorectal cancer (CRC) remain unclear. Here, hundreds of significantly expressed circRNAs, and thousands of lncRNAs as well as mRNAs were identified. By qRT-PCR, one abnormal circRNA, lncRNA, and three mRNAs were verified in 24 pairs of tissues and blood samples, respectively. Then, by GO analysis, we found that the gene expression profile of linear counterparts of upregulated circRNAs in human CRC tissues preferred positive regulation of GTPase activity, cellular protein metabolic process, and protein binding, while that of downregulated circRNAs of CRC preferred positive regulation of cellular metabolic process, acetyl-CoA metabolic process, and protein kinase C activity. Moreover, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that p53 signaling pathway was an important pathway in upregulated protein-coding genes, whereas cyclic guanosine monophosphate-protein kinase G (cGMP–PKG) signaling pathway was the top enriched KEGG pathway for downregulated transcripts. Furthermore, lncRNA–mRNA co-expression analysis demonstrated that downregulated lncRNA uc001tma.3 was negatively with CDC45 and positively with ELOVL4, BVES, FLNA, and HSPB8, while upregulated lncRNA NR_110882 was positively with FZD2. In addition, lncRNA–transcription factor (TF) co-expression analysis showed that the most relevant TFs were forkhead box protein A1 (FOXA1), transcription initiation factor TFIID submint 7 (TAF7), and adenovirus early region 1A(E1A)-associated protein p300 (EP300). Our findings offer a fresh view on circRNAs and lncRNAs and provide the foundation for further study on the potential roles of circRNAs and lncRNAs in colorectal cancer.