Major Determinants of Nucleosome Positioning

Interpretable deep residual network uncovers nucleosome positioning and associated features

Nucleosomes represent elementary building units of eukaryotic chromosomes and consist of DNA wrapped around a histone octamer flanked by linker DNA segments. Nucleosomes are central in epigenetic pathways and their genomic positioning is associated with regulation of gene expression, DNA replication, DNA methylation and DNA repair, among other functions. Building on prior discoveries that DNA sequences noticeably affect nucleosome positioning, our objective is to identify nucleosome positions and related features across entire genome. Here, we introduce an interpretable framework based on the concepts of deep residual networks (NuPoSe). Trained on high-coverage human experimental MNase-seq data, NuPoSe is able to learn sequence and structural patterns associated with nucleosome organization in human genome. NuPoSe can be also applied to unseen data from different organisms and cell types. Our findings point to 43 informative features, most of them constitute tri-nucleotides, di-nucleotides and one tetra-nucleotide. Most features are significantly associated with the nucleosomal structural characteristics, namely, periodicity of nucleosomal DNA and its location with respect to a histone octamer. Importantly, we show that features derived from the 27 bp linker DNA flanking nucleosomes contribute up to 10% to the quality of the prediction model. This, along with the comprehensive training sets, deep-learning architecture, and feature selection method, may contribute to the NuPoSe's 80-89% classification accuracy on different independent datasets.

nucMACC: An MNase-seq pipeline to identify structurally altered nucleosomes in the genome

Micrococcal nuclease sequencing is the state-of-the-art method for determining chromatin structure and nucleosome positioning. Data analysis is complex due to the AT-dependent sequence bias of the endonuclease and the requirement for high sequencing depth. Here, we present the nucleosome-based MNase accessibility (nucMACC) pipeline unveiling the regulatory chromatin landscape by measuring nucleosome accessibility and stability. The nucMACC pipeline represents a systematic and genome-wide approach for detecting unstable ("fragile") nucleosomes. We have characterized the regulatory nucleosome landscape in Drosophila melanogaster, Saccharomyces cerevisiae, and mammals. Two functionally distinct sets of promoters were characterized, one associated with an unstable nucleosome and the other being nucleosome depleted. We show that unstable nucleosomes present intermediate states of nucleosome remodeling, preparing inducible genes for transcriptional activation in response to stimuli or stress. The presence of unstable nucleosomes correlates with RNA polymerase II proximal pausing. The nucMACC pipeline offers unparalleled precision and depth in nucleosome research and is a valuable tool for future nucleosome studies.

Examining chromatin heterogeneity through PacBio long-read sequencing of M.EcoGII methylated genomes: an m6A detection efficiency and calling bias correcting pipeline

Recent studies have combined DNA methyltransferase footprinting of genomic DNA in nuclei with long-read sequencing, resulting in detailed chromatin maps for multi-kilobase stretches of genomic DNA from one cell. Theoretically, nucleosome footprints and nucleosome-depleted regions can be identified using M.EcoGII, which methylates adenines in any sequence context, providing a high-resolution map of accessible regions in each DNA molecule. Here, we report PacBio long-read sequence data for budding yeast nuclei treated with M.EcoGII and a bioinformatic pipeline which corrects for three key challenges undermining this promising method. First, detection of m6A in individual DNA molecules by the PacBio software is inefficient, resulting in false footprints predicted by random gaps of seemingly unmethylated adenines. Second, there is a strong bias against m6A base calling as AT content increases. Third, occasional methylation occurs within nucleosomes, breaking up their footprints. After correcting for these issues, our pipeline calculates a correlation coefficient-based score indicating the extent of chromatin heterogeneity within the cell population for every gene. Although the population average is consistent with that derived using other techniques, we observe a wide range of heterogeneity in nucleosome positions at the single-molecule level, probably reflecting cellular chromatin dynamics.

General regulatory factors exert differential effects on nucleosome sliding activity of the ISW1a complex

BACKGROUND

Chromatin dynamics is deeply involved in processes that require access to DNA, such as transcriptional regulation. Among the factors involved in chromatin dynamics at gene regulatory regions are general regulatory factors (GRFs). These factors contribute to establishment and maintenance of nucleosome-depleted regions (NDRs). These regions are populated by nucleosomes through histone deposition and nucleosome sliding, the latter catalyzed by a number of ATP-dependent chromatin remodeling complexes, including ISW1a. It has been observed that GRFs can act as barriers against nucleosome sliding towards NDRs. However, the relative ability of the different GRFs to hinder sliding activity is currently unknown.

RESULTS

Considering this, we performed a comparative analysis for the main GRFs, with focus in their ability to modulate nucleosome sliding mediated by ISW1a. Among the GRFs tested in nucleosome remodeling assays, Rap1 was the only factor displaying the ability to hinder the activity of ISW1a. This effect requires location of the Rap1 cognate sequence on linker that becomes entry DNA in the nucleosome remodeling process. In addition, Rap1 was able to hinder nucleosome assembly in octamer transfer assays. Concurrently, Rap1 displayed the highest affinity for and longest dwell time from its target sequence, compared to the other GRFs tested. Consistently, through bioinformatics analyses of publicly available genome-wide data, we found that nucleosome occupancy and histone deposition in vivo are inversely correlated with the affinity of Rap1 for its target sequences in the genome.

CONCLUSIONS

Our findings point to DNA binding affinity, residence time and location at particular translational positions relative to the nucleosome core as the key features of GRFs underlying their roles played in nucleosome sliding and assembly.

Genome organization across scales: mechanistic insights from in vitro reconstitution studies

