Nucleosome Positioning Research Articles

Bi-directional nucleosome sliding by the Chd1 chromatin remodeler integrates intrinsic sequence-dependent and ATP-dependent nucleosome positioning.

Chromatin remodelers use a helicase-type ATPase motor to shift DNA around the histone core. Although not directly reading out the DNA sequence, some chromatin remodelers exhibit a sequence-dependent bias in nucleosome positioning, which presumably reflects properties of the DNA duplex. Here, we show how nucleosome positioning by the Chd1 remodeler is influenced by local DNA perturbations throughout the nucleosome footprint. Using site-specific DNA cleavage coupled with next-generation sequencing, we show that nucleosomes shifted by Chd1 can preferentially localize DNA perturbations - poly(dA:dT) tracts, DNA mismatches, and single-nucleotide insertions - about a helical turn outside the Chd1 motor domain binding site, super helix location 2 (SHL2). This phenomenon occurs with both the Widom 601 positioning sequence and the natural+1 nucleosome sequence from the Saccharomycescerevisiae SWH1 gene. Our modeling indicates that localization of DNA perturbations about a helical turn outward from SHL2 results from back-and-forth sliding due to remodeler action on both sides of the nucleosome. Our results also show that barrier effects from DNA perturbations can be extended by the strong phasing of nucleosome positioning sequences.

Proteomics contributions to epigenetic drug discovery.

The combined activity of epigenetic features, which include histone post-translational modifications, DNA methylation, and nucleosome positioning, regulates gene expression independently from changes in the DNA sequence, defining how the shared genetic information of an organism is used to generate different cell phenotypes. Alterations in epigenetic processes have been linked with a multitude of diseases, including cancer, fueling interest in the discovery of drugs targeting the proteins responsible for writing, erasing, or reading histone and DNA modifications. Mass spectrometry (MS)-based proteomics has emerged as a versatile tool that can assist drug discovery pipelines from target validation, through target deconvolution, to monitoring drug efficacy in vivo. Here, we provide an overview of the contributions of MS-based proteomics to epigenetic drug discovery, describing the main approaches that can be used to support different drug discovery pipelines and highlighting how they contributed to the development and characterization of epigenetic drugs.

Short Sequence Aligner Benchmarking for Chromatin Research.

Much of today's molecular science revolves around next-generation sequencing. Frequently, the first step in analyzing such data is aligning sequencing reads to a reference genome. This step is often taken for granted, but any analysis downstream of the alignment will be affected by the aligner's ability to correctly map sequences. In most cases, for research into chromatin structure and nucleosome positioning, ATAC-seq, ChIP-seq, and MNase-seq experiments use short read lengths. How well aligners manage these reads is critical. Most aligner programs will output mapped reads and unmapped reads. However, from a biological point of view, reads will fall into one of three categories: correctly mapped, incorrectly mapped, and unmapped. While increased sequencing depth can often compensate for unmapped reads, incorrectly and correctly mapped reads appear algorithmically identical but can produce biologically significant alterations in the results. For this reason, we are benchmarking various alignment programs to determine their propensity to incorrectly map short reads. As short-read alignment is an important step in ATAC-seq, ChIP-seq, and MNase-seq experiments, caution should be taken in mapping reads to ensure that the most accurate conclusions can be made from the data generated. Our analysis is intended to help investigators new to the field pick the alignment program best suited for their experimental conditions. In general, the aligners we tested performed well. BWA, Bowtie2, and Chromap were all exceptionally accurate, and we recommend using them. Furthermore, we show that longer read lengths do in fact lead to more accurate mappings.

Nucleosome repositioning in chronic lymphocytic leukemia.

The location of nucleosomes in the human genome determines the primary chromatin structure and regulates access to regulatory regions. However, genome-wide information on deregulated nucleosome occupancy and its implications in primary cancer cells is scarce. Here, we conducted a genome-wide comparison of high-resolution nucleosome maps in peripheral blood B cells from patients with chronic lymphocytic leukemia (CLL) and healthy individuals at single-base-pair resolution. Our investigation uncovered significant changes of nucleosome positioning in CLL. Globally, the spacing between nucleosomes-the nucleosome repeat length (NRL)-is shortened in CLL. This effect is stronger in the more aggressive IGHV-unmutated CLL subtype than in the IGHV-mutated CLL subtype. Changes in nucleosome occupancy at specific sites are linked to active chromatin remodeling and reduced DNA methylation. Nucleosomes lost or gained in CLL marks differential binding of 3D chromatin organizers such as CTCF as well as immune response-related transcription factors and delineated mechanisms of epigenetic deregulation. The principal component analysis of nucleosome occupancy in cancer-specific regions allowed the classification of samples between cancer subtypes and normal controls. Furthermore, patients could be better assigned to CLL subtypes according to differential nucleosome occupancy than based on DNA methylation or gene expression. Thus, nucleosome positioning constitutes a novel readout to dissect molecular mechanisms of disease progression and to stratify patients. Furthermore, we anticipate that the global nucleosome repositioning detected in our study, such as changes in the NRL, can be exploited for liquid biopsy applications based on cell-free DNA to stratify patients and monitor disease progression.

