Regulation of breast cancer metastasis by Runx2 and estrogen signaling: the role of SNAI2
© Chimge et al.; licensee BioMed Central Ltd. 2011
Received: 6 July 2011
Accepted: 9 December 2011
Published: 9 December 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
Functional Annotation Clustering of Genes Independently Regulated by Runx2 and E2 in MCF7/Rx2dox Cells
Annotation Cluster 1
Enrichment Score: 4.50
GO:0048514~blood vessel morphogenesis
GO:0001568~blood vessel development
Annotation Cluster 2
Enrichment Score: 4.48
GO:0051674~localization of cell
Genes Independently Regulated by Runx2 and E2 in MCF7/Rx2dox Cells
Functional link between Runx2, E2 and SNAI2 in BCa
E2 antagonizes Runx2-induced EMT and invasiveness in vitro
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
WNT and TGFβ signaling are required for Runx2-mediated stimulation of SNAI2 expression
Negative association of SNAI2 with metastasis-free survival in BCa patients
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.
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.
- 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.PubMedView ArticleGoogle Scholar
- Nguyen DX, Bos PD, Massague J: Metastasis: from dissemination to organ-specific colonization. Nat Rev Cancer. 2009, 9: 274-284. 10.1038/nrc2622.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Nieto MA: The snail superfamily of zinc-finger transcription factors. Nat Rev Mol Cell Biol. 2002, 3: 155-166.PubMedView ArticleGoogle Scholar
- Roodman GD: Mechanisms of bone metastasis. N Eng J Med. 2004, 350: 1655-1664. 10.1056/NEJMra030831.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Yager JD, Davidson NE: Estrogen carcinogenesis in breast cancer. N Engl J Med. 2006, 354: 270-282. 10.1056/NEJMra050776.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.