Published online Jun 26, 2021. doi: 10.12998/wjcc.v9.i18.4520
Peer-review started: December 11, 2020
First decision: January 17, 2021
Revised: January 26, 2021
Accepted: February 26, 2021
Article in press: February 26, 2021
Published online: June 26, 2021
Processing time: 182 Days and 0.5 Hours
The high morbidity and mortality of colorectal cancer (CRC) have posed great threats to human health. Circular RNA (CircRNA) and microRNA (miRNA), acting as competing endogenous RNAs (ceRNAs), have been found to play vital roles in carcinogenesis. However, the biological function of ceRNAs in CRC pathogenesis and prognosis remains largely unexplored.
To identify the CRC-specific circRNA-miRNA-mRNA regulatory network and uncover the subnetwork associated with its prognosis.
CircRNAs, miRNAs and mRNAs differentially expressed (DE) in CRC tissues were selected by expression file analysis in the Gene Expression Omnibus (GEO) database, and the downstream target molecules of circRNAs and miRNAs were predicted. Then, the intersection of differentially expressed RNA molecules with the predicted targets was determined to obtain a ceRNA network. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were conducted to elucidate the possible mechanism of pathogenesis. A survival analysis using the gene profiles and clinical information in The Cancer Genome Atlas (TCGA) database was performed to identify the mRNAs associated with the clinical outcome of CRC patients and construct a prognostic subnetwork.
We downloaded three datasets (GSE126095, GSE41655 and GSE41657) of large-scale CRC samples from the GEO database. There were 55 DEcircRNAs, 114 DEmiRNAs and 267 DEmRNAs in CRC tissues compared with normal tissues. After intersecting these molecules with predicted targets, 19 circRNAs, 13 miRNAs and 28 mRNAs were chosen to develop a circRNA-miRNA-mRNA network. GO and KEGG functional enrichment analyses indicated that the retinol metabolic process, leukocyte chemotaxis, extracellular matrix remodeling, endoplasmic reticulum stress, alcohol dehydrogenase activity, gastric acid secretion, nitrogen metabolism and NOD-like receptor signaling pathway might participate in the tumorigenesis of CRC. After verifying the identified mRNA effect in the TCGA database, we finally recognized 3 mRNAs (CA2, ITLN1 and LRRC19) that were significantly associated with the overall survival of CRC patients and constructed a ceRNA subnetwork including 5 circRNAs (hsa_circ_0080210, hsa_circ_0007158, hsa_circ_0000375, hsa_circ_0018909 and hsa_circ_0011536) and 3 miRNAs (hsa-miR-601, hsa-miR-671-5p and hsa-miR-765), which could contain innovative and noninvasive indicators for the early screening and prognostic prediction of CRC.
We proposed a circRNA-miRNA-mRNA regulatory network closely associated with the progression and clinical outcome of CRC that might include promising biomarkers for carcinogenesis and therapeutic targets.
Core Tip: The biological functions of circRNA and miRNA interactions and their potential as noninvasive biomarkers have not been well elucidated in colorectal cancer (CRC). In this study, we constructed a circRNA-miRNA-mRNA regulatory network with 19 circRNAs, 13 miRNAs and 28 mRNAs. GO and KEGG analyses indicated several signaling pathways probably involved in tumorigenesis. After being combined with survival analysis, a prognostic subnetwork was constructed including 5 circRNAs, 3 miRNAs and 3 mRNAs, which may represent novel diagnostic and prognostic candidate biomarkers, as well as therapeutic targets of CRC.