Epigenomic signatures of cis-regulatory elements in the developing mouse and pig forelimb

Cis-regulatory elements (CREs) orchestrate the spatiotemporal regulation of key transcriptional programs. These genomic regions-including promoters, enhancers, and insulators-play a crucial role during embryonic development, and their functional diversification is thought to contribute significantly to the morphological evolution of animal body plans. We performed chromatin immunoprecipitation for various histone modifications (H3K4me3, H3K27ac, and H3K4me1) from equivalent developmental stages of mouse (E11.5) and pig (day 23; D23) to identify active regulatory regions during forelimb development in both species. The overlap of these epigenomic signatures with the pattern of open chromatin allowed us to classify these putative regulatory regions into different chromatin states in mouse and pig limb primordia. Our profiling of the regulatory genome in mouse and pig limb buds offers a valuable resource in the evo-devo field for exploring mechanisms underlying the morphological evolution of the tetrapod limb.

Systematic study of hybrid triplex topology and stability suggests a general triplex-mediated regulatory mechanism

By combining in silico, biophysical, and in vitro experiments, we decipher the topology, physical, and potential biological properties of hybrid-parallel nucleic acids triplexes, an elusive structure at the basis of life. We found that hybrid triplex topology follows a stability order: r(Py)-d(Pu)·r(Py) > r(Py)-d(Pu)·d(Py) > d(Py)-d(Pu)·d(Py) > d(Py)-d(Pu)·r(Py). The r(Py)-d(Pu)·d(Py) triplex is expected to be preferred in the cell as it avoids the need to open the duplex reducing the torsional stress required for triplex formation in the r(Py)-d(Pu)·r(Py) topology. Upon a massive collection of melting data, we have created the first predictor for hybrid triplex stability. Leveraging this predictor, we conducted a comprehensive scan to assess the likelihood of the human genome and transcriptome to engage in triplex formation. Our findings unveil a remarkable inclination-of both the human genome and transcriptome-to generate hybrid triplex formation, particularly within untranslated (UTRs) and regulatory regions, thereby corroborating the existence of a triplex-mediated regulatory mechanism. Furthermore, we found a correlation between nucleosome linkers and Triplex-forming sequence (TFS) which agree with a putative role of triplexes in arranging chromatin structure.

H3K56 acetylation regulates chromatin maturation following DNA replication

Following DNA replication, the newly reassembled chromatin is disorganized and must mature to its steady state to maintain both genome and epigenome integrity. However, the regulatory mechanisms governing this critical process remain poorly understood. Here, we show that histone H3K56 acetylation (H3K56ac), a mark on newly-synthesized H3, facilitates the remodeling of disorganized nucleosomes in nascent chromatin, and its removal at the subsequent G2/M phase of the cell cycle marks the completion of chromatin maturation. In vitro, H3K56ac enhances the activity of ISWI chromatin remodelers, including yeast ISW1 and its human equivalent SNF2h. In vivo, a deficiency of H3K56ac in nascent chromatin results in the formation of closely packed di-nucleosomes and/or tetra-nucleosomes. In contrast, abnormally high H3K56ac levels disrupt chromatin maturation, leading to genome instability. These findings establish a central role of H3K56ac in chromatin maturation and reveal a mechanism regulating this critical aspect of chromosome replication.

An integrated machine-learning model to predict nucleosome architecture

We demonstrate that nucleosomes placed in the gene body can be accurately located from signal decay theory assuming two emitters located at the beginning and at the end of genes. These generated wave signals can be in phase (leading to well defined nucleosome arrays) or in antiphase (leading to fuzzy nucleosome architectures). We found that the first (+1) and the last (-last) nucleosomes are contiguous to regions signaled by transcription factor binding sites and unusual DNA physical properties that hinder nucleosome wrapping. Based on these analyses, we developed a method that combines Machine Learning and signal transmission theory able to predict the basal locations of the nucleosomes with an accuracy similar to that of experimental MNase-seq based methods.

Comparative 3D genome analysis between neural retina and retinal pigment epithelium reveals differential cis-regulatory interactions at retinal disease loci

BACKGROUND

Vision depends on the interplay between photoreceptor cells of the neural retina and the underlying retinal pigment epithelium (RPE). Most genes involved in inherited retinal diseases display specific spatiotemporal expression within these interconnected retinal components through the local recruitment of cis-regulatory elements (CREs) in 3D nuclear space.

RESULTS

To understand the role of differential chromatin architecture in establishing tissue-specific expression at inherited retinal disease loci, we mapped genome-wide chromatin interactions using in situ Hi-C and H3K4me3 HiChIP on neural retina and RPE/choroid from human adult donor eyes. We observed chromatin looping between active promoters and 32,425 and 8060 candidate CREs in the neural retina and RPE/choroid, respectively. A comparative 3D genome analysis between these two retinal tissues revealed that 56% of 290 known inherited retinal disease genes were marked by differential chromatin interactions. One of these was ABCA4, which is implicated in the most common autosomal recessive inherited retinal disease. We zoomed in on retina- and RPE-specific cis-regulatory interactions at the ABCA4 locus using high-resolution UMI-4C. Integration with bulk and single-cell epigenomic datasets and in vivo enhancer assays in zebrafish revealed tissue-specific CREs interacting with ABCA4.

CONCLUSIONS

Through comparative 3D genome mapping, based on genome-wide, promoter-centric, and locus-specific assays of human neural retina and RPE, we have shown that gene regulation at key inherited retinal disease loci is likely mediated by tissue-specific chromatin interactions. These findings do not only provide insight into tissue-specific regulatory landscapes at retinal disease loci, but also delineate the search space for non-coding genomic variation underlying unsolved inherited retinal diseases.

Evolutionary Dynamics of Chromatin Structure and Duplicate Gene Expression in Diploid and Allopolyploid Cotton

