Skip to main content

Pho dynamically interacts with Spt5 to facilitate transcriptional switches at the hsp70 locus

Abstract

Background

Numerous target genes of the Polycomb group (PcG) are transiently activated by a stimulus and subsequently repressed. However, mechanisms by which PcG proteins regulate such target genes remain elusive.

Results

We employed the heat shock-responsive hsp70 locus in Drosophila to study the chromatin dynamics of PRC1 and its interplay with known regulators of the locus before, during and after heat shock. We detected mutually exclusive binding patterns for HSF and PRC1 at the hsp70 locus. We found that Pleiohomeotic (Pho), a DNA-binding PcG member, dynamically interacts with Spt5, an elongation factor. The dynamic interaction switch between Pho and Spt5 is triggered by the recruitment of HSF to chromatin. Mutation in the protein–protein interaction domain (REPO domain) of Pho interferes with the dynamics of its interaction with Spt5. The transcriptional kinetics of the heat shock response is negatively affected by a mutation in the REPO domain of Pho.

Conclusions

We propose that a dynamic interaction switch between PcG proteins and an elongation factor enables stress-inducible genes to efficiently switch between ON/OFF states in the presence/absence of the activating stimulus.

Background

Polycomb group (PcG) proteins are a highly conserved class of epigenetic regulators that play an essential role during development and adulthood in many organisms [1]. In any given cell type, while PcG proteins maintain stable silencing at a certain subset of genes, other PcG target genes are capable of being activated in the presence of the appropriate stimulus. Thus, by conferring both stability and plasticity upon the transcriptome, PcG proteins bestow competence upon the cells to respond to varied stimuli over their lifetime, in order to maintain cellular homoeostasis [2]. However, the underlying mechanisms that enable PcG-regulated loci to dynamically switch between ON and OFF states remain elusive.

In Drosophila, the targeting of PcG proteins to chromatin occurs mainly through the recognition of cis-regulatory elements, known as Polycomb group response elements (PREs), by Pleiohomeotic (Pho) and other DNA-binding proteins [3, 4]. Pho interacts with dSfmbt to form the Pho repressive complex (PhoRC), which acts as a recruitment platform for the Polycomb repressive complex 1 and 2 (PRC1 and PRC2) [5, 6]. Subsequently, PRC1 and PRC2 generate a chromatin structure that is incompatible with gene expression, through their enzymatic and chromatin compaction activities [7, 8]. The transcriptionally silent state is then faithfully inherited through cell division, which involves PRE-dependent mechanisms [9, 10]. Although the mechanisms by which PcG proteins establish and maintain stable silencing at their targets have been subject to intense investigation, little is known about PcG-mediated regulation of target genes with dynamic expression patterns.

In this study, we employed the heat shock-responsive hsp70 gene in Drosophila as a model to study dynamic gene regulation by PcG proteins (Additional file 1: Figure S1a). The hsp70 locus is an attractive inducible gene model as it is targeted by PRC1 [11]. It is regulated at the level of transcriptional elongation by promoter-proximal RNA polymerase II pausing [12]. At normal growth temperature (25 °C), RNA polymerase II is maintained in the paused state by the DSIF complex (composed of Spt4–Spt5) and the NELF complex [13, 14]. The heat shock stimulus leads to the recruitment of the transcription factor, HSF, to chromatin [15]. This is followed by the recruitment of P-TEFb, which phosphorylates Spt5 and converts it to function as a positive elongation factor [16, 17]. This leads to rapid up-regulation of the hsp70 locus in response to the heat stress [18].

While the transcriptional response and the chromatin dynamics at the hsp70 locus have been well characterised in the acute phase of the heat shock response, the effect of the withdrawal of the activating stimulus on the transcription and chromatin dynamics remains poorly understood [18]. By following the chromatin dynamics of PRC1 and Pho during the recovery phase of the heat shock, we sought to dissect the role of these two factors in the re-establishment of silencing at the hsp70 locus. Pho was previously found to co-localise with RNA polymerase II upon heat shock [19]. Furthermore, Pho interacts with Spt5 [20], and Spt5 is involved in both, promoter-proximal pausing and transcriptional elongation [21]. We thus asked whether the conversion of Spt5 from a pausing factor into an elongation factor upon gene activation and its interaction with Pho could mechanistically explain the spreading of Pho in the body of active genes. Genetic analyses revealed the necessity of Pho gene function in the re-establishment of silencing at the hsp70 locus in the absence of the activating stimulus [19]. We tested the hypothesis whether Pho could dynamically interact with repressors (PRC1 via dSfmbt) and activators (Spt5) during the heat shock response and the significance of its dynamic protein–protein interactions on transcriptional regulation of the hsp70 locus.

By combining multiple read-outs, we identified a role for the dynamic protein–protein interaction between Pho and Spt5 in the transcriptional regulation of the hsp70 locus. Our results highlight how a dynamic interaction switch between PcG proteins and an elongation factor could facilitate PcG-regulated genes, which are responsive to inductive signals, to efficiently switch between ON/OFF states.

Results

Antagonistic binding patterns of HSF and PRC1 at the hsp70 locus

To characterise the transcriptional and chromatin dynamics of the heat shock response, we focused on the binding kinetics of RNA polymerase II (Pol II), HSF and Polyhomeotic (Ph, a component of PRC1) at the hsp70 locus using Drosophila S2 tissue culture cells. Our analysis revealed a rapid increase in hsp70 transcripts immediately upon heat shock, followed by a steady decrease in mRNA levels as the locus is re-repressed, during the recovery phase in the absence of the heat shock stimulus (Fig. 1b). We investigated the binding dynamics of two forms of Pol II, the unphosphorylated C-terminal domain (CTD) form and the hyperphosphorylated (S2P) form of Pol II, at the promoter-proximal region and within the gene body of the hsp70 locus, respectively. We detected a rapid increase in both forms of Pol II at the promoter-proximal region and within the gene body of the hsp70 locus upon heat shock (Fig. 1c, d). This was accompanied by a decrease in Pol II binding levels at the Act42A gene (Additional file 2: Figure S2a, b), which is in accordance with the previously described phenomenon of the global decrease in Pol II levels genome-wide in Drosophila upon heat shock [22]. However, levels of bound Pol II at the hsp70 locus steadily decreased as the cells recovered from the heat shock, with the locus returning to the paused state 90 min into the recovery phase (Fig. 1c, d). In contrast, Pol II binding at the Act42A gene increased over the recovery phase, indicating that the cells resumed a normal transcriptional programme within 90 min of the heat shock stimulus (Additional file 2: Figure S2a, b).

Fig. 1
figure 1

Antagonistic binding patterns of HSF and PRC1 at the hsp70 locus. a Experimental scheme used throughout the study to measure the dynamics of mRNA/nascent RNA production, chromatin binding patterns of proteins and protein–protein interactions during the heat shock response. S2 DRSC cells were heat shocked at 37 °C for 15 min and thereafter allowed to recover from the heat shock at 25 °C for up to 90 min. At time points indicated in the scheme, cells were harvested and used for the following experimental procedures: mRNA quantification, ChIP-qPCR, co-IP and nascent RNA quantification. b Quantitative RT-PCR (qRT-PCR) detection of hsp70 transcript levels during the heat shock response detailed in a. Distance of the location of the primers used from the hsp70 transcription start site (TSS) is + 687 bp and has been depicted in the form of a cartoon below the data figure. The control line represents cells that were maintained at 25 °C for the entire duration of the time course detailed in a. cf ChIP-qPCR measurements of occupancy levels of RNA polymerase II CTD, S2P form of RNA polymerase II, heat shock binding factor (HSF) and Polyhomeotic (Ph, a component of PRC1) at the hsp70 locus during the heat shock response detailed in a. Distances of the location of the primers used from the hsp70 TSS are as follows: for c and f (+ 78 bp), for d (+ 802 bp) and for e (− 85 bp) has been depicted in the form of a cartoon below the data figure. The control line represents cells that were maintained at 25 °C for the entire duration of the time course detailed in a. Data information: In bf, data are presented as mean ± SEM for n = 2 (except for f where n = 3). Statistical significance was determined by performing a paired two-tailed t test, **P ≤ 0.01, *P ≤ 0.05 and ns = non-significant

HSF binding is critical to release paused Pol II into productive elongation at the hsp70 locus [23]. In agreement, we detected a rapid increase in HSF binding upon heat shock followed by steady decrease as the cells recovered from the heat shock (Fig. 1e). Thus, the chromatin binding dynamic of Pol II and HSF during the heat shock response was similar, emphasising a role for HSF binding at chromatin for the transcriptional activation of the hsp70 locus [23].

