Genome-wide mapping of FOXM1 binding reveals co-binding with estrogen receptor alpha in breast cancer cells
Genome Biology volume 14, Article number: R6 (2013)
The forkhead transcription factor FOXM1 is a key regulator of the cell cycle. It is frequently over-expressed in cancer and is emerging as an important therapeutic target. In breast cancer FOXM1 expression is linked with estrogen receptor (ERα) activity and resistance to endocrine therapies, with high levels correlated with poor prognosis. However, the precise role of FOXM1 in ER positive breast cancer is not yet fully understood.
The study utilizes chromatin immunoprecipitation followed by high-throughput sequencing to map FOXM1 binding in both ERα-positive and -negative breast cancer cell lines. The comparison between binding site distributions in the two cell lines uncovered a previously undescribed relationship between binding of FOXM1 and ERα. Further molecular analyses demonstrated that these two factors can bind simultaneously at genomic sites and furthermore that FOXM1 regulates the transcriptional activity of ERα via interaction with the coactivator CARM1. Inhibition of FOXM1 activity using the natural product thiostrepton revealed down-regulation of a set of FOXM1-regulated genes that are correlated with patient outcome in clinical breast cancer samples.
These findings reveal a novel role for FOXM1 in ERα transcriptional activity in breast cancer and uncover a FOXM1-regulated gene signature associated with ER-positive breast cancer patient prognosis.
The forkhead transcription factor FOXM1 is a key regulator of the cell cycle [1, 2] critical for the G1 to S phase transition and G2 to M progression . Expression of FOXM1 is essential for mitotic spindle assembly and correct chromosome segregation with depletion leading to mitotic catastrophe and cell cycle arrest . FOXM1 is also known to regulate the expression of genes involved in angiogenesis , metastasis  and response to oxidative stress and DNA damage [7, 8]. Overexpression of FOXM1 has been reported in many types of cancer  and is correlated with poor prognosis [10, 11]. Aberrant FOXM1 expression is an early event in oncogenesis , possibly acting as an initiating factor  and has been associated with genomic instability .
Breast cancer is one of the leading causes of cancer mortality in women and numerous studies have shown a correlation between FOXM1 expression and breast cancer progression [4, 14, 15], suggesting that FOXM1 is a potential prognostic breast tumor marker . FOXM1 expression in breast cancer was found to correlate with levels of YWHAZ, a member of the 14-3-3 family of proteins  and also with HER2 status [15, 16]. Meta-analysis of gene expression data from breast cancer patient studies identified FOXM1 as one of 117 genes comprising a gene expression signature predictive of survival . FOXM1 over-expression has also been linked with drug resistance in breast cancer chemotherapy [19, 20] and therefore poor clinical prognosis.
Approximately 70% of breast cancers are estrogen receptor (ERα)-positive and there is increasing evidence to suggest that ERα and FOXM1 act as co-regulators. FOXM1 and ERα regulate the expression of each other in a positive cross-regulatory loop [21, 22]. FOXM1 has previously been identified as an ERα-responsive gene  and has been suggested to act as a prognostic marker in endocrine-positive cancers . Furthermore, resistance to anti-estrogen treatment has been correlated with increased FOXM1 expression .
We investigated the relationship between FOXM1 and ERα in breast cancer by mapping global FOXM1 binding in an ERα-positive cell line (MCF7) and an ERα-negative cell line (MDA-MB-231) using chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq). We show that there are cell-line dependent patterns of FOXM1 binding. We identify a common set of FOXM1 binding sites in the promoter regions of cell cycle-regulating genes but additionally in MCF7 cells; the majority of binding is located in intronic and intragenic regions with a high concordance to ERα binding, similar to the distribution of FOXA1 , another forkhead factor. These data suggest a distinct role for FOXM1 in different cellular contexts.
FOXM1 binding overlaps with ER binding genome-wide
As FOXM1 has been implicated as an important transcription factor in breast cancer, we mapped FOXM1 binding genome-wide using ChIP-seq in asynchronous MCF7 cells to determine the regulatory regions bound by FOXM1. Four biological replicates were conducted in MCF7 cells, resulting in 21,029 FOXM1 discrete binding events detected in at least two replicates. FOXM1 is known to be a key regulator of the cell cycle by regulating the transcription of genes required for G1/S and G2/M phase transition , and indeed we find binding in the promoter regions of many cell cycle regulating genes (Figure 1a). However, analysis of the genomic distribution of FOXM1 binding events (Figure 1b) reveals that the majority of binding is found in intronic and intergenic regions, showing a similar pattern to that of nuclear hormone receptors . We carried out functional analysis of the genes within 50 kb of a binding peak and identified 34 over-represented categories with a false discovery rate (FDR) < 0.01 (Figure 1c; Table S1 in Additional file 1). As expected, the top enriched biological processes were cell proliferation, M phase regulation and mitosis. Response to steroid hormone stimulus was another enriched category. Motif enrichment analysis was conducted on the 200 bp region around the peak center of the binding peaks using the motif-based sequence analysis tool AME . Amongst the top significantly enriched motifs present in regions of FOXM1 binding were members of the Forkhead, leucine zipper, GATA and nuclear hormone receptor families (examples shown in Figure 1d; additional motifs identified are shown in Table S2 in Additional file 1).
As the ER binding motif was significantly enriched in the FOXM1 binding data, we overlapped the FOXM1 peaks with a previously published dataset for ERα binding  in asynchronous MCF7 cells. Considering overlapping peaks as those that share at least one base pair, we find a high degree of overlap between FOXM1 and ERα sites, with approximately 80% of FOXM1 binding occurring at ERα binding regions (Figure 2a, b). As FOXA1 has been shown to act as an ERα pioneer factor and is known to be required for ERα chromatin binding, we overlapped FOXM1 binding data with published genome-wide FOXA1 binding data . FOXM1 and FOXA1 showed a high degree of concordance, with approximately 71% of FOXM1 binding events overlapping with FOXA1 binding events (Figure 2b).
We looked at the enriched motifs present in the regions of joint FOXM1/ERα compared to sites where only FOXM1 is located (Figure 2c). There is a high degree of overlap in the motifs identified (Table S2 in Additional file 1); in both cases the forkhead motif (FOXA1) is the most significantly enriched motif and also the leucine zipper motifs for FOS and AP1, which have previously been associated with ER and FOXA1 binding sites  and recently with potential FOXM1 binding regions . The FOXM1-only regions did not contain ESR1 motifs or indeed any motifs for nuclear hormone receptor family members present in the joint binding regions. We identified a number of motifs present only in the regions of FOXM1 binding (Table S3 and Figure S1A in Additional file 1); interestingly, these included the NFY CCAAT motif, which has recently been associated with FOXM1 binding in promoter regions of cell cycle-regulated genes in both U2OS  and HeLa cells . The CTCF motif was enriched in FOXM1 only regions. CTCF has been linked with both ER and FOXA1 binding sites and a recent ChIP-seq study  demonstrated a set of co-localized CTCF/ER or CTCF/FOXA1 binding sites in MCF7 cells; however, the majority were CTCF/FOXA1 joint sites, with the CTCF/ER sites located predominantly in promoter regions.
To confirm that FOXM1 and ERα can simultaneously co-occupy the same chromatin location, re-ChIP experiments followed by quantitative PCR (qPCR) were performed focusing on regions of co-binding and, as controls, regions with only ERα binding (example regions are shown in Figure 2d). ERα ChIP was performed first using asynchronous MCF7 cells followed by FOXM1 ChIP or IgG as a negative control. Results confirmed (Figure 2e) that these transcription factors simultaneously co-occupy the same genomic locations. Using co-immunoprecipitation Figure 2f), we observed a direct interaction between FOXM1 and ERα, suggesting that FOXM1 binds as part of the ER transcriptional complex .
FOXM1 binding at co-bound sites is dependent on ERα
To investigate the relationship between FOXM1 and ERα at sites of co-binding, we looked at the recruitment of FOXM1 following treatment of MCF7 cells with fulvestrant, an estrogen receptor antagonist  that forms an unstable complex with ER, thereby increasing the rate of degradation. Following a 3 hour treatment with fulvestrant, there was a significant reduction in ERα protein but FOXM1 protein and transcript levels were unchanged (Figure 3a; Figure S2 in Additional file 1). We performed ChIP followed by qPCR for regions where FOXM1 and ERα both bind (determined by analysis of ChIP-seq data; Figure 2d) and regions where only FOXM1 binds. This showed (Figure 3b) a significant reduction in FOXM1 binding (P < 0.05) at regions where ERα also binds but not at ERα-negative regions such as the promoter of PLK1. We also looked at the mRNA levels of the target genes of the binding regions (Figure S2 in Additional file 1) and observed a significant reduction in expression for genes that are co-bound by FOXM1 and ERα (P < 0.05). This result suggests that ERα binding is required for FOXM1 recruitment at co-bound regions and confirms the re-ChIP results showing co-occupancy of FOXM1 and ERα at these same genomic loci. In contrast, treatment with the proteasome inhibitor MG132, which is known to reduce FOXM1 expression , reduces FOXM1 binding at both ERα co-bound and also ERα-negative regions (Figure S3 in Additional file 1). We assessed the effect of small interfering RNA (siRNA)-mediated knockdown of FOXM1 on ERα binding by ChIP-qPCR. The level of ERα protein was not significantly altered by FOXM1 knockdown as shown by western blot (Figure 3c) and we did not see any significant differences in ERα binding at known binding sites (Figure 3d), suggesting that FOXM1, unlike FOXA1, is not necessary for ERα to bind to the chromatin. Previous studies have shown that FOXM1 expression is required for ERα-induced proliferation  and indeed we observed that depletion of FOXM1 significantly reduced the expression of a number of ERα-regulated genes (Figure 3e), suggesting a functional role for FOXM1 at these co-bound sites.
Interaction of FOXM1 with ERα and CARM1 co-activator
Since FOXM1 was not required for ERα to bind to DNA, we hypothesized that FOXM1 may be involved in tethering co-factors to the ERα complex. We used an initial proteomic screen (data not shown) to identify FOXM1 interacting proteins and checked for proteins known to be part of the ERα transcriptional complex . Using this method we found a possible interaction with co-activator-associated arginine methyltransferase (CARM1), one of the key co-factors involved in ERα signaling [37, 38]. We demonstrate (Figure 3f) using co-immunoprecipitation that there is a direct interaction between these proteins and propose that this might be important for regulating transcriptional activity at these sites. CARM1 is known to act as a transcriptional regulator by methylation of arginine residues in histones and co-activator proteins such as p160 and CBP/p300 . We therefore investigated whether FOXM1 binding at ERα co-bound sites regulated the methylation activity of CARM1. ChIP-qPCR for H3 arginine 17 modifications in MCF7 cells treated with siRNA for FOXM1 showed a significant decrease following depletion of FOXM1 protein at the sites of co-binding (Figure 3g). This suggests that FOXM1 may assist in recruitment of co-factors, such as CARM1, to ERα binding sites (Figure 3h).
FOXM1 binding is altered in ERα-negative cell lines
To explore the hypothesis that FOXM1 is recruited to DNA by ERα in ER-positive breast cancer cells, we performed ChIP-seq to map FOXM1 binding in the ERα-negative cell line MDA-MB-231 and compared the binding events to ERα-positive MCF7 (Figure 4a). We observed little overlap in FOXM1 binding events between the two cell lines, with only 14% of MCF7 FOXM1 binding occurring at the same genomic location in MDA-MB-231 cells. Approximately 52% of the common binding regions in the two cell lines were found at sites where ERα is not bound in MCF7 cells (example shown in Figure 4b), implying a subset of biologically reproducible ERα-independent FOXM1 binding events. FOXM1 binding in the different cell lines showed different genomic distributions (Figure 4c) with a significant enrichment of binding in promoter regions in the MDA-MB-231 cells when compared to MCF7 cells (approximately 31% compared to approximately 9% in MCF7, compared to approximately 2% for the genome average distribution) and 5' UTR (approximately 13% compared to approximately 2% in MCF7 and approximately 0.4% for the genome average distribution) whilst binding in the intergenic and intronic regions was reduced (approximately 50% compared to approximately 85% in MCF7 and 92% for the genome average distribution). Most of the shared FOXM1 binding events are in promoter regions (Figure S4 in Additional file 1).
Motif analysis of MDA-MB-231 binding peaks showed little overlap with those identified in MCF7 binding peaks in regions of joint binding with ERα; however, comparison with the top enriched motifs in the FOXM1 binding regions in MCF7 cells revealed a similar pattern of cis-elements present (Figure 4d; Table S4 in Additional file 1). Interestingly, although present, the FKH motif was not amongst the top 50 enriched motifs in this cell line. The most significantly enriched motifs were for members of the Ets and leucine zipper families of transcription factors; also present was the NFY CCAAT-binding motifs, which were recently shown to be present in sites of FOXM1 binding in U2OS cells . This suggests that while there is a core set of common FOXM1 binding sites, there are also cell line-specific binding patterns. The lower enrichment of FKH motifs in the MDA-MB-231 cells suggests recruitment by other transcription factors. We found enrichment of the helix-loop-helix transcription factor family motif in the MDA-MB-231 cell line and have observed a high level of concordance with c-MYC binding in this cell line (data not shown).
ChIP-qPCR of FOXM1 binding sites was carried out in additional ERα-positive (namely ZR-751 and T47D) and -negative (namely MDA-MB-453 and MCF10A) cell lines (Figure 4e, f) looking at regions identified in MCF7 cells as FOXM1-only sites or sites of ERα co-binding. We found that FOXM1 binding in all cell lines tested was enriched in the promoter regions of cell cycle-regulating proteins whilst in the ERα-positive cells (ZR-751 and T47D) binding was also detected in regions known to be ERα binding sites in MCF7 (GREB1, RARα and TFF1) but not in the ERα-negative cells (MDA-MB-453 and MCF10A). The pattern of FOXM1 binding was confirmed using two additional antibodies with specificity to different regions of FOXM1 (Figure S5 in Additional file 1) in MCF7 cells.
Inhibition of FOXM1 binding in MCF7 cells modulates binding at specific sites
To study the effect of inhibition of FOXM1 binding on gene expression and cellular phenotype, we treated MCF7 cells with the thiazole antibiotic thiostrepton (Figure 5a). This molecule inhibits FOXM1 activity [40, 41] and we have previously shown that thiostrepton interacts directly with FOXM1 inhibiting DNA binding  at certain genomic loci. In this study we aimed to extend this analysis genome-wide. Thiostrepton is known to reduce FOXM1 expression and as such we treated cells for 4 hours, a timepoint where there is no significant change in FOXM1 protein level  and performed ChIP-seq to map global changes in FOXM1 binding. ChIP-seq was performed in quadruplicate with at least 25 million sequencing reads per replicate (Table S5 in Additional file 1). FOXM1 binding events were identified using MACS  and any FOXM1 binding events present in at least two out of the four replicates were considered. Differential binding analysis (DBA; see Materials and methods) was used on this dataset of 21,494 FOXM1 binding events to identify binding events that were statistically differentially bound in either the control or thiostrepton-treated condition. Using an FDR < 0.05, a total of 1,902 differentially bound peaks were identified from the replicate sets of control/thiostrepton-treated samples (Figure 5b; Additional file 2), with the majority of FOXM1 binding events showing decreased binding after treatment with thiostrepton (1,446, representing approximately 76% of all differential FOXM1 binding events) (Figure 5c), consistent with our previous report that thiostrepton inhibits FOXM1 DNA binding . Examples of differentially bound regions are shown in Figure 5d. Pathway analysis using GREAT (Genomic Region Enrichment Analysis Tool)  on regions associated with differentially bound peaks (Figure 5b) showed that regions with increased FOXM1 binding after thiostrepton treatment were enriched for apoptosis signaling pathways, consistent with the known induction of apoptosis with thiostrepton . Interestingly, in the FOXM1 binding regions with decreased binding after thiostrepton treatment, the top enriched processes were related to the steroid hormone-signaling pathway as well as fibroblast growth factor receptor (FGFR) signaling, a key pathway involved in epithelial to mesenchymal transition and cell proliferation [46, 47]. This supports the qPCR data (Figure 3e) showing decreased expression of ERα-regulated genes following depletion of FOXM1, suggesting a functional role for FOXM1 in ERα activity at co-binding sites. Cis-regulatory element annotation system (CEAS) analysis of the DBA peaks (Figure S6 in Additional file 1) showed that the increased binding was enriched in promoter regions. Interestingly, motif analysis of the DBA peak regions (Figure 5e; Table S6 in Additional file 1) showed that three out of the four enriched motifs were for FKH factors in the down-regulated regions whilst there was no enrichment of FKH motifs in the up-regulated regions.
To test whether the differentially bound FOXM1 regions affect gene expression, we performed gene expression analysis on MCF7 cells after treatment with dimethylsulphoxide (DMSO) or thiostrepton for 6 hours. Using an FDR < 0.01 and log ratio ≥ 0.5, there were 2,322 differentially expressed genes (Figure 6a), with 1,098 down-regulated and 1,224 up-regulated genes. Gene Ontology analysis of the function annotations for the differentially expressed genes showed that out of the top ten categories based on P-value, cell cycle regulation and apoptosis pathways were over represented (Table S7 in Additional file 1 lists the top 10 up/down-regulated genes and Additional file 3 lists all differentially expressed genes). Interestingly the ERα signaling pathway was also highly enriched, again supporting the importance of FOXM1 as a key factor in ERα signaling.
We integrated the differentially expressed genes and differentially bound peaks (identified from the DBA) using either a 10 or 50 kb window. For genes within 10 kb of the nearest peak there was a correlation of 0.48 (95% confidence intervals 0.29 to 0.63) between the thiostrepton regulated FOXM1 binding events and the thiostrepton regulated genes (Figure 6c). Extending the window to 50 kb resulted in a correlation of 0.41 (95% confidence intervals 0.32 to 0.49). Using the 50 kb window we identified 383 peaks associated with genes (Additional file 4). To validate these results we selected 17 down-regulated and 5 up-regulated genes (Table S8 in Additional file 1) and performed independent ChIP experiments with MCF7 cells treated with DMSO/thiostrepton for 4 hours followed by qPCR for the peak region from the DBA (Figure 6d). We found that there was a significant reduction (P < 0.05) in binding in all of the down-regulated regions corresponding with the reduced expression in the microarray dataset. However, for the up-regulated regions we found that the level of enrichment was too low to verify reproducibly.
Inhibition of FOXM1 regulates a gene-signature correlated with prognosis
Gene Set Enrichment Analysis (GSEA)  was used to compare thiostrepton-regulated genes containing a FOXM1 binding peak (FDR < 0.05) within 50 kb of the transcription start site (TSS) (405 unique genes) with curated gene sets in the Molecular Signature Database (MSigDB). We analyzed the up- and down-regulated genes as separate lists. Using the down-regulated gene set (198 identified genes) the top 10 significantly overlapping sets were all related to breast cancer and ERα signaling (Table S9 and Figure S7 in Additional file 1). There was a significant correlation between the thiostrepton down-regulated gene set and genes shown in previous studies to be estrogen-regulated in both ERα-positive cell lines [49, 50] and tumor samples [51, 52] and gene sets related to resistance to endocrine therapy [53–55]. Similar analysis using the genes up-regulated by thiostrepton (111 identified genes) revealed that the ten most significantly overlapped gene sets were from diverse studies (Table S10 in Additional file 1), none of which related to breast cancer or ERα signaling.
As FOXM1 itself is a poor prognosis marker in ERα breast cancer (Figure S8 in Additional file 1) we tested whether the genes regulated by thiostrepton and with a FOXM1 binding site within 50 kb of the TSS were correlated with clinical outcome in breast cancer. In this case we again considered the up- and down-regulated genes separately and tested whether there was any correlation with patient survival using a publically available data set of 286 breast cancer patients with associated survival data . We restricted analysis to ERα-positive patients (209) and used this set to test for enrichment of genes correlated with prognosis in the thiostrepton-regulated gene sets. In the case of the up-regulated genes there was no evidence of overrepresentation (P = 0.217). However, for the down-regulated genes there was a significant correlation (P < 10-6) with 38 genes identified (Table S11 in Additional file 1) whose expression is related to prognosis (high expression correlated to reduced time to relapse; Figure 7a, b; Figure S9 in Additional file 1). Survival analysis using this 38 gene-signature showed a significant correlation (P < 0.0034) with patient survival (Figure 7b), with high expression correlating with reduced time to relapse; these results were also confirmed using another independent dataset (P = 0.0094)  (Figure 7c; Figure S10 in Additional file 1). We looked for known functional interactions within the genes identified using STRING (Search Tool for the Retrieval of Interacting Genes/Proteins)  and found high confidence interactions linking 20 proteins relating to mitotic spindle formation, chromosome segregation and cell cycle regulation (Figure S11 in Additional file 1).
This confirms the therapeutic potential of targeting FOXM1, as thiostrepton treatment reduces the expression of genes whose high expression correlates with poor patient outcome in breast cancer.
The forkhead factor FOXM1 was first linked with cancer in 2002  and has subsequently been shown to be involved in cancer initiation , progression  and prognosis  in diverse tissue types. FOXM1 is a known regulator of cell proliferation and mitotic spindle assembly genes but has also been linked with regulating epithelial to mesenchymal transition, angiogenesis, metastasis and more recently with induction of DNA methylation .
In this study we focus on the role of FOXM1 in breast cancer with the aim of mapping global FOXM1 binding sites in both ER-positive and -negative cancer cell lines and relating this to biological function. We show that FOXM1 binding occurs in the promoter regions of many cell cycle-related genes and also genes known to regulate mitotic spindle assembly; BUB1, ASPM, MAD2L1 and members of the kinesin family. We also confirmed binding in target genes such as CAV1  and MMP family members; MMP2 and 8, VEGF and FGF linking with the role of FOXM1 in metastasis.
The high degree of overlap between FOXM1 binding and ERα binding in MCF7 is of significant interest as previous studies have shown co-regulation of these factors [21, 22] and importantly that FOXM1 expression is essential for ERα-induced proliferation. We show that FOXM1 and ERα bind to the chromatin at the same place at the same time; however, the presence of FOXM1 at ERα co-bound chromatin sites is not required for ERα recruitment but is necessary for the loading of other cofactors to the ER transcriptional complex, such as CARM1. CARM1 is known to activate transcription by methylating histone H3 arginine residues  and constitutes a key part of nuclear hormone receptor transcriptional complexes acting in conjunction with other coactivators, such as CBP/p300, SRC-2 and GRIP1 [33, 64]. FOXM1 has previously been shown to interact directly with CARM1 binding factors such as p300/CBP  but, to our knowledge, this is the first time that a direct interaction has been shown between FOXM1 and CARM1.
The differences in FOXM1 binding with ERα status in breast cancer cells contribute another dimension to the already complex mechanisms known to regulate FOXM1 transactivation. FOXM1 activity is regulated by multiple protein-protein interactions; during the cell cycle the temporal control required for precise G2/M transition is regulated both by phosphorylation steps involving different Cyclin/Cdks [66–68] and by inhibition by B55a, a subunit of protein phosphatase 2A . Transcriptional activity is modulated by recruitment of coactivators and interaction with other transcription factors . The data from our study suggest another mode of interaction of FOXM1, in this case co-binding with ERα in ERα-positive breast cancer cells.
A similar role for FOXM1 in transcriptional regulation by recruitment of co-factors at specific genomic locations was proposed in a recent study by Carr et al. . They showed that FOXM1 represses GATA-3 expression in the mammary gland by recruiting the methyltransferase DNMT3b to binding sites within the GATA-3 promoter, thereby leading to methylation-induced gene silencing.
Inhibition of FOXM1 binding in MCF7 cells using thiostrepton confirmed the importance of this transcription factor in the regulation of ERα signaling pathways. Analysis of the regions where binding was significantly reduced showed that these were cis-regulatory regions involved in ERα pathways and indeed the gene expression changes confirmed that ERα signaling pathways were affected. Somewhat surprisingly, FOXM1 binding was increased at a number of genomic sites following thiostrepton treatment; however, these sites did not contain consensus FKH binding motifs and may represent a redistribution of FOXM1 in a similar manner to that described for the androgen receptor following FOXA1 depletion . In the absence of FOXA1, androgen receptor binding was delocalized around the genome, generating, in addition to FOXA1 dependent-sites, both a subset of FOXA1-inhibited sites and a subset of FOXA1-independent binding sites where androgen receptor recruitment may occur via the presence of transcription factors other than FOXA1 with pioneer activity. Furthermore, a recent publication showed that FOXM1 has an atypical DNA binding mechanism , with binding at some genomic sites occurring via recruitment and binding to the MMB transcriptional activator complex. This mechanism not only explains how binding is increased at certain sites after inhibition of direct DNA binding but also fits with the binding pattern we observed in MDA-MB-231 cells compared to MCF7 cells, where the FKH motif was poorly enriched compared to other motifs such as the ETS family members. In fact we cannot rule out that some of the enriched FKH consensus sites in MCF7 cells actually represent FOXA1 binding motifs with FOXM1 tethered via other transcriptional co-factors. Thus, FOXM1, as a key regulator of the cell cycle, may show cell line-dependent patterns of DNA binding due to recruitment by different transcriptional complexes responsible for driving the proliferation in a particular cellular context.
As FOXM1 expression has previously been correlated with prognosis in breast cancer we used our dataset to identify a novel FOXM1-regulated gene set significantly correlated with ER-positive breast cancer prognosis and drug resistance, thus confirming previous studies showing that FOXM1 over-expression is associated with resistance to drug treatment [19, 21]. This gene-signature of 38 FOXM1-regulated genes down-regulated by thiostrepton treatment was predictive of prognosis in ER-positive breast cancer patient datasets. Within this set are many well described FOXM1 target genes, such as those involved in mitotic spindle formation; however, there are a number of interesting novel FOXM1 target genes, such as ABCC5, a transporter protein associated with multi-drug resistance that may provide useful insights for future studies into the role of FOXM1 in breast cancer.
We have demonstrated that FOXM1 shows distinct patterns of binding depending on ERα status in breast cancer cells, but within an ER-positive context FOXM1 plays an important role in ERα signaling pathways. Specifically, FOXM1 regulates a gene signature that correlates with poor prognosis in breast cancer patients, supporting the therapeutic potential in targeting FOXM1 in ER-positive breast cancer.
Materials and methods
Human MCF7, MDA-MB-231, ZR-751, T47D, MDA-MB-453 and MCF10A cell lines were obtained from the ECACC (European Collection of Animal Cell Cultures) and grown in DMEM or RPMI (ZR-751, T47D) supplemented with 10% FBS and the MCF10A in MEGB (Lonza Biologics, Slough, Berkshire, UK). Thiostrepton (Sigma-Aldrich Company Ltd, Gillingham, Dorset, UK) and MG132 (Sigma) were made up as stock solutions of 10 mM in DMSO and used at a final concentration of 10 μM (thiostrepton) and 3 μM (MG132) and fulvestrant (ICI 182780) was added at a final concentration of 10 nM.
Western blots were performed using antibodies for anti-FOXM1 (sc-502) and anti-ER (sc-543) from Santa Cruz Biotechnology (Dallas, Texas, USA) and anti-β-actin (ab6276) from Abcam (Cambridge, Cambridgeshire, UK). Further details are provided in the Supplementary materials and methods in Additional file 1.
Quantitative real-time PCR analysis
MCF7 cells were treated with thiostrepton/MG132/fulvestrant and RNA was collected after the indicated timepoints and qPCR was performed using Power Sybr mix (ABI, Warrington, Cheshire, UK) on a CFX96 Real-time thermal cycler (Bio-Rad, Hemel Hempstead, Hertfordshire, UK). Further details and primer sequences are given in Additional file 1.
Small interfering siRNA
MCF7 cells were transfected with Dharmacon siGenome pools for FOXM1 or siCONTROL at 15 nM using lipofectamine 2000. Knockdown of FOXM1 protein levels was confirmed by western blotting after 48 h incubation.
MCF7 cells were treated with either thiostrepton or DMSO for 6 h. Each treatment was carried out in replicates of six. Gene expression analysis was carried out on Illumina Human HT12 version 4 bead arrays. Details of analysis are given in the Supplementary materials and methods in Additional file 1. Data are available through the NCBI's Gene Expression Omnibus  using GEO Series accession number GSE40767.
ChIP experiments were performed as previously described . Antibodies used were anti-FOXM1 (Santa Cruz sc-502, Genetex GTX1000276, Genetex GTX102170 [GeneTex, Irvine, California, USA]) and anti-ER (Santa Cruz sc-543), anti-histone H3 (Abcam ab1791) and anti-histone H3 (asym-dimethyl Arg17; Novus NB21-1132 (Novus Biologicals, Cambridge, Cambridgeshire, UK)). Experimental details and primer sequences are available in Additional file 1.
To detect FOXM1 binding sites affected by thiostrepton versus DMSO, ChIP-Seq experiments were performed in four biological replicates in MCF7 cells and in two biological replicates in MDA-MB-231 (Table S5 in Additional file 1). ChIP DNA was processed for Illumina sequencing as previously described . Details of analysis are in Additional file 3. Data are available through the NCBI's Gene Expression Omnibus  using GEO Series accession number GSE40767.
Analysis of differential binding
To identify regions of differential FOXM1 binding between the control (DMSO) and thiostrepton-treated samples, a general linear model was fitted to each putative binding site to test for the difference in read count between treatments (see Supplementary materials and methods in Additional file 1 for detection and definition of binding sites). Model fitting and testing was performed using the Bioconductor library edgeR  using the function estimateGLMTagwiseDisp for estimating the dispersion parameter of the negative binomial distribution and glmFit and glmLRT for fitting and testing the difference of treatment of each binding site . The heatmaps in Figures 2 and 5 were prepared with Java Treeview .
Motif analysis and genomic distribution of binding events
The CEAS  function in Cistrome  was used to functionally annotate binding sites. Known transcription factor motifs significantly enriched in the binding sites were identified with AME . The -log p-value of significance was used to scale the word clouds using the R package wordcloud .
Gene Ontology pathway analysis
Gene Otology pathway enrichment was performed using GeneGo metacore (MetaCore from Thomson Reuters (New York, USA)) and visualized with REViGO (Reduce and Visualize Gene Ontology) .
Re-ChIP was performed as described  using the same antibodies as for ChIP and normal rabbit IgG (2729) from Cell Signaling Technology (Danvers, Massachusetts, USA).
Experiments were performed using the nuclear co-immunoprecipitation kit from Active motif (Carlsbad, California, USA) with pull-down using FOXM1 (Santa Cruz sc-502) following the manufacturer's protocol with immunoprecipitation carried out using either the low or high buffers provided supplemented with 1X protease inhibitor cocktail and 1 mM dithiothreitol (DTT). Detection was performed by western blotting using ERα (Novacastra, Milton Keynes, Buckingham, UK) and CARM1 (Santa Cruz sc-5421) with LiCor IRDye secondary antibodies; 800LT goat anti-mouse, 680LT goat anti-rabbit and 680LT donkey anti-goat. Further details are in Supplementary materials and methods in Additional file 1.
The clinical and gene expression dataset GSE2034  was analyzed to investigate whether genes regulated by thiostrepton are involved in poor cancer prognosis; only the 209 ER-positive samples were considered. Details of the analysis are in Supplementary materials and methods in Additional file 1). Gene Set Enrichment Analysis (GSEA)  was performed using the web-based tool  to compare gene lists with the Molecular Signatures Database (MSigDB).
All statistical analyses not described above were performed with GraphPad prism software or R . The tests for difference between means were performed using the two-tailed Student's t-test. If not otherwise stated, P-value < 0.05 was considered statistically significant. Error bars represent standard deviations.
cis-regulatory element annotation system
differential binding analysis
false discovery rate
Genomic Region Enrichment Analysis Tool
quantitative polymerase chain reaction
Reduce and Visualize Gene Ontology
small interfering RNA
transcription start site
Korver W, Roose J, Clevers H: The winged-helix transcription factor Trident is expressed in cycling cells. Nucleic Acids Res. 1997, 25: 1715-1719. 10.1093/nar/25.9.1715.
Costa RH: FoxM1 dances with mitosis. Nat Cell Biol. 2005, 7: 108-110. 10.1038/ncb0205-108.
Wang IC, Chen YJ, Hughes D, Petrovic V, Major ML, Park HJ, Tan Y, Ackerson T, Costa RH: Forkhead box M1 regulates the transcriptional network of genes essential for mitotic progression and genes encoding the SCF (Skp2-Cks1) ubiquitin ligase. Mol Cell Biol. 2005, 25: 10875-10894. 10.1128/MCB.25.24.10875-10894.2005.
Wonsey DR, Follettie MT: Loss of the forkhead transcription factor FoxM1 causes centrosome amplification and mitotic catastrophe. Cancer Res. 2005, 65: 5181-5189. 10.1158/0008-5472.CAN-04-4059.
Wang Z, Banerjee S, Kong D, Li Y, Sarkar FH: Down-regulation of Forkhead Box M1 transcription factor leads to the inhibition of invasion and angiogenesis of pancreatic cancer cells. Cancer Res. 2007, 67: 8293-8300. 10.1158/0008-5472.CAN-07-1265.
Raychaudhuri P, Park HJ: FoxM1: a master regulator of tumor metastasis. Cancer Res. 2011, 71: 4329-4333. 10.1158/0008-5472.CAN-11-0640.
Park HJ, Carr JR, Wang Z, Nogueira V, Hay N, Tyner AL, Lau LF, Costa RH, Raychaudhuri P: FoxM1, a critical regulator of oxidative stress during oncogenesis. EMBO J. 2009, 28: 2908-2918. 10.1038/emboj.2009.239.
Alvarez-Fernandez M, Halim VA, Krenning L, Aprelia M, Mohammed S, Heck AJ, Medema RH: Recovery from a DNA-damage-induced G2 arrest requires Cdk-dependent activation of FoxM1. EMBO Rep. 2010, 11: 452-458. 10.1038/embor.2010.46.
Pilarsky C, Wenzig M, Specht T, Saeger HD, Grutzmann R: Identification and validation of commonly overexpressed genes in solid tumors by comparison of microarray data. Neoplasia. 2004, 6: 744-750. 10.1593/neo.04277.
Chu XY, Zhu ZM, Chen LB, Wang JH, Su QS, Yang JR, Lin Y, Xue LJ, Liu XB, Mo XB: FOXM1 expression correlates with tumor invasion and a poor prognosis of colorectal cancer. Acta Histochem. 2012, 114: 755-762. 10.1016/j.acthis.2012.01.002.
Priller M, Poschl J, Abrao L, von Bueren AO, Cho YJ, Rutkowski S, Kretzschmar HA, Schuller U: Expression of FoxM1 is required for the proliferation of medulloblastoma cells and indicates worse survival of patients. Clin Cancer Res. 2011, 17: 6791-6801. 10.1158/1078-0432.CCR-11-1214.
Teh MT, Gemenetzidis E, Chaplin T, Young BD, Philpott MP: Upregulation of FOXM1 induces genomic instability in human epidermal keratinocytes. Mol Cancer. 2010, 9: 45-10.1186/1476-4598-9-45.
Gemenetzidis E, Elena-Costea D, Parkinson EK, Waseem A, Wan H, Teh MT: Induction of human epithelial stem/progenitor expansion by FOXM1. Cancer Res. 2010, 70: 9515-9526. 10.1158/0008-5472.CAN-10-2173.
Park YY, Jung SY, Jennings N, Rodriguez-Aguayo C, Peng G, Lee SR, Kim SB, Kim K, Leem SH, Lin SY, Lopez-Berestein G, Sood A, Lee JS: FOXM1 mediates Dox resistance in breast cancer by enhancing DNA repair. Carcinogenesis. 2012, 33: 1843-1853. 10.1093/carcin/bgs167.
Francis RE, Myatt SS, Krol J, Hartman J, Peck B, McGovern UB, Wang J, Guest SK, Filipovic A, Gojis O, Palmieri C, Peston D, Shousha S, Yu Q, Sicinski P, Coombes RC, Lam EW: FoxM1 is a downstream target and marker of HER2 overexpression in breast cancer. Int J Oncol. 2009, 35: 57-68.
Bektas N, Haaf A, Veeck J, Wild PJ, Luscher-Firzlaff J, Hartmann A, Knuchel R, Dahl E: Tight correlation between expression of the Forkhead transcription factor FOXM1 and HER2 in human breast cancer. BMC Cancer. 2008, 8: 42-10.1186/1471-2407-8-42.
Bergamaschi A, Christensen BL, Katzenellenbogen BS: Reversal of endocrine resistance in breast cancer: interrelationships among 14-3-3zeta, FOXM1, and a gene signature associated with mitosis. Breast Cancer Res. 2011, 13: R70-10.1186/bcr2913.
Abba MC, Lacunza E, Butti M, Aldaz CM: Breast cancer biomarker discovery in the functional genomic age: a systematic review of 42 gene expression signatures. Biomarker Insights. 2010, 5: 103-118.
Kwok JM, Peck B, Monteiro LJ, Schwenen HD, Millour J, Coombes RC, Myatt SS, Lam EW: FOXM1 confers acquired cisplatin resistance in breast cancer cells. Mol Cancer Res. 2010, 8: 24-34. 10.1158/1541-7786.MCR-09-0432.
Carr JR, Park HJ, Wang Z, Kiefer MM, Raychaudhuri P: FoxM1 mediates resistance to herceptin and paclitaxel. Cancer Res. 2010, 70: 5054-5063. 10.1158/0008-5472.CAN-10-0545.
Millour J, Constantinidou D, Stavropoulou AV, Wilson MS, Myatt SS, Kwok JM, Sivanandan K, Coombes RC, Medema RH, Hartman J, Lykkesfeldt AE, Lam EW: FOXM1 is a transcriptional target of ERalpha and has a critical role in breast cancer endocrine sensitivity and resistance. Oncogene. 2010, 29: 2983-2995. 10.1038/onc.2010.47.
Madureira PA, Varshochi R, Constantinidou D, Francis RE, Coombes RC, Yao KM, Lam EW: The Forkhead box M1 protein regulates the transcription of the estrogen receptor alpha in breast cancer cells. J Biol Chem. 2006, 281: 25167-25176. 10.1074/jbc.M603906200.
Weisz A, Basile W, Scafoglio C, Altucci L, Bresciani F, Facchiano A, Sismondi P, Cicatiello L, De Bortoli M: Molecular identification of ERalpha-positive breast cancer cells by the expression profile of an intrinsic set of estrogen regulated genes. J Cell Physiol. 2004, 200: 440-450. 10.1002/jcp.20039.
Yau C, Wang Y, Zhang Y, Foekens JA, Benz CC: Young age, increased tumor proliferation and FOXM1 expression predict early metastatic relapse only for endocrine-dependent breast cancers. Breast Cancer Res Treat. 2011, 126: 803-810. 10.1007/s10549-011-1345-1.
Hurtado A, Holmes KA, Ross-Innes CS, Schmidt D, Carroll JS: FOXA1 is a key determinant of estrogen receptor function and endocrine response. Nat Genet. 2011, 43: 27-33. 10.1038/ng.730.
Tang Q, Chen Y, Meyer C, Geistlinger T, Lupien M, Wang Q, Liu T, Zhang Y, Brown M, Liu XS: A comprehensive view of nuclear receptor cancer cistromes. Cancer Res. 2011, 71: 6940-6947. 10.1158/0008-5472.CAN-11-2091.
McLeay RC, Bailey TL: Motif Enrichment Analysis: a unified framework and an evaluation on ChIP data. BMC Bioinformatics. 2010, 11: 165-10.1186/1471-2105-11-165.
Ross-Innes CS, Stark R, Teschendorff AE, Holmes KA, Ali HR, Dunning MJ, Brown GD, Gojis O, Ellis IO, Green AR, Ali S, Chin SF, Palmieri C, Caldas C, Carroll JS: Differential oestrogen receptor binding is associated with clinical outcome in breast cancer. Nature. 2012, 481: 389-393.
Joseph R, Orlov YL, Huss M, Sun W, Kong SL, Ukil L, Pan YF, Li G, Lim M, Thomsen JS, Ruan Y, Clarke ND, Prabhakar S, Cheung E, Liu ET: Integrative model of genomic factors for determining binding site selection by estrogen receptor-alpha. Mol Systems Biol. 2010, 6: 456-
Sadasivam S, Duan S, DeCaprio JA: The MuvB complex sequentially recruits B-Myb and FoxM1 to promote mitotic gene expression. Genes Dev. 2012, 26: 474-489. 10.1101/gad.181933.111.
Chen X, Muller GA, Quaas M, Fischer M, Han N, Stutchbury B, Sharrocks AD, Engeland K: The forkhead transcription factor FOXM1 controls cell cycle-dependent gene expression through an atypical chromatin binding mechanism. Mol Cell Biol. 2012, 33: 227-236.
Ross-Innes CS, Brown GD, Carroll JS: A co-ordinated interaction between CTCF and ER in breast cancer cells. BMC Genomics. 2011, 12: 593-10.1186/1471-2164-12-593.
Klinge CM, Jernigan SC, Mattingly KA, Risinger KE, Zhang J: Estrogen response element-dependent regulation of transcriptional activation of estrogen receptors alpha and beta by coactivators and corepressors. J Mol Endocrinol. 2004, 33: 387-410. 10.1677/jme.1.01541.
Osborne CK, Wakeling A, Nicholson RI: Fulvestrant: an oestrogen receptor antagonist with a novel mechanism of action. Br J Cancer. 2004, 90 (Suppl 1): S2-6.
Bhat UG, Halasi M, Gartel AL: FoxM1 is a general target for proteasome inhibitors. PloS One. 2009, 4: e6593-10.1371/journal.pone.0006593.
Metivier R, Penot G, Hubner MR, Reid G, Brand H, Kos M, Gannon F: Estrogen receptor-alpha directs ordered, cyclical, and combinatorial recruitment of cofactors on a natural target promoter. Cell. 2003, 115: 751-763. 10.1016/S0092-8674(03)00934-6.
Al-Dhaheri M, Wu J, Skliris GP, Li J, Higashimato K, Wang Y, White KP, Lambert P, Zhu Y, Murphy L, Xu W: CARM1 is an important determinant of ERalpha-dependent breast cancer cell differentiation and proliferation in breast cancer cells. Cancer Res. 2011, 71: 2118-2128. 10.1158/0008-5472.CAN-10-2426.
Naeem H, Cheng D, Zhao Q, Underhill C, Tini M, Bedford MT, Torchia J: The activity and stability of the transcriptional coactivator p/CIP/SRC-3 are regulated by CARM1-dependent methylation. Mol Cell Biol. 2007, 27: 120-134. 10.1128/MCB.00815-06.
Ceschin DG, Walia M, Wenk SS, Duboe C, Gaudon C, Xiao Y, Fauquier L, Sankar M, Vandel L, Gronemeyer H: Methylation specifies distinct estrogen-induced binding site repertoires of CBP to chromatin. Genes Dev. 2011, 25: 1132-1146. 10.1101/gad.619211.
Radhakrishnan SK, Bhat UG, Hughes DE, Wang IC, Costa RH, Gartel AL: Identification of a chemical inhibitor of the oncogenic transcription factor forkhead box M1. Cancer Res. 2006, 66: 9731-9735. 10.1158/0008-5472.CAN-06-1576.
Kwok JM, Myatt SS, Marson CM, Coombes RC, Constantinidou D, Lam EW: Thiostrepton selectively targets breast cancer cells through inhibition of forkhead box M1 expression. Mol Cancer Ther. 2008, 7: 2022-2032. 10.1158/1535-7163.MCT-08-0188.
Hegde NS, Sanders DA, Rodriguez R, Balasubramanian S: The transcription factor FOXM1 is a cellular target of the natural product thiostrepton. Nat Chem. 2011, 3: 725-731. 10.1038/nchem.1114.
Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, Nusbaum C, Myers RM, Brown M, Li W, Liu XS: Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008, 9: R137-10.1186/gb-2008-9-9-r137.
McLean CY, Bristor D, Hiller M, Clarke SL, Schaar BT, Lowe CB, Wenger AM, Bejerano G: GREAT improves functional interpretation of cis-regulatory regions. Nat Biotechnol. 2010, 28: 495-501. 10.1038/nbt.1630.
Bhat UG, Halasi M, Gartel AL: Thiazole antibiotics target FoxM1 and induce apoptosis in human cancer cells. PloS One. 2009, 4: e5592-10.1371/journal.pone.0005592.
Tomlinson DC, Baxter EW, Loadman PM, Hull MA, Knowles MA: FGFR1-induced epithelial to mesenchymal transition through MAPK/PLCgamma/COX-2-mediated mechanisms. PloS One. 2012, 7: e38972-10.1371/journal.pone.0038972.
Koziczak M, Holbro T, Hynes NE: Blocking of FGFR signaling inhibits breast cancer cell proliferation through downregulation of D-type cyclins. Oncogene. 2004, 23: 3501-3508. 10.1038/sj.onc.1207331.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005, 102: 15545-15550. 10.1073/pnas.0506580102.
Gozgit JM, Pentecost BT, Marconi SA, Ricketts-Loriaux RS, Otis CN, Arcaro KF: PLD1 is overexpressed in an ER-negative MCF-7 cell line variant and a subset of phospho-Akt-negative breast carcinomas. Br J Cancer. 2007, 97: 809-817. 10.1038/sj.bjc.6603926.
Frasor J, Stossi F, Danes JM, Komm B, Lyttle CR, Katzenellenbogen BS: Selective estrogen receptor modulators: discrimination of agonistic versus antagonistic activities by gene expression profiling in breast cancer cells. Cancer Res. 2004, 64: 1522-1533. 10.1158/0008-5472.CAN-03-3326.
Doane AS, Danso M, Lal P, Donaton M, Zhang L, Hudis C, Gerald WL: An estrogen receptor-negative breast cancer subset characterized by a hormonally regulated transcriptional program and response to androgen. Oncogene. 2006, 25: 3994-4008. 10.1038/sj.onc.1209415.
Smid M, Wang Y, Zhang Y, Sieuwerts AM, Yu J, Klijn JG, Foekens JA, Martens JW: Subtypes of breast cancer show preferential site of relapse. Cancer Res. 2008, 68: 3108-3114. 10.1158/0008-5472.CAN-07-5644.
Massarweh S, Osborne CK, Creighton CJ, Qin L, Tsimelzon A, Huang S, Weiss H, Rimawi M, Schiff R: Tamoxifen resistance in breast tumors is driven by growth factor receptor signaling with repression of classic estrogen receptor genomic function. Cancer Res. 2008, 68: 826-833. 10.1158/0008-5472.CAN-07-2707.
Riggins RB, Lan JP, Zhu Y, Klimach U, Zwart A, Cavalli LR, Haddad BR, Chen L, Gong T, Xuan J, Ethier SP, Clarke R: ERRgamma mediates tamoxifen resistance in novel models of invasive lobular breast cancer. Cancer Res. 2008, 68: 8908-8917. 10.1158/0008-5472.CAN-08-2669.
Creighton CJ, Massarweh S, Huang S, Tsimelzon A, Hilsenbeck SG, Osborne CK, Shou J, Malorni L, Schiff R: Development of resistance to targeted therapies transforms the clinically associated molecular profile subtype of breast tumor xenografts. Cancer Res. 2008, 68: 7493-7501. 10.1158/0008-5472.CAN-08-1404.
Wang Y, Klijn JG, Zhang Y, Sieuwerts AM, Look MP, Yang F, Talantov D, Timmermans M, Meijer-van Gelder ME, Yu J, Jatkoe T, Berns EM, Atkins D, Foekens JA: Gene-expression profiles to predict distant metastasis of lymph-node-negative primary breast cancer. Lancet. 2005, 365: 671-679.
Loi S, Haibe-Kains B, Desmedt C, Lallemand F, Tutt AM, Gillet C, Ellis P, Harris A, Bergh J, Foekens JA, Klijn JG, Larsimont D, Buyse M, Bontempi G, Delorenzi M, Piccart MJ, Sotiriou C: Definition of clinically distinct molecular subtypes in estrogen receptor-positive breast carcinomas through genomic grade. J Clin Oncol. 2007, 25: 1239-1246. 10.1200/JCO.2006.07.1522.
Snel B, Lehmann G, Bork P, Huynen MA: STRING: a web-server to retrieve and display the repeatedly occurring neighbourhood of a gene. Nucleic Acids Res. 2000, 28: 3442-3444. 10.1093/nar/28.18.3442.
Teh MT, Wong ST, Neill GW, Ghali LR, Philpott MP, Quinn AG: FOXM1 is a downstream target of Gli1 in basal cell carcinomas. Cancer Res. 2002, 62: 4773-4780.
Kalinichenko VV, Major ML, Wang X, Petrovic V, Kuechle J, Yoder HM, Dennewitz MB, Shin B, Datta A, Raychaudhuri P, Costa RH: Foxm1b transcription factor is essential for development of hepatocellular carcinomas and is negatively regulated by the p19ARF tumor suppressor. Genes Dev. 2004, 18: 830-850. 10.1101/gad.1200704.
Teh MT, Gemenetzidis E, Patel D, Tariq R, Nadir A, Bahta AW, Waseem A, Hutchison IL: FOXM1 induces a global methylation signature that mimics the cancer epigenome in head and neck squamous cell carcinoma. PloS One. 2012, 7: e34329-10.1371/journal.pone.0034329.
Huang C, Qiu Z, Wang L, Peng Z, Jia Z, Logsdon CD, Le X, Wei D, Huang S, Xie K: A novel FoxM1-caveolin signaling pathway promotes pancreatic cancer invasion and metastasis. Cancer Res. 2012, 72: 655-665. 10.1158/0008-5472.CAN-11-3102.
Bauer UM, Daujat S, Nielsen SJ, Nightingale K, Kouzarides T: Methylation at arginine 17 of histone H3 is linked to gene activation. EMBO Rep. 2002, 3: 39-44. 10.1093/embo-reports/kvf013.
Frietze S, Lupien M, Silver PA, Brown M: CARM1 regulates estrogen-stimulated breast cancer growth through up-regulation of E2F1. Cancer Res. 2008, 68: 301-306. 10.1158/0008-5472.CAN-07-1983.
Major ML, Lepe R, Costa RH: Forkhead box M1B transcriptional activity requires binding of Cdk-cyclin complexes for phosphorylation-dependent recruitment of p300/CBP coactivators. Mol Cell Biol. 2004, 24: 2649-2661. 10.1128/MCB.24.7.2649-2661.2004.
Wierstra I, Alves J: Transcription factor FOXM1c is repressed by RB and activated by cyclin D1/Cdk4. Biol Chem. 2006, 387: 949-962.
Wierstra I, Alves J: FOXM1c is activated by cyclin E/Cdk2, cyclin A/Cdk2, and cyclin A/Cdk1, but repressed by GSK-3alpha. Biochem Biophys Res Commun. 2006, 348: 99-108. 10.1016/j.bbrc.2006.07.008.
Anders L, Ke N, Hydbring P, Choi YJ, Widlund HR, Chick JM, Zhai H, Vidal M, Gygi SP, Braun P, Sicinski P: A systematic screen for CDK4/6 substrates links FOXM1 phosphorylation to senescence suppression in cancer cells. Cancer Cell. 2011, 20: 620-634. 10.1016/j.ccr.2011.10.001.
Alvarez-Fernandez M, Halim VA, Aprelia M, Laoukili J, Mohammed S, Medema RH: Protein phosphatase 2A (B55alpha) prevents premature activation of forkhead transcription factor FoxM1 by antagonizing cyclin A/cyclin-dependent kinase-mediated phosphorylation. J Biol Chem. 2011, 286: 33029-33036. 10.1074/jbc.M111.253724.
Petrovic V, Costa RH, Lau LF, Raychaudhuri P, Tyner AL: Negative regulation of the oncogenic transcription factor FoxM1 by thiazolidinediones and mithramycin. Cancer Biol Ther. 2010, 9: 1008-1016. 10.4161/cbt.9.12.11710.
Carr JR, Kiefer MM, Park HJ, Li J, Wang Z, Fontasarosa J, DeWaal D, Kopanja D, Benevolenskaya EV, Guzman G, Raychaudhuri P: FOXM1 regulates mammary luminal cell fate. Cell Reports. 2012, 1-15.
Sahu B, Laakso M, Ovaska K, Mirtti T, Lundin J, Rannikko A, Sankila A, Turunen JP, Lundin M, Konsti J, Vesterinen T, Nordling S, Kallioniemi O, Hautaniemi S, Janne OA: Dual role of FoxA1 in androgen receptor binding to chromatin, androgen signalling and prostate cancer. EMBO J. 2011, 30: 3962-3976. 10.1038/emboj.2011.328.
Gene Expression Omnibus. [http://www.ncbi.nlm.nih.gov/geo/]
Schmidt D, Wilson MD, Spyrou C, Brown GD, Hadfield J, Odom DT: ChIP-seq: using high-throughput sequencing to discover protein-DNA interactions. Methods. 2009, 48: 240-248. 10.1016/j.ymeth.2009.03.001.
Robinson MD, McCarthy DJ, Smyth GK: edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010, 26: 139-140. 10.1093/bioinformatics/btp616.
McCarthy DJ, Chen Y, Smyth GK: Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation. Nucleic Acids Res. 2012, 40: 4288-4297. 10.1093/nar/gks042.
Saldanha AJ: Java Treeview--extensible visualization of microarray data. Bioinformatics. 2004, 20: 3246-3248. 10.1093/bioinformatics/bth349.
Shin H, Liu T, Manrai AK, Liu XS: CEAS: cis-regulatory element annotation system. Bioinformatics. 2009, 25: 2605-2606. 10.1093/bioinformatics/btp479.
Supek F, Bosnjak M, Skunca N, Smuc T: REVIGO summarizes and visualizes long lists of gene ontology terms. PLoS One. 2011, 6: e21800-10.1371/journal.pone.0021800.
Bararia D, Trivedi AK, Zada AA, Greif PA, Mulaw MA, Christopeit M, Hiddemann W, Bohlander SK, Behre G: Proteomic identification of the MYST domain histone acetyltransferase TIP60 (HTATIP) as a co-activator of the myeloid transcription factor C/EBPalpha. Leukemia. 2008, 22: 800-807. 10.1038/sj.leu.2405101.
Gene Set Enrichment Analysis. [http://www.broadinstitute.org/gsea/]
The R Project for Statistical Computing. [http://www.R-project.org]
We would like to acknowledge the Genomics and bioinformatics core at the CRUK Research institute for the Illumina sequencing and microarray analysis and in particular Dr Suraj Menon and Dr Sarah Vowler for statistical advice. We would also like to acknowledge the assistance and advice of Gordon Brown for ChIP-seq analysis, Rosalind Launchbury for technical advice and Dr Nagaratna Hedge for helpful suggestions. We acknowledge the support of The University of Cambridge and Cancer Research UK. The Balasubramanian Laboratory is supported by core funding from Cancer Research UK.
The authors declare that they have no competing interests.
DS, CSR, SB and JSC designed all experiments. DS and CSR performed the experimental work. DB performed analysis and processing of ChIP-seq data and statistical analysis. All authors contributed to the interpretation of the data. The manuscript was written by DS, CSR and DB with help from other authors. All authors read and approved the final manuscript.
Deborah A Sanders, Caryn S Ross-Innes contributed equally to this work.
Electronic supplementary material
Additional file 2: FOXM1 binding peaks in MCF7 cells. Excel spreadsheet containing the location of FOXM1 binding peaks identified using MACS in replicate samples of MCF7 cells treated with DMSO or thiostrepton. Differential bound peaks (FDR < 0.05) were identified using edgeR. (XLSX 2 MB)
Additional file 3: Differentially expressed gene list from microarray analysis. Excel spreadsheet with details of differentially expressed gene (FDR < 0.01) from microarray analysis for DMSO versus thiostrepton-treated MCF7 cells. (XLSX 1 MB)
Additional file 4: FOXM1 binding peaks from DBA versus expression changes from microarray. Excel spreadsheet containing LogFC values for differential expression of genes with TSSs located within ± 50 kb of differentially bound FOXM1 identified from DBA in MCF7 cells treated with DMSO or thiostrepton. (XLSX 6 MB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Sanders, D.A., Ross-Innes, C.S., Beraldi, D. et al. Genome-wide mapping of FOXM1 binding reveals co-binding with estrogen receptor alpha in breast cancer cells. Genome Biol 14, R6 (2013). https://doi.org/10.1186/gb-2013-14-1-r6