Circular RNA-Associated Competing Endogenous RNA Network and Prognostic Nomogram for Patients With Colorectal Cancer

Background: Genetic characteristics remain underutilized for establishing prognostic models for colorectal cancer (CRC). We explored the underlying regulatory mechanisms of circular RNAs (circRNAs) that act as competing endogenous RNAs (ceRNAs) and constructed a gene-based nomogram to predict overall survival (OS) in patients with CRC. Methods: We obtained circRNA expression profiling data from the Gene Expression Omnibus (GEO) database. MicroRNA (miRNA) and mRNA expression profiles, with associated clinical data, were obtained from The Cancer Genome Atlas (TCGA). A ceRNA network was established using Cytoscape. Interactions between differential genes were analyzed, and hub genes were identified using the cytoHubba application. The R package clusterProfiler was used to evaluate the Gene Ontology (GO) annotations of the differentially expressed mRNAs and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Database-extracted patients were randomized into a training and validation cohorts. A prognostic model was developed using the training set based on multivariate Cox analyses and was then assessed in the validation set. The accuracy of the model was evaluated using discrimination and calibration plots. Results: Thirteen circRNAs, 62 miRNAs, and 301 mRNAs were used to construct the ceRNA network; 10 hub genes were identified via the PPI network. Next, a circRNA- miRNA hub of gene-regulatory modules was established based on four differentially expressed circRNAs, eight differentially expressed miRNAs, and nine differentially expressed mRNAs (DEmRNAs). GO and KEGG pathway analyses indicated the possible association of DEmRNAs with CRC onset and progression. Multivariate analyses revealed that age, tumor stage, and CXCR5 expression were independent risk factors for OS. A CXCR5-based model was developed to predict the OS of patients with CRC in our training set. Our nomogram showed relatively good accuracy, with C-indices of 0.757 and 0.702 in the training and validation sets, respectively. The areas under the curve of the nomograms predicting 3- and 5-years OS were 0.749 and 0.805 in the training set and 0.706 and 0.779 in the validation set, respectively. Conclusions: Our data suggested that the hsa_circ_00001666/has-mir-1229/CXCR5 axis plays an important role in the pathogenesis of CRC, thereby identifying a potential therapeutic target. The proposed CXCR5-based nomogram may also assist surgeons in devising personalized treatments for patients with this disease.