Virus hijacking of host epigenetic machinery to impair immune response.

Epigenetic modifications, such as DNA hypermethylation, histone acetylation/methylation, or nucleosome positioning, result in differential gene expression. These modifications can have an impact on various pathways, including host antiviral immune responses. In this review, we summarize the current understanding of epigenetic modifications induced by viruses to counteract host antiviral immune responses, which are crucial for establishing and maintaining infection of viruses. Finally, we provide insights into the potential use of epigenetic modulators in combating viral infections and virus-induced diseases.

Analysis of the complex between MBD2 and the histone deacetylase core of NuRD reveals key interactions critical for gene silencing.

The nucleosome remodeling and deacetylase (NuRD) complex modifies nucleosome positioning and chromatin compaction to regulate gene expression. The methyl-CpG-binding domain proteins 2 and 3 (MBD2 and MBD3) play a critical role in complex formation; however, the molecular details of how they interact with other NuRD components have yet to be fully elucidated. We previously showed that an intrinsically disordered region (IDR) of MBD2 is necessary and sufficient to bind to the histone deacetylase core of NuRD. Building on that work, we have measured the inherent structural propensity of the MBD2-IDR using solvent and site-specific paramagnetic relaxation enhancement measurements. We then used the AlphaFold2 machine learning software to generate a model of the complex between MBD2 and the histone deacetylase core of NuRD. This model is remarkably consistent with our previous studies, including the current paramagnetic relaxation enhancement data. The latter suggests that the free MBD2-IDR samples conformations similar to the bound structure. We tested this model of the complex extensively by mutating key contact residues and measuring binding using an intracellular bioluminescent resonance energy transfer assay. Furthermore, we identified protein contacts that, when mutated, disrupted gene silencing by NuRD in a cell model of fetal hemoglobin regulation. Hence, this work provides insights into the formation of NuRD and highlights critical binding pockets that may be targeted to block gene silencing for therapy. Importantly, we show that AlphaFold2 can generate a credible model of a large complex that involves an IDR that folds upon binding.

Centromere plasticity with evolutionary conservation and divergence uncovered by wheat 10+ genomes.

Centromeres are the chromosomal regions that play a crucial role in maintaining genomic stability. The underling highly repetitive DNA sequences can evolve quickly in most eukaryotes, and promote karyotype evolution. Despite their variability, it is not fully understood how these widely variable sequences ensure the homeostasis of centromere function. In this study, we investigated the genetics and epigenetics of centromeres in a population of wheat lines from global breeding programs. We captured a high degree of sequences, positioning, and epigenetic variations in the large and complex wheat centromeres. We found that the most CENH3-associated repeats are Cereba element of retrotransposons and exhibit phylogenetic homogenization across different wheat lines, but the less-associated repeat sequences diverge on their own way in each wheat line, implying specific mechanisms for selecting certain repeat types as functional core centromeres. Furthermore, we observed that CENH3 nucleosome structures display looser wrapping of DNA termini on complex centromeric repeats, including the repositioned centromeres. We also found that strict CENH3 nucleosome positioning and intrinsic DNA features play a role in determining centromere identity among different lines. Specific non-B form DNAs were substantially associated with CENH3 nucleosomes for the repositioned centromeres. These findings suggest that multiple mechanisms were involved in the adaptation of CENH3 nucleosomes that can stabilize centromeres. Ultimately, we proposed a remarkable epigenetic plasticity of centromere chromatin within the diverse genomic context, and the high robustness is crucial for maintaining centromere function and genome stability in wheat 10+ lines as a result of past breeding selections.

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.

Screening and identification of potential biomarkers for pancreatic cancer: An integrated bioinformatics analysis

