Why Do Nucleosomes Unwrap Asymmetrically?

Chromatin remodelers: a concise introduction for biophysicists

Chromatin remodelers are molecular motors that act on nucleosomes: they move them along DNA or (dis-)assemble them. Despite the fact that they perform essential regulatory functions in cells-their deregulation can contribute to the development of cancers and lead to cell death-chromatin remodelers have only received meager attention in the biophysics community so far. In this short text, we attempt to present the key features of this interesting class of enzymes obtained with different experimental and theoretical methods, thereby providing a concise introduction for biophysicists to further stimulate interest in their properties.

Real-Time Multistep Asymmetrical Disassembly of Nucleosomes and Chromatosomes Visualized by High-Speed Atomic Force Microscopy

During replication, expression, and repair of the eukaryotic genome, cellular machinery must access the DNA wrapped around histone proteins forming nucleosomes. These octameric protein·DNA complexes are modular, dynamic, and flexible and unwrap or disassemble either spontaneously or by the action of molecular motors. Thus, the mechanism of formation and regulation of subnucleosomal intermediates has gained attention genome-wide because it controls DNA accessibility. Here, we imaged nucleosomes and their more compacted structure with the linker histone H1 (chromatosomes) using high-speed atomic force microscopy to visualize simultaneously the changes in the DNA and the histone core during their disassembly when deposited on mica. Furthermore, we trained a neural network and developed an automatic algorithm to track molecular structural changes in real time. Our results show that nucleosome disassembly is a sequential process involving asymmetrical stepwise dimer ejection events. The presence of H1 restricts DNA unwrapping, significantly increases the nucleosomal lifetime, and affects the pathway in which heterodimer asymmetrical dissociation occurs. We observe that tetrasomes are resilient to disassembly and that the tetramer core (H3·H4)2 can diffuse along the nucleosome positioning sequence. Tetrasome mobility might be critical to the proper assembly of nucleosomes and can be relevant during nucleosomal transcription, as tetrasomes survive RNA polymerase passage. These findings are relevant to understanding nucleosome intrinsic dynamics and their modification by DNA-processing enzymes.

Asymmetrical nucleosomal DNA signatures regulate transcriptional directionality

Despite the symmetrical structure of nucleosomes, in vitro studies have shown that transcription proceeds with different efficiency depending on the orientation of the DNA sequence around them. However, it is unclear whether this functional asymmetry is present in vivo and whether it could regulate transcriptional directionality. Here, we report that the proximal and distal halves of nucleosomal DNA contribute differentially to nucleosome stability in the genome. In +1 nucleosomes, this asymmetry facilitates or hinders transcription depending on the orientation of its underlying DNA, and this difference is associated with an asymmetrical interaction between DNA and histones. These properties are encoded in the DNA signature of +1 nucleosomes, since its incorporation in the two orientations into downstream nucleosomes renders them asymmetrically accessible to MNase and inverts the balance between sense and antisense transcription. Altogether, our results show that nucleosomal DNA endows nucleosomes with asymmetrical properties that modulate the directionality of transcription.

Why Are Nucleosome Breathing Dynamics Asymmetric?

In eukaryotic cells, DNA is bound to nucleosomes, but DNA segments occasionally unbind in the process known as nucleosome breathing. Although DNA can unwrap simultaneously from both ends of the nucleosome (symmetric breathing), experiments indicate that DNA prefers to dissociate from only one end (asymmetric breathing). However, the molecular origin of the asymmetry is not understood. We developed a new theoretical approach that gives microscopic explanations of asymmetric breathing. It is based on a stochastic description that leads to a comprehensive evaluation of dynamics by using effective free-energy landscapes. It is shown that asymmetric breathing follows the kinetically preferred pathways. In addition, it is also found that asymmetric breathing leads to a faster target search by transcription factors. Theoretical predictions, supported by computer simulations, agree with experiments. It is proposed that nature utilizes the symmetry of nucleosome breathing to achieve a better dynamic accessibility of chromatin for more efficient genetic regulation.

Stochastic nucleosome disassembly mediated by remodelers and histone fragmentation