Recent ChIP-seq studies in Drosophila and mammalian cell lines have revealed PRC1 binding at genes regulated by paused Pol II, leading to the hypothesis that PRC1 could regulate paused genes [11, 24, 25]. We, therefore, assessed the binding dynamics of Ph during the heat shock response. In keeping with PRC1’s role as a repressor, Ph occupancy at the hsp70 locus decreased significantly upon activation. However, as the locus recovered from the heat shock and occupancy of Pol II and HSF decreased, Ph binding to chromatin was steadily restored to pre-activation levels (Fig. 1f).

Taken together, these results reveal an antagonism between HSF and Ph (PRC1) for occupancy at the hsp70 locus.

Co-localisation of Pleiohomeotic (Pho) and Pol II at the heat shock locus is independent of transcriptional elongation

We confirmed the previously identified co-localisation of Pho with Pol II at active genes by co-staining these factors on polytene chromosomes from larval salivary glands. Pho and Pol II co-localised at the majority of bands under basal conditions (Fig. 2a). Upon heat shock, Pol II binding was most prominent at the heat shock puffs (Fig. 2b). In addition, Pol II signal intensity decreased outside the heat shock loci, in accordance with a genome-wide decrease in Pol II occupancy upon heat shock (Fig. 2b). Interestingly, we observed a striking co-localisation of Pho and Pol II at all activated heat shock genes (Fig. 2b).

Fig. 2
figure 2

Co-localisation of Pleiohomeotic (Pho) and RNA polymerase II at the heat shock loci upon activation is independent of transcriptional elongation. ad Double polytene immunostaining for RNA polymerase II (green) and Pho (magenta) under control, HS (heat shock), FP (flavopiridol) control and FP HS (flavopiridol + heat shock), respectively. The white block arrows point at the double heat shock puffs (87A/C) that were observed under HS and FP HS conditions. Magnification used 63×. e, f ChIP-qPCR measurements of occupancy levels for HSF and Ph, respectively, at the hsp70 locus in S2 DRSC cells. Distances of the location of the primers used from the hsp70 TSS are as follows: for e (− 85 bp) and for f (+ 78 bp) has been depicted in the form of a cartoon below the data figure. The control bar represents cells that were maintained at 25 °C for 15 min, the HS bar represents cells that were maintained at 37 °C for 15 min, the FP control bar represents cells that were treated with 500 nM flavopiridol for 40 min and then maintained at 25 °C for 15 min and the FP HS bar represents cells that were treated with 500 nM flavopiridol for 40 min and then maintained at 37 °C for 15 min. Data information: In e, f, data are presented as mean ± SEM for n = 2

To determine whether this co-localisation is dependent upon transcription, we blocked transcriptional elongation by treatment with flavopiridol, a P-TEFb inhibitor, prior to heat shock [26]. The striking overlap in binding patterns of Pho and Pol II was not perturbed by flavopiridol treatment, suggesting that the interaction between these two factors could occur independently of transcriptional elongation (Fig. 2c). Furthermore, inhibition of P-TEFb followed by heat shock resulted in a marked decrease in the size of the heat shock puffs (Fig. 2d) and strongly reduced hsp70 transcript levels (Additional file 3: Figure S3d). However, co-localisation of Pho and Pol II was not perturbed by flavopiridol (Fig. 2d). This result indicated that the co-localisation between Pho and Pol II was intact during the paused state as well as during the elongation phase of the transcriptional cycle. It is conceivable that Pho might co-localise with the paused form of Pol II. However, we found it intriguing that a DNA-binding protein like Pho was also capable of co-localising with actively elongating Pol II at chromatin. In order to understand whether the interaction between Pho and elongating Pol II was dependent upon the DNA-binding activity of Pho, we employed native ChIP (N-ChIP) qPCR to measure the occupancy levels of Pho over the time course detailed in Fig. 1a. We employed N-ChIP over formaldehyde cross-linked ChIP to measure the occupancy levels of Pho since N-ChIP only detects proteins that are associated with chromatin due to their direct interaction with DNA [27]. Indeed, while Pho was bound at the hsp70 locus under control conditions, its occupancy decreased upon heat shock. However, as the cells recovered from the heat shock, the binding of Pho to the hsp70 locus via DNA was restored (Additional file 3: Figure S3a). This suggested that the association of Pho to chromatin upon heat shock cannot be explained by its DNA-binding ability and could possibly be due to a protein–protein interaction with RNA Pol II or an elongation factor.

Previous observations established that flavopiridol treatment leads to considerable decrease in mature hsp70 transcript production, even though the density of Pol II in the gene body remains normal [28]. We have now shown that chromatin binding levels of Ph, a component of PRC1, anti-correlate with mature hsp70 transcript levels (Fig. 1f). Thus, we hypothesised that HSF recruitment in the presence of flavopiridol could be used to uncouple its effect on the occupancy of Ph at chromatin from the production of mature hsp70 transcripts.

We, therefore, examined the consequence of flavopiridol treatment on the binding of PRC1 at the hsp70 locus. We treated cells with flavopiridol prior to heat shock (hereafter referred to as FP + HS) followed by ChIP-qPCR. In agreement with previous observations, flavopiridol treatment led to a significant decrease in S2P Pol II levels in the hsp70 gene body upon heat shock (Additional file 3: Figure S3c), whereas unphosphorylated Pol II and HSF binding levels were not perturbed (Additional file 3: Figure S3b and Fig. 2e). Ph binding at the hsp70 locus decreased under FP + HS conditions, once again demonstrating that HSF and PRC1 binding patterns are mutually exclusive at the hsp70 locus (Fig. 2e, f).

These data revealed that the co-localisation between Pho and Pol II was intact during the paused as well as during the elongation phase of the transcriptional cycle. Interestingly, while Pho continued to localise at the hsp70 locus upon heat shock, Ph (PRC1) was evicted from the locus. We, thus, concluded that HSF binding could be capable of disrupting the ability of Pho to recruit PRC1 at the hsp70 locus.

Pho dynamically interacts with Spt5 and dSfmbt

Since we determined that the co-localisation of Pho and Pol II at the hsp70 locus upon activation was independent of the DNA-binding ability of Pho, we hypothesised that this observation could possibly be explained by the reported genetic and biochemical interaction between Pho and Spt5 [20]. We also determined by ChIP-qPCR that chromatin occupancy of Spt5 does not change dramatically over the time course detailed in Fig. 1a (Additional file 4: Figure S4a). In order to gain mechanistic insight into the significance of the interaction between Pho and Spt5, we sought to identify the amino acid residues in Pho that mediate the protein–protein interaction with Spt5. Pho exhibits a modular structure with its protein–protein interaction domains localised in the N-terminus and its DNA-binding property ascribed to the Zn-finger domains in the C-terminus [5, 29, 30]. To identify putative Spt5-interacting regions of Pho, we searched for evolutionary conserved functional domain(s) by performing a multiple protein alignment between Pho and YY1, the mouse and human homologue of Pho. The Zn-finger domains in Pho and YY1 showed a high degree of conservation in the multiple protein alignment (Fig. 3a). However, a second stretch of amino acids, located in the N-terminus of Pho, also showed conservation with YY1 in the multiple protein alignment (Fig. 3a). This region mapped to a domain known as the Pho spacer region [5], and its homologous region within YY1 is known as the Recruitment of Polycomb (REPO) as it was shown to be necessary and sufficient for PcG-mediated repression in Drosophila [31].

Fig. 3
figure 3

Pho interacts dynamically with Spt5 and dSfmbt. a Multiple protein alignment between Pho and YY1, the mouse and human homologue of Pho. The REPO domain is surrounded by a blue rectangle, while the Zn-finger domains are surrounded by an orange rectangle. The asterisk mark on the bottom line ‘consensus’ denotes the amino acid residues that are fully conserved between the three species. b Representation of the N- and C-termini of Pho and the mutant constructs generated for this study. c, d co-IP assays of S2 DRSC cells transiently transfected with plasmids expressing FLAG-tagged dSfmbt or Spt5 (FLAG-dSfmbt, FLAG-Spt5) and HA-tagged Pho, PhoV164D, Pho ΔREPO, EGFP (HA-Pho, HA-PhoV164D, HA-Pho REPO del, HA-EGFP) as indicated. Cell lysates were used for pull-downs using an anti-FLAG antibody and were later probed by Western blot using an anti-HA antibody. MW = molecular weight in kDa. e co-IP assays of S2 DRSC cells transiently transfected with plasmids expressing FLAG-tagged dSfmbt or Spt5 (FLAG-dSfmbt, FLAG-Spt5) and HA-tagged Pho or Spt4 (HA-Pho, HA-Spt4) as indicated. S2 DRSC cells were either maintained at normal growth temperature (25 °C) or heat shocked at 37 °C for 15 min, cell lysates were used for pull-downs using an anti-FLAG antibody and were later probed by Western blot using an anti-HA antibody. MW = molecular weight in kDa