Eukaryotic genomes are compacted and organized into distinct three-dimensional (3D) structures, which range from small-scale nucleosome arrays to large-scale chromatin domains. These chromatin structures play an important role in the regulation of transcription and other nuclear processes. The molecular mechanisms that drive the formation of chromatin structures across scales and the relationship between chromatin structure and function remain incompletely understood. Because the processes involved are complex and interconnected, it is often challenging to dissect the underlying principles in the nuclear environment. Therefore, in vitro reconstitution systems provide a valuable approach to gain insight into the molecular mechanisms by which chromatin structures are formed and to determine the cause-consequence relationships between the processes involved. In this review, we give an overview of in vitro approaches that have been used to study chromatin structures across scales and how they have increased our understanding of the formation and function of these structures. We start by discussing in vitro studies that have given insight into the mechanisms of nucleosome positioning. Next, we discuss recent efforts to reconstitute larger-scale chromatin domains and loops and the resulting insights into the principles of genome organization. We conclude with an outlook on potential future applications of chromatin reconstitution systems and how they may contribute to answering open questions concerning chromatin architecture.

Yeast heterochromatin stably silences only weak regulatory elements by altering burst duration

Transcriptional silencing in Saccharomyces cerevisiae involves the generation of a chromatin state that stably represses transcription. Using multiple reporter assays, a diverse set of upstream activating sequence enhancers and core promoters were investigated for their susceptibility to silencing. We show that heterochromatin stably silences only weak and stress-induced regulatory elements but is unable to stably repress housekeeping gene regulatory elements, and the partial repression of these elements did not result in bistable expression states. Permutation analysis of enhancers and promoters indicates that both elements are targets of repression. Chromatin remodelers help specific regulatory elements to resist repression, most probably by altering nucleosome mobility and changing transcription burst duration. The strong enhancers/promoters can be repressed if silencer-bound Sir1 is increased. Together, our data suggest that the heterochromatic locus has been optimized to stably silence the weak mating-type gene regulatory elements but not strong housekeeping gene regulatory sequences.

HIV Expression in Infected T Cell Clones

The principal barrier to an HIV-1 cure is the persistence of infected cells harboring replication-competent proviruses despite antiretroviral therapy (ART). HIV-1 transcriptional suppression, referred to as viral latency, is foremost among persistence determinants, as it allows infected cells to evade the cytopathic effects of virion production and killing by cytotoxic T lymphocytes (CTL) and other immune factors. HIV-1 persistence is also governed by cellular proliferation, an innate and essential capacity of CD4+ T cells that both sustains cell populations over time and enables a robust directed response to immunological threats. However, when HIV-1 infects CD4+ T cells, this capacity for proliferation can enable surreptitious HIV-1 propagation without the deleterious effects of viral gene expression in latently infected cells. Over time on ART, the HIV-1 reservoir is shaped by both persistence determinants, with selective forces most often favoring clonally expanded infected cell populations harboring transcriptionally quiescent proviruses. Moreover, if HIV latency is incomplete or sporadically reversed in clonal infected cell populations that are replenished faster than they are depleted, such populations could both persist indefinitely and contribute to low-level persistent viremia during ART and viremic rebound if treatment is withdrawn. In this review, select genetic, epigenetic, cellular, and immunological determinants of viral transcriptional suppression and clonal expansion of HIV-1 reservoir T cells, interdependencies among these determinants, and implications for HIV-1 persistence will be presented and discussed.

A single fiber view of the nucleosome organization in eukaryotic chromatin

Eukaryotic cells are thought to arrange nucleosomes into extended arrays with evenly spaced nucleosomes phased at genomic landmarks. Here we tested to what extent this stereotypic organization describes the nucleosome organization in Saccharomyces cerevisiae using Fiber-Seq, a long-read sequencing technique that maps entire nucleosome arrays on individual chromatin fibers in a high throughput manner. With each fiber coming from a different cell, Fiber-Seq uncovers cell-to-cell heterogeneity. The long reads reveal the nucleosome architecture even over repetitive DNA such as the ribosomal DNA repeats. The absolute nucleosome occupancy, a parameter that is difficult to obtain with conventional sequencing approaches, is a direct readout of Fiber-Seq. We document substantial deviations from the stereotypical nucleosome organization with unexpectedly long linker DNAs between nucleosomes, gene bodies missing entire nucleosomes, cell-to-cell heterogeneity in nucleosome occupancy, heterogeneous phasing of arrays and irregular nucleosome spacing. Nucleosome array structures are indistinguishable throughout the gene body and with respect to the direction of transcription arguing against transcription promoting array formation. Acute nucleosome depletion destroyed most of the array organization indicating that nucleosome remodelers cannot efficiently pack nucleosomes under those conditions. Given that nucleosomes are cis-regulatory elements, the cell-to-cell heterogeneity uncovered by Fiber-Seq provides much needed information to understand chromatin structure and function.

Epigenetic remodelling under hypoxia

Hypoxia is intrinsic to tumours and contributes to malignancy and metastasis while hindering the efficiency of existing treatments. Epigenetic mechanisms play a crucial role in the regulation of hypoxic cancer cell programs, both in the initial phases of sensing the decrease in oxygen levels and during adaptation to chronic lack of oxygen. During the latter, the epigenetic regulation of tumour biology intersects with hypoxia-sensitive transcription factors in a complex network of gene regulation that also involves metabolic reprogramming. Here, we review the current literature on the epigenetic control of gene programs in hypoxic cancer cells. We highlight common themes and features of such epigenetic remodelling and discuss their relevance for the development of therapeutic strategies.

Examining chromatin heterogeneity through PacBio long-read sequencing of M.EcoGII methylated genomes: an m6A detection efficiency and calling bias correcting pipeline

Recent studies have combined DNA methyltransferase footprinting of genomic DNA in nuclei with long-read sequencing, resulting in detailed chromatin maps for multi-kilobase stretches of genomic DNA from one cell. Theoretically, nucleosome footprints and nucleosome-depleted regions can be identified using M.EcoGII, which methylates adenines in any sequence context, providing a high-resolution map of accessible regions in each DNA molecule. Here we report PacBio long-read sequence data for budding yeast nuclei treated with M.EcoGII and a bioinformatic pipeline which corrects for three key challenges undermining this promising method. First, detection of m6A in individual DNA molecules by the PacBio software is inefficient, resulting in false footprints predicted by random gaps of seemingly unmethylated adenines. Second, there is a strong bias against m6A base calling as AT content increases. Third, occasional methylation occurs within nucleosomes, breaking up their footprints. After correcting for these issues, our pipeline calculates a correlation coefficient-based score indicating the extent of chromatin heterogeneity within the cell population for every gene. Although the population average is consistent with that derived using other techniques, we observe a wide range of heterogeneity in nucleosome positions at the single-molecule level, probably reflecting cellular chromatin dynamics.

