Background: Liver metastases are the major cause of colorectal cancer (CRC)-related deaths. However, there is no reliable clinical predictor for CRC progression to liver metastasis. In this study, we investigated possible predictors (miRNAs and biomarkers) for clinical application.
Methodology: The Gene Expression Omnibus (GEO) datasets GSE49355, GSE41258 and GSE81558 for genes and GSE54088 and GSE56350 for miRNAs were used to identify common differentially expressed genes (DEGs) and miRNAs between primary CRC tissues and liver metastases. The identified miRNAs and their targets from the DEGs were verified in datasets comprising gene, miRNA and miRNA exosome profiles of CRC patients with no distant metastases (M0) and distant metastases (M1); the interaction networks and pathways were also mapped.
Results: There were 49 upregulated and 13 downregulated DEGs and 16 downregulated and 14 upregulated miRNAs; between the DEGs and miRNA targets, there were five upregulated and four downregulated genes. MiR-20a was strongly correlated with the status of liver metastasis. MiR-20a, miR499a, and miR-576-5p were highly correlated with the metastatic outcomes. MiR-20a was significantly highly expressed in the M1 group. In an analysis of the miRNA target genes, we found that CDH2, KNG1, and MMP2 were correlated with CRC metastasis. We demonstrated a new possible pathway for CRC metastasis: miR-576-5p/F9, miR20a/MMP2, CTSK, MMP3, and miR449a/P2RY14. The regulation of IGF transport and uptake by IGFBPs, extracellular matrix organization, signal transduction and the immune system were the enriched pathways.
Conclusion: This model can predict CRC to liver metastases and the pathways involved, which can be clinically applicable.