- Research article
- Open Access
Regulation of breast cancer metastasis by Runx2 and estrogen signaling: the role of SNAI2
Breast Cancer Research volume 13, Article number: R127 (2011)
In contrast to its role in breast cancer (BCa) initiation, estrogen signaling has a protective effect in later stages, where estrogen receptor (ER)α loss associates with aggressive metastatic disease. We asked whether the beneficial effect of estrogen signaling in late-stage BCa is attributable to the recently reported estrogen-mediated antagonism of the pro-metastatic transcription factor Runx2.
MCF7/Rx2dox breast cancer cells were engineered with a lentivirus expressing Runx2 in response to doxycycline (dox). Cells treated with dox and/or estradiol (E2) were subjected to genome-wide expression profiling, RT-qPCR analysis of specific genes, and Matrigel™ invasion assays. Knockdown of genes of interest was performed using lentiviruses expressing appropriate shRNAs, either constitutively or in response to dox. Gene expression in BCa tumors was investigated using a cohort of 557 patients compiled from publicly available datasets. Association of gene expression with clinical metastasis was assessed by dichotomizing patients into those expressing genes of interest at either high or low levels, and comparing the respective Kaplan-Meier curves of metastasis-free survival.
Runx2 induced epithelial-mesenchymal transition (EMT) evidenced by acquisition of a fibroblastic morphology, decreased expression of E-cadherin, increased expression of vimentin and invasiveness. Runx2 stimulated SNAI2 expression in a WNT- and transforming growth factor (TGF)β-dependent manner, and knockdown of SNAI2 abrogated the pro-metastatic activities of Runx2. E2 antagonized the pro-metastatic activities of Runx2, including SNAI2 upregulation. In primary BCa tumors, Runx2 activity, SNAI2 expression, and metastasis were positively correlated, and SNAI2 expression was negatively correlated with ERα. However, the negative correlation between SNAI2 and ERα in bone-seeking BCa cells was weaker than the respective negative correlation in tumors seeking lung. Furthermore, the absence of ERα in primary tumors was associated with lung- and brain- but not with bone metastasis, and tumor biopsies from bone metastatic sites displayed the unusual combination of high Runx2/SNAI2 and high ERα expression.
E2 antagonizes Runx2-induced EMT and invasiveness of BCa cells, partly through attenuating expression of SNAI2, a Runx2 target required for mediating its pro-metastatic property. That ERα loss promotes non-osseous metastasis by unleashing Runx2/SNAI2 is supported by the negative correlation observed in corresponding tumors. Unknown mechanisms in bone-seeking BCa allow high Runx2/SNAI2 expression despite high ERα level
Metastasis of primary tumors to distant sites is a complex process that involves a sequence of interdependent events including intravasation, survival within the circulation, extravasation, and colonization. Epithelial-mesenchymal transition (EMT) has been strongly implicated in metastasis [1, 2]. During EMT, epithelial cells dissociate from each other, in part due to loss of E-cadherin expression, upregulate mesenchymal markers, acquire a fibroblast-like morphology, reorganize their cytoskeleton, and become more motile and invasive [1, 2]. Several transcription factors, including members of the SNAI family have been shown to promote EMT and thus tumor dissemination [3–5]. Bone metastasis is a frequent complication of breast cancer (BCa), with distinct gene signatures defining bone-seeking tumors [6–9].
Prolonged exposure to estradiol (E2) is associated with an increased risk of BCa [10–14]. The mechanisms through which estrogens contribute to BCa initiation and progression are complex and implicate estrogen receptor (ER)-mediated genomic and nongenomic signaling as well as the action of genotoxic estrogen metabolites . In contrast to E2-mediated carcinogenesis, the presence of ERα is a favorable prognostic marker associated with less invasive tumors, and those negative for ERα are more aggressive . A randomized clinical trial of postmenopausal women receiving equine estrogen treatment revealed a decrease in BCa incidence  and low-dose estradiol treatment has been proposed as a treatment modality for advanced ERα-positive BCa that does not respond to aromatase inhibition . Indeed, introduction of ERα into ERα-negative BCa cells attenuated their pro-cancerous properties in vitro [18–21]. Thus, anti-estrogen therapy for BCa patients, while antagonizing the oncogenic properties of estrogens, may inadvertently result in loss of their anti-metastatic property. Better understanding of mechanisms underlying the beneficial role of estrogen signaling in advanced disease may inform the development of novel therapeutic approaches and improved treatment plans for BCa patients.
Runx2 is a lineage-specific transcription factor with crucial roles in both bone biology and carcinogenesis [22–24]. During development Runx2 is involved in the process of osteogenesis. Targeted disruption of Runx2 in mice leads to failure of osteoblast differentiation and bone formation [25, 26], and Runx2 haploinsufficiency in humans results in the skeletal disorder cleidocranial dysplasia, with a similar phenotype observed in Runx2 haploinsufficient mice . Although Runx proteins have tumor suppressor properties , recent studies assigned a role for Runx2 in promoting breast and prostate cancer metastasis [27–32]. Thus, Runx2 and E2 signaling play dual roles in BCa, with each functioning to either promote or suppress tumor progression. The mechanisms underlying these contrasting manifestations in cancer are poorly understood.
We previously showed that in the presence of ligand, ERα physically binds Runx2 and inhibits expression of several Runx2 target genes ; our recent study revealed that in MCF7 BCa cells E2 independently regulated about half of the Runx2-responsive genes . We therefore hypothesized that gene(s) stimulated by Runx2 and inhibited by E2 may contribute to manifestations of the pro-metastatic or tumor suppressor functions of Runx2 and E2, respectively. Using a combination of tissue culture modeling and bioinformatics analysis of gene expression in BCa biopsies, we present evidence suggesting that SNAI2/SLUG plays a role in mediating the pro- and anti-metastatic effects of Runx2 and E2, respectively.
Materials and methods
Cell culture assays and reagents
MCF7 and T47D BCa cells were obtained from the American Type Culture Collection (Rockville, MD, USA). MCF7/Rx2dox cells, engineered to express Runx2 in response to doxycycline (dox) , were maintained in Dulbecco's modified Eagle's medium (DMEM) containing 10% fetal bovine serum (FBS) and 50 μg/ml Hygromycin B (GIBCO, Carlsbad, CA). Data presented herein were obtained from cells cultured in phenol red-free DMEM containing 10% charcoal-stripped FBS (Hyclone, South Logan, UT), 0.5 μg/ml dox (Calbiochem, San Diego, CA) and/or 10 nM E2 (Sigma, St Louis, MO). For invasion assays, cells further engineered to constitutively express luciferase were placed in Matrigel™-containing inserts (BD Biosciences, Bedford, MA) and invasion was assessed based on luciferase activity in cells that had crossed the membrane . The Wnt inhibitor ICG-001 was synthesized as previously described . Mouse monoclonal anti-SNAI2 and anti-Runx2 antibodies were from Millipore (Billerica, MA) and Invitrogen (Grand Island, NY), respectively. Anti-TGF-β type I receptor (TGFBRI) blocking antibody (sc-398) and negative control antibody (sc-2027) were from Santa Cruz Biotechnology Inc. (Santa Cruz, CA).
Real-time quantitative RT-qPCR analysis
Total RNA was extracted using Aurum™ Total RNA Mini-Kit (BioRad, Hercules, CA) and cDNA was synthesized with iScript Reverse Transcription Kit (BioRad). RT-qPCR was carried out in triplicate using iQ™ SYBR Green Supermix (BioRad) and a CFX96 qPCR machine (BioRad). All transcript levels were normalized to that of GAPDH. Primers used for PCR are listed in Additional file 1.
T47D/shRx2dox cells, conditionally expressing shRNA to silence Runx2 in response to dox, were previously described . For SNAI2 silencing, specific Mission-shRNA lentiviral plasmids (SHCLNG) and a negative control shRNA plasmid (SHC002) were purchased from Sigma and packaged as described  using HEK293T cells and the helper plasmids pMD.G1 and pCMV.R8.91. Viral particles were used for transduction of MCF7/Rx2dox cells followed by selection with 1 μg/ml Puromycin (Invitrogen, Carlsbad, CA). Of five shRNAs tested, two were found to reduce SNAI2 expression by > 75% and the corresponding cells were propagated as described above. shRNA sequences used in this study are listed in Additional file 1.
For primary analysis of gene expression profiles in clinical samples, we compiled microarray data from 557 BCa patients by combining three Gene Expression Omnibus (GEO) datasets, GSE2034, GSE2603 and GSE12276 [37–39]. These datasets were generated using Affymetrix chips HG-U133A (GSE2034 and GSE2603) and HG-U133plus2 (GSE12276). The Affymetrix CEL files were downloaded into Partek, normalized using the RMA (Robust Multiarray Averaging) method and filtered for common probes using interplatform comparison. The data were standardized by centering and scaling and patients were dichotomized into low and high expressing groups with regard to the indicated genes of interest. Dichotomizing of the tumors according to the ESR1 Affymetrix probe '205225_at' was consistent with the histological definition of their ER status in 329 of 368 cases (89.4%) for which pathological annotation was available, and with the molecular definition of their ER status in 177 of 189 cases (93.7%) where it was based on different processing of the microarray data . Clinical features, including hormonal status, metastasis-free survival time and sites of metastasis, have been compiled from the supplementary information of the published studies [9, 37]. For each patient, metastasis-free survival was defined as the time interval between surgery and the diagnosis of metastasis. Only tumors with unambiguous information about sites of metastasis were included in the cohort. For metastatic tissue analyses we used the microarray data from 58 archival human breast carcinoma metastasis specimens deposited into GEO under accession number GSE14020 . Because all the data were from publicly available resources, we did not require ethical approval to carry out this study.
Survival analyses and hierarchical clustering were performed using Partek Genomics Suite 6.6 and the correlation study was performed using Graphpad. Functional annotation of gene groups was performed using both the web-based Database for Annotation, Visualization, and Integrated Discovery (DAVID) system  and the Ingenuity Pathways Analysis (IPA™) package.
Metastasis-related genes in BCa cells are stimulated by Runx2 and inhibited by E2
A recent study with prostate cancer cells and a dox-inducible lentiviral system identified a Runx2-regulated metastasis-related gene network . We used the same system to conditionally express Runx2 in the ER-positive/Runx2-negative MCF7 BCa cells, and determined the genes that were responsive to Runx2 (N = 1,081) or E2 (N = 2,435) . Interestingly, almost half of the Runx2 responsive genes were also E2-responsive (N = 530) and these 'common' genes are listed in Additional file 2 along with their level of responsiveness to Runx2 and E2. Functional computational analysis linked these 530 'common' genes to angiogenesis and cell movement (Table 1), suggesting their importance in tumor metastasis. A literature survey for the top ten 'common' genes (Table 2) identified three -- SNAI2, S100A9 and CXCL12, which have been implicated in cell motility and BCa metastasis [4, 41, 42]. Moreover, SNAI2 and S100A9 were downregulated by E2. These genes may contribute to the pro- and anti-invasive properties of Runx2 and E2, respectively.
Functional link between Runx2, E2 and SNAI2 in BCa
Of the top 'common' genes, SNAI2 was of particular interest not only because its expression was the second most stimulated by Runx2 (6.6-fold) but also because it was strongly inhibited by E2 (-3.4-fold) (Table 2). Furthermore, the SNAI family of transcription factors is known to play a role in EMT [3–5]. By Western blot analysis, SNAI2 protein was detected in MCF7/Rx2dox cells only after the induction of Runx2 (Figure 1A). Additional support for the regulation of SNAI2 by Runx2 was provided by RT-qPCR analysis of other cell types, namely C4-2B/Rx2dox prostate cancer cells, in which SNAI2 was stimulated by dox-induced Runx2 (Figure 1B), and T47D/shRx2dox BCa cells, in which SNAI2 expression was inhibited after dox-induced shRNA-mediated silencing of endogenous Runx2 (Figure 1C). We then sought evidence for the positive and negative regulation of SNAI2 by Runx2 and ERα, respectively, in BCa tumors. We compiled a cohort of 557 BCa patients from the publicly available datasets GSE2034, GSE2603 and GSE12276 [37–39] and initially investigated the correlation between SNAI2 expression and a Runx2 metagene, the latter defined as the average normalized expression of all genes that Runx2 stimulated by ≥ 2-fold in the MCF7/Rx2dox system (except SNAI2 itself). As shown in Figure 1D, there was a positive correlation between expression of SNAI2 and the Runx2 metagene across the 557 BCa patient cohort. On the other hand, SNAI2 expression in the clinical specimens was negatively correlated with ERα/ESR1: tumors with low ESR1 expression (mostly histologically ER-negative) had significantly more SNAI2 transcripts than those with high ESR1 expression (Figure 1E), and the latter (mostly ER-positive tumors) exhibited a negative correlation between ESR1 and SNAI2 expression (Figure 1F). These data suggest that SNAI2 is stimulated by Runx2 and inhibited by E2 not only in the MCF7/Rx2dox model (Table 2), but also in other cell lines, as well as in BCa tumors.
E2 antagonizes Runx2-induced EMT and invasiveness in vitro
Because SNAI2 has been implicated in EMT and in the invasive metastatic behavior of breast and other cancer cells [3, 4] we employed the MCF7/Rx2dox culture system to test the hypotheses that Runx2 induces and E2 antagonizes EMT and invasiveness, and that SNAI2 plays a role in these processes. We initially measured by RT-qPCR the expression levels of EMT and motility-related genes in MCF7/Rx2dox cells treated with dox, E2, or both together. Consistent with the microarray results (Table 2), SNAI2 expression was increased by 9.7 fold and repressed by 4.6 fold in response to 48-hour treatment with dox and E2, respectively (Figure 2A). The mRNA for another well-established marker of motility and invasion, S100A4 [43, 44], was also stimulated by Runx2. Although E2 did not affect basal S100A4 expression, it completely blocked the Runx2-mediated induction (Figure 2A). Despite the stimulation of SNAI2 and S100A4 by Runx2, these changes in gene expression were not accompanied by either noticeable morphological changes (data not shown) or a decrease in the epithelial marker E-cadherin (CDH1) or an increase in the mesenchymal marker vimentin (VIM) (Figure 2A). Taken together these results suggest that induction of Runx2 for two days at best initiated partial EMT in MCF7 cells.
We further assessed expression of EMT and motility-related genes seven days after Runx2 induction in MCF7/Rx2dox cells. As shown in Figure 2B, this prolonged induction resulted in a more complete EMT. Not only were SNAI2 and S100A4 stimulated, CDH1 was downregulated and VIM was upregulated. Finally, microscopic examination of cultures seven days after Runx2 induction demonstrated a fibroblast-like and scattered morphology as compared to control cultures (Figure 2C). Thus, the partial EMT indicated by increased SNAI2 and S100A4 expression on day 2 (Figure 2A) developed into a more comprehensive EMT phenotype by day 7. More importantly, E2 significantly antagonized most of the Runx2-induced EMT characteristics. It blunted the Runx2-mediated effects on SNAI2 and CDH1 expression, and attenuated the Runx2-mediated effects on the expression of VIM and S100A4 (Figure 2B), as well as its effects on cell morphology (Figure 2C).
Finally, we measured the effects of Runx2 and E2 on BCa cell invasiveness. MCF7/Rx2dox cells were transduced with lentiviruses constitutively expressing firefly luciferase and placed on top of Matrigel™ membranes within BD Biocoat™ inserts. Luciferase activity associated with cells that had crossed the membrane, a measure of invasiveness , was two-fold higher in cultures treated with dox to induce Runx2 as compared to controls, while co-treatment with E2 attenuated the Runx2-mediated invasiveness (Figure 2D). E2 alone had little, if any, effect on MCF7/Rx2dox cell invasiveness through Matrigel™.
SNAI2 plays a critical role in Runx2-mediated MCF7 cell invasion
We next tested the hypothesis that SNAI2 was required for Runx2-mediated EMT and invasion in MCF7/Rx2dox BCa cells. To this end, MCF7/Rx2dox BCa cells were transduced with lentiviruses expressing two independent shRNAs, sh1-SNAI2 and sh2-SNAI2, which target distinct regions of the SNAI2 transcript. Both shRNAs reduced SNAI2 mRNA levels by > 75% compared to cells expressing a control shRNA (Figure 3A). Cells were then treated with dox and/or E2, followed by RT-qPCR analysis of EMT markers and Matrigel™ invasion assays. Remarkably, SNAI2 knockdown by RNA interference abolished the Runx2-mediated decrease in CDH1 mRNA and increase in VIM mRNA, implicating SNAI2 in Runx2-mediated EMT (Figure 3B). Although SNAI2 silencing did not antagonize the stimulation of S100A4 by Runx2 (Figure 3B), it completely abrogated the ability of Runx2 to increase MCF7 cell invasion through Matrigel™ (Figure 3C).
WNT and TGFβ signaling are required for Runx2-mediated stimulation of SNAI2 expression
EMT is regulated by various signal transduction pathways including WNT and TGFβ [45, 46], and IPA™ analysis of the 530 'common' genes for representation of developmental and growth factor signaling pathways indicated enrichment of WNT/β-catenin pathway-related genes (Figure 4A). Both the WNT and TGFβ pathways have been also implicated in the regulation of SNAI2 expression [47, 48]. Given that SNAI2 was required for Runx2-mediated EMT and invasiveness, we asked whether Runx2-mediated stimulation of SNAI2 expression involved WNT and TGFβ signaling pathways. To this end, MCF7/Rx2dox cells were treated with dox, E2 or both, along with inhibitors of either WNT or TGFβ signaling. WNT signaling was inhibited by using ICG-001, a small molecule that interferes with the recruitment of CAAT-enhancer binding protein (CBP) to β-catenin/LEF/TCF complexes . For TGFβ blockade, we used neutralizing antibodies against TGFBR1. As shown in Figure 4B, both ICG-001 and the anti-TGFBR1 antibody attenuated Runx2-induced SNAI2 expression. Interestingly, attenuation of the Runx2 response with either ICG-001 or the anti-TGFBR1 antibodies was no longer observed in cells co-treated with both dox and E2 (Figure 4B). Co-blockade of the WNT and TGFβ pathways was as effective as blocking each separately (data not shown). Thus, inhibition of either WNT or TGFβ signaling, or the presence of E2, suppressed Runx2-stimulated SNAI2 expression in MCF7 BCa cells. Neither WNT nor TGFβ interference affected Runx2-stimulated SNAI2 expression in the presence of E2. These results suggest that the WNT and TGFβ blockade may have targeted the same pathways by which E2 antagonizes Runx2 in BCa cells.
Negative association of SNAI2 with metastasis-free survival in BCa patients
We next determined the association between SNAI2 expression level and clinical BCa progression by constructing Kaplan-Meier curves for overall or site-specific metastasis-free survival of tumors expressing either high or low levels of SNAI2. Of the 557 patients compiled from the publicly available datasets GSE2034, GSE2603 and GSE12276 [37–39], 288 (52%) relapsed, with 184, 28 and 76 developing bone, brain and lung metastasis, respectively [38, 39]. In cases of metastasis to multiple organs, each of the corresponding sites was included in the analysis. Overall, tumors expressing SNAI2 at high levels were more likely to relapse than those with low SNAI2 expression (Figure 5A-D). Accordingly, tumors which did not metastasize were mostly SNAI2low (Figure 5M). The association between SNAI2 expression and metastasis was particularly strong for bone (P = 4.04E-010) (Figure 5D). These results suggest that high expression of SNAI2 in BCa primes tumors for metastasis in general and bone metastasis in particular.
Anomalous high expression of SNAI2, Runx2 and ERα in bone-seeking BCa tumors
Suppression of SNAI2 expression in BCa cells by E2 (Figure 2A) may be related to the protective effect of ERα against metastasis. Indeed, consistent with previous reports [15, 17–21], tumors with high ESR1 expression were less likely to metastasize in our cohort, particularly during the initial 20 to 40 months after surgery (Figure 5E). Analysis of the association between ESR1 expression and metastasis to specific sites, however, revealed an interesting exception. Whereas high ESR1 expression was negatively correlated with lung and brain metastasis (Figure 5F-G), it positively correlated with bone metastasis (Figure 5H). Thus, unleashing SNAI2 gene expression by loss of ERα is a plausible mechanism promoting metastasis in general; however, high expression of SNAI2 in cells that metastasize to bone occurs despite high ESR1 expression. Furthermore, SNAI2 expression was associated with lung and brain metastasis only for tumors expressing low levels of ESR1 (Figure 5J-K), whereas bone metastasis was associated with SNAI2 for tumors expressing ESR1 at either low or high levels (Figure 5L), suggesting weaker control of SNAI2 by ERα in bone-seeking tumors. Consequently, the combination of low ESR1 and high SNAI2 expression was predominant within lung and brain-seeking tumors (Figure 5O-P), whereas bone-seeking tumors were mostly those expressing both ESR1 and SNAI2 at high levels (Figure 5Q). Consistent with these observations, the negative correlation between ESR1 and SNAI2 expression was weaker in bone compared to lung-seeking tumors (Figure 5R-U; brain-seeking tumors were not analyzed because only eight of them were ER-positive). Finally, we investigated Runx2, SNAI2 and ESR1 gene expression in 58 BCa tumors derived from metastatic sites: 16 from bone, 19 from brain, 18 from lung and 5 from liver . The standardized signal intensities plot showed preferential upregulation of Runx2, SNAI2 and ESR1 in the bone-derived biopsies (Figure 5V). The unusual concomitant high expression of ESR1 and Runx2/SNAI2 may reflect dominant regulation of Runx2/SNAI2 by E2-resistant signaling pathways that promote BCa bone metastasis.
Apart from promoting BCa progression, estrogen signaling has paradoxical anti-metastatic properties [15, 17, 19–21]. Potentially contributing to this, E2 antagonized the transcription factor Runx2, whose role in metastasis is being increasingly recognized [22, 24, 30–32]. Specifically, Runx2 promoted EMT and invasion of BCa cells in vitro, and this was antagonized by E2 (Figure 2). At the center stage of the opposing effects of Runx2 and E2 signaling on EMT and invasion was the transcription factor SNAI2. Consistent with previous reports [30, 48, 49], SNAI2 was stimulated by Runx2 and inhibited by E2 (Figure 1A and 2A); Runx2 no longer enhanced EMT and invasion after SNAI2 knockdown; SNAI2 expression in BCa biopsies positively correlated with a metagene that reports on Runx2 activity and negatively correlated with ERα mRNA (Figure 1B), and expression of SNAI2 was associated with BCa metastasis in general and bone metastasis in particular (Figure 5A-D).
SNAI2 likely plays different roles during various stages of BCa progression, and its ability to transcriptionally repress E-cadherin and induce EMT is well documented [3–5, 49–51]. The role of SNAI2 downstream of Runx2 in BCa is reminiscent of the role of its homologue, SNAI1, in EMT and metastasis of breast, ovarian, colon, lung and squamous cell carcinomas . However, SNAI1 is not stimulated by Runx2 in either MCF7 BCa  or C4-2B PCa cells . Instead, it is SNAI2 that is strongly stimulated by Runx2 in BCa and PCa cells (Figure 2). Furthermore, SNAI2 and not SNAI1 was repressed by E2 in our study (data not shown), and SNAI2 and not SNAI1 exhibited an inverse correlation with ERα and E-cadherin in human BCa tumors . Like Runx2, Twist1 has recently been shown to induce SNAI2 expression, and SNAI2 was essential for Twist1-mediated EMT of human mammary epithelial cells . Thus, depending on context, either SNAI1 or SNAI2 regulate BCa metastasis, with the latter mediating the pro-metastatic activity of Runx2 and Twist1 as well as the anti-metastatic activity of E2.
It remains to be investigated whether Runx2 directly regulates SNAI2 transcription. On the one hand, there are several Runx-binding motifs upstream of the SNAI2 transcription start site (Additional file 3: Supplemental Figure 1A). However, Runx2 ChIP-seq analysis in C4-2B cells, where SNAI2 expression is also stimulated by Runx2 (Figure 1B), did not suggest occupancy of Runx2 at these SNAI2 upstream Runx motifs . Instead, we believe that Runx2 indirectly stimulates SNAI2 expression via modulation of ETS signaling. This hypothesis is based on the observation that Runx2 down-regulated expression of SPDEF (Additional file 3: Supplemental Figure 1B and ), a transcriptional inhibitor belonging to the ETS family, which suppresses tumor progression and metastasis  in part through direct inhibition of SNAI2 expression . Indeed, our Runx2 ChIP-seq analysis in C4-2B/Rx2dox cells  demonstrated strong occupancy of the SPDEF transcription start site by Runx2 (Additional file 3: Supplemental Figure 1C). We are currently investigating how such Runx2 occupancy represses SPDEF transcription, and which of the many ETS sites upstream of the SNAI2 transcription start site (Additional file 3: Supplemental Figure 1A) mediate its repression by SPDEF and/or its activation by other members of the ETS family. Regardless of the precise mechanism by which Runx2 controls SNAI2, expression of the two is tightly correlated in BCa tumors (Figure 1D). In fact, much of the reported association between Runx2 expression and BCa metastasis  is attributable to SNAI2. Indeed, the results of our independent analysis validating the association of Runx2 with metastasis (Additional file 4) are remarkably similar to those demonstrating the association of SNAI2 with metastasis (Figure 5).
The anti-Runx2/SNAI2, anti-metastatic properties of E2 signaling in BCa cells may be relevant for early in vivo metastatic events (Figure 5E), and possibly to targeting BCa metastasis to non-osseous tissues (Figure 5F-G). However, they appear less relevant for bone metastasis. Whereas SNAI2 expression is associated with bone-seeking primary BCa tumors and BCa bone metastasis even more strongly than with lung or brain metastasis, ERα is negatively associated with only non-osseous metastasis (Figure 5). In fact, consistent with previous reports [9, 55], the correlation between ERα expression and bone metastasis was positive, not negative in our cohort (Figure 5H); BCa tumors that metastasized to bone were mostly ER-positive (Figure 5Q) and BCa bone metastases co-expressed high levels of Runx2/SNAI2 and ERα (Figure 5V). The bone-seeking property that ERα bestows on BCa cells may be counteracted during the first two years after surgery by the general anti-metastatic property of estrogens, resulting in a total neutral effect, but at the end it is the former that prevails (Figure 5H). Possibly, activation of bone-seeking pathways in primary BCa cells, such as Src1 [8, 9], results in the stimulation of SNAI2 and other Runx2-regulated genes despite high levels of ERα.
Runx2 stimulates EMT and invasiveness in MCF7 BCa cell cultures, adding to the mounting evidence for its role in metastasis. Our data suggest that E2 attenuates BCa metastasis by antagonizing Runx2. Indeed, EMT and invasiveness of MCF7 cells were severely compromised by E2 specifically after induction of Runx2. The linkage between the pro- and anti-metastatic properties of Runx2 and E2, respectively, is attributable in part to SNAI2, which is upregulated by Runx2 and downregulated by E2. Furthermore, we observe a strong positive association between Runx2 and SNAI2 expression, a negative association between ESR1 and SNAI2, and a positive correlation between SNAI2 expression and metastasis in BCa tumors. These observations are consistent with the hypothesis that unopposed stimulation of Runx2 target genes such as SNAI2 in ER-negative tumors contributes to their aggressive metastatic phenotype. However, the relationship between ER expression and metastasis is more complicated in bone-seeking tumors. Here, the negative correlation between SNAI2 and ER is weak, and the anti-metastatic property of ER is likely masked by an opposite property, which ultimately prevails. Mimicking ER signaling to specifically antagonize Runx2 and SNAI2 offers a research avenue towards the development of novel therapeutic approaches for the management of BCa patients who fail first line therapy.
Database for Annotation: Visualization and Integrated Discovery
estrogen receptor alpha
Ingenuity Pathway Analysis
Robust Multiarray Averaging
snail homolog 2
SAM pointed domain-containing ets transcription factor.
Thiery JP, Acloque H, Huang RY, Nieto MA: Epithelial-mesenchymal transitions in development and disease. Cell. 2009, 139: 871-890. 10.1016/j.cell.2009.11.007.
Nguyen DX, Bos PD, Massague J: Metastasis: from dissemination to organ-specific colonization. Nat Rev Cancer. 2009, 9: 274-284. 10.1038/nrc2622.
Foubert E, De Craene B, Berx G: Key signalling nodes in mammary gland development and cancer. The Snail1-Twist1 conspiracy in malignant breast cancer progression. Breast Cancer Res. 2010, 12: 206-10.1186/bcr2585.
Peinado H, Olmeda D, Cano A: Snail, Zeb and bHLH factors in tumour progression: an alliance against the epithelial phenotype?. Nat Rev Cancer. 2007, 7: 415-428. 10.1038/nrc2131.
Nieto MA: The snail superfamily of zinc-finger transcription factors. Nat Rev Mol Cell Biol. 2002, 3: 155-166.
Roodman GD: Mechanisms of bone metastasis. N Eng J Med. 2004, 350: 1655-1664. 10.1056/NEJMra030831.
Jemal A, Siegel R, Ward E, Hao Y, Xu J, Thun MJ: Cancer statistics, 2009. CA CancerJ Clin. 2009, 59: 225-249. 10.3322/caac.20006.
Kang Y, Siegel PM, Shu W, Drobnjak M, Kakonen SM, Cordon-Cardo C, Guise TA, Massague J: A multigenic program mediating breast cancer metastasis to bone. Cancer Cell. 2003, 3: 537-549. 10.1016/S1535-6108(03)00132-6.
Zhang XH, Wang Q, Gerald W, Hudis CA, Norton L, Smid M, Foekens JA, Massague J: Latent bone metastasis in breast cancer tied to Src-dependent survival signals. Cancer Cell. 2009, 16: 67-78. 10.1016/j.ccr.2009.05.017.
Henderson IC: Risk factors for breast cancer development. Cancer. 1993, 71: 2127-2140. 10.1002/1097-0142(19930315)71:6+<2127::AID-CNCR2820711602>3.0.CO;2-2.
Ma H, Bernstein L, Pike MC, Ursin G: Reproductive factors and breast cancer risk according to joint estrogen and progesterone receptor status: a meta-analysis of epidemiological studies. Breast Cancer Res. 2006, 8: R43-10.1186/bcr1525.
Yager JD, Davidson NE: Estrogen carcinogenesis in breast cancer. N Engl J Med. 2006, 354: 270-282. 10.1056/NEJMra050776.
Chlebowski RT, Hendrix SL, Langer RD, Stefanick ML, Gass M, Lane D, Rodabough RJ, Gilligan MA, Cyr MG, Thomson CA, Khandekar J, Petrovitch H, McTiernan A, WHI Investigators: Influence of estrogen plus progestin on breast cancer and mammography in healthy postmenopausal women: the Women's Health Initiative Randomized Trial. JAMA. 2003, 289: 3243-3253. 10.1001/jama.289.24.3243.
Pike MC, Krailo MD, Henderson BE, Casagrande JT, Hoel DG: 'Hormonal' risk factors, 'breast tissue age' and the age-incidence of breast cancer. Nature. 1983, 303: 767-770. 10.1038/303767a0.
Sotiriou C, Piccart MJ: Taking gene-expression profiling to the clinic: when will molecular signatures become relevant to patient care?. Nat Rev Cancer. 2007, 7: 545-553. 10.1038/nrc2173.
Anderson GL, Limacher M, Assaf AR, Bassford T, Beresford SA, Black H, Bonds D, Brunner R, Brzyski R, Caan B, Chlebowski R, Curb D, Gass M, Hays J, Heiss G, Hendrix S, Howard BV, Hsia J, Hubbell A, Jackson R, Johnson KC, Judd H, Kotchen JM, Kuller L, LaCroix AZ, Lane D, Langer RD, Lasser N, Lewis CE, Manson J, et al: Effects of conjugated equine estrogen in postmenopausal women with hysterectomy: the Women's Health Initiative randomized controlled trial. JAMA. 2004, 291: 1701-1712.
Ellis MJ, Gao F, Dehdashti F, Jeffe DB, Marcom PK, Carey LA, Dickler MN, Silverman P, Fleming GF, Kommareddy A, Jamalabadi-Majidi S, Crowder R, Siegel BA: Lower-dose vs high-dose oral estradiol therapy of hormone receptor-positive, aromatase inhibitor-resistant advanced breast cancer: a phase 2 randomized study. JAMA. 2009, 302: 774-780. 10.1001/jama.2009.1204.
Maynadier M, Nirde P, Ramirez JM, Cathiard AM, Platet N, Chambon M, Garcia M: Role of estrogens and their receptors in adhesion and invasiveness of breast cancer cells. Adv Exp Med Biol. 2008, 617: 485-491. 10.1007/978-0-387-69080-3_48.
Platet N, Cunat S, Chalbos D, Rochefort H, Garcia M: Unliganded and liganded estrogen receptors protect against cancer invasion via different mechanisms. Mol Endocrinol. 2000, 14: 999-1009. 10.1210/me.14.7.999.
Maynadier M, Ramirez JM, Cathiard AM, Platet N, Gras D, Gleizes M, Sheikh MS, Nirde P, Garcia M: Unliganded estrogen receptor alpha inhibits breast cancer cell growth through interaction with a cyclin-dependent kinase inhibitor (p21(WAF1)). FASEB J. 2008, 22: 671-681.
Platet N, Cathiard AM, Gleizes M, Garcia M: Estrogens and their receptors in breast cancer progression: a dual role in cancer proliferation and invasion. Crit Rev Oncol Hematol. 2004, 51: 55-67. 10.1016/j.critrevonc.2004.02.001.
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.
Frenkel B, Hong A, Baniwal SK, Coetzee GA, Ohlsson C, Khalid O, Gabet Y: Regulation of adult bone turnover by sex steroids. J Cell Physiol. 2010, 224: 305-310. 10.1002/jcp.22159.
Blyth K, Vaillant F, Jenkins A, McDonald L, Pringle MA, Huser C, Stein T, Neil J, Cameron ER: Runx2 in normal tissues and cancer cells: A developing story. Blood Cells Mol Dis. 2010, 45: 117-123. 10.1016/j.bcmd.2010.05.007.
Otto F, Thornell AP, Crompton T, Denzel A, Gilmour KC, Rosewell IR, Stamp GW, Beddington RS, Mundlos S, Olsen BR, Selby PB, Owen MJ: Cbfa1, a candidate gene for cleidocranial dysplasia syndrome, is essential for osteoblast differentiation and bone development. Cell. 1997, 89: 765-771. 10.1016/S0092-8674(00)80259-7.
Komori T, Yagi H, Nomura S, Yamaguchi A, Sasaki K, Deguchi K, Shimizu Y, Bronson RT, Gao YH, Inada M, Sato M, Okamoto R, Kitamura Y, Yoshiki S, Kishimoto T: Targeted disruption of Cbfa1 results in a complete lack of bone formation owing to maturational arrest of osteoblasts. Cell. 1997, 89: 755-764. 10.1016/S0092-8674(00)80258-5.
Javed A, Barnes GL, Pratap J, Antkowiak T, Gerstenfeld LC, van Wijnen AJ, Stein JL, Lian JB, Stein GS: Impaired intranuclear trafficking of Runx2 (AML3/CBFA1) transcription factors in breast cancer cells inhibits osteolysis in vivo. Proc Natl Acad Sci USA. 2005, 102: 1454-1459. 10.1073/pnas.0409121102.
Pratap J, Imbalzano KM, Underwood JM, Cohet N, Gokul K, Akech J, van Wijnen AJ, Stein JL, Imbalzano AN, Nickerson JA, Lian JB, Stein GS: Ectopic runx2 expression in mammary epithelial cells disrupts formation of normal acini structure: implications for breast cancer progression. Cancer Res. 2009, 69: 6807-6814. 10.1158/0008-5472.CAN-09-1471.
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.
Baniwal SK, Khalid O, Gabet Y, Shah RR, Purcell DJ, Mav D, Kohn-Gabet AE, Shi Y, Coetzee GA, Frenkel B: Runx2 transcriptome of prostate cancer cells: insights into invasiveness and bone metastasis. Mol Cancer. 2010, 9: 258-10.1186/1476-4598-9-258.
Akech J, Wixted JJ, Bedard K, van der Deen M, Hussain S, Guise TA, van Wijnen AJ, Stein JL, Languino LR, Altieri DC, Pratap J, Keller E, Stein GS, Lian JB: Runx2 association with progression of prostate cancer in patients: mechanisms mediating bone osteolysis and osteoblastic metastatic lesions. Oncogene. 2010, 29: 811-821. 10.1038/onc.2009.389.
Pratap J, Lian JB, Stein GS: Metastatic bone disease: role of transcription factors and future targets. Bone. 2011, 48: 30-36. 10.1016/j.bone.2010.05.035.
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.
Chimge NO, Baniwal SK, Luo J, Coetzee S, Khalid O, Berman B, Tripathy D, Ellis MJ, Frenkel B: Opposing effects of Runx2 and estradiol on breast cancer cell proliferation: in vitro identification of reciprocally-regulated gene signature related to clinical letrozole responsiveness. Clin Cancer Res.
Emami KH, Nguyen C, Ma H, Kim DH, Jeong KW, Eguchi M, Moon RT, Teo JL, Kim HY, Moon SH, Ha JR, Kahn M: A small molecule inhibitor of beta-catenin/CREB-binding protein transcription [corrected]. Proc Nat Acad Sci USA. 2004, 101: 12682-12687. 10.1073/pnas.0404875101.
Baniwal SK, Little GH, Chimge NO, Frenkel B: Runx2 controls a feed-forward loop between androgen and prolactin-induced protein (PIP) in stimulating T47D cell proliferation. J Cell Physiol.
Bos PD, Zhang XH, Nadal C, Shu W, Gomis RR, Nguyen DX, Minn AJ, van de Vijver MJ, Gerald WL, Foekens JA, Massagué J: Genes that mediate breast cancer metastasis to the brain. Nature. 2009, 459: 1005-1009. 10.1038/nature08021.
Minn AJ, Gupta GP, Siegel PM, Bos PD, Shu W, Giri DD, Viale A, Olshen AB, Gerald WL, Massague J: Genes that mediate breast cancer metastasis to lung. Nature. 2005, 436: 518-524. 10.1038/nature03799.
Wang Y, Klijn JG, Zhang Y, Sieuwerts AM, Look MP, Yang F, Talantov D, Timmermans M, Meijer-van Gelder ME, Yu J, Jatkoe T, Berns EM, Atkins D, Foekens JA: Gene-expression profiles to predict distant metastasis of lymph-node-negative primary breast cancer. Lancet. 2005, 365: 671-679.
Huang da W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4: 44-57.
Williams SA, Harata-Lee Y, Comerford I, Anderson RL, Smyth MJ, McColl SR: Multiple functions of CXCL12 in a syngeneic model of breast cancer. Mol Cancer. 2010, 9: 250-10.1186/1476-4598-9-250.
Ghavami S, Chitayat S, Hashemi M, Eshraghi M, Chazin WJ, Halayko AJ, Kerkhoff C: S100A8/A9: a Janus-faced molecule in cancer therapy and tumorgenesis. Eur J Pharmacol. 2009, 625: 73-83. 10.1016/j.ejphar.2009.08.044.
Ma X, Yang Y, Wang Y, An G, Lv G: Small interfering RNA-directed knockdown of S100A4 decreases proliferation and invasiveness of osteosarcoma cells. Cancer Lett. 299: 171-181.
Kwak JM, Lee HJ, Kim SH, Kim HK, Mok YJ, Park YT, Choi JS, Moon HY: Expression of protein S100A4 is a predictor of recurrence in colorectal cancer. World J Gastroenterol. 2010, 16: 3897-3904. 10.3748/wjg.v16.i31.3897.
Miettinen PJ, Ebner R, Lopez AR, Derynck R: TGF-beta induced transdifferentiation of mammary epithelial cells to mesenchymal cells: involvement of type I receptors. J Cell Biol. 1994, 127: 2021-2036. 10.1083/jcb.127.6.2021.
Moustakas A, Heldin CH: Signaling networks guiding epithelial-mesenchymal transitions during embryogenesis and cancer progression. Cancer Sci. 2007, 98: 1512-1520. 10.1111/j.1349-7006.2007.00550.x.
Vincent T, Neve EP, Johnson JR, Kukalev A, Rojo F, Albanell J, Pietras K, Virtanen I, Philipson L, Leopold PL, Crystal RG, de Herreros AG, Moustakas A, Pettersson RF, Fuxe J: A SNAIL1-SMAD3/4 transcriptional repressor complex promotes TGF-beta mediated epithelial-mesenchymal transition. Nat Cell Biol. 2009, 11: 943-950. 10.1038/ncb1905.
Lambertini E, Franceschetti T, Torreggiani E, Penolazzi L, Pastore A, Pelucchi S, Gambari R, Piva R: SLUG: a new target of lymphoid enhancer factor-1 in human osteoblasts. BMC Mol Biol. 2010, 11: 13-10.1186/1471-2199-11-13.
Ye Y, Xiao Y, Wang W, Yearsley K, Gao JX, Shetuni B, Barsky SH: ERalpha signaling through slug regulates E-cadherin and EMT. Oncogene. 2010, 29: 1451-1462. 10.1038/onc.2009.433.
Casas E, Kim J, Bendesky A, Ohno-Machado L, Wolfe CJ, Yang J: Snail2 is an essential mediator of Twist1-induced epithelial mesenchymal transition and metastasis. Cancer Res. 2011, 71: 245-254. 10.1158/0008-5472.CAN-10-2330.
Bolos V, Peinado H, Perez-Moreno MA, Fraga MF, Esteller M, Cano A: The transcription factor Slug represses E-cadherin expression and induces epithelial to mesenchymal transitions: a comparison with Snail and E47 repressors. J Cell Sci. 2003, 116: 499-511. 10.1242/jcs.00224.
Little GH, Noushmehr H, Baniwal SK, Berman BP, Coetzee GA, Frenkel B: Genome-wide Runx2 occupancy in prostate cancer cells suggests a role in regulating secretion. Nucleic Acids Res.
Steffan JJ, Koul HK: Prostate derived ETS factor (PDEF): a putative tumor metastasis suppressor. Cancer Lett. 2011, 310: 109-117. 10.1016/j.canlet.2011.06.011.
Findlay VJ, Turner DP, Yordy JS, McCarragher B, Shriver MR, Szalai G, Watson PM, Larue AC, Moussa O, Watson DK: Prostate-derived ETS factor regulates epithelial-to-mesenchymal transition through both SLUG-dependent and independent mechanisms. Genes Cancer. 2011, 2: 120-129. 10.1177/1947601911410424.
Hess KR, Pusztai L, Buzdar AU, Hortobagyi GN: Estrogen receptors and distinct patterns of breast cancer relapse. Breast Cancer Res Treat. 2003, 78: 105-118. 10.1023/A:1022166517963.
We thank USC's Dr. Susan Groshen at the Norris Cancer Center Biostatistics Services Core for guidance with data analysis, Ms. Meng Li at the Bioinformatics Service Program, Norris Medical Library for helpful discussions, and Ms. Yunfan Shi for excellent technical support. This work was funded by NIH grants RO1 DK071122 and DK071122S1 to BF, who holds the J. Harold and Edna L. LaBriola Chair in Genetic Orthopedic Research at USC. ZB is funded by NIH grant RO1 HL089445 and holds the Ralph Edgington Chair in Medicine, DT holds the Priscilla and Art Ulene Chair in Women's Cancer.
The authors declare that they have no competing interests.
NOC, SKB and GHL performed the experiments. NOC and YBC performed the bioinformatics analyses. NOC, MK, ZB, DT and BF designed the study and wrote the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Oligonucleotide sequences used in this study. List of primer pair sequences for qRT-PCR and shRNA sequences. (XLS 22 KB)
Additional file 2: The complete list of the 'common' genes. The complete list of the 'common' genes, independently regulated by Runx2 and E2 in MCF7/Rx2dox cells are depicted along with respective fold changes and p values. (XLS 88 KB)
Additional file 3: Potential Mechanisms Linking Runx2 to SNAI2. Based on the observation that Runx2 down-regulated expression of SPDEF and strong occupancy of the SPDEF transcription start site by Runx2 we hypothesize that Runx2 indirectly stimulates SNAI2 expression in part via modulation of ETS signaling. (PPT 145 KB)
Additional file 4: Association of Runx2 activity with metastasis. Runx2 activity was associated with both osseous and non-osseous metastases. Bone-seeking BCa tumors exhibit an unusual combination of high RUNX2 activity and high ERα levels. (PPTX 446 KB)
About this article
Cite this article
Chimge, NO., Baniwal, S.K., Little, G.H. et al. Regulation of breast cancer metastasis by Runx2 and estrogen signaling: the role of SNAI2. Breast Cancer Res 13, R127 (2011). https://doi.org/10.1186/bcr3073
- Bone Metastasis
- Estrogen Signaling
- Cleidocranial Dysplasia
- SNAI2 Expression
- Runx2 Target Gene