Pancreatic cancer is one of the highly invasive and the seventh most common cause of death among cancers worldwide. To identify essential genes and the involved mechanisms in pancreatic cancer, we used bioinformatics analysis to identify potential biomarkers for pancreatic cancer management. Gene expression profiles of pancreatic cancer patients and normal tissues were screened and downloaded from The Cancer Genome Atlas (TCGA) bioinformatics database. The Differentially expressed genes (DEGs) were identified among gene expression signatures of normal and pancreatic cancer, using R software. Then, enrichment analysis of the DEGs, including Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, was performed by an interactive and collaborative HTML5 gene list enrichment analysis tool (enrichr) and ToppGene. The protein-protein interaction (PPI) network was also constructed using the Search Tool for the Retrieval of Interacting Genes (STRING) database and ToppGenet web based tool followed by identifying hub genes of the top 100 DEGs in pancreatic cancer using Cytoscape software. Over 2000 DEGs with variable log2 fold (LFC) were identified among 34,706 genes. Principal component analysis showed that the top 20 DEGs, including H1–4, H1–5, H4C3, H4C2, RN7SL2, RN7SL3, RN7SL4P, RN7SKP80, SCARNA12, SCARNA10, SCARNA5, SCARNA7, SCARNA6, SCARNA21, SCARNA9, SCARNA13, SNORA73B, SNORA53, SNORA54 might distinguish pancreatic cancer from normal tissue. GO analysis showed that the top DEGs have more enriched in the negative regulation of gene silencing, negative regulation of chromatin organization, negative regulation of chromatin silencing, nucleosome positioning, regulation of chromatin silencing, and nucleosomal DNA binding. KEGG analysis identified an association between pancreatic cancer and systemic lupus erythematosus, alcoholism, neutrophil extracellular trap formation, and viral carcinogenesis. In PPI network analysis, we found that the different types of histone-encoding genes are involved as hub genes in the carcinogenesis of pancreatic cancer. In conclusion, our bioinformatics analysis identified genes that were significantly related to the prognosis of pancreatic cancer patients. These genes and pathways could serve as new potential prognostic markers and be used to develop treatments for pancreatic cancer patients.

Structure and Dynamics of Compact Dinucleosomes: Analysis by Electron Microscopy and spFRET.

Formation of compact dinucleosomes (CODIs) occurs after collision between adjacent nucleosomes at active regulatory DNA regions. Although CODIs are likely dynamic structures, their structural heterogeneity and dynamics were not systematically addressed. Here, single-particle Förster resonance energy transfer (spFRET) and electron microscopy were employed to study the structure and dynamics of CODIs. spFRET microscopy in solution and in gel revealed considerable uncoiling of nucleosomal DNA from the histone octamer in a fraction of CODIs, suggesting that at least one of the nucleosomes is destabilized in the presence of the adjacent closely positioned nucleosome. Accordingly, electron microscopy analysis suggests that up to 30 bp of nucleosomal DNA are involved in transient uncoiling/recoiling on the octamer. The more open and dynamic nucleosome structure in CODIs cannot be stabilized by histone chaperone Spt6. The data suggest that proper internucleosomal spacing is an important determinant of chromatin stability and support the possibility that CODIs could be intermediates of chromatin disruption.

Nucleosome retention by histone chaperones and remodelers occludes pervasive DNA-protein binding.

DNA packaging within chromatin depends on histone chaperones and remodelers that form and position nucleosomes. Cells express multiple such chromatin regulators with overlapping in-vitro activities. Defining specific in-vivo activities requires monitoring histone dynamics during regulator depletion, which has been technically challenging. We have recently generated histone-exchange sensors in Saccharomyces cerevisiae, which we now use to define the contributions of 15 regulators to histone dynamics genome-wide. While replication-independent exchange in unperturbed cells maps to promoters, regulator depletions primarily affected gene bodies. Depletion of Spt6, Spt16 or Chd1 sharply increased nucleosome replacement sequentially at the beginning, middle or end of highly expressed gene bodies. They further triggered re-localization of chaperones to affected gene body regions, which compensated for nucleosome loss during transcription complex passage, but concurred with extensive TF binding in gene bodies. We provide a unified quantitative screen highlighting regulator roles in retaining nucleosome binding during transcription and preserving genomic packaging.

The regulation of chromatin configuration at AGAMOUS locus by LFR-SYD-containing complex is critical for reproductive organ development in Arabidopsis.

