Skip to main content
  • Research article
  • Open access
  • Published:

Global analysis of estrogen receptor beta binding to breast cancer cell genome reveals an extensive interplay with estrogen receptor alpha for target gene regulation

Abstract

Background

Estrogen receptors alpha (ERα) and beta (ERβ) are transcription factors (TFs) that mediate estrogen signaling and define the hormone-responsive phenotype of breast cancer (BC). The two receptors can be found co-expressed and play specific, often opposite, roles, with ERβ being able to modulate the effects of ERα on gene transcription and cell proliferation. ERβ is frequently lost in BC, where its presence generally correlates with a better prognosis of the disease. The identification of the genomic targets of ERβ in hormone-responsive BC cells is thus a critical step to elucidate the roles of this receptor in estrogen signaling and tumor cell biology.

Results

Expression of full-length ERβ in hormone-responsive, ERα-positive MCF-7 cells resulted in a marked reduction in cell proliferation in response to estrogen and marked effects on the cell transcriptome. By ChIP-Seq we identified 9702 ERβ and 6024 ERα binding sites in estrogen-stimulated cells, comprising sites occupied by either ERβ, ERα or both ER subtypes. A search for TF binding matrices revealed that the majority of the binding sites identified comprise one or more Estrogen Response Element and the remaining show binding matrixes for other TFs known to mediate ER interaction with chromatin by tethering, including AP2, E2F and SP1. Of 921 genes differentially regulated by estrogen in ERβ+ vs ERβ- cells, 424 showed one or more ERβ site within 10 kb. These putative primary ERβ target genes control cell proliferation, death, differentiation, motility and adhesion, signal transduction and transcription, key cellular processes that might explain the biological and clinical phenotype of tumors expressing this ER subtype. ERβ binding in close proximity of several miRNA genes and in the mitochondrial genome, suggests the possible involvement of this receptor in small non-coding RNA biogenesis and mitochondrial genome functions.

Conclusions

Results indicate that the vast majority of the genomic targets of ERβ can bind also ERα, suggesting that the overall action of ERβ on the genome of hormone-responsive BC cells depends mainly on the relative concentration of both ERs in the cell.

Background

Estrogens are key regulators of cell growth and differentiation in the mammary gland [1, 2], where they are involved in the pathogenesis and clinical outcome of breast cancer (BC) [3]. These steroid hormones exert their effects in normal and transformed mammary epithelial cells by binding to specific receptors, ERα and ERβ, that mediate estrogen signaling by functioning as ligand-dependent transcription factors. Ligand-activated ERs drive gene cascades comprising primary genes, whose transcription is directly controlled by the hormone through physical interaction of ERs with regulatory sites in the genome (genomic pathway) and/or with signal transduction effectors (non genomic pathway), as well as downstream genes whose activity depends upon the functions encoded by the primary responders [1, 4].

ERs are able to bind DNA at specific sites in the genome and thereby control gene activity by recruiting transcriptional mediators and co-regulators, as well as a host of other nuclear proteins with different roles in ER-mediated control of gene activity [5, 6]. The two ERs show 55% identity in their estrogen-binding domains (LBDs) and approximately 97% similarity in the DNA-binding domains (DBDs) [7]. Reflecting the high degree of similarity in their DBDs, both receptors interact with the same conserved estrogen response element (ERE) (5'-GGTCAnnnTGACC-3') as either homodimers or alpha/beta heterodimers [8, 9]. ERβ, however, holds low trans-acting capability on ERE-containing estrogen target genes and alpha/beta heterodimers are less efficient than ERα homodimers in promoting target genes activity [10]. The different behaviour of ERα/ERβ heterodimers respect to ERα homodimers on transcriptional regulation of ERE-containing genes might be explained by different co-factor recruitment, as ERβ could prevent efficient co-activator binding to the ERα moiety of the heterodimer, conversely inducing recruitment of co-repressors and/or driving assembly of co-regulatory complexes other than those involving ERα only. [8, 11, 12].

Although the two receptors are quite similar in sequence and structure, in BC ERβ has considerably different biological effects than ERα [1, 13, 14]. Furthermore, the two ERs show a remarkably different expression pattern in BCs, with higher ERα and lower ERβ levels observed in malignant cells compared to normal mammary epithelial or benign tumor cells [15, 16]. Furthermore, while ERα induces a mitogenic response to estrogen, when expressed alone the β subtype is not only unable to induce the same mitogenic response, but it reduces basal, hormone-independent cell proliferation [17-18 and R. Tarallo et al., unpublished]. Finally, ERβ was shown to change dramatically ERα-positive BC cell behaviour in vivo, as its expression in the cell prevents tumorigenicity in mouse xenograft models by reducing tumor growth and angiogenesis [19, 20].

Gene expression studies performed in BC cell lines expressing endogenous ERα and recombinant ERβ [21–23] revealed multiple signaling pathways involving the α and/or β receptor subtypes [1]. The two ERs appear thus to share many target genes, although each of them may affect specific downstream targets. For this reason, inhibition of hormone-responsive BC cell growth by ERβ might be due to direct interference with ERα activity on growth-promoting pathways as well as to the activity of ERβ-specific target genes [24].

Recently, next-generation sequencing technologies combining chromatin immunoprecipitation (ChIP) either with genomic DNA hybridization to microarrays (ChIP-on-chip) or massively parallel sequencing (ChIP-Seq, ChIP-PET), opened new venues for our understanding of physical and functional associations between transcription factors and chromatin in vivo. These analytical strategies led to genome-wide mapping of ERα-binding regions in intact chromatin of cultured cell lines [25–28], revealing important new information relative to ERα interaction with the genome. Carroll et al.[25], for example, using ChIP-on-chip demonstrated that the Forkhead factor FoxA1 plays an important role as pioneering factor for ERα binding to chromatin in BC cells, while Cicatiello et al.[26] identified novel gene regulation cascades mediating estrogen actions in hormone-responsive BC cells. In contrast, although several studies focused on ERβ interaction with the genome [29–32], a thorough characterization of this important aspect of ERβ biology in BC cells, essential to clarify the mechanisms mediating its control of estrogen-dependent gene pathways and the hormone-responsive phenotype, is still missing. For this reason, we performed a comprehensive analysis of ERβ and ERα target sites in the genome of MCF-7 cells engineered to express both receptors to comparable levels, by integrating global mapping of in vivo ER binding to the genome by ChIP-Seq with comparative gene expression profiling in ERβ+/ERα+ vs ERβ-/ERα+ cells during early stimulation with 17-beta-estradiol (E2), followed by in silico analyses of the ERβ binding regions and responsive genes identified.

Results and Discussion

Establishment and characterization of ERβ-expressing MCF-7 cells

Stabilized human BC cell lines expressing endogenous ERα and ERβ at comparable levels are not available. For this reason, we first generated and characterized cell lines derived from ERα-positive MCF-7 cells expressing full-length human ERβ (ERβ1) at levels similar to those of endogenous ERα. This strategy was adopted to prevent artefacts due to ERβ over-expression in the cell and to mimic what observed in primary breast tumors, where very high expression of this receptor has never been observed. As suitable antibodies for efficient immunoprecipitation of chromatin-bound ERβ are not available, the expressed proteins were tagged on either their C- (Ct-ERβ) or N-terminus (Nt-ERβ) with the TAP epitope. This approach allows to track tagged ERs in different cell compartments and to efficiently immunoprecipitate and purify them in vitro by Tandem Affinity Purification, to identify their molecular partners [5, 33], and ex vivo in chromatin immunoprecipitation assays (see below). Preliminary tests were performed to verify whether the presence of the TAP moiety could influence intracellular redistribution of ERβ in response to 17β-estradiol (E2) and its ability either to trans-activate an estrogen-responsive reporter gene or to interfere with ERα activity on reporter gene transcription and BC cell proliferation. To this end, ER expression and nuclear translocation in response to E2 was determined in wt MCF-7, Ct-ERβ and Nt-ERβ cells by subcellular fractionation followed by SDS-PAGE and immunoblotting (Figure 1A). In absence of hormone a larger fraction of both ERs was found in the cytosol in all cases. Following estrogen stimulation, both receptors migrated to the nucleus, a crucial event to trigger target gene transcription via the genomic pathway of the estrogen signaling cascade. An antibody against α-tubulin was used as control, and the absence of this protein in the nuclear fractions indicates that they were indeed free from cytosolic contaminants. The ERβ-expressing clones selected for this study showed a ERβ/ERα ratio <2, as verified by immunoblotting analysis of the proteins in whole cell extracts and quantitative rtPCR of the corresponding RNAs [5, and data not shown]. To control that the presence of the TAP tag did not interfere with ERβ activity, ER-negative SKBR-3 BC cells were transiently transfected with expression vectors encoding wt ERβ, Ct-ERβ, Nt-ERβ, ERα (HEG0) or 'empty' vector (pSG5), as controls, and ERE-tk-luc [34], a reporter gene where luciferase expression is driven by an estrogen-responsive minimal promoter. Exposure of transiently transfected cells to E2 induced reporter gene activation in the presence of ERα, ERβ, Ct-ERβ or Nt-ERβ, with the activity of both tagged ERβ proteins slightly (15-20%) lower than that of wt ERβ (Figure 1B, left). We then tested whether the two recombinant forms of ERβ were able to interfere with target gene activation by the endogenous ERα resident in MCF-7 cells. To this end, wt, Ct-ERβ+ and Nt-ERβ+ cells were transfected with ERE-tk-luc and the response of the reporter gene to E2 was determined. As shown in Figure 1B (right), ERβ-expressing cells showed in all cases a marked (50-60%) reduction in reporter gene response to the hormone, when compared to wt MCF-7 cells, indicating that both tagged ERβs are able to interfere with the activity of endogenous ERα. Results show that cell lines stably expressing Ct-ERβ and Nt-ERβ display a marked reduction in proliferative response to the hormone, when compared to wt MCF-7 cells (Figure 1C), in agreement with the known effects of ERβ in ERα-positive cells [23, 35–37]. Furthermore, comparative RNA expression profiling in exponentially growing Ct-ERβ and Nt-ERβ vs wt MCF-7 cells revealed extensive overlapping effects of the two tagged ERβ proteins on the activity of estrogen target genes [O. Paris et al., manuscript in preparation and data not shown]. Taken together, these observations confirmed that both tag-ERβ expressing cell lines generated for this study show a well defined phenotype, with respect to the known activities of this ER subtype in BC cells, and are thus suitable to investigate the genomic bases of ERβ actions in this cell type. As we could not exclude that the presence of the TAP tag at either the N- or C- term of ERβ may specifically influence its activity on cellular targets or pathways different from those investigated above, all experiments reported in this study were performed in both Ct-ERβ and Nt-ERβ cells and the data were combined for analysis, with the aim to focus on the most significant and reproducible actions of ERβ independently from position of the tag in the receptor moiety.

Figure 1
figure 1

