Colorectal cancer consensus molecular subtypes (CMSs) are widely accepted and constitutes the basis for patient stratification to improve clinical practice. We aimed to find whether miRNAs could reproduce molecular subtypes, and to identify miRNA targets associated to the High-stroma/CMS4 subtype. The expression of 939 miRNAs was analyzed in tumors classified in CMS. TALASSO was used to find gene-miRNA interactions. A miR-mRNA regulatory network was constructed using Cytoscape. Candidate gene-miR interactions were validated in 293T cells. Hierarchical-Clustering identified three miRNA tumor subtypes (miR-LS; miR-MI; and miR-HS) which were significantly associated (p < 0.001) to the reported mRNA subtypes. miR-LS correlated with the low-stroma/CMS2; miR-MI with the mucinous-MSI/CMS1 and miR-HS with high-stroma/CMS4. MicroRNA tumor subtypes and association to CMSs were validated with TCGA datasets. TALASSO identified 1462 interactions (p < 0.05) out of 21,615 found between 176 miRs and 788 genes. Based on the regulatory network, 88 miR-mRNA interactions were selected as candidates. This network was functionally validated for the pair miR-30b/SLC6A6. We found that miR-30b overexpression silenced 3′-UTR-SLC6A6-driven luciferase expression in 293T-cells; mutation of the target sequence in the 3′-UTR-SLC6A6 prevented the miR-30b inhibitory effect. In conclusion CRC subtype classification using a miR-signature might facilitate a real-time analysis of the disease course and treatment response.

BACKGROUND: MicroRNAs (miRs) are frequently altered in colorectal cancer (CRC) and can be used as prognostic factors. OBJECTIVE: To confirm in stage III CRC patients a reported miR signature that was associated to the presence of metastatic disease. To correlate miR expression with microsatellite instability (MSI) and mutations in RAS and BRAF. METHODS: miR-21, miR-135a, miR-206, miR-335 and miR-Let-7a expression was analyzed by RT-qPCR in 150 patients out of the 329 patients used to analyze MSI and RAS and BRAF mutations. Association with disease free survival (DFS) and overall survival (OS) was analyzed. Data was confirmed by a multivariate analysis. RESULTS: MiR-21 high expression (p= 0.034) and miR-335 low expression (p= 0.0061) were significantly associated with MSI-H. A positive trend (p= 0.0624) between miR-135a high expression and RAS mutations was found. Lower miR-21 expression levels are associated with DFS (HR = 2.654, 95% CI: 1.066–6.605, p= 0.036) and a trend with OS (HR = 2.419, 95% CI: 0.749–7.815, p= 0.140). MiR-21 high expression significantly improves DFS of the poor prognosis group (T4 or N2) (p= 0.03). CONCLUSIONS: Association of increased expression of miR-21 and better prognosis in the poor prognostic group may be of interest and could be explored in future prospective clinical trials.