Switch defective/sucrose non-fermentable (SWI/SNF) chromatin remodeling complexes are evolutionarily conserved, multi-subunit machinery that play vital roles in the regulation of gene expression by controlling nucleosome positioning and occupancy. However, little is known about the subunit composition of SPLAYED (SYD)-containing SWI/SNF complexes in plants. Here, we show that the Arabidopsis thaliana Leaf and Flower Related (LFR) is a subunit of SYD-containing SWI/SNF complexes. LFR interacts directly with multiple SWI/SNF subunits, including the catalytic ATPase subunit SYD, in vitro and in vivo. Phenotypic analyses of lfr-2 mutant flowers revealed that LFR is important for proper filament and pistil development, resembling the function of SYD. Transcriptome profiling revealed that LFR and SYD shared a subset of co-regulated genes. We further demonstrate that the LFR and SYD interdependently activate the transcription of AGAMOUS (AG), a C-class floral organ identity gene, by regulating the occupation of nucleosome, chromatin loop, histone modification, and Pol II enrichment on the AG locus. Furthermore, the chromosome conformation capture (3C) assay revealed that the gene loop at AG locus is negatively correlated with the AG expression level, and LFR-SYD was functional to demolish the AG chromatin loop to promote its transcription. Collectively, these results provide insight into the molecular mechanism of the Arabidopsis SYD-SWI/SNF complex in the control of higher chromatin conformation of the floral identity gene essential to plant reproductive organ development.

Epigenetic mechanisms to propagate histone acetylation by p300/CBP

Histone acetylation is important for the activation of gene transcription but little is known about its direct read/write mechanisms. Here, we report cryogenic electron microscopy structures in which a p300/CREB-binding protein (CBP) multidomain monomer recognizes histone H4 N-terminal tail (NT) acetylation (ac) in a nucleosome and acetylates non-H4 histone NTs within the same nucleosome. p300/CBP not only recognized H4NTac via the bromodomain pocket responsible for reading, but also interacted with the DNA minor grooves via the outside of that pocket. This directed the catalytic center of p300/CBP to one of the non-H4 histone NTs. The primary target that p300 writes by reading H4NTac was H2BNT, and H2BNTac promoted H2A-H2B dissociation from the nucleosome. We propose a model in which p300/CBP replicates histone N-terminal tail acetylation within the H3-H4 tetramer to inherit epigenetic storage, and transcribes it from the H3-H4 tetramer to the H2B-H2A dimers to activate context-dependent gene transcription through local nucleosome destabilization.

The role of transcript regions and amino acid choice in nucleosome positioning.

