# Linker histones are fine-scale chromatin architects modulating developmental decisions in Arabidopsis

Kinga Rutowicz, Maciej Lirski, Benoît Mermaz, Gianluca Teano, Jasmin Schubert, Imen Mestiri, Magdalena Kroteń, Tohnyui Ndinyanka Fabrice, Simon Fritz, Stefan Grob, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Kinga Rutowicz, Maciej Lirski, Benoît Mermaz, Gianluca Teano, Jasmin Schubert, et al. Linker histones are fine-scale chromatin architects modulating developmental decisions in Arabidopsis. Genome Biology, 2019, 20, pp.157. 10.1186/s13059-019-1767-3. hal-02378703

HAL Id: hal-02378703

https://hal.science/hal-02378703

Submitted on 3 Jun 2021

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>
RESEARCH

Open Access

<figure><img src="image_3.png" /></figure>

Check for
updates

# Linker histones are fine-scale chromatin architects modulating developmental decisions in *Arabidopsis*

Kinga Rutowicz<fnref n="1" />†, Maciej Lirski<fnref n="2" />†, Benoît Mermaz<fnref n="1" />,³, Gianluca Teano⁴, Jasmin Schubert<fnref n="1" />, Imen Mestiri⁴, Magdalena A. Kroten⁵, Tohnyui Ndinyanka Fabrice<fnref n="1" />, Simon Fritz<fnref n="1" />, Stefan Grob<fnref n="1" />, Christoph Ringli<fnref n="1" />, Lusik Cherkezyan⁶, Fredy Barneche⁴, Andrzej Jerzmanowski<fnref n="2" />,⁷* and Célia Baroux<fnref n="1" />*

## Abstract

**Background:** Chromatin provides a tunable platform for gene expression control. Besides the well-studied core nucleosome, H1 linker histones are abundant chromatin components with intrinsic potential to influence chromatin function. Well studied in animals, little is known about the evolution of H1 function in other eukaryotic lineages for instance plants. Notably, in the model plant *Arabidopsis*, while H1 is known to influence heterochromatin and DNA methylation, its contribution to transcription, molecular, and cytological chromatin organization remains elusive.

**Results:** We provide a multi-scale functional study of *Arabidopsis* linker histones. We show that H1-deficient plants are viable yet show phenotypes in seed dormancy, flowering time, lateral root, and stomata formation—complemented by either or both of the major variants. H1 depletion also impairs pluripotent callus formation. Fine-scale chromatin analyses combined with transcriptome and nucleosome profiling reveal distinct roles of H1 on hetero- and euchromatin: H1 is necessary to form heterochromatic domains yet dispensable for silencing of most transposable elements; H1 depletion affects nucleosome density distribution and mobility in euchromatin, spatial arrangement of nanodomains, histone acetylation, and methylation. These drastic changes affect moderately the transcription but reveal a subset of H1-sensitive genes.

**Conclusions:** H1 variants have a profound impact on the molecular and spatial (nuclear) chromatin organization in *Arabidopsis* with distinct roles in euchromatin and heterochromatin and a dual causality on gene expression. Phenotypical analyses further suggest the novel possibility that H1-mediated chromatin organization may contribute to the epigenetic control of developmental and cellular transitions.

**Keywords:** Linker histones, H1, Chromatin, Heterochromatin, Histone methylation, Developmental transitions, Nucleosome occupancy

## Background

Linker histones (H1) belong to the major constituents of plant and animal chromatin besides the core nucleosomal histones. H1 variants (collectively referred to as H1 thereafter) appeared early during evolution as lysine-rich, proto-linker histones found in the ancestral

eukaryotes such as protists but not Archaea [1]. In contrast to core nucleosomal constituents, H1 form a highly divergent class of histones [1]. H1 typically possesses a tripartite structure conserved across all eukaryotes: it consists of a short and flexible N-terminal tail, a structured globular domain (GH1) which interacts with a nucleosome dyad, and a structurally disordered, lysine-rich (highly basic) C-terminal tail. The C-terminal tail of H1 variants varies in length and composition among isotypes and organisms, conferring their various chromatin compaction potential by interacting with internucleosomal DNA and drawing adjacent nucleosomes together

<footnote>* Correspondence: andyj@ibb.waw.pl; cbaroux@botinst.uzh.ch</footnote>

<footnote>†Kinga Rutowicz and Maciej Lirski contributed equally to this work.</footnote>

<footnote n="2">²Institute of Biochemistry and Biophysics, Polish Academy of Sciences, Pawinskiego 5a, 02-106 Warsaw, Poland</footnote>

<footnote n="1">¹Institute of Plant and Microbial Biology, Zürich-Basel Plant Science Center, University of Zürich, Zürich, Switzerland</footnote>

<footnote>Full list of author information is available at the end of the article</footnote>
[2, 3]. In animals, several H1 isotypes can coexist in the same cell type playing both redundant and specific roles in chromatin organization [4]. H1 proteins constitute a highly mobile fraction of the animal chromatin; their apparent constitutive presence results from a steady-state level of dynamic binding [5]. H1-mediated higher-order chromatin organization influences RNA polymerase II accessibility, hence gene expression, DNA replication, chromosome segregation, and DNA repair [4, 6–9]. In *Drosophila*, *Neurospora crassa*, and mammalian cells, H1 deposition is part of a crosstalk with the epigenetic landscape, notably DNA methylation [10–12] and histone H3 methylation [12–14]. Given the profound impact of H1 proteins on chromatin organization, it is surprising that H1 depletion is tolerated in some organisms: this is the case of *Tetrahymena*, yeast, fungi, worms [15–19]. However, in mouse and *Drosophila*, it impairs viability [10, 20]. The intrinsic role of H1 in chromatin organization and yet the variable impact of its depletion create an apparent paradox which is difficult to address given the high functional redundancy between variants, particularly in mammalian genomes.

In the plant kingdom, H1 gene variants can be traced back to the earliest land plants [21]. The flowering plant *Arabidopsis thaliana* possesses only three canonical H1 variants, and the reduced number of isotypes makes it suitable for analyzing H1 functions [21, 22]. Two of them, H1.1 and H1.2, are ubiquitous in differentiated shoot and root tissues and are expressed throughout vegetative development. Whether they play distinct or similar roles on chromatin organization is not known. Despite their apparent ubiquitous expression, both variants are evicted at the somatic-to-reproductive cell fate transition and undetectable in the reproductive lineage [23–25]. Unlike in animals where germline-specific H1 variants have been described [26], no gametophyte-specific variants have yet been characterized in plants. H1.1 protein abundance is nevertheless restored in male gametes (Baroux, unpublished and [27]) where it could thus play a sperm cell-specific function, yet to be determined. H1.3, by contrast, is not a routine variant but is specifically incorporated in response to abiotic stresses: this stress-inducible variant is thought to act as a pioneer factor priming transcriptional reprogramming that leads to physiological adaptations [28]. The genomic distribution of H1 variants is widely present along the *Arabidopsis* genome, spanning both heterochromatin and euchromatin chromosomal domains [28, 29]. Specifically, H1.1 and H1.2 are enriched at the 3' and 5' ends of transposable elements (TEs) and over gene bodies where its profile is anti-correlated to the transcription and H3K4me3 levels [28]. They are subjected to numerous post-translational modifications with undescribed functions [30]. H1 variants further influence DNA

methylation patterns in *Arabidopsis*, primarily but not exclusively, in heterochromatin and in all sequence contexts (i.e., CG, CHG, CHH) [22, 28, 31, 32]. H1-containing nucleosomes are thought to create a barrier to DNA methyltransferases, a structural conflict that can be resolved by the activity of SWI/SNF types of chromatin remodelers [32, 33].

Overall, however, the functional influence of H1 on chromatin organization and function in plants remains elusive, and our understanding of H1 function on cell specification and development is limited. Whereas H1 downregulation by RNAi knockdown in *Arabidopsis* produced pleiotropic phenotypes [22], where possible off-target effects on H1-related proteins cannot be excluded [21], double and triple knockout lines are viable and provide a unique opportunity to investigate the role of H1 during development [28, 31].

Here, we report that H1 has a profound role not only in heterochromatin but also in euchromatin in *Arabidopsis*, which was so far overlooked. It regulates chromatin compaction and organization at the nuclear level, being required for heterochromatin formation and influencing on the spatial distribution of nanoscale, compact domains in euchromatin. Mutant analyses reveal that H1 restricts nucleosome mobility and influences H3K9 acetylation, H3K27, and H3K4 methylation levels. We also show that H1 regulates nucleosome distribution and provides distinct nucleosome density over expressed genes correlating with their expression levels. Interestingly, these nucleosomal states and large-scale chromatin structures mediated by H1 are not epistatic to transcriptional regulation for a vast majority of loci. A small fraction of genes and TEs is, however, H1-sensitive suggesting that chromatin level regulation primes over other transcriptional controls for these loci. Finally, the analysis of mutant and complemented mutant lines unveils a new role of H1-mediated organization during development.

## Results

### H1-deficient *Arabidopsis* lines show deregulation of several developmental transitions

T-DNA insertion mutant alleles for the three *H1* genes were previously introgressed producing double and triple mutants which are viable [25, 28, 31, 32]. This is in contrast to RNAi lines inducing pleiotropic aberrations [22], where possible off-target effects on H1-related genes [21] cannot be excluded. *H1* triple (T-DNA) mutant plants, thereafter called *3h1*, were previously characterized at the molecular level and do not show detectable levels of *H1* transcripts nor H1 proteins [25]. *3h1* appeared at first normal with regular shoot and root apparatus indicating that organogenesis coordinated by the apical meristems are not affected. We investigated the development of *3h1*-mutant plants in further details and
noticed subtle phenotypes at developmental transitions. We noticed for example that a large fraction of *3h1* seeds underwent a prolonged dormancy after harvest (Fig. 1a). In addition, *3h1* plants showed precocious in-duction of flowering (Fig. 1b, Additional file 1: Figure S1), a phenotype which was rescued by the introduction of tagged H1 variants in *3h1* plants (Additional file 1: Figure S1C). We then rationalized that H1 might be im-plicated in other developmental transitions. To test this hypothesis, we analyzed the formation of lateral roots, root hairs, and epidermal stomata as traits reflecting the cell identity establishment concurrent to organogenesis. In *Arabidopsis*, the specification of lateral root primordia from pericycle founder cells and the differentiation of root hairs from epidermal cells follow a regular pattern modulated by developmental and environmental cues and subjected to epigenetic regulation [34, 35]. Firstly, compared to wild type, *3h1* seedlings produced a higher number of lateral roots per root length unit (Fig. 1c) and a higher density of root hairs (Fig. 1d). Again, both de-fects were rescued upon complementation with H1.1 and H1.2 transgenes (Fig. 1c, d; Additional file 1: Figure S2A). These observations suggest that H1 is required for proper spatial and temporal control of lateral root and root hair lineage initiation. In addition, the unicellular fate of root hairs was occasionally compromised in *3h1* but not in wild-type plants, with the appearance of mul-tiple nuclei and cell boundaries (Fig. 1e, Additional file 1: Figure S2B). Interestingly, this trait was previously corre-lated with an unstable epigenetic root hair cell fate in mutants impaired in *Polycomb*-group Repressive Com-plex2 (PRC2)-based chromatin regulation [36].

Secondly, we noticed that stomata patterning in the epidermis of cotyledons was altered in *3h1* with a higher occurrence of high-degree (tertiary and quaternary) clus-ters, associated with complex arrangements, collated sto-mata, or atypical division patterns in early stages. These features were not found in the wild type and were re-stored upon the introduction of H1.1 and H1.2 trans-genes (Fig. 1f, Additional file 1: Figure S2C). These observations suggest a loose control by H1 of stomatal spacing presumably involving occasional re-initiation events [37]. Finally, we also tested how *3h1* tissues re-spond to reprogramming in in vitro culture. We mea-sured a decreased efficiency in callus development compared to the wild type (Fig. 1g), a feature mostly at-tributed to H1.3 in our complementation experiments (Additional file 1: Figure S3). It was shown before that for efficient leaf-to-callus transition, genome-wide repro-gramming of H3K27me3 is a critical step, and in PRC2-mutants callus, development is defective [38].

In summary, we identified defects in seed dormancy control, flowering time control, lateral root initiation, and the fate of root hair as well as guard cells which

were previously overlooked in H1-depleted plants. This collectively indicates that linker histones are re-quired for a tight control of developmental and cellu-lar transitions.

### H1 variants are necessary to form compact heterochromatin domains but are dispensable for TE silencing and peripheral positioning of chromocenters

H1.1 and H1.2 variants are largely ubiquitously expressed in plant tissues except in the reproductive lineage [24, 25], while H1.3 is usually not expressed in most cell types [28]. Transient depletion of H1 proteins at the somatic-to-reproductive transition precedes, and probably causes, drastic chromatin changes at the nu-clear level [24, 25]. In order to analyze the consequence of H1 depletion on chromatin organization at the cyto-logical scale, we quantified the following parameters in isolated leaf nuclei from *3h1* and compared them to that in wild-type nuclei: number of chromocenters, the rela-tive heterochromatin fraction (RHF) measuring the frac-tion of chromatin condensed in conspicuous foci or domains [39], and nuclear size informing on the global compaction level of chromatin. These analyses showed that *3h1* nuclei have a stark reduction in heterochroma-tin content with a RHF of less than 5%, were signifi-cantly larger than the wild type and failed to form the typical 6–8 heterochromatic chromocenters (CCs) nor-mally seen in most wild-type somatic nuclei (Fig. 2a, b; Additional file 2: Table S1).

