Identification and verification of key cancer genes associated with prognosis of colorectal cancer based on bioinformatics analysis.

Objectives: The biomarkers targeting colorectal cancer (CRC) prognosis are short of high accuracy and sensitivity in clinic. Through bioinformatics analysis, we aim to identify and confirm a series of key genes referred to the diagnosis and prognosis of CRC.

Methods: GSE31905, GSE35279, and GSE41657 were selected as complete RNA sequencing data sets of CRC and colorectal mucosa (CRM) tissues from the NCBI-GEO database, and the differentially expressed genes (DEGs) were analyzed. The common DEGs in these 3 data sets were obtained by Venn map, and enriched by STRING network system and Cytoscape software. The Kaplan-Meier plotter website was used to verify the correlation between the enriched genes and the prognosis of CRC.

Results: For the whole RNA sequencing data sets of CRC and normal intestinal mucosa samples, the DEGs of CRC and CRM in the 3 data sets (|logFC|>2 and <0.05) were screened by GEO2R tool in NCBI-GEO database. By using Venn graph analysis software, the intersection of up-regulated/down-regulated genes in 3 GSE datasets was obtained, and a total 105 up-regulated genes and 140 down-regulated genes were found in the 3 samples. The up-regulated/down-regulated genes were introduced into the STRING network system to obtain the interacting genes. The interacting gene sets were introduced into Cytoscape software, and 61 up-regulated genes were found by Molecular Complex Detection (MCODE) plug-in. Through the Kaplan-Meier plotter website, we found that EPHB2, KLK8, DIAPH3, STC2, OXTR, MMP7, MET, KRT85, KRT6B, KRT23, and KLK10 genes were highly expressed in CRC, and were related to the prognosis.

Conclusions: The above 11 genes verified by bioinformatics retrieval and analysis can predict the poor prognosis of CRC to a certain extent, and they provide a possible target for the diagnosis and treatment of CRC.