Nucleosome spacing and chromatin higher-order folding

Genome wide nucleosome landscape shapes 3D chromatin organization

The hierarchical chromatin organization begins with formation of nucleosomes, which fold into chromatin domains punctuated by boundaries and ultimately chromosomes. In a hierarchal organization, lower levels shape higher levels. However, the dependence of higher-order 3D chromatin organization on the nucleosome-level organization has not been studied in cells. We investigated the relationship between nucleosome-level organization and higher-order chromatin organization by perturbing nucleosomes across the genome by deleting Imitation SWItch (ISWI) and Chromodomain Helicase DNA-binding (CHD1) chromatin remodeling factors in budding yeast. We find that changes in nucleosome-level properties are accompanied by changes in 3D chromatin organization. Short-range chromatin contacts up to a few kilo-base pairs decrease, chromatin domains weaken, and boundary strength decreases. Boundary strength scales with accessibility and moderately with width of nucleosome-depleted region. Change in nucleosome positioning seems to alter the stiffness of chromatin, which can affect formation of chromatin contacts. Our results suggest a biomechanical "bottom-up" mechanism by which nucleosome distribution across genome shapes 3D chromatin organization.

Cryoelectron tomography reveals the multiplex anatomy of condensed native chromatin and its unfolding by histone citrullination

Nucleosome chains fold and self-associate to form higher-order structures whose internal organization is unknown. Here, cryoelectron tomography (cryo-ET) of native human chromatin reveals intrinsic folding motifs such as (1) non-uniform nucleosome stacking, (2) intermittent parallel and perpendicular orientations of adjacent nucleosome planes, and (3) a regressive nucleosome chain path, which deviates from the direct zigzag topology seen in reconstituted nucleosomal arrays. By examining the self-associated structures, we observed prominent nucleosome stacking in cis and anti-parallel nucleosome interactions, which are consistent with partial nucleosome interdigitation in trans. Histone citrullination strongly inhibits nucleosome stacking and self-association with a modest effect on chromatin folding, whereas the reconstituted arrays undergo a dramatic unfolding into open zigzag chains induced by histone citrullination. This study sheds light on the internal structure of compact chromatin nanoparticles and suggests a mechanism for how epigenetic changes in chromatin folding are retained across both open and condensed forms.

Human macroH2A1 drives nucleosome dephasing and genome instability in histone-humanized yeast

In addition to replicative histones, eukaryotic genomes encode a repertoire of non-replicative variant histones providing additional layers of structural and epigenetic regulation. Here, we systematically replaced individual replicative human histones with non-replicative human variant histones using a histone replacement system in yeast. Variants H2A.J, TsH2B, and H3.5 complemented for their respective replicative counterparts. However, macroH2A1 failed to complement and its expression was toxic in yeast, negatively interacting with native yeast histones and kinetochore genes. To isolate yeast with "macroH2A1 chromatin" we decoupled the effects of its macro and histone fold domains, which revealed that both domains sufficed to override native yeast nucleosome positioning. Furthermore, both modified constructs of macroH2A1 exhibited lower nucleosome occupancy that correlated with decreased short-range chromatin interactions (<20 Kb), disrupted centromeric clustering, and increased chromosome instability. While supporting viability, macroH2A1 dramatically alters chromatin organization in yeast, leading to genome instability and massive fitness defects.

Electronic Mapping of a Bacterial Genome with Dual Solid-State Nanopores and Active Single-Molecule Control