*Arabidopsis* CCs are largely composed of centromeric and pericentromeric transposable element (TE) repeats of the five chromosome pairs, and a subset of two to four CCs comprising ribosomal (rDNA) repeats is asso-ciated with the nucleolus [39, 41]. To understand where these repeats were located in the *3h1* nuclei in the ab-sence of conspicuous chromocenters, we carried out 3D fluorescence in situ hybridization (FISH) on intact nuclei embedded in acrylamide gel pads [42]. Centromeric re-peats were remarkably dispersed in *3h1*-mutant nuclei, while rDNA repeats were localized in the remaining compact CCs, as in the wild-type nuclei (Fig. 2c). H1 is thus essential for maintaining the structural, compact domains at the (peri-) centromere regions but is dis-pensable for the heterochromatinization of rDNA repeat loci. Interestingly, the decondensed centromeric repeats in the *3h1* nuclei remain located at the periphery as de-scribed in the wild-type nuclei [43] (Fig. 2d). Thus, H1-mediated CC compaction occurs downstream of the spatial positioning of the centromeric regions. In addition, high-resolution imaging indicated the presence of nanoscopic bodies of condensed chromatin in the *3h1* nuclei (Fig. 2e). These were not detected in the wild-type nuclei and therefore could correspond to dispersed
<figure><img src="image_4.png" /><figcaption>**B** Flowering time</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**C** Lateral Roots</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**D** Root hairs</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**E** Root hair fate</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**F** Stomatal lineage</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**G** Callus production</figcaption></figure>

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /><figcaption>**Fig. 1** (See legend on next page.)</figcaption></figure>
(See figure on previous page.)

**Fig. 1** H1 depletion relaxes the epigenetic control of several developmental and cellular transitions. *3h1* mutant plants show a relaxed control of seed dormancy (**a**), flowering time (**b**), lateral root formation (**c**), root hair density (**d**), root hair fate (**e**), stomatal spacing (**f**), and are impaired in callus production in vitro. *3h1* shows, compared to wild-type **a** a prolonged dormancy, i.e., lower germination rate 1 day post-harvest but alleviated 3 weeks post-harvest, **b** early flowering measured by the number of rosette leaves at bolting (*n* = number of plants analyzed in this replicate; more experiment replicates are shown in Additional file 1: Figure S1), **c** increased number of lateral roots (eight DAG seedlings), **d** increased root hair density (scale bar, 200 µm), **e** occasional multi-cellular root hairs showing visible cell walls (arrows) in *Renaissance* staining (i, iii) compared to in wild-type and additional nuclei in DAPI counterstaining (ii: arrows, red; additional examples and comments are presented in Additional file 1: Figure S2; scale bar, 10 µm), **f** stomatal complexes with reduced spacing (adaxial cotyledon epidermis; see Additional file 1: Figure S2 for additional examples and quantifications; scale bar, 20 µm), and **g** decreased callus size produced from excised cotyledons in vitro (scale bar, 1 cm). Wild-type segregants (wt) were compared with triple mutant tissues/seedlings (*3h1*) and, whenever indicated, with complemented lines expressing H1.1 and H1.2 variants only (*3h1*; H1) or all three H1 variants (*3h1*; H1*). Statistical tests (**a**, **b** Welch *t* test; **c** Fisher exact test) were performed against wt replicates, \*\*\* *p* < 0.001; ns, not significant

heterochromatin regions that were not assembled into chromocenter structures.

The three *Arabidopsis* H1 variants vary in protein structure and expression pattern and may not equally contribute to heterochromatin organization [28]. We quantified heterochromatin content in whole-mount roots and cotyledons in wild-type, mutant, and complemented mutant seedlings. We found an intermediate reduction in heterochromatin content in the double *h1.1h1.2* mutant compared to that in the triple mutant in roots, an effect restored by tagged H1.1 or H1.2 variants (Additional file 1: Figure S4A, B). We explain this by a possible compensatory expression of *H1.3* when H1.1 and H1.2 are depleted, as inferred from reporter analysis in *3h1* roots showing ectopic H1.3-GFP levels (Additional file 1: Figure S4D). However, in cotyledons, which are leaf structures of embryonic origin, CC formation and compaction of centromeric and pericentromeric repeats seem controlled by the H1.1 and H1.2 variants only (Additional file 1: Figure S4C). In wild-type plants, the relative levels of H1.1 and H1.2 variants change along the meristematic-elongation differentiation transition in roots (Additional file 1: Figure S5) suggesting a specialization of these variants, however with a tissue-dependent, functional redundancy revealed in our mutant complementation analyses.

To investigate whether heterochromatin dispersion in *3h1* correlates with looser chromatin organization at the molecular level, we generated chromatin accessibility profiles using micrococcal nuclease sequencing (MNase-seq). We focused the analysis on nucleosomal coverage and distribution in the regions characteristic for heterochromatin and indexed as chromatin states CS8 and CS9 in the established nomenclature [40, 44]. In *Arabidopsis*, CS8 and CS9 are enriched in DNA methylation, H3.1 histone variant, and H3K27me1 and H3K9me2 modifications. In our analyses, the typical CS8 and CS9 regions have a consistent 12–15% reduction in nucleosomal density in the *3h1* nuclei (Fig. 2f, Additional file 1: Figure S6). In addition, nucleosome distribution is more variable in *3h1* heterochromatin as shown by the

broader distribution of MNase-protected regions compared to wild type: notably, we found a higher frequency of both short (< 150 nt) and unusually long (> 300 nt) fragments, while the average nucleosome repeat length (NRL) was globally shorter by 10 nt in the *3h1* mutant (Fig. 2g). We concluded from these analyses that H1 constrains nucleosomal spacing and enhances regularity in nucleosome distribution along the heterochromatin regions in *Arabidopsis*. Variability in nucleosomal spacing and distribution in heterochromatic regions might be responsible for the unstable chromocenters in the *3h1* nuclei. Interestingly, the absence of microscopically visible chromocenters does not seem to impair the deposition of their corresponding epigenetic silencing marks which remain abundant and widely redistributed (Fig. 2h-i) most likely following the dispersion of cognate heterochromatin regions.

We then asked whether this structural and spatial dispersion of heterochromatin in *3h1* is associated with reduced silencing of transposable elements (TEs) which are typically confined to chromocenters. Strikingly, RNA-seq profiling of wild-type and *3h1* seedlings did not show changes for the vast majority of TE transcript levels. This demonstrates that H1-mediated heterochromatin condensation into larger compact CCs is dispensable for silencing of the majority of heterochromatic TEs (Additional file 3: Table S2 and Additional file 4: Table S3). Nevertheless, a moderate number of 450 elements (*p* value < 0.05 and fold change > 2) representing about 1.5% of TAIR10-annotated TEs were significantly reactivated in *3h1* plants. A third of them corresponded to LTR/gypsy elements that are mainly distributed within the pericentromeric regions (Additional file 1: Figure S7, Additional file 3: Table S2, Additional file 4: Table S3).

Collectively, our observations indicate that linker histones contribute to regular nucleosomal distribution and density over heterochromatin regions and are absolutely required for the formation and/or maintenance of conspicuous chromocenter domains.
<figure><img src="image_12.png" /><figcaption>**Fig. 2** (See legend on next page.)</figcaption></figure>
(See figure on previous page.)

**Fig. 2** Loss of H1 variants leads to global chromatin decondensation but is dispensable for heterochromatin identity. Cytogenetic (a-e, h-i) and nucleosome profile (f, g) analyses of chromatin organization in 3h1 triple mutant and wild-type segregant (wt) seedlings. **a** H1 depletion induces a significant reduction of the relative heterochromatin fraction (RHF) and in the number of chromocenters (CCs) as well as an increase in nuclear size (area). ***two-sided t test, *p* < 0.001; error bars, standard error to the mean (SEM). Cytological analyses on isolated, spread leaf nuclei: **b** Typical wt and 3h1 nuclei as used in **a**, stained with DAPI. Scale bar, 2 µm. **c** H1 depletion induces a spatial dispersion of the centromeric repeats (CEN, purple) but not of the 45S rDNA, nucleolar organization region repeats (NOR, green) as shown by fluorescence in situ hybridization (FISH). **d** 3D image segmentation of the CEN signals shows that the preferentially peripheral localization of CEN repeats is unaffected in 3h1 nuclei despite their lack of condensation. **e** High-resolution imaging and deconvolution-based reconstruction of 3h1 and wt nuclei. Nanoscopic bodies of condensed chromatin are dispersed throughout the nucleus in 3h1 instead of conspicuous chromocenters as in wt. **f** Nucleosome occupancy is lower in 3h1 heterochromatin, as defined by the chromatin states (CS) 8 and 9 [40]. **g** Distribution of nucleosomal repeat lengths (NRLs) in wt and 3h1, chi-square test, ***$p$* < 0.0001. **h, i** The heterochromatic marks H3K9me1 and H3K27me1 are not reduced but redistributed in 3h1 nuclei. Scale bar, 2 µm. Isolated leaf nuclei were flow-sorted according to their 2C DNA content (a-e, h, i)

## H1 variants enable a regular spatial distribution of nanoscale chromatin domains and regulate nucleosomal density and mobility in euchromatin

As shown by genome-wide profiling, *Arabidopsis* H1 variants are abundant throughout the genome and, besides heterochromatin, are present in the euchromatin regions [28, 29]. Immunolocalization with an antibody specifically targeting plant H1 variants allowed visualizing the native distribution of H1 in situ, revealing discrete regions in euchromatin and encompassing heterochromatic CCs (Fig. 3a). Furthermore, H1.1-GFP was found to be distributed homogenously in euchromatin with, however, the occurrence of enrichment foci relative to nucleosome distribution visualized with an RFP-tagged H2B variant (Fig. 3b). We thus asked whether H1 depletion also impacts the structural organization of euchromatin regions. To resolve nanoscale level of organization, we measured chromatin density patterns on ultrathin transmission electron microscopy (TEM) preparations (Fig. 3c). For this, we used a spatial pattern analysis approach that was previously validated to capture relevant structural features of chromatin organization in cancerogenous animal cells [45] and that we adapted for images of the *Arabidopsis* nuclei [46]. In brief, a spatial autocorrelation function (ACF) of chromatin staining spatial distribution is calculated inside multiple regions of interests (ROIs, Fig. 3d, Additional file 1: Figure S8A) within the euchromatin region of each nucleus and is used to infer the distribution of structured signal intensities at given length scales (Fig. 3e). In wild-type nuclei, the spatial autocorrelation fit (ACF) is shallow, indicating high autocorrelation (i.e., highly regular arrangements) of nanodomains within length scales of ~40 nm. By contrast, the exponential decrease of the ACF fit of 3h1 euchromatin was dampen compared to the wild-type nuclei, particularly within 30–60 nm length scales (gray range, graph Fig. 3e). This means that chromatin densities also distribute over large length scales in the mutant. This characteristic is also quantified by measuring the shape of the density ACF (D value in box plot Fig. 3e). This showed a higher

dispersion of length scales in 3h1 compared to the wild-type nuclei (Fig. 3g). In other words, 3h1 euchromatin shows a significant loss of spatial homogeneity with chromatin nanodomains (high-density patches) spatially distributed along variable intervals (length scales) compared to a highly regular distribution in wild-type chromatin. This trend was reversed in mutants complemented by a tagged H1.1 variant (Additional file 1: Figure S8B) and confirmed in an independent quantification made on fluorescently immunolabeled nucleosomes captured by super resolution microscopy (Additional file 1: Figure S8C).

We then investigated the euchromatin structure at the molecular level and analyzed both nucleosome occupancy profiles and density in our MNase-seq experiments along with defined genomic elements. Interestingly, in 3h1, the shape of nucleosome occupancy profiles over the most prevalent chromatin state (CS) [40] is comparable to that of the wild-type samples (Additional file 1: Figure S6). However, nucleosomal density is clearly affected in 3h1 chromatin. Comparative density profiles for each chromatin state showed enhanced or diminished average density levels relative to the element boundaries, as for instance, the CS1- and H3K27me3-marked CS4 states, respectively (Fig. 3g, Additional file 1: Figure S6). This is interpreted as a partial loss of structural differentiation of the chromatin states in the absence of H1. As a comparison, we analyzed the impact of the loss-of-function of the CAF-1 histone chaperone contributing nucleosome assembly on nucleosome density distribution across these CS as done for [47]. The analysis showed some alteration in global levels for some but not all CS states but without modifying the profile amplitudes as H1 depletion does (Additional file 1: Figure S6).

Modular alterations in nucleosome density profiles in 3h1 may possibly arise from decreased structural constraints facilitating nucleosome redistribution compared to in wild type. To address the question of whether nucleosome turnover is increased in the absence of linker histone, we carried out fluorescence recovery after photobleaching (FRAP) analyses on the root epidermal nuclei expressing an RFP-tagged H2B reporter [48] in a wild-type or 3h1 background. The analysis showed a ~2.5 times
<figure><img src="image_13.png" /><figcaption>**Fig. 3** H1 depletion has a strong impact on euchromatin organization with increased dispersion of nanoscopic domains, altered distribution of nucleosome coverage, and increased mobility. **a, b** H1 is abundant in euchromatin distributed as discrete foci partially colocalizing with H2B. **a** H1 immunostaining and propidium iodide (PI) counterstaining as in Fig. 1 (fixed, spread nuclei; scale bar, 2 µm). **b** Representative distribution of GFP- and RFP-tagged H1.1 and H2B imaging in the nuclei from fresh root tissue (epidermis, single plane, confocal laser scanning microscopy). The inset shows a close up (1 × 1 µm) in euchromatin displaying H1-enriched regions (dominant green signals) at the nanoscale, but below optical resolution. Scale bar, 1 µm. **c–f** Ultrastructural analysis of euchromatin organization in wt vs 3h1. **c** Typical TEM image of the nuclei stained with uranylacetate on 7-nm cryosection (root epidermis, see the “Methods” section). Scale bar, 1 µm. **d** Representative region of interest (ROI) in euchromatin of wt and 3h1 nuclei used for spatial autocorrelation function (ACF) analyses. Scale bar, 500 nm. **e, f** Spatial chromatin density analyses show decreased regularity in the spatial chromatin distribution pattern in 3h1 revealed by a less shallow ACF curve within length scales of 20–60 nm (gray zone, graph, **e**) and higher dispersion of length scales as shown by a bigger range of the estimate **D** characterizing the spatial autocorrelation fit (**f**). These differences in 3h1 are restored upon complementation with an H1.1 expressing construct. ***Unpaired t test, *p* < 0.001, see also Additional file 1: Figure S8. **g** Nucleosome coverage but not qualitative distribution is altered in H1-depleted euchromatin. Antagonist effects are seen for regions of chromatin states CS1, 3, and 7 (CS1 only is shown here) and CS4 (CS according to [40]), see also Additional file 1: Figure S6 for nucleosome occupancy in 3h1 and wt over regions from all chromatin states. **h** H2B-RFP fluorescence recovery is ~2.5-fold faster in 3h1 compared to in wild-type as measured in FRAP experiments, see the “Methods” section. **i** Histone acetylation levels are lower in the 3h1 leaf nuclei compared to wt; t test, *p* < 0.001, see also Additional file 2: Table S1. Scale bar, 2 µm</figcaption></figure>