Functional characterization of ERβ-expressing MCF-7 cells. (A) Nuclear translocation of ERα and ERβ shown by western blot analysis on cytosolic (c) and nuclear (n) protein extracts, prepared from wt MCF-7, Nt-ERβ and Ct-ERβ cells after treatment with either 17β-estradiol (10-8M, +E2) or vehicle alone (-E2) for 45 minutes. The amount of α-tubulin was also analyzed to verify the absence of cytosolic contaminants in the nuclear fractions. (B) The transcriptional activity of ERα, Nt-ERβ and Ct-ERβ was measured by transient transfection in SKBR3 cells (left) and the ability of tagged ERβ to interfere with ERα activity was assessed by comparing estrogen effects in wt, in Nt-ERβ and Ct-ERβ MCF-7 cells (right); in all cases transiently transfected ERE-tk-luc was used as reporter gene. (C) Proliferation of wt MCF-7, Nt-ERβ and Ct-ERβ cells was measured by stimulating hormone-starved cells with 10-8M E2, followed by cell counting with a colorimetric assay at the indicated times.

Effects of ERβ on the estrogen-responsive MCF-7 cell transcriptome

Expression of ERβ is known to cause significant changes in the genomic response to estrogen in target cells. To identify the genes whose estrogen regulation in hormone-responsive BC cells is perturbed by ERβ, we performed gene expression profiling with microarrays in estrogen-starved, quiescent wt and TAP-ERβ+ (Ct-ERβ and Nt-ERβ) MCF-7 cells following stimulation with 10-8M E2. Total RNA was extracted from the three cell lines either before or after 2, 4 or 8 hrs stimulation, fluorescently labelled and analyzed on whole-genome microarrays. Samples extracted from the two ERβ-expressing cell lines were pooled before analysis, to reduce the impact of clone-specific differences and to focus on the most significant effects of ERβ, independent from tag position in the protein. Results obtained in these samples were then compared with those obtained under the same conditions in wt MCF-7 cells. This analysis yielded 921 transcripts differentially regulated by the hormone in ERβ+ vs ERβ- cells (Figure 2), including 234 mRNA whose regulation was detected only in wt MCF-7 cells, 516 regulated only in ERβ+ cells (see Venn diagram in Figure 2), 154 showing a similar pattern of response in both cell types (up- or down-regulated in all cases, although at different levels) and 17 showing opposite responses to the hormone in ERβ+ vs ERβ- cells (14 transcripts repressed in wt cells but activated in ERβ+ cells and 3 showing an opposite behaviour). The full list of these differentially regulated transcripts is reported, with relevant information, in Additional Table S1. Taken together, these results indicate that the presence of ERβ greatly influences the response of the MCF-7 cell genome to estrogen, by interfering with ERα-mediated hormonal regulation of 405 genes (Figure 2, left and central panels) and promoting de novo regulation of 516 genes (Figure 2, right panel). It should be noted that these analyses were performed with data obtained after 8 hrs of hormonal stimulation, a timing that allowed us to focus on early response genes, positioned upstream in the composite transcriptional cascade set in motion by the hormone in this cell type and more likely to include primary genomic targets of ERs [26]. It is thus possible that this analysis missed ERβ-responsive genes showing significant changes in expression only at later times after hormonal stimulation. However, analysis of the global effects of ERβ on gene expression in these same cells, performed as described above in cultures exponentially growing under continuous hormonal stimulus, suggests that the number of regulated transcripts identified here is rather close to the total number of genes targeted by this ER subtype in MCF-7 cells [O. Paris et al., manuscript in preparation].

Figure 2
figure 2

Gene expression differences in absence or presence of ERβ. Top: Heatmap summarizing the effects of ERβ expression of the estrogen responsive transcriptome of MCF-7 cells, showing changes in expression (log2 of the fold-change) of 921 transcripts after cell exposure to 10-8M E2 for the indicated times. Transcripts are grouped as follows: regulated only in wt MCF-7 cells (1), in both cell lines (2) and only in TAP-ERβ cells (3). Bottom: Venn diagram showing the numbers of differentially regulated by E2 in wt MCF-7 only (1), both cell lines (2) or ERβ expressing cells only (3).

Global mapping of ERα and ERβ binding to MCF-7 cell genome

The widespread effects of ERβ on MCF-7 cell transcriptome are likely to result from multiple effects of this receptor in the cells, including direct regulation of primary response genes via genomic or non genomic mechanisms, and indirect gene regulation events mediated by the products of primary genes. The primary ERβ target genes are most likely to comprise also master regulators of complex cellular responses to the receptor, mediating its effects on the biological and clinical phenotype of BC cells. To identify such primary genomic targets and investigate the mechanisms that allow their regulation by ERβ, a global analysis of in vivo binding of this receptor to the genome was carried out in TAP-ERβ cells by chromatin immunoprecipitation coupled to massively parallel sequencing (ChIP-Seq) [38], that allows detailed mapping of in vivo TF binding to the genome. In parallel, we studied ERα binding to the genome under the same conditions, to allow comparative analyses between the two ER subtypes. Replicate chromatin samples were prepared from both Ct-ERβ and Nt-ERβ cells before and after stimulation for 45 minutes with 10-8M E2 and DNA-bound proteins were immunoprecipitated either with antibodies against the N- and C-terminus of ERα, or with IgGs binding with high affinity the TAP moiety of tagged ERβ (see Methods). Preliminary testing on several known ERβ binding sites, including the promoter-near region of pS2/TFF1 gene [26], confirmed that the method selected to immunoprecipitate chromatin-bound Ct-ERβ and Nt-ERβ was efficient and specific (data not shown). The resulting DNAs were used to generate ChIP-Seq libraries for ERα and ERβ, respectively, that were then sequenced with the Illumina Genome Analyzer. The sequence tags obtained were then aligned to the human genome sequence and peaks enriched in the libraries generated after E2-treatment were identified using MACS (Model-based Analysis of ChIP-Seq) [39]. This led to the identification of 9702 binding sites for ERβ and 6024 sites for ERα, of which 4506 were shared by both receptors (Figure 3A), with an average False Discovery Rate (FDR) of 3%. The full list of these binding sites is available, with relevant information, in Additional Files 1 and 2. Interestingly, about half (4862) ERβ binding sites identified map within transcription units, mainly (3942 sites) in intronic regions. This distribution is maintained also in 424 ERβ-regulated transcription units (see below), where 966 ERβ binding sites located in the gene or within 10 kbps from it are distributed as follows: 154 in promoter regions, 51 in exons, 471 within one or more introns and the remaining either upstream of promoters (156) or downstream of the gene (134). In both cases the ERβ binding sites within genes did not show any preference with respect to exon or intron position nor for know intragenic regulatory elements (splice sites, polyadenylation sites, etc). It should be mentioned that the number of ERβ binding sites identified is significantly higher that those mapped in MCF-7 cells by ChIP-on-chips [30, 31], possibly for technical differences due to ERβ expression levels in the different MCF-7 cell-derived clones used, in immunoprecipitation efficiency and/or in DNA analysis. Since only Zhao et al. [31] performed an unbiased search for ERβ binding with whole-genome chips, we could confront our results only with those reported in that study. This showed that 86% of high confidence ERβ sites described in that study appear also in our dataset. The binding sites identified here were then subjected to sequence analysis, searching first for the presence of EREs (Estrogen Receptor Elements), the characteristic ER binding signature (Figure 3B). This analysis revealed that in all three cases (i.e. ERβ, ERα and ERβ+ERα) a high percentage of sites displayed one or more imperfectly palindromic ERE (ERE+), with a slightly higher positivity in ERα sites (58.89 vs 53.51%). As ERs have been shown to bind both in vitro and in vivo to PuGGTCA hemi-palindromes (hEREs), we searched the sequence of the remaining (ERE-) sites for perfect matches to this sequence. Results showed that almost half of them indeed contained one or more hEREs. The percentages of sites not carrying a known ER-binding element (ERE- and hERE-) were similar for both receptors (ERα: 22.34%, ERβ: 28.38%). We observed that ERα and ERβ binding sites were often found in close proximity to each other, a confounding factor when attempting to discern and analyze separately ER subtype-specific sites and target genes. This could be due to the limits of the ChIP-Seq technology or of the algorithm used for peak selection. To overcome this problem, and allow the identification of potential ER subtype-specific sites, we used a cartographic approach to group nearby binding sites that might be the result, at least in part, of shortfalls of the mapping methods applied. Each binding peak was thus elongated in both directions by 1000 bp and the overlapping ones obtained were merged into 8536 ERβ and 5371 ERα 'extended' binding regions. These regions were intersected to define ERα only, ERβ only or ERα+ERβ binding regions. In this way we could identify 1271 ERα-only and 4541 ERβ-only binding sites, comprised in these regions, none of which showed nearby binding of the other receptor. These were named: ER subtype 'prevalent' sites. The binding peak sequences included in each of the three regions obtained (ERα only, ERβ only or ERα+ERβ) were then re-analyzed for the presence of ERE or hERE elements. In this way we could observe that sites within the ERα+ERβ regions showed now a much higher percentage of ERE+ sequences (62.90%), respect to those present in the ERα-only or ERβ-only regions (45.63% and 44.62%, respectively, Figure 3C). Since all three types of sites showed almost identical proportions of hERE+, this result suggests that perfectly or imperfectly palindromic EREs are preferential binding sequences for ERα-ERβ heterodimers, while ERα and ERβ homodimers appear to be more flexible in DNA recognition. ERE+ sequences were then analyzed in more detail with MEME [40], to investigate if the three classes of sites identified showed any difference in the relative base composition of their respective ERE signatures. For each list of sequences, the most significant position-specific probability matrix generated by MEME was compared to the matrices present in the JASPAR transcription factor binding profile database [41], using the STAMP tool-kit for DNA motif comparison [42]. As shown in Figure 3C (left panels), this analysis revealed that the ERE matrices derived from the three types of binding regions identified (ERα selective, ERβ selective and ERα+ERβ) are identical and, as a consequence, that ERβ does not appear to display ERE variant selectivity.

Figure 3
figure 3

Sequence analysis of ERα, ERβ and ERα+ERβ binding sites. (A) Venn diagram showing a summary of ERα and ERβ binding sites identified in TAP-ERβ cells by ChIP-Seq. (B) Classification of ERα and β binding sites based on the presence of a perfectly or imperfectly palindromic Estrogen Response Element (ERE, green), an ERE hemipalindrome (hERE, blue) or no ERE (none, red). (C) ERE motif matrices identified in each of the three ER binding regions indicated (left), classification of the binding sites belonging to each region according to the presence of ERE (center) and grid summarizing the results of TFBS matrix enrichment (overrepresentation) analyses performed on the binding sites groups indicated (right). Z-Score cut-off was 3.0 and only TFBSs showing an over-representation score ≥4.0 in at least one of the ERE- (none) binding site groups. Light grey cells indicate a Z-Score <3.0 while dark grey cells indicate absence of the matrix.

We then examined the ERE- sequences to search for enriched binding motifs for other transcription factors that might play a role in ER binding to chromatin in the absence of canonical EREs. ERs are known to be able to bind chromatin indirectly, by physically interacting with DNA-bound TFs (tethering), including SP1 [43, 44] or the AP1complex [31, 45, 46], for gene trans-regulation. TFBS enrichment respect to the whole genome was calculated thus in ERE- sites with using the RegionMiner tool [47] and only statistically significant (Z-score ≥3.0) and highly enriched (over-representation ≥4-fold) matrices were further considered. The results are summarized in the right panel of Figure 3C, showing for comparison the over-representation values scored in hERE+ and ERE+ sites by the same matrices selected in the ERE- sites. These numerical values, together with relevant information, are reported in Additional Table S2 [Additional file 3]. The enriched matrices found in the ERE- set of ERβ prevalent sites include V$SP1.01, V$SP1.02, V$SP1.03, V$SP2.01 and V$GC.01, all belonging to the GC-Box family targeted by SP1 and GCFC1 (GC-rich sequence DNA-binding) factors, V$CTCF.01 and V$CTCF.02, binding site matrices for the CCCTC-binding factor CTCF, that is a known transcriptional repressor of c-myc [48], V$NRF1.01, binding NRF1 (nuclear respiratory factor 1), a transcription factor that regulates the expression of nuclear-encoded mitochondrial genes [49], V$ZF5.01, for the POZ domain zinc finger protein ZF5, and V$ZNF9.01, recognized by the zinc finger proteins ZNF148, 202, 219 and 281. The majority of these TFs bind GC- and C-rich sequences that are structurally related to each other, suggesting the possibility that a significant portion of the sequence elements listed above might indeed be recognized by one or a limited number of TFs. On the other hand, the V$GAGA.01 matrix was specifically enriched only in the ERα prevalent ERE- binding sites. This sequence binds a transcription factor known to influence chromatin structure in Drosophila[50] and to bind throughout the genome [51], but nothing is known about physical or functional interactions of this factor with ERs or other nuclear receptors. The results of this analysis point to TFs that could act as partners of ERβ for binding to chromatin in the absence of canonical EREs. Interestingly, the majority of these same matrices were found enriched also in the ERβ binding sites comprising hEREs or EREs (Figure 3C), suggesting that one or more such DNA matrices might cooperate with ERβ for either DNA binding or gene trans-regulation. We performed a direct search for conserved motifs in the ERE- binding sites of ERβ with MEME [40], and the most significant position-specific probability matrices were compared to those present in the JASPAR TF binding profile database [41]. The results failed to provide any conclusive information, as each of several sequence motifs obtained with this analysis was found only in a small fraction of the binding sites analyzed.

Identification of primary ERβ target genes

To identify the genes directly controlled by ERβ binding to the genome, and analyze the interplay between the two ERs in gene regulation, we combined the ChIP-Seq data with those relative to estrogen responsive genes differentially regulated by E2 in ERβ+ vs ERβ- MCF-7 cells under the same experimental conditions (Figure 2). Three sets of ER binding regions (defined as described in Methods) were used for this analysis: ERβ (8872, Set 1) and ERα (5558, Set 2) binding regions identified in TAP-ERβ cells, and all ERα binding regions identified so far in wt MCF-7 (17888, Set 3). The regions from Set 1 (ERβ in ERβ+ cells) were intersected with those from Set 2 (ERα in ERβ+ cells), to define which of them was binding both receptor subtypes ('heterodimer ERα+ERβ': 4186), only ERβ ('homodimer ERβ': 4686) or only ERα ('homodimer ERα': 1372) in TAP-ERβ cells. The 'homodimer ERβ' and the 'heterodimer ERα+ERβ' groups were further filtered against Set 3 binding regions (ERα sites detected in ERβ- cells), to identify the genomic sites recognized by ERβ, with or without ERα, but never by ERα alone. This allowed us to classify the sites comprised in these regions as follows: Class 1, including 2126 sites occupied by ERα in wt MCF-7 cells and by ERβ in TAP-ERβ cells (ERβ vs ERα competition); Class 2, showing 4340 sites where ERα can bind in wt MCF-7 cells and both receptors are detected in TAP-ERβ cells (ERβ+ERα); Class 3, with 2707 sites binding only ERβ and never, under any condition, ERα (ERβ specific); Class 4, comprising 529 sites where both receptors are found in TAP-ERβ cells but none in wt MCF-7 cells (ERβ+ERα specific); Class 5, including 617 sites where ERα binds only in TAP-ERβ but never in wt cells (ERα displacement); Class 6, composed of 773 sites that bind only ERα both in wt and TAP-ERβ cells (ERα specific). When combined with the results of the sequence analyses described above, this classification reveals that ERβ-specific cis-acting regulatory elements are unlikely to exist in the genome, as all evidence point to the fact that the two ER subtypes can recognize identical sequence features.

To identify among the genes differentially regulated by estrogen in ERβ+ vs ERβ- cells those representing direct targets for transcriptional regulation by DNA-bound ERβ in our cell model, we extracted from the list in Additional Table S1 [Additional file 3] the genes bearing one or more ERβ binding sites inside or within 10 kb of the TU, and termed them 'primary', to indicate that they are most likely to respond directly to the signal conveyed by the E2-activated receptor [26]. Of these 424 genes - listed in Additional Table S3 [Additional file 3], whose kinetics of response to E2 in wt and TAP-ERβ cells is shown in Figure 4, 52 show one ERβ site of Class 1 (ERβ vs ERα competition), 90 a site of Class 2 (ERβ+ERα), 71 a site of Class 3 (ERβ specific) and only 9 a Class 4 site (ERβ+ERα specific), while 202 showed multiple ERβ sites belonging to different classes and were thus classified accordingly (grey in Figure 4). In the right panels of Figure 4 are reported examples for each of the primary gene classes described above, showing the location of the receptor binding sites respect to the promoter and structural gene. It is worth mentioning that when the gene expression data from wt MCF-7 cells stimulated with E2 for 8hrs (Figure 2) were combined with information concerning ERα binding regions identified in wt MCF-7 cells under comparable conditions (Set 3 described above), 228 putative primary ERα target genes were identified -listed in Additional Table S4 [Additional file 3], 71% of which (163 genes) showed ERβ binding in hormone-stimulated ERβ+ cells. This result supports the view that the two ER subtypes tend to interact with the same targets in BC cells genome.

Figure 4
figure 4

Putative ERβ primary targets. Left: Heatmap summarizing the effects of estrogen stimulation on 424 mRNAs encoded by genes showing one (left) or more (right) ERβ binding sites within 10 kb of the TU (primary response genes) transcriptome of MCF-7 cells, as changes in expression (log2 of the fold-change) after cell exposure to 10-8M E2 for the indicated times. Transcripts are grouped as follows: regulated only in wt MCF-7 cells (1), in both cell lines (2) and only in TAP-ERβ cells (3). Vertical bars to the right of each heatmap indicate the class of ERβ binding site present, as indicated in the legend. When a regulated gene showed multiple ERβ binding sites belonging to different classes it was included in a separate group, classified as 'combination of ERβ sites' (grey bar).Right: Genome Browser view of genomic loci representative of the different ERα and ERβ binding site categories identified. ChIP-Seq ERα and ChIP-Seq ERβ indicate sites identified in this study, ChIP-Seq ERα1, ChIP-Seq ERα2 and ChIP-on-chip indicates sites identified in MCF-7 cells by Cicatiello et al. [26], Fullwood et al. [27], and Hurtado et al. [28], respectively.

A functional analysis of the primary ERβ target genes identified here with Ontologizer [52] showed that most primary ERβ responsive genes are involved in key cellular processes, including control of cell proliferation, survival and differentiation status as well as cell motility, migration and adhesion, and can all greatly influence BC cell phenotype and response to estrogen - listed in Additional Table S5 [Additional file 4]. When GO analysis was performed on the ERα target gene set from wt MCF-7 cells (228 genes), results showed that the genes controlled directly by this ER subtype appear to be involved in the same cellular processes described above for ERβ -compare results reported in Additional Tables S5 and S6 [Additional files 4 and 5], providing a further indication of the significant overlapping between gene pathways targeted by ERβ and ERα in BC cells. Focusing on the genes differentially regulated by E2 in ERβ- vs ERβ+ cells known for their involvement in cell proliferation, we observed that many of them encode for transcription factors and other key proteins controlling large gene networks of cell division cycle and cell survival and, in general, cell proliferation and differentiation pathways. These include, in particular, CDK-6 (cyclin-dependent kinase 6), CEBPA (CCAAT/enhancer binding protein, alpha), DAB2 (Disabled homolog 2, mitogen-responsive phosphoprotein), HES-1 (Hairy and enhancer of split homologue 1), IGFBP-4 (Insulin-like growth factor binding protein 4), IRS-1 and -2 (Insulin receptor substrates 1 and 2), JAK-2 (Janus kinase 2), JunB, MITF (Microphthalmia-associated transcription factor), MYC, SLIT-2 (Slit homolog 2, Drosophila), SMARCA-2 (SWI/SNF related, matrix associated, actin dependent regulator of chromatin, subfamily a, member 2), SOX-9 (Sex determining region Y-box 9), TGFB-2 (Transforming growth factor beta 2), TGFBI/LCD-1 (transforming growth factor, beta-induced, 68 kDa), TGM-2 (Transglutaminase 2), TNS-3 (Tensin 3) and WISP-2 (WNT1 inducible signaling pathway protein 2). Interestingly, the role of all these genes in tumor cell proliferation and differentiation is known and an involvement in hormone-mediated BC cell responses to ERα has been reported for most of them, suggesting that discovery of an ERβ direct effect on transcription of these genes provides a new molecular framework to elucidate the anti-proliferative and differentiative effects of this receptor subtype in hormone-responsive cells.

Among the RNAs encoded in the genome, microRNAs (miRNAs) have emerged as master regulators of gene expression for their ability to influence mRNA concentration and activity by post-transcriptional mechanisms. Recent results highlighted the role of miRNA in BC cells response to estrogen [26, 53–58] and, in addition, several lines of evidence indicate extensive miRNA deregulation in BC, including differential expression of miRNAs in normal vs transformed mammary epithelial cells [59–61]. For these reasons, we focused our attention on the TUs encoding pre-miRNAs, to test the possibility that ERβ binding sites might be located in close proximity of these genes. Results show that 52 miRNA-encoding loci (isolated or in genomic clusters) show one or more ERβ site within 10 kb from the pre-miR sequence of the host gene - listed in Additional Table S7 [Additional file 6]. Distribution of these sites among the ER binding Classes described above was comparable to what observed for primary genes. Interestingly, in several cases ERβ binds both up- and down-stream of the pre-miR locus, further suggesting that receptor docking might exert multiple regulatory actions on miRNA biogenesis.

