Skip to main content


Springer Nature is making SARS-CoV-2 and COVID-19 research free. View research | View latest news | Sign up for updates

The prognostic significance of RUNX2 and miR-10a/10b and their inter-relationship in breast cancer



The major cancer related mortality is caused by metastasis and invasion. It is important to identify genes regulating metastasis and invasion in order to curtail metastatic spread of cancer cells.


This study investigated the association between RUNX2 and miR-10a/miR-10b and the risk of breast cancer relapse. Expression levels of RUNX2 and miR-10a/b in108 pairs of tumor and non-tumor tissue of breast cancer were assayed by quantitative PCR analysis and evaluated for their prognostic implications.


The median expression levels of RUNX2 and miR-10b in tumor tissue normalized using adjacent non-tumor tissue were significantly higher in relapsed patients than in relapse-free patients. Higher expression of these three genes were significantly correlated with the hazard ratio for breast cancer recurrence (RUNX2: 3.02, 95% CI = 1.50 ~ 6.07; miR-10a: 2.31, 95% CI = 1.00 ~ 5.32; miR-10b: 3.96, 95% CI = 1.21 ~ 12.98). The joint effect of higher expression of all three genes was associated with a hazard ratio of 12.37 (95% CI = 1.62 ~ 94.55) for relapse. In a breast cancer cell line, RUNX2 silencing reduced the expression of miR-10a/b and also impaired cell motility, while RUNX2 overexpression elicited opposite effects.


These findings indicate that higher expression of RUNX2 and miR-10a/b was associated with adverse outcome of breast cancer. Expression levels of RUNX2 and miR-10a/b individually or jointly are potential prognostic factors for predicting breast cancer recurrence. Data from in vitro studies support the notion that RUNX2 promoted cell motility by upregulating miR-10a/b.


Breast cancer is the most common malignant disease in women. Approximately 10 ~ 15% breast cancer patients have an aggressive disease and develop distant metastases [1] to bone, lung, liver and brain [2]. Up to 70% of bone metastases occur in advanced breast cancer patients [3]. Moreover, five year survivability for breast cancer patients with bone metastasis is only 20% after its detection [4].

Runt-related transcription factor 2 (RUNX2) has been shown to play an important role in osteogenesis and development of osteosarcoma [5]. Runx2 was also reported to be highly expressed in breast cancer with poor clinical outcomes [6]. Recently, McDonald et al. showed that high RUNX2 expression is significantly associated with estrogen receptor (ER)/progesterone receptor (PR)/HER2-negative breast cancers and that patients with high RUNX2 expression have a poorer survival rate than those with negative or low expression [7]. Moreover, in non-small cell lung cancer-patients, higher RUNX2 expression was significantly correlated with tumor progression and metastasis [8]. In epithelial ovarian cancer, various genes involved in tumor invasion and metastasis were suppressed upon RUNX2 knockdown [9]. Studies of breast cancer revealed regulation of several genes involved in bone invasion, such as MMPs, VEGF, OP, and BSP, by RUNX2, suggesting that this master transcription factor might contribute to bone metastasis in breast tumor [10]-[13]. This is consistent with the report that RUNX2 silencing reduced cell motility of metastatic breast cancer cell line, MDA-MB-231. On the other hand, RUNX2 overexpression increased cell migration ability in non-metastatic MCF7 breast cancer cell line [14].

MicroRNAs (miRNAs), a group of ~22 nucleotides endogenous and evolutionarily conserved single-stranded small non-coding RNAs, are crucial post-transcriptional regulators of a variety of biological processes, including the initiation, progression and metastases of cancer [15]-[18]. As reported in several studies, the miRNA-10 (miR-10) family, including miR-10a and miR-10b which are identical except for the 12th nucleotide [19], play an important role in tumorigenesis and progression [20],[21]. MiR-10a was reported to be downregulated in chronic myeloid leukaemia and acute myeloid leukaemia, and upregulated in colon cancer and hepatocellular carcinomas [20]. Mir-10a was also reported to gain in copy number in melanoma and breast cancer [22] and overexpression of miR-10a promoted cell migration and invasion of hepatoma cancer cell lines [23] and cervical cancer cell lines in vitro[24]. On the other hand, miR-10b was upregulated in pancreatic cancer and B-cell chronic lymphocytic leukemia [20] In addition, miR-10b was highly expressed in breast tumor with poor clinical outcomes [25] and facilitated cell migration and invasion in breast cancer [26].

These findings suggest that RUNX2 and miR-10a/b play important role in progression and metastases in breast cancer, but the association between RUNX2 and miR-10a/b, if any, is unknown. In this study, we try to decipher the relationship between RUNX2 and miR-10a/b in clinical breast cancer samples as well as in cell lines. We demonstrated that expression of RUNX2 significantly correlated with miR-10a/b in ER negative and triple negative breast cancers and the expression levels of RUNX2 and miR-10a/b individually or jointly were significant prognostic factors for predicting breast cancer recurrence. Furthermore, RUNX2 silencing in MDA-MB-231 cells downregulated miR-10a/b transcription and clearly impeded cell motility. These results indicated that RUNX2 plays an important role in regulating breast cancer progression.


Study patients and tissues

Ninety five of the 108 breast cancer patients examined in this study had clinicopathologically confirmed primary ductal carcinoma of the breast, and the remaining 13 patients had non-ductal carcinoma of the breast. All of them were diagnosed at the Tri-Service General Hospital, Taipei, Taiwan between October 1994 and February 2013. Patients’ clinical information, including cancer stage, tumor grade, estrogen receptor (ER), progesterone receptor (PR), HER-2/neu, recurrence and survival status, were also noted. Recurrent breast tumors were subjected to pathological confirmation to exclude the possibility of second primary tumors. Moreover, the cause of death was verified from death certificate; patients whose death was clearly documented to be due to breast cancer were considered to have died of breast cancer, whereas other causes of deaths were considered censored events. Tumor and adjacent non-tumor breast tissues from patients were obtained at the time of initial surgery, and were fully encoded to protect patient confidentiality. Clinical specimens were utilized under a protocol approved by the Institutional Review Board of the Human Subjects Research Ethics Committee of Academia Sinica and Tri-Service General Hospital, Taipei, Taiwan.