faster recovery in H1-depleted chromatin compared to in the wild-type nuclei (Fig. 3h). Specifically, while the diffu- sional component of the recovery curve (until ~10s) ap- peared similar in both wild-type and mutant chromatins, H2B-RFP recovery rapidly reaches the slow phase

characteristic of limited binding-unbinding events in wild-
type chromatin [49]. The steep and continuous recovery
in 3h1 indicates that H1 depletion unleashed constraints
in chromatin organization permitting rapid turnover of
H2B-RFP. Furthermore, this chromatin property
measured in *3h1*-mutant differentiated root cells resem-
bled that of the property of wild-type chromatin in meri-
stematic cells (pluripotent and undifferentiated) cells
(Additional file 1: Figure S9). This increased mobility did
not, however, coincide with higher levels of histone acetyl-
ation as found in meristematic cells [50]. By contrast, a
moderate decrease was measured in immunostaining
(Fig. 3i, Additional file 2: Table S1) which was further con-
firmed by immunoblotting on whole seedling chromatin
extracts (Additional file 1: Figure S16A,B).

We concluded that H1-depleted cells harbor a relaxed
and highly mobile chromatin with a low degree of struc-
tural differentiation between chromatin states. These
analyses thus uncover a role for *Arabidopsis* H1 in eu-
chromatin which has not been described previously.

H1 loss-of-function allows identifying both H1-sensitive and H1-independent genes

Next, we asked whether euchromatin relaxation induced
by H1 depletion was reflected at the molecular level.
Combined analysis of MNase-seq and RNA-seq profiles
allowed us to infer correlations between nucleosome oc-
cupancy and relative transcription level grouped into
quantiles in wild-type and *3h1* seedlings. As previously
reported [51], nucleosomal coverage inversely correlates
with gene expression levels in wild-type tissues (Fig. 4a).
Low nucleosomal density over highly expressed genes is
commonly interpreted as a result of, if not a require-
ment for, high RNA polymerase accessibility linked with
high transcriptional rates. In *3h1*, we observed a notable
reduction of the structural differentiation among tran-
scription levels with a general tendency for higher nu-
cleosomal density, specifically downstream of the
transcriptional start sites (TSS) (Fig. 4a).

Surprisingly though, this higher nucleosomal occu-
pancy over genes was not correlated with a massive de-
regulation of transcription. Indeed, comparison of *3h1*
transcriptomes determined that about 701 genes are
misregulated as compared to wild-type plants (*p* value <
0.05 and fold change > 2, Fig. 4b, Additional file 1: Fig-
ure S10, Additional file 5: Table S4). Taken together,
these observations have two important implications: (i)
H1-mediated nucleosome occupancy is a structural sig-
nature correlating with gene expression levels and (ii)
for a minor group of ~ 700 H1-sensitive genes, H1-me-
diated nucleosome occupancy is epistatic to transcrip-
tional regulation. This invites revisiting the idea that
nucleosome density influences transcription.

The class of H1-sensitive genes is characterized by 43
down- and 658 upregulated loci in *3h1* seedlings (Fig. 4b,
Additional file 5: Table S4, Additional file 1: Figure S10).
Analyzing functional categories represented among
downregulated genes clearly showed a collective role in
light-related metabolism (Additional file 4: Table S5,

Additional file 5: Table S4). The group of upregulated
genes, however, did not show any specific enrichment in
GO terms (not shown) nor a dramatic overrepresenta-
tion of specific chromatin states (Additional file 1: Figure
S11). Nevertheless, it is noteworthy that they displayed a
high periodicity in nucleosome positioning within 800
bp downstream the transcriptional start site (TSS) indi-
cating a strong phasing for this class of genes (Fig. 4c).
Highly phased nucleosomes profiles are thought to be a
feature of highly expressed genes [54], whereas, intri-
guingly, this class of H1-sensitive genes is weakly
expressed in wild-type plants (Additional file 1: Figure
S12).

Collectively, our data indicate that H1 variants provide
structural attributes defined by nucleosome density over
gene bodies that differentiate transcriptional states with
distinct expression levels. In a counterintuitive manner,
H1 depletion only affects the expression of a relatively
small set of (701) expressed genes. This observation
poses questions particularly when considering the
strong, genome-wide alterations in nucleosome occu-
pancy and euchromatin mobility in the *3h1* mutant.

H1 depletion impairs the distribution of H3K27me3 and
to a lesser extent of H3K4me3

We previously described that natural eviction of H1 in
spore mother cells (SMC) precedes a breadth of global
chromatin changes at the structural and epigenetic levels
during the somatic-to-reproductive transition [25].
These include heterochromatin decondensation, histone
hyperacetylation, elevation of H3K4me3, and decrease of
H3K27me3 levels [25]. These changes are further associ-
ated with a transient decrease of DNA methylation levels
in the CHH but not the CG sequence context [23]. To
assess whether H1 depletion is functionally linked to
such chromatin changes, we introgressed the mCG and
mCHH DynaMET reporters [23] into the triple *3h1* mu-
tant to analyze the cytological distribution and abun-
dance of DNA methylation. Albeit cytological imaging
would not capture locus-specific alterations as those de-
scribed previously in H1-depleted tissues by molecular
profiling [22, 32], the global levels and distribution pat-
terns of methylated DNA in *3h1*-mutant root nuclei
were similar to that in wild-type nuclei (Additional file 1:
Figure S13). We then performed quantitative immuno-
staining for the canonical marks H3K4me3 and
H3K27me3 previously found to be dynamically redistrib-
uted in H1-depleted SMCs [25]. We found that
H3K4me3 levels were moderately but reproducibly lower
in the *3h1* nuclei compared to in the wild-type nuclei
(Fig. 4d, Additional file 2: Table S1, Additional file 1: Fig-
ure S14). More dramatically, H3K27me3 levels were
considerably lowered in H1-depleted nuclei compared to
in wild-type nuclei (Fig. 4e, Additional file 2: Table S1)
<figure><img src="image_14.png" /><figcaption>**Fig. 4** H1 is necessary to secure transcriptional state-specific nucleosomal and epigenetic profiles yet influence only a moderate gene fraction. **a** Nucleosome distribution profiles clearly define distinct gene classes according to expression levels in wild type but no longer in *3h1*. Quintiles 5 to 1 represent categories of genes with expression levels ranked from the highest to the lowest level, respectively, as previously described [28]. **b** H1 depletion induces moderate changes in the transcriptional profile, yet a subset of 701 genes (*p* < 0.05 and fold change > 2) is misregulated. The volcano plot was cropped around the denser part of the dataset. The full plot is presented in Additional file 1: Figure S10. **c** Upregulated loci show a characteristic nucleosome occupancy with high periodicity and a higher coverage in *3h1* downstream the TSS. TSS, transcription start site. **d, e** Decreased abundance of H3K27 and H3K4 trimethylation (green, immunosignals (Ab); red, propidium iodide (PI) counterstaining) in *3h1* measured by quantitative immunostaining on isolated leaf nuclei. The graphs show the relative abundance of these respective histone modifications (Ab/PI ratio relative to wt). Scale bar, 2 µm. **f** Genes which are upregulated (*p* < 0.05) in *3h1* share a significant overlap with H3K27me3 targets defined by [52] (*p* = 0.0007, Fisher exact test) but remain distinct from those affected by the *clf-29* mutation [53]. **g** A twofold reduction of H3K27me3 levels upon H1 depletion is detected by Western blot on chromatin extracts from seedlings yet is less dramatic than in a loss of PRC2 function mutant, *clf-29*. Left panel: the blots show H3K27me3 and total histone H4 levels (loading control) detected on the same membrane. A dilution series of wild-type samples, together with the absence of signal saturation in luminescence imaging, allows showing proportionality of detection signals with protein amounts. The two original gel blots are presented in Additional file 1: Figure S15. Right panel: the histogram shows mean H3K27me3 relative to H4 signals from three individual biological replicates. The wild-type level has been arbitrarily set to 1. *3h1*, triple H1 mutant; *2h1*, *h1.1h1.2* double mutant</figcaption></figure>

while H3K27me2 was visibly unaffected (Additional file 1: Figure S14). A two-fold reduction of H3K27me3 and H3K4me3 global levels were further confirmed by immunoblotting on whole seedling chromatin extracts (Fig. 4g, Additional file 1: Figure S15 and Figure S16C, D). Of note, expression of genes encoding PRC2 subunits were not significantly altered in *3h1* plants

(Additional file 4: Table S6). Consequently, our observations indicate that linker histones are required for targeting, maintenance, or both, of H3K27me3. This raised the possibility that gene upregulation in *3h1* may indirectly result from downstream consequences of a modified epigenome landscape. Yet only 10% of the upregulated genes overlapped with known H3K27me3
genomic targets. Among those, a small share (4%) corresponds to genes misregulated in the PRC2 mutant *clf* (*curlyleaf*), compromised in histone H3-Lys27 methyl-transferase activity [53] (Fig. 4f). Thus, loss of H3K27me3 does not appear to represent a major cause of gene upregulation in *3h1* seedlings. This indicates that altered H3K4me3 and H3K27me3 landscapes may rather be a consequence of altered chromatin structure affecting targeting or spreading of histone modifications, or both.

## Discussion

H1, due to its inherent properties in promoting a higher order of chromatin folding, has long been considered a typical structural chromatin protein [55]. Accordingly in plants, H1 was shown to induce ectopic heterochromatinization when expressed in a heterologous system [56]. In addition, like in animals, plant H1 variants interplay with the DNA methylation machinery, particularly affecting TEs. H1 also affects DNA methylation at the genic regions, although only a few have been experimentally interrogated so far [22, 31–33]. Despite their expected fundamental role in chromatin accessibility, very little was known about the specific impact of H1 at the structural, epigenetic, and gene expression level in plants. In this study, we unveiled distinct effects of H1 function on nucleosome occupancy and transcriptional activity over heterochromatin and protein-coding genes, respectively (model presented in Fig. 5). We further unveiled that H1 is strictly required for heterochromatin formation but not for the establishment of heterochromatic hallmarks, that H1 provides spatial regularity in nanodomain chromatin compaction and distribution which coincides with a role in the maintenance of H3K27me3 in euchromatin, and that *in fine* H1-mediated chromatin organization secures a proper control of specific developmental transitions.

## H1 as a heterochromatin architect largely uncoupled from transcriptional silencing

We showed that H1 plays an unexpected role in heterochromatin regulation in *Arabidopsis*. This function is largely but not entirely uncoupled from transcriptional silencing, as shown by massive relaxation of peripheral heterochromatin and minimal reactivation of TE expression in plants lacking H1. H1-dependent chromocenter formation also appears to occur downstream the epigenetic marking of heterochromatin, notably DNA, H3K9, and H3K27 methylation, which remain abundant in *3h1* nuclei. Consistently, heterochromatin dispersion despite strong H3K9me2 marking was also reported in *2h1* leaf nuclei (i.e., lacking H1.1 and H1.2 variants) in a study published during our manuscript revision [58]. In the same study, the absence of H1 in the vegetative cell of

pollen grains correlates with the activation of about, and only, hundred TEs [58]. Collectively, these observations argue against the common view that heterochromatin compaction in chromocenter domains is directly responsible for the general transcriptional repression of TEs. Transposable elements upregulated in *3h1* define a functional category of loci where H1-mediated chromatin organization primes over other controls of transcriptional silencing. Those TEs are mostly pericentromeric and enriched in LINE, gypsy, and copia elements. DNA methylation at pericentromeric TEs is independent of the RNA-dependent DNA methylation (RdDM) pathway but requires the DNA methyltransferase CHROMOMETHYLTRANSFERASE2 (CMT2) and the chromatin SWI/SNF remodeler DECREASE in DNA METHYLATION1 (DDM1) [32]. At pericentromeric TE loci, H1 is thought to modulate (but not hinder) DNA methylation in the CHG context by reducing access to CMT2, a configuration resolved by DDM1 [32]. This explains the global gain of pericentromeric TE CHG methylation in the absence of H1 [32]. This global effect does not explain the derepression of H1-sensitive TEs. The observation that H1 can interact with CMT3 [59] suggests also a scenario in which H1-bound nucleosomes are favored for CMT3-mediated DNA methylation, thereby reinforcing silencing at these TE loci. Thus in *Arabidopsis*, H1 may play a dual role in restricting methylation through antagonizing DDM1 function [32] but also reinforcing methylation in a CMT3-dependent manner. Such a model, however, remains to be experimentally validated. Alternatively, and possibly not exclusively to other mechanisms, the absence of H1 may expose additional TEs to the DNA demethylase DEMETER; this was suggested for TE loci activated in the vegetative cells of pollen grains where the main H1 variants are naturally depleted [58]. In animals, interestingly, H1 depletion also leads to partial TE reactivation but following distinct mechanisms: independently of common epigenetic silencing marks in mammalian stem cells [60] or through an H3K9me2-dependant manner in *Drosophila* [14, 61].