We construct and analyze monomeric and multimeric models of the stochastic disassembly of a single nucleosome. Our monomeric model predicts the time needed for a number of histone-DNA contacts to spontaneously break, leading to dissociation of a non-fragmented histone from DNA. The dissociation process can be facilitated by DNA binding proteins or processing molecular motors that compete with histones for histone-DNA contact sites. Eigenvalue analysis of the corresponding master equation allows us to evaluate histone detachment times under both spontaneous detachment and protein-facilitated processes. We find that competitive DNA binding of remodeling proteins can significantly reduce the typical detachment time but only if these remodelers have DNA-binding affinities comparable to those of histone-DNA contact sites. In the presence of processive motors, the histone detachment rate is shown to be proportional to the product of the histone single-bond dissociation constant and the speed of motor protein procession. Our simple intact-histone model is then extended to allow for multimeric nucleosome kinetics that reveal additional pathways of disassembly. In addition to a dependence of complete disassembly times on subunit-DNA contact energies, we show how histone subunit concentrations in bulk solutions can mediate the disassembly process by rescuing partially disassembled nucleosomes. Moreover, our kinetic model predicts that remodeler binding can also bias certain pathways of nucleosome disassembly, with higher remodeler binding rates favoring intact-histone detachment.

Molecular dynamics simulations for the study of chromatin biology

The organization of Eukaryotic DNA into chromatin has profound implications for the processing of genetic information. In the past years, molecular dynamics (MD) simulations proved to be a powerful tool to investigate the mechanistic basis of chromatin biology. We review recent all-atom and coarse-grained MD studies revealing how the structure and dynamics of chromatin underlie its biological functions. We describe the latest method developments; the structural fluctuations of nucleosomes and the various factors affecting them; the organization of chromatin fibers, with particular emphasis on its liquid-like character; the interactions and dynamics of transcription factors on chromatin; and how chromatin organization is modulated by molecular motors acting on DNA.

Stochastic models of nucleosome dynamics reveal regulatory rules of stimulus-induced epigenome remodeling

The genomic positions of nucleosomes are a defining feature of the cell's epigenomic state, but signal-dependent transcription factors (SDTFs), upon activation, bind to specific genomic locations and modify nucleosome positioning. Here we leverage SDTFs as perturbation probes to learn about nucleosome dynamics in living cells. We develop Markov models of nucleosome dynamics and fit them to time course sequencing data of DNA accessibility. We find that (1) the dynamics of DNA unwrapping are significantly slower in cells than reported from cell-free experiments, (2) only models with cooperativity in wrapping and unwrapping fit the available data, (3) SDTF activity produces the highest eviction probability when its binding site is adjacent to but not on the nucleosome dyad, and (4) oscillatory SDTF activity results in high location variability. Our work uncovers the regulatory rules governing SDTF-induced nucleosome dynamics in live cells, which can predict chromatin accessibility alterations during inflammation at single-nucleosome resolution.

Chemical map-based prediction of nucleosome positioning using the Bioconductor package nuCpos

Background

Assessing the nucleosome-forming potential of specific DNA sequences is important for understanding complex chromatin organization. Methods for predicting nucleosome positioning include bioinformatics and biophysical approaches. An advantage of bioinformatics methods, which are based on in vivo nucleosome maps, is the use of natural sequences that may contain previously unknown elements involved in nucleosome positioning in vivo. The accuracy of such prediction attempts reflects the genomic coordinate resolution of the nucleosome maps applied. Nucleosome maps are constructed using micrococcal nuclease digestion followed by high-throughput sequencing (MNase-seq). However, as MNase has a strong preference for A/T-rich sequences, MNase-seq may not be appropriate for this purpose. In addition to MNase-seq-based maps, base pair-resolution chemical maps of in vivo nucleosomes from three different species (budding and fission yeasts, and mice) are currently available. However, these chemical maps have yet to be integrated into publicly available computational methods.

Results

We developed a Bioconductor package (named nuCpos) to demonstrate the superiority of chemical maps in predicting nucleosome positioning. The accuracy of chemical map-based prediction in rotational settings was higher than that of the previously developed MNase-seq-based approach. With our method, predicted nucleosome occupancy reasonably matched in vivo observations and was not affected by A/T nucleotide frequency. Effects of genetic alterations on nucleosome positioning that had been observed in living yeast cells could also be predicted. nuCpos calculates individual histone binding affinity (HBA) scores for given 147-bp sequences to examine their suitability for nucleosome formation. We also established local HBA as a new parameter to predict nucleosome formation, which was calculated for 13 overlapping nucleosomal DNA subsequences. HBA and local HBA scores for various sequences agreed well with previous in vitro and in vivo studies. Furthermore, our results suggest that nucleosomal subsegments that are disfavored in different rotational settings contribute to the defined positioning of nucleosomes.

