Snail-induced epithelial-to-mesenchymal transition of MCF-7 breast cancer cells: systems analysis of molecular changes and their effect on radiation and drug sensitivity
BMC Cancer volume 16, Article number: 236 (2016)
Epithelial-to-mesenchymal transition (EMT) has been associated with the acquisition of metastatic potential and the resistance of cancer cells to therapeutic treatments. MCF-7 breast cancer cells engineered to constitutively express the zinc-finger transcriptional repressor gene Snail (MCF-7-Snail cells) have been previously shown to display morphological and molecular changes characteristic of EMT. We report here the results of a comprehensive systems level molecular analysis of changes in global patterns of gene expression and levels of glutathione and reactive oxygen species (ROS) in MCF-7-Snail cells and the consequence of these changes on the sensitivity of cells to radiation treatment and therapeutic drugs.
Snail-induced changes in global patterns of gene expression were identified by microarray profiling using the Affymetrix platform (U133 Plus 2.0). The resulting data were processed and analyzed by a variety of system level analytical methods. Levels of ROS and glutathione (GSH) were determined by fluorescent and luminescence assays, and nuclear levels of NF-κB protein were determined by an ELISA based method. The sensitivity of cells to ionizing radiation and anticancer drugs was determined using a resazurin-based cell cytotoxicity assay.
Constitutive ectopic expression of Snail in epithelial-like, luminal A-type MCF-7 cells induced significant changes in the expression of >7600 genes including gene and miRNA regulators of EMT. Mesenchymal-like MCF-7-Snail cells acquired molecular profiles characteristic of triple-negative, claudin-low breast cancer cells, and displayed increased sensitivity to radiation treatment, and increased, decreased or no change in sensitivity to a variety of anticancer drugs. Elevated ROS levels in MCF-7-Snail cells were unexpectedly not positively correlated with NF-κB activity.
Ectopic expression of Snail in MCF-7 cells resulted in morphological and molecular changes previously associated with EMT. The results underscore the complexity and cell-type dependent nature of the EMT process and indicate that EMT is not necessarily predictive of decreased resistance to radiation and drug-based therapies.
Breast cancer is the most common female malignancy worldwide with an estimated 1.67 million new cases in 2012 . Despite significant recent progress in the diagnosis and treatment of this biologically and clinically heterogeneous disease, breast cancer remains the most frequent cause of cancer death among women in less developed regions of the world and the second-leading cause of cancer death among women in developed nations [1, 2]. As is the case with most cancers, breast cancer-related deaths are primarily due to metastasis. Metastatic breast cancer (MBC) is present in ~6 % of patients at the time of initial diagnosis and eventually develops in 20–50 % of all breast cancer patients . Since MBC is currently an incurable condition with median survival time of only 0.5–2.2 years, depending on subtype , it continues to be a challenging problem in both basic and clinical cancer research.
Epithelial-to-mesenchymal transition (EMT) is an essential process in normal embryonic development [4, 5] and has been associated with the acquisition of metastatic potential [6, 7] and the resistance of breast and other types of cancers to ionizing radiation  and anticancer drugs (reviewed in ). One of the genes frequently associated with EMT is the zinc-finger transcriptional repressor Snail (SNAI1) . Snail, together with Slug (SNAI2) and Smuc (SNAI3), comprises the Snail family of transcription factors . Previous studies indicate that both Snail and Slug may contribute to the progression of breast and other types of cancer by the down regulation of E-cadherin (CDH1) and other genes associated with the epithelial phenotype and the up regulation of genes associated with the mesenchymal phenotype (reviewed in [10, 12]).
In this study, we were interested in characterizing, on a molecular systems level, the role of Snail in breast cancer EMT and the consequence of this transition on the sensitivity of breast cancer cells to a variety of therapeutic treatments. Toward this end, we performed system level analyses of differences in global patterns of gene expression and therapeutic response profiles between two cell lines derived from the well-studied epithelial breast cancer cell line MCF-7 (Michigan Cancer Foundation-7) . MCF-7-Snail is a derivative of MCF-7 that has been stably transfected with a variant (Snail-6SA) of Snail and displays a mesenchymal-like morphology. Snail-6SA is a more stable protein than wild-type Snail and it has been shown to display constitutive activity and ability to induce EMT [14, 15]. MCF-7-Control is a derivative of MCF-7 that has been transfected with an empty vector and displays the same epithelial morphology as the parental MCF-7 cell line .
We report here that MCF-7-Snail cells display significant changes in the expression of several master regulators of EMT, including various zinc-finger and basic helix-loop-helix transcription factors, as well as members of the miR-200 family of microRNAs. While MCF-7-Control cells display molecular profiles characteristic of the luminal A (ER-positive, PR-positive, HER2-negative) breast cancer subtype, MCF-7-Snail cells were found to display molecular profiles characteristic of the aggressive triple-negative (ER-negative, PR-negative, HER2-negative), claudin-low breast cancer subtype. In addition, we found that relative to the MCF-7-Control, MCF-7-Snail cells display a higher level of cellular ROS, lower levels of GSH and NF-κB (nuclear factor kappa-light-chain-enhancer of activated B cells) activity, increased sensitivity to ionizing radiation and increased, decreased or no change in sensitivity to several anti-cancer drugs. Our results underscore the complexity of the EMT process in breast cancer cells and its consequence on cancer therapies.
MCF-7-Snail and MCF-7-Control cells, developed as previously described , were kindly provided by Dr. Valerie Odero-Marah (Clark Atlanta University). Transfected MCF-7-Snail and MCF-7-Control cells were selected from several clones to display the highest expression of Snail or the highest phenotypic similarity (doubling time) to the parental MCF-7 cells, respectively. Over-expression of Snail in MCF-7-Snail cells has been demonstrated using the western blot analysis . Cells were routinely maintained in RPMI 1640 medium supplemented with 10 % FBS (Atlanta Biologicals, Lawrenceville, GA), 1 % antibiotic-antimycotic solution (Mediatech-Cellgro, Manassas, VA) and 400 μg/mL G418 (Geneticin, GIBCO) at 37 °C in a humidified atmosphere with 5 % CO2 and sub-cultured when they reach ~80 % confluence. In all experiments, cells were no more than four passages from the originally received MCF-7-Snail and MCF-7-Control cells.
Expression analysis by microarray
MCF-7-Snail and MCF-7-Control cells (three replicates per cell line) were grown in the above-described medium and processed for microarray analysis using the Human Genome U133 Plus 2.0 Array (Affymetrix, Santa Clara, CA, USA). The resulting data were acquired as CEL files and processed with Expression Console software Build 18.104.22.168 (Affymetrix, Santa Clara, CA, USA) using the Affymetrix default analysis setting for PLIER and MAS 5.0 algorithms with annotation file HG-U133 Plus_2, Release 34 from 10/24/2013 (www.affymetrix.com). A detailed description of the microarray experiment and the resulting data are available in the Gene Expression Omnibus repository (GEO, http://www.ncbi.nlm.nih.gov/geo/) under the accession number GSE58252.
Differential expression analysis
Expression signals were converted to PLIER+16 and log2-transformed. Probe sets that displayed absent detection calls (MAS5.0 algorithm) across all chips were removed and log2 PLIER+16 values were used to identify genes differentially expressed between MCF-7-Snail and MCF-7-Control cells using the Significance Analysis of Microarrays (SAM) version 4.01 . Genes were reported as differentially expressed at FDR = 2.12 % and absolute fold change (FC) ≥1.5. Differential gene expression was interpreted in the context of EMT and resistance to anticancer drugs using manually curated lists of 71 genes relevant to EMT and 53 genes relevant to anticancer drug resistance (these genes and their Affymetrix probe set IDs are listed in Additional file 1). The threshold for the expression signal intensities that allows identification of genes as highly likely “not expressed” was calculated by the “funnel-shaped procedure” described by Saviozzi et al.  and used to support lack of expression of selected genes (Additional file 2: Figure S1).
Pathway enrichment analysis
Probesets corresponding to differentially expressed genes were employed for enrichment analysis using the MetaCore suite 6.18 build 65,505 (Thomson Reuters, New York, NY, USA). Briefly, significantly perturbed pathways and process networks were identified by mapping differentially expressed genes onto manually curated GeneGO canonical pathway maps and cell process network models .
For each protein from the list of differentially expressed genes between MCF-7-Snail and MCF-7-Control cells, one step interaction neighbors from the global human interactome were identified using the MetaCore “interactome by protein function” tool (MetaCore suite 6.18 build 65,505; Thomson Reuters) and the local interactome was built by adding them to the protein interaction network built from genes differentially expressed between MCF-7-Snail and MCF-7-Control cells. Observed connectivity of each protein (network object) from this local interactome was compared to its expected connectivity based on the global human interactome and relative connectivity (connectivity ratio) was calculated to identify over-connected or under-connected network objects. Statistical significance of differences between observed and expected connectivities was evaluated using the hypergeometric test and multiplicity was controlled by the FDR procedure . The list of over-connected network objects at FDR = 0.01 was reported.
Transcriptional network building
To elucidate complex relationships among the regulators of EMT in our dataset, a custom transcriptional network was built using the results of the differential expression analysis, previously reported associations between genes and EMT, as well as previously reported information on transcriptional regulation and influence on expression between selected network objects. Differentially expressed (i) transcription factors that were previously reported as major regulators of EMT, (ii) microRNA-200 family members, and (iii) epithelial or mesenchymal phenotype-associated genes coding for adherence junctions, tight junctions and intermediate filaments were employed to build the transcriptional network using the knowledge-based system MapEditor (MetaCore suite 6.18 build 65,505; Thomson Reuters). Relative expression data for network objects were color coded (red: up-regulation; blue: down-regulation in MCF-7-Snail relative to MCF-7-Control cells) and mapped on the transcriptional network. Network objects (genes) were connected in the network if their transcription regulation relationships were previously documented and included in the MetaCore knowledge base.
Gene Set Enrichment Analysis (GSEA)
To identify gene sets significantly enriched in a given phenotype (MCF-7-Snail or MCF-7-Control), GSEA  was performed on the data processed by PLIER+16 without any pre-filtering of probe sets, using categorical phenotype labels, gene set permutation type and signal-to-noise metrics. The following gene sets were employed in the analysis: C2: Curated Gene Sets (4722 gene sets) and C6: Oncogenic Signatures (189 gene sets) from the Molecular Signatures Database (http://www.broadinstitute.org/gsea/msigdb/collections.jsp).
In all enrichment analyses, the statistical significance of enrichment was evaluated using p-values calculated based on hypergeometric distribution and corrected for multiplicity using the false discovery rate (FDR) procedure. Unless stated otherwise, pathways, process networks or gene sets were considered to be significantly enriched, if their q-values were ≤ FDR threshold, for which the expected number of false positive entities was ≤1.
MicroRNA expression analysis by qPCR
Relative expression of miRNA-429, miR-200a, miR-200b and miR-141 in MCF-7-Snail vs MCF-7-Control cells was determined by qPCR using specific TaqMan miRNA assays for miRNA-429, miR-200a, miR-200b and miR-141, and non-coding small nuclear RNA RNU6B as an endogenous reference (Applied Biosystems/Life Technologies, Carlsbad, CA). Total cell RNA was isolated using the mirVana miRNA Isolation Kit (Ambion, Foster City, CA, USA) and cDNAs were prepared using the miRNA-specific stem-loop RT primers and TaqMan MicroRNA Reverse Transcription Kit following the manufacturer’s recommendation. Thereafter, cDNA was amplified using the TaqMan Universal Master Mix II with UNG in the CFX96 Real Time PCR Detection System (BioRad, Hercules, CA) following the manufacturer’s recommendation. Expressions of individual miRNAs in MCF-7-Snail relative to MCF-7-Control cells was calculated from the threshold cycles using the REST 2009 Software (Qiagen, Valencia, CA, USA)  and expressed as means, and the 95 % confidence intervals calculated by bootstrapping technique without normality or symmetrical distribution assumptions. P-values determined by a randomization test represent the probability that the observed difference in expression between MCF-7-Snail and MCF-7-Control cells is due to chance.
Determination of radiation sensitivity
One hundred thousand cells were plated in 2.5 mL of RPMI 1640 medium supplemented with 10 % FBS in 35 mm tissue culture dishes (Corning Incorporated, Corning, NY, USA). After 24 h, the cultures were irradiated in an RS-2000 X-ray irradiator (Rad Source Technologies, Suwanee, GA) at 160 kV and 25 mA on an aluminum specimen shelf four at dose rate ~ 311 cGy/min and single dose levels 2 Gy (39 s), 4 Gy (77 s) and 8 Gy (154 s). Control medium was irradiated at 4 Gy. After the irradiation, cells were allowed to grow for 72 h at 37 °C in a humidified atmosphere with 5 % CO2. For quantification of viable cells, 200 μL of Tox-8 reagent were added to each dish and incubated for 2.5 h at 37 °C in a humidified atmosphere with 5 % CO2. Thereafter, the specimens were transferred to a 96-well plate (200 μL/well) and viable cells were quantified via fluorescence at 560 nm excitation and 590 nm emission. The results were expressed as % of non-irradiated control.
Determination of cell cycle distribution
Cells plated in parallel with cells used in the radiation sensitivity experiment were cultured for 24 h, harvested by trypsinization, fixed and stained for DNA analysis by flow cytometry as previously described . Cell cycle distribution was determined by deconvolution of DNA content histograms, after discrimination of doublets and other cellular aggregates by FlowJo 7.6.5 software (Tree Star, Inc., Ashland, OR, USA) using the Dean-Jet-Fox Model. For each cell line, the flow cytometry DNA analysis was performed on three independent cell cultures and the results are presented as means from these three experiments.
Determination of intracellular level of ROS
MCF-7-Snail and MCF-7-Control cells in the RPMI-1640 medium supplemented with 10 % FBS (20,000 cells/mL) were plated into 96-well black-walled plate (100 μL/well) and incubated for 48 h at 37 °C in a humidified atmosphere with 5 % CO2. Thereafter, the medium was removed and 10 μM solution of 2′,7′-dichlorodihydrofluorescein diacetate (H2DCF-DA, Molecular Probes, Inc., Eugene, OR) in PBS was added to each well (100 μL/well). H2DCF-DA is a general oxidative stress indicator that can detect several types of ROS including hydrogen peroxide, hydroxyl radicals and peroxynitrite . Cells were incubated for additional 30 min at 37 °C in a humidified atmosphere with 5 % CO2 and the fluorescence of the ROS-sensitive dye was measured by a Synergy 4 microplate reader (Biotek, Winooski, VT) with filter set 485/20 nm (excitation), 528/20 nm (emission) and 510 nm full-size mirror. Fluorescence intensity corresponding to the ROS signal was normalized to the quantity of viable cells per well as determined by TOX-8 assay and expressed as mean±SD.
Determination of cellular glutathione
MCF-7-Snail and MCF-7-Control cells in DMEM medium (glutathione-free) supplemented with 10 % FBS (20,000 cells/mL) were plated onto the tissue culture-treated 96-well white-walled plate (100 μL/well) and incubated for 48 h at 37 °C in a humidified atmosphere with 5 % CO2.
Reduced glutathione (GSH) and total cellular glutathione (GSH+GSSG) in MCF-7-Snail and MCF-7-Control cells were quantified using the GSH-Glo Glutathione Assay (Promega, Madison, WI, USA). In this assay, the luciferin derivative Luc-NT is converted in the presence of GSH and glutathione S-transferase (GST) to luciferin that generates a luminescent signal in a coupled reaction catalyzed by firefly luciferase. The assay was performed following the manufacturer’s instructions for adherent cell cultures.
Total cellular glutathione was determined after reduction of GSSG to GSH with tris(2-carboxyethyl) phosphine (TCEP, final concentration 1 mM). The luminescence signal after subtraction of blanks (net RLU) was normalized to the number of viable cells determined by resazurin-based cell viability assay TOX-8 (Sigma–Aldrich, St. Louis, MO). All experiments were performed in triplicate.
Determination of the level of nuclear NF-κB
MCF-7-Snail and MCF-7-Control cells each in three replicated cultures were grown in full growth medium to ~80 % confluence, harvested by scraping and processed to obtain nuclear protein extracts using the CelLytic NuCLEAR Extraction Kit (Sigma-Aldrich) following the manufacturer’s protocol. Protein concentration in nuclear extracts was determined using the Pierce 660 nm Protein Assay (Thermo Scientific, Rockford, IL). NF-κB (p50 subunit) was determined in nuclear protein extracts by an ELISA-based assay using the NF-κB (human p50) Transcription Factor Assay Kit (Cayman Chemical Company, Ann Arbor, MI, USA) in 96-well assay format following the manufacturer’s protocol. After developing plates, the stop solution was added and signals corresponding to the p50 protein levels were read as A450 - A570. Concentration of nuclear NF-κB was expressed as A450 - A570 corrected for non-specific binding signal and normalized to protein concentration in nuclear extracts. The results were expressed as means±SDs.
Determination of drug sensitivity
Sensitivity of MCF-7-Snail and MCF-7-Control cells to the cytotoxic effects of selected conventional anticancer drugs was evaluated using the resazurin-based in vitro toxicology assay kit TOX-8 (Sigma–Aldrich) as previously described . Aliquots of cell suspensions (100 μL/well) were plated onto 96-well black-walled plates at 30,000 cells/mL in RPMI 1640 medium supplemented with 10 % FBS, 1 % antibiotic-antimycotic solution and 400 μg/mL G418. Tested compounds were diluted from the following stock solutions: vincristine (VCR) - 0.4 mM in DMSO; doxorubicin (DOX) – 2 mM in DMSO; methotrexate (MTX) – 1 mM in DMSO; gemcitabine (GEM) – 10 mM in H2O; mitomycin C (MMC) – 10 mM in DMSO; 5-fluorouracil (5-FU) – 16.5 mM in H2O; cisplatin (CPT) – 1.7 mM in 0.9 % NaCl/H2O. Tested compounds dissolved in growth medium at a concentration twice the desired final concentration were added in quadruplicates at 100 μL volumes per well. Incubation of cells with drugs or control medium proceeded for 72 h. After that, 20 μL of the TOX-8 reagent were added to each well and incubated for the next 3 h. The increase of fluorescence was measured at a wavelength of 590 nm using an excitation wavelength of 560 nm. The emission of control wells (no drug treatment) after the subtraction of a blank was taken as 100 % and the results for treatments were expressed as a percentage of the control. The experiment was performed four times and GI50 values (concentrations of tested agents that inhibited growth of cell cultures after 72-h incubation to 50 % of the untreated control) were determined by non-linear regression of log-transformed data using a normalized response-variable slope model (GraphPad Prism 5.01; GraphPad Software, Inc.) and expressed as mean ± SD.
Unless stated otherwise, the statistical significance of differences between means of continuous data was evaluated by Welch-corrected t-test and considered significant for two-tailed p-values <0.05. In the analysis of microarray data, multiplicity of statistical tests was corrected by the FDR approach and the discoveries were considered significant if their rank and q-value would allow not more than one false discovery. When multiple comparisons were performed in other experiments, p-values were adjusted by the Šidák-Holm approach and considered significant for padj <0.05. Unless stated otherwise, results of continuous variables were expressed as means±SD.
MCF-7-Snail cells display morphological and molecular changes characteristic of EMT
MCF-7-Snail cells display an elongated morphology characteristic of mesenchymal-like breast cancer cells (Fig. 1a-b). In contrast MCF-7-Control cells, like their parental MCF-7 cells, display the classic “cobble-stone” morphology characteristic of breast cancer epithelial cells (Fig. 1c-d). Consistent with these morphological differences, MCF-7-Snail and MCF-7-Control cells are remarkably different in their respective patterns of gene expression (Additional file 2: Figure S2).
Among 38,226 probe sets included in the differential expression analysis, over 12,000 probe sets corresponding to 7602 genes were found to display statistically significant differences in expression (4242 up-regulated; 3291 down-regulated; 69 discordant) between MCF-7-Snail and MCF-7-Control cells (FDR = 2.12 % and |FC| ≥1.5 (Additional file 2: Figure S3, Additional files 3 and 4). Among the genes significantly differentially expressed are many that have been previously implicated in EMT, including transcription factors Slug, Zeb1, Zeb2, Twist1  and TCF4  (Fig. 2). Among the differentially expressed genes 69 genes displayed discordant expression with some probe sets detecting their up regulation and other probe sets detecting down regulation (Additional file 4). The over-expression of SNAI1 (Snail) and SNAI2 (Slug) genes in MCF-7-Snail relative to MCF-7-Control cells was confirmed by qPCR (Additional file 2: Figure S4 and Additional file 2: Method S1. We have previously reported the up-regulation of mesenchymal markers CDH2 and VIM and down-regulation of epithelial marker CDH1 in MCF-7-Snail cells by RT-PCR .
The miR-200 family of microRNAs are significantly down-regulated in MCF-7-Snail cells
The regulatory role of microRNAs in EMT is an area of growing interest for both developmental and cancer biologists . Members of the miR-200 family of microRNAs are of particular interest because they have been previously shown to target genes that play central roles in EMT (e.g., Zeb1, Zeb2, Slug) [29, 30]. In addition, more recent studies have demonstrated that miR-429 and other members of the miR-200 family are down regulated in ovarian cancer mesenchymal-like cells and that ectopic over-expression of these miRNAs in these cells is sufficient to induce mesenchymal-to-epithelial transition (MET) [31–33]. In light of these findings and because Snail-induced repression of miR-200 family miRNAs has recently been implicated with EMT in embryonic stem cells (ESC) , we examined levels of miR-200 members in MCF-7-Snail cells relative to controls.
The results of comparative qPCR expression analyses of four miRNA-200 family members in MCF-7-Snail and MCF-7-Control cells indicate that levels of miR-200 family microRNAs are consistently and significantly reduced in the mesenchymal-like MCF-7-Snail cells (Fig. 3). Since a number of EMT associated genes are known or predicted targets of miR-200 family members (Fig. 4), these findings suggest that down regulation of members of the miR-200 family of microRNAs may contribute to the regulatory changes associated with EMT in MCF-7-Snail cells.
Systems analysis provides evidence of a complex regulatory interplay between EMT-associated genes and miR-200 family miRNAs in MCF-7-Snail cells
Transcriptional network analysis of genes and miR-200 miRNAs differentially expressed between MCF-7-Snail and MCF-7-Control cells suggests a complex regulatory relationship among EMT-associated transcription factors, miR-200 family members and various cytoskeletal and junction proteins previously associated with the epithalial/mesenchymal phenotype (Fig. 4).
Pathway Enrichment Analysis of the 7634 network objects that were recognized by MetaCore suite among 7602 differentially expressed genes between MCF-7-Snail and MCF-Control cells, identified 18 significantly enriched pathways for both up- and down-regulated genes (Table 1, FDR = 0.06089). Likewise, six significantly enriched pathways were identified for up-regulated genes only (Additional file 2: Table S1, FDR = 0.1677) and 69 pathways for down-regulated genes only (Additional file 2: Table S2, FDR = 0.01514). Mapping up- and down-regulated genes onto GeneGO Process Network Maps identified 30 significantly enriched networks (Table 2, FDR = 0.03104). As expected, the enriched pathways and networks include those related to EMT (e.g., Fig. 5, Tables 1 and 2). Additionally implicated were other related cellular processes including the Wnt-signaling, Hedgehog signaling, estrogen receptor-mediated signaling, NOTCH-signaling, ERBB-signaling, the endoplasmic reticulum stress pathway, and reactive oxygen species (ROS)-associated processes (Tables 1 and 2).
To further explore genes differentially expressed between MCF-7-Snail and MCF-7-Control cells, we examined their paired (binary) protein interactions with those in the MetaCore human protein interaction network . This interactome analysis identified 164 significantly over-connected, but no significantly under-connected human interactome proteins (Additional file 5). The list of over-connected proteins includes: (i) Snail, Slug and Twist1 consistent with their previously recognized role in EMT ; (ii) GSK3B, RYK, β-catenin (CTNNB1) and TCF7L2 (TCF4) consistent with a role in the Wnt signaling pathway in Snail-induced EMT ; (iii) ESR1, PGR (PR) and androgen receptor suggesting a role of the sex-hormone-receptor-mediated signaling pathway; (iv) cancer-associated transcription factors p53 (TP53) and c-Myc (MYC), and (iv) the pluripotency-associated transcription factors KLF-4, Oct3/4, Nanog, Sox2 [35, 36].
We next focused on transcription factor sub-networks represented among genes significantly differentially expressed between MCF-7-Snail and MCF-7-Control cells. This analysis identified 31 transcription factor-centered networks (Additional file 6) confirming a likely role of several transcription factors previously implicated in EMT (MYC, Oct3/4, ESR1, p53, Nanog, Sox2 and TCF-4), in addition to others (CREB1, SP1, NF-κB, and ETS1). For example, among the 294 genes differentially expressed between MCF-7-Snail and MCF-7-Control cells known to be transcriptionally activated by MYC, 240 (~82 %) were up regulated in MCF-7-Snail cells (Additional file 7). Of the 162 genes differentially expressed between MCF-7-Snail and MCF-7-Control cells known be transcriptionally repressed by MYC, 107 (~66 %) were found to be down regulated in MCF-7-Snail cells (Additional file 7). These results suggest a significant role and up-regulated activity of MYC in Snail-induced EMT in breast cancer.
MCF-7-Snail cells display molecular profiles characteristic of the triple-negative breast cancer subtype
Estrogen receptor 1 (ESR1), progesterone receptor (PGR) and ERBB2 genes were found to be significantly down-regulated in MCF-7-Snail relative to MCF-Control cells (Fig. 6a). In addition, the intensities of ESR1 and PGR transcripts were below computed threshold levels, indicating that these genes are not significantly expressed in MCF-7-Snail cells (Additional file 2: Figure S1B).
Further evidence of the involvement of ESR1 in the molecular changes underlying the phenotypic differences between MCF-7-Snail and MCF-7-Control cells, is the observation that ESR1 is significantly over-connected with genes differentially expressed between the two cell lines (Additional file 5). In addition, of the 262 genes differentially expressed between MCF-7-Snail and MCF-7-Control cells known to be transcriptionally activated by ESR1, 172 (~66 %) were found to be down regulated in MCF-7-Snail cells (Additional file 8). Collectively, these findings consistently indicate that ESR1-mediated signaling is significantly reduced in MCF-7-Snail relative to MCF-Control cells and we have confirmed the down-regulation of ESR1 in MCF-7-Snail cells by qPCR (Additional file 2: Figure S5a and Additional file 2: Method S2).
Although the expression intensity of the ERBB2 gene that encodes HER-2/neu protein appears to exceed the calculated expression threshold (Additional file 2: Figure S1B), we conclude that the ERBB2 gene is not significantly expressed in MCF-7-Snail cells for the following two reasons: (i) the probe set that detected ERBB2 expression in our dataset (216836_s_at) is not specific and may also detect other transcripts, and (ii) differential expression analysis identified the ERBB2 gene as down regulated upon ectopic expression of Snail in MCF-7 cells, which are known to be HER-2/neu negative . This conclusion is further supported by the quantification of total HER-2/neu protein in MCF-7-Snail and MCF-7-Control cells using the ELISA (Additional file 2: Figure S6 and Additional file 2: Method S3).
Taken together our results indicate that while MCF-7-Control cells display molecular profiles characteristic of the luminal A subtype of breast cancer (ER-positive, PR-positive and HER-2/neu-negative) [37, 38], MCF-7-Snail have acquired profiles characteristic of the triple-negative (ER-negative, PR-negative, HER-2/neu-negative), breast cancer subtype.
Gene Set Enrichment Analysis (GSEA) further supports this conclusion in that the enrichment of gene sets previously associated with the ductal-invasive, non-luminal, mesenchymal and/or metaplastic phenotype (e.g., activation of KRAS, LEF1, EGFR and suppression of PTEN functions) are also significantly enriched in MCF-7-Snail cells. Conversely, GSEA of MCF-7-Control cells demonstrated enrichment of gene sets previously associated with the estrogen receptor positive, luminal type, and epithelial differentiated phenotype (Additional file 9).
Triple-negative breast cancers have been sub-classified into specific molecular subtypes that include basal-like breast cancers (CK5/6-positive and/or EGFR-positive)  and claudin-low breast cancers exhibiting low expression of the CDH1 and claudin 3, 4 and 7 genes . Cytokeratins 5 and 6 were not found among genes differentially expressed between MCF-7-Snail and MCF-7-Control cells. Consistent with the previously reported absence of the expression of CK5/6 in MCF-7 cells , our result suggests that these genes are not significantly expressed in either MCF-7-Snail or MCF-7-Control cells. Although the EGFR gene appears to be expressed over the calculated threshold, and up-regulated in MCF-7-Snail cells (Additional file 2: Figure S5b, Additional files 3 and 4), its signal intensity is low (Additional file 2: Figure S1B). Moreover, the ITGB4 gene previously associated with the basal-like subtype  is down regulated in MCF-7-Snail cells (Additional files 3 and 4).
Overall, our gene expression results do not support the basal-like phenotype for MCF-7-Snail cells. In contrast, genes whose down regulation is known to be associated with the claudin-low subtype (CLDN3, CLDN4, CLDN7, CDH1, and ITGB4)  are all significantly down regulated in MCF-7-Snail cells (Fig. 6b, Additional file 2: Figure S1B, Additional file 3 and Additional file 4). Collectively, our findings indicate that ectopic expression of Snail has induced the transformation of MCF-7 cells from the luminal A-like cells to the claudin-low triple-negative breast cancer subtype.
Levels of reactive oxygen species (ROS) are significantly elevated in MCF-7-Snail cells
Reactive oxygen species (ROS) are known to induce a variety of cellular responses contributing to the development and progression of breast and other types of cancer . Since our systems analyses indicated that ROS-associated cellular process networks are enriched by genes differentially expressed between MCF-7-Snail and MCF-7-Control cells (Table 2), we explored the possibility that the ectopic expression of Snail in MCF-7 cells might be associated with changes in levels of ROS contributing to the process of EMT and other Snail-mediated phenotypic changes in MCF-7-Snail cells. The results, presented in Fig. 7, are consistent with this hypothesis and demonstrate significantly higher levels of ROS in mesenchymal-like MCF-7-Snail cells relative to the MCF-7-Controls. These results are also consistent with recently reported findings in mesenchymal-like prostate cancer cells .
Elevated levels of ROS in MCF-7-Snail cells are associated with decreased levels of cellular glutathione and NF-κB activity
Increased intracellular concentrations of ROS in cancer cells have been previously associated with decreased levels of antioxidant enzymes and/or glutathione . We compared intracellular levels of glutathione in MCF-7-Snail and MCF-7-Control cells using a luminescence-based assay. The results presented in Fig. 8a demonstrate lower levels of both reduced (GSH) and total (GSH+GSSG) glutathione levels in MCF-7-Snail cells relative to MCF-7-Controls.
ROS have been shown to induce EMT in mouse mammary epithelial cells via NF-κB-mediated activation of Snail . In addition, NF-κB was reported to play an essential role in the induction and maintenance of EMT in breast cancer . Given the previously observed associations between the ROS and NF-κB activity in breast and other cancer types , we compared the status of NF-κB activity in MCF-7-Snail and MCF-7-Control cells. Our microarray data indicate down regulation of the NF-κB1 gene in MCF-7-Snail cells (Additional files 3 and 4). In addition, we found that the NF-κB transcription network is enriched by genes differentially expressed between MCF-7-Snail and MCF-7-Control cells (Additional file 6). NF-κB transcript levels do not necessarily correlate with NF-κB activity since the protein is activated by proteasome degradation of IκB, and the subsequent translocation of NF-κB dimers from the cytoplasm to the nucleus . Thus, measurement of NF-κB protein levels in cell nuclei is necessary for biologically accurate estimates of NF-κB activity. We determined the level of NF-κB (p50/p65- the most common member of NF-κB/Rel family) by its quantification in nuclear protein extracts using an ELISA-based method. The results, presented in Fig. 8b, demonstrate a significant decrease in nuclear protein levels of NF-κB in MCF-7-Snail cells relative to controls indicating decreased nuclear translocation and activation.
The consequence of ectopic expression of Snail on drug sensitivity is variable
Since EMT has been reported by several groups to enhance the resistance of cancer cells to anticancer drugs (reviewed in ), we evaluated the relative sensitivity of MCF-7-Snail and MCF-7-Control cells to a variety of cancer drugs previously and/or currently employed in breast cancer treatment (vincristine, doxorubicin , mitomycin C, methotrexate , gemcitabine, cisplatin  and 5-fluorouracil ). The sensitivity assay was designed to compare numbers of viable cells in drug-treated cell cultures relative to untreated controls independent of the specific mechanisms of cell cytotoxicity. GI50 values for each drug and cell type (Fig. 9) were determined from dose response curves (Additional file 2: Figure S7).
The results indicate that mesenchymal-like MCF-7-Snail cells are significantly more sensitive to doxorubicin, methotrexate, cisplatin, and 5-fluorouracil relative to epithelial-like MCF-7-Control cells. Conversely, MCF-7-Snail cells display significantly reduced sensitivity to gemcitabine and mitomycin C. While the difference in GI50 values between the two cell types is statistically significant for MMC, the fold change (~1.6 ×) is relatively low and, thus, of questionable biological significance. MCF-7-Snail and MCF-7-Control cells did not display significant differences in sensitivity to vincristine.
Differences in sensitivity of MCF-7-Snail and MCF-7-Control cells to some anticancer drugs tested in this work can be interpreted in the context of changes in the expression of genes known to contribute to resistance against specific anticancer drugs (Fig. 10). For example, increased sensitivity of MCF-7-Snail cells to cisplatin can be attributed to the combined effect of (i) down regulation of cisplatin-sequestering metallothioneins MT1E, MT1F and MT1G , (ii) down regulation of the ABCC3 transporter , down regulation of the nucleotide excision repair enzyme ERCC1 , as well as (iii) down regulation of the anti-apoptotic protein XIAP , and (iv) decreased concentration of GSH (Fig. 8a). In fact, the lower concentration of GSH in MCF-7-Snail cells likely suppressed the effect of up regulation of GSTP1 in MCF-7-Snail cells that would otherwise contribute to the resistance against cisplatin . Likewise, increased sensitivity of MCF-7-Snail cells to doxorubicin could be attributed to the down regulation of HSPB1  and metallothioneins , while increased sensitivity to 5-FU can be attributed to the down regulation of DPYD in MCF-7-Snail cells . On the other hand, lower sensitivity of MCF-7-Snail cells to gemcitabine can be attributed to the up regulation of its target, RRM1 . Overall, these results indicate that EMT is not always associated with increased resistance to anticancer drugs.
MCF 7-Snail cells display increased radiosensitivity relative to MCF-7-Controls
Radiotherapy is an important treatment modality for breast cancer . The efficacy of radiotherapy, however, can be decreased by intrinsic or acquired resistance  possibly associated with EMT [8, 65]. To explore this relationship in our experimental system, we exposed MCF-7-Snail and MCF-7-Control cells to a single dose of X-ray radiation at 2–8 Gy and cells were allowed to replicate for 72 h at which point the number of metabolically active viable cells was determined. The results demonstrate a significantly lower proportion of viable MCF-7-Snail cells at 72-hours post-irradiation (Fig. 11) indicating higher radiation sensitivity of the mesenchymal-like MCF-7-Snail cells relative to epithelial-like MCF-7-Control cells. Consistent with the response to radiation treatment, we also observed prominent multinucleated cells in cultures of irradiated MCF-7-Snail cells (Fig. 11) indicative of radiation-induced deregulation of mitosis and deficient separation of nuclei during cytokinesis [66, 67].
Since radiosensitivity has previously been shown to be dependent on the phase of the cell cycle , we examined the distribution of MCF7-Snail and MCF7-Control cells in specific phases of the cell cycle before irradiation. Cultures of MCF-7-Snail and MCF-7-Control cells plated at densities identical to those used in the radiation sensitivity experiments displayed a significantly higher proportion of MCF-7-Snail cells in the G2/M phase relative to MCF-7-Control cells (Holm-Šidák corrected t-test p = 0.0059; CI95 = 2.483–14.29 %). Differences in the proportion of the two cell types in G0/G1 and S phases of the cell cycle were not significant (Fig. 12 and Additional file 2: Figure S8). Since cells are known to display their greatest sensitivity to irradiation during mitosis and G2 phases of the cell cycle , the higher proportion of MCF-7-Snail cells in the G2/M phase of the cell cycle may contribute to their higher radiosensitivity relative to MCF-7-Control cells.
However, different slopes of semi-log radiosensitivity plots at a medium-high dose range (Additional file 2: Figure S9) suggest that the observed difference in the radiosensitivity between MCF-7-Snail and MCF-7-Control cells cannot be attributed only to the cell cycle differences. In an initial effort to identify molecular differences potentially responsible for the observed differences in radiation sensitivity of these two cell types, we focused on genes previously associated with radiation resistance in breast and other cancers and found that many of them were significantly down regulated in MCF-7-Snail cells (Additional file 4), e.g. IGF1R (insulin-like growth factor one receptor) , REG4 (Regenerating islet-derived protein four) , RALBP1 (RalA-binding protein one)  and ERCC1 (excision repair cross-complementation group one) ) suggesting their potential involvement in the increased sensitivity of MCF-7-Snail cells to radiation treatment.
An abundance of prior evience indicates that breast cancer progression and the development of MBC is intimately associated with EMT [74–76]. In addition, EMT has also been associated with the resistance of breast and other cancers to a variety of anticancer drugs  and ionizing radiation [8, 65]. While dysregulation of the transcription factor Snail has been previously associated with EMT in breast and other cancers [10, 12], its effect on system-wide molecular changes and consequent changes in response to radiation and drug therapies is incompletely understood. In an effort to better understand the role of Snail in breast cancer development, we examined an isogenic pair of breast cancer cell lines MCF-7-Snail and MCF-7-Control. MCF-7-Snail cell line has been previously engineered to ectopically express Snail . Detailed anayses of MCF-7-Snail cells provides a unique opportunity to evaluate biologically and clinically significant consequences of Snail-induced EMT in breast cancer.
Previous studies have demonstrated that MCF-7-Snail cells display reduced expression of the epithelial marker E-cadherin (CDH1), increased expression of the mesenchymal markers VIM and FN1 and enhanced migratory capabilities relative to controls . Consistent with these earlier findings, our microarray analyses of MCF-7-Snail cells demonstrate concerted changes in the expression of a number of genes and pathways previously implicated in EMT . Among the genes significantly up regulated/activated in response to ectopic Snail expression are several key transcription factors known to be critical to EMT development (reviewed in ). For example, the initial stages of EMT are believed to require Snail-induced repression of E-cadherin, while subsequent maintenance may require cooperation of other key transcription factors, including Slug, E47, Zeb2 and Twist . Consistent with these expectations, we observed significant down regulation of E-cadherin and up regulation of Slug, E47, Zeb2, and Twist 1 in MCF-7-Snail relative to MCF-7-Control cells.
Snail and Slug are believed to work cooperatively in EMT [77–79]. Our results suggest that the mechanism for Snail-mediated up regulation of Slug is through the down regulation of estrogen receptor 1 (ESR1). Snail is a known transcriptional repressor of ESR1  that, in turn, is a documented transcriptional repressor of Slug . Consistent with this interactive model, we observed a significant up regulation of Slug and concurrent down regulation of ESR1 in MCF-7-Snail cells relative to controls. Interestingly, this down regulation of ESR1 was found be correlated with other molecular changes indicative of Snail-induced transformation of MCF-7 cells from the luminal A to the clinically more aggressive triple-negative, claudin-low sub-type. This finding is consistent with the potential of Snail-targeted therapy for the treatment of metastatic breast cancer .
Interactions between Snail and Slug may also contribute to our observed down regulation in the expression of the miR-200 family of microRNAs in MCF-7-Snail cells. In agreement with prior studies implicating members of the miR-200 family of microRNAs with EMT/MET (mesenchymal-to-epithelial transition) [28, 32, 33], we observed a significant down regulation of members of the miR-200 family in MCF-7-Snail cells. Our gene expression analyses suggest several possible mechanisms by which ectopic expression of Snail in MCF-7 cells may be contributing to repression of miR-200 family of miRNAs.
Although Snail has been previously implicated as an activator of Zeb1  and Zeb2 , both of which are well-established transcriptional repressors of miR-200 family members , Snail-activation of Zeb1/2 is believed not to be direct [83, 86]. In contrast, Slug is known to be a direct transcriptional activator of Zeb1  and can itself directly repress transcription of miR-200 family members . Thus, Snail-induced activation of Slug in MCF-7-Snail cells may contribute to the down regulation of miR-200 family microRNAs.
Our results also demonstrate a significant down regulation of the transcription factor Krüppel-like factor 5 (KLF5) in MCF-7-Snail cells (Additional files 3 and 4). KLF5 is a known transcriptional activator of miR-200 family members, as well as several other microRNAs previously implicated in EMT (e.g., miR-205 and Let-7a/7b/7c/7d/7e/7g) . It has been previously suggested that KLF5-Smads-p300 complexes activate while ZEB1/2-Smads-p300 complexes repress the transcription of miR-200, and that KLF5 and ZEB1/2 may physically interact or compete to execute opposing functions in miR-200 regulation .
A third possible mechanism of miR-200 family down regulation suggested by our results involves the nuclear receptor PELP1 and the histone deacetylase HDAC2. Both of these genes were found to be significantly up regulated in MCF-7-Snail cells (Additional files 3 and 4). It has recently been shown that PELP1, in cooperation with histone deacetylase HDAC2, can transcriptionally repress mir-200 family microRNAs . Interestingly, we also found that MYC, a known transcriptional activator of PELP1 , is functionally activated in MCF-7-Snail cells. None of the above possible modes of miR-200 family repression are mutually exclusive and suggest a redundant and highly interactive nature of Snail regulatory controls.
Previous studies have implicated Snail-induced activation of TGF-β with EMT in breast cancer cells [91, 92]. While our pathway enrichment analysis identified “TGF-beta-dependent induction of EMT via SMADs” as a significantly enriched pathway (Table 1), this does not necessarily imply TGF-β activation. In fact, our gene expression analysis indicates down-regulation of TGFB1-3, SMAD2 and SMAD3 genes in MCF-7-Snail cells (Additional files 3 and 4) which argues against activation of TGF-β mediated signaling via receptor-regulated Smads . In addition, we failed to observe up regulation of known downstreatm targets of TGF-β (e.g., CTGF, ligand TGFB2 and receptor TGFBR2) in MCF-7-Snail cells, again suggesting the absence of TGF-β-signaling. Consistent with these findings, previous studies of claudin-low breast cancer subtypes indicate that Snail-induced EMT need not be dependent on TGF-β-mediated autocrine signaling .
Several previous studies have noted an interactive relationship between changes in the expression of key regulatory genes and intracellular levels of ROS in EMT [44, 46, 94]. This regulatory relationship coupled with their association with oxidative damage and chronic inflammation has identified ROS as a major contributing factor in the progression of breast and other cancers [43, 94]. Our finding that ROS levels are significantly elevated in MCF-7-Snail cells is consistent with these and other results indicating that elevated levels of ROS contribute to EMT and the maintenance of the mesenchymal phenotype in metastatic cancer cells.
Several studies have previously correlated increasing levels of ROS with activation of the transcription factor NF-κB, possibly as part of the stress response pathway . For example, ROS were shown to activate NF-κB that subsequently induced the expression of Snail in MMP3-mediated EMT of mammary epithelial cells . In contrast to expectations, we observed a significant decrease in the expression of NF-κB1 in MCF-7-Snail cells relative to controls. We also observed a decrease in nuclear protein levels of NF-κB in MCF-7-Snail cells indicating reduced NF-κB activation in these cells. Collectively, our results indicate that Snail-induced increases in intracellular levels of ROS are not coupled with NF-κB activation in MCF-7-Snail cells. The reduced activation of NF-κB in MCF-7-Snail cells may be associated with glutathione depletion, consistent with a previously report that demonstrated IκB kinase-dependent and independent down-regulation of NF-κB activity upon glutathione depletion .
In addition to a systems analysis of the molecular consequences of Snail-induced EMT in MCF-7 cells, we were also interested in documenting the effect of these molecular changes on drug and radiation sensitivity. A number of prior studies have suggested that EMT is associated with increased resistance of cancer cells to anticancer drugs (reviewed in ). However, our findings suggest that this generalization may not be universally correct. While we found that MCF-7-Snail cells are significantly more resistant to gemcitabine and mitomycin C (MMC), they were significantly more sensitive to doxorubicin, methotrexate, cisplatin, and 5-fluorouracil. Similar inconsistencies between EMT and drug sensitivities have been previously reported [96–98]. For example, in a survey of the sensitivity of 54 adherent human cancer cell lines from the NCI-60 cell panel, statistically significant positive correlations between expression levels of the epithelial biomarker E-cadherin and drug sensitivity was found for only ten out of 118 drugs . In another study, non-small cell lung cancer (NSCLC) cells designated as mesenchymal-like based on their molecular profile, were found to be more resistant than epithelial-like NSCLC cells to erlotinib and the PI3K/AKT/mTOR inhibitors GDC0941 and 8-aminoadenosine, but more sensitive to pemetrexed, paclitaxel and docetaxel . Similarly, two mesenchymal-like pancreatic cancer cell lines were found to display greater sensitivity to paclitaxel than epithelial-like pancreatic cancer cell lines . Collectively, these findings indicate that epithelial-mesenchymal cellular phenotypes, per se, are not reliable predictors of relative drug sensitivity. While such inconsistencies may simply be attributable to other differences between cell lines that are not consistently associated with epithelial or mesenchymal phenotype, it may also be a reflection of a continuum of potential EMT-mediated effects on cellular phenotypes .
Another reflection of the inconsistency of EMT-mediated effects on cellular phenotypes is our finding that mesenchymal-like MCF-7-Snail cells are more sensitive to radiation treatment than their parental epithelial-like cells. In contrast to our results, a number of previous studies have associated EMT and the mesenchymal cell phenotype with increased radioresistance [8, 65]. For example, epithelial ovarian cancer cells, in which the ectopic expression of Snail or Slug resulted in the acquisition of stem cell-like properties, were found to be associated with radioresistance . Likewise, mesenchymal-like breast cancer cells that displayed higher resistance to radiation than their parental epithelial-like cells, were also found to exhibit a stem cell-like phenotype including a CD44+/CD24-/low molecular signature . One possible explanation of the apparent discrepancy between our results and the accumulated evidence that mesenchymal-like cancer cells are generally associated with radiation resistance (reviewed in ) is that it is the stem cell-like properties of these cells and not EMT per se that is responsible for their radiation resistance. While stem cells typically display a mesenchymal-like phenotype, not all EMT-induced mesenchymal-like cells may display features characteristic of stem cells. For example, it has been previously shown that cancer stem cells typically contain lower levels of ROS than their more differentiated progeny and that this difference is critical for the maintenance of stem cell phenotype . In contrast, we found that MCF-7-Snail cells contain higher levels of ROS than MCF-7-Control cells. In addition, MCF-7-Snail cells display down regulation of Sox2 and Nanog (Additional files 3 and 4) that have been previously shown to contribute to the maintenance of pluripotency and the self-renewal properties of human embryonic stem cells . In addition, we found that other stem cell biomarkers are also down regulated or unchanged in expression in MCF-7-Snail cells relative to controls (e.g., UTF1 , FBXO15 , ALDH1A3 , CD44, CD133 and ITGB1 ). Furthermore, previous studies have also reported down-regulation of several genes associated with stem cell phenotype early (ABCG2 , CD44, ALDH1A3) or late (SOX9 ) after induction of Snail expression in immortalized human mammary epithelial MCF10A cells . Collectively, these findings indicate that MCF-7-Snail cells are not stem cell-like cells and that EMT per se does not necessarily lead to the acquisition of resistance to radiation and anticancer drugs.
Constitutive ectopic expression of Snail in the epithelial-like, luminal A-type, breast cancer cell line MCF-7 results in significant changes in the expression >7600 genes including master gene and miRNA regulators of EMT. Ectopic expression of Snail induced MCF-7 cells to undergo EMT and to acquire features characteristic of triple-negative, claudin-low breast cancer cells but not of breast cancer stem-like cells. Snail-induced EMT of MCF-7 cells resulted in increased sensitivity to radiation treatment but increased, decreased or no change in sensitivity to a variety of anticancer drugs indicating that EMT is not necessarily predictive of decreased responsiveness to therapeutic treatments. These results underscore the complexity and cell-context dependent nature of EMT-mediated changes in breast cancer cells.
Availability of supporting data
A detailed description of the microarray experiment and the resulting data are available in the Gene Expression Omnibus repository (GEO, http://www.ncbi.nlm.nih.gov/geo/) under the accession number GSE58252.
cancer stem-like cells
glutathione (reduced form)
glutathione (oxidized disulfide form)
metastatic breast cancer
non-small cell lung cancer
reactive oxygen species
Ferlay J, Soerjomataram I, Ervik M, Dikshit R, Eser S, Mathers C, et al. GLOBOCAN 2012 v1.0, Cancer Incidence and Mortality Worldwide: IARC CancerBase No. 11 [Internet]. 2013. http://globocan.iarc.fr. Acessed 07/01/2014.
Lu J, Steeg PS, Price JE, Krishnamurthy S, Mani SA, Reuben J, et al. Breast cancer metastasis: challenges and opportunities. Cancer Res. 2009;69:4951–3.
Kennecke H, Yerushalmi R, Woods R, Cheang MCU, Voduc D, Speers CH, et al. Metastatic behavior of breast cancer subtypes. J Clin Oncol. 2010;28:3271–7.
Bolender DL, Markwald RR. Epithelial-mesenchymal transformation in chick atrioventricular cushion morphogenesis. Scan Electron Microsc. 1979;3:313–21.
Baum B, Settleman J, Quinlan MP. Transitions between epithelial and mesenchymal states in development and disease. Semin Cell Dev Biol. 2008;19:294–308.
Yang J, Mani SA, Donaher JL, Ramaswamy S, Itzykson RA, Come C, et al. Twist, a master regulator of morphogenesis, plays an essential role in tumor metastasis. Cell. 2004;117:927–39.
Thiery JP, Acloque H, Huang RY, Nieto MA. Epithelial-mesenchymal transitions in development and disease. Cell. 2009;139:871–90.
Theys J, Jutten B, Habets R, Paesmans K, Groot AJ, Lambin P, et al. E-Cadherin loss associated with EMT promotes radioresistance in human tumor cells. Radiother Oncol. 2011;99:392–7.
Foroni C, Broggini M, Generali D, Damia G. Epithelial-mesenchymal transition and breast cancer: role, molecular mechanisms and clinical impact. Cancer Treat Rev. 2012;38:689–97.
Lamouille S, Xu J, Derynck R. Molecular mechanisms of epithelial-mesenchymal transition. Nat Rev Mol Cell Biol. 2014;15:178–96.
Nieto MA. The snail superfamily of zinc-finger transcription factors. Nat Rev Mol Cell Biol. 2002;3:155–66.
de Herreros AG, Peiro S, Nassour M, Savagner P. Snail family regulation and epithelial mesenchymal transitions in breast cancer progression. J Mammary Gland Biol Neoplasia. 2010;15:135–47.
Soule HD, Vazquez J, Long A, Albert S, Brennan M. Human cell line from a pleural effusion derived from a breast carcinoma. J Natl Cancer Inst. 1973;51:1409–16.
Zhou BP, Deng J, Xia W, Xu J, Li YM, Gunduz M, et al. Dual regulation of Snail by GSK-3beta-mediated phosphorylation in control of epithelial-mesenchymal transition. Nat Cell Biol. 2004;6:931–40.
Javaid S, Zhang J, Anderssen E, Black JC, Wittner BS, Tajima K, et al. Dynamic chromatin modification sustains epithelial-mesenchymal transition following inducible expression of Snail-1. Cell Rep. 2013;5:1679–89.
Smith BN, Burton LJ, Henderson V, Randle DD, Morton DJ, Smith BA, et al. Snail promotes epithelial mesenchymal transition in breast cancer cells in part via activation of nuclear ERK2. PLoS One. 2014;9:e104987.
Tusher VG, Tibshirani R, Chu G. Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci U S A. 2001;98:5116–21.
Saviozzi S, Iazzetti G, Caserta E, Guffanti A, Calogero RA. Microarray data analysis and mining. Methods Mol Med. 2004;94:67–90.
Shi W, Bessarabova M, Dosymbekov D, Dezso Z, Nikolskaya T, Dudoladova M, et al. Functional analysis of multiple genomic signatures demonstrates that classification algorithms choose phenotype-related genes. Pharmacogenomics J. 2010;10:310–23.
Bessarabova M, Ishkin A, JeBailey L, Nikolskaya T, Nikolsky Y. Knowledge-based analysis of proteomics data. BMC Bioinformatics. 2012;13 Suppl 16:S13.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102:15545–50.
Pfaffl MW, Horgan GW, Dempfle L. Relative expression software tool (REST (c)) for group-wise comparison and statistical analysis of relative expression results in real-time PCR. Nucleic Acids Res. 2002;30(9):e36.
Mezencev R, Updegrove T, Kutschy P, Repovska M, McDonald JF. Camalexin induces apoptosis in T-leukemia Jurkat cells by increased concentration of reactive oxygen species and activation of caspase-8 and caspase-9. J Nat Med. 2011;65:488–99.
Lebel CP, Ischiropoulos H, Bondy SC. Evaluation of the probe 2′,7′-dichlorofluorescin as an indicator of reactive oxygen species formation and oxidative stress. Chem Res Toxicol. 1992;5:227–31.
Zheng H, Kang Y. Multilayer control of the EMT master regulators. Oncogene. 2014;33:1755–63.
Medici D, Hay ED, Olsen BR. Snail and Slug promote epithelial-mesenchymal transition through beta-catenin-T-cell factor-4-dependent expression of transforming growth factor-beta3. Mol Biol Cell. 2008;19:4875–87.
McGrail DJ, Mezencev R, Kieu QM, McDonald JF, Dawson MR. SNAIL-induced epithelial-to-mesenchymal transition produces concerted biophysical changes from altered cytoskeletal gene expression. FASEB J. 2015;29:1280–9.
Zhang J, Ma L. MicroRNA control of epithelial-mesenchymal transition and metastasis. Cancer Metastasis Rev. 2012;31:653–62.
Korpal M, Lee ES, Hu G, Kang Y. The miR-200 family inhibits epithelial-mesenchymal transition and cancer cell migration by direct targeting of E-cadherin transcriptional repressors ZEB1 and ZEB2. J Biol Chem. 2008;283:14910–4.
Liu YN, Yin JJ, Abou-Kheir W, Hynes PG, Casey OM, Fang L, et al. MiR-1 and miR-200 inhibit EMT via Slug-dependent and tumorigenesis via Slug-independent mechanisms. Oncogene. 2013;32:296–306.
Chen J, Wang L, Matyunina LV, Hill CG, McDonald JF. Overexpression of miR-429 induces mesenchymal-to-epithelial transition (MET) in metastatic ovarian cancer cells. Gynecol Oncol. 2011;121:200–5.
Wang L, Mezencev R, Svajdler M, Benigno BB, McDonald JF. Ectopic over-expression of miR-429 induces mesenchymal-to-epithelial transition (MET) and increased drug sensitivity in metastasizing ovarian cancer cells. Gynecol Oncol. 2014;134:96–103.
Park SM, Gaur AB, Lengyel E, Peter ME. The miR-200 family determines the epithelial phenotype of cancer cells by targeting the E-cadherin repressors ZEB1 and ZEB2. Genes Dev. 2008;22:894–907.
Gill JG, Langer EM, Lindsley RC, Cai M, Murphy TL, Kyba M, et al. Snail and the microRNA-200 family act in opposition to regulate epithelial-to-mesenchymal transition and germ layer fate restriction in differentiating ESCs. Stem Cells. 2011;29:764–76.
Ling GQ, Chen DB, Wang BQ, Zhang LS. Expression of the pluripotency markers Oct3/4, Nanog and Sox2 in human breast cancer cell lines. Oncol Lett. 2012;4:1264–8.
Takahashi K, Yamanaka S. Induction of pluripotent stem cells from mouse embryonic and adult fibroblast cultures by defined factors. Cell. 2006;126:663–76.
Subik K, Lee JF, Baxter L, Strzepek T, Costello D, Crowley P, et al. The expression patterns of ER, PR, HER2, CK5/6, EGFR, Ki-67 and AR by immunohistochemical analysis in breast cancer cell lines. Breast Cancer (Auckl). 2010;4:35–41.
Holliday DL, Speirs V. Choosing the right cell line for breast cancer research. Breast Cancer Res. 2011;13:215.
Nielsen TO, Hsu FD, Jensen K, Cheang M, Karaca G, Hu Z, et al. Immunohistochemical and clinical characterization of the basal-like subtype of invasive breast carcinoma. Clin Cancer Res. 2004;10:5367–74.
Herschkowitz JI, Simin K, Weigman VJ, Mikaelian I, Usary J, Hu Z, et al. Identification of conserved gene expression features between murine mammary carcinoma models and human breast tumors. Genome Biol. 2007;8:R76.
Trask DK, Band V, Zajchowski DA, Yaswen P, Suh T, Sager R. Keratins as markers that distinguish normal and tumor-derived mammary epithelial cells. Proc Natl Acad Sci U S A. 1990;87:2319–23.
Perou CM, Sorlie T, Eisen MB, van de Rijn M, Jeffrey SS, Rees CA, et al. Molecular portraits of human breast tumours. Nature. 2000;406:747–52.
Waris G, Ahsan H. Reactive oxygen species: role in the development of cancer and various chronic conditions. J Carcinog. 2006;5:14.
Barnett P, Arnold RS, Mezencev R, Chung LW, Zayzafoon M, Odero-Marah V. Snail-mediated regulation of reactive oxygen species in ARCaP human prostate cancer cells. Biochem Biophys Res Commun. 2011;404:34–9.
Liu RM, Pravia KAG. Oxidative stress and glutathione in TGF-beta-mediated fibrogenesis. Free Radic Biol Med. 2010;48:1–15.
Cichon MA, Radisky DC. ROS-induced epithelial-mesenchymal transition in mammary epithelial cells is mediated by NF-kB-dependent activation of Snail. Oncotarget. 2014;5:2827–38.
Huber MA, Beug H, Wirth T. Epithelial-mesenchymal transition: NF-kappaB takes center stage. Cell Cycle. 2004;3:1477–80.
Morgan MJ, Liu ZG. Crosstalk of reactive oxygen species and NF-kappaB signaling. Cell Res. 2011;21:103–15.
Yamamoto Y, Gaynor RB. Therapeutic potential of inhibition of the NF-kappaB pathway in the treatment of inflammation and cancer. J Clin Invest. 2001;107:135–42.
Rainey JM, Jones SE, Salmon SE. Combination chemotherapy for advanced breast-cancer utilizing vincristine, adriamycin, and cyclophosphamide (Vac). Cancer. 1979;43:66–71.
Powles TJ, Jones AL, Judson IR, Hardy JR, Ashley SE. A randomized trial comparing combination chemotherapy using mitomycin-C, mitoxantrone and methotrexate (3 m) with vincristine, anthracycline and cyclophosphamide (Vac) in advanced breast-cancer. Br J Cancer. 1991;64:406–10.
Julka PK, Chacko RT, Nag S, Parshad R, Nair A, Koppiker CB, et al. A phase 2 study of sequential neoadjuvant chemotherapy with gemcitabine and doxorubicin followed by gemcitabine and cisplatin in patients with large or locally advanced operable breast cancer: results from long-term follow-up. Breast Cancer. 2013;20:357–62.
Livi L, Saieva C, Borghesi S, Paoletti L, Meattini I, Rampini A, et al. Concurrent cyclophosphamide, methotrexate, and 5-fluorouracil chemotherapy and radiotherapy for early breast carcinoma. Int J Radiat Oncol Biol Phys. 2008;71:705–9.
Satoh M, Cherian MG, Imura N, Shimizu H. Modulation of resistance to anticancer drugs by inhibition of metallothionein synthesis. Cancer Res. 1994;54:5255–7.
Zhao YB, Lu HL, Yan A, Yang YM, Meng QW, Sun LC, et al. ABCC3 as a marker for multidrug resistance in non-small cell lung cancer. Sci Rep. 2013;3:3120.
Li QD, Yu JJ, Mu CJ, Yunmbam MK, Slavsky D, Cross CL, et al. Association between the level of ERCC-1 expression and the repair of cisplatin-induced DNA damage in human ovarian cancer cells. Anticancer Res. 2000;20:645–52.
Ma JJ, Chen BL, Xin XY. XIAP gene downregulation by small interfering RNA inhibits proliferation, induces apoptosis, and reverses the cisplatin resistance of ovarian carcinoma. Eur J Obstet Gynecol Reprod Biol. 2009;146:222–6.
Goto S, Iida T, Cho S, Oka M, Kohno S, Kondo T. Overexpression of glutathione S-transferase pi enhances the adduct formation of cisplatin with glutathione in human cancer cells. Free Radic Res. 1999;31:549–58.
Hansen RK, Parra I, Lemieux P, Oesterreich S, Hilsenbeck SG, Fuqua SAW. Hsp27 overexpression inhibits doxorubicin-induced apoptosis in human breast cancer cells. Breast Cancer Res Treat. 1999;56:187–96.
Yap XL, Tan HY, Huang JX, Lai YY, Yip GWC, Tan PH, et al. Over-expression of metallothionein predicts chemoresistance in breast cancer. J Pathol. 2009;217:563–70.
Li LH, Dong H, Zhao F, Tang J, Chen X, Ding J, et al. The upregulation of dihydropyrimidine dehydrogenase in liver is involved in acquired resistance to 5-fluorouracil. Eur J Cancer. 2013;49:1752–60.
Davidson JD, Ma LD, Flagella M, Geeganage S, Gelbert LM, Slapak CA. An increase in the expression of ribonucleotide reductase large subunit 1 is associated with gemcitabine resistance in non-small cell lung cancer cell lines. Cancer Res. 2004;64:3761–6.
Langlands FE, Horgan K, Dodwell DD, Smith L. Breast cancer subtypes: response to radiotherapy and potential radiosensitisation. Br J Radiol. 2013;86:20120601.
Li JY, Li YY, Jin W, Yang Q, Shao ZM, Tian XS. ABT-737 reverses the acquired radioresistance of breast cancer cells by targeting Bcl-2 and Bcl-xL. J Exp Clin Cancer Res. 2012;31:102.
Kurrey NK, Jalgaonkar SP, Joglekar AV, Ghanate AD, Chaskar PD, Doiphode RY, et al. Snail and slug mediate radioresistance and chemoresistance by antagonizing p53-mediated apoptosis and acquiring a stem-like phenotype in ovarian cancer cells. Stem Cells. 2009;27:2059–68.
Vakifahmetoglu H, Olsson M, Zhivotovsky B. Death through a tragedy: mitotic catastrophe. Cell Death Differ. 2008;15:1153–62.
Kroemer G, Galluzzi L, Vandenabeele P, Abrams J, Alnemri ES, Baehrecke EH, et al. Classification of cell death: recommendations of the Nomenclature Committee on Cell Death 2009. Cell Death Differ. 2009;16:3–11.
Biade S, Stobbe CC, Chapman JD. The intrinsic radiosensitivity of some human tumor cells throughout their cell cycles. Radiat Res. 1997;147:416–21.
Pawlik TM, Keyomarsi K. Role of cell cycle in mediating sensitivity to radiotherapy. Int J Radiat Oncol Biol Phys. 2004;59:928–42.
Jameel JK, Rao VS, Cawkwell L, Drew PJ. Radioresistance in carcinoma of the breast. Breast. 2004;13:452–60.
Kobunai T, Watanabe T, Fukusato T. REG4, NEIL2, and BIRC5 gene expression correlates with gamma-radiation sensitivity in patients with rectal cancer receiving radiotherapy. Anticancer Res. 2011;31:4147–53.
Awasthi S, Singhal SS, Yadav S, Singhal J, Drake K, Nadkar A, et al. RLIP76 is a major determinant of radiation sensitivity. Cancer Res. 2005;65:6022–8.
Kawashima A, Nakayama M, Kakuta Y, Abe T, Hatano K, Mukai M, et al. Excision repair cross-complementing group 1 may predict the efficacy of chemoradiation therapy for muscle-invasive bladder cancer. Clin Cancer Res. 2011;17:2561–9.
Tomaskovic-Crook E, Thompson EW, Thiery JP. Epithelial to mesenchymal transition and breast cancer. Breast Cancer Res. 2009;11:213.
Teschendorff AE, Journee M, Absil PA, Sepulchre R, Caldas C. Elucidating the altered transcriptional programs in breast cancer using independent component analysis. PLoS Comput Biol. 2007;3:e161.
Mironchik Y, Winnard Jr PT, Vesuna F, Kato Y, Wildes F, Pathak AP, et al. Twist overexpression induces in vivo angiogenesis and correlates with chromosomal instability in breast cancer. Cancer Res. 2005;65:10801–9.
Lee K, Nelson CM. New insights into the regulation of epithelial-mesenchymal transition and tissue fibrosis. Int Rev Cell Mol Biol. 2012;294:171–221.
Lee K, Gjorevski N, Boghaert E, Radisky DC, Nelson CM. Snail1, Snail2, and E47 promote mammary epithelial branching morphogenesis. EMBO J. 2011;30:2662–74.
Aybar MJ, Nieto MA, Mayor R. Snail precedes slug in the genetic cascade required for the specification and migration of the Xenopus neural crest. Development. 2003;130:483–94.
Dhasarathy A, Kajita M, Wade PA. The transcription factor snail mediates epithelial to mesenchymal transitions by repression of estrogen receptor-alpha. Mol Endocrinol. 2007;21:2907–18.
Ye Y, Xiao Y, Wang W, Yearsley K, Gao JX, Barsky SH. ERalpha suppresses slug expression directly by transcriptional repression. Biochem J. 2008;416:179–87.
Azmi AS, Bollig-Fischer A, Bao B, Park BJ, Lee SH, Yong-Song G, et al. Systems analysis reveals a transcriptional reversal of the mesenchymal phenotype induced by SNAIL-inhibitor GN-25. BMC Syst Biol. 2013;7:85.
Guaita S, Puig I, Franci C, Garrido M, Dominguez D, Batlle E, et al. Snail induction of epithelial to mesenchymal transition in tumor cells is accompanied by MUC1 repression and ZEB1 expression. J Biol Chem. 2002;277:39209–16.
Taki M, Verschueren K, Yokoyama K, Nagayama M, Kamata N. Involvement of Ets-1 transcription factor in inducing matrix metalloproteinase-2 expression by epithelial-mesenchymal transition in human squamous carcinoma cells. Int J Oncol. 2006;28:487–96.
Bracken CP, Gregory PA, Kolesnikoff N, Bert AG, Wang J, Shannon MF, et al. A double-negative feedback loop between ZEB1-SIP1 and the microRNA-200 family regulates epithelial-mesenchymal transition. Cancer Res. 2008;68:7846–54.
Wels C, Joshi S, Koefinger P, Bergler H, Schaider H. Transcriptional activation of ZEB1 by Slug leads to cooperative regulation of the epithelial-mesenchymal transition-like phenotype in melanoma. J Invest Dermatol. 2011;131:1877–85.
Ding X, Park SI, McCauley LK, Wang CY. Signaling between transforming growth factor beta (TGF-beta) and transcription factor SNAI2 represses expression of microRNA miR-203 to promote epithelial-mesenchymal transition and tumor metastasis. J Biol Chem. 2013;288:10241–53.
Zhang B, Zhang Z, Xia S, Xing C, Ci X, Li X, et al. KLF5 activates microRNA 200 transcription to maintain epithelial characteristics and prevent induced epithelial-mesenchymal transition in epithelial cells. Mol Cell Biol. 2013;33:4919–35.
Roy SS, Gonugunta VK, Bandyopadhyay A, Rao MK, Goodall GJ, Sun LZ, et al. Significance of PELP1/HDAC2/miR-200 regulatory network in EMT and metastasis of breast cancer. Oncogene. 2014;33:3707–16.
Zeller KI, Zhao X, Lee CW, Chiu KP, Yao F, Yustein JT, et al. Global mapping of c-Myc binding sites and target gene networks in human B cells. Proc Natl Acad Sci U S A. 2006;103:17834–9.
Dhasarathy A, Phadke D, Mav D, Shah RR, Wade PA. The transcription factors Snail and Slug activate the transforming growth factor-beta signaling pathway in breast cancer. PLoS One. 2011;6:e26514.
Taube JH, Herschkowitz JI, Komurov K, Zhou AY, Gupta S, Yang J, et al. Core epithelial-to-mesenchymal transition interactome gene-expression signature is associated with claudin-low and metaplastic breast cancer subtypes. Proc Natl Acad Sci U S A. 2010;107:15449–54.
Shirakihara T, Saitoh M, Miyazono K. Differential regulation of epithelial and mesenchymal markers by deltaEF1 proteins in epithelial mesenchymal transition induced by TGF-beta. Mol Biol Cell. 2007;18:3533–44.
Radisky DC, Levy DD, Littlepage LE, Liu H, Nelson CM, Fata JE, et al. Rac1b and reactive oxygen species mediate MMP-3-induced EMT and genomic instability. Nature. 2005;436:123–7.
Lou H, Kaplowitz N. Glutathione depletion down-regulates tumor necrosis factor alpha-induced NF-kappaB activity via IkappaB kinase-dependent and -independent mechanisms. J Biol Chem. 2007;282:29470–81.
Reinhold WC, Reimers MA, Lorenzi P, Ho J, Shankavaram UT, Ziegler MS, et al. Multifactorial regulation of E-Cadherin expression: an integrative study. Mol Cancer Ther. 2010;9:1–16.
Byers LA, Diao LX, Wang J, Saintigny P, Girard L, Peyton M, et al. An epithelial-mesenchymal transition gene signature predicts resistance to EGFR and PI3K Inhibitors and Identifies axl as a therapeutic target for overcoming egfr inhibitor resistance. Clin Cancer Res. 2013;19:279–90.
Arumugam T, Ramachandran V, Fournier KF, Wang H, Marquis L, Abbruzzese JL, et al. Epithelial to mesenchymal transition contributes to drug resistance in pancreatic cancer. Cancer Res. 2009;69:5820–8.
Soon L, Tachtsidis A, Fok S, Williams ED, Newgreen DF, Thompson EW. The Continuum of Epithelial Mesenchymal Transition – Implication of Hybrid States for Migration and Survival in Development and Cancer. In: Lyden D, Welch D, Psaila B, editors. Cancer metastasis: biologic basis and therapeutics. New York: Cambridge University Press; 2011.
Santisteban M, Reiman JM, Asiedu MK, Behrens MD, Nassar A, Kalli KR, et al. Immune-induced epithelial to mesenchymal transition in vivo generates breast cancer stem cells. Cancer Res. 2009;69:2887–95.
Pajonk F, Vlashi E, McBride WH. Radiation resistance of cancer stem cells: the 4 R’s of radiobiology revisited. Stem Cells. 2010;28:639–48.
Kobayashi CI, Suda T. Regulation of reactive oxygen species in stem cells and cancer stem cells. J Cell Physiol. 2012;227:421–30.
Fong H, Hohenstein KA, Donovan PJ. Regulation of self-renewal and pluripotency by Sox2 in human embryonic stem cells. Stem Cells. 2008;26:1931–8.
Campbell PA, Perez-Iratxeta C, Andrade-Navarro MA, Rudnicki MA. Oct4 targets regulatory nodes to modulate stem cell function. PLoS One. 2007;2:e553.
Okita K, Ichisaka T, Yamanaka S. Generation of germline-competent induced pluripotent stem cells. Nature. 2007;448:313–7.
Marcato P, Dean CA, Pan D, Araslanova R, Gillis M, Joshi M, et al. Aldehyde dehydrogenase activity of breast cancer stem cells is primarily due to isoform ALDH1A3 and its expression is predictive of metastasis. Stem Cells. 2011;29:32–45.
Liu Y, Nenutil R, Appleyard MV, Murray K, Boylan M, Thompson AM, et al. Lack of correlation of stem cell markers in breast cancer stem cells. Br J Cancer. 2014;110:2063–71.
Ding XW, Wu JH, Jiang CP. ABCG2: a potential marker of stem cells and novel target in stem cell and cancer therapy. Life Sci. 2010;86:631–7.
Guo W, Keckesova Z, Donaher JL, Shibue T, Tischler V, Reinhardt F, et al. Slug and Sox9 cooperatively determine the mammary stem cell state. Cell. 2012;148:1015–28.
This work was supported by an endowment from The Deborah Nash Foundation.
The authors declare that they have no competing interest.
RM and JM conceived the study and wrote the manuscript. RM performed the experiments, analyzed and interpreted the data. LM processed cells for RNA and generated the microarray data files. NJ took part in drug sensitivity tests. All authors read and approved the final manuscript.
List of 71 genes relevant to EMT and 53 genes relevant to anticancer drug resistance assembled by manual curation from appended references. (XLSX 54 kb)
Supplemental Figures S1-S9, supplemental tables S1 and S2 and supplemental methods S1-S3. (PDF 528 kb)
List of 12,104 probe sets corresponding to differentially expressed transcripts between MCF-7-Snail and MCF-7-Control cells identified by Significance Analysis of Microarray (FDR = 2.12 % and |FC| ≥1.5). (XLSX 2199 kb)
List of 7602 unique annotated genes differentially expressed between MCF-7-Snail and MCF-7-Control cells. (XLSX 178 kb)
List of 164 network objects from the MetaCore human protein interactome that are statistically significantly over-connected to the genes differentially expressed between MCF-7-Snail and MCF-7-Control cells. (XLS 71 kb)
List of 31 significant networks centered on transcription factors built from genes differentially expressed between MCF-7-Snail and MCF-7-Control cells. (XLS 39 kb)
List of differentially expressed genes between MCF-7-Snail and MCF-7-Control cells known to MetaCore knowledgebase to be directly transcriptionally regulated by c-Myc. (XLS 610 kb)
List of differentially expressed genes between MCF-7-Snail and MCF-7-Control cells known to MetaCore knowledgebase to be directly transcriptionally regulated by ESR-1. (XLS 339 kb)
List of gene sets from C2: Curated Gene Sets and C6: Oncogenic Signatures gene sets (http://www.broadinstitute.org/gsea/msigdb/collections.jsp) significantly enriched in MCF-7-Snail or MCF-7-Control phenotypes. (XLSX 24 kb)
About this article
Cite this article
Mezencev, R., Matyunina, L.V., Jabbari, N. et al. Snail-induced epithelial-to-mesenchymal transition of MCF-7 breast cancer cells: systems analysis of molecular changes and their effect on radiation and drug sensitivity. BMC Cancer 16, 236 (2016). https://doi.org/10.1186/s12885-016-2274-5