- Open Access
Local compartment changes and regulatory landscape alterations in histone H1-depleted cells
Genome Biology volume 16, Article number: 289 (2015)
Linker histone H1 is a core chromatin component that binds to nucleosome core particles and the linker DNA between nucleosomes. It has been implicated in chromatin compaction and gene regulation and is anticipated to play a role in higher-order genome structure. Here we have used a combination of genome-wide approaches including DNA methylation, histone modification and DNase I hypersensitivity profiling as well as Hi-C to investigate the impact of reduced cellular levels of histone H1 in embryonic stem cells on chromatin folding and function.
We find that depletion of histone H1 changes the epigenetic signature of thousands of potential regulatory sites across the genome. Many of them show cooperative loss or gain of multiple chromatin marks. Epigenetic alterations cluster to gene-dense topologically associating domains (TADs) that already showed a high density of corresponding chromatin features. Genome organization at the three-dimensional level is largely intact, but we find changes in the structural segmentation of chromosomes specifically for the epigenetically most modified TADs.
Our data show that cells require normal histone H1 levels to expose their proper regulatory landscape. Reducing the levels of histone H1 results in massive epigenetic changes and altered topological organization particularly at the most active chromosomal domains. Changes in TAD configuration coincide with epigenetic landscape changes but not with transcriptional output changes, supporting the emerging concept that transcriptional control and nuclear positioning of TADs are not causally related but independently controlled by the locally associated trans-acting factors.
DNA in the eukaryotic nucleus is packaged into arrays of nucleosome core particles that are the basic unit of chromatin [1, 2]. Each nucleosome consists of an octamer of four core histones (H2A, H2B, H3 and H4) around which about 145 bp of DNA is wrapped. Chromatin also contains a fifth histone, the linker histone, usually referred to as H1. H1 binds to nucleosome core particles near the DNA entry/exit position and to the linker DNA between core particles, stabilizing the association of core particles and DNA and facilitating folding of oligonucleosome arrays into compact structures. Mice and humans express 11 H1 subtypes, including H1a to H1e found at varying levels in most cell types, a replacement subtype (H1(0)) generally associated with terminal differentiation and quiescent states, four germ cell-specific H1s (H1t, H1T2, H1LS1, and H1oo) and a less well-studied subtype (H1x) [3–5]. In addition to their different developmental expression patterns and abundance, the amino acid sequences of these H1 subtypes differ significantly. Despite these differences, we [6–8] and others [9–11] found that elimination of any one of several H1 subtypes, and even some pairs of subtypes, does not noticeably affect mouse development. The absence of phenotypes in these mice appears to be due to up-regulation of the remaining subtypes, resulting in maintenance of a normal H1 to nucleosome core particle stoichiometry. However, elimination of three H1 subtypes (H1c, H1d, and H1e) together led to a 50 % reduction in the ratio of H1 to core particles and embryonic lethality. Embryonic stem (ES) cells derived from H1c, H1d, H1e null embryos are viable and also exhibit a 50 % reduction in H1:core histone stoichiometry. They show a decrease in the average spacing between nucleosome core particles of about 15 bp, from ~189 bp in normal ES cells to ~174 bp in the triple H1 knock-out (TKO) ES cells . These TKO ES cells also showed decreased local chromatin compaction and selective changes in gene expression. Importantly, up-regulation of certain imprinted and X chromosome-linked genes were prominent and found to be due to H1-dependent alterations in DNA methylation and H3 histone methylation at the regulatory regions of the affected gene [12–14]. Reduced H1 levels were also found to enable CTCF to bind to normally occluded DNA sequences at some imprinted gene loci . CTCF is a central factor in setting up local chromatin loops and defining structural domains across mammalian chromosomes [15, 16]. Collectively, this suggests that histone H1 may also have an important function in shaping higher order genome structures in vivo, either directly through its capacity to compact DNA or indirectly by controlling the DNA accessibility of chromatin architectural proteins.
To investigate the role of H1 in genome-wide, higher order chromatin structure, we studied the regulatory landscape and overall genome conformation in the H1-depleted TKO ES cells . ES cells show several unique features of nuclear organization compared with somatic cells. For example, ES cells display hypermobility of chromatin proteins, including the core histones and histone H1, indicative of their loosened binding to DNA . Restricting the dynamic state of these core chromatin components compromises the differentiation capacity of ES cells, suggesting that this feature is essential for ES cell identity [17, 18]. ES cells also have an unusually low H1 to nucleosome core stoichiometry: whereas this ratio is typically 0.75 or more in differentiated cells, in wild-type ES cells it is only about 0.5 . In the H1 TKO ES cells, this ratio is further reduced to one histone H1 molecule per four nucleosomes . ES cells also display a distinctively disorganized three-dimensional (3D) genome with particularly inactive chromosomal regions that fail to cluster as efficiently as seen in somatic cells . Here we applied Hi-C  and other genome-wide approaches for mapping epigenetic features to compare wild-type and H1-depleted ES cells to better understand how the H1 linker histone impacts the regulatory and 3D landscape of the genome and its consequences on transcription.
Results and discussion
Clustered DNA methylation changes in histone H1-depleted ES cells
Local DNA demethylation was previously observed in the H1-depleted TKO ES cells particularly at the imprinting control regions of H19-Igf2 and Gtl2-Dlk1 loci  and at the Rhox 5 promoter on the X chromosome . The activity of H1 in promoting DNA methylation at the imprinting control regions was later attributed to H1’s function in recruiting the DNA methyltransferases DNMT1 and DNMT3B . No global changes in DNA methylation were observed at that time, based on digestion with methylation-sensitive restriction enzymes . To study genome-wide methylation changes in more detail, here we used the HELP-tagging assay that enables high-throughput identification of sequences neighboring restriction sites of the methylation-sensitive HpaII enzyme . Approximately 15,000 sites displayed differential methylation across the genome. Consistent with a function of H1 in recruiting DNA methyltransferases , more than two-thirds of these sites appeared hypo-methylated in TKO cells, while ~30 % were hypermethylated. Hypermethylation often was not as pronounced as hypomethylation though, and when more stringent criteria were applied for selecting differentially methylated regions (DMRs), 4315 hypomethylated sites were found versus only 308 hypermethylated sites (93 % versus 7 % of all DMRs). The differentially methylated sites were not uniformly distributed across the genome but appeared significantly clustered (Fig. 1a, b), with the sex chromosomes being surprisingly protected against methylation changes (Figure S1 in Additional file 1). To further delineate the clustered distribution, we investigated the density of differentially methylated sites in topologically associating domains (TADs). We considered TADs as chromosomal units of interest as they mark genomic segments within which sequences preferentially contact each other. As such, TADs are believed to be the structural and functional genomic units that encompass the genes and their cognate regulatory sites [15, 24]. We categorized TADs according to gene content and created five bins with equal numbers of TADs; each bin had a similar average TAD size, but the gene-poorest bin encompassed only 3–4 % of all genes while the gene-richest bin carried over 50 % (Fig. 1c). Overall, DNA methylation density followed gene density (Fig. 1c). Taking this distribution into account we found that changes in DNA methylation significantly clustered in the most gene-rich TADs, no matter whether we considered hyper- or hypo-methylation events (Fig. 1d). Thus, bins in gene-dense TADs showed more changes in methylation than bins with equal CpG density in gene-poor TADs (Figure S2 in Additional file 1). Collectively, the data show that DNA methylation changes preferentially occur in chromosomal segments with a high gene density.
To investigate whether the DNA methylation changes associated with H1 depletion localized to any specific type of chromatin, we used published ES cell ChIP-seq profiles and a partitioning of the genome into different chromatin states by means of ChromHMM, a segmentation algorithm for the identification of chromatin states based on the presence of combinations of chromatin modifications . We found that particularly predicted enhancers were strongly overrepresented among the DMRs (Fig. 1e). CpG-rich gene promoters, on the other hand, appeared underrepresented among the hypomethylated sites (Fig. 1e).
In summary, a reduction in the amount of histone H1 protein results in abundant changes in DNA methylation, with some sites gaining but most sites losing methylation. Methylation changes accumulate at the most gene-dense TADs. They occur particularly at enhancer sequences, further indicating that histone H1 plays a role in controlling the DNA methylation status at potential regulatory sequences. In contrast, CpG-rich promoter sequences appear to stably maintain their methylation status in TKO cells, indicating that they control methylation levels in a histone H1-independent manner. A previous observation that H1c and H1d variants are depleted around the transcription start site of active promoters may be in agreement with this idea .
Altered genomic regulatory landscape in H1 TKO cells
To further characterize the consequences of histone H1 depletion on the regulatory chromatin landscape of ES cells, we determined the distribution of DNase hypersensitivity sites (DHSs), H3K4me1, H3K4me3, H3K27me3, and H3K9me3 sites across the genome of H1 TKO ES cell lines and their wild-type counterparts. Overall, the histone modification distributions looked similar between the wild-type, TKO cells and those published by ENCODE for another ES cell line (Fig. 2a; Figure S3 in Additional file 1). Roughly identical numbers of DHSs were scored in the two conditions (281,934 in TKO versus 293,319 in wild type, medians over triplicates). When ranked according to their delta DHS signal, de novo formed DHSs were clearly appreciable, but there was little evidence for complete loss of DHSs in H1 TKO cells (Fig. 2b). Chromatin-associated histone H1 was previously found to interfere with the binding of histone methyltransferase SET7/9, thereby preventing methylation of H3K4 . Statistical analysis of ChIP-seq differential enrichment (see “Materials and methods”) revealed no change in the overall numbers of H3K4me1 sites across the genome but did show a high number of sites that gained (6536) or lost (7319) mono-methylation, indicating dynamic changes in this enhancer mark (Figure S4a in Additional file 1). When looking at the H3K4me3 ChIP-seq results we found four times more sites with increased than with decreased levels of trimethyation (2043 versus 495) (Figure S4b in Additional file 1). Surprisingly, compared with the abundant changes observed for the active marks H3K4me1 and H3K4me3, changes in the repressive histone modifications H3K9me3 and H3K27me3 were almost negligible. When using the same fold change cutoffs as used to identify sites with an altered H3K4 methylation status, only a few dozen sites showed dynamic changes in the repressive H3K9me3 and H3K27me3 marks (not shown). A largely unaltered distribution of these repressive marks seems surprising given our earlier finding that histone H1 physically recruits the heterochromatin-specific histone H3 lysine 9 methyltransferase Su(var)3-9 in Drosophila  but is in agreement with our observations that the intranuclear distribution of histone marks H3K27me3/H3K9me2 and heterochromatin-associated factors such as HP1a, HP1b, and MeCP2 appeared normal by immunofluorescence .
We next wished to understand the relationship between these epigenetic changes. Since differences in DHSs were clearest for the 2123 newly formed DHSs, we focused on those DHSs and asked whether their formation coincided with other epigenetic changes. Interestingly, these sites were statistically significantly enriched (Figure S5 in Additional file 1) for the binding motifs of a number of pluripotency factors, including Klf4 (three-fold enrichment, as judged by HOMER ), but also Oct4 (two-fold) and Sox2 (two-fold). This suggests that histone H1 normally serves to occlude these sites, which may be in agreement with the earlier observation that wild-type H1 levels are necessary for normal ES cell differentiation and the concomitant repression of Oct4 expression . Nearly one-third of the new DHSs also showed a gain in either H3K4me1 (p < 1e-6, significance of overlap in hypergeometric test) or H3K4me3 (p < 1e-6) or both, whereas loss of these marks was very infrequently observed at new DHSs (Figure 2c; Additional file 2). More than 10 % (256/2123) of the new DHSs also revealed loss of DNA methylation, while the opposite, hypermethylation, was rarely found at these sites (19 times) (Figure S6 in Additional file 1). When focusing on H3K4me3 sites, those with increased H3k4me3 levels often (>25 %) also showed a gain (and seldom a loss) in H3K4me1, while sites losing H3K4me3 frequently showed concomitant loss of H3K4me1 (also >25 %) (Fig. 2c). Finally, when considering the differentially methylated CpGs, sites with reduced methylation in TKO cells were often enriched for H3K4me1 marks and, vice versa, hypermethylated sites frequently lost H3K4me1 (Fig. 2e).
In summary, the depletion of histone H1 has little impact on the genome-wide distribution of the repressive histone modifications H3K27me3 and H3K9me3, but changes the active chromatin H3K4me1 and H3K4me3 signatures of thousands of sites across the genome. Many of them show concomitant loss or gain of multiple chromatin marks associated with regulatory activity. Nearly invariably these combinatorial changes are either all positively or all negatively contributing to an active chromatin signature, implying that they can cooperate to strengthen or dampen the regulatory potential of a site.
Epigenetic changes accumulate in gene-dense TADs
We then wished to understand where in the genome these epigenetic changes take place. For this we again considered TADs as the genomic units of interest and we intersected the various datasets with the previously defined five classes of TADs. Not surprisingly, the general distribution of DHSs and H3K4me1 and H3K4me3 sites in wild-type (and TKO) cells closely followed that of genes, with all of these marks specifically accumulating in the most gene-dense TADs (Figure S7a in Additional file 1). The same was true for H3K27me3, while H3K9me3 distributed more equally across TADs with different gene density, albeit slightly accumulated in both the most gene-poor as well as the most gene-rich TADs (Figure S7b in Additional file 1). When correcting for their overall distribution there was no obvious enrichment of sites losing DHS or H3K4me3 signal in any of the TAD bins (Fig. 3a). This suggests that sites showing loss of hypersensitivity or loss of the promoter mark H3K4me3 are distributed proportional to the overall genomic localization of DHSs and H3K4me3 sites. In contrast, sites losing H3K4 monomethylation in H1 TKO cells were significantly depleted (chi-squared test p < 10−6) from the most gene-dense TADs and appeared to accumulate in the gene-poorest TADs (Fig. 3b). This could indicate that normal histone H1 levels are needed for proper maintenance of H3K4me1 levels in inactive chromatin surroundings. Alternatively, H3K4me1 sites in active chromatin surroundings are relatively protected against demethylation. To further investigate if epigenetic alterations occurred at specific genomic locations, we looked at de novo acquired active chromatin marks. We defined de novo DHSs as those that were exclusively identified in TKO cells but also lacked threshold H3K4me1 or H3K4me3 levels in wild-type cells. Similarly, we defined de novo formed H3K4me1 sites as those scoring positive for this mark only in TKO cells and also lacking significant H3K4me3 and DHS signal in wild-type cells. These non-bookmarked sites are ubiquitously present and their conversion into active sites may, therefore, a priori take place anywhere in the genome. However, new DHSs and new H3K4me1 sites both preferentially accumulated again in the most gene-dense TADs (Fig. 3c, d). Thus, despite being a generic chromatin component present throughout the genome, depletion of H1 results in a preferential gain of the active H3K4me1 and H3K4me3 chromatin marks within the most gene-dense TADs. These TADs are already dense in such regulatory chromatin signatures and this, we speculate, may create sensitized chromatin that is exquisitely susceptible to further epigenetic changes. An alternative, not mutually exclusive, explanation is that these TADs form nuclear compartments where the corresponding chromatin modifying enzymes accumulate to cooperatively establish and maintain the required dense landscape of regulatory sites. Perturbing the integrity of the chromatin template, as occurs when H1 histone levels are reduced, may further stimulate local mass action and increase the chance to modify neighboring chromatin sites. Precedence for local cooperative action exists: it was recently shown in Drosophila that clustered low affinity binding sites better accumulate PcG proteins than their more isolated counterparts elsewhere in the genome .
Genes with altered expression are proportionally distributed across the genome
To investigate how the altered regulatory chromatin landscape functionally translates into gene expression changes, we assayed the genome-wide transcriptome. RNA-seq confirmed previous observations obtained by microarray analysis. Transcription of the vast majority of genes is not affected in TKOs and the cells still clearly harbor an ES cell identity (Fig. 4a). We again found a small subset of genes (<200), apparently randomly distributed across all chromosomes, that showed a significant twofold or higher change in expression, the majority (>75 %) of which showed reduced levels of transcription (Fig. 4b). Among those were the previously described Hox genes , while the most prominently upregulated genes included a series of paternally imprinted genes  (Fig. 4c). The slight overrepresentation of X-linked genes that was previously apparent among 29 dysregulated genes  was no longer appreciable in this larger set of differentially expressed genes. Previous detailed characterization of two of the most strongly upregulated loci in TKO cells, the paternally imprinted Gtl2 locus and the H19 locus, revealed hypomethylation of their imprinting control regions . To investigate whether loss of DNA methylation generally underlies transcriptome changes we compared the genomic distribution of up- and down-regulated genes and differentially methylated sites at the level of TADs. To maximally exploit the benefit of an integrative analysis, we considered a less stringent set of 598 differentially expressed genes. We ranked TADs based on the number of DNA de-methylated sites and computed the fractions of differentially regulated genes. Figure 4d shows that indeed TADs with most changes in DNA methylation co-segregated with those most enriched for differentially expressed genes. However, given the non-uniform genomic distribution of differentially methylated sites over gene-dense TADs (Fig. 1d), we considered the overall distribution of genes to be a confounding factor here. To investigate this in more detail we ranked TADs according to gene content. Indeed, this categorization highly correlated with the distribution of differentially expressed genes (Fig. 4e), implying that, from a genomic distribution point of view, they are a proportional and apparently random collection of genes. Possibly in agreement with this, a gene ontology enrichment analysis on the set of differentially expressed genes did not reveal any specific gene ontology categories to be highly enriched. For the sites with changes in DNase I hypersensitivity, the analysis at TAD level is not really appropriate as they are too scarce in individual TADs, so instead we computed the percentages of genes where a significant change in DNase I hypersensitivity occurred within 2500 bp up- or downstream of the transcription start site of the gene. For up-regulated genes, we saw that 6.5 % and 3.2 % had a significant loss and gain, respectively, in DHSs. For down-regulated genes these percentages were 4.6 % and 2.4 %. Hence, the vast majority of differentially expressed genes did not show significant changes in DHSs at their promoters. In cases where changes were found, gain and loss of hypersensitivity was uncorrelated to up- and down-regulated gene activity. The same was true when comparing the distribution of differentially expressed genes to that of TADs with increased H3K4me1 signal (differentially expressed genes strictly followed the overall gene distribution; Fig. 4f), whereas no correlation was observed when comparing differentially expressed genes with TADs with decreased H3K4me1 signal in TKO cells (Fig. 4g).
In summary, while new epigenetic features de novo acquired in H1 TKO cells preferentially appeared in TADs that in wild-type cells already carried the highest density of these marks, these same TADs, when normalized for gene content, were not enriched for deregulated genes. The mere consideration of epigenetic changes is therefore not sufficient to predict gene expression changes. Rather than being dependent on the overall density of epigenetic marks in the TAD, the transcriptional output of individual genes is, therefore, more likely controlled by the regulatory status of only a limited collection of regulatory modules that, possibly via chromatin looping, act on the target gene.
Higher-order topological changes follow epigenetic but not transcriptional changes
We finally wished to understand how histone H1 depletion and the accompanying epigenome and transcriptome changes affect the overall 3D organization of the genome. To this end we performed replicate Hi-C experiments, each with a different, frequently cutting, restriction enzyme (NlaIII and DpnII) , in both TKO and matched wild-type ES cells. Each dataset contained between 26 and 42 million valid Hi-C read pairs, adding up to 53 milliom (wild-type) and 76 million (TKO) valid Hi-C read pairs per cell type. All Hi-C libraries showed an equal high ratio of intra- over inter-chromosomal contacts (~75 %), indicative of good quality Hi-C libraries . We normalized and processed Hi-C data by binning reads per 100 kb chromosomal segments to generate contact heatmaps as described before . Visual inspection of the heatmaps suggested that chromosomes folded very similar between wild-type and TKO cells (Fig. 5a). Principal component analysis of Hi-C data was previously used to uncover an A and B compartment where active and inactive chromatin regions, respectively, preferentially cluster. When applied to our datasets it showed that chromosomal domain organization and overall 3D genome structure is indeed very similar between the two cell types (Fig. 5a). In fact, the contact profiles that we generated for wild-type and histone H1-depleted ES cells were more similar to each other than any of the two was to a previously published wild-type ES cell contact profile  or pro-B-cell contact profile  (Figure S8 in Additional file 1). Thus, a 50 % depletion of linker histone H1 can be tolerated without profound changes in the overall 3D genome.
We previously reported that ES cells harbor a unique 3D genome, hallmarked in general by a more random higher order topology with particularly the inactive chromatin compartment being spatially unorganized. Another feature of the pluripotent 3D genome is the specific clustering of genomic regions dense in binding of pluripotency factors . This preferential clustering of dense pluripotency factor-associated regions was also appreciable in both our wild-type and H1 TKO cells (Fig. 5b) and confirmed that the TKO cells harbored an ES cell identity (Fig. 5c). Compared with their wild-type counterparts, the TKO cells showed a slightly increased capacity for chromosomal regions to contact each other over distance, but this effect clearly was not as pronounced as seen in, for example, differentiated pro-B cells (Fig. 5c).
Realizing that the overall topology of chromosomes is unaltered in H1-depleted cells but that domains gain some capacity to reach out and contact other domains elsewhere on their chromosome, we used the recently published TAD cross-boundary ratio to search for more subtle topological changes. The TAD cross-boundary ratio divides the intra-domain contacts over the inter-domain contacts  and as such can serve as a domain score. It is the most gene-dense TADs that are most enriched for DHSs and for active chromatin marks which show the highest such domain score (Fig. 5d). Active TADs therefore seem the structurally most isolated chromosomal entities. We calculated this value for each TAD in wild-type and H1 TKO cells and compared them. In agreement with the observation that in histone H1 depleted cells chromosomal sites more easily engage in contacts over very large distances, TADs generally had a lower domain score in TKO versus wild-type cells (Fig. 5e). Histone H1 therefore seems to contribute to the topological segmentation of chromosomes.
We then wished to identify the TADs that were most sensitive to topological changes upon histone H1 depletion. For this we computed the difference in domain score between wild-type and TKO cells and ranked TADs accordingly. Despite almost all TADs showing a reduction in the domain score and despite this difference being relatively modest, we nevertheless found that the degree of structural changes significantly correlated (chi-squared test p < < 10−6) with the amount of epigenetic changes observed per TAD (Fig. 5f). This was true no matter which of the investigated marks was considered. Thus, the top 20 % topologically most altered TADs were those that also carried most sites with altered hypersensitivity, most sites with modified levels of H3K4me1 and/or H3K4me3 and most differentially methylated CpGs. In contrast, the TADs most resistant to topological changes were those showing the least epigenetic changes upon H1 depletion. Interestingly, such correlations were not found with gene density, TAD size or differential gene expression (Figure S9 in Additional file 1). It is, therefore, not necessarily the TADs with highest gene content nor the larger TADs and also not the TADs with most striking changes in transcriptional output that are most sensitive to topological changes. Rather, alterations in the epigenetic landscape appear to best correlate with topological alterations of TADs. Thus, while gene expression clearly correlates with the nuclear positioning of TADs relative to each other and to, for example, the nuclear periphery, our current data provide further evidence that gene expression and higher-order chromosome topology are not causally related [36–38]. Rather, they may be independently controlled by the locally associated trans-acting factors.
Our data show that cells require normal histone H1 levels to expose their proper regulatory landscape. Reducing the levels of histone H1 results in massive epigenetic changes and altered topological organization, particularly at the most active chromosomal domains. Changes in TAD configuration coincide with epigenetic landscape changes but not with transcriptional output changes, supporting the emerging concept that transcriptional control and nuclear positioning of TADs are not causally related but independently controlled by the locally associated trans-acting factors.
Materials and methods
Wild-type and H1 TKO ES cells  were grown on irradiated mouse embryonic fibroblasts in Dulbecco's modified Eagle's medium (high glucose, Gibco) with 15 % fetal bovine serum, 1× non-essential amino acids (NEAA; Gibco), 1× penicillin–streptomycin (Gibco), 1:1000 b-mercaptoethanol (Invitrogen), 1× L-glutamine (Gibco) and 1000 U/ml leukaemia inhibitory factor (Gibco).
Cells were trypsinized and plated on uncoated plates for 30 minutes at 37 °C to get rid of the feeder cells. Then a 3C template was generated as previous described . In brief, 10 million cells were cross-linked by 2 % formaldehyde, then digested with DpnII or NlaIII, and ligated to form 3C circles. Purified 3C products were then further sheared to 600–800 bp. Sheared DNA (1 μg)was used to generate high-throughput sequencing-ready sample by using TruSeq DNA sample prep kit (Illumina), following the standard commercial protocol. The Hi-C library was sequenced with Illumina paired-end sequencing.
Chromatin immunoprecipitation (ChIP) was performed as described previously [12, 13] with ChIP grade H3K4me1- and H3K4me3-specific antibodies purchased from Abcam with a few modifications. Fixed cells were lysed in a buffer containing 10 mM Na-butyrate and isolated chromatin was sonicated to 500–800 bp with a Covaris S2 sonicator at 4 °C. ChIP-seq library preparation and sequencing were performed by the Epigenomics Core Facility at the Albert Einstein College of Medicine using on an Illumina 2500 HiSeq instrument. We generated duplicate ChIP-seq libraries (both duplicate input and immunoprecipitation samples, with antibodies against H3K4me1 and H3K4me3) for both conditions (wild type and TKO). For the H3K9me3 and H3K27me3 ChIP-seq experiments, 40 million cells of both conditions (wild type and TKO) were cross-linked, washed and lysed as described previously . The obtained nuclei were dissolved in 80 μl sonication buffer, transferred to microtubes and sonicated for 12 cycles of 60 seconds using microtubes in the Covaris S series with the following settings: intensity 3, duty cycle 20 %, 200 cycles/bursts. Supernatant was cleared  and added to DynaI protein G beads that were pre-incubated with ab6002 for H3K27me3 and ab8898 for H3K9me3 from Abcam. After immunoprecipitation, beads were washed and DNA was eluted, reverse cross-linked and further purified as described before . ChIP-seq libraries were made according to the Illumina Truseq DNA library protocol, and sequencing was performed at the Utrecht Sequencing Facility on a NextSeq500. Reads from all different libraries where aligned to the reference genome (NCBI37/mm9) with bowtie2  with default settings and the --qc-filter switch. Duplicates were marked using Picard (http://broadinstitute.github.io/picard/) and were removed from the data for subsequent analyses. Regions significantly enriched for H3K4me1, H3K4me3, H3K27me3, and H3K9me3 compared with matched input samples were identified using the MACS2 peak caller  with default settings. For the H3K9me3 and H3K27me3 histone marks, the parameter - - broad was set. Analysis of differential ChIP enrichment was done using diffReps  with parameters -me gt --pval 0.001 --frag 150.
Total RNAs were prepared from ES cells adapted to gelatinized dishes using RiboPure RNA purification kits (Ambion). Paired-end library construction was performed using Tru-seq kits (Illumina). Resulting libraries were run on the Hi-seq 2000 (Salk Institute), generating 2× 100-bp paired-end reads. We aligned reads of two replicate wild-type ES cell RNA-seq libraries and three replicate H1 TKO ES cell RNA-seq libraries to the reference genome (NCBI37/mm9) with TopHat  and used Cufflinks and CuffDiff  for differential expression analysis of RNA-seq expression for a non-redundant collection of 20,876 known RefSeq transcripts. We considered genes with a marginal p value smaller than 0.05 and an absolute log2 fold-change bigger than 1 to be differentially expressed (598 genes).
Genome-wide DNA methylation analysis using the HELP-tagging assay
Genomic DNA was isolated and digested with HpaII and MspI, and Illumina library preparation was performed exactly as described previously . Library sequencing was performed in the Epigenomics Core Facility at the Albert Einstein College of Medicine. We computed the HELP angle as described in  and used it as a measure for the percentage of methylated cytosines. We performed binomial tests for differential methylation and this resulted in 15,492 differentially methylated sites with a p value smaller than 1e-6.
DNase I hypersensitivity assay
DNase I hypersensitivity assay was essentially carried out as described in . In brief, nuclei were extracted in lysis buffer (15 mM Tris‐HCl, 15 mM NaCl, 60 mM KCl, 1 mM EDTA, 0.5 mM EGTA, 0.5 mM spermidine) by incubating for 10 minutes on ice. Then, nuclei were incubated for 3 minutes at 37 °C in the same lysis buffer with 1 mM CaCl2 and with limiting concentrations of the DNA endonuclease deoxyribonuclease I (DNase I). Reactions were stopped by adding stop buffer (50 mM Tris‐HCl, 100 mM NaCl, 0.1 % SDS, 100 mM EDTA, 1 mM spermidine, 0.5 spermine, pH 8.0) and purified fragments were recovered by sucrose ultracentrifugation, end‐repaired and ligated with adapters, followed by sequencing on the Illumina sequencing platform. From an initial collection of 33 different DNase-seq libraries (17 wild type, 16 TKO), we filtered out three high quality replicates in each condition according to their SPOT score. Reads were aligned to the reference genome (NCBI37/mm9) and we considered 36-bp reads that aligned uniquely and contained no more than two mismatches as properly mapped reads. We used the Hotspot  algorithm to identify DHSs in all six samples separately. For differential DNase-seq analysis between wild type and TKO, we used the PoissonSeq R package . We compared DNase-seq coverage in 89,875 different regions (with sufficient coverage in either condition) and this resulted in a set of 4166 regions with statistically significant difference in coverage after applying a multiple testing procedure (false discovery rate 5 %).
All statistical analyses were performed under R/Bioconductor  using custom R scripts. Manipulation with and computation of statistics on genomic intervals and domains was done using the GenomicRanges package . Analysis of ChIP-seq data and DNase I hypersensitivity data and generation of heatmaps were done using the compEpiTools package (http://genomics.iit.it/groups/computational-epigenomics.html).
The authors state that no ethical approval was required for this study.
Availability of supporting data
All the raw and processed data for the experiments described in this paper have been submitted to NCBI Gene Expression Omnibus (GEO) under accession number GSE75426.
DNase I hypersensitivity site
Differentially methylated region
Topologically associating domain
Van Holde KE. Chromatin. New York: Springer; 1989.
Wolffe A. Chromatin: structure and function. 3rd ed. San Diego, CA; London: Academic; 1998.
Brown DT. Histone variants: are they functionally heterogeneous? Genome Biol. 2001;2:REVIEWS0006.
Happel N, Doenecke D. Histone H1 and its isoforms: contribution to chromatin structure and function. Gene. 2009;431:1–12.
Izzo A, Kamieniarz K, Schneider R. The histone H1 family: specific members, specific functions? Biol Chem. 2008;389:333–43.
Sirotkin AM, Edelmann W, Cheng G, Klein-Szanto A, Kucherlapati R, Skoultchi AI. Mice develop normally without the H1(0) linker histone. Proc Natl Acad Sci U S A. 1995;92:6434–8.
Fan Y, Sirotkin A, Russell RG, Ayala J, Skoultchi AI. Individual somatic H1 subtypes are dispensable for mouse development even in mice lacking the H1(0) replacement subtype. Mol Cell Biol. 2001;21:7933–43.
Lin Q, Sirotkin A, Skoultchi AI. Normal spermatogenesis in mice lacking the testis-specific linker histone H1t. Mol Cell Biol. 2000;20:2122–8.
Drabent B, Saftig P, Bode C, Doenecke D. Spermatogenesis proceeds normally in mice without linker histone H1t. Histochem Cell Biol. 2000;113:433–42.
Rabini S, Franke K, Saftig P, Bode C, Doenecke D, Drabent B. Spermatogenesis in mice is not affected by histone H1.1 deficiency. Exp Cell Res. 2000;255:114–24.
Fantz DA, Hatfield WR, Horvath G, Kistler MK, Kistler WS. Mice with a targeted disruption of the H1t gene are fertile and undergo normal changes in structural chromosomal proteins during spermiogenesis. Biol Reprod. 2001;64:425–31.
Fan Y, Nikitina T, Zhao J, Fleury TJ, Bhattacharyya R, Bouhassira EE, et al. Histone H1 depletion in mammals alters global chromatin structure but causes specific changes in gene regulation. Cell. 2005;123:1199–212.
Yang SM, Kim BJ, Norwood Toro L, Skoultchi AI. H1 linker histone promotes epigenetic silencing by regulating both DNA methylation and histone H3 methylation. Proc Natl Acad Sci U S A. 2013;110:1708–13.
Maclean JA, Bettegowda A, Kim BJ, Lou CH, Yang SM, Bhardwaj A, et al. The rhox homeobox gene cluster is imprinted and selectively targeted for regulation by histone h1 and DNA methylation. Mol Cell Biol. 2011;31:1275–87.
Dixon JR, Selvaraj S, Yue F, Kim A, Li Y, Shen Y, et al. Topological domains in mammalian genomes identified by analysis of chromatin interactions. Nature. 2012;485:376–80.
Splinter E, Heath H, Kooren J, Palstra RJ, Klous P, Grosveld F, et al. CTCF mediates long-range chromatin looping and local histone modification in the beta-globin locus. Genes Dev. 2006;20:2349–54.
Meshorer E, Yellajoshula D, George E, Scambler PJ, Brown DT, Misteli T. Hyperdynamic plasticity of chromatin proteins in pluripotent embryonic stem cells. Dev Cell. 2006;10:105–16.
Melcer S, Hezroni H, Rand E, Nissim-Rafinia M, Skoultchi A, Stewart CL, et al. Histone modifications and lamin A regulate chromatin protein dynamics in early embryonic stem cell differentiation. Nat Commun. 2012;3:910.
Woodcock CL, Skoultchi AI, Fan Y. Role of linker histone in chromatin structure and function: H1 stoichiometry and nucleosome repeat length. Chromosome Res. 2006;14:17–25.
Fan Y, Nikitina T, Morin-Kensicki EM, Zhao J, Magnuson TR, Woodcock CL, et al. H1 linker histones are essential for mouse development and affect nucleosome spacing in vivo. Mol Cell Biol. 2003;23:4559–72.
de Wit E, Bouwman BA, Zhu Y, Klous P, Splinter E, Verstegen MJ, et al. The pluripotent genome in three dimensions is shaped around pluripotency factors. Nature. 2013;501:227–31.
Lieberman-Aiden E, van Berkum NL, Williams L, Imakaev M, Ragoczy T, Telling A, et al. Comprehensive mapping of long-range interactions reveals folding principles of the human genome. Science. 2009;326:289–93.
Suzuki M, Jing Q, Lia D, Pascual M, McLellan A, Greally JM. Optimized design and data analysis of tag-based cytosine methylation assays. Genome Biol. 2010;11:R36.
Nora EP, Lajoie BR, Schulz EG, Giorgetti L, Okamoto I, Servant N, et al. Spatial partitioning of the regulatory landscape of the X-inactivation centre. Nature. 2012;485:381–5.
Ernst J, Kellis M. ChromHMM: automating chromatin-state discovery and characterization. Nat Methods. 2012;9:215–6.
Cao K, Lailler N, Zhang Y, Kumar A, Uppal K, Liu Z, et al. High-resolution mapping of h1 linker histone variants in embryonic stem cells. PLoS Genet. 2013;9:e1003417.
Lu X, Wontakal SN, Kavi H, Kim BJ, Guzzardo PM, Emelyanov AV, et al. Drosophila H1 regulates the genetic activity of heterochromatin by recruitment of Su(var)3-9. Science. 2013;340:78–81.
Heinz S, Benner C, Spann N, Bertolino E, Lin YC, Laslo P, et al. Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. Mol Cell. 2010;38:576–89.
Zhang Y, Cooke M, Panjwani S, Cao K, Krauth B, Ho PY, et al. Histone h1 depletion impairs embryonic stem cell differentiation. PLoS Genet. 2012;8:e1002691.
Schuettengruber B, Oded Elkayam N, Sexton T, Entrevan M, Stern S, Thomas A, et al. Cooperativity, specificity, and evolutionary stability of Polycomb targeting in Drosophila. Cell Rep. 2014;9:219–33.
Zhang Y, Liu Z, Medrzycki M, Cao K, Fan Y. Reduction of Hox gene expression by histone H1 depletion. PLoS One. 2012;7:e38829.
Sexton T, Yaffe E, Kenigsberg E, Bantignies F, Leblanc B, Hoichman M, et al. Three-dimensional folding and functional organization principles of the Drosophila genome. Cell. 2012;148:458–72.
Rao SS, Huntley MH, Durand NC, Stamenova EK, Bochkov ID, Robinson JT, et al. A 3D map of the human genome at kilobase resolution reveals principles of chromatin looping. Cell. 2014;159:1665–80.
Zhang Y, McCord RP, Ho YJ, Lajoie BR, Hildebrand DG, Simon AC, et al. Spatial organization of the mouse genome and its role in recurrent chromosomal translocations. Cell. 2012;148:908–21.
Chandra T, Ewels PA, Schoenfelder S, Furlan-Magaril M, Wingett SW, Kirschner K, et al. Global reorganization of the nuclear landscape in senescent cells. Cell Rep. 2015;10:471–83.
Noordermeer D, Branco MR, Splinter E, Klous P, van Ijcken W, Swagemakers S, et al. Transcription and chromatin organization of a housekeeping gene cluster containing an integrated beta-globin locus control region. PLoS Genet. 2008;4:e1000016.
Palstra RJ, Simonis M, Klous P, Brasset E, Eijkelkamp B, de Laat W. Maintenance of long-range DNA interactions after inhibition of ongoing RNA polymerase II transcription. PLoS One. 2008;3:e1661.
Therizols P, Illingworth RS, Courilleau C, Boyle S, Wood AJ, Bickmore WA. Chromatin decondensation is sufficient to alter nuclear organization in embryonic stem cells. Science. 2014;346:1238–42.
Splinter E, de Wit E, van de Werken HJ, Klous P, de Laat W. Determining long-range chromatin interactions for selected genomic sites using 4C-seq technology: from fixation to computation. Methods. 2012;58:221–30.
Vermunt MW, Reinink P, Korving J, de Bruijn E, Creyghton PM, Basak O, et al. Large-scale identification of coregulated enhancer networks in the adult human brain. Cell Rep. 2014;9:767–79.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–9.
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:R137.
Shen L, Shao NY, Liu X, Maze I, Feng J, Nestler EJ. diffReps: detecting differential chromatin modification sites from ChIP-seq data with biological replicates. PLoS One. 2013;8:e65598.
Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25:1105–11.
Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012;7:562–78.
Vierstra J, Rynes E, Sandstrom R, Zhang M, Canfield T, Hansen RS, et al. Mouse regulatory DNA landscapes reveal global principles of cis-regulatory evolution. Science. 2014;346:1007–12.
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:75–82.
Li J, Witten DM, Johnstone IM, Tibshirani R. Normalization, testing, and false discovery rate estimation for RNA-sequencing data. Biostatistics. 2012;13:523–38.
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5:R80.
Lawrence M, Huber W, Pages H, Aboyoun P, Carlson M, Gentleman R, et al. Software for computing and annotating genomic ranges. PLoS Comput Biol. 2013;9:e1003118.
We thank Britta Bouwman for help with the figures, Michal Mokry and Edwin Cuppen for help with high-throughput sequencing, Stieneke van den Brink for ES cell culture, and Masako Suzuki, Varun Gupta, and John Greally for guidance in preparing HELP-tagging libraries and data analysis.
This work was supported by National Institutes of Health Grant CA079057 and GM116143 to A.I.S. and an NWO/CW TOP grant (714.012.002), an NWO VICI grant 724.012.003, an EU grant 2010-259743 (MODHEP) and a European Research Council Starting Grant (209700, ‘4C’) to W.d.L.
The authors declare that they have no competing interests.
GG and BAB computationally analyzed the data; YZ, BJK, SMY, MJAMV, HP, MWV, MPC and PJW performed the experiments; TSM generated and analyzed RNA-seq data; JAS generated and analyzed Dnase I data; GG, AIS, and WdL conceived the study, participated in its design and drafted the manuscript. All authors read and approved the final manuscript.
Geert Geeven, Yun Zhu, Byung Ju Kim and Boris A. Bartholdy contributed equally to this work.
DNA methylation changes per chromosome in histone H1-depleted ES cells. Figure S2. Hypomethylation in TKO cells occurs preferentially in gene-dense TADs regardless of GC content. Figure S3. ChIP-seq of four different histone marks in wild-type and TKO mouse ES cells. Figure S4. Large-scale changes in H3K4me1 and H3K4me3 sites upon depletion of histone H1. Figure S5. Analysis of DNA binding motifs at de novo formed DHSs in TKO cells. Figure S6. Over 10 % of de novo formed DHSs also show loss in CpG methylation in H1 TKO cells. Figure S7. Active chromatin marks accumulate in the most gene-dense TADs. Figure S8. Higher-order genome topology is very similar between wild-type and histone H1 TKO cells. Figure S9. Changes in compartment organization are not related to gene content, TAD size or differential expression. (PDF 4637 kb)
The number of loci at which changes are observed for various combinations of epigenetic marks. (XLS 13 kb)
About this article
Cite this article
Geeven, G., Zhu, Y., Kim, B.J. et al. Local compartment changes and regulatory landscape alterations in histone H1-depleted cells. Genome Biol 16, 289 (2015) doi:10.1186/s13059-015-0857-0
- Histone H1
- Chromosome conformation