Improving genome-wide mapping of nucleosomes in Trypanosome cruzi

In Trypanosoma cruzi DNA is packaged into chromatin by octamers of histone proteins that form nucleosomes. Transcription of protein coding genes in trypanosomes is constitutive producing polycistronic units and gene expression is primarily regulated post-transcriptionally. However, chromatin organization influences DNA dependent processes. Hence, determining nucleosome position is of uppermost importance to understand the peculiarities found in trypanosomes. To map nucleosomes genome-wide in several organisms, digestion of chromatin with micrococcal nuclease followed by deep sequencing has been applied. Nonetheless, the special requirements for cell manipulation and the uniqueness of the chromatin organization in trypanosomes entails a customized analytical approach. In this work, we adjusted this broadly used method to the hybrid reference strain, CL Brener. Particularly, we implemented an exhaustive and thorough computational workflow to overcome the difficulties imposed by this complex genome. We tested the performance of two aligners, Bowtie2 and HISAT2, and discuss their advantages and caveats. Specifically, we highlight the relevance of using the whole genome as a reference instead of the commonly used Esmeraldo-like haplotype to avoid spurious alignments. Additionally, we show that using the whole genome refines the average nucleosome representation, but also the quality of mapping for every region represented. Moreover, we show that the average nucleosome organization around trans-splicing acceptor site described before, is not just an average since the same chromatin pattern is detected for most of the represented regions. In addition, we extended the study to a non-hybrid strain applying the experimental and analytical approach to Sylvio-X10 strain. Furthermore, we provide a source code for the construction of 2D plots and heatmaps which are easy to adapt to any T. cruzi strain.

Yeast Heterochromatin Only Stably Silences Weak Regulatory Elements by Altering Burst Duration

The interplay between nucleosomes and transcription factors leads to programs of gene expression. Transcriptional silencing involves the generation of a chromatin state that represses transcription and is faithfully propagated through DNA replication and cell division. Using multiple reporter assays, including directly visualizing transcription in single cells, we investigated a diverse set of UAS enhancers and core promoters for their susceptibility to heterochromatic gene silencing. These results show that heterochromatin only stably silences weak and stress induced regulatory elements but is unable to stably repress housekeeping gene regulatory elements and the partial repression did not result in bistable expression states. Permutation analysis of different UAS enhancers and core promoters indicate that both elements function together to determine the susceptibility of regulatory sequences to repression. Specific histone modifiers and chromatin remodellers function in an enhancer specific manner to aid these elements to resist repression suggesting that Sir proteins likely function in part by reducing nucleosome mobility. We also show that the strong housekeeping regulatory elements can be repressed if silencer bound Sir1 is increased, suggesting that Sir1 is a limiting component in silencing. Together, our data suggest that the heterochromatic locus has been optimized to stably silence the weak mating type gene regulatory elements but not strong housekeeping gene regulatory sequences which could help explain why these genes are often found at the boundaries of silenced domains.

Measuring occupancies of the nucleosome and nucleosome-interacting factors in vivo in Saccharomyces cerevisiae genome-wide

Nucleosomes are the repeating units of chromatin. The presence of nucleosomes poses a major impediment to all DNA-dependent processes. As a result, access to DNA in chromatin is dynamically regulated by many factors, including ATP-dependent chromatin remodeling complexes. Digestion of chromatin by micrococcal nuclease (MNase) followed by chromatin immunoprecipitation (ChIP) and sequencing can be leveraged to determine nucleosome occupancy, positioning, and the ability of chromatin interacting factors to alter chromatin accessibility. Here we describe the procedure for performing MNase and MNase ChIP-seq in detail.

Role of Nucleosome Sliding in the Protein Target Search for Covered DNA Sites

Associations of transcription factors (TFs) with specific sites on DNA initiate major cellular processes. But DNA in eukaryotic cells is covered by nucleosomes which prevent TFs from binding. However, nucleosome structures on DNA are not static and exhibit breathing and sliding. We develop a theoretical framework to investigate the effect of nucleosome sliding on a protein target search. By analysis of a discrete-state stochastic model of nucleosome sliding, search dynamics are explicitly evaluated. It is found that for long sliding lengths the target search dynamics are faster for normal TFs that cannot enter the nucleosomal DNA. But for more realistic short sliding lengths, the so-called pioneer TFs, which can invade nucleosomal DNA, locate specific sites faster. It is also suggested that nucleosome breathing, which is a faster process, has a stronger effect on protein search dynamics than that of nucleosome sliding. Theoretical arguments to explain these observations are presented.

Transcription and FACT facilitate the restoration of replication-coupled chromatin assembly defects

Genome duplication occurs through the coordinated action of DNA replication and nucleosome assembly at replication forks. Defective nucleosome assembly causes DNA lesions by fork breakage that need to be repaired. In addition, it causes a loss of chromatin integrity. These chromatin alterations can be restored, even though the mechanisms are unknown. Here, we show that the process of chromatin restoration can deal with highly severe chromatin defects induced by the absence of the chaperones CAF1 and Rtt106 or a strong reduction in the pool of available histones, and that this process can be followed by analyzing the topoisomer distribution of the 2µ plasmid. Using this assay, we demonstrate that chromatin restoration is slow and independent of checkpoint activation, whereas it requires the action of transcription and the FACT complex. Therefore, cells are able to "repair" not only DNA lesions but also chromatin alterations associated with defective nucleosome assembly.

Dynamic nucleosome remodeling mediated by YY1 underlies early mouse development

Nucleosome positioning can alter the accessibility of DNA-binding proteins to their cognate DNA elements, and thus its precise control is essential for cell identity and function. Mammalian preimplantation embryos undergo temporal changes in gene expression and cell potency, suggesting the involvement of dynamic epigenetic control during this developmental phase. However, the dynamics of nucleosome organization during early development are poorly understood. In this study, using a low-input MNase-seq method, we show that nucleosome positioning is globally obscure in zygotes but becomes well defined during subsequent development. Down-regulation of the chromatin assembly in embryonic stem cells can partially reverse nucleosome organization into a zygote-like pattern, suggesting a possible link between the chromatin assembly pathway and fuzzy nucleosomes in zygotes. We also reveal that YY1, a zinc finger-containing transcription factor expressed upon zygotic genome activation, regulates the de novo formation of well-positioned nucleosome arrays at the regulatory elements through identifying YY1-binding sites in eight-cell embryos. The YY1-binding regions acquire H3K27ac enrichment around the eight-cell and morula stages, and YY1 depletion impairs the morula-to-blastocyst transition. Thus, our study delineates the remodeling of nucleosome organization and its underlying mechanism during early mouse development.

Context dependent prediction in DNA sequence using neural networks

One way to better understand the structure in DNA is by learning to predict the sequence. Here, we trained a model to predict the missing base at any given position, given its left and right flanking contexts. Our best-performing model was a neural network that obtained an accuracy close to 54% on the human genome, which is 2% points better than modelling the data using a Markov model. In likelihood-ratio tests, the neural network performed significantly better than any of the alternative models by a large margin. We report on where the accuracy was obtained, first observing that the performance appeared to be uniform over the chromosomes. The models performed best in repetitive sequences, as expected, although their performance far from random in the more difficult coding sections, the proportions being ~70:40%. We further explored the sources of the accuracy, Fourier transforming the predictions revealed weak but clear periodic signals. In the human genome the characteristic periods hinted at connections to nucleosome positioning. We found similar periodic signals in GC/AT content in the human genome, which to the best of our knowledge have not been reported before. On other large genomes similarly high accuracy was found, while lower predictive accuracy was observed on smaller genomes. Only in the mouse genome did we see periodic signals in the same range as in the human genome, though weaker and of a different type. This indicates that the sources of these signals are other or more than nucleosome arrangement. Interestingly, applying a model trained on the mouse genome to the human genome resulted in a performance far below that of the human model, except in the difficult coding regions. Despite the clear outcomes of the likelihood-ratio tests, there is currently a limited superiority of the neural network methods over the Markov model. We expect, however, that there is great potential for better modelling DNA using different neural network architectures.

Production of a viral surface protein in Nannochloropsis oceanica for fish vaccination against infectious pancreatic necrosis virus

Abstract

Nannochloropsis oceanica is a unicellular oleaginous microalga of emerging biotechnological interest with a sequenced, annotated genome, available transcriptomic and proteomic data, and well-established basic molecular tools for genetic engineering. To establish N. oceanica as a eukaryotic host for recombinant protein synthesis and develop molecular technology for vaccine production, we chose the viral surface protein 2 (VP2) of a pathogenic fish virus that causes infectious pancreatic necrosis as a model vaccine. Upon stable nuclear transformation of N. oceanica strain CCMP1779 with the codon-optimized VP2 gene, a Venus reporter fusion served to evaluate the strength of different endogenous promoters in transformant populations by qPCR and flow cytometry. The highest VP2 yields were achieved for the elongation factor promoter, with enhancer effects by its N-terminal leader sequence. Individual transformants differed in their production capability of reporter-free VP2 by orders of magnitude. When subjecting the best candidates to kinetic analyses of growth and VP2 production in photobioreactors, recombinant protein integrity was maintained until the early stationary growth phase, and a high yield of 4.4% VP2 of total soluble protein was achieved. The maximum yield correlated with multiple integrations of the expression vector into the nuclear genome. The results demonstrate that N. oceanica was successfully engineered to constitute a robust platform for high-level production of a model subunit vaccine. The molecular methodology established here can likely be adapted in a straightforward manner to the production of further vaccines in the same host, allowing their distribution to fish, vertebrates, or humans via a microalgae-containing diet.

Key points

• We engineered N. oceanica for recombinant protein production.

• The antigenic surface protein 2 of IPN virus could indeed be expressed in the host.

• A high yield of 4.4% VP2 of total soluble protein was achieved in N. oceanica.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00253-022-12106-7.

A study of strong nucleosomes in the human genome

Micrococcal nuclease (MNase) is widely used to map nucleosomes. However, nucleosomes are highly dynamic and susceptible to experimental conditions, resulting in extreme variability across nucleosome maps, which complicates the generation of accurate nucleosome organization data. We mapped nucleosomes from different individuals using improved MNase-seq. The improvements included setting different digestion levels (low, medium, high) and naked DNA correction to remove the noise caused by experimental manipulation and comparing maps to obtain the accurate position and occupancy of strong nucleosomes (SNs) in the whole genome. In addition, the characteristics of SNs were further excavated. SNs were enriched in Alu elements and near the centromere of Chr12. SNs contain some specific sequences, and the GC content of SNs is different from that of dynamic nucleosomes. The findings suggest that nucleosome location in the genome and the DNA sequence may affect nucleosome stability.

•

Naked DNA correction improved the accuracy of nucleosome map in partial digestion

•

Level of MNase digestion has effects on nucleosome organization

•

A type of strong nucleosomes (SNs) exist across different nucleosome maps

•

Nucleosome stability may be related to its location and the DNA sequence

Aberrant nucleosome organization in mouse SCNT embryos revealed by ULI-MNase-seq

Somatic cell nuclear transfer (SCNT) can reprogram terminally differentiated somatic cells into totipotent embryos, but with multiple defects. The nucleosome positioning, as an important epigenetic regulator for gene expression, is largely unexplored during SCNT embryonic development. Here, we mapped genome-wide nucleosome profiles in mouse SCNT embryos using ultra-low-input MNase-seq (ULI-MNase-seq). We found that the nucleosome-depleted regions (NDRs) around promoters underwent dramatic reestablishment, which is consistent with the cell cycle. Dynamics of nucleosome position in SCNT embryos were delayed compared to fertilized embryos. Subsequently, we found that the aberrant gene expression levels in inner cell mass (ICM) were positively correlated with promoter NDRs in donor cells, which indicated that the memory of nucleosome occupancy in donor cells was a potential barrier for SCNT-mediated reprogramming. We further confirmed that the histone acetylation level of donor cells was associated with the memory of promoter NDRs. Our study provides insight into nucleosome reconfiguration during SCNT preimplantation embryonic development.

Dynamic nucleosome organization after fertilization reveals regulatory factors for mouse zygotic genome activation

Chromatin remodeling is essential for epigenome reprogramming after fertilization. However, the underlying mechanisms of chromatin remodeling remain to be explored. Here, we investigated the dynamic changes in nucleosome occupancy and positioning in pronucleus-stage zygotes using ultra low-input MNase-seq. We observed distinct features of inheritance and reconstruction of nucleosome positioning in both paternal and maternal genomes. Genome-wide de novo nucleosome occupancy in the paternal genome was observed as early as 1 h after the injection of sperm into ooplasm. The nucleosome positioning pattern was continually rebuilt to form nucleosome-depleted regions (NDRs) at promoters and transcription factor (TF) binding sites with differential dynamics in paternal and maternal genomes. NDRs formed more quickly on the promoters of genes involved in zygotic genome activation (ZGA), and this formation is closely linked to histone acetylation, but not transcription elongation or DNA replication. Importantly, we found that NDR establishment on the binding motifs of specific TFs might be associated with their potential pioneer functions in ZGA. Further investigations suggested that the predicted factors MLX and RFX1 played important roles in regulating minor and major ZGA, respectively. Our data not only elucidate the nucleosome positioning dynamics in both male and female pronuclei following fertilization, but also provide an efficient method for identifying key transcription regulators during development.

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.

Mapping of Nonhomologous End Joining-Mediated Integration Facilitates Genome-Scale Trackable Mutagenesis in Yarrowia lipolytica

Genome-scale mutagenesis, phenotypic screening, and tracking the causal mutations is a powerful approach for genetic analysis. However, classic mutagenesis approaches require extensive effort to identify causal mutations. It is desirable to demonstrate a powerful approach for rapid trackable mutagenesis. Here, we mapped the distribution of nonhomologous end joining (NHEJ)-mediated integration for the first time and demonstrated that it can be used for constructing the genome-scale trackable mutagenesis library in Yarrowia lipolytica. The sequencing of 9.15 × 10⁵ insertions showed that NHEJ-mediated integration inserted DNA randomly across the chromosomes, and the transcriptional regulatory regions exhibited integration preference. The insertions were located in both nucleosome-occupancy regions and nucleosome-free regions. Using NHEJ-mediated integration to construct the genome-scale mutagenesis library, the new targets that improved β-carotene biosynthesis and acetic acid tolerance were identified rapidly. This mutagenesis approach is readily applicable to other organisms with strong NHEJ preference and will contribute to cell factory construction.

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.

Formation of artificial chromosomes in Caenorhabditis elegans and analyses of their segregation in mitosis, DNA sequence composition and holocentromere organization

To investigate how exogenous DNA concatemerizes to form episomal artificial chromosomes (ACs), acquire equal segregation ability and maintain stable holocentromeres, we injected DNA sequences with different features, including sequences that are repetitive or complex, and sequences with different AT-contents, into the gonad of Caenorhabditis elegans to form ACs in embryos, and monitored AC mitotic segregation. We demonstrated that AT-poor sequences (26% AT-content) delayed the acquisition of segregation competency of newly formed ACs. We also co-injected fragmented Saccharomyces cerevisiae genomic DNA, differentially expressed fluorescent markers and ubiquitously expressed selectable marker to construct a less repetitive, more complex AC. We sequenced the whole genome of a strain which propagates this AC through multiple generations, and de novo assembled the AC sequences. We discovered CENP-AHCP-3 domains/peaks are distributed along the AC, as in endogenous chromosomes, suggesting a holocentric architecture. We found that CENP-AHCP-3 binds to the unexpressed marker genes and many fragmented yeast sequences, but is excluded in the yeast extremely high-AT-content centromeric and mitochondrial DNA (> 83% AT-content) on the AC. We identified A-rich motifs in CENP-AHCP-3 domains/peaks on the AC and on endogenous chromosomes, which have some similarity with each other and similarity to some non-germline transcription factor binding sites.

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.

Epigenetic Mechanisms of HIV-1 Persistence

Eradicating HIV-1 in infected individuals will not be possible without addressing the persistence of the virus in its multiple reservoirs. In this context, the molecular characterization of HIV-1 persistence is key for the development of rationalized therapeutic interventions. HIV-1 gene expression relies on the redundant and cooperative recruitment of cellular epigenetic machineries to cis-regulatory proviral regions. Furthermore, the complex repertoire of HIV-1 repression mechanisms varies depending on the nature of the viral reservoir, although, so far, few studies have addressed the specific regulatory mechanisms of HIV-1 persistence in other reservoirs than the well-studied latently infected CD4⁺ T cells. Here, we present an exhaustive and updated picture of the heterochromatinization of the HIV-1 promoter in its different reservoirs. We highlight the complexity, heterogeneity and dynamics of the epigenetic mechanisms of HIV-1 persistence, while discussing the importance of further understanding HIV-1 gene regulation for the rational design of novel HIV-1 cure strategies.

The yeast ISW1b ATP-dependent chromatin remodeler is critical for nucleosome spacing and dinucleosome resolution

