RUNX2 facilitates aggressiveness and chemoresistance of triple negative breast cancer cells via activating MMP1

Breast cancer remains the most common malignancy in women and constantly threatens the lives of patients worldwide. State-of-the-art renewal has indicated the involvement of RUNX-associated transcription factor 2 (RUNX2) in tumorigenesis and cancer progression, yet the detailed information during breast cancer is largely obscure. Herein, we took advantage of breast cancer cell lines and in vivo tumorigenicity test as well as multifaceted phenotypic analyses (e.g., RNA-sequencing, ChIP and qRT-PCR assay) to verify the pathogenic mechanism of RUNX2 in triple negative breast cancer aggressiveness and chemoresistance. Strikingly, the proliferation, migration, invasion and chemoresistance of resistant cell lines in triple negative breast cancer was effectively suppressed by RUNX2 silencing, and the in vivo tumorigenicity was significantly weakened as well. Furthermore, with the aid of transcriptomic and bioinformatic analyses, we found MMP1 was highly expressed in triple negative breast cancer (TNBC) and showed a strong correlation with the poor prognosis of the patients, which was consistent with the expression pattern of RUNX2. Finally, by conducting ChIP and qRT-PCR assessment, we verified that RUNX2 functioned via directly binding to the specific motifs in the promoter of MMP1 and thus activating the transcriptional process. Collectively, our data demonstrated the facilitating effect of RUNX2 during triple negative breast cancer progression by directly orchestrating the expression of MMP1, which supplied overwhelming new references for RUNX2-MMP1 axis serving as a novel candidate for breast cancer diagnosis and treatment.

Introduction

Breast cancer is one of the most common malignancies and the leading cause of cancer related death in women, which also has drawn public-health concerns worldwide (1). Despite the dramatic progress in clinical practice, surgery and chemotherapy remain the major treatment options for breast cancer administration (2). As a particularly aggressive subtype, triple negative breast cancer (TNBC) has been regarded as the most challenging breast cancer without hormone receptors and human epidermal growth factor receptor 2 (HER2) gene amplification, which commonly results in greater metastatic potential, higher rates of relapse, together with shorter overall survival compared with other subtypes (3). Worse still, due to the high metastatic nature and tumor heterogeneity as well as drug resistance, breast cancer is still intractable for clinicians (4). Thus, there’s an urgent need to identify novel biomarkers and dissect the underlying molecular mechanisms of metastasis and drug resistance, which will provide novel therapeutic targets for drug development and the clinical treatment decisions of breast cancer patients.

RUNX2 is a transcription factor that has been reported essential for skeletal development, chondrocyte maturation and osteoblast differentiation (5–7). To date, literatures in the field have also suggested the involvement of RUNX2 in various malignant progression. For instance, He et al. and Li et al. recently indicated the aberrant phosphorylation and expression of RUNX2 in breast cancer invasion and gastric cancer metastasis, respectively (8, 9). Of note, Vishal and the colleagues further reviewed the potential dysfunction and regulation of RUNX2 in breast cancer mediated bone metastasis, yet the regulatory mechanism during the metastatic process and drug resistance in triple negative breast cancer was insufficiently disclosed (10).

For the purpose, we conducted multifaceted analyses both in vitro and in vivo to verify the dysfunction and the concomitant mechanism of RUNX2 in modulating proliferation, migration, invasion, and chemoresistance of breast cancer cells. Notably, with the aid of transcriptomic analysis and ChIP assay, we found that RUNX2 was adequate to regulating MMP1 expression by directly binding to the specific motifs in the promoter region. Compared to normal breast epithelium, triple negative breast cancer cells and tissues revealed higher level of MMP1 expression, which also showed perfect positive correlation with poor prognosis in breast cancer patients. Taken together, our findings highlighted the pathogenic mechanism of RUNX2-MMP1 axis in breast cancer, which further provided new insights into metastasis and chemoresistance of triple negative breast cancer.

Materials and methods

Cell lines

MDA-MB-231 and SUM-149 were purchased from the American Type Culture Collection (ATCC, USA), and MDA-MB-231 was cultured in high glucose DMEM medium (PM150210, Procell) supplemented with 10% fetal bovine serum (164210, Procell) and 1% penicillin and streptomycin (PB180120, Procell). SUM-149 was cultured in Ham’s F12 medium (L410KJ, BasalMedia) supplemented with 5% fetal bovine serum, 5µg/ml of insulin (UNITED LABORATORIES), and 1µg/ml of hydrocortisone (R011855, RHAWN). MDA-MB-231-Re epirubicin-resistant cell lines were established as previously described (11). All cells were cultured in 5% CO₂, at 37°C.

Lentiviral stable cell lines

MDA-MB-231-Re cells were seeded in six-well plates and incubated for 24 h. Lentivirus containing shRNA (shRUNX2, NC; MOI=10) was added when cell reached 20-35% confluency. Then, the supernatant was replaced with fresh medium after 10 h, followed by the observation of GFP fluorescence under a microscope at 72 h and stable line selection by puromycin addition. MDA-MB-231-Re cells with shRNA NC infection were used as a negative control.

Similarly, the RUNX2 over-expression (denoted as RUNX2) and the negative control (denoted as Vector) SUM-149 cell lines were generated by utilizing the aforementioned lentiviral strategy.

SiRNA transfection