Chromocenters containing ribosomal DNA repeats are not affected by the loss of H1 variants in *Arabidopsis*, further indicating the existence of an H1-independent control of such heterochromatic structures. It would be interesting to investigate whether other structural proteins such as H1-related, GH1-containing proteins [21] contribute NOR condensation or if hypoacetylation of the NOR nucleosomes by HDA6 [62] is self-sufficient.

In addition, dispersed heterochromatin regions still locate peripherally in the *3h1* nuclei indicating that not only structural compaction but also peripheral location is uncoupled from heterochromatin silencing in *Arabidopsis*. Interestingly, this role in chromocenter formation is not a shared, inherent feature of linker histones
<figure><img src="image_15.png" /></figure>

among eukaryotes. Indeed, depletion of the H1c and
H1d mammalian variants in mouse embryonic stem cells
leads to chromocenter clustering but not decondensa-
tion [60]. In *Drosophila*, H1 also promotes the conden-
sation of pericentromeric heterochromatin, but this
observation is specific to salivary gland cells which are
singular due to their polytenic chromosomes [20]. In
both *Drosophila* and mouse somatic cells, chromocenter
condensation seems instead to be controlled by the D1
and HMGA1 multiAT hook proteins, respectively [63].
Therefore, our study uncovers an influence of H1 in

chromocenter condensation which may differ from the
animal system.

H1 as a fine-scale architect of euchromatin organization
from nucleosomes to nanodomains

We have shown that global chromatin decondensation in
H1-depleted nuclei correlates with a greater heterogeneity
in the spatial distribution of chromatin domains in the nu-
cleoplasm. Yet, compact nanodomains as those visualized
by high- and ultrastructural microscopy imaging still form
in the absence of H1. This is consistent with the intrinsic
property of chromatin fibers to fold along various configurations even in the absence of linker histones: in vitro assembly assays have shown that nucleosomal arrays can be formed in the absence of H1, yet with a lesser degree of molecular organization, forming “ladder” or “puddle” type of arrangements [64]. Thus, our observations in vivo might relate to this intrinsic property of H1. Furthermore, the increased spatial dispersion of compact nanoscopic domains free of H1 in *3h1* nuclei is highly reminiscent of euchromatin organization in tumorigenic nuclei losing fractal property of organization [45].

At the molecular level, we observed that *Arabidopsis* H1 influences nucleosomal spacing, resulting in deregulated density over defined genomic regions. Notably, H1 depletion affects the structural differentiation, in terms of nucleosomal density, of chromatin states. This effect is distinct from that of a loss-of-function of a general chromatin factor regulating nucleosome assembly such as the histone chaperone CAF-1 (Additional file 1: Figure S6). In addition and strikingly, besides a shift of NRL distribution peak from 170–180 to 160–170 bp similar to the situation in animals [10], *3h1* chromatin harbors a higher fraction of both short NRLs (< 160 bp) or unexpectedly long (> 200 bp). This indicates a relaxed control of nucleosome distribution, with a more permissive (i.e. less constraint) nucleosome positioning [64]. H1 incorporation is thought to reduce the sterical occupation of linker DNA and to favor a compact zig-zag chromatin folding [65]. The question arises, how NRL longer than 197 bp can be achieved. One possibility is that such NRL corresponds to the MNase inaccessible regions spanning two closely associated nucleosomes; another possibility is suggested by the occurrence of the so-called stretched nucleosomes where DNA-histone contacts are lost in the absence of H1 [66]. Ultimately, variations in NRL and associated linker DNA length can contribute to chromatin fiber polymorphisms with in turn an influence on the formation of chromatin domains [67–69]. The relaxed NRL distribution hence variable nucleosome positioning landscape seen in *3h1* is thus likely responsible for the observed local variations in chromatin compaction and spatial heterogeneity.

How this fine-scale topology connects with chromatin function, notably gene expression, requires a biophysical appraisal of chromatin organization. Emerging models linking molecular-to-spatial (macromolecular) levels of chromatin dynamics are very attractive to explain the emergence of transcriptional dynamics during cellular differentiation [70]. Notably, local nucleosomal arrays appear to connect with molecular crowding, a dynamic property directly influencing the access and residence of the transcription machinery to chromosomal regions [57, 70]. In such a model, spatial heterogeneity of chromatin nanodomains seems to favor the occurrence of

transcriptional bursts such as those observed in an animal model of cellular differentiation [57, 70]. In addition, specific disturbances in spatial chromatin organization in cancerogenous cells have been correlated with transcriptional heterogeneity linked to the paired effect of increased accessibility and increased local compaction [57]. It is tempting to link these biophysical, molecular models of the nanoscale level of chromatin organization with the role of linker histones. Consistent with it, pluri-potent mammalian cells in which H1 levels are low display highly dispersed, small nucleosome clutches favoring RNA PolII redistribution [9]. Collectively, in light of these models and observations, we propose that the molecular and spatial distribution heterogeneity of H1-depleted chromatin in *Arabidopsis* resembles that of animal un-/de-differentiated cells. This proposal invites future studies on the role of *Arabidopsis* H1 variants in regulating transcriptional dynamics notably during cellular differentiation or in response to environmental cues.

## H1 defines a dual causality between chromatin organization and gene expression

Our cytological analyses clearly showed a profound impact of H1 on hetero- and euchromatin organizations at the molecular and spatial distribution level. But the daunting question remains whether this influences chromatin function and notably gene expression. Our study uncovers an unexpected, complex relationship between chromatin organization, probed by nucleosomal coverage, and gene expression. For instance, we found a lower nucleosome coverage in *3h1* along a large fraction of genes (excluding heterochromatic regions, discussed earlier), yet this reduced nucleosome occupancy did not trigger a significant upregulation in *3h1* plants. This was true for the genomic regions associated with H3K27me3 (e.g., chromatin state CS4, [40]), hence PRC2-target loci. By contrast, the genomic regions normally associated with transcriptional competence (e.g., CS1 enriched in H3K4me3 and H3.3) gained nucleosomal coverage in *3h1*, still without detectable consequences on most transcript levels. Thus, for a large fraction of the genome, both in eu- and heterochromatin contexts, H1-mediated chromatin organization appears to act downstream of the transcriptional controls. The observation that the structural differentiation of loci relative to the gene expression strength is blurred in *3h1* reinforces the conclusion that gene body nucleosome density is not epistatic to transcription.

These findings challenge a common view that transcription is influenced by nucleosome density, coming from the observation that nucleosome occupancy is inversely correlated to gene expression level ([51] and this study). Whether H1-mediated chromatin organization serves another purpose, for instance, facilitating the
rewiring of transcriptional programs when cells are en-
vironmentally challenged or securing transcriptional ro-
bustness is an exciting possibility to be investigated. If
they hold true, these concepts would provide an explan-
ation for the developmental phenotypes in *3h1* plants
with an apparent relaxed control of several developmen-
tal transitions. These deficiencies can be interpreted as
instable transcriptional states releasing seed dormancy,
flowering, lateral root organ initiation, root hair, and sto-
matal cell fate earlier than in the wild-type counterparts,
although the severity of such developmental phenotypes
might be attenuated by the buffering capacity of devel-
opmental processes inherent to plants [71].

The dual role of H1 is unveiled by the several hundred
of H1-sensitive loci for which chromatin organization
has a prime contribution with respect to transcriptional
control. Gene upregulation in *3h1* correlated with a not-
able nucleosome enrichment in the gene body but not at
the upstream elements. Yet, whether this average in-
creased coverage effectively associates with high tran-
scriptional rates or whether it is contributed by a
fraction of cells harboring silent loci cannot be resolved
in bulk tissue profiling. Possibly for these loci, increased
nucleosome mobility provides a chromatin template fa-
vorable for transcription and/or unfavorable to propagate
repressive epigenetic marks, as suggested by compromised
H3K27me3 maintenance in *3h1*.

A new connection between H1 and epigenetic
reprogramming

To a certain extent, chromatin organization and proper-
ties in *Arabidopsis* *3h1* mutant nuclei are reminiscent of
a pluripotent chromatin state described in different plant
cell types and studies: (i) chromatin decondensation and
reduced H3K27me3 abundance resemble the chromatin
status of H1-depleted spore mother cells (SMC) at the
somatic-to-reproductive fate transition [25, 72]; (ii) het-
erochromatin decondensation and redistribution of epi-
genetic marks such as those observed in H1-depleted
nuclei are hallmarks of in vitro de-differentiated cells
and at flowering transition [73, 74]; and (iii) high chro-
matin mobility in H1-depleted chromatin is reminiscent
of the chromatin state in meristematic cells [50].

These collective observations allow to propose the hy-
pothesis that H1 depletion contributes to reprogram-
ming of the chromatin landscape at the structural and
epigenetic level in plants cells undergoing fate reorienta-
tion as in the cases mentioned above. In support of this
hypothesis, which of course calls for further investiga-
tions, the stark decrease in H3K27me3 levels in H1-de-
pleted chromatin such as in *3h1* plant cells or in wild-
type SMC [25, 72]. H3K27me3 is notoriously associated
with the memory of epigenetic states and contributes

gene expression reprogramming during cell pluripo-
tency acquisition and differentiation in animal and
plants [38, 75–78]. Interestingly, H1 depletion in
mammalian cells, while moderately affecting global
transcriptome patterns as here in *Arabidopsis*, affects
more particularly pluripotency genes [79].

In addition, and consistently with our proposal, the devel-
opmental phenotypes observed in *3h1*—including delayed
release of seed dormancy, precocious flowering, deregulated
lateral root and root hair initiations, and the control in sto-
matal spacing—are reminiscent of phenotypes arising in
PRC2 loss-of-function mutants with altered H3K27me3
levels [38, 80–85].

The influence of H1-mediated chromatin structure on
H3K27me3 levels is a conserved feature in plants (as
shown by our study) and in animals [10] although this
may be cell type-dependent [13]. An explanation for this
may be a substrate preference of the PRC2 complex to
H1-containing nucleosomes as shown for the animal
EZH2 subunit in vitro [86] leading presumably to a
more efficient propagation of H3K27me3 in chromatin
states structured by H1 in vivo. Conversely, H1-depleted
cells may inefficiently propagate H3K27me3 landscapes,
which over numerous divisions may lead to drastic loss,
as measured in (terminally differentiated) plant tissues
(this study) and in mouse tissues [10]. Involving H1 in
H3K27me3 maintenance rather than the establishment
is a hypothesis supported by the observation that H1-de-
pleted plants are not impaired in basic growth, nor or-
ganogenesis or body patterning processes, but, instead,
show altered controls of developmental and cellular
transitions normally modulated by PRC2.

Conclusions

This work explored for the first time the functional rela-
tionship between fine-scale chromatin organization at
the molecular and spatial levels, transcriptional control
(Fig. 5), and phenotypic impact in *Arabidopsis*. Our ana-
lyses of the impact of H1 depletion in *Arabidopsis* unveil
a remarkable evolutionary convergence of H1 function
in plant and animal kingdoms despite a diversification of
the H1 variants family [87]. In *Drosophila* and mamma-
lian cells, H1 depletion generates large-scale chromatin
alteration including blurring of topological domains, re-
laxation of nucleosome distribution periodicity, and im-
paired propagation of the epigenetic marks [10, 88–90].
In *Arabidopsis*, the misregulation of a small proportion
of genes in the absence of H1 is intriguing when con-
sidering the massive structural chromatin alterations.
While it suggests conserved functional principles for H1,
it also unveils a dual, causal relationship between tran-
scriptional controls and H1-mediated chromatin organ-
isation. Our findings prompt revisiting the common
view that chromatin compaction and nucleosome
density influence gene expression in a major way. This
idea may only hold true for a subset of H1-sensitive
genes and TEs. But generally, transcriptional regulation
remains epistatic to H1-mediated chromatin
organization for a majority of heterochromatin and eu-
chromatin loci. Based on these findings, we propose that
chromatin organization in the absence of H1 loses versa-
tility and is less prone to rapid reprogramming at the
structural, epigenetic and transcriptional level. This
property is not challenged under controlled growth con-
ditions hence requires testing in modular and vari-
able environments. We suggest a working model, in
which linker histones act as fine-scale chromatin archi-
tects with a dual function in (i) the facilitation of tran-
scriptional reprogramming (in development or under
environmental cues) and (ii) securing stable transcrip-
tional landscapes buffering against variability. Our find-
ings echo with recent models in mammalian cells
implicating higher-order folding chromatin topology as
an independent route influencing transcriptional dynam-
ics [57] and in facilitating transcriptional robustness in
differentiation [70].

**Methods**

**Plant materials and growth conditions**