Isw1 and Chd1 are ATP-dependent nucleosome-spacing enzymes required to establish regular arrays of phased nucleosomes near transcription start sites of yeast genes. Cells lacking both Isw1 and Chd1 have extremely disrupted chromatin, with weak phasing, irregular spacing and a propensity to form close-packed dinucleosomes. The Isw1 ATPase subunit occurs in two different remodeling complexes: ISW1a (composed of Isw1 and Ioc3) and ISW1b (composed of Isw1, Ioc2 and Ioc4). The Ioc4 subunit of ISW1b binds preferentially to the H3-K36me3 mark. Here we show that ISW1b is primarily responsible for setting nucleosome spacing and resolving close-packed dinucleosomes, whereas ISW1a plays only a minor role. ISW1b and Chd1 make additive contributions to dinucleosome resolution, such that neither enzyme is capable of resolving all dinucleosomes on its own. Loss of the Set2 H3-K36 methyltransferase partly phenocopies loss of Ioc4, resulting in increased dinucleosome levels with only a weak effect on nucleosome spacing, suggesting that Set2-mediated H3-K36 trimethylation contributes to ISW1b-mediated dinucleosome separation. The H4 tail domain is required for normal nucleosome spacing but not for dinucleosome resolution. We conclude that the nucleosome spacing and dinucleosome resolving activities of ISW1b and Chd1 are critical for normal global chromatin organisation.

Long non-coding RNAs: the tentacles of chromatin remodeler complexes

Chromatin remodeler complexes regulate gene transcription, DNA replication and DNA repair by changing both nucleosome position and post-translational modifications. The chromatin remodeler complexes are categorized into four families: the SWI/SNF, INO80/SWR1, ISWI and CHD family. In this review, we describe the subunits of these chromatin remodeler complexes, in particular, the recently identified members of the ISWI family and novelties of the CHD family. Long non-coding (lnc) RNAs regulate gene expression through different epigenetic mechanisms, including interaction with chromatin remodelers. For example, interaction of lncBRM with BRM inhibits the SWI/SNF complex associated with a differentiated phenotype and favors assembly of a stem cell-related SWI/SNF complex. Today, over 50 lncRNAs have been shown to affect chromatin remodeler complexes and we here discuss the mechanisms involved.

On the role of transcription in positioning nucleosomes

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

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

CpG Islands Shape the Epigenome Landscape

Epigenetic modifications and nucleosome positioning play an important role in modulating gene expression. However, how the patterns of epigenetic modifications and nucleosome positioning are established around promoters is not well understood. Here, we have addressed these questions in a series of genome-wide experiments coupled to a novel bioinformatic analysis approach. Our data reveal a clear correlation between CpG density, promoter activity and accumulation of active or repressive histone marks. CGI boundaries define the chromatin promoter regions that will be epigenetically modified. CpG-rich promoters are targeted by histone modifications and histone variants, while CpG-poor promoters are regulated by DNA methylation. CGIs boundaries, but not transcriptional activity, are essential determinants of H2A.Z positioning in vicinity of the promoters, suggesting that the presence of H2A.Z is not related to transcriptional control. Accordingly, H2A.Z depletion has no impact on gene expression of arrested mouse embryonic fibroblasts. Therefore, the underlying DNA sequence, the promoter CpG density and, to a lesser extent, transcriptional activity, are key factors implicated in promoter chromatin architecture.

CORENup: a combination of convolutional and recurrent deep neural networks for nucleosome positioning identification

BACKGROUND

Nucleosomes wrap the DNA into the nucleus of the Eukaryote cell and regulate its transcription phase. Several studies indicate that nucleosomes are determined by the combined effects of several factors, including DNA sequence organization. Interestingly, the identification of nucleosomes on a genomic scale has been successfully performed by computational methods using DNA sequence as input data.

RESULTS

In this work, we propose CORENup, a deep learning model for nucleosome identification. CORENup processes a DNA sequence as input using one-hot representation and combines in a parallel fashion a fully convolutional neural network and a recurrent layer. These two parallel levels are devoted to catching both non periodic and periodic DNA string features. A dense layer is devoted to their combination to give a final classification.

CONCLUSIONS

Results computed on public data sets of different organisms show that CORENup is a state of the art methodology for nucleosome positioning identification based on a Deep Neural Network architecture. The comparisons have been carried out using two groups of datasets, currently adopted by the best performing methods, and CORENup has shown top performance both in terms of classification metrics and elapsed computation time.

Nucleosome movement analysis based on second-order information entropy and density functional theory

Dynamics of +1 and -1 nucleosomes near TSS of yeast chromosome 2 were analyzed by using second-order information entropy and density functional theory method. Second-order information entropy can measure the interaction intensity between nucleosome sequences and nucleosome histones based on the intensity of base association. In addition, density functional theory method can be used to obtain the global interaction intensity between nucleosome sequences and nucleosome histones based on energy state size and active or non-active state of binucleoside pairs. Our results showed asymmetry of interaction intensity on both sides of the nucleosome central site, and that +1 nucleosomes tend to move toward the 5'-end and -1 nucleosomes tend to move toward the 3'-end. Under the dynamic balance of nucleosome movement, in roder to shut down gene transcription, +1 and -1 nucleosomes will cover TSS. If the dynamic balance is destroyed, +1 and -1 nucleosomes stay away from each other to expose TSS to restart gene transcription. The movement trend of +1 and -1 nucleosomes coincides with the biological mechanism of gene transcription and non-transcription, and the nucleosome sequences contain the dynamic information of nucleosome movement, which provides effective technical support for the study of gene transcription regulation mechanism.

Interplay among ATP-Dependent Chromatin Remodelers Determines Chromatin Organisation in Yeast