We tested the hypothesis that the observed distribution of sites near the pre-miRNA loci occurred at random by applying a bootstrap approach. To this end, we repeatedly sampled 1000 times the same number of sites of the real ERβ binding set, with the same length distribution, a similar distribution among chromosomes but randomly selected coordinates. We then counted the number of pre-miRNA loci and the number of randomly generated sites found within 10 kb of each other and compared their distributions with that of the experimentally detected ones. The number of randomly generated sites within 10 kb of a pre-miRNA never reached the value detected experimentally, while the number of miRNA loci scoring an artificial site in close proximity was equal or above what measured only in 7.6% of the cases. These results can be explained also by the observation that in several cases ERβ binds both up- and down-stream of the pre-miR locus, a result that support a functional significance of this observation. In fact, the ratio between the number of ERβ binding sites and the number of pre-miRNA loci within 10 kb of each other is 1.5, while this varied between 0.5 and 1 for the randomly generated sets (data not shown). We thus conclude that although some of the ERβ sites detected in or near pre-miRNA loci may be the result of a random, non functional event, these is likely to represent rare events, as random distribution never reaches the enrichment level observed experimentally. Indeed, preliminary miRNA profiling analyses carried out in wt MCF-7 and TAP-ERβ cells indicate that mature miRNA expression undergo extensive deregulation in the presence of ERβ [O. Paris et al., manuscript in preparation and data not shown], suggesting that at least some of the sites identified here might indeed be involved in ERβ-mediated regulation of miRNA biosynthesis in BC cells.

ERβ binding to the mitochondrial genome

Mitochondrial DNA (mtDNA) is usually overlooked in whole-genome ChIP-seq analyses, since identification of enriched peaks is more difficult here due to a much higher noise, consequence of the high and variable number of mtDNA copies in the cell. ERβ has been shown to localize in the mitochondria in different cell types [62] including human BC cells [63, 64], and a role for estrogen in mitochondrial function, with implications on cell growth, has been established [65, 66]. We thus analyzed the sequence reads that aligned on the mtDNA sequence with the same method used for whole genome data analysis, but applying a supplementary fold-intensity filter (see Methods) to deal with the higher background noise. This analysis revealed one ERβ binding site in proximity of the mtDNA D-loop, but no ERα binding sites either in wt MCF-7 or in TAP-ERβ cells (Figure 5A). Independent ChIP analysis confirmed this results, showing ERβ binding to this same mtDNA region upon activation by E2 or its selective agonist 2,3-bis(4-hydroxyphenyl) propionitrile (DPN) and lack of ERα binding (Figure 5B). Furthermore, the presence of ERβ in mitochondria of TAP-ERβ cells was confirmed biochemically, by western blotting (Figure 5C), an analysis that revealed also the presence of ERα in the organelle. When we analyzed the sequence of the ERβ mitochondrial site with MatInspector, we observed the presence of the matrix V$GATA1.06, bound by GATA1, a factor whose activity is strongly repressed by ERα [67], and V$HMGA.01, bound by HMGA1, a non-histone chromosomal protein that is highly overexpressed in cancer cells [68] and has been shown to interact with ERα and to enhance its binding to DNA [69].

Figure 5
figure 5

Mitochondrial ER-beta binding sites. (A) Genome Browser view of the ERβ binding site identified in mitochondrial genome. (B) Validation of ERβ binding to mitochondrial DNA by ChIP. Results shown are representative of duplicate analyses. E2: 10-8M 17β-estradiol; PPT: 10-8M 1,3,5-tris(4-hydroxyphenyl)-4-propyl-1H-pyrazole (selective ERα agonist); DPN: 10-8M 2,3-bis(4-hydroxyphenyl) propionitrile (ERβ agonist). (C) Western blot analysis of ERβ and/or ERα in purified mitochondria from Ct-ERβ and Ct-ERα [33] cells. Cyt: cytosol depleted of mitochondria, H: whole cell homogenate; Mito: purified mitochondrial fraction. ATPase is a mitochondrial marker and α-tubulin was included to determine the level cytosolic contaminants in 'Mito' samples.

Conclusions

The results of this study indicate that in vivo ERβ can interact with hormone-responsive BC cell chromatin either alone or complexed with ERα, but in all cases the two receptors share the same genomic targets. An observation that is in agreement with the conclusions of previous studies based on analysis of ERα and ERβ heterodimerization and binding to the ERE [9, 70–72]. When both ER subtypes are expressed in the same cell, the main action of ERβ in the genome is thus achieved in combination with ERα, by either heterodimerization or competition for binding to the same target sites in chromatin. Based on these observations, we propose that in hormone-responsive BC the final cellular response to estrogen is likely to depend upon the relative concentration of the two ERs in the cell, their activation status and DNA binding kinetics and the presence of other factors influencing their respective functions.

Methods

Plasmid preparation, cell lines, cell culture and stable transfections

Different BC cell lines were used: MCF-7 TET Off (ER-alpha positive; ATCC, Cat No. HTB-22) and SKBR3 (ER-alpha negative; ATCC, Cat No. HTB-30). MCF-7 TET Off cells (described here as wt or ERβ-) were used to produce stable clones expressing ERβ tagged with TAP-tag respectively at the C-term and at the N-term (Ct-ERβ and Nt-ERβ) or C-tagged ERα (Ct-ERα), as previously described [5, 33]. All were grown in Dulbecco's modified Eagle's medium (DMEM), supplemented with 10% fetal bovine serum (FBS) (both from Invitrogen), 100 U/ml penicillin, 100 mg/ml streptomycin, 250 ng/ml Amfotericin-B, 50 μg/ml G418 (normal growing condition). For estrogen starvation, cells were plated at 40% confluence in steroid-free medium (phenol red-free DMEM medium, with 5% fetal bovine serum pre-treated with dextran-coated charcoal and antibiotics) and maintained for 5 days with replacement of fresh medium before stimulation with 10-8M 17β-estradiol (E2).

Preparation of cell extracts, mitochondria isolation and immunoblotting analyses

Cells starved in 100-mm dishes were stimulated for 45 minutes, harvested in cold PBS and collected by centrifugation at 1000 × g. The cell pellets were then resuspended in three volumes of Hypotonic Buffer (HB) (20 mM HEPES pH 7.4, 5 mM Sodium Fluoride, 10 μM Sodium Molybdate, 0.1 mM EDTA, 1 mM dithiothreitol, 1 mM protease inhibitors, 1 mM Phenylmethyl-Sulfonyl Fluoride). Cells were then incubated on ice for 15 min and 0.5% Nonidet P-40 followed by spinning 30 sec at 4°C at 13000 × g. Supernatants were recovered and clarified at 13,000 × g for 15 min at 4°C while pellets were resuspended in hypotonic buffer, stratified on 25% sucrose-HB solution and centrifuged at 6000 × g for 15 min at 4°C. The resulting pellets were then resuspended in one volume of Nuclear Lysis Buffer [73] containing 800 mM NaCl, incubated for 30 min at 4°C with gentle shaking and centrifuged for 15 min at 4°C at 13000 × g. The supernatant fraction was recovered.

Mitochondria were isolated from 20 × 106 Ct-ERβ or Ct-ERα cells (in 150 mm culture dishes) as described [74], with minor modifications. All steps during mitochondria isolation were performed at 4°C, cells were washed twice in PBS, scraped, and centrifuged at 290 × g, 5 minutes. The samples were resuspended in buffer A (250 mM Sucrose, 50 mM Tris-HCl, 2 mM EGTA), homogenized in Glass-Teflon Potter homogenizer and centrifuged at 600 × g for 3 min, then the supernatants were re-centrifuged at the same speed. Mitochondria were pelleted by centrifugation at 10400 × g for 10 min, resuspended in buffer A and centrifuged again at 5300 × g for 10 min, in order to eliminate microsomal and cytosolic contamination. The samples were dissolved in buffer A and centrifuged at 1500 × g for 4 min and then pelleted again at 9600 × g for 10 minutes. The final pellet was resuspended in Buffer B (50 mM Tris pH 7.5, 150 mM NaCl, 1% Triton X-100, 0.1% SDS, 1% Na Deoxycholate, 1 mM PMSF, 1X Protease Inhibitor) and incubated on ice for 20 min to extract mitochondrial proteins. A small portion of sample was collected after homogenization and processed to obtain the samples corresponding to homogenate and cytosol. An equivalent protein amount was fractionated by SDS-PAGE on Mini Protean Precast polyacrylamide gels 4-20% from Biorad.

Homogenate, cytosol, nuclear or mitochondrial extracts from equivalent cell number were fractionated by SDS-PAGE. Immunoblot analysis was performed using the following primary antibodies: Ct-ERα (HC-20; sc-543) from Santa Cruz Biotechnology, Inc., TAP tag (CAB1001) from Open Biosystems, α-Tubulin (T 6199) from Sigma-Aldrich, ATPase B (ab14730) from Abcam. Peroxidase-labelled anti-rabbit or -mouse immunoglobulin antisera were used according to the manufacturer's instructions (Amersham Italia).

Transient Transfections and Luciferase Assay

Wild type MCF-7, Ct-ERβ and Nt-ERβ clones were starved for 5 days in estrogen-free medium. Then 5 × 105 cells/dish were seeded in 60 mm culture dishes and transfected by using 25 μg/dish polyethylenimine (Polysciences, Inc.) with 2.5 μg/dish DNA, including 300 ng ERE-tk-Luc [75], 500 ng pSG-Δ2-NLS-LacZ vector [76], co-transfected as an internal control for transfer efficiency, and carrier DNA (Bluescribe M13+). At 48 hrs after transfection, cells were treated for 24 hrs with either vehicle (EtOH) or E2 (10-8M).

SKBR3 cells were grown to 60-70% confluence and transfected with Lipofectamine 2000 reagent (Invitrogen) and OPTI-MEM (Invitrogen) according to the manufacturer's instructions. The plasmids used were pSG5-ERβ, encoding full-length ERβ (ERβ1), pSG-HEGO, encoding the full-length ERα and the corresponding pSG5 empty vector (Stratagene), pUSE-C-TAP-ERβ and pUSE-N-TAP-ERβ, encoding full-length ERβ tagged, respectively, at the C-term and at the N-term [5], ERE-tk-Luc and pSG-Δ2-NLS-LacZ. Six hours after transfection, the medium was changed and 24 hrs later cells were stimulated as described above for 24 hrs. At the end of treatment, cells were washed with cold PBS and lysed in 100 μl lysis buffer (Promega). Luciferase activity was measured in extracts using the Luciferase Assay Reagent (Promega), according to the manufacturer's instructions, and values were expressed as relative light units normalized to the β-galactosidase activity or to the protein concentrations measured using the Bradford technique. For each condition, average luciferase activity was calculated from the data obtained from three independent dishes.

Cell Cycle Analysis

Estrogen-deprivation was always controlled by cell cycle analysis as follows. Cells (1.5 × 105 cells/dish) were starved in 60 mm culture dishes, stimulated for 27 hrs and collected in PBS containing 50 μg/ml propidium iodide, 0.1% (v/v) sodium citrate, 0.1% (v/v) Nonidet P-40. Cell samples were incubated in the dark for at least 15 min at room temperature, or overnight at 4°C, and analyzed by a FACScalibur flow cytometer using the CellQuest software package (BD Biosciences), according to standard protocols suggested by the manufacturer [77, 78]. Data analysis was performed with Modfit software (Verity Software, Topsham). Values were plotted as increasing of S-phase respect to unstimulated controls. Results showed were obtained from two independent experiments.

