Profiling gene promoter occupancy of Sox2 in two phenotypically distinct breast cancer cell subsets using chromatin immunoprecipitation and genome-wide promoter microarrays
Breast Cancer Research volume 16, Article number: 470 (2014)
Aberrant expression of the embryonic stem cell marker Sox2 has been reported in breast cancer (BC). We previously identified two phenotypically distinct BC cell subsets separated based on their differential response to a Sox2 transcription activity reporter, namely the reporter-unresponsive (RU) and the more tumorigenic reporter-responsive (RR) cells. We hypothesized that Sox2, as a transcription factor, contributes to their phenotypic differences by mediating differential gene expression in these two cell subsets.
We used chromatin immunoprecipitation and a human genome-wide promoter microarray (ChIP-chip) to determine the promoter occupancies of Sox2 in the MCF7 RU and RR breast cancer cell populations. We validated our findings with conventional chromatin immunoprecipitation, quantitative reverse transcription polymerase chain reaction (qPCR), and western blotting using cell lines, and also performed qPCR using patient RU and RR samples.
We found a largely mutually exclusive profile of gene promoters bound by Sox2 between RU and RR cells derived from MCF7 (1830 and 456 genes, respectively, with only 62 overlapping genes). Sox2 was bound to stem cell- and cancer-associated genes in RR cells. Using quantitative RT-PCR, we confirmed that 15 such genes, including PROM1 (CD133), BMI1, GPR49 (LGR5), and MUC15, were expressed significantly higher in RR cells. Using siRNA knockdown or enforced expression of Sox2, we found that Sox2 directly contributes to the higher expression of these genes in RR cells. Mucin-15, a novel Sox2 downstream target in BC, contributes to the mammosphere formation of BC cells. Parallel findings were observed in the RU and RR cells derived from patient samples.
In conclusion, our data supports the model that the Sox2 induces differential gene expression in the two distinct cell subsets in BC, and contributes to their phenotypic differences.
Sex-determining region Y (SRY)-box binding protein-2 (Sox2) is a transcription factor essential to the maintenance of the pluripotent stem cell state in embryonic stem cells (ESCs) and induced pluripotent stem cells -. In human ESCs, Sox2 governs their pluripotency by binding to the promoters of its target genes and transcriptionally regulating their expressions both positively and negatively . A previous study of Sox2 promoter occupancy in human ESCs using chromatin immunoprecipitation promoter microarray chip analysis (ChIP-chip) has revealed target genes positively regulated by Sox2 (including SOX2, OCT4, NANOG and MYC) . In normal adult tissues, Sox2 is largely restricted to somatic stem cells; specifically, Sox2 expression has been detected in the stem/progenitor cells of the brain, stomach, colon, and anus . In normal mammary glands, Sox2 is largely restricted to the stem cell populations -.
In recent years, Sox2 has been discovered to be aberrantly expressed in cancer cells, including those of the lungs, brain, ovaries, bone, colon, skin, and breasts -. In many of these studies, Sox2 was found in the cancer stem cell population ,,-, supporting the hypothesis that cancer stemness is related to the aberrant expression of ESC proteins. It has been demonstrated that Sox2 promotes key tumorigenic properties in cancer cells, including enhanced proliferation, invasion, migration, colony formation, non-adherent stem cell-associated sphere formations in vitro, and tumorigenicity in vivo ,,-. Further, Sox2 has been shown to correlate with a worse prognosis in cancer patients, including those with breast cancer (BC) ,,-. Up to 30% of BC, including all four major molecular subtypes, have been reported to express Sox2 ,. In a relatively small number of in vitro studies, Sox2 has been directly implicated in promoting cell proliferation, mammosphere formation, invasion and tumorigenesis in BC ,,.
We recently identified and characterized two distinct cell subsets of BC, separated based on their differential responsiveness to a Sox2 transcription activity reporter . Using two estrogen receptor-positive (ER+) cell lines, MCF7 and ZR751, we found that the vast majority of these cells, despite robust levels of Sox2, were reporter unresponsive (labeled as RU cells), while a relatively small cell subset were reporter responsive (labeled as RR cells) . Importantly, RU and RR cells are phenotypically distinct, with RR cells showing a higher expression of the stem cell marker CD49f and exhibiting a higher tumorigenic potential . In view of the fact that Sox2 is a transcription factor, we hypothesized that Sox2 mediates differential gene expressions in RU and RR cells, thereby contributing to their phenotypic differences. To test this hypothesis, we analyzed and compared the global promoter occupancy of Sox2 in RU and RR cells using ChIP-chip. As detailed below, we found that the Sox2 gene promoter occupancy between RU and RR cells are mutually exclusive. Importantly, we identified a number of stem cell- or cancer-associated genes that were more highly expressed in RR cells.
Cell lines and materials
MCF7 and ZR751 parental cells were purchased from American Type Culture Collection (ATCC, Rockville, MD, USA). MCF7 and ZR751 parental cells, unsorted cells, RU (previously referred to as GFP Neg), and RR (previously referred to as GFP Pos) cells were cultured and derived as previously described . Triptolide was purchased from Sigma-Aldrich (T3652, Sigma-Aldrich Canada, Oakville, ON, Canada).
Soxtranscription activity reporter
The commercially available Sox2 transcription activity reporter is driven by a minimal CMV promoter followed by three tandem repeats of the Sox2 regulatory region 2 (SRR2), a sequence containing a Sox2 consensus sequence demonstrated to be bound by Sox2 in mouse and human embryonic stem cells .
ChIP (chromatin immunoprecipitation)-chip and ChIP-PCR
ChIP-chip was performed based on a previously described ChIP-PCR protocol . The starting material was scaled up four times, such that starting materials were four 15-cm plates of both MCF7 RU and RR cells, and four identical immunoprecipitations were performed for each condition (MCF7 RU and RR, IgG and Sox2 IPs). The resulting DNA was further purified using the QIAquick PCR Purification Kit (Qiagen Canada, Toronto, ON, Canada), lyophilized, and reconstituted in 10 μL of UltraPure DNase/RNase-free distilled water (Life Technologies, Burlington, ON, Canada). The DNA was subsequently amplified twice using the Sigma GenomePlex Complete Whole Genome Amplification Kit (#WGA2, Sigma-Aldrich Canada) using a published adapted protocol . ChIP-PCR was performed as previously described . ChIP input DNA was run on an agarose gel to check for microarray optimized DNA fragments of 200 to 1200 bp (Additional file 1: Figure S1A ). DNA samples were sent in two replicates to Roche Nimblegen ChIP-chip Microarray Services for quality assessment, and full service ChIP-chip microarray service and analysis. Briefly, DNA samples were hybridized to the Roche Nimblegen Human ChIP-chip 3x720K RefSeq Promoter array, with promoter tiling ranging from -3,200 to +800 relative to the transcription start site. The ChIP-chip microarray data have been submitted to the public repository Gene Expression Omnibus [GEO: GSE61703]. Primers for ChIP-PCR were designed to flank the promoter peaks identified by ChIP-chip analysis for each gene.
RNA extraction, cDNA synthesis, quantitative reverse transcription PCR (qPCR)
Total RNA extraction was performed with the Qiagen RNeasy Kit (Qiagen Canada) according to the manufacturer’s protocol: 1 μg of RNA was reverse transcribed using Oligo dT and Superscript II (Life Technologies) according to the manufacturer’s protocol. 1 μL of the resulting cDNA mixture was added to the Platinum SYBR Green qPCR SuperMix-UDG with Rox (Life Technologies) and amplified with target gene-specific primers. Please see Additional file 2: Table S1 for list of PrimerBank primer sequences ,. All genes of interest are normalized to glyceraldehyde-3-phosphate dehydrogenase (GAPDH) transcript expression levels except for the Triptolide experiments where 18S rRNA was used as the housekeeping gene for its superior stability.
Sox2 siRNAs (SMARTpool: ON-TARGETplus SOX2 siRNA, Dharmacon, Thermo Scientific, Waltham, MA, USA) or scrambled (Scr) siRNAs (ON-TARGETplus Non-targeting Pool, #477C20, Dharmacon, ThermoScientific) at 40 pmol per rxn (20 nM final concentration) and 5 μL of Lipofectamine RNAiMAX (Life Technologies) were added to 0.5 mL of OptiMEM media (Life Technologies) and reverse transfected to 800,000 cells in normal culture medium in a 6-well plate format. Cells were incubated with siRNAs for 72 hours before harvesting. Muc15 siRNA (#SI04331166, Qiagen Canada, and SMARTpool: ON-TARGETplus MUC15 siRNA, Dharmacon) was transfected in the same manner at 80 pmol and 200 pmol per rxn respectively (40 nM and 100 nM final concentration).
Western blot analyses were performed as previously described . All antibodies were diluted in 5% BSA in TBST: Sox2 (1:500, #2683-1, Epitomics, Burlingame, CA, USA), FlagM2 (1:1000, #F1804, Sigma-Aldrich), Muc15 (1:500, #ab98045, Abcam, Cambridge, UK), and vinculin (1:1000, #4650, Cell Signaling Technologies, Danvers, MA, USA). Vinculin was used as loading control for all western blots.
We transfected 3 μg of pcDNA-Flag-EV or pcDNA-Flag-Sox2 with 5 μL of Lipofectamine 2000 (Life Technologies) in 0.5 mL of OptiMEM media (Life Technologies) to 1.2 million MCF7 cells seeded the day before. Cells were incubated for 72 hours before harvesting.
Mammospheres were generated as previously described . Mammospheres were collected by centrifugation at 300 × g for 5 minutes and trypsinized before subjecting to trypan blue exclusion assay of mammosphere-derived cells.
Primary patient breast tumor cells isolation, lentiviral infections, fluorescence-activated cell sorting (FACS) purification
Patient material and clinical information were collected with full written consent from the patients and with approval by the University of Alberta Human Research Ethics Board, approved project ID Pro00044942. Fresh breast tumors were collected in cold 100% FBS and harvested within hours. We isolated breast tumor cells from fresh breast tumor tissues with no exposure to radiation therapy or chemotherapy. We harvested purified primary BC cells first by mechanical dissociation and then by using the Cancer Cell Isolation Kit (Panomics Solutions, Affymetrix, Santa Clara, CA, USA) as per manufacturer’s protocol. Cells were cultured in 10% RPMI medium for 48 hours before virus infection. We generated a new dual green fluorescent protein (GFP)/red fluorescent protein (RFP) lentiviral Sox2 reporter by replacing the puromycin resistance gene in the Sox2 reporter with the red fluorescent protein (RFP) gene. Isolated tumor cells were infected with our modified lentiviral Sox2 GFP-RFP dual-color reporter, SRR2-mCMV-GFP-EF1-RFP, twice 24 hours apart. RFP+ cells were gated to include only successfully infected primary breast tumor cells in subsequent analyses and experiments. Using flow cytometry, we analyzed and collected RFP+/GFP- (RU) and RFP+/GFP+ (RR) cells.
The paired Student’s t-test was used for statistical analysis of experiments throughout: *P <0.05; **P <0.01.
The Sox2-bound gene promoter regions are largely mutually exclusive between RU and RR cells
Using ChIP-chip, we queried the global promoter occupancy profile of Sox2 in the two phenotypically distinct cell subsets, namely RU and RR cells. Using a stringent threshold (a promoter array peak signal >2.0, compared to the input DNA signal) and a false discovery rate of <0.05, we found that Sox2 was bound to the promoter regions of 1,830 genes in RU cells and 456 genes in RR cells, with an overlap of only 62 genes between the two cell subsets (illustrated in Figure 1A). The complete RU and RR gene lists can be found in Additional file 3: Table S2. ChIP-chip gene promoter analyses are detailed in Additional file 1 (Supplementary materials and methods).
To understand the possible biological effects exerted by Sox2 in BC cells, we annotated the functions of the identified genes using the Protein Analysis THrough Evolutionary Relationships (PANTHER) Protein Class classification system software . As shown in Figure 1B, the biological functions associated with the identified genes are largely similar between the RU and RR cells, with the functions falling most frequently into the categories of hydrolases, nucleic acid binding, and receptors.
The RR gene list comprise markers associated with cancer stem cells
As we have previously shown that RR cells exhibit more tumorigenic and stem-like properties than RU cells  we hypothesized that the ChIP-chip gene list derived from RR cells will contain genes that are known to be associated with cancer stem cells. To test this hypothesis, we searched our RR gene lists for reported cancer stem cell markers, based on those described in two recent publications ,. We found that Sox2 was bound to the gene promoters of three established stem cell markers in solid tumors, including CD133 (PROM1), Lgr5 (GPR49), and Bmi-1 (BMI1). Importantly, these three genes were not on the RU gene list.
When we examined the remaining 453 genes identified in RR cells, we identified 12 additional genes that have been previously implicated in cancer initiation and/or progression (Table 1). These genes include FZD4 (the Wnt pathway) ,, PLAU (encoding metastasis-promoting protein urokinase plasminogen activator)  and ELF5 (a normal mammary stem/progenitor cell gene) . None of these 12 genes were found in the RU gene list and the majority of these genes (8 of 15) had a very high microarray signal >2.5. Interestingly, ANTXR1, also found in our RR gene list, which encodes anthrax toxin receptor-1, has just been recently reported as a stem cell gene important to the tumorigenesis of BC ,. Again, this gene was not found in the RU gene list.
Validation of the ChIP-chip data using ChIP-PCR
We then aimed to validate the observation that the gene promoters bound by Sox2 in RU and RR cells are largely non-overlapping. To do so, we employed ChIP-PCR and used two genes from the RR gene list that show relatively high microarray signals and robust mRNA expression in BC cells, namely GPR49 and MUC15 ,. The ChIP-PCR primers for these two genes were designed to flank the exact promoter locations specified by the ChIP-chip microarray probes. As shown in Figure 1C, in RR cells, we detected more robust Sox2 binding at both the GPR49 and MUC15 gene promoters than in the RU cells that showed barely detectable to no binding. These ChIP-PCR results support the validity of the ChIP-chip findings. To further validate our ChIP-PCR findings, we also pursued ChIP-qPCR analyses of six gene promoters of interest with high peak scores from the RR gene list, and validated that Sox2 was significantly more frequently bound to these promoters in the RR cells when compared to the RU cells (Figure 1C).
To further test if Sox2 binds to different sets of gene promoters between RU and RR cells, we performed ChIP-PCR to detect the binding of Sox2 to CCND1 (Cyclin D1) promoter, a direct Sox2 gene target previously shown by us and others ,. We found the interaction between Sox2 and the CCND1 gene promoter, but only in RR (data not shown, previously reported by us . We also validated our ChIP DNA by looking at several Sox2 target genes found in human ESCs previously described in the literature, including BCL2 and CDH1  As shown in Additional file 1: Figure S1B, we found that Sox2 showed significantly greater binding at the promoters of BCL2 and CDH1 in RR cells than in RU cells. Of note, CCDN1, BCL2, and CDH1 were not found in our ChIP-chip gene list, likely due to our very stringent analysis criteria, which were used to identify only the most frequently bound DNA sequences in BC cells.
RR cells express elevated levels of target genes compared to RU cells
We next asked if the differential Sox2 gene promoter occupancy between RU and RR cells correlates to significant differences in gene expression between these two cell subsets. Using qRT-PCR, we measured and compared the expression levels of the 15 genes of interest described in Table 1. As compared to RU cells, RR cells expressed significantly higher (2- to 5-fold) gene transcript levels of 14 out of these 15 genes (Figure 2). These results support our hypothesis that Sox2 mediates differential gene expression between RU and RR cells.
Overexpression of Soxup-regulates target genes in RR cells but not RU cells
To demonstrate the direct role of Sox2 in contributing to the differential gene expression between RU and RR, we examined if enforced expression of Sox2 in MCF7 cells results in significant alterations of their expressions. For the purpose of this study, we chose 7 of the 15 genes, based on their relatively high ChIP-chip peak scores, including PLXNA2, FZD4, MUC15, PLAU, ELF5, GPR49 and PROM1. As shown in Figure 3A, with transient transfection of Sox2 into RR cells, all seven genes examined showed a significant increase in their transcript levels in RR cells (3- to 7-fold); conversely, RU cells showed no significant alterations of any of these seven genes.
SiRNA knockdown of Soxdownregulates target genes
Next, we examined if siRNA knockdown of Sox2 also can modulate the expression of the seven target genes tested. As shown in Figure 3B, the efficiency of the knockdown was demonstrated by western blotting and quantitative RT-PCR. We found that Sox2 siRNAs significantly downregulated these target genes in RR cells. Surprisingly, the same treatment also significantly downregulated the expression of these seven genes in RU cells. Similar findings were also observed in MCF7 parental cells, which predominantly comprise RU cells (Additional file 1: Figure S2). As Sox2 did not induce an increase in the expression of Sox2 target genes in RU cells (Figure 3A), we hypothesized that the downregulation of Sox2 target genes in RU cells induced by Sox2 siRNA was mediated in a transcription-independent manner. If this is the case, the gene transcripts in RR cells are expected to be more sensitive to transcription inhibition than those in RU cells. In keeping with this concept, the addition of the transcription inhibitor, Triptolide, significantly decreased the transcript level of PROM1 (CD133) in RR cells but paradoxically increased that in RU cells (Figure 3C).
Mucin-15, a novel Sox2 target, contributes to mammosphere formation
To further support the concept that Sox2 contributes to tumorigenesis and stemness in BC by upregulating these stem cell- or cancer cell-associated genes, we examined the oncogenic effects of Mucin-15 (Muc15), which has not been previously shown to be a Sox2 downstream target. While Muc15 has been shown to play a key role in increasing invasiveness and tumorigenic capacity in colon cancer , it has not been linked to BC. As shown in Figures 2 and 4A, Muc15 was more highly expressed at the mRNA and protein levels in RR cells, as compared to RU cells. Furthermore, as shown above, overexpression or knockdown of Sox2 significantly modulated the expression of Muc15. As shown in Figure 4B, knockdown of Muc15 using siRNA significantly decreased the number of mammospheres formed from MCF7 unsorted cells, which comprise natural proportions of RU and RR subsets. Furthermore, using trypan blue exclusion assay, we found that siRNA knockdown of Muc15 significantly reduced the number of viable cells derived from the mammospheres (Figure 4B). The same experiment was repeated using four pooled unique siRNA sequences and we observed the same results, with Muc15 knockdown verified (Additional file 1: Figure S3).
RR cells derived from primary patient breast tumors exhibit elevated tumorigenic properties and expression of target genes
Lastly, we examined if BC cells derived from patient samples displayed similar findings to the MCF7 cells. Due to the relatively small number of tumor cells available, and the relatively low proportions of RR cells, we modified our Sox2 reporter such that it carried two signals, with the expression of RFP indicating successful infection with the viral vector, and the GFP signal indicating Sox2 reporter activity (detailed in main Methods section). Only cells expressing RFP but not GFP were regarded as RU cells, whereas those lacking both RFP and GFP were excluded from the analysis. Results from 19 primary BC tumors are summarized in Table 2. All 19 samples contained a detectable subset of RR cells, and the size of this population ranged from 0.3% to 23.8%. Interestingly, estrogen receptor-negative tumors (n = 3) had a significantly lower proportion of RR cells (p = 0.001). Functional studies were performed in eight samples in total. As shown in Figure 5A, RR cells were more efficient in forming colonies on methylcellulose agar in four out of four patient cells sampled. Importantly, as we gated our cells using RFP and GFP expression, we demonstrate that the RFP+ GFP- cells were healthy in culture (Additional file 1: Figure S4A-B). Under a fluorescence microscope, the cells were confirmed to be RFP+ (data not shown). As shown in Additional file 1: Figure S4B-C, RU and RR cells derived from patient samples had a similar Sox2 protein expression level in the nuclei, suggesting that the differences observed are not simply due to a lack of Sox2 protein or Sox2 nuclear localization in RU cells. Using fresh primary patient samples, we went on to test if RU and RR cells also differ in the expression of Sox2 downstream targets. Due to the relatively small number of primary samples available for testing, we chose three genes, including PROM1 (CD133), GPR49 (LGR5), and MUC15, based on the fact that the expression of these genes were amongst the most responsive to modulation of Sox2 (Figure 3A and B). As shown in Figure 5B, in a total of seven fresh primary patient samples, we detected higher expression of these three genes in patient RR cells as compared to their RU counterparts, although statistical analysis was not possible for all due to limitations in patient materials. Further, some patient samples did not contain enough RNA for analysis for all genes. Nevertheless, the overall findings from patient samples appear to mirror those in MCF7 cells.
We recently identified two Sox2-expressing, phenotypically distinct cell subsets in BC cells, separated based on their differential response to a Sox2 transcription activity reporter, with RR cells showing higher tumorigenicity and more stem-like features relative to RU cells . In the same study, we also found that these phenotypic differences are dependent on Sox2, as siRNA knockdown of Sox2 abrogates many of these phenotypic differences . Because Sox2 is a transcription factor, we hypothesized that Sox2 contributes to the phenotypic differences between RU and RR cells by mediating differential gene expression. To test this hypothesis, we compared the Sox2 gene occupancy in RU cells with that of RR cells. Importantly, we found a largely mutually exclusive Sox2 promoter occupancy between these two cell subsets. Furthermore, there were a number of cancer- or stem-cell associated genes that are only found in the RR gene lists. Experiments using enforced expression or siRNA knockdown of Sox2 support the direct role of Sox2 in regulating these genes. The biological significance of our findings is supported by our results generated from the use of patient samples. Taken together, we believe that the overall findings lend support to our hypothesis.
Although aberrant Sox2 expression is well-documented in cancer, its mechanism of action in the regulation of downstream targets is incompletely understood. Currently, with the exception of CCND1 (encoding Cyclin D1) , no other gene has been identified as a direct downstream target of Sox2 in BC. Nevertheless, a few Sox2 downstream gene targets have been reported in other cancer types, including PROM1 (encoding CD133) in human lung cancer cells  and ITGA6 (encoding CD49f) in human mesenchymal stem cells . Regarding the functional importance of Sox2 in cancer, an exciting finding from our ChIP-chip study is that Sox2 was bound to the promoters of many cancer- and stem cell-associated genes in RR cells. This finding correlates well with the prevailing concept that the expression of embryonic stem cell markers in cancer cells results in stem-like features, which are often associated with an aggressive clinical course and treatment resistance ,. We believe that our finding of Sox2 regulating an array of cancer- and stem cell-associated genes provides a mechanistic explanation as to how Sox2 enhances stemness and tumorigenesis in cancer cell subsets. The importance of stem cell markers in identifying cancer stem cells, including Frizzled-4, Lgr5, and CD133 have previously been demonstrated ,,, and here our data suggest that their expressions may be dependent on common precursor protein Sox2. Furthermore, we have demonstrated that CD133 and Lgr5 mRNA transcripts were also upregulated in primary tumor-derived RR cells from patients.
Importantly, it should be noted that the identification of our list of 15 novel Sox2 targets were hand-picked by us using a manual search approach. As we were most interested in better understanding how Sox2 contributes to BC and/or BC stem cell biology, we chose genes with published roles in that context. As a result, we have discovered that Sox2 does regulate an intriguing list of genes in the RR cells, but this does not exclude the possibility that other important cancer and/or stem cell genes exist in our ChIP-chip lists. Additionally, as the ChIP-chip assay is limited by the detection of hybridization of our Sox2-bound DNA samples to the microarray, our list certainly does not exhaust all the possible promoter interactions of Sox2 in BC cells. Importantly, we have done motif analyses on the ChIP-chip data, and have confirmed that Sox2 motifs previously published by others are enriched in our Sox2 ChIP DNA from both subsets .
We hypothesize that Sox2 in RU and RR cells are biochemically distinct, allowing for differential transcription activation ability at unique promoter regions. The RU cells exhibit no transcription activity as reported by our Sox2 reporter, and here we have shown that Sox2 overexpression did not transactivate the RR ChIP-chip promoters. Moreover, we have shown by ChIP-chip and conventional ChIP-PCR that Sox2 does not occupy the same promoters in RU and RR cells. These results suggest that Sox2 does not interact with these promoters in RU cells the same way as in RR cells. Conversely, Sox2 in RU cells binds to its own large cohort of gene promoters. This suggests multiple possibilities for the role of Sox2 at the RU gene promoters: 1) Sox2 could be suppressing gene expression of these genes as we have recently reported ; 2) Sox2 is transcriptionally active in RU cell gene promoters but did not transactivate luciferase or GFP expression from the reporter due to discrepancies between the reporter and gene promoters in Sox2 consensus binding sequences and/or adjacent sequences that can recruit other co-factors; 3) Sox2 occupancy at these promoters serves as a positive or negative facilitator to other transcriptional co-factor binding and/or activation; and 4) Sox2 is non-functional at these promoter regions due to an absent co-factor or post-translational modification that is present in RR cells.
While we found that Sox2 is directly involved in regulating the expression of its target genes in RR cells, the finding that siRNA knockdown of Sox2 decreased gene transcript expression in RU cells is a rather unexpected finding. From our previous studies, we found that Sox2 exists in the cytoplasm , and it can potentially carry out functions related to post-transcriptional modifications and/or translational modulations. One possible explanation is that Sox2 regulates the expression of these genes by non-transcriptional mechanisms. It is possible that Sox2 can prolong the integrity and half-life of specific gene transcripts, or it functions as a translation factor. To examine the contributions of transcriptional and non-transcriptional mechanisms in RU and RR cells, we used transcription inhibitor Triptolide. In RR cells, we found that Sox2 target PROM1 transcripts were sensitive to the treatment and the mRNA levels decreased with increasing concentrations, supporting the hypothesis that Sox2 is transcriptionally activating PROM1. In RU cells, we did not observe decreased PROM1 transcript levels with transcription inhibitor treatment, suggesting that Sox2 in RU cells may have a distinct regulatory mechanism for Sox2 target PROM1.
We have focused on Muc15 in our studies as it is a new putative onco-protein, consistently highly expressed in RR cells, and responsive to Sox2 regulation. In particular, Muc15 is of interest to us as it is relatively unknown in the cancer biology of any tissue. Muc15 is a highly glycosylated extracellular mucin protein previously reported to be expressed in normal epithelial cells, including the breast, but elevated in tumor cell populations ,,-. In this report, we are the first to identify very high Muc15 expression in BC cells. Importantly, we detected increased Muc15 mRNA transcript and protein levels in RR cells compared to RU cells in cell lines and primary patient samples. Muc15 was previously demonstrated to promote oncogenesis in colon cancer cells in vitro and in vivo . Thus, our ChIP-chip study is a good resource for novel putative therapeutic BC targets.
We showed that the patient RU and RR cells have distinct phenotypes as demonstrated in an anchorage-independent methylcellulose colony formation assay and underlying biology as determined by qPCR. Importantly, we have confirmed that RU cells from patients, although reporter unresponsive, do express nuclear Sox2 as detected by immunohistochemistry techniques. We have also uncovered potential mechanisms underlying the more tumorigenic RR cells as the patient RR cells exhibited higher expression levels of Sox2 target genes, PROM1, GPR49, and MUC15 transcripts. Thus, the response of BC cells to the Sox2 transcription activity reporter has distinguished primary patient and cultured cell lines cancer-cell subpopulations with distinct phenotypic and molecular features.
Taken together, we have shown that Sox2 behaves heterogeneously in breast tumor cell populations. Sox2 is strongly bound to a subset of cancer and stem cell gene promoters and can upregulate the corresponding gene transcripts in RR cells but not in RU cells. Importantly, we have identified a novel Sox2 target Muc15 that is important for mammosphere formation, and is also upregulated in the tumorigenic RR cells derived from primary patient breast tissue samples. In summary, we depict in a schematic diagram where Sox2 in RR cells interacts with DNA, and/or transcriptionally activates promoters differently compared to Sox2 in RU cells (Figure 6).
KJ conceived and designed the research plan, optimized and performed the chromatin immunoprecipitation, analyzed the microarray results, performed all the cell line work and experiments, assisted with some of the patient sample work and experiments, analyzed data, and wrote the manuscript. PW, NG, KG, FW, and XY performed the majority of the patient sample work and experiments. AA provided assistance in performing and analyzing the qPCR experiments. GB and TPM collected the patient tumours, provided clinical data, and gave intellectual input. BSA provided intellectual input and critical reading of the manuscript. RL conceived and designed the research plan and wrote the manuscript. All authors read and approved the final manuscript.
bovine serum albumin
chromatin immunoprecipitation promoter microarray chip analysis
fluorescence-activated cell sorting
fetal bovine serum
green fluorescent protein
Protein Analysis THrough Evolutionary Relationships
quantitative reverse transcription polymerase chain reaction
red fluorescent protein
sex-determining region Y (SRY)-box binding protein-2
Sox2 regulatory region 2
Takahashi K, Tanabe K, Ohnuki M, Narita M, Ichisaka T, Tomoda K, Yamanaka S: Induction of pluripotent stem cells from adult human fibroblasts by defined factors. Cell. 2007, 131: 861-872. 10.1016/j.cell.2007.11.019.
Boyer LA, Lee TI, Cole MF, Johnstone SE, Levine SS, Zucker JP, Guenther MG, Kumar RM, Murray HL, Jenner RG, Gifford DK, Melton DA, Jaenisch R, Young RA: Core transcriptional regulatory circuitry in human embryonic stem cells. Cell. 2005, 122: 947-956. 10.1016/j.cell.2005.08.020.
Sarkar A, Hochedlinger K: The sox family of transcription factors: versatile regulators of stem and progenitor cell fate. Cell Stem Cell. 2013, 12: 15-30. 10.1016/j.stem.2012.12.007.
Arnold K, Sarkar A, Yram MA, Polo JM, Bronson R, Sengupta S, Seandel M, Geijsen N, Hochedlinger K: Sox2(+) adult stem and progenitor cells are important for tissue regeneration and survival of mice. Cell Stem Cell. 2011, 9: 317-329. 10.1016/j.stem.2011.09.001.
Roy S, Gascard P, Dumont N, Zhao J, Pan D, Petrie S, Margeta M, Tlsty TD: Rare somatic cells from human breast tissue exhibit extensive lineage plasticity. Proc Natl Acad Sci U S A. 2013, 110: 4598-4603. 10.1073/pnas.1218682110.
Wang Y, Dong J, Li D, Lai L, Siwko S, Li Y, Liu M: Lgr4 regulates mammary gland development and stem cell activity through the pluripotency transcription factor Sox2. Stem Cells. 2013, 31: 1921-1931. 10.1002/stem.1438.
Leis O, Eguiara A, Lopez-Arribillaga E, Alberdi MJ, Hernandez-Garcia S, Elorriaga K, Pandiella A, Rezola R, Martin AG: Sox2 expression in breast tumors and activation in breast cancer stem cells. Oncogene. 2012, 31: 1354-1365. 10.1038/onc.2011.338.
Chen Y, Shi L, Zhang L, Li R, Liang J, Yu W, Sun L, Yang X, Wang Y, Zhang Y, Shang Y: The molecular mechanism governing the oncogenic potential of SOX2 in breast cancer. J Biol Chem. 2008, 283: 17969-17978. 10.1074/jbc.M802917200.
Rudin CM, Durinck S, Stawiski EW, Poirier JT, Modrusan Z, Shames DS, Bergbower EA, Guan Y, Shin J, Guillory J, Rivers CS, Foo CK, Bhatt D, Stinson J, Gnad F, Haverty PM, Gentleman R, Chaudhuri S, Janakiraman V, Jaiswal BS, Parikh C, Yuan W, Zhang Z, Koeppen H, Wu TD, Stern HM, Yauch RL, Huffman KE, Paskulin DD, Illei PB, et al: Comprehensive genomic analysis identifies SOX2 as a frequently amplified gene in small-cell lung cancer. Nat Genet. 2012, 44: 1111-1116. 10.1038/ng.2405.
Annovazzi L, Mellai M, Caldera V, Valente G, Schiffer D: SOX2 expression and amplification in gliomas and glioma cell lines. Cancer Genomics Proteomics. 2011, 8: 139-147.
Zhang J, Chang DY, Mercado-Uribe I, Liu J: Sex-determining region Y-box 2 expression predicts poor prognosis in human ovarian carcinoma. Hum Pathol. 2012, 43: 1405-1412. 10.1016/j.humpath.2011.10.016.
Basu-Roy U, Seo E, Ramanathapuram L, Rapp TB, Perry JA, Orkin SH, Mansukhani A, Basilico C: Sox2 maintains self renewal of tumor-initiating cells in osteosarcomas. Oncogene. 2012, 31: 2270-2282. 10.1038/onc.2011.405.
Neumann J, Bahr F, Horst D, Kriegl L, Engel J, Luque RM, Gerhard M, Kirchner T, Jung A: SOX2 expression correlates with lymph-node metastases and distant spread in right-sided colon cancer. BMC Cancer. 2011, 11: 518-10.1186/1471-2407-11-518.
Laga AC, Zhan Q, Weishaupt C, Ma J, Frank MH, Murphy GF: SOX2 and nestin expression in human melanoma: an immunohistochemical and experimental study. Exp Dermatol. 2011, 20: 339-345. 10.1111/j.1600-0625.2011.01247.x.
Lengerke C, Fehm T, Kurth R, Neubauer H, Scheble V, Muller F, Schneider F, Petersen K, Wallwiener D, Kanz L, Fend F, Perner S, Bareiss PM, Staebler A: Expression of the embryonic stem cell marker SOX2 in early-stage breast carcinoma. BMC Cancer. 2011, 28: 11:42-
Hagerstrand D, He X, Bradic Lindh M, Hoefs S, Hesselager G, Ostman A, Nister M: Identification of a SOX2-dependent subset of tumor- and sphere-forming glioblastoma cells with a distinct tyrosine kinase inhibitor sensitivity profile. Neuro Oncol. 2011, 13: 1178-1191. 10.1093/neuonc/nor113.
Nakatsugawa M, Takahashi A, Hirohashi Y, Torigoe T, Inoda S, Murase M, Asanuma H, Tamura Y, Morita R, Michifuri Y, Kondo T, Hasegawa T, Takahashi H, Sato N: SOX2 is overexpressed in stem-like cells of human lung adenocarcinoma and augments the tumorigenicity. Lab Invest. 2011, 91: 1796-1804. 10.1038/labinvest.2011.140.
Wu F, Zhang J, Wang P, Ye X, Jung K, Bone KM, Pearson JD, Ingham RJ, McMullen TP, Ma Y, Lai R: Identification of two novel phenotypically distinct breast cancer cell subsets based on Sox2 transcription activity. Cell Signal. 2012, 24: 1989-1998. 10.1016/j.cellsig.2012.07.008.
Zhang Y, Eades G, Yao Y, Li Q, Zhou Q: Estrogen receptor alpha signaling regulates breast tumor-initiating cells by down-regulating miR-140 which targets the transcription factor SOX2. J Biol Chem. 2012, 287: 41514-41522. 10.1074/jbc.M112.404871.
Singh S, Trevino J, Bora-Singhal N, Coppola D, Haura E, Altiok S, Chellappan SP: EGFR/Src/Akt signaling modulates Sox2 expression and self-renewal of stem-like side-population cells in non-small cell lung cancer. Mol Cancer. 2012, 11: 73-10.1186/1476-4598-11-73.
Lin CW, Liao MY, Lin WW, Wang YP, Lu TY, Wu HC: Epithelial cell adhesion molecule regulates tumor initiation and tumorigenesis via activating reprogramming factors and epithelial-mesenchymal transition gene expression in colon cancer. J Biol Chem. 2012, 287: 39449-39459. 10.1074/jbc.M112.386235.
Chen YL, Wang SY, Liu RS, Wang HE, Chen JC, Chiou SH, Chang CA, Lin LT, Tan DT, Lee YJ: Remnant living cells that escape cell loss in late-stage tumors exhibit cancer stem cell-like characteristics. Cell Death Dis. 2012, 3: e399-10.1038/cddis.2012.136.
Han X, Fang X, Lou X, Hua D, Ding W, Foltz G, Hood L, Yuan Y, Lin B: Silencing SOX2 induced mesenchymal-epithelial transition and its expression predicts liver and lymph node metastasis of CRC patients. PLoS One. 2012, 7: e41335-10.1371/journal.pone.0041335.
Lu X, Deng Q, Li H, Suo Z: Altered characteristics of cancer stem/initiating cells in a breast cancer cell line treated with persistent 5-FU chemotherapy. Exp Ther Med. 2011, 2: 821-826.
Wang X, Liang Y, Chen Q, Xu HM, Ge N, Luo RZ, Shao JY, He Z, Zeng YX, Kang T, Yun JP, Xie F: Prognostic significance of SOX2 expression in nasopharyngeal carcinoma. Cancer Invest. 2012, 30: 79-85. 10.3109/07357907.2011.630049.
Ben-Porath I, Thomson MW, Carey VJ, Ge R, Bell GW, Regev A, Weinberg RA: An embryonic stem cell-like gene expression signature in poorly differentiated aggressive human tumors. Nat Genet. 2008, 40: 499-507. 10.1038/ng.127.
Gomez-Mateo Mdel C, Piqueras M, Pahlman S, Noguera R, Navarro S: Prognostic value of SOX2 expression in neuroblastoma. Genes Chromosomes Cancer. 2011, 50: 374-377. 10.1002/gcc.20859.
Sholl LM, Barletta JA, Yeap BY, Chirieac LR, Hornick JL: Sox2 protein expression is an independent poor prognostic indicator in stage I lung adenocarcinoma. Am J Surg Pathol. 2010, 34: 1193-1198. 10.1097/PAS.0b013e3181e5e024.
Simoes BM, Piva M, Iriondo O, Comaills V, Lopez-Ruiz JA, Zabalza I, Mieza JA, Acinas O, Vivanco MD: Effects of estrogen on the proportion of stem cells in the breast. Breast Cancer Res Treat. 2011, 129: 23-35. 10.1007/s10549-010-1169-4.
Tomioka M, Nishimoto M, Miyagi S, Katayanagi T, Fukui N, Niwa H, Muramatsu M, Okuda A: Identification of Sox-2 regulatory region which is under the control of Oct-3/4-Sox-2 complex. Nucleic Acids Res. 2002, 30: 3202-3213. 10.1093/nar/gkf435.
O'Geen H, Nicolet CM, Blahnik K, Green R, Farnham PJ: Comparison of sample preparation methods for ChIP-chip assays. Biotechniques. 2006, 41: 577-580. 10.2144/000112268.
Spandidos A, Wang X, Wang H, Seed B: PrimerBank: a resource of human and mouse PCR primer pairs for gene expression detection and quantification. Nucleic Acids Res. 2010, 38: D792-D799. 10.1093/nar/gkp1005.
Wang X, Spandidos A, Wang H, Seed B: PrimerBank: a PCR primer database for quantitative gene expression analysis, 2012 update. Nucleic Acids Res. 2012, 40: D1144-D1149. 10.1093/nar/gkr1013.
Zhang J, Wang P, Wu F, Li M, Sharon D, Ingham RJ, Hitt M, McMullen TP, Lai R: Aberrant expression of the transcriptional factor Twist1 promotes invasiveness in ALK-positive anaplastic large cell lymphoma. Cell Signal. 2012, 24: 852-858. 10.1016/j.cellsig.2011.11.020.
Mi H, Muruganujan A, Thomas PD: PANTHER in 2013: modeling the evolution of gene function, and other gene attributes, in the context of phylogenetic trees. Nucleic Acids Res. 2013, 41: D377-D386. 10.1093/nar/gks1118.
Medema JP: Cancer stem cells: the challenges ahead. Nat Cell Biol. 2013, 15: 338-344. 10.1038/ncb2717.
Alison MR, Islam S, Wright NA: Stem cells in cancer: instigators and propagators?. J Cell Sci. 2010, 123: 2357-2368. 10.1242/jcs.054296.
Gupta S, Iljin K, Sara H, Mpindi JP, Mirtti T, Vainio P, Rantala J, Alanen K, Nees M, Kallioniemi O: FZD4 as a mediator of ERG oncogene-induced WNT signaling and epithelial-to-mesenchymal transition in human prostate cancer cells. Cancer Res. 2010, 70: 6735-6745. 10.1158/0008-5472.CAN-10-0244.
Jin X, Jeon HY, Joo KM, Kim JK, Jin J, Kim SH, Kang BG, Beck S, Lee SJ, Kim JK, Park AK, Park WY, Choi YJ, Nam DH, Kim H: Frizzled 4 regulates stemness and invasiveness of migrating glioma cells established by serial intracranial transplantation. Cancer Res. 2011, 71: 3066-3075. 10.1158/0008-5472.CAN-10-1495.
Nielsen TO, Andrews HN, Cheang M, Kucab JE, Hsu FD, Ragaz J, Gilks CB, Makretsov N, Bajdik CD, Brookes C, Neckers LM, Evdokimova V, Huntsman DG, Dunn SE: Expression of the insulin-like growth factor I receptor and urokinase plasminogen activator in breast cancer is associated with poor survival: potential for intervention with 17-allylamino geldanamycin. Cancer Res. 2004, 64: 286-291. 10.1158/0008-5472.CAN-03-1242.
Oakes SR, Naylor MJ, Asselin-Labat ML, Blazek KD, Gardiner-Garden M, Hilton HN, Kazlauskas M, Pritchard MA, Chodosh LA, Pfeffer PL, Lindeman GJ, Visvader JE, Ormandy CJ: The Ets transcription factor Elf5 specifies mammary alveolar cell fate. Genes Dev. 2008, 22: 581-586. 10.1101/gad.1614608.
Chen D, Bhat-Nakshatri P, Goswami C, Badve S, Nakshatri H: ANTXR1, a stem cell-enriched functional biomarker, connects collagen signaling to cancer stem-like cells and metastasis in breast cancer. Cancer Res. 2013, 73: 5821-5833. 10.1158/0008-5472.CAN-13-1080.
Chaudhary A, Hilton MB, Seaman S, Haines DC, Stevenson S, Lemotte PK, Tschantz WR, Zhang XM, Saha S, Fleming T, St Croix B: TEM8/ANTXR1 blockade inhibits pathological angiogenesis and potentiates tumoricidal responses against multiple cancer types. Cancer Cell. 2012, 21: 212-226. 10.1016/j.ccr.2012.01.004.
Pallesen LT, Berglund L, Rasmussen LK, Petersen TE, Rasmussen JT: Isolation and characterization of MUC15, a novel cell membrane-associated mucin. Eur J Biochem. 2002, 269: 2755-2763. 10.1046/j.1432-1033.2002.02949.x.
Oskarsson T, Acharyya S, Zhang XH, Vanharanta S, Tavazoie SF, Morris PG, Downey RJ, Manova-Todorova K, Brogi E, Massague J: Breast cancer cells produce tenascin C as a metastatic niche component to colonize the lungs. Nat Med. 2011, 17: 867-874. 10.1038/nm.2379.
Medeiros RB, Papenfuss KJ, Hoium B, Coley K, Jadrich J, Goh SK, Elayaperumal A, Herrera JE, Resnik E, Ni HT: Novel sequential ChIP and simplified basic ChIP protocols for promoter co-occupancy and target gene identification in human embryonic stem cells. BMC Biotechnol. 2009, 9: 59-10.1186/1472-6750-9-59.
Huang J, Che MI, Huang YT, Shyu MK, Huang YM, Wu YM, Lin WC, Huang PH, Liang JT, Lee PH, Huang MC: Overexpression of MUC15 activates extracellular signal-regulated kinase 1/2 and promotes the oncogenic potential of human colon cancer cells. Carcinogenesis. 2009, 30: 1452-1458. 10.1093/carcin/bgp137.
Iida H, Suzuki M, Goitsuka R, Ueno H: Hypoxia induces CD133 expression in human lung cancer cells by up-regulation of OCT3/4 and SOX2. Int J Oncol. 2012, 40: 71-79.
Yu KR, Yang SR, Jung JW, Kim H, Ko K, Han DW, Park SB, Choi SW, Kang SK, Scholer H, Kang KS: CD49f enhances multipotency and maintains stemness through the direct regulation of OCT4 and SOX2. Stem Cells. 2012, 30: 876-887. 10.1002/stem.1052.
Apostolou P, Toloudi M, Chatziioannou M, Ioannou E, Papasotiriou I: Cancer stem cells stemness transcription factors expression correlates with breast cancer disease stage. Curr Stem Cell Res Ther. 2012, 7: 415-419. 10.2174/157488812804484639.
Sato T, Clevers H: Growing self-organizing mini-guts from a single intestinal stem cell: mechanism and applications. Science. 2013, 340: 1190-1194. 10.1126/science.1234852.
Singh SK, Hawkins C, Clarke ID, Squire JA, Bayani J, Hide T, Henkelman RM, Cusimano MD, Dirks PB: Identification of human brain tumor initiating cells. Nature. 2004, 432: 396-401. 10.1038/nature03128.
Wu F, Ye X, Wang P, Jung K, Wu C, Douglas D, Kneteman N, Bigras G, Ma Y, Lai R: Sox2 suppresses the invasiveness of breast cancer cells via a mechanism that is dependent on Twist1 and the status of Sox2 transcription activity. BMC Cancer. 2013, 13: 317-10.1186/1471-2407-13-317.
Wang RY, Chen L, Chen HY, Hu L, Li L, Sun HY, Jiang F, Zhao J, Liu GM, Tang J, Chen CY, Yang YC, Chang YX, Liu H, Zhang J, Yang Y, Huang G, Shen F, Wu MC, Zhou WP, Wang HY: MUC15 inhibits dimerization of EGFR and PI3K-AKT signaling and is associated with aggressive hepatocellular carcinomas in patients. Gastroenterology. 2013, 145: 1436-1448. 10.1053/j.gastro.2013.08.009.
Nam KH, Noh TW, Chung SH, Lee SH, Lee MK, Hong SW, Chung WY, Lee EJ, Park CS: Expression of the membrane mucins MUC4 and MUC15, potential markers of malignancy and prognosis, in papillary thyroid carcinoma. Thyroid. 2011, 21: 745-750. 10.1089/thy.2010.0339.
Riker AI, Enkemann SA, Fodstad O, Liu S, Ren S, Morris C, Xi Y, Howell P, Metge B, Samant RS, Shevde LA, Li W, Eschrich S, Daud A, Ju J, Matta J: The gene expression profiles of primary and metastatic melanoma yields a transition point of tumor progression and metastasis. BMC Med Genomics. 2008, 1: 13-8794-1-13-10.1186/1755-8794-1-13.
This work was financially supported by grants from the Canadian Institutes of Health Research (CIHR), Alberta Cancer Foundation, and Canadian Breast Cancer Foundation awarded to RL. KJ is a recipient of the Izaak Walton Killam Memorial Scholarship, CIHR Vanier Canada Graduate Scholarships, and the University of Alberta President’s Doctoral Prize of Distinction. FW was awarded the Alberta Cancer Foundation Cancer Research Postdoctoral Fellowship. AA is funded by the Taibah University, Al Madina Al Munawarrah, Saudi Arabia. We acknowledge the work of Dr Theodore Perkins and Gareth Palidwor, University of Ottawa for their work in the ChIP-chip motif analyses. We thank Kathleen Bone for proofreading.
The authors declare that they have no competing interests.
Electronic supplementary material
Additional file 2: Table S1.: Primer sequences used for chromatin immunoprecipitation (ChIP)-PCR and quantitative (q)-PCR. All ChIP-PCR primers were designed using Primer 3. All q-PCR primers were designed using PrimerBank, Massachusetts General Hospital, Harvard Medical School (http://pga.mgh.harvard.edu/primerbank/). (XLSX 11 KB)
Additional file 3: Table S2.: MCF7 reporter unresponsive (RU) and reporter responsive (RR) gene lists with chromatin immunoprecipitation (ChIP)-chip promoter microarray peak score >2.0-fold compared to input DNA (1830 + 456 genes). (XLSX 141 KB)
Authors’ original submitted files for images
About this article
Cite this article
Jung, K., Wang, P., Gupta, N. et al. Profiling gene promoter occupancy of Sox2 in two phenotypically distinct breast cancer cell subsets using chromatin immunoprecipitation and genome-wide promoter microarrays. Breast Cancer Res 16, 470 (2014). https://doi.org/10.1186/s13058-014-0470-2
- Breast Cancer
- Green Fluorescent Protein
- Cancer Stem Cell
- Stem Cell Marker