Cellular DNA is packaged into chromatin, which is composed of regularly-spaced nucleosomes with occasional gaps corresponding to active regulatory elements, such as promoters and enhancers, called nucleosome-depleted regions (NDRs). This chromatin organisation is primarily determined by the activities of a set of ATP-dependent remodeling enzymes that are capable of moving nucleosomes along DNA, or of evicting nucleosomes altogether. In yeast, the nucleosome-spacing enzymes are ISW1 (Imitation SWitch protein 1), Chromodomain-Helicase-DNA-binding (CHD)1, ISW2 (Imitation SWitch protein 2) and INOsitol-requiring 80 (INO80); the nucleosome eviction enzymes are the SWItching/Sucrose Non-Fermenting (SWI/SNF) family, the Remodeling the Structure of Chromatin (RSC) complexes and INO80. We discuss the contributions of each set of enzymes to chromatin organisation. ISW1 and CHD1 are the major spacing enzymes; loss of both enzymes results in major chromatin disruption, partly due to the appearance of close-packed di-nucleosomes. ISW1 and CHD1 compete to set nucleosome spacing on most genes. ISW1 is dominant, setting wild type spacing, whereas CHD1 sets short spacing and may dominate on highly-transcribed genes. We propose that the competing remodelers regulate spacing, which in turn controls the binding of linker histone (H1) and therefore the degree of chromatin folding. Thus, genes with long spacing bind more H1, resulting in increased chromatin compaction. RSC, SWI/SNF and INO80 are involved in NDR formation, either directly by nucleosome eviction or repositioning, or indirectly by affecting the size of the complex that resides in the NDR. The nature of this complex is controversial: some suggest that it is a RSC-bound “fragile nucleosome”, whereas we propose that it is a non-histone transcription complex. In either case, this complex appears to serve as a barrier to nucleosome formation, resulting in the formation of phased nucleosomal arrays on both sides.

Chromatin replication and epigenetic cell memory

Propagation of the chromatin landscape across cell divisions is central to epigenetic cell memory. Mechanistic analysis of the interplay between DNA replication, the cell cycle, and the epigenome has provided insights into replication-coupled chromatin assembly and post-replicative chromatin maintenance. These breakthroughs are critical for defining how proliferation impacts the epigenome during cell identity changes in development and disease. Here we review these findings in the broader context of epigenetic inheritance across mitotic cell division.

Epigenome Regulation by Dynamic Nucleosome Unwrapping

Gene regulation in eukaryotes requires the controlled access of sequence-specific transcription factors (TFs) to their sites in a chromatin landscape dominated by nucleosomes. Nucleosomes are refractory to TF binding, and often must be removed from regulatory regions. Recent genomic studies together with in vitro measurements suggest that the nucleosome barrier to TF binding is modulated by dynamic nucleosome unwrapping governed by ATP-dependent chromatin remodelers. Genome-wide occupancy and the regulation of subnucleosomal intermediates have gained recent attention with the application of high-resolution approaches for precision mapping of protein-DNA interactions. We summarize here recent findings on nucleosome substructures and TF binding dynamics, and highlight how unwrapped nucleosomal intermediates provide a novel signature of active chromatin.

Creating 2D Occupancy Plots Using plot2DO

Chromatin organization and epigenetic marks play a critical role in stem cell pluripotency and differentiation. Chromatin digestion by micrococcal nuclease (MNase) followed by high-throughput sequencing (MNase-seq) is the most widely used genome-wide method for studying nucleosome organization, that is, the first level of DNA packaging into chromatin. Combined with chromatin immunoprecipitation (ChIP), MNase-ChIP-seq represents a high-resolution method for investigating both chromatin organization and the distribution of epigenetic marks and histone variants. The plot2DO package presented here is a flexible tool for evaluating the quality of MNase-seq and MNase-ChIP-seq data, and for visualizing the distribution of nucleosomes near the functional regions of the genome. The plot2DO package is open-source software, and it is freely available from https://github.com/rchereji/plot2DO under the MIT license.

Differential nucleosome spacing in neurons and glia

Eukaryotic chromosomes are composed of chromatin, in which regularly spaced nucleosomes containing ∼147 bp of DNA are separated by linker DNA. Most eukaryotic cells have a characteristic average nucleosome spacing of ∼190 bp, corresponding to a ∼45 bp linker. However, cortical neurons have a shorter average spacing of ∼165 bp. The significance of this atypical global chromatin organization is unclear. We have compared the chromatin structures of purified mouse dorsal root ganglia (DRG) neurons, cortical oligodendrocyte precursor cells (OPCs) and cortical astrocytes. DRG neurons have short average spacing (∼165 bp), whereas OPCs (∼182 bp) and astrocytes (∼183 bp) have longer spacing. We measured nucleosome positions by MNase-seq and gene expression by RNA-seq. Most genes in all three cell types have a promoter chromatin organization typical of active genes: a nucleosome-depleted region at the promoter flanked by regularly spaced nucleosomes phased relative to the transcription start site. In DRG neurons, the spacing of phased nucleosomes downstream of promoters (∼182 bp) is longer than expected from the genomic average for DRG neurons, whereas phased nucleosome spacing in OPCs and astrocytes is similar to the global average for these cells (∼183 bp). Thus, the atypical nucleosome spacing of neuronal chromatin does not extend to promoter-proximal regions.

Good News for Nuclear Transgene Expression in Chlamydomonas

Chlamydomonas reinhardtii is a well-established model system for basic research questions ranging from photosynthesis and organelle biogenesis, to the biology of cilia and basal bodies, to channelrhodopsins and photoreceptors. More recently, Chlamydomonas has also been recognized as a suitable host for the production of high-value chemicals and high-value recombinant proteins. However, basic and applied research have suffered from the inefficient expression of nuclear transgenes. The combined efforts of the Chlamydomonas community over the past decades have provided insights into the mechanisms underlying this phenomenon and have resulted in mutant strains defective in some silencing mechanisms. Moreover, many insights have been gained into the parameters that affect nuclear transgene expression, like promoters, introns, codon usage, or terminators. Here I critically review these insights and try to integrate them into design suggestions for the construction of nuclear transgenes that are to be expressed at high levels.

Training-free measures based on algorithmic probability identify high nucleosome occupancy in DNA sequences

We introduce and study a set of training-free methods of an information-theoretic and algorithmic complexity nature that we apply to DNA sequences to identify their potential to identify nucleosomal binding sites. We test the measures on well-studied genomic sequences of different sizes drawn from different sources. The measures reveal the known in vivo versus in vitro predictive discrepancies and uncover their potential to pinpoint high and low nucleosome occupancy. We explore different possible signals within and beyond the nucleosome length and find that the complexity indices are informative of nucleosome occupancy. We found that, while it is clear that the gold standard Kaplan model is driven by GC content (by design) and by k-mer training; for high occupancy, entropy and complexity-based scores are also informative and can complement the Kaplan model.