Cell Proliferation Assay

Hormone-starved cells (3000/well) were seeded in 96-well dishes. After 12 hrs, medium was changed to include the indicated compounds. After appropriate stimulation, cells were washed in phosphate-buffered saline (PBS) and fixed with 12.5% glutaraldehyde for 20 min at room temperature, followed by washing with distilled water, incubation with 0.05% methylene blue for 30 min, rinsing and incubation with 0.33 M HCl for 18 hrs. Absorption was measured at 620 nm.

RNA purification

Total RNA was extracted from wt MCF-7, Ct-ERβ and Nt-ERβ clones, generated as described above, using the standard RNA Extraction method with TRIzol (Invitrogen) method, as described previously [79, 80]. Cells were estrogen-deprived and total RNA was extracted before or at the indicated times after stimulation with 10-8M 17β-estradiol (+E2) or ethanol vehicle.

In each case RNA derived from two independent experiments performed in duplicate was used. Before use, RNA concentration in each sample was assayed with a ND-1000 spectrophotometer (NanoDrop) and its quality assessed with the Agilent 2100 Bioanalyzer with Agilent RNA 6000 nanokit (Agilent Technologies).

Microarray analyses

Total RNA extracted from Ct-ERβ and Nt-ERβ cells were pooled. For mRNA expression profiling, 500 ng total RNA were reverse transcribed, as described previously [26, 81] and used for synthesis of cDNA and biotinylated cRNA according to the Illumina TotalPrep RNA Amplification Kit (Ambion, Cat. n. IL1791) protocol. For each sample, 750 ng of cRNA were hybridized for 18 hrs at 55°C on Illumina HumanHT-12 v3.0 BeadChips, containing 48,804 probes (Illumina Inc.), according to the manufacturer's protocol and subsequently scanned with the Illumina BeadArray Reader 500. Data analyses were performed with GenomeStudio software version 2009 (Illumina Inc.), by comparing all values obtained at each time point against the 0 hrs values. Data was normalized with the quantile normalization algorithm, and genes were considered as detected if the detection p-value was lower than 0.01. Statistical significance was calculated with the Illumina DiffScore, a proprietary algorithm that uses the bead standard deviation to build an error model. Only genes with a DiffScore ≤-40 and ≥40, corresponding to a p-value of 0.0001, were considered as statistical significant.

Chromatin Immunoprecipitation

Ct-ERβ and Nt-ERβ cells were hormone-deprived for 4 days and chromatin was extracted in several replicates before (-E2) and after stimulation for 45 minutes with 10-8M 17β-estradiol (+E2) or, where indicated, with the 10-8M selective ERα agonist 1,3,5-tris(4-hydroxyphenyl)-4-propyl-1H-pyrazole (PPT) or 10-8M ERβ agonist 2,3-bis(4-hydroxyphenyl) propionitrile (DPN), from Tocris Cookson.

Chromatin was prepared with the Millipore/Upstate Chromatin Immunoprecipitation (ChIP) Assay Kit (Millipore) according to the instruction provided by the producer, using a variation of the protocol described at the Upstate website. For each assay, a total of 5 × 106 cells were fixed with 1% formaldeyde for 10 min at room temperature, the reaction was then stopped by adding glycine at final concentration of 0.125 M. Fixed cells were washed twice with ice-cold PBS, harvested by scraping, centrifuged and the cell pellets were re-suspended in SDS lysis buffer. Samples were sonicated with a Diagenode Bioruptor (Diagenode) for 12 cycles of 30 sec at high power, centrifuged at 12500 xg for 15 minutes and diluted 8-fold in ChIP dilution buffer. After removing an aliquot (whole-cell extract input), the chromatin sample was divided in three aliquots, that were incubated at 4°C overnight with antibodies against either the C-term (HC-20, from Santa Cruz Biotechnology) or N-term (anti-Estrogen Receptor 18-32, from SigmaAldrich) of human ERα or with IgG Sepharose 6 fast Flow (GE Healthcare Bio-Science AB) for TAP-ERβ [5]. As control, an aliquot of the same chromatins were processed in the same way but Abs or IgGs were omitted form the reaction. The samples were then precipitated by binding to protein-A Agarose/Salmon Sperm DNA beads (Millipore), for ERα, or as such for to ERβ. The beads were washed sequentially with 'low-salt immune complex wash buffer, 'high salt immune complex wash buffer, 'LiCl immune complex wash buffer' and TE buffer, before elution in Elution buffer by ON incubation at 65°C and treatment with Proteinase K. DNA was purified from immunoprecipitated (IPP) chromatin by extraction with phenol:chloroform:isoamyl alcohol (25:24:1) and ethanol precipitation according to standard procedures. DNA pellets were dissolved in nuclease-free water and kept frozen before further use.

Primers for ChIP-QPCR validation of the mitochondrial genome ERβ binding site were the following: 5'-GATCACAGGTCTATCACCCTATTAACC (forward) and 5'-CAGCGTCTCGCAATGCTATC (reverse).

Samples preparation for ChIP-Seq

DNAs from Ct-ERβ and Nt-ERβ cells treated with E2 were pooled together to generate an ERβ +E2 sample and the same was done for DNAs from hormone-starved cells (ERβ -E2 sample). Similarly, IPP DNAs obtained with anti-C-term and anti-N-term ERα Abs were pooled together to generate ERα +E2 and ERα -E2 samples. About 20 ng of ChIP DNA was purified using the MinElute PCR Purification Kit (QIAGEN, Italy), with a recovery of 55-70%, as assessed with the Quant-IT DNA Assay Kit-High Sensitivity and a Qubit Fluorometer (Invitrogen). Preparation of IPP DNA libraries for massively parallel sequencing was performed from 10 ng purified DNA according to the Illumina ChIP-Seq sample preparation kit protocol (Illumina Inc.). Libraries were sequenced with the Illumina Cluster Station and Genome Analyzer II according to manufacturer's instructions.

ChIP-Seq data analysis

The sequence tags generated by massively parallel sequencing were aligned on the human genome (hg18) with the software ELAND, allowing up to 2 mismatches. Enriched regions from ERα +E2 and ERβ +E2 samples were compared with the same from ERα -E2 and ERβ -E2 samples, respectively. The enriched ChIP-Seq peaks were identified using MACS (Model-based Analysis of ChIP-Seq) version 1.3.7.1 [39], with standard parameters (p-value cut-off of 1e-5, mfold of 32). For mtDNA we further filtered out sites with tag density below 0.5 (N/l; N = number of tags, l = length of site)

Computational searches for ERE sequences

The ERE binding motif was searched in binding sites with the MatInspector application [82], a part of GenomatixSuite software (Genomatix Software GmbH, Germany). The matrices ER.01, ER.02, ER.03 and ER.04 (Genomatix Matrix Library 8.2), were used with a core similarity threshold of 0.75 and a matrix similarity threshold of Optimal -0.02. The sequences bearing a match of any of the four matrices were termed ERE+ sequences. On the remaining sequences, the hemi-palindromic ERE motif was searched with the same application, by defining a custom half ERE matrix PuGGTCA (hERE). The remaining sequences were classified as ERE- and were scanned for other TF binding sites motifs contained in the Genomatix Matrix Library using the standard parameters, as described previously [26].

TFBS over-representation analysis

ERE- sequences were analyzed to search for known TF binding sites (TFBSs) that were enriched (over-represented) with respect to their relative frequency in the whole human genome. This analysis was performed with the RegionMiner application of the Genomatix software suite [47]. The software automatically searches for all TFBS matches present in the submitted sequence list and calculates the over-representation value of the actual number of matches over the expected value based on its frequency in the reference set (genome or promoters) for each matrix. It reports also the significance of the over-representation ratio, expressed as Z-scores [83]. Enrichment values with a Z-score <3.0 were not considered further. A filter was applied also on the over-representation values, depending upon the range of values set for each analysis, to highlight only the stronger associations. Results are shown as heatmaps representing over-representation values, generated with MeV software [84].

Classification of ERβ binding sites and identification of primary ERβ responsive genes by combining ChIP-Seq and expression profiling data

Three sets of ER binding sites were taken in consideration. The first and second comprised, respectively, the ERβ and ERα sites mapped in TAP-ERβ cells and the third (named 'MCF-7 ERα') included all ERα binding sites identified in wt MCF-7 cells by ChIP-Seq [26, 27] and/or by ChIP-on-chip [28]. First of all, the ER binding sites from ChIP-Seq analyses were elongated in both directions by 1000 bp. Subsequently, using UCSC Table Browser [85], for each of the first two sets the extended sites overlapping with each other were merged in ERβ or ERα binding 'regions', respectively. For the third set (wt MCF-7 ERα), the extended ChIP-Seq sites and the ChIP-on-chip sites overlapping with each other in the genome were all merged to generate unique MCF-7 ERα binding 'regions'. To identify putative primary ERβ responsive genes, the TUs corresponding to RNAs differentially regulated by E2 in ERβ+ vs ERβ- cells that showed one or more ERβ binding region inside or within 10 kb were extracted using UCSC Table Browser, as described previously [26].

Gene Ontology analysis

To identify Biological Process GO terms statistically overrepresented in our regulated gene lists, we used Ontologizer 2.0, a tool for GO term enrichment analysis of genes derived from an experiment [52]. We identified enriched terms in primary ERβ or ERα target genes against all genes expressed (detected) in the cell lines under study, set as background of the analysis, with a p value threshold of 0.01.

Microarray and ChIP-Seq data accession numbers

The microarray and ChIP-Seq data have been deposited in the Array Express database ( HYPERLINK "http://www.ebi.ac.uk/arrayexpress" http://www.ebi.ac.uk/arrayexpress) with Accession Numbers E-TABM-1051 and E-MTAB-345, respectively.

References

  1. Heldring N, Pike A, Andersson S, Matthews J, Cheng G, Hartman J, Tujague M, Strom A, Treuter E, Warner M, Gustafsson JA: Estrogen receptors: how do they signal and what are their targets. Physiol Rev. 2007, 87: 905-931. 10.1152/physrev.00026.2006.

    Article  CAS  PubMed  Google Scholar 

  2. Hall JM, Couse JF, Korach KS: The multifaceted mechanisms of estradiol and estrogen receptor signaling. J Biol Chem. 2001, 276: 36869-36872. 10.1074/jbc.R100029200.

    Article  CAS  PubMed  Google Scholar 

  3. Bai Z, Gust R: Breast cancer, estrogen receptor and ligands. Arch Pharm (Weinheim). 2009, 342: 133-149. 10.1002/ardp.200800174.

    Article  CAS  Google Scholar 

  4. Ordonez-Moran P, Munoz A: Nuclear receptors: genomic and non-genomic effects converge. Cell Cycle. 2009, 8: 1675-1680.

    Article  CAS  PubMed  Google Scholar 

  5. Nassa G, Tarallo R, Ambrosino C, Bamundo A, Ferraro L, Paris O, Ravo M, Guzzi PH, Cannataro M, Baumann M, Nyman TA, Nola E, Weisz A: A large set of estrogen receptor beta-interacting proteins identified by tandem affinity purification in hormoneresponsive human breast cancer cell nuclei. Proteomics.

  6. Tarallo R, Bamundo A, Nassa G, Nola E, Paris O, Ambrosino C, Facchiano A, Baumann M, Nyman TA, Weisz A: Identification of proteins associated with ligand-activated estrogen receptor alpha in human breast cancer cell nuclei by tandem affinity purification and nanoLC-MS/MS. Proteomics.

  7. Matthews J, Gustafsson JA: Estrogen signaling: a subtle balance between ER alpha and ER beta. Mol Interv. 2003, 3: 281-292. 10.1124/mi.3.5.281.

    Article  CAS  PubMed  Google Scholar 

  8. 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.

    Article  CAS  PubMed  Google Scholar 

  9. Cowley SM, Hoare S, Mosselman S, Parker MG: Estrogen receptors alpha and beta form heterodimers on DNA. J Biol Chem. 1997, 272: 19858-19862. 10.1074/jbc.272.32.19858.

    Article  CAS  PubMed  Google Scholar 

  10. Zhao C, Matthews J, Tujague M, Wan J, Strom A, Toresson G, Lam EW, Cheng G, Gustafsson JA, Dahlman-Wright K: Estrogen receptor beta2 negatively regulates the transactivation of estrogen receptor alpha in human breast cancer cells. Cancer Res. 2007, 67: 3955-3962. 10.1158/0008-5472.CAN-06-3505.

    Article  CAS  PubMed  Google Scholar 

  11. Nassa G, Tarallo R, Guzzi PH, Ferraro L, Cirillo F, Ravo M, Nola E, Baumann M, Nyman TA, Cannataro M, Ambrosino C, Weisz A: Comparative analysis of nuclear estrogen receptor alpha and beta interactomes in breast cancer cells. Mol BioSyst. 2010,

    Google Scholar 

  12. Matthews J, Wihlén B, Tujague M, Wan J, Ström A, Gustafsson JA: Estrogen receptor (ER) beta modulates ERalpha-mediated transcriptional activation by altering the recruitment of c-Fos and c-Jun to estrogen-responsive promoters. Mol Endocrinol. 2006, 20: 534-43. 10.1210/me.2005-0140.

    Article  CAS  PubMed  Google Scholar 

  13. Katzenellenbogen BS, Montano MM, Ediger TR, Sun J, Ekena K, Lazennec G, Martini PG, McInerney EM, Delage-Mourroux R, Weis K, Katzenellenbogen JA: Estrogen receptors: selective ligands, partners, and distinctive pharmacology. Recent Prog Horm Res. 2000, 55: 163-93.

    CAS  PubMed  Google Scholar 

  14. Pettersson K, Delaunay F, Gustafsson JA: Estrogen receptor beta acts as a dominant regulator of estrogen signaling. Oncogene. 2000, 19: 4970-4978. 10.1038/sj.onc.1203828.

    Article  CAS  PubMed  Google Scholar 

  15. Sugiura H, Toyama T, Hara Y, Zhang Z, Kobayashi S, Fujii Y, Iwase H, Yamashita H: Expression of estrogen receptor beta wild-type and its variant ERbetacx/beta2 is correlated with better prognosis in breast cancer. Jpn J Clin Oncol. 2007, 37: 820-828. 10.1093/jjco/hym114.

    Article  PubMed  Google Scholar 

  16. Saji S, Hirose M, Toi M: Clinical significance of estrogen receptor β in breast cancer. Cancer Chemother Pharmacol. 2005, 56 (Suppl 1): 21-26. 10.1007/s00280-005-0107-3.

    Article  PubMed  Google Scholar 

  17. Lindberg MK: Estrogen Receptor (ER)-beta Reduces ERalpha-Regulated Gene Transcription, Supporting a "Ying Yang" Relationship between ERalpha and ERbeta in Mice. Mol Endocrinol. 2003, 17: 203-208. 10.1210/me.2002-0206.

    Article  CAS  PubMed  Google Scholar 

  18. Lazennec G, Bresson D, Lucas A, Chauveau C, Vignon F: ER beta Inhibits Proliferation and Invasion of Breast Cancer Cells. Endocrinology. 2001, 142: 4120-4130. 10.1210/en.142.9.4120.

    CAS  PubMed  PubMed Central  Google Scholar 

  19. Behrens D, Gill JH, Fichtner I: Loss of tumourigenicity of stably ERbeta-transfected MCF-7 breast cancer cells. Mol Cell Endocrinol. 2007, 274: 19-29. 10.1016/j.mce.2007.05.012.

    Article  CAS  PubMed  Google Scholar 

  20. Paruthiyil S, Parmar H, Kerekatte V, Cunha GR, Firestone GL, Leitman DC: Estrogen Receptor Inhibits Human Breast Cancer Cell Proliferation and Tumor Formation by Causing a G2 Cell Cycle Arrest. Cancer Res. 2004, 64: 423-428. 10.1158/0008-5472.CAN-03-2446.

    Article  CAS  PubMed  Google Scholar 

  21. Chang EC, Frasor J, Komm B, Katzenellenbogen BS: Impact of estrogen receptor beta on gene networks regulated by estrogen receptor alpha in breast cancer cells. Endocrinology. 2006, 147: 4831-4842. 10.1210/en.2006-0563.

    Article  CAS  PubMed  Google Scholar 

  22. Lin CY, Ström A, Li Kong S, Kietz S, Thomsen JS, Tee JB, Vega VB, Miller LD, Smeds J, Bergh J, Gustafsson JA, Liu ET: Inhibitory effects of estrogen receptor beta on specific hormone-responsive gene expression and association with disease outcome in primary breast cancer. Breast Cancer Res. 2007, 9: R25-10.1186/bcr1667.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Williams C, Edvardsson K, Lewandowski SA, Ström A, Gustafsson J: A genome-wide study of the repressive effects of estrogen receptor beta on estrogen receptor alpha signaling in breast cancer cells. Oncogene. 2008, 27: 1019-1032. 10.1038/sj.onc.1210712.

    Article  CAS  PubMed  Google Scholar 

  24. Fritzemeier KH, Hillisch A, Elger W, Kaufmann U, Kollenkirchen U, Kosemund D, Lindenthal B, Müller G, Muhn P, Nubbemeyer R, Peters O, Siebel P, Hegele-Hartung C: Biological effects of ERalpha- and ERbeta-selective estrogens. Ernst Schering Res Found Workshop. 2004, 46: 127-150.

    CAS  Google Scholar 

  25. Carroll JS, Liu XS, Brodsky AS, Li W, Meyer CA, Szary AJ, Eeckhoute J, Shao W, Hestermann EV, Geistlinger TR, Fox EA, Silver PA, Brown M: Chromosome-wide mapping of estrogen receptor binding reveals long-range regulation requiring the forkhead protein FoxA1. Cell. 2005, 122: 33-43. 10.1016/j.cell.2005.05.008.

    Article  CAS  PubMed  Google Scholar 

  26. Cicatiello L, Mutarelli M, Grober OM, Paris O, Ferraro L, Ravo M, Tarallo R, Luo S, Schroth GP, Seifert M, Zinser C, Chiusano ML, Traini A, De Bortoli M, Weisz A: Estrogen Receptor alpha Controls a Gene Network in Luminal-Like Breast Cancer Cells Comprising Multiple Transcription Factors and MicroRNAs. Am J Pathol. 2010, 176: 2113-2130. 10.2353/ajpath.2010.090837.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Fullwood MJ, Liu MH, Pan YF, Liu J, Xu H, Mohamed YB, Orlov YL, Velkov S, Ho A, Mei PH, Chew EG, Huang PY, Welboren WJ, Han Y, Ooi HS, Ariyaratne PN, Vega VB, Luo Y, Tan PY, Choy PY, Wansa KD, Zhao B, Lim KS, Leow SC, Yow JS, Joseph R, Li H, Desai KV, Thomsen JS, Lee YK, Karuturi RK, Herve T, Bourque G, Stunnenberg HG, Ruan X, Cacheux-Rataboul V, Sung WK, Liu ET, Wei CL, Cheung E, Ruan Y: An oestrogen-receptor-alpha-bound human chromatin interactome. Nature. 2009, 462: 58-64. 10.1038/nature08497.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Hurtado A, Holmes KA, Geistlinger TR, Hutcheson IR, Nicholson RI, Brown M, Jiang J, Howat WJ, Ali S, Carroll JS: Regulation of ERBB2 by oestrogen receptor-PAX2 determines response to tamoxifen. Nature. 2008, 456: 663-666. 10.1038/nature07483.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Vivar OI, Zhao X, Saunier EF, Griffin C, Mayba OS, Tagliaferri M, Cohen I, Speed TP, Leitman DC: Estrogen receptor beta binds to and regulates three distinct classes of target genes. J Biol Chem. 2010, 285: 22059-22066. 10.1074/jbc.M110.114116.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Charn TH, Liu ET, Chang EC, Lee YK, Katzenellenbogen JA, Katzenellenbogen BS: Genome-wide dynamics of chromatin binding of estrogen receptors alpha and beta: mutual restriction and competitive site selection. Mol Endocrinol. 2010, 24: 47-59. 10.1210/me.2009-0252.

    Article  CAS  PubMed  Google Scholar 

  31. Zhao C, Gao H, Liu Y, Papoutsi Z, Jaffrey S, Gustafsson JA, Dahlman-Wright K: Genome-wide mapping of estrogen receptor-beta-binding regions reveals extensive cross-talk with transcription factor activator protein-1. Cancer Res. 2010, 70: 5174-83. 10.1158/0008-5472.CAN-09-4407.

    Article  CAS  PubMed  Google Scholar 

  32. Liu Y, Gao H, Marstrand TT, Ström A, Valen E, Sandelin A, Gustafsson JA, Dahlman-Wright K: The genome landscape of ERalpha- and ERbeta-binding DNA regions. Proc Natl Acad Sci USA. 2008, 105: 2604-2609. 10.1073/pnas.0712085105.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Ambrosino C, Tarallo R, Bamundo A, Cuomo D, Franci G, Nassa G, Paris O, Ravo M, Giovane A, Zambrano N, Lepikhova T, Jänne OA, Baumann M, Nyman TA, Cicatiello L, Weisz A: Identification of a hormone-regulated dynamic nuclear actin network associated with estrogen receptor alpha in human breast cancer cell nuclei. Mol Cell Proteomics. 2010, 9: 1352-1367. 10.1074/mcp.M900519-MCP200.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Cicatiello L, Addeo R, Sasso A, Altucci L, Petrizzi VB, Borgo R, Cancemi M, Caporali S, Caristi S, Scafoglio C, Teti D, Bresciani F, Perillo B, Weisz A: Estrogens and progesterone promote persistent CCND1 gene activation during G1 by inducing transcriptional derepression via c-Jun/c-Fos/estrogen receptor (progesterone receptor) complex assembly to a distal regulatory element and recruitment of cyclin D1 to its own gene promoter. Mol Cell Biol. 2004, 24: 7260-7274. 10.1128/MCB.24.16.7260-7274.2004.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. 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.

    Article  CAS  PubMed  Google Scholar 

  36. Frasor J, Danes JM, Komm B, Chang KC, Lyttle CR, Katzenellenbogen BS: Profiling of estrogen up- and down-regulated gene expression in human breast cancer cells: insights into gene networks and pathways underlying estrogenic control of proliferation and cell phenotype. Endocrinology. 2003, 144: 4562-4574. 10.1210/en.2003-0567.

    Article  CAS  PubMed  Google Scholar 

  37. Buterin T, Koch C, Naegeli H: Convergent transcriptional profiles induced by endogenous estrogen and distinct xenoestrogens in breast cancer cells. Carcinogenesis. 2006, 27: 1567-1578. 10.1093/carcin/bgi339.

    Article  CAS  PubMed  Google Scholar 

  38. Mardis ER: ChIP-seq: welcome to the new frontier. Nature methods. 2007, 4: 613-614. 10.1038/nmeth0807-613.

    Article  CAS  PubMed  Google Scholar 

  39. Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, Nussbaum C, Myers RM, Brown M, Li W, Liu XS: Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008, 9 (9): R137-10.1186/gb-2008-9-9-r137.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Bailey TL, Elkan C: Fitting a mixture model by expectation maximization to discover motifs in biopolymers. Proc Int Conf Intell Syst Mol Biol. 1994, 2: 28-36.

    CAS  PubMed  Google Scholar 

  41. Portales-Casamar E, Thongjuea S, Kwon AT, Arenillas D, Zhao X, Valen E, Yusuf D, Lenhard B, Wasserman WW, Sandelin A: JASPAR 2010: the greatly expanded open-access database of transcription factor binding profiles. Nucleic Acids Res. 2010, 38: D105-10. 10.1093/nar/gkp950.

    Article  CAS  PubMed  Google Scholar 

  42. Mahony S, Auron PE, Benos PV: DNA familial binding profiles made easy: comparison of various motif alignment and clustering strategies. PLoS Comput Biol. 2007, 3: e61-10.1371/journal.pcbi.0030061.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Safe S, Kim K: Non-classical genomic estrogen receptor (ER)/specificity protein and ER/activating protein-1 signaling pathways. J Mol Endocrinol. 2008, 41: 263-275. 10.1677/JME-08-0103.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Vyhlidal C, Samudio I, Kladde MP, Safe S: Transcriptional activation of transforming growth factor alpha by estradiol: requirement for both a GC-rich site and an estrogen response element half-site. J Mol Endocrinol. 2000, 24: 329-338. 10.1677/jme.0.0240329.

    Article  CAS  PubMed  Google Scholar 

  45. Webb P, Nguyen P, Valentine C, Lopez GN, Kwok GR, McInerney E, Katzenellenbogen BS, Enmark E, Gustafsson JA, Nilsson S, Kushner PJ: The estrogen receptor enhances AP-1 activity by two distinct mechanisms with different requirements for receptor transactivation functions. Mol Endocrinol. 1999, 13: 1672-1685. 10.1210/me.13.10.1672.

    Article  CAS  PubMed  Google Scholar 

  46. Leung Y, Mak P, Hassan S, Ho S: Estrogen receptor (ER)-beta isoforms: a key to understanding ER-beta signaling. Proc Natl Acad Sci USA. 2006, 103: 13162-13167. 10.1073/pnas.0605676103.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Weigelt K, Moehle C, Stempfl T, Weber B, Langmann T: An integrated workflow for analysis of ChIP-chip data. Biotechniques. 2008, 45: 131-132. 10.2144/000112819. 134, 136 passim

    Article  CAS  PubMed  Google Scholar 

  48. Aulmann S, Bläker H, Penzel R, Rieker RJ, Otto HF, Sinn HP: CTCF gene mutations in invasive ductal breast cancer. Breast Cancer Res Treat. 2003, 80: 347-352. 10.1023/A:1024930404629.

    Article  CAS  PubMed  Google Scholar 

  49. Mattingly KA, Ivanova MM, Riggs KA, Wickramasinghe NS, Barch MJ, Klinge CM: Estradiol stimulates transcription of nuclear respiratory factor-1 and increases mitochondrial biogenesis. Mol Endocrinol. 2008, 22: 609-622. 10.1210/me.2007-0029.

    Article  CAS  PubMed  Google Scholar 

  50. Granok H, Leibovitch B, Shaffer C: Ga-ga over GAGA factor. Curr Biol. 1995, 5: 238-241. 10.1016/S0960-9822(95)00048-0.

    Article  CAS  PubMed  Google Scholar 

  51. Van Steensel B, Delrow J, Bussemaker HJ: Genomewide analysis of Drosophila GAGA factor target genes reveals context-dependent DNA binding. Proc Natl Acad Sci USA. 2003, 100: 2580-2585. 10.1073/pnas.0438000100.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Bauer S, Grossmann S, Vingron M, Robinson PN: Ontologizer 2.0--a multifunctional tool for GO term enrichment analysis and data exploration. Bioinformatics. 2008, 24: 1650-1651. 10.1093/bioinformatics/btn250.

    Article  CAS  PubMed  Google Scholar 

  53. Al-Nakhle H, Burns PA, Cummings M, Hanby AM, Hughes TA, Satheesha S, Shaaban AM, Smith L, Speirs V: Estrogen receptor {beta}1 expression is regulated by miR-92 in breast cancer. Cancer Res. 2010, 70: 4778-4784. 10.1158/0008-5472.CAN-09-4104.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Bhat-Nakshatri P, Wang G, Collins NR, Thomson MJ, Geistlinger TR, Carroll JS, Brown M, Hammond S, Srour EF, Liu Y, Nakshatri H: Estradiol-regulated microRNAs control estradiol response in breast cancer cells. Nucleic Acids Res. 2009, 37: 4850-4861. 10.1093/nar/gkp500.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Castellano L, Giamas G, Jacob J, Coombes RC, Lucchesi W, Thiruchelvam P, Barton G, Jiao LR, Wait R, Waxman J, Hannon GJ, Stebbing J: The estrogen receptor-alpha-induced microRNA signature regulates itself and its transcriptional response. Proc Natl Acad Sci USA. 2009, 106: 15732-15737. 10.1073/pnas.0906947106.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Di Leva G, Gasparini P, Piovan C, Ngankeu A, Garofalo M, Taccioli C, Iorio MV, Li M, Volinia S, Alder H, Nakamura T, Nuovo G, Liu Y, Nephew KP, Croce CM: MicroRNA Cluster 221-222 and Estrogen Receptor {alpha} Interactions in Breast Cancer. J Natl Cancer Inst. 2010, 102: 706-721. 10.1093/jnci/djq102.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Maillot G, Lacroix-Triki M, Pierredon S, Gratadou L, Schmidt S, Bénès V, Roché H, Dalenc F, Auboeuf D, Millevoi S, Vagner S: Widespread estrogen-dependent repression of micrornas involved in breast tumor cell growth. Cancer Res. 2009, 69: 8332-8340. 10.1158/0008-5472.CAN-09-2206.

    Article  CAS  PubMed  Google Scholar 

  58. Wickramasinghe NS, Manavalan TT, Dougherty SM, Riggs KA, Li Y, Klinge CM: Estradiol downregulates miR-21 expression and increases miR-21 target gene expression in MCF-7 breast cancer cells. Nucleic Acids Res. 2009, 37: 2584-2595. 10.1093/nar/gkp117.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Adams BD, Guttilla IK, White BA: Involvement of microRNAs in breast cancer. Semin Reprod Med. 2008, 26: 522-536. 10.1055/s-0028-1096132.

    Article  CAS  PubMed  Google Scholar 

  60. Iorio MV, Ferracin M, Liu CG, Veronese A, Spizzo R, Sabbioni S, Magri E, Pedriali M, Fabbri M, Campiglio M, Ménard S, Palazzo JP, Rosenberg A, Musiani P, Volinia S, Nenci I, Calin GA, Querzoli P, Negrini M, Croce CM: MicroRNA gene expression deregulation in human breast cancer. Cancer Res. 2005, 65: 7065-7070. 10.1158/0008-5472.CAN-05-1783.

    Article  CAS  PubMed  Google Scholar 

  61. Iorio MV, Casalini P, Tagliabue E, Ménard S, Croce CM: MicroRNA profiling as a tool to understand prognosis, therapy response and resistance in breast cancer. Eur J Cancer. 2008, 44: 2753-2759. 10.1016/j.ejca.2008.09.037.

    Article  CAS  PubMed  Google Scholar 

  62. Yang SH, Liu R, Perez EJ, Wen Y, Stevens SM, Valencia T, Brun-Zinkernagel AM, Prokai L, Will Y, Dykens J, Koulen P, Simpkins JW: Mitochondrial localization of estrogen receptor beta. Proc Natl Acad Sci USA. 2004, 101: 4130-4135. 10.1073/pnas.0306948101.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Chen JQ, Eshete M, Alworth WL, Yager JD: Binding of MCF-7 cell mitochondrial proteins and recombinant human estrogen receptors alpha and beta to human mitochondrial DNA estrogen response elements. J Cell Biochem. 2004, 93: 358-373. 10.1002/jcb.20178.

    Article  CAS  PubMed  Google Scholar 

  64. Chen J, Russo PA, Cooke C, Russo IH, Russo J: ERbeta shifts from mitochondria to nucleus during estrogen-induced neoplastic transformation of human breast epithelial cells and is involved in estrogen-induced synthesis of mitochondrial respiratory chain proteins. Biochim Biophys Acta. 2007, 1773: 1732-1746. 10.1016/j.bbamcr.2007.05.008.

    Article  CAS  PubMed  Google Scholar 

  65. Felty Q, Singh KP, Roy D: Estrogen-induced G1/S transition of G0-arrested estrogen-dependent breast cancer cells is regulated by mitochondrial oxidant signaling. Oncogene. 2005, 24: 4883-4893. 10.1038/sj.onc.1208667.

    Article  CAS  PubMed  Google Scholar 

  66. Simpkins JW, Yang S, Sarkar SN, Pearce V: Estrogen actions on mitochondria--physiological and pathological implications. Mol Cell Endocrinol. 2008, 290: 51-59. 10.1016/j.mce.2008.04.013.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  67. Blobel GA, Sieff CA, Orkin SH: Ligand-dependent repression of the erythroid transcription factor GATA-1 by the estrogen receptor. Mol Cell Biol. 1995, 6: 3147-3153.

    Article  Google Scholar 

  68. Mao L, Wertzler KJ, Maloney SC, Wang Z, Magnuson NS, Reeves R: HMGA1 levels influence mitochondrial function and mitochondrial DNA repair efficiency. Mol Cell Biol. 2009, 29: 5426-5440. 10.1128/MCB.00105-09.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Massaad-Massade L, Tacine R, Dulauroy S, Reeves R, Barouki R: The functional interaction between HMGA1 and the estrogen receptor requires either the N- or the C-terminal domain of the receptor. FEBS Lett. 2004, 559: 89-95. 10.1016/S0014-5793(04)00032-8.

    Article  CAS  PubMed  Google Scholar 

  70. Pace P, Taylor J, Suntharalingam S, Coombes RC, Ali S: Human estrogen receptor beta binds DNA in a manner similar to and dimerizes with estrogen receptor alpha. J Biol Chem. 1997, 272: 25832-25838. 10.1074/jbc.272.41.25832.

    Article  CAS  PubMed  Google Scholar 

  71. Pettersson K, Grandien K, Kuiper GG, Gustafsson JA: Mouse estrogen receptor beta forms estrogen response element-binding heterodimers with estrogen receptor alpha. Mol Endocrinol. 1997, 11: 1486-1496. 10.1210/me.11.10.1486.

    CAS  PubMed  Google Scholar 

  72. Li X, Huang J, Yi P, Bambara RA, Hilf R, Muyan M: Single-chain estrogen receptors (ERs) reveal that the ERalpha/beta heterodimer emulates functions of the ERalpha dimer in genomic estrogen signaling pathways. Mol Cell Biol. 2004, 24: 7681-7694. 10.1128/MCB.24.17.7681-7694.2004.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Andrews NC, Faller DV: A rapid micropreparation technique for extraction of DNA-binding proteins from limiting numbers of mammalian cells. Nucleic Acids Res. 1991, 19: 2499-10.1093/nar/19.9.2499.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. Frezza C, Cipolat S, Scorrano L: Organelle isolation: functional mitochondria from mouse liver, muscle and cultured fibroblasts. Nat Protoc. 2007, 2: 287-295. 10.1038/nprot.2006.478.

    Article  CAS  PubMed  Google Scholar 

  75. Caporali S, Imai M, Altucci L, Cancemi M, Caristi S, Cicatiello L, Matarese F, Penta R, Sarkar DK, Bresciani F, Weisz A: Distinct signaling pathways mediate stimulation of cell cycle progression and prevention of apoptotic cell death by estrogen in rat pituitary tumor PR1 cells. Mol Biol Cell. 2003, 14: 5051-5059. 10.1091/mbc.E03-05-0303.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  76. Bermont L, Lamielle-Musard F, Chezy E, Weisz A, Adessi GL: 17beta-estradiol inhibits forskolin-induced vascular endothelial growth factor promoter in MCF-7 breast adenocarcinoma cells. J Steroid Biochem Mol Biol. 2001, 78: 343-349. 10.1016/S0960-0760(01)00111-X.

    Article  CAS  PubMed  Google Scholar 

  77. Licznar A, Caporali S, Lucas A, Weisz A, Vignon F, Lazennec G: Identification of genes involved in growth inhibition of breast cancer cells transduced with estrogen receptor. FEBS Lett. 2003, 553: 445-450. 10.1016/S0014-5793(03)01090-1.

    Article  CAS  PubMed  Google Scholar 

  78. Bonapace IM, Addeo R, Altucci L, Cicatiello L, Bifulco M, Laezza C, Salzano S, Sica V, Bresciani F, Weisz A: 17 beta-Estradiol overcomes a G1 block induced by HMG-CoA reductase inhibitors and fosters cell cycle progression without inducing ERK-1 and -2 MAP kinases activation. Oncogene. 1996, 12: 753-763.

    CAS  PubMed  Google Scholar 

  79. Cicatiello L, Scafoglio C, Altucci L, Cancemi M, Natoli G, Facchiano A, Iazzetti G, Calogero R, Biglia N, De Bortoli M, Sfiligoi C, Sismondi P, Bresciani F, Weisz A: A genomic view of estrogen actions in human breast cancer cells by expression profiling of the hormone-responsive transcriptome. J Mol Endocrinol. 2004, 32: 719-775. 10.1677/jme.0.0320719.

    Article  CAS  PubMed  Google Scholar 

  80. Scafoglio C, Ambrosino C, Cicatiello L, Altucci L, Ardovino M, Bontempo P, Medici N, Molinari AM, Nebbioso A, Facchiano A, Calogero RA, Elkon R, Menini N, Ponzone R, Biglia N, Sismondi P, De Bortoli M, Weisz A: Comparative gene expression profiling reveals partially overlapping but distinct genomic actions of different antiestrogens in human breast cancer cells. J Cell Biochem. 2006, 98: 1163-1184. 10.1002/jcb.20820.

    Article  CAS  PubMed  Google Scholar 

  81. Ravo M, Mutarelli M, Ferraro L, Grober OM, Paris O, Tarallo R, Vigilante A, Cimino D, De Bortoli M, Nola E, Cicatiello L, Weisz A: Quantitative expression profiling of highly degraded RNA from formalin-fixed, paraffin-embedded breast tumor biopsies by oligonucleotide microarrays. Lab Invest. 2008, 88: 430-440. 10.1038/labinvest.2008.11.

    Article  CAS  PubMed  Google Scholar 

  82. Cartharius K, Frech K, Grote K, Klocke B, Haltmeier M, Klingenhoff A, Frisch M, Bayerlein M, Werner T: MatInspector and beyond: promoter analysis based on transcription factor binding sites. Bioinformatics. 2005, 21: 2933-2942. 10.1093/bioinformatics/bti473.

    Article  CAS  PubMed  Google Scholar 

  83. Ho Sui SJ, Fulton DL, Arenillas DJ, Kwon AT, Wasserman WW: oPOSSUM: identification of over-represented transcription factor binding sites in co-expressed genes. Nucleic Acids Res. 2005, 33: 3154-64. 10.1093/nar/gki624.

    Article  PubMed  PubMed Central  Google Scholar 

  84. Saeed AI, Sharov V, White J, Li J, Liang W, Bhagabati N, Braisted J, Klapa M, Currier T, Thiagarajan M, Sturn A, Snuffin M, Rezantsev A, Popov D, Ryltsov A, Kostukovich E, Borisovsky I, Liu Z, Vinsavich A, Trush V, Quackenbush J: TM4: a free, open-source system for microarray data management and analysis. Biotechniques. 2003, 34: 374-378.

    CAS  PubMed  Google Scholar 

  85. Kuhn RM, Karolchik D, Zweig AS, Wang T, Smith KE, Rosenbloom KR, Rhead B, Raney BJ, Pohl A, Pheasant M, Meyer L, Hsu F, Hinrichs AS, Harte RA, Giardine B, Fujita P, Diekhans M, Dreszer T, Clawson H, Barber GP, Haussler D, Kent WJ: The UCSC Genome Browser Database: update 2009. Nucleic Acids Res. 2009, 37 (Database issue): D755-761. 10.1093/nar/gkn875.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Work supported by the European Union (CRESCENDO I.P., contract n.er LSHM-CT2005-018652), the Italian Association for Cancer Research (grant IG-8586) and the Italian Ministry for Education, University and Research (grant PRIN 2008CJ4SYW_004), Fondazione per il Sud (grant 2009-PdP-22) to AW. MR, GN and RT are recipient of fellowships from Fondazione per il Sud. The authors thank Dr Rosario Casale for technical assistance and Francesca Rizzo for revision of the manuscript.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Alessandro Weisz.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

OMVG, MM, GG, MR, LC, MRDF, LF, GN, MFP, OP and RT participated in the conception and design of the study, performed in vivo and in vitro experimental work (OP, GN, MFP, RT), statistical and in silico data analyses (OMVG, MM, GG, MRDF) and participated in drafting the manuscript. OP, LC and MR prepared and tested the sequencing libraries. LC, LF and MR carried out the microarray experiments and participated in data analyses. SL, GPS and VB participated in the conception and design of the study and performed massively parallel sequencing. AW coordinated the project, participated in conception and design of the study and participated in drafting and finalization of the manuscript. All authors read and approved the final manuscript.

Oli MV Grober, Margherita Mutarelli contributed equally to this work.

Electronic supplementary material

12864_2010_3175_MOESM1_ESM.BED

Additional File 1:Grober_ER-alpha_Binding_Sites ERα binding sites. The UCSC genome BED formatted file lists the chromosome, start coordinate, stop coordinate and identifier of the ERα binding sites. (BED 222 KB)

12864_2010_3175_MOESM2_ESM.BED

Additional File 2:Grober_ER-beta_Binding_Sites ERβ binding sites. The UCSC genome BED formatted lists the chromosome, start coordinate, stop coordinate and identifier of the ERβ binding sites. (BED 359 KB)

12864_2010_3175_MOESM3_ESM.XLS

Additional File 3:Sheet 1: Additional Table S1 Differentially Regulated Genes. Overview of the 921 genes differentially regulated by E2 in ERβ+ vs ERβ- cells, containing the following additional information: gene set membership, symbol, Entrez ID, gene name and expression values (fold-change). Sheet 2: Additional Table S2 TFBS enrichment matrix. The worksheet shows the over-representation values for the TF binding matrices from ERE- binding sites. Sheet 3: Additional Table S3 Primary ERβ target genes, showing ERβ binding sites within 10 kb of the TU. Overview of the 424 putative ERβ primary gene targets containing the following additional information: gene set membership, category of ERβ binding sites, Symbol, Entrez ID, Gene Name, TU Coordinates, ERα Binding Sites Coordinates and ERβ Binding Sites Coordinates. Sheet 4: Additional Table S4 Primary ERα target genes, showing ERα binding sites within 10 kb of the TU in ERβ- cells. Overview of the 228 putative ERα primary gene targets containing the following additional information: Symbol, Entrez ID, Gene Name, TU Coordinates. (XLS 440 KB)

12864_2010_3175_MOESM4_ESM.DOC

Additional File 4:GO analysis of primary ERβ target genes. Containing the following information: Biological process, Gene Ontology term, Name, Count in total GO population, Count in selected genes, % genes and p-value. (DOC 112 KB)

12864_2010_3175_MOESM5_ESM.DOC

Additional File 5:GO analysis of primary ERα target genes. Containing the following information: Biological process, Gene Ontology term, Name, Count in total GO population, Count in selected genes, % genes and p-value. (DOC 88 KB)

12864_2010_3175_MOESM6_ESM.DOC

Additional File 6:ERβ binding sites in proximity of miRNA loci. Containing the following information: miRNA name, ID of ER-β binding site upstream, Distance from the closest ERβ binding site upstream, ID of ER-β binding site downstream and Distance from the closest ERβ binding site downstream. (DOC 118 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Grober, O.M., Mutarelli, M., Giurato, G. et al. Global analysis of estrogen receptor beta binding to breast cancer cell genome reveals an extensive interplay with estrogen receptor alpha for target gene regulation. BMC Genomics 12, 36 (2011). https://doi.org/10.1186/1471-2164-12-36

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-12-36

Keywords