Expression of epithelial to mesenchymal transition-related markers in lymph node metastases as a surrogate for primary tumor metastatic potential in breast cancer
- Aleksandra Markiewicz†1, 2,
- Tomasz Ahrends†1,
- Marzena Wełnicka-Jaśkiewicz3,
- Barbara Seroczyńska4,
- Jarosław Skokowski4, 5,
- Janusz Jaśkiewicz5,
- Jolanta Szade6,
- Wojciech Biernat6 and
- Anna J Żaczek1Email author
© Markiewicz et al.; licensee BioMed Central Ltd. 2012
Received: 3 August 2012
Accepted: 14 November 2012
Published: 19 November 2012
Breast cancers are phenotypically and genotypically heterogeneous tumors containing multiple cancer cell populations with various metastatic potential. Aggressive tumor cell subpopulations might more easily be captured in lymph nodes metastases (LNM) than in primary tumors (PT). We evaluated mRNA and protein levels of master EMT regulators: TWIST1, SNAIL and SLUG, protein levels of EMT-related markers: E-cadherin, vimentin, and expression of classical breast cancer receptors: HER2, ER and PgR in PT and corresponding LNM. The results were correlated with clinicopathological data and patients outcomes.
Formalin-fixed paraffin-embedded samples from PT and matched LNM from 42 stage II-III breast cancer patients were examined. Expression of TWIST1, SNAIL and SLUG was measured by reverse-transcription quantitative PCR. Protein expression was examined by immunohistochemistry on tissue microarrays. Kaplan-Meier curves for disease-free survival (DFS) and overall survival (OS) were compared using F-Cox test. Hazard ratios (HRs) with 95% confidence intervals (95% CI) were computed using Cox regression analysis.
On average, mRNA expression of TWIST1, SNAIL and SLUG was significantly higher in LNM compared to PT (P < 0.00001 for all). Gene and protein levels of TWIST1, SNAIL and SLUG were highly discordant between PT and matched LNM. Increased mRNA expression of TWIST1 and SNAIL in LNM was associated with shorter OS (P = 0.04 and P = 0.02, respectively) and DFS (P = 0.02 and P = 0.01, respectively), whereas their expression in PT had no prognostic impact. Negative-to-positive switch of SNAIL protein correlated with decreased OS and DFS (HR = 4.6; 1.1-18.7; P = 0.03 and HR = 3.8; 1.0-48.7; P = 0.05, respectively).
LNM are enriched in cells with more aggressive phenotype, marked by elevated levels of EMT regulators. High expression of TWIST1 and SNAIL in LNM, as well as negative-to-positive conversion of SNAIL confer worse prognosis, confirming the correlation of EMT with aggressive disease behavior. Thus, molecular profiling of LNM may be used as surrogate marker for aggressiveness and metastatic potential of PT.
KeywordsBreast cancer Primary tumor Lymph node metastasis Gene expression Immunohistochemistry Biomarkers Epithelial-mesenchymal transition
Distant metastasis remain the main cause of cancer-related death in breast cancer. Metastatic lymph node involvement is still the most powerful prognostic factor for relapse and death, but lymph node dissection does not affect patients survival . It is still the matter of debate if positive lymph nodes are able to metastasize . Currently, lymph node metastases (LNM) are considered rather a manifestation of the widespread metastatic process and a marker of aggressive phenotype of the primary tumor (PT) than the “bridge-heads” for further metastatic spread . It has been confirmed by the clinical observation of poorer survival after relapse in node-positive patients compared to node-negative ones . Experimental models provide further evidence that development of LNM indicates the increased potential of PT to disseminate aggressive cells and produce metastasis promoting growth factors , according to the recently proposed stromal progression model . In this model mutual regulatory interactions between stroma and tumor cells play equally important roles in tumor progression as genetic and epigenetic changes. Those interactions contribute to the process of epithelial-mesenchymal transition (EMT), similarly to the bone marrow mesenchymal stem cells actively recruited by tumor cells to the surrounding stroma .
EMT has been found to be crucial in cancer dissemination, endowing cells with metastatic and cancer stem cell properties [7, 8]. It is characterized by downregulation of epithelial markers (e.g. cytokeratin 8, 18, 19, E-cadherin, claudins, occludins) and upregulation of mesenchymal markers (e.g. vimentin, N-cadherin) , what results in numerous phenotypic changes such as the loss of cell-cell adhesion and cell polarity, and the acquisition of migratory and invasive properties [10, 11]. TWIST1, SNAIL and SLUG are transcription factors among those governing EMT (EMT-TFs) . Increased expression of EMT-TFs examined in PT has been associated with poor prognostic clinico-pathological features and outcome in breast and other cancers [12–14] as well as multidrug resistance . However, to the best of our knowledge there are no data available on the status of TWIST1, SNAIL and SLUG in LNM.
Phenotype of the aggressive cancer cells subpopulations resulting from EMT might not be easily captured in the PT since they form a minority of cells within tumor bulk [16–18]. Moreover, expression of genes associated with EMT is transient and space-limited . We hypothesized that the analysis of pre-selected subpopulations of cancer cells found in LNM could be more informative in terms of aggressiveness than the analysis of the PT bulk. The primary objective of the present study was to evaluate possible changes occurring in classical and EMT-related marker status between PT and corresponding synchronous axillary LNM before initiation of the therapy and to relate them to clinical outcome. Additionally, we aimed to investigate the feasibility of the quantitative PCR-based gene expression profiling of low-level transcripts TWIST1, SNAIL and SLUG in formalin-fixed, paraffin-embedded (FFPE) tissues compared to matched frozen counterparts.
Material and methods
Studied material included 44 tissue specimens from patients with operable breast cancer and lymph node involvement who were treated between 2006–2008 at the Medical University of Gdansk Hospital. Patients were treated with surgery by modified radical mastectomy or local tumor resection, with axillary node dissection followed by postoperative breast irradiation. Adjuvant therapy with chemotherapy and/or hormone therapy was given in standard care settings based on the nodal and hormone receptor status. Availability of PT and matched LNM was mandatory. Patients with no evidence of lymph node involvement or earlier chemotherapy were not eligible for this study. Non-cancer control breast tissue samples were acquired during mastectomy ensuring the greatest possible distance to the main tumor mass, and sections of non-involved lymph nodes were collected. The study was conducted in accordance with the Declaration of Helsinki and approved by the Ethics Committee of the Medical University of Gdansk. All patients signed informed consent forms.
RNA extraction from formalin-fixed paraffin-embedded (FFPE) tissue
Tissue specimens were fixed in 10% (v/v) neutral-buffered formalin for up to 24 h, dehydrated in 70% ethanol and embedded in paraffin. FFPE tissue blocks were stored at room temperature for up to 6 years. The percentage of tumor cells in each FFPE specimen was evaluated by hematoxylin and eosin staining reviewed by a certified pathologist. Only the tissue section with confirmed presence of invasive carcinoma and tumor cells content over 50% were included. 2–4 slices of 10 μm thickness were cut using a microtome and placed in 1.5 ml centrifuge tubes. Tissues were de-paraffinized by treatment with xylene and 100% ethanol. Total RNA was isolated using RNeasy FFPE Kit (Qiagen, Germany) according to the manufacturer's protocol, including on-column DNase I treatment.
RNA extraction from fresh-frozen (FF) tissue
After collection, tissue samples were immediately frozen in liquid nitrogen and stored at −80°C for further analysis. 20–30 mg tissue sections were homogenized with zircon beads in MagNA Lyzer (Roche, Germany) for 40 s. Total RNA was isolated using RNeasy Mini Kit (Qiagen, Germany) according to the manufacturer’s protocol, including on-column DNase I treatment.
RNA analysis and reverse transcription
For all samples RNA concentration and purity was determined using the Nano-Drop ND-1000 spectrophotometer (Thermo Scientific, USA). Qualitative analysis of RNA was performed by microcapillary electrophoresis using the Agilent 2100 Bioanalyzer (Expert software version B.02.08) with an RNA Nano Chip (Agilent Technologies, UK). For each sample, whenever possible, 1 μg of RNA was used as the template in the reverse transcription reaction using Transcriptor First Strand cDNA Synthesis Kit (Roche, Germany) in a 20 μl volume with random hexamer primers, according to the manufacturer's protocol.
qPCR and gene expression analysis
Gene expression levels were determined by RT-qPCR on the CFX96 Real-Time System (Bio-Rad, USA) with specific TaqMan Gene Expression Assays (Applied Biosystems/Life Technologies, USA) containing validated primers and probe set spanning exon-exon boundaries for detection of TWIST1 (Hs00361186_m1), SNAIL (Hs00195591_m1) and SLUG (Hs00950344_m1). Relative expression values of each gene were calculated by the delta-delta-Ct method normalized to the reference gene GAPDH (glyceraldehyde 3-phosphate dehydrogenase, Hs99999905_m1), and non-cancer breast tissue as a calibrator with the use of qBasePLUS software (Biogazelle, ver. 2.0).
Stability of four reference genes: β-actin (Hs99999903_m1), HPRT (Hs99999909_m1), GAPDH and YWHAZ (Hs03044281_g1) was evaluated in 11 breast cancer samples. GAPDH was chosen as a reference gene based on its highest expression stability measure (M) calculated in GeNorm.
The qPCR cycling profile was programmed as follows: predenaturation at 95°C for 10 minutes, followed by amplification phase, which included denaturation at 95°C for 15 seconds, annealing and extension at 60°C for 60 seconds for 45 cycles. 40 ng of cDNA in 4 μl was added per reaction and mixed with 10 μl of TaqMan Universal PCR Master Mix (Applied Biosystems/Life Technologies, USA), 5 μl of water and 1 μl of specific primer and probe mix (20x concentrated). To verify that the qPCR signals derived from cDNA, not genomic DNA, for each gene tested a control without reverse transcriptase in the RT reaction (no RT control) was included. Each sample was analyzed in duplicate, and the average Ct value of duplicates was used as a quantitative value. The specificity of the polymerase chain reaction was confirmed by gel electrophoresis using 1.5% agarose gel containing Gel-Red (Gentaur, Belgium) and viewed under ultraviolet illuminator Gel Doc (Bio-Rad, USA). Positive result was defined when a relative gene expression was higher than median expression in all tumor samples.
Efficiency of the reaction
Sensitivity: limit of detection (LOD) and limit of quantification (LOQ)
Intra- and interassay variation
Detailed description of the standardization methodology is described in the Additional file 1.
Immunohistochemistry (IHC) on tissue microarrays (TMA)
TMAs were prepared as described before . Protein expression was examined using antibodies against TWIST1 [ab50581, Abcam], SNAIL [ab85936, Abcam], SLUG [ab38551, Abcam], E-cadherin [Nch 38, Dako], vimentin [V9, Dako], HER2 [HercepTest Kit, Dako], ER [1D5, Dako], PgR [636, Dako] and peroxidase-based detection system (Novolink Polymer Detection System, Novocastra, UK) in accordance with the manufacturer’s guidelines. Antigen retrieval was carried out by heat induced epitope retrieval at pH 6. Positive and negative controls for each marker were used according to the supplier’s data sheet. The material was analyzed using a transmission light microscope (Olympus BX 41) with 400x magnification. Two cores from each tumor were assessed individually. IHC analysis was performed by two independent observers blinded to the clinical data and patients outcomes. Discordant results were reviewed to achieve agreement. The same protocol of staining and scoring was used for both PT and LNM. There was half a year difference in staining of PT and LNM.
ER and PgR were scored according to classical Allred system with cut-point 3 for positive result, while HER2 - according to HercepTest criteria, with 3+ score defining positive result. For vimentin and E-cadherin positive result was considered when 10% or more cells stained positively [22, 23]. For TWIST1, SNAIL and SLUG only nuclear staining was considered with a 10% cut off value of positivity [22–24]. Results were considered concordant if PT and LNM were both positive or both negative, other combinations were considered discordant and denoted the conversion rate.
To fulfill the primary objective of determining the probability of conversion in biomarker status, 28 paired samples were required to detect a discordance rate of 20% with 80% power using a one-sided alpha of 5%.
Concordance between PT and LNM was measured by estimating Cohen’s kappa coefficient (κ) with Medcalc software, version 18.104.22.168 (MedCalc Software, Belgium). The level of agreement based on κ values was assessed using the Landis and Koch criteria: 0.00-0.20, slight agreement; 0.21-0.40, fair agreement; 0.41-0.60, moderate agreement; 0.61-0.80, substantial agreement; and 0.81-1.00, almost perfect agreement .
Categorical variables were compared by Fisher’s exact test, and continuous variables were compared by the Spearman’s rank order test. Kaplan-Meier curves for disease-free survival (DFS) and overall survival (OS) were compared using F-Cox test. DFS was defined as the time from surgery to an event (local or distant relapse, second malignancy or death, whichever came first) or censoring. A censoring was defined as lost to follow-up or alive without relapse at the end of follow-up. Hazard ratios (HRs) with 95% confidence intervals (95% CI) were computed using Cox regression analysis. Significance was defined as P ≤ 0.05. STATISTICA software version 10.0 for Windows was used for all statistical analyses.
Number of cases N = 44
Biomarkers’ conversion rate between paired PT and LNM
Positive in PT
Positive in LNM
Conversion rate PT → LNM
(−) → (+) N (%)
(+) → (−) N (%)
kappa coefficient (95% CI)
IHC staining was successful for paired samples of PT and LNM for TWIST1, SNAIL and SLUG in 88%, 86% and 88% of cases (Figure 1B). Cores were missing, folded, or contained no invasive cancer in 12%, 14%, and 12%, respectively. Exemplary photographs of immunohistochemical staining of lymph nodes metastases are presented in Additional file 3: Figure S3 and in Additional file 4: Figure S4.
The conversion rates in TWIST1, SNAIL and SLUG protein status between PT and paired LNM were 46%, 53% and 28%, respectively (Table 2). Negative-to-positive conversion of all TF occurred more frequently than positive-to-negative one (Table 2).
IHC staining of TWIST1 protein was concordant with expression levels of TWIST1 mRNA measured by RT-qPCR in 83% (P = 0.002). No correlation was observed for other TFs. The conversion rates of E-cadherin and vimentin were much lower (0% and 7%, respectively), while the HER2, ER and PgR conversion occurred in 3%, 18% and 29% of cases, respectively.
Expression of TWIST1, SNAIL and SLUG on mRNA and protein level in PT and LNM was correlated with the number of involved lymph nodes. Increased number of involved lymph nodes (more than 3) correlated with elevated expression of TWIST1 protein in LNM (P = 0.02) and showed a trend towards increased expression of TWIST1 mRNA in LNM (P = 0.07) and SNAIL protein in LNM (P = 0.07) (in Additional file 5: Table S2).
RT-qPCR standardization results
Comparison of RNA quantity, quality and RT-qPCR performance for FF and FFPE derived templates
RNA yield [ng/1 mg tissue]*
210 ± 217
227 ± 159
median 260/280 ratio
9.5 ± 0.6
2.1 ± 0.5
efficiency of qPCR
29.9 ± 2.1
35.4 ± 2.7
interassay variation CV
intraassay variation CV
Despite the years of routine use of lymph node dissection in breast cancer management, inspired by Halstedt , our understanding of the role of lymph nodes in the metastatic process is still marginal. The analysis of pre-selected subpopulations of cancer cells found in lymph nodes metastases could provide insights into biological background of cancer progression. Our work explored patterns of conversion in classical and EMT-related biomarker status between primary breast tumors and corresponding synchronous axillary lymph node metastases to determine whether phenotypic variability is associated with different clinical outcome.
Numerous studies have shown discordant expression of classical molecular markers, ER, PgR and HER2 between the PT and both lymph nodes and distant metastases [28–32]. Discordance in ER and PgR status has been demonstrated within the range of 10-32% and 34%–41% cases, respectively [28, 33–35] and 3%–24% for HER2 status [28, 32, 34–36] when PT were compared with metastatic relapse. Generally, lower discordance rates have been observed in LNM, as showed in individual studies [30, 31, 37] and confirmed in recent meta-analysis on HER2 status in primary and metastatic cancer, including 26 studies of 2520 subjects . However, there are also reports showing the opposite, as exemplified by the study of Aitken et al., who found different breast/node status of at least one receptor (ER, PgR or HER2) in almost 47% of cases .
Receptor discordance has been frequently associated with poor survival [32, 34, 35], what has been attributed to the inappropriate use of hormone and targeted therapy prescribed based on the characteristics of PT  or selection of tumors with more unstable phenotype and therefore more aggressive behavior . Adverse impact of receptor discordance was abolished, when treatment was modified according to the results of metastatic biopsy . In general, the cancer management scheme in 14-20% patients was changed on the basis of metastasis biopsy [28, 33, 39]. The growing body of evidence supports reassessment of ER, PgR and HER2 at the time of relapse diagnosis to tailor the most effective treatment for each patient at all times .
Until now, main research focus has been put on the discordance of ER, PgR and HER2 status between the PT and metastasis, as discussed above. Some studies examined cell proliferation, differentiation and apoptosis markers [30, 37, 41]. However, to the best of our knowledge there is no data available on the status of EMT regulators: TWIST1, SNAIL and SLUG in LNM compared to PT. The clinical outcome in relation to these biomarker status in lymph nodes has not yet been reported.
Our results show frequent conversion of TWIST1, SNAIL and SLUG status between PT and LNM, occurring both at mRNA and protein level (within the range of 28-53% of cases). Expression levels of EMT-TFs were significantly higher in LNM compared to PT. This could be explained by the higher frequency of pre-selected aggressive subpopulations of cancer cells resulting from EMT present in LNM than in the PT, where these aggressive cells constitute a minority of cells [16, 17] and therefore might not be easily captured. It seems that indeed molecular profile of LNM might be a surrogate marker of aggressive phenotype of the PT, as postulated .
We have found that both the elevated level of TFs in LNM and their negative-to-positive switch were associated with poor clinical outcome. Increased TWIST1 and SNAIL expression in LNM correlated with shorter OS (P = 0.04 and P = 0.02, respectively) and DFS (P = 0.02 and P = 0.01, respectively), whereas their expression in PT had no prognostic impact. Negative-to-positive conversion of SNAIL status also correlated with worse survival compared to unchanged status (OS: HR = 4.6; 1.1-18.7; P = 0.03; DFS: HR = 3.8; 1.0-48.7; P = 0.05). Numerous studies show poor prognostic impact of increased expression of TWIST1, SNAIL and SLUG in PT [12–14, 42–44]. The contradictory results have also been reported . But no study have examined their clinical outcome in lymph nodes until now.
We also described the protein status of EMT markers – E-cadherin and vimentin in PT and LNM. We have observed no or little change in E-cadherin and vimentin status between PT and LNM. Our results and the results of other groups  show that E-cadherin is at least as frequently expressed in LNM as in primary tumors and it is not downregulated in LNM, what would be expected if EMT was involved. However, this observation does not exclude the possibility that E-cadherin undergoes transient downregulation during EMT and is re-expressed at the stage of circulating tumor cells to facilitate adhesion and metastases formation . Similar to other groups we observed that vimentin is rarely expressed in LNM  or is reduced in comparison to the PT . However, expression of vimentin in metastases was also shown to be heterogeneous, metastasis size and side dependent , what supports the notion that vimentin expression is transient and environment-controlled. As the observed change in EMT-TFs status between PT and LNM is a new finding, a question arises if it results from true biological variation or from inconsistent measurement.
There are many hypotheses that could explain biomarkers change between PT and corresponding metastasis at biological level. Clonal selection, with subsequent clonal expansion during tumor progression, has been proposed as the mechanism inducing the differences in the genetic composition of primary and metastatic breast cancer [50, 51]. Clonal selection may be related to intra-tumor heterogeneity [52, 53] and/or to various selective pressures such as the immune surveillance of the host, stromal or growth factor interactions, nutritional deficiencies, hypoxia, and therapy. Indeed, Niikura reported recently a significantly higher discordance rate in HER2 status among women who received chemotherapy than among those who did not . To the biological variation between primary and metastatic tumor could also contribute independent evolution of an early stem cell clones in both sites, instead of a linear progression from the PT to metastasis .
The discordance between PT and LNM may also be caused by inconsistent measurement resulting from numerous technical issues . For example, recent meta-analysis revealed 15 pre-analytical variables capable of impacting IHC, including fixation delay, fixative type, time in fixative, reagents and conditions of dehydration, clearing, and paraffin impregnation, and conditions of slide drying and storage . Analytical procedures of antigen retrieval, immunostaining and interpretation of results add additional variance to the final result . To ensure maximal consistency, in our study all specimens were from single cancer hospital, examined by the same two observers with the same protocol of staining and scoring for PT and LNM. Additionally, independent observers were blinded to the clinical data and patient outcome.
The results should be, however, interpreted with caution because there is a possibility of selection bias due to retrospective character of the study. Another limitation might be the small sample size. But even with that number of patients the study had still sufficient power to detect reliably the 20% difference in biomarker status, and all of the discussed differences are above that value. It has been claimed that to reduce reporting errors the use of confirmatory test is recommended . To strengthen our study we confirmed increase of EMT-TFs at both mRNA and protein level with RT-qPCR and IHC, respectively. Taking into consideration discussed biological and technical issues, it seems that current study demonstrates true biological variation in TWIST1, SNAIL and SLUG in PT and LNM.
To make gene expression most reliable we have carefully validated the method we used. Since gene expression profiling in lymph nodes could only be examined in formalin-fixed, paraffin-embedded (FFPE) material, we have undertaken the methodological substudy to assess the reliability of obtained results of RT-qPCR. We aimed to investigate the feasibility of RT-qPCR-based gene expression profiling of low-level transcripts TWIST1, SNAIL and SLUG in FFPE tissues compared to matched frozen counterparts.
Expression of genes associated with EMT is transient and space-limited , what makes them specially difficult to study in clinical setting. Technical difficulties that interfere with reliable gene expression analysis can be particularly prominent in FFPE samples, in which RNA is degraded and chemically modified, resulting in lower RNA Integrity Number (RIN) values and rendering it inaccessible for amplification.
In order to examine to what extent tissue fixation influences RT-qPCR based methods we have measured sensitivity of the method in use by defining limits of detection and quantification (LOD and LOQ). No sample presented relative gene expression value of TWIST1 lower than calculated LOD and LOQ, which confirms the reliability of our data, as experimental results less than the theoretically possible LOD should never be reported.
The overall profile of TWIST1, SNAIL and SLUG expression was similar to that generated with well-preserved RNA from matched FF tissue. We presume however, that decreased functionality of RNA (due to degradation or modification) is the underlying factor responsible for increased intra- and interassay variation performed on FFPE samples. This process did not affect the sensitivity of the method as the assay measures relative quantity of the gene expression, which does not change due to similar degree of reference gene and examined genes degradation.
Thus, despite low RIN values RNA from FFPE samples may work fine in quantitative PCR in terms of specificity, sensitivity and reproducibility, even for low-level transcripts, as exemplified by TWIST1, SNAIL and SLUG.
In conclusion, this is the first study to date demonstrating that samples derived from LNM are enriched in cells with more aggressive phenotype marked by elevated levels of EMT regulators. High expression of these markers correlates with increased number of involved lymph nodes and decreased survival of breast cancer patients. We postulate that examination of molecular profile of lymph node metastasis could provide information about aggressiveness potential of the PT and the characteristics of seeded cells. It therefore has the potential to be developed into a powerful prognostic biomarker in breast cancer, ensuring an accurate prognosis and treatment selection. A prospective study of larger cohort of patients would be necessary to confirm the clinical significance of the changes in EMT-related molecular markers.
Formalin-fixed, paraffin-embedded specimen
Human epidermal growth factor receptor 2, LNM, lymph nodes metastases
RNA integrity number
This research was supported by the grants from Ministry of Science and Higher Education: IP2010 050370 and IP2011 064571.
- Veronesi U, Marubini E, Mariani L, Valagussa P, Zucali R: The dissection of internal mammary nodes does not improve the survival of breast cancer patients. 30-year results of a randomised trial. Eur J Cancer. 1999, 35: 1320-1325. 10.1016/S0959-8049(99)00133-1.View ArticlePubMedGoogle Scholar
- Engel J, Emeny RT, Holzel D: Positive lymph nodes do not metastasize. Cancer Metastasis Rev. 2012, 31: 235-246. 10.1007/s10555-011-9343-7.View ArticlePubMedGoogle Scholar
- Sleeman JP, Cady B, Pantel K: The connectivity of lymphogenous and hematogenous tumor cell dissemination: biological insights and clinical implications. Clin Exp Metastasis. 2012, 29: 737-746. 10.1007/s10585-012-9489-x.View ArticlePubMedGoogle Scholar
- Jatoi I, Hilsenbeck SG, Clark GM, Osborne CK: Significance of axillary lymph node metastasis in primary breast cancer. J Clin Oncol. 1999, 17: 2334-PubMedGoogle Scholar
- Sleeman JP, Christofori G, Fodde R, Collard JG, Berx G, Decraene C, Ruegg C: Concepts of metastasis in flux: the stromal progression model. Semin Cancer Biol. 2012, 22: 174-186. 10.1016/j.semcancer.2012.02.007.View ArticlePubMedGoogle Scholar
- El-Haibi CP, Bell GW, Zhang J, Collmann AY, Wood D, Scherber CM, Csizmadia E, Mariani O, Zhu C, Campagne A: Critical role for lysyl oxidase in mesenchymal stem cell-driven breast cancer malignancy. Proc Natl Acad Sci USA. 2012, 109: 17460-17465. 10.1073/pnas.1206653109.View ArticlePubMedPubMed CentralGoogle Scholar
- Mani SA, Guo W, Liao MJ, Eaton EN, Ayyanan A, Zhou AY, Brooks M, Reinhard F, Zhang CC, Shipitsin M: The epithelial-mesenchymal transition generates cells with properties of stem cells. Cell. 2008, 133: 704-715. 10.1016/j.cell.2008.03.027.View ArticlePubMedPubMed CentralGoogle Scholar
- Polyak K, Weinberg RA: Transitions between epithelial and mesenchymal states: acquisition of malignant and stem cell traits. Nat Rev Cancer. 2009, 9: 265-273. 10.1038/nrc2620.View ArticlePubMedGoogle Scholar
- Thiery JP, Sleeman JP: Complex networks orchestrate epithelial-mesenchymal transitions. Nat Rev Mol Cell Biol. 2006, 7: 131-142. 10.1038/nrm1835.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Książkiewicz M, Markiewicz A, Żaczek AJ: Epithelial-mesenchymal transition: a hallmark in metastasis formation linking circulating tumor cells and cancer stem cells. Pathobiology. 2012, 79: 195-208. 10.1159/000337106.View ArticlePubMedGoogle Scholar
- van Nes J, de Kruijf E, Putter H, Faratian D, Munro A, Campbell F, Smit V, Liefers G-J, Kuppen P, van de Velde C, Bartlett J: Co-expression of SNAIL and TWIST determines prognosis in estrogen receptor–positive early breast cancer patients. Breast Cancer Res Treat. 2012, 133: 49-59. 10.1007/s10549-011-1684-y.View ArticlePubMedGoogle Scholar
- Elloul S, Bukholt Elstrand M, Nesland JM, Tropé CG, Kvalheim G, Goldberg I, Reich R, Davidson B: Snail, Slug, and Smad-interacting protein 1 as novel parameters of disease aggressiveness in metastatic ovarian and breast carcinoma. Cancer. 2005, 103: 1631-1643. 10.1002/cncr.20946.View ArticlePubMedGoogle Scholar
- Martin T, Goyal A, Watkins G, Jiang W: Expression of the transcription factors snail, slug, and twist and their clinical significance in human breast cancer. Ann Surg Oncol. 2005, 12: 488-496. 10.1245/ASO.2005.04.010.View ArticlePubMedGoogle Scholar
- Saxena M, Stephens MA, Pathak H, Rangarajan A: Transcription factors that mediate epithelial-mesenchymal transition lead to multidrug resistance by upregulating ABC transporters. Cell Death and Dis. 2011, 2: e179-10.1038/cddis.2011.61.View ArticleGoogle Scholar
- Al-Hajj M, Wicha MS, Benito-Hernandez A, Morrison SJ, Clarke MF: Prospective identification of tumorigenic breast cancer cells. Proc Natl Acad Sci USA. 2003, 100: 3983-3988. 10.1073/pnas.0530291100.View ArticlePubMedPubMed CentralGoogle Scholar
- Ding L, Ellis MJ, Li S, Larson DE, Chen K, Wallis JW, Harris CC, McLellan MD, Fulton RS, Fulton LL: Genome remodelling in a basal-like breast cancer metastasis and xenograft. Nature. 2010, 464: 999-1005. 10.1038/nature08989.View ArticlePubMedPubMed CentralGoogle Scholar
- Grunewald T, Herbst S, Heinze J, Burdach S: Understanding tumor heterogeneity as functional compartments - superorganisms revisited. J Transl Med. 2011, 9: 79-10.1186/1479-5876-9-79.View ArticlePubMedPubMed CentralGoogle Scholar
- De Wever O, Pauwels P, De Craene B, Sabbah M, Emami S, Redeuilh G, Gespach C, Bracke M, Berx G: Molecular and pathological signatures of epithelial-mesenchymal transitions at the cancer invasion front. Histochem Cell Biol. 2008, 130: 481-494. 10.1007/s00418-008-0464-1.View ArticlePubMedPubMed CentralGoogle Scholar
- Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, Mueller R, Nolan T, Pfaffl MW, Shipley GL: The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 2009, 55: 611-622. 10.1373/clinchem.2008.112797.View ArticlePubMedGoogle Scholar
- Supernat A, Markiewicz A, Welnicka-Jaskiewicz M, Seroczynska B, Skokowski J, Sejda A, Szade J, Czapiewski P, Biernat W, Zaczek A: CD73 Expression as a Potential Marker of Good Prognosis in Breast Carcinoma. Appl Immunohistochem Mol Morphol. 2012, 20: 103-107. 10.1097/PAI.0b013e3182311d82. 110.1097/PAI.1090b1013e3182311d3182382.View ArticlePubMedGoogle Scholar
- Prasad C, Rath G, Mathur S, Bhatnagar D, Parshad R, Ralhan R: Expression analysis of E-cadherin, slug and GSK3beta in invasive ductal carcinoma of breast. BMC Cancer. 2009, 9: 325-10.1186/1471-2407-9-325.View ArticlePubMedPubMed CentralGoogle Scholar
- Kim MA, Lee HS, Lee HE, Kim JH, Yang H-K, Kim WH: Prognostic importance of epithelial–mesenchymal transition-related protein expression in gastric carcinoma. Histopathology. 2009, 54: 442-451. 10.1111/j.1365-2559.2009.03247.x.View ArticlePubMedGoogle Scholar
- Montserrat N, Gallardo A, Escuin D, Catasus L, Prat J, Gutiérrez-Avignó FJ, Peiró G, Barnadas A, Lerma E: Repression of E-cadherin by SNAIL, ZEB1, and TWIST in invasive ductal carcinomas of the breast: a cooperative effort?. Hum Pathol. 2011, 42: 103-110. 10.1016/j.humpath.2010.05.019.View ArticlePubMedGoogle Scholar
- Landis JR, Koch GG: The measurement of observer agreement for categorical data. Biometrics. 1977, 33: 159-174. 10.2307/2529310.View ArticlePubMedGoogle Scholar
- Altman DG, McShane LM, Sauerbrei W, Taube SE: Reporting recommendations for tumor marker prognostic studies (REMARK): explanation and elaboration. BMC Med. 2012, 10: 51-10.1186/1741-7015-10-51.View ArticlePubMedPubMed CentralGoogle Scholar
- Halstedt WS: The Results of radical operations for the cure of carcinoma of the breast. Ann Surg. 1907, 46: 1-19. 10.1097/00000658-190707000-00001.View ArticleGoogle Scholar
- Amir E, Miller N, Geddie W, Freedman O, Kassam F, Simmons C, Oldfield M, Dranitsaris G, Tomlinson G, Laupacis A: Prospective study evaluating the impact of tissue confirmation of metastatic disease in patients with breast cancer. J Clin Oncol. 2012, 30: 587-592. 10.1200/JCO.2010.33.5232.View ArticlePubMedGoogle Scholar
- Aitken SJ, Thomas JS, Langdon SP, Harrison DJ, Faratian D: Quantitative analysis of changes in ER, PR and HER2 expression in primary breast cancer and paired nodal metastases. Ann Oncol. 2010, 21: 1254-1261. 10.1093/annonc/mdp427.View ArticlePubMedGoogle Scholar
- Jensen JD, Knoop A, Ewertz M, Laenkholm AV: ER, HER2, and TOP2A expression in primary tumor, synchronous axillary nodes, and asynchronous metastases in breast cancer. Breast Cancer Res Treat. 2012, 132: 511-521. 10.1007/s10549-011-1610-3.View ArticlePubMedGoogle Scholar
- Falck AK, Ferno M, Bendahl PO, Ryden L: Does analysis of biomarkers in tumor cells in lymph node metastases give additional prognostic information in primary breast cancer?. World J Surg. 2010, 34: 1434-1441. 10.1007/s00268-010-0499-z.View ArticlePubMedGoogle Scholar
- Niikura N, Liu J, Hayashi N, Mittendorf EA, Gong Y, Palla SL, Tokuda Y, Gonzalez-Angulo AM, Hortobagyi GN, Ueno NT: Loss of human epidermal growth factor receptor 2 (HER2) expression in metastatic sites of HER2-overexpressing primary breast tumors. J Clin Oncol. 2012, 30: 593-599. 10.1200/JCO.2010.33.8889.View ArticlePubMedPubMed CentralGoogle Scholar
- Thompson A, Jordan L, Quinlan P, Anderson E, Skene A, Dewar J, Purdie C, Group tBRiTS: Prospective comparison of switches in biomarker status between primary and recurrent breast cancer: the Breast Recurrence In Tissues Study (BRITS). Breast Cancer Res. 2010, 12: R92-10.1186/bcr2771.View ArticlePubMedPubMed CentralGoogle Scholar
- Lindström LS, Karlsson E, Wilking UM, Johansson U, Hartman J, Lidbrink EK, Hatschek T, Skoog L, Bergh J: Clinically used breast cancer markers such as estrogen receptor, progesterone receptor, and human epidermal growth factor receptor 2 are unstable throughout tumor progression. J Clin Oncol. 2012, 30: 2601-2608. 10.1200/JCO.2011.37.2482.View ArticlePubMedGoogle Scholar
- Liedtke C, Broglio K, Moulder S, Hsu L, Kau S-W, Symmans WF, Albarracin C, Meric-Bernstam F, Woodward W, Theriault RL: Prognostic impact of discordance between triple-receptor measurements in primary and recurrent breast cancer. Ann Oncol. 2009, 20: 1953-1958. 10.1093/annonc/mdp263.View ArticlePubMedPubMed CentralGoogle Scholar
- Macfarlane R, Seal M, Speers C, Woods R, Masoudi H, Aparicio S, Chia SK: Molecular alterations between the primary breast cancer and the subsequent locoregional/metastatic tumor. Oncologist. 2012, 17: 172-178. 10.1634/theoncologist.2011-0127.View ArticlePubMedPubMed CentralGoogle Scholar
- D'Andrea MR, Limiti MR, Bari M, Zambenedetti P, Montagutti A, Ricci F, Pappagallo GL, Sartori D, Vinante O, Mingazzini PL: Correlation between genetic and biological aspects in primary non-metastatic breast cancers and corresponding synchronous axillary lymph node metastasis. Breast Cancer Res Treat. 2007, 101: 279-284. 10.1007/s10549-006-9300-2.View ArticlePubMedGoogle Scholar
- Houssami N, Macaskill P, Balleine RL, Bilous M, Pegram MD: HER2 discordance between primary breast cancer and its paired metastasis: tumor biology or test artefact? Insights through meta-analysis. Breast Cancer Res Treat. 2011, 129: 659-674. 10.1007/s10549-011-1632-x.View ArticlePubMedGoogle Scholar
- Simmons C, Miller N, Geddie W, Gianfelice D, Oldfield M, Dranitsaris G, Clemons MJ: Does confirmatory tumor biopsy alter the management of breast cancer patients with distant metastases?. Ann Oncol. 2009, 20: 1499-1504. 10.1093/annonc/mdp028.View ArticlePubMedPubMed CentralGoogle Scholar
- Dawood S, Gonzalez-Angulo AM: To biopsy or not to biopsy: is that the only question?. Oncologist. 2012, 17: 151-153. 10.1634/theoncologist.2011-0432.View ArticlePubMedPubMed CentralGoogle Scholar
- Nishimura R, Osako T, Okumura Y, Tashima R, Toyozumi Y, Arima N: Changes in the ER, PgR, HER2, p53 and Ki-67 biological markers between primary and recurrent breast cancer: discordance rates and prognosis. World J Surg Oncol. 2011, 9: 131-10.1186/1477-7819-9-131.View ArticlePubMedPubMed CentralGoogle Scholar
- Moody SE, Perez D, Pan T-c, Sarkisian CJ, Portocarrero CP, Sterner CJ, Notorfrancesco KL, Cardiff RD, Chodosh LA: The transcriptional repressor Snail promotes mammary tumor recurrence. Cancer Cell. 2005, 8: 197-209. 10.1016/j.ccr.2005.07.009.View ArticlePubMedGoogle Scholar
- Toyama T, Zhang Z, Iwase H, Yamashita H, Ando Y, Hamaguchi M, Mizutani M, Kondo N, Fujita T, Fujii Y, Iwata H: Low expression of the snail gene is a good prognostic factor in node-negative invasive ductal carcinomas. Jpn J Clin Oncol. 2006, 36: 357-363. 10.1093/jjco/hyl038.View ArticlePubMedGoogle Scholar
- Soini Y, Tuhkanen H, Sironen R, Virtanen I, Kataja V, Auvinen P, Mannermaa A, Kosma V-M: Transcription factors zeb1, twist and snai1 in breast carcinoma. BMC Cancer. 2011, 11: 73-10.1186/1471-2407-11-73.View ArticlePubMedPubMed CentralGoogle Scholar
- Kowalski PJ, Rubin MA, Kleer CG: E-cadherin expression in primary carcinomas of the breast and its distant metastases. Breast Cancer Res. 2003, 5: R217-222. 10.1186/bcr651.View ArticlePubMedPubMed CentralGoogle Scholar
- Kleer CG, van Golen KL, Braun T, Merajver SD: Persistent E-cadherin expression in inflammatory breast cancer. Mod Pathol, Inc. 2001, 14: 458-464. 10.1038/modpathol.3880334.View ArticleGoogle Scholar
- Kokkinos MI, Wafai R, Wong MK, Newgreen DF, Thompson EW, Waltham M: Vimentin and epithelial-mesenchymal transition in human breast cancer–observations in vitro and in vivo. Cells Tissues Organs. 2007, 185: 191-203. 10.1159/000101320.View ArticlePubMedGoogle Scholar
- Gilles C, Polette M, Piette J, Delvigne AC, Thompson EW, Foidart JM, Birembaut P: Vimentin expression in cervical carcinomas: association with invasive and migratory potential. J Pathol. 1996, 180: 175-180. 10.1002/(SICI)1096-9896(199610)180:2<175::AID-PATH630>3.0.CO;2-G.View ArticlePubMedGoogle Scholar
- Bonnomet A, Syne L, Brysse A, Feyereisen E, Thompson EW, Noel A, Foidart JM, Birembaut P, Polette M, Gilles C: A dynamic in vivo model of epithelial-to-mesenchymal transitions in circulating tumor cells and metastases of breast cancer. Oncogene. 2012, 31: 3741-3753. 10.1038/onc.2011.540.View ArticlePubMedGoogle Scholar
- Nowell PC: The clonal evolution of tumor cell populations. Science. 1976, 194: 23-28. 10.1126/science.959840.View ArticlePubMedGoogle Scholar
- Greaves M, Maley CC: Clonal evolution in cancer. Nature. 2012, 481: 306-313. 10.1038/nature10762.View ArticlePubMedPubMed CentralGoogle Scholar
- Gerlinger M, Rowan AJ, Horswell S, Larkin J, Endesfelder D, Gronroos E, Martinez P, Matthews N, Stewart A, Tarpey P: Intratumor heterogeneity and branched evolution revealed by multiregion sequencing. N Engl J Med. 2012, 366: 883-892. 10.1056/NEJMoa1113205.View ArticlePubMedPubMed CentralGoogle Scholar
- Shah SP, Roth A, Goya R, Oloumi A, Ha G, Zhao Y, Turashvili G, Ding J, Tse K, Haffari G: The clonal and mutational evolution spectrum of primary triple-negative breast cancers. Nature. 2012, 486: 395-399.PubMedGoogle Scholar
- Klein CA: Parallel progression of primary tumours and metastases. Nat Rev Cancer. 2009, 9: 302-312. 10.1038/nrc2627.View ArticlePubMedGoogle Scholar
- Pusztai L, Viale G, Kelly CM, Hudis CA: Estrogen and HER-2 receptor discordance between primary breast cancer and metastasis. Oncologist. 2010, 15: 1164-1168. 10.1634/theoncologist.2010-0059.View ArticlePubMedPubMed CentralGoogle Scholar
- Engel KB, Moore HM: Effects of preanalytical variables on the detection of proteins by immunohistochemistry in formalin-fixed, paraffin-embedded tissue. Arch Pathol Lab Med. 2011, 135: 537-543.PubMedGoogle Scholar
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2 − ΔΔCT method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.View ArticlePubMedGoogle 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.