Chemical map of Schizosaccharomyces pombe reveals species-specific features in nucleosome positioning

A fission yeast CENP-B homolog Abp1 prevents RNAi-mediated heterochromatin formation at ribosomal DNA repeats

In response to nutritional starvation, living cells sensitively regulate the production rates of molecules required for survival. Under glucose starvation, a facultative heterochromatinization of ribosomal DNA is considered to regulate ribosomal RNA production. However, the molecular mechanism is still unclear. Here, we report a novel function of CENP-B homolog Abp1 in forming facultative heterochromatin at ribosomal DNA repeats. We find that the loss of Abp1 induces an ectopic nucleosome assembly at rDNA repeats. Interestingly, the loss of Abp1 induces two mutually exclusive changes at ribosomal DNA repeats: an excess accumulation of methylation of histone H3 at lysine 9, a hallmark of heterochromatin, and an active RNA polymerase II transcription. This excess heterochromatin represses ribosomal RNA expression and requires RNA interference machinery for its formation. Furthermore, we show that the excess heterochromatin does not affect cellular viability under glucose starvation but prevents the return to the proliferation cycle in recovering glucose-rich conditions. Since glucose starvation rapidly induces partial Abp1 disassociation from ribosomal DNA repeats, we propose that Abp1 regulates activity of RNA polymerase II transcription that is paradoxically required for RNA interference-mediated heterochromatin formation and controls an appropriate level of heterochromatinization at ribosomal DNA repeats under glucose starvation.

Genome-Wide Nucleosome Mapping by H3Q85C-Directed Chemical Cleavage in Saccharomyces cerevisiae

Mapping the position of nucleosomes in vivo is key to our understanding of chromatin-based processes such as gene transcription, DNA replication, and repair. Methods based on micrococcal nuclease (MNase) digestion of chromatin are widely used to map nucleosomes, but these methods suffer from some limitations. More recently, nucleosome mapping methods that take advantage of our ability to convert cysteine residues-carefully inserted in histone proteins-into nucleases provide alternatives to MNase-based assays. Here, we provide a detailed protocol for the mapping of nucleosomes via H3Q85C-directed DNA cleavage in Saccharomyces cerevisiae.

Nucleosome Spacing Can Fine-Tune Higher Order Chromatin Assembly

Cellular chromatin displays heterogeneous structure and dynamics, properties that control diverse nuclear processes. Models invoke phase separation of conformational ensembles of chromatin fibers as a mechanism regulating chromatin organization in vivo. Here we combine biochemistry and molecular dynamics simulations to examine, at single base-pair resolution, how nucleosome spacing controls chromatin phase separation. We show that as DNA linkers extend from 25 bp to 30 bp, as examplars of 10N+5 and 10N (integer N) bp lengths, chromatin condensates become less thermodynamically stable and nucleosome mobility increases. Simulations reveal that this is due to trade-offs between inter- and intramolecular nucleosome stacking, favored by rigid 10N+5 and 10N bp linkers, respectively. A remodeler can induce or inhibit phase separation by moving nucleosomes, changing the balance between intra- and intermolecular stacking. The intrinsic phase separation capacity of chromatin enables fine tuning of compaction and dynamics, likely contributing to heterogeneous chromatin organization in vivo.

Physical modeling of nucleosome clustering in euchromatin resulting from interactions between epigenetic reader proteins

Euchromatin is an accessible phase of genetic material containing genes that encode proteins with increased expression levels. The structure of euchromatin in vitro has been described as a 30-nm fiber formed from ordered nucleosome arrays. However, recent advances in microscopy have revealed an in vivo euchromatin architecture that is much more disordered, characterized by variable-length linker DNA and sporadic nucleosome clusters. In this work, we develop a theoretical model to elucidate factors contributing to the disordered in vivo architecture of euchromatin. We begin by developing a 1D model of nucleosome positioning that captures the interactions between bound epigenetic reader proteins to predict the distribution of DNA linker lengths between adjacent nucleosomes. We then use the predicted linker lengths to construct 3D chromatin configurations consistent with the physical properties of DNA within the nucleosome array, and we evaluate the distribution of nucleosome cluster sizes in those configurations. Our model reproduces experimental cluster-size distributions, which are dramatically influenced by the local pattern of epigenetic marks and the concentration of reader proteins. Based on our model, we attribute the disordered arrangement of euchromatin to the heterogeneous binding of reader proteins and subsequent short-range interactions between bound reader proteins on adjacent nucleosomes. By replicating experimental results with our physics-based model, we propose a mechanism for euchromatin organization in the nucleus that impacts gene regulation and the maintenance of epigenetic marks.

Energy-driven genome regulation by ATP-dependent chromatin remodellers

The packaging of DNA into chromatin in eukaryotes regulates gene transcription, DNA replication and DNA repair. ATP-dependent chromatin remodelling enzymes (re)arrange nucleosomes at the first level of chromatin organization. Their Snf2-type motor ATPases alter histone-DNA interactions through a common DNA translocation mechanism. Whether remodeller activities mainly catalyse nucleosome dynamics or accurately co-determine nucleosome organization remained unclear. In this Review, we discuss the emerging mechanisms of chromatin remodelling: dynamic remodeller architectures and their interactions, the inner workings of the ATPase cycle, allosteric regulation and pathological dysregulation. Recent mechanistic insights argue for a decisive role of remodellers in the energy-driven self-organization of chromatin, which enables both stability and plasticity of genome regulation - for example, during development and stress. Different remodellers, such as members of the SWI/SNF, ISWI, CHD and INO80 families, process (epi)genetic information through specific mechanisms into distinct functional outputs. Combinatorial assembly of remodellers and their interplay with histone modifications, histone variants, DNA sequence or DNA-bound transcription factors regulate nucleosome mobilization or eviction or histone exchange. Such input-output relationships determine specific nucleosome positions and compositions with distinct DNA accessibilities and mediate differential genome regulation. Finally, remodeller genes are often mutated in diseases characterized by genome dysregulation, notably in cancer, and we discuss their physiological relevance.

The structure of the NuA4-Tip60 complex reveals the mechanism and importance of long-range chromatin modification