We present an electronic mapping of a bacterial genome using solid-state nanopore technology. A dual-nanopore architecture and active control logic are used to produce single-molecule data that enables estimation of distances between physical tags installed at sequence motifs within double-stranded DNA. Previously developed "DNA flossing" control logic generates multiple scans of each captured DNA. We extended this logic in two ways: first, to automate "zooming out" on each molecule to progressively increase the number of tags scanned during flossing, and second, to automate recapture of a molecule that exited flossing to enable interrogation of the same and/or different regions of the molecule. Custom analysis methods were developed to produce consensus alignments from each multiscan event. The combined multiscanning and multicapture method was applied to the challenge of mapping from a heterogeneous mixture of single-molecule fragments that make up the Escherichia coli (E. coli) chromosome. Coverage of 3.1× across 2355 resolvable sites of the E. coli genome was achieved after 5.6 h of recording time. The recapture method showed a 38% increase in the merged-event alignment length compared to single-scan alignments. The observed intertag resolution was 150 bp in engineered DNA molecules and 166 bp natively within fragments of E. coli DNA, with detection of 133 intersite intervals shorter than 200 bp in the E. coli reference map. We present results on estimating distances in repetitive regions of the E. coli genome. With an appropriately designed array, higher throughput implementations could enable human-sized genome and epigenome mapping applications.

Slow nucleosome dynamics set the transcriptional speed limit and induce RNA polymerase II traffic jams and bursts

Nucleosomes are recognized as key regulators of transcription. However, the relationship between slow nucleosome unwrapping dynamics and bulk transcriptional properties has not been thoroughly explored. Here, an agent-based model that we call the dynamic defect Totally Asymmetric Simple Exclusion Process (ddTASEP) was constructed to investigate the effects of nucleosome-induced pausing on transcriptional dynamics. Pausing due to slow nucleosome dynamics induced RNAPII convoy formation, which would cooperatively prevent nucleosome rebinding leading to bursts of transcription. The mean first passage time (MFPT) and the variance of first passage time (VFPT) were analytically expressed in terms of the nucleosome rate constants, allowing for the direct quantification of the effects of nucleosome-induced pausing on pioneering polymerase dynamics. The mean first passage elongation rate γ(h_c, h_o) is inversely proportional to the MFPT and can be considered to be a new axis of the ddTASEP phase diagram, orthogonal to the classical αβ-plane (where α and β are the initiation and termination rates). Subsequently, we showed that, for β = 1, there is a novel jamming transition in the αγ-plane that separates the ddTASEP dynamics into initiation-limited and nucleosome pausing-limited regions. We propose analytical estimates for the RNAPII density ρ, average elongation rate v, and transcription flux J and verified them numerically. We demonstrate that the intra-burst RNAPII waiting times t_in follow the time-headway distribution of a max flux TASEP and that the average inter-burst interval tIBI¯ correlates with the index of dispersion D_e. In the limit γ→0, the average burst size reaches a maximum set by the closing rate h_c. When α≪1, the burst sizes are geometrically distributed, allowing large bursts even while the average burst size NB¯ is small. Last, preliminary results on the relative effects of static and dynamic defects are presented to show that dynamic defects can induce equal or greater pausing than static bottle necks.

To perform specific functions, cells must express specific genes by copying the information in DNA into RNA via transcription. Structural proteins called nucleosomes are spaced every 200 base pairs along the length of a strand of DNA and play a crucial function in the regulation of gene activity by tightly binding DNA strands and condensing them into heterochromatin, preventing transcription by RNA polymerase II (RNAPII). Even on active genes where nucleosomes are loosely attached to DNA strands, the wrapping and unwrapping of nucleosomes pause transcription as RNAPII passes by. Previous mathematical models of transcription have compared this biological process to traffic on a one lane highway without obstructions. In contrast, our proposed model simulates transcription like traffic in a grid system where nucleosomes can be thought of as pedestrians or other vehicles crossing the road at regularly spaced intersections. Just as side street traffic and pedestrian crossings can cause cars to form convoys and cause jams limiting the max speed in an area, nucleosomes can cause RNAPII to form convoys that lead to bursts of mRNA production and limit the average polymerase flux through the gene.

Histone Tail Conformations: A Fuzzy Affair with DNA