Theory of Active Chromatin Remodeling

Nucleosome positioning controls the accessible regions of chromatin and plays essential roles in DNA-templated processes. ATP driven remodeling enzymes are known to be crucial for its establishment in vivo, but their nonequilibrium nature has hindered the development of a unified theoretical framework for nucleosome positioning. Using a perturbation theory, we show that the effect of these enzymes can be well approximated by effective equilibrium models with rescaled temperatures and interactions. Numerical simulations support the accuracy of the theory in predicting both kinetic and steady-state quantities, including the effective temperature and the radial distribution function, in biologically relevant regimes. The energy landscape view emerging from our study provides an intuitive understanding for the impact of remodeling enzymes in either reinforcing or overwriting intrinsic signals for nucleosome positioning, and may help improve the accuracy of computational models for its prediction in silico.

Quantitative MNase-seq accurately maps nucleosome occupancy levels

Micrococcal nuclease (MNase) is widely used to map nucleosomes. However, its aggressive endo-/exo-nuclease activities make MNase-seq unreliable for determining nucleosome occupancies, because cleavages within linker regions produce oligo- and mono-nucleosomes, whereas cleavages within nucleosomes destroy them. Here, we introduce a theoretical framework for predicting nucleosome occupancies and an experimental protocol with appropriate spike-in normalization that confirms our theory and provides accurate occupancy levels over an MNase digestion time course. As with human cells, we observe no overall differences in nucleosome occupancies between Drosophila euchromatin and heterochromatin, which implies that heterochromatic compaction does not reduce MNase accessibility of linker DNA.

Accessibility of promoter DNA is not the primary determinant of chromatin-mediated gene regulation

DNA accessibility is thought to be of major importance in regulating gene expression. We test this hypothesis using a restriction enzyme as a probe of chromatin structure and as a proxy for transcription factors. We measured the digestion rate and the fraction of accessible DNA at almost all genomic AluI sites in budding yeast and mouse liver nuclei. Hepatocyte DNA is more accessible than yeast DNA, consistent with longer linkers between nucleosomes, suggesting that nucleosome spacing is a major determinant of accessibility. DNA accessibility varies from cell to cell, such that essentially no sites are accessible or inaccessible in every cell. AluI sites in inactive mouse promoters are accessible in some cells, implying that transcription factors could bind without activating the gene. Euchromatin and heterochromatin have very similar accessibilities, suggesting that transcription factors can penetrate heterochromatin. Thus, DNA accessibility is not likely to be the primary determinant of gene regulation.

Old cogs, new tricks: the evolution of gene expression in a chromatin context

Sophisticated gene-regulatory mechanisms probably evolved in prokaryotes billions of years before the emergence of modern eukaryotes, which inherited the same basic enzymatic machineries. However, the epigenomic landscapes of eukaryotes are dominated by nucleosomes, which have acquired roles in genome packaging, mitotic condensation and silencing parasitic genomic elements. Although the molecular mechanisms by which nucleosomes are displaced and modified have been described, just how transcription factors, histone variants and modifications and chromatin regulators act on nucleosomes to regulate transcription is the subject of considerable ongoing study. We explore the extent to which these transcriptional regulatory components function in the context of the evolutionarily ancient role of chromatin as a barrier to processes acting on DNA and how chromatin proteins have diversified to carry out evolutionarily recent functions that accompanied the emergence of differentiation and development in multicellular eukaryotes.

How the Genome Folds: The Biophysics of Four-Dimensional Chromatin Organization

The genetic information that instructs transcription and other cellular functions is carried by the chromosomes, polymers of DNA in complex with histones and other proteins. These polymers are folded inside nuclei five orders of magnitude smaller than their linear length, and many facets of this folding correlate with or are causally related to transcription and other cellular functions. Recent advances in sequencing and imaging-based techniques have enabled new views into several layers of chromatin organization. These experimental findings are accompanied by computational modeling efforts based on polymer physics that can provide mechanistic insights and quantitative predictions. Here, we review current knowledge of the main levels of chromatin organization, from the scale of nucleosomes to the entire nucleus, our current understanding of their underlying biophysical and molecular mechanisms, and some of their functional implications.

ADF1 and BEAF-32 chromatin proteins affect nucleosome positioning and DNA decompaction in 61C7/C8 interband region of Drosophila melanogaster polytene chromosomes

The formation of interphase chromosomes is a multi-level process in which DNA is compacted several thousandfold by association with histones and non-histone proteins. The first step of compaction includes the formation of nucleosomes – the basic repeating units of chromatin. Further packaging occurs due to DNA binding to histone H1 and non-histone proteins involved in enhancer-promoter and insulator interactions. Under these conditions, the genome retains its functionality due to the dynamic and uneven DNA compaction along the chromatin fiber. Since the DNA compaction level affects the transcription activity of a certain genomic region, it is important to understand the interplay between the factors acting at different levels of the packaging process. Drosophila polytene chromosomes are an excellent model system for studying the molecular mechanisms that determine DNA compaction degree. The unevenness of DNA packaging along the chromatin fiber is easily observed along these chromosomes due to their large size and specific banding pattern. The purpose of this study was to figure out the role of two non-histone regulatory proteins, ADF1 and BEAF-32, in the DNA packaging process from nucleosome positioning to the establishment of the final chromosome structure. We studied the impact of mutations that affect ADF1 and BEAF-32 binding sites on the formation of 61C7/C8 interband – one of the decompacted regions of Drosophila polytene chromosomes. We show that such mutations led to the collapse of an interband, which was accompanied with increased nucleosome stability. We also find that ADF1 and BEAF-32 binding sites are essential for the rescue of lethality caused by the null allele of bantam microRNA gene located in the region 61C7/C8.