Histone acetylation regulates most DNA transactions and is dynamically controlled by highly conserved enzymes. The only essential histone acetyltransferase (HAT) in yeast, Esa1, is part of the 1-MDa NuA4 complex, which plays pivotal roles in both transcription and DNA-damage repair. NuA4 has the unique capacity to acetylate histone targets located several nucleosomes away from its recruitment site. Neither the molecular mechanism of this activity nor its physiological importance are known. Here we report the structure of the Pichia pastoris NuA4 complex, with its core resolved at 3.4-Å resolution. Three subunits, Epl1, Eaf1 and Swc4, intertwine to form a stable platform that coordinates all other modules. The HAT module is firmly anchored into the core while retaining the ability to stretch out over a long distance. We provide structural, biochemical and genetic evidence that an unfolded linker region of the Epl1 subunit is critical for this long-range activity. Specifically, shortening the Epl1 linker causes severe growth defects and reduced H4 acetylation levels over broad chromatin regions in fission yeast. Our work lays the foundations for a mechanistic understanding of NuA4's regulatory role and elucidates how its essential long-range activity is attained.

Sequence dependencies and mutation rates of localized mutational processes in cancer

BACKGROUND

Cancer mutations accumulate through replication errors and DNA damage coupled with incomplete repair. Individual mutational processes often show nucleotide sequence and functional region preferences. As a result, some sequence contexts mutate at much higher rates than others, with additional variation found between functional regions. Mutational hotspots, with recurrent mutations across cancer samples, represent genomic positions with elevated mutation rates, often caused by highly localized mutational processes.

METHODS

We count the 11-mer genomic sequences across the genome, and using the PCAWG set of 2583 pan-cancer whole genomes, we associate 11-mers with mutational signatures, hotspots of single nucleotide variants, and specific genomic regions. We evaluate the mutation rates of individual and combined sets of 11-mers and derive mutational sequence motifs.

RESULTS

We show that hotspots generally identify highly mutable sequence contexts. Using these, we show that some mutational signatures are enriched in hotspot sequence contexts, corresponding to well-defined sequence preferences for the underlying localized mutational processes. This includes signature 17b (of unknown etiology) and signatures 62 (POLE deficiency), 7a (UV), and 72 (linked to lymphomas). In some cases, the mutation rate and sequence preference increase further when focusing on certain genomic regions, such as signature 62 in transcribed regions, where the mutation rate is increased up to 9-folds over cancer type and mutational signature average.

CONCLUSIONS

We summarize our findings in a catalog of localized mutational processes, their sequence preferences, and their estimated mutation rates.

An essential Noc3p dimerization cycle mediates ORC double-hexamer formation in replication licensing

Replication licensing, a prerequisite of DNA replication, helps to ensure once-per-cell-cycle genome duplication. Some DNA replication-initiation proteins are sequentially loaded onto replication origins to form pre-replicative complexes (pre-RCs). ORC and Noc3p bind replication origins throughout the cell cycle, providing a platform for pre-RC assembly. We previously reported that cell cycle-dependent ORC dimerization is essential for the chromatin loading of the symmetric MCM double-hexamers. Here, we used Saccharomyces cerevisiae separation-of-function NOC3 mutants to confirm the separable roles of Noc3p in DNA replication and ribosome biogenesis. We also show that an essential and cell cycle-dependent Noc3p dimerization cycle regulates the ORC dimerization cycle. Noc3p dimerizes at the M-to-G1 transition and de-dimerizes in S-phase. The Noc3p dimerization cycle coupled with the ORC dimerization cycle enables replication licensing, protects nascent sister replication origins after replication initiation, and prevents re-replication. This study has revealed a new mechanism of replication licensing and elucidated the molecular mechanism of Noc3p as a mediator of ORC dimerization in pre-RC formation.

Chromatin localization of nucleophosmin organizes ribosome biogenesis

Ribosome biogenesis takes place in the nucleolus, a nuclear membrane-less organelle. Although well studied, it remains unknown how nascent ribosomal subunits separate from the central chromatin compartment and move to the outer granular component, where maturation occurs. We find that the Schizosaccharomyces pombe nucleophosmin-like protein Fkbp39 localizes to rDNA sites encoding the 60S subunit rRNA, and this localization contributes to its specific association with nascent 60S subunits. Fkbp39 dissociates from chromatin to bind nascent 60S subunits, causing the latter to partition away from chromatin and from nascent 40S subunits through liquid-liquid phase separation. In vivo, Fkbp39 binding directs the translocation of nascent 60S subunits toward the nucleophosmin-rich granular component. This process increases the efficiency of 60S subunit assembly, facilitating the incorporation of 60S RNA domain III. Thus, chromatin localization determines the specificity of nucleophosmin in sorting nascent ribosomal subunits and coordinates their movement into specialized assembly compartments within the nucleolus.

Deciphering the mechanical code of the genome and epigenome

Diverse DNA-deforming processes are impacted by the local mechanical and structural properties of DNA, which in turn depend on local sequence and epigenetic modifications. Deciphering this mechanical code (that is, this dependence) has been challenging due to the lack of high-throughput experimental methods. Here we present a comprehensive characterization of the mechanical code. Utilizing high-throughput measurements of DNA bendability via loop-seq, we quantitatively established how the occurrence and spatial distribution of dinucleotides, tetranucleotides and methylated CpG impact DNA bendability. We used our measurements to develop a physical model for the sequence and methylation dependence of DNA bendability. We validated the model by performing loop-seq on mouse genomic sequences around transcription start sites and CTCF-binding sites. We applied our model to test the predictions of all-atom molecular dynamics simulations and to demonstrate that sequence and epigenetic modifications can mechanically encode regulatory information in diverse contexts.

Chromatin accessibility: methods, mechanisms, and biological insights

Access to DNA is a prerequisite to the execution of essential cellular processes that include transcription, replication, chromosomal segregation, and DNA repair. How the proteins that regulate these processes function in the context of chromatin and its dynamic architectures is an intensive field of study. Over the past decade, genome-wide assays and new imaging approaches have enabled a greater understanding of how access to the genome is regulated by nucleosomes and associated proteins. Additional mechanisms that may control DNA accessibility in vivo include chromatin compaction and phase separation - processes that are beginning to be understood. Here, we review the ongoing development of accessibility measurements, we summarize the different molecular and structural mechanisms that shape the accessibility landscape, and we detail the many important biological functions that are linked to chromatin accessibility.