The core histone tails are critical in chromatin structure and signaling. Studies over the past several decades have provided a wealth of information on the histone tails and their interaction with chromatin factors. However, the conformation of the histone tails in a chromatin relevant context has remained elusive. Only recently has enough evidence emerged to start to build a structural model of the tails in the context of nucleosomes and nucleosome arrays. Here, we review these studies and propose that the histone tails adopt a high-affinity fuzzy complex with DNA, characterized by robust but dynamic association. Furthermore, we discuss how these DNA-bound conformational ensembles promote distinct chromatin structure and signaling, and that their fuzzy nature is important in transitioning between functional states.

On the role of transcription in positioning nucleosomes

Nucleosome positioning is crucial for the genome’s function. Though the role of DNA sequence in positioning nucleosomes is well understood, a detailed mechanistic understanding on the impact of transcription remains lacking. Using numerical simulations, we investigated the dependence of nucleosome density profiles on transcription level across multiple species. We found that the low nucleosome affinity of yeast, but not mouse, promoters contributes to the formation of phased nucleosomes arrays for inactive genes. For the active genes, a heterogeneous distribution of +1 nucleosomes, caused by a tug-of-war between two types of remodeling enzymes, is essential for reproducing their density profiles. In particular, while positioning enzymes are known to remodel the +1 nucleosome and align it toward the transcription start site (TSS), spacer enzymes that use a pair of nucleosomes as their substrate can shift the nucleosome array away from the TSS. Competition between these enzymes results in two types of nucleosome density profiles with well- and ill-positioned +1 nucleosome. Finally, we showed that Pol II assisted histone exchange, if occurring at a fast speed, can abolish the impact of remodeling enzymes. By elucidating the role of individual factors, our study reconciles the seemingly conflicting results on the overall impact of transcription in positioning nucleosomes across species.

Nucleosome positioning plays a key role in the genome’s function by regulating the accessibility of protein binding sites as well as higher-order chromatin organization. Though significant progress has been made towards studying the role of DNA sequence in positioning the nucleosomes, our understanding on the impact of transcription lags behind. Our study uses kinetic simulations to explore the role of DNA sequence specificity, transcription factor binding, enzyme remodeling, and Pol II elongation in positioning nucleosomes. It suggests that the differences in nucleosome density profiles observed at various transcription levels in yeast and mouse embryonic stem cells can be understood from a tug-of-war between two types of remodeling enzymes.

Elucidating the influence of linker histone variants on chromatosome dynamics and energetics

Linker histones are epigenetic regulators that bind to nucleosomes and alter chromatin structures and dynamics. Biophysical studies have revealed two binding modes in the linker histone/nucleosome complex, the chromatosome, where the linker histone is either centered on or askew from the dyad axis. Each has been posited to have distinct effects on chromatin, however the molecular and thermodynamic mechanisms that drive them and their dependence on linker histone compositions remain poorly understood. We present molecular dynamics simulations of chromatosomes with the globular domain of two linker histone variants, generic H1 (genGH1) and H1.0 (GH1.0), to determine how their differences influence chromatosome structures, energetics and dynamics. Results show that both unbound linker histones adopt a single compact conformation. Upon binding, DNA flexibility is reduced, resulting in increased chromatosome compaction. While both variants enthalpically favor on-dyad binding, energetic benefits are significantly higher for GH1.0, suggesting that GH1.0 is more capable than genGH1 of overcoming the large entropic reduction required for on-dyad binding which helps rationalize experiments that have consistently demonstrated GH1.0 in on-dyad states but that show genGH1 in both locations. These simulations highlight the thermodynamic basis for different linker histone binding motifs, and details their physical and chemical effects on chromatosomes.

Rouse model with transient intramolecular contacts on a timescale of seconds recapitulates folding and fluctuation of yeast chromosomes

