- Open Access
RUNX1 promotes tumour metastasis by activating the Wnt/β-catenin signalling pathway and EMT in colorectal cancer
Journal of Experimental & Clinical Cancer Researchvolume 38, Article number: 334 (2019)
Runt-related transcription factor 1 (RUNX1) plays the roles of an oncogene and an anti-oncogene in epithelial tumours, and abnormally elevated RUNX1 has been suggested to contribute to the carcinogenesis of colorectal cancer (CRC). However, the mechanism remains unclear.
The expression of RUNX1 in CRC and normal tissues was detected by real-time quantitative PCR and Western blotting. The effect of RUNX1 on CRC migration and invasion was conducted by functional experiments in vitro and in vivo. Chromatin Immunoprecipitation assay verified the direct regulation of RUNX1 on the promoter of the KIT, which leads to the activation of Wnt/β-catenin signaling.
RUNX1 expression is upregulated in CRC tissues. Upregulated RUNX1 promotes cell metastasis and epithelial to mesenchymal transition (EMT) of CRC both in vitro and in vivo. Furthermore, RUNX1 can activate Wnt/β-catenin signalling in CRC cells by directly interacting with β-catenin and targeting the promoter and enhancer regions of KIT to promote KIT transcription. These observations demonstrate that RUNX1 upregulation is a common event in CRC specimens and is closely correlated with cancer metastasis and that RUNX1 promotes EMT of CRC cells by activating Wnt/β-catenin signalling. Moreover, RUNX1 is regulated by Wnt/β-catenin.
Our findings first demonstrate that RUNX1 promotes CRC metastasis by activating the Wnt/β-catenin signalling pathway and EMT.
Colorectal cancer (CRC) is one of the most common human malignancies and is a leading cause of cancer-related death worldwide . Wnt/β-catenin signalling plays a vital role during development and in maintaining homeostasis in multiple tissues throughout the body, but Wnt/β-catenin signalling also plays an essential role in the initiation and early progression, as well as in later stages of invasion and metastasis, of CRC [2,3,4,5]. The process known as epithelial to mesenchymal transition (EMT) is a key mechanism in the invasiveness and metastatic drive in most cases of carcinoma . Moreover, dysregulation of the Wnt/β-catenin signalling pathway has been shown to play a significant role in the EMT required for colorectal tumour metastasis [7, 8].
Runt-related transcription factor 1 (RUNX1), also called acute myeloid leukaemia 1, is a member of RUNX family of transcription factors (RUNX1, RUNX2 and RUNX3), and this family is composed of evolutionarily conserved transcription factors that function as critical lineage determinants in various tissues [9,10,11]. RUNX1 is one of the most frequently mutated genes in a variety of haematological malignancies and has been proposed to play tumour suppressor roles in leukaemia; however, more recent studies suggest that wild-type RUNX1 is required for the growth and survival of certain types of leukaemia cells [12, 13]. Additionally, RUNX1 plays the role of an oncogene and an anti-oncogene in epithelial tumours, and its oncogenic effect has been given increasing attention [10, 14,15,16,17,18,19].
In previous studies, RUNX1 and the Wnt/β-catenin signalling pathway were shown to be closely related, and we found that their relationship involved inhibition  or interdependence [17, 20,21,22,23]. The role of RUNX1 in CRC has been examined in several studies [24, 25], and its possible mechanism involves inducing EMT by regulating the Wnt signalling pathway.
In the present study, we determined that the expression of RUNX1 was significantly upregulated in CRC tissues compared to its expression in adjacent tissues. We showed that RUNX1 is a critical activator of CRC cell metastasis both in vitro and in vivo. Moreover, RUNX1 activated the Wnt/β-catenin signalling pathway to promote EMT and tumour metastasis.
Materials and methods
A total of 73 patients who underwent radical operation for CRC at Nanfang Hospital of Southern Medical University were included in the study after obtaining informed consent. A diagnosis of CRC was histopathologically confirmed for each patient sample. Cancer tissues and matched normal tissues were stored at -80C° for quantitative real-time PCR (qRT-PCR) and Western blotting analyses. Additionally, none of the patients had received radiotherapy or chemotherapy before surgery. The protocols used in this study were approved by Nanfang Hospital’s Protection of Human Subjects Committee.
Cell culture, plasmid construction, lentiviral construction and cell transfection
Human embryonic kidney cell line (293 T), human normal colon epithelial cell line (FHC) and six human CRC cell lines (SW480, SW620, RKO, HCT116 HT29 and LOVO) were purchased from the Cell Bank of Type Culture Collection (CBTCC, Chinese Academy of Sciences, Shanghai, China) and were cultured in Dulbecco’s modified Eagle medium (DMEM) (Gibco, Carlsbad, CA) supplemented with 10% fetal bovine serum (FBS; Gibco, Carlsbad, CA). Cells were maintained at 37 °C in a humidified 5% CO2 atmosphere.
Plasmid construction, lentiviral construction and cell transfections
RUNX1 overexpression and knockdown were performed using a lentiviral packaging system. To construct overexpressing exogenous and RNA-interfered endogenous RUNX1 cell lines, full-length RUNX1 (NM_001754) was cloned into the expression vector pLenti-EF1a-EGFP-P2A-Puro-CMV (Obio Technology, Shanghai, China) and transfected into HCT116 and RKO cell lines according to the manufacturer’s instructions. Knockdown of endogenous RUNX1 was mediated by designed shRNAs (Cyagen, Guangzhou, China) that were transfected into SW480 and RKO cell lines according to the manufacturer’s instructions. The RUNX1 shRNA sequence was sense 5′-CCAGGTTGCAAGATTTAAT-3′ and 5′-GGCAGAAACTAGATGATCA-3′, and the scramble sequence was sense 5′-CCTAAGGTTAAGTCGCCCTCG-3′. Transduced cells were selected in medium containing puromycin (#EZ2811D376, BioFroxx, Germany) (2 μg/ml) and maintained in medium containing puromycin (1 μg/ml).
RNA isolation and qRT-PCR
Total RNA was extracted from cells or tissues with TRIzol reagent (TaKaRa, China). qRT-PCR was performed using the PrimeScript RT Reagent Kit (#RR035A, TaKaRa, China) and SYBR Premix Ex Taq (#RR820A, TaKaRa, Dalian, China) following the manufacturer’s instructions. Our results were normalized to the expression of glyceraldehyde-3-phosphate dehydrogenase (GAPDH) or U6. The specific primers used are listed in Additional file 5: Table S1. The qRT-PCR results were analyzed to obtain the Ct values of the amplified products, and data were analyzed by the 2-ΔΔCt method.
Western blot and immunohistochemistry (IHC) analysis
We performed Western blotting according to the methods of a previous study. Protein lysates were prepared, subjected to sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE), transferred onto polyvinylidene difluoride (PVDF) membranes and blotted according to standard methods using the following antibodies: RUNX1 (#4336, CST), KIT (#3074, CST), CD44 (#3570, CST), cyclin D1 (#2978, CST), c-Jun (#9165, CST), LEF1 (#2230, CST), Met (#8198, CST), c-Myc (#5605, CST), TCF1/TCF7 (#2230, CST), β-catenin (#8480, CST), E-cadherin (#3195, CST), N-cadherin (#13116, CST), Vimentin (#5741, CST), TCF8/ZEB1 (#3396, CST), AXIN1 (#2087, CST) and GAPDH (60004–1-Ig, Proteintech).
Immunohistochemistry was performed following the manufacturer’s instructions (PV-6001, ZSGB-BIO, Beijing, China) using the E-cadherin antibody (60335–1-Ig, Proteintech), N-cadherin (66219–1-Ig, Proteintech), and Vimentin (60330–1-Ig, Proteintech). One independent pathologist used software ImageJ to calculate gray values for pathological scoring.
Cell Transwell migration and invasion assays
Approximately 2 × 105 cells were suspended in 100 μl serum-free medium and seeded in the upper 8-μm-pore Transwell chambers (3422, Corning, USA), and the lower chambers were filled with 500 μL of 20% FBS medium. For the cell invasion assay, the Transwell chambers were Matrigel-coated (354234, Corning, USA). Then, the cells were incubated at 37 °C for 24–48 h to allow for migration or invasion. For quantification, the cells were fixed with 4% paraformaldehyde for 20 min, stained with hematoxylin for 20 min at room temperature, and counted in five randomly chosen fields (× 200) under a microscope.
Cell wound healing assay
Cells were seeded on six-well culture plates at a density of approximately 1 × 106 cells per well and incubated for 24 h (80–90% confluence). Scratch wounds were produced by a 10 μl plastic pipette tip, and then cells were cultured in DMEM with 2% FBS. Wound margins were photographed, and migration was monitored after 48 h of wound formation. Cell motility was quantified by measuring the distance between the advancing margins of cells in three randomly selected microscopic fields (× 200) at each time point.
For immunofluorescence, cells were seeded on cover slips. Overexpression or shRNA knockdown treatment was performed after 24 h. After the indicated treatment, the cells were cultured for 48 h, fixed with 4% formaldehyde for 10 min at room temperature, and washed three times with wash buffer (0.02% Tween 20/PBS). Then, the cells were permeabilized with 0.5% Triton X-100/PBS for 10 min at room temperature. The cells were washed three times with wash buffer (5 min each time) and then incubated with 1.5% bovine serum albumin (BSA)/phosphate-buffered saline (PBS) solution (blocking solution) for 30 min at room temperature. The E-cadherin antibody (60335–1-Ig, Proteintech), N-cadherin (66219–1-Ig, Proteintech), Vimentin (60330–1-Ig, Proteintech), and β-catenin (66379–1-Ig, Proteintech) were incubated in blocking solution at 4 °C overnight. Rhodamine phalloidin (1:1500, Cytoskeleton, PHDR1), which was used for detecting F-actin, was incubated in blocking solution at room temperature in the dark for 60 min. After washing, the cells were incubated with Alexa594-conjugated secondary antibodies for 60 min at room temperature in a protected environment in the dark (Life Technologies, A-21235, 1:500 in blocking buffer) followed by counterstaining with DAPI (Thermo Fisher). Samples were mounted with ProLong Gold antifade reagent (Life Technologies) and imaged on a confocal microscope (× 200).
ChIP assays were performed using a kit (#17–10085, Merck, German), and all the experimental procedures were performed according to the manufacturer’s instructions. Briefly, cells (1 × 107) in a 15-cm culture dish were treated with 1% formaldehyde to cross-link chromatin-associated proteins to DNA. The cell lysates were subjected to ultrasound for 9–10 sets of 10-s pulses at 40% output to shear the DNA into fragments between 200 and 1000 bps. Equal cell lysates were incubated with 10 μl of DYKDDDDK Tag (bound to the same epitope as the Sigma Anti-Flag M2 antibody) (#14793, CST), anti-IgG antibody (Merck) as a negative control and anti-RNA polymerase II as a positive control. All the above chromatin supernatants were incubated with 20 μL magnetic protein A/G beads overnight at 4 °C with rotation. On the second day, the protein-DNA complexes were reversed and purified for pure DNA. The human c-Kit binding sites were amplified with qRT-PCR and PCR. The specific primers used are listed in Additional file 6: Table S2.
Orthotopic injection metastatic mouse model
Female athymic 4 to 5-week-old BALB/c nude mice were purchased from the Laboratory Animal Services Centre of Guangdong Province and were maintained in a specific pathogen-free facility. For the orthotopic injection metastatic mouse model assay, the nude mice of the experimental group were anesthetized with 1% pentobarbital, and 5 × 106 cells were injected under the ileocecal serosa. The general growth conditions in the nude mice were observed after surgery to determine whether there were signs of cachexia, such as wasting and arched backs. Then, 40 days after the operation, all mice were sacrificed, and the intestines, liver, spleen and lungs were removed for photographing. After formalin fixation, paraffin sections were generated and stained with hematoxylin and eosin (H&E).
The SPSS 22.0 (IBM; Chicago, IL, USA) and Microsoft Excel 2016 (Microsoft, Redmond, WA, USA) statistical analysis software programs were used for statistical analysis of the experimental data. The significance of differences between groups was estimated by Student’s t-test. Additionally, multiple group comparisons were analyzed with one-way ANOVA. Bivariate correlations between study variables were calculated by Spearman’s rank correlation coefficients. * P < 0.05, ** P < 0.01, and *** P < 0.001 were considered significant.
RUNX1 is upregulated in CRC
An analysis of The Cancer Genome Atlas (TCGA) database using the online Gene Expression Profiling Interactive Analysis (GEPIA) tool (http://gepia.cancer-pku.cn/) indicated that RUNX1 is highly expressed in most types of digestive system tumours (Fig. 1a). To screen CRC gene expression profiles, we selected two qualified gene expression microarray datasets (TCGA and GSE106582) (Fig. 1b) and determined that RUNX1 was differentially expressed in colorectal tumour tissues compared to its expression in non-neoplastic tissues (p < 0.001) in two independent samples and in paired samples (Fig. 1c). An analysis of the correlation between CRC-related gene set enrichment and RUNX1 expression showed statistically significant (NES = 1.54, p = 0.006) results in TCGA database, and Gene Set Enrichment Analysis (GSEA) (http://software.broadinstitute.org/gsea/index.jsp) suggested a significant correlation between RUNX1 expression and these genes (Fig. 1d). Notably, Kaplan-Meier survival analysis revealed that patients with high RUNX1 expression levels had poorer disease-free survival and overall survival than patients with low RUNX1 expression levels. Similar results were found when the GEPIA web server was used to analyze TCGA database (Fig. 1e) as well as in the PROGgeneV2 - Pan Cancer Prognostics Database using the GEO database (http://genomics.jefferson.edu/proggene/) (Additional file 1: Figure S1A and B).
Furthermore, to determine whether RUNX1 protein expression is also increased in CRC, qRT-PCR and Western blotting were performed in 61 pairs of fresh samples and in 12 pairs of fresh samples, including cancers and matched adjacent normal tissues or independent normal tissues (Fig. 1f and g). In this analysis, the RUNX1 levels were significantly upregulated in most CRC tissues. In the patient clinical information statistics table (Additional file 7: Table S3), we found that patients with high expression of RUNX1 were worse in AJCC stage, T stage and M stage, when the patients were divided into high and low RUNX1 expression groups by the median expression of RUNX1 mRNA. Western blotting and qRT-PCR analysis revealed varying levels of RUNX1 expression in six CRC cell lines: SW480, SW620, HT29, HCT116, RKO and LoVo. RUNX1 expression was highest in SW480 and SW620 cells; lowest in HCT116, HT29 and LoVo cells; and moderately expressed in RKO cells (Additional file 1: Figure S1C).
RUNX1 promotes CRC cell metastasis in vitro and in vivo
To investigate the biological function of RUNX1 in CRC cells in vitro and in vivo, we established stable RUNX1-overexpressing CRC cell lines with HCT116 and RKO cells and RUNX1-silenced cell lines using a lentiviral vector carrying two specific shRNAs in SW480 and RKO cells. Forty-eight hours after transfection, qRT-PCR and Western blotting analysis showed a significant increase in RUNX1 mRNA and protein expression in HCT116-RUNX1 and RKO-RUNX1 cells and a significant decrease in RUNX1 expression in SW480-shRUNX1_1, SW480-shRUNX1_2, RKO-shRUNX1_1 and RKO-shRUNX1_2 cells compared to the expression levels in control cells (Fig. 2a).
Markedly faster wound closure was observed in HCT116-RUNX1 and RKO-RUNX1 cells, but SW480-shRUNX1_1, SW480-shRUNX1_2, RKO-shRUNX1_1 and RKO-shRUNX1_2 cells exhibited slower wound closure, as revealed by a wound healing assay (Fig. 2c and Additional file 2: Figure S2B). Moreover, the results from Transwell assays revealed that overexpression of RUNX1 significantly increased the migration and invasion abilities of CRC cells. Next, our results showed that the migration and invasion abilities were significantly decreased in cells with RUNX1 downregulation (Fig. 2b and Additional file 2: Figure S2A).
Furthermore, we also achieved consistent results for in vivo metastasis assays. As shown by our results, tumours formed by cells overexpressing RUNX1 showed stronger metastatic ability, and more tumours were formed in the spleen (Additional file 2: Figure S2C) in the nude mice cecum tumour orthotopic model transfected with RUNX1-overexpressing CRC HCT116-RUNX1 cells than in the control group transfected with the HCT116-vector (Fig. 2d).
RUNX1 promotes EMT
GSEA analysis using the TCGA and GSE17538 datasets showed that RUNX1 expression is associated with regulation of the actin cytoskeleton and focal adhesion (Fig. 3a). We retrieved genes with more than 100 RUNX1 transcription factor binding sites from the Gene Transcription Regulation Database (GTRD) (http://gtrd.biouml.org/) and used the results to conduct a Gene Ontology (GO) enrichment analysis using DAVID Bioinformatics Resources 6.8 (https://david.ncifcrf.gov/). The GO-Cellular Components enrichment analysis results showed gene enrichment in the cellular components of the cytoskeleton, cell junction, actin cytoskeleton and cell-cell junction (Fig. 3b). Next, immunofluorescence labelling of F-actin was performed using phalloidin. Immunofluorescence analysis showed that the cell lines with upregulated RUNX1 were characterized by a spindle-shaped morphology and the rearrangement of F-actin fibres (Fig. 3c). The spindle shape change and loss of cell-cell contact are characteristics of the cell morphology changes that occur during EMT. Through mRNA expression correlation analysis on website GEPIA, we found that the expression of RUNX1 is positively correlated with multiple molecules of the EMT and is statistically significant, such as N-cadherin (CDH2), Vimentin (VIM), Snail and ZEB1 (Additional file 3: Figure S3A). Western blotting was used to analyze EMT markers induced by RUNX1 in the HCT116, SW480 and RKO cell lines. E-cadherin was downregulated in HCT116-RUNX1 cells and upregulated in SW480-shRUNX1 cells compared to its expression in control cells. The upregulation of ZEB1, N-cadherin and vimentin was also induced by RUNX1 in HCT116-RUNX1 cells and inhibited in SW480-shRUNX1 cells (Fig. 3d). Consistent with these results, the immunofluorescence assay results revealed decreased expression of E-cadherin, whereas the expression levels of N-cadherin and vimentin were decreased in HCT116-RUNX1 cells, and the opposite result was observed in SW480-shRUNX1_2 cells (Fig. 3e). Consistent with these findings, decreased E-cadherin and increased N-cadherin and vimentin expression levels were shown by immunohistochemistry in situ in the RUNX1-overexpressing cecum tumours of nude mice (Fig. 3f). Through mRNA expression correlation analysis on website GEPIA, we found that the expression of RUNX1 is positively correlated with multiple molecules of the MMP family and is statistically significant, such as MMP3, MMP7, MMP9 and MMP14 (Additional file 3: Figure S3C). Subsequently, we also used qPCR experiments to detect elevated expression of the above molecules in the HCT116 cell line overexpressing RUNX1 (Additional file 4: Figure S4A). These results strongly suggest that RUNX1 is involved in enhancing the metastatic capacity of CRC.
RUNX1 enhances Wnt/β-catenin pathway activation in CRC
GSEA analysis using both TCGA and the GSE17538 database showed that RUNX1 expression is associated with the Wnt signalling pathway (Fig. 4a). Next, we found a positive correlation between the expression of RUNX1 and CTNNB1, also called β-catenin, in the gene expression microarray dataset (GSE17538) (Fig. 4b) and in CRC cell lines (Fig. 4c). Through mRNA expression correlation analysis on website GEPIA, we found that the expression of RUNX1 is positively correlated with multiple molecules of the Wnt/β-catenin signalling pathway and is statistically significant, such as CTNNB1, WNT5a, CD44 and LEF1 (Additional file 3: Figure S3B). Furthermore, we verified the expression of β-catenin by Western blot and immunofluorescence analysis in lentivirus-infected cell lines, and the results showed that β-catenin expression increased after RUNX1 was overexpressed and decreased after RUNX1 was downregulated (Fig. 4d and g). Western blot analysis showed that β-catenin was upregulated in the nucleus with RUNX1 overexpression and downregulated in the nucleus with RUNX1 knockdown (Fig. 4e). A co-immunoprecipitation assay indicated that RUNX1 directly interacts with β-catenin in SW480 cells (Fig. 4f). We suggest that RUNX1 and β-catenin possibly enter the nucleus together to perform their functions. We used Western blotting to measure downstream targets, such as Met, C-Jun, C-Myc, TCF1/7, CD44, cyclin D1 and LEF1 (Fig. 4h). Our results showed that RUNX1 increased the expression of proteins, including Met, C-Jun, C-Myc, TCF1/7, CD44, cyclin D1 and LEF1, that are associated with the Wnt/β-catenin signalling pathway.
RUNX1 targets KIT, and KIT knock down suppresses Wnt/β-catenin signalling
Using genes with more than 100 RUNX1 transcription factor binding sites, a GO-Biological Process enrichment analysis was performed using DAVID Bioinformatics Resources 6.8. The three most significant biological process involved positive regulation of GTPase activity, intracellular signal transduction and signal transduction (Fig. 5a). Through a Wayne diagram using the website (http://bioinfogp.cnb.csic.es/tools/venny/index.html), we determined the intersection of these three gene sets, which contained eight molecules (Fig. 5b), including KIT, PLCB1, BCR, ARHGEF7, KALRN, MY09B, ABR and CHN2. Using the protein interaction analysis website STRING (https://version-10-5.string-db.org/cgi/input.pl), we found that KIT may be a bridge between RUNX1 and the Wnt/β-catenin signalling pathway (Fig. 5c). Moreover, Western blot results showed that the expression of KIT increased and decreased along with the expression of RUNX1 (Fig. 5d). After siRNA was used to interfere with KIT expression, the expression of several proteins, including C-Myc and TCF1, that activate the Wnt/β-catenin signalling pathway was decreased (Fig. 5e) RUNX1 transcription factor binding sites have been reported in the promoter and enhancer regions of KIT . Therefore, we used ChIP-qPCR and ChIP-PCR to examine the binding sites at + 700 bp and + 30 kb from the transcription start site KIT, and the results showed that RUNX1 can bind at those two sites (Fig. 5f).
RUNX1 is affected by the Wnt/β-catenin signalling pathway
XAV-939, a Wnt/β-catenin signalling pathway specific inhibitor, promotes β-catenin degradation. Moreover, using Transwell assays, we found that when cells were treated with XAV-939, the migration and invasion ability of HCT116 and RKO RUNX1-overexpressing cell lines was decreased (Fig. 6a). At the same time, XAV-939 can inhibit the expression changes of EMT-related molecules induced by RUNX1 (Fig. 6b). We used qRT-PCR and Western blotting to examine RUNX1 expression after treatment with siRNA of β-catenin for 48 h, and the results showed that the expression of RUNX1 decreased when the β-catenin expression was disturbed. (Fig. 6c).
The roles of RUNX1 have been described in many different cancers; RUNX1 is overexpressed in various epithelial tumours, especially during tumour initiation . Some studies on ovarian carcinoma have suggested that RUNX1 contributes to cell proliferation, migration and invasion [27, 28]. In studies of oesophageal cancer, RUNX1a, a transcript of RUNX1, has been shown to promote cell proliferation and tumour growth . Interestingly, RUNX1 functions as both an oncogene and a tumour suppressor in breast cancer . RUNX1 loss promotes oestrogen receptor ER+ breast cancer epithelial cell growth and stem cell markers, but this growth promotion does not occur in ER- breast epithelial cell lines [14, 17, 18]. RUNX1 also plays dual role in gastric cancer [31,32,33]. Although the oncogenic role of RUNX1 in CRC has not been thoroughly examined, some studies have laterally confirmed its function [24, 34].
An unresolved question is whether RUNX1 functions to promote or suppress tumour metastasis in CRC. This report is the first direct investigation of the function and mechanism of action of RUNX1 in CRC. First, we found that RUNX1 expression was upregulated in CRC by detecting its levels in paired samples. Next, we identified RUNX1 as an oncogenic transcription factor that promotes CRC cell metastasis both in vitro and in vivo, and downregulation of RUNX1 had the opposite effect. Moreover, we found that RUNX1 promotes EMT in CRC. Mechanistically, RUNX1 promotes CRC metastasis by directly interacting with β-catenin and activating KIT transcription to promote β-catenin nucleation and the subsequent activation of the Wnt signalling pathway.
In this study, we found that RUNX1 expression was significantly increased in CRC tissues, which is consistent with findings from TCGA and the Gene Expression Omnibus (GEO) database. Furthermore, the expression of RUNX1 was positively associated with CRC metastasis. Therefore, it can be concluded that RUNX1 acts as an oncogene or risk factor in CRC tumorigenesis, which further confirms previous research [24, 25, 34]. The results of our in vitro and in vivo experiments showed that RUNX1 increased CRC the migration and invasion capacities of cell lines, and this effect is related to EMT.
It is noteworthy that some research has verified that the loss of RUNX1 expression in ER+ mammary epithelial cells increases β-catenin, deregulates mitosis and stimulates cell proliferation and expression of stem cell markers . However, in complementary experiments, we also found that RUNX1 could affect the expression of AXIN1 protein, but no clear positive correlation between the two molecules was observed (Additional file 4: Figure S4B). Objectively, other results suggest a positive correlation between RUNX and the Wnt/β-catenin signalling pathway [20,21,22,23]. The RUNX1 gene can be a Wnt4 signalling target, and RUNX1 and Wnt4 are mutually interdependent in their expression in ovary cells . A study involving leukaemia reported that Wnt/β-catenin signalling increases the expression of the ETO and RUNX1 genes in human haematopoietic progenitors , and the distal P1-RUNX1 promoter is a direct transcriptional target of Wnt/β-catenin signalling . All of the abovementioned results are consistent with those of this study. Treatment with the Wnt/β-catenin pathway inhibitor XAV-939 resulted in decreases in RUNX1 expression and cell migration ability. Canonical Wnt signalling was investigated in RUNX1-deficient bone marrow stem cells, and decreases in β-catenin, LEF1, TCF1 and Wnt10b expression were found in these cells . Our results also showed the corresponding conclusions that RUNX1 is positively correlated with CTNNB1 expression, and the Wnt/β-catenin pathway was activated when RUNX1 was overexpressed and inhibited when RUNX1 was downregulated.
In many cancers, Wnt/β-catenin signalling is constitutively active and promotes EMT . Nuclear β-catenin binds to members of the TCF/LEF family of transcription factors to promote EMT. LEF1 is another key transcription factor that can directly induce EMT by repressing E-cadherin [36, 37]. In this study, we found that RUNX1 promotes β-catenin expression and activates Wnt signalling, and some transcription factors, such as LEF1, which can repress E-cadherin, are also activated. Previous studies have shown that MMPs are associated with the occurrence of EMT , while MMPs are regulated by Wnt/β-catenin signalling pathways [39,40,41,42].
As a transcription factor, RUNX1 may regulate the transcription of target genes in many ways. Thirty percent of all RUNX1 binding sites are intergenic, indicating its diverse roles in promoter and enhancer regulation and suggesting additional functions of RUNX1 . We validated two strong RUNX1 binding regions by ChIP-qPCR and ChIP-PCR in the first intron of KIT: one at + 700 bp and another at + 30 kb from the transcription start site . Moreover, our results showed that the expression of KIT increased with the overexpression of RUNX1 and decreased with the reduced expression of RUNX1. KIT signalling may play a growth-stimulatory role in colon cancer . Furthermore, KIT activation may be associated with tumour aggressiveness via the activation of Wnt/β-catenin signalling . We observed that the Wnt pathway was inhibited when KIT expression was disturbed by siRNA transfection, which may be the mechanism by which RUNX1 regulates the Wnt/β-catenin signalling pathway.
In conclusion, the current study illustrates that RUNX1 functions as an oncogene to facilitate metastasis and EMT in CRC by directly interacting with β-catenin and activating KIT transcription to enhance the Wnt/β-catenin signalling pathway (Fig. 6d). These findings enhance our understanding of CRC metastasis. In addition, RUNX1 is downstream of the Wnt pathway and is regulated by Wnt/β-catenin. After further research, RUNX1 may become a vital prognostic biomarker and an effective target for anti-metastasis therapies for CRC.
In summary, our results highlight a crucial role for RUNX1 in the regulation of tumour metastasis in CRC by activating the Wnt/β-catenin signalling pathway and EMT, and RUNX1 might be regarded as a potential prognostic marker and as an effective therapeutic target for CRC.
Availability of data and materials
All data analysed during this study are included in this manuscript.
Liver hepatocellular carcinoma
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2017. CA Cancer J Clin. 2017;67:7–30.
Qi J, Yu Y, Akilli Ozturk O, Holland JD, Besser D, Fritzmann J, Wulf-Goldenberg A, Eckert K, Fichtner I, Birchmeier W. New Wnt/beta-catenin target genes promote experimental metastasis and migration of colorectal cancer cells through different signals. Gut. 2016;65:1690–701.
Zheng K, Zhou X, Yu J, Li Q, Wang H, Li M, Shao Z, Zhang F, Luo Y, Shen Z, et al. Epigenetic silencing of miR-490-3p promotes development of an aggressive colorectal cancer phenotype through activation of the Wnt/beta-catenin signaling pathway. Cancer Lett. 2016;376:178–87.
Liu N, Jiang F, Han XY, Li M, Chen WJ, Liu QC, Liao CX, Lv YF. MiRNA-155 promotes the invasion of colorectal cancer SW-480 cells through regulating the Wnt/beta-catenin. Eur Rev Med Pharmacol Sci. 2018;22:101–9.
Li J, He M, Xu W, Huang S. LINC01354 interacting with hnRNP-D contributes to the proliferation and metastasis in colorectal cancer through activating Wnt/beta-catenin signaling pathway. J Exp Clin Cancer Res. 2019;38:161.
Hu JL, Wang W, Lan XL, Zeng ZC, Liang YS, Yan YR, Song FY, Wang FF, Zhu XH, Liao WJ, et al. CAFs secreted exosomes promote metastasis and chemotherapy resistance by enhancing cell stemness and epithelial-mesenchymal transition in colorectal cancer. Mol Cancer. 2019;18:91.
Yang S, Liu Y, Li MY, Ng CSH, Yang SL, Wang S, Zou C, Dong Y, Du J, Long X, et al. FOXP3 promotes tumor growth and metastasis by activating Wnt/beta-catenin signaling pathway and EMT in non-small cell lung cancer. Mol Cancer. 2017;16:124.
Gu Y, Wang Q, Guo K, Qin W, Liao W, Wang S, Ding Y, Lin J. TUSC3 promotes colorectal cancer progression and epithelial-mesenchymal transition (EMT) through WNT/beta-catenin and MAPK signalling. J Pathol. 2016;239:60–71.
Levanon D, Groner Y. Structure and regulated expression of mammalian RUNX genes. Oncogene. 2004;23:4211–9.
Ito Y, Bae SC, Chuang LS. The RUNX family: developmental regulators in cancer. Nat Rev Cancer. 2015;15:81–95.
Ito Y. Oncogenic potential of the RUNX gene family: 'overview'. Oncogene. 2004;23:4198–208.
Choi A, Illendula A, Pulikkan JA, Roderick JE, Tesell J, Yu J, Hermance N, Zhu LJ, Castilla LH, Bushweller JH, Kelliher MA. RUNX1 is required for oncogenic Myb and Myc enhancer activity in T-cell acute lymphoblastic leukemia. Blood. 2017;130:1722–33.
Sood R, Kamikubo Y, Liu P. Role of RUNX1 in hematological malignancies. Blood. 2017;129:2070–82.
Hong D, Fritz AJ, Finstad KH, Fitzgerald MP, Weinheimer A, Viens AL, Ramsey J, Stein JL, Lian JB, Stein GS. Suppression of breast Cancer stem cells and tumor growth by the RUNX1 transcription factor. Mol Cancer Res. 2018;16:1952–64.
Han R, Wei J, Zhang H, Su X, Chu X, Chen Y, Gong Y, Wang X, Shi J, Chen J. Influence of TS (rs34743033) and RUNX1 (rs2014300) gene polymorphisms on survival outcomes of fluorouracil-based chemotherapy in Chinese advanced gastric cancer patients. Cancer Manag Res. 2018;10:1429–37.
Deltcheva E, Nimmo R. RUNX transcription factors at the interface of stem cells and cancer. Biochem J. 2017;474:1755–68.
Chimge NO, Ahmed-Alnassar S, Frenkel B. Relationship between RUNX1 and AXIN1 in ER-negative versus ER-positive breast Cancer. Cell Cycle. 2017;16:312–8.
Chimge NO, Little GH, Baniwal SK, Adisetiyo H, Xie Y, Zhang T, O'Laughlin A, Liu ZY, Ulrich P, Martin A, et al. RUNX1 prevents oestrogen-mediated AXIN1 suppression and beta-catenin activation in ER-positive breast cancer. Nat Commun. 2016;7:10751.
Scheitz CJ, Lee TS, McDermitt DJ, Tumbar T. Defining a tissue stem cell-driven Runx1/Stat3 signalling axis in epithelial cancer. EMBO J. 2012;31:4124–39.
Naillat F, Yan W, Karjalainen R, Liakhovitskaia A, Samoylenko A, Xu Q, Sun Z, Shen B, Medvinsky A, Quaggin S, Vainio SJ. Identification of the genes regulated by Wnt-4, a critical signal for commitment of the ovary. Exp Cell Res. 2015;332:163–78.
Ugarte GD, Vargas MF, Medina MA, Leon P, Necunir D, Elorza AA, Gutierrez SE, Moon RT, Loyola A, De Ferrari GV. Wnt signaling induces transcription, spatial proximity, and translocation of fusion gene partners in human hematopoietic cells. Blood. 2015;126:1785–9.
Medina MA, Ugarte GD, Vargas MF, Avila ME, Necunir D, Elorza AA, Gutierrez SE, De Ferrari GV. Alternative RUNX1 promoter regulation by Wnt/beta-catenin signaling in leukemia cells and human hematopoietic progenitors. J Cell Physiol. 2016;231:1460–7.
Luo Y, Zhang Y, Miao G, Zhang Y, Liu Y, Huang Y. Runx1 regulates osteogenic differentiation of BMSCs by inhibiting adipogenesis through Wnt/beta-catenin pathway. Arch Oral Biol. 2019;97:176–84.
Zhou Y, Zhang X, Zhang J, Fang J, Ge Z, Li X. LRG1 promotes proliferation and inhibits apoptosis in colorectal cancer cells via RUNX1 activation. PLoS One. 2017;12:e0175122.
Voora D, Rao AK, Jalagadugula GS, Myers R, Harris E, Ortel TL, Ginsburg GS. Systems pharmacogenomics finds RUNX1 is an aspirin-responsive transcription factor linked to cardiovascular disease and Colon Cancer. EBioMedicine. 2016;11:157–64.
Debaize L, Jakobczyk H, Avner S, Gaudichon J, Rio AG, Serandour AA, Dorsheimer L, Chalmel F, Carroll JS, Zornig M, et al. Interplay between transcription regulators RUNX1 and FUBP1 activates an enhancer of the oncogene c-KIT and amplifies cell proliferation. Nucleic Acids Res. 2018;46:11214–28.
Keita M, Bachvarova M, Morin C, Plante M, Gregoire J, Renaud MC, Sebastianelli A, Trinh XB, Bachvarov D. The RUNX1 transcription factor is expressed in serous epithelial ovarian carcinoma and contributes to cell proliferation, migration and invasion. Cell Cycle. 2013;12:972–86.
Eldholm V, Haugen A, Zienolddiny S. CTCF mediates the TERT enhancer-promoter interactions in lung cancer cells: identification of a novel enhancer region involved in the regulation of TERT gene. Int J Cancer. 2014;134:2305–13.
Wu H, Zheng J, Deng J, Zhang L, Li N, Li W, Li F, Lu J, Zhou Y. LincRNA-uc002yug.2 involves in alternative splicing of RUNX1 and serves as a predictor for esophageal cancer and prognosis. Oncogene. 2015;34:4723–34.
Barutcu AR, Hong D, Lajoie BR, McCord RP, van Wijnen AJ, Lian JB, Stein JL, Dekker J, Imbalzano AN, Stein GS. RUNX1 contributes to higher-order chromatin organization and gene regulation in breast cancer cells. Biochim Biophys Acta. 2016;1859:1389–97.
Mitsuda Y, Morita K, Kashiwazaki G, Taniguchi J, Bando T, Obara M, Hirata M, Kataoka TR, Muto M, Kaneda Y, et al. RUNX1 positively regulates the ErbB2/HER2 signaling pathway through modulating SOS1 expression in gastric cancer cells. Sci Rep. 2018;8:6423.
Zhuang M, Gao W, Xu J, Wang P, Shu Y. The long non-coding RNA H19-derived miR-675 modulates human gastric cancer cell proliferation by targeting tumor suppressor RUNX1. Biochem Biophys Res Commun. 2014;448:315–22.
Liu G, Xiang T, Wu QF, Wang WX. Long noncoding RNA H19-derived miR-675 enhances proliferation and invasion via RUNX1 in gastric Cancer cells. Oncol Res. 2016;23:99–107.
Slattery ML, Lundgreen A, Herrick JS, Caan BJ, Potter JD, Wolff RK. Associations between genetic variation in RUNX1, RUNX2, RUNX3, MAPK1 and eIF4E and riskof colon and rectal cancer: additional support for a TGF-beta-signaling pathway. Carcinogenesis. 2011;32:318–26.
Liang Z, Lu L, Mao J, Li X, Qian H, Xu W. Curcumin reversed chronic tobacco smoke exposure induced urocystic EMT and acquisition of cancer stem cells properties via Wnt/beta-catenin. Cell Death Dis. 2017;8:e3066.
Jamora C, DasGupta R, Kocieniewski P, Fuchs E. Links between signal transduction, transcription and adhesion in epithelial bud development. Nature. 2003;422:317–22.
Gonzalez DM, Medici D. Signaling mechanisms of the epithelial-mesenchymal transition. Sci Signal. 2014;7:re8.
Ding C, Luo J, Li L, Li S, Yang L, Pan H, Liu Q, Qin H, Chen C, Feng J. Gab2 facilitates epithelial-to-mesenchymal transition via the MEK/ERK/MMP signaling in colorectal cancer. J Exp Clin Cancer Res. 2016;35:5.
Blavier L, Lazaryev A, Shi XH, Dorey FJ, Shackleford GM, DeClerck YA. Stromelysin-1 (MMP-3) is a target and a regulator of Wnt1-induced epithelial-mesenchymal transition (EMT). Cancer Biol Ther. 2010;10:198–208.
Yang X, Du X, Sun L, Zhao X, Zhu J, Li G, Tian J, Li X, Wang Z. SULT2B1b promotes epithelial-mesenchymal transition through activation of the beta-catenin/MMP7 pathway in hepatocytes. Biochem Biophys Res Commun. 2019;510:495–500.
Yang Z, Li K, Liang Q, Zheng G, Zhang S, Lao X, Liang Y, Liao G. Elevated hydrostatic pressure promotes ameloblastoma cell invasion through upregulation of MMP-2 and MMP-9 expression via Wnt/beta-catenin signalling. J Oral Pathol Med. 2018;47:836–46.
Wan G, Liu Y, Zhu J, Guo L, Li C, Yang Y, Gu X, Deng LL, Lu C. SLFN5 suppresses cancer cell migration and invasion by inhibiting MT1-MMP expression via AKT/GSK-3beta/beta-catenin pathway. Cell Signal. 2019;59:1–12.
Chen EC, Karl TA, Kalisky T, Gupta SK, O'Brien CA, Longacre TA, van de Rijn M, Quake SR, Clarke MF, Rothenberg ME. KIT signaling promotes growth of Colon xenograft tumors in mice and is up-regulated in a subset of human Colon cancers. Gastroenterology. 2015;149:705–17 e702.
Chau WK, Ip CK, Mak AS, Lai HC, Wong AS. c-kit mediates chemoresistance and tumor-initiating capacity of ovarian cancer cells through activation of Wnt/beta-catenin-ATP-binding cassette G2 signaling. Oncogene. 2013;32:2767–81.
This work was supported by the Guangdong Gastrointestinal Disease Research Center (No.2017B020209003), the National Natural Science Funds of China (No.81772964), China’s Postdoctoral Science Fund (2018M643129), and the Foundation for the President of Nanfang Hospital of Southern Medical University (2018C027).
Ethics approval and consent to participate
All the patients were informed of sample collection and usage. The tissue samples were collected and used in accordance with approval by the Institutional Ethical Committee Board (Nanfang Hospital, Guangzhou, China). Use of animal was approved by the Nanfang hospital animal ethic committee.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. A. Relapse free survival with low/high RUNX1 gene expression was analyzed using a colorectal cancer data set of GSE14333/17536/31595/39582. B. Overall survival with low/high RUNX1 gene expression was analyzed using a colorectal cancer data set of GSE17536/39582. C. Expression of RUNX1 mRNA and protein detected by qPCR and western blot in colorectal cancer cell lines. (TIF 2047 kb)
Figure S2. A. Migration and invasion ability determined by Transwell assay in HCT116, RKO and SW480 cells with RUNX1 overexpressed or silencing. B. Migration ability detected by Wound healing in HCT116, RKO and SW480 cells with RUNX1 overexpressed or silencing. C. External whole-body、colon and ileocecus images of mice by orthotopic injection of HCT116/RUNX1 and HCT116/Vector cells. HE staining of orthotopic colorectal cancer were also exhibited. (TIF 11573 kb)
Figure S3. A. Positive correlation between RUNX1 expression and EMT targeted genes. B. Positive correlation between RUNX1 expression and CTNNB1、WNT5A、CD44 and LEF1 were observed. C. Relationship between RUNX1 expression and MMP3、MMP7、MMP9 and MMP14 were also found positively in GEPIA (GEPIA: http://gepia.cancer-pku.cn/). (TIF 3666 kb)
Figure S4. A. The mRNA expression level of MMP3、MMP7、MMP9 and MMP14 were higher in HCT116/RUNX1 group than that in HCT116/vector group. B. The protein level of AXIN1 detected by western blot in HCT116 and SW480 cells with RUNX1 overexpression or silencing. (TIF 431 kb)
Table S1. Primer Sequences Used for Real-time PCR (5′ to 3′). (DOC 31 kb)
Table S2. Chip primer sequenced used for quantitative real-time PCR. (DOC 29 kb)
Table S3. Relationship between RUNX1 mRNA expression and the clinicpathological features of CRC patients. (DOCX 14 kb)