Statistical analysis

Histogram data were given as mean ± standard deviation (SD). Comparison between the tumor and non-tumor breast tissues was performed using Student’s t-test. We examined whether there was any correlation between the expression levels of the miR-10a/b and RUNX2 genes, individually or jointly, in primary cancer tissue and the clinicopathological features of the tumor. The goal was to evaluate the usefulness of these three genes as prognostic biomarkers for breast cancer. The expression levels of the three genes in tumor or non-tumor tissues of different clinicopathological feathers were described by mean ± SD and median. The results of statistical analyses were examined against those of Student’s t-test and Wilcoxon rank-sum test. The estimated areas under the curves (AUC) for assessing the accuracy of prediction were calculated based on the receiver operating characteristic (ROC) curves, where an optimal diagnostic cutoff point was determined to differentiate over expression of individual genes. Based on these cutoffs, quantitative measurements of the expression levels of individual genes were converted into binary measurements to help examine whether over expression of these genes could be prospectively associated with breast cancer progression using the relapse-free survival (RFS) and overall survival (OS) as outcomes of interest. RFS was measured as the time from surgery to recurrence or the end of the study, and OS was defined as the time from surgery to death or the end of the last follow-up for a patient. Survival curves were plotted by using the Kaplan–Meier method and 5-year RFS-rate and 10-year OS-rate point estimates with 95% confidence interval (CI) were reported. The significance of these associations was assessed using the log-rank test. The Cox regression model was used to compute hazard ratios. P-value less than 0.05 were considered statistically significant. All data analyses were performed using the SAS statistical software for Windows, version 9.3 (SAS institute, Cary, NC, USA).

Cell cultures

MCF7 breast cancer cell line was cultured in Modified Eagle Medium (MEM) supplemented with 10% fetal bovine serum (FBS). MDA-MB-231 and MB157 breast cancer cell lines, U2OS osteosarcoma cell line and 293 T embryonic kidney were cultured in Dulbecco’s Modified Eagle Medium (DMEM) supplemented with 10% FBS. BT483 breast cancer cell line was cultured in Roswell Park Memorial Institute medium (RPMI) supplemented with 10% FBS. All cell lines were obtained from American Type Culture Collection (ATCC, Manassas, VA, USA).

RUNX2 plasmid, RUNX2 siRNA and miRNA Oligo synthesis and transfection

Total RNA from the osteosarcoma cell line U2OS was extracted using Trizol reagent (Invitrogen, Carlsbad, CA, USA), and 1 μg of RNA was reverse transcribed using the SuperScript III RT–PCR kit (Invitrogen). For cDNA cloning, full-length RUNX2wt was amplified from cDNA of the U2OS cell line using the PCR primer pair (RUNX2wt-forward: 5′-GCTAGCATGGACTACAAAGACGATGACGACAAGGTGATGCGTATTCCCGTAGATCCGAG-3′) and (RUNX2wt-reverse:5′-CTCGAGATATGGTCGCCAAACAGATTCATCC-3′) in an amplification mixture containing 1X Pfu polymerase buffer (Stratagene, La Jolla, CA, USA), 100 μM 4dNTPs, 1.5 mM MgCl2 and 1U Pfu polymerase (Stratagene). The samples were amplified for 35 cycles of 95°C for 1 min 30 sec, 60°C for 35 sec, and 72°C for 2 min each, preceded by a 5 min denaturation at 95°C and followed by a 10 min extension at 72°C. The amplified cDNAs of RUNX2wt were then cloned into pcDNA3.1 (+) vector (Invitrogen) by ligation of the NheIXhoI fragments to generate the pcDNA-RUNX2 expression constructs. The sequence of pcDNA-RUNX2 was confirmed by DNA sequencing.