In Vitro Mapping of Nucleosome Positions at Base-Pair Resolution Using Ortho-Phenanthroline

The positions of nucleosomes along genomic DNA play a role in defining patterns of gene expression and chromatin organization. Determination of nucleosome positions in vivo and in vitro, as revealed by the locations of histones on DNA, has provided insight into mechanisms of nucleosome sliding, spacing, assembly, and disassembly. Here, we describe methods for the in vitro determination of histone-DNA contacts at base-pair (bp) resolution. The protocol involves the labeling of histones with ortho-phenanthroline (OP), site-specific cleavage of nucleosomal DNA, and processing and analysis of the resulting DNA fragments. This methodology provides an efficient and high-resolution means for studying kinetics and behavior of enzymes that alter nucleosome structure and/or positioning, and can be used to identify preferred distributions of nucleosomes on natural DNA sequences. © 2022 Wiley Periodicals LLC. Basic Protocol 1: Cysteine-specific chemical modification of folded histones with ortho-phenanthroline (OP) Basic Protocol 2: Nucleosome sliding assay adapted for OP mapping of histone-DNA contacts Basic Protocol 3: OP-mediated cleavage, processing, and analysis of DNA fragments using a sequencing gel Support Protocol 1: Preparation of dideoxy sequencing ladders Support Protocol 2: Preparation and running of a denaturing DNA sequencing gel.

GC content strongly influences the role of poly(dA) in the intrinsic nucleosome positioning in Saccharomyces cerevisiae

The nucleosome is the basic structural element of genomic DNA packaging and plays a role in transcription, replication, and recombination. Poly(dA) tracts are considered major sequence determinants of nucleosome positioning, although their role is not well understood. Here, we show that the homopolymeric character and the low GC content of poly(dA)s play different roles in nucleosome formation. We found that the inherent low GC content of poly(dA) alone can account for the deep and anisotropic nucleosome depletion at structurally and functionally important regions of promoters and origins of replication. We also show that the level of nucleosome occupancy at poly(dA) is strongly related to the local nucleotide background and its high frequency of occurrence in Saccharomyces cerevisiae does not appear merely to be associated with its intrinsic nucleosome-excluding properties. In addition, we show that the GC content alone can predict more than 60% of the in vitro nucleosome map, providing further evidence that the intrinsic nucleosome positioning is more greatly determined by GC content than poly(dA) stretches. Our results are consistent with a model in which poly(dA) stretches act at two distinct levels: first, by its low GC content, which intrinsically contributes to hinder nucleosome formation, and second, by its contiguous runs of dA that selectively drive the recruitment of non-histone proteins with structural and functional roles.

Proximal-end bias from in-vitro reconstituted nucleosomes and the result on downstream data analysis

The most basic level of eukaryotic gene regulation is the presence or absence of nucleosomes on DNA regulatory elements. In an effort to elucidate in vivo nucleosome patterns, in vitro studies are frequently used. In vitro, short DNA fragments are more favorable for nucleosome formation, increasing the likelihood of nucleosome occupancy. This may in part result from the fact that nucleosomes prefer to form on the terminal ends of linear DNA. This phenomenon has the potential to bias in vitro reconstituted nucleosomes and skew results. If the ends of DNA fragments are known, the reads falling close to the ends are typically discarded. In this study we confirm the phenomenon of end bias of in vitro nucleosomes. We describe a method in which nearly identical libraries, with different known ends, are used to recover nucleosomes which form towards the terminal ends of fragmented DNA. Finally, we illustrate that although nucleosomes prefer to form on DNA ends, it does not appear to skew results or the interpretation thereof.

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.

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.

An Essential and Cell-Cycle-Dependent ORC Dimerization Cycle Regulates Eukaryotic Chromosomal DNA Replication

Eukaryotic DNA replication licensing is a prerequisite for, and plays a role in, regulating genome duplication that occurs exactly once per cell cycle. ORC (origin recognition complex) binds to and marks replication origins throughout the cell cycle and loads other replication-initiation proteins onto replication origins to form pre-replicative complexes (pre-RCs), completing replication licensing. However, how an asymmetric single-heterohexameric ORC structure loads the symmetric MCM (minichromosome maintenance) double hexamers is controversial, and importantly, it remains unknown when and how ORC proteins associate with the newly replicated origins to protect them from invasion by histones. Here, we report an essential and cell-cycle-dependent ORC "dimerization cycle" that plays three fundamental roles in the regulation of DNA replication: providing a symmetric platform to load the symmetric pre-RCs, marking and protecting the nascent sister replication origins for the next licensing, and playing a crucial role to prevent origin re-licensing within the same cell cycle.

DNA mechanics and its biological impact

Almost all nucleoprotein interactions and DNA manipulation events involve mechanical deformations of DNA. Extraordinary progresses in single-molecule, structural, and computational methods have characterized the average mechanical properties of DNA, such as bendability and torsional rigidity, in high resolution. Further, the advent of sequencing technology has permitted measuring, in high-throughput, how such mechanical properties vary with sequence and epigenetic modifications along genomes. We review these recent technological advancements, and discuss how they have contributed to the emerging idea that variations in the mechanical properties of DNA play a fundamental role in regulating, genome-wide, diverse processes involved in chromatin organization.

A deformation energy model reveals sequence-dependent property of nucleosome positioning

We present a deformation energy model for predicting nucleosome positioning, in which a position-dependent structural parameter set derived from crystal structures of nucleosomes was used to calculate the DNA deformation energy. The model is successful in predicting nucleosome occupancy genome-wide in budding yeast, nucleosome free energy, and rotational positioning of nucleosomes. Our model also indicates that the genomic regions underlying the MNase-sensitive nucleosomes in budding yeast have high deformation energy and, consequently, low nucleosome-forming ability, while the MNase-sensitive non-histone particles are characterized by much lower DNA deformation energy and high nucleosome preference. In addition, we also revealed that remodelers, SNF2 and RSC8, are likely to act in chromatin remodeling by binding to broad nucleosome-depleted regions that are intrinsically favorable for nucleosome positioning. Our data support the important role of position-dependent physical properties of DNA in nucleosome positioning.