The *Arabidopsis thaliana* plants used in all experiments
were in the Col-0 background unless it is specified
otherwise. The *h1.1h1.2h1.3* (*3h1*) mutant was described
before [25, 28]. The snapshots from RNAseq reads of
H1 variant expressions in *3h1* and wt are presented in
Additional file 1: Figure S17. The mutant showed no de-
tectable levels of H1 in immunostaining experiments
[25] and in Western blot (www.agrisera.com/en/artiklar/
h1-histone-h1.html). Complemented mutant lines were
generated by transforming *3h1* via floral dip method
[91] with H1-tagged variants (prom. H1.1::H1.1-RFP,
prom. H1.2::H1.1-(G/C) FP, prom. H1.3::H1.3-GFP) de-
scribed previously [25, 28]. The *3h1* was complemented
with either two main (H1.1, H1.2) or all three H1 vari-
ants to generate the following lines: *3h1-comp*<sup>1,2=</sup>
* h1.1h1.2h1.3; H1.1-RFP; H1.2-GFP* (line #KR276), *3h1-
comp*<sup>1,2,3</sup> = *h1.1h1.2h1.3; prom. H1.1::H1.1-RFP; prom.
H1.2::H1.2-CFP; prom. H1.3::H1.3-GFP* (lines #KR264
and #KR265). For FRAP experiments, the UBQ10::H2B-
RFP [48] was crossed with *3h1*, and in the subsequent
generations, by genotyping the *3h1/UBQ10::H2B-RFP*
and wt segregants were identified.

Seeds were surface sterilized and rinsed in sterile water
before transferring onto germination medium (0.5× MS
medium, 0.8% agar). They were placed on the medium
using toothpicks to ensure a uniform distribution, strati-
fied 2–4 days at 4 °C, and transferred into a plant growth
incubator (Percival, Germany) with long-day photo-
period (16 h, 22 °C day/8 h, 18 °C night) and light flux

around 120 µM s⁻¹ m⁻² for routine experiments. Growth
of calli and scoring of lateral root production were tested
under continuous light (light flux around 100 µM s⁻¹ m⁻²,
Aralab FitoClima 1200). When the flowering stage was
necessary, the 10-day-old seedlings were transferred into
the soil and grown at 19–21 °C with a 16-h day/8-h night
photoperiod.

**Chromatin analyses and immunostaining**

The nuclei area, heterochromatin (RHF, CCs), and im-
munostaining analyses were carried out essentially as de-
scribed [92] with minor modifications. The nuclei were
isolated from rosette leaves of 3–4-week-old seedlings;
per extraction, five leaves were fixed during 20 min
under vacuum in a fresh 4% formaldehyde solution prior
to isolation and resuspension of the nuclei in a final vol-
ume of 1 mL nuclei isolation buffer (NIB). DAPI was
added at a concentration of 0.1 mg/mL for flow sorting
according to DNA content. Diploid (2C) nuclei have
been flow-sorted using a BD FACSAria IIIu flow cyt-
ometer with a 450/50-nm filter (405 nm laser), equipped
with a 100-µm nozzle and 25 Psi pressure. The nuclei
were collected in 200 µL of NIB before spreading on
Superfrost plus slides (1000 nuclei per slide) and stored
at 4 °C until use. Mutant and wild-type plants were
grown and processed for nuclei isolation and immuno-
staining in parallel.

For heterochromatin analysis, the slides were rinsed in
SSC2X then PBS before staining with DAPI 1 µg/mL in
Vectashield (Vector Laboratory). For immunostaining,
the protocol essentially followed previously described
steps [92]. As primary antibodies, rabbit anti-Histone
H3 (Abcam; ab1791), anti-Histone H1 (Agrisera;
as111801), anti-H3K27me3 (Active Motif; 39155), anti-
H3K27me1 (Abcam; ab113671), antiH3K4me3 (Abcam;
ab8580), anti-H3K9ac (Abcam; ab10812), and anti-
H3K9me1 (Abcam8896) were used at a dilution of 1:200
and incubated at 37 °C for 1 h. As secondary antibody,
Alexa Fluor 488-conjugated goat anti-rabbit IgG (Mo-
lecular Probes; A-11008) was used at a dilution of 1:
1000 and incubated for 2 h at 37 °C. The nuclei were
counterstained for DNA with propidium iodide (PI).

Sixteen-bit images were acquired using a Leica TCS
SP5 confocal laser scanning microscope (CLSM) (Leica
Microsystems, GmBH, Germany) using a ×63 GLY lens
(NA 1.4) for heterochromatin and immunostaining ana-
lyses. Exposure times, illumination intensities, zoom fac-
tor, scanning speed, and pinhole were kept identical for
the image series in an experiment. For RHF measure-
ments, signal intensities were recorded in manually
drawn ROIs capturing chromocenters and normalized
over the whole nucleus intensity using Fiji [93]. For im-
munostaining, the signal intensities for antibodies were
normalized against PI levels. Graphs were plotted in
Excel, and the data were statistically assessed using a Student *t* test (unpaired, unequal variance) for comparing wild-type and mutant samples.

## Fluorescence in situ hybridization and 3D image processing
### FISH analysis of leaf nuclei

The nuclei were isolated from leaves of 35-day-old rosettes grown under a 16-h/8-h photoperiod. Nuclei extraction and embedding in acrylamide gel pads on the slide were done as described [42]. Centromeric and 45S rDNA repeats were detected by FISH using pAL1 and pTA9 to generate DNA probes, respectively [39]. FISH was done as described [42] with the following labeling kits and fluorescent immunolabeling reagents: DIG-Nick (Sigma Aldrich, 11745816910), mouse IgG anti-DIG (1:250, Sigma Aldrich, 11333062910), goat IgG anti-mouse IgG~Alexa 488 (1:200, Life Technologies, A11001); Biotin-Nick translation kit (Sigma Aldrich, 11745824910), Biotinylated Anti-Avidin D (1:250, Vector Labs, BA-0300), and Texas Red Avidin D (1:1000, Vector Labs, A-2006). The nuclei were counterstained for DNA with DAPI in Vectashield (Vector Laboratory). FISH signals in 3D nuclei were imaged using Stimulated emission depletion (STED) microscopy (Leica SP8R WL 3xSTED, Leica microsystems, Germany).

### FISH analysis of cotyledon nuclei

The nuclei were isolated from dissected cotyledons of 5-day-old seedlings grown under a 16-h/8-h photoperiod. Nuclei extraction, fixation, and hybridization with pAL1-derived and F28D6-derived (180-bprepeats) probes [39] were performed as previously described [94]. The slides were washed and mounted in Vectashield with 2 µg/µL DAPI, and image acquisition was performed as in [94].

## Nuclei isolation for MNase-seq

The nuclei were isolated from 3-week-old seedlings frozen in liquid nitrogen as previously described [95] with the following modifications: after resuspending in HBB, the nuclei were applied to a layer of HBB with 40% Percoll (GE Healthcare); centrifuged at 1000g, 6 min; resuspended in HBB; applied to 40/75% Percoll gradient; centrifuged at 400g, 40 min; collected; and washed three times with HBC. The integrity of the extracted nuclei was monitored using DAPI staining and fluorescence microscopy. The quantity of nuclei was measured by qPCR with primers targeting nuclear DNA.

Digestion was performed by incubating the nuclei suspended in DB buffer (16 mM Tris-HCl pH = 7.6, 50 mM NaCl, 2.5 mM CaCl₂, 0.01 mM PMSF, 1x Complete EDTA-free Protease Inhibitors (Roche)) with 1.5 µL (final concentration 0.3 U/µL) of micrococcal nuclease (Thermo Fisher), and 2 µL (final concentration 0.2 U/µL) of RNase A (Thermo Fisher) at 8 °C for 90 min with

gentle mixing. The reaction was stopped by adding an equal volume of 2x Lysis buffer with EDTA (100 mM Tris-HCl pH = 8, 200 mM NaCl, 50 mM EDTA, 1% SDS). The samples were lysed by incubation at 37 °C for 60 min with shaking (1000 rpm). DNA was purified using phenol-chloroform extraction, precipitated with isopropanol and sodium acetate, and resuspended in water.

DNA was size selected by electrophoresis on 2% agarose gel with 1x TAE buffer with SYBR Gold (Invitrogen) stain. The mononucleosomal band was excised, frozen, and squeezed by 3 cycles of spinning and rehydration on centrifuge column. DNA was purified and concentrated using Agencourt AMPure XP beads (Beckman Coulter). Barcoded libraries were synthesized from 100 ng of mononucleosomal DNA using Ion Xpress™ Plus gDNA Fragment Library Preparation Kit and Ion Xpress™ Barcode Adapters. DNA was end-repaired prior to adapter ligation and size selection, and amplification steps were omitted. The resulting libraries were quantified with Ion Library Quantitation Kit, pooled, and used to prepare the template by clonal PCR with Ion PI™ Template OT2 200 Kit v3 on Ion OneTouch™ 2 System. Sequencing was performed on Ion PI™ chip v2 and Ion Proton™ sequencer using Ion PI™ Sequencing 200 Kit v2 (all Ion Torrent kits and software are trademarks of Thermo Fisher).

## FRAP imaging and data analyses

A *promUBQ10::H2B-RFP* marker [96] was introgressed in *3h1*-mutant plants by crossing. Both wild-type and triple-mutant segregants were analyzed. Measurements were done on the root tips of 2-week-old seedlings grown as previously described. One sample was prepared at a time: the root was excised and delicately mounted (i.e., without squashing) in 0.5x MS between the slide and coverslip (precleaned with EtOH), sealed with transparent nail polish, and let to equilibrate upside down for 10 min on the microscope platform before measurements. The imaging chamber was set at a constant temperature of 20 °C (higher/fluctuating temperatures induce nuclei juggling). Bleaching and imaging were done using an APO PL × 40 oil immersion objective, NA 1.3, over a single plane capturing an optical section of ~ 2 µm encompassing a single nucleus (pinhole opening to 5 AU) with a 256 × 256 pixels image format, threefold zoom factor. Bleaching was performed in eucromatin within ROI of 1 µm diameter using five or more pulses until near-total bleach was obtained (Argon laser at 80% power, 100% transmission in 488 nm), and post-bleach images were recorded using 5–7% laser transmission for excitation, in a series of 10 time points, 1-s interval, followed by 10 time points, 60-s interval. For analyzing the fluorescence recovery, images were first corrected for nuclear drifts occurring
during acquisition, using a rigid registration approach
in Fiji ([93], plugin/registration/stack reg/rigid trans-
formation). When a single image captured several nu-
clei, single nuclei were cropped for registration and
analysis. Fluorescence measurements were done on the
bleach ROI, a control ROI near and outside the nu-
cleus, and over the whole nucleus. Calculation of fluor-
escence recovery was done as described in [50, 97]
whereby the initial intensity was normalized at 1 for
each image before average calculation.

TEM sample preparation, imaging, and image analysis

Seventy-nanometer tissue sections were prepared
from 2-week-old seedling roots, using a high-pres-
sure freezing/freeze substitution and uranyl acetate
staining approach as described in details previously
[46]. The sections from the elongation zone were se-
lected for the analysis (i.e., meristematic zone was
avoided), and nuclei pictures were consistently re-
corded from the epidermal layer at the 24,500-fold
magnification yielding a resolution of 1 pixel = 2 nm
in our setup. For the analysis, the square regions of
interests (ROIs) of similar size (ca 800 × 800 ± 200
pixels) were captured in the euchromatin regions
(i.e., excluding strongly staining chromocenters) for
the analysis. We used a spatial pattern analysis ap-
proach that was previously validated to capture rele-
vant structural features of chromatin organization in
cancerogenous animal cells [45]. Spatial autocorrel-
ation analysis delivers a mathematical model of chro-
matin density distribution for each ROI with respect
to the physical length scales within which signal pat-
terns (i.e., local objects of similar intensities) are re-
peated in a regular pattern (periodicity) [45]. We
used a user-friendly graphical interface developed in
Matlab for batch processing of multiple ROIs avail-
able at https://github.com/barouxlab/ChromDensity-
Nano and described in details previously [46].

Analysis of developmental transitions

Flowering time

Plants for flowering experiments were grown in the greenhouse or growth chamber under the long daylight regime. To avoid positional effect, different genotypes were always randomly arranged over the growth area. The number of rosette leaves was counted when the inflorescence was about 0.5 cm long. Three replicate experiments were conducted: replicate 1 was based on 40 to 50 plants depending on genotype as detailed in Additional file 1: Figure S1A; replicates 2 and 3 were based on 10 and 12 plants per genotype, respectively, as detailed in Additional file 1: Figure S1C,D.

Root length and lateral root scoring

Seedlings were grown vertically on square Petri dishes under a continuous light regime. The plates were scanned 8 days after germination to score for the number of lateral roots. Root (main and lateral) lengths were scored using manual vector tracing in Fiji, reported at scale [93]. For microscopic observations of lateral root primordia, 5-day-old seedlings grown under continuous light were fixed in 70% ethanol, rinsed once in sterile water, and mounted in water on microscope slides (five roots aligned/slide covered with 40 × 22 mm coverglass). Primordia were scored according to published developmental scale [98]. Graphs were plotted in R.

Root hair observations

Increased root hair number was visible under growth
conditions as described above for lateral root scoring.
The phenotype was however increased using MS
medium supplemented with 5% sucrose, and this
medium was used to take the micrographs in Fig. 1.

Stomata patterning

Fresh epidermal peals of 14-day-old cotyledons were
mounted in water. Images of the adaxial surface were re-
corded with DIC microscopy, and stomatal clusters were
scored following as described [99].

Seed dormancy

The experiment was designed as described previously
[100] with minor modifications. Plants were grown in a
growth chamber under long daylight regime with con-
trolled humidity. Freshly harvested seeds were collected
and stored under constant conditions. Around 180 seeds
per plant were placed on a wet filter paper in a Petri dish
and incubated in the growth incubator at 22°C under
long day light regime. After 3 days, the number of emer-
ging radicles was counted. For the time point “day 1,”
the seeds were used 1 day after harvesting. For the time
point “3 weeks,” the seeds from the same batch were
used 3 weeks after harvesting.

Callus induction

Cotyledons from 7-day-old seedlings grown under a 16-
h/8-h photoperiod were excised, transferred onto callus
induction medium (CIM, Gamborg B5, 0.05% MES, 2%
glucose, 0.1 mg/L kinetin, 0.5 mg/L 2,4-D), and let to de-
velop for 5 weeks under a 16-h/8-h photoperiod. Callus
size (area) was determined from images using manually
drawn contours in Fiji [93]. Graphs were plotted in R.

Immunoblot analyses

Seeds were surface sterilized in 70% ethanol 0.05% SDS for 3 min and rinsed with 90% ethanol before drying and plating on MS medium supplemented with 0.5% sucrose
and 0.9% agar. Eight-day-old seedlings were used for chromatin extraction protocol as described previously [101]. Forty micrograms of protein samples, as estimated by Pierce BCA Protein Assay Kit (Thermo Scientific), was loaded on 14% LiDs Tris-Tricine gels and blotted onto Immobilon-P membranes (Millipore) before immunodetection using antibodies recognizing either unmodified histone H4 (Millipore #05-858) or H3K27me3 (Millipore #07-449), H3K4me3 (Millipore #07-745R), and H3K9ac (Millipore 06-942, lot: 31636). Chemiluminescent signals were detected and quantified using a LAS4000 luminescence imager (Fuji) in three biological replicates.

### RNA-seq and differential gene expression analyses

RNA was isolated using modified TRIzol method [102] in three biological replicates for wt and *3h1* from 3-week-old seedlings pooled and frozen in liquid nitrogen. Ribosomal RNA was removed using RiboMinus Plant Kit (Thermo Fisher), and ERCC RNA Spike-In Mix 1 (Thermo Fisher) was added. Libraries were prepared with Ion Total RNA-Seq Kit v2 and Ion Xpress RNA-Seq Barcode 1-16 Kit according to the user guide. Sequencing template was generated with Ion PI™ Template OT2 200 Kit v3 on Ion OneTouch™ 2 System. Sequencing was performed on Ion PI™ chip v2 and Ion Proton™ sequencer using Ion PI™ Sequencing 200 Kit v2 (all Ion Torrent kits and software are trademarks of Thermo Fisher).

Base calling and adapter trimming were performed automatically by Torrent Suite software. Residual rRNA and ERCC reads were identified and filtered out using bbsplit and filterbyname scripts from BBTools suite [103]. Reads were aligned to TAIR10 genome using TMAP 5.0.13. with soft clipping from both ends and returning all the mappings with the best score. Other settings were set according to Torrent Suite defaults. Unaligned reads were aligned with BBMap (Brian Bushnell).

Quantitation to ARAPORT11 transcripts was performed in Partek Flow (Partek Inc.) using Partek E/M (expectation-maximization) algorithm. This allows inclusion of multi-mapping reads into quantitation. The algorithm quantitates unique alignments first and then iteratively divides the value of each multi-mapping read between possible loci based on the quantitation value in their close vicinity. The quantitation was performed on gene level (one value per gene). The read counts were normalized by the trimmed mean of M values (TMM) method [104] followed by the UQ (upper quartile) method [105]. Differential expression analysis was performed with the Partek gene-specific analysis (GSA) algorithm. It is a sophisticated version of global linear model, fitting model of variables, and theoretical distribution of expression values to each gene independently. The model was allowed to use normal and log-normal

distributions to fit to the data. For log-normal distribution, the “shrinkage” method was allowed. This method, published as limma Voom [106] is specifically designed for read-based quantitative studies with low replication. It works by calculating the average trend of gene values' standard errors in relation to average gene expression levels. Next, the errors observed for each gene are normalized towards the trend. This approach was developed to reduce the occurrence of false-positive and false-negative rate of differentially expressed gene calls, especially in low replication studies and genes with low expression, where high or low dispersion of quantitated values can result from sampling effects. Pearson correlation between replicate samples is shown in Additional file 1: Figure S10B.

To divide *Arabidopsis* genes into quintiles of expression, the quantified values from Partek E/M algorithm (described above) were normalized using reads per kilobase of million (RPKM) mapped reads method [107]. The expression values were averaged across replicates. The genes with zero values of expression in each genotype were sorted into a separate group. The remaining genes were divided into quintiles based on the expression levels separately for wt and *3h1* datasets.

Boxplots were generated with tables of RPKM values using the “boxplot” function in R language.

To calculate the distribution of genes upregulated in *3h1* versus all genes in *Arabidopsis thaliana* genome across chromatin states, genomic coordinates for chromatin state (CS) locations across the genome were downloaded from published data [40]. The genomic coordinates for genes upregulated in *3h1* were taken from TAIR 9 to be consistent with CS coordinates. Then, for each gene, the percentage of overlapped chromatin states was calculated, and for the final graph, the summary of all analyzed genes was presented.

### MNase-seq data analysis

Base calling and adapter trimming were performed automatically by Torrent Suite software. Reads were aligned with TMAP 5.0.13. Soft clipping was turned off, end repair was allowed, and all alignments for multi-mapping reads were reported. Other settings were set according to the Torrent Suite defaults. Multi-mapping read positions were resolved using MMR [108] with default settings.

Peak calling was performed on reads reaching terminal adapter with length range between 147 and 220 nt using a set of custom-made Python scripts. First read centers were piled up, and then the highest coverage positions were selected using greedy algorithm. Ends of the longest read used to define a position were used as peak boundaries. A peak was called only if its boundaries were not overlapping those of the neighboring peak. NRLs were defined by calculating peak-to-peak distances from peak calling results. The frequency of distances
was calculated as a percentage of all measurements and binned into three groups. Histograms were plotted in Microsoft Excel.

Quantile-normalized wiggle occupancy files were generated with DANPOS2 [109] using the dpos function with default settings. To avoid shifting of read positions (automatic procedure for single-end reads), the program was fed with “fake 75-nt paired-end” bed files, generated from both ends of alignments of fully sequenced reads. To create fake fragments, the genomic interval of each mapped read (single-end) was divided into two 75-nt intervals on opposite strands that were used to feed the Danpos algorithm as if those were paired-end data. This did not change the informative, real read data content but enabled to circumvent a limitation of Danpos imposing read shifting for single-end reads, which is not appropriate for single-end sequencing of full fragments after MNase digestion. Using these “fake fragments,” we could thus switch off the read shifting function in Danpos for further data processing.

Wig files were converted to BigWig format using UCSC wigToBigWig [110] and used in deepTools [111] for plotting.

Nucleosome occupancy was calculated as the number of reads, adjusted to a length of 80 nt summed over bins of 10 nt, smoothed by 20-nt running window, and quintile normalized.

## Filtering out *Ler* residual sequences

Despite series of five backcrosses after introduction of *h1.3*-mutant allele from *Ler* background into our *h1.1h1.2h1.3* (Col-0) line, some residual *Ler* sequences were still present, mainly neighboring the *H1.3* gene. To avoid interference from those sequences in our analyses, we identified their precise genomic coordinates using SNP and coverage analyses by comparing to sequenced genome of parent *Ler* *h1.3* line. We used those coordinates to generate bed files and filter out all reads overlapping residual *Ler* sequences using bedtools intersect [112].

## Additional files

**Additional file 1:** Supplementary Figures S1–S17. (DOCX 7482 kb)

**Additional file 2:** Table S1. Quantifications of chromatin cytology corresponding to Figs. 2, 3 and 4. (DOCX 17 kb)

**Additional file 3:** Table S2. Transposable element (TE) expression in *3h1*. Available as an Excel table. (XLSX 64 kb)

**Additional file 4:** Table S3, Table S5, Table S6. Classes of TEs upregulated in *3h1*; Gene Ontology (GO) analysis of genes which are misregulated in *3h1* mutant; Expression of histone modifying enzymes in *3h1*. (PDF 903 kb)

**Additional file 5:** Table S4. Gene expression in *3h1*. Available as an Excel table. (XLSX 4453 kb)

**Additional file 6:** Review history. (DOCX 57 kb)

## Abbreviations

*3h1:* Triple H1 knockout mutant *h1.1*, *h1.2*, *h1.3*; Col: Columbia (*Arabidopsis* accession); CS: Chromatin state; FISH: Fluorescence in situ hybridization; H1: Histone H1; H3: Histone H3; H3K27me3: Histone H3 lysine 27 trimethylation; *Ler*: Landsberg *erecta* (*Arabidopsis* accession); PRC2: Polycomb-group Repressive Complex 2; TE: Transposable Element

## Acknowledgements

We are indebted to Ueli Grossniklaus (UG; University of Zurich) for the valuable advice and insightful discussions, for sharing laboratory facilities, and for the financial support; to Chris Bowler and Vincent Colot (Ecole Normale Supérieure Paris) for the constant support and helpful discussions to FB; to Valeria Gagliardini, Christof Eichenberger, Arturo Bolaños, Daniela Guthörl, and Peter Kopf (University of Zurich) for the general lab support; to Daniel Prata (University of Zurich) for the assistance in plant growth and selection; to Stefan Wyder (University Research Priority Program “Evolution in Action”, University of Zurich) for the advice on the statistical tests; to Lea Wolff (FB Lab, Ecole Normale Supérieure Paris -IBENS CNRS 8197 Paris) for the assistance in Western blot experiments; to Dominique Bergman (Stanford University) for the advice on stomata analysis; to Matthieu Ingouff (IRD Montpellier) for sharing the DynaMET lines; to Alexis Maizel and Daniel von Wangenheim (University of Heidelberg) for the pUBQ10:H2B-RFP line; to Mariamawit Ashenafi (University of Zurich) for the technical assistance in FISH and immunostaining; and to Marek Kalinowski for the scripts to analyze chromatin state distributions. We thank Urs Ziegler, Jana Dohner, Andres Kaech, Claudia Dumrese, and Moritz Kirschman (University of Zurich) for the service and support at the Cytometry Facility and the Center for Microscopy and Image Analysis for nuclei sorting, STED, and GSD imaging, TEM preparation, FRAP imaging, and advice on data analysis. We also thank Marta Kobłowska (Institute of Biochemistry and Biophysics PAS) for sharing the facilities of Laboratory of Microarray Analysis for RNA-seq and MNase-seq experiments and Maciej Kotliński (University of Warsaw) for the conceptual support in designing and analyzing for RNA-seq and MNase-seq experiments.

## Review history

The review history is available as Additional file 6.

## Authors' contributions

KR and ML contributed equally to this work. ML and KR were responsible for the profiling of the experiments and data analysis (MNase-seq and RNAseq). SG was responsible for the data analysis supplement (S7B). BM was responsible for the cytological work (immunostaining, quantifications, RHF analyses). CB was responsible for the RHF analyses. IM was responsible for the cytology supplement (S4B). GT, MAK, and FB were responsible for the Western blot analysis. TNF and CB were responsible for the TEM preparations and images. CB and LC (script improvement) were responsible for the image processing and analysis. KR (callus, dormancy, flowering, root), JS (flowering, root), SF (root measurements), and CB (root hair, stomata) were responsible for the plant phenotype analysis. CB and KR were responsible for the reporter line imaging. CB was responsible for the FRAP experiment. CB, AJ, and FB acquired funding for the study. CR supervised and funded TNF. CB, AJ, FB, KR, and ML designed the experiments and wrote the manuscript. All authors read and approved the final manuscript.

## Funding

This work was supported by the University of Zurich, grants of the Swiss National Science Foundation to CB (SNF 310013A_149974) and UG (310030B_160336), grants of SystemsX.ch through an iPhD project to CB (2014/235) and the RTD project MecanX to UG and CR, fellowships from the "Forschungskredit" of the University of Zurich and the COFUND Program PLANT FELLOWS (2010-267243) of the European Union's 7th Framework Program to KR, support from the Polish National Science Centre to ML (2011/01/N/NZ2/04849), a grant of the Ministry of Sciences and Higher Education to AJ (MNIsw/PO4A/03928), the "Investissements d'Avenir" program launched by the French government and implemented by the "Agence Nationale de la Recherche (ANR)" to FB (ANR-10-LABX-54 MÉMOLIFE, ANR-10-IDEX-0001-02 PSL\*Research University, ANR-18-CE13-0004-01 Chromalight), support to FB lab by the EPIPLANT Groupement de Recherche, a PhD fellowship from the Université Paris-Sud Doctoral School in Plant Sciences to GT, and COST Action CA16212 INDEPTH.
**Availability of data and materials**

The data for MNase-seq and RNA-seq discussed in this publication have been deposited in NCBI's Gene Expression Omnibus [113] and are accessible through GEO Series accession number GSE113558 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE113558) [114].

The image data used for cytological quantifications and FRAP analyses are available at https://figshare.com/articles/3h1 vs WT_nuclei/8223977 [115]. The Matlab-script to analyze chromatin density distribution is freely available at https://www.github.com/barouxlab/ChromDensityNano.

**Ethics approval and consent to participate**

Not applicable.

**Consent for publication**

Not applicable.

**Competing interests**

The authors declare that they have no competing interests.

**Author details**

¹Institute of Plant and Microbial Biology, Zürich-Basel Plant Science Center, University of Zürich, Zürich, Switzerland. ²Institute of Biochemistry and Biophysics, Polish Academy of Sciences, Pawinskiego 5a, 02-106 Warsaw, Poland. ³Institute of Molecular, Cellular & Developmental Biology, Yale University, 352a Osborn memorial laboratories, New Haven CT, 06511, USA.

⁴Département de Biologie, IBENS, Ecole Normale Supérieure, CNRS, INSERM, PSL Research University, 46 rue d'Ulm, F-75005 Paris, France. ⁵College of Inter-Faculty Individual Studies in Mathematics and Natural Sciences, University of Warsaw, 02-089 Warsaw, Poland. ⁶Department of Biomedical Engineering, Northwestern University, Evanston, IL 60208, USA. ⁷Faculty of Biology, University of Warsaw, Pawinskiego 5a, 02-106 Warsaw, Poland.

Received: 10 January 2019 Accepted: 21 July 2019

Published online: 07 August 2019

**References**

1. Kasinsky HE, Lewis JD, Dacks JB, Ausio J. Origin of H1 linker histones. FASEB J. 2001;15(1):34-42.

2. Bednar J, Garcia-Saez I, Boopathi R, Cutter AR, Papai G, Reymer A, Syed SH, Lone IN, Tonchev O, Crucifix C, Menoni H, Papin C, Skoufias DA, Kurumizaka H, Lavery R, Hamiche A, Hayes JJ, Schultz P, Angelov D, Petosa C, Dimitrov S. Structure and dynamics of a 197 bp nucleosome in complex with linker histone H1. Mol Cell. 2017;66(3):384-97 e388.

3. Zhou BR, Feng H, Kato H, Dai L, Yang Y, Zhou Y, Bai Y. Structural insights into the histone H1-nucleosome complex. Proc Natl Acad Sci U S A. 2013; 110(48):19390-5.

4. Fyodorov DV, Zhou BR, Skoultschi AI, Bai Y. Emerging roles of linker histones in regulating chromatin structure and function. Nat Rev Mol Cell Biol. 2018; 19(3):192-206.

5. Bustin M, Catez F, Lim JH. The dynamics of histone H1 function in chromatin. Mol Cell. 2005;17(5):617-20.

6. Almeida R, Fernandez-Justel JM, Santa-Maria C, Cadoret JC, Cano-Aroca L, Lombrana R, Herranz G, Agresti A, Gomez M. Chromatin conformation regulates the coordination between DNA replication and transcription. Nat Commun. 2018;9(1):1590.

7. Christophorou MA, Castelo-Branco G, Halley-Stott RP, Oliveira CS, Loos R, Radzisheuskaya A, Mowen KA, Bertone P, Silva JC, Zernicka-Goetz M, Nielsen ML, Gurdon JB, Kouzarides T. Citrullination regulates pluripotency and histone H1 binding to chromatin. Nature. 2014;507(7490):104-8.

8. Hergeth SP, Schneider R. The H1 linker histones: multifunctional proteins beyond the nucleosomal core particle. EMBO Rep. 2015;16(11):1439-53.

9. Ricci MA, Manzo C, Garcia-Parajo MF, Lakadamyali M, Cosma MP. Chromatin fibers are formed by heterogeneous groups of nucleosomes in vivo. Cell. 2015;160(6):1145-58.

10. Fan Y, Nikitina T, Zhao J, Fleury TJ, Bhattacharyya R, Bouhassira EE, Stein A, Woodcock CL, Skoultschi AI. Histone H1 depletion in mammals alters global chromatin structure but causes specific changes in gene regulation. Cell. 2005;123(7):1199-212.

11. Seymour M, Ji L, Santos AM, Kamei M, Sasaki T, Basenko EY, Schmitz RJ, Zhang X, Lewis ZA. Histone H1 limits DNA methylation in Neurospora crassa. G3 (Bethesda). 2016;6(7):1879-89.

12. Yang SM, Kim BJ, Norwood Toro L, Skoultschi 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(5):1708-13.

13. Geeven G, Zhu Y, Kim BJ, Bartholdy BA, Yang SM, Macfarlan TS, Gifford WD, Pfaff SL, Verstegen MJ, Pinto H, Vermunt MW, Creyghton MP, Wijchers PJ, Stamatoynnopoulos JA, Skoultschi AI, de Laat W. Local compartment changes and regulatory landscape alterations in histone H1-depleted cells. Genome Biol. 2015;16:289.

14. Lu X, Wontakal SN, Kavi H, Kim BJ, Guzzardo PM, Emelyanov AV, Xu N, Hannon GJ, Zavadil J, Fyodorov DV, Skoultschi AI. Drosophila H1 regulates the genetic activity of heterochromatin by recruitment of Su (var)3-9. Science. 2013;340(6128):78-81.

15. Ausio J. Are linker histones (histone H1) dispensable for survival? Bioessays. 2000;22(10):873-7.

16. Jedrusik MA, Schulze E. A single histone H1 isoform (H1.1) is essential for chromatin silencing and germline development in Caenorhabditis elegans. Development. 2001;128(7):1069-80.

17. Patterton HG, Landel CC, Landsman D, Peterson CL, Simpson RT. The biochemical and phenotypic characterization of Hho1p, the putative linker histone H1 of Saccharomyces cerevisiae. J Biol Chem. 1998; 273(13):7268-76.

18. Shen X, Yu L, Weir JW, Gorovsky MA. Linker histones are not essential and affect chromatin condensation in vivo. Cell. 1995;82(1):47-56.

19. Ushinsky SC, Bussey H, Ahmed AA, Wang Y, Friesen J, Williams BA, Storms RK. Histone H1 in Saccharomyces cerevisiae. Yeast. 1997;13(2):151-61.

20. Lu X, Wontakal SN, Emelyanov AV, Morcillo P, Konev AY, Fyodorov DV, Skoultschi AI. Linker histone H1 is essential for Drosophila development, the establishment of pericentric heterochromatin, and a normal polytene chromosome structure. Genes Dev. 2009;23(4):452-65.

21. Kotlinski M, Knizewski L, Muszewska A, Rutowicz K, Lirski M, Schmidt A, Baroux C, Ginalski K, Jerzmanowski A. Phylogeny-based systematization of Arabidopsis proteins with histone H1 globular domain. Plant Physiol. 2017; 174(1):27-34.

22. Wierzbicki AT, Jerzmanowski A. Suppression of histone H1 genes in Arabidopsis results in heritable developmental defects and stochastic changes in DNA methylation. Genetics. 2005;169(2):997-1008.

23. Ingouff M, Selles B, Michaud C, Vu TM, Berger F, Schorn AJ, Autran D, Van Durme M, Nowack MK, Martienssen RA, Grimanelli D. Live-cell analysis of DNA methylation during sexual reproduction in Arabidopsis reveals context and sex-specific dynamics controlled by noncanonical RdDM. Genes Dev. 2017;31(1):72-83.

24. She W, Baroux C. Chromatin dynamics in pollen mother cells underpin a common scenario at the somatic-to-reproductive fate transition of both the male and female lineages in Arabidopsis. Front Plant Sci. 2015;6:294.

25. She W, Grimanelli D, Rutowicz K, Whitehead MW, Puzio M, Kotlinski M, Jerzmanowski A, Baroux C. Chromatin reprogramming during the somatic-to-reproductive cell fate transition in plants. Development. 2013;140(19): 4008-19.

26. Perez-Montero S, Carbonell A, Azorin F. Germline-specific H1 variants: the "sex'y" linker histones. Chromosoma. 2016;125(1):1-13.

27. Hsieh PH, He S, Buttress T, Gao H, Couchman M, Fischer RL, Zilberman D, Feng X. Arabidopsis male sexual lineage exhibits more robust maintenance of CG methylation than somatic tissues. Proc Natl Acad Sci U S A. 2016; 113(52):15132-7.

28. Rutowicz K, Puzio M, Halibart-Puzio J, Lirski M, Kotlinski M, Kroten MA, Knizewski L, Lange B, Muszewska A, Sniegowska-Swierk K, Koscielniak J, Iwanicka-Nowicka R, Buza K, Janowiak F, Zmuda K, Joesaar I, Laskowska-Kaszub K, Fogtman A, Kollist H, Zielenkiewicz P, Tiuryn J, Siedlecki P, Swiezewski S, Ginalski K, Koblowska M, Archacki R, Wilczynski B, Rapacz M, Jerzmanowski A. A specialized histone H1 variant is required for adaptive responses to complex abiotic stress and related DNA methylation in Arabidopsis. Plant Physiol. 2015;169(3):2080-101.

29. Wollmann H, Stroud H, Yelagandula R, Tarutani Y, Jiang D, Jing L, Jamge B, Takeuchi H, Holec S, Nie X, Kakutani T, Jacobsen SE, Berger F. The histone H3 variant H3.3 regulates gene body DNA methylation in *Arabidopsis thaliana*. Genome Biol. 2017;18(1):94.

30. Kotlinski M, Rutowicz K, Knizewski L, Palusinski A, Oledzki J, Fogtman A, Rubel T, Koblowska M, Dadlez M, Ginalski K, Jerzmanowski A. Histone H1 variants in Arabidopsis are subject to numerous post-translational modifications, both conserved and previously unknown in histones, suggesting complex functions of H1 in plants. PLoS One. 2016;11(1):e0147908.
31. Rea M, Zheng W, Chen M, Braud C, Bhangu D, Rognan TN, Xiao W. Histone H1 affects gene imprinting and DNA methylation in Arabidopsis. Plant J. 2012;71(5):776-86.

32. Zemach A, Kim MY, Hsieh PH, Coleman-Derr D, Eshed-Williams L, Thao K, Harmer SL, Zilberman D. The Arabidopsis nucleosome remodeler DDM1 allows DNA methyltransferases to access H1-containing heterochromatin. Cell. 2013;153(1):193-205.

33. Lyons DB, Zilberman D. DDM1 and Lsh remodelers allow methylation of DNA wrapped in nucleosomes. Elife. 2017;6:e30674.

34. Salazar-Henao JE, Velez-Bermudez IC, Schmidt W. The regulation and plasticity of root hair patterning and morphogenesis. Development. 2016; 143(11):1848-58.

35. Van Norman JM, Xuan W, Beeckman T, Benfey PN. To branch or not to branch: the role of pre-patterning in lateral root formation. Development. 2013;140(21):4301-10.

36. Ikeuchi M, Iwase A, Rymen B, Harashima H, Shibata M, Ohnuma M, Breuer C, Morao AK, de Lucas M, De Veylder L, Goodrich J, Brady SM, Roudier F, Sugimoto K. PRC2 represses dedifferentiation of mature somatic cells in Arabidopsis. Nat Plants. 2015;1:15089.

37. Lau OS, Bergmann DC. Stomatal development: a plant's perspective on cell polarity, cell fate transitions and intercellular communication. Development. 2012;139(20):3683-92.

38. He C, Chen X, Huang H, Xu L. Reprogramming of H3K27me3 is critical for acquisition of pluripotency from cultured Arabidopsis tissues. PLoS Genet. 2012;8(8):e1002911.

39. Fransz P, De Jong JH, Lysak M, Castiglione MR, Schubert I. Interphase chromosomes in Arabidopsis are organized as well defined chromocenters from which euchromatin loops emanate. Proc Natl Acad Sci U S A. 2002;99(22):14584-9.

40. Sequeira-Mendes J, Araguez I, Peiro R, Mendez-Giraldez R, Zhang X, Jacobsen SE, Bastolla U, Gutierrez C. The functional topography of the Arabidopsis genome is organized in a reduced number of linear motifs of chromatin states. Plant Cell. 2014;26(6):2351-66.

41. Soppe WJ, Jasencakova Z, Houben A, Kakutani T, Meister A, Huang MS, Jacobsen SE, Schubert I, Fransz PF. DNA methylation controls histone H3 lysine 9 methylation and heterochromatin assembly in Arabidopsis. EMBO J. 2002;21(23):6549-59.

42. Ashenafi MS, Baroux C. Automated 3D gene position analysis using a customized Imaris plugin: XTFISHInsideNucleus. Methods Mol Biol. 2018; 1675:591-613.

43. Andrey P, Kieu K, Kress C, Lehmann G, Tirichine L, Liu Z, Biot E, Adenot PG, Hue-Beauvais C, Houba-Herin N, Duranthon V, Devinoy E, Beaujean N, Gaudin V, Maurin Y, Debey P. Statistical analysis of 3D images detects regular spatial distributions of centromeres and chromocenters in animal and plant nuclei. PLoS Comput Biol. 2010; 6(7):e1000853.

44. Vergara Z, Gutierrez C. Emerging roles of chromatin in the maintenance of genome organization and function in plants. Genome Biol. 2017;18(1):96.

45. Cherkezian L, Stypula-Cyrus Y, Subramanian H, White C, Dela Cruz M, Wali RK, Goldberg MJ, Bianchi LK, Roy HK, Backman V. Nanoscale changes in chromatin organization represent the initial steps of tumorigenesis: a transmission electron microscopy study. BMC Cancer. 2014;14:189.

46. Fabrice T, Cherkezeyan NL, Ringli C, Baroux C. Transmission electron microscopy imaging to analyse chromatin density distribution at the nanoscale level. Methods Mol Biol. 2018;1675:633-51. https://doi.org/10.1 007/978-1-4939-7318-7_34

47. Munoz-Viana R, Wildhaber T, Trejo-Arellano MS, Mozgova I, Hennig L. Arabidopsis chromatin assembly factor 1 is required for occupancy and position of a subset of nucleosomes. Plant J. 2017;92(3):363-74.

48. Lucas M, Kenobi K, von Wangenheim D, Vobeta U, Swarup K, De Smet I, Van Damme D, Lawrence T, Peret B, Moscardi E, Barbeau D, Godin C, Salt D, Guyomarc'h S, Stelzer EH, Maizel A, Laplaze L, Bennett MJ. Lateral root morphogenesis is dependent on the mechanical properties of the overlaying tissues. Proc Natl Acad Sci U S A. 2013;110(13):5229-34.

49. Carrero G, Crawford E, Th'ng J, de Vries G, Hendzel MJ. Quantification of protein-protein and protein-DNA interactions in vivo, using fluorescence recovery after photobleaching. Methods Enzymol. 2004;375:415-42.

50. Rosa S, Ntoukakis V, Ohmido N, Pendle A, Abranches R, Shaw P. Cell differentiation and development in Arabidopsis are associated with changes in histone dynamics at the single-cell level. Plant Cell. 2014;26(12):4821-33.

51. Li G, Liu S, Wang J, He J, Huang H, Zhang Y, Xu L. ISWI proteins participate in the genome-wide nucleosome distribution in Arabidopsis. Plant J. 2014; 78(4):706-14.

52. Zhang X, Clarenz O, Cokus S, Bernatavichute YV, Pellegrini M, Goodrich J, Jacobsen SE. Whole-genome analysis of histone H3 lysine 27 trimethylation in Arabidopsis. PLoS Biol. 2007;5(5):e129.

53. Wang H, Liu C, Cheng J, Liu J, Zhang L, He C, Shen WH, Jin H, Xu L, Zhang Y. Arabidopsis flower and embryo developmental genes are repressed in seedlings by different combinations of Polycomb group proteins in association with distinct sets of cis-regulatory elements. PLoS Genet. 2016;12(1):e1005771.

54. Pass DA, Sornay E, Marchbank A, Crawford MR, Paszkiewicz K, Kent NA, Murray JAH. Genome-wide chromatin mapping with size resolution reveals a dynamic sub-nucleosomal landscape in Arabidopsis. PLoS Genet. 2017; 13(9):e1006988.

55. Jerzmanowski A. SWI/SNF chromatin remodeling and linker histones in plants. Biochim Biophys Acta. 2007;1769(5-6):330-45.

56. Prymakowska-Bosak M, Przewloka MR, Iwkiewicz J, Egierszdroff S, Kuras M, Chaubet N, Gigot C, Spiker S, Jerzmanowski A. Histone H1 overexpressed to high level in tobacco affects certain developmental programs but has limited effect on basal cellular functions. Proc Natl Acad Sci U S A. 1996; 93(19):10250-5.

57. Almassalha LM, Tiwari A, Ruhoff PT, Stypula-Cyrus Y, Cherkezian L, Matsuda H, Dela Cruz MA, Chandler JE, White C, Maneval C, Subramanian H, Szleifer I, Roy HK, Backman V. The global relationship between chromatin physical topology, fractal structure, and gene expression. Sci Rep. 2017;7:41061.

58. He S, Vickers M, Zhang J, Feng X. Natural depletion of H1 in sex cells causes DNA demethylation, heterochromatin decondensation and transposon activation. Elife. 2019;8:e42530.

59. Du J, Zhong X, Bernatavichute YV, Stroud H, Feng S, Caro E, Vashisht AA, Terragni J, Chin HG, Tu A, Hetzel J, Wohlschlegel JA, Pradhan S, Patel DJ, Jacobsen SE. Dual binding of chromomethylase domains to H3K9me2-containing nucleosomes directs DNA methylation in plants. Cell. 2012; 151(1):167-80.

60. Cao K, Lailler N, Zhang Y, Kumar A, Uppal K, Liu Z, Lee EK, Wu H, Medrzycki M, Pan C, Ho PY, Cooper GP Jr, Dong X, Bock C, Bouhassira EE, Fan Y. High-resolution mapping of h1 linker histone variants in embryonic stem cells. PLoS Genet. 2013;9(4):e1003417.

61. Iwasaki YW, Murano K, Ishizu H, Shibuya A, Iyoda Y, Siomi MC, Siomi H, Saito K. Piwi modulates chromatin accessibility by regulating multiple factors including histone H1 to repress transposons. Mol Cell. 2016;63(3):408-19.

62. Probst AV, Fagard M, Proux F, Mourrain P, Boutet S, Earley K, Lawrence RJ, Pikaard CS, Murfett J, Furner I, Vaucheret H, Mittelsten Scheid O. Arabidopsis histone deacetylase HDA6 is required for maintenance of transcriptional gene silencing and determines nuclear organization of rDNA repeats. Plant Cell. 2004;16(4):1021-34.

63. Jagannathan M, Cummings R, Yamashita YM. A conserved function for pericentromeric satellite DNA. Elife. 2018;7:e34122.

64. Beshnova DA, Cherstvy AG, Vainshtein Y, Teif VB. Regulation of the nucleosome repeat length in vivo by the DNA sequence, protein concentrations and long-range interactions. PLoS Comput Biol. 2014;10(7):e1003698.

65. Woodcock CL, Ghosh RP. Chromatin higher-order structure and dynamics. Cold Spring Harb Perspect Biol. 2010;2(5):a000596.

66. Usachenko SI, Gavin IM, Bavykin SG. Alterations in nucleosome core structure in linker histone-depleted chromatin. J Biol Chem. 1996;271(7):3831-6.

67. Boule JB, Mozziconacci J, Lavelle C. The polymorphisms of the chromatin fiber. J Phys Condens Matter. 2015;27(3):033101.

68. Correll SJ, Schubert MH, Grigoryev SA. Short nucleosome repeats impose rotational modulations on chromatin fibre folding. EMBO J. 2012;31(10): 2416-26.

69. Routh A, Sandin S, Rhodes D. Nucleosome repeat length and linker histone stoichiometry determine chromatin fiber structure. Proc Natl Acad Sci U S A. 2008;105(26):8872-7.

70. Golkaram M, Jang J, Hellander S, Kosik KS, Petzold LR. The role of chromatin density in cell population heterogeneity during stem cell differentiation. Sci Rep. 2017;7(1):13307.

71. Lempe J, Lachowiec J, Sullivan AM, Queitsch C. Molecular mechanisms of robustness in plants. Curr Opin Plant Biol. 2013;16(1):62-9.

72. She W, Baroux C. Chromatin dynamics during plant sexual reproduction. Front Plant Sci. 2014;5:354.

73. Tessadori F, Chupeau MC, Chupeau Y, Knip M, Germann S, van Driel R, Fransz P, Gaudin V. Large-scale dissociation and sequential reassembly of
pericentric heterochromatin in dedifferentiated Arabidopsis cells. J Cell Sci.
2007a;120(Pt 7):1200-8.

74. Tessadori F, Schulkes RK, van Driel R, Fransz P. Light-regulated large-scale reorganization of chromatin during the floral transition in Arabidopsis. Plant J. 2007b;50(5):848-57.

75. Lafos M, Kroll P, Hohenstatt ML, Thorpe FL, Clarenz O, Schubert D. Dynamic regulation of H3K27 trimethylation during Arabidopsis differentiation. PLoS Genet. 2011;7(4):e1002040.

76. Liu X, Zhou S, Wang W, Ye Y, Zhao Y, Xu Q, Zhou C, Tan F, Cheng S, Zhou DX. Regulation of histone methylation and reprogramming of gene expression in the rice inflorescence meristem. Plant Cell. 2015;27(5):1428-44.

77. Mozgova I, Kohler C, Hennig L. Keeping the gate closed: functions of the polycomb repressive complex PRC2 in development. Plant J. 2015;83(1):121-32.

78. You Y, Sawikowska A, Neumann M, Pose D, Capovilla G, Langenecker T, Neher RA, Krajewski P, Schmid M. Temporal dynamics of gene expression and histone marks at the Arabidopsis shoot meristem during flowering. Nat Commun. 2017;8:15120.

79. Zhang Y, Cooke M, Panjwani S, Cao K, Krauth B, Ho PY, Medrzycki M, Berhe DT, Pan C, McDevitt TC, Fan Y. Histone H1 depletion impairs embryonic stem cell differentiation. PLoS Genet. 2012;8(5):e1002691.

80. Bouyer D, Roudier F, Heese M, Andersen ED, Gey D, Nowack MK, Goodrich J, Renou JP, Grini PE, Colot V, Schnittger A. Polycomb repressive complex 2 controls the embryo-to-seedling phase transition. PLoS Genet. 2011;7(3):e1002014.

81. Gu X, Xu T, He Y. A histone H3 lysine-27 methyltransferase complex represses lateral root formation in Arabidopsis thaliana. Mol Plant. 2014;7(6):977-88.

82. Lee E, Lucas JR, Goodrich J, Sack FD. Arabidopsis guard cell integrity involves the epigenetic stabilization of the FLP and FAMA transcription factor genes. Plant J. 2014;78(4):566-77.

83. Molitor AM, Bu Z, Yu Y, Shen WH. Arabidopsis AL PHD-PRC1 complexes promote seed germination through H3K4me3-to-H3K27me3 chromatin state switch in repression of seed developmental genes. PLoS Genet. 2014; 10(1):e1004091.

84. Wood CC, Robertson M, Tanner G, Peacock WJ, Dennis ES, Helliwell CA. The Arabidopsis thaliana vernalization response requires a polycomb-like protein complex that also includes VERNALIZATION INSENSITIVE 3. Proc Natl Acad Sci U S A. 2006;103(39):14631-6.

85. Zhu D, Rosa S, Dean C. Nuclear organization changes and the epigenetic silencing of FLC during vernalization. J Mol Biol. 2015;427(3):659-69.

86. Martin C, Cao R, Zhang Y. Substrate preferences of the EZH2 histone methyltransferase complex. J Biol Chem. 2006;281(13):8365-70.

87. Jerzmanowski A, Przewloka M, Grasser KD. Linker histones and HMG1 proteins of higher plants. Plant Biol. 2000;2(6):586-97.

88. Hellauer K, Sirard E, Turcotte B. Decreased expression of specific genes in yeast cells lacking histone H1. J Biol Chem. 2001;276(17):13587-92.

89. Sancho M, Diani E, Beato M, Jordan A. Depletion of human histone H1 variants uncovers specific roles in gene expression and cell growth. PLoS Genet. 2008;4(10):e1000227.

90. Shen X, Gorovsky MA. Linker histone H1 regulates specific gene expression but not global transcription in vivo. Cell. 1996;86(3):475-83.

91. Clough SJ, Bent AF. Floral dip: a simplified method for Agrobacterium-mediated transformation of Arabidopsis thaliana. Plant J. 1998;16(6):735-43.

92. Pavlova P, Tessadori F, de Jong HJ, Fransz P. Immunocytological analysis of chromatin in isolated nuclei. Methods Mol Biol. 2010;655:413-32.

93. Schindelin J, Arganda-Carreras I, Frise E, Kaynig V, Longair M, Pietzsch T, Preibisch S, Rueden C, Saalfeld S, Schmid B, Tinevez JY, White DJ, Hartenstein V, Eliceiri K, Tomancak P, Cardona A. Fiji: an open-source platform for biological-image analysis. Nat Methods. 2012;9(7):676-82.

94. Bourbousse C, Mestiri I, Zabulon G, Bourge M, Formiggini F, Koini MA, Brown SC, Fransz P, Bowler C, Barneche F. Light signaling controls nuclear architecture reorganization during seedling establishment. Proc Natl Acad Sci U S A. 2015;112(21):E2836-44.

95. Chodavarapu RK, Feng S, Bernatavichute YV, Chen PY, Stroud H, Yu Y, Hetzel JA, Kuo F, Kim J, Cokus SJ, Casero D, Bernal M, Huijser P, Clark AT, Kramer U, Merchant SS, Zhang X, Jacobsen SE, Pellegrini M. Relationship between nucleosome positioning and DNA methylation. Nature. 2010;466(7304):388-92.

96. Maizel A, von Wangenheim D, Federici F, Haseloff J, Stelzer EH. High-resolution live imaging of plant growth in near physiological bright conditions using light sheet fluorescence microscopy. Plant J. 2011;68(2):377-85.

97. Phair RD, Gorski SA, Misteli T. Measurement of dynamic protein binding to chromatin in vivo, using photobleaching microscopy. Methods Enzymol. 2004;375:393-414.

98. Malamy JE, Benfey PN. Organization and cell differentiation in lateral roots of Arabidopsis thaliana. Development. 1997;124(1):33-44.

99. Kutter C, Schob H, Stadler M, Meins F Jr, Si-Ammour A. MicroRNA-mediated regulation of stomatal development in Arabidopsis. Plant Cell. 2007;19(8):2417-29.

100. Nakabayashi K, Bartsch M, Xiang Y, Miatton E, Pellengahr S, Yano R, Seo M, Soppe WJ. The time required for dormancy release in Arabidopsis is determined by DELAY OF GERMINATION1 protein levels in freshly harvested seeds. Plant Cell. 2012;24(7):2826-38.

101. Bowler C, Benvenuto G, Laflamme P, Molino D, Probst AV, Tariq M, Paszkowski J. Chromatin techniques for plant cells. Plant J. 2004;39(5):776-89.

102. Chomczynski P, Sacchi N. Single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction. Anal Biochem. 1987; 162(1):156-9.

103. Bushnell B. BBMap: a fast, accurate, splice-aware aligner; 2014.

104. Robinson MD, Oshlack A. A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol. 2010;11(3):R25.

105. Bullard JH, Purdom E, Hansen KD, Dudoit S. Evaluation of statistical methods for normalization and differential expression in mRNA-Seq experiments. BMC Bioinformatics. 2010;11:94.

106. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, Smyth GK. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.

107. Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5(7):621-8.

108. Kahles A, Behr J, Ratsch G. MMR: a tool for read multi-mapper resolution. Bioinformatics. 2016;32(5):770-2.

109. Chen K, Xi Y, Pan X, Li Z, Kaestner K, Tyler J, Dent S, He X, Li W. DANPOS: dynamic analysis of nucleosome position and occupancy by sequencing. Genome Res. 2013;23(2):341-51.

110. Kent WJ, Zweig AS, Barber G, Hinrichs AS, Karolchik D. BigWig and BigBed: enabling browsing of large distributed datasets. Bioinformatics. 2010;26(17):2204-7.

111. Ramirez F, Ryan DP, Gruning B, Bhardwaj V, Kilpert F, Richter AS, Heyne S, Dundar F, Manke T. deepTools2: a next generation web server for deep-sequencing data analysis. Nucleic Acids Res. 2016;44(W1):W160-5.

112. Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26(6):841-2.

113. Edgar R, Domrachev M, Lash AE. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002; 30(1):207-10.

114. Rutowicz K, Lirski M, Merma B, Teano G, Schubert J, Mestiri I, Kroteń MA, Fabrice TN, Fritz S, Grob S, et al., Linker histones are fine-scale chromatin architects modulating developmental decisions in Arabidopsis data sets. NCBI. 2019. https://www.ncbi.nlm.nih.gov/geo/ query/acc.cgi?acc=GSE113558.

115. Rutowicz K, Lirski M, Merma B, Teano G, Schubert J, Mestiri I, Kroteń MA, Fabrice TN, Fritz S, Grob S, et al. Linker histones are fine-scale chromatin architects modulating developmental decisions in Arabidopsis data sets. Figshare. 2019. https://figshare.com/articles/3h1 vs_WT_nuclei/8223977.

Publisher's Note

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

Ready to submit your research? Choose BMC and benefit from:

• fast, convenient online submission

* thorough peer review by experienced researchers in your field

* rapid publication on acceptance

• support for research data, including large and complex data types

* gold Open Access which fosters wider collaboration and increased citations

• maximum visibility for your research: over 100M website views per year

At BMC, research is always in progress.

Learn more biomedcentral.com/submissions

<figure><img src="image_16.png" /></figure>