For RUNX2 interference with the specific siRNAs (siRUNX2#1, siRUNX2#2) (GenePharma), MDA-MB-231 cells were cultured in 60-mm dishes for 24 h and the medium was replaced with fresh medium containing 200 pmol siRNA and 10µl Lipo6000™ Transfection Reagent (C0526, Beyotime) for 5 h when cells reached 30-50% confluence. Then, the siRNA-Transfected MDA-MB-231 cells was cultured in fresh medium, and the efficiency of interference for RUNX2 (siRUNX2#1, siRUNX2#2) and negative mimic (NC) was verified 48 h later by utilizing the qRT-PCR assay.

Cell counting kit-8 assay

To verify the cytotoxicity of the indicated anti-cancer drug at various concentrations, the cells were seeded in 96-well plates at a density of 5000 cells/well in 100 µl culture medium and incubated at 37°C, 5% CO₂ for 24 h. Then, the cells were incubated with anti-cancer drug epirubicin (Pfizer) for 24 h, and 10 µl/well CCK-8 solution (GK10001, GlpBio) was added for another 2 h. Finally, the absorbance value was read at 450 nm by using a microplate reader (Thermo).

RNA isolation and quantitative real-time PCR

Total RNAs were isolated by using TRIzol Reagent (CW0580S, CWBIO) according to the manufacturer’s instructions. Then, cDNAs were synthesized with the Reverse Transcriptase kit (E047-01B, Novoprotein) as we previously reported (12). After that, qRT-PCR analysis was conducted by using the SYBR Green Pro Taq HS premixed qPCR kit (AG11701-S, Accurate Biotechnology) and the ABI PRISM 7900 (Applied Biosystems). In details, an initial denaturation step of 95°C for 30s, followed by 40 cycles of 95°C for 5s and 60°C for 30s. The results were calculated by using 2^-ΔΔCt method. The primer sequences were shown in Supplementary Table S1. GAPDH was used as internal reference.

Protein extraction and western blotting assay

Total protein was extracted with RIPA lysis buffer from the indicated cells, and the concentration was measured with the BCA Protein Assay Kit (P0012, Beyotime). The samples were subjected to 10% sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and proteins were transferred to a PVDF membrane (IPVH00010, Millipore). After that, the membranes were incubated with primary antibody overnight at 4°C and secondary antibody for 1 h at room temperature, respectively. Finally, the protein bands were developed by using an enhanced chemiluminescence detection reagent (D045, Bridgen) and imaged by using an Image Lab (BIO-RAD). Primary antibodies against RUNX2 (12556, CST, dilution 1:1000), MMP1 (DF6325, Affinity, dilution 1:1000), E-CADHERIN (20874-1-AP, Proteintech, dilution 1:5000), VIMENTIN (10366-1-AP, Proteintech, dilution 1:5000), and β-ACTIN (AF7018, Affinity, dilution 1:5000) were used. Three independent experiments were conducted.

Transwell assay

Cells preconditioned in DMEM basal medium (without 10% FBS) for 12 h were seeded in the upper chamber of the 24-well plate transwell insert (353097, Corning) at a density of 10⁵ cells/well, and 600 µl of DMEM complete medium (with 10% FBS) were added to the lower chamber. Then, the cells were fixed with 10% formalin (YULU) for 20 min after incubation for 24- 48 h, following by staining with 0.1% crystal violet (C8470, Solarbio) for 20 min. After wiping the non-migratory cells in the inner side of the upper chamber, images of migrated cells were captured from three randomly selected fields under the microscope, and the number of cells in each field was calculated.

The invasion assay was similar to the migration assay, where the upper chamber was coated with matrix gel (356234, Corning).

Plate cloning assay

Breast cancer cells were seeded in 6-well plates at 1000 cells/well, and cultured for 14 days at 37°C, 5% CO₂. The cells were fixed with 10% formalin for 20 min, followed by staining with 0.1% crystal violet for 20 min. Typical colonies (foci > 100 µm) were counted. Three independent experiments were conducted.

Soft agar colony formation assay

To evaluate the malignant proliferative capacity of the indicated breast cancer cells, we conducted the soft agar colony formation assay. Briefly, the bottom layer was composed of 1.2% low melting point agarose (BS144, Biosharp) for preventing cell adhesion, and the upper layer was made up of 0.7% low melting point agarose containing 2000 cells/well. After 14 days, the morphology and count of colonies were recorded under the microscope.

Xenograft assay and histological analysis

5-week-old female BALB/c nude mice were applied for xenograft assay according to the Guide for the Care and Use of Laboratory Animals and were approved by the Experimental Animal Ethics Committee of Anhui Medical University (Approval No.: LLSC20210808). 5×10⁶ NC or shRUNX2-transduced MDA-MB-231-Re cells were subcutaneously injected into 5-week-old female BALB/c nude mice (N=4) in each group. Tumor long diameter (L) and short diameter (W) were measured every 3 days. The volume of the xenograft was calculated according to the formula: volume (mm³(=L×W²×0.5. Finally, the morphology and weight of the xenograft in nude mice in the indicated groups were recorded at day 33 of the xenograft assay. Besides, mice were sacrificed and tumors were resected for immunohistochemical staining.

RNA-sequencing and bioinformatic analyses

For preparation of RNA-SEQ samples, total RNAs were isolated by using the TRIzol Reagent according to the manufacturer’s instructions. Then, the RNAs were sequenced by Personal Bio Company with the Illumina mRNA-Seq sample preparation kit and the Illumina Hiseq sequencing platform (Illumina, Inc.). The multifaceted bioinformatic analyses, including HeatMap, volcano plots, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were accomplished as we recently reported with several modifications (13).

Chromatin immunoprecipitation and qRT-PCR assay

Cells were crosslinked with 37% formaldehyde and lysed in SDS buffer. Then, the DNA was fragmented by using sonication. For RUNX2 binding sites on MMP1 promoter, ChIP assay was performed as we previously reported (14).

The eluted DNA fragments were turned to qRT-PCR assay, and the values generated from immunoprecipitated samples were normalized to those of the corresponding input samples. Besides, qRT-PCR products were stored at -20°C in a refrigerator until subsequent agarose gel electrophoresis (AGE). Primers sequences were shown in Supplementary Table S2.

Bioinformatic analysis

For the purpose of exploring the correlation of MMP1 in breast cancer, we analyzed the expression level of MMP1 according to the TCGA datasets. In details, we verified the MMP1 expression level in subtypes of breast cancer by the UALCAN database (http://ualcan.path.uab.edu). To assess the potential clinical significance of MMP1 in breast cancer, we took advantage of the Kaplan-Meier Plotter database (http://kmplot.com/analysis/index.php). The potential binding sites of RUNX2 in the promoter region of MMP1 were predicted via utilizing the JASPAR database (http://jaspar.genereg.net/).

Statistical analysis

Student’s t-test was used to analyze the differences between the two groups. All statistical analyses were performed with the GraphPad Prism (GraphPad Software) and SPSS 20.0 software as we described before (14). Only when P<0.05 was considered statistically significant differences. All data were shown as mean ± SD (n=3 independent experiments). *,P<0.05; **, P<0.01; ***, P<0.001; ****, P<0.0001.

Results

MDA-MB-231-Re cells manifested increased resistance and facilitated EMT compared with parental cells

To generate the triple negative breast cancer resistant cell line MDA-MB-231-Re cells, we took advantage of the chronic-induced MDA-MB-231 cells for 8 months by step-up therapy (11). According to the CCK-8 values of the half of inhibition rate (IC₅₀), we found the drug resistance of MDA-MB-231-Re cells was higher than that of MDA-MB-231 parental cells (15.10 nmol/l vs. 2.035 nmol/l). In details, the Epirubicin drug resistance index of MDA-MB-231-Re cells was 7.42-fold over that of the MDA-MB-231 cells (Figure 1A). Meanwhile, we noticed the visible differences in cytomorphology between the MDA-MB-231-Re cells and MDA-MB-231 parental cells (Figure 1B). Distinguish from the epithelial-like shape in the MDA-MB-231 group, MDA-MB-231-Re cells showed typical mesenchymal-like morphology instead (Figure 1B). To investigate whether MDA-MB-231-Re cells could facilitate EMT, qRT-PCR and western-blotting analyses were performed to detect the expression pattern of EMT-associated biomarkers in the MDA-MB-231-Re cells and MDA-MB-231 parental cells. As show in Figures 1C, D, the expression of epithelial-related biomarker E-CADHERIN was reduced, while the mesenchymal-related biomarker VIMENTIN was up-regulated. Taken together, the triple negative breast cancer resistant MDA-MB-231-Re cells revealed enhanced drug-resistance against Epirubicin and facilitated EMT compared to the parental MDA-MB-231 cells.

FIGURE 1

Figure 1 Drug resistance and cytomorphology of the drug-resistant triple negative breast cancer cell lines. (A) Cell viability detection after 24 h of Epirubicin treatment at various concentrations by CCK-8 assay. (B) Cell morphology was observed under a microscope. Scale bar =100 µm. (C) qRT-PCR assay for analyzing the expression of E-CADHERIN and VIMENTIN in MDA-MB-231-Re cells or MDA-MB-231-Pa cells. (D) Protein levels of E-CADHERIN, VIMENTIN and β-ACTIN were determined by western-blotting assay. All data were shown as mean ± SD (N=3 independent experiments). **P <0.01.

RUNX2 was essential for migration and invasion of triple negative breast cancer cells

To assess the influence of RUNX2 in the migration and invasion of triple negative breast cancer, we transduced lentivirus-mediated shRUNX2 into MDA-MB-231-Re cells (denoted as shRUNX2) to knockdown the endogenous expression of RUNX2. By conducting qRT-PCR and western-blotting analyses, we found the expression of RUNX2 both at the transcription and translation level was consistently impaired in the shRUNX2 group compared with the negative control group (denoted as NC) (Figures 2A, B). Strikingly, MDA-MB-231-Re cells with stable RUNX2 knockdown restored the epithelial-like cytomorphology, which was distinguish from that in the NC group (Figure 2C).

FIGURE 2

Figure 2 Influence of shRUNX2 on the migration and invasion of drug-resistant cells in triple negative breast cancer. (A) qRT-PCR analysis of RUNX2 expression in MDA-MB-231-Re cells with stable RUNX2 knockdown (shRUNX2) or Scramble-transduced negative control (NC). (B) Protein levels of RUNX2 and β-ACTIN were determined by western-blotting assay. (C) Cytomorphology of the shRUNX2 and NC MDA-MB-231-Re cells was observed under the microscope. Scale bar =100 µm. (D) qRT-PCR assay for analyzing the expression of E-CADHERIN and VIMENTIN in the shRUNX2 and NC MDA-MB-231-Re cells. (E) protein levels of E-CADHERIN, VIMENTIN and β-ACTIN were determined by western-blotting assay. (F) The migration ability of the shRUNX2 and NC MDA-MB-231-Re cells assessed by Transwell assay. Scale bar =100 um. (G) The invasion ability of the shRUNX2 and NC MDA-MB-231-Re cells assessed by Transwell assay. Scale bar =100 µm. All data were shown as mean ± SD (N=3 independent experiments). *P <0.05; **P <0.01; ****P<0.0001.

Considering the switching action in orchestrating the cytomorphology, we supposed that RUNX2 might play a pivotal role in the progression of epithelial-mesenchymal transition (EMT). Therefore, qRT-PCR and western-blotting analyses were performed to verify the expression pattern of EMT-associated biomarkers in the shRUNX2 group and NC group, respectively. Intuitively, we found that RUNX2 knockdown increased the expression of E-CADHERIN, while reduced the expression of VIMENTIN (Figures 2D, E). Consistently, we took advantage of the transwell assay to evaluate the influence of RUNX2 knockdown upon MDA-MB-231-Re cells. Notably, the migration and invasion ability of the MDA-MB-231-Re cells with RUNX2 knockdown was consistently impaired (Figures 2F, G). Collectively, RUNX2 knockdown was adequate to inhibit EMT and reduce the in vitro migration and invasion in the triple negative breast cancer resistant MDA-MB-231-Re cells.

Knockdown of RUNX2 inhibited the proliferation, drug resistance, and anchorage-independent growth in triple negative breast cancer cells

Having clarified the inhibitory effect of RUNX2 knockdown upon cytomorphology as well as migration and invasion, we were next curious about the potential impact upon proliferation, drug resistance and anchorage-independent growth in MDA-MB-231-Re cells. With the aid of plate cloning assay, we found that fewer clones were formed in the shRUNX2 group compared to that in the NC group (Figure 3A). CCK-8 assay showed that the IC₅₀ of MDA-MB-231-Re cells with shRUNX2 transfection was much lower than that in the NC group (4.875 nmol/l vs. 17.74nmol/l) as well, which indicated the weakened drug resistance according to the epirubicin resistance index (RI=27.48%) (Figure 3B). In consistence, soft agar colony formation assay further revealed the reduced anchorage-independent growth capacity of shRUNX2 cells compared to NC cells (Figure 3C). Taken together, the proliferation, drug resistance, and anchorage-independent growth in triple negative breast cancer MDA-MB-231-Re cells were consistently impaired by RUNX2 knockdown.

FIGURE 3

Figure 3 Influence of RUNX2 on the proliferation and resistance of drug-resistant cells in breast cancer. (A) Proliferative capacity of the shRUNX2 and NC MDA-MB-231-Re cells was determined by plate cloning assay. (B) Resistance to the shRUNX2 and NC MDA-MB-231-Re cells was determined by the CCK-8 assay. (C) The anchorage-independent growth capacity of the shRUNX2 and NC MDA-MB-231-Re cells was determined by a soft agar colony formation assay. Scale bar =100 µm, 50 µm. All data were shown as mean ± SD (N=3 independent experiments). **P<0.01; ***P<0.001.

Knockdown of RUNX2 efficaciously inhibited tumor growth in nude mouse

To further verify the oncogenic potential of RUNX2 in vivo, the aforementioned NC cells and shRUNX2 cells were subcutaneously injected into BALB/c nude mice for xenograft tumor formation. Distinguish from those formed in the NC group, tumor sizes in the shRUNX2 group were significantly decreased at day 33 after transplantation (Figures 4A, B). Furthermore, as shown by the spatio-temporal growth curves, the tumor volume and tumor growth rate of shRUNX2 cells were less than that of the NC cells (Figure 4C). Similarly, the tumor weight in the shRUNX2 group was also sharply declined when compared with the NC group (Figure 4D). Overall, our data indicated the crucial role of RUNX2 in tumorigenesis during triple negative breast cancer.

FIGURE 4

Figure 4 RUNX2 promoted the progression of breast cancer cells in vivo. (A) Image of in vivo subcutaneous neoplasia in BALB/c nude mice formed by the transplanted shRUNX2 and NC MDA-MB-231-Re cells. (B) Image of the tumor nodules formed by the transplanted shRUNX2 and NC MDA-MB-231-Re cells. (C) Tumor volume was measured in each group on the indicated days. (D) Tumors were removed after the last measurement on the 33rd day and weighed. All data were shown as mean ± SD (N=3 independent experiments). *P<0.05.

shRUNX2 cells revealed discrepancy in gene expression profiling with NC cells

Having verified the cellular characteristics of shRUNX2 cells and NC cells for triple negative breast cancer, we next turned to dissect the underlying molecular mechanisms. Therewith, we took advantage of the RNA-SEQ analysis to explore the potent similarities and differences at the transcriptomic level. As shown by the principal component analysis (PCA), the three independent duplicate sample (rep1, 2, 3) in the shRUNX2 group (in pink) or the NC group (in green) were clustered, respectively (Figure 5A). Hierarchical cluster analysis further showed the discrepancy in the profiling of differentially expressed genes (DEGs) between the indicated two groups (Figure 5B). Furthermore, the volcano and MA plots intuitively revealed the distribution of the DEGs including 896 up-regulated genes and 217 down-regulated genes (|log2Fold Change|>1, P <0.05) (Figures 5C, D).

FIGURE 5

Figure 5 Gene expression profiling and DEGs in shRUNX2 and NC MDA-MB-231-Re cells. (A) Principal component analysis (PCA) of the three replications of shRUNX2 and NC MDA-MB- -231-Re cells. (B) Hierarchical cluster analysis of the DEGs in shRUNX2 and NC MDA-MB-231-Re cells. (C, D) Volcano plots (C) and MA plots (D) of the total genes expressed in shRUNX2 and NC MDA-MB-231-Re cells.

MMP1 expression was significantly downregulated in shRUNX2 cells

Aiming to clarify the biological significance of the DEGs, we performed top gene ontology (topGO) analysis and observed the enrichment of certain datasets such as extracellular matrix (cell components, short for CC), protein binding (molecular function, short for MF), and regulation of cell population proliferation (biological processes, short for BP) (Figure 6A). Meanwhile, EMT-associated subsets were enriched by conducting GO enrichment analysis, including extracellular region, cell periphery, cell communication, regulation of cell population proliferation (Figure 6B).

FIGURE 6

Figure 6 Function and signal prediction of the DEGs in TNBC. (A, B) Bar chart of the top 10 GO terms (A) and bubble plots of the top 20 GO terms (B) in each GO classification based on the enriched DEGs. (C, D) Bar chart (C) and bubble diagram (D) of the top 20 KEGG pathways based on the enriched DEGs. (E) qRT-PCR analysis of candidate gene expression in MDA-MB-231-Re cells in shRUNX2 and NC MDA-MB-231-Re cells. All data were shown as mean ± SD (N=3 independent samples). **P<0.01; ****P<0.0001; ns, not significant.

Subsequently, we took advantage of the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database to further dissect the underlying variations accompanied with the DEGs. Notably, a certain number of signaling pathways associated with proliferation and tumorigenesis were enrich based on the -log10 (P value) such as PI3K-Akt, MAPK and Jak-STAT signaling pathways (Figure 6C). In consistence, KEGG pathway enrichment assay based on FDR values (FDR<0.001) further confirmed the specific enrichment of pathways in cancer, Jak-STAT, MAPK and PI3K-Akt signaling pathways, together with cell adhesion molecules(CAMs) (Figure 6D).

Therewith, we detected the mRNA expression pattern of extracellular matrix-, cell adhesion-, cell invasion and metastasis-associated DEGs in triple negative breast cancer MDA-MB-231-Re cells with stable NC or shRUNX2 transfection. Of the indicated candidates (e.g., TGF-β, MDR1, NTN4, NOTCH1), we noticed that the expression of MMP1 with the most downregulation in shRUNX2 cells when compared with the NC group (Figure 6E). Collectively, our data indicated the biological significance of DEGs were mainly related with EMT process (e.g., migration, invasion) and resistance. MMP1 might play a critical role in mediating the tumor-promoting effect of RUNX2 in triple negative breast cancer cells.

MMP1 expression manifested positive correlation with poor prognosis of breast cancer

For the purpose of further exploring the correlation of MMP1 in breast cancer, we detected the MMP1 expression level in 1098 breast cancer samples and 113 normal breast tissues according to the TCGA datasets. Notably, we found breast cancer tissues with higher level of MMP1 expression over that in normal breast tissue, which was confirmed by the UALCAN database (http://ualcan.path.uab.edu) (Figures 7A, B). Furthermore, compared with other types of breast cancer, we noticed that triple negative breast cancer with the highest level of MMP1 expression (Figure 7C).

FIGURE 7

Figure 7 Increased expression of MMP1 indicated a poor prognosis in breast cancer. (A) Relative expression level of MMP1 in normal breast tissue samples (n=113) and breast cancer samples (n=1098) according to the TCGA database. (B) Relative expression level of MMP1 in normal breast tissue samples (n=114) and breast cancer samples (n=1097) according to the ualcan database. (C) Expression level of MMP1 in various types of breast cancer samples according to the ualcan database. (D-F) The relation of MMP1 expression levels according to the kmplot database and overall survival (D), relapse free survival (E) and distant metastasis free survival (F) in breast cancer patients. **P<0.01; ***P<0.001; ****P<0.0001.

Meanwhile, we turn to assess the potential clinical significance in breast cancer. Based on the Kaplan-Meier Plotter database (http://kmplot.com/analysis/index.php), we found that MMP1 expression showed positive correlation with the overall survival, relapse free survival and distant metastasis free survival of patients with breast cancer (Figures 7D–F). Taken together, our data revealed the positive correlation of MMP1expression with poor prognosis of breast cancer.

RUNX2 directly targeted MMP1 in triple negative breast cancer cells

To further illustrate the molecular mechanism of MMP1 and RUNX2 in TNBC, we initially examined the expression pattern of RUNX2 and MMP1 in the SUM-149 and MDA-MB-231 triple negative breast cancer cell line and the MCF-7 non-triple negative breast cancer cell line. Consistently, we found that the mRNA expression levels of RUNX2 and MMP1 revealed similar expression pattern, and in particular, those in triple negative breast cancer cell line were significantly higher than in non-triple negative breast cancer cell line (Figure 8A). Furthermore, by conducting qRT-PCR and western-blotting assays, we found the expression levels of MMP1 and RUNX2 both at transcription and translation level were higher in MDA-MB-231-Re cells compared with those in MDA-MB-231 cells (Figures 8B, C). Then, with the aid of shRNA- and siRNA- mediated knockdown of RUNX2, we observed the sharp decrease of MMP1 expression in MDA-MB-231-Re and MDA-MB-231 cells compared to the NC groups (Figure 8D). Simultaneously, the expression of MMP1 was significantly upregulated in SUM-149 breast cancer cells with RUNX2 overexpression (Figure 8E). Besides, the expression levels of RUNX2 and MMP1 in tumor sections from both groups of xenograft tumor were detected by immunohistochemical staining, which showed that they were obviously lower in tumors of the shRUNX2 group compared with those of the NC group (Figure 8F).

FIGURE 8

Figure 8 MMP1 and RUNX2 showed positively correlation in relative mRNA and protein expression. (A) qRT-PCR analysis of the expression of RUNX2 and MMP1 in various breast cancer cell lines. (B, C) qRT-PCR and western-blotting analyses of RUNX2 (B) and MMP1 (C) expression in MDA-MB-231-Re cells or MDA-MB-231-Pa cells. (D, E) qRT-PCR analysis of MMP1 expression in RUNX2 knockdown MDA-MB-231-Re and MDA-MB-231 cells (D) as well as the SUM-149 with RUNX2 overexpression (E). (F) RUNX2 and MMP1 expression of tumor sections from each group were determined by immunohistochemistry. Scale bar =50 µm. *P <0.05; **P <0.01; ***P<0.001; ****P<0.0001.

Having preliminarily verified the potential regulatory relationship between RUNX2 and MMP1, we assumed whether RUNX2 functions via direct activation of MMP1 expression. Therefore, we predicted the potential binding sites of RUNX2 in the promoter region of MMP1 according to the JASPAR database (http://jaspar.genereg.Net/) (Figures 9A, B; Table 1). Of the seven candidates, the top 3 binding sites (site 1, 2, 3) with the highest score were selected, for further identification (Figure 9C). Then, with the aid of chromatin immunoprecipitation (ChIP) and qRT-PCR analyses together with agarose gel electrophoresis, we found only site 2 was specifically enriched by RUNX2 in MDA-MB-231-Re cells with high level of RUNX2 expression (Figures 9D, E). This was finally verified by agarose gel electrophoresis (Figure 9E). Overall, these data collectively indicated the direct modulation of MMP1 by RUNX2 via specific binding motif within the promoter region.

FIGURE 9

Figure 9 RUNX2 binds directly to the MMP1 promoter region. (A, B) The binding motif prediction of the RUNX2 according to the jaspar database. (C) Illustration of the three potential RUNX2 binding sites in the MMP1 promoter region. (D) ChIP-qRT-PCR analysis of the binding sites of RUNX2 in the MMP1 promoter region in MDA-MB-231-Re cells. (E) Agarose gel electrophoresis analysis of RUNX2 binding to the MMP1 promoter region in MDA-MB-231-Re cells. All data were shown as mean ± SD (N=3 independent samples). ***P<0.001; ns, not significant.

TABLE 1

Table 1 RUNX2 binding motif prediction in MMP1 promoter.

Discussion

The heterogeneous TNBC, characterized by the absence of progesterone receptor (PR) and HER2 expression, commonly results in poor prognosis due to the chemoresistance and infeasibility of targeted therapeutic approaches (3, 15). Therefore, despite chemotherapy serve as the “one size fits all” approach, yet the current treatment paradigm for TNBC administration is challenging and inopportune largely due to the deficiency of dissecting the clinical response and underlying molecular characteristics. Of note, state-of-the-art renewal has indicated the oncogenic role of RUNX2 in multiple cancer types (16–18), which the expression pattern and regulatory mechanism during chemoresistance in triple negative breast cancer are largely indistinct. For the purpose, we verified that RUNX2 knockdown was adequate to inhibit the in vitro proliferation, migration, invasion, and chemoresistance of resistant cells in triple negative breast cancer as well as the tumor growth in BALB/c nude mice. Furthermore, with the aid of transcriptomic analysis and ChIP-qRT-PCR assay, we identified that RUNX2/MMP1 axis functioned a crucial role in the malignant progression and poor prognosis in triple negative breast cancer (Figure 10).

FIGURE 10

Figure 10 Illustration of the working model of RUNX2-MMP1 axis in TNBC. In chemotherapy drug-resistant TNBC, RUNX2 expression was upregulated against Epirubicin treatment. Then, RUNX2 was adequate to directly bind to the specific motif in the promoter region of MMP1 and thus activate MMP1 expression. Finally, the upregulated MMP1 functioned critically in the degradation of ECMs and the resultant modulation of proliferation, metastasis (e.g., migration, invasion) and drug resistance of TNBC.

RUNX2 has been recognized as the basic transcription factor for various physiological development (e.g., skeletal and osteogenic development) (6), and the concomitant pathological processes including cancer development and tumor bone metastasis (19). For instance, Guo et al. demonstrated the oncogenic role of RUNX2 during gastric carcinogenesis via YAP1 signaling (20). Wang and colleagues verified the promoting effect of Cbf-β in colorectal cancer progression in a RUNX2-dependent manner (21). Previously, we demonstrated that RUNX2 was adequate to activate TGF-β signaling pathway, and thus performed a pivotal role in TNBC Epirubicin-resistance by orchestrating the stemness, EMT and apoptosis of CD44⁺/CD24^- BCSCs (11). Notably, current literature has suggested the critical role of RUNX2 for malignant progression of breast cancer by modulating the CD44 ⁺/CD24^- breast cancer stem cells (22, 23). Instead, we not only verified the tumorigenic role of RUNX2 in TNBC progression and cancer chemoresistance both in vitro and in vivo, but also revealed the regulatory mechanism of RUNX2 via directly targeting the specific binding site in the promoter region of MMP1.

Matrix metalloproteinases (MMPs) are zinc-dependent endopeptidases, which are modulated by tissue inhibitors of metalloproteinases (TIMPs) and involved in the degradation of a certain number of proteins in the extracellular matrix (ECM) (24). To date, MMPs have been reported crucially in many processes, including embryonic development and tissue repair, tumor growth and metastasis (25). For instance, Liu et al. and Yang et al. ascertained the MMP1-dependent effect in hepatocellular carcinoma (HCC) progression via interacting with HuR and the metastasis of colorectal cancer via orchestrating the DDX3/YY1/MMP1/PI3K-AKT axis, respectively (26, 27). It’s noteworthy that MMP1 upregulation by promoter hypomethylation has been reported with enhanced tamoxifen resistance in breast cancer (28). In this study, by conducting multifaceted cellular and molecular analyses, we for the first time identified a specific binding motif in the promoter of MMP1 for the direct transmission of the tumorigenic effect of RUNX2, which would supply novel target candidates for drug development and clinical intervention decision against the proliferation, metastasis (e.g., migration, invasion) and drug resistance of TNBC (29, 30).

Overall, our data indicated the oncogenic effect and concomitant mechanism of the RUNX2-MMP1 axis in the modulation of the proliferation, migration, invasion and drug resistance in triple negative breast cancer drug-resistant cells in vitro, together with tumor growth of breast cancer cells in vivo, which would supply overwhelming new references for dissecting the pathogenesis of TNBC and facilitating the development of clinical remedies and novel therapeutic targets for triple negative breast cancer.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Ethics statement

The animal study was reviewed and approved by the Ethics Committee of the Experimental Animal Ethics Committee of Anhui Medical University.

Author contributions

WS and FFL conceived and planned the study. WS, XDX, LW, and FXL performed the experiments. XDX, WW, XJX, and WL helped to prepare and analyze the samples. WS and XDX analyzed the data. WS and FFL wrote the manuscript. FFL, LZ, and DH completed critical revisions and proofread the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the National Natural Science Foundation of China (82173377, 81302319, 82260031), the Major Research Project of Education Department of Anhui Province (KJ2018ZD018), the Non-profit Central Research Institute Fund of Chinese Academy of Medical Sciences (2019PT320005), the 2021 Central-Guided Local Science and Technology Development Fund (ZYYDDFFZZZJ-1), the Anhui Provincial Natural Science Foundation of China (No.2008085MH299), the Key research and development Program of Anhui Province, China (No. 2022e07020012), the Fundamental Research Funds for the Central Universities (No. WK9110000086), Postdoctoral Research Funding of Anhui Province in 2019 (No. 2019B371).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2022.996080/full#supplementary-material

References

1. Ginsburg O, Bray F, Coleman MP, Vanderpuye V, Eniu A, Kotha SR, et al. The global burden of women's cancers: A grand challenge in global health. Lancet (2017) 389(10071):847–60. doi: 10.1016/S0140-6736(16)31392-7

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Van de Wiel M, Dockx Y, Van den Wyngaert T, Stroobants S, Tjalma WAA, Huizing MT. Neoadjuvant systemic therapy in breast cancer: Challenges and uncertainties. Eur J Obstet Gynecol Reprod Biol (2017) 210:144–56. doi: 10.1016/j.ejogrb.2016.12.014

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Garrido-Castro AC, Lin NU, Polyak K. Insights into molecular classifications of triple-negative breast cancer: Improving patient selection for treatment. Cancer Discov (2019) 9(2):176–98. doi: 10.1158/2159-8290.CD-18-1177

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Ferlay J, Colombet M, Soerjomataram I, Mathers C, Parkin DM, Pineros M, et al. Estimating the global cancer incidence and mortality in 2018: Globocan sources and methods. Int J Cancer (2019) 144(8):1941–53. doi: 10.1002/ijc.31937

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Komori T. Whole aspect of Runx2 functions in skeletal development. Int J Mol Sci (2022) 23(10):5776. doi: 10.3390/ijms23105776

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Komori T. Regulation of proliferation, differentiation and functions of osteoblasts by Runx2. Int J Mol Sci (2019) 20(7):1694. doi: 10.3390/ijms20071694

CrossRef Full Text | Google Scholar

7. Komori T. Regulation of skeletal development by the runx family of transcription factors. J Cell Biochem (2005) 95(3):445–53. doi: 10.1002/jcb.20420

PubMed Abstract | CrossRef Full Text | Google Scholar

8. He F, Matsumoto Y, Asano Y, Yamamura Y, Katsuyama T, La Rose J, et al. Runx2 phosphorylation by tyrosine kinase abl promotes breast cancer invasion. Front Oncol (2021) 11:665273. doi: 10.3389/fonc.2021.665273

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Li Y, Sun R, Zhao X, Sun B. Runx2 promotes malignant progression in gastric cancer by regulating Col1a1. Cancer biomark (2021) 31(3):227–38. doi: 10.3233/CBM-200472

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Vishal M, Swetha R, Thejaswini G, Arumugam B, Selvamurugan N. Role of Runx2 in breast cancer-mediated bone metastasis. Int J Biol Macromol (2017) 99:608–14. doi: 10.1016/j.ijbiomac.2017.03.021

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Xu X, Zhang L, He X, Zhang P, Sun C, Xu X, et al. Tgf-beta plays a vital role in triple-negative breast cancer (Tnbc) drug-resistance through regulating stemness, emt and apoptosis. Biochem Biophys Res Commun (2018) 502(1):160–5. doi: 10.1016/j.bbrc.2018.05.139

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Wei Y, Zhang L, Chi Y, Ren X, Gao Y, Song B, et al. High-efficient generation of vcam-1(+) mesenchymal stem cells with multidimensional superiorities in signatures and efficacy on aplastic anaemia mice. Cell Prolif (2020) 53(8):e12862. doi: 10.1111/cpr.12862

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Huo J, Zhang L, Ren X, Li C, Li X, Dong P, et al. Multifaceted characterization of the signatures and efficacy of mesenchymal Stem/Stromal cells in acquired aplastic anemia. Stem Cell Res Ther (2020) 11(1):59. doi: 10.1186/s13287-020-1577-2

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Zhang L, Wang H, Liu C, Wu Q, Su P, Wu D, et al. Msx2 initiates and accelerates mesenchymal Stem/Stromal cell specification of hpscs by regulating Twist1 and prame. Stem Cell Rep (2018) 11(2):497–513. doi: 10.1016/j.stemcr.2018.06.019

CrossRef Full Text | Google Scholar

15. Lyons TG. Targeted therapies for triple-negative breast cancer. Curr Treat Options Oncol (2019) 20(11):82. doi: 10.1007/s11864-019-0682-x

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Zhang D, Song Y, Li D, Liu X, Pan Y, Ding L, et al. Cancer-associated fibroblasts promote tumor progression by lncrna-mediated Runx2/Gdf10 signaling in oral squamous cell carcinoma. Mol Oncol (2022) 16(3):780–94. doi: 10.1002/1878-0261.12935

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Zhang C, Xu X, Trotter TN, Gowda PS, Lu Y, Suto MJ, et al. Runx2 deficiency in osteoblasts promotes myeloma resistance to bortezomib by increasing tsp-1-Dependent Tgfbeta1 activation and suppressing immunity in bone marrow. Mol Cancer Ther (2022) 21(2):347–58. doi: 10.1158/1535-7163.MCT-21-0310

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Hong D, Knelson EH, Li Y, Durmaz YT, Gao W, Walton E, et al. Plasticity in the absence of notch uncovers a Runx2-dependent pathway in small cell lung cancer. Cancer Res (2022) 82(2):248–63. doi: 10.1158/0008-5472.CAN-21-1991

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Yin X, Teng X, Ma T, Yang T, Zhang J, Huo M, et al. Runx2 recruits the Nurd(Mta1)/Crl4b complex to promote breast cancer progression and bone metastasis. Cell Death Differ (2022) 29(11):2203–17. doi: 10.1038/s41418-022-01010-2

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Guo Z, Zhou K, Wang Q, Huang Y, Ji J, Peng Y, et al. The transcription factor Runx2 fuels Yap1 signaling and gastric cancer tumorigenesis. Cancer Sci (2021) 112(9):3533–44. doi: 10.1111/cas.15045

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Wang C, Shi Z, Zhang Y, Li M, Zhu J, Huang Z, et al. Cbfbeta promotes colorectal cancer progression through transcriptionally activating opn, Fam129a, and Upp1 in a Runx2-dependent manner. Cell Death Differ (2021) 28(11):3176–92. doi: 10.1038/s41418-021-00810-2

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Zhang P, Liu L, Zhang L, He X, Xu X, Lu Y, et al. Runx2 is required for activity of Cd44(+)/Cd24(-/Low) breast cancer stem cell in breast cancer development. Am J Transl Res (2020) 12(5):2305–18.

PubMed Abstract | Google Scholar

23. Zhang L, Liu L, Xu X, He X, Wang G, Fan C, et al. Mir-205/Runx2 axis negatively regulates Cd44(+)/Cd24(-) breast cancer stem cell activity. Am J Cancer Res (2020) 10(6):1871–87.

PubMed Abstract | Google Scholar

24. Wang X, Khalil RA. Matrix metalloproteinases, vascular remodeling, and vascular disease. Adv Pharmacol (2018) 81:241–330. doi: 10.1016/bs.apha.2017.08.002

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Najafi M, Farhood B, Mortezaee K. Extracellular matrix (Ecm) stiffness and degradation as cancer drivers. J Cell Biochem (2019) 120(3):2782–90. doi: 10.1002/jcb.27681

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Liu H, Lan T, Li H, Xu L, Chen X, Liao H, et al. Circular rna Circdlc1 inhibits Mmp1-mediated liver cancer progression Via interaction with hur. Theranostics (2021) 11(3):1396–411. doi: 10.7150/thno.53227

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Yang P, Li J, Peng C, Tan Y, Chen R, Peng W, et al. Tcons_00012883 promotes proliferation and metastasis Via Ddx3/Yy1/Mmp1/Pi3k-akt axis in colorectal cancer. Clin Transl Med (2020) 10(6):e211. doi: 10.1002/ctm2.211

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Kim HW, Park JE, Baek M, Kim H, Ji HW, Yun SH, et al. Matrix metalloproteinase-1 (Mmp1) upregulation through promoter hypomethylation enhances tamoxifen resistance in breast cancer. Cancers (Basel) (2022) 14(5):1232. doi: 10.3390/cancers14051232

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Xiao D, Bi R, Liu X, Mei J, Jiang N, Zhu S. Notch signaling regulates mmp-13 expression Via Runx2 in chondrocytes. Sci Rep (2019) 9(1):15596. doi: 10.1038/s41598-019-52125-5

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Kim B, Kim H, Jung S, Moon A, Noh DY, Lee ZH, et al. A ctgf-Runx2-Rankl axis in breast and prostate cancer cells promotes tumor progression in bone. J Bone Miner Res (2020) 35(1):155–66. doi: 10.1002/jbmr.3869

PubMed Abstract | CrossRef Full Text | Google Scholar