DNA folding and dynamics along with major nuclear functions are determined by chromosome structural properties, which remain, thus far, elusive in vivo. Here, we combine polymer modeling and single particle tracking experiments to determine the physico-chemical parameters of chromatin in vitro and in living yeast. We find that the motion of reconstituted chromatin fibers can be recapitulated by the Rouse model using mechanical parameters of nucleosome arrays deduced from structural simulations. Conversely, we report that the Rouse model shows some inconsistencies to analyze the motion and structural properties inferred from yeast chromosomes determined with chromosome conformation capture techniques (specifically, Hi-C). We hence introduce the Rouse model with Transient Internal Contacts (RouseTIC), in which random association and dissociation occurs along the chromosome contour. The parametrization of this model by fitting motion and Hi-C data allows us to measure the kinetic parameters of the contact formation reaction. Chromosome contacts appear to be transient; associated to a lifetime of seconds and characterized by an attractive energy of -0.3 to -0.5 kBT. We suggest attributing this energy to the occurrence of histone tail-DNA contacts and notice that its amplitude sets chromosomes in 'theta' conditions, in which they are poised for compartmentalization and phase separation.

Nucleosome positioning and spacing: from genome-wide maps to single arrays

The positioning of nucleosomes relative to DNA and their neighboring nucleosomes represents a fundamental layer of chromatin organization. Changes in nucleosome positioning and spacing affect the accessibility of DNA to regulatory factors and the formation of higher order chromatin structures. Sequencing of mononucleosomal fragments allowed mapping nucleosome positions on a genome-wide level in many organisms. This revealed that successions of evenly spaced and well-positioned nucleosomes-so called phased nucleosome arrays-occur at the 5' end of many active genes and in the vicinity of transcription factor and other protein binding sites. Phased arrays arise from the interplay of barrier elements on the DNA, which position adjacent nucleosomes, and the nucleosome spacing activity of ATP-dependent chromatin remodelers. A shortcoming of classic mononucleosomal mapping experiments is that they only reveal nucleosome spacing and array regularity at select sites in the genome with well-positioned nucleosomes. However, new technological approaches elucidate nucleosome array structure throughout the genome and with single-cell resolution. In the future, it will be interesting to see whether changes in nucleosome array regularity and spacing contribute to the formation of higher order chromatin structures and the spatial organization of the genome in vivo.

Chromatin fiber structural motifs as regulatory hubs of genome function?

Nucleosomes cover eukaryotic genomes like beads on a string and play a central role in regulating genome function. Isolated strings of nucleosomes have the potential to compact and form higher order chromatin structures, such as the well-characterized 30-nm fiber. However, despite tremendous advances in observing chromatin fibers in situ it has not been possible to confirm that regularly ordered fibers represent a prevalent structural level in the folding of chromosomes. Instead, it appears that folding at a larger scale than the nucleosome involves a variety of random structures with fractal characteristics. Nevertheless, recent progress provides evidence for the existence of structural motifs in chromatin fibers, potentially localized to strategic sites in the genome. Here we review the current understanding of chromatin fiber folding and the emerging roles that oligonucleosomal motifs play in the regulation of genome function.

Probing Chromatin Structure with Magnetic Tweezers

Magnetic tweezers form a unique tool to study the topology and mechanical properties of chromatin fibers. Chromatin is a complex of DNA and proteins that folds the DNA in such a way that meter-long stretches of DNA fit into the micron-sized cell nucleus. Moreover, it regulates accessibility of the genome to the cellular replication, transcription, and repair machinery. However, the structure and mechanisms that govern chromatin folding remain poorly understood, despite recent spectacular improvements in high-resolution imaging techniques. Single-molecule force spectroscopy techniques can directly measure both the extension of individual chromatin fragments with nanometer accuracy and the forces involved in the (un)folding of single chromatin fibers. Here, we report detailed methods that allow one to successfully prepare in vitro reconstituted chromatin fibers for use in magnetic tweezers-based force spectroscopy. The higher-order structure of different chromatin fibers can be inferred from fitting a statistical mechanics model to the force-extension data. These methods for quantifying chromatin folding can be extended to study many other processes involving chromatin, such as the epigenetic regulation of transcription.

