Transcription-coupled recruitment of human CHD1 and CHD2 influences chromatin accessibility and histone H3 and H3.3 occupancy at active chromatin regions
Epigenetics & Chromatin volume 8, Article number: 4 (2015)
CHD1 and CHD2 chromatin remodeling enzymes play important roles in development, cancer and differentiation. At a molecular level, the mechanisms are not fully understood but include transcriptional regulation, nucleosome organization and turnover.
Here we show human CHD1 and CHD2 enzymes co-occupy active chromatin regions associated with transcription start sites (TSS), enhancer like regions and active tRNA genes. We demonstrate that their recruitment is transcription-coupled. CHD1 and CHD2 show distinct binding profiles across active TSS regions. Depletion of CHD1 influences chromatin accessibility at TSS and enhancer-like chromatin regions. CHD2 depletion causes increased histone H3 and reduced histone variant H3.3 occupancy.
We conclude that transcription-coupled recruitment of CHD1 and CHD2 occurs at transcribed gene TSSs and at intragenic and intergenic enhancer-like sites. The recruitment of CHD1 and CHD2 regulates the architecture of active chromatin regions through chromatin accessibility and nucleosome disassembly.
ATP-dependent chromatin remodeling enzymes disrupt DNA-histone interactions facilitating nucleosome mobilization, including the disassembly, eviction, sliding and spacing of nucleosomes . The ability of ATP-dependent chromatin remodeling enzymes to manipulate nucleosomes also facilitates diverse chromatin-associated processes including transcription, DNA repair, incorporation of histone variants, DNA methylation, the covalent modification of histone tails and gene activation or silencing . All ATP-dependent chromatin remodeling enzymes contain a SNF2-like helicase domain . The chromodomain helicase DNA binding protein (CHD) subfamily of ATP-dependent chromatin remodeling enzymes are characterized by the presence of dual chromatin organization modifier (chromo) domains N-terminal to the SNF2 domain . Human CHD1 and CHD2 are grouped as subclass I CHD remodelers based on both their modular organization, with N-terminal dual chromodomains and on their sequence similarity within the SNF2 helicase domain in relation to CHD1 homologues in other organisms [3, 5]. In contrast, other human CHD family remodelers are grouped into class II and class III CHD proteins on the basis of additional regulatory domains such as plant homeobox domains as well as protein sequence differences in the SNF2 helicase domain [3, 5].
Mammalian CHD1 has important biological roles in stem cell function, prostate cancer and transcription, while CHD2 has emerged as an important regulator of genetic stability, development, kidney function, muscle cell and hematopoietic differentiation [6–15]. In model organisms, studies on CHD1 homologues have demonstrated various mechanistic roles in nucleosome disassembly, histone turnover, nucleosome spacing, transcription through chromatin, transcription-coupled chromatin assembly and deposition of the histone H3 variant H3.3 during development [10, 16–19]. Relatively little is known about CHD2 function on a molecular level with the exception of muscle cell differentiation where it was demonstrated that CHD2 is an essential cofactor for the transcription factor MyoD in H3.3 deposition at muscle cell gene promoters .
Interestingly, CHD family remodelers function as chromatin readers in some cases by recognizing specific histone modifications . For example human CHD1 has specific affinity for histone 3 (H3) tails when di- or tri-methylated at lysine 4 (H3K4me2 and H3K4me3) . The role of H3K4me2/3 in the physical recruitment of CHD1 per se is not clear although both are present at transcribed gene promoters . For example, some studies have demonstrated that CHD1 recruitment is at least partially independent of the chromodomains [21, 22], while other studies have linked CHD1 to chromatin associated with transcriptional elongation [23–25]. In the model organisms, budding yeast and Drosophila, CHD1 homologues do not show any specific affinity for H3K4me2/3 and bind unmodified H3 [20, 21].
The fission yeast CHD1 homologue Hrp1 was shown to co-purify with the fission yeast mediator complex . Consistent with this finding, it was later demonstrated that mediator-dependent assembly of the pre-initiation complex in vitro is associated with CHD1 recruitment in mammalian cells . Surprisingly, the latter study demonstrated CHD1 recruitment to pre-initiation complexes on both naked DNA and chromatin-based templates, including chromatin templates with and without H3K4me3 . Unlike CHD1, the chromodomains of CHD2 do not have specific affinity H3K4me3 . In mouse and human cells, the myogenic transcription factor MyoD is sufficient for the recruitment of CHD2 to muscle lineage gene promoters, which may be mediated through a direct interaction between MyoD and CHD2 . It was also suggested that CHD2 does not bind to housekeeping genes, but operates as a muscle specific factor in muscle lineage activation. Thus, exactly how chromatin remodeling enzymes such as CHD1 and CHD2 are recruited to target sites remains an open question.
To examine how subclass I CHD remodelers, CHD1 and CHD2 are recruited and function in human cells, we first carefully examined the occupancy of each in relation to chromatin states utilizing data from the Encyclopedia of DNA Elements (ENCODE) consortium (http://genome.ucsc.edu/ENCODE/). Closer inspection of the data demonstrates that neither H3K4me3 nor micrococcal nuclease (MNase) sensitive DNA at the NDR is correlated to CHD1 or CHD2. We demonstrate that both CHD1 and CHD2 are recruited in association with the RNA polymerase II (Pol II) machinery to active chromatin regions. Depletion of CHD1 in human cells reduces the accessibility of active regions as measured by DNase sensitivity. We also show that at active chromatin regions H3 occupancy is increased following knockdown of CHD2 suggesting a role in promoting nucleosome disassembly at active regions genome-wide. Knockdown of CHD2 also leads to a reduction in the relative enrichment of H3.3 further supporting a role for CHD2 in nucleosome turnover. We speculate that noncoding transcription at regulatory elements such as promoters, enhancers and active tRNA loci functions partially to modulate the chromatin environment through the recruitment and action of transcription coupled factors such as the remodeling enzymes CHD1 and CHD2.
Results and discussion
Recruitment of CHD1 and CHD2 to promoters is linked to the Pol II machinery
We analyzed CHD1 and CHD2 occupancy by calculating the enrichment of CHD1 and CHD2 from ChIP-seq data in relation to the chromatin states delineated by Ernst et al. . The chromatin states established by Ernst et al. were based on ChIP-seq data of combinatorial patterns of eight histone modifications and the insulator binding protein CTCF . The promoter and enhancer chromatin states are summarized in Table 1. In relation to promoter chromatin states; H3K4me2/3 and high H3K9/K27ac define the active promoter state. The weak promoter state is lower in H3K9/27 ac compared to the active promoter state, and H3K27me3 defines the inactive promoter. Enrichment of CHD1 and CHD2 was strongest at the active promoter chromatin state (Figure 1A). In all three cell types examined for both CHD1 and CHD2, the enrichment observed at active promoters was significantly higher than at weak or inactive promoters (Kruskall-Wallis, k = 3, P <0.0001). Previous studies have demonstrated in model organisms and human cells that CHD1 remodelers are recruited to promoters in a gene expression-dependent manner [29, 30]. Importantly, the weak promoter state from the ENCODE project is also defined by H3K4me2/me3 . This suggested that H3K4 methylation is therefore not sufficient for recruiting CHD1 to promoters in vivo since occupancy is strongly reduced or absent at weak promoters, which are marked by H3K4me2/3.
We hypothesized that the presence and absence of the Pol II transcription machinery occupancy may explain the difference in occupancy observed between the active and weak promoter states respectively. Consistent with this notion, Pol II was strongly enriched at active promoters only, whereas the H3K4me2/3 levels were comparable at both promoter states (Figure 1B). This was also apparent when inspecting individual gene promoters. In K562 cells at the GALK1 transcription start site (TSS), (Figure 1C; TSS 1), MNase-sequencing data and ChIP-seq data showed low nucleosome occupancy and high levels of H3K4me2/3, yet Pol II, CHD1 and CHD2 were not bound. In contrast, the TSS and gene body of the highly transcribed H3F3B (Figure 1C; TSS 2) gene were bound by Pol II and strongly enriched for the CHD1 and CHD2 remodelers. The UNK gene TSS (Figure 1C; TSS 3) showed lower levels of Pol II, CHD1 and CHD2 occupancy but abundant H3K4me2/3. Such examples clearly demonstrated that Pol II levels associate with CHD1 and CHD2 occupancy, indicating that transcription-coupled recruitment of CHD1 and CHD2 is important in targeting the enzymes to active promoters. Morettini et al. demonstrated that the chromodomains of Drosophila melanogaster CHD1 are critical for CHD1 function but not recruitment and concluded that ChIP-seq studies were needed to confirm this . Our analysis of human CHD1 provides support for this conclusion. Studies on CHD7, which is also reported to bind H3K4me2/3, revealed that CHD7 occupancy occurs on a fraction of H3K4me2/3 sites [31, 32]. These studies and the CHD7 data available in the ENCODE project suggest that CHD7, like CHD1, is recruited to a fraction of H3K4me2/3 sites independently of the chromodomain-H3K4me2/3 interaction since at sites such as weak promoters marked by H3K4me3, CHD7 is not recruited (see Additional file 1).
CHD remodeling enzymes also contain nonsequence-specific DNA binding domains, which play important roles in the regulation of chromatin remodeling. Zentner et al. hypothesized that this may drive recruitment of the CHD family of remodeling enzymes to active promoters, where nucleosome disassembly provides naked DNA . In human cells, the majority of promoters have CpG islands (CGI), which consist of stretches of abundant CpG base pairs. Such promoters are MNase sensitive and are thought to have constitutively lower nucleosome occupancies relative to non-CGI promoters independently of transcription rates . Thus, in mammalian cells, nucleosome occupancy and transcription can be partially uncoupled since non-CGI promoters show low nucleosome occupancy only when highly transcribed, while CGI promoters are constitutively MNase sensitive . Grouping of CGI promoters into quintiles of ascending Pol II occupancy clearly demonstrated that CHD1 and CHD2 recruitment is correlated to Pol II occupancy at these sites (Figure 2). Furthermore, despite showing low nucleosome occupancy in the absence of transcription, CGI TSSs do not recruit CHD1 or CHD2; that is, the accessible DNA at these nucleosome-depleted regions (NDRs) is not sufficient for the recruitment.
Recruitment of CHD1 and CHD2 at enhancer-like chromatin and active tRNA genes
Another chromatin state strongly associated with CHD1 or CHD2 enrichment was enhancer chromatin state 4 (Figure 3A). The enhancer chromatin states defined by Ernst et al. were identified through combinations of H3K4-methylation including the enhancer-specific H3K4me1 and histone acetylation (Table 1, Figure 3B). We found that enhancer state 4, which is distinguished by higher H3K4me3 and H3K9/27 ac than other enhancer states, was associated with Pol II, CHD1 and CHD2 enrichment. Other enhancer chromatin states 5 to 7, defined by lower levels of H3K4 methylation and H3K9/27 ac, showed weaker enrichment for Pol II, CHD1 and CHD2. This suggested that transcription of enhancers is coupled to chromatin remodeling by CHD1 and CHD2. Furthermore it has been demonstrated that a fraction of Pol II binding sites overlap with a distinct set of RNA polymerase III binding sites including active tRNA genes [34, 35]. Such active tRNA genes are bound by BRF1, the transcription factor IIIB subunit . Consistent with this, we found CHD1 and CHD2 were enriched at BRF1 but not at the alternate Pol III transcription factor, BRF2 binding sites in K562 cells (see Additional file 2). In this scenario, BRF2 binding sites function as control Pol III sites that do not overlap with Pol II binding. This suggests CHD1 and CHD2 may function at a subset of Pol III genes where Pol II transcription occurs, such as tRNAs.
CHD1 and CHD2 show overlapping occupancy but discrete binding profiles at transcription start sites
To validate the finding that CHD1 and CHD2 both prefer active promoter and enhancer chromatin states we identified ChIP-seq peaks using model-based analysis for ChIP-seq (MACS)  with a stringent cut-off p value of 1 × 10-15. For both data sets, peaks were merged if they fell within approximately one nucleosome length from another (200 bp). This produced 9,827 CHD1 peaks and 7,736 CHD2 peaks. For CHD1 and CHD2 respectively, 94% and 78% of high confidence peaks fell within a region associated with the active promoter or enhancer chromatin state 4 or both, based on chromatin states defined by Ernst et al. (Figure 4A, Table 1). Since chromatin states are continuous, a given peak may fall across more than one chromatin state.
The majority of CHD1 and CHD2 sites overlapped (60%) (Figure 4B). Despite this, the exact binding profiles of CHD1 and CHD2 at transcription start sites of genes in the top expression quintile differed (Figure 4C). CHD1 was enriched in a manner mimicking its preferred H3K4me2/3 nucleosomal substrates, while CHD2 bound across the NDR (Figure 4C). For CHD1, we propose that transcription-coupled recruitment with the Pol II machinery at the pre-initiation complex stage functions as the primary recruitment mechanism. In support of this, it is well established that CHD1 interacts with multiple Pol II-associated complexes [9, 22, 23, 38, 39]. Following recruitment, CHD1 then presumably interacts preferentially with H3K4me2/3 marked nucleosomes. In contrast, CHD2 does not have any specific affinity for H3K4me2/3 and binds strongest at the NDR. An alternative possibility is that differences in cross-linking preferences could affect the results obtained from conventional ChIP-seq. Future studies could address this by using higher resolution techniques, such as those employed by Skene et al., for murine Chd1 or using ChIP-exo to gain a more exact chromatin remodeler-DNA footprint [29, 40].
To test if transcription was linked to CHD1 and CHD2 recruitment, K562 cells were incubated with the Pol II CTD kinase inhibitor 5,6-dichloro-1-β-D-ribofuranosyl-1H-benzimidazole (DRB) to block elongation and synchronize transcription at the initiation stage (Figure 5A). This led to an accumulation of Pol II and a reduced H3 occupancy at the active NPM1 gene promoter as measured by ChIP-qPCR (Figure 5B, Students t-test, ***P <0.001). The inactive PRSS1 gene promoter was used as a control. Consistent with the proposed model of Pol II-driven recruitment, the binding of both CHD1 and CHD2 to the NPM1 gene promoter was significantly increased (Figure 5B, Students t-test ***P <0.001). Thus, we conclude that CHD1 and CHD2 are recruited with the Pol II machinery. The Henikoff laboratory recently demonstrated that a catalytically inactive, ATPase mutant of murine Chd1 reduces Pol II release from promoters, resulting in increased Pol II occupancy at TSSs. Consistent with our data (Figure 5B) increased Pol II occupancy in cells expressing mutant Chd1 was associated with increased Chd1 occupancy directly across the TSS . Our findings are also consistent with a recently published study showing that budding yeast Chd1 more closely matches Pol II phosphorylated at Serine 5 across the 5’ regions of genes than it does with the elongating form phosphorylated at Serine 2 . This is the case for both human CHD1 and CHD2, which bind across the 5’ region of genes similar to total Pol II ChIP-seq data while serine 2 phosphorylated Pol II shows stronger enrichment at the 3’ end of genes (see Additional file 3).
CHD1 affects chromatin accessibility at transcription start, intragenic and intergenic sites
To begin to dissect the functions of CHD1 and CHD2 and to confirm that they target active loci such as promoters and enhancers, we performed siRNA knockdowns in K562 cells combined with DNase accessibility assays. Effective knockdown of CHD1 and CHD2 was demonstrated at both mRNA and protein levels (Additional file 4).
DNase hypersensitivity sites (DHS) define regions of accessible chromatin associated with TSSs of expressed genes, while intergenic and intragenic DHS represent potential enhancer regions associated with transcription factor binding sites . We analyzed ENCODE DHS associated with the TSSs of highly expressed genes in K562 cells (n = 15) as well as intergenic and intragenic DHS sites (n = 15) and expressed tRNA genes (n = 6), which overlap DHS sites. Depletion of CHD1 with siRNA reduced the accessibility of all DHS types tested, while CHD2 knockdown did not (Figure 6A-B, **P <0.01, ***P >0.001, one-way ANOVA). Inaccessible centromeric repeats were used as a negative control and remained inaccessible (Figure 6A). The decrease in accessibility of DHS sites in CHD1 siRNA treated K562 cells was reduced to less than 10% to more than 90% of control (Figure 6B). Yet, since all DHS sites tested remained accessible to DNase to some degree, other factors must operate independently of CHD1 to promote accessible chromatin. Furthermore since there are on average more than 200,000 DHS loci per cell type , and we only tested 36 DHS sites, this study is not comprehensive. It does suggest, however, that CHD1 co-operates with other factors to promote accessible chromatin at regulatory regions.
DHS are embedded within chromatin marked by H3K4me2/3 at TSS and at intergenic and intragenic DHS sites including enhancers [43–46]. One possible function of CHD1 in light of the role in promoting chromatin accessibility could be to facilitate the binding of other transcription factors and complexes. This has been demonstrated in the case of post-initiation factors such as the splicing machinery, whereby CHD1 promotes efficient recruitment . Interestingly accessible chromatin is associated with both nucleosome-free and nucleosome containing regions and one of the strongest correlations between chromatin accessibility and histone modifications is the positioning and strength of the H3K4me3 nucleosomes, which may act to promote CHD1-mediated mobilization or relaxation of nucleosomes at regulatory regions .
CHD1 and CHD2 differentially affect H3 and H3.3 occupancy
One possible explanation for the effect of CHD1 on chromatin accessibility could be through the disassembly of nucleosomes. In fission yeast, the CHD1 homologues Hrp1 and Hrp3 function to reduce H3 occupancy at promoter regions [10, 19]. In addition, the H3 variant H3.3 may act to promote chromatin accessibility in some circumstances, and CHD1 homologues in model organisms and mammalian CHD2 have both been linked to H3.3 deposition [12, 16].
To analyze H3 occupancy and H3.3 occupancy by ChIP-qPCR, we first tested the specificity of the total H3 and H3.3 antibodies with recombinant human H3.1 and H3.3 proteins. The anti-H3 antibody detected both H3.1 and H3.3, while the H3.3 specific antibody detected only recombinant H3.3 protein by western blotting (see Additional file 5). We then tested the performance of these antibodies by ChIP-qPCR. When normalizing to total H3 occupancy, which is reduced at active regions compared to inactive regions, the relative enrichment of H3.3 was higher at the active GAPDH promoter compared to inactive MYOD promoter (see Additional file 5). Similarly, at the NPM1 gene, the H3.3/H3 total ratio was higher at the promoter compared to an intron from the gene body (see Additional file 4).
We selected TSS, intragenic, intergenic and tRNA-associated DHS sites, as well as a number of non-DHS control loci. for analysis by ChIP-qPCR (Figure 7). At TSSs, CHD1 knockdown lead to a nonsignificant trend toward increased H3 occupancy, while CHD2 knockdown lead to a striking increase in total H3 occupancy (Figure 7A, ***P > 0.001, one-way ANOVA). Double CHD1 and CHD2 knockdowns showed an additive effect with an elevated H3 occupancy compared to the single CHD2 knockdown (Figure 7A, 3.8% input compared to 4.3% input). In contrast the relative H3.3 enrichment observed in CHD2 and double CHD1 and CHD2 knockdowns was significantly lower than that of control and CHD1 knockdowns (Figure 7A, ***P <0.001, **P <0.01, one-way ANOVA). Similar results were obtained at intragenic and intergenic DHS sites with CHD2 knockdown giving significantly elevated H3 occupancy and reduced H3.3 enrichment (Figure 7B, *P <0.05, ****P <0.0001, one-way ANOVA). Similar trends were also observed at active tRNA loci associated DHS for both H3 and H3.3 (H3 ***P <0.001, H3.3, n.s. one-way ANOVA).
Surprisingly, DNase accessibility, H3 occupancy and H3.3 enrichment do not appear to be directly related in this study since, CHD1 affects accessibility strongly with weaker affects on H3 occupancy, while CHD2 affects H3 occupancy and H3.3 enrichment without affecting accessibility. Nucleosome occupancy is not the only factor that regulates chromatin accessibility. Accessible chromatin sites at intergenic loci for example show higher nucleosome occupancy and smaller nucleosome-free regions than transcription start sites but are still highly accessible [48, 49]. Accessible regulatory regions are composed of both nucleosome-free transcription factor binding sites and adjacent accessible nucleosomes [48, 49]. Hormone-inducible transcription factors bind to nucleosome-containing, accessible sites and are associated with chromatin remodeling but not necessarily nucleosome eviction [43, 50]. For example, nucleosomes are present at accessible intergenic progesterone receptor binding sites; the ensuing progesterone receptor activation accessibility is increased while total H3 occupancy remains unaltered . It is also possible that CHD1 depletion leads to abnormalities in higher order chromatin structure in a similar manner to what was demonstrated in Drosophila polytene chromosomes, following CHD1 deletion .
It is perhaps less surprising that the strong increase in H3 occupancy does not restrict chromatin accessibility since previous studies in fission yeast revealed that hrp1 and hrp3 deletions increase promoter H3 occupancy but do not affect MNase sensitivity of promoter regions [10, 19]. One possibility is that CHD2 reduces interactions of non-nucleosomal H3 with naked DNA. It also suggests that in the absence of CHD2, other factors prevent a mature and access restricting nucleosome from forming. Reduced H3.3 occupancy following CHD2 knockdown is likely a result of reduced CHD2-dependent nucleosome disassembly since it has been clearly demonstrated that H3.3 deposition correlates with histone turnover at active loci .
Implications for the biological functions of CHD1 and CHD2
Our analysis of ChIP-seq data for the CHD1 and CHD2 in human cells helps to reinterpret previously suggested models of CHD1 and CHD2 recruitment. We show that the recruitment is transcription-coupled and not driven by H3K4me2/3 or binding to NDR regions. The primary mechanisms of CHD1 recruitment to promoters and active chromatin regions must be driven by interactions between CHD1 and the mediator and must facilitate chromatin transcription (FACT); polymerase associated factor (PAF); and Spt-, Ada-, and Gcn5 acetyltransferase (SAGA) complexes [22, 38, 39]. The chromodomains of CHD1 instead act in an important post-recruitment step in the activation and regulation of remodeling. Consistent with this notion, chromodomain-nucleosome interactions of CHD1 regulate DNA-dependent activation of the ATPase motor .
An important role of CHD1 is to facilitate Pol II transcription through chromatin, regulating nucleosome turnover and transcription coupled chromatin assembly [9, 18, 29]. The affinity of human CHD1 for H3K4me2/3 may target the remodeling function of CHD1 to the 5’ region of active genes. For example it was recently suggested that mammalian CHD1 plays a key role in promoting Pol II escapes from active promoters by alleviating the barrier of the +1 nucleosome .
CHD2 has been linked to genetic stability, DNA repair and differentiation. In muscle cell differentiation, the myogenic transcription factor MyoD requires CHD2 for H3.3 deposition into muscle specific gene promoters. The strongest effects were likely seen at muscle-specific genes in this study since these regions are heavily transcribed in the model used . We observe that the H3.3 enrichment is not specific to muscle cell differentiation and speculate that CHD2 is an essential enzyme for the disassembly of H3 nucleosomes, which is a pre-requisite for the deposition of new H3.3 during differentiation. In other cell types it is possible that other transcription factors promote the recruitment of CHD2 in an analogous manner to MyoD. The effect of CHD1 and CHD2 in particular, on H3.3 occupancy suggest they may also play a role in the reassembly of H3.3 nucleosomes following transcription.
The transcription coupled recruitment of both CHD1 and CHD2 suggest that transcription could also play a role in noncoding regions of the genome such as promoters, enhancers, and active tRNA loci. Previous reports have suggested such functions for Pol II transcription at tRNA genes, leading to the hypothesis that Pol III transcription is facilitated by the Pol II associated machinery [34, 35]. Active enhancers are also identified by Pol II transcription . In such instances, we speculate transcription may be a mechanism for modulating chromatin structure. It was demonstrated that another transcription-coupled, SNF2 chromatin-remodeling enzyme Brg1 (SMARCA4), was required for enhancer activation of the MYC gene in leukemia . Another key question is how multiple remodeling factors determine the activity of each particular locus. A recent study examining distinct remodeling enzyme families demonstrated overlapping binding of multiple remodelers to active regions . Despite redundancy at the level of genome-wide occupancy, inactive mutants of the murine chromatin remodelers Brg1, Chd4 and hSNF2 all regulated chromatin accessibility at a subset of loci .
We conclude that CHD1 and CHD2 are recruited along with the transcription machinery to promoters and enhancer sites associated with transcription. For CHD1, this establishes a clear model of recruitment via interactions with Pol II-associated complexes. This provides a basis for interpretation of previous studies and models of CHD1 recruitment. At active regions, CHD1 and CHD2 cooperate to regulate the chromatin architecture. Thus, transcription coupled chromatin remodeling promotes chromatin accessibility and nucleosome disassembly at transcribed regulatory regions.
Analysis of ENCODE data
Data from the ENCODE project were accessed and downloaded through the UCSC web portal (http://genome.ucsc.edu/ENCODE/). The files used in this study are listed in Additional file 6. ChIP-seq and chromatin state data were analyzed using Seqmonk (http://www.bioinformatics.babraham.ac.uk/projects/seqmonk/). All ChIP-seq data, including CHD1 and CHD2 occupancy, was measured as the fold enrichment over input. Where indicated in the text, differences between chromatin states were tested for significance using the Kruskall-Wallis test with Dunn’s post-test correction. High confidence CHD1 and CHD2 peaks were calculated using model-based analysis for ChIP-seq (MACS) with a stringency cut off of 1 × 10-15 and binding sites within 200 bp or less were merged . MACS was used from Seqmonk version 0.27.0. For probe trend plots calculating the enrichment of H3k4me2/3 and Pol II over active and weak promoters and CHD1 or CHD2 occupancy across TSS regions, read counts were normalized to the total number of reads and the presented as the average signal per probe.
Cell culture and drug treatment
K562 cells were cultured in RPMI 1640 medium containing glutamine (Life Technologies, Stockholm, Sweden), which was supplemented with 10% fetal bovine serum, 10 U/mL penicillin and 10 U/mL streptomycin (Life Technologies). To arrest transcriptional elongation K562 cells were treated overnight (12 h) with 100 μM 5,6-dichloro-1-β-D-ribofuranosyl-1H-benzimidazole (DRB) (Sigma-Aldrich, Stockholm, Sweden) or vehicle (DMSO).
siRNA knockdown and validation
siRNA targeting CHD1 (s2974 and s2975) and CHD2 (s2979), as well as nontargeting negative controls, were obtained from Life Technologies (Stockholm, Sweden. K562 cells were transfected with the Neon electroporation system (Life Technologies, Stockholm, Sweden; 100 ul kit). The transfection protocol was performed as per manufacturer’s instructions with the following electroporation parameters: three pulses, 1,450 v, and 10 ms pulse width. The final total concentration of siRNA per sample was 25 nM in 8 mL of antibiotic-free RPMI 1640 medium with 10% serum per 3 × 106 cells.
RNA from 5 × 106 K562 cells was isolated using the RNeasy plus kit (QIAGEN, Sollentuna, Sweden). A total of 1 μg of RNA was reverse-transcribed using Superscript™ III First Strand Synthesis Supermix (Life Technologies, Stockholm, Sweden) according to the manufacturer’s instructions. Resultant cDNA was diluted 1:20, and qPCR analysis of CHD1 and CHD2 transcript abundance was normalized to GAPDH using Taqman assays (Life Technologies, Stockholm, Sweden; CHD1 assay ID - Hs00154405_m1; CHD2 assay ID - Hs00172280_m1; GAPDH assay ID - Hs02758991_g1).
For western blotting, 5 × 106 K562 were lyzed using single detergent lysis buffer (50 mM Tris, 150 mM NaCl, and 1% Triton X-100, pH 8.0) and protein concentration was calculated using the BCA protein assay (Fisher Scientific, Stockholm, Sweden). A total of 40 μg of protein per sample was denatured and resolved on a 4 to 12% Bis-Tris gel with MOPS running buffer (Novex, NuPAGE, precast gels, Life Technologies, Stockholm, Sweden). Gels were then wet-transferred to a 0.2-μM PVDF membrane (Millipore, Solna, Sweden) at a constant voltage of 100 V for 120 min. Membranes were air-dried, reactivated for 5 min in methanol and then blocked for 4 h at room temperature in 5% milk in PBS with 0.1% Tween. Membranes were incubated overnight in primary antibodies at 1:1,000 dilution in blocking buffer then washed three times for 10 min each in PBS 0.1% Tween and incubated for 1 h at room temperature in secondary antibody (GE Healthcare, Stockholm, Sweden). Finally, the PBS-Tween washes were repeated, and the membranes were developed with chemiluminescent ECL substrates. Antibodies used were as follows: CHD1 (#4351, Cell Signalling Technology, Danvers, MA, USA), CHD2 (#4170 Cell Signalling, Danvers, MA, USA), and H3 (ab1791, Abcam, Cambridge U.K.), all at a concentration of 1:1,000.
To measure chromatin accessibility, we used the EpiQ™ chromatin analysis kit (BIORAD, Stockholm, Sweden) according to the manufacturer’s instructions at 28 hours post transfection with siRNA. Briefly, 0.25 × 106 K562 cells were used per reaction and the resultant DNA was diluted to 1 in 400 in qPCR reactions. Accessibility was calculated as; 2(Ct(undigested)-Ct(digested)). We then normalized each loci to the percentage of accessibility in control siRNA treated cells. Primers were designed to DNase hypersensitivity peaks in K562 cells based on ENCODE data. For TSS, DHS regions we focused on expressed genes, whereas for intragenic and intergenic regions, we selected DHS sites observed when looking for accessible and inaccessible TSS regions. For tRNA loci, we chose loci from expressed tRNA genes that did not overlap with an annotated TSS but did overlap with a DHS. Primers are listed in Additional file 7.
K562 cells were fixed in 1% paraformaldehyde for 8 minutes, and ChIP was performed using the iDeal ChIP-seq kit (Diagenode, Liege, Belgium) according to the manufacturer’s instructions. Chromatin extracts from one million cells were sonicated to yield 100 to 500 bp fragments using a Bioruptor™ Plus (Diagenode) on high power for two rounds of 12 cycles with 30s on 30s off at 4°C. Antibodies used in ChIP were as follows; CHD1 - 3 μl per reaction of #4351 (Cell Signalling Technology, Danvers, MA, USA), CHD2 - 3 μl per reaction of ab68301 (Abcam, Cambridge U.K.), H3 - 2 μl per reaction of ab1791 (Abcam, Cambridge U.K.), and H3.3 - 3 μl per reaction of 09–838 (Millipore, Solna, Sweden). Primers for qPCR were designed across DNase hypersensitivity sites at TSS regions, intergenic or intragenic regions, or active tRNA loci and are listed in Additional file 7: Table S2.
chromodomain helicase DNA binding protein
chromatin immunoprecipitation high throughput sequencing
deoxyribonucleic acid I
DNase hypersensitive site
Encyclopedia of DNA Elements
embryonic stem cell
model-based analysis for ChIP-seq
- Pol II:
RNA polymerase II
transcription start site.
Narlikar GJ, Sundaramoorthy R, Owen-Hughes T: Mechanisms and functions of ATP-dependent chromatin-remodeling enzymes. Cell 2013,154(3):490–503. 10.1016/j.cell.2013.07.011
Hargreaves DC, Crabtree GR: ATP-dependent chromatin remodeling: genetics, genomics and mechanisms. Cell Res 2011,21(3):396–420. 10.1038/cr.2011.32
Flaus A, Martin DM, Barton GJ, Owen-Hughes T: Identification of multiple distinct Snf2 subfamilies with conserved structural motifs. Nucleic Acids Res 2006,34(10):2887–905. 10.1093/nar/gkl295
Murawska M, Brehm A: CHD chromatin remodelers and the transcription cycle. Transcription 2011,2(6):244–53. 10.4161/trns.2.6.17840
Marfella CG, Imbalzano AN: The Chd family of chromatin remodelers. Mutat Res 2007,618(1–2):30–40.
Burkhardt L, Fuchs S, Krohn A, Masser S, Mader M, Kluth M, et al.: CHD1 is a 5q21 tumor suppressor required for ERG rearrangement in prostate cancer. Cancer Res 2013,73(9):2795–805. 10.1158/0008-5472.CAN-12-1342
Gaspar-Maia A, Alajem A, Polesso F, Sridharan R, Mason MJ, Heidersbach A, et al.: Chd1 regulates open chromatin and pluripotency of embryonic stem cells. Nature 2009,460(7257):863–8.
Huang S, Gulzar ZG, Salari K, Lapointe J, Brooks JD, Pollack JR: Recurrent deletion of CHD1 in prostate cancer with relevance to cell invasiveness. Oncogene 2012,31(37):4164–70. 10.1038/onc.2011.590
Lin JJ, Lehmann LW, Bonora G, Sridharan R, Vashisht AA, Tran N, et al.: Mediator coordinates PIC assembly with recruitment of CHD1. Genes Dev 2011,25(20):2198–209. 10.1101/gad.17554711
Pointner J, Persson J, Prasad P, Norman-Axelsson U, Strålfors A, Khorosjutina O, et al.: CHD1 remodelers regulate nucleosome spacing in vitro and align nucleosomal arrays over gene coding regions in S. pombe. EMBO J 2012,31(23):4388–403. 10.1038/emboj.2012.289
Baca SC, Prandi D, Lawrence MS, Mosquera JM, Romanel A, Drier Y, et al.: Punctuated evolution of prostate cancer genomes. Cell 2013,153(3):666–77. 10.1016/j.cell.2013.03.021
Harada A, Okada S, Konno D, Odawara J, Yoshimi T, Yoshimura S, et al.: Chd2 interacts with H3.3 to determine myogenic cell fate. EMBO J 2012,31(13):2994–3007. 10.1038/emboj.2012.136
Nagarajan P, Onami TM, Rajagopalan S, Kania S, Donnell R, Venkatachalam S: Role of chromodomain helicase DNA-binding protein 2 in DNA damage response signaling and tumorigenesis. Oncogene 2009,28(8):1053–62. 10.1038/onc.2008.440
Marfella CG, Ohkawa Y, Coles AH, Garlick DS, Jones SN, Imbalzano AN: Mutation of the SNF2 family member Chd2 affects mouse development and survival. J Cell Physiol 2006,209(1):162–71. 10.1002/jcp.20718
Marfella CG, Henninger N, LeBlanc SE, Krishnan N, Garlick DS, Holzman LB, et al.: A mutation in the mouse Chd2 chromatin remodeling enzyme results in a complex renal phenotype. Kidney Blood Press Res 2008,31(6):421–32. 10.1159/000190788
Konev AY, Tribus M, Park SY, Podhraski V, Lim CY, Emelyanov AV, et al.: CHD1 motor protein is required for deposition of histone variant H3.3 into chromatin in vivo. Science 2007,317(5841):1087–90. 10.1126/science.1145339
Lusser A, Urwin DL, Kadonaga JT: Distinct activities of CHD1 and ACF in ATP-dependent chromatin assembly. Nat Struct Mol Biol 2005,12(2):160–6. 10.1038/nsmb884
Radman-Livaja M, Quan TK, Valenzuela L, Armstrong JA, van Welsem T, Kim T, et al.: A key role for Chd1 in histone H3 dynamics at the 3' ends of long genes in yeast. PLoS Genet 2012,8(7):e1002811. 10.1371/journal.pgen.1002811
Walfridsson J, Khorosjutina O, Matikainen P, Gustafsson CM, Ekwall K: A genome-wide role for CHD remodelling factors and Nap1 in nucleosome disassembly. EMBO J 2007,26(12):2868–79. 10.1038/sj.emboj.7601728
Sims RJ, Chen CF, Santos-Rosa H, Kouzarides T, Patel SS, Reinberg D: Human but not yeast CHD1 binds directly and selectively to histone H3 methylated at lysine 4 via its tandem chromodomains. J Biol Chem 2005,280(51):41789–92. 10.1074/jbc.C500395200
Morettini S, Tribus M, Zeilner A, Sebald J, Campo-Fernandez B, Scheran G, et al.: The chromodomains of CHD1 are critical for enzymatic activity but less important for chromatin localization. Nucleic Acids Res 2011,39(8):3103–15. 10.1093/nar/gkq1298
Kelley DE, Stokes DG, Perry RP: CHD1 interacts with SSRP1 and depends on both its chromodomain and its ATPase/helicase-like domain for proper association with chromatin. Chromosoma 1999,108(1):10–25. 10.1007/s004120050347
Simic R, Lindstrom DL, Tran HG, Roinick KL, Costa PJ, Johnson AD, et al.: Chromatin remodeling protein Chd1 interacts with transcription elongation factors and localizes to transcribed genes. EMBO J 2003,22(8):1846–56. 10.1093/emboj/cdg179
Shema-Yaacoby E, Nikolov M, Haj-Yahya M, Siman P, Allemand E, Yamaguchi Y, et al.: Systematic identification of proteins binding to chromatin-embedded ubiquitylated H2B reveals recruitment of SWI/SNF to regulate transcription. Cell Rep 2013,4(3):601–8. 10.1016/j.celrep.2013.07.014
Smolle M, Venkatesh S, Gogol MM, Li H, Zhang Y, Florens L, et al.: Chromatin remodelers Isw1 and Chd1 maintain chromatin structure during transcription by preventing histone exchange. Nat Struct Mol Biol 2012,19(9):884–92. 10.1038/nsmb.2312
Khorosjutina O, Wanrooij PH, Walfridsson J, Szilagyi Z, Zhu X, Baraznenok V, et al.: A chromatin-remodeling protein is a component of fission yeast mediator. J Biol Chem 2010,285(39):29729–37. 10.1074/jbc.M110.153858
Flanagan JF, Blus BJ, Kim D, Clines KL, Rastinejad F, Khorasanizadeh S: Molecular implications of evolutionary differences in CHD double chromodomains. J Mol Biol 2007,369(2):334–42. 10.1016/j.jmb.2007.03.024
Ernst J, Kheradpour P, Mikkelsen TS, Shoresh N, Ward LD, Epstein CB, et al.: Mapping and analysis of chromatin state dynamics in nine human cell types. Nature 2011,473(7345):43–9. 10.1038/nature09906
Skene PJ, Hernandez AE, Groudine M, Henikoff S: The nucleosomal barrier to promoter escape by RNA polymerase II is overcome by the chromatin remodeler Chd1. Elife 2014, 3:e02042.
Zentner GE, Tsukiyama T, Henikoff S: ISWI and CHD chromatin remodelers bind promoters but act in gene bodies. PLoS Genet 2013,9(2):e1003317. 10.1371/journal.pgen.1003317
Schnetz MP, Bartels CF, Shastri K, Balasubramanian D, Zentner GE, Balaji R, et al.: Genomic distribution of CHD7 on chromatin tracks H3K4 methylation patterns. Genome Res 2009,19(4):590–601. 10.1101/gr.086983.108
Schnetz MP, Handoko L, Akhtar-Zaidi B, Bartels CF, Pereira CF, Fisher AG, et al.: CHD7 targets active gene enhancer elements to modulate ES cell-specific gene expression. PLoS Genet 2010,6(7):e1001023. 10.1371/journal.pgen.1001023
Vavouri T, Lehner B: Human genes with CpG island promoters have a distinct transcription-associated chromatin organization. Genome Biol 2012,13(11):R110. 10.1186/gb-2012-13-11-r110
Barski A, Chepelev I, Liko D, Cuddapah S, Fleming AB, Birch J, et al.: Pol II and its associated epigenetic marks are present at Pol III-transcribed noncoding RNA genes. Nat Struct Mol Biol 2010,17(5):629–34. 10.1038/nsmb.1806
Listerman I, Bledau AS, Grishina I, Neugebauer KM: Extragenic accumulation of RNA polymerase II enhances transcription by RNA polymerase III. PLoS Genet 2007,3(11):e212. 10.1371/journal.pgen.0030212
Canella D, Praz V, Reina JH, Cousin P, Hernandez N: Defining the RNA polymerase III transcriptome: genome-wide localization of the RNA polymerase III transcription machinery in human cells. Genome Res 2010,20(6):710–21. 10.1101/gr.101337.109
Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, et al.: Model-based analysis of ChIP-Seq (MACS). Genome Biol 2008,9(9):R137. 10.1186/gb-2008-9-9-r137
Dey P, Ponnusamy MP, Deb S, Batra SK: Human RNA polymerase II-association factor 1 (hPaf1/PD2) regulates histone methylation and chromatin remodeling in pancreatic cancer. PLoS One 2011,6(10):e26926. 10.1371/journal.pone.0026926
Pray-Grant MG, Daniel JA, Schieltz D, Yates JR, Grant PA: Chd1 chromodomain links histone H3 methylation with SAGA- and SLIK-dependent acetylation. Nature 2005,433(7024):434–8. 10.1038/nature03242
Yen K, Vinayachandran V, Batta K, Koerber RT, Pugh BF: Genome-wide nucleosome specificity and directionality of chromatin remodelers. Cell 2012,149(7):1461–73. 10.1016/j.cell.2012.04.036
Park D, Shivram H, Iyer VR: Chd1 co-localizes with early transcription elongation factors independently of H3K36 methylation and releases stalled RNA polymerase II at introns. Epigenetics Chromatin 2014, 7:32. doi:10.1186/1756–8935–7-32 10.1186/1756-8935-7-32
Thurman RE, Rynes E, Humbert R, Vierstra J, Maurano MT, Haugen E, et al.: The accessible chromatin landscape of the human genome. Nature 2012,489(7414):75–82. 10.1038/nature11232
He HH, Meyer CA, Chen MW, Jordan VC, Brown M, Liu XS: Differential DNase I hypersensitivity reveals factor-dependent chromatin dynamics. Genome Res 2012,22(6):1015–25. 10.1101/gr.133280.111
Pekowska A, Benoukraf T, Zacarias-Cabeza J, Belhocine M, Koch F, Holota H, et al.: H3K4 tri-methylation provides an epigenetic signature of active enhancers. EMBO J 2011,30(20):4198–210. 10.1038/emboj.2011.295
Kowalczyk MS, Hughes JR, Garrick D, Lynch MD, Sharpe JA, Sloane-Stanley JA, et al.: Intragenic enhancers act as alternative promoters. Mol Cell 2012,45(4):447–58. 10.1016/j.molcel.2011.12.021
Djebali S, Davis CA, Merkel A, Dobin A, Lassmann T, Mortazavi A, et al.: Landscape of transcription in human cells. Nature 2012,489(7414):101–8. 10.1038/nature11233
Sims RJ, Millhouse S, Chen CF, Lewis BA, Erdjument-Bromage H, Tempst P, et al.: Recognition of trimethylated histone H3 lysine 4 facilitates the recruitment of transcription postinitiation factors and pre-mRNA splicing. Mol Cell 2007,28(4):665–76. 10.1016/j.molcel.2007.11.010
Vierstra J, Wang H, John S, Sandstrom R, Stamatoyannopoulos JA: Coupling transcription factor occupancy to nucleosome architecture with DNase-FLASH. Nat Methods 2014,11(1):66–72.
Buenrostro JD, Giresi PG, Zaba LC, Chang HY, Greenleaf WJ: Transposition of native chromatin for fast and sensitive epigenomic profiling of open chromatin, DNA-binding proteins and nucleosome position. Nat Methods 2013,10(12):1213–8. 10.1038/nmeth.2688
Ballaré C, Castellano G, Gaveglia L, Althammer S, González-Vallinas J, Eyras E, et al.: Nucleosome-driven transcription factor binding and gene regulation. Mol Cell 2013,49(1):67–79.
Bugga L, McDaniel IE, Engie L, Armstrong JA: The Drosophila melanogaster CHD1 chromatin remodeling factor modulates global chromosome structure and counteracts HP1a and H3K9me2. PLoS One 2013,8(3):e59496. 10.1371/journal.pone.0059496
Ahmad K, Henikoff S: The histone variant H3.3 marks active chromatin by replication-independent nucleosome assembly. Mol Cell 2002,9(6):1191–200. 10.1016/S1097-2765(02)00542-7
Hauk G, McKnight JN, Nodelman IM, Bowman GD: The chromodomains of the Chd1 chromatin remodeler regulate DNA access to the ATPase motor. Mol Cell 2010,39(5):711–23. 10.1016/j.molcel.2010.08.012
Andersson R, Gebhard C, Miguel-Escalada I, Hoof I, Bornholdt J, Boyd M, et al.: An atlas of active enhancers across human cell types and tissues. Nature 2014,507(7493):455–61. 10.1038/nature12787
Shi J, Whyte WA, Zepeda-Mendoza CJ, Milazzo JP, Shen C, Roe JS, et al.: Role of SWI/SNF in acute leukemia maintenance and enhancer-mediated Myc regulation. Genes Dev 2013,27(24):2648–62. 10.1101/gad.232710.113
Morris SA, Baek S, Sung MH, John S, Wiench M, Johnson TA, et al.: Overlapping chromatin-remodeling systems collaborate genome wide at dynamic chromatin transitions. Nat Struct Mol Biol 2014,21(1):73–81.
We would like to thank JP Svensson, P Prasad and O Khorosjutina for fruitful discussions during this study and about the manuscript. The Ekwall laboratory is supported by grants from the Swedish Cancer Society (grant number CAN 2012/238), the Swedish Research Council, (grant numbers VR-M 2579 and VR-NT 4448) and the Knut and Alice Wallenberg foundation.
The authors declare that they have no competing interests.
LS performed the ENCODE data analysis, siRNA knockdown studies, designed the studies and wrote the manuscript. LC performed chromatin immunoprecipitation, siRNA transfections and RT-qPCR. MR performed the ENCODE data analysis. AL and KE designed the study, aided in the interpretation of the data and assisted in writing the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: CHD7 occupancy at promoter chromatin states in K562 and H1 embryonic stem cells (ESCs). The occupancy of CHD7 was calculated as the fold enrichment over input. (PDF 1 MB)
Additional file 2: Recruitment of CHD1 and CHD2 to active tRNA loci. (A) The occupancy of CHD1 and CHD2 across BRF1 and BRF2 binding sites in K562 cells was calculated as the fold enrichment over input. (B) A selected, expressed tRNA gene in an intergenic region showing BRF1, CHD1, CHD2 and Pol II occupancy in K562 cells. (PDF 1 MB)
Additional file 3: Average enrichment of Pol II, Pol II ser2, CHD1 and CHD2 at transcribed genes. Transcribed genes were identified as in the top quintile of Pol II ser 2 occupancy and analyzed for Pol II, CHD1 and CHD2 enrichment across the gene plus 2.5 kb up and downstream. (PDF 361 KB)
Additional file 4: CHD1 and CHD2 siRNA knockdown. (A) Transcript abundance of CHD1 and CHD2 relative to GAPDH in K562 cells 48 h following siRNA transfection. (B) Western blot analysis of CHD1 and CHD2 expression 72 h following transfection. (PDF 2 MB)
Additional file 5: Analysis of H3.3 antibody specificity by western blot and ChIP-qPCR. (A) The reactivity of total H3 and H3.3 specific antibodies toward 10 to 250 ng of recombinant human H3.1 and H3.3 was analyzed by western blot. (B) The relative enrichment of H3.3 to total H3 was examined by ChIP qPCR at active GAPDH and inactive MYOD gene TSS and the active NPM1 promoter and NPM1 intron. (PDF 1 MB)
About this article
Cite this article
Siggens, L., Cordeddu, L., Rönnerblad, M. et al. Transcription-coupled recruitment of human CHD1 and CHD2 influences chromatin accessibility and histone H3 and H3.3 occupancy at active chromatin regions. Epigenetics & Chromatin 8, 4 (2015). https://doi.org/10.1186/1756-8935-8-4