- Research article
- Open Access
Direct repression of MYB by ZEB1 suppresses proliferation and epithelial gene expression during epithelial-to-mesenchymal transition of breast cancer cells
Breast Cancer Research volume 15, Article number: R113 (2013)
Epithelial-to-mesenchymal transition (EMT) promotes cell migration and is important in metastasis. Cellular proliferation is often downregulated during EMT, and the reverse transition (MET) in metastases appears to be required for restoration of proliferation in secondary tumors. We studied the interplay between EMT and proliferation control by MYB in breast cancer cells.
MYB, ZEB1, and CDH1 expression levels were manipulated by lentiviral small-hairpin RNA (shRNA)-mediated knockdown/overexpression, and verified with Western blotting, immunocytochemistry, and qRT-PCR. Proliferation was assessed with bromodeoxyuridine pulse labeling and flow cytometry, and sulforhodamine B assays. EMT was induced with epidermal growth factor for 9 days or by exposure to hypoxia (1% oxygen) for up to 5 days, and assessed with qRT-PCR, cell morphology, and colony morphology. Protein expression in human breast cancers was assessed with immunohistochemistry. ZEB1-MYB promoter binding and repression were determined with Chromatin Immunoprecipitation Assay and a luciferase reporter assay, respectively. Student paired t tests, Mann–Whitney, and repeated measures two-way ANOVA tests determined statistical significance (P < 0.05).
Parental PMC42-ET cells displayed higher expression of ZEB1 and lower expression of MYB than did the PMC42-LA epithelial variant. Knockdown of ZEB1 in PMC42-ET and MDA-MB-231 cells caused increased expression of MYB and a transition to a more epithelial phenotype, which in PMC42-ET cells was coupled with increased proliferation. Indeed, we observed an inverse relation between MYB and ZEB1 expression in two in vitro EMT cell models, in matched human breast tumors and lymph node metastases, and in human breast cancer cell lines. Knockdown of MYB in PMC42-LA cells (MYBsh-LA) led to morphologic changes and protein expression consistent with an EMT. ZEB1 expression was raised in MYBsh-LA cells and significantly repressed in MYB-overexpressing MDA-MB-231 cells, which also showed reduced random migration and a shift from mesenchymal to epithelial colony morphology in two dimensional monolayer cultures. Finally, we detected binding of ZEB1 to MYB promoter in PMC42-ET cells, and ZEB1 overexpression repressed MYB promoter activity.
This work identifies ZEB1 as a transcriptional repressor of MYB and suggests a reciprocal MYB-ZEB1 repressive relation, providing a mechanism through which proliferation and the epithelial phenotype may be coordinately modulated in breast cancer cells.
Epithelial-to-mesenchymal transition (EMT), well described in development , enables carcinoma cells to invade local tissues and metastasize to distant sites . EMT causes cell-cell detachment and basement membrane degradation, permitting cell migration aided by actin cytoskeletal rearrangements. EMT triggers myriad intracellular and extracellular signals, which combine to generate motile cells and provide protection against pro-death signals from the host and anticancer therapies, on the journey to secondary sites and while in the systemic circulation (reviewed in ).
ZEB1 (zinc-finger E-box-binding homeobox 1) is a dual zinc-finger, DNA-binding transcription factor, recognizing bipartite E-boxes (CACCTG, CAGGTG) and/or Z-boxes (CAGGTA) [4, 5]. ZEB1 as with ZEB2, Snail1 and 2, Twist1 and 2, TCF3 and 4, FoxC2, Goosecoid, KLF8 and Id1 orchestrate EMT transcriptional and morphologic changes (reviewed in ). In EMT, ZEB1 is a direct transcriptional repressor of E-cadherin  plakophilin3 , Crumbs3, HUGL2, and Pals1 [9, 10]. ZEB1 may also promote metastasis, as shown in a xenograft mouse model  and significantly higher ZEB1 expression is seen in human breast cancer cell lines of the more mesenchymal/invasive basal B subgroup [11–13].
The transcription factor MYB is an oncogene in human leukemias, and in epithelial cancers of the colon and breast (reviewed in [14, 15]). MYB promotes proliferation and inhibits differentiation . We have shown that MYB drives proliferation and suppresses apoptosis and differentiation in estrogen receptor (ER)-positive breast cancer cells in response to estrogen [16, 17], and that it is essential for mammary carcinogenesis in xenograft and transgenic models .
Mutual regulatory relations have been defined for MYB and ZEB1 in the hematopoietic system. MYB and Ets-1 synergize to overcome transcriptional repression of MYB by ZEB1 , and MYB has been shown to regulate ZEB1 expression in the developing inner ear . Conversely, ZEB1 maintains tight regulatory control over MYB during T-cell differentiation . However, the mechanism of this relation has not been defined, and it has not been reported in a solid tumor (cell) context.
A number of transcriptional repressors of CDH1 have been demonstrated to impede cell-cycle progression directly (reviewed in ). Colon cancer cells undergoing an EMT at the invasive front coincide with the region where ZEB1 is expressed  and display a downregulation of proliferation . Conversely, miR-200 family members, which target ZEB mRNA for degradation , have been shown to have a pro-proliferative role [25, 26], thus promoting the growth of breast cancer cell metastases . However, a pro-proliferative role has also been described for ZEB1, because in some contexts, it represses the cell-cycle inhibitors p21 and p73 [28, 29]. The current study sought to determine the ZEB1/MYB/proliferation interplay in the epidermal growth factor (EGF)-responsive PMC42 model of breast cancer EMT.
The PMC42 model system  comprises the parental cell line PMC42-ET (ET) and its more epithelial variant PMC42-LA (LA). Both lines exhibit EMT in response to EGF [30, 31], with marked differences in EMT-marker protein expression and arrangement . Here we have identified an inverse relation between ZEB1 and MYB throughout these cell states, and also in the breast cancer cell lines MDA-MB-231 and MDA-MB-468. We showed that ZEB1 is a key player in promoting the mesenchymal phenotype and regulating the proliferative rate in ET cells through the direct transcriptional repression of MYB. Release of MYB repression promotes an epithelial phenotype, or MET, in which proliferation is restored. MYB and CDH1 gene expression were found to correlate in human breast cancer cell lines and in primary human breast tumors and metastases. Collectively our data link the transcriptional regulation of MYB by ZEB1 to the proliferative state of cells during EMT-MET in breast cancer, and also indicate a contribution of MYB to the epithelial phenotype.
PMC42-ET (ET) cells were derived from a breast cancer pleural effusion by Dr. Robert Whitehead, Ludwig Institute for Cancer Research, Melbourne, Australia, with appropriate institutional ethics clearance (Institutional Review Board of the Peter MacCallum Hospital, Melbourne) and informed consent of the patient. The PMC42-LA (LA) subline was derived further from the parental PMC42- ET cells by Dr. Leigh Ackland, Deakin University, Melbourne, Australia [30, 33–35]. These cells, and all modified (shRNA) derivatives, were maintained at 37°C, 5% CO2 in RPMI1640 containing 10% fetal calf serum (FCS; Thermo Fisher Scientific, Waltham, MA, USA). T47D cells and MCF-7 cells (used in EGF and hypoxic assays) were obtained from the American Type Culture Collection (ATCC, Rockville, MD, USA). MCF-7 cells used for constitutive MYB expression were provided by the Michigan Cancer Foundation to the Garvan Institute, Sydney, Australia. MDA-MB-468 cells originally from the ATCC were transferred from the Lombardi Cancer Center, Washington, DC, USA. All were maintained in Dulbecco Modified Eagle Medium (DMEM; Thermo Fisher Scientific). For induction of EMT via hypoxia, MDA-MB-468 cells were cultured at 1% O2, 5% CO2 in a hypoxia workstation (INVIVO2400; Ruskinn Technology Ltd, Bridgend, UK). For EGF-induced EMT, approximately 30% confluent LA cells were treated with EGF (Collaborative Biomedical Products, Bedford, MA, USA) at 10 ng/ml. After trialing incubation times and concentrations we concluded that treatment of PMC42-LA cells with EGF for 9 days generated a robust and reproducible EMT (as indicated by QRT-PCR and morphology; data not shown). Three dimensional (3D) culture conditions were created by using the Matrigel “on-top” assay .
Creation of modified cell lines
Generation of shRNA-expressing ET, LA, MDA-MB-468, or MDA-MB-231 cells was achieved by using the Lenti-X HTX packaging system (Clonetech Laboratories, Mountain View, CA, USA) and following the manufacturer’s instructions; shRNA sequences (Table 1) were encoded within a GFP-expressing lentiviral vector (pGIPZ; Thermo Scientific Open Biosystems, distributed by Millenium Science, Surrey Hills, VIC, Australia). Successfully transduced cells were selected by cell sorting for GFP; thus a transduced pool rather than individual clones was used in these studies. Stable MYBs expressing MDA-MB-231 and MCF-7 cells were generated by lentiviral transduction, whereby MYB and green fluorescent protein (GFP) were co-expressed by using the pLV411G vector [37, 38]. Successfully transduced cells were selected by cell sorting for GFP.
Immunofluorescence microscopy, multiplex tandem polymerase chain reaction (MT-PCR), and quantitative real-time PCR (QRT-PCR)
These were performed as previously described . Details of antibodies and QRT-PCR primers used are provided in Tables 2 and 3, respectively. Gene-symbol abbreviations used in the current study are shown in Table 4.
Immunobead-negative selection and FACS analysis for EpCAM
LA, SCRsh-ET (SCR, scrambled control shRNA) and ZEB1sh-ET were trypsinized and counted, and approximately 106 cells were added to 20 μl of beads (Collection Epithelial Enrich, Dynal (Invitrogen) to a final volume of 5 ml and incubated at 4°C with rotation for 2 hours. The bead-cell mixture was then placed in a DynaMag-15 magnet (Invitrogen) to separate beads from solution. The number of unbound cells was then determined and represented as a proportion of the original population. Anti-human EpCAM and CDH1 antibodies (Biolegend, San Diego, CA, USA) were used for fluorescence-activated cell sorting (FACS). Cells were trypsinized and counted, and approximately 106 cells in PBS were incubated with 20 μl of antibody for 30 minutes in the dark at 4°C. Cells were then washed twice with 1 ml of phosphate-buffered saline (PBS) and resuspended in 300 μl of PBS, transferred to a 5-ml polystyrene tube and fluorescence measured by using the LSR Fortessa II (BD Biosciences). Data were analyzed by using FlowJo Software (version 7.6.4).
BrdU (bromodeoxyuridine) pulse labeling and flow cytometry, sulforhodamine B (SRB) colorimetric assay, Western blotting
Clonogenic assay and colony morphology analysis
Cells were trypsinized to a single-cell suspension, and 2,000 cells plated in 10 ml of media in an 8-cm-diameter culture dish. After 9 to 10 days, the cells were fixed in 10% neutral buffered formalin, stained with crystal violet (0.5% wt/vol), and washed with PBS. Plates were dried overnight face down. Images were taken of each plate by using the VersaDoc system, and a ‘‘thresholded” image was generated by the imaging program Image J when particular upper and lower pixel-intensity values were set (in the case of the images in Additional file 1: Figure S1C, and Additional file 2: Figure S2B, it was set to automatically detect the threshold, usually at pixel intensities between 129 and 255). This effectively segmented the image into colonies and background. Colonies were then counted by using this imaging program. Colony morphology was assessed by counting a minimum of 100 colonies/plate and classifying them based on the phenotypes displayed in Figure 1B(ii).
Monolayer wound-healing assay
MYB-transfected MDA-MB-231 cells and sh-MYB-knocked-down LA cells were plated in a 24-well plate set up in triplicate and incubated at 37°C for 24 hours to allow the formation of a confluent monolayer. The cells were then wounded by using a P200 pipette tip. The wounded monolayers were washed with complete media to remove the detached cells. Images of the wound were taken at 0, 24, and 48 hours. Wound areas at each time point were analyzed and quantitated by using ImageJ software.
Cells were trypsinized to form a single-cell suspension and resuspended into invasion medium (serum-free RPMI supplemented with 0.1% BSA). Cell counts and viability were estimated by using the Countess Automated Cell Counter (Invitrogen). Then 1 × 105 cells in 100 μl of invasion medium were added into the upper compartment of each migration chamber (6.5-mm Transwell with 8.0-μm pore polycarbonate membrane insert, Corning). The lower chamber contained 600 μl of chemoattractant medium (RPMI supplemented with 10% FBS and 0.1% BSA). The plates were then incubated at 37°C for 4 hours. The membranes were then washed twice in 1× PBS, fixed in cold methanol for 2 minutes, stained with Diff-Quick (1 minute with eosin, followed by 1 minute with buffered thiazole), and washed with water to remove any excess stains. The nonmigrated cells that remained on the upper surface of the membrane were removed by carefully wiping with a cotton bud. The membranes were left to air dry and subsequently mounted on slides. Five random high-power fields per membrane were imaged and counted by using ImageJ software.
Live cell-migration analyses
Live-cell imaging was performed on cells plated at 1 × 105 cells/dish in 3-cm TC glass-bottom dishes (FluoroDish; World Precision Instruments, Inc., Sarasota, FL, USA) by using magnification × 10 phase-contrast objective on an Olympus IX70 microscope attached to a Spot Monochrome camera (2.1.1). Images were acquired for 120 to 150 frames at one frame per 2 minutes, by using Image-Pro Plus 4.5 (MediaCybernetics, Silver Spring, MD, USA). Analysis was via Image-Pro–Analyser 6.1 (MediaCybernetics).
Antibodies are detailed in Table 2. Cells were cultured on Superfrost plus slides (Thermo Scientific, Scoresby, VIC, Australia) overnight then fixed in ice-cold methanol for 20 minutes at −20°C. After 3× PBS 5-minute washes, slides were immunostained by following the protocol described in the Immunohistochemistry section later.
Human breast primary and matched sentinel lymph-node sections
Formalin-fixed, paraffin-embedded surgical specimens were obtained from St. Vincent’s Pathology, Melbourne, with approval of the St. Vincent’s Hospital Human Research Ethics Committee (006–09). Specific patient consent was not required for this study, on the basis that (i) the tissues were collected by the hospital for hospital procedures, (ii) individual patients could not be identified, and (iii) this project did not affect tissue donors’ disease or treatment.
Antibodies are detailed in Table 2. Paraffin-embedded blocks were sectioned, mounted on Superfrost plus slides, de-waxed, and then brought to water. Slides were incubated in primary antibody overnight at 4°C (all antibodies were diluted in 0.1% bovine serum albumin (BSA)/PBS), washed in 0.1% Tween20/1 × PBS, and then incubated with a biotinylated secondary antibody. Slides were then washed in 1 × PBS and incubated with streptavidin-linked HRP tertiary antibody. Staining was visualized with peroxidase-sensitive Sigmafast 3,3′-diaminobenzidine tablets (DAB; Sigma), resuspended to 0.5 mg/ml in 50 mM Tris, pH 7.6, containing 150 μl H2O2. Slides were counterstained with Mayer’s hematoxylin solution (Amber Scientific, Midvale, WA) and mounted in DPX (BDH, Poole, England).
Chromatin immunoprecipitation (ChIP) assay
SCRsh-ET and ZEB1sh-ET cells were grown in T175 flasks to 70% to 80% confluence (approximately 2 to × 107 cells). Cells were washed 3 times with PBS and incubated 20 minutes at room temperature with 1% formaldehyde. Cross-linking was stopped by the addition of 0.1 M glycine for 5 minutes at room temperature. Cells were washed 3 times with PBS and resuspended in 1.0 ml of lysis buffer (50 mM Tris–HCl, pH 8.1, 1.0% SDS, 10.0 mM EDTA, and protease inhibitor cocktail). Chromatin was sonicated (Bioruptor UCD-200, Diagenode) for 10 cycles of 20-second sonication by using 40% efficiency with 2-minute incubation at 4°C between the sonication steps to achieve an average chromatin length of 1 kbp. Chromatin extracts were centrifuged for 14,000 rpm for 30 minutes at 4°C, and the supernatant was diluted 1:10 with IP dilution buffer (20 mM Tris–HCl, pH 8.1, 1.0% Triton X-100, 2.0 mM EDTA, 150 mM NaCl). Diluted chromatin extracts were pre-cleared with ProteinA/G-agarose (Santa Cruz) for 1 hour at 4°C and subsequently used for immunoprecipitation with 2 μg of ZEB1 antibody (E-20; Santa Cruz) or 2 μg of control anti-goat IgG antibody overnight at 4°C. Immunoprecipitants were collected by the addition of 100 μl ProteinA/G-agarose (Santa Cruz) precoated with BSA and salmon sperm DNA. After incubation for 1 hour at 4°C, the agarose was pelleted and washed for 15 minutes with each of the following buffers: low-salt buffer (20 mM Tris–HCl, pH 8.1, 0.1% SDS, 1.0% Triton X-100, 2.0 mM EDTA, and 150 mM NaCl), high-salt buffer (20 mM Tris–HCl, pH 8.1, 0.1% SDS, 1.0% Triton X-100, 2.0 mM EDTA, 500 mM NaCl), LiCl buffer (10 mM Tris HCl, pH 8.1, 1 mM EDTA, 1.0% deoxycholate, 1.0% NP-40, 0.25 M LiCl), followed by three washes in 1× TE. Complexes were eluted with 500 μl elution buffer (1.0% SDS, 0.1 M NaHCO3). Eluates were adjusted with NaCl (0.2 M) and cross-links reversed by heating the samples at 65°C overnight. Eluates were adjusted with EDTA (10 mM) and 50 mM Tris HCl, pH 6.5) and subsequently treated with Proteinase K (2 μl of 10 mg/ml) for 1 hour at 45°C. DNA was recovered by phenol/chloroform/isoamyl alcohol extraction and ethanol precipitation and resuspended in 50 μl 1× TE. Samples were quantified with SYBR-green real-time PCR analysis (Applied Biosystems). Details of primers used for Q-PCR analysis are found in Table 3. Results were expressed as the percentage of total DNA immunoprecipitated with ZEB1 antibody or control anti-goat IgG relative to unprecipitated input samples that were used as reference for all QRT-PCR reactions.
Plasmid DNA, transfection, and reporter gene assay
pMYB prom + intron1 CAT (chloramphenicol acetyltransferase) and pMYB prom CAT were described in a previous study , and pMYB prom + intron1 CAT b (used in this study) was in the previous study named 5′-3′. In pMYB prom + intron1 CAT construct, the promoter, exon 1, and intron 1 sequences of MYB gene were cloned upstream of the CAT. Both sets of E/Z boxes identified in Figure 2A were included in this construct. pMYB prom + intron1 CAT construct b was a shortened version, containing only the first set of Z/E boxes upstream of CAT. pMYB prom contained the promoter region only, encompassing only the first Z-box. This was then mutated by using the QuickChange Site-Directed Mutagenesis Kit (Stratagene/Integrated Sciences), such that CAGGTA was mutated to GTGGTA, and named MYB prom MUT CAT. Expression of exogenous ZEB1 from the expression construct (pcDNA3.1-ZEB1)  was confirmed with Western blot by using a goat anti-ZEB1 antibody C-20, sc-10570; Santa Cruz). A β-galactosidase expression vector was used as an internal control in the reporter assay. Approximately 7 μg of total plasmid DNA was transfected into HEK293T cells by using Fugene HD. In brief, 1.6 × 106 cells/T25 tissue-culture flask was seeded, and transfections were carried out according to manufacturer protocol. After 72 hours of transfection, cells were lyzed by using 1 ml of lysis buffer/T25 flask, and 200 μl of extract was used to quantitate CAT reporter gene activity by using a CAT-ELISA enzyme immunoassay kit (Roche Applied Science, Castle Hill, NSW, Australia).
Student paired t tests were performed by using Microsoft Excel version 2003 (18.8335.8333) SP3. Mann–Whitney and analysis of variance (ANOVA) tests were performed by using Graph Pad Prism (GraphPad Software, San Diego, CA, USA). Where “n” is used (as in n = 3), this refers to the number of independent, biologic replicate experiments performed.
Gene-expression analysis in the PMC42 EMT model
The parental ET cells display mesenchymal characteristics compared with the more-epithelial LA derivative, providing a unique opportunity to investigate the epithelial –mesenchymal axis in isogenic breast cancer cell lines [6, 31, 32]. In investigating the cellular localization and expression of three significant EMT regulators, SNAI1, SNAI2, and ZEB1  between LA and ET cells, ZEB1 nuclear staining and mRNA abundance had the greatest differential and favored the ET cells (Figure 3A, B). By contrast, MYB expression was significantly higher in LA cells (Figure 3B). ET cells also expressed lower levels of E-cadherin (CDH1) and higher levels of vimentin (VIM) than did LA cells (Figure 3B), as previously shown immunocytochemically . Somewhat surprisingly, lower levels of Twist1 (TWIST1) were seen in ET cells compared with LA (Figure 3B).
ZEB1sh-ET cells are epithelially shifted, express higher MYB, and are more proliferative
We previously showed that transient knockdown of ZEB1 with siRNA caused re-expression of CDH1 in ET cells . Here we established stable ZEB1 knocked-down ET cell lines with four shRNA variants and a”:scrambled” nontargeting hairpin control (SCRsh-ET). The level of ZEB1 knockdown was proportional to CDH1 reexpression (Figure 4A, parts i, ii; Further characterization of ZEB1sh-ET cells is shown in Additional file 3). Variant number 4 (hereafter designated ZEB1sh-ET) displayed the greatest ZEB1 knockdown and commensurate CDH1 protein reexpression. Examination of the expression of EMT-related genes revealed that ZEB1 knockdown in ET cells caused MET-like gene-expression changes, including increased expression of CDH1 and CD24, and reduced expression of MMP1, a known target of ZEB1 for induction  (Figure 4A, part iii). Further proepithelial changes due to ZEB1 knockdown are shown in Additional files 3 and 4. These experiments suggested that suppression of ZEB1 caused gene-expression changes similar to the paired LA-versus-ET model of MET. Similarly, expression of MYB protein (Figure 4B) and mRNA (Additional file 3D) was higher in ZEB1sh-ET cells. Pax-2 was also increased (Figure 4A, part iii), but no change was seen in a variety of other markers, including CDH4, CDH13, ZEB2, SNAI1, SNAI2, TWIST1, MMP9, MMP13, and MMP14 (Additional file 1A).
ZEB1 knockdown in ET cells was associated with an increased rate of proliferation, as shown by BrdU (bromodeoxyuridine) pulse labeling and flow cytometry (Figure 4C, part i). Over the 15-hour period after arrest, we observed a significantly higher proportion of ZEB1sh-ET cells in S phase compared with SCRsh-ET cells, indicative of an increase in proliferation. The ZEB1sh-ET cells also proliferated faster than SCRsh-ET cells in the sulforodamine B (SRB) growth assay , (Figure 4C, part ii, further repeats shown in Additional file 1B). ZEB1sh-ET cells also grew larger colonies from a single cell compared with the SCRsh-ET cells (Additional file 1C), an increase that was negated when MYB was knocked down (Additional file 2A and B). Fewer colonies were seen after ZEB1 knockdown, suggesting that ZEB1 may control other models of regulation including survival, apoptosis, anoikis, all of which have been linked to EMT. In MDA-MB-231 cells, ZEB1 knockdown also led to increased CDH1 at the cell membrane and MYB expression (Figure 4D, part i), the latter being shown to be nuclear (Figure 4D, part ii). Consistent with previous work , we observed that these ZEB1-knockdown cells, when plated in 3D, exhibited diminished sprouting and retarded wound closure (data not shown).
MYB and ZEB1 expression are inversely correlated in various in vivo and in vitrobiologic settings
We examined the ZEB1/MYB relation more broadly in EMT and in human breast tumors. Hypoxia has been reported to induce EMT in MDA-MB-468 cells [48, 49]. In a time course of 5 days of hypoxia, only ZEB1 was significantly induced of the eight CDH1 repressor genes examined (Figure 5A, part i, immunocytochemistry for E-cadherin and vimentin: Additional file 5A). This was observed at day 5, when MYB expression was also significantly reduced. MYB repression also accompanied hypoxia-induced EMT-like expression changes in MCF-7 cells (see Additional file 6A).
After 9-day treatment of PMC42-LA cells with 10 ng/ml EGF, VIM and ZEB1 expression were upregulated, and MYB expression was downregulated (Figure 5A, part ii, morphologic changes: Additional file 5B). EGF-induced EMT also led to the repression of MYB in the cell line T47D (Additional file 6B).
Given that tumor aggressiveness has been linked with EMT (, reviewed in ) we examined the expression of MYB in a publically available microarray dataset series of human breast tumors  and observed that MYB expression was generally high in Luminal A and Normal-like tumors but dramatically lower in Basal and Her2 + tumors (Figure 5B). This finding is consistent with previous work that also showed that MYB expression in luminal tumors correlated with a good prognosis . Analysis of publically available expression data  from 51 human breast cancer cell lines that clustered into Luminal, BasalA, and BasalB subgroups [6, 11] also revealed significantly higher MYB expression in Luminal and Basal A subgroups compared with Basal B (Additional file 5C). This pattern was the inverse to that of ZEB1, but correlated with CDH1.
In our study, serial sections from paraffin-embedded matched tissue sets of primary tumors and involved lymph nodes from 11 individual patients were immunostained with the various markers listed in Figure 5C, part i. Histologic examination of another primary tumor-lymph node set from this group of 11 is shown in Additional file 5D. MYB was found to be nuclear and generally tumoral, and associated with epithelial markers CDH1 and cytokeratin along with ERα, particularly in primary tumor samples (invasive ductal breast carcinomas). MYB and associated markers were inversely expressed with regard to VIM. This is further supported by a positive correlation of MYB and CDH1 in breast cancer bone metastases, with both markers staining inversely to VIM (Figure 5C, part ii).
The manipulation of MYB expression influences epithelial versus mesenchymal state
Given that ZEB1 knockdown led to MYB reexpression and epithelialization, we investigated whether these changes were causally linked. Knockdown of MYB mRNA (Figure 1A, i) and protein (Figure 1A, ii) in the more-epithelial LA cells (MYBsh-LA) reduced CDH1 and induced VIM protein expression (Figure 1A, ii). To investigate this further, we created a stable MYB-overexpressing MDA-MB-231 cell line (Figure 1B, part i), and found a shift from mesenchymal, scattered colonies to more-compact epithelial colonies in these cells compared with vector control cells, when plated very sparsely, such that colonies would be seeded from single cells (Figure 1B, part ii). Despite these morphologic changes, preliminary experiments showed no effect of MYB manipulation on monolayer wound closure in either the MYB-transfected MDA-MB-231 cells or sh-MYB-knocked-down LA cells (n = 1; data not shown), and although a trend toward reduced Transwell migration was seen after MYB transfection of MDA-MB-231 cells, this was not significant (n = 3; data not shown). However, live cell-imaging analysis of monolayer cells plated relatively sparsely revealed an altered morphology and migration in the MYB-transfected MDA-MB-231 cells, which grew mostly as connected, flatter cells in epithelial groups compared with a carpet of smaller, fibroblast-like, bipolar, and round cells in the control cultures, and showed a slower, more-tortuous migration pattern than the vector control cells over a 4-hour period (see Additional file 7). In the 9-day EGF-EMT treatment regimen used previously, MYB knockdown appeared further to enhance the repression of CDH1 (see Additional file 8A, n = 1). In MYB-overexpressing MCF-7 cells, the epithelial phenotype was enhanced, as evidenced by significant reductions in VIM and SNAI1 expression (Additional file 8B).
Potential reciprocal repression of ZEB1 by MYB
QRT-PCR analyses of the MYBsh-LA cells revealed an upregulation of ZEB1 expression in variant 2 (Figure 6(i)), the variant with the best MYB knockdown (Figure 1A, part i), suggesting that MYB represses ZEB1 expression in PMC42-LA. Similarly, ZEB1 was significantly repressed in MYB-transfected MDA-MB-231 cells (Figure 6(ii)), despite no detectable change in CDH1 expression (data not shown).
MYB is a direct ZEB1 target
These experiments indicated that MYB and ZEB1 are inversely expressed in LA versus ET cells and in MDA-MB-231. We therefore reasoned that ZEB1 might directly repress MYB by binding to Z- or E-box consensus sequences within MYB. We identified four sequences within MYB that conformed to either the Z- or E-box consensus [4, 5] (indicated as Z-box 1, E-box 1, E-box 2, and Z-box 2; Figure 2A part i and Figure 2B). Z-box 1 (ACAGGTAA) encompassed nt position +219 to +266, such that this sequence spanned the MYB exon 1 and intron 1 junction (Figure 2A, part I, and Figure 2B). E-box 1 (GCAGGTGG) encompassed nt position +294 to +301 and was positioned within intron 1 (Figure 2A, part I, and Figure 2B). E-box 2 (CCAGGTGC) encompassed nt position +550 to +558 and was positioned within intron 1 (Figure 2A, part I, and Figure 2B). Z-box 2 (GCAGGTAG) encompassed nt position +687 to +694 and was also positioned within intron 1 (Figure 2A, part i, and Figure 2B). Each of the Z- and E-box sites was conserved between human, mouse, and dog DNA sequences (Additional file 9A).
We next examined via a CAT reporter-gene assay whether ZEB1 has a direct functional effect on the MYB via the region containing Z-box 1, E-box 1, E-box 2 and Z-box 2 (Figure 2A, parts i and iii) . In this assay, we co-transfected an MYB promoter-intron 1-CAT construct containing Z-box 1, E-box 1, E-box 2, and Z-box 2 with a ZEB1 expression vector (pcDNA3.1-ZEB1) into HEK 293 cells and assessed CAT reporter activity (Figure 2A, parts i and iii). Protein expression from the ZEB1 construct was confirmed by Western blot (Figure 2A, part ii). CAT reporter-gene assays showed that the expression of exogenous ZEB1 in HEK 293 cells repressed the activity of the MYB promoter-intron 1-CAT reporter (Figure 2A, part iii). These data suggest that ZEB1 can negatively regulate MYB. This finding is consistent with the previous observations that have shown that ZEB1 acts as a repressor of gene activity [8–10].
To determine whether the Z-box 1, E-box 1, E-box 2, and Z-box 2 sequences contributed to the ZEB1-mediated repression of CAT activity, we next assessed a series of CAT-reporter constructs in which the Z- and E-box sequences were removed or mutated (MYB prom + intron 1 CATb, deletes E-box 2 and Z-box 2; MYB prom CAT, deletes E-box 1; MYB prom MUT CAT, mutates Z-box 1). CAT reporter-gene assays showed that ZEB1 repressed the activity of all of the deletion and mutant CAT reporter constructs including MYB prom MUT CAT (Figure 2A, parts i and iii). The MYB gene can positively autoregulate itself via sites in its proximal promoter . Co-transfection of an MYB expression vector with MYB prom MUT CAT confirmed the activity of this reporter.
Collectively these data suggest that, in the context of the CAT reporter construct, the repressive activity of ZEB1 in HEK 293 cells was not dependent on the presence of the Z-box 1, E-box 1, E-box 2, or Z-box 2. We did not observe additional Z- or E-box consensus sites within the MYB promoter, exon1, or intron 1. This suggests that ZEB regulation of MYB may occur via ZEB1 binding to alternative nonconsensus site(s) within intron 1 or via alternative upstream sequences within exon 1 or the MYB promoter.
We next examined whether endogenous ZEB1 directly bound MYB in vivo by chromatin immunoprecipitation (ChIP) analysis. We undertook ChIP of SCRsh-ET and ZEB1sh-ET cells with an antibody directed against ZEB1, and experiments were performed at a standard chromatin fragment resolution of 1 kbp (see Methods). QPCR analysis of the immunoprecipitated chromatin used an optimized primer pair located at the 5′ end of intron 1 (Figure 2B). This approach allowed the analysis of ZEB1 occupancy of MYB across the region encompassing the 5′ end of intron 1 and the downstream proximal promoter and exon 1 (Figure 2B). ChIP data in Figure 2C suggest that in SCRsh-ET cells, MYB was enriched for ZEB1 binding. In contrast in ZEB1sh-ET cells, in which ZEB protein expression was decreased (shown by Western blot in Figure 4A, part i), ZEB1 binding was reduced (Figure 2C). Additional repeats of this experiment are shown in Additional file 9B. These data suggest that endogenous ZEB1 directly bound MYB in vivo in a region encompassing the 5′ end of intron 1 and the flanking proximal promoter and exon 1.
Collectively, the CAT reporter and ChIP analyses suggest that MYB is a target of ZEB1 and that ZEB1 exerts its repressive activity via nonconsensus sites located within the region spanning the 5′ end of intron 1 and the flanking proximal promoter of MYB.
Breast cancer metastases in soft tissue and bone have been reported to reexpress CDH1 [55, 56]. We have proposed that the restoration of proliferation required for a secondary tumor to form may be linked to the reversion of EMT, or MET, and we suggest that MYB plays a central role here in promoting tumor growth and the epithelial phenotype. Furthermore, we have established a functional and reciprocal relation between ZEB1 and MYB that is prevalent in breast cancer systems. Direct repression of MYB by the EMT driver ZEB1 builds onto the paradigm through which tumor growth, invasion, and metastasis are integrated (Figure 7). Resistance to apoptosis , reduced proliferation [5, 58], and tumor dormancy [59, 60] have all been associated with EMT in the past, but the potential role of MYB in this process is novel.
Examining the consequences of MYB knockdown in LA cells (Figure 1A), and the overexpression of MYB expression in MDA-MB-231 cells (Figure 1B) and MCF-7 cells (see Additional file 8B) has shown (a) that the manipulation of MYB expression has a direct consequence on epithelial versus mesenchymal state, and (b) MYB expression is concordant with the epithelial phenotype in terms of marker expression, culture and colony morphology, and random migration in the monolayer. We have also found increased MYB expression on CDH1 knockdown in MDA-MB-468 cells (not shown).
How could MYB expression determine epithelial state? One possibility that we have considered is through promotion of KLF4 expression. KLF4 was noted as a novel gene upregulated by MYB via occupation of the KLF4 promoter in MCF-7 cells treated with estradiol . KLF4 is a metastasis-suppressor gene that drives CDH1 expression, preventing EMT . Thus in an actively growing, ER-positive primary tumor, MYB may transduce estrogen signaling into both proliferation  and promotion of the epithelial phenotype. Further studies will be needed to determine whether this MYB-KLF4 partnership occurs in secondary tumors, whether this may be a means through which MYB could mediate the restoration of CDH1 at these sites, and the functional importance of the epithelial aspect regulated by MYB.
We present evidence to suggest that MYB may repress ZEB1 and hence relieve ZEB1-mediated CDH1 transcriptional repression, restoring the epithelial phenotype. This is suggested by the upregulation of ZEB1 on the knockdown of MYB in LA cells (Figure 6(i)) and the repression of ZEB1 in MYB-overexpressing MDA-MB-231 cells (Figure 6(ii)), and also by previous work in hematopoietic cells . A predicted MYB binding site was identified in intron 1 of ZEB1, the site which is mutated in mice with multiple developmental malformations, some of which were suggested to be related to dysfunctional EMT or MET . Thus, similar to the ZEB1-miR200 reciprocal regulatory relation , our data suggest a mutually antagonistic negative-feedback loop between ZEB1 and MYB. ZEB1 and MYB may therefore act as biologic switches, translating molecular changes in the tumor microenvironment into alterations in cell state. These may be epithelial or mesenchymal, as illustrated in Figure 7.
We have shown that ZEB1 binds to the MYB gene in a region encompassing the promoter and intron 1 (Figure 2B). Furthermore, CAT-reporter experiments suggest that the overexpression of ZEB1 represses MYB promoter-intron 1 activity (Figure 2C). This finding is supported by ET-LA cell comparative studies (Figure 3), the characterization of ZEB1sh-ET cells (Figure 4), and a potential role in the EMT of MDA-MB-468 (hypoxia) and LA cells (EGF) (Figure 5A), leading us to conclude that ZEB1 directly targets and represses MYB. Our deletion-mapping studies suggest that ZEB1 can operate through alternative regions, such as the MYB promoter. Although we did not identify exact consensus binding sites for ZEB1 within the promoter, this suggests that ZEB1 may work through divergent binding sites or indirectly by impinging on the activity of a transactivator, such as MYB. Given that MYB is a gene that positively autoregulates itself via tandem MYB binding sites located in the proximal promoter region , and that ZEB1 has been shown to inhibit the transcriptional activity of MYB , we propose that ZEB1 represses MYB gene expression by intercepting the positive-feedback cycle of MYB. ZEB1 downregulates miR34a in a similar manner by directly repressing the miR34a-positive transactivating protein, ΔNp63, resulting in the acquisition of invasive tumor cell capabilities .
Likewise, our studies outline the consequences of ZEB1-mediated gene repression, that is, the suppression of epithelial gene expression and cellular proliferation in breast cancer cells.
The role of MYB in EMT is controversial, and most likely is context dependent. MYB has been positively implicated in the induction of an EMT of avian embryonic neural crest cells  and recently in the induction of Snail2 expression in EMT of embryonic kidney, colon carcinoma, chronic myeloid leukemia-blast crisis, and neuroblastoma cells . Knopfova et al. recently showed increased migration and Matrigel invasion, but not collagen invasion, in MYB-transfected-MDA-MB-231 cells, which contrasts our live-cell imaging results for random monolayer migration. In further contrast to our results, it has been shown that TGF-β-induced EMT of ER-positive breast cancer cells is dependent on MYB upregulation, in which putative miR-200 family binding sites were identified in the MYB mRNA, and the downregulation of these miRs was pinpointed as the mechanism driving this EMT . We have observed a significantly higher level of expression of miR-200 family in LA cells compared with ET (unpublished observation), but increased MYB expression in the LA cells that overexpress all miR-200 family members. In addition, in MYB-overexpressing MCF-7 or MDA-MB-231 cells, we did not observe an induction of Snail2 (Additional file 8B, C) nor did we find any association between MYB induction and EMT, but found the opposite (Figure 5A; Additional file 6A, B). We also found no association between MYB and Snail2 in Luminal versus Basal breast cancer subtypes (see Additional file 10). These dramatic differences could be due to cellular context, because much of our work was performed in the PMC42 model system, which clusters with Basal A/Basal breast cancer cell lines (T. Blick, E. Tomascovic-Crook, unpublished data). However, our results were supported by similar findings in luminal, ER+ breast cancer cells such as MCF-7 and T47D. It seems more likely that our findings are specific to ZEB1, which we find to be expressed later during the EMT, after SNAI1 and 2 induction in the EGF and EGF/staurosporine-induced EMT in our system . Of the CDH1 repressor genes, ZEB1 has been found by others to be the more-dominant EMT driver, important in stabilizing and maintaining the mesenchymal phenotype in breast cancer systems [5, 9, 68]. It is possible that MYB contributes to early EMT responses by driving Snail 1 and 2 expression while repressing ZEB1, until stronger driving forces activate ZEB1, suppress MYB, and thus suppress SNAI1 and 2. Further work is required to delineate these potentially distinct EMT states and the capricious role that MYB may play.
Comparison of the EMT gene-expression profiles in control versus ZEB1 knockdown PMC42-ET cells (Figure 4A, part iii) may provide insights into specific gene-expression changes that may occur during MET at the secondary tumor site. For example, PAX6 was significantly upregulated in ZEB1sh-ET cells. PAX6 has been described as a tumor suppressor, suppressing invasiveness of glioblastoma cells and MMP2 expression , both of which are hallmarks of EMT . PAX6 has also been implicated in promoting the proliferation of breast cancer cells , and we have reported that PAX6 is a direct MYB target gene in the neurogenic zones of the brain . Thus the upregulation of PAX6 may also have contributed to the enhanced proliferative rate observed for ZEB1sh-ET. CDH5 expression was significantly downregulated in ZEB1sh-ET cells (Figure 4A, iii). CDH5 has been shown to be induced in an EMT of mammary tumor cells, influencing the levels of Smad2 phosphorylation and upregulating the expression of TGF-β target genes . CD24 was found to be significantly increased in the ZEB1sh-ET cells, correlating with increased CD24 receptor expression, as determined by FACS (data not shown), consistent with our previous findings that EMT traits in human breast cancer cell lines correlate with the CD44(hi/+)/CD24(lo/-) stem-cell signature [12, 73]. MMP1 expression was significantly reduced in the ZEB1sh-ET cells. Therefore, the ZEB1sh-ET cells reexpress markers of the epithelial phenotype.
We have shown that ZEB1 transcriptionally represses MYB to achieve two outcomes: to downregulate cell proliferation and to stabilize the mesenchymal phenotype. Both outcomes are likely to be beneficial to the survival of the cancer cell in the systemic circulation and at sites of dissemination. EMT has been linked to the acquisition of cell dormancy and reduced chemotherapeutic efficiency (reviewed in ). In addition, a sustained mesenchymal state may enable escape from immune-related apoptotic signals via expression of factors such as SNAI1, which may act to suppress these signals and the immune system itself ; reviewed in . An ongoing debate in relation to the treatment of metastatic breast cancer is whether to customize therapies that reactivate proliferation in dormant, circulating tumor cells and hence improve chemotherapeutic efficiency, or to maintain dormancy in these rogue cells [77, 78]. Each approach is equally problematic, but gene-targeted modulation along the ZEB1-MYB axis may provide an avenue to tackle this foremost clinical dilemma.
Note added in proof
As we noted earlier, the Raschella group initially reported that MYB expression was upregulated in ER-positive human breast cancer cell lines after TGF-β treatment because of transcriptional activation, release of miR-200 family suppression, and protein stabilization, and was instrumental in promoting EMT effects . Since our manuscript was submitted, the Raschella group published a follow-up article  showing that MYB also activates the expression of the miR-200 family, which blocks ZEB1 expression. These new findings now support our key conclusion, that MYB favors the epithelial rather than mesenchymal state and opposes EMT. The authors also support our proposal of contextual issues surrounding MYB regulation of EMT processes in response to signals such as TGF-β, because in the combined presence of MYB and ZEB1, the ZEB1 suppression of miR200 is dominant. This may explain the discrepancies seen in our PMC42 model system compared with their original study. MYB seems poised for only transient expression in the epithelial setting because of the MYB-miR-200 feedback loop. The effects of TGF-β are also known to be very much context dependent.
HJH was supported by a National Breast Cancer Foundation Postdoctoral Training Fellowship; and EWT and HJH were supported also by the Cancer Council of Victoria; NPADG was supported by The Australian Government’s Endeavour Awards Scholarship Program; BvD and ETLS were supported by the National Breast Cancer Foundation (Australia) and Cancer Australia; EWT, DFN, TJG and CP are active members of the EMPathy Breast Cancer Network, a National Breast Cancer Foundation (Australia)-funded National Collaborative Research Program. This study benefited from support by the Victorian Government’s Operational Infrastructure Support Program to St. Vincent’s Institute and Murdoch Children’s Research Institute.
Analysis of variance
American Type Culture Collection
Bovine serum albumin
Dulbecco modified eagle medium
Epidermal growth factor
Fluorescence-activated cell sorting
Fetal calf serum
- LN MET:
Lymph node metastatic
Multiplex tandem polymerase chain reaction
Quantitative reverse-transcription polymerase chain reaction
Roswell Park Memorial Institute medium
Short-hairpin ribonucleic acid
Nieto MA: The ins and outs of the epithelial to mesenchymal transition in health and disease. Annu Rev Cell Dev Biol. 2011, 27: 347-376. 10.1146/annurev-cellbio-092910-154036.
Said NA, Williams ED: Growth factors in induction of epithelial-mesenchymal transition and metastasis. Cells Tissues Organs. 2011, 193: 85-97. 10.1159/000320360.
Gunasinghe NPAD, Thompson EW, Hugo HJ: Mesenchymal-epithelial transition (MET) as a mechanism for metastatic colonisation in breast cancer. Ca Met Rev. 2012, 31: 469-478. 10.1007/s10555-012-9377-5.
Burk U, Schubert J, Wellner U, Schmalhofer O, Vincan E, Spaderna S, Brabletz T: A reciprocal repression between ZEB1 and members of the miR-200 family promotes EMT and invasion in cancer cells. EMBO Rep. 2008, 9: 582-589. 10.1038/embor.2008.74.
Peinado H, Olmeda D, Cano A: Snail, Zeb and bHLH factors in tumour progression: an alliance against the epithelial phenotype?. Nat Rev Cancer. 2007, 7: 415-428. 10.1038/nrc2131.
Hugo H, Ackland ML, Blick T, Lawrence MG, Clements JA, Williams ED, Thompson EW: Epithelial–mesenchymal and mesenchymal–epithelial transitions in carcinoma progression. J Cell Physiol. 2007, 213: 374-383. 10.1002/jcp.21223.
Eger A, Aigner K, Sonderegger S, Dampier B, Oehler S, Schreiber M, Berx G, Cano A, Beug H, Foisner R: DeltaEF1 is a transcriptional repressor of E-cadherin and regulates epithelial plasticity in breast cancer cells. Oncogene. 2005, 24: 2375-2385. 10.1038/sj.onc.1208429.
Aigner K, Descovich L, Mikula M, Sultan A, Dampier B, Bonne S, van Roy F, Mikulits W, Schreiber M, Brabletz T, Sommergruber W, Schweifer N, Wernitznig A, Beug H, Foisner R, Eger A: The transcription factor ZEB1 (deltaEF1) represses Plakophilin 3 during human cancer progression. FEBS Lett. 2007, 581: 1617-1624. 10.1016/j.febslet.2007.03.026.
Aigner K, Dampier B, Descovich L, Mikula M, Sultan A, Schreiber M, Mikulits W, Brabletz T, Strand D, Obrist P, Sommergruber W, Schweifer N, Wernitznig A, Beug H, Foisner R, Eger A: The transcription factor ZEB1 (deltaEF1) promotes tumour cell dedifferentiation by repressing master regulators of epithelial polarity. Oncogene. 2007, 26: 6979-6988. 10.1038/sj.onc.1210508.
Spaderna S, Schmalhofer O, Wahlbuhl M, Dimmler A, Bauer K, Sultan A, Hlubek F, Jung A, Strand D, Eger A, Kirchner T, Behrens J, Brabletz T: The transcriptional repressor ZEB1 promotes metastasis and loss of cell polarity in cancer. Cancer Res. 2008, 68: 537-544. 10.1158/0008-5472.CAN-07-5682.
Blick T, Widodo E, Hugo H, Waltham M, Lenburg ME, Neve RM, Thompson EW: Epithelial mesenchymal transition traits in human breast cancer cell lines. Clin Exp Metastasis. 2008, 25: 629-642. 10.1007/s10585-008-9170-6.
Blick T, Hugo H, Widodo E, Waltham M, Pinto C, Mani SA, Weinberg RA, Neve RM, Lenburg ME, Thompson EW: Epithelial mesenchymal transition traits in human breast cancer cell lines parallel the CD44(hi/)CD24 (lo/-) stem cell phenotype in human breast cancer. J Mammary Gland Biol Neoplasia. 2010, 15: 235-252. 10.1007/s10911-010-9175-z.
Neve RM, Chin K, Fridlyand J, Yeh J, Baehner FL, Fevr T, Clark L, Bayani N, Coppe JP, Tong F, Speed T, Spellman PT, DeVries S, Lapuk A, Wang NJ, Kuo WL, Stilwell JL, Pinkel D, Albertson DG, Waldman FM, McCormick F, Dickson RB, Johnson MD, Lippman M, Ethier S, Gazdar A, Gray JW: A collection of breast cancer cell lines for the study of functionally distinct cancer subtypes. Cancer Cell. 2006, 10: 515-527. 10.1016/j.ccr.2006.10.008.
Ramsay RG, Gonda TJ: MYB function in normal and cancer cells. Nat Rev Cancer. 2008, 8: 523-534. 10.1038/nrc2439.
Zhou Y, Ness SA: Myb proteins: angels and demons in normal and transformed cells. Front Biosci. 2011, 16: 1109-1131. 10.2741/3738.
Drabsch Y, Hugo H, Zhang R, Dowhan DH, Miao YR, Gewirtz AM, Barry SC, Ramsay RG, Gonda TJ: Mechanism of and requirement for estrogen-regulated MYB expression in estrogen-receptor-positive breast cancer cells. Proc Natl Acad Sci USA. 2007, 104: 13762-13767. 10.1073/pnas.0700104104.
Drabsch Y, Robert RG, Gonda TJ: MYB suppresses differentiation and apoptosis of human breast cancer cells. Breast Cancer Res. 2010, 12: R55-10.1186/bcr2614.
Miao RY, Drabsch Y, Cross RS, Cheasley D, Carpinteri S, Pereira L, Malaterre J, Gonda TJ, Anderson RL, Ramsay RG: MYB is essential for mammary tumorigenesis. Cancer Res. 2011, 71: 7029-7037. 10.1158/0008-5472.CAN-11-1015.
Postigo AA, Sheppard AM, Mucenski ML, Dean DC: c-Myb and Ets proteins synergize to overcome transcriptional repression by ZEB. EMBO J. 1997, 16: 3924-3934. 10.1093/emboj/16.13.3924.
Kurima K, Hertzano R, Gavrilova O, Monahan K, Shpargel KB, Nadaraja G, Kawashima Y, Lee KY, Ito T, Higashi Y, Eisenman DJ, Strome SE, Griffith AJ: A noncoding point mutation of Zeb1 causes multiple developmental malformations and obesity in Twirler mice. PLoS Genet. 2011, 7: e1002307-10.1371/journal.pgen.1002307.
Postigo AA, Dean DC: Independent repressor domains in ZEB regulate muscle and T-cell differentiation. Mol Cell Biol. 1999, 19: 7961-7971.
Browne G, Sayan AE, Tulchinsky E: ZEB proteins link cell motility with cell cycle control and cell survival in cancer. Cell Cycle. 2010, 9: 886-891. 10.4161/cc.9.5.10839.
Spaderna S, Schmalhofer O, Hlubek F, Berx G, Eger A, Merkel S, Jung A, Kirchner T, Brabletz T: A transient, EMT-linked loss of basement membranes indicates metastasis and poor survival in colorectal cancer. Gastroenterology. 2006, 131: 830-840. 10.1053/j.gastro.2006.06.016.
Rubio CA: Further studies on the arrest of cell proliferation in tumor cells at the invading front of colonic adenocarcinoma. J Gastroenterol Hepatol. 2007, 22: 1877-1881. 10.1111/j.1440-1746.2007.04839.x.
Kent OA, Mullendore M, Wentzel EA, Lopez-Romero P, Tan AC, Alvarez H, West K, Ochs MF, Hidalgo M, Arking DE, Maitra A, Mendell JT: A resource for analysis of microRNA expression and function in pancreatic ductal adenocarcinoma cells. Cancer Biol Ther. 2009, 8: 2013-2024. 10.4161/cbt.8.21.9685.
Xia W, Li J, Chen L, Huang B, Li S, Yang G, Ding H, Wang F, Liu N, Zhao Q, Fang T, Song T, Wang T, Shao N: MicroRNA-200b regulates cyclin D1 expression and promotes S-phase entry by targeting RND3 in HeLa cells. Mol Cell Biochem. 2010, 344: 261-266. 10.1007/s11010-010-0550-2.
Dykxhoorn DM, Wu Y, Xie H, Yu F, Lal A, Petrocca F, Martinvalet D, Song E, Lim B, Lieberman J: miR-200 enhances mouse breast cancer cell colonization to form distant metastases. PLoS One. 2009, 4: e7181-10.1371/journal.pone.0007181.
Hu F, Wang C, Du J, Sun W, Yan J, Mi D, Zhang J, Qiao Y, Zhu T, Yang S: DeltaEF1 promotes breast cancer cell proliferation through down-regulating p21 expression. Biochim Biophys Acta. 1802, 2010: 301-312.
Fontemaggi G, Gurtner A, Strano S, Higashi Y, Sacchi A, Piaggio G, Blandino G: The transcriptional repressor ZEB regulates p73 expression at the crossroad between proliferation and differentiation. Mol Cell Biol. 2001, 21: 8461-8470. 10.1128/MCB.21.24.8461-8470.2001.
Ackland ML, Michalczyk A, Whitehead RH: PMC42, a novel model for the differentiated human breast. Exp Cell Res. 2001, 263: 14-22. 10.1006/excr.2000.5106.
Ackland ML, Newgreen DF, Fridman M, Waltham MC, Arvanitis A, Minichiello J, Price JT, Thompson EW: Epidermal growth factor-induced epithelio-mesenchymal transition in human breast carcinoma cells. Lab Invest. 2003, 83: 435-448. 10.1097/01.LAB.0000059927.97515.FD.
Hugo HJ, Kokkinos MI, Blick T, Ackland ML, Thompson EW, Newgreen DF: Defining the E-cadherin repressor interactome in epithelial-mesenchymal transition: the PMC42 model as a case study. Cells Tissues Organs. 2011, 193: 23-40. 10.1159/000320174.
Whitehead RH, Bertoncello I, Webber LM, Pedersen JS: A new human breast carcinoma cell line (PMC42) with stem cell characteristics. I. Morphologic characterization. J Natl Cancer Inst. 1983, 70: 649-661.
Whitehead RH, Monaghan P, Webber LM, Bertoncello I, Vitali AA: A new human breast carcinoma cell line (PMC42) with stem cell characteristics. II. Characterization of cells growing as organoids. J Natl Cancer Inst. 1983, 71: 1193-1203.
Whitehead RH, Quirk SJ, Vitali AA, Funder JW, Sutherland RL, Murphy LC: A new human breast carcinoma cell line (PMC42) with stem cell characteristics. III. Hormone receptor status and responsiveness. J Natl Cancer Inst. 1984, 73: 643-648.
Lee GY, Kenny PA, Lee EH, Bissell MJ: Three-dimensional culture models of normal and malignant breast epithelial cells. Nat Methods. 2007, 4: 359-365. 10.1038/nmeth1015.
Barry SC, Harder B, Brzezinski M, Flint LY, Seppen J, Osborne WR: Lentivirus vectors encoding both central polypurine tract and posttranscriptional regulatory element provide enhanced transduction and transgene expression. Hum Gene Ther. 2001, 12: 1103-1108. 10.1089/104303401750214311.
Skalamera D, Ranall MV, Wilson BM, Leo P, Purdon AS, Hyde C, Nourbakhsh E, Grimmond SM, Barry SC, Gabrielli B, Gonda TJ: A high-throughput platform for lentiviral overexpression screening of the human ORFeome. PLoS One. 2011, 6: e20057-10.1371/journal.pone.0020057.
Hugo HJ, Wafai R, Blick T, Thompson EW, Newgreen DF: Staurosporine augments EGF-mediated EMT in PMC42-LA cells through actin depolymerisation, focal contact size reduction and Snail1 induction: a model for cross-modulation. BMC Cancer. 2009, 9: 235-10.1186/1471-2407-9-235.
Kaproth-Joslin KA, Li X, Reks SE, Kelley GG: Phospholipase C delta 1 regulates cell proliferation and cell-cycle progression from G1- to S-phase by control of cyclin E-CDK2 activity. Biochem J. 2008, 415: 439-448. 10.1042/BJ20080233.
Franken NA, Rodermond HM, Stap J, Haveman J, van Bree C: Clonogenic assay of cells in vitro. Nat Protoc. 2006, 1: 2315-2319. 10.1038/nprot.2006.339.
Hammacher A, Thompson EW, Williams ED: Interleukin-6 is a potent inducer of S100P, which is up-regulated in androgen-refractory and metastatic prostate cancer. Int J Biochem Cell Biol. 2005, 37: 442-450. 10.1016/j.biocel.2004.07.011.
Voigt W: Sulforhodamine B assay and chemosensitivity. Methods Mol Med. 2005, 110: 39-48.
Hugo H, Cures A, Suraweera N, Drabsch Y, Purcell D, Mantamadiotis T, Phillips W, Dobrovic A, Zupi G, Gonda TJ, Iacopetta B, Ramsay RG: Mutations in the MYB intron I regulatory sequence increase transcription in colon cancers. Genes Chromosomes Cancer. 2006, 45: 1143-1154. 10.1002/gcc.20378.
Escriva M, Peiro S, Herranz N, Villagrasa P, Dave N, Montserrat-Sentis B, Murray SA, Franci C, Gridley T, Virtanen I, García de Herreros A: Repression of PTEN phosphatase by Snail1 transcriptional factor during gamma radiation-induced apoptosis. Mol Cell Biol. 2008, 28: 1528-1540. 10.1128/MCB.02061-07.
Hu F, Wang C, Guo S, Sun W, Mi D, Gao Y, Zhang J, Zhu T, Yang S: deltaEF1 promotes osteolytic metastasis of MDA-MB-231 breast cancer cells by regulating MMP-1 expression. Biochim Biophys Acta. 1809, 2011: 200-210.
Vichai V, Kirtikara K: Sulforhodamine B colorimetric assay for cytotoxicity screening. Nat Protoc. 2006, 1: 1112-1116. 10.1038/nprot.2006.179.
Jo M, Lester RD, Montel V, Eastman B, Takimoto S, Gonias SL: Reversibility of epithelial-mesenchymal transition (EMT) induced in breast cancer cells by activation of urokinase receptor-dependent cell signaling. J Biol Chem. 2009, 284: 22825-22833. 10.1074/jbc.M109.023960.
Lester RD, Jo M, Montel V, Takimoto S, Gonias SL: uPAR induces epithelial-mesenchymal transition in hypoxic breast cancer cells. J Cell Biol. 2007, 178: 425-436. 10.1083/jcb.200701092.
Hu Z, Fan C, Oh DS, Marron JS, He X, Qaqish BF, Livasy C, Carey LA, Reynolds E, Dressler L, Nobel A, Parker J, Ewend MG, Sawyer LR, Wu J, Liu Y, Nanda R, Tretiakova M, Ruiz Orrico A, Dreher D, Palazzo JP, Perreard L, Nelson E, Mone M, Hansen H, Mullins M, Quackenbush JF, Ellis MJ, Olopade OI, Bernard PS, Perou CM: The molecular portraits of breast tumors are conserved across microarray platforms. BMC Genomics. 2006, 7: 96-10.1186/1471-2164-7-96.
Aokage K, Ishii G, Ohtaki Y, Yamaguchi Y, Hishida T, Yoshida J, Nishimura M, Nagai K, Ochiai A: Dynamic molecular changes associated with epithelial-mesenchymal transition and subsequent mesenchymal-epithelial transition in the early phase of metastatic tumor formation. Int J Cancer. 2011, 128: 1585-1595. 10.1002/ijc.25500.
Gavert N, Ben-Ze'ev A: Epithelial-mesenchymal transition and the invasive potential of tumors. Trends Mol Med. 2008, 14: 199-209. 10.1016/j.molmed.2008.03.004.
Thorner AR, Parker JS, Hoadley KA, Perou CM: Potential tumor suppressor role for the c-Myb oncogene in luminal breast cancer. PLoS One. 2010, 5: e13073-10.1371/journal.pone.0013073.
Nicolaides NC, Gualdi R, Casadevall C, Manzella L, Calabretta B: Positive autoregulation of c-myb expression via Myb binding sites in the 5′ flanking region of the human c-myb gene. Mol Cell Biol. 1991, 11: 6166-6176.
Kowalski PJ, Rubin MA, Kleer CG: E-cadherin expression in primary carcinomas of the breast and its distant metastases. Breast Cancer Res. 2003, 5: R217-R222. 10.1186/bcr651.
Chao YL, Shepard CR, Wells A: Breast carcinoma cells re-express E-cadherin during mesenchymal to epithelial reverting transition. Mol Cancer. 2010, 9: 179-10.1186/1476-4598-9-179.
Valdes F, Alvarez AM, Locascio A, Vega S, Herrera B, Fernandez M, Benito M, Nieto MA, Fabregat I: The epithelial mesenchymal transition confers resistance to the apoptotic effects of transforming growth factor beta in fetal rat hepatocytes. Mol Cancer Res. 2002, 1: 68-78.
Barrallo-Gimeno A, Nieto MA: The Snail genes as inducers of cell movement and survival: implications in development and cancer. Development. 2005, 132: 3151-3161. 10.1242/dev.01907.
Valastyan S, Weinberg RA: Tumor metastasis: molecular insights and evolving paradigms. Cell. 2011, 147: 275-292. 10.1016/j.cell.2011.09.024.
Yang J, Weinberg RA: Epithelial-mesenchymal transition: at the crossroads of development and tumor metastasis. Dev Cell. 2008, 14: 818-829. 10.1016/j.devcel.2008.05.009.
Quintana AM, Liu F, O'Rourke JP, Ness SA: Identification and regulation of c-Myb target genes in MCF-7 cells. BMC Cancer. 2011, 11: 30-10.1186/1471-2407-11-30.
Yori JL, Johnson E, Zhou G, Jain MK, Keri RA: Kruppel-like factor 4 inhibits epithelial-to-mesenchymal transition through regulation of E-cadherin gene expression. J Biol Chem. 2011, 285: 16854-16863.
Ahn YH, Gibbons DL, Chakravarti D, Creighton CJ, Rizvi ZH, Adams HP, Pertsemlidis A, Gregory PA, Wright JA, Goodall GJ, Flores ER, Kurie JM: ZEB1 drives prometastatic actin cytoskeletal remodeling by downregulating miR-34a expression. J Clin Invest. 2012, 122: 3170-3183. 10.1172/JCI63608.
Karafiat V, Dvorakova M, Krejci E, Kralova J, Pajer P, Snajdr P, Mandikova S, Bartunek P, Grim M, Dvorak M: Transcription factor c-Myb is involved in the regulation of the epithelial-mesenchymal transition in the avian neural crest. Cell Mol Life Sci. 2005, 62: 2516-2525. 10.1007/s00018-005-5297-7.
Tanno B, Sesti F, Cesi V, Bossi G, Ferrari-Amorotti G, Bussolari R, Tirindelli D, Calabretta B, Raschella G: Expression of Slug is regulated by c-Myb and is required for invasion and bone marrow homing of cancer cells of different origin. J Biol Chem. 2011, 285: 29434-29445.
Knopfova L, Beneš P, Pekarčíková L, Hermanová M, Masařík M, Pernicová Z, Souček K, Šmarda J: c-Myb regulates matrix metalloproteinases 1/9, and cathepsin D: implications for matrix-dependent breast cancer cell invasion and metastasis. Mol Cancer. 2012, 11: 15-10.1186/1476-4598-11-15.
Cesi V, Casciati A, Sesti F, Tanno B, Calabretta B, Raschella G: TGF-beta-induced c-Myb affects the expression of EMT-associated genes and promotes invasion of ER + breast cancer cells. Cell Cycle. 2011, 10: 4149-4161. 10.4161/cc.10.23.18346.
Taube JH, Herschkowitz JI, Komurov K, Zhou AY, Gupta S, Yang J, Hartwell K, Onder TT, Gupta PB, Evans KW, Hollier BG, Ram PT, Lander ES, Rosen JM, Weinberg RA, Mani SA: Core epithelial-to-mesenchymal transition interactome gene-expression signature is associated with claudin-low and metaplastic breast cancer subtypes. Proc Natl Acad Sci USA. 2010, 107: 15449-15454. 10.1073/pnas.1004900107.
Mayes DA, Hu Y, Teng Y, Siegel E, Wu X, Panda K, Tan F, Yung WK, Zhou YH: PAX6 suppresses the invasiveness of glioblastoma cells and the expression of the matrix metalloproteinase-2 gene. Cancer Res. 2006, 66: 9809-9817. 10.1158/0008-5472.CAN-05-3877.
Zong X, Yang H, Yu Y, Zou D, Ling Z, He X, Meng X: Possible role of Pax-6 in promoting breast cancer cell proliferation and tumorigenesis. BMB Rep. 2011, 44: 595-600. 10.5483/BMBRep.2011.44.9.595.
Malaterre J, Mantamadiotis T, Dworkin S, Lightowler S, Yang Q, Ransome MI, Turnley AM, Nichols NR, Emambokus NR, Frampton J, Ramsay RG: c-Myb is required for neural progenitor cell proliferation and maintenance of the neural stem cell niche in adult brain. Stem Cells. 2008, 26: 173-181. 10.1634/stemcells.2007-0293.
Labelle M, Schnittler HJ, Aust DE, Friedrich K, Baretton G, Vestweber D, Breier G: Vascular endothelial cadherin promotes breast cancer progression via transforming growth factor beta signaling. Cancer Res. 2008, 68: 1388-1397. 10.1158/0008-5472.CAN-07-2706.
Mani SA, Guo W, Liao MJ, Eaton EN, Ayyanan A, Zhou AY, Brooks M, Reinhard F, Zhang CC, Shipitsin M, Ramsay RG: The epithelial-mesenchymal transition generates cells with properties of stem cells. Cell. 2008, 133: 704-715. 10.1016/j.cell.2008.03.027.
Dave B, Mittal V, Tan NM, Chang JC: Epithelial-mesenchymal transition, cancer stem cells and treatment resistance. Breast Cancer Res. 2012, 14: 202-
Kudo-Saito C, Shirako H, Takeuchi T, Kawakami Y: Cancer metastasis is accelerated through immunosuppression during Snail-induced EMT of cancer cells. Cancer Cell. 2009, 15: 195-206. 10.1016/j.ccr.2009.01.023.
Wu Y, Zhou BP: Snail: more than EMT. Cell Adh Migr. 2010, 4: 199-203. 10.4161/cam.4.2.10943.
Goss PE, Chambers AF: Does tumour dormancy offer a therapeutic target?. Nat Rev Cancer. 2010, 10: 871-877. 10.1038/nrc2933.
Perret GY, Uzzan B: An anticancer strategic dilemma: to kill or to contain: the choice of the pharmaceutical industry in 2009. Fundam Clin Pharmacol. 2011, 25: 283-295. 10.1111/j.1472-8206.2010.00849.x.
Pieraccioli M, Imbastari F, Antonov A, Melino G, Raschella G: Activation of miR200 by c-Myb depends on ZEB1 expression and miR200 promoter methylation. Cell Cycle. 2013, 12: 1-12. 10.4161/cc.23067.
We thank the various sources of financial support listed. We are appreciative of the technical advice given by Mr. Tony Blick and his reviewing of the manuscript. We thank Antonio García de Herreros, Dept Ciencies Experimentals i de la Salut, Universitat Pompeu Fabra, Barcelona, Spain, for providing us with the pcDNA3-ZEB1 expression construct.
The authors declare that they have no competing interests.
HJH conceived the study design and performed all experiments not attributed to others. LP performed chromatin immunoprecipitation studies, Western blot analysis (Figure 4B), and provided support for the CAT reporter analyses; RS and BS carried out cell-cycle analyses and interpretation; YD and TJG created the MYB-modified MCF-7 cell line; CP performed the EpCAM bead-binding assay and EpCAM FACS analyses; CP, ETLS, BvD, and DFN performed cell-migration assays and/or analyses. DG contributed to the experiments with MDA-MB-468 cells; PH assisted in histopathologic interpretations and provided archival sections of breast cancers and matched lymph node metastases; HJH RGR, TJG, DFN, and EWT participated in the study design, coordination, and data interpretation. All authors participated in drafting the manuscript, and read and approved the final version.
Donald F Newgreen and Erik W Thompson contributed equally to this work.
Electronic supplementary material
Additional file 1: (A) The expression (MT-PCR) of other EMT-related genes in SCRsh-ET versus ZEB1sh-ET, not shown in Figure 4 (A, part iii). (B) Additional SRB experiments supporting the SRB data shown in Figure 4C, part ii. (C) Clonogenic assay of SCRsh-ET versus ZEB1sh-ET showing that ZEB1 knockdown led to larger overall colony size, consistent with a higher rate of proliferation. Colony number and size were averaged from thresholded images of six individual 10-cm-diameter plates (n = 6) by using the Image J image-analysis program (*P < 0.05; error bars represent SEM. Statistics were determined by using Student paired t test). (PDF 667 KB)
Additional file 2: MYB-targeted knockdown in ZEB1sh-ET cells abrogates the effect of ZEB1KD on colony size. (A) (i) Validation of MYB knockdown in MYB shRNA expressing SCRsh-ET and ZEB1sh-ET cells by QRT-PCR (fold expression relative to SCRsh-ET/SCRsh-ET shown, n = 1, error bars represent SD) and (ii) Western blotting analysis. (B) clonogenicity assay revealing ZEB1sh/MYBsh-ET cells form smaller colonies, indicating a reduction in proliferative rate (n = 3 independent biologic replicates, error bars represent SEM). (PDF 96 KB)
Additional file 3: (A) DAB immunocytochemistry staining of ZEB1sh-ET shows reduced nuclear ZEB1, magnification 600×, scale bar, 100 μm; (B) The EpCAM receptor was reexpressed in ZEB1sh-ET, as shown by (i) greater retention to EpCAM coated beads (inferred by reduced number of unbound cells in the flow-through population); and (ii) FACS for EpCAM. CDH1 protein reexpression shown in Figure 4A was further confirmed by FACS to be membrane bound. The results from three independent experiments (biologic replicates) are shown. (C) (i) Confocal microscopy of PMC42-ET spheroids grown on Matrigel revealed CDH1 was expressed at the cell membrane in ZEB1sh-ET, magnification 400×, scale bar, 30 μM; (ii) ZEB1sh-ET spheroids were generally larger than SCRsh-ET, consistent with previous work in MDA-MB-231, in which ZEB1 was knocked down . PMC42-ET WT cell organoids were not enlarged and were comparable in size to SCR shRNA-ET control (not shown). 100× magnification; scale bar, 200 μM. (iii) ZEB1sh-ET cells grew as tighter clusters in 2D culture (magnification 200×; scale bar, 50 μM) and exhibited slower epithelial-like movement in this dimension (see Additional file 4). (D) MYB was expressed higher in ZEB1sh-ET cells, regardless of the dimension of culture. (PDF 42 MB)
Additional file 4: PMC42 cells transfected with siRNA: scrambled (A) or Zeb1 (B). Stills from movies at an interval of 3 hours are shown, and paths of movement of seven cells in each field are projected in yellow on the first time frame of each. Traces for PMC42ET-shSCR control (147.3 ± 52.1 μm/4 hours) and PMC42ET-shZeb1 cells (60.9 ±17.2 μm/4 hours; P < 0.0001). Note that the more-epithelial appearance of the Zeb1 siRNA-transfected cells correlates with very short paths. Tracking point was the center of the cell nucleus. (PDF 2 MB)
Additional file 5: (A) Immunofluorescence of MDA-MB-468 cells seeded in Terasaki plates exposed to hypoxic conditions (1% O2) for various times, as indicated. CDH1 sublocalization to the cytoplasm from the membrane was seen, along with an increase in the expression of VIM protein in a time-dependent manner on exposure to hypoxia; magnification 600×; scale bar, 50 μm. (B) Phase-contrast images of untreated or 10 ng/ml EGF-treated PMC42-LA, depicting the acquisition of mesenchymal features: scattering and an elongation of cell shape, magnification 100×; scale bar, 200 μM. (C) Comparison of expression levels of MYB, ZEB1, and CDH1 in Luminal + BasalA versus BasalB subgroups of human breast cancer cells from the Neve dataset . Significance (*) set at P < 0.05, Mann–Whitney statistical test. (D) An additional matched human breast tumor primary (Primary T) and sentinel lymph node metastases (LN MET) set to that shown in Figure 5; magnification, 400×; scale bar, 150 μm. (PDF 55 MB)
Additional file 6: Additional cell models of EMT in which MYB was repressed (A, B). (A) Expression (QRT-PCR) analyses of MCF-7 cells exposed to hypoxia for up to 72 hours. (B) T47D cells treated for up to 120 hours with 10 ng/ml EGF. (PDF 88 KB)
Additional file 7: Representative images of the (A) vector control and (B) MYB-transfected MDA-MB-231 cells under phase-contract microscopy, showing tracking of path length over a 4–hour period. For each, two groups of seven closely spaced cells were selected for tracking. Tracking point was the center of the cell nucleus. Yellow denotes representative 4-hour traces for MDA-MB-231-vector control (147.7 ± 30.6 μm/4 hours; 4-hour displacement = 42% of path length, that is, relatively straight) and MYB-transfected MDA-MB-231 cells (91.2 ±19.01 μm/4 hours; 4-hour displacement = 25% of path length; that is, quite tortuous). (PDF 2 MB)
Additional file 8: (A) Knockdown of MYB in PMC42-LA cells (which were characterized in Figure 1) led to a further reduction of CDH1 expression after induction of EMT by EGF. Data shown are fold change, relative to untreated SCRsh-LA; n = 1, error bars represent standard deviation. (B) Stable overexpression of MYB in MCF-7 cells (i) validated by immunocytochemistry (DAB), magnification 100×, scale bar, 200 μm; and (ii) fold changes in mesenchymal gene expression compared with WT MCF-7. Data shown are an average of three independent experiments; error bars represent SEM, significance (P < 0.05) indicated by *, as determined by Student paired t test. (C) Stable MYB overexpression in MDA-MB-231 cells does not result in Snail2 expression upregulation: MT-PCR data expressed as fold change (Vector control MDA-MB-231 = 1); n = 1; error bars represent standard deviation. (PDF 21 MB)
Additional file 9: (A) E and Z-boxes number1 within the MYB promoter captured in the ChIP analysis are highly conserved between Homo sapiens, Mus muscularis, and Canis familiaris, as are E and Z-boxes number 2 downstream. (B) Additional graphs/replications of ChIP analyses summarized in Figure 2: QRT-PCR of two independent replications of ChIP assay examining the region of the MYB promoter defined in Figure 1. (C) QRT-PCR positive control for ChIP assay: amplification of the region of the CDH1 promoter containing E-boxes at which ZEB1 has been demonstrated to bind. (D) QRT-PCR negative control: amplification of a nonrelated sequence within MYB gene in intron 10. (PDF 681 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Hugo, H.J., Pereira, L., Suryadinata, R. et al. Direct repression of MYB by ZEB1 suppresses proliferation and epithelial gene expression during epithelial-to-mesenchymal transition of breast cancer cells. Breast Cancer Res 15, R113 (2013). https://doi.org/10.1186/bcr3580
- Human Breast Cancer Cell Line
- Epithelial Phenotype
- ZEB1 Expression
- Mesenchymal State
- Epithelial Gene Expression