DNA Methylation Landscape Reflects the Spatial Organization of Chromatin in Different Cells

The relationship between DNA methylation and chromatin structure is still largely unknown. By analyzing a large set of published sequencing data, we observed a long-range power law correlation of DNA methylation with cell class-specific scaling exponents in the range of tens of kilobases. We showed that such cell class-specific scaling exponents are caused by different patchiness of DNA methylation in different cells. By modeling the chromatin structure using high-resolution chromosome conformation capture data and mapping the methylation level onto the modeled structure, we demonstrated that the patchiness of DNA methylation is related to chromatin structure. The scaling exponents of the power law correlation are thus a display of the spatial organization of chromatin. Besides the long-range correlation, we also showed that the local correlation of DNA methylation is associated with nucleosome positioning. The local correlation of partially methylated domains is different from that of nonpartially methylated domains, suggesting that their chromatin structures differ at the scale of several hundred base pairs (covering a few nucleosomes). Our study provides a novel, to our knowledge, view of the spatial organization of chromatin structure from a perspective of DNA methylation, in which both long-range and local correlations of DNA methylation along the genome reflect the spatial organization of chromatin.

Chromatin remodeler CHD1 promotes XPC-to-TFIIH handover of nucleosomal UV lesions in nucleotide excision repair

Ultraviolet (UV) light induces mutagenic cyclobutane pyrimidine dimers (CPDs) in nucleosomal DNA that is tightly wrapped around histone octamers. How global-genome nucleotide excision repair (GG-NER) processes CPDs despite that this chromatin arrangement is poorly understood. An increased chromatin association of CHD1 (chromodomain helicase DNA-binding 1) upon UV irradiation indicated possible roles of this chromatin remodeler in the UV damage response. Immunoprecipitation of chromatin fragments revealed that CHD1 co-localizes in part with GG-NER factors. Chromatin fractionation showed that the UV-dependent recruitment of CHD1 occurs to UV lesions in histone-assembled nucleosomal DNA and that this CHD1 relocation requires the lesion sensor XPC (xeroderma pigmentosum group C). In situ immunofluorescence analyses further demonstrate that CHD1 facilitates substrate handover from XPC to the downstream TFIIH (transcription factor IIH). Consequently, CHD1 depletion slows down CPD excision and sensitizes cells to UV-induced cytotoxicity. The finding of a CHD1-driven lesion handover between sequentially acting GG-NER factors on nucleosomal histone octamers suggests that chromatin provides a recognition scaffold enabling the detection of a subset of CPDs.

Single-molecule force spectroscopy on histone H4 tail-cross-linked chromatin reveals fiber folding

The eukaryotic genome is highly compacted into a protein-DNA complex called chromatin. The cell controls access of transcriptional regulators to chromosomal DNA via several mechanisms that act on chromatin-associated proteins and provide a rich spectrum of epigenetic regulation. Elucidating the mechanisms that fold chromatin fibers into higher-order structures is therefore key to understanding the epigenetic regulation of DNA accessibility. Here, using histone H4-V21C and histone H2A-E64C mutations, we employed single-molecule force spectroscopy to measure the unfolding of individual chromatin fibers that are reversibly cross-linked through the histone H4 tail. Fibers with covalently linked nucleosomes featured the same folding characteristics as fibers containing wild-type histones but exhibited increased stability against stretching forces. By stabilizing the secondary structure of chromatin, we confirmed a nucleosome repeat length (NRL)-dependent folding. Consistent with previous crystallographic and cryo-EM studies, the obtained force-extension curves on arrays with 167-bp NRLs best supported an underlying structure consisting of zig-zag, two-start fibers. For arrays with 197-bp NRLs, we previously inferred solenoidal folding, which was further corroborated by force-extension curves of the cross-linked fibers. The different unfolding pathways exhibited by these two types of arrays and reported here extend our understanding of chromatin structure and its potential roles in gene regulation. Importantly, these findings imply that chromatin compaction by nucleosome stacking protects nucleosomal DNA from external forces up to 4 piconewtons.