Conclusions

Our results demonstrate that chemical map-based statistical models are beneficial for studying nucleosomal DNA features. Studies employing nuCpos software can enhance understanding of chromatin regulation and the interpretation of genetic alterations and facilitate the design of artificial sequences.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04240-2.

High-throughput AFM analysis reveals unwrapping pathways of H3 and CENP-A nucleosomes

Nucleosomes, the fundamental units of chromatin, regulate readout and expression of eukaryotic genomes. Single-molecule experiments have revealed force-induced nucleosome accessibility, but a high-resolution unwrapping landscape in the absence of external forces is currently lacking. Here, we introduce a high-throughput pipeline for the analysis of nucleosome conformations based on atomic force microscopy and automated, multi-parameter image analysis. Our data set of ∼10 000 nucleosomes reveals multiple unwrapping states corresponding to steps of 5 bp DNA. For canonical H3 nucleosomes, we observe that dissociation from one side impedes unwrapping from the other side, but in contrast to force-induced unwrapping, we find only a weak sequence-dependent asymmetry. Notably, centromeric CENP-A nucleosomes do not unwrap anti-cooperatively, in stark contrast to H3 nucleosomes. Finally, our results reconcile previous conflicting findings about the differences in height between H3 and CENP-A nucleosomes. We expect our approach to enable critical insights into epigenetic regulation of nucleosome structure and stability and to facilitate future high-throughput AFM studies that involve heterogeneous nucleoprotein complexes.

Phase-plate cryo-EM structure of the Widom 601 CENP-A nucleosome core particle reveals differential flexibility of the DNA ends

The histone H3 variant CENP-A marks centromeres epigenetically and is essential for mitotic fidelity. Previous crystallographic studies of the CENP-A nucleosome core particle (NCP) reconstituted with a human α-satellite DNA derivative revealed both DNA ends to be highly flexible, a feature important for CENP-A mitotic functions. However, recent cryo-EM studies of CENP-A NCP complexes comprising primarily Widom 601 DNA reported well-ordered DNA ends. Here, we report the cryo-EM structure of the CENP-A 601 NCP determined by Volta phase-plate imaging. The data reveal that one (‘left’) 601 DNA end is well ordered whereas the other (‘right’) end is flexible and partly detached from the histone core, suggesting sequence-dependent dynamics of the DNA termini. Indeed, a molecular dynamics simulation of the CENP-A 601 NCP confirmed the distinct dynamics of the two DNA extremities. Reprocessing the image data using two-fold symmetry yielded a cryo-EM map in which both DNA ends appeared well ordered, indicating that such an artefact may inadvertently arise if NCP asymmetry is lost during image processing. These findings enhance our understanding of the dynamic features that discriminate CENP-A from H3 nucleosomes by revealing that DNA end flexibility can be fine-tuned in a sequence-dependent manner.

Translational nucleosome positioning: A computational study

About three-quarters of eukaryotic DNA is wrapped into nucleosomes; DNA spools with a protein core. The affinity of a given DNA stretch to be incorporated into a nucleosome is known to depend on the base-pair sequence-dependent geometry and elasticity of the DNA double helix. This causes the rotational and translational positioning of nucleosomes. In this study we ask the question whether the latter can be predicted by a simple coarse-grained DNA model with sequence-dependent elasticity, the rigid base-pair model. Whereas this model is known to be rather robust in predicting rotational nucleosome positioning, we show that the translational positioning is a rather subtle effect that is dominated by the guanine-cytosine content dependence of entropy rather than energy. A correct qualitative prediction within the rigid base-pair framework can only be achieved by assuming that DNA elasticity effectively changes on complexation into the nucleosome complex. With that extra assumption we arrive at a model which gives an excellent quantitative agreement to experimental in vitro nucleosome maps, under the additional assumption that nucleosomes equilibrate their positions only locally.

Nucleosomal DNA Dynamics Mediate Oct4 Pioneer Factor Binding

