Identification of Hub Genes for Colorectal Cancer with Liver Metastasis Using miRNA-mRNA Network

Liver metastasis is an important cause of death in patients with colorectal cancer (CRC). Increasing evidence indicates that microRNAs (miRNAs) are involved in the pathogenesis of colorectal cancer liver metastasis (CRLM). This study is aimed at exploring the potential miRNA-mRNA regulatory network.

From the GEO database, we downloaded the microarray datasets GSE56350 and GSE73178. GEO2R was used to conduct differentially expressed miRNAs (DEMs) between CRC and CRLM using the GEO2R tool. Then, GO and KEGG pathway analysis for differentially expressed genes (DEGs) performed via DAVID. A protein-protein interaction (PPI) network was constructed by the STRING and identified by Cytoscape. Hub genes were identified by miRNA-mRNA network. Finally, the expression of the hub gene expression was assessed in the GSE81558.

The four DEMs (hsa-miR-204-5p, hsa-miR-122-5p, hsa-miR-95-3p, and hsa-miR-552-3p) were identified as common DEMs in GSE56350 and GSE73178 datasets. The SP1 was likely to adjust the upregulated DEMs; however, the YY1 could regulate both the upregulated and downregulated DEMs. A total of 3925 genes (3447 upregulated DEM genes and 478 downregulated DEM genes) were screened. These predicted genes were mainly linked to Platinum drug resistance, Cellular senescence, and ErbB signaling pathway. Through the gene network construction, most of the hub genes were found to be modulated by hsa-miR-204-5p, hsa-miR-122-5p, hsa-miR-95-3p, and hsa-miR-552-3p. Among the top 20 hub genes, the expression of CREB1, RHOA, and EGFR was significantly different in the GSE81558 dataset.

In this study, miRNA-mRNA networks in CRLM were screened between CRC patients and CRLM patients to provide a new method to predict for the pathogenesis and development of CRC.