Identification of Transcriptional Markers and microRNA–mRNA Regulatory Networks in Colon Cancer by Integrative Analysis of mRNA and microRNA Expression Profiles in Colon Tumor Stroma

The aberrant expression of microRNAs (miRNAs) and genes in tumor microenvironment (TME) has been associated with the pathogenesis of colon cancer. An integrative exploration of transcriptional markers (gene signatures) and miRNA–mRNA regulatory networks in colon tumor stroma (CTS) remains lacking. Using two datasets of mRNA and miRNA expression profiling in CTS, we identified differentially expressed miRNAs (DEmiRs) and differentially expressed genes (DEGs) between CTS and normal stroma. Furthermore, we identified the transcriptional markers which were both gene targets of DEmiRs and hub genes in the protein–protein interaction (PPI) network of DEGs. Moreover, we investigated the associations between the transcriptional markers and tumor immunity in colon cancer. We identified 17 upregulated and seven downregulated DEmiRs in CTS relative to normal stroma based on a miRNA expression profiling dataset. Pathway analysis revealed that the downregulated DEmiRs were significantly involved in 25 KEGG pathways (such as TGF-β, Wnt, cell adhesion molecules, and cytokine–cytokine receptor interaction), and the upregulated DEmiRs were involved in 10 pathways (such as extracellular matrix (ECM)-receptor interaction and proteoglycans in cancer). Moreover, we identified 460 DEGs in CTS versus normal stroma by a meta-analysis of two gene expression profiling datasets. Among them, eight upregulated DEGs were both hub genes in the PPI network of DEGs and target genes of the downregulated DEmiRs. We found that three of the eight DEGs were negative prognostic factors consistently in two colon cancer cohorts, including COL5A2, EDNRA, and OLR1. The identification of transcriptional markers and miRNA–mRNA regulatory networks in CTS may provide insights into the mechanism of tumor immune microenvironment regulation in colon cancer.