International Journal of Surgical OncologyVolume 2012 (2012), Article ID 984346, 9 pagesdoi:10.1155/2012/984346

Research Article

Cell Polarity, Epithelial-Mesenchymal Transition, and Cell-Fate Decision Gene Expression in Ductal Carcinoma In Situ

Danila Coradini,1,2 Patrizia Boracchi,¹ Federico Ambrogi,¹ Elia Biganzoli,¹ and Saro Oriana²

¹Department of Work Medicine “Clinica del Lavoro L. Devoto”, Section of Medical Statistics and Biometry “G.A. Maccacaro”, University of Milano, 20133 Milan, Italy²Senology Center, Casa di Cura Ambrosiana, Istituto Sacra Famiglia, Cesano Boscone, 20090 Milano, Italy

Received 29 November 2011; Revised 17 January 2012; Accepted 25 January 2012

Academic Editor: Lucio Fortunato

Copyright © 2012 Danila Coradini et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Loss of epithelial cell identity and acquisition of mesenchymal features are early events in the neoplastic transformation of mammary cells. We investigated the pattern of expression of a selected panel of genes associated with cell polarity and apical junction complex or involved in TGF-β-mediated epithelial-mesenchymal transition and cell-fate decision in a series of DCIS and corresponding patient-matched normal tissue. Additionally, we compared DCIS gene profile with that of atypical ductal hyperplasia (ADH) from the same patient. Statistical analysis identified a “core” of genes differentially expressed in both precursors with respect to the corresponding normal tissue mainly associated with a terminally differentiated luminal estrogen-dependent phenotype, in agreement with the model according to which ER-positive invasive breast cancer derives from ER-positive progenitor cells, and with an autocrine production of estrogens through androgens conversion. Although preliminary, present findings provide transcriptomic confirmation that, at least for the panel of genes considered in present study, ADH and DCIS are part of a tumorigenic multistep process and strongly arise the necessity for the regulation, maybe using aromatase inhibitors, of the intratumoral and/or circulating concentration of biologically active androgens in DCIS patients to timely hamper abnormal estrogens production and block estrogen-induced cell proliferation.

1. Introduction

Ductal carcinoma in situ (DCIS), also known as intraductal carcinoma, is the most common type of noninvasive breast cancer in women [1]. From the mid-1970s, the incidence of DCIS has sharply increased, primarily because of the adoption of radiographic screening for invasive carcinoma. Currently, it accounts for approximately 25% of newly diagnosed breast cancer cases [1].

DCIS is a nonobligate precursor to invasive breast cancer, and for this reason some members of the 2009 US National Institutes of Health DCIS consensus conference proposed to remove the word “carcinoma” from the term DCIS [2]. Nevertheless, experimental studies have shown the presence of carcinoma precursor cells in DCIS lesions [3–5], and clinical evidence indicates that approximately 50% of cases will progress to invasive breast cancer if untreated [6, 7]. As a result, the malignant nature of DCIS remains debated, primarily because of the limited knowledge of DCIS arising and development. In fact, while several genome-wide studies have compared the gene profiles of DCIS and invasive breast cancer, very few studies have investigated and recognized the molecular alteration that characterize DCIS with respect to normal tissue [8–11].

DCIS is defined as an abnormal proliferation of transformed mammary epithelial cells within the closed environment of a duct, likely in response to microenvironment alterations including hypoxia and nutrient deprivation [4, 12]. Among the processes early affected during mammary cells transformation, those involved in the establishment and maintenance of epithelial cell identity and tissue specificity are of particular relevance. In fact, epithelial mammary cells are characterized by an asymmetric distribution of cytoplasmic and membrane proteins, termed apicobasolateral cell polarity, essential for a correct cell-cell adhesion and the formation of an epithelial sheet. As a result of an epithelial to mesenchymal transition (EMT), during neoplastic transformation, cell polarity and epithelial morphology are early lost polarized and immotile epithelial cells acquire a fibroblast-like morphology and increased cell motility [13].

EMT process can be induced by a variety of signaling pathways among which the main and best characterized is that involving transforming-growth factor-β (TGF-β). TGF-β is a multifunctional cytokine and a powerful tumor suppressor that governs many aspects of mammary epithelial cells physiology and homeostasis [14]. Under abnormal microenvironment conditions; however, some mammary epithelial cells may acquire resistance to TGF-β, circumvent its cytostatic effect and tumor suppressive activity, and activate EMT [15, 16].