Nucleosome Positioning and Spacing: From Mechanism to Function

Eukaryotes associate their genomes with histone proteins, forming nucleosome particles. Nucleosomes regulate and protect the genetic information. They often assemble into evenly spaced arrays of nucleosomes. These regular nucleosome arrays cover significant portions of the genome, in particular over genes. The presence of these evenly spaced nucleosome arrays is highly conserved throughout the entire eukaryotic domain. Here, we review the mechanisms behind the establishment of this primary structure of chromatin with special emphasis on the biogenesis of evenly spaced nucleosome arrays. We highlight the roles that transcription, nucleosome remodelers, DNA sequence, and histone density play towards the formation of evenly spaced nucleosome arrays and summarize our current understanding of their cellular functions. We end with key unanswered questions that remain to be explored to obtain an in-depth understanding of the biogenesis and function of the nucleosome landscape.

Prospect of reprogramming replication licensing for cancer drug development

Eukaryotic chromosomal DNA replication is preceded by replication licensing which involves the identification of the origin of replication by origin recognition complex (ORC). The ORC loads pre-replication complexes (pre-RCs) through a series of tightly regulated mechanisms where the ORC interacts with Cdc6 to recruit cdt1-MCM2-7 complexes to the origin of replication. In eukaryotes, adherence to regulatory mechanisms of the replication program is required to ensure that all daughter cells carry the exact copy of genetic material as the parent cell. Failure of which leads to the development of genome instability syndromes like cancer, diabetes, etc. In an event of such occurrence, preventing cells from carrying the defaulted genetic material and passing it to other cells hinges on the regulation of chromosomal DNA replication. Thus, understanding the mechanisms underpinning chromosomal DNA replication and particularly replication licensing can expose druggable enzymes, effector molecules, and secondary messengers that can be targeted for diagnosis and therapeutic purposes. Effectively drugging these molecular markers to reprogram pre-replication events can be used to control the fate of chromosomal DNA replication for the treatment of genome instability disorders and in this case, cancer. This review discusses available knowledge of replication licensing in the contest of molecular drug discovery for the treatment of cancer.

The origin and evolution of a distinct mechanism of transcription initiation in yeasts

The molecular process of transcription by RNA Polymerase II is highly conserved among eukaryotes (“classic model”). A distinct way of locating transcription start sites (TSSs) has been identified in a budding yeast Saccharomyces cerevisiae (“scanning model”). Herein, we applied genomic approaches to elucidate the origin of the scanning model and its underlying genetic mechanisms. We first identified TSSs at single-nucleotide resolution for 12 yeast species using the nAnT-iCAGE technique, which significantly improved the annotations of these genomes by providing accurate 5′ boundaries for protein-coding genes. We then inferred the initiation mechanism of each species based on its TSS maps and genome sequences. We discovered that the scanning model likely originated after the split of Yarrowia lipolytica and the other budding yeasts. Species that use the scanning model showed an adenine-rich region immediately upstream of the TSS that might facilitate TSS selection. Both initiation mechanisms share a strong preference for pyrimidine–purine dinucleotides surrounding the TSS. Our results suggest that the purine is required to accurately recruit the first nucleotide, thereby increasing the chances of a messenger RNA of being capped during mRNA maturation, which is critical for efficient translation initiation during protein biosynthesis. Based on our findings, we propose a model for TSS selection in the scanning-model species, as well as a model for the stepwise process responsible for the origin and evolution of the scanning model.

Nucleosome Clutches are Regulated by Chromatin Internal Parameters

Nucleosomes cluster together when chromatin folds in the cell to form heterogeneous groups termed "clutches". These structural units add another level of chromatin regulation, for example during cell differentiation. Yet, the mechanisms that regulate their size and compaction remain obscure. Here, using our chromatin mesoscale model, we dissect clutch patterns in fibers with different combinations of nucleosome positions, linker histone density, and acetylation levels to investigate their role in clutch regulation. First, we isolate the effect of each chromatin parameter by studying systems with regular nucleosome spacing; second, we design systems with naturally-occurring linker lengths that fold onto specific clutch patterns; third, we model gene-encoding fibers to understand how these combined factors contribute to gene structure. Our results show how these chromatin parameters act together to produce different-sized nucleosome clutches. The length of nucleosome free regions (NFRs) profoundly affects clutch size, while the length of linker DNA has a moderate effect. In general, higher linker histone densities produce larger clutches by a chromatin compaction mechanism, while higher acetylation levels produce smaller clutches by a chromatin unfolding mechanism. We also show that it is possible to design fibers with naturally-occurring DNA linkers and NFRs that fold onto specific clutch patterns. Finally, in gene-encoding systems, a complex combination of variables dictates a gene-specific clutch pattern. Together, these results shed light into the mechanisms that regulate nucleosome clutches and suggest a new epigenetic mechanism by which chromatin parameters regulate transcriptional activity via the three-dimensional folded state of the genome at a nucleosome level.

Bioinformatical dissection of fission yeast DNA replication origins

Replication origins in eukaryotes form a base for assembly of the pre-replication complex (pre-RC), thereby serving as an initiation site of DNA replication. Characteristics of replication origin vary among species. In fission yeast Schizosaccharomyces pombe, DNA of high AT content is a distinct feature of replication origins; however, it remains to be understood what the general molecular architecture of fission yeast origin is. Here, we performed ChIP-seq mapping of Orc4 and Mcm2, two representative components of the pre-RC, and described the characteristics of their binding sites. The analysis revealed that fission yeast efficient origins are associated with two similar but independent features: a ≥15 bp-long motif with stretches of As and an AT-rich region of a few hundred bp. The A-rich motif was correlated with chromosomal binding of Orc, a DNA-binding component in the pre-RC, whereas the AT-rich region was associated with efficient binding of the DNA replicative helicase Mcm. These two features, in combination with the third feature, a transcription-poor region of approximately 1 kb, enabled to distinguish efficient replication origins from the rest of chromosome arms with high accuracy. This study, hence, provides a model that describes how multiple functional elements specify DNA replication origins in fission yeast genome.