We hypothesised that the REPO domain of Pho could be involved in mediating its protein–protein interaction with Spt5. To test this hypothesis, we generated two mutant constructs with the aim of disrupting the structure of the REPO domain [5] and analysed the effect of the mutation on the protein–protein interactions of Pho by co-immunoprecipitation (co-IP). The first construct, PhoV164D, contained a point mutation resulting in a valine (a hydrophobic amino acid residue) to aspartic acid (a hydrophilic amino acid residue) substitution to interfere with its ability to contact hydrophobic surfaces, while the second construct, ΔREPO, carried a 66-nucleotide deletion in the REPO domain [5] (Fig. 3b).

In agreement with previous findings, both mutant Pho constructs were incapable of interacting with dSfmbt [5] (Fig. 3c). However, both mutant Pho constructs also significantly impaired the interaction between Pho and Spt5, indicating that the REPO domain represents an interaction surface with multiple protein interaction partners (Fig. 3d). Next, we sought to explain the residual interaction between Pho and Spt5 in the mutant background using different N-terminal and C-terminal Pho constructs (1-351a.a N-Pho, 1-200a.a N-Pho, 173-351a.a N-Pho, 351-520a.a C-Pho). Our analysis revealed that in addition to the REPO domain containing N-terminal fragments of Pho, the C-terminus of Pho could also interact with Spt5, thus explaining the residual interaction in the absence of a functional REPO domain (Additional file 4: Figure S4c). In contrast, the Pho–dSfmbt interaction was solely dependent on the REPO domain (Additional file 4: Figure S4b). In conclusion, the REPO domain facilitated the interaction of Pho with dSfmbt (a PcG protein) as well as with Spt5 (a transcription elongation factor). These data led us to hypothesise that the protein–protein interactions of Pho might be dynamic under ON/OFF transcriptional states.

To test this hypothesis, we performed co-IPs of Pho–Spt5 and Pho–dSfmbt under control and heat shock conditions. Upon heat shock, we observed a considerable increase in the interaction between Pho and Spt5 and a concomitant reduction in the interaction between Pho and dSfmbt, which was stronger in non-stimulated cells (Fig. 3e and Additional file 4: Figure S4d). As a negative control, we tested whether the interaction between Spt4 and Spt5 is also dynamic upon heat shock. However, the interaction between Spt4 and Spt5 appeared stable under control and heat shock conditions (Fig. 3e). Thus, our co-IP experiments support the view that Pho can dynamically interact with Spt5 and dSfmbt depending upon the transcriptional state of the hsp70 locus.

Next, we tested whether the dynamic protein–protein interactions of Pho recapitulated the dynamic chromatin binding patterns of activators and repressors during the heat shock response. To this end, we analysed Pho–Spt5 and Pho–dSfmbt interactions during the heat shock response. Interestingly, the Pho–Spt5 interaction followed the Pol II and HSF (activator) binding dynamic. In contrast, the interaction between Pho and dSfmbt closely resembled the dynamic exhibited by Ph (PRC1, repressor) (Fig. 4a–d). This result suggested the possibility of a relationship between the chromatin binding dynamic of HSF and the dynamic protein–protein interaction of Pho with an activator (Spt5).

Fig. 4
figure 4

Pho preferentially interacts with Spt5 and dSfmbt during the active and silencing phase, respectively, of the heat shock response. a, c co-IP assays of S2 DRSC cells transiently transfected with plasmids expressing FLAG-tagged dSfmbt or Spt5 (FLAG-dSfmbt, FLAG-Spt5) and HA-tagged Pho (HA-Pho). The S2 DRSC cells were subjected to a heat shock and then allowed to recover from the heat shock according to the experimental scheme depicted in Fig. 1a. At the indicated time points in the experimental scheme, cell lysates were prepared and were used for pull-downs using an anti-FLAG antibody. They were later probed by Western blot using an anti-HA antibody. b, d Quantification of the co-IP assays in A and C. The control line represents cells that were maintained at 25 °C for the entire duration of the time course detailed in Fig. 1a. Details of the procedure used for quantification can be found in Materials and Methods. Data information: In b and d, data are presented as mean ± SEM for n = 2

HSF recruitment to chromatin is necessary and sufficient for the dynamic interaction between Pho and Spt5

Since our data suggested that HSF binding kinetics at the hsp70 locus could influence the protein–protein interaction between Pho and Spt5, we asked whether HSF recruitment to chromatin would be sufficient and necessary for this dynamic protein–protein interaction. We, thus, relied on heat shock-independent ways to recruit HSF to chromatin.

HSF is canonically chaperoned by Hsp90, which prevents it from trimerising and assuming its active conformation [32]. However, Hsp90 inhibition by radicicol triggers HSF to acquire its active confirmation through trimerisation, thus binding to its cognate site in the hsp70 locus and causing transcriptional activation. We confirmed both recruitment of HSF to the hsp70 locus by ChIP-qPCR and transcriptional activation of the locus upon radicicol treatment (Fig. 5a, b). Interestingly, the interaction between Pho and Spt5 was enhanced by Hsp90 inhibition (Fig. 5c, d). Therefore, we set out to decouple this increase in interaction from the transcriptional activation of the hsp70 locus by recruiting HSF under conditions that caused little or no transcription from the hsp70 locus.

Fig. 5
figure 5

HSF recruitment to chromatin is necessary and sufficient for the dynamic protein–protein interaction between Pho and Spt5. a ChIP-qPCR measurements for the occupancy levels of HSF at the hsp70 locus in S2 DRSC cells upon HS (heat shock), radicicol, FP + HS (flavopiridol + heat shock) and Na Sal (Sodium Salicylate) treatment conditions, respectively. In case of radicicol treatment, the control condition was DMSO treatment. For FP + HS, the control condition was FP control and for Na Sal, the control condition was nuclease-free water. Distance of the location of the primers used from the hsp70 TSS is − 85 bp and has been depicted in the form of a cartoon below the data figure. b qRT-PCR measurements of hsp70 transcript levels in S2 DRSC cells under HS, radicicol, FP + HS and Na Sal treatment conditions, respectively. Distance of the location of the primers used from the hsp70 TSS is + 1492 bp and has been depicted in the form of a cartoon below the data figure. c co-IP assays of S2 DRSC cells transiently transfected with plasmids expressing FLAG-tagged Spt5 (FLAG-Spt5) and HA-tagged Pho (HA-Pho). The S2 DRSC cells were subjected to HS, radicicol, FP + HS and Na Sal, respectively; cell lysates were used for pull-downs with an anti-FLAG antibody and later probed by Western blot with an anti-HA antibody. d Quantification of the co-IP assays in c. Details of the procedure used for quantification can be found in Materials and Methods. e co-IP assays of S2 DRSC cells transiently transfected with plasmids expressing FLAG-tagged Spt5 (FLAG-Spt5) and HA-tagged Pho (HA-Pho). Cells were either treated with HSF or LacZ RNAi for 4 days, then either maintained at control conditions or heat shocked, cell lysates were prepared and used for pull-downs using an anti-FLAG antibody and later probed by Western blot using an anti-HA antibody. f qRT-PCR measurements of hsp70 transcript levels in S2 DRSC cells under HSF and LacZ RNAi conditions. Distance of the location of the primers used from the hsp70 TSS is + 1492 bp and has been depicted in the form of a cartoon below the data figure. Data information: In a and b, data are presented as mean ± SEM for n = 2, while in e, data are presented as mean ± SEM for n = 3

Upon flavopiridol treatment followed by heat shock (FP + HS), HSF was recruited to the hsp70 locus, but inhibition of P-TEFb led to a significant decrease in hsp70 transcript levels (Fig. 5a, b). However, under FP + HS conditions, the interaction between Pho and Spt5 was still enhanced (Fig. 5c, d). Treatment with sodium salicylate has also been shown to recruit HSF to chromatin without affecting transcription from the hsp70 locus [33]. While HSF was robustly recruited to chromatin upon sodium salicylate treatment, hsp70 transcript levels remained unchanged when compared to the control condition (Fig. 5a, b). Upon sodium salicylate treatment, the interaction between Pho and Spt5 was, also, increased (Fig. 5c, d). We could, thus, uncouple the increase in interaction between Pho and Spt5 upon HSF recruitment from transcriptional activation. Next, we tested whether HSF binding was also necessary for this interaction switch between Pho and Spt5 by RNAi-mediated knockdown of HSF. By knocking down HSF, the dynamic nature of the Pho–Spt5 protein–protein interaction was undetectable when compared to the control LacZ RNAi condition (Fig. 5e, f).

Thus, by performing gain-of-function and loss-of-function experiments, we could establish that recruitment of HSF to chromatin is necessary and sufficient for the dynamic protein–protein interaction switch between Pho and Spt5.