Transcription factor (TF) proteins bind to DNA to regulate gene expression. Normally, accessibility to DNA is required for their function. However, in the nucleus, the DNA is often inaccessible, wrapped around histone proteins in nucleosomes forming the chromatin. Pioneer TFs are thought to induce chromatin opening by recognizing their DNA binding sites on nucleosomes. For example, Oct4, a master regulator and inducer of stem cell pluripotency, binds to DNA in nucleosomes in a sequence-specific manner. Here, we reveal the structural dynamics of nucleosomes that mediate Oct4 binding from molecular dynamics simulations. Nucleosome flexibility and the amplitude of nucleosome motions such as breathing and twisting are enhanced in nucleosomes with multiple TF binding sites. Moreover, the regions around the binding sites display higher local structural flexibility. Probing different structures of Oct4-nucleosome complexes, we show that alternative configurations in which Oct4 recognizes partial binding sites display stable TF-DNA interactions similar to those observed in complexes with free DNA and compatible with the DNA curvature and DNA-histone interactions. Therefore, we propose a structural basis for nucleosome recognition by a pioneer TF that is essential for understanding how chromatin is unraveled during cell fate conversions.

Ensembles of Breathing Nucleosomes: A Computational Study

About three-fourths of the human DNA molecules are wrapped into nucleosomes, protein spools with DNA. Nucleosomes are highly dynamic, transiently exposing their DNA through spontaneous unspooling. Recent experiments allowed to observe the DNA of an ensemble of such breathing nucleosomes through x-ray diffraction with contrast matching between the solvent and the protein core. In this study, we calculate such an ensemble through a Monte Carlo simulation of a coarse-grained nucleosome model with sequence-dependent DNA mechanics. Our analysis gives detailed insights into the sequence dependence of nucleosome breathing observed in the experiment and allows us to determine the adsorption energy of the DNA bound to the protein core as a function of the ionic strength. Moreover, we predict the breathing behavior of other potentially interesting sequences and compare the findings to earlier related experiments.

Quantitative Modeling of Nucleosome Unwrapping from Both Ends

In eukaryotic cells, DNA is packaged into chromatin where nucleosomes are the basic packaging unit. Important cellular processes including gene expression, DNA replication, and DNA repair require nucleosomal DNA to be unwrapped so that functional proteins can access their target sites, which otherwise are sterically occluded. A key question in this process is what the unwrapped conformations individual nucleosomes adopt within chromatin are. Here, we develop a concurrent nucleosome unwrapping model to address this question. We hypothesize that for a given end-to-end distance of the nucleosomal DNA, the nucleosomal DNA stochastically unwraps from the histone core from both ends independently and that this combination of unwrapping from both sides results in a significant increase in the average distance between the DNA extending from both sides of the nucleosomes. We test our model on recently published experiments using a DNA origami nanocaliper that quantifies nucleosome unwrapping and achieve good agreement between experiment and model prediction. We then investigate the DNA origami caliper distribution when attached to a hexasome (a nucleosome lacking an H2A/H2B dimer). A significant shift in the caliper angle distribution caused by the asymmetric structural features of the hexasome seen experimentally is consistent with the model. Our modeling approach may be more broadly useful to the interpretation of other studies of nucleosome dynamics, chromatin dynamics, and regulatory processes involving nucleosome unwrapping, as well as more generally to optimization of future DNA origami designs to probe mechanical properties of biomolecules.

Challenges for base excision repair enzymes: Acquiring access to damaged DNA in chromatin

Repair of damaged DNA plays a crucial role in maintaining genomic integrity and normal cell function. The base excision repair (BER) pathway is primarily responsible for removing modified nucleobases that would otherwise cause deleterious and mutagenic consequences and lead to disease. The BER process is initiated by a DNA glycosylase, which recognizes and excises the target nucleobase lesion, and is completed via downstream enzymes acting in a well-coordinated manner. A majority of our current understanding about how BER enzymes function comes from in vitro studies using free duplex DNA as a simplified model. In eukaryotes, however, BER is challenged by the packaging of genomic DNA into chromatin. The fundamental structural repeating unit of chromatin is the nucleosome, which presents a more complex substrate context than free duplex DNA for repair. In this chapter, we discuss how BER enzymes, particularly glycosylases, engage in the context of packaged DNA with insights obtained from both in vivo and in vitro studies. Furthermore, we review factors and mechanisms that can modify chromatin architecture and/or influence DNA accessibility to BER machinery, such as the geometric location of the damage site, nucleosomal DNA unwrapping, histone post-translational modifications, histone variant incorporation, and chromatin remodeling.