Nucleosome Dynamics: a new tool for the dynamic analysis of nucleosome positioning

We present Nucleosome Dynamics, a suite of programs integrated into a virtual research environment and created to define nucleosome architecture and dynamics from noisy experimental data. The package allows both the definition of nucleosome architectures and the detection of changes in nucleosomal organization due to changes in cellular conditions. Results are displayed in the context of genomic information thanks to different visualizers and browsers, allowing the user a holistic, multidimensional view of the genome/transcriptome. The package shows good performance for both locating equilibrium nucleosome architecture and nucleosome dynamics and provides abundant useful information in several test cases, where experimental data on nucleosome position (and for some cases expression level) have been collected for cells under different external conditions (cell cycle phase, yeast metabolic cycle progression, changes in nutrients or difference in MNase digestion level). Nucleosome Dynamics is a free software and is provided under several distribution models.

Absolute nucleosome occupancy map for the Saccharomyces cerevisiae genome

Mapping of nucleosomes, the basic DNA packaging unit in eukaryotes, is fundamental for understanding genome regulation because nucleosomes modulate DNA access by their positioning along the genome. A cell-population nucleosome map requires two observables: nucleosome positions along the DNA ("Where?") and nucleosome occupancies across the population ("In how many cells?"). All available genome-wide nucleosome mapping techniques are yield methods because they score either nucleosomal (e.g., MNase-seq, chemical cleavage-seq) or nonnucleosomal (e.g., ATAC-seq) DNA but lose track of the total DNA population for each genomic region. Therefore, they only provide nucleosome positions and maybe compare relative occupancies between positions, but cannot measure absolute nucleosome occupancy, which is the fraction of all DNA molecules occupied at a given position and time by a nucleosome. Here, we established two orthogonal and thereby cross-validating approaches to measure absolute nucleosome occupancy across the Saccharomyces cerevisiae genome via restriction enzymes and DNA methyltransferases. The resulting high-resolution (9-bp) map shows uniform absolute occupancies. Most nucleosome positions are occupied in most cells: 97% of all nucleosomes called by chemical cleavage-seq have a mean absolute occupancy of 90 ± 6% (±SD). Depending on nucleosome position calling procedures, there are 57,000 to 60,000 nucleosomes per yeast cell. The few low absolute occupancy nucleosomes do not correlate with highly transcribed gene bodies, but correlate with increased presence of the nucleosome-evicting chromatin structure remodeling (RSC) complex, and are enriched upstream of highly transcribed or regulated genes. Our work provides a quantitative method and reference frame in absolute terms for future chromatin studies.

The nucleosome position-encoding WW/SS sequence pattern is depleted in mammalian genes relative to other eukaryotes

Nucleosomal DNA sequences generally follow a well-known pattern with ∼10-bp periodic WW (where W is A or T) dinucleotides that oscillate in phase with each other and out of phase with SS (where S is G or C) dinucleotides. However, nucleosomes with other DNA patterns have not been systematically analyzed. Here, we focus on an opposite pattern, namely anti-WW/SS pattern, in which WW dinucleotides preferentially occur at DNA sites that bend into major grooves and SS (where S is G or C) dinucleotides are often found at sites that bend into minor grooves. Nucleosomes with the anti-WW/SS pattern are widespread and exhibit a species- and context-specific distribution in eukaryotic genomes. Unlike non-mammals (yeast, nematode and fly), there is a positive correlation between the enrichment of anti-WW/SS nucleosomes and RNA Pol II transcriptional levels in mammals (mouse and human). Interestingly, such enrichment is not due to underlying DNA sequence. In addition, chromatin remodeling complexes have an impact on the abundance but not on the distribution of anti-WW/SS nucleosomes in yeast. Our data reveal distinct roles of cis- and trans-acting factors in the rotational positioning of nucleosomes between non-mammals and mammals. Implications of the anti-WW/SS sequence pattern for RNA Pol II transcription are discussed.

Biochemical analysis of nucleosome targeting by Tn5 transposase

Tn5 transposase is a bacterial enzyme that integrates a DNA fragment into genomic DNA, and is used as a tool for detecting nucleosome-free regions of genomic DNA in eukaryotes. However, in chromatin, the DNA targeting by Tn5 transposase has remained unclear. In the present study, we reconstituted well-positioned 601 dinucleosomes, in which two nucleosomes are connected with a linker DNA, and studied the DNA integration sites in the dinucleosomes by Tn5 transposase in vitro. We found that Tn5 transposase preferentially targets near the entry-exit DNA regions within the nucleosome. Tn5 transposase minimally cleaved the dinucleosome without a linker DNA, indicating that the linker DNA between two nucleosomes is important for the Tn5 transposase activity. In the presence of a 30 base-pair linker DNA, Tn5 transposase targets the middle of the linker DNA, in addition to the entry-exit sites of the nucleosome. Intriguingly, this Tn5-targeting characteristic is conserved in a dinucleosome substrate with a different DNA sequence from the 601 sequence. Therefore, the Tn5-targeting preference in the nucleosomal templates reported here provides important information for the interpretation of Tn5 transposase-based genomics methods, such as ATAC-seq.

Bridging chromatin structure and function over a range of experimental spatial and temporal scales by molecular modeling

Chromatin structure, dynamics, and function are being intensely investigated by a variety of methods, including microscopy, X-ray diffraction, nuclear magnetic resonance, biochemical crosslinking, chromosome conformation capture, and computation. A range of experimental techniques combined with modeling is clearly valuable to help interpret experimental data and, importantly, generate configurations and mechanisms related to the 3D organization and function of the genome. Contact maps, in particular, as obtained by a variety of chromosome conformation capture methods, are of increasing interest due to their implications on genome structure and regulation on many levels. In this perspective, using seven examples from our group's studies, we illustrate how molecular modeling can help interpret such experimental data. Specifically, we show how computed contact maps related to experimental systems can be used to explain structures of nucleosomes, chromatin higher-order folding, domain segregation mechanisms, gene organization, and the effect on chromatin structure of external and internal fiber parameters, such as nucleosome positioning, presence of nucleosome free regions, histone posttranslational modifications, and linker histone binding. We argue that such computations on multiple spatial and temporal scales will be increasingly important for the integration of genomic, epigenomic, and biophysical data on chromatin structure and related cellular processes.