Epigenetic characteristics of the mitotic chromosome in 1D and 3D

While chromatin characteristics in interphase are widely studied, characteristics of mitotic chromatin and their inheritance through mitosis are still poorly understood. During mitosis, chromatin undergoes dramatic changes: transcription stalls, chromatin-binding factors leave the chromatin, histone modifications change and chromatin becomes highly condensed. Many key insights into mitotic chromosome state and conformation have come from extensive microscopy studies over the last century. Over the last decade, the development of 3C-based techniques has enabled the study of higher order chromosome organization during mitosis in a genome-wide manner. During mitosis, chromosomes lose their cell type-specific and locus-dependent chromatin organization that characterizes interphase chromatin and fold into randomly positioned loop arrays. Upon exit of mitosis, cells are capable of quickly rearranging the chromosome conformation to form the cell type-specific interphase organization again. The information that enables this rearrangement after mitotic exit is thought to be encoded at least in part in mitotic bookmarks, e.g. histone modifications and variants, histone remodelers, chromatin factors, and non-coding RNA. Here we give an overview of the chromosomal organization and epigenetic characteristics of interphase and mitotic chromatin in vertebrates. Second, we describe different ways in which mitotic bookmarking enables epigenetic memory of the features of interphase chromatin through mitosis. And third, we explore the role of epigenetic modifications and mitotic bookmarking in cell differentiation.

Regulation of nucleosome positioning by a CHD Type III chromatin remodeler and its relationship to developmental gene expression in Dictyostelium

Nucleosome placement and repositioning can direct transcription of individual genes; however, the precise interactions of these events are complex and largely unresolved at the whole-genome level. The Chromodomain-Helicase-DNA binding (CHD) Type III proteins are a subfamily of SWI2/SNF2 proteins that control nucleosome positioning and are associated with several complex human disorders, including CHARGE syndrome and autism. Type III CHDs are required for multicellular development of animals and Dictyostelium but are absent in plants and yeast. These CHDs can mediate nucleosome translocation in vitro, but their in vivo mechanism is unknown. Here, we use genome-wide analysis of nucleosome positioning and transcription profiling to investigate the in vivo relationship between nucleosome positioning and gene expression during development of wild-type (WT) Dictyostelium and mutant cells lacking ChdC, a Type III CHD protein ortholog. We demonstrate major nucleosome positional changes associated with developmental gene regulation in WT. Loss of chdC caused an increase of intragenic nucleosome spacing and misregulation of gene expression, affecting ∼50% of the genes that are repositioned during WT development. These analyses demonstrate active nucleosome repositioning during Dictyostelium multicellular development, establish an in vivo function of CHD Type III chromatin remodeling proteins in this process, and reveal the detailed relationship between nucleosome positioning and gene regulation, as cells transition between developmental states.

Correlation among DNA Linker Length, Linker Histone Concentration, and Histone Tails in Chromatin

Eukaryotic cells condense their genetic material in the nucleus in the form of chromatin, a macromolecular complex made of DNA and multiple proteins. The structure of chromatin is intimately connected to the regulation of all eukaryotic organisms, from amoebas to humans, but its organization remains largely unknown. The nucleosome repeat length (NRL) and the concentration of linker histones (ρLH) are two structural parameters that vary among cell types and cell cycles; the NRL is the number of DNA basepairs wound around each nucleosome core plus the number of basepairs linking successive nucleosomes. Recent studies have found a linear empirical relationship between the variation of these two properties for different cells, but its underlying mechanism remains elusive. Here we apply our established mesoscale chromatin model to explore the mechanisms responsible for this relationship, by investigating chromatin fibers as a function of NRL and ρLH combinations. We find that a threshold of linker histone concentration triggers the compaction of chromatin into well-formed 30-nm fibers; this critical value increases linearly with NRL, except for long NRLs, where the fibers remain disorganized. Remarkably, the interaction patterns between core histone tails and chromatin elements are highly sensitive to the NRL and ρLH combination, suggesting a molecular mechanism that could have a key role in regulating the structural state of the fibers in the cell. An estimate of the minimized work and volume associated with storage of chromatin fibers in the nucleus further suggests factors that could spontaneously regulate the NRL as a function of linker histone concentration. Both the tail interaction map and DNA packing considerations support the empirical NRL/ρLH relationship and offer a framework to interpret experiments for different chromatin conditions in the cell.

