Letter to the Editor | Open | Published:
EphB2 stem-related and EphA2 progression-related miRNA-based networks in progressive stages of CRC evolution: clinical significance and potential miRNA drivers
Molecular Cancervolume 17, Article number: 169 (2018)
EphB2 and EphA2 control stemness and differentiation in the intestinal mucosa, but the way they cooperate with the complex mechanisms underlying tumor heterogeneity and how they affect the therapeutic outcome in colorectal cancer (CRC) patients, remain unclear. MicroRNA (miRNA) expression profiling along with pathway analysis provide comprehensive information on the dysregulation of multiple crucial pathways in CRC.
Through a network-based approach founded on the characterization of progressive miRNAomes centered on EphA2/EphB2 signaling during tumor development in the AOM/DSS murine model, we found a miRNA-dependent orchestration of EphB2-specific stem-like properties in earlier phases of colorectal tumorigenesis and the EphA2-specific control of tumor progression in the latest CRC phases. Furthermore, two transcriptional signatures that are specifically dependent on the EphA2/EphB2 signaling pathways were identified, namely EphA2, miR-423-5p, CREB1, ADAMTS14, and EphB2, miR-31-5p, mir-31-3p, CRK, CXCL12, ARPC5, SRC.
EphA2- and EphB2-related signatures were validated for their expression and clinical value in 1663 CRC patients. In multivariate analysis, both signatures were predictive of survival and tumor progression.
The early dysregulation of miRs-31, as observed in the murine samples, was also confirmed on 49 human tissue samples including preneoplastic lesions and tumors. In light of these findings, miRs-31 emerged as novel potential drivers of CRC initiation.
Our study evidenced a miRNA-dependent orchestration of EphB2 stem-related networks at the onset and EphA2-related cancer-progression networks in advanced stages of CRC evolution, suggesting new predictive biomarkers and potential therapeutic targets.
EphA2 and EphB2 are tyrosine kinase receptors that are involved in complex biochemical mechanisms underlying tumor heterogeneity primarily by guiding the positioning of intestinal cells with distinct stemness or differentiation properties . EphA2high cells are localized in the upper differentiated crypt region, while EphB2 is expressed in a decreasing gradient from the crypt base toward the upper cell compartment and corresponds to an intestinal stem cell (ISC) marker . At the onset of colorectal cancer (CRC), early events trigger continuous stem-like self-renewing state, giving rise to adenoma, while additional driver pathway alterations confer invasive behavior in advanced carcinoma . EphA2 and EphB2 undergo progressive dysregulation in carcinogenesis that may resemble these different stages, so representing attractive druggable targets .
MicroRNA (miRNA) redundancy and the capacity of individual miRNAs to simultaneously regulate large cohorts of genes have shaped by evolution combinatorial miRNA-target networks that profoundly affect cellular properties, including the promotion of tumorigenic processes . A network-based approach founded on the characterization of progressive miRNAomes centered on EphA2/EphB2 signaling during tumor development may unveil important regulatory networks and oncogenic targets.
Given that EphA2 and EphB2 have opposite distribution and roles along the intestinal crypt, we hypothesized that these receptors may influence alternative cell behaviors on initial or advanced CRC stages. Through an integrative translational approach (Additional File 1: Figure S1) and a miRNAome-guided pathway analysis, we defined two transcriptional signatures that are associated with EphB2 cells/early CRC phases and EphA2 cells/late CRC phases, with significant prognostic value.
Results and discussion
Distinct miRNAs characterize progressive stages of CRC development
Given the poorly foreseeable multiphase pattern of human CRC (hCRC) development, we used the well-characterized azoxymethane (AOM)/dextran sodium sulfate (DSS) murine model of sporadic CRC to obtain a predictable neoplastic evolution through the colorectal “Aberrant crypt foci (ACF)-microadenoma-adenoma-carcinoma” sequence, where each phase was expected to most likely progress to the following one [6, 7] (Fig. 1a; Additional file 1: Figure S2).
Of 641 miRNAs analyzed by TaqMan Low-Density Array (TLDA), 121 miRNAs were differentially expressed in the four murine CRC phases relative to normal mucosa (P < 0.05) (Fig. 1b, c, and d; Additional file 1: Figures S3 and S4). We found that the AOM/DSS model not only recapitulates several genetic modifications of hCRC , but also showed the aberrant expression of miRNAs belonging to hCRC-related pathways. In particular, miR-135b, targeting APC/β-catenin pathway , was up-regulated in murine adenoma (P < 0.05); let-7 family and miR-143, acting as tumor-suppressors in KRAS-driven CRC , were down-regulated already in ACF (let-7a, c, d, f P < 0.05) and in microadenoma (miR-143 P < 0.005); miR-9 was down-regulated in all the stages (P < 0.05), correlating with its tumor-suppressor function in hCRC by mediating invasion and metastasis  (Fig. 1c, d). More importantly, miR-31-5p and miR-31-3p were the most up-regulated in the earliest phases (P < 0.05 in ACF) (Fig. 1c, d), suggesting that these miRNAs may be involved in the CRC initiation.
Distinct EphA2/EphB2-specific miRNAomes orchestrate cell proliferation/stemness/migration in progressive murine CRC stages
Based on cytofluorimetric, immunohistochemical and molecular characterization (Additional file 2: Methods and Materials) of EphA2high/low and EphB2high/low cells isolated from murine adenocarcinoma and normal colon (Additional file 1: Figure S5), we ascertained that EphA2 and EphB2 mark murine tumor cell subpopulations bearing either a differentiated or a stem-like phenotype, respectively.
With this assumption, we performed a TLDA-based miRNA expression analysis to characterize miRNAomes related to EphA2high/low and EphB2high/low CRC cells (Additional file 3: Tables S1, S2).
A list of 180 miRNAs was associated with EphA2 and CRC. Moreover, according to the EphB2high CRC cells well-accepted identity of putative CSCs , we found they were enriched with miRNAs dysregulated in colon CSCs, such as miR-137, miR-34a, miR-215, miR-328 and miR-203 (P < 0.0001) .
Figure 1e describes a functional enrichment of the EphA2/EphB2 differential miRNAomes. All the functions analyzed were principally under control of miRNA targets listed in Additional file 3: Table S3.
When compared with the miRNAomes of murine ACF, microadenoma, adenoma and adenocarcinoma, the EphA2 and EphB2 tumor cells’ miRNAomes resulted in the alternate control of specific functions in different CRC stages (Fig. 1f; Additional file 3: Table S4).
It is worth noting that functions concerning cell morphology/cell differentiation and the inactivation of cell proliferation and cell survival/death were exclusively associated with EphB2 and enriched in ACF. In contrast, cell survival/death and cell proliferation, cell migration and metastasis and the reduction of cell adhesion were significantly associated with EphA2 and enriched in adenoma and adenocarcinoma.
These results indicate that while ACF might be enriched with EphB2high tumor-initiating cells, which exhibit driver-gene/miRNA influencing CRC initiation via stem-related functions, the latest stages of murine CRC show an enrichment with EphA2high cells orchestrating functions primarily related to tumor malignancy.
Establishment of two CRC EphA2/EphB2-specific signatures
Considering the EphA2/EphB2 macro-pathway (Fig. 2a), we compared Eph pathway-related miRNAs included in our preclinical results to those dysregulated in CRC patients (The Cancer Genome Atlas - colon adenocarcinoma, TCGA-COAD). Eight miRNAs were concordantly dysregulated in both murine and human CRC (Fig. 2b). Strikingly, among them, miR-31-5p and miR-31-3p were up-regulated in all four murine CRC stages and in the EphB2high tumor cells. Conversely, miR-423-5p was down-regulated in the advanced murine CRC phases and in the EphA2high tumor cells. These miRNAs displayed inverse regulation in the two cell subpopulations (Fig. 2b), supporting the strong association of miRs-31 and miR-423-5p not only with specific CRC phases but also with specific Eph receptor-expressing cells.
The pathway analysis of 617 genes target of miRs-31 and miR-423-5p (Additional file 3: Tables S5 and S6) resulted to significantly enrich CRC-related signaling pathways (P value < 0.05) (Additional file 1: Figures S6 and S7).
On that basis, we defined two EphA2-/EphB2-related signatures by the combination of miRs-31 and miR-423-5p with their experimentally validated target genes belonging to the EphA2/EphB2 macro-pathway. Therefore, the EphB2 signature comprehended miRs-31, ROHA, CRKL, CRK, ITGA5, CXCL12, ARPC5, RASA1 and SRC, while the EphA2 signature was composed of miR-423-5p, ARPC1A, PAK1, GNG12, CREB1 and ADAMTS14.
Up-regulation of miRs-31 occurs at very early stages of hCRC development
Given the importance of the early dysregulation of both miR-31-5p and miR-31-3p in the murine CRC, for their possible involvement in driving tumor initiation, we focused on the validation of their expression levels in the earliest stages of hCRC. We used a selected set of 49 human samples including low dysplasia (LD)-adenoma (n = 14), severe dysplasia (SD)-adenoma (n = 10), adenocarcinoma (n = 13) and normal mucosa (n = 12), that are not available in public datasets. qPCR analysis confirmed a significant overexpression of both miRs-31 in all the samples with an upward trend from LD-adenoma (preneoplastic) to SD-adenoma and adenocarcinoma, relative to normal mucosa (Student’s t-test, P < 0.05) (Fig. 2c). This is in line with recent demonstrations of miR-31-5p acting as a master modulator of the ISCs niche signaling during normal homeostasis, regeneration and tumorigenesis  or as potential driver of lung cancer initiation via KRAS pathway .
Prognostic value of EphA2/EphB2-specific signatures
Both EphA2 and EphB2 signatures were validated in 1663 CRC patients from nine cohorts of TCGA-COAD and Gene Expression Omnibus (GEO) databases, through a differential expression analysis between tumor and normal samples (Additional file 3: Table S7; Table 1a). EPHA2 was overexpressed in different datasets (P < 0.02) and miR-423-5p was concomitantly down-regulated (P < 0.01), whereas CREB1 and ADAMTS14 were overexpressed (P < 0.05). However, no significant changes in the expression levels of both EPHA2 and miR-423-5p were observed when they were analyzed in correlation with tumor staging progression (Additional file 1: Figure S8). EPHB2 was significantly overexpressed in TCGA-COAD (P = 0.03), GSE35834 cohort (P = 0.0008) and in patients with colorectal adenoma of GSE4183 cohort (P < 0.001); miR-31-5p was also overexpressed (P < 0.03) and its targets CRK, CXCL12, ARPC5 and SRC were down-regulated in three out of six cohorts. It is worth noting that in the analysis of the GSE4183 cohort, the expression level of EphB2 decreased in the transition from adenoma to carcinoma. In line with this observation, we found that EPHB2 expression levels gradually decreased with the progression of the TNM/AJCC stage and with metastasis (Fig. 2d). These results support the model that during CRC progression, stem EphB2-dependent functions are confined to a very small cell fraction of the tumor mass [7, 14], which likely represents the tumorigenic cell reservoir. Moreover, a gradual increase of miR-31-5p expression levels was observed with the progression of the TNM/AJCC stage (Fig. 2d), suggesting that it could be associated not only with particular miRNA patterns of stem-like EphB2high cells during tumorigenesis, but also with miRNA patterns of different cell populations replenishing the tumor mass in progressive CRC stages. Our findings are in line with functional studies showing that miR-31 has pleiotropic activity and is typically overexpressed with high expression correlating with advanced CRC disease .
In addition, we evaluated the collective prognostic effect of the EphA2 and EphB2 signatures (multivariate Cox-regression analysis) (Table 1b). Regarding the EphA2 fingerprint, ADAMTS14 was a risk factor for disease-free survival (DFS) (P < 0.001) and metastasis or recurrence within 3 years (MR3) (P < 0.05). Concerning the EphB2 fingerprint, SRC was a risk factor for overall survival (OS) (P < 0.01), CXCL12 was significant for DFS (P < 0.01), CRK, CXCL12 and ARPC5 were significant for MR3 (P < 0.05). Furthermore, miR-31-5p was deemed a risk factor in TCGA-COAD (P < 0.001). Our results confirm the involvement of these molecules in a complex network affected by EphA2/EphB2 signaling pathways and their crosstalk.
EphA2 and EphB2 contribute, in a pleiotropic manner, to the CRC pathogenesis. We demonstrated that the early dysregulation of EphB2 at the CRC onset correlates with specific stem-like properties and that the late overexpression of EphA2 triggers tumor progression signals in the invasive CRC phase. These properties are orchestrated by distinct miRNAomes and are associated with miR-31-5p/miR-31-3p, miR-423-5p and their EphA2/EphB2-related target genes. These signatures have significant clinical value and may unveil clinical predictive biomarkers and potential therapeutic markers.
Aberrant crypt foci
Cancer stem cells
Dextran sodium sulfate
Gene Expression Omnibus
Intestinal stem cells
Metastasis or recurrence within 3 years
The Cancer Genome Atlas - colon adenocarcinoma
TaqMan Low-Density Array
Kania A, Klein R. Mechanisms of ephrin-Eph signalling in development, physiology and disease. Nat Rev Mol Cell Biol. 2016;17:240–56.
Merlos-Suárez A, Barriga FM, Jung P, Iglesias M, Céspedes MV, Rossell D, et al. The intestinal stem cell signature identifies colorectal cancer stem cells and predicts disease relapse. Cell Stem Cell. 2011;8:511–24.
Tauriello DVF, Calon A, Lonardo E, Batlle E. Determinants of metastatic competency in colorectal cancer. Mol Oncol. 2017;11:97–119.
Lodola A, Giorgio C, Incerti M, Zanotti I, Tognolini M. Targeting Eph/ephrin system in cancer therapy. Eur J Med Chem. 2017;142:152–62.
Bracken CP, Scott HS, Goodall GJ. A network-biology perspective of microRNA function and dysfunction in cancer. Nat Rev Genetics. 2016;17:719–32.
De Robertis M, Massi E, Poeta ML, Carotti S, Morini S, Cecchetelli L, et al. The AOM/DSS murine model for the study of colon carcinogenesis: from pathways to diagnosis and therapy studies. J Carcinog. 2011;10:9.
De Robertis M, Loiacono L, Fusilli C, Poeta ML, Mazza T, Sanchez M, et al. Dysregulation of EGFR pathway in EphA2 cell subpopulation significantly associates with poor prognosis in colorectal Cancer. Clin Cancer Res. 2017;23:159–70.
Nagel R, le Sage C, Diosdado B, van der Waal M, Oude Vrielink JA, Bolijn A, et al. Regulation of the adenomatous polyposis coli gene by the miR-135 family in colorectal cancer. Cancer Res. 2008;68:5795–802.
Kent OA, Chivukula RR, Mullendore M, et al. Repression of the miR-143/145 cluster by oncogenic Ras initiates a tumor-promoting feed-forward pathway. Genes Dev. 2010;24:2754–9.
Park YR, Lee ST, Kim SL, Liu YC, Lee MR, Shin JH, et al. MicroRNA-9 suppresses cell migration and invasion through downregulation of TM4SF1 in colorectal cancer. Int J Oncol. 2016;48:2135–43.
De Robertis M, Poeta ML, Signori E, Fazio VM. Current understanding and clinical utility of miRNAs regulation of colon cancer stem cells. Semin Cancer Biol. 2018. https://doi.org/10.1016/j.semcancer.2018.08.008. [Epub ahead of print].
Tian Y, Ma X, Lv C, Sheng X, Li X, Zhao R, et al. Stress responsive miR-31 is a major modulator of mouse intestinal stem cells during regeneration and tumorigenesis. Elife. 2017;6:e29538. https://doi.org/10.7554/eLife.29538.
Edmonds MD, Boyd KL, Moyo T, Mitra R, Duszynski R, Arrate MP, et al. MicroRNA-31 initiates lung tumorigenesis and promotes mutant KRAS-driven lung cancer. J Clin Invest. 2016;126:349–64.
Zeuner A, Todaro M, Stassi G, De Maria R. Colorectal cancer stem cells: from the crypt to the clinic. Cell Stem Cell. 2014;15:692–705.
Schee K, Boye K, Abrahamsen TW, Fodstad O, Flatmark K. Clinical relevance of microRNA miR-21, miR-31, miR-92a, miR-101, miR-106a and miR-145 in colorectal cancer. BMC Cancer. 2012;12:505.
This work was partially supported by grants from Italian Ministry of Education, University and Research (MIUR) national research program: MIUR RBAP10KJC5 to ALV, MIUR PON02_00576_3329762/3 AMIDERHA to VMF, and MIUR PON03PE_00232_Distretto H-Bio to VMF and GP, and Italian Association for Cancer Research: AIRC IG 14368 to ALV and AIRC MFAG 10520 to MLP. This work was also partially supported by the Italian Ministry of Health, IRCCS Ricerca Corrente RC1803LO36 to VMF. Bioinformatics analyses carried out by GP were performed using the infrastructures of the Italian node of ELIXIR funded by MIUR and Consiglio Nazionale delle Ricerche (CNR).
Availability of data and materials
The authors declare that the data supporting the findings of this study are available within the article and its supplementary information files.
Ethics approval and consent to participate
Use and treatment of mice was approved by the Italian Ministry of Health with decree no. 336/2013-B.
The collection of clinical samples was approved by the Ethics Committee of the IRCCS “Regina Elena” National Cancer Institute and all patients signed an informed consent in accordance with the WMA Declaration of Helsinki 2013.
Consent for publication
All authors declare that they have no competing interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.