Chromatin accessibility and the regulatory epigenome

Physical access to DNA is a highly dynamic property of chromatin that plays an essential role in establishing and maintaining cellular identity. The organization of accessible chromatin across the genome reflects a network of permissible physical interactions through which enhancers, promoters, insulators and chromatin-binding factors cooperatively regulate gene expression. This landscape of accessibility changes dynamically in response to both external stimuli and developmental cues, and emerging evidence suggests that homeostatic maintenance of accessibility is itself dynamically regulated through a competitive interplay between chromatin-binding factors and nucleosomes. In this Review, we examine how the accessible genome is measured and explore the role of transcription factors in initiating accessibility remodelling; our goal is to illustrate how chromatin accessibility defines regulatory elements within the genome and how these epigenetic features are dynamically established to control gene expression.

Genome-wide search of nucleosome patterns using visual analytics

MOTIVATION

The Burrows-Wheeler transform (BWT) is widely used for the fast alignment of high-throughput sequence data. This method also has potential applications in other areas of bioinformatics, and it can be specially useful for the fast searching of patterns on coverage data from different sources.

RESULTS

We present a nucleosome pattern search method that converts levels of nucleosomal occupancy to a sequence-like format to which BWT searches can be applied. The method is embedded in a nucleosome map browser, 'Nucleosee', an interactive visual tool specifically designed to enhance BWT searches, giving them context and making them suitable for visual discourse analysis of the results. The proposed method is fast, flexible and sufficiently generic for the exploration of data in a broad and interactive way.

AVAILABILITY AND IMPLEMENTATION

The proposed algorithm and visual browser are available for testing at http://cpg3.der.usal.es/nucleosee. The source code and installation packages are also available at https://github.com/rodrigoSantamaria/nucleosee.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

Comprehensive profiling of the fission yeast transcription start site activity during stress and media response

Fission yeast, Schizosaccharomyces pombe, is an attractive model organism for transcriptional and chromatin biology research. Such research is contingent on accurate annotation of transcription start sites (TSSs). However, comprehensive genome-wide maps of TSSs and their usage across commonly applied laboratory conditions and treatments for S. pombe are lacking. To this end, we profiled TSS activity genome-wide in S. pombe cultures exposed to heat shock, nitrogen starvation, hydrogen peroxide and two commonly applied media, YES and EMM2, using Cap Analysis of Gene Expression (CAGE). CAGE-based annotation of TSSs is substantially more accurate than existing PomBase annotation; on average, CAGE TSSs fall 50-75 bp downstream of PomBase TSSs and co-localize with nucleosome boundaries. In contrast to higher eukaryotes, dispersed TSS distributions are not common in S. pombe. Our data recapitulate known S. pombe stress expression response patterns and identify stress- and media-responsive alternative TSSs. Notably, alteration of growth medium induces changes of similar magnitude as some stressors. We show a link between nucleosome occupancy and genetic variation, and that the proximal promoter region is genetically diverse between S. pombe strains. Our detailed TSS map constitutes a central resource for S. pombe gene regulation research.

A unified computational framework for modeling genome-wide nucleosome landscape

Nucleosomes form the fundamental building blocks of eukaryotic chromatin, and previous attempts to understand the principles governing their genome-wide distribution have spurred much interest and debate in biology. In particular, the precise role of DNA sequence in shaping local chromatin structure has been controversial. This paper rigorously quantifies the contribution of hitherto-debated sequence features-including G+C content, 10.5 bp periodicity, and poly(dA:dT) tracts-to three distinct aspects of genome-wide nucleosome landscape: occupancy, translational positioning and rotational positioning. Our computational framework simultaneously learns nucleosome number and nucleosome-positioning energy from genome-wide nucleosome maps. In contrast to other previous studies, our model can predict both in vitro and in vivo nucleosome maps in Saccharomyces cerevisiae. We find that although G+C content is the primary determinant of MNase-derived nucleosome occupancy, MNase digestion biases may substantially influence this GC dependence. By contrast, poly(dA:dT) tracts are seen to deter nucleosome formation, regardless of the experimental method used. We further show that the 10.5 bp nucleotide periodicity facilitates rotational but not translational positioning. Applying our method to in vivo nucleosome maps demonstrates that, for a subset of genes, the regularly-spaced nucleosome arrays observed around transcription start sites can be partially recapitulated by DNA sequence alone. Finally, in vivo nucleosome occupancy derived from MNase-seq experiments around transcription termination sites can be mostly explained by the genomic sequence. Implications of these results and potential extensions of the proposed computational framework are discussed.

DMS-Seq for In Vivo Genome-wide Mapping of Protein-DNA Interactions and Nucleosome Centers

Protein-DNA interactions provide the basis for chromatin structure and gene regulation. Comprehensive identification of protein-occupied sites is thus vital to an in-depth understanding of genome function. Dimethyl sulfate (DMS) is a chemical probe that has long been used to detect footprints of DNA-bound proteins in vitro and in vivo. Here, we describe a genomic footprinting method, dimethyl sulfate sequencing (DMS-seq), which exploits the cell-permeable nature of DMS to obviate the need for nuclear isolation. This feature makes DMS-seq simple in practice and removes the potential risk of protein re-localization during nuclear isolation. DMS-seq successfully detects transcription factors bound to cis-regulatory elements and non-canonical chromatin particles in nucleosome-free regions. Furthermore, an unexpected preference of DMS confers on DMS-seq a unique potential to directly detect nucleosome centers without using genetic manipulation. We expect that DMS-seq will serve as a characteristic method for genome-wide interrogation of in vivo protein-DNA interactions.

Precise genome-wide mapping of single nucleosomes and linkers in vivo