In order to confirm that the protein–protein interaction between Pho and Spt5 was independent of the DNA-binding ability of Pho, we performed a co-IP between Pho and Spt5 under DNase I conditions. Indeed, we could confirm that the interaction as well as the dynamic switch between Pho and Spt5 was independent of DNA (Additional file 5: Figure S5a). We further went on to confirm the DNA-independent nature of the Pho–Spt5 interaction by performing a co-IP between N-Pho (lacking the DNA-binding Zn-finger domains) and Spt5 (Additional file 5: Figure S5b). However, we did observe dependence for the dynamic interaction switch between Pho and Spt5 on RNA binding (Additional file 5: Figure S5a). This observation is in line with our conclusion that HSF is necessary and sufficient for the dynamic interaction switch between Pho and Spt5 since HSF binding can trigger release of paused RNA Pol II leading to the production of nascent RNA transcripts [23]. Furthermore, the KOW4-5 domain of Spt5 has been shown to be capable of binding RNA and thereby influence the position of paused RNA Pol II [34].

Absence of a functional REPO domain disrupts the transcriptional kinetics of the hsp70 locus during the heat shock response

We next tested the effect of mutating the REPO domain on the dynamic protein–protein interaction between Pho and Spt5 during the heat shock response. To this end, we performed a co-IP for Pho V164D–Spt5 over the time course detailed in Fig. 1a and compared it to the Pho WT–Spt5 protein–protein interaction dynamic.

To our surprise, mutant and wild-type Pho exhibited comparable interaction strengths with Spt5 upon heat shock, indicating that the dynamic interaction switch upon heat shock is not affected by the V164D point mutation in the REPO domain (Fig. 6a, b). However, while the interaction between Pho WT–Spt5 diminished as the hsp70 locus was repressed in the absence of the activating stimulus, the increase in interaction between Pho V164D–Spt5 remained sustained during the recovery phase (Fig. 6a, b). This could be explained by the fact that in the absence of a functional REPO domain, dSfmbt could no longer compete with Spt5 for interaction with Pho during the recovery phase of the heat shock response (Fig. 6a, b). Similar to the Pho WT–Spt5 interaction dynamic, the Pho V164D–Spt5 interaction switch was independent of DNA but dependent upon RNA (Additional file 6: Figure S6a).

Fig. 6
figure 6

Absence of a functional REPO domain disrupts the transcriptional kinetics of the hsp70 locus during the heat shock response. a co-IP assays of S2 DRSC cells transiently transfected with plasmids expressing FLAG-tagged Spt5 (FLAG-Spt5) and HA-tagged PhoV164D (HA-PhoV164D). The S2 DRSC cells were subjected to a heat shock and then allowed to recover from the heat shock according to the experimental scheme depicted in Fig. 1a. At the indicated time points in the experimental scheme, cell lysates were prepared and were used for pull-downs using an anti-FLAG antibody. They were later probed by Western blot using an anti-HA antibody. b Quantification of the co-IP assay in a and its comparison with the quantification of the Pho–Spt5 co-IP assay across the same time course. The control lines represent cells that were maintained at 25 °C for the entire duration of the time course detailed in Fig. 1a. Details of the procedure used for quantification can be found in Materials and Methods. c qRT-PCR measurements of nascent hsp70 transcripts produced over the course of the heat shock response from S2 DRSC cells upon over-expression of WT Pho and mutant Pho (PhoV164D), respectively. Distance of location of the primers from the hsp70 TSS is + 687 bp and has been depicted in the form of a cartoon below the data figure. Data information: In b, data are presented as mean ± SEM for n = 2, while in c, data are presented as mean ± SEM for n = 3

To directly assess the effect of the PhoV164D mutant on transcription from the hsp70 locus upon heat shock, we expressed the PhoV164D mutant by placing it under the control of the copper-inducible metallothionein promoter. As a control, PhoWT was expressed under the control of the same promoter. The expression levels of the copper-inducible constructs were compared to endogenous Pho after 72 h of induction (Additional file 6: Figure S6b). After expression of the WT and mutant constructs for 72 h, respectively, the transcriptional kinetics of the heat shock response was assessed by measuring nascent RNA production from the hsp70 locus. The expression of the PhoV164D mutant did not affect the up-regulation of the hsp70 locus in response to a heat shock as compared to the control (Fig. 6c). However, upon withdrawing the heat shock stimulus, while nascent RNA production in the control cells sharply decreased, this was not observed in the case of PhoV164D mutant (Fig. 6c). Thus, the absence of a functional REPO domain negatively affected the recovery kinetics of the hsp70 locus during the heat shock response.

Discussion

Our study aimed at gaining mechanistic understanding about how PcG proteins regulate their target genes, which are transiently activated by an extracellular stimulus. To this end, we used the heat shock-responsive hsp70 gene in Drosophila as it exhibits a dynamic pattern of expression during the well-characterised heat shock response and it is also a target for PRC1. Furthermore, mechanistic analysis supports a role for PcG proteins in the re-establishment of silencing at the hsp70 locus in the absence of the heat shock stimulus.

Our results highlight antagonism between HSF and PRC1 for binding at the hsp70 locus, suggesting a role for PRC1 in repressing the locus in the absence of the heat shock stimulus. Furthermore, we could uncouple the eviction of PRC1 from the hsp70 locus upon heat shock from production of mRNA by treating cells with flavopiridol, a P-TEFb inhibitor. We demonstrated that Pho could preferentially interact with Spt5 and dSfmbt under ON/OFF conditions, respectively, at the hsp70 locus. Our data also conclusively showed that the recruitment of HSF to chromatin is sufficient and necessary for the dynamic interaction between Pho and Spt5. The REPO domain of Pho is central for the interaction switch, since a mutation in the REPO domain not only disrupted the kinetics of the protein–protein interaction between Pho and Spt5 but also disrupted the kinetics of nascent RNA production from the hsp70 locus during the heat shock response.

Several reports have described a role for Pho in both the activation and the repression of its target genes [19, 35, 36]. In order to dissect the relationship between binding patterns of PRC1, Pho and Trx and their respective effect on transcriptional activity, ChIP-chip was employed to identify their binding sites in the ANT-C and BX-C gene clusters in two different Drosophila cell lines. At silent loci, the Pho signal localised at PREs. However, in case of active loci, the Pho signal spread over the gene body [19]. Another study demonstrated the requirement of a Pho-binding site and a functional Pho protein for a pairing-sensitive silencing element in the Drosophila even skipped (eve) locus for propagation of the silenced as well as the activated state in different cellular contexts, respectively [35]. Yet another reported example is the Ubx gene in Ubx ON and Ubx OFF cells in the developing Drosophila larvae; the levels of Pho remained unchanged at the bx PRE in Ubx ON cells, while the occupancy of PRC1 and PRC2 was reduced by twofold [36].

Our results demonstrate the ability of Pho to switch its protein interaction partners in a HSF-dependent manner. Our current observations at the hsp70 locus lead us to propose that, in the absence of transcription factor binding, Pho interacts with dSfmbt, forming a recruitment platform for PRC1 and PRC2, and stably maintains the silent state [6, 9, 37]. However, upon activation, TF binding could lead to an increase in interaction affinity between Pho and Spt5, thereby relieving Polycomb-mediated silencing and enabling a switch to the activated state. It is thus plausible that dynamic protein–protein interaction switches between Polycomb group members and Spt5 could potentially play a role in enabling rapid switches in expression states at Pho target genes.

Our study also assessed the contribution of the dynamic switch in Pho protein–protein interactions on the transcriptional output of the hsp70 locus during the heat shock response. It was notable that Pho preferentially interacts with Spt5 and dSfmbt when the hsp70 locus is active and silent, respectively. A point mutation in the REPO domain of Pho resulted in ablation of its interaction with dSfmbt. While the PhoV164D mutant was capable of strongly interacting with Spt5, immediately upon heat shock, it was incapable of switching back to its interaction with dSfmbt in the absence of the heat shock stimulus. This result underscores the importance of a functional REPO domain in mediating the switch in the protein interaction partners of Pho depending upon the binding kinetics of HSF at the hsp70 locus. Nascent RNA analysis also suggested that the dynamic protein–protein interactions of Pho could contribute to the transcriptional regulation of the hsp70 locus. Expression of PhoV164D mutant negatively affected the transcriptional kinetics of the hsp70 locus in response to a heat shock when compared to the control. We speculate that in the absence of a functional REPO domain, the dynamic nature of the protein–protein interactions of Pho is lost, thereby preventing cells from efficiently switching between ON/OFF states in the presence/absence of the activating stimulus.

Taken together, our study highlights how dynamic protein interaction switches between PcG members and Spt5 could facilitate the transcriptional regulation of a Polycomb-regulated inducible gene like hsp70. Additional open questions remain regarding the complex relationship between PcG proteins and paused RNA polymerase II. Nonetheless, this study provides a new perspective on how the plasticity of Polycomb gene silencing is mediated by dynamic protein–protein interaction switches between activators and repressors.

Conclusions

