Research article | Open | Published:
Homeobox transcription factor muscle segment homeobox 2 (Msx2) correlates with good prognosis in breast cancer patients and induces apoptosis in vitro
Breast Cancer Researchvolume 12, Article number: R59 (2010)
The homeobox-containing transcription factor muscle segment homeobox 2 (Msx2) plays an important role in mammary gland development. However, the clinical implications of Msx2 expression in breast cancer are unclear. The aims of this study were to investigate the potential clinical value of Msx2 as a breast cancer biomarker and to clarify its functional role in vitro.
Msx2 gene expression was first examined in a well-validated breast cancer transcriptomic dataset of 295 patients. Msx2 protein expression was then evaluated by immunohistochemistry in a tissue microarray (TMA) containing 281 invasive breast tumours. Finally, to assess the functional role of Msx2 in vitro, Msx2 was ectopically expressed in a highly invasive breast tumour cell line (MDA-MB-231) and an immortalised breast cell line (MCF10a), and these cell lines were examined for changes in growth rate, cell death and cell signalling.
Examination of Msx2 mRNA expression in a breast cancer transcriptomic dataset demonstrated that increased levels of Msx2 were associated with good prognosis (P = 0.011). Evaluation of Msx2 protein expression on a TMA revealed that Msx2 was detectable in both tumour cell nuclei and cytoplasm. Cytoplasmic Msx2 expression was associated with low grade tumours (P = 0.012) and Ki67 negativity (P = 0.018). Nuclear Msx2 correlated with low-grade tumours (P = 0.015), estrogen receptor positivity (P = 0.038), low Ki67 (P = 0.005) and high cyclin D1 expression (P = 0.037). Increased cytoplasmic Msx2 expression was associated with a prolonged breast cancer-specific survival (P = 0.049), recurrence-free survival (P = 0.029) and overall survival (P = 0.019). Ectopic expression of Msx2 in breast cell lines resulted in radically decreased cell viability mediated by induction of cell death via apoptosis. Further analysis of Msx2-expressing cells revealed increased levels of p21 and phosphorylated extracellular signal-regulated kinase (ERK) and decreased levels of Survivin and the 'split ends' (SPEN) protein family member RBM15.
We conclude that increased Msx2 expression results in improved outcome for breast cancer patients, possibly by increasing the likelihood of tumour cell death by apoptosis.
Homeobox genes are important during embryonic development, where they function to control cell fate and positioning, thereby regulating the morphological development of several organs, including skeletal structures, the heart, teeth, eyes and mammary glands [1–6]. The homeobox family of transcription factors were originally isolated from Drosophila and contain a common 61-amino acid domain, known as the homeodomain, which can directly bind DNA and regulate gene transcription . Mutations in the muscle segment homeobox 2 (Msx2) homeodomain which cause loss or gain of Msx2 DNA binding activity can both result in cranial defects [8, 9]. Msx2 function can also be affected by subcellular localisation and protein-protein interactions .
Within the mammary gland, homeobox genes are thought to be involved in assimilating systemic signals into the precise local interactions required for correct morphogenesis [7, 10]. During puberty, the core epithelial structure of the mammary tree is established through the invasion of the terminal end buds from the nipple into the surrounding fat pad . Apoptosis plays an important role in this process, and the level of apoptosis in the mouse pubertal mammary gland is higher than at any other stage of development . In the mouse, Msx2 is expressed during pubertal development and early pregnancy, downregulated during late pregnancy and lactation, and reexpressed during involution [3, 13]. Msx2 expression is stimulated by both estradiol and progesterone, and the role of progesterone in promoting branching morphogenesis in the mouse mammary gland is thought to be mediated partially through Msx2 . However, one of the few studies on Msx2 expression in human breast tissues reported a complicated regulation of Msx2 by steroid hormones: Msx2 could be either increased or decreased by steroid hormone treatment, depending on the estrogen/progesterone receptor (ER/PR) status and whether the breast tissue sample was normal or malignant .
Given its widespread regulatory role in growth and development, it would not be surprising to find that Msx2 is involved in tumourigenic processes. Indeed, studies in various cell line models have suggested that increased expression of Msx2 can induce neoplastic transformation and epithelial-to-mesenchymal transition (EMT) [16, 17]. However, Barnes et al. found that Msx2 repressed the activity of the prometastatic bone sialoprotein (BSP) in breast cancer cell lines , casting doubt on a proinvasive role for Msx2, as BSP is associated with poor outcome in breast cancer . A role has also been proposed for Msx2 in the regulation of cell death. A recent study of Sonic hedgehog-1 (Shh)-knockout mice  found that expression of Msx1 and Msx2 was increased, and this was partially responsible for the massive apoptosis and severe limb defects seen in the Shh-null phenotype. Furthermore, ectopic expression of Msx2 can induce apoptosis in pluripotent murine embryonic carcinoma cells , and cranial neural crest-derived cells . This effect on cell death can also be induced by expression of the bone morphogenetic protein (BMP) family: Msx2 and p21 are induced following BMP4 treatment, and both molecules are necessary for BMP4-mediated cell death to occur in several cell types [21, 22].
Large-scale clinical studies of Msx2 expression in cancer are few in number. A study of 32 pancreatic adenocarcinomas demonstrated an association between Msx2 expression and high tumour grade and vascular invasion . The only study of Msx2 in breast tumours to date involved four invasive ductal carcinomas and found that Msx2 expression was increased in infiltrating compared to noninfiltrating cells . The prognostic influence of Msx2 expression has not been investigated in a large cohort of breast cancer patients. The aim of this study was to investigate Msx2 expression as a prognostic biomarker in breast cancer. We also sought to clarify the in vitro functional role of Msx2. Our findings indicate that increased expression of both Msx2 mRNA and protein are associated with improved outcome in breast cancer. We also found that ectopic expression of Msx2 in breast cell lines leads to induction of apoptosis and radically decreased cell viability. These data suggest that increased Msx2 results in improved outcome for breast cancer patients, possibly by increasing the likelihood of tumour cell death by apoptosis.
Materials and methods
Statistical analysis of DNA microarray data
Relevant gene expression and clinical data relating to 295 patients with breast cancer  were downloaded from Rosetta Inpharmatics Inc. . The log ratios of gene expression values were used without modification and classified using a previously published method . Tumour samples were classified by first separating into quartiles according to mRNA expression. Adjacent groups with significant overlap of Kaplan-Meier survival curves were combined, and the survival curves of the resulting two groups were compared using the log-rank test. The χ2 test and Fisher's exact test were used for relating mRNA levels to clinicopathological variables.
The tissue microarray (TMA) used in this study was derived from a reference cohort of 512 consecutive invasive breast cancer cases diagnosed at Malmö University Hospital, Sweden, between 1988 and 1992, and has been previously described [26, 27]. From the original cohort of 512 patients, samples were available from 281 patients for analysis of Msx2 protein expression (this reduced number was primarily due to core loss). Patient and tumour characteristics of the available and missing cohorts are outlined in Additional file 1. The study has been approved by the Ethics Committee at Lund University.
The TMA was constructed as previously described . Briefly, two 0.6-mm tissue cores were extracted from each donor block using an automated tissue arrayer (MTA-27; Beecher Inc., Sun Prairie, WI, USA) and placed into a recipient block. To confirm Western blot analysis results, a cell pellet array (CPA) was constructed as previously described  using formalin-fixed, paraffin-embedded (FFPE) breast cell lines.
All cell lines (MCF-7, T47 D, BT474, ZR75-1, MDA-MB-231, and SKBR3) were purchased from the European Collection of Cell Cultures (Wiltshire, UK), except for the MCF10a cell line, which was a gift from Dr. Geert Berx, University of Ghent, Belgium, and the Hs578t isogenic cell line series (Hs578t P and i8) , which was a gift from Dr. Susan McDonnell, University College Dublin. All cell lines were maintained as previously described .
Western blot analysis
For analysis of cell lines, protein was extracted as previously described . Lysates were separated by reducing sodium dodecyl sulphate-polymerase gel electrophoresis (SDS-PAGE), transferred to polyvinylidene fluoride (PVDF), and immunoblotted using antibodies against Msx2 (clone 2E12, 1:1,000; Abcam, Cambridge, UK), extracellular signal-regulated kinase (ERK) (1:1,000; Santa Cruz Biotechnology, Santa Cruz, CA, USA), p-ERK (1:1,000; Santa Cruz Biotechnology), p21 (1:1,000; BD Biosciences, Franklin Lakes, NJ, USA), Survivin (1:100; Santa Cruz Biotechnology), cyclin D1 (1:1,000; Santa Cruz Biotechnology), DLX5 (clone 3B11; 1:1,000; Abnova, Taipei, Taiwan), Twist (1:1,000; Santa Cruz Biotechnology), and Smad4 (1:1,000; Cell Signaling, Danvers, MA, USA). Membranes were stripped and reprobed with anti-β-actin (1:5,000 dilution; Abcam) as a loading control.
TMA and CPA sections (4 μm) were rehydrated in descending gradient alcohols. Heat-mediated antigen retrieval was performed using 10 mM sodium citrate buffer (pH 6) in a PT module (LabVision; Thermo-Fisher Scientific, Fremont, CA, USA) for 15 min at 95 °C, followed by immunohistochemistry (IHC) in a Lab Vision Autostainer 360 (LabVision, Fremont, CA, USA) for Msx2 (clone 2E12; 1:25; Abcam) or in the Ventana Benchmark system (Ventana Medical Systems Inc., Oro Valley, AZ) using prediluted antibodies to estrogen receptor (clone 6F11; Ventana Medical Systems Inc.), progesterone receptor (clone 16; Ventana Medical Systems Inc.) and Her2 (Pathway, Clone CB-11, Ventana Medical Systems) or in the Dako Techmate 500 system (Dako, Glostrup, Denmark) for Ki-67 (M7240, 1:200; Dako) and cyclin D1 (clone DSC-6, 1:100; Dako). A mouse IgG2a isotype control (Abcam) was used to evaluate Msx2 antibody specificity.
Slides were scanned at ×20 magnification using a ScanScopeXT slide scanner (Aperio Technologies, Vista, CA, USA). Tumour samples were evaluated by two independent observers (RH and FL) and scored for Msx2 staining intensity in both the cytoplasmic and nuclear compartments on a scale of 0 to 3, where 0 is negative, 1 is weakly positive, 2 is medium positive and 3 is strongly positive. Nuclear staining for Msx2 was found to be relatively homogeneous within tumour samples, and thus the percentage of nuclear staining was not included in the scoring system. The mean value of both scores for each patient was used for statistical analysis. ER, PR, HER2 and Ki-67 were assessed as previously described .
To control for the subjectivity inherent in the manual scoring process, we utilised a co-localisation image analysis algorithm (Aperio Technologies) to further examine Msx2 localisation and expression. This algorithm classifies each pixel as either blue (negative nuclear), blue and brown (positive nuclear), or brown (positive cytoplasmic). The average positive pixel intensities of 3,3'-diaminobenzidine (DAB) staining in the cytoplasm and nuclei were used for statistical analysis, divided at the 50th percentile.
Statistical analysis of TMA data
The χ2 test and Fisher's exact test were used to evaluate associations between Msx2 expression and clinicopathological characteristics. Kaplan-Meier plots were used for survival analysis, and the log-rank test was used to compare curves separated according to Msx2 expression. Cox proportional hazards regression was used to estimate hazard ratios (HR). All calculations were carried out using SPSS version 15 (SPSS Inc., Chicago, IL, USA).
Lentiviral overexpression of Msx2
Human embryonic kidney (HEK)-293t cells, at ~60-70% confluence, were transfected using a calcium phosphate transfection method with an LLCIEP vector (Trono Laboratory, Lausanne, Switzerland) with no insert (EV) or containing full-length Msx2 with a V5 tag (MSX2), together with packaging and envelope plasmids (PSPAX2 and PMD2G; Trono Laboratory). The media were refreshed after 6-8 hr and, after a further 48 hr, the viral supernatant was removed and filtered through a 0.45-μm low-protein binding filter (Millipore, Billerica, MA, USA). This was added to ~50% confluent MDA-MB-231 or MCF10a cells at a 1:3 dilution with fresh media, along with 8 μg/ml polybrene. MCF10a cells were centrifuged at 300 × g for 1 hr at 30°C following addition of the viral supernatant. Media were refreshed after 24 hr. The transduction efficiency was estimated at >95% using immunofluorescence microscopy for Msx2 expression.
Cells growing on chamber slides (Nunc, Roskilde, Denmark) were fixed with 3.7% paraformaldehyde for 20 min, permeabilised with 0.5% Triton-X 100 for 10 min, and labelled immunofluorescently using an anti-Msx2 mouse monoclonal antibody (Abcam), followed by incubation with a rhodamine-coupled secondary antibody (Abcam) and 4',6-diamidino-2-phenylindole (DAPI) (Sigma-Aldrich, Dorset, UK). Images were obtained using a Zeiss LSM-510-Meta confocal microscope.
MTT proliferation assay
Cellular proliferation was measured using an MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide) colorimetric assay over a period of 5 days. Cells were seeded at a density of 2,000 cells/well in 96-well plates, with one plate measured every day for 5 days. Cells were incubated with MTT reagent (5 mg/ml) for 4 hr, solubilised in dimethyl sulfoxide, and the absorbance at 570 nm was measured.
Colony formation assay
Cells were seeded at a density of 500 cells/well in a six-well plate. Plates were incubated for 3 weeks with media changed weekly. Cells were fixed for 15 min in 10% neutral-buffered formalin and stained with 0.25% crystal violet solution. Colonies were counted by eye, with three wells counted per cell line in each replicate experiment.
The Caspase-Glo 3/7 kit (Promega, Madison, WI, USA) was used to determine levels of apoptosis in cell lines. Cells were seeded in two 96-well plates at 5,000 cells/well and allowed to grow for 2 days. Media were removed from one plate, replaced with Caspase-Glo working solution in a 1:1 dilution to normal media, left at room temperature for 1 hr, and read using a luminescent plate reader (GloMax Multi Detection System, Promega). The cells on the second plate were trypsinized and counted to normalize the luminescent signal to cell number.
Cell cycle analysis
Cells at ~60-70% confluence were harvested and washed in ice-cold phosphate-buffered saline (PBS) before fixation for 1 hr in ice-cold 70% ethanol in PBS. Cells were collected by centrifugation, washed in PBS and recentrifuged before resuspension in 1 ml of ice-cold PBS containing 0.4 mg/ml propidium iodide (PI). Cells were left for 30 min on ice prior to flow cytometry analysis on a Coulter Epic XL flow cytometer (Beckman Coulter Inc., Brea, CA, USA), to measure the PI uptake and thus the DNA content of the cells.
Msx2 mRNA expression is associated with good prognosis in a DNA microarray dataset
Expression levels of Msx2 mRNA were analysed in a publicly available DNA microarray dataset derived from 295 primary invasive breast tumours . Using a previously described method , 74 tumours were classified as expressing low levels of Msx2 mRNA, and 221 tumours were classified as expressing high levels of Msx2 mRNA, based on the log ratios of gene expression values. Survival analysis of these tumours revealed that increased expression of Msx2 mRNA was associated with good prognosis (P = 0.011) (Figure 1). Univariate Cox regression analysis of Msx2 mRNA expression as a continuous variable supported this association with favourable outcome (HR 0.47; 95% CI 0.24-0.95; P = 0.035). Associations between Msx2 mRNA expression and commonly used clinical biomarkers were also examined (Table 1). Increased Msx2 mRNA expression correlated with ER-positive (P < 0.001) and low-grade (P = 0.003) tumours; moreover, it was associated with good prognosis according to the 70-gene signature defined by van't Veer et al.  (P < 0.001). In relation to breast cancer molecular subtypes described by Sorlie et al. , Msx2 mRNA expression was particularly low in the basal subtype and high in the luminal B and Her2 subtypes of breast cancer (P = 0.001).
Msx2 protein expression in primary breast tumours
The specificity of the Msx2 antibody was validated first by Western blot analysis on a panel of breast cancer cell lines (Figure 2a) and second by IHC in FFPE breast cancer cell lines (Figure 2b). High Msx2 expression was detected by Western blot analysis in MCF7, T47 D, SKBR3 and ZR75-1 cell lines, with low expression seen in BT474 cells and no detectable expression in MDA-MB231, Hs578t, Hs578t-i8 and MCF10a cell lines. This expression pattern appears for the most part to correlate with ER status, apart from the high expression of Msx2 seen in the ER-negative cell line SK-BR3. MCF7 and MCF10a cell lines were used as positive and negative controls, respectively, for all further IHC procedures. Msx2 protein expression was assessed using IHC in a breast cancer TMA constructed from a cohort of 512 patients. Owing to core loss, it was possible to evaluate Msx2 protein expression in 281 tumours (55%) of the 512 tumours represented on the TMA. To evaluate our study for any potential selection bias, baseline clinicopathological characteristics from both the evaluated or 'Msx2 known' cohort (n = 281) and the unevaluated or 'Msx2 unknown' cohort (n = 231) are presented in Additional file 1. No difference was seen in the expression of any variable between either cohort.
Tumour-specific nuclear and cytoplasmic expression of Msx2 protein was evident, with levels of Msx2 in these subcellular compartments being scored separately (Figure 2c). Staining in each compartment was scored on a scale from 0 to 3 on the basis of staining intensity. To divide the tumours into two relatively equal-sized groups, high cytoplasmic Msx2 was defined as a staining intensity ≥2, and high nuclear Msx2 was defined as a staining intensity ≥1. On the basis of this analysis, 168 tumours (59.8%) were classified as expressing low levels of cytoplasmic Msx2, and 113 tumours (40.2%) were classified as expressing high cytoplasmic Msx2. Similarly, 122 tumours (43.4%) were classified as expressing low nuclear Msx2, and 159 tumours (56.6%) were classified as expressing high nuclear Msx2.
Associations between nuclear and cytoplasmic Msx2 protein expression and a number of well-established clinicopathological variables were then investigated (Table 2). Cytoplasmic Msx2 expression was associated with low-grade tumours (P = 0.012) and low expression of Ki67 (proliferation-related Ki-67 antigen) (P = 0.018). Nuclear Msx2 expression was associated with low-grade (P = 0.015), ER-positive (P = 0.038) tumours, low Ki67 expression (P = 0.005), and increased nuclear cyclin D1 expression (P = 0.037).
Increased Msx2 protein expression is associated with prolonged patient survival
The relationship between Msx2 protein expression and survival was then examined. In agreement with our previous findings, cytoplasmic Msx2 was associated with longer breast cancer-specific survival (BCSS) (P = 0.049), recurrence-free survival (RFS) (P = 0.029), and overall survival (OS) (P = 0.019) (Figure 3). In contrast, nuclear Msx2 expression was not associated with outcome at any endpoint: BCSS (P = 0.220), RFS (P = 0.439) or OS (P = 0.123). Univariate Cox regression analysis confirmed that increased cytoplasmic Msx2 expression was associated with an extended BCSS (HR 0.54; 95% CI 0.29-1.01; P = 0.053), RFS (HR 0.57; 95% CI 0.35-0.92; P = 0.021) and OS (HR 0.69; 95% CI 0.49-0.97; P = 0.031). To compare the prognostic impact of Msx2 with well-established clinical variables, a multivariate Cox regression analysis was carried out using OS as an endpoint (Table 3). Cytoplasmic Msx2 expression was an independent predictor of prolonged OS (HR 0.58; 95% CI 0.36-0.93; P = 0.023), along with tumour grade (HR 1.78; 95% CI 1.06-2.99; P = 0.028) and nodal status (HR 3.06; 95% CI 1.95-4.80; P < 0.001). Univariate Cox regression analysis of nuclear Msx2 showed no association with outcome (HR 0.78; 95% CI 0.56-1.07; P = 0.124).
To control for subjectivity in the manual scoring process, we used a co-localisation algorithm (Aperio) to impartially examine Msx2 expression and localisation. We found a strong correlation between automated and manual scores (Spearman's ρ = 0.745; P < 0.001). Using the automated data, cytoplasmic Msx2 expression was associated with overall survival (P = 0.008), whereas nuclear Msx2 expression was not (P = 0.097), in agreement with manual scoring data (Additional file 2).
Msx2 induces cell death by apoptosis in both transformed and immortalized breast cell lines
Having shown that Msx2 expression correlates with good prognosis in breast tumours, we proceeded to examine the role of Msx2 in vitro. Cells suitable for Msx2 overexpression were determined by Western blot analysis (Figure 2a). MDA-MB-231, Hs578t, and MCF10a cells all had undetectable levels of Msx2 protein. We chose a highly invasive breast cancer cell line (MDA-MB-231) and an immortalized mammary epithelial cell line (MCF10a) for our study. Ectopic Msx2 was overexpressed in these cell lines, and expression was verified by Western blot analysis, which showed that both the MDA-MB-231-Msx2 (designated MDA-Msx2) and MCF10a-Msx2 cell lines displayed high levels of Msx2 protein compared to the corresponding empty vector (EV) controls (Figure 4a). We also examined these cell lines by immunofluorescence microscopy using a fluorescently tagged secondary antibody (Figure 4b). This revealed high levels of Msx2 in both MDA-Msx2 and MCF10a-Msx2 cell lines, the expression of which was predominantly nuclear and perinuclear.
We also observed that Msx2-expressing cell lines appeared to be undergoing cell death, whereas the EV-transduced cells grew unperturbed (Figure 4c); thus, the effect of Msx2 overexpression on cell growth, anchorage-independent growth, and cell death was examined (Figure 5). A significant reduction in growth rates, as measured by MTT assay over a period of 5 days, was observed between EV-transduced cells and their Msx2-overexpressing counterparts, namely, MDA-Msx2 (P = 0.005) and MCF10a-Msx2 (P = 0.023) (Figure 5a). Clonogenic assays also revealed a dramatic decrease in cell viability in both MSX2-expressing cell lines (P < 0.001 for both) (Figure 5b). The mechanism of cellular death was then investigated using a luminescence-based assay (Caspase-Glo 3/7), which detects apoptosis by measuring levels of caspase-3 and -7 within cells. This revealed a significant increase in caspase activity in both MDA-Msx2 (P < 0.001) and MCF10a-Msx2 (P = 0.003) cell lines relative to controls, indicating that the mode of cell death was apoptosis (Figure 5c). The proportion of cells in sub-G1, as measured by flow cytometric analysis, was also significantly higher in MDA-Msx2 (P < 0.001) and MCF10a-Msx2 (P < 0.001) cell lines compared to controls (Figure 5d). Taken together, these data indicate that increased levels of Msx2 in breast cells can lead to induction of cell death by apoptosis.
Msx2 overexpression in cell line models leads to an alteration in a number of signalling pathways
To elucidate the mechanism by which Msx2 might induce apoptosis, the expression of several proteins involved in cell signalling was examined in the Msx2-overexpressing and control cells (Figure 4a). We first examined ERK activation and found that ectopic Msx2 expression induced activation of ERK by phosphorylation at Tyr204 in both the MDA-Msx2 and MCF10a-Msx2 cell lines compared to EV controls. We proceeded to look at a number of proteins which regulate cell cycle and apoptosis and found that the Msx2-overexpressing cell lines had notably increased levels of the cyclin-dependent kinase inhibitor p21 (Cip1) and decreased levels of the apoptotic inhibitor Survivin (BIRC5). Levels of cyclin D1 protein were unaffected by ectopic Msx2 expression in these cell lines.
We next focused on proteins involved in transcriptional regulation which may interact with Msx2 and affect its activity as a transcription factor. Levels of the homeobox transcription factor, DLX5, were unaffected by Msx2 overexpression, as was the transcriptional regulator Smad4. Msx2 is known to interact with the SPEN (a derivation of 'split ends' in Drosophila) family of transcriptional regulators [7, 32], and we found that levels of RBM15 (RNA-binding motif protein 15), a member of the SPEN family, were decreased in Msx2-overexpressing cells. Finally, we looked at expression of Twist, a transcription factor known to promote EMT, which has previously been shown to be induced by Msx2 . Twist was expressed in the transformed MDA-MB-231 cell line, but not in the immortalized MCF10a cell line. We saw no change in Twist expression between control and Msx2-overexpressing cells.
Msx2 mRNA expression correlates with expression of cell cycle proteins in tumour samples
To correlate our in vitro cell line data with the clinical situation, we examined the expression of a number of cell survival, cell cycle and Msx2-interacting genes within the previously described Van de Vijver dataset  (Table 4). In relation to cell survival, high levels of Msx2 mRNA were associated with low mRNA levels of the apoptotic inhibitor Survivin (P < 0.001). With regard to cell cycle, increased Msx2 mRNA expression was associated with high cyclin D1 (P < 0.001) and low levels of the cyclins A1 (P = 0.037), A2 (P = 0.007), B2 (P = 0.014), and E1 (P = 0.003). We also examined the Msx2-interacting protein, Dlx5, which was associated with poor prognosis (P = 0.032) but did not significantly correlate with Msx2 expression (P = 0.077). In addition, we found that high Msx2 expression was associated with increased levels of BMP4 (P = 0.017), which has previously been shown to induce Msx2 expression . BMP4 alone was associated with good prognosis (P < 0.001). Expression data on p21 or RBM15 was not available within this dataset.
The objective of this study was to clarify the role of Msx2 in breast cancer. Previous studies have postulated a role for Msx2 in the advancement of the invasive phenotype in breast cancer . However, it has become clear from our analysis of two independent patient cohorts that increased Msx2 mRNA and protein expression is associated with improved outcome in breast cancer. Interestingly, in the TMA-based study, cytoplasmic as opposed to nuclear MSX2 expression was more significant in relation to improved survival, suggesting that the DNA-binding activity of Msx2 may not be the main reason for the beneficial effects of Msx2 overexpression. However, some factors which were not associated with cytoplasmic Msx2, such as ER status and cyclin D1, were associated with nuclear Msx2, indicating that Msx2 may have some additional activity upon translocation to the nucleus.
Our hypothesis that Msx2 is a good prognostic marker in breast cancer was further supported by the subsequent in vitro investigations, when we overexpressed Msx2 in both an immortalized mammary epithelial cell line (MCF10a) and an invasive breast cancer cell line (MDA-MB-231). In both cell lines, a drastic reduction in cell viability, apparently due to induction of apoptosis, was observed. Investigation of how this apoptosis might occur through dissection of signalling pathways showed that ERK activation was increased in Msx2-overexpressing cells. Although ERK activation is known to stimulate cell growth in many situations, it can also induce apoptosis in certain cases, depending on the downstream signals activated . We also found an increase in p21 expression concurrent with Msx2 in these cell lines. Previous studies of BMP4-induced apoptosis in various cell types [21, 22] have shown that both Msx2 and p21 are induced following BMP4 treatment. Inhibition of either molecule alone was sufficient to block cell death, suggesting that both are needed for apoptosis to occur and that they may be part of the same apoptotic cascade. Indeed, the BMP pathway may be involved in modulating the effect of Msx2 expression in breast cancer: analysis of BMP4 expression in 295 tumours from a transcriptomic dataset revealed that Msx2 is associated with high levels of BMP4 in breast tumours, and BMP4 itself was associated with improved patient outcome within this dataset.
In-depth analysis of cell cycle machinery within Msx2-overexpressing cells showed that levels of the apoptosis inhibitor Survivin were decreased, possibly leaving these cells susceptible to apoptosis. Survivin is a significant predictor of poor prognosis in breast cancer . Although we saw no change in levels of cyclin D1 following Msx2 overexpression in vitro, when we examined the van de Vijver dataset, we found that Msx2 was associated with high cyclin D1 and low cyclin A1, A2, B2 and E1. High cyclin D1 is needed for cells to pass the G1 checkpoint, and it was previously thought to be associated with increased proliferation and poor clinical outcome. However, further studies have revealed that cyclin D1 can be associated with improved outcome in breast cancer, possibly due to the blocking effects of p27, which is often increased in parallel with cyclin D1 [35, 36]. In this transcriptomic dataset, cyclin D1 was associated with good prognosis and also with p27 expression (P < 0.001), indicating a possible mechanism whereby high expression of Msx2 can be linked to cell cycle arrest. The fact that Msx2 mRNA expression in these tumours correlates with markers of cell cycle arrest and apoptosis suggests that the apoptosis which we observed in vitro may represent a real phenomenon in breast tumours in vivo. However, further mechanistic analyses are needed to confirm this hypothesis.
Further analysis of the Msx2-overexpressing cell lines revealed that levels of the SPEN family member RBM15 were downregulated concurrently with Msx2 overexpression. Although very little is known about the physiological role of RBM15, it is closely related to SHARP, another member of the SPEN family which is known to interact with Msx2 . This family of proteins are characterized by N-terminal RNA recognition motifs and a conserved SPOC (SPEN paralog and ortholog C-terminal) domain . The conserved SPOC domain in both SHARP and RBM15 is known to interact with SMRT (silencing mediator of retinoid and thyroid hormone receptors) and NCoR (nuclear receptor corepressor) corepressors to mediate transcriptional repression . RBM15 is also involved in regulating activation of Notch signalling , nuclear export of mRNA , and hematopoietic development and cell fate . Deletion of RBM15 in mice is embryonic lethal, and conditional deletion of RBM15 in adult mice results in a reduction in levels of white blood cells . Although the mechanism of action of RBM15 is unclear, the downregulation of this protein in Msx2-overexpressing cell lines indicates that RBM15 may interact directly or indirectly with Msx2 and could mediate, in part, the effects of Msx2 overexpression.
Perhaps it is the background level of Msx2 interacting proteins and cofactors that determines the effect of increased levels of Msx2 in breast tumours. A link with the ER signalling pathway is also likely, as previous studies have shown that Msx2 expression can be induced in breast cancer cell lines and breast explants following estradiol or progesterone treatment [14, 15]. We observed a correlation between Msx2 expression and ER status in the majority of the cell lines that were examined in this study (Figure 2a). The only previous study of Msx2 overexpression in the breast found an association between Msx2 overexpression and induction of EMT and cellular invasion , with no evidence of a reduction in cell viability or induction of apoptosis. However, the earlier study used mouse mammary epithelial cells in which to overexpress Msx2, which would have a different genetic background to the breast cell lines used in this study and may therefore have altered downstream responses to Msx2 overexpression. This variation in the cellular response to Msx2 indicates that the context of Msx2 overexpression, as well as its subcellular localisation and intensity, is an important factor to consider when examining the role of Msx2 in breast tumours. In addition, the intracellular milieu of transcriptional regulators and cofactors may play an important role in regulating the downstream effects of Msx2 overexpression.
The varying effect of Msx2 expression in different tumour types also suggests a cell type-specific effect. Indeed, Msx2 has been associated with high tumour grade in pancreatic cancer . Studies in breast cancer have suggested a link to EMT through induction of Cripto-1 and the c-Src pathway , although this is contradicted by studies showing that Msx2 downregulates the prometastatic factor BSP . The reasons for the context-specific effects of Msx2 overexpression are yet to be explained and may be a result of both changes in cell lineage and protein localisation. However, the results from our clinical study make it clear that Msx2 is a marker of good prognosis in human breast cancer, illustrating the principle that data obtained from in vitro laboratory work must be supplemented by clinical investigation to determine the real prognostic value of a biomarker.
These findings indicate that increased expression of both Msx2 mRNA and protein are associated with improved patient outcome in breast cancer. Furthermore, we have shown that ectopic expression of Msx2 in breast cell lines leads to the induction of apoptosis and radically decreased cell viability. This study suggests that increased Msx2 results in improved outcome for breast cancer patients, possibly by increasing the likelihood of tumour cell death by apoptosis.
breast cancer-specific survival
bone morphogenetic protein
cell pellet array
distal-less homeobox 5
epithelial to mesenchymal transition
extracellular signal-regulated kinase
human epidermal growth factor 2
proliferation-related Ki-67 antigen
muscle segment homeobox 2
RNA-binding motif protein 15
a derivation of 'split ends' (Drosophila)
SPEN paralog and ortholog C-terminal
vascular endothelial growth factor.
Chen YH, Ishii M, Sun J, Sucov HM, Maxson RE: Msx1 and Msx2 regulate survival of secondary heart field precursors and post-migratory proliferation of cardiac neural crest in the outflow tract. Dev Biol. 2007, 308: 421-437. 10.1016/j.ydbio.2007.05.037.
Lallemand Y, Bensoussan V, Cloment CS, Robert B: Msx genes are important apoptosis effectors downstream of the Shh/Gli3 pathway in the limb. Dev Biol. 2009, 331: 189-198. 10.1016/j.ydbio.2009.04.038.
Phippard DJ, Weber-Hall SJ, Sharpe PT, Naylor MS, Jayatalake H, Maas R, Woo I, Roberts-Clark D, Francis-West PH, Liu YH, Maxson R, Hill RE, Dale TC: Regulation of Msx-1, Msx-2, Bmp-2 and Bmp-4 during foetal and postnatal mammary gland development. Development. 1996, 122: 2729-2737.
Satokata I, Ma L, Ohshima H, Bei M, Woo I, Nishizawa K, Maeda T, Takano Y, Uchiyama M, Heaney S, Peters H, Tang Z, Maxson R, Maas R: Msx2 deficiency in mice causes pleiotropic defects in bone growth and ectodermal organ formation. Nat Genet. 2000, 24: 391-395. 10.1038/74231.
Takahashi K, Nuckolls GH, Tanaka O, Semba I, Takahashi I, Dashner R, Shum L, Slavkin HC: Adenovirus-mediated ectopic expression of Msx2 in even-numbered rhombomeres induces apoptotic elimination of cranial neural crest cells in ovo. Development. 1998, 125: 1627-1635.
Wu LY, Li M, Hinton DR, Guo L, Jiang S, Wang JT, Zeng A, Xie JB, Snead M, Shuler C, Maxson RE, Liu YH: Microphthalmia resulting from MSX2-induced apoptosis in the optic vesicle. Invest Ophthalmol Vis Sci. 2003, 44: 2404-2412. 10.1167/iovs.02-0317.
Satoh K, Ginsburg E, Vonderhaar BK: Msx-1 and Msx-2 in mammary gland development. J Mammary Gland Biol Neoplasia. 2004, 9: 195-205. 10.1023/B:JOMG.0000037162.84758.b5.
Ma L, Golden S, Wu L, Maxson R: The molecular basis of Boston-type craniosynostosis: the Pro148-->His mutation in the N-terminal arm of the MSX2 homeodomain stabilizes DNA binding without altering nucleotide sequence preferences. Hum Mol Genet. 1996, 5: 1915-1920. 10.1093/hmg/5.12.1915.
Wilkie AO, Tang Z, Elanko N, Walsh S, Twigg SR, Hurst JA, Wall SA, Chrzanowska KH, Maxson RE: Functional haploinsufficiency of the human homeobox gene MSX2 causes defects in skull ossification. Nat Genet. 2000, 24: 387-390. 10.1038/74224.
Abate-Shen C: Deregulated homeobox gene expression in cancer: cause or consequence?. Nat Rev Cancer. 2002, 2: 777-785. 10.1038/nrc907.
Lanigan F, O'Connor D, Martin F, Gallagher WM: Molecular links between mammary gland development and breast cancer. Cell Mol Life Sci. 2007, 64: 3159-3184. 10.1007/s00018-007-7386-2.
Humphreys RC, Krajewska M, Krnacik S, Jaeger R, Weiher H, Krajewski S, Reed JC, Rosen JM: Apoptosis in the terminal endbud of the murine mammary gland: a mechanism of ductal morphogenesis. Development. 1996, 122: 4013-4022.
McBryan J, Howlin J, Kenny PA, Shioda T, Martin F: ERalpha-CITED1 co-regulated genes expressed during pubertal mammary gland development: implications for breast cancer prognosis. Oncogene. 2007, 26: 6406-6419. 10.1038/sj.onc.1210468.
Satoh K, Hovey RC, Malewski T, Warri A, Goldhar AS, Ginsburg E, Saito K, Lydon JP, Vonderhaar BK: Progesterone enhances branching morphogenesis in the mouse mammary gland by increased expression of Msx2. Oncogene. 2007, 26: 7526-7534. 10.1038/sj.onc.1210555.
Malewski T, Milewicz T, Krzysiek J, Gregoraszczuk EL, Augustowska K: Regulation of Msx2 gene expression by steroid hormones in human nonmalignant and malignant breast cancer explants cultured in vitro. Cancer Invest. 2005, 23: 222-228. 10.1081/CNV-200055958.
di Bari MG, Ginsburg E, Plant J, Strizzi L, Salomon DS, Vonderhaar BK: Msx2 induces epithelial-mesenchymal transition in mouse mammary epithelial cells through upregulation of Cripto-1. J Cell Physiol. 2009, 219: 659-666. 10.1002/jcp.21712.
Satoh K, Hamada S, Kimura K, Kanno A, Hirota M, Umino J, Fujibuchi W, Masamune A, Tanaka N, Miura K, Egawa S, Motoi F, Unno M, Vonderhaar BK, Shimosegawa T: Up-regulation of MSX2 enhances the malignant phenotype and is associated with twist 1 expression in human pancreatic cancer cells. Am J Pathol. 2008, 172: 926-939. 10.2353/ajpath.2008.070346.
Barnes GL, Javed A, Waller SM, Kamal MH, Hebert KE, Hassan MQ, Bellahcene A, Van Wijnen AJ, Young MF, Lian JB, Stein GS, Gerstenfeld LC: Osteoblast-related transcription factors Runx2 (Cbfa1/AML3) and MSX2 mediate the expression of bone sialoprotein in human metastatic breast cancer cells. Cancer Res. 2003, 63: 2631-2637.
Bellahcene A, Menard S, Bufalino R, Moreau L, Castronovo V: Expression of bone sialoprotein in primary human breast cancer is associated with poor survival. Int J Cancer. 1996, 69: 350-353. 10.1002/(SICI)1097-0215(19960822)69:4<350::AID-IJC19>3.0.CO;2-9.
Marazzi G, Wang Y, Sassoon D: Msx2 is a transcriptional regulator in the BMP4-mediated programmed cell death pathway. Dev Biol. 1997, 186: 127-138. 10.1006/dbio.1997.8576.
Gomes WA, Kessler JA: Msx-2 and p21 mediate the pro-apoptotic but not the anti-proliferative effects of BMP4 on cultured sympathetic neuroblasts. Dev Biol. 2001, 237: 212-221. 10.1006/dbio.2001.0344.
Israsena N, Kessler JA: Msx2 and p21(CIP1/WAF1) mediate the proapoptotic effects of bone morphogenetic protein-4 on ventricular zone progenitor cells. J Neurosci Res. 2002, 69: 803-809. 10.1002/jnr.10362.
van de Vijver MJ, He YD, van't Veer LJ, Dai H, Hart AA, Voskuil DW, Schreiber GJ, Peterse JL, Roberts C, Marton MJ, Parrish M, Atsma D, Witteveen A, Glas A, Delahaye L, van der Velde T, Bartelink H, Rodenhuis S, Rutgers ET, Friend SH, Bernards R: A gene-expression signature as a predictor of survival in breast cancer. N Engl J Med. 2002, 347: 1999-2009. 10.1056/NEJMoa021967.
Merck & Co., Inc: [http://www.rii.com]
O'Brien SL, Fagan A, Fox EJ, Millikan RC, Culhane AC, Brennan DJ, McCann AH, Hegarty S, Moyna S, Duffy MJ, Higgins DG, Jirstrom K, Landberg G, Gallagher WM: CENP-F expression is associated with poor prognosis and chromosomal instability in patients with primary breast cancer. Int J Cancer. 2007, 120: 1434-1443. 10.1002/ijc.22413.
Helczynska K, Larsson AM, Holmquist Mengelbier L, Bridges E, Fredlund E, Borgquist S, Landberg G, Pahlman S, Jirstrom K: Hypoxia-inducible factor-2alpha correlates to distant recurrence and poor outcome in invasive breast cancer. Cancer Res. 2008, 68: 9212-9220. 10.1158/0008-5472.CAN-08-1135.
Lanigan F, McKiernan E, Brennan DJ, Hegarty S, Millikan RC, McBryan J, Jirstrom K, Landberg G, Martin F, Duffy MJ, Gallagher WM: Increased claudin-4 expression is associated with poor prognosis and high tumour grade in breast cancer. Int J Cancer. 2009, 124: 2088-2097. 10.1002/ijc.24159.
Hughes L, Malone C, Chumsri S, Burger AM, McDonnell S: Characterisation of breast cancer cell lines and establishment of a novel isogenic subclone to study migration, invasion and tumourigenicity. Clin Exp Metastasis. 2008, 25: 549-557. 10.1007/s10585-008-9169-z.
Borgquist S, Holm C, Stendahl M, Anagnostaki L, Landberg G, Jirstrom K: Oestrogen receptors alpha and beta show different associations to clinicopathological parameters and their co-expression might predict a better response to endocrine treatment in breast cancer. J Clin Pathol. 2008, 61: 197-203. 10.1136/jcp.2006.040378.
van 't Veer LJ, Dai H, van de Vijver MJ, He YD, Hart AA, Mao M, Peterse HL, van der Kooy K, Marton MJ, Witteveen AT, Schreiber GJ, Kerkhoven RM, Roberts C, Linsley PS, Bernards R, Friend SH: Gene expression profiling predicts clinical outcome of breast cancer. Nature. 2002, 415: 530-536. 10.1038/415530a.
Sorlie T, Perou CM, Tibshirani R, Aas T, Geisler S, Johnsen H, Hastie T, Eisen MB, van de Rijn M, Jeffrey SS, Thorsen T, Quist H, Matese JC, Brown PO, Botstein D, Eystein Lonning P, Borresen-Dale AL: Gene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications. Proc Natl Acad Sci USA. 2001, 98: 10869-10874. 10.1073/pnas.191367098.
Newberry EP, Latifi T, Towler DA: The RRM domain of MINT, a novel Msx2 binding protein, recognizes and regulates the rat osteocalcin promoter. Biochemistry. 1999, 38: 10678-10690. 10.1021/bi990967j.
Mebratu Y, Tesfaigzi Y: How ERK1/2 activation controls cell proliferation and cell death: Is subcellular localization the answer?. Cell Cycle. 2009, 8: 1168-1175. 10.4161/cc.8.8.8147.
Brennan DJ, Rexhepaj E, O'Brien SL, McSherry E, O'Connor DP, Fagan A, Culhane AC, Higgins DG, Jirstrom K, Millikan RC, Landberg G, Duffy MJ, Hewitt SM, Gallagher WM: Altered cytoplasmic-to-nuclear ratio of survivin is a prognostic indicator in breast cancer. Clin Cancer Res. 2008, 14: 2681-2689. 10.1158/1078-0432.CCR-07-1760.
Barnes DM, Gillett CE: Cyclin D1 in breast cancer. Breast Cancer Res Treat. 1998, 52: 1-15. 10.1023/A:1006103831990.
Bostrom P, Soderstrom M, Palokangas T, Vahlberg T, Collan Y, Carpen O, Hirsimaki P: Analysis of cyclins A, B1, D1 and E in breast cancer in relation to tumour grade and other prognostic factors. BMC Res Notes. 2009, 2: 140-10.1186/1756-0500-2-140.
Sanchez-Pulido L, Rojas AM, van Wely KH, Martinez AC, Valencia A: SPOC: a widely distributed domain associated with cancer, apoptosis and transcription. BMC Bioinformatics. 2004, 5: 91-10.1186/1471-2105-5-91.
Hiriart E, Gruffat H, Buisson M, Mikaelian I, Keppler S, Meresse P, Mercher T, Bernard OA, Sergeant A, Manet E: Interaction of the Epstein-Barr virus mRNA export factor EB2 with human Spen proteins SHARP, OTT1, and a novel member of the family, OTT3, links Spen proteins with splicing regulation and mRNA export. J Biol Chem. 2005, 280: 36935-36945. 10.1074/jbc.M501725200.
Ma X, Renda MJ, Wang L, Cheng EC, Niu C, Morris SW, Chi AS, Krause DS: Rbm15 modulates Notch-induced transcriptional activation and affects myeloid differentiation. Mol Cell Biol. 2007, 27: 3056-3064. 10.1128/MCB.01339-06.
Lindtner S, Zolotukhin AS, Uranishi H, Bear J, Kulkarni V, Smulevitch S, Samiotaki M, Panayotou G, Felber BK, Pavlakis GN: RNA-binding motif protein 15 binds to the RNA transport element RTE and provides a direct link to the NXF1 export pathway. J Biol Chem. 2006, 281: 36915-36928. 10.1074/jbc.M608745200.
Raffel GD, Mercher T, Shigematsu H, Williams IR, Cullen DE, Akashi K, Bernard OA, Gilliland DG: Ott1(Rbm15) has pleiotropic roles in hematopoietic development. Proc Natl Acad Sci USA. 2007, 104: 6001-6006. 10.1073/pnas.0609041104.
We thank Dr. Darran O'Connor and Dr. Sharon McGee for helpful discussions on virus experiments. Funding is acknowledged from Enterprise Ireland, the Irish Research Council for Science Engineering and Technology under the EMBARK initiative (for support of Ms. Fiona Lanigan's and Ms. Gabriela Gremel's postgraduate studies), Cancer Research Ireland (for part-support of Dr. Brennan's postgraduate studies), Science Foundation Ireland, the British Association for Cancer Research and the Health Research Board of Ireland, the latter under the auspices of the 'Breast Cancer Metastasis: Biomarkers and Functional Mediators' research programme. This study has also been supported by grants from the Swedish Cancer Society, Swegene/Wallenberg Consortium North, Gunnar, Arvid and Elisabeth Nilsson Cancer Foundation, Per-Eric and Ulla Schyberg Foundation, Lund University Research Funds and Malmö University Hospital Research and Cancer Funds. Finally, the cross-national component of the project was facilitated by the Marie Curie Transfer of Knowledge Industry-Academia Partnership research programme, TargetBreast . The UCD Conway Institute is funded by the Programme for Third Level Institutions (PRTLI), as administered by the Higher Education Authority (HEA) of Ireland.
The authors declare that they have no competing interests.
FL participated in the design of the study, carried out the transcriptomic dataset mining, antibody validation, TMA staining, lentiviral overexpression, and functional analysis; performed the statistical analysis, and drafted the manuscript. GG participated in the design of the lentiviral study, and designed and validated the lentiviral construct used in this study. RH and FL scored the immunohistochemical staining on the TMA. DB, FM and WG participated in the design of the study and edited the manuscript. KJ constructed the TMAs and collated the clinical information. All authors read and approved the final manuscript.