Recent studies have demonstrated that EMT may generate cells with stemness-like properties, especially in the transitioning mammary epithelial cell compartment [17, 18]. Therefore, a interrelationship among EMT (TGF-β mediated), disruption of the mechanisms deputed to cell polarity and adhesion control, and acquisition of stemness-like features can be assumed already in DCIS [19, 20].

Taking advantage from the only microarray dataset publicly accessible at the ArrayExpress web site, we investigated the pattern of expression of a selected panel of genes involved in TGF-β-mediated EMT or associated with epithelial cells identity (i.e., cell polarity and apical junction complex) and cell-fate decision in a series of DCIS and corresponding patient-matched histologically normal (HN) epithelium [11]. As the whole-gene expression profile of patient-matched atypical ductal hyperplasia (ADH) was also available, we further compared DCIS and ADH profiles to verify the hypothesis according to which breast cancer progression is a multistep process involving a continuum of changes from normal phenotype through hyperplastic lesions, carcinoma in situ, and invasive carcinoma [21].

2. Materials and Methods

2.1. Materials

As reported in the original paper [11], patient-matched samples (HN, ADH, and DCIS) were isolated via laser capture microdissection from surgical specimens of 12 preoperative untreated patients with an ER-positive (immunohistochemically evaluated) sporadic breast cancer. Gene expression was determined by using the Affymetrix Human Genome HG-U133A GeneChip; corresponding microarray dataset was publicly available at the ArrayExpress web site (http://www.ebi.ac.uk/arrayexpress/) with the Accession number E-GEOD-16873.

2.2. Gene Selection

To select the panel of genes specifically involved in the TGF-β-activated EMT, cell polarity and apical junction complex, and cell-fate decision, we combined Gene Ontology (http://www.geneontology.org) and PubMed (http://www.ncbi.nlm.nih.gov) information. In addition, since cancer cells often increase their autocrine production of TGF-β to activate angiogenesis in response to oxygen and nutrients deprivation [14], we also included some genes coding for angiogenesis-inducing factors. On the whole, a set of 199 genes evolutionarily conserved in Homo sapiens was established (see Supplementary Table  1 available on line at doi 10.1155/2012/984346). However, because 27 genes had no corresponding probe-sets on the HG-U133A GeneChip, the gene set was actually composed of 172 elements, 33 of which involved in EMT activation [13–15, 21–23], 75 involved in cell polarity and apical junction assembly [24–26], 36 involved in cell fate-decisions and in the maintenance of a self-renewal state in tumorigenic adult tissues [27–29], 28 involved in hormone steroid signaling [30–32], angiogenesis activation [33, 34] or used as luminal and basal markers [35–38]. These 172 genes corresponded to 339 Affymetrix probe-sets, as verified by GeneAnnot system v2.0 (http://bioinfo2.weizmann.ac.il/geneannot/), that additionally provided us information about the quality of each probe-set in terms of sensitivity and specificity score [39] (see Supplementary).

2.3. Statistical Analysis

As some genes are recognized by more than a single probe set, each of which characterized by an individual specificity and sensitivity that differently contribute to gene expression value, a gene expression mean value was calculated after weighting each probe-set for its own sensitivity and specificity score. Specifically, each expression value (already log 2 transformed in the original dataset) was multiplied for the semi sum of sensitivity and specificity scores of the corresponding probe set. Given the patient-matched samples study design, all statistical analyses were performed considering a regression model for repeated measures with random effect, and the differential gene expression was evaluated by t-test on regression coefficients. To correct for multiple testing, the false discovery rate (FDR) was used [40].

To evidence latent variables accounting for genes correlations, a factor analysis was applied [41] in the following three comparisons: DCIS and paired ADH, DCIS and paired HN, and ADH and paired HN. The number of retained factors was selected according to the scree test [42]. To facilitate the interpretation of the factors, varimax rotation was applied. Loading values lower than 0.3 were not considered.

All analyses were performed using open source software R 2.11.1 packages HDMD (http://www.R-project.org).

3. Results and Discussion

Genes found differentially expressed () between DCIS and NH or ADH and NH are reported in Table 1. Specifically, 47 of the 172 selected genes were found differentially expressed between DCIS and NH (11 with an estimated ) and 28 were found differentially expressed between ADH and NH (only one with an estimated ). Notably, 24 of the 28 genes found differentially expressed between ADH and NH were found differentially expressed (in a similar manner) also between DCIS and NH. The persistence of this “core” of genes, dysregulated in a similar manner in both invasive breast cancer precursors, seems to support the hypothesis of ADH as the direct precursor of DCIS. In fact, in agreement with the proposed multistep process, DCIS showed an increased number of genes differentially expressed with respect to ADH.

Table 1: Differentially expressed genes between ADH or DCIS and histologically normal (HN) tissue (ordered according to value).

ADH versus HN
DCIS versus HN

Gene symbol
value
Variation
Gene symbol
value
Variation

JAM3
0.000045
↓
JAM3
0.000016
↓

JAG2
0.000337
↓
EGFR
0.000103
↓

CD24
0.000364
↑
SNAI2
0.000127
↓

SNAI2
0.000594
↓
CLDN5
0.000136
↓

EGFR
0.001897
↓
JAM2
0.000150
↓

FOXC1
0.002222
↓
FOXC1
0.000167
↓

EGF
0.003163
↓
CD24
0.000277
↑

ID2
0.003547
↑
JAG2
0.000297
↓

JAM2
0.005052
↓
KRT5
0.000420
↓

TJP3
0.005070
↑
KRT14
0.000425
↓

TGFBR3
0.005537
↓
GATA3
0.000490
↑

KRT17
0.005624
↓
TGFBR3
0.000710
↓

GATA3
0.005795
↑
TJP3
0.001046
↑

TP53
0.006737
↑
KRT17
0.001402
↓

CDH4
0.006916
↓
SOX4
0.001942
↑

AKT1
0.009044
↑
AKT1
0.002867
↑

CLDN7
0.012026
↑
CLDN8
0.003090
↓

EPCAM
0.015207
↑
CDC42
0.003643
↑

ABCG2
0.019210
↓
EGF
0.003752
↓

KRT5
0.019854
↓
EPCAM
0.004605
↑

KRT14
0.020117
↓
CLDN7
0.004749
↑

CLDN11
0.023005
↓
ESR1
0.007188
↑

PARD3
0.030202
↓
DLG1
0.007249
↑

SOX4
0.036127
↑
KRT19
0.007637
↑

CDC42
0.036404
↑
FOXA1
0.007844
↑

TIAM1
0.038034
↑
KRT18
0.010091
↑

ESR1
0.038527
↑
CDH3
0.010638
↓

PVR
0.039104
↓
TIAM1
0.010722
↑

TGFBR2
0.013243
↓

RHOA
0.013540
↑

BRCA1
0.014098
↑

ID4
0.017030
↓

ID2
0.018115
↑

NOTCH4
0.018550
↓

PVRL2
0.019914
↑

AKT3
0.020275
↓

ACTN1
0.020619
↓

PROM1
0.022238
↓

CDH4
0.022933
↓

DLG3
0.023798
↑

F11R
0.031657
↑

CTNNA1
0.033234
↑

MTA2
0.033999
↓

ABCG2
0.037966
↓

MPP5
0.039023
↑

HIF1A
0.039489
↑

CDKN1A
0.047769
↑

In bold, genes with an estimated FDR < 0.01.

With respect to normal tissue, both DCIS and ADH showed the overexpression of ESR1, coding for the estrogen receptor; CD24, coding for a mucin-like cell-adhesion molecule positively associated with a terminally differentiated luminal phenotype [43, 44]; GATA3, coding for a transcription factor involved in mammary gland morphogenesis [35, 36]; CLDN7, EPCAM, and TJP3, coding for tight junction components; CDC42 and TIAM1, coding for two small GTPase family members involved in cell polarity and apical junction complex formation. Concomitantly, both breast cancer precursors showed the underexpression of EGFR gene, in which expression is generally negatively associated with ESR1 expression, KRT5, KRT14, and KRT17, all coding for cytokeratins associated with a basal phenotype. On the whole, this pattern of expression clearly indicates that DCIS and ADH are both characterized by a terminally differentiated luminal phenotype. Since all specimens were derived from patients with an ER-positive ductal carcinoma, it is conceivable the hypothesis that the establishment of an estrogen-dependent phenotype, in response to estrogens present in the microenvironment, should be a very early event in the tumorigenic process. Such a finding is in agreement with the model for breast cancer development proposed by Dontu et al. [45], according to which ER-positive cancers should derive from transiently amplifying ERpositive progenitor cells. Escaped from proliferation control as a consequence of genetic and epigenetic alterations in genes involved in cell-fate decision, these ER-positive progenitor cells should generate cells constitutively expressing estrogen receptor. Once established in ADH, this terminally differentiated luminal phenotype seems to consolidate in DCIS as demonstrated by the presence, among the genes exclusively expressed in DCIS, of KRT18, KRT19 (coding for some luminal-associated cytokeratins) and DLG1, DLG3, and MPP5 (coding for some cell polarity complex components).

With respect to both tumor precursors, histologically normal tissue expressed genes coding for some transcription factors involved in EMT (SNAI2 and TGFBR3) and cell-fate decision (FOXC1 and JAG2), and for stemness-associated features (ABCG2). This finding is more evident considering the 47 genes differentially expressed between DCIS and normal tissue among which we found some other genes involved in the TGF-β-mediated EMT (AKT3, ID4 and TGFBR2), and cell-fate decision and self-renewal (NOTCH4 and PROM1).

Such an apparently paradoxical finding, that is, the positive association between histologically normal tissues and EMT- and stemness-related genes (more likely expected in transformed tissues) is not surprising since we observed a similar behavior in an independent cases series composed of primary breast cancers and corresponding patient-matched normal tissue (submitted paper) and in normal pleura with respect to pleural mesothelioma [46]. In agreement with the physiological remodeling of the mammary gland, such a finding should support the persistence of resident stem/progenitor cells in normal tissue [47, 48].

Factor analysis, applied to investigate the latent variables intrinsically associated with the selected 172 genes, corroborated these findings and highlighted some other interesting interrelations. In agreement with the results provided by t-paired test, factor analysis indicated that in both precursors (Figure 1(a) for DCIS and Figure 2(a) for ADH), the first factor (F1) was principally characterized by genes associated with an estrogen-dependent epithelial phenotype. In fact, within the genes with a positive loading value on F1, we found those coding for hormone steroid receptors (AR and ESR1) and transcriptional coactivator for steroid and nuclear hormone receptors (NCOA1), for tight (EPCAM, MAGI1, MAGI2, MPDZ, TJP1, TJP2 and TJP3) and adherens (CTNND1, PFN1 and PVRL2) junction components, and for small GTPase family members involved in epithelial cell polarization processes (CDC42 and RHOA). In addition, we found some genes associated with cell-fate decision: BMI-1, coding for a member of Polycomb group required to maintain the transcriptionally repressive state of many genes throughout embryo development [49] and adult tissues differentiation including mammary gland [50]; FDZ4 and FDZ6, two members of the frizzled gene family involved in β-catenin-signaling transduction and intercellular transmission of polarity information in differentiated tissues [51].

(a)

(b)

Figure 1: Factor analysis in DCIS subgroup. Schematic representation of genes with a loading value < characterizing the first (F1) and the second factor (F2). Solid color indicates a positive loading value whereas dashed color indicates a negative loading value. Color correspondence: light yellow, tight junction components; dark yellow, adherens junction components; light blue, polarity complexes components; dark blue, angiogenesis; orange, cell-fate decision; light green, luminal markers and hormone steroid; dark green, basal markers; pink, epithelial-mesenchymal transition, grey, GTPase family members.

(a)

(b)

Figure 2: Factor analysis in ADH subgroup. Schematic representation of genes with a loading value < characterizing the first (F1) and the second factor (F2). Solid color indicates a positive loading value whereas dashed color indicates a negative loading value. Color correspondence: light yellow, tight junction components; dark yellow, adherens junction components; light blue, polarity complexes components; dark blue, angiogenesis; orange, cell-fate decision; light green, luminal markers and hormone steroid; dark green, basal markers; pink, epithelial-mesenchymal transition, grey, GTPase family members.

Furthermore, in agreement with the notion that ERα, and TGF-β-signaling pathways are major regulators during mammary gland development [52], we found genes coding for TGF-β receptor (TGFBR3) or for proteins involved in canonical (ID1, SMAD2 and SMAD4) and noncanonical (PTEN, RHOA and ROCK1) TGF-β pathways [15, 16].

Finally, probably associated with the adaptation of DCIS and ADH to hypoxic stress caused by the unbalance between cell proliferation and oxygen supply, we also found JAM2, JAM3, VEGFA, VEGFB, and VEGFC, all genes involved in endothelial cell proliferation.

The second factor (F2) identified by factor analysis (Figure 1(b) for DCIS and Figure 2(b) for ADH) was conversely characterized by the positive loading value of several genes coding for mesenchymal markers (EGF, EGFR, KRT5, KRT6B, KRT14 and KRT17), and the negative loading value of GATA3, the gene coding for the transcription factor driving the luminal morphogenesis of the mammary gland [35, 36]. In addition, F2 was characterized by the presence of some genes coding for proteins playing a critical role in cell-fate decision and cell-renewal (ALDH1A3, FOXC1, NOTCH2, PROM1, and SOX9) [53].

Notably, when applied to ADH subgroup, factor analysis identified CYP19A1 as included in the panel of genes characterizing the first factor. This gene encodes for cytochrome P450 (better known as aromatase), the enzyme that catalyzes the conversion from circulating and rostenedione to estrone or testosterone to estradiol, and its presence should support the hypothesis of a very early activation of an autocrine production of estrogen. Furthermore, the concomitant presence of CYP19A1, AR, ESR1, and NCOA1 should provide a transcriptomic confirmation for the clinical evidence that high androgens level may have detrimental effect on breast carcinogenesis and progression due to a persistent local estrogen production that incessantly stimulates epithelial cell proliferation [54–56].

When applied to DCIS subgroup, factor analysis seems to indicate a consolidation of such an estrogen dependence as suggested by the additional presence of PGR, the gene coding for progesterone receptor and which expression is under estrogenic control.

4. Conclusions

Elucidating the initial steps of breast tumorigenesis is of paramount importance to allow an even early diagnosis and consequently an adequate treatment strategy aimed to prevent the malignant transformation of preneoplastic alterations. That is of particular importance for DCIS because of its high incidence [1] and facility in progressing to invasive breast cancer if untreated [6, 7].

Experimental evidence till now accumulated has clearly indicated that at the basis of the neoplastic transformation of mammary epithelial cells there are loss of apicobasal epithelial cell identity and acquisition of a functional mesenchymal morphology. Therefore, we investigated the pattern of expression of a selected panel of genes associated with epithelial cells identity (i.e., cell polarity and apical junction complex) or involved in TGF-β-mediated EMT and cell-fate decision in a series of DCIS and corresponding patient-matched normal tissue. In addition, we compared DCIS profile with that of patient-matched ADH to investigate the hypothesis according to which breast cancer progression is a multistep process involving a continuum of changes from normal phenotype through hyperplastic lesions, carcinoma in situ, and invasive carcinoma [21].

Statistical analysis seems to support this hypothesis because it identified a “core” of genes, mainly associated with a terminally differentiated luminal phenotype, and differentially expressed in both precursors with respect to the corresponding normal tissue. Notably, these alterations in gene expression did not result in a progressive mesenchymal transition but rather in a terminally differentiated luminal phenotype, in agreement with the model according to which ER-positive invasive breast cancer derives from ER-positive progenitor cells [45]. The constitutive expression of ER should make ADH and DCIS forming cells able to exploit the proliferative stimulus induced by estrogens whereas the establishment of an autocrine production of estrogens, through androgens conversion, should provide an additional selective advantage. The detrimental effect of such continuous estrogen stimulation should be corroborated by the observation that all patients included in the present study developed an invasive ER-positive ductal carcinoma.

Experimental evidence supporting the hypothesis that androgens conversion may be involved in DCIS development, and progression has been provided by a recent study in which the intratumoral concentration of estradiol and 5α-dihydrotestosterone (DHT), and the expression of some sex steroid-producing enzymes, including aromatase, has been evaluated in DCIS specimens [57]. The study clearly demonstrated that aromatase expression level was significantly higher in DCIS with respect to nonneoplastic tissue suggesting the self-sustaining process adopted by DCIS.

Taken together all these findings provide transcriptomic confirmation that, at least for the panel of genes considered in present study, ADH and DCIS are part of a tumorigenic multistep process and strongly arise the necessity for the regulation, maybe using aromatase inhibitors, of the intratumoral and/or circulating concentration of biologically active androgens in DCIS patients to timely hamper abnormal estrogens production and block estrogen-induced cell proliferation [58].

There is no doubt that present in silico study suffers for the limitation common to the majority of studies involving gene expression profile, that is, the lack of validation, at protein level, of the modulations observed at mRNA level. In fact, it is well known that mRNA transcript levels do not always reflect protein expression. However, the immunohistochemical data provided in The Human Protein Atlas largely confirms the differential gene expression that we observed between histologically normal and cancerous tissue. For instance, when we considered the 11 genes with an estimated (Table 1, in bold), we found that protein expression of normal breast (glandular cell) and breast cancer tissue mainly paralleled the mRNA differential expression we observed in our dataset, making our preliminary findings more reliable and worthy of further investigation.

Author’s Contribution

D. Coradini and P. Boracchi contributed equally to this paper.

References

B. A. Virnig, T. M. Tuttle, T. Shamliyan, and R. L. Kane, “Ductal carcinoma in Situ of the breast: a systematic review of incidence, treatment, and outcomes,” Journal of the National Cancer Institute, vol. 102, no. 3, pp. 170–178, 2010.

C. J. Allegra, D. R. Aberle, P. Ganschow et al., “NIH state-of-the-science conference statement: diagnosis and management of ductal carcinoma in situ (DCIS),” NIH Consensus and State-of-the-Science Statements, vol. 26, no. 2, pp. 1–27, 2009.

N. P. Castro, C. A. B. T. Osório, C. Torres et al., “Evidence that molecular changes in cells occur before morphological alterations during the progression of breast ductal carcinoma,” Breast Cancer Research, vol. 10, no. 5, article no. R87, 2008.

V. Espina, B. D. Mariani, R. I. Gallagher et al., “Malignant precursor cells pre-exist in human breast DCIS and require autophagy for survival,” PLoS One, vol. 5, no. 4, Article ID e10240, 2010.

D. C. Sgroi, “Preinvasive breast cancer,” Annual Review of Pathology, vol. 5, pp. 193–221, 2010.

H. M. Kuerer, C. T. Albarracin, W. T. Yang et al., “Ductal carcinoma in situ: state of the science and roadmap to advance the field,” Journal of Clinical Oncology, vol. 27, no. 2, pp. 279–288, 2009.

L. C. Collins, R. M. Tamimi, H. J. Baer, J. L. Connolly, G. A. Colditz, and S. J. Schnitt, “Outcome of patients with ductal carcinoma in situ untreated after diagnostic biopsy: results from the nurses' health study,” Cancer, vol. 103, no. 9, pp. 1778–1784, 2005.

A. Vincent-Salomon, C. Lucchesi, N. Gruel et al., “Integrated genomic and transcriptomic analysis of ductal carcinoma in situ of the breast,” Clinical Cancer Research, vol. 14, no. 7, pp. 1956–1965, 2008.

R. M. Tamimi, H. J. Baer, J. Marotti et al., “Comparison of molecular phenotypes of ductal carcinoma in situ and invasive breast cancer,” Breast Cancer Research, vol. 10, no. 4, article no. R67, 2008.

D. C. Allred, Y. Wu, S. Mao et al., “Ductal carcinoma in situ and the emergence of diversity during breast cancer evolution,” Clinical Cancer Research, vol. 14, no. 2, pp. 370–378, 2008.

L. A. Emery, A. Tripathi, C. King et al., “Early dysregulation of cell adhesion and extracellular matrix pathways in breast cancer progression,” American Journal of Pathology, vol. 175, no. 3, pp. 1292–1302, 2009.

M. L. Gauthier, H. K. Berman, C. Miller et al., “Abrogated response to cellular stress identifies DCIS associated with subsequent tumor events and defines basal-like breast tumors,” Cancer Cell, vol. 12, no. 5, pp. 479–491, 2007.

E. Tomaskovic-Crook, E. W. Thompson, and J. P. Thiery, “Epithelial to mesenchymal transition and breast cancer,” Breast Cancer Research, vol. 11, no. 6, article no. 213, 2009.

R. Derynck, R. J. Akhurst, and A. Balmain, “TGF-β signaling in tumor suppression and cancer progression,” Nature Genetics, vol. 29, no. 2, pp. 117–129, 2001.

J. Xu, S. Lamouille, and R. Derynck, “TGF-β-induced epithelial to mesenchymal transition,” Cell Research, vol. 19, no. 2, pp. 156–172, 2009.

C. H. Heldin, M. Landström, and A. Moustakas, “Mechanism of TGF-β signaling to growth arrest, apoptosis, and epithelial-mesenchymal transition,” Current Opinion in Cell Biology, vol. 21, no. 2, pp. 166–176, 2009.

S. A. Mani, W. Guo, M. J. Liao et al., “The epithelial-mesenchymal transition generates cells with properties of stem cells,” Cell, vol. 133, no. 4, pp. 704–715, 2008.

A. P. Morel, M. Lièvre, C. Thomas, G. Hinkal, S. Ansieau, and A. Puisieux, “Generation of breast cancer stem cells through epithelial-mesenchymal transition,” PLoS One, vol. 3, no. 8, Article ID e2888, 2008.

C. J. Creighton, J. C. Chang, and J. M. Rosen, “Epithelial-mesenchymal transition (EMT) in tumor-initiating cells and its clinical implications in breast cancer,” Journal of Mammary Gland Biology and Neoplasia, vol. 15, no. 2, pp. 253–260, 2010.

G. Moreno-Bueno, F. Portillo, and A. Cano, “Transcriptional regulation of cell polarity in EMT and cancer,” Oncogene, vol. 27, no. 55, pp. 6958–6969, 2008.

P. T. Simpson, J. S. Reis-Filho, T. Gale, and S. R. Lakhani, “Molecular evolution of breat cancer,” Journal of Pathology, vol. 205, no. 2, pp. 248–254, 2005.

D. S. Micalizzi and H. L. Ford, “Epithelial-mesenchymal transition in development and cancer,” Future Oncology, vol. 5, no. 8, pp. 1129–1143, 2009.

T. Blick, E. Widodo, H. Hugo et al., “Epithelial mesenchymal transition traits in human breast cancer cell lines,” Clinical and Experimental Metastasis, vol. 25, no. 6, pp. 629–642, 2008.

A. Wodarz and I. Näthke, “Cell polarity in development and cancer,” Nature Cell Biology, vol. 9, no. 9, pp. 1016–1024, 2007.

H. A. J. Müller, “Genetic control of epithelial cell polarity: lessons from Drosophila,” Developmental Dynamics, vol. 218, no. 1, pp. 52–67, 2000.

F. Martin-Belmonte and K. Mostov, “Regulation of cell polarity during epithelial morphogenesis,” Current Opinion in Cell Biology, vol. 20, no. 2, pp. 227–234, 2008.

M. Liu, M. C. Casimiro, C. Wang et al., “p21CIP1 attenuates Ras- and c-Myc-dependent breast tumor epithelial mesenchymal transition and cancer stem cell-like gene expression in vivo,” Proceedings of the National Academy of Sciences of the United States of America, vol. 106, no. 45, pp. 19035–19039, 2009.

A. Lugli, G. Iezzi, I. Hostettler et al., “Prognostic impact of the expression of putative cancer stem cell markers CD133, CD166, CD44s, EpCAM, and ALDH1 in colorectal cancer,” British Journal of Cancer, vol. 103, no. 3, pp. 382–390, 2010.

J. C. Lawson, G. L. Blatch, and A. L. Edkins, “Cancer stem cells in breast cancer and metastasis,” Breast Cancer Research and Treatment, vol. 118, no. 2, pp. 241–254, 2009.

B. Manavathi, K. Singh, and R. Kumar, “MTA family of coregulators in nuclear receptor biology and pathology,” Nuclear Receptor Signaling, vol. 5, article e010, 2007.

N. Fujita, M. Kajita, P. Taysavang, and P. A. Wade, “Hormonal regulation of metastasis-associated protein 3 transcription in breast cancer cells,” Molecular Endocrinology, vol. 18, no. 12, pp. 2937–2949, 2004.

J. Xu, R. C. Wu, and B. W. O'Malley, “Normal and cancer-related functions of the p160 steroid receptor co-activator (SRC) family,” Nature Reviews Cancer, vol. 9, no. 9, pp. 615–630, 2009.

A. E. Greijer, P. van der Groep, D. Kemming et al., “Up-regualtion of gene expression by hypoxia is mediated predominantly by hypoxia-inducible factor I (HIF-I),” Journal of Pathology, vol. 206, no. 3, pp. 291–304, 2005.

S. F. Schoppmann, D. Tamandl, L. Roberts et al., “HER2/neu expression correlates with vascular endothelial growth factor-C and lymphangiogenesis in lymph node-positive breast cancer,” Annals of Oncology, vol. 21, no. 5, pp. 955–960, 2009.

H. Kouros-Mehr, J. W. Kim, S. K. Bechis, and Z. Werb, “GATA-3 and the regulation of the mammary luminal cell fate,” Current Opinion in Cell Biology, vol. 20, no. 2, pp. 164–170, 2008.

M. L. Asselin-Labat, K. D. Sutherland, H. Barker et al., “Gata-3 is an essential regulator of mammary-gland morphogenesis and luminal-cell differentiation,” Nature Cell Biology, vol. 9, no. 2, pp. 201–209, 2007.

D. M. Abd El-Rehim, S. E. Pinder, C. E. Paish et al., “Expression of luminal and basal cytokeratins in human breast carcinoma,” Journal of Pathology, vol. 203, no. 2, pp. 661–671, 2004.

K. A. Hoadley, V. J. Weigman, C. Fan et al., “EGFR associated expression profiles vary with breast tumor subtype,” BMC Genomics, vol. 8, article no. 258, 2007.

F. Ferrari, S. Bortoluzzi, A. Coppe et al., “Novel definition files for human GeneChips based on GeneAnnot,” BMC Bioinformatics, vol. 8, article no. 446, 2007.

Y. Benjamini and Y. Hochberg, “Controlling the false discovery rate. A practical and powerful approach to multiple testing,” Journal of Royal Statistical Society B, vol. 57, pp. 284–300, 1995.

F. Husson, S. Lê, and J. Pagès, Exploratory Multivariate Analysis by Example Using R, Computer Science and Data Analysis Series, CRC Press, Boca Raton, Fla, USA, 2010.

J. D. Jobson, Applied Multivariate Data Analysis. Volume 2. Categorical and Multivariate Methods, Springer, Berlin, Germany, 1992.

K. E. Sleeman, H. Kendrick, A. Ashworth, C. M. Isacke, and M. J. Smalley, “CD24 staining of mouse mammary gland cells defines luminal epithelial, myoepithelial/basal and non-epithelial cells,” Breast Cancer Research, vol. 8, no. 1, article no. R7, 2005.

D. Ponti, A. Costa, N. Zaffaroni et al., “Isolation and in vitro propagation of tumorigenic breast cancer cells with stem/progenitor cell properties,” Cancer Research, vol. 65, no. 13, pp. 5506–5511, 2005.

G. Dontu, D. El-Ashry, and M. S. Wicha, “Breast cancer, stem/progenitor cells and the estrogen receptor,” Trends in Endocrinology and Metabolism, vol. 15, no. 5, pp. 193–197, 2004.

C. Casarsa, N. Bassani, F. Ambrogi et al., “Epithelial-to-mesenchymal transition, cell polarity and stemness-associated features in malignant pleural mesothelioma,” Cancer Letters, vol. 302, no. 2, pp. 136–143, 2011.

O. W. Petersen and K. Polyak, “Stem cells in the human breast,” Cold Spring Harbor perspectives in biology, vol. 2, no. 5, Article ID a003160, 2010.

B. Tiede and Y. Kang, “From milk to malignancy: the role of mammary stem cells in development, pregnancy and breast cancer,” Cell Research, vol. 21, no. 2, pp. 245–257, 2011.

A. P. Bracken, N. Dietrich, D. Pasini, K. H. Hansen, and K. Helin, “Genome-wide mapping of polycomb target genes unravels their roles in cell fate transitions,” Genes and Development, vol. 20, no. 9, pp. 1123–1136, 2006.

A. M. Pietersen, B. Evers, A. A. Prasad et al., “Bmi1 regulates stem cells and proliferation and differentiation of committed cells in mammary epithelium,” Current Biology, vol. 18, no. 14, pp. 1094–1099, 2008.

L. Tickenbrock, S. Hehn, B. Sargin et al., “Activation of Wnt signalling in acute myeloid leukemia by induction of Frizzled-4,” International Journal of Oncology, vol. 33, no. 6, pp. 1215–1221, 2008.

A. M. Band and M. Laiho, “Crosstalk of TGF-β and estrogen receptor signaling in breast cancer,” Journal of Mammary Gland Biology and Neoplasia, vol. 16, no. 2, pp. 109–115, 2011.

G. Dontu, K. W. Jackson, E. McNicholas, M. J. Kawamura, W. M. Abdallah, and M. S. Wicha, “Role of Notch signaling in cell-fate determination of human mammary stem/progenitor cells,” Breast Cancer Research, vol. 6, no. 6, pp. R605–615, 2004.

T. Suzuki, Y. Miki, N. Ohuchi, and H. Sasano, “Intratumoral estrogen production in breast carcinoma: significance of aromatase,” Breast Cancer, vol. 15, no. 4, pp. 270–277, 2008.

G. Hudelist, P. Wülfing, C. Kersting et al., “Expression of aromatase and estrogen sulfotransferase in preinvasive and invasive breast cancer,” Journal of Cancer Research and Clinical Oncology, vol. 134, no. 1, pp. 67–73, 2008.

E. S. Díaz-Cruz, Y. Sugimoto, G. I. Gallicano, R. W. Brueggemeier, and P. A. Furth, “Comparison of increased aromatase versus ERα in the generation of mammary hyperplasia and cancer,” Cancer Research, vol. 71, no. 16, pp. 5477–5487, 2011.

R. Shibuya, T. Suzuki, Y. Miki et al., “Intratumoral concentration of sex steroids and expression of sex steroid-producing enzymes in ductal carcinoma in situ of human breast,” Endocrine-Related Cancer, vol. 15, no. 1, pp. 113–124, 2008.

P. E. Lønning, “The potency and clinical efficacy of aromatase inhibitors across the breast cancer continuum,” Annals of Oncology, vol. 22, no. 3, pp. 503–514, 2011.