We developed a chemical cleavage method that releases single nucleosome dyad-containing fragments, allowing us to precisely map both single nucleosomes and linkers with high accuracy genome-wide in yeast. Our single nucleosome positioning data reveal that nucleosomes occupy preferred positions that differ by integral multiples of the DNA helical repeat. By comparing nucleosome dyad positioning maps to existing genomic and transcriptomic data, we evaluated the contributions of sequence, transcription, and histones H1 and H2A.Z in defining the chromatin landscape. We present a biophysical model that neglects DNA sequence and shows that steric occlusion suffices to explain the salient features of nucleosome positioning.

Parallel mapping with site-directed hydroxyl radicals and micrococcal nuclease reveals structural features of positioned nucleosomes in vivo

Micrococcal nuclease (MNase) has been widely used for analyses of nucleosome locations in many organisms. However, due to its sequence preference, the interpretations of the positions and occupancies of nucleosomes using MNase have remained controversial. Next-generation sequencing (NGS) has also been utilized for analyses of MNase-digests, but some technical biases are commonly present in the NGS experiments. Here, we established a gel-based method to map nucleosome positions in Saccharomyces cerevisiae, using isolated nuclei as the substrate for the histone H4 S47C-site-directed chemical cleavage in parallel with MNase digestion. The parallel mapping allowed us to compare the chemically and enzymatically cleaved sites by indirect end-labeling and primer extension mapping, and thus we could determine the nucleosome positions and the sizes of the nucleosome-free regions (or nucleosome-depleted regions) more accurately, as compared to nucleosome mapping by MNase alone. The analysis also revealed that the structural features of the nucleosomes flanked by the nucleosome-free region were different from those within regularly arrayed nucleosomes, showing that the structures and dynamics of individual nucleosomes strongly depend on their locations. Moreover, we demonstrated that the parallel mapping results were generally consistent with the previous genome-wide chemical mapping and MNase-Seq results. Thus, the gel-based parallel mapping will be useful for the analysis of a specific locus under various conditions.

A Reply to "MNase-Sensitive Complexes in Yeast: Nucleosomes and Non-histone Barriers," by Chereji et al

In this issue of Molecular Cell, Chereji et al. (2017) present new data on MNase-sensitive particles previously identified upstream of transcription start sites at many promoters in budding yeast, and they argue, based upon negative histone-ChIP results, that they are non-nucleosomal signals generated by transcription factors (TFs). We show instead, based upon functional experiments where the relevant TFs are rapidly depleted, that this explanation does not hold, and we argue instead that histone ChIP and chemical cleavage assays have a limited capacity to capture these highly dynamic, MNase-sensitive "fragile" nucleosomes.

Evolutionary analysis of nucleosome positioning sequences based on New Symmetric Relative Entropy

New Symmetric Relative Entropy (NSRE) was applied innovatively to analyze the nucleosome sequences in S. cerevisiae, S. pombe and Drosophila. NSRE distributions could well reflect the characteristic differences of nucleosome sequences among three organisms, and the differences indicate a concerted evolution in the sequence usage of nucleosome. Further analysis about the nucleosomes around TSS shows that the constitutive property of +1/-1 nucleosomes in S. cerevisiae is different from that in S. pombe and Drosophila, which indicates that S. cerevisiae has a different transcription regulation mechanism based on nucleosome. However, in either case, the nucleosome dyad region is conserved and always has a higher NSRE. Base composition analysis shows that this conservative property in nucleosome dyad region is mainly determined by base A and T, and the dependence degrees on base A and T are consistent in three organisms.

Capturing Structural Heterogeneity in Chromatin Fibers

Chromatin fiber organization is implicated in processes such as transcription, DNA repair and chromosome segregation, but how nucleosomes interact to form higher-order structure remains poorly understood. We solved two crystal structures of tetranucleosomes with approximately 11-bp DNA linker length at 5.8 and 6.7 Å resolution. Minimal intramolecular nucleosome-nucleosome interactions result in a fiber model resembling a flat ribbon that is compatible with a two-start helical architecture, and that exposes histone and DNA surfaces to the environment. The differences in the two structures combined with electron microscopy reveal heterogeneous structural states, and we used site-specific chemical crosslinking to assess the diversity of nucleosome-nucleosome interactions through identification of structure-sensitive crosslink sites that provide a means to characterize fibers in solution. The chromatin fiber architectures observed here provide a basis for understanding heterogeneous chromatin higher-order structures as they occur in a genomic context.

Nucleosomal organization and DNA base composition patterns

Nucleosomes are the basic units of chromatin. They compact the genome inside the nucleus and regulate the access of proteins to DNA. In the yeast genome, most nucleosomes occupy well-defined positions, which are maintained under many different physiological situations and genetic backgrounds. Although several short sequence elements have been described that favor or reduce the affinity between histones and DNA, the extent to which the DNA sequence affects nucleosome positioning in the genomic context remains unclear. Recent analyses indicate that the base composition pattern of mononucleosomal DNA differs among species, and that the same sequence elements have a different impact on nucleosome positioning in different genomes despite the high level of phylogenetic conservation of histones. These studies have also shown that the DNA sequence contributes to nucleosome positioning to the point that it is possible to design synthetic DNA molecules capable of generating regular and species-specific nucleosomal patterns in vivo.

Towards quantitative analysis of gene regulation by enhancers

Enhancers are regulatory DNA sequences that can activate transcription over large distances. Recent studies have revealed the widespread role of distant activation in eukaryotic gene regulation and in the development of various human diseases, including cancer. Here we review recent progress in the field, focusing on new experimental and computational approaches that quantify the role of chromatin structure and dynamics during enhancer-promoter interactions in vitro and in vivo.

Genome-wide Mapping of the Nucleosome Landscape by Micrococcal Nuclease and Chemical Mapping

Nucleosomes regulate the transcription output of the genome by occluding the underlying DNA sequences from DNA-binding proteins that must act on it. Knowledge of the precise locations of nucleosomes in the genome is thus essential towards understanding how transcription is regulated. Current nucleosome-mapping strategies involve digesting chromatin with nucleases or chemical cleavage followed by high-throughput sequencing. In this review, we compare the traditional micrococcal nuclease (MNase)-based approach with a chemical cleavage strategy, with discussion on the important insights each has uncovered about the role of nucleosomes in shaping transcriptional processes.