RUNX2 siRNA (Cat. # 1027020) and siRNA control (Cat. # 1027280) were guaranteed products by manufacturers and purchased from Qiagen (Valencia, CA, USA). Precursor miRNA oligos (pre-miR control # 1: Cat. # AM17110, pre-has-miR-10a: Cat. # PM10787 and pre-has-miR-10b: Cat. # PM11108) and anti-miRNA oligos (anti-miR control # 1: Cat. # AM17010, anti-has-miR-10a: Cat. # AM10787 and anti-has-miR-10b-5p: Cat. # AM11108) were guaranteed products by manufacturers and purchased from Ambion (Austin, TX, USA).

For plasmid transfection, 4 g of pcDNA-RUNX2 construct plasmid per well with 4 μl Lipofectamine 2000 (Invitrogen) was transfected into cell lines. For siRNA transfection, 25 nM RUNX2 siRNA per well with 2 μl Lipofectamine RNAiMax (Invitrogen) was transfected into cell lines. For miRNA oligos transfection, 25 nM miRNA oligos per well with 2 μl RNAiMax were transfected into cell lines. For anti-miRNA oligo transfection, 15 nM anti-miRNA oligos per well with 2 μl RNAiMax were transfected into cell lines. The experiments were repeated three times using six-well plates.

Real-time PCR quantification of miRNA-10a/miR-10b and RUNX2

Total RNAs were extracted from tumor and non-tumor breast tissues for each patient or breast cancer cell lines using Trizol reagent according to the manufacturer’s specifications. 4 ng of total RNAs were used for quantification of miRNAs expression, including miRNA-10a/miR-10b and RNU6B (U6) RNA, by TaqMan RT-qPCR kit (miRNA-10a: Cat. # 000387, miRNA-10b: Cat. # 002218, and RNU6B: Cat. # 001093; Applied Biosystems, Foster City, CA, USA) according to the manufacturer’s instructions. For the quantification of mRNA of RUNX2 and GAPDH, 1 g of total RNAs was reverse transcribed into cDNA using the SuperScript III kit and the specific mRNAs were detected by Power SYBR Green PCR Master Mix (Applied Biosystems) according to the manufacturer’s instructions. The primer sequences of RUNX2 and GAPDH qPCR were: RUNX2 5′-AGCAAGGTTCAACGATCTGAGATT-3′ (forward) and 5′-AAGACGGTTATGGTCAAGGTGAAA-3′ (reverse); GAPDH 5′-CTGCTCCTCCTGTT CGACAGT-3′ (forward) and 5′-ACCTTCCCCATGGTGTCTGA-3′ (reverse). The gene expression was normalized based on threshold cycle (Ct) values of common internal control for miRNA quantification assays, U6 RNA, and RUNX2 quantification assays, GAPDH. The level of RUNX2 and miR-10a/b expression were measured using the 2-Ct method [27]. The results were presented as fold change of each gene in the tumor tissue relative to the non-tumor tissue. The RT-qPCR was performed on the 7300 Real-Time PCR System (Applied Biosystems). The procedures were repeated three times for all samples, and mean values were obtained for statistical analysis.

Western blotting

Cell lysates were extracted in 0.1 M NaCl, 0.01 M TriszHCl (pH 7.6), 1 MM EDTA (pH 8.0), 1 Mg/ml aprotinin, and 100 mg/ml PMSF, while protein concentration was determined by a Bio-Rad protein assay. Protein (50 μg) was boiled at 95°C in sodium dodecyl sulfate (SDS) sample buffer for 5 min, electrophoresed on 7–12% SDS-PAGE gels, and transferred to polyvinyldifluoridine membranes. These were then incubated overnight at 4°C with (i) rabbit monoclonal antibody against human RUNX2 (Cell Signaling Technology, Danvers, MA, USA; #8486; 1:250), (ii) goat polyclonal antibody against human HOXA1 (R&D Systems Inc., Minneapolis, MN, USA; AF5014; 1:500), (iii) goat polyclonal antibody against human HOXD10 (Santa Cruz Biotechnology, Dallas, TX, USA; sc-33005; 1:500), or (vi) rabbit polyclonal antibody against human Actin (Santa Cruz Biotechnology; sc-1615-R; 1:3000). Membranes were washed three times for 15 min with phosphate buffered saline (PBS) containing 0.1% Tween20 (PBST), incubated with anti-goat or anti-rabbit secondary antibody (both Santa Cruz Biotechnology; 1:5000) at room temperature for 60 min, and then washed three times for 15 min with PBST. The signals of protein bands were visualized using ECF western blotting kit (Amersham Biosciences, Piscataway, NJ, USA) and detected by Typhoon 9400 imager (Amersham Biosciences).

Cell migration and invasion assays

Cell migration and invasion were assayed using a transwell chamber (Millipore, Temecula, CA, USA) with or without Matrigel (BD, Franklin Lakes, NJ, USA). 1 Ml serum free medium was added into a 24-well plate, and a transwell chamber was placed into the well in both assays. The chamber was coated with 50 μl Matrigel and incubated for 2 hours at 37°C for the invasion assay. After transfection for 72 hours, in migration assay: cells were trypsinized and seeded into chambers at the density of 1 × 104 cells for MDA-MB-231 and 3 × 104 cells for MCF7 per well; and in invasion assay: cells were 2 × 104 cells for MDA-MB-231 and 6 × 104 cells for MCF7. Eight hours later for migration and 24 hours later for invasion in MDA-MB-231, and 12 hours later for migration and 32 hours later for invasion in MCF7, migrated cells remaining in the basal chamber were fixed and stained with 10% ethanol in 1% crystal violet for 30 min and washed three times with PBS. Non-migrated cells were removed by cotton swabs. Images of the stained cells were obtained under a phase-contrast microscope (Olympus, Tokyo, Japan). To quantify the number of invading cells, crystal violet-stained cells in visual fields were counted. The experiments were repeated three times and mean values were obtained for statistical analysis. Detailed materials and methods can be found in Additional file 1.


Clinical features of breast cancer patients

It has been shown that RUNX2 plays an important role in cell migration and invasion in breast cancer and higher expression is associated with poor outcome [6],[8],[14]. Meanwhile, others have reported the contribution of miR-10a/b to cell migration and invasion in breast cancer and miR-10b overexpression correlates with poor prognosis [23]-[26]. However, the inter-relationship between RUNX2 and miR-10a/b and their joint impact on prognosis in breast cancer has yet to be investigated. We examined one hundred and eight patients with histologically proven breast cancers to delineate the relationship between the expression of RUNX2 and miR-10a/b, and clinical outcomes. Their clinical and pathological features were summarized in Table 1. The mean age was 59.2 × 13.2 years and the median follow up for relapsed patients and overall survival was 54.6 and 63.2 months, respectively. The distribution of cancer stages was stage 0 (2.80%), stage I (31.78%), stage II (42.99%) and stage III (22.43%), respectively, and of the molecular subtypes was Luminal A (ER + or PR+, Her2-) (33.33%), Luminal B (ER + or PR+, Her2+) (35.19%), Her2 positive (ER-, PR-, Her2+) (15.74%) and triple negative (ER-, PR-, Her2-) (15.74%), respectively. These clinical and pathological characteristics of the breast cancer subjects across the strata of age, hormone receptor status, and Her2, were similar to those reported in other breast cancer clinics in Taiwan [28]-[30].

Table 1 Clinical and pathological characteristics of breast cancer patients

The relationship between expression of RUNX2 and miR-10a/b genes in breast cancer tumors and clinical parameters

The expression levels of RUNX2 and miR-10a/b were determined in tumor tissues and adjacent non-tumor tissues by qPCR. As shown in Table 2, the median expression levels of RUNX2 and miR-10b were significantly higher in tumor tissue than in non-tumor tissues by 2.16 and 1.73 fold, respectively (p < 0.0001, p < 0.0001). There is also a trend for higher expression of miR-10a in tumor tissues than in non-tumor part, but the difference did not reach statistical significance (p = 0.1318). Patients with ER- or triple negative subtypes displayed higher-level expression of these three genes than those with ER + or non-triple negative subtypes. Analysis of clinical outcome revealed that patients who relapsed or died of breast cancer showed higher expression of these three genes than non-relapse or alive patients. On the other hand, no significant association was found between the expression levels of these three genes and gender, histology, stage, regional lymph node involvement or status of PR and Her2, although there was a trend for higher expression in patients with higher histology grade (Table 2; Additional file 1: Table S1).

Table 2 Correlation of expression of miR-10a, miR-10b and RUNX2 with clinical parameters of breast cancer patients

Overexpression of RUNX2 and miR-10a/b genes in tumors as predictors of breast cancer outcome

To evaluate the prognostic significance of expression levels of RUNX2 and miR-10a/b, we used the receiver operating characteristic (ROC) curves based on the expression levels of RUNX2 and miR-10a/b in individual tumor/non-tumor ratios to determine the cutoff values defining “overexpression” of the three genes. According to these cutoff values, the areas under the curves (AUC) values for individual or various combinations of these three genes were found to be significantly lower in non-relapsed patients than in relapsed patients. However, when comparing patients with fatal outcome to those alive, only AUC values for miR-10b and combination of these three genes were significantly lower in alive patients than in those who died (Additional file 1: Table S2). Accordingly, we used the cutoff values to define the status of “overexpression” in subsequent logistic and COX regression analyses to ascertain correlations between gene expression and prognostic features. Under these cutoff values, the frequencies of overexpression of RUNX2 and miR-10a/b genes were significantly higher in patients who relapsed, and that of miR-10b overexpression was higher in patients who died (all p < 0.05). Using the same dichotomized values for logistic regression analysis, the odds ratios (ORs) for relapse were higher in patients with increased mRNA or miRNA levels of these 3 individual genes, while the OR for death was only associated with miR-10b overexpression (all p < 0.05). Patients with no overexpression of these three genes served as the reference group. In terms of the joint effects of overexpression of all three genes on prognosis, the ORs were respectively 20.32 (95% confidence interval (CI) = 1.50 ~ 179.92) and 17.05 (95% CI = 1.76 ~ 165.01) for relapse and death (P < 0.05). There was a trend toward an incremental additive effect on significantly elevated risk of relapse (OR = 2.66; 95% CI = 1.50 ~ 4.75), death (OR = 2.90; 95% CI = 1.22 ~ 6.89) and regional lymph node involvement (OR = 1.51: 95% CI = 1.03 ~ 2.21), respectively, as patients carried a greater number of overexpressed genes (P-trend < 0.05; Table 3; Additional file 1: Table S3).

Table 3 The relationship between overexpression of individual gene or joint effects of miR-10a, miR-10b, and RUNX2 with prognosis in breast cancer patients

The correlations between overexpression of RUNX2 and miR-10a/b and poor clinical outcome

The correlations between overexpression of RUNX2 and miR-10a/b, individually or jointly, and clinical outcomes were examined by the Kaplan–Meier method of patients. Both RFS (Figure 1) and OS (Additional file 1: Figure S1) were significantly lower in patients with overexpression of RUNX2 or miR-10b individually (Figure 1A, C; Additional file 1: Figure S1 A, C). Compared to the patients with no overexpression of these three genes, those with a greater number of overexpressed genes (1, 2, or 3) in tumors showed lower rates of RFS (log-rank p < 0.0001) and OS (log-rank p = 0.0002) (Figure 1D and Additional file 1: Figure S1D). Multivariate Cox regression analysis showed a significant correlation between higher hazard ratios (HRs) for relapse and greater gene signals for all three genes individually and a correlation between higher HR for death and expression of miR-10b (Table 4; Additional file 1: Table S4). As to the joint effects of overexpression of all three genes compared with no overexpression of any of them on prognosis, the HRs were 12.37 (95% CI = 1.62 ~ 94.55) and 14.59 (95% CI = 1.66 ~ 128.24) for relapse and death, respectively (P = 0.015, 0.016). The multiple gene additive model displayed a statistically significant HR of 2.25 for RFS (95% CI = 1.44 ~ 3.52, P-trend = 0.0004) and 2.88 for OS (95% CI = 1.26 ~ 6.61, P-trend = 0.0125).

Figure 1

Prognostic significance of increased expression levels of RUNX2, miR-10a or miR-10b in breast cancer patients. Kaplan-Meier statistical analyses were conducted to examine the association between relapse-free survival and the expression of RUNX2 (A), miR-10a (B) or miR-10b (C), and the number of overexpressed genes (D) in all patients.

Table 4 Cox regression model analyses of the gene expression levels in association with RFS of breast cancer patients

Suppression of RUNX2 decreased miR-10a/b expression and impaired cell migration and invasion ability

In view of our findings of significant correlations among RUNX2 or miR-10a/b and relapse or death, we next investigated the inter-relationship of these 3 genes in breast cancer cell lines in vitro. Using qPCR analysis, we determined the expression levels of RUNX2 and miR-10a/b in 4 breast cancer cell lines. As shown in Figure 2A, MDA-MB-231, the invasive estrogen receptor negative (ER-) cell line, showed the highest expression levels of all these three molecules. Since MDA-MB-231 cells expressed very little Twist [31],[32], which is reported to upregulate miR-10b [26] and since RUNX2 is involved in transcriptional activities of Twist [33], we explored the possible regulation of miR-10a/b by RUNX2. Silencing of RUNX2 by siRNA transfection of MDA-MB-231 cell line to 24% (p < 0.0001) of control resulted in a significant decrease in expression of miR-10a and miR-10b to 52% (p = 0.0349) and 51% (p < 0.0001), respectively (Figure 2B). Consistent with this, HOXA1 and HOXD10 are putative targets for miR-10a/b in Garzon et al. 2006 [34], suppression or overexpression of miR-10a/b increased or decreased HOXA1/HOXD10 protein expression, respectively, in breast cancer cell lines explained in Additional file 1: Figure S2, escalated in the silenced RUNX2 cells by 3.86 and 2.33 fold, respectively (Figure 2C). In line with the increased HOXA1 and HOXD10, migration and invasion ability of MDA-MB-231 cells were significantly reduced to 21% (p = 0.0035) and 29% (p = 0.0029), respectively. (Figure 2D, E). Thus, suppression of RUNX2 reduced miR-10a/b expression, and thereby impairing breast cancer cell motility.

Figure 2

RUNX2 silencing reduced the expression of miR-10a and miR-10b and inhibited cell migration and invasion. (A) Expression level of miR-10a and miR-10b in breast cancer cell lines were analyzed by qPCR analysis and normalized to U6 snRNA expression. (B) MDA-MB-231 cells were transfected with 25 μM control siRNA or RUNX2 siRNA for 72 hours. RUNX2, miR-10a or miR-10b expression levels were determined by qPCR analysis. (C) Expression levels of RUNX2, HOXA1 and HOXD10 proteins were assessed by Western blot analysis and normalized to actin. (D) Cell migration and invasion ability were assessed by Transwell assay or Matrigel assay, respectively as described in Material and Methods. Black line: 20 μl, Red line: 100 μl (E) The quantification result was based on three independent experiments. *: p < 0.05, t-test.

Overexpression of RUNX2 induced transcription of miR-10a/b and promoted cell migration and invasion

In order to clarify the relationship between RUNX2 and miR-10a/b, RUNX2 was overexpressed in high and low RUNX2 expressing breast cancer cell lines, MDA-MB-231 and MCF7 (Figure 3A, B). RUNX2 increased the transcription of miR-10a by 1.61 fold (p = 0.0079) and that of miR-10b by 1.60 fold (p < 0.0001) in MCF7 cells and increased only the expression of miR-10b by 1.52 fold (p = 0.0413) in MDA-MB-231 cells. Consistently, protein levels of HOXA1 and HOXD10 markedly declined to 0.15 and 0.09 fold, respectively, in RUNX2 overexpressing MCF7 cells, but only slightly decreased to 0.91 and 0.71 fold in MDA-MB-231 cells (Figure 3C). Functional assay showed that cell migration and invasion were significantly enhanced in RUNX2 overexpressing MDA-MB-231 (1.74 fold, p = 0.0315 and 2.51 fold, p = 0.0016, for migration and invasion, respectively) and MCF7 cells (2.96 fold, p = 0.0143 and 4.95 fold, p = 0.0018, respectively) (Figure 3D, E, F, G). Thus, overexpression of RUNX2 upregulated miR-10a/b expression, thereby promoting breast cancer cell motility.

Figure 3

Overexpression of RUNX2 induced expression of miR-10a and miR-10b and promoted cell migration and invasion. (A, B) MDA-MB-231 or MCF7 cells were transfected with 4 μg control pcDNA3.1 vector or pcDNA-RUNX2 plasmid for 72 hours. RUNX2, miR-10a or miR-10b expression levels were analyzed by qPCR. (C) RUNX2, HOXA1 and HOXD10 protein levels were assessed by Western blot analysis and normalized to actin. (D, E) Cell migration and invasion ability were assessed by Transwell assay or Matrigel assay, respectively. Black line: 20 μl, Red line: 100 μl (F, G). The quantification result was based on three independent experiments. *: p < 0.05, t-test.

To sum up, these findings suggest that RUNX2 plays an important role in controlling breast cancer cell motility through regulation of miR-10a/b expression, and overexpression of these three genes has poor prognosis.


Metastasis is the major cause of mortality in cancer and affects the clinical progression of cancer patients [35]. In this study, primary breast tumor tissues were found to have higher expression levels of RUNX2 and miR-10a/b than adjacent non-tumor tissue, and overexpression in any of these three genes was significantly associated with poor clinical outcome.

Recently, Ma et al. has shown that the transcription factor Twist regulates the expression of miR-10b, which inhibits HOXD10 protein synthesis, thereby permitting the production of RHOC proteins and activating migration and invasion of breast cancer cells. Analyses of clinical specimens also indicate that miR-10b is highly expressed in not only metastatic breast cancer [26], but also early stage breast cancer [25]. However, as suggested in another study, in early-stage breast cancer patients, the expression of miR-10b appears to show no significant association with the development of distant metastases, recurrence-free survival, distant-relapse-free survival, or breast-cancer-specific survival [36]. In this study, we showed that RUNX2 upregulated miR-10a/b and promoted breast cancer cell migration and invasion. Furthermore, we demonstrated that miR-10a played a role in cell motility, similar to miR-10b [26], thus providing an additional mechanism in regulating breast cancer metastasis and invasion.

We found that, the expression of RUNX2 and miR-10a/b in ER- and triple negative breast cancer was higher than in ER + breast cancer. This observation may be related to the intricate interactions among estrogen, estrogen receptor alpha (ERα) and RUNX2, which is of crucial importance in osteoporosis and breast cancer. ERα can form a complex with RUNX2 to affect the capability of RUNX2 in transcribing downstream genes [37]. RUNX2 regulates the expression of ERα by interacting with its F promoter, one of the multiple promoters of the human ERα gene [33]. Also to be noted is that estrogen may also regulate the expression and activation of RUNX2 through various signal transduction pathways, notably the estrogen-transforming growth factor-beta (TGFβ) pathway [38], the Wnt pathway [39], the Fas/Fas Ligand pathway [40], and the nuclear factor-kappaB (NFκB) pathway [41]. Estrogen, ER and RUNX2 work to mutually regulate their expression and activation. This may help to explain the lower expression of RUNX2 in ER + breast cancer than ER- tumors.

We illustrated the regulation of the expression of miR-10a/b by RUNX2 in breast cancer cells. Several putative RUNX2 transcription factor binding sites were found at the proximal sequences preceding pre-miR-10a/b, based on website software PROMO analysis ( [42]-[44]. Several different reporter constructs that carried putative RUNX2 binding sites within 15 k/24 k upstream of pre-miR-10a/b, respectively, were generated (Additional file 1: Table S5). We then compared the ability of RUNX2 to promote the transcription of reporter constructs carrying the RUNX2 sites. Unfortunately, RUNX2 failed to enhance the transcription of the miR-10a/b promoter regions that we cloned (Additional file 1: Figure S3). Studies to delineate the transcriptional regulation of miR-10a/b by RUNX2 will be further pursued in the future. Moreover, the RUNX2 cofactors, including STAT1, Twist, YAP or ERα protein [37],[45], have been reported to play important roles in RUNX2 transcription ability and regulation of downstream gene expression. These RUNX2 cofactors will be incorporated into future studies to further clarify how RUNX2 regulates miR-10a/b transcription.


Through in vitro experiments, our study demonstrates the impact of RUNX2 on breast cancer migration and invasion and its regulation of miR-10a/b expression. Importantly, our clinical-laboratory correlative analyses indicate that breast cancer patients with higher expression of RUNX2 and miR-10a/b – either individually or jointly – tend to suffer a greater risk of recurrence or death. Therefore, the three genes – RUNX2, miR-10a, and miR-10b – are deemed valuable markers for prognostication of breast cancer patients.

Authors’ contributions

CHC, TCF, YCL, ACS, WHL and ALY carried out the conception and experimental design; JCY and GSL carried out the acquisition of clinical sample and clinical information; CHC and TCF carried out the performed the experiments; CHC, ACS, WHL and ALY carried out the analysis and interpretation of data; CHC, TCF, ACS, WHL and ALY carried out the writing, review and/or revision of the manuscript. All authors read and approved the final manuscript.

Additional file



Runt-related transcription factor 2


miRNA, microRNA


Estrogen receptor


Progesterone receptor


Standard deviation


Receiver operating characteristic


Areas under the curves


Relapse-free survival


Overall survival


Confidence interval


Modified eagle medium


Fetal bovine serum


Dulbecco’s modified eagle medium


Roswell Park Memorial Institute medium


American Type Culture Collection




Threshold cycle


Sodium dodecyl sulfate


Phosphate buffered saline


PBS containing 0.1% Tween20


Odds ratio


Hazard ratio


Estrogen receptor alpha


Transforming growth factor-beta


Nuclear factor-kappaB


  1. 1.

    Weigelt B, Peterse JL, van `t Veer LJ: Breast cancer metastasis: markers and models. Nat Rev Cancer. 2005, 5: 591-602. 10.1038/nrc1670.

  2. 2.

    Lee YT: Breast carcinoma: pattern of metastasis at autopsy. J Surg Oncol. 1983, 23: 175-180. 10.1002/jso.2930230311.

  3. 3.

    Coleman RE, Rubens RD: The clinical course of bone metastases from breast cancer. Br J Cancer. 1987, 55: 61-66. 10.1038/bjc.1987.13.

  4. 4.

    Coleman RE: Metastatic bone disease: clinical features, pathophysiology and treatment strategies. Cancer Treat Rev. 2001, 27: 165-176. 10.1053/ctrv.2000.0210.

  5. 5.

    Ducy P, Zhang R, Geoffroy V, Ridall AL, Karsenty G: Osf2/Cbfa1: a transcriptional activator of osteoblast differentiation. Cell. 1997, 89: 747-754. 10.1016/S0092-8674(00)80257-3.

  6. 6.

    Onodera Y, Miki Y, Suzuki T, Takagi K, Akahira J, Sakyu T, Watanabe M, Inoue S, Ishida T, Ohuchi N, Sasano H: Runx2 in human breast carcinoma: its potential roles in cancer progression. Cancer Sci. 2010, 101: 2670-2675. 10.1111/j.1349-7006.2010.01742.x.

  7. 7.

    McDonald L, Ferrari N, Terry A, Bell M, Mohammed ZM, Orange C, Jenkins A, Muller WJ, Gusterson BA, Neil JC, Edwards J, Morris JS, Cameron ER, Blyth K: RUNX2 correlates with subtype-specific breast cancer in a human tissue microarray, and ectopic expression of Runx2 perturbs differentiation in the mouse mammary gland. Dis Model Mech. 2014, 7: 525-534. 10.1242/dmm.015040.

  8. 8.

    Li H, Zhou RJ, Zhang GQ, Xu JP: Clinical significance of RUNX2 expression in patients with nonsmall cell lung cancer: a 5-year follow-up study. Tumour Biol. 2013, 34: 1807-1812. 10.1007/s13277-013-0720-4.

  9. 9.

    Li W, Xu S, Lin S, Zhao W: Overexpression of runt-related transcription factor-2 is associated with advanced tumor progression and poor prognosis in epithelial ovarian cancer. J Biomed Biotechnol. 2012, 2012: 456534-

  10. 10.

    Pratap J, Lian JB, Javed A, Barnes GL, van Wijnen AJ, Stein JL, Stein GS: Regulatory roles of Runx2 in metastatic tumor and cancer cell interactions with bone. Cancer Metastasis Rev. 2006, 25: 589-600. 10.1007/s10555-006-9032-0.

  11. 11.

    Pratap J, Wixted JJ, Gaur T, Zaidi SK, Dobson J, Gokul KD, Hussain S, van Wijnen AJ, Stein JL, Stein GS, Lian JB: Runx2 transcriptional activation of Indian Hedgehog and a downstream bone metastatic pathway in breast cancer cells. Cancer Res. 2008, 68: 7795-7802. 10.1158/0008-5472.CAN-08-1078.

  12. 12.

    Chen YC, Sosnoski DM, Mastro AM: Breast cancer metastasis to the bone: mechanisms of bone loss. Breast Cancer Res. 2010, 12: 215-10.1186/bcr2781.

  13. 13.

    Pratap J, Javed A, Languino LR, van Wijnen AJ, Stein JL, Stein GS, Lian JB: The Runx2 osteogenic transcription factor regulates matrix metalloproteinase 9 in bone metastatic cancer cells and controls cell invasion. Mol Cell Biol. 2005, 25: 8581-8591. 10.1128/MCB.25.19.8581-8591.2005.

  14. 14.

    Leong DT, Lim J, Goh X, Pratap J, Pereira BP, Kwok HS, Nathan SS, Dobson JR, Lian JB, Ito Y, Voorhoeve PM, Stein GS, Salto-Tellez M, Cool SM, van Wijnen AJ: Cancer-related ectopic expression of the bone-related transcription factor RUNX2 in non-osseous metastatic tumor cells is linked to cell proliferation and motility. Breast Cancer Res. 2010, 12: R89-10.1186/bcr2762.

  15. 15.

    Lee RC, Feinbaum RL, Ambros V: The C elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell. 1993, 75: 843-854. 10.1016/0092-8674(93)90529-Y.

  16. 16.

    Lee RC, Ambros V: An extensive class of small RNAs in Caenorhabditis elegans. Science. 2001, 294: 862-864. 10.1126/science.1065329.

  17. 17.

    Browne G, Taipaleenmaki H, Stein GS, Stein JL, Lian JB: MicroRNAs in the control of metastatic bone disease. Trends Endocrinol Metab. 2014, 25: 320-327. 10.1016/j.tem.2014.03.014.

  18. 18.

    Lian JB, Stein GS, van Wijnen AJ, Stein JL, Hassan MQ, Gaur T, Zhang Y: MicroRNA control of bone formation and homeostasis. Nat Rev Endocrinol. 2012, 8: 212-227. 10.1038/nrendo.2011.234.

  19. 19.

    Griffiths-Jones S, Grocock RJ, van Dongen S, Bateman A, Enright AJ: miRBase: microRNA sequences, targets and gene nomenclature. Nucleic Acids Res. 2006, 34: D140-D144. 10.1093/nar/gkj112.

  20. 20.

    Lund AH: miR-10 in development and cancer. Cell Death Differ. 2010, 17: 209-214. 10.1038/cdd.2009.58.

  21. 21.

    Tehler D, Hoyland-Kroghsbo NM, Lund AH: The miR-10 microRNA precursor family. RNA Biol. 2011, 8: 728-734. 10.4161/rna.8.5.16324.

  22. 22.

    Zhang L, Huang J, Yang N, Greshock J, Megraw MS, Giannakakis A, Liang S, Naylor TL, Barchetti A, Ward MR, Yao G, Medina A, O’Brien-Jenkins A, Katsaros D, Hatzigeorgiou A, Gimotty PA, Weber BL, Coukos G: microRNAs exhibit high frequency genomic alterations in human cancer. Proc Natl Acad Sci U S A. 2006, 103: 9136-9141. 10.1073/pnas.0508889103.

  23. 23.

    Yan Y, Luo YC, Wan HY, Wang J, Zhang PP, Liu M, Li X, Li S, Tang H: MicroRNA-10a is involved in the metastatic process by regulating Eph tyrosine kinase receptor A4-mediated epithelial-mesenchymal transition and adhesion in hepatoma cells. Hepatology. 2013, 57: 667-677. 10.1002/hep.26071.

  24. 24.

    Long MJ, Wu FX, Li P, Liu M, Li X, Tang H: MicroRNA-10a targets CHL1 and promotes cell growth, migration and invasion in human cervical cancer cells. Cancer Lett. 2012, 324: 186-196. 10.1016/j.canlet.2012.05.022.

  25. 25.

    Markou A, Yousef GM, Stathopoulos E, Georgoulias V, Lianidou E: Prognostic significance of metastasis-related microRNAs in early breast cancer patients with a long follow-up. Clin Chem. 2014, 60: 197-205. 10.1373/clinchem.2013.210542.

  26. 26.

    Ma L, Teruya-Feldstein J, Weinberg RA: Tumour invasion and metastasis initiated by microRNA-10b in breast cancer. Nature. 2007, 449: 682-688. 10.1038/nature06174.

  27. 27.

    Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(−delta delta C(T)) method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.

  28. 28.

    Yu JC, Ding SL, Chang CH, Kuo SH, Chen ST, Hsu GC, Hsu HM, Hou MF, Jung LY, Cheng CW, Wu PE, Shen CY: Genetic susceptibility to the development and progression of breast cancer associated with polymorphism of cell cycle and ubiquitin ligase genes. Carcinogenesis. 2009, 30: 1562-1570. 10.1093/carcin/bgp173.

  29. 29.

    Ding SL, Yu JC, Chen ST, Hsu GC, Hsu HM, Ho JY, Lin YH, Chang CC, Fann CS, Cheng CW, Wu PE, Shen CY: Diverse associations between ESR1 polymorphism and breast cancer development and progression. Clin Cancer Res. 2010, 16: 3473-3484. 10.1158/1078-0432.CCR-09-3092.

  30. 30.

    Lin CH, Liau JY, Lu YS, Huang CS, Lee WC, Kuo KT, Shen YC, Kuo SH, Lan C, Liu JM, Kuo WH, Chang KJ, Cheng AL: Molecular subtypes of breast cancer emerging in young women in Taiwan: evidence for more than just westernization as a reason for the disease in Asia. Cancer Epidemiol Biomarkers Prev. 2009, 18: 1807-1814. 10.1158/1055-9965.EPI-09-0096.

  31. 31.

    Yang J, Mani SA, Donaher JL, Ramaswamy S, Itzykson RA, Come C, Savagner P, Gitelman I, Richardson A, Weinberg RA: Twist, a master regulator of morphogenesis, plays an essential role in tumor metastasis. Cell. 2004, 117: 927-939. 10.1016/j.cell.2004.06.006.

  32. 32.

    Horikawa T, Yang J, Kondo S, Yoshizaki T, Joab I, Furukawa M, Pagano JS: Twist and epithelial-mesenchymal transition are induced by the EBV oncoprotein latent membrane protein 1 and are associated with metastatic nasopharyngeal carcinoma. Cancer Res. 2007, 67: 1970-1978. 10.1158/0008-5472.CAN-06-3933.

  33. 33.

    Lambertini E, Penolazzi L, Tavanti E, Schincaglia GP, Zennaro M, Gambari R, Piva R: Human estrogen receptor alpha gene is a target of Runx2 transcription factor in osteoblasts. Exp Cell Res. 2007, 313: 1548-1560. 10.1016/j.yexcr.2007.02.002.

  34. 34.

    Garzon R, Pichiorri F, Palumbo T, Iuliano R, Cimmino A, Aqeilan R, Volinia S, Bhatt D, Alder H, Marcucci G, Calin GA, Liu CG, Bloomfield CD, Andreeff M, Croce CM: MicroRNA fingerprints during human megakaryocytopoiesis. Proc Natl Acad Sci U S A. 2006, 103: 5078-5083. 10.1073/pnas.0600587103.

  35. 35.

    Chiang AC, Massague J: Molecular basis of metastasis. N Engl J Med. 2008, 359: 2814-2823. 10.1056/NEJMra0805239.

  36. 36.

    Gee HE, Camps C, Buffa FM, Colella S, Sheldon H, Gleadle JM, Ragoussis J, Harris AL: MicroRNA-10b and breast cancer metastasis. Nature. 2008, 455: E8-E9. 10.1038/nature07362. author reply E9

  37. 37.

    Khalid O, Baniwal SK, Purcell DJ, Leclerc N, Gabet Y, Stallcup MR, Coetzee GA, Frenkel B: Modulation of Runx2 activity by estrogen receptor-alpha: implications for osteoporosis and breast cancer. Endocrinology. 2008, 149: 5984-5995. 10.1210/en.2008-0680.

  38. 38.

    Hawse JR, Subramaniam M, Ingle JN, Oursler MJ, Rajamannan NM, Spelsberg TC: Estrogen-TGFbeta cross-talk in bone and other cell types: role of TIEG, Runx2, and other transcription factors. J Cell Biochem. 2008, 103: 383-392. 10.1002/jcb.21425.

  39. 39.

    Gaur T, Lengner CJ, Hovhannisyan H, Bhat RA, Bodine PV, Komm BS, Javed A, van Wijnen AJ, Stein JL, Stein GS, Lian JB: Canonical WNT signaling promotes osteogenesis by directly stimulating Runx2 gene expression. J Biol Chem. 2005, 280: 33132-33140. 10.1074/jbc.M500608200.

  40. 40.

    Kovacic N, Lukic IK, Grcevic D, Katavic V, Croucher P, Marusic A: The Fas/Fas ligand system inhibits differentiation of murine osteoblasts but has a limited role in osteoblast and osteoclast apoptosis. J Immunol. 2007, 178: 3379-3389. 10.4049/jimmunol.178.6.3379.

  41. 41.

    Chen JR, Haley RL, Hidestrand M, Shankar K, Liu X, Lumpkin CK, Simpson PM, Badger TM, Ronis MJ: Estradiol protects against ethanol-induced bone loss by inhibiting up-regulation of receptor activator of nuclear factor-kappaB ligand in osteoblasts. J Pharmacol Exp Ther. 2006, 319: 1182-1190. 10.1124/jpet.106.109454.

  42. 42.

    Messeguer X, Escudero R, Farre D, Nunez O, Martinez J, Alba MM: PROMO: detection of known transcription regulatory elements using species-tailored searches. Bioinformatics. 2002, 18: 333-334. 10.1093/bioinformatics/18.2.333.

  43. 43.

    Farre D, Roset R, Huerta M, Adsuara JE, Rosello L, Alba MM, Messeguer X: Identification of patterns in biological sequences at the ALGGEN server: PROMO and MALGEN. Nucleic Acids Res. 2003, 31: 3651-3653. 10.1093/nar/gkg605.

  44. 44.

    Zheng Q, Zhou G, Morello R, Chen Y, Garcia-Rojas X, Lee B: Type X collagen gene regulation by Runx2 contributes directly to its hypertrophic chondrocyte-specific expression in vivo. J Cell Biol. 2003, 162: 833-842. 10.1083/jcb.200211089.

  45. 45.

    Schroeder TM, Jensen ED, Westendorf JJ: Runx2: a master organizer of gene transcription in developing and maturing osteoblasts. Birth Defects Res C Embryo Today. 2005, 75: 213-225. 10.1002/bdrc.20043.

Download references


We sincerely thank Drs. Yao-Ming Chang and Te-Chin Chu for calculating bioinformatics analyses, study nurse Hsiang-Min Kung for gaining clinical sample and clinical information and excellent technical assistance from Li-Wen Lo for using of Typhoon Variable Mode Imager, TopCount NXT, and Real-Time PCR System at the Division of Medical Biology, Genomics Research Center, Academia Sinica, Taipei, Taiwan. This work was supported by Academia Sinica, Chang Gung Medical Foundation and National Science Council grants (AL Yu, NSC97-2323-B-001 and NSC 100-2321-B-001-037).

Author information

Correspondence to Wen-Hsiung Li or Alice Lin-Tsing Yu.

Additional information

Competing interests

The authors declare that they have no competing interests.

Electronic supplementary material

Additional file 1: Table S1.: Correlative analysis of expression of miR-10a, miR-10b and RUNX2 with clinical parameters of breast cancer patients. Table S2. ROC analysis of expression of miR-10a, miR-10b, RUNX2 for predicting outcome in breast cancer patients. Table S3. The relationship between overexpression of individual gene or joint effects of miR-10a, miR-10b, and RUNX2 with prognosis in breast cancer patients. Table S4. Cox regression model analyses of the gene expression regarding OS of breast cancer patients. Table S5. Predicted RUNX2 binding sites at proximal sequences preceding pre- miR-10a and pre-miR-10b. Figure S1. Prognostic significance of increased expression levels of RUNX2, miR-10a or miR-10b in breast cancer patients. Figure S2. Target gene expression of suppression or overexpression miR-10a and miR-10b. Figure S3. Luciferase reporter assay of putative RUNX2 target genes. (PDF 1 MB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Chang, C., Fan, T., Yu, J. et al. The prognostic significance of RUNX2 and miR-10a/10b and their inter-relationship in breast cancer. J Transl Med 12, 257 (2014).

Download citation


  • RUNX2
  • miR-10a
  • miR-10b
  • Breast cancer prognosis