New insights into two distinct nucleosome distributions: comparison of cross-platform positioning datasets in the yeast genome

Coupling of replisome movement with nucleosome dynamics can contribute to the parent-daughter information transfer

Positioning of nucleosomes along the genomic DNA is crucial for many cellular processes that include gene regulation and higher order packaging of chromatin. The question of how nucleosome-positioning information from a parent chromatin gets transferred to the daughter chromatin is highly intriguing. Accounting for experimentally known coupling between replisome movement and nucleosome dynamics, we propose a model that can obtain de novo nucleosome assembly similar to what is observed in recent experiments. Simulating nucleosome dynamics during replication, we argue that short pausing of the replication fork, associated with nucleosome disassembly, can be a event crucial for communicating nucleosome positioning information from parent to daughter. We show that the interplay of timescales between nucleosome disassembly (τ_p) at the replication fork and nucleosome sliding behind the fork (τ_s) can give rise to a rich ‘phase diagram’ having different inherited patterns of nucleosome organization. Our model predicts that only when τ_p ≥ τ_s the daughter chromatin can inherit nucleosome positioning of the parent.

NucTools: analysis of chromatin feature occupancy profiles from high-throughput sequencing data

Background

Biomedical applications of high-throughput sequencing methods generate a vast amount of data in which numerous chromatin features are mapped along the genome. The results are frequently analysed by creating binary data sets that link the presence/absence of a given feature to specific genomic loci. However, the nucleosome occupancy or chromatin accessibility landscape is essentially continuous. It is currently a challenge in the field to cope with continuous distributions of deep sequencing chromatin readouts and to integrate the different types of discrete chromatin features to reveal linkages between them.

Results

Here we introduce the NucTools suite of Perl scripts as well as MATLAB- and R-based visualization programs for a nucleosome-centred downstream analysis of deep sequencing data. NucTools accounts for the continuous distribution of nucleosome occupancy. It allows calculations of nucleosome occupancy profiles averaged over several replicates, comparisons of nucleosome occupancy landscapes between different experimental conditions, and the estimation of the changes of integral chromatin properties such as the nucleosome repeat length. Furthermore, NucTools facilitates the annotation of nucleosome occupancy with other chromatin features like binding of transcription factors or architectural proteins, and epigenetic marks like histone modifications or DNA methylation. The applications of NucTools are demonstrated for the comparison of several datasets for nucleosome occupancy in mouse embryonic stem cells (ESCs) and mouse embryonic fibroblasts (MEFs).

Conclusions

The typical workflows of data processing and integrative analysis with NucTools reveal information on the interplay of nucleosome positioning with other features such as for example binding of a transcription factor CTCF, regions with stable and unstable nucleosomes, and domains of large organized chromatin K9me2 modifications (LOCKs). As potential limitations and problems we discuss how inter-replicate variability of MNase-seq experiments can be addressed.

Electronic supplementary material

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

Nucleosome architecture throughout the cell cycle

Nucleosomes provide additional regulatory mechanisms to transcription and DNA replication by mediating the access of proteins to DNA. During the cell cycle chromatin undergoes several conformational changes, however the functional significance of these changes to cellular processes are largely unexplored. Here, we present the first comprehensive genome-wide study of nucleosome plasticity at single base-pair resolution along the cell cycle in Saccharomyces cerevisiae. We determined nucleosome organization with a specific focus on two regulatory regions: transcription start sites (TSSs) and replication origins (ORIs). During the cell cycle, nucleosomes around TSSs display rearrangements in a cyclic manner. In contrast to gap (G1 and G2) phases, nucleosomes have a fuzzier organization during S and M phases, Moreover, the choreography of nucleosome rearrangements correlate with changes in gene expression during the cell cycle, indicating a strong association between nucleosomes and cell cycle-dependent gene functionality. On the other hand, nucleosomes are more dynamic around ORIs along the cell cycle, albeit with tighter regulation in early firing origins, implying the functional role of nucleosomes on replication origins. Our study provides a dynamic picture of nucleosome organization throughout the cell cycle and highlights the subsequent impact on transcription and replication activity.

Nucleosome positioning: resources and tools online

Nucleosome positioning is an important process required for proper genome packing and its accessibility to execute the genetic program in a cell-specific, timely manner. In the recent years hundreds of papers have been devoted to the bioinformatics, physics and biology of nucleosome positioning. The purpose of this review is to cover a practical aspect of this field, namely, to provide a guide to the multitude of nucleosome positioning resources available online. These include almost 300 experimental datasets of genome-wide nucleosome occupancy profiles determined in different cell types and more than 40 computational tools for the analysis of experimental nucleosome positioning data and prediction of intrinsic nucleosome formation probabilities from the DNA sequence. A manually curated, up to date list of these resources will be maintained at http://generegulation.info.

Comparative analysis of methods for genome-wide nucleosome cartography

Nucleosomes contribute to compacting the genome into the nucleus and regulate the physical access of regulatory proteins to DNA either directly or through the epigenetic modifications of the histone tails. Precise mapping of nucleosome positioning across the genome is, therefore, essential to understanding the genome regulation. In recent years, several experimental protocols have been developed for this purpose that include the enzymatic digestion, chemical cleavage or immunoprecipitation of chromatin followed by next-generation sequencing of the resulting DNA fragments. Here, we compare the performance and resolution of these methods from the initial biochemical steps through the alignment of the millions of short-sequence reads to a reference genome to the final computational analysis to generate genome-wide maps of nucleosome occupancy. Because of the lack of a unified protocol to process data sets obtained through the different approaches, we have developed a new computational tool (NUCwave), which facilitates their analysis, comparison and assessment and will enable researchers to choose the most suitable method for any particular purpose. NUCwave is freely available at http://nucleosome.usal.es/nucwave along with a step-by-step protocol for its use.

Nucleosome positioning and kinetics near transcription-start-site barriers are controlled by interplay between active remodeling and DNA sequence

We investigate how DNA sequence, ATP-dependent chromatin remodeling and nucleosome-depleted ‘barriers’ co-operate to determine the kinetics of nucleosome organization, in a stochastic model of nucleosome positioning and dynamics. We find that ‘statistical’ positioning of nucleosomes against ‘barriers’, hypothesized to control chromatin structure near transcription start sites, requires active remodeling and therefore cannot be described using equilibrium statistical mechanics. We show that, unlike steady-state occupancy, DNA site exposure kinetics near a barrier is dominated by DNA sequence rather than by proximity to the barrier itself. The timescale for formation of positioning patterns near barriers is proportional to the timescale for active nucleosome eviction. We also show that there are strong gene-to-gene variations in nucleosome positioning near barriers, which are eliminated by averaging over many genes. Our results suggest that measurement of nucleosome kinetics can reveal information about sequence-dependent regulation that is not apparent in steady-state nucleosome occupancy.

Impact of methylation on the physical properties of DNA

There is increasing evidence for the presence of an alternative code imprinted in the genome that might contribute to gene expression regulation through an indirect reading mechanism. In mammals, components of this coarse-grained regulatory mechanism include chromatin structure and epigenetic signatures, where d(CpG) nucleotide steps are key players. We report a comprehensive experimental and theoretical study of d(CpG) steps that provides a detailed description of their physical characteristics and the impact of cytosine methylation on these properties. We observed that methylation changes the physical properties of d(CpG) steps, having a dramatic effect on enriched CpG segments, such as CpG islands. We demonstrate that methylation reduces the affinity of DNA to assemble into nucleosomes, and can affect nucleosome positioning around transcription start sites. Overall, our results suggest a mechanism by which the basic physical properties of the DNA fiber can explain parts of the cellular epigenetic regulatory mechanisms.

Nucleosome positioning in a model of active chromatin remodeling enzymes

Accounting for enzyme-mediated active sliding, disassembly, and sequence-dependent positioning of nucleosomes, we simulate nucleosome occupancy over cell-cycle-scale times using a stochastic kinetic model. We show that ATP-dependent active nucleosome sliding and nucleosome removal processes are essential to obtain in vivo-like nucleosome positioning. While active sliding leads to dense nucleosome filling, sliding events alone cannot ensure sequence-dependent nucleosome positioning: Active nucleosome removal is the crucial remodeling event that drives positioning. We also show that remodeling activity changes nucleosome dynamics from glassy to liquid-like, and that remodeling dramatically influences exposure dynamics of promoter regions.

Structural constraints revealed in consistent nucleosome positions in the genome of S. cerevisiae

Background

Recent advances in the field of high-throughput genomics have rendered possible the performance of genome-scale studies to define the nucleosomal landscapes of eukaryote genomes. Such analyses are aimed towards providing a better understanding of the process of nucleosome positioning, for which several models have been suggested. Nevertheless, questions regarding the sequence constraints of nucleosomal DNA and how they may have been shaped through evolution remain open. In this paper, we analyze in detail different experimental nucleosome datasets with the aim of providing a hypothesis for the emergence of nucleosome-forming sequences.

Results

We compared the complete sets of nucleosome positions for the budding yeast (Saccharomyces cerevisiae) as defined in the output of two independent experiments with the use of two different experimental techniques. We found that < 10% of the experimentally defined nucleosome positions were consistently positioned in both datasets. This subset of well-positioned nucleosomes, when compared with the bulk, was shown to have particular properties at both sequence and structural levels. Consistently positioned nucleosomes were also shown to occur preferentially in pairs of dinucleosomes, and to be surprisingly less conserved compared with their adjacent nucleosome-free linkers.

Conclusion

Our findings may be combined into a hypothesis for the emergence of a weak nucleosome-positioning code. According to this hypothesis, consistent nucleosomes may be partly guided by nearby nucleosome-free regions through statistical positioning. Once established, a set of well-positioned consistent nucleosomes may impose secondary constraints that further shape the structure of the underlying DNA. We were able to capture these constraints through the application of a recently introduced structural property that is related to the symmetry of DNA curvature. Furthermore, we found that both consistently positioned nucleosomes and their adjacent nucleosome-free regions show an increased tendency towards conservation of this structural feature.