Abstract
Expression of oestrogen receptor (ESR1) determines whether a breast cancer patient receives endocrine therapy, but does not guarantee patient response. The molecular factors that define endocrine response in ESR1-positive breast cancer patients remain poorly understood. Here we characterize the DNA methylome of endocrine sensitivity and demonstrate the potential impact of differential DNA methylation on endocrine response in breast cancer. We show that DNA hypermethylation occurs predominantly at oestrogen-responsive enhancers and is associated with reduced ESR1 binding and decreased gene expression of key regulators of ESR1 activity, thus providing a novel mechanism by which endocrine response is abated in ESR1-positive breast cancers. Conversely, we delineate that ESR1-responsive enhancer hypomethylation is critical in transition from normal mammary epithelial cells to endocrine-responsive ESR1-positive cancer. Cumulatively, these novel insights highlight the potential of ESR1-responsive enhancer methylation to both predict ESR1-positive disease and stratify ESR1-positive breast cancer patients as responders to endocrine therapy.
Similar content being viewed by others
Introduction
The steroid hormone oestrogen activates the oestrogen receptor (ESR1) to mediate a variety of functions that are central to the normal development and maintenance of multiple tissues1. The unique transcriptional response to oestrogen in each tissue-specific cell subtype is, in part, regulated by the epigenome2. Differential DNA methylation and chromatin remodelling serve to dictate accessibility to functional, oestrogen-responsive regions of the genome, and thus define endocrine response3,4. Inappropriate activation of the ESR1 signalling network in mammary epithelial cells initiates neoplastic transformation and drives ESR1-positive breast cancer1. Patients with this disease commonly receive adjuvant endocrine therapy, which serves to inhibit ESR1 signalling1,5. Although endocrine therapy reduces the risk of disease recurrence, a third of patients acquire drug resistance and experience disease relapse6. Thus, endocrine sensitivity of both normal breast cells and breast cancer cells is dynamic, raising the hypothesis that global epigenetic reprogramming of oestrogen-responsive regions of the genome can modulate endocrine sensitivity and contributes to the onset of ESR1-positive breast cancer and the acquisition of endocrine resistance.
While recent studies have provided excellent proof of principle that the DNA methylation profile of mammary epithelial cells is altered in early carcinogenesis7, and further modified in cell models of endocrine-resistant breast cancer8,9, they do not address how these changes could directly affect endocrine sensitivity. Here we identify DNA methylation as a key determinant of endocrine response in breast cancer. We show that differential DNA hypermethylation occurs predominantly at oestrogen-responsive enhancer, not promoter regions, and is associated with reduced ESR1 binding and decreased gene expression of key regulators of ESR1 activity. In addition, we demonstrate that the methylation status of these regulatory regions is associated with endocrine resistance in human disease, thus providing a novel mechanism by which endocrine response is abated in ESR1-positive breast cancers.
Results
Methylation of enhancer loci in endocrine-resistant cells
To interrogate DNA methylation remodelling as a critical component of acquired endocrine resistance, we performed methylation profiling in duplicate using the Infinium HumanMethylation 450 beadchip, on ESR1-positive hormone sensitive MCF7 cells, and three different well-characterized endocrine-resistant MCF7-derived cell lines; tamoxifen-resistant (TAMR)10, fulvestrant-resistant (FASR)11 and oestrogen deprivation-resistant (MCF7X)12 cells. Density plots showing the correlation between the DNA methylation profile of parent MCF7 cells and individual endocrine-resistant cell lines indicate that the MCF7X and TAMR cells, which are both ESR1 positive10,12, predominantly gained DNA methylation as indicated by the increased density of points above the trend line. In contrast, FASR cells, which are ESR1 negative11, exhibited both hyper and hypomethylation events relative to parent MCF7 cells as indicated by a symmetrical density distribution (Fig. 1a–c). We first sought to identify the common differential DNA methylation events present in each of the three uniquely derived endocrine-resistant cell models by carrying out paired analyses (that is, each endocrine-resistant cell line versus MCF7 parent control) and overlapping the data (Fig. 1d). We found that across the individual resistant cell lines, 14,749 CpG probes were commonly hypermethylated (false discovery rate, FDR<0.01), whereas only 192 probes exhibited shared hypomethylation (FDR<0.01; Fig. 1d).
(a–c) A colorimetric density plot showing correlation between the HM450 methylation profile of the endocrine-resistant MCF7X (a), TAMR (b) and FASR (c) cells and the parent (endocrine-sensitive) MCF7 cells. The plots show that while the methylation profile of the endocrine-resistant cell lines is strongly correlated with the parent MCF7 cells (MCF7X, r2=0.895; TAMR, r2=0.91; FASR, r2=0.848; Pearson’s coefficient), both the MCF7X and TAMR cells predominantly gain DNA methylation, whereas the FASR cells exhibit both hyper- and hypomethylation events relative to parent MCF7 cells. (d) A Venn diagram showing the overlap of HM450 methylation probes that are more heavily methylated in multiple endocrine-resistant cells compared with the parent MCF7 cells (FDR<0.01). (e) A bar plot showing the association of differentially methylated HM450 probes that were common to all endocrine-resistant cell lines (compared with the parent MCF7 cells) across functional/regulatory regions of the genome as determined by MCF7 ChromHMM annotation13. The height of the bars represents the level of enrichment measured as a ratio between the frequency of hypermethylated (dark blue) or hypomethylated (light blue) probes overlapping a functional element over the expected frequency if such overlaps were to occur at random in the genome. Statistically significant enrichments (P value<<0.0001; hypergeometric test) are marked with an asterisk. The numbers of commonly hyper/hypomethylated probes located within each specific region are presented in the respective column.
To comprehensively characterize the functional genomic location of differential methylation observed in the endocrine-resistant cell models, we used ChromHMM segmentation of the MCF7 genome (previously described in Taberlay et al.13; Fig. 1e). Strikingly, significant enrichment of commonly hypermethylated probes was exclusively observed in enhancer regions of the genome (n=3,932 probes, P<<0.0001; hypergeometric test; Fig. 1e). We next sought to determine whether the enhancer regions identified as being more heavily methylated in all endocrine resistance models were regulated by the ESR1 in the parental MCF7 cells. Using reprocessed, publically available MCF7 ESR1 (ref. 14), GATA3 (ref. 15) and FOXA1 ChIP-Seq data16 (two transcription factors closely associated with ESR1 activity), we found that enhancer-specific CpG-hypermethylated probes were enriched in ESR1-binding sites by approximately sixfold, FOXA1-binding sites by fivefold and GATA3-binding sites by eightfold (P<<0.0001; hypergeometric test; Fig. 2a). The greatest number of hypermethylated enhancer probes were found to overlap ESR1-binding sites (n=801), which represents ∼20% of all hypermethylated probes in enhancer regions. Significantly, 47% (379 out of 801) of the hypermethylated enhancer probes that were located within an ESR1-binding site were also located within a FOXA1 and/or GATA3-binding site (Fig. 2b), which is particularly noteworthy since these transcription factors cooperatively modulate ESR1-transcriptional networks by forming a functional enhanceosome17.
(a) A bar plot showing the association of HM450 probes that were more heavily methylated in endocrine-resistant cell models (compared with MCF7 cells) and also specifically located in enhancer regions, across ESR1-, FOXA1- and GATA3-binding sites in MCF7 cells. The height of the bars represents the enrichment measured as a ratio between the frequency of hypermethylated probes in enhancers overlapping a transcription factor binding site over the expected frequency if such overlaps were to occur at random across the genome (*P value<<0.0001; hypergeometric test). The numbers of commonly hyper/hypomethylated probes located within each specific region are presented in the columns. (b) A Venn diagram showing the overlap of enhancer-specific HM450 methylation probes that are more heavily methylated in multiple endocrine-resistant cell models (compared with MCF7 cells) across ESR1-, FOXA1- and GATA3-binding sites. (c) A box plot showing the log-fold change (logFC) in ESR1 binding signal at ESR1-enhancer sites that contain at least one commonly hypermethylated probe (yellow box) and all other ESR1-enhancer sites that overlap a HM450 probe (grey box) in TAMR cells compared with the parent MCF7 cells. The mean logFC in ESR1 binding at hypermethylated ER-enhancer sites is −2.29 and the mean logFC of all other ESR1-enhancer sites is −0.52 (*P<<0.0001; t-test). (The whiskers of the box plot extend to the most extreme data point, which is no more than 1.5 × interquartile range from the box). (d) IGV screen shots to illustrate the loss of ESR1 binding in TAMR cells compared with the parent MCF7 cells in enhancer regions that overlap methylation probes that are more heavily methylated in the endocrine-resistant cell models. The MCF7 ChromHMM regions are colour coded as follows—blue, enhancer; yellow, transcribed; green, promoter; light blue, CTCF; and burgundy, transcribed. The HM450 β values are shown for the MCF7 (green), MCF7X (burgundy), TAMR (orange) and FASR cells (red) and are representative of biological duplicates. ESR1 ChIP data (blue) is presented in duplicate for both MCF7 and TAMR cells. The ESR1 enhancers that overlap the regions of endocrine-resistant-specific hypermethylation are highlighted by the blue boxes.
Enhancer DNA hypermethylation and diminished ESR1 binding
Having defined a subset of ESR1-binding sites that overlap enhancer regions which contain hypermethylated loci in multiple models of endocrine resistance (see Methods section; n=856 sites, Supplementary Data 1), we sought to determine whether DNA methylation affected the intensity of ESR1 binding at these sites. Using MCF7 and TAMR ESR1 ChIP data14, we compared the change in ESR1 binding signal intensity at ESR1-enhancer sites that contained (a) hypermethylated probe(s) to that of all other ESR1-enhancer sites (Fig. 2c). At methylated ESR1-enhancer sites, there was a 2.29-log-fold reduction in ESR1 binding in TAMR compared with MCF7 cells. In contrast, at all other ESR1-enhancer-binding sites, there was a 0.52-log-fold reduction in ESR1 binding in TAMR compared with MCF7 cells. Thus, increased methylation at ESR1-enhancer sites is associated with reduction in ESR1 binding (P<<0.0001; t-test; Fig. 2c). Four illustrative examples show the loss of ESR1 binding in the TAMR cells at enhancer regions that are more heavily methylated in the endocrine-resistant versus the parent MCF7 (Fig. 2d). The examples include enhancer regions located within the gene body of death-associated protein 6 (DAXX), golgi to ER traffic protein 4 homologue (GET4; a member of the BAG6-UBL4A-GET4 DNA damage response/cell death complex18), ESR1 itself and nuclear receptor co-repressor 2 (NCOR2; Fig. 2d).
Enhancer DNA hypermethylation and related gene expression
Since the vast majority of ESR1-enhancer-binding sites identified as hypermethylated in the endocrine-resistant cell lines compared with the parent MCF7 cells were intragenic (that is, 617 out of 856, 72% with at least partial overlap; Supplementary Data 1), we next sought to determine if the DNA methylation of these regions correlated with the expression of the genes in which they were located (or closest TSS if intergenic) in human breast cancer. Using RNA-seq and HM450 methylation data derived from TCGA breast cohort19 (n=459 patients), we determined that out of the 856 ESR1-enhancer-binding sites of interest, hypermethylation of 328 sites (that is, 38% of ESR1-enhancer sites) correlated with the reduced expression of the genes with which they were most closely associated (Spearman’s correlation coefficient; P<0.001; Supplementary Data 2). The 328 ESR1-enhancer-binding sites represented 291 unique genes (including those presented in Fig. 2d; Supplementary Data 3). Gene set enrichment analysis revealed that these genes were over-represented in gene sets upregulated by ESR1 activation, downregulated in the acquisition of endocrine resistance and gene sets lowly expressed in basal versus luminal disease, thus suggesting that such genes were critical drivers of oestrogen-driven tumours (Supplementary Fig. 1a). Interestingly, using unsupervised clustering analysis, this gene set (n=291) stratifies ESR1-positive and ESR1-negative breast cancer patients (Supplementary Fig. 1b). Cumulatively, this indicates that the methylation events occurring throughout the acquisition of endocrine resistance are serving to facilitate an oestrogen-independent phenotype reflective of a breast cancer subtype that is refractory to endocrine therapy.
ESR1-enhancer methylation defines breast cancer subtype
We next sought to determine whether ESR1-enhancer hypermethylation was indicative of breast cancer subtype. We assessed the median methylation of all hypermethylated ESR1-enhancer probes (n=801) in TCGA normal (n=97), luminal A (n=301), luminal B (n=52) and ESR1-negative (n=105) patient HM450 data (Fig. 3a). In normal breast tissue (which is reported to be ∼7% ESR1 positive20), the median methylation of the ESR1-enhancer sites was highest, while median DNA methylation was significantly reduced in luminal A disease (P<<0.0001; Mann–Whitney U-test), which is indicative of its endocrine-responsive state. Interestingly, median ESR1-enhancer methylation was greater in luminal B patients compared with luminal A patients (P=0.017; Mann–Whitney U-test), who are almost twice as likely to acquire endocrine resistance21. In ESR1-negative disease, median methylation was higher than in luminal disease (versus luminal A, P<<0.0001; versus luminal B, P<<0.0001; Mann–Whitney U-test; Fig. 3a). A heatmap highlights the hypomethylated status of the ESR1-enhancer sites in luminal A disease relative to normal breast tissue and the other breast cancer subtypes (Fig. 3b). This trend is clearly illustrated at the DAXX enhancer region in which each CpG within the ESR1-binding site was hypomethylated in luminal A disease compared with normal tissue and luminal B and ESR1-negative cancer (Fig. 3c). Critically, no such variability was apparent at the DAXX promoter region (1,000 bp upstream and 100 bp downstream of the transcription start site; Fig. 3c), suggesting a significant regulatory effect of increased methylation at the enhancer locus.
(a) A box plot showing the median methylation of all HM450 probes that overlap an enhancer region, an ESR1-binding site and demonstrate hypermethylation in endocrine-resistant versus parental MCF7 cells (n=801 probes), in normal breast tissue (green; n=97), luminal A (light blue; n=301), luminal B (dark blue; n=52) and ESR1-negative (red; n=105) breast cancer (data obtained from TCGA breast cancer cohort; *P<0.05, **P<<0.0001; Mann–Whitney U-test). (The whiskers of the box plot extend to the most extreme data point, which is no more than 1.5 × interquartile range from the box). (b) A heatmap showing the methylation profile of 801 ESR1-enhancer-specific HM450 probes that are more heavily methylated in endocrine-resistant versus parent MCF7 cells in normal breast tissue (green; n=97), luminal A (light blue; n=301), luminal B (dark blue; n=52) and ESR1-negative (red; n=105) breast cancer. Columns are patient samples and rows are HM450 probes. The level of methylation is represented by a colour scale—blue for low levels and red for high levels of methylation. (c) Box plots showing distribution of methylation β values in normal n=97 (green), luminal A (light blue; n=301), luminal B (dark blue; n=52) and ESR1-negative (red; n=105) breast cancer samples across HM450 probes overlapping the ESR1-binding site located within the DAXX enhancer (Chr6: 33288112-33288670; left panel) and the DAXX promoter region (1,000 bp upstream and 100 bp downstream of the transcription start site; Chr6: 33290693-33291793; right panel). (The whiskers of the box plots extend to the most extreme data point, which is no more than 1.5 × interquartile range from the box).
ESR1-enhancer hypermethylation predicts endocrine failure
Given that ESR1-enhancer hypermethylation is prevalent in acquired endocrine resistance in vitro (Figs 1e and 2a-d) and in molecular subclassifications of breast cancer that are intrinsically less responsive to endocrine therapy (Fig. 3a–c), we next sought to determine the methylation status of a panel of these loci in ESR1-positive (luminal A) breast cancer samples from patients with different outcomes. Primary samples were sourced from patients who received endocrine therapy for 5 years and either experienced relapse-free survival (RFS; >14 years) or those who had relapsed (<6 years), defined as no relapse-free survival (n/RFS). Matched local relapse samples were also compared with the primary n/RFS patient samples. All patients received the same endocrine therapy (tamoxifen; anonymized patient data is given in Supplementary Data 4). Using a multiplex bisulphite-PCR resequencing methodology specifically devised for formalin-fixed, paraffin-embedded (FFPE)-derived DNA22, the methylation of multiple CpG sites across a panel of nine oestrogen-responsive enhancer regions was interrogated (technical duplicate correlates for all amplicons investigated are shown in Supplementary Fig. 2). These enhancer regions included those located within DAXX, MSI2, NCOR2, RXRA and C8orf46 (Fig. 4a–e) and enhancer regions located within GATA3, ITPK1, ESR1 and GET4 (Supplementary Fig. 3a–d). The assay was repeated with DNA extracted from biological duplicates of the endocrine-resistant cell lines and the parent MCF7 cells to ensure its viability (Supplementary Fig. 4a–i; technical duplicate correlates for all amplicons investigated are shown in Supplementary Fig. 5). The average methylation levels detected at all enhancer loci were significantly higher in the recurrent tumours compared with the matched primary (n/RFS) tumours (DAXX, P<0.0001; ESR1, P<0.0005; RXRA, P<0.005; GET4, NCOR2, GATA3, MSI2, P<0.01; C8orf46, ITPK1, P<0.05; t-test), confirming that DNA methylation at ESR1-responsive enhancers is acquired in resistant disease (Fig. 4 and Supplementary Fig. 3). The difference in DNA methylation between RFS and n/RFS primary tumours was less considerable, although a statistically significant difference was observed for DAXX, P<0.0001; RXRA, P<0.01; C8orf46, P=0.01; NCOR2 and MSI2 (P<0.05; t-test) enhancer regions (Fig. 4).
(a–e) (Left panel) A scatter plot showing the methylation of individual CpG sites across the ESR1-enhancer region of interest ((a)-DAXX—Chr6: 33288296-33288372; (b)-MSI2—Chr17: 55371693-55371786; (c)-NCOR2—Chr12: 124844786-124844883; (d)-RXRA—Chr9: 137252867-137252967; (e)-C8orf46—Chr8: 67425069-67425134) in three primary luminal A breast cancers from patients that received adjuvant endocrine therapy and exhibited RFS (green), three primary luminal A breast cancers from patients that relapsed following adjuvant endocrine therapy, defined as no n/RFS (blue) and their matched local relapse (red). Each dot represents the % methylation at an individual CpG site for a single patient and the lines represent the average methylation for the region in primary RFS (green), primary n/RFS (blue) and matched recurrent tumours (red). (Right panel) Box plots showing the distribution of methylation values across the ESR1-enhancer region depicted in the left panel for RFS (green), prognosis/RFS (blue) and matched recurrent tumours (red); P values correspond to t-test comparison between RFS versus n/RFS, and n/RFS versus relapse tumours. (The whiskers of the box plots extend to the most extreme data point, which is no more than 1.5 × interquartile range from the box).
Discussion
Our results support a model whereby ESR1-responsive enhancer DNA methylation is a fundamental unifying characteristic that defines endocrine sensitivity in breast cancer. Interestingly, previous studies interrogating DNA methylation changes in endocrine-resistant cell models have predominantly reported ESR1-regulated promoter methylation8,9,23,24,25,26. Our study is the first to combine in depth MCF7 ChromHMM annotation and genome-wide methylation data from multiple resistance models to more comprehensively characterize global differential methylation across diverse genomic regions. We show for the first time that the methylation status of enhancers is associated with the inhibition of ESR1 binding in vitro and with the reduced expression of critical regulators and effectors of ESR1 activity in human disease. The identification of ESR1-responsive enhanceosome hypermethylation is both novel and considerably pertinent in the context of endocrine resistance, since genome-wide positional analyses defining the set of cis-regulatory elements that recruit ESR1 in breast cancer cells have revealed its predominant recruitment to enhancers as opposed to promoter regions3,27,28,29,30. Enhancers are more common than promoters in the mammalian genome31 and can regulate gene transcription from tens to thousands of kilobases away by promoting communication with target promoters through chromatin looping32,33. In our study, the majority of ESR1-regulated enhancer regions identified as hypermethylated in the resistant cells were located within gene bodies. Strikingly, hypermethylation of these enhancer regions was frequently correlated with reduced expression of the host gene, which is in line with recent studies that have shown that over half of all enhancer regions are located within a gene body and that the activation of these enhancers can indeed affect the transcription of the host gene34,35. Examples of genes whose expression inversely correlated with ESR1-enhancer DNA methylation include DAXX and GET4, which have been previously associated with roles in apoptosis18,36. It is conceivable that the loss of expression of genes associated with pro-apoptotic functions facilitates the progression of endocrine resistance by reducing the efficacy of apoptotic signalling pathways activated by endocrine therapies37.
Importantly, the ESR1-responsive enhancer hypermethylation events identified in the endocrine-resistant cell lines were also differentially methylated in endocrine-sensitive and endocrine-resistant breast cancer patient samples. Therefore, it is feasible that ESR1-responsive enhancer methylation status is reflective of endocrine dependence and could potentially be used to stratify patients as responders to endocrine therapy. For example, NCOR2, a gene whose expression has previously been associated with metastasis-free survival in 620 lymph node-negative patients with ESR1-positive breast cancer38, was shown to negatively correlate with ESR1-enhancer methylation. In the present study, NCOR2 enhancer methylation was significantly higher in the poor (non-relapse-free) prognosis patients, compared with the good (relapse-free) prognosis primary luminal A breast cancer patients. Critically, however, in matched recurrent tumours, enhancer DNA methylation was further increased, supporting the hypothesis that the endocrine-resistant methylation profile is acquired, rather than pre-existing, limiting its potential prognostic value. Intriguingly, it could be a combination of both acquired and intrinsic methylation differences that give rise to endocrine-resistant disease. One possible explanation is that sparse, or ‘seeding’ methylation at ESR1-responsisive enhancer sites in primary tumours could reflect a propensity to gain extensive methylation that spreads as resistance develops, which then becomes firmly established in recurrent disease (as discussed in ref. 39). Further characterization of ESR1-responsive enhancer methylation in endocrine-resistant disease will hereafter be an important area of future investigation, as will be the assessment of its potential predictive and prognostic application in breast cancer.
Methods
Cell culture and HumanMethylation450K array
MCF7 breast cancer cells and the corresponding endocrine-resistant sub-cell lines were kindly given to our laboratory by Dr Julia Gee (Cardiff University, UK). In brief, MCF7 cells were maintained in RPMI-1640-based medium containing 5% (v/v) fetal calf serum (FCS). TAMR MCF7 cells were generated by the long-term culture of MCF7 cells in phenol red-free RPMI medium containing 5% charcoal stripped FCS and 4-OH-tamoxifen (1 × 10−7 M; TAM). FASR MCF7 cells were generated by the long-term culture of MCF7 cells in phenol red-free RPMI medium containing 5% charcoal stripped FCS and fulvestrant (1 × 10−7 M; FAS). Long-term oestrogen-deprived MCF7 (MCF7X) cells were generated by the long-term culture of MCF7 cells in phenol red-free RPMI medium containing 5% charcoal stripped FCS. Endocrine-resistant sub-lines were established and characterized following 6 months endocrine challenge/oestrogen deprivation exposure10,11,12. All cell lines were authenticated by short-tandem repeat profiling (Cell Bank, Australia) and cultured for <6 months after authentication. Genomic DNA was extracted using the Qiagen DNeasy Blood and Tissue kit according to the manufacturer’s instructions. HumanMethylation450K arrays were carried out by the Australian Genome Research Facility (AGRF; Melbourne, Australia).
HM450 analysis
Two biological replicates per condition—MCF7, TAMR, MCF7X or FASR—were profiled on Illumina’s HumanMethylation450K array. Raw HM450 data was preprocessed and background normalized with the Biconductor minfi package40 using preprocess Illumina(..., bg.correct=TRUE, normalize=’controls’, reference=1); resulting M values were used for statistical analyses and β values for heatmap visualizations and clustering. Differential methylation analysis of the preprocessed data was performed using the Bioconductor limma package.
Genomic segmentation and annotation
The ChromHMM segmentation of the MCF7 genome was obtained from Taberlay et al.13. Enhancer (‘Enhancer’ and ‘Enhancer+CTCF’) and Promoter categories (‘Promoter’, ‘Promoter+CTCF’ and ‘Poised Promoter’) were collapsed into a single ‘Enhancer’ and ‘Promoter’ state respectively for the purposes of our analysis. RefSeq transcript annotations were obtained from UCSC genome browser41,42.
ChIP-seq data acquisition and analysis
ESR1 ChiP-seq data for ESR1 in MCF7 and TAMR14 was utilized in this study. Reads were mapped to genome build HG19 (GRCh37) with bowtie and mismatched (>3 mismatched bases), multiple mapping and duplicate reads were excluded from downstream analysis. ESR1 enrichment peaks were identified with the HOMER software suite43 using the findPeaks utility (-style factor -fragLength 200 -size 300 -F 0 -L 0 -C 0 -poisson 1e-06) on each experiment separately. We merged the resulting peaks to produce a ground set of 120,735 regions for subsequent analysis. Active ESR1 regions were identified in MCF7 by comparing the distribution of reads overlapping the ground set of ESR1 regions in the three MCF7 ESR1 experiments (GSM798423, GSM798424 and GSM798425) and MCF7 input experiment (GSM798440) with edegR44. This yielded 54,265 active ESR1 regions in MCF7 (FDR<0.05). A similar strategy was applied to TAMR data to yield 49,511 ESR1 regions in TAMR cells. Regions of differential ESR1 binding were identified by comparing the distribution of sequence reads in MCF7 and TAMR across the ground set of ESR1 regions using edgeR and potential variation in copy number was accounted for using DiffBind14. This analysis resulted in 24,711 regions with statistical significant gain (FDR 5%) and 32,343 regions with statistical significant loss (FDR 5%) of ESR1 binding in TAMR cells as compared to MCF7 cells. ESR1 peaks overlapping HM450 probes were assigned to the nearest RefSeq transcript (<20 kb distance) for the purposes of gene expression analysis. Raw MCF7 GATA3 and FOXA1 ChIP-Seq data was obtained from Theodorou et al.15 and Hurtado et al.16, respectively. Data were processed in the same manner as outlined for ESR1 ChIP-seq above.
Gene set enrichment analysis
GSEA was performed against the Molecular Signatures Database v4.0 (MSigDB)45 C2 Collection. Enrichment was assessed by hypergeometric testing as implemented in the R stats package.
TCGA data acquisition
DNA methylation analysis utilized clinical data available through the TCGA Breast Invasive Carcinoma cohort19. Raw HM450 methylation data (level 1) were obtained from the TCGA data portal (normal samples=97, ESR1-positive tumours=353 and ESR1-negative tumours=105). ESR1-positive tumours were further divided into luminal A (lumA=301) and luminal B (lumB=52) populations using progesterone receptor (PR) expression, such that lumA were ESR1+/PR+ and lumB were ESR1+/PR−. Processed RNA-Seq expression data (level 3) were obtained from TCGA data portal (588 ESR1 positive tumours with 73 matched normals and 174 ESR1 negative samples with 19 matched normals).
Multiplex bisulfite-PCR resequencing of clinical FFPE DNA
Bisulfite DNA conversions were performed using a manual protocol. For each conversion, ∼100 ng was bisulfite converted at a time. Conversion took place at 80 °C for 45 min in the presence of 0.3 M NaOH, 3.75 mM quinone and 2.32 M sodium metabisulfite, as per the method of Clark et al.46. The multiplex bisulfite-PCR reaction was performed as follows22. In brief, Promega HotStart GoTaq with Flexi buffer (M5005) was used with the following components at the indicated concentrations: 5 × green (1 × ), CES 5 × , (0.5 × , N.B. refer to ref. 47 for CES recipe), MgCl2 (4.5 mM), dNTP’s (200 μM each), primers (forward and reverse at 100 mM), Hot Start Taq (0.025 U μl−1), DNA (2 ng μl−1 ). All primers used are listed in Supplementary Data 5. Cycling conditions were: 94 °C, 5 min; 12 cycles of (95 °C, 20 s; 60, 1 min); 12 cycles of (94 °C, 20 s; 65 °C, 1 min 30 s); 65 °C, 3 min, 10 hold. Agencourt XP beads were using to clean-up and concentrate the multiplex reaction for subsequent barcoding (that is, addition of Illumina p5/p7 sequences and sample-specific DNA barcodes). The barcoding PCR used the following reagents at the indicated final concentrations in a 100-μl reaction: 1 × GoTaq Green Flexi buffer; 0.25 × CES; 4.5 mM MgCl2; 200 μM dNTPs; 0.05 U μl−1 HotStart Taq; 25 μl of pooled template after Agencourt XP bead clean-up; and 20 μl MiSeq (Fluidigm PN FLD-100-3771). Cycling conditions were: 94 °C, 5 min; 9 cycles of (97 °C, 15 s; 60 °C, 30 s; 72 °C, 2 min); 72 °C, 2 min; 6 °C, 5 min. MiSeq sequencing was performed used the MiSeq Reagent Kit v2, 300 cycle; PN MS-102-2002. Bioinformatic analysis started with adaptor trimming using Trim galore (options: --length 100). Mapping used the Bismark methylation mapping programme48 running Bowtie2 (ref. 49) (options: --bowtie2 -N 1 -L 15 --bam -p 2 --score L,-0.6,-0.6 --non_directional; bismark_methylation_extractor -s -merge_non_CpG –comprehensive --cytosine_report). To reduce computational overhead, mapping took place against only those genomic regions which were being investigated, plus an additional 100 bp–1 kb of flanking sequence.
Clinical sample acquisition and DNA extraction
FFPE breast cancer samples were obtained from the St George Hospital, Kogarah, Australia (Ethics approval reference from St George Hospital Human Research Ethics Committee is 96/84). The deidentified haematoxylin–eosin-stained sections were reviewed by a pathologist and representative tumour areas were marked and blocks were cored accordingly. Genomic DNA was extracted using the Qiagen AllPrep DNA/RNA FFPE kit according to the manufacturer’s instructions.
Additional information
Accession codes: Cell line HumanMethylation450K array data is available online at GEO (GSE69118).
How to cite this article: Stone, A. et al. DNA methylation of oestrogen-regulated enhancers defines endocrine sensitivity in breast cancer. Nat. Commun. 6:7758 doi: 10.1038/ncomms8758 (2015).
Accession codes
References
Musgrove, E. A. & Sutherland, R. L. Biological determinants of endocrine resistance in breast cancer. Nat. Rev. Cancer 9, 631–643 (2009).
Biddie, S. C., John, S. & Hager, G. L. Genome-wide mechanisms of nuclear receptor action. Trends Endocrinol. Metab. 21, 3–9 (2010).
Lupien, M. et al. FoxA1 translates epigenetic signatures into enhancer-driven lineage-specific transcription. Cell 132, 958–970 (2008).
Joseph, R. et al. Integrative model of genomic factors for determining binding site selection by estrogen receptor‐α. Mol. Syst. Biol. 6, 456 (2010).
Ignatiadis, M. & Sotiriou, C. Luminal breast cancer: from biology to treatment. Nat. Rev. Clin. Oncol. 10, 494–506 (2013).
EBCTCG. Effects of chemotherapy and hormonal therapy for early breast cancer on recurrence and 15-year survival: an overview of the randomised trials. Lancet 365, 1687–1717 (2005).
Locke, W. et al. Coordinated epigenetic remodelling of transcriptional networks occurs during early breast carcinogenesis. Clin. Epigenet. 7, 1–15 (2015).
Lin, X. et al. Integrative analyses of gene expression and DNA methylation profiles in breast cancer cell line models of tamoxifen-resistance indicate a potential role of cells with stem-like properties. Breast Cancer Res. 15, R119 (2013).
Williams, K. E., Anderton, D. L., Lee, M. P., Pentecost, B. T. & Arcaro, K. F. High-density array analysis of DNA methylation in tamoxifen-resistant breast cancer cell lines. Epigenetics 9, 297–307 (2014).
Knowlden, J. M. et al. Elevated levels of epidermal growth factor receptor/c-erbB2 heterodimers mediate an autocrine growth regulatory pathway in tamoxifen-resistant MCF-7 cells. Endocrinology 144, 1032–1044 (2003).
McClelland, R. A. et al. Enhanced epidermal growth factor receptor signaling in MCF7 breast cancer cells after long-term culture in the presence of the pure antiestrogen ICI 182,780 (Faslodex). Endocrinology 142, 2776–2788 (2001).
Staka, C. M., Nicholson, R. I. & Gee, J. M. Acquired resistance to oestrogen deprivation: role for growth factor signalling kinases/oestrogen receptor cross-talk revealed in new MCF-7X model. Endocr. Relat. Cancer 12, S85–S97 (2005).
Taberlay, P. C., Statham, A. L., Kelly, T. K., Clark, S. J. & Jones, P. A. Reconfiguration of nucleosome depleted regions at distal regulatory elements accompanies DNA methylation of enhancers and insulators in cancer. Genome Res. 24, 1421–1432 (2014).
Ross-Innes, C. S. et al. Differential oestrogen receptor binding is associated with clinical outcome in breast cancer. Nature 481, 389–393 (2012).
Theodorou, V., Stark, R., Menon, S. & Carroll, J. S. GATA3 acts upstream of FOXA1 in mediating ESR1 binding by shaping enhancer accessibility. Genome Res. 23, 12–22 (2013).
Hurtado, A., Holmes, K. A., Ross-Innes, C. S., Schmidt, D. & Carroll, J. S. FOXA1 is a key determinant of estrogen receptor function and endocrine response. Nat. Genet. 43, 27–33 (2011).
Kong, S. L., Li, G., Loh, S. L., Sung, W. K. & Liu, E. T. Cellular reprogramming by the conjoint action of ERα, FOXA1, and GATA3 to a ligand‐inducible growth state. Mol. Syst. Biol. 7, 526 (2011).
Krenciute, G. et al. Nuclear BAG6-UBL4A-GET4 complex mediates DNA damage signaling and cell death. J. Biol. Chem. 288, 20547–20557 (2013).
TCGA. Comprehensive molecular portraits of human breast tumours. Nature 490, 61–70 (2012).
Petersen, O. W., Høyer, P. E. & van Deurs, B. Frequency and distribution of estrogen receptor-positive cells in normal, nonlactating human breast tissue. Cancer Res. 47, 5748–5751 (1987).
Kennecke, H. et al. Metastatic behavior of breast cancer subtypes. J. Clin. Oncol. 28, 3271–3277 (2010).
Korbie, D. et al. Multiplex bisulfite PCR resequencing of clinical FFPE DNA. Clin. Epigenet. 7, 28 (2015).
Stone, A. et al. BCL-2 Hypermethylation is a potential biomarker of sensitivity to antimitotic chemotherapy in endocrine-resistant breast cancer. Mol. Cancer Ther. 12, 1874–1885 (2013).
Stone, A. et al. Tamoxifen-induced epigenetic silencing of oestrogen-regulated genes in anti-hormone resistant breast cancer. PLoS ONE 7, 10 (2012).
Miller, T. W. et al. Hyperactivation of phosphatidylinositol-3 kinase promotes escape from hormone dependence in estrogen receptor–positive human breast cancer. J. Clin. Invest. 120, 2406–2413 (2010).
Pathiraja, T. N. et al. Epigenetic reprogramming of HOXC10 in endocrine-resistant breast cancer. Sci. Transl. Med. 6, 229ra241 (2014).
Carroll, J. S. et al. Chromosome-wide mapping of estrogen receptor binding reveals long-range regulation requiring the forkhead protein FoxA1. Cell 122, 33–43 (2005).
Carroll, J. S. et al. Genome-wide analysis of estrogen receptor binding sites. Nat. Genet. 38, 1289–1297 (2006).
Lin, C.-Y. et al. Whole-genome cartography of estrogen receptor α binding sites. PLoS Genet. 3, e87 (2007).
Liu, Y. et al. The genome landscape of ERα- and ERβ-binding DNA regions. Proc. Natl Acad. Sci. USA 105, 2604–2609 (2008).
Bulger, M. & Groudine, M. Enhancers: the abundance and function of regulatory sequences beyond promoters. Dev. Biol. 339, 250–257 (2010).
Bulger, M. & Groudine, M. Functional and mechanistic diversity of distal transcription enhancers. Cell 144, 327–339 (2011).
Fullwood, M. J. et al. An oestrogen-receptor-[agr]-bound human chromatin interactome. Nature 462, 58–64 (2009).
Heintzman, N. D. et al. Histone modifications at human enhancers reflect global cell-type-specific gene expression. Nature 459, 108–112 (2009).
Kowalczyk Monika, S. et al. Intragenic enhancers act as alternative promoters. Mol. Cell 45, 447–458 (2012).
Yang, X., Khosravi-Far, R., Chang, H. Y. & Baltimore, D. Daxx, a novel Fas-binding protein that activates JNK and apoptosis. Cell 89, 1067–1076 (1997).
Thiantanawat, A., Long, B. J. & Brodie, A. M. Signaling pathways of apoptosis activated by aromatase inhibitors and antiestrogens. Cancer Res. 63, 8037–8050 (2003).
van Agthoven, T. et al. CITED2 and NCOR2 in anti-oestrogen resistance and progression of breast cancer. Br. J. Cancer 101, 1824–1832 (2009).
Stone, A., Valdes-Mora, F. & Clark, S. J. Exploring and exploiting the aberrant DNA methylation profile of endocrine-resistant breast cancer. Epigenomics 5, 595–598 (2013).
Aryee, M. J. et al. Minfi: a flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinformatics 30, 1363–1369 (2014).
Kent, W. J. et al. The Human Genome Browser at UCSC. Genome Res. 12, 996–1006 (2002).
Meyer, L. R. et al. The UCSC Genome Browser database: extensions and updates 2013. Nucleic Acids Res. 41, D64–D69 (2013).
Heinz, S. et al. Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. Mol. Cell 38, 576–589 (2010).
Robinson, M. D., McCarthy, D. J. & Smyth, G. K. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139–140 (2010).
Subramanian, A. et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl Acad. Sci. USA 102, 15545–15550 (2005).
Clark, S. J., Statham, A., Stirzaker, C., Molloy, P. L. & Frommer, M. DNA methylation: bisulphite modification and analysis. Nat. Protoc. 1, 2353–2364 (2006).
Ralser, M. et al. An efficient and economic enhancer mix for PCR. Biochem. Biophys. Res. Commun. 347, 747–751 (2006).
Krueger, F. & Andrews, S. R. Bismark: a flexible aligner and methylation caller for Bisulfite-Seq applications. Bioinformatics 27, 1571–1572 (2011).
Langmead, B. & Salzberg, S. L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357–359 (2012).
Acknowledgements
We thank Dr Jason Carroll and members of his lab for providing raw ChIP-Seq data used in this study, and Dr Brigid O’Gorman for her careful review of the manuscript. We also acknowledge the input of the late Professor Robert L. Sutherland. This work is supported by the National Breast Cancer Foundation (NBCF) programme and project grants and National Health and Medical Research Council (NHMRC 1029579) project grant and NHMRC Fellowship (S.J.C.). In addition, J.M.W.G. received support from Breast Cancer Campaign UK and E.A.M. received support from The Cancer Institute NSW (11/CDF/3-26) and Cancer Research UK (C596/A18076).
Author information
Authors and Affiliations
Contributions
A.S., E.Z., S.J.C. were responsible for the concept and designA.S., E.Z., D.K., W.J.L. and S.J.C. developed the methods A.S., E.Z., W.J.L., D.K., E.A.M. and S.J.C. were responsible for acquisition of data; A.S., E.K.A.M. and P.G. provided clinical samples (preparation of DNA); E.Z., W.J.L., A.S., D.K., R.P., S.J.C. and C.S. analysed and interpreted the data (for example, statistical analysis, biostatistics and computational analysis); all authors, with the exception of P.G., were responsible for writing, figures and review of the manuscript; E.A.M, C.S. and S.J.C supervised the study..
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Supplementary information
Supplementary Figures
Supplementary Figures 1-5 (PDF 909 kb)
Supplementary Data 1
ESR1 binding sites that contain hypermethylated loci in multiple models of endocrine resistant breast cancer. (XLS 1792 kb)
Supplementary Data 2
Correlation of DNA hypermethylation at ESR1 enhancer loci and differentail gene expression in TCGA breast cohort. (XLS 174 kb)
Supplementary Data 3
Genes whose expression negatively correlate with increased ESR1 enhancer DNA hypermethylation in multiple models of endocrine resistant breast cancer (n = 291 unique genes). (XLS 96 kb)
Supplementary Data 4
Clinical parameters of patient samples used. (XLS 19 kb)
Supplementary Data 5
Primer Sequences. (XLS 19 kb)
Rights and permissions
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
About this article
Cite this article
Stone, A., Zotenko, E., Locke, W. et al. DNA methylation of oestrogen-regulated enhancers defines endocrine sensitivity in breast cancer. Nat Commun 6, 7758 (2015). https://doi.org/10.1038/ncomms8758
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/ncomms8758