In this study, we sought to understand the mechanistic basis of dynamic gene regulation by Polycomb group proteins at stress-inducible loci like hsp70. We identified a role for the dynamic protein–protein interaction switch between Pho and Spt5 in the transcriptional regulation of the hsp70 locus. Crucially, we show that the recruitment of HSF to chromatin and a functional REPO domain in Pho is necessary for the dynamic interaction switch between Pho and Spt5. Thus, our study reveals how a dynamic interaction switch between PcG proteins and an elongation factor could enable inducible genes regulated by Polycomb to switch between active and repressed transcriptional states.

Materials and methods

Cell culture

S2 DRSC cells (DGRC, stock #181) were cultured in Schneider’s medium (Sigma-Aldrich) + 10% heat-inactivated FBS (PAN Biotech GmbH) at 25 °C.

Cells were instantaneously heat shocked by adding Schneider’s medium, which was pre-warmed to 48 °C, and then, the temperature was maintained at 37 °C for 15 min by incubating the cells in a water bath as previously described [18]. The cells were then placed into a water bath at 25 °C to enable recovery from the heat shock. In case of treatment with 40 μM radicicol, 500 nm flavopiridol and 10 μM sodium salicylate, the treatment was performed for 1 h, 40 and 10 min, respectively, at 25 °C.

Chemicals and antibodies

The following chemicals were used in this study: flavopiridol (F3055, Sigma-Aldrich), radicicol (R-370, Alomone), sodium salicylate (S3007, Sigma-Aldrich), copper sulphate (102790, Merck).

The following primary antibodies were used in this study: α-mouse RNA polymerase II CTD repeat (Abcam, ab817), α-mouse RNA polymerase II CTD repeat phosphor S5 (Abcam, ab5408), α-rabbit RNA polymerase II CTD repeat phospho S2 (Abcam, ab5095), α-rat RNA polymerase II CTD repeat phosphor S2 (Active Motif, Cat#61083), α-rabbit HSF (gift from John Lis), α-rabbit Polyhomeotic (Paro lab stock), α-rabbit Pleiohomeotic (gift from Judith Kassis), α-rabbit Pleiohomeotic (gift from Jurg Muller), α-mouse FLAG (Sigma, F1804), α-rat HA (Roche, 11 867 423 001).

The following secondary antibodies were used in this study: α-mouse IgG HRP-linked whole antibody (GE Healthcare, NA931), α-rabbit IgG HRP-linked whole antibody (GE Healthcare, NA934V), α-rat IgG HRP-linked whole antibody (GE Healthcare, NA935V), Goat α-mouse IgG Alexa Fluor 488 (Thermo Fisher Scientific, A-11001), Goat α-rabbit IgG Alexa Fluor 568 (Thermo Fisher Scientific, A-11011).

Gene cloning

cDNA clones for Spt5 (LD10265) and dSfmbt (LD14884) were obtained from the Drosophila Genome Resource Centre (DGRC). cDNA for Pho and Spt4 was prepared from S2 DRSC cells using the First Strand cDNA synthesis kit (Thermo Fisher Scientific, K1612).

In order to express the constructs in S2 DRSC cells, the cDNA clones were amplified using appropriate primers and cloned in pENTR/D-TOPO vector using the pENTR™/D-TOPO™ cloning kit (Thermo Fisher Scientific, K240020). The entry clones were subsequently cloned into pAFW, pAHW, pAWH, pMT/FHW (Drosophila Gateway Vector Collection by DGRC) by using the Gateway™ LR clonase™ II enzyme mix (Thermo Fisher Scientific, 11791020). All plasmids were sequenced by Sanger Sequencing (Microsynth AG). All primer sequences used in this study are listed in Additional file 7: Table S1.

Cell Transfection

In case of transient co-transfection for co-IP experiments, 1 μg of each plasmid was used for transfection, and procedure was performed according to the Effectene transfection reagent protocol (Qiagen, 301425).

In case of stable transfection for generation of cell lines, 2 μg of pMT/FHW HygRes [38] was transfected into cells according to the Effectene transfection reagent protocol and cells were selected with Hygromycin B (Thermo Fisher Scientific, 10687010) 3 days post-transfection until the emergence of drug-resistant colonies was observed.

RNAi

Sequences used for the production of double-stranded RNA against HSF and LacZ have been previously described [23]. The procedure of dsRNA production was performed according to the manufacturer’s protocol (MEGAscript™ RNAi kit, Thermo Fisher Scientific, AM1626).

Cells were treated with dsRNA according to the protocol by DRSC (Cell RNAi, 6-well) with the following exceptions. One million cells were plated in 1 ml of serum-free medium onto 6-well plates. Ten micrograms of dsRNA was added to the well, and this was followed by a 45-min incubation period. Thereafter, 1 ml of Schneider’s medium + 20% heat-inactivated FBS was added to the wells, and the cells were incubated for a period of 96 h.

Chromatin immunoprecipitation

Approximately 25 × 106 cells were resuspended in PBS and fixed with 1% formaldehyde (Sigma, F8775) for 10 min. The reaction was quenched with 125 mM glycine for 5 min. Cells were washed with ice-cold PBS and centrifuged at 500g for 5 min. The cells were resuspended in buffer A (100 mM Tris ph8, 10 mM DTT) and incubated at 4 °C on a rotor for 15 min. Thereafter, the cells were incubated on a thermomixer at 30 °C at 300 rpm. The cells were briefly vortexed for 5 s and centrifuged at 500g for 5 min. The cells were resuspended in buffer B (10 mM HEPES ph8, 10 mM EDTA, 10 mM EGTA, 0.25% Triton X-100) and incubated at 4 °C on a rotor for 5 min. The cells were centrifuged at 500g for 5 min. The cells were resuspended in buffer C (10 mM HEPES ph8, 10 mM EDTA, 0.5 mM EGTA, 200 mM NaCl) and incubated at 4 °C on a rotor for 5 min. The cells were centrifuged at 500g for 5 min. The cells were resuspended in buffer D (50 mM HEPES ph8, 10 mM EDTA, 0.5 mM EGTA, 0.1% SDS) and sonicated using S220 ultrasonicator (Covaris) to obtain chromatin fragments approximately 500 bp long.

The chromatin was adjusted to RIPA conditions and centrifuged at 4 °C for 10 min at 13,000 rpm. One hundred microlitres of the chromatin was aliquoted as an input fraction. Appropriate amount of the antibody (amount of each antibody used per IP is stated in Additional file 8: Table S2) was added to the chromatin solution and incubated overnight on a rotor at 4 °C.

In order to prepare the Dynabeads™ Protein A/G (Thermo Fisher Scientific) for purification of the immunocomplexes, 40 μl of the magnetic beads was washed with 1 ml RIPA. The ChIP sample was centrifuged at 4 °C for 5 min at 13,000 rpm, transferred to the tube containing the magnetic beads and incubated at 4 °C on a rotor for 2–3 h. The immunocomplexes along with the magnetic beads were washed five times with RIPA buffer (10 mM Tris ph8, 1 mM EDTA, 1% Triton X-100, 0.1% SDS, 0.1% Na-deoxycholate), once with LiCL buffer (10 mM Tris ph8, 250 mM LiCl, 1 mM EDTA, 0.5% NP40, 0.5% Na-deoxycholate) and twice with TE (10 mM Tris ph8, 1 mM EDTA). The immunocomplexes were eluted thrice with the elution buffer (1% SDS, 0.1 M NaHCO3) at 65 °C for 15 min with gentle agitation (800 rpm). The ChIP sample was reverse cross-linked with 200 mM NaCl and incubated at 65 °C overnight. The input and ChIP samples were treated with RNase (Roche, 11 119 915 001) and incubated at 37 °C for 90 min. This was followed by treatment with Proteinase K (Axon Lab AG, A4392.0001) and incubated at 55 °C for 90 min. The DNA was purified using QIAquick PCR purification kit (QIAGEN, 28106). All buffers were supplemented with protease inhibitors (cOmplete EDTA-free, Roche, 4693132001).

Polytene immunostaining

Drosophila salivary glands were dissected from third-instar larvae, permeabilised in solution 1 (1xPBS, 1% Triton X-100) and then fixed for 10 min (50% acetic acid, 3.7% p-formaldehyde). The squashed salivary glands were prepared as described previously [39]. The slides were washed with PBS and blocked for 60 min with the blocking solution (3% BSA, 0.2% NP40, 0.2% Tween 20, 10% non-fat dry milk powder). Dilutions of the primary antibody (Pho, 1:50 and RNA pol II, 1:100) were made with the blocking solution, and the slides were incubated with the primary antibody overnight in a humid chamber at 4 °C. The slides were washed with PBS, dilutions of the secondary antibody (1:200) were made with the blocking solution and the slides were incubated with the secondary antibodies for 60 min in a humid chamber at room temperature. The slides were washed with wash buffer ‘300’ and wash buffer ‘400’ (1xPBS, 300/400 mM NaCl, 0.2% NP40, 0.2% Tween 20). The slides were then stained with 1 μg/ml DAPI (Roche, 236 276), mounted using 40 μl Fluoromount G (Southern Biotech, 0100-01) and analysed using a Leica DMI6000 B fluorescent microscope.

Native ChIP (N-ChIP)

The protocol for N-ChIP had been previously described in [27]. Approximately 2 × 108 cells were harvested by scraping, washed with ice-cold PBS and centrifuged at 1200g at 4 °C for 4 min. The cells were resuspended in TM2+ buffer (10 mM Tris pH 7.5, 2 mM MgCl2, 0.5 mM PMSF). Nuclei were released by adding 0.6% NP-40 to the cell suspension and incubating for 3 min and 30 s on ice. During the incubation, the cell suspension was gently vortexed at intervals of 1 min. The nuclei were centrifuged at 4 °C for 10 min at 1000g and washed once with TM2+. The nuclei were resuspended in TM2+ with protease inhibitors (cOmplete EDTA-free, Roche, 4693132001). MNase treatment was performed as follows: the nuclei were preheated at 37 °C for 3 min. 1 mM CaCl2 was added to the nuclei followed by addition of 80U MNase (M0247S, NEB). The digestion with MNase was carried out for precisely 6 min, and the reaction was stopped by addition of 2 mM EGTA. The nuclei were incubated on ice for 2 min after the MNase treatment and centrifuged at 4 °C for 10 min at 1000g. Nuclei were washed with TM2+ containing protease inhibitors and 2 mM EGTA and centrifuged at 4 °C for 10 min at 1000g. The nuclei were resuspended in 80T buffer (70 mM NaCl, 10 mM Tris pH 7.5, 2 mM MgCl2, 2 mM EGTA, 0.1% Triton X-100, 0.5 mM PMSF, protease inhibitors), cavitated 10 times through a 26G × ½” needle (4710004512, Henke Sass Wolf), and chromatin was centrifuged at 4 °C for 10 min at 1000g. Ten percentage of the chromatin was aliquoted as the input fraction. Anti-Pho antibody (rabbit polyclonal, a gift from Jurg Muller) was used for each IP at 1:100 dilution and incubated at 4 °C on a rotating wheel overnight. Forty microlitres of Dynabeads™ Protein G magnetic beads was added to the IP and incubated at 4 °C for 2 h. Beads were washed twice with 80T buffer, and beads were resuspended in 80T buffer. Following treatment with RNase A (Thermo Fisher Scientific, EN0531) and Proteinase K (Axon Lab AG, A4392.0001), DNA was extracted by phenol–chloroform and ethanol precipitation.

Co-immunoprecipitation (co-IP)

The co-IP protocol was adapted from [40]. Cells were harvested, washed with ice-cold PBS and centrifuged at 4 °C for 5 min at 500g. The cells were resuspended and lysed in co-IP buffer (50 mM Tris ph7.5, 150 mM NaCl, 1% NP40, 1 mM EGTA, 20 mM NaF) for 30 min at 4 °C. The cells were centrifuged at 4 °C for 20 min at 12,000 rpm. Prior to treating the cell lysate with nucleases, 0.4 mM MgCl2 was supplemented to the cell lysate. The cell lysate was incubated with RNase A (Thermo Fisher Scientific, EN0531) or DNase I (Thermo Fisher Scientific, AM2222) for 10 min at room temperature. The reaction with the nuclease was terminated by addition of 5 mM EDTA to the cell lysate. The supernatant was transferred to a siliconised tube. Ten percentage was aliquoted as the input fraction. 2.5 μg of α-mouse FLAG was added to the IP fraction and incubated at 4 °C for 2–3 h. Forty microlitres of Dynabeads™ Protein G magnetic beads was added to the IP fraction and incubated at 4 °C for 2–3 h. The immunocomplexes were washed thrice with co-IP wash buffer (10 mM Tris ph7.5, 1 mM EDTA, 1 mM EGTA, 150 mM NaCl, 1% Triton X-100). The immunocomplexes were eluted with co-IP buffer, denatured with a mixture containing NuPage® 4x LDS sample buffer (Thermo Fisher Scientific, NP0008) and 200 mM DTT upon heating to 70 °C with mild agitation at 800 rpm for 15 min. Western blots were visualised using Amersham ECL Western Blotting Detection reagent (GE Healthcare, RPN2106) and Amersham Hyperfilm ECL (GE Healthcare, 28906836). All buffers were supplemented with protease inhibitors (cOmplete EDTA-free, Roche, 4693132001) and phosphatase inhibitors (PhosSTOP phosphatase inhibitor cocktail tablets, Roche, 4906837001). For the quantification of HA-Pho co-IP, the intensity of the co-IP band was normalised against the intensity of the corresponding FLAG IP band using ImageJ software.

Native elongating transcript isolation and quantification

The protocol used for the isolation of native elongating transcripts had been previously described [41]. Cells were grown in medium containing 0.5 mM CuSO4 for 72 h to induce the expression of stably integrated constructs. Cells were harvested and centrifuged at 2000 rpm for 5 min at 4 °C. The cells were resuspended in PBS and centrifuged at 2000 rpm for 5 min at 4 °C. The cells were resuspended in cytoplasmic buffer (10 mM Tris pH 7, 150 mM NaCl, 0.15% NP-40) and incubated on ice for 5 min. The cell lysate was layered on top of NET-seq sucrose buffer (10 mM Tris pH 7, 150 mM NaCl, 0.34 M sucrose) and centrifuged at 16,000g for 10 min at 4 °C. Nuclei were washed twice with nuclei wash buffer (1 mM EDTA, 0.05% Triton X-100, PBS) and centrifuged at 2000g for 5 min at 4 °C. The nuclei were resuspended in glycerol buffer (20 mM Tris pH 8, 75 mM NaCl, 0.5 mM EDTA, 50% glycerol, 0.85 mM DTT) and harsh nuclei lysis buffer (20 mM HEPES ph 8, 300 mM NaCl, 1% NP-40, 1 M urea, 0.2 mM EDTA, 1 mM DTT). This was followed by 4 rounds of vortexing for 2 s each with an incubation period of 1 min on ice between each round of vortexing. Samples were incubated on ice for 2 min and centrifuged at 18,500g for 2 min at 4 °C. The supernatant represented the nucleoplasmic fraction while the pellet contained the chromatin fraction. The chromatin fraction was washed once using a mix of the glycerol and harsh nuclei lysis buffer. The pellet was resuspended in 500 μl of TRIzol. All the buffers were supplemented with protease inhibitors (cOmplete EDTA-free, Roche, 4693132001) and actinomycin D (Sigma, A1410).

Nascent RNA was isolated from the chromatin fraction following manufacturer’s protocol (Thermo Fisher Scientific, 15596018). The RNA was precipitated overnight using ethanol. Mature RNAs were depleted from the sample by performing two rounds of Oligo dT-based depletion (Dynabeads™ Oligo(dT)25, Thermo Fischer Scientific, 61002). The RNA was DNase-treated using Turbo DNA-free™ kit (Thermo Fisher Scientific, AM1907) and reverse transcribed using First Strand cDNA synthesis kit (Thermo Fisher Scientific, K1612). cDNA was diluted (1:50), mixed with FastStart Essential DNA Green Master (Roche, 06402712001) and assayed using Light Cycler® 96 (Roche). ΔΔCT method was used to measure transcript levels relative to control.

Quantitative real-time PCR

For measurement of mRNA levels using qRT-PCR, cells were lysed using TRIzol™ reagent (Thermo Fisher Scientific, 15596018), and RNA was extracted using Direct-zol™ RNA MiniPrep (Zymo Research, R2052). The RNA was DNase-treated using Turbo DNA-free™ kit (Thermo Fisher Scientific, AM1907) and reverse transcribed using First Strand cDNA synthesis kit (Thermo Fisher Scientific, K1612). cDNA was diluted (1:50), mixed with FastStart Essential DNA Green Master (Roche, 06402712001) and assayed using Light Cycler® 96 (Roche). ΔΔCT method was used to measure transcript levels relative to control.

For measurement of % input values for the ChIP and N-ChIP samples, the input DNA samples were diluted in fivefold dilution series (1:10 up to 1:1250) and the ChIP sample was diluted (1:5). The standard curve obtained from the input DNA dilutions was used to determine the amount of DNA in the ChIP, which was calculated in terms of % input.

Abbreviations

PRC1:

Polycomb repressive complex 1

hsp70 :

Heat shock protein 70

HSF:

Heat shock factor

Pho:

Pleiohomeotic

PRE:

Polycomb response element

DSIF:

DRB sensitivity-inducing factor

NELF:

Negative elongation factor

P-TEFb:

Positive transcription elongation factor

dSfmbt:

Drosophila Scm-related gene containing four MBT domains

Spt5:

Suppressor of TY’s

CTD:

C-terminal domain

S2P:

Serine 2 phosphorylated

ChIP:

Chromatin immunoprecipitation

N-ChIP:

Native chromatin immunoprecipitation

co-IP:

Co-immunoprecipitation

YY1:

Ying yang 1

REPO:

Recruitment of Polycomb domain

References

  1. Prezioso C, Orlando V. Polycomb proteins in mammalian cell differentiation and plasticity. FEBS Lett. 2011;585:2067–77.

    Article  CAS  PubMed  Google Scholar 

  2. Schwartz YB, Pirrotta V. A new world of Polycombs: unexpected partnerships and emerging functions. Nat Rev Genet. 2013;14:853–64.

    Article  CAS  PubMed  Google Scholar 

  3. Brown JL, Mucci D, Whiteley M, Dirksen ML, Kassis JA. The Drosophila Polycomb group gene pleiohomeotic encodes a DNA binding protein with homology to the transcription factor YY1. Mol Cell. 1998;1:1057–64.

    Article  CAS  PubMed  Google Scholar 

  4. Müller J, Kassis JA. Polycomb response elements and targeting of Polycomb group proteins in Drosophila. Curr Opin Genet Dev. 2006;16:476–84.

    Article  PubMed  Google Scholar 

  5. Alfieri C, Gambetta MC, Matos R, Glatt S, Sehr P, Fraterman S, et al. Structural basis for targeting the chromatin repressor Sfmbt to Polycomb response elements. Genes Dev. 2013;27:2367–79.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Frey F, Sheahan T, Finkl K, Stoehr G, Mann M, Benda C, et al. Molecular basis of PRC1 targeting to Polycomb response elements by PhoRC. Genes Dev. 2016;30:1116–27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Francis NJ, Kingston RE, Woodcock CL. Chromatin compaction by a polycomb group protein complex. Science. 2004;306:1574–7.

    Article  CAS  PubMed  Google Scholar 

  8. Pengelly AR, Copur Ö, Jäckle H, Herzig A, Müller J. A histone mutant reproduces the phenotype caused by loss of histone-modifying factor Polycomb. Science. 2013;339:698–9.

    Article  CAS  PubMed  Google Scholar 

  9. Laprell F, Finkl K, Müller J. Propagation of Polycomb-repressed chromatin requires sequence-specific recruitment to DNA. Science. 2017;356:85–8.

    Article  CAS  PubMed  Google Scholar 

  10. Coleman RT, Struhl G. Causal role for inheritance of H3K27me3 in maintaining the OFF state of a Drosophila HOX gene. Science. 2017;356:eaai8236.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Enderle D, Beisel C, Stadler MB, Gerstung M, Athri P, Paro R. Polycomb preferentially targets stalled promoters of coding and noncoding transcripts. Genome Res. 2011;21:216–26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Rougvie AE, Lis JT. The RNA polymerase II molecule at the 5′ end of the uninduced hsp70 gene of D. melanogaster is transcriptionally engaged. Cell. 1988;54:795–804.

    Article  CAS  PubMed  Google Scholar 

  13. Missra A, Gilmour DS. Interactions between DSIF (DRB sensitivity inducing factor), NELF (negative elongation factor), and the Drosophila RNA polymerase II transcription elongation. Proc Natl Acad Sci USA. 2010;107:11301–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Wu CH, Yamaguchi Y, Benjamin LR. NELF and DSIF cause promoter proximal pausing on the hsp70 promoter in Drosophila. Genes Dev. 2003;17:1402–14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Westwood JT, Clos J, Wu C. Stress-induced oligomerization and chromosomal relocation of heat-shock factor. Nature. 1991;353:822.

    Article  CAS  PubMed  Google Scholar 

  16. Lis JT, Mason P, Peng J, Price DH, Werner J. P-TEFb kinase recruitment and function at heat shock loci. Genes Dev. 2000;14:792–803.

    CAS  PubMed  PubMed Central  Google Scholar 

  17. Ivanov D, Kwak YT, Guo J, Gaynor RB. Domains in the SPT5 protein that modulate its transcriptional regulatory properties. Mol Cell Biol. 2000;20:2970–83.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Boehm AK, Saunders A, Werner J, Lis JT. Transcription factor and polymerase recruitment, modification, and movement on dhsp70 in vivo in the minutes following heat shock. Mol Cell Biol. 2003;23:7628–37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Beisel C, Buness A, Roustan-Espinosa IM, Koch B, Schmitt S, Haas SA, et al. Comparing active and repressed expression states of genes controlled by the Polycomb/Trithorax group proteins. Proc Natl Acad Sci USA. 2007;104:16615–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Harvey R, Schuster E, Jennings BH. Pleiohomeotic interacts with the core transcription elongation factor Spt5 to regulate gene expression in Drosophila. PLoS ONE. 2013;8:e70184.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Andrulis ED, Guzmán E, Döring P, Werner J, Lis JT. High-resolution localization of Drosophila Spt5 and Spt6 at heat shock genes in vivo: roles in promoter proximal pausing and transcription elongation. Genes Dev. 2000;14:2635–49.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Teves SS, Henikoff S. Heat shock reduces stalled RNA polymerase II and nucleosome turnover genome-wide. Genes Dev. 2011;25:2387–97.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Duarte FM, Fuda NJ, Mahat DB, Core LJ, Guertin MJ, Lis JT. Transcription factors GAF and HSF act at distinct regulatory steps to modulate stress-induced gene activation. Genes Dev. 2016;30:1731–46.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Boyer LA, Plath K, Zeitlinger J, Brambrink T, Medeiros LA, Lee TI, et al. Polycomb complexes repress developmental regulators in murine embryonic stem cells. Nature. 2006;441:349–53.

    Article  CAS  PubMed  Google Scholar 

  25. Gaertner B, Johnston J, Chen K, Wallaschek N, Paulson A, Garruss AS, et al. Poised RNA polymerase II changes over developmental time and prepares genes for future expression. Cell Rep. 2012;2:1670–83.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Chao SH. Flavopiridol inactivates P-TEFb and blocks most RNA polymerase II transcription in vivo. J Biol Chem. 2001;276:31793–9.

    Article  CAS  PubMed  Google Scholar 

  27. Kasinathan S, Orsi GA, Zentner GE, Ahmad K, Henikoff S. High-resolution mapping of transcription factor binding sites on native chromatin. Nat Methods. 2013;11:203–9.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Ni Z, Schwartz BE, Werner J, Suarez J-R, Lis JT. Coordination of transcription, RNA processing, and surveillance by P-TEFb kinase on heat shock genes. Mol Cell. 2004;13:55–65.

    Article  CAS  PubMed  Google Scholar 

  29. Mohd-Sarip A, Venturini F, Chalkley GE, Verrijzer CP. Pleiohomeotic can link polycomb to DNA and mediate transcriptional repression. Mol Cell Biol. 2002;22:7473–83.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Klymenko T, Papp B, Fischle W, Köcher T, Schelder M, Fritsch C, et al. A Polycomb group protein complex with sequence-specific DNA-binding and selective methyl-lysine-binding activities. Genes Dev. 2006;20:1110–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Wilkinson FH, Park K, Atchison ML. Polycomb recruitment to DNA in vivo by the YY1 REPO domain. Proc Natl Acad Sci USA. 2006;103:19296–301.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Zou J, Guo Y, Guettouche T, Smith DF, Voellmy R. Repression of heat shock transcription factor HSF1 activation by HSP90 (HSP90 complex) that forms a stress-sensitive complex with HSF1. Cell. 1998;94:471–80.

    Article  CAS  PubMed  Google Scholar 

  33. Winegarden NA, Wong KS, Sopta M, Westwood JT. Sodium salicylate decreases intracellular ATP, induces both heat shock factor binding and chromosomal puffing, but does not induce hsp 70 gene transcription in Drosophila. J Biol Chem. 1996;271:26971–80.

    Article  CAS  PubMed  Google Scholar 

  34. Qiu Y, Gilmour DS. Identification of regions in the Spt5 subunit of DRB sensitivity-inducing factor (DSIF) that are involved in promoter-proximal pausing. J Biol Chem. 2017;292:5555–70.

    Article  CAS  PubMed  Google Scholar 

  35. Fujioka M, Yusibova GL, Zhou J, Jaynes JB. The DNA-binding Polycomb-group protein Pleiohomeotic maintains both active and repressed transcriptional states through a single site. Development. 2008;135:4131–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Papp B. Histone trimethylation and the maintenance of transcriptional ON and OFF states by trxG and PcG proteins. Genes Dev. 2006;20:2041–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Wang L, Brown JL, Cao R, Zhang Y, Kassis JA, Jones RS. Hierarchical recruitment of polycomb group silencing complexes. Mol Cell. 2004;14:637–46.

    Article  CAS  PubMed  Google Scholar 

  38. Kockmann T, Gerstung M, Schlumpf T, Xhinzhou Z, Hess D, Beerenwinkel N, et al. The BET protein FSH functionally interacts with ASH1 to orchestrate global gene activity in Drosophila. Genome Biol. 2013;14:R18.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Paro R. Mapping protein distributions on polytene chromosomes by immunostaining. Cold Spring Harbor Protoc. 2008;2008:pdb.prot4714-4.

    Article  Google Scholar 

  40. Genevet A, Wehr MC, Brain R, Thompson BJ, Tapon N. Kibra is a regulator of the Salvador/Warts/Hippo signaling network. Dev Cell. 2010;18:300–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Mayer A, Churchman LS. Genome-wide profiling of RNA polymerasetranscription at nucleotide resolution in humancells with native elongating transcript sequencing. Nat Protoc. 2016;11:813–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Authors’ contributions

The study was conceived by AP and RP. All the experiments and analyses were performed by AP. The results of the experiments were interpreted by AP and RP. The manuscript was written by AP and RP. Both authors read and approved the final manuscript.

Acknowledgements

We would like to thank Jorge Beira, Federico Comoglio and Cem Sievers for their critical comments on the manuscript. We would also like to acknowledge Fabian Rudolf for the discussion and input during the generation of the Pho mutants and Christian Beisel for the discussion and input on performing N-ChIP. We would like to thank John Lis, Judith Kassis and Jurg Muller for sharing antibodies. Allwyn Pereira is a member of the Molecular Life Sciences Programme, Life Science Graduate School Zurich.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

All data generated or analysed during this study are included in this published article [and its additional files].

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.

Funding

This work was supported by a grant from the Swiss National Science Foundation to Renato Paro and the ETH Zurich.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Renato Paro.

Additional files

13072_2017_166_MOESM1_ESM.jpg

Additional file 1: Figure S1. Schematic representation of the proteins involved in the silencing, activation and re-silencing of the hsp70 locus. a At optimal growth temperature, Pho, a DNA-binding PcG member, binds to promoter region of the hsp70 locus. Pho interacts with dSfmbt, which together form a recruitment platform for PRC1 to the hsp70 locus. In addition, RNA polymerase II is maintained in the paused state by NELF and Spt5, which act as pausing factors. Upon heat shock, HSF, along with P-TEFb, is recruited to chromatin and releases RNA polymerase II from the paused state. P-TEFb modifies Spt5 and converts into an elongation factor. It also modifies the CTD of RNA polymerase II to enable productive elongation. However, upon removal of the heat shock stimulus, the locus should eventually return to its paused state. Thus, hsp70 is an ideal model gene to study the eviction and recruitment of PRC1 upon activation and re-silencing, respectively. The colour code for the protein names is as follows: silencing in red, pausing factors in orange and activators in green.

13072_2017_166_MOESM2_ESM.jpg

Additional file 2: Figure S2. RNA polymerase II binding dynamics at Act42A locus during the heat shock response. a, b ChIP-qPCR measurements of occupancy levels of RNA polymerase II CTD and S2P form of RNA polymerase II at the Act42A locus during the heat shock response detailed in Fig. 1a. The control line represents cells that were maintained at 25 °C for the entire duration of the time course detailed in Fig. 1a. Distance of the location of the primers used from the Act42A TSS is as follows: for a and b (+ 104 bp) and has been depicted in the form of a cartoon below the data figure. Data information: In (a, b), data are presented as mean ± SEM for n = 2.

13072_2017_166_MOESM3_ESM.jpg

Additional file 3: Figure S3. Effect of flavopiridol on RNA polymerase II occupancy at the hsp70 locus upon heat shock. a N-ChIP-qPCR measurement of the occupancy level of Pho at the hsp70 locus in S2 DRSC cells during the heat shock response detailed in Fig. 1a. Distance of the location of the primers used from the hsp70 TSS is − 85 bp and has been depicted in the form of a cartoon below the data figure. b, c ChIP-qPCR measurements of occupancy levels of RNA polymerase II CTD and S2P form of RNA polymerase II, respectively, at the hsp70 locus in S2 DRSC cells. Distances of the location of the primers used from the hsp70 TSS are as follows: for a (+ 78 bp) and for b (+ 802 bp) and has been depicted in the form of a cartoon below the data figure. d qRT-PCR measurements of hsp70 transcript levels in third-instar Drosophila larvae under the conditions used for double polytene immunostaining. Distance of the location of the primers used from the hsp70 TSS is + 687 bp and has been depicted in the form of a cartoon below the data figure. For b–d, the control bar represents cells that were maintained at 25 °C for 15 min, the HS bar represents cells that were maintained at 37 °C for 15 min, the FP control bar represents cells that were treated with 500 nM flavopiridol for 40 min and then maintained at 25 °C for 15 min and the FP HS bar represents cells that were treated with 500 nM flavopiridol for 40 min and then maintained at 37 °C for 15 min. Data information: In (a–d), data are presented as mean ± SEM (n = 2).

13072_2017_166_MOESM4_ESM.jpg

Additional file 4: Figure S4. Chromatin binding dynamics of Spt5 and the dissection of the protein–protein interaction domains of Pho. a ChIP-qPCR measurements of occupancy levels of FLAG-Spt5 at the hsp70 locus over the time course detailed in Fig. 1a. The control line represents cells that were maintained at 25 °C for the entire duration of the time course detailed in Fig. 1a. The cartoon at the bottom of the figure represents the distance of the location of the primers used from the hsp70 TSS. b–c co-IP assays of S2 DRSC cells transiently transfected with plasmids expressing FLAG-tagged dSfmbt or Spt5 (FLAG-dSfmbt, FLAG-Spt5) and HA-tagged Pho, N-Pho (a.a 1-351), 1-200N-Pho (a.a 1-200), 173-351N-Pho (a.a 173-351), C-Pho (a.a 352-520). Cell lysates were used for pull-downs using an anti-FLAG antibody and were later probed by Western blot using an anti-HA antibody. MW = molecular weight in kDa. d co-IP assays for S2 DRSC cells transiently transfected with plasmids expressing FLAG-tagged Pho and HA-tagged Spt5 or HA-tagged dSfmbt. S2 DRSC cells were either maintained at normal growth temperature (25 °C) or heat shocked at 37 °C for 15 min. Cell lysates were used for pull-downs using an anti-FLAG antibody and were later probed by Western blot using an anti-HA antibody. MW = molecular weight in kDa.

13072_2017_166_MOESM5_ESM.jpg

Additional file 5: Figure S5. The dynamic interaction switch between Pho and Spt5 is independent of DNA but dependent upon RNA. a co-IP assays of S2 DRSC cells transiently transfected with plasmids expressing FLAG-tagged Spt5 (FLAG-Spt5) and HA-tagged Pho (HA-Pho). S2 DRSC cells were either maintained at normal growth temperature (25 °C) or heat shocked at 37 °C for 15 min. Prior to performing the pull-down, the cell lysates were treated with either RNase A or DNase I. Thereafter, the treated cell lysates were used for pull-downs using an anti-FLAG antibody and later probed by Western blot using an anti-HA antibody. b co-IP assays of S2 DRSC cells transiently transfected with plasmids expressing FLAG-tagged Spt5 (FLAG-Spt5) and HA-tagged N-Pho (HA-N-Pho). S2 DRSC cells were either maintained at normal growth temperature (25 °C) or heat shocked at 37 °C for 15 min. Cell lysates were prepared and were used for pull-downs using an anti-FLAG antibody. They were later probed by Western blot using an anti-HA antibody.

13072_2017_166_MOESM6_ESM.jpg

Additional file 6: Figure S6. The dynamic interaction switch between Pho V164D–Spt5 is independent of DNA but dependent upon RNA. a co-IP assays of S2 DRSC cells transiently transfected with plasmids expressing FLAG-tagged Spt5 (FLAG-Spt5) and HA-tagged Pho V164D (HA-Pho V164D). S2 DRSC cells were either maintained at normal growth temperature (25 °C) or heat shocked at 37 °C for 15 min. Prior to performing the pull-down, the cell lysates were treated with either RNase A or DNase I. Thereafter, the treated cell lysates were used for pull-downs using an anti-FLAG antibody and later probed by Western blot using an anti-HA antibody. b A western blot depicting the levels of the copper-inducible constructs (Pho WT/Pho V164D) after 72 h of induction to the endogenous Pho. MW = molecular weight in kDa.

Additional file 7: Table S1. List of the primer sequences used in the study.

Additional file 8: Table S2. Amount of antibody used per IP in ChIP.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Pereira, A., Paro, R. Pho dynamically interacts with Spt5 to facilitate transcriptional switches at the hsp70 locus. Epigenetics & Chromatin 10, 57 (2017). https://doi.org/10.1186/s13072-017-0166-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13072-017-0166-9

Keywords