Eukaryotic DNA is organized and compacted in a string of nucleosomes, DNA-wrapped protein cylinders. The positions of nucleosomes along DNA are not random but show well-known base pair sequence preferences that result from the sequence-dependent elastic and geometric properties of the DNA double helix. Here, we focus on DNA around transcription start sites, which are known to typically attract nucleosomes in multicellular life forms through their high GC content. We aim to understand how these GC signals, as observed in genome-wide averages, are produced and encoded through different genomic regions (mainly 5' UTRs, coding exons, and introns). Our study uses a bioinformatics approach to decompose the genome-wide GC signal into between-region and within-region signals. We find large differences in GC signal contributions between vertebrates and plants and, remarkably, even between closely related species. Introns contribute most to the GC signal in vertebrates, while in plants the exons dominate. Further, we find signal strengths stronger on DNA than on mRNA, suggesting a biological function of GC signals along the DNA itself, as is the case for nucleosome positioning. Finally, we make the surprising discovery that both the choice of synonymous codons and amino acids contribute to the nucleosome positioning signal.

Cooperation between bHLH transcription factors and histones for DNA access

The basic helix–loop–helix (bHLH) family of transcription factors recognizes DNA motifs known as E-boxes (CANNTG) and includes 108 members1. Here we investigate how chromatinized E-boxes are engaged by two structurally diverse bHLH proteins: the proto-oncogene MYC-MAX and the circadian transcription factor CLOCK-BMAL1 (refs. 2,3). Both transcription factors bind to E-boxes preferentially near the nucleosomal entry–exit sites. Structural studies with engineered or native nucleosome sequences show that MYC-MAX or CLOCK-BMAL1 triggers the release of DNA from histones to gain access. Atop the H2A–H2B acidic patch4, the CLOCK-BMAL1 Per-Arnt-Sim (PAS) dimerization domains engage the histone octamer disc. Binding of tandem E-boxes5–7 at endogenous DNA sequences occurs through direct interactions between two CLOCK-BMAL1 protomers and histones and is important for circadian cycling. At internal E-boxes, the MYC-MAX leucine zipper can also interact with histones H2B and H3, and its binding is indirectly enhanced by OCT4 elsewhere on the nucleosome. The nucleosomal E-box position and the type of bHLH dimerization domain jointly determine the histone contact, the affinity and the degree of competition and cooperativity with other nucleosome-bound factors.

Molecular Mechanism of Nucleosome Recognition by the Pioneer Transcription Factor Sox.

Pioneer transcription factors (PTFs) have the remarkable ability to directly bind to chromatin to stimulate vital cellular processes. In this work, we dissect the universal binding mode of Sox PTF by combining extensive molecular simulations and physiochemistry approaches, along with DNA footprinting techniques. As a result, we show that when Sox consensus DNA is located at the solvent-facing DNA strand, Sox binds to the compact nucleosome without imposing any significant conformational changes. We also reveal that the base-specific Sox:DNA interactions (base reading) and Sox-induced DNA changes (shape reading) are concurrently required for sequence-specific nucleosomal DNA recognition. Among three different nucleosome positions located on the positive DNA arm, a sequence-specific reading mechanism is solely satisfied at the superhelical location 2 (SHL2). While SHL2 acts transparently for solvent-facing Sox binding, among the other two positions, SHL4 permits only shape reading. The final position, SHL0 (dyad), on the other hand, allows no reading mechanism. These findings demonstrate that Sox-based nucleosome recognition is essentially guided by intrinsic nucleosome properties, permitting varying degrees of DNA recognition.

High-yield ligation-free assembly of DNA constructs with nucleosome positioning sequence repeats for single-molecule manipulation assays

Force and torque spectroscopy have provided unprecedented insights into the mechanical properties, conformational transitions, and dynamics of DNA and DNA-protein complexes, notably nucleosomes. Reliable single-molecule manipulation measurements require, however, specific and stable attachment chemistries to tether the molecules of interest. Here, we present a functionalization strategy for DNA that enables high-yield production of constructs for torsionally constrained and very stable attachment. The method is based on two subsequent PCRs: first ∼380bp long DNA strands are generated that contain multiple labels, which are used as "megaprimers" in a second PCR to generate ∼kbp long double-stranded DNA constructs with multiple labels at the respective ends. To achieve high-force stability, we use dibenzocyclooctyne-based click chemistry for covalent attachment to the surface and biotin-streptavidin coupling to the bead. The resulting tethers are torsionally constrained and extremely stable under load, with an average lifetime of 70± 3h at 45 pN. The high yield of the approach enables nucleosome reconstitution by salt dialysis on the functionalized DNA, and we demonstrate proof-of-concept measurements on nucleosome assembly statistics and inner turn unwrapping under force. We anticipate that our approach will facilitate a range of studies of DNA interactions and nucleoprotein complexes under forces and torques.

Using Synthetic DNA Libraries to Investigate Chromatin and Gene Regulation.

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

Grand challenge in chromatin epigenomics: everything, everywhere, all at once

Our understanding of the regulation and functions of histone modifications has come a long way since they were first reported in the mid-1960s. So too has our understanding of the importance of DNA methylation, histone variants, nucleosome locations and arrangements, and progressively higher order structures that impact when and where DNA-templated processes take place. Recent advances have even allowed the first ever complete sequencing and epigenomic profiles of individual chromosomes from telomere to telomere, including highly repetitive regions that were previously refractory to analysis. The regulatory power of chromatin organization for gene transcription, DNA replication, recombination and repair is undisputable. Still, an ongoing challenge is to understand the full spectrum of changes (everything) that impact processes in cells and tissues (everywhere) and how each change impacts others (all at once).

Development of multifunctional synthetic nucleosomes to interrogate chromatin-mediated protein interactions

Various proteins bind to chromatin to regulate DNA and its associated processes such as replication, transcription, and damage repair. The identification and characterization of these chromatin-associating proteins remain a challenge, as their interactions with chromatin often occur within the context of the local nucleosome or chromatin structure, which makes conventional peptide-based strategies unsuitable. Here, we developed a simple and robust protein labeling chemistry to prepare synthetic multifunctional nucleosomes that carry a photoreactive group, a biorthogonal handle, and a disulfide moiety to examine chromatin-protein interactions in a nucleosomal context. Using the prepared protein- and nucleosome-based photoaffinity probes, we examined a number of protein-protein and protein-nucleosome interactions. In particular, we (i) mapped the binding sites for the HMGN2-nucleosome interaction, (ii) provided the evidence for transition between the active and poised states of DOT1L in recognizing H3K79 within the nucleosome, and (iii) identified OARD1 and LAP2α as nucleosome acidic patch-associating proteins. This study provides powerful and versatile chemical tools for interrogating chromatin-associating proteins.