Critical role of histone tail entropy in nucleosome unwinding

The nucleosome is the fundamental packaging unit for the genome. It must remain tightly wound to ensure genome stability while simultaneously being flexible enough to keep the DNA molecule accessible for genome function. The set of physicochemical interactions responsible for the delicate balance between these naturally opposed processes have not been determined due to challenges in resolving partially unwound nucleosome configurations at atomic resolution. Using a near atomistic protein-DNA model and advanced sampling techniques, we calculate the free energy cost of nucleosome DNA unwinding. Our simulations identify a large energetic barrier that decouples the outer and the inner DNA unwinding into two separate processes, occurring on different time scales. This dynamical decoupling allows the exposure of outer DNA at a modest cost to ensure accessibility while keeping the inner DNA and the histone core intact to maintain stability. We also reveal that this energetic barrier arises from a delayed loss of contacts between disordered histone tails and the DNA and is, surprisingly, largely offset by an entropic contribution from these tails. Implications of this enthalpy entropy compensation for the regulation of nucleosome stability and genome function are discussed.

Histone H2A Variants Enhance the Initiation of Base Excision Repair in Nucleosomes

Substituting histone variants for their canonical counterparts can profoundly alter chromatin structure, thereby impacting multiple biological processes. Here, we investigate the influence of histone variants from the H2A family on the excision of uracil (U) by the base excision repair (BER) enzymes uracil DNA glycosylase (UDG) and single-strand selective monofunctional uracil DNA glycosylase. Using a DNA population with globally distributed U:G base pairs, enhanced excision is observed in H2A.Z and macroH2A-containing nucleosome core particles (NCPs). The U with reduced solution accessibility exhibit limited UDG activity in canonical NCPs but are more readily excised in variant NCPs, reflecting the ability of these variants to facilitate excision at sites that are otherwise poorly repaired. We also find that U with the largest increase in the level of excision in variant NCPs are clustered in regions with differential structural features between the variants and canonical H2A. Within 35-40 bp of the DNA terminus in macroH2A NCPs, the activities of both glycosylases are comparable to that on the free duplex. We show that this high level of activity results from two distinct species within the macroH2A NCP ensemble: octasomes and hexasomes. These observations reveal potential functions for H2A variants in promoting BER and preventing mutagenesis within the context of chromatin.

Shortest paths through synonymous genomes

The elasticity of the DNA double helix varies with the underlying base pair sequence. This allows one to put mechanical cues into sequences that in turn influence the packaging of DNA into nucleosomes, DNA-wrapped protein cylinders. Nucleosomes dictate a broad range of biological processes, ranging from gene regulation, recombination, and replication to chromosome condensation. Here we map base pair sequences onto graphs and use shortest paths algorithms to determine which DNA stretches are easiest or hardest to bend inside a nucleosome. We further demonstrate how genetic and mechanical information can be multiplexed by studying paths through graphs of synonymous codons. Using this method we find that nucleosomes can be placed by mechanical cues nearly everywhere on the genome of baker's yeast (Saccharomyces cerevisiae).

Local DNA Sequence Controls Asymmetry of DNA Unwrapping from Nucleosome Core Particles

DNA is tightly wrapped around histone proteins in nucleosome core particles (NCPs) yet must become accessible for processing in the cell. This accessibility, a key component of transcription regulation, is influenced by the properties of both the histone proteins and the DNA itself. Small angle x-ray scattering with contrast variation is used to examine how sequence variations affect DNA unwrapping from NCPs at different salt concentrations. Salt destabilizes NCPs, populating multiple unwrapped states as many possible unwrapping pathways are explored by the complexes. We apply coarse-grained Monte Carlo methods to generate realistic sequence-dependent unwrapped structures for the nucleosomal DNA with thermal variations. An ensemble optimization method is employed to determine the composition of the overall ensemble as electrostatic interactions are weakened. Interesting DNA-sequence-dependent differences are revealed in the unwrapping paths and equilibrium constants. These differences are correlated with specific features within the nucleic acid sequences.

The role of DNA sequence in nucleosome breathing

Roughly 3/4 of human genomes are sequestered by nucleosomes, DNA spools with a protein core, dictating a broad range of biological processes, ranging from gene regulation, recombination, and replication, to chromosome condensation. Nucleosomes are dynamical structures and temporarily expose wrapped DNA through spontaneous unspooling from either end, a process called site exposure or nucleosome breathing. Here we ask how this process is influenced by the mechanical properties of the wrapped DNA, which is known to depend on the underlying base pair sequence. Using a coarse-grained nucleosome model we calculate the accessibility profiles for site exposure. We find that the process is very sensitive to sequence effects, so that evolution could potentially tune the accessibility of nucleosomal DNA and would only need a small number of mutations to do so.

Physics behind the mechanical nucleosome positioning code

The positions along DNA molecules of nucleosomes, the most abundant DNA-protein complexes in cells, are influenced by the sequence-dependent DNA mechanics and geometry. This leads to the "nucleosome positioning code", a preference of nucleosomes for certain sequence motives. Here we introduce a simplified model of the nucleosome where a coarse-grained DNA molecule is frozen into an idealized superhelical shape. We calculate the exact sequence preferences of our nucleosome model and find it to reproduce qualitatively all the main features known to influence nucleosome positions. Moreover, using well-controlled approximations to this model allows us to come to a detailed understanding of the physics behind the sequence preferences of nucleosomes.

Designing nucleosomal force sensors

About three quarters of our DNA is wrapped into nucleosomes: DNA spools with a protein core. It is well known that the affinity of a given DNA stretch to be incorporated into a nucleosome depends on the geometry and elasticity of the basepair sequence involved, causing the positioning of nucleosomes. Here we show that DNA elasticity can have a much deeper effect on nucleosomes than just their positioning: it affects their "identities". Employing a recently developed computational algorithm, the mutation Monte Carlo method, we design nucleosomes with surprising physical characteristics. Unlike any other nucleosomes studied so far, these nucleosomes are short-lived when put under mechanical tension whereas other physical properties are largely unaffected. This suggests that the nucleosome, the most abundant DNA-protein complex in our cells, might more properly be considered a class of complexes with a wide array of physical properties, and raises the possibility that evolution has shaped various nucleosome species according to their genomic context.

Benchmarking and refining probability-based models for nucleosome-DNA interaction

Background

In investigations of nucleosome positioning preferences, a model that assigns an affinity to a given sequence is necessary to make predictions. One important class of models, which treats a nucleosome sequence as a Markov chain, has been applied with success when informed with experimentally measured nucleosomal sequence preferences.

Results

We find that we can also use such models as a fast approximative scheme for computationally expensive biophysical models, vastly increasing their reach. Employing these models in this way also allows us to benchmark them for the first time. Doing so for the approximative in silico models indirectly tells us about the accuracy we can expect of them when applied to real data.

Conclusion

We find that models presented in the literature should perform well, but this performance depends on factors such as the order of the Markov model, the preprocessing of the probability distributions on which the model is based, and the size and quality of the sequence ensemble from which those distributions are calculated.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-017-1569-0) contains supplementary material, which is available to authorized users.

Genomes of Multicellular Organisms Have Evolved to Attract Nucleosomes to Promoter Regions

Sequences that influence nucleosome positioning in promoter regions, and their relation to gene regulation, have been the topic of much research over the last decade. In yeast, significant nucleosome-depleted regions are found, which facilitate transcription. With the arrival of nucleosome positioning maps for the human genome, it was discovered that in our genome, unlike in that of yeast, promoters encode for high nucleosome occupancy. In this work, we look at the genomes of a range of different organisms, to provide a catalog of nucleosome positioning signals in promoters across the tree of life. We utilize a computational model of the nucleosome, based on crystallographic analyses of the structure and elasticity of the nucleosome, to predict the nucleosome positioning signals in promoter regions. To be able to apply our model to large genomic datasets, we introduce an approximative scheme that makes use of the limited range of correlations in nucleosomal sequence preferences to create a computationally efficient approximation of the full biophysical model. Our predictions show that a clear distinction between unicellular and multicellular life is visible in the intrinsically encoded nucleosome affinity. Furthermore, the strength of the nucleosome positioning signals correlates with the complexity of the organism. We conclude that encoding for high nucleosome occupancy, as in the human genome, is in fact a universal feature of multicellular life.