A comparison of nucleosome organization in Drosophila cell lines

Changes in the distribution of nucleosomes along the genome influence chromatin structure and impact gene expression by modulating the accessibility of DNA to transcriptional machinery. However, the role of genome-wide nucleosome positioning in gene expression and in maintaining differentiated cell states remains poorly understood. Drosophila melanogaster cell lines represent distinct tissue types and exhibit cell-type specific gene expression profiles. They thus could provide a useful tool for investigating cell-type specific nucleosome organization of an organism's genome. To evaluate this possibility, we compared genome-wide nucleosome positioning and occupancy in five different Drosophila tissue-specific cell lines, and in reconstituted chromatin, and then tested for correlations between nucleosome positioning, transcription factor binding motifs, and gene expression. Nucleosomes in all cell lines were positioned in accordance with previously known DNA-nucleosome interactions, with helically repeating A/T di-nucleotide pairs arranged within nucleosomal DNAs and AT-rich pentamers generally excluded from nucleosomal DNA. Nucleosome organization in all cell lines differed markedly from in vitro reconstituted chromatin, with highly expressed genes showing strong nucleosome organization around transcriptional start sites. Importantly, comparative analysis identified genomic regions that exhibited cell line-specific nucleosome enrichment or depletion. Further analysis of these regions identified 91 out of 16,384 possible heptamer sequences that showed differential nucleosomal occupation between cell lines, and 49 of the heptamers matched one or more known transcription factor binding sites. These results demonstrate that there is differential nucleosome positioning between these Drosophila cell lines and therefore identify a system that could be used to investigate the functional significance of differential nucleosomal positioning in cell type specification.

Kilobase Pair Chromatin Fiber Contacts Promoted by Living-System-Like DNA Linker Length Distributions and Nucleosome Depletion

Nucleosome placement, or DNA linker length patterns, are believed to yield specific spatial features in chromatin fibers, but details are unknown. Here we examine by mesoscale modeling how kilobase (kb) range contacts and fiber looping depend on linker lengths ranging from 18 to 45 bp, with values modeled after living systems, including nucleosome free regions (NFRs) and gene encoding segments. We also compare artificial constructs with alternating versus randomly distributed linker lengths in the range of 18-72 bp. We show that nonuniform distributions with NFRs enhance flexibility and encourage kb-range contacts. NFRs between neighboring gene segments diminish short-range contacts between flanking nucleosomes, while enhancing kb-range contacts via hierarchical looping. We also demonstrate that variances in linker lengths enhance such contacts. In particular, moderate sized variations in fiber linker lengths (∼27 bp) encourage long-range contacts in randomly distributed linker length fibers. Our work underscores the importance of linker length patterns, alongside bound proteins, in biological regulation. Contacts formed by kb-range chromatin folding are crucial to gene activity. Because we find that special linker length distributions in living systems promote kb contacts, our work suggests ways to manipulate these patterns for regulation of gene activity.

Establishing nucleosome architecture and stability at promoters: Roles of pioneer transcription factors and the RSC chromatin remodeler

Improvements in deep sequencing, together with methods to rapidly deplete essential transcription factors (TFs) and chromatin remodelers, have recently led to a more detailed picture of promoter nucleosome architecture in yeast and its relationship to transcriptional regulation. These studies revealed that ∼40% of all budding yeast protein-coding genes possess a unique promoter structure, where we propose that an unusually unstable nucleosome forms immediately upstream of the transcription start site (TSS). This "fragile" nucleosome (FN) promoter architecture relies on the combined action of the essential RSC (Remodels Structure of Chromatin) nucleosome remodeler and pioneer transcription factors (PTFs). FNs are associated with genes whose expression is high, coupled to cell growth, and characterized by low cell-to-cell variability (noise), suggesting that they may promote these features. Recent studies in metazoans suggest that the presence of dynamic nucleosomes upstream of the TSS at highly expressed genes may be conserved throughout evolution.

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.

In vitro reconstitution and biochemical analyses of the Schizosaccharomyces pombe nucleosome

Schizosaccharomyces pombe, which has a small genome but shares many physiological functions with higher eukaryotes, is a useful single-cell, model eukaryotic organism. In particular, many features concerning chromatin structure and dynamics, including heterochromatin, centromeres, telomeres, and DNA replication origins, are well conserved between S. pombe and higher eukaryotes. However, the S. pombe nucleosome, the fundamental structural unit of chromatin, has not been reconstituted in vitro. In the present study, we established the method to purify S. pombe histones H2A, H2B, H3, and H4, and successfully reconstituted the S. pombe nucleosome in vitro. Our thermal stability assay and micrococcal nuclease treatment assay revealed that the S. pombe nucleosome is markedly unstable and its DNA ends are quite accessible, as compared to the canonical human nucleosome. These findings are important to understand the mechanisms of epigenetic genomic DNA regulation in fission yeast.

Insights into Nucleosome Organization in Mouse Embryonic Stem Cells through Chemical Mapping

Nucleosome organization influences gene activity by controlling DNA accessibility to transcription machinery. Here, we develop a chemical biology approach to determine mammalian nucleosome positions genome-wide. We uncovered surprising features of nucleosome organization in mouse embryonic stem cells. In contrast to the prevailing model, we observe that for nearly all mouse genes, a class of fragile nucleosomes occupies previously designated nucleosome-depleted regions around transcription start sites and transcription termination sites. We show that nucleosomes occupy DNA targets for a subset of DNA-binding proteins, including CCCTC-binding factor (CTCF) and pluripotency factors. Furthermore, we provide evidence that promoter-proximal nucleosomes, with the +1 nucleosome in particular, contribute to the pausing of RNA polymerase II. Lastly, we find a characteristic preference for nucleosomes at exon-intron junctions. Taken together, we establish an accurate method for defining the nucleosome landscape and provide a valuable resource for studying nucleosome-mediated gene regulation in mammalian cells.