Nucleosome spacing generated by ISWI and CHD1 remodelers is constant regardless of nucleosome density

Arrays of regularly spaced nucleosomes are a hallmark of chromatin, but it remains unclear how they are generated. Recent genome-wide studies, in vitro and in vivo, showed constant nucleosome spacing even if the histone concentration was experimentally reduced. This counters the long-held assumption that nucleosome density determines spacing and calls for factors keeping spacing constant regardless of nucleosome density. We call this a clamping activity. Here, we show in a purified system that ISWI- and CHD1-type nucleosome remodelers have a clamping activity such that they not only generate regularly spaced nucleosome arrays but also generate constant spacing regardless of nucleosome density. This points to a functionally attractive nucleosome interaction that could be mediated either directly by nucleosome-nucleosome contacts or indirectly through the remodelers. Mutant Drosophila melanogaster ISWI without the Hand-Sant-Slide (HSS) domain had no detectable spacing activity even though it is known to remodel and slide nucleosomes. This suggests that the role of ISWI remodelers in generating constant spacing is not just to mediate nucleosome sliding; they actively contribute to the attractive interaction. Additional factors are necessary to set physiological spacing in absolute terms.

Forced unraveling of chromatin fibers with nonuniform linker DNA lengths

The chromatin fiber undergoes significant structural changes during the cell's life cycle to modulate DNA accessibility. Detailed mechanisms of such structural transformations of chromatin fibers as affected by various internal and external conditions such as the ionic conditions of the medium, the linker DNA length, and the presence of linker histones, constitute an open challenge. Here we utilize Monte Carlo (MC) simulations of a coarse grained model of chromatin with nonuniform linker DNA lengths as found in vivo to help explain some aspects of this challenge. We investigate the unfolding mechanisms of chromatin fibers with alternating linker lengths of 26-62 bp and 44-79 bp using a series of end-to-end stretching trajectories with and without linker histones and compare results to uniform-linker-length fibers. We find that linker histones increase overall resistance of nonuniform fibers and lead to fiber unfolding with superbeads-on-a-string cluster transitions. Chromatin fibers with nonuniform linker DNA lengths display a more complex, multi-step yet smoother process of unfolding compared to their uniform counterparts, likely due to the existence of a more continuous range of nucleosome-nucleosome interactions. This finding echoes the theme that some heterogeneity in fiber component is biologically advantageous.

Evolution of histone 2A for chromatin compaction in eukaryotes

During eukaryotic evolution, genome size has increased disproportionately to nuclear volume, necessitating greater degrees of chromatin compaction in higher eukaryotes, which have evolved several mechanisms for genome compaction. However, it is unknown whether histones themselves have evolved to regulate chromatin compaction. Analysis of histone sequences from 160 eukaryotes revealed that the H2A N-terminus has systematically acquired arginines as genomes expanded. Insertion of arginines into their evolutionarily conserved position in H2A of a small-genome organism increased linear compaction by as much as 40%, while their absence markedly diminished compaction in cells with large genomes. This effect was recapitulated in vitro with nucleosomal arrays using unmodified histones, indicating that the H2A N-terminus directly modulates the chromatin fiber likely through intra- and inter-nucleosomal arginine-DNA contacts to enable tighter nucleosomal packing. Our findings reveal a novel evolutionary mechanism for regulation of chromatin compaction and may explain the frequent mutations of the H2A N-terminus in cancer.