Polyploidy is a prominent mechanism of plant speciation and adaptation, yet the mechanistic understandings of duplicated gene regulation remain elusive. Chromatin structure dynamics are suggested to govern gene regulatory control. Here, we characterized genome-wide nucleosome organization and chromatin accessibility in allotetraploid cotton, Gossypium hirsutum (AADD, 2n = 4X = 52), relative to its two diploid parents (AA or DD genome) and their synthetic diploid hybrid (AD), using DNS-seq. The larger A-genome exhibited wider average nucleosome spacing in diploids, and this intergenomic difference diminished in the allopolyploid but not hybrid. Allopolyploidization also exhibited increased accessibility at promoters genome-wide and synchronized cis-regulatory motifs between subgenomes. A prominent cis-acting control was inferred for chromatin dynamics and demonstrated by transposable element removal from promoters. Linking accessibility to gene expression patterns, we found distinct regulatory effects for hybridization and later allopolyploid stages, including nuanced establishment of homoeolog expression bias and expression level dominance. Histone gene expression and nucleosome organization are coordinated through chromatin accessibility. Our study demonstrates the capability to track high-resolution chromatin structure dynamics and reveals their role in the evolution of cis-regulatory landscapes and duplicate gene expression in polyploids, illuminating regulatory ties to subgenomic asymmetry and dominance.

Skeletal muscle differentiation induces wide-ranging nucleosome repositioning in muscle gene promoters

In a previous report, we demonstrated that Cbx1, PurB and Sp3 inhibited cardiac muscle differentiation by increasing nucleosome density around cardiac muscle gene promoters. Since cardiac and skeletal muscle express many of the same proteins, we asked if Cbx1, PurB and Sp3 similarly regulated skeletal muscle differentiation. In a C2C12 model of skeletal muscle differentiation, Cbx1 and PurB knockdown increased myotube formation. In contrast, Sp3 knockdown inhibited myotube formation, suggesting that Sp3 played opposing roles in cardiac muscle and skeletal muscle differentiation. Consistent with this finding, Sp3 knockdown also inhibited various muscle-specific genes. The Cbx1, PurB and Sp3 proteins are believed to influence gene-expression in part by altering nucleosome position. Importantly, we developed a statistical approach to determine if changes in nucleosome positioning were significant and applied it to understanding the architecture of muscle-specific genes. Through this novel statistical approach, we found that during myogenic differentiation, skeletal muscle-specific genes undergo a set of unique nucleosome changes which differ significantly from those shown in commonly expressed muscle genes. While Sp3 binding was associated with nucleosome loss, there appeared no correlation with the aforementioned nucleosome changes. In summary, we have identified a novel role for Sp3 in skeletal muscle differentiation and through the application of quantifiable MNase-seq have discovered unique fingerprints of nucleosome changes for various classes of muscle genes during myogenic differentiation.

Automated Composition Assessment of Natural Extracts: Untargeted Mass Spectrometry-Based Metabolite Profiling Integrating Semiquantitative Detection

Recent developments in mass spectrometry-based metabolite profiling allow unprecedented qualitative coverage of complex biological extract composition. However, the electrospray ionization used in metabolite profiling generates multiple artifactual signals for a single analyte. This leads to thousands of signals per analysis without satisfactory means of filtering those corresponding to abundant constituents. Generic approaches are therefore needed for the qualitative and quantitative annotation of a broad range of relevant constituents. For this, we used an analytical platform combining liquid chromatography-mass spectrometry (LC-MS) with Charged Aerosol Detection (CAD). We established a generic metabolite profiling for the concomitant recording of qualitative MS data and semiquantitative CAD profiles. The MS features (recorded in high-resolution tandem MS) are grouped and annotated using state-of-the-art tools. To efficiently attribute features to their corresponding extracted and integrated CAD peaks, a custom signal pretreatment and peak-shape comparison workflow is built. This strategy allows us to automatically contextualize features at both major and minor metabolome levels, together with a detailed reporting of their annotation including relevant orthogonal information (taxonomy, retention time). Signals not attributed to CAD peaks are considered minor metabolites. Results are illustrated on an ethanolic extract of Swertia chirayita (Roxb.) H. Karst., a bitter plant of industrial interest, exhibiting the typical complexity of plant extracts as a proof of concept. This generic qualitative and quantitative approach paves the way to automatically assess the composition of single natural extracts of interest or broader collections, thus facilitating new ingredient registrations or natural-extracts-based drug discovery campaigns.

Genome-wide chromatin interaction map for Trypanosoma cruzi

Trypanosomes are eukaryotic, unicellular parasites, such as Trypanosoma brucei, which causes sleeping sickness, and Trypanosoma cruzi, which causes Chagas disease. Genomes of these parasites comprise core regions and species-specific disruptive regions that encode multigene families of surface glycoproteins. Few transcriptional regulators have been identified in these parasites, and the role of spatial organization of the genome in gene expression is unclear. Here we mapped genome-wide chromatin interactions in T. cruzi using chromosome conformation capture (Hi-C), and we show that the core and disruptive regions form three-dimensional chromatin compartments named C and D. These chromatin compartments differ in levels of DNA methylation, nucleosome positioning and chromatin interactions, affecting genome expression dynamics. Our data reveal that the trypanosome genome is organized into chromatin-folding domains and transcription is affected by the local chromatin structure. We propose a model in which epigenetic mechanisms affect gene expression in trypanosomes.

Biased eviction of variant histone H3 nucleosomes triggers biofilm growth in Candida albicans

IMPORTANCE

Candida albicans lives as a commensal in most healthy humans but can cause superficial skin infections to life-threatening systemic infections. C. albicans also forms biofilms on biotic and abiotic surfaces. Biofilm cells are difficult to treat and highly resistant to antifungals. A specific set of genes is differentially regulated in biofilm cells as compared to free-floating planktonic cells of C. albicans. In this study, we addressed how a variant histone H3VCTG, a previously identified negative regulator of biofilm formation, modulates gene expression changes. By providing compelling evidence, we show that biased eviction of H3VCTG nucleosomes at the promoters of biofilm-relevant genes facilitates the accessibility of both transcription activators and repressors to modulate gene expression. Our study is a comprehensive investigation of genome-wide nucleosome occupancy in both planktonic and biofilm states, which reveals transition to an open chromatin landscape during biofilm mode of growth in C. albicans, a medically relevant pathogen.

Changes in adenoviral chromatin organization precede early gene activation upon infection

Within the virion, adenovirus DNA associates with the virus-encoded, protamine-like structural protein pVII. Whether this association is organized, and how genome packaging changes during infection and subsequent transcriptional activation is currently unclear. Here, we combined RNA-seq, MNase-seq, ChIP-seq, and single genome imaging during early adenovirus infection to unveil the structure- and time-resolved dynamics of viral chromatin changes as well as their correlation with gene transcription. Our MNase mapping data indicates that the adenoviral genome is arranged in precisely positioned nucleoprotein particles with nucleosome-like characteristics, that we term adenosomes. We identified 238 adenosomes that are positioned by a DNA sequence code and protect about 60-70 bp of DNA. The incoming adenoviral genome is more accessible at early gene loci that undergo additional chromatin de-condensation upon infection. Histone H3.3 containing nucleosomes specifically replaces pVII at distinct genomic sites and at the transcription start sites of early genes. Acetylation of H3.3 is predominant at the transcription start sites and precedes transcriptional activation. Based on our results, we propose a central role for the viral pVII nucleoprotein architecture, which is required for the dynamic structural changes during early infection, including the regulation of nucleosome assembly prior to transcription initiation. Our study thus may aid the rational development of recombinant adenoviral vectors exhibiting sustained expression in gene therapy.

Histone H3 serine-57 is a CHK1 substrate whose phosphorylation affects DNA repair

Histone post-translational modifications promote a chromatin environment that controls transcription, DNA replication and repair, but surprisingly few phosphorylations have been documented. We report the discovery of histone H3 serine-57 phosphorylation (H3S57ph) and show that it is implicated in different DNA repair pathways from fungi to vertebrates. We identified CHK1 as a major human H3S57 kinase, and disrupting or constitutively mimicking H3S57ph had opposing effects on rate of recovery from replication stress, 53BP1 chromatin binding, and dependency on RAD52. In fission yeast, mutation of all H3 alleles to S57A abrogated DNA repair by both non-homologous end-joining and homologous recombination, while cells with phospho-mimicking S57D alleles were partly compromised for both repair pathways, presented aberrant Rad52 foci and were strongly sensitised to replication stress. Mechanistically, H3S57ph loosens DNA-histone contacts, increasing nucleosome mobility, and interacts with H3K56. Our results suggest that dynamic phosphorylation of H3S57 is required for DNA repair and recovery from replication stress, opening avenues for investigating the role of this modification in other DNA-related processes.

The esBAF and ISWI nucleosome remodeling complexes influence occupancy of overlapping dinucleosomes and fragile nucleosomes in murine embryonic stem cells

BACKGROUND

Nucleosome remodeling factors regulate the occupancy and positioning of nucleosomes genome-wide through ATP-driven DNA translocation. While many nucleosomes are consistently well-positioned, some nucleosomes and alternative nucleosome structures are more sensitive to nuclease digestion or are transitory. Fragile nucleosomes are nucleosome structures that are sensitive to nuclease digestion and may be composed of either six or eight histone proteins, making these either hexasomes or octasomes. Overlapping dinucleosomes are composed of two merged nucleosomes, lacking one H2A:H2B dimer, creating a 14-mer wrapped by ~ 250 bp of DNA. In vitro studies of nucleosome remodeling suggest that the collision of adjacent nucleosomes by sliding stimulates formation of overlapping dinucleosomes.

RESULTS

To better understand how nucleosome remodeling factors regulate alternative nucleosome structures, we depleted murine embryonic stem cells of the transcripts encoding remodeler ATPases BRG1 or SNF2H, then performed MNase-seq. We used high- and low-MNase digestion to assess the effects of nucleosome remodeling factors on nuclease-sensitive or "fragile" nucleosome occupancy. In parallel we gel-extracted MNase-digested fragments to enrich for overlapping dinucleosomes. We recapitulate prior identification of fragile nucleosomes and overlapping dinucleosomes near transcription start sites, and identify enrichment of these features around gene-distal DNaseI hypersensitive sites, CTCF binding sites, and pluripotency factor binding sites. We find that BRG1 stimulates occupancy of fragile nucleosomes but restricts occupancy of overlapping dinucleosomes.

CONCLUSIONS

Overlapping dinucleosomes and fragile nucleosomes are prevalent within the ES cell genome, occurring at hotspots of gene regulation beyond their characterized existence at promoters. Although neither structure is fully dependent on either nucleosome remodeling factor, both fragile nucleosomes and overlapping dinucleosomes are affected by knockdown of BRG1, suggesting a role for the complex in creating or removing these structures.

Ageing-associated changes in transcriptional elongation influence longevity

Physiological homeostasis becomes compromised during ageing, as a result of impairment of cellular processes, including transcription and RNA splicing1-4. However, the molecular mechanisms leading to the loss of transcriptional fidelity are so far elusive, as are ways of preventing it. Here we profiled and analysed genome-wide, ageing-related changes in transcriptional processes across different organisms: nematodes, fruitflies, mice, rats and humans. The average transcriptional elongation speed (RNA polymerase II speed) increased with age in all five species. Along with these changes in elongation speed, we observed changes in splicing, including a reduction of unspliced transcripts and the formation of more circular RNAs. Two lifespan-extending interventions, dietary restriction and lowered insulin-IGF signalling, both reversed most of these ageing-related changes. Genetic variants in RNA polymerase II that reduced its speed in worms5 and flies6 increased their lifespan. Similarly, reducing the speed of RNA polymerase II by overexpressing histone components, to counter age-associated changes in nucleosome positioning, also extended lifespan in flies and the division potential of human cells. Our findings uncover fundamental molecular mechanisms underlying animal ageing and lifespan-extending interventions, and point to possible preventive measures.

MiOS, an integrated imaging and computational strategy to model gene folding with nucleosome resolution

The linear sequence of DNA provides invaluable information about genes and their regulatory elements along chromosomes. However, to fully understand gene function and regulation, we need to dissect how genes physically fold in the three-dimensional nuclear space. Here we describe immuno-OligoSTORM, an imaging strategy that reveals the distribution of nucleosomes within specific genes in super-resolution, through the simultaneous visualization of DNA and histones. We combine immuno-OligoSTORM with restraint-based and coarse-grained modeling approaches to integrate super-resolution imaging data with Hi-C contact frequencies and deconvoluted micrococcal nuclease-sequencing information. The resulting method, called Modeling immuno-OligoSTORM, allows quantitative modeling of genes with nucleosome resolution and provides information about chromatin accessibility for regulatory factors, such as RNA polymerase II. With Modeling immuno-OligoSTORM, we explore intercellular variability, transcriptional-dependent gene conformation, and folding of housekeeping and pluripotency-related genes in human pluripotent and differentiated cells, thereby obtaining the highest degree of data integration achieved so far to our knowledge.

Early in an SV40 infection, histone modifications correlate with the presence or absence of RNAPII and direction of transcription

Since epigenetic regulation seemed likely to be involved in SV40 early transcription following infection, we have analyzed the organization of nucleosomes carrying histone modifications (acetyl-H3, acetyl-H4, H3K9me1, H3K9me3, H3K4me1, H3K4me3, H3K27me3, H4K20me1) at 30 min and 2 h post infection in SV40 minichromosomes prepared in the absence or presence of the transcription inhibitor dichloro-1-beta-d-ribofuranosyl benzimidazole. The former condition was used to determine how SV40 chromatin structure changed during early transcription, and the latter was used to determine the role of active transcription. The location of RNAPII was used as a marker to identify where histone modifications were most likely to be involved in regulation. Acetyl-H3 acted like epigenetic memory by being present at sites subsequently bound by RNAPII, while H3K9me1 and H3K27me3 were reorganized to the late side of the SV40 regulatory region apparently to repress late transcription. The organization of acetyl-H3 and H3K9me1 but not H3K27me3 required active transcription.

NucleoMap: A computational tool for identifying nucleosomes in ultra-high resolution contact maps

Although poorly positioned nucleosomes are ubiquitous in the eukaryotic genome, they are difficult to identify with existing nucleosome identification methods. Recently available enhanced high-throughput chromatin conformation capture techniques such as Micro-C, DNase Hi-C, and Hi-CO characterize nucleosome-level chromatin proximity, probing the positions of mono-nucleosomes and the spacing between nucleosome pairs at the same time, enabling nucleosome profiling in poorly positioned regions. Here we develop a novel computational approach, NucleoMap, to identify nucleosome positioning from ultra-high resolution chromatin contact maps. By integrating nucleosome read density, contact distances, and binding preferences, NucleoMap precisely locates nucleosomes in both prokaryotic and eukaryotic genomes and outperforms existing nucleosome identification methods in both precision and recall. We rigorously characterize genome-wide association in eukaryotes between the spatial organization of mono-nucleosomes and their corresponding histone modifications, protein binding activities, and higher-order chromatin functions. We also find evidence of two tetra-nucleosome folding structures in human embryonic stem cells and analyze their association with multiple structural and functional regions. Based on the identified nucleosomes, nucleosome contact maps are constructed, reflecting the inter-nucleosome distances and preserving the contact distance profiles in original contact maps.

Inter-nucleosomal potentials from nucleosomal positioning data

No systematic method exists to derive inter-nucleosomal potentials between nucleosomes along a chromosome consistently across a given genome. Such potentials can yield information on nucleosomal ordering, thermal as well as mechanical properties of chromosomes. Thus, indirectly, they shed light on a possible mechanical genomic code along a chromosome. To develop a method yielding effective inter-nucleosomal potentials between nucleosomes, a generalized Lennard-Jones potential for the parameterization is developed based on nucleosomal positioning data. This approach eliminates some of the problems that the underlying nucleosomal positioning data have, rendering the extraction difficult on the individual nucleosomal level. Furthermore, patterns on which to base a classification along a chromosome appear on larger domains, such as hetero- and euchromatin. An intuitive selection strategy for the noisy optimization problem is employed to derive effective exponents for the generalized potential. The method is tested on the Candida albicans genome. Applying k-means clustering based on potential parameters and thermodynamic compressibilities, a genome-wide clustering of nucleosome sequences is obtained for C. albicans. This clustering shows that a chromosome beyond the classical dichotomic categories of hetero- and euchromatin is more feature-rich.

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

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

Mapping nucleosome and chromatin architectures: A survey of computational methods

With ever-growing genomic sequencing data, the data variabilities and the underlying biases of the sequencing technologies pose significant computational challenges ranging from the need for accurately detecting the nucleosome positioning or chromatin interaction to the need for developing normalization methods to eliminate systematic biases. This review mainly surveys the computational methods for mapping the higher-resolution nucleosome and higher-order chromatin architectures. While a detailed discussion of the underlying algorithms is beyond the scope of our survey, we have discussed the methods and tools that can detect the nucleosomes in the genome, then demonstrated the computational methods for identifying 3D chromatin domains and interactions. We further illustrated computational approaches for integrating multi-omics data with Hi-C data and the advance of single-cell (sc)Hi-C data analysis. Our survey provides a comprehensive and valuable resource for biomedical scientists interested in studying nucleosome organization and chromatin structures as well as for computational scientists who are interested in improving upon them.

DeNOPA: decoding nucleosome positions sensitively with sparse ATAC-seq data

As the basal bricks, the dynamics and arrangement of nucleosomes orchestrate the higher architecture of chromatin in a fundamental way, thereby affecting almost all nuclear biology processes. Thanks to its rather simple protocol, assay for transposase-accessible chromatin using sequencing (ATAC)-seq has been rapidly adopted as a major tool for chromatin-accessible profiling at both bulk and single-cell levels; however, to picture the arrangement of nucleosomes per se remains a challenge with ATAC-seq. In the present work, we introduce a novel ATAC-seq analysis toolkit, named decoding nucleosome organization profile based on ATAC-seq data (deNOPA), to predict nucleosome positions. Assessments showed that deNOPA outperformed state-of-the-art tools with ultra-sparse ATAC-seq data, e.g. no more than 0.5 fragment per base pair. The remarkable performance of deNOPA was fueled by the short fragment reads, which compose nearly half of sequenced reads in the ATAC-seq libraries and are commonly discarded by state-of-the-art nucleosome positioning tools. However, we found that the short fragment reads enrich information on nucleosome positions and that the linker regions were predicted by reads from both short and long fragments using Gaussian smoothing. Last, using deNOPA, we showed that the dynamics of nucleosome organization may not directly couple with chromatin accessibility in the cis-regulatory regions when human cells respond to heat shock stimulation. Our deNOPA provides a powerful tool with which to analyze the dynamics of chromatin at nucleosome position level with ultra-sparse ATAC-seq data.

Genome-wide mapping reveals R-loops associated with centromeric repeats in maize

R-loops are stable chromatin structures comprising a DNA:RNA hybrid and a displaced single-stranded DNA. R-loops have been implicated in gene expression and chromatin structure, as well as in replication blocks and genome instability. Here, we conducted a genome-wide identification of R-loops and identified more than 700,000 R-loop peaks in the maize (Zea mays) genome. We found that sense R-loops were mainly enriched in promoters and transcription termination sites and relatively less enriched in gene bodies, which is different from the main gene-body localization of sense R-loops in Arabidopsis and Oryza sativa. At the chromosome scale, maize R-loops were enriched in pericentromeric heterochromatin regions, and a significant portion of R-loops were derived from transposable elements. In centromeres, R-loops preferentially formed within the binding regions of the centromere-specific histone CENH3, and centromeric retrotransposons were strongly associated with R-loop formation. Furthermore, centromeric retrotransposon R-loops were observed by applying the single-molecule imaging technique of atomic force microscopy. These findings elucidate the fundamental character of R-loops in the maize genome and reveal the potential role of R-loops in centromeres.

Impact of DNA methylation on 3D genome structure

Determining the effect of DNA methylation on chromatin structure and function in higher organisms is challenging due to the extreme complexity of epigenetic regulation. We studied a simpler model system, budding yeast, that lacks DNA methylation machinery making it a perfect model system to study the intrinsic role of DNA methylation in chromatin structure and function. We expressed the murine DNA methyltransferases in Saccharomyces cerevisiae and analyzed the correlation between DNA methylation, nucleosome positioning, gene expression and 3D genome organization. Despite lacking the machinery for positioning and reading methylation marks, induced DNA methylation follows a conserved pattern with low methylation levels at the 5’ end of the gene increasing gradually toward the 3’ end, with concentration of methylated DNA in linkers and nucleosome free regions, and with actively expressed genes showing low and high levels of methylation at transcription start and terminating sites respectively, mimicking the patterns seen in mammals. We also see that DNA methylation increases chromatin condensation in peri-centromeric regions, decreases overall DNA flexibility, and favors the heterochromatin state. Taken together, these results demonstrate that methylation intrinsically modulates chromatin structure and function even in the absence of cellular machinery evolved to recognize and process the methylation signal.

Multi-layered epigenetic regulation in higher eukaryotes makes it challenging to disentangle the individual effects of modifications on chromatin structure and function. Here, the authors expressed mammalian DNA methyltransferases in yeast, which have no DNA methylation, to show that methylation has intrinsic effects on chromatin structure.

Role of H1 and DNA methylation in selective regulation of transposable elements during heat stress

Heat-stressed Arabidopsis plants release heterochromatin-associated transposable element (TE) silencing, yet it is not accompanied by major reductions of epigenetic repressive modifications. In this study, we explored the functional role of histone H1 in repressing heterochromatic TEs in response to heat stress. We generated and analyzed RNA and bisulfite-sequencing data of wild-type and h1 mutant seedlings before and after heat stress. Loss of H1 caused activation of pericentromeric Gypsy elements upon heat treatment, despite these elements remaining highly methylated. By contrast, nonpericentromeric Copia elements became activated concomitantly with loss of DNA methylation. The same Copia elements became activated in heat-treated chromomethylase 2 (cmt2) mutants, indicating that H1 represses Copia elements through maintaining DNA methylation under heat. We discovered that H1 is required for TE repression in response to heat stress, but its functional role differs depending on TE location. Strikingly, H1-deficient plants treated with the DNA methyltransferase inhibitor zebularine were highly tolerant to heat stress, suggesting that both H1 and DNA methylation redundantly suppress the plant response to heat stress.

Genome-wide prediction of topoisomerase IIβ binding by architectural factors and chromatin accessibility

DNA topoisomerase II-β (TOP2B) is fundamental to remove topological problems linked to DNA metabolism and 3D chromatin architecture, but its cut-and-reseal catalytic mechanism can accidentally cause DNA double-strand breaks (DSBs) that can seriously compromise genome integrity. Understanding the factors that determine the genome-wide distribution of TOP2B is therefore not only essential for a complete knowledge of genome dynamics and organization, but also for the implications of TOP2-induced DSBs in the origin of oncogenic translocations and other types of chromosomal rearrangements. Here, we conduct a machine-learning approach for the prediction of TOP2B binding using publicly available sequencing data. We achieve highly accurate predictions, with accessible chromatin and architectural factors being the most informative features. Strikingly, TOP2B is sufficiently explained by only three features: DNase I hypersensitivity, CTCF and cohesin binding, for which genome-wide data are widely available. Based on this, we develop a predictive model for TOP2B genome-wide binding that can be used across cell lines and species, and generate virtual probability tracks that accurately mirror experimental ChIP-seq data. Our results deepen our knowledge on how the accessibility and 3D organization of chromatin determine TOP2B function, and constitute a proof of principle regarding the in silico prediction of sequence-independent chromatin-binding factors.

Differential SP1 interactions in SV40 chromatin from virions and minichromosomes

SP1 binding in SV40 chromatin in vitro and in vivo was characterized in order to better understand its role during the initiation of early transcription. We observed that chromatin from disrupted virions, but not minichromosomes, was efficiently bound by HIS-tagged SP1 in vitro, while the opposite was true for the presence of endogenous SP1 introduced in vivo. Using ChIP-Seq to compare the location of SP1 to nucleosomes carrying modified histones, we found that SP1 could occupy its whole binding site in virion chromatin but only the early side of its binding site in most of the minichromosomes carrying modified histones due to the presence of overlapping nucleosomes. The results suggest that during the initiation of an SV40 infection, SP1 binds to an open region in SV40 virion chromatin but quickly triggers chromatin reorganization and its own removal.

Q-Nuc: a bioinformatics pipeline for the quantitative analysis of nucleosomal profiles

Nucleosomal profiling is an effective method to determine the positioning and occupancy of nucleosomes, which is essential to understand their roles in genomic processes. However, the positional randomness across the genome and its relationship with nucleosome occupancy remains poorly understood. Here we present a computational method that segments the profile into nucleosomal domains and quantifies their randomness and relative occupancy level. Applying this method to published data, we find on average ~ 3-fold differences in the degree of positional randomness between regions typically considered "well-ordered", as well as an unexpected predominance of only two types of domains of positional randomness in yeast cells. Further, we find that occupancy levels between domains actually differ maximally by ~ 2-3-fold in both cells, which has not been described before. We also developed a procedure by which one can estimate the sequencing depth that is required to identify nucleosomal positions even when regional positional randomness is high. Overall, we have developed a pipeline to quantitatively characterize domain-level features of nucleosome randomness and occupancy genome-wide, enabling the identification of otherwise unknown features in nucleosomal organization.

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.

DNA stretching induces Cas9 off-target activity

CRISPR/Cas9 is a powerful genome-editing tool, but spurious off-target edits present a barrier to therapeutic applications. To understand how CRISPR/Cas9 discriminates between on-targets and off-targets, we have developed a single-molecule assay combining optical tweezers with fluorescence to monitor binding to λ-DNA. At low forces, the Streptococcus pyogenes Cas9 complex binds and cleaves DNA specifically. At higher forces, numerous off-target binding events appear repeatedly at the same off-target sites in a guide-RNA-sequence-dependent manner, driven by the mechanical distortion of the DNA. Using single-molecule Förster resonance energy transfer (smFRET) and cleavage assays, we show that DNA bubbles induce off-target binding and cleavage at these sites, even with ten mismatches, as well as at previously identified in vivo off-targets. We propose that duplex DNA destabilization during cellular processes (for example, transcription, replication, etc.) can expose these cryptic off-target sites to Cas9 activity, highlighting the need for improved off-target prediction algorithms.

Application of MNase-Seq in the Global Mapping of Nucleosome Positioning in Plants

The precise positioning of nucleosomes along the underlying DNA is critical for a variety of biological processes, especially in regulating transcription. The interplay between nucleosomes and transcription factors for accessing the underlying DNA sequences is one of the key determinants that affect transcriptional regulation. Moreover, nucleosomes with various packing statuses confer distinct functions in regulating gene expressions in response to various internal or external signals. Therefore, global mapping of nucleosome positions is one informative way to elucidate the relationship between patterns of nucleosome positioning/occupancy and transcriptional regulations. MNase digestion coupled with high-throughput sequencing (MNase-seq) has been utilized widely for global mapping of nucleosome positioning in eukaryotes that have a sequenced genome. We have developed a robust MNase-seq procedure in plants. It mainly includes plant nuclei isolation, treatment of purified nuclei with MNase, gel recovery of MNase-trimmed mononucleosomal DNA with an approximate size of 150 bp, MNase-seq library preparation followed by Illumina sequencing, and data analysis. MNase-seq has already been successfully applied to identify genome-wide nucleosome positioning in model plants, rice, and Arabidopsis thaliana.

H3K23me1 is an evolutionarily conserved histone modification associated with CG DNA methylation in Arabidopsis

Amino-terminal tails of histones are targets for diverse post-translational modifications whose combinatorial action may constitute a code that will be read and interpreted by cellular proteins to define particular transcriptional states. Here, we describe monomethylation of histone H3 lysine 23 (H3K23me1) as a histone modification not previously described in plants. H3K23me1 is an evolutionarily conserved mark in diverse species of flowering plants. Chromatin immunoprecipitation followed by high-throughput sequencing in Arabidopsis thaliana showed that H3K23me1 was highly enriched in pericentromeric regions and depleted from chromosome arms. In transposable elements it co-localized with CG, CHG and CHH DNA methylation as well as with the heterochromatic histone mark H3K9me2. Transposable elements are often rich in H3K23me1 but different families vary in their enrichment: LTR-Gypsy elements are most enriched and RC/Helitron elements are least enriched. The histone methyltransferase KRYPTONITE and normal DNA methylation were required for normal levels of H3K23me1 on transposable elements. Immunostaining experiments confirmed the pericentromeric localization and also showed mild enrichment in less condensed regions. Accordingly, gene bodies of protein-coding genes had intermediate H3K23me1 levels, which coexisted with CG DNA methylation. Enrichment of H3K23me1 along gene bodies did not correlate with transcription levels. Together, this work establishes H3K23me1 as a so far undescribed component of the plant histone code.

Nucleosomes influence multiple steps during replication initiation

Eukaryotic replication origin licensing, activation and timing are influenced by chromatin but a mechanistic understanding is lacking. Using reconstituted nucleosomal DNA replication assays, we assessed the impact of nucleosomes on replication initiation. To generate distinct nucleosomal landscapes, different chromatin-remodeling enzymes (CREs) were used to remodel nucleosomes on origin-DNA templates. Nucleosomal organization influenced two steps of replication initiation: origin licensing and helicase activation. Origin licensing assays showed that local nucleosome positioning enhanced origin specificity and modulated helicase loading by influencing ORC DNA binding. Interestingly, SWI/SNF- and RSC-remodeled nucleosomes were permissive for origin licensing but showed reduced helicase activation. Specific CREs rescued replication of these templates if added prior to helicase activation, indicating a permissive chromatin state must be established during origin licensing to allow efficient origin activation. Our studies show nucleosomes directly modulate origin licensing and activation through distinct mechanisms and provide insights into the regulation of replication initiation by chromatin.

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.

Using informative Multinomial-Dirichlet prior in a t-mixture with reversible jump estimation of nucleosome positions for genome-wide profiling

Genome-wide mapping of nucleosomes has revealed a great deal about the relationships between chromatin structure and control of gene expression. Recent next generation CHIP-chip and CHIP-Seq technologies have accelerated our understanding of basic principles of chromatin organization. These technologies have taught us that nucleosomes play a crucial role in gene regulation by allowing physical access to transcription factors. Recent methods and experimental advancements allow the determination of nucleosome positions for a given genome area. However, most of these methods estimate the number of nucleosomes either by an EM algorithm using a BIC criterion or an effective heuristic strategy. Here, we introduce a Bayesian method for identifying nucleosome positions. The proposed model is based on a Multinomial-Dirichlet classification and a hierarchical mixture distributions. The number and the positions of nucleosomes are estimated using a reversible jump Markov chain Monte Carlo simulation technique. We compare the performance of our method on simulated data and MNase-Seq data from Saccharomyces cerevisiae against PING and NOrMAL methods.

Recombination Rate Heterogeneity within Arabidopsis Disease Resistance Genes

Meiotic crossover frequency varies extensively along chromosomes and is typically concentrated in hotspots. As recombination increases genetic diversity, hotspots are predicted to occur at immunity genes, where variation may be beneficial. A major component of plant immunity is recognition of pathogen Avirulence (Avr) effectors by resistance (R) genes that encode NBS-LRR domain proteins. Therefore, we sought to test whether NBS-LRR genes would overlap with meiotic crossover hotspots using experimental genetics in Arabidopsis thaliana. NBS-LRR genes tend to physically cluster in plant genomes; for example, in Arabidopsis most are located in large clusters on the south arms of chromosomes 1 and 5. We experimentally mapped 1,439 crossovers within these clusters and observed NBS-LRR gene associated hotspots, which were also detected as historical hotspots via analysis of linkage disequilibrium. However, we also observed NBS-LRR gene coldspots, which in some cases correlate with structural heterozygosity. To study recombination at the fine-scale we used high-throughput sequencing to analyze ~1,000 crossovers within the RESISTANCE TO ALBUGO CANDIDA1 (RAC1) R gene hotspot. This revealed elevated intragenic crossovers, overlapping nucleosome-occupied exons that encode the TIR, NBS and LRR domains. The highest RAC1 recombination frequency was promoter-proximal and overlapped CTT-repeat DNA sequence motifs, which have previously been associated with plant crossover hotspots. Additionally, we show a significant influence of natural genetic variation on NBS-LRR cluster recombination rates, using crosses between Arabidopsis ecotypes. In conclusion, we show that a subset of NBS-LRR genes are strong hotspots, whereas others are coldspots. This reveals a complex recombination landscape in Arabidopsis NBS-LRR genes, which we propose results from varying coevolutionary pressures exerted by host-pathogen relationships, and is influenced by structural heterozygosity.

DChIPRep, an R/Bioconductor package for differential enrichment analysis in chromatin studies

The genome-wide study of epigenetic states requires the integrative analysis of histone modification ChIP-seq data. Here, we introduce an easy-to-use analytic framework to compare profiles of enrichment in histone modifications around classes of genomic elements, e.g. transcription start sites (TSS). Our framework is available via the user-friendly R/Bioconductor package DChIPRep. DChIPRep uses biological replicate information as well as chromatin Input data to allow for a rigorous assessment of differential enrichment. DChIPRep is available for download through the Bioconductor project at http://bioconductor.org/packages/DChIPRep. Contact. DChIPRep@gmail.com.

Nucleosome architecture throughout the cell cycle

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

Gene Expression and Chromatin Modifications Associated with Maize Centromeres

Centromeres are defined by the presence of CENH3, a variant of histone H3. Centromeres in most plant species contain exclusively highly repetitive DNA sequences, which has hindered research on structure and function of centromeric chromatin. Several maize centromeres have been nearly completely sequenced, providing a sequence-based platform for genomic and epigenomic research of plant centromeres. Here we report a high resolution map of CENH3 nucleosomes in the maize genome. Although CENH3 nucleosomes are spaced ∼190 bp on average, CENH3 nucleosomes that occupied CentC, a 156-bp centromeric satellite repeat, showed clear positioning aligning with CentC monomers. Maize centromeres contain alternating CENH3-enriched and CENH3-depleted subdomains, which account for 87% and 13% of the centromeres, respectively. A number of annotated genes were identified in the centromeres, including 11 active genes that were located exclusively in CENH3-depleted subdomains. The euchromatic histone modification marks, including H3K4me3, H3K36me3 and H3K9ac, detected in maize centromeres were associated mainly with the active genes. Interestingly, maize centromeres also have lower levels of the heterochromatin histone modification mark H3K27me2 relative to pericentromeric regions. We conclude that neither H3K27me2 nor the three euchromatic histone modifications are likely to serve as functionally important epigenetic marks of centromere identity in maize.

Nucleosome positioning: resources and tools online

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

Zelda overcomes the high intrinsic nucleosome barrier at enhancers during Drosophila zygotic genome activation

The Drosophila genome activator Vielfaltig (Vfl), also known as Zelda (Zld), is thought to prime enhancers for activation by patterning transcription factors (TFs). Such priming is accompanied by increased chromatin accessibility, but the mechanisms by which this occurs are poorly understood. Here, we analyze the effect of Zld on genome-wide nucleosome occupancy and binding of the patterning TF Dorsal (Dl). Our results show that early enhancers are characterized by an intrinsically high nucleosome barrier. Zld tackles this nucleosome barrier through local depletion of nucleosomes with the effect being dependent on the number and position of Zld motifs. Without Zld, Dl binding decreases at enhancers and redistributes to open regions devoid of enhancer activity. We propose that Zld primes enhancers by lowering the high nucleosome barrier just enough to assist TFs in accessing their binding motifs and promoting spatially controlled enhancer activation if the right patterning TFs are present. We envision that genome activators in general will utilize this mechanism to activate the zygotic genome in a robust and precise manner.

DNA sequence templates adjacent nucleosome and ORC sites at gene amplification origins in Drosophila

Eukaryotic origins of DNA replication are bound by the origin recognition complex (ORC), which scaffolds assembly of a pre-replicative complex (pre-RC) that is then activated to initiate replication. Both pre-RC assembly and activation are strongly influenced by developmental changes to the epigenome, but molecular mechanisms remain incompletely defined. We have been examining the activation of origins responsible for developmental gene amplification in Drosophila. At a specific time in oogenesis, somatic follicle cells transition from genomic replication to a locus-specific replication from six amplicon origins. Previous evidence indicated that these amplicon origins are activated by nucleosome acetylation, but how this affects origin chromatin is unknown. Here, we examine nucleosome position in follicle cells using micrococcal nuclease digestion with Ilumina sequencing. The results indicate that ORC binding sites and other essential origin sequences are nucleosome-depleted regions (NDRs). Nucleosome position at the amplicons was highly similar among developmental stages during which ORC is or is not bound, indicating that being an NDR is not sufficient to specify ORC binding. Importantly, the data suggest that nucleosomes and ORC have opposite preferences for DNA sequence and structure. We propose that nucleosome hyperacetylation promotes pre-RC assembly onto adjacent DNA sequences that are disfavored by nucleosomes but favored by ORC.

Nucleosome positioning in yeasts: methods, maps, and mechanisms

Eukaryotic nuclear DNA is packaged into nucleosomes. During the past decade, genome-wide nucleosome mapping across species revealed the high degree of order in nucleosome positioning. There is a conserved stereotypical nucleosome organization around transcription start sites (TSSs) with a nucleosome-depleted region (NDR) upstream of the TSS and a TSS-aligned regular array of evenly spaced nucleosomes downstream over the gene body. As nucleosomes largely impede access to DNA and thereby provide an important level of genome regulation, it is of general interest to understand the mechanisms generating nucleosome positioning and especially the stereotypical NDR-array pattern. We focus here on the most advanced models, unicellular yeasts, and review the progress in mapping nucleosomes and which nucleosome positioning mechanisms are discussed. There are four mechanistic aspects: How are NDRs generated? How are individual nucleosomes positioned, especially those flanking the NDRs? How are nucleosomes evenly spaced leading to regular arrays? How are regular arrays aligned at TSSs? The main candidates for nucleosome positioning determinants are intrinsic DNA binding preferences of the histone octamer, specific DNA binding factors, nucleosome remodeling enzymes, transcription, and statistical positioning. We summarize the state of the art in an integrative model where nucleosomes are positioned by a combination of all these candidate determinants. We highlight the predominance of active mechanisms involving nucleosome remodeling enzymes which may be recruited by DNA binding factors and the transcription machinery. While this mechanistic framework emerged clearly during recent years, the involved factors and their mechanisms are still poorly understood and require future efforts combining in vivo and in vitro approaches.

Comparative analysis of methods for genome-wide nucleosome cartography

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

PuFFIN--a parameter-free method to build nucleosome maps from paired-end reads

Background

We introduce a novel method, called PuFFIN, that takes advantage of paired-end short reads to build genome-wide nucleosome maps with larger numbers of detected nucleosomes and higher accuracy than existing tools. In contrast to other approaches that require users to optimize several parameters according to their data (e.g., the maximum allowed nucleosome overlap or legal ranges for the fragment sizes) our algorithm can accurately determine a genome-wide set of non-overlapping nucleosomes without any user-defined parameter. This feature makes PuFFIN significantly easier to use and prevents users from choosing the "wrong" parameters and obtain sub-optimal nucleosome maps.

Results

PuFFIN builds genome-wide nucleosome maps using a multi-scale (or multi-resolution) approach. Our algorithm relies on a set of nucleosome "landscape" functions at different resolution levels: each function represents the likelihood of each genomic location to be occupied by a nucleosome for a particular value of the smoothing parameter. After a set of candidate nucleosomes is computed for each function, PuFFIN produces a consensus set that satisfies non-overlapping constraints and maximizes the number of nucleosomes.

Conclusions

We report comprehensive experimental results that compares PuFFIN with recently published tools (NOrMAL, TEMPLATE FILTERING, and NucPosSimulator) on several synthetic datasets as well as real data for S. cerevisiae and P. falciparum. Experimental results show that our approach produces more accurate nucleosome maps with a higher number of non-overlapping nucleosomes than other tools.

Time resolved DNA occupancy dynamics during the respiratory oscillation uncover a global reset point in the yeast growth program

The structural dynamics of chromatin have been implicated in the regulation of fundamental eukaryotic processes, such as DNA transcription, replication and repair. Although previous studies have revealed that the chromatin landscape, nucleosome remodeling and histone modification events are intimately tied into cellular energetics and redox state, few studies undertake defined time-resolved measurements of these state variables. Here, we use metabolically synchronous, continuously-grown yeast cultures to measure DNA occupancy and track global patterns with respect to the metabolic state of the culture. Combined with transcriptome analyses and ChIP-qPCR experiments, these paint an intriguing picture where genome-wide nucleosome focusing occurs during the recovery of energy charge, followed by clearance of the promoter regions and global transcriptional slow-down, thus indicating a nucleosome-mediated “reset point” for the cycle. The reset begins at the end of the catabolic and stress-response transcriptional programs and ends prior to the start of the anabolic and cell-growth transcriptional program, and the histones on genes from both the catabolic and anabolic superclusters are deacetylated.

Binding of the transcription factor Atf1 to promoters serves as a barrier to phase nucleosome arrays and avoid cryptic transcription

Schizosaccharomyces pombe displays a large transcriptional response common to several stress conditions, regulated primarily by the transcription factor Atf1. Atf1-dependent promoters contain especially broad nucleosome depleted regions (NDRs) prior to stress imposition. We show here that basal binding of Atf1 to these promoters competes with histones to create wider NDRs at stress genes. Moreover, deletion of atf1 results in nucleosome disorganization specifically at stress coding regions and derepresses antisense transcription. Our data indicate that the transcription factor binding to promoters acts as an effective barrier to fix the +1 nucleosome and phase downstream nucleosome arrays to prevent cryptic transcription.

Genome-wide nucleosome positioning is orchestrated by genomic regions associated with DNase I hypersensitivity in rice

Nucleosome positioning dictates the DNA accessibility for regulatory proteins, and thus is critical for gene expression and regulation. It has been well documented that only a subset of nucleosomes are reproducibly positioned in eukaryotic genomes. The most prominent example of phased nucleosomes is the context of genes, where phased nucleosomes flank the transcriptional starts sites (TSSs). It is unclear, however, what factors determine nucleosome positioning in regions that are not close to genes. We mapped both nucleosome positioning and DNase I hypersensitive site (DHS) datasets across the rice genome. We discovered that DHSs located in a variety of contexts, both genic and intergenic, were flanked by strongly phased nucleosome arrays. Phased nucleosomes were also found to flank DHSs in the human genome. Our results suggest the barrier model may represent a general feature of nucleosome organization in eukaryote genomes. Specifically, regions bound with regulatory proteins, including intergenic regions, can serve as barriers that organize phased nucleosome arrays on both sides. Our results also suggest that rice DHSs often span a single, phased nucleosome, similar to the H2A.Z-containing nucleosomes observed in DHSs in the human genome.

The fundamental unit of chromatin is the nucleosome, which consists of 147 bp of DNA wrapped around a histone octamer containing four core histones (H3, H4, H2A, and H2B). Nucleosome positioning in the genome affects the DNA accessibility for regulatory proteins, and thus is critical for gene expression and regulation. Genomic regions associated with regulatory proteins are associated with a pronounced sensitivity to DNase I digestion, and are thus called DNase I hypersensitive sites (DHSs). It is well known that only a subset of nucleosomes are reproducibly positioned in eukaryotic genomes. However, it is less clear what factors determine genome-wide nucleosome positioning, especially in intergenic regions. We mapped both nucleosome positioning and DHS datasets across the rice genome. We discovered that DHSs located in a variety of contexts, both genic and intergenic, were flanked by strongly phased nucleosome arrays. We confirmed the same association of DHSs with phased nucleosomes in the human genome. We conclude that genomic loci associated with a diverse set of regulatory proteins are major determinants of nucleosome phasing, and this is true in both genic and intergenic regions.