- Research article
- Open Access
Systematic identification of Ctr9 regulome in ERα-positive breast cancer
BMC Genomicsvolume 17, Article number: 902 (2016)
We had previously identified Ctr9, the key scaffold subunit of the human RNA polymerase II (RNAPII) associated factor complex (PAFc), as a key factor regulating a massive ERα target gene expression and ERα-positive breast cancer growth. Furthermore, we have shown that knockdown of Ctr9 reduces ERα protein stability and decreases the occupancy of ERα and RNAPII at a few ERα-target loci. However, it remains to be determined whether Ctr9 controls ERα-target gene expression by regulating the global chromatin occupancy of ERα and RNAPII in the presence of estrogen.
In this study, we determined the genome-wide ERα and RNAPII occupancy in response to estrogen treatment and/or Ctr9 knockdown by performing chromatin immunoprecipitation coupled with high-throughput sequencing (ChIP-seq). We found that loss of Ctr9 dramatically decreases the global occupancy of ERα and RNAPII, highlighting the significance of Ctr9 in regulating estrogen signaling in ERα-positive breast cancer cells. Combining this resource with previously published genomic data sets, we identified a unique subset of ERα and Ctr9 target genes, and further delineated the independent function of Ctr9 from other subunits in PAFc when regulating transcription.
Our data demonstrated that Ctr9, independent of other PAFc subunits, controls ERα-target gene expression by regulating global chromatin occupancies of ERα and RNAPII.
Estrogen receptor α (ERα) is the paramount transcription factor in ERα-positive breast cancer, which constitutes 70% of all human breast cancers. Emerging evidence has shown that the execution of ERα signaling cascade is through multiple pathways to regulate the expression of ERα target genes in the presence or absence of ligand, i.e. 17β-estradiol (estrogen; E2) [1, 2]. The genomic ERα signaling pathway involves ligand binding, receptor dimerization and recognition of cognate estrogen response element (ERE) which consists of a consensus sequence containing the inverted repeat GGTCAxxxTGACC (where x is any nucleotide), located in either the distal or promoter regions of ERα-target genes . The ERα-dependent transcriptional program is delicately regulated by a large number of cofactors including pioneering factors, chromatin modifiers as well as basal transcriptional machineries . Perturbations of ERα-dependent transcriptional networks by the aberrant expression of ERα cofactors have been strongly linked to breast tumorigenesis . Thus, comprehensively understanding their mechanisms of action could lead to therapeutic interventions for the treatment of ERα-positive breast cancer.
Human RNA polymerase II (RNAPII) associated factor complex (PAFc) is composed of five subunits including PAF1, Ctr9, Cdc73 (also known as parafibromin), Leo1, and eukaryotic specific Ski8 (also known as WDR61) . PAFc is evolutionarily conserved from yeast to Drosophila and humans and it participates in the step-wise regulation of transcription from transcriptional initiation, elongation, to termination [6–8]. The PAFc regulated transcriptional elongation is coupled with histone modifications such as histone H2B mono-ubiquitination (H2Bub1), H3K4 tri-methylation (H3K4me3), and H3K36 tri-methylation (H3K36me3) via protein-protein interaction of PAFc subunits with the respective histone modifying enzymes . Additionally, PAFc plays a role in mRNA processing and maturation, through the maintenance of proper poly (A) tail length . Emerging evidence has demonstrated the diverse functions of human PAFc as well as its individual subunits in various cellular milieus. For instance, the human PAFc is required for the maintenance of embryonic stem cell identity , neuronal migration in mammalian brain , antiviral and pro-inflammatory response , and oncogenesis regulation [14, 15]. Our studies have demonstrated the involvement of Ctr9, a key scaffold subunit of human PAFc, in ERα-positive breast cancer progression and ERα target gene expression . Specifically, using loss-of-function approach, we observed that depletion of Ctr9 led to apparent morphological change, decrease of proliferation, reduced colony formation, and impaired ERα-target gene expression in ERα-positive breast cancer cells . Moreover, Ctr9 regulates ERα protein stability and facilitates RNAPII recruitment and transcription-coupled histone modifications such as H2Bub1, H3K27ac, and H3K36me3 at the selected ERα targets . However, the global control of Ctr9 on ERα occupancy and RNAPII recruitment awaits investigation to deeply understand the functional significance of Ctr9 in ERα transcription network in breast cancer.
Profiling the genome-wide RNAPII and ERα chromatin occupancy affected by Ctr9 is advantageous over gene expression study to elucidate the global effects of Ctr9 on ERα-dependent transcriptional program and to identify the target genes regulated by Ctr9. In agreement with several published gene expression studies in MCF7 cells [17–28], our gene expression profiling identified less than 2000 E2-responsive genes [16, 29]. This is in contrast to the large scale ERα ChIP profiling where much more ERα binding sites across the genome were identified [18, 19, 21, 23, 25, 30–32]. For example, the Myles Brown lab identified 3665 ERα binding sites  and the Hendrik G Stunnenberg lab identified 10,205 ERα binding sites  in response to E2 treatment, by ChIP-chip and ChIP-seq, respectively. The discordance in E2-responsive transcription and ERα binding events is probably due to the prevalence of enhancers over promoters in regulating transcription and also suggests that some ERα binding events do not fulfill gene regulation. In addition, the mRNA level of a target gene may not faithfully reflect the transcriptional event due to the mRNA degradation. Therefore profiling the genome-wide RNAPII occupancy would provide a direct readout of transcriptional activity and thus yield mechanistic insights into the transcriptional regulation by estrogen and Ctr9 in ERα-positive breast cancer.
In this study, we performed chromatin immunoprecipitation followed by massive parallel sequencing (ChIP-seq) and demonstrated that Ctr9 knockdown (KD) elicits global decrease of the genome-wide occupancy of ERα and RNAPII, consistent with the function of Ctr9 in regulating ERα protein stability and ERα target gene expression. Integrative analyses of the ChIP-seq and microarray gene expression data sets identified the primary ERα and Ctr9 target genes, which have prognostic value of breast cancer outcome. Comparative analysis of RNAPII genome-wide occupancy alteration resulted from Ctr9 KD and PAF1 KD provided the first line of evidence in support of the PAFc independent function of Ctr9 on transcription regulation.
Ctr9 knockdown triggers genome-wide decrease of ERα occupancy
To investigate the estrogen and Ctr9 mediated transcription regulation in ERα-positive breast cancer cells, we used the MCF7-tet-on-shCtr9 inducible knockdown cell line previously developed in our lab as a model system (Fig. 1a) [16, 33]. ChIP-seq was first performed using an anti-ERα antibody in MCF7-tet-on-shCtr9 cells pre-treated with vehicle or doxycycline (Dox) followed by DMSO control or E2 induction. Classical ERα target genes, such as TFF1 and GREB1, showed significant enrichment of ERα over a narrow range in their enhancer and promoter regions in response to E2 treatment (Fig. 1b and c), These results are highly consistent with several previous reports [18, 25, 32]. Consistent with our previous ChIP-qPCR results , we found that Ctr9 KD by Dox treatment significantly reduced ERα binding at the enhancer and promoter regions of TFF1 and GREB1 genes (Fig. 1b and c).
In accordance with others’ findings that substantial ERα binding is induced by estrogen treatment [25, 32], we identified a total of 5264 significant ERα binding sites upon E2 induction (Fig. 2a). We visually inspected the genomic locations of these 5264 ERα binding sites using Cistrome , and found that the majority of the binding sites are located within introns (42.1%) or distal regions (44.4%) and only a small percentage (7.0%) in promoter regions (Additional file 1: Figure S1) in agreement with previous reports [23, 25, 32], validating the reliability of our ChIP-seq datasets. Comparison between the four treatment conditions revealed that, in the presence of E2, Ctr9 KD decreased ERα binding at 3330 sites and only increased ERα binding at 26 sites (Fig. 2a), suggesting that Ctr9 is in general a positive regulator of genome-wide ERα occupancy. Aggregate plots further confirmed that loss of Ctr9 dramatically decreased the genome-wide ERα binding (Fig. 2b). As a further support, the heatmaps of ERα ChIP-seq density showed that Ctr9 depletion significantly reduced ERα occupancy at the majority of ERα binding sites (Fig. 2c). To obtain a more intuitive view of the effect of Ctr9 KD on global ERα binding, we performed a Venn diagram analysis for 5264 E2-induced ERα binding sites in control cells, 2848 E2-induced ERα binding sites in Ctr9 KD cells, and 3330 ERα binding sites decreased by Ctr9 KD (Fig. 2d). Among the E2-induced ERα binding sites in control cells, 1630 was significantly reduced by Ctr9 KD, although E2 treatment can still induce the binding of ERα at these sites (Fig. 2d). Much stronger effect was observed for 1572 ERα binding sites where Ctr9 KD almost completely abrogated the ERα binding (Fig. 2d). Therefore, we termed the 1630 sites and 1572 sites less and more sensitive to Ctr9 KD, respectively. ERE motif analyses identified that 65 and 55% of these two groups of ERα cistrome contains at least half ERE, respectively. We had previously shown that Ctr9 KD decreased ERα protein level thus expected reduced ERα occupancy on all targets. These results instead support the notion that Ctr9 elicits target-specific regulation towards ERα target genes, underscoring the complexity of Ctr9 function in the regulation of ERα-mediated transcriptional program beyond its conventional roles in the PAFc to regulate transcription initiation and elongation.
Ctr9 knockdown alters the repertoire of ERα cistrome as illustrated by motif analyses
Next, we further investigated the impact of Ctr9 KD on ERα cistrome. We interrogated the sequence of the aforementioned three groups of ERα binding sites (Fig. 2d) for overrepresentation of DNA motifs using the MEME suite . As expected, the most prevalent motifs at the ERα binding sites in all three groups are EREs (Fig. 2e, No. 1). Moreover, among the top five DNA motifs for E2-induced ERα cistrome, we found significant enrichment of DNA sequences featuring progesterone responsive element (PRE) (Fig. 2e, No. 3) and FOXA1 binding motif (Fig. 2e, No. 5). These results are consistent with the previous reports that progesterone receptor associates with ERα and modulates ERα action in breast cancer , and FOXA1 is the pioneer factor for shaping the genome-wide distribution of ERα , respectively. Importantly, PREs are no longer enriched in the E2-induced ERα cistrome in the context of Ctr9 KD, this is consistent with our previous finding that Ctr9 KD significantly decreases the expression of progesterone receptor . Furthermore, all three groups of ERα cistrome displayed overlapping yet different DNA motif enrichment, suggesting that Ctr9 might affect or cooperate with other factors to regulate ERα action and estrogen signaling in ERα-positive breast cancer.
Ctr9 positively regulates RNA polymerase II binding in MCF7 cells
Given that the vast majority of ERα binding sites are located at great distances from the annotated genes (Additional file 1: Figure S1), identification of ERα direct target genes simply based on the ERα binding sites has been problematic. To identify genes in response to E2 and/or Ctr9 KD in a more intuitive manner, we performed ChIP-seq for RNAPII in MCF7-tet-on-shCtr9 cells under four different conditions (Vehicle, DMSO; Vehicle, E2; Dox, DMSO; Dox, E2). RNAPII occupancy throughout the gene body represents a direct readout of the transcriptional activity and thus provides more intuitional evidence for estrogen and Ctr9 mediated transcription regulation. Classical ERα target genes, for example TFF1 and GREB1, showed a dramatic increase of RNAPII occupancy over their enhancer, promoter, and gene body regions upon E2 treatment (Fig. 3a). Ctr9 KD by Dox treatment significantly reduced RNAPII occupancy at the aforementioned regions of TFF1 and GREB1 genes (Fig. 3a), in line with our previous ChIP-qPCR results and microarray gene expression results . Genomic location analysis of E2-induced RNAPII binding sites showed that unlike ERα, RNAPII binding displayed more dispersed occupancy across the genome (Additional file 1: Figure S2).
At a global scale, E2 treatment and Ctr9 KD triggered redistribution of genome-wide RNAPII binding pattern (Fig. 3b). To directly identify genes responding to E2 treatment and/or Ctr9 KD, we first defined RNAPII-target gene using the criteria that one or more RNAPII binding sites can be found within its gene body (TSS, TTS, 5′-UTR, exon, intron, and 3′-UTR). Comparisons between the different conditions showed that, in response to E2 treatment, 2353 genes showed an increase and 562 genes showed a decrease in RNAPII occupancy (Fig. 3c). In contrast, E2 stimulation only increased RNAPII occupancy over 894 genes in Ctr9 KD cells (Fig. 3c), suggesting that Ctr9 positively regulates E2-induced RNAPII binding. The positive regulatory role of Ctr9 was further supported by the decrease of RNAPII occupancy over thousands of genes upon Ctr9 KD (Fig. 3c). Aggregate plot also confirmed that loss of Ctr9 significantly decreased the genome-wide RNAPII occupancy particularly at the transcription start sites (TSS) (Fig. 3d). We previously performed microarray gene expression analysis in MCF7-tet-on-shCtr9 cells treated with vehicle or Dox followed by DMSO control or E2 induction to identify E2 and Ctr9 regulated transcriptome [16, 29]. To gain additional insights into the regulation of RNAPII occupancy by Ctr9, we further examined the RNAPII occupancy for Ctr9 KD decreased and increased genes selected from our microarray data set (Fig. 3e). Aggregate plots showed that RNAPII occupancy was drastically decreased by Ctr9 KD for Ctr9 KD decreased genes, while increased particularly within the gene body by Ctr9 KD for Ctr9 KD increased genes, suggesting the context dependent role of Ctr9 in regulating RNAPII occupancy in MCF7 cells. To directly assess the impact of Ctr9 KD on estrogen signaling, we compared the genes with RNAPII binding increased by E2 treatment and genes with RNAPII binding decreased by Ctr9 KD in the presence of E2 (Fig. 3f). The Venn diagram analysis showed that about half of E2-responsive genes (1156 genes) displayed RNAPII binding decrease upon Ctr9 KD (Fig. 3f). Besides, Ctr9 KD also decreased RNAPII occupancy at 1583 non-E2-responsive genes (Fig. 3f), suggesting that Ctr9-dependent RNAPII chromatin occupancy is not restricted to the estrogen regulated genes in breast cancer.
Integrative analyses of the combined ChIP-seq and microarray data reveal the ERα and Ctr9 signature which has clinical significance
To directly identify ERα and Ctr9 targets, we performed integrative analyses of the combination of ChIP-seq and microarray data. Comparing genes with RNAPII binding increased by E2 and genes with increased expression upon E2 treatment revealed an overlap of 279 genes, termed E2-activated genes (Fig. 4a and Additional file 2: Table S1). Further comparing with genes with ERα binding increased by E2 treatment identified 88 common genes among the three groups including some well-known ERα target genes such as TFF1 and GREB1 (Fig. 4a and Additional file 2: Table S1). Thus, we defined the 88-gene set as primary ERα target signature. Note that many genes with increased expression upon E2 treatment failed to display elevated RNAPII binding (Fig. 4a), probably because ChIP-seq with short E2 treatment (45 min) predominantly identifies early responsive target genes whereas gene expression profiling at 4 h after E2 treatment identifies the late responsive and some indirect targets. As ERα signature, the identified 88 gene set is highly expressed in ERα-positive human breast tumors (Fig. 4b) and the high expression level of the ERα signature is associated with poor overall survival in ERα-positive, lymph node-negative breast cancer patients (Fig. 4c).
Given that ERα-positive breast cancer depends on estrogen for growth and target gene transcription regulation, we next sought to identify primary Ctr9 targets in the presence of E2 by comparing genes with RNAPII binding decreased by Ctr9 KD and genes with decreased expression upon Ctr9 KD (Fig. 4d). The Venn diagram analysis identified 240 common genes between these two data sets. We defined the 240 gene set as Ctr9 signature (Fig. 4d and Additional file 2: Table S1). To determine if the Ctr9 signature has a broader clinical significance, we mined breast cancer outcome-linked gene expression data using the Gene Expression-Based Outcome for Breast Cancer Online (GOBO) tool . In line with our previous finding that Ctr9 is highly expressed in ERα-positive breast cancer , the Ctr9 signature is also enriched in ERα-positive breast tumors (Fig. 4e). Kaplan-Meier analysis demonstrated that high expression levels of the Ctr9 signature strongly associates with poor overall survival in ERα-positive, lymph node-negative breast cancer patients (Fig. 4f), further highlighting the functional significance of Ctr9 in human breast cancer. Furthermore, comparing the Ctr9 signature with E2-responsive genes revealed an overlap of 59 genes (Fig. 4g and Additional file 2: Table S1). Taken together, these integrative analyses provided a deeper understanding of estrogen and Ctr9 regulated transcriptome in ERα-positive breast cancer cells.
Distinct roles of Ctr9 in regulating RNAPII occupancy independent of PAFc
Our previous study implies that Ctr9 might function independently of the other subunits of PAFc in regulating target gene expression and breast cancer progression . To further distinguish the functions of Ctr9 from the other PAFc subunits, we analyzed recently published RNAPII ChIP-seq data obtained in wild-type and PAF1 KD MCF7 cells cultured in complete medium . Similar to Ctr9 KD, PAF1 KD resulted in decrease of RNAPII occupancy over 5135 genes and increase over only 190 genes (Fig. 5a), suggesting that PAF1 is a positive regulator of RNAPII genomic binding in MCF7 cells. Aggregate plot and heatmap of ChIP-seq density also confirmed that loss of PAF1 significantly decreases genome-wide RNAPII occupancy (Fig. 5b and c). To examine whether PAF1 and Ctr9 exhibit the same regulation on genome-wide RNAPII occupancy, we compared the genes that elicited decreased RNAPII binding by knocking down PAF1 and Ctr9; the results showed an overlap of 1325 genes (Fig. 5d). Interestingly, a large number of PAF1 regulated genes and Ctr9 regulated genes are mutually exclusive (Fig. 5d), indicating that Ctr9 functions independently from PAF1. Similarly, comparison of genes that elicited increased RNAPII binding by knocking down PAF1 and Ctr9 revealed only 12 common genes (Fig. 5e). Taken the Ctr9 target gene IGFBP4 and GREB1 as examples, Ctr9 KD dramatically decreased the occupancy of RNAPII across the gene body whereas PAF1 KD elicited no effect on RNAPII occupancy at IGFBP4 gene and increased RNAPII occupancy at GREB1 gene, respectively (Fig. 5f). To further validate our ChIP-seq findings, we performed validation experiments by ChIP-qPCR to examine the RNAPII occupancy at the transcription start site (TSS) and transcribed region (TR) of GREB1 and IGFBP4 genes in MCF7 cells stably expressing control shRNA or shRNAs targeting Ctr9, Cdc73, Ski8, and PAF1 (Fig. 6). The results showed that Ctr9 KD with two independent shRNAs (shCtr9 #3 and shCtr9 #5) significantly reduced the RNAPII occupancy at TSS and TR regions of GREB1 and IGFBP4 genes, whereas KD of the other PAFc subunits elicited little effect on the RNAPII occupancy at the same regions (Fig. 6). Altogether, these findings suggest that Ctr9 may have a unique role in regulating RNAPII chromatin association independent of PAFc.
Estrogen signaling in ERα-positive breast cancer is tightly regulated by ERα itself as well as its associated co-regulator proteins. It is essential to understand the ERα target gene network and the transcription regulation mechanisms. In this study, we have applied ChIP-seq to determine the genome-wide alterations of ERα and RNAPII occupancy in response to E2 induction and/or KD of Ctr9. The major finding is that loss of Ctr9 led to dramatic decrease of genome-wide occupancy of both ERα and RNAPII, supporting a positive regulatory role of Ctr9 in mediating estrogen-dependent transcriptional regulation. Integrative analyses of the combined ChIP-seq and microarray gene expression data led to the identification of 88 primary ERα targets and 240 primary Ctr9 targets which both have significant clinical relevance. We acknowledge that ChIP-seq for Ctr9 in MCF7 cells could provide a more precise identification of direct Ctr9 target genes. However, Ctr9 ChIP-seq has been technically challenging. We performed the Ctr9 ChIP-seq using the commercially available anti-Ctr9 antibody. However, ChIP-seq signals were very weak with less than 30 peaks identified (data not shown). This is consistent with the report by Fei Xavier Chen, et al.  in which this anti-Ctr9 antibody also failed in ChIP-seq. Therefore, ChIP-grade antibodies for Ctr9 need to be developed for future ChIP-seq experiments.
Emerging evidence has shown that some transcriptional cofactors can exert a global impact on ERα target gene expression and/or genome-wide ERα binding [18, 36, 39–41]. For example, loss-of-function approach coupled with gene expression microarray profiling demonstrated that BRD4, KDM3A, and Cohesin are required for the expression of a large portion of ERα target genes [18, 40–42]. GATA3 and FOXA1 have been implied as the pioneer factors for priming the genome-wide distribution of ERα [36, 39]. Although Ctr9 acts as a master regulator of estrogen signaling by stabilizing ERα protein , controlling genome-wide ERα and RNAPII binding, and affecting the target gene expression, the importance of Ctr9 in the estrogen signaling network orchestrated by the family of transcriptional factors and cofactors is of great interest and awaits future investigation. In this regard, our genomic analyses performed in this study provided valuable resources for possible comparison with the other factors regulating ERα transcriptome. Moreover, given that ERα protein level is regulated by Ctr9, we cannot distinguish whether the global change in ERα cistrome upon Ctr9 KD is specific to loss of Ctr9 or is due to reduced ERα protein itself, which awaits further clarification in the future.
Our previous study has also shown that Ctr9 modulates MCF7 cell morphology, cell proliferation, and ERα target gene expression . However, we could not distinguish whether Ctr9 functions independently or in the context of PAFc. Here, by comparing the alterations in RNAPII occupancy resulted from PAF1 KD and Ctr9 KD, we provided the first evidence showing that PAF1 and Ctr9, despite both positively regulate genome-wide RNAPII binding, might regulate RNAPII binding through different mechanisms in MCF7 cells. Moreover, ChIP-seq analyses in human acute myeloid leukemia THP1 cells showed that PAF1 occupancy generally overlaps with the promoter-proximal RNAPII peaks whereas Ctr9 occupancy does not , reinforcing the PAFc independent functions of Ctr9.
In this study, we profiled the genome-wide binding patterns of ERα and RNAPII in response to estrogen stimulation and/or Ctr9 KD. We showed that Ctr9 acts as a master regulator of the gene transcription program in ERα-positive breast cancer cells by affecting the global occupancy of ERα and RNAPII. Combination of ChIP-seq and the previous microarray gene expression data revealed primary target genes of ERα and Ctr9. Moreover, comparing RNAPII binding alterations between Ctr9 KD and PAF1 KD identified a unique subset of Ctr9 target genes that implicate PAFc-independent function of Ctr9. Altogether, our study provides a molecular basis to further understand the transcriptional control in ERα-positive breast cancers.
MCF7-tet-on-shCtr9 cells were generated previously  and maintained in Dulbecco’s modified medium (DMEM) (Life Technologies) supplemented with 10% fetal bovine serum (FBS) (Life Technologies). Cells were cultured at 37 °C in a humidified atmosphere containing 5% CO2. Doxycycline (Dox) was purchased from Clontech (Mountain View, CA) and used at a final concentration of 500 ng/ml. Dimethyl sulfoxide (DMSO) and 17 β-estradiol (E2) were purchased from Sigma (St. Louis, MO). Prior to chromatin immunoprecipitation (ChIP) assays, MCF7-tet-on-shCtr9 cells were cultured in DMEM supplemented with 10% FBS in the absence or presence of 500 ng/mL Dox for 4 days, followed by continuing cultured in phenol red-free DMEM supplemented with 5% charcoal-dextran-treated FBS in the absence or presence of 500 ng/mL Dox for another 3 days. Cells were then treated with DMSO or 10 nM E2 for 45 min. Antibodies used for ChIP-seq in this study were anti-ERα (sc-543, HC-20) purchased from Santa Cruz Biotechnology and anti-RNA polymerase II (#05-623, clone CTD4H8) purchased from Millipore.
Chromatin immunoprecipitation (ChIP)
ChIP assays were performed as described previously . In detail, cells were grown in 15-cm dish as described above and cross-linked for 10 min at room temperature by the addition of 37% formaldehyde stock solution to the cell culture media to achieve a final concentration of 1%. Crosslink was quenched for 5 min at room temperature by the addition of glycine to a final concentration of 0.125 M followed by two washes with ice-cold PBS. Cells were scraped, collected by centrifugation and subjected to cell lysis in lysis buffer 1 (10 mM HEPES pH 7.0, 10 mM EDTA, 0.5 mM EGTA, 0.25% Triton X-100, supplemented with 0.5 mM PMSF before use) with rotation for 10 min at 4 °C. The crude nuclear pellets were collected by centrifugation (1500 rpm, 4 °C, 4 min), washed with lysis buffer 2 (10 mM HEPES pH 7.0, 200 mM NaCl, 1 mM EDTA, 0.5 mM EGTA, supplemented with 0.5 mM PMSF before use) with rotation for 10 min at 4 °C. Nuclear pellets were collected by centrifugation (1500 rpm, 4 °C, 4 min), resuspended in nuclear lysis buffer [50 mM Tris–HCl pH 8.1, 10 mM EDTA, 1% SDS, supplemented with 1 mM PMSF and 1 × protease inhibitor cocktail (Sigma-Aldrich) before use] and incubated on ice for 10 min. The chromatin was sheared by sonication in ice-water bath at 4 °C using a Branson Sonifier 450 with a microtip (40% amplitude, 3 s on, 10 s off, 5 min total pulse time). Sonicated nuclear lysates were cleared by centrifugation (15,000 rpm, 10 °C, 15 min) and nuclear protein concentration was determined using the BioRad Protein Assay (BioRad). Equal amounts of total nuclear proteins were used for ChIP and supplemented with corresponding volume of nuclear lysis buffer to achieve same final volumes between different samples, and then diluted 1:10 with dilution buffer (20 mM Tris–HCl pH 8.1, 150 mM NaCl, 2 mM EDTA, 1% Triton X-100, supplemented with 1 × protease inhibitor cocktail before use). Five percent of the material was removed and saved as input, and the rest was incubated overnight at 4 °C with the antibody of interest or a normal IgG control.
The following day, the immune complexes were collected by adding Dynabeads Protein A or Dynabeads M-280 Sheep anti-Mouse IgG (Life Technologies) (beads were pre-washed with ChIP dilution buffer three times before use) and incubating with rotation for 2 h at 4 °C. The immunoprecipitated material was subsequently washed once with low salt wash buffer (20 mM Tris–HCl pH 8.1, 150 mM NaCl, 2 mM EDTA, 0.1% SDS, 1% Triton X-100), once with high salt wash buffer (20 mM Tris–HCl pH 8.1, 500 mM NaCl, 2 mM EDTA, 0.1% SDS, 1% Triton X-100), once with LiCl wash buffer (10 mM Tris–HCl pH 8.1, 0.25 M LiCl, 1 mM EDTA, 1% NP-40, 1% deoxycholate), and twice with TE buffer (10 mM Tris–HCl pH 8.0, 1 mM EDTA pH 8.0). Each wash was done with rotation for 5 min at 4 °C followed by separation on magnetic stand. The immunoprecipitated material was eluted twice at room temperature in freshly prepared elution buffer (1% SDS, 0.1 M NaHCO3) by shaking on a vortexer for 20 min. The eluted material and the input material were then digested with proteinase K (200 μg/ml final concentration) for 2 h at 55 °C, and the crosslinks were reversed overnight by incubating at 65 °C in a hybridization oven. DNA was purified using a Qiagen PCR Purification Kit according to the manufacturer’s protocol.
ChIP-seq library preparation
Prior to ChIP-seq library preparation, ChIPed DNA samples were submitted to the Sequencing Facility at University of Wisconsin Madison Biotechnology Center for determining DNA concentration and size distribution, using Qubit Fluorometer (Thermo Fisher Scientific) and Agilent High Sensitivity DNA Kit (Agilent Technologies), respectively. 10 ng of ChIPed DNA from each condition was used to generate ChIP-seq library using the Ovation Ultralow System V2 1–16 Kit (NuGEN Technologies), according to the manufacturer’s protocol. Briefly, the DNA was end-repaired and subsequently ligated to Illumina sequencing adaptors. The ligated DNA was purified using the Agencourt RNAClean XP bead (Beckman Coulter). A subsequent PCR amplification step (15 cycles) added additional linker sequence to the purified fragments to prepare them for annealing to the Genome Analyzer flow-cell. Following PCR amplification, the library was size selected to a narrow range of fragment sizes by separation on a 2% agarose gel (120 V, 1.5 h) and a band between 300–500 bp was excised. The library was purified from the excised agarose using the Qiagen MiniElute PCR Purification Kit, according to the manufacturer’s protocol. Purified libraries were then submitted to the Sequencing Facility at University of Wisconsin Madison Biotechnology Center for quality control to determine the size, purity, and concentration of the final ChIP-seq libraries. Qualified libraries were submitted to New York Genome Center (NYGC) for deep sequencing using an Illumina HiSeq 2000 per the manufacturer’s instructions.
ChIP-seq data analyses
ChIP-seq raw data were obtained from NYGC, and the reads were mapped to the human reference genome (UCSC hg19) using Bowtie (version 2.1.0)  using default parameters. ChIP-seq raw peaks were visualized in Integrative Genomics Viewer (IGV) (version 2.3.40) . Peak calling was done by Model-based Analysis of ChIP-seq (MACS) (version 1.4.2)  using default parameters and by taking one treatment group as control and comparing the other treatment group to the corresponding control. Coverage was determined by normalizing the total number of mapped reads per hundred million. MEME suite (version 4.10.1) was used for motif discovery. The ±100 bp around a peak summit was used for motif discovery. Annotated genes with its 50 kb upstream that overlapped with ERα peaks were defined as ERα target genes. Annotated genes with their gene body overlapped with RNAPII peaks were defined as RNAPII target genes. For plots of ChIP reads over gene bodies, each gene was divided into 20 intervals (5% each interval) separately for the body of the genes, 2 kb upstream of the TSS, and 2 kb downstream of the TTS. Normalized ChIP read density for each interval was plotted. Standardized z-scores were used in heatmap visualization.
Primers used in ChIP-qPCR experiments
GREB1 TSS F: 5′-GCCAAATGGAAGAAGGACAG-3′
GREB1 TSS R: 5′-ACCACCTACCTCCAGTCACC -3′
GREB1 TR F: 5′-CCCAAGCCTTCTCTGGGTTC-3′
GREB1 TR R: 5′-AGCAGACGAGAAGTAGGGCT-3′
IGFBP4 TSS F: 5′-GGTGTGGCCAGGGACCGGTATAAA-3′
IGFBP4 TSS R: 5′- AAGTTAGCAGGCGTGCCGGA-3′
IGFBP4 TR F: 5′-TTGGGGAGGATGAGGGAGTG-3′
IGFBP4 TR R: 5′-CCTCAACGCCAAAGTCCCTCTA-3′
chromatin immunoprecipitation coupled with high-throughput sequencing
estrogen response element
RNAPII associated factor complex
RNA polymerase II
Marino M, Galluzzo P, Ascenzi P. Estrogen signaling multiple pathways to impact gene transcription. Curr Genomics. 2006;7(8):497–508.
Hall JM, Couse JF, Korach KS. The multifaceted mechanisms of estradiol and estrogen receptor signaling. J Biol Chem. 2001;276(40):36869–72.
Lonard DM, Lanz RB, O’Malley BW. Nuclear receptor coregulators and human disease. Endocr Rev. 2007;28(5):575–87.
Wang L, Yu Y, Chow DC, Yan F, Hsu CC, Stossi F, Mancini MA, Palzkill T, Liao L, Zhou S, et al. Characterization of a steroid receptor coactivator small molecule stimulator that overstimulates cancer cells and leads to cell stress and death. Cancer Cell. 2015;28(2):240–52.
Chaudhary K, Deb S, Moniaux N, Ponnusamy MP, Batra SK. Human RNA polymerase II-associated factor complex: dysregulation in cancer. Oncogene. 2007;26(54):7499–507.
Crisucci EM, Arndt KM. The Roles of the Paf1 Complex and Associated Histone Modifications in Regulating Gene Expression. Genet Res Int. 2011;2011:15. doi:10.4061/2011/707641.
Jaehning JA. The Paf1 complex: platform or player in RNA polymerase II transcription? Biochim Biophys Acta. 2010;1799(5–6):379–88.
Kim J, Guermah M, Roeder RG. The human PAF1 complex acts in chromatin transcription elongation both independently and cooperatively with SII/TFIIS. Cell. 2010;140(4):491–503.
Tomson BN, Arndt KM. The many roles of the conserved eukaryotic Paf1 complex in regulating transcription, histone modifications, and disease states. Biochim Biophys Acta. 2013;1829(1):116–26.
Mueller CL, Porter SE, Hoffman MG, Jaehning JA. The Paf1 complex has functions independent of actively transcribing RNA polymerase II. Mol Cell. 2004;14(4):447–56.
Ding L, Paszkowski-Rogacz M, Nitzsche A, Slabicki MM, Heninger AK, de Vries I, Kittler R, Junqueira M, Shevchenko A, Schulz H, et al. A genome-scale RNAi screen for Oct4 modulators defines a role of the Paf1 complex for embryonic stem cell identity. Cell Stem Cell. 2009;4(5):403–15.
Zhang C, Mejia LA, Huang J, Valnegri P, Bennett EJ, Anckar J, Jahani-Asl A, Gallardo G, Ikeuchi Y, Yamada T, et al. The X-linked intellectual disability protein PHF6 associates with the PAF1 complex and regulates neuronal migration in the mammalian brain. Neuron. 2013;78(6):986–93.
Marazzi I, Ho JS, Kim J, Manicassamy B, Dewell S, Albrecht RA, Seibert CW, Schaefer U, Jeffrey KL, Prinjha RK, et al. Suppression of the antiviral response by an influenza histone mimic. Nature. 2012;483(7390):428–33.
Muntean AG, Tan J, Sitwala K, Huang Y, Bronstein J, Connelly JA, Basrur V, Elenitoba-Johnson KS, Hess JL. The PAF complex synergizes with MLL fusion proteins at HOX loci to promote leukemogenesis. Cancer Cell. 2010;17(6):609–21.
Tan J, Muntean AG, Hess JL. PAFc, a key player in MLL-rearranged leukemogenesis. Oncotarget. 2010;1(6):461–5.
Zeng H, Xu W. Ctr9, a key subunit of PAFc, affects global estrogen signaling and drives ERα-positive breast tumorigenesis. Genes Dev. 2015;29(20):2153–67.
Xue X, Yang YA, Zhang A, Fong KW, Kim J, Song B, Li S, Zhao JC, Yu J. LncRNA HOTAIR enhances ER signaling and confers tamoxifen resistance in breast cancer. Oncogene. 2016;35:2746–55.
Nagarajan S, Hossan T, Alawi M, Najafova Z, Indenbirken D, Bedi U, Taipaleenmäki H, Ben-Batalla I, Scheller M, Loges S, et al. Bromodomain protein BRD4 is required for estrogen receptor-dependent enhancer activation and gene transcription. Cell Rep. 2014;8(2):460–9.
Mohammed H, D’Santos C, Serandour AA, Ali HR, Brown GD, Atkins A, Rueda OM, Holmes KA, Theodorou V, Robinson JL, et al. Endogenous purification reveals GREB1 as a key estrogen receptor regulatory factor. Cell Rep. 2013;3(2):342–9.
Kininis M, Isaacs GD, Core LJ, Hah N, Kraus WL. Postrecruitment regulation of RNA polymerase II directs rapid signaling responses at the promoters of estrogen target genes. Mol Cell Biol. 2009;29(5):1123–33.
Kininis M, Chen BS, Diehl AG, Isaacs GD, Zhang T, Siepel AC, Clark AG, Kraus WL. Genomic analyses of transcription factor binding, histone acetylation, and gene expression reveal mechanistically distinct classes of estrogen-regulated promoters. Mol Cell Biol. 2007;27(14):5090–104.
Kwon YS, Garcia-Bassets I, Hutt KR, Cheng CS, Jin M, Liu D, Benner C, Wang D, Ye Z, Bibikova M, et al. Sensitive ChIP-DSL technology reveals an extensive estrogen receptor alpha-binding program on human gene promoters. Proc Natl Acad Sci U S A. 2007;104(12):4852–7.
Lin CY, Vega VB, Thomsen JS, Zhang T, Kong SL, Xie M, Chiu KP, Lipovich L, Barnett DH, Stossi F, et al. Whole-genome cartography of estrogen receptor alpha binding sites. PLoS Genet. 2007;3(6):e87.
Stender JD, Frasor J, Komm B, Chang KC, Kraus WL, Katzenellenbogen BS. Estrogen-regulated gene networks in human breast cancer cells: involvement of E2F1 in the regulation of cell proliferation. Mol Endocrinol. 2007;21(9):2112–23.
Carroll JS, Meyer CA, Song J, Li W, Geistlinger TR, Eeckhoute J, Brodsky AS, Keeton EK, Fertuck KC, Hall GF, et al. Genome-wide analysis of estrogen receptor binding sites. Nat Genet. 2006;38(11):1289–97.
Coser KR, Chesnes J, Hur J, Ray S, Isselbacher KJ, Shioda T. Global analysis of ligand sensitivity of estrogen inducible and suppressible genes in MCF7/BUS breast cancer cells by DNA microarray. Proc Natl Acad Sci U S A. 2003;100(24):13994–9.
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(10):4562–74.
Charpentier AH, Bednarek AK, Daniel RL, Hawkins KA, Laflin KJ, Gaddis S, MacLeod MC, Aldaz CM. Effects of estrogen on global gene expression: identification of novel targets of estrogen action. Cancer Res. 2000;60(21):5977–83.
Zeng H, Xu W. Gene expression profiling of Ctr9-regulated transcriptome in ERα-positive breast cancer. Genom Data. 2016;7:103–4.
Franco HL, Nagari A, Kraus WL. TNFα signaling exposes latent estrogen receptor binding sites to alter the breast cancer cell transcriptome. Mol Cell. 2015;58(1):21–34.
Mohammed H, Russell IA, Stark R, Rueda OM, Hickey TE, Tarulli GA, Serandour AA, Birrell SN, Bruna A, Saadi A, et al. Progesterone receptor modulates ERα action in breast cancer. Nature. 2015;523(7560):313–7.
Welboren WJ, van Driel MA, Janssen-Megens EM, van Heeringen SJ, Sweep FC, Span PN, Stunnenberg HG. ChIP-Seq of ERalpha and RNA polymerase II defines genes differentially responding to ligands. EMBO J. 2009;28(10):1418–28.
Wu J, Xu W. Histone H3R17me2a mark recruits human RNA polymerase-associated factor 1 complex to activate transcription. Proc Natl Acad Sci U S A. 2012;109(15):5675–80.
Liu T, Ortiz JA, Taing L, Meyer CA, Lee B, Zhang Y, Shin H, Wong SS, Ma J, Lei Y, et al. Cistrome: an integrative platform for transcriptional regulation studies. Genome Biol. 2011;12(8):R83.
Bailey TL, Boden M, Buske FA, Frith M, Grant CE, Clementi L, Ren J, Li WW, Noble WS. MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res. 2009;37:W202–8.
Hurtado A, Holmes KA, Ross-Innes CS, Schmidt D, Carroll JS. FOXA1 is a key determinant of estrogen receptor function and endocrine response. Nat Genet. 2011;43(1):27–33.
Ringnér M, Fredlund E, Häkkinen J, Borg Å, Staaf J. GOBO: gene expression-based outcome for breast cancer online. PLoS One. 2011;6(3):e17911.
Chen FX, Woodfin AR, Gardini A, Rickels RA, Marshall SA, Smith ER, Shiekhattar R, Shilatifard A. PAF1, a molecular regulator of promoter-proximal pausing by RNA polymerase II. Cell. 2015;162(5):1003–15.
Theodorou V, Stark R, Menon S, Carroll JS. GATA3 acts upstream of FOXA1 in mediating ESR1 binding by shaping enhancer accessibility. Genome Res. 2013;23(1):12–22.
Wade MA, Jones D, Wilson L, Stockley J, Coffey K, Robson CN, Gaughan L. The histone demethylase enzyme KDM3A is a key estrogen receptor regulator in breast cancer. Nucleic Acids Res. 2015;43(1):196–207.
Antony J, Dasgupta T, Rhodes JM, McEwan MV, Print CG, O’Sullivan JM, Horsfield JA. Cohesin modulates transcription of estrogen-responsive genes. Biochim Biophys Acta. 2015;1849(3):257–69.
Prenzel T, Kramer F, Bedi U, Nagarajan S, Beissbarth T, Johnsen SA. Cohesin is required for expression of the estrogen receptor-alpha (ESR1) gene. Epigenetics Chromatin. 2012;5(1):13.
Yu M, Yang W, Ni T, Tang Z, Nakadai T, Zhu J, Roeder RG. RNA polymerase II-associated factor 1 regulates the release and phosphorylation of paused RNA polymerase II. Science. 2015;350(6266):1383–6.
Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10(3):R25.
Thorvaldsdóttir H, Robinson JT, Mesirov JP. Integrative genomics viewer (IGV): high-performance genomics data visualization and exploration. Brief Bioinform. 2013;14(2):178–92.
Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, Nusbaum C, Myers RM, Brown M, Li W, et al. Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008;9(9):R137.
This project was supported by a Department of Defense ERA of Hope Award W81XWYH-11-1-0237 to W.X. Work in X.Z.’s laboratory was supported by startup funds from UW-Madison, an NSF CAREER award (MCB-1552455), Alexander von Humboldt foundation (Alfred Toepfer Faculty Fellow), and USDA & National Institute of Food and Agriculture grant (Hatch 1002874). P.A. was supported by NIH grant CA22443. P.A. is an investigator of the Howard Hughes Medical Institute and the Morgridge Institute for Research.
Availability of data and materials
The data sets supporting the conclusions of this article (ChIP-seq data of ERα and RNAPII) are available in the NCBI’s GEO database with the accession number GSE80728 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE80728). Publically available RNAPII ChIP-seq data sets in wild-type or PAF1 knockdown MCF7 cells were downloaded from GSE70408 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE70408).
HZ, PA and WX conceived and designed the research. HZ and NTC performed the experiments. HZ, LL and XZ performed the bioinformatics analyses and edited the manuscript. HZ and WX wrote the manuscript. All authors provided intellectual input and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate