Thermodynamic modeling of genome-wide nucleosome depleted regions in yeast

Native nucleosome-positioning elements for the investigation of nucleosome repositioning

Nucleosome repositioning is essential for establishing nucleosome-depleted regions (NDRs) to initiate transcription. This process has been extensively studied using structural, biochemical, and single-molecule approaches, which require homogenously positioned nucleosomes. This is often achieved using the Widom 601 sequence, a highly efficient nucleosome positioning element (NPE) selected for its unusually strong binding to the H3-H4 histone tetramer. Due to the artificial nature of 601, native NPEs are needed to explore the role of DNA sequence in nucleosome repositioning. Here, we characterize the position distributions and nucleosome formation free energy for a set of yeast native nucleosomes (YNNs) from Saccharomyces cerevisiae. We show these native NPEs can be used in biochemical studies of nucleosome repositioning by transcription factors (TFs) and the chromatin remodeler Chd1. TFs could directly reposition a fraction of nucleosomes containing native NPEs, but not 601-containing nucleosomes. In contrast, partial unwrapping was similar for 601 and native NPE sequences, and the rate of ATP-dependent remodeling by Chd1 was within the range of the fast and slow directions of the 601 nucleosomes. This set of native NPEs provides an alternative to the 601 NPE that can be used for probing the repositioning of nucleosomes that contain native DNA sequences.

Thermodynamics of nucleosome breathing and positioning

Nucleosomes are fundamental units of chromatin in which a length of genomic DNA is wrapped around a histone octamer spool in a left-handed superhelix. Large-scale nucleosome maps show a wide distribution of DNA wrapping lengths, which in some cases are tens of base pairs (bp) shorter than the 147 bp canonical wrapping length observed in nucleosome crystal structures. Here, we develop a thermodynamic model that assumes a constant free energy cost of unwrapping a nucleosomal bp. Our model also incorporates linker DNA-short DNA segments between neighboring nucleosomes imposed by the folding of nucleosome arrays into chromatin fibers and other higher-order chromatin structures. We use this model to study nucleosome positioning and occupancy in the presence of nucleosome "breathing"-partial unwrapping and rewrapping of nucleosomal DNA due to interactions with the neighboring particles. We find that, as the unwrapping cost per bp and the chemical potential are varied, the nucleosome arrays are characterized by three distinct states, with low, intermediate, and high densities. The transition between the latter two states proceeds through an equiprobable state in which all nucleosome wrapping lengths are equally likely. We study the equiprobable state theoretically using a mean-field approach, obtaining an excellent agreement with numerical simulations. Finally, we use our model to reproduce S. cerevisiae nucleosome occupancy profiles observed in the vicinity of transcription start sites, as well as genome-wide distributions of nucleosome wrapping lengths. Overall, our results highlight the key role of partial nucleosome unwrapping in shaping the genome-wide patterns of nucleosome positioning and occupancy.

DNA methylation-based high-resolution mapping of long-distance chromosomal interactions in nucleosome-depleted regions

3C-based methods have significantly advanced our understanding of 3D genome organization. However, it remains a formidable task to precisely capture long-range chromosomal interactions between individual loci, such as those between promoters and distal enhancers. Here, we present Methyltransferase Targeting-based chromosome Architecture Capture (MTAC), a method that maps the contacts between a target site (viewpoint) and the rest of the genome in budding yeast with high resolution and sensitivity. MTAC detects hundreds of intra- and inter-chromosomal interactions within nucleosome-depleted regions (NDRs) that cannot be captured by 4C, Hi-C, or Micro-C. By applying MTAC to various viewpoints, we find that (1) most long-distance chromosomal interactions detected by MTAC reflect tethering by the nuclear pore complexes (NPCs), (2) genes co-regulated by methionine assemble into inter-chromosomal clusters near NPCs upon activation, (3) mediated by condensin, the mating locus forms a highly specific interaction with the recombination enhancer (RE) in a mating-type specific manner, and (4) correlation of MTAC signals among NDRs reveal spatial mixing and segregation of the genome. Overall, these results demonstrate MTAC as a powerful tool to resolve fine-scale long-distance chromosomal interactions and provide insights into the 3D genome organization.

Using Synthetic DNA Libraries to Investigate Chromatin and Gene Regulation

Despite the recent explosion in genome-wide studies in chromatin and gene regulation, we are still far from extracting a set of genetic rules that can predict the function of the regulatory genome. One major reason for this deficiency is that gene regulation is a multi-layered process that involves an enormous variable space, which cannot be fully explored using native genomes. This problem can be partially solved by introducing synthetic DNA libraries into cells, a method that can test the regulatory roles of thousands to millions of sequences with limited variables. Here, we review recent applications of this method to study transcription factor (TF) binding, nucleosome positioning, and transcriptional activity. We discuss the design principles, experimental procedures, and major findings from these studies and compare the pros and cons of different approaches.

Investigating pioneer factor activity and its coordination with chromatin remodelers using integrated synthetic oligo assay

Chromatin accessibility is regulated by pioneer factors (PFs) and chromatin remodelers (CRs). Here, we present a protocol, based on integrated synthetic oligonucleotide libraries in yeast, to systematically interrogate the nucleosome-displacing activities of PFs and their coordination with CRs. We describe steps for designing oligonucleotide sequences, constructing yeast libraries, measuring nucleosome configurations, and data analyses. This approach potentially can be adapted for use in higher eukaryotes to investigate the activities of many types of chromatin-associated factors. For complete details on the use and execution of this protocol, please refer to Yan et al.,1 and Chen et al.2.

Functional interplay between Mediator and RSC chromatin remodeling complex controls nucleosome-depleted region maintenance at promoters

Chromatin organization is crucial for transcriptional regulation in eukaryotes. Mediator is an essential and conserved co-activator thought to act in concert with chromatin regulators. However, it remains largely unknown how their functions are coordinated. Here, we provide evidence in the yeast Saccharomyces cerevisiae that Mediator establishes physical contact with RSC (Remodels the Structure of Chromatin), a conserved and essential chromatin remodeling complex that is crucial for nucleosome-depleted region (NDR) formation. We determine the role of Mediator-RSC interaction in their chromatin binding, nucleosome occupancy, and transcription on a genomic scale. Mediator and RSC co-localize on wide NDRs of promoter regions, and specific Mediator mutations affect nucleosome eviction and TSS-associated +1 nucleosome stability. This work shows that Mediator contributes to RSC remodeling function to shape NDRs and maintain chromatin organization on promoter regions. It will help in our understanding of transcriptional regulation in the chromatin context relevant for severe diseases.

Predicting nucleosome positioning using statistical equilibrium models in budding yeast

We present a protocol using thermodynamic models to predict nucleosome positioning with transcription factors (TFs) and chromatin remodelers. We describe step-by-step approaches to annotate genome-wide nucleosome-depleted regions (NDRs), compute nucleosome and TF occupancy, optimize parameters, and evaluate model performance. These models identify nucleosome-displacing TFs in the budding yeast genome and predict the locations and sizes of NDRs solely based on DNA sequence and TF motifs. The protocol can be applied to all organisms with prior knowledge of TF motifs. For complete details on the use and execution of this protocol, please refer to Kharerin and Bai (2021).1.

Partitioned usage of chromatin remodelers by nucleosome-displacing factors

Nucleosome-displacing-factors (NDFs) in yeast, similar to pioneer factors in higher eukaryotes, can open closed chromatin and generate nucleosome-depleted regions (NDRs). NDRs in yeast are also affected by ATP-dependent chromatin remodelers (CRs). However, how NDFs and CRs coordinate in nucleosome invasion and NDR formation is still unclear. Here, we design a high-throughput method to systematically study the interplay between NDFs and CRs. By combining an integrated synthetic oligonucleotide library with DNA methyltransferase-based, single-molecule nucleosome mapping, we measure the impact of CRs on NDRs generated by individual NDFs. We find that CRs are dispensable for nucleosome invasion by NDFs, and they function downstream of NDF binding to modulate the NDR length. A few CRs show high specificity toward certain NDFs; however, in most cases, CRs are recruited in a factor-nonspecific and NDR length-dependent manner. Overall, our study provides a framework to investigate how NDFs and CRs cooperate to regulate chromatin opening.

Chromatin accessibility in yeast is regulated by nucleosome-displacing-factors (NDFs) and chromatin remodelers (CRs). Chen et al. show that NDFs first invade into nucleosomes and then recruit CRs to modulate the NDR length. NDF-specific and NDR length-dependent recruitment of CRs allow partitioned usage of CRs by NDFs.

Nucleosome-directed replication origin licensing independent of a consensus DNA sequence

The numerous enzymes and cofactors involved in eukaryotic DNA replication are conserved from yeast to human, and the budding yeast Saccharomyces cerevisiae (S.c.) has been a useful model organism for these studies. However, there is a gap in our knowledge of why replication origins in higher eukaryotes do not use a consensus DNA sequence as found in S.c. Using in vitro reconstitution and single-molecule visualization, we show here that S.c. origin recognition complex (ORC) stably binds nucleosomes and that ORC-nucleosome complexes have the intrinsic ability to load the replicative helicase MCM double hexamers onto adjacent nucleosome-free DNA regardless of sequence. Furthermore, we find that Xenopus laevis nucleosomes can substitute for yeast ones in engaging with ORC. Combined with re-analyses of genome-wide ORC binding data, our results lead us to propose that the yeast origin recognition machinery contains the cryptic capacity to bind nucleosomes near a nucleosome-free region and license origins, and that this nucleosome-directed origin licensing paradigm generalizes to all eukaryotes.

Superstructure Detection in Nucleosome Distribution Shows Common Pattern within a Chromosome and within the Genome

Nucleosome positioning plays an important role in crucial biological processes such as replication, transcription, and gene regulation. It has been widely used to predict the genome’s function and chromatin organisation. So far, the studies of patterns in nucleosome positioning have been limited to transcription start sites, CTCFs binding sites, and some promoter and loci regions. The genome-wide organisational pattern remains unknown. We have developed a theoretical model to coarse-grain nucleosome positioning data in order to obtain patterns in their distribution. Using hierarchical clustering on the auto-correlation function of this coarse-grained nucleosome positioning data, a genome-wide clustering is obtained for Candida albicans. The clustering shows the existence beyond hetero- and eu-chromatin inside the chromosomes. These non-trivial clusterings correspond to different nucleosome distributions and gene densities governing differential gene expression patterns. Moreover, these distribution patterns inside the chromosome appeared to be conserved throughout the genome and within species. The pipeline of the coarse grain nucleosome positioning sequence to identify underlying genomic organisation used in our study is novel, and the classifications obtained are unique and consistent.

The Active Mechanism of Nucleosome Depletion by Poly(dA:dT) Tracts In Vivo

Poly(dA:dT) tracts cause nucleosome depletion in many species, e.g., at promoters and replication origins. Their intrinsic biophysical sequence properties make them stiff and unfavorable for nucleosome assembly, as probed by in vitro nucleosome reconstitution. The mere correlation between nucleosome depletion over poly(dA:dT) tracts in in vitro reconstituted and in in vivo chromatin inspired an intrinsic nucleosome exclusion mechanism in vivo that is based only on DNA and histone properties. However, we compile here published and new evidence that this correlation does not reflect mechanistic causation. (1) Nucleosome depletion over poly(dA:dT) in vivo is not universal, e.g., very weak in S. pombe. (2) The energy penalty for incorporating poly(dA:dT) tracts into nucleosomes is modest (<10%) relative to ATP hydrolysis energy abundantly invested by chromatin remodelers. (3) Nucleosome depletion over poly(dA:dT) is much stronger in vivo than in vitro if monitored without MNase and (4) actively maintained in vivo. (5) S. cerevisiae promoters evolved a strand-biased poly(dA) versus poly(dT) distribution. (6) Nucleosome depletion over poly(dA) is directional in vivo. (7) The ATP dependent chromatin remodeler RSC preferentially and directionally displaces nucleosomes towards 5′ of poly(dA). Especially distribution strand bias and displacement directionality would not be expected for an intrinsic mechanism. Together, this argues for an in vivo mechanism where active and species-specific read out of intrinsic sequence properties, e.g., by remodelers, shapes nucleosome organization.