Nucleosome-depleted Regions Research Articles

Measuring DNA mechanics on the genome scale.

Mechanical deformations of DNA such as bending are ubiquitous and implicated in diverse cellular functions1. However, the lack of high-throughput tools to directly measure the mechanical properties of DNA limits our understanding of whether and how DNA sequences modulate DNA mechanics and associated chromatin transactions genome-wide. We developed an assay called loop-seq to measure the intrinsic cyclizability of DNA – a proxy for DNA bendability – in high throughput. We measured the intrinsic cyclizabilities of 270,806 50 bp DNA fragments that span the entire length of S. cerevisiae chromosome V and other genomic regions, and also include random sequences. We discovered sequence-encoded regions of unusually low bendability upstream of Transcription Start Sites (TSSs). These regions disfavor the sharp DNA bending required for nucleosome formation and are co-centric with known Nucleosome Depleted Regions (NDRs). We show biochemically that low bendability of linker DNA located about 40 bp away from a nucleosome edge inhibits nucleosome sliding into the linker by the chromatin remodeler INO80. The observation explains how INO80 can create promoter-proximal nucleosomal arrays in the absence of any other factors2 by reading the DNA mechanical landscape. We show that chromosome wide, nucleosomes are characterized by high DNA bendability near dyads and low bendability near the linkers. This contrast increases for nucleosomes deeper into gene bodies, suggesting that DNA mechanics plays a previously unappreciated role in organizing nucleosomes far from the TSS, where nucleosome remodelers predominate. Importantly, random substitution of synonymous codons does not preserve this contrast, suggesting that the evolution of codon choice has been impacted by selective pressure to preserve sequence-encoded mechanical modulations along genes. We also provide evidence that transcription through the TSS-proximal nucleosomes is impacted by local DNA mechanics. Overall, this first genome-scale map of DNA mechanics hints at a ‘mechanical code’ with broad functional implications.

Chromatin accessibility dynamics of Chlamydia-infected epithelial cells

Chlamydia are Gram-negative, obligate intracellular bacterial pathogens responsible for a broad spectrum of human and animal diseases. In humans, Chlamydia trachomatis is the most prevalent bacterial sexually transmitted infection worldwide and is the causative agent of trachoma (infectious blindness) in disadvantaged populations. Over the course of its developmental cycle, Chlamydia extensively remodels its intracellular niche and parasitises the host cell for nutrients, with substantial resulting changes to the host cell transcriptome and proteome. However, little information is available on the impact of chlamydial infection on the host cell epigenome and global gene regulation. Regions of open eukaryotic chromatin correspond to nucleosome-depleted regions, which in turn are associated with regulatory functions and transcription factor binding. We applied formaldehyde-assisted isolation of regulatory elements enrichment followed by sequencing (FAIRE-Seq) to generate temporal chromatin maps of C. trachomatis-infected human epithelial cells in vitro over the chlamydial developmental cycle. We detected both conserved and distinct temporal changes to genome-wide chromatin accessibility associated with C. trachomatis infection. The observed differentially accessible chromatin regions include temporally-enriched sets of transcription factors, which may help shape the host cell response to infection. These regions and motifs were linked to genomic features and genes associated with immune responses, re-direction of host cell nutrients, intracellular signalling, cell–cell adhesion, extracellular matrix, metabolism and apoptosis. This work provides another perspective to the complex response to chlamydial infection, and will inform further studies of transcriptional regulation and the epigenome in Chlamydia-infected human cells and tissues.

Distinctive regulatory architectures of germline-active and somatic genes in C. elegans

RNA profiling has provided increasingly detailed knowledge of gene expression patterns, yet the different regulatory architectures that drive them are not well understood. To address this, we profiled and compared transcriptional and regulatory element activities across five tissues of Caenorhabditis elegans, covering ∼90% of cells. We find that the majority of promoters and enhancers have tissue-specific accessibility, and we discover regulatory grammars associated with ubiquitous, germline, and somatic tissue–specific gene expression patterns. In addition, we find that germline-active and soma-specific promoters have distinct features. Germline-active promoters have well-positioned +1 and −1 nucleosomes associated with a periodic 10-bp WW signal (W = A/T). Somatic tissue–specific promoters lack positioned nucleosomes and this signal, have wide nucleosome-depleted regions, and are more enriched for core promoter elements, which largely differ between tissues. We observe the 10-bp periodic WW signal at ubiquitous promoters in other animals, suggesting it is an ancient conserved signal. Our results show fundamental differences in regulatory architectures of germline and somatic tissue–specific genes, uncover regulatory rules for generating diverse gene expression patterns, and provide a tissue-specific resource for future studies.

Nucleosome Positioning Regulates the Establishment, Stability, and Inheritance of Heterochromatin in Saccharomyces cerevisiae

Heterochromatic domains are complex structures composed of nucleosome arrays that are bound by silencing factors. This composition raises the possibility that certain configurations of nucleosome arrays facilitate heterochromatic silencing. We tested this possibility in Saccharomyces cerevisiae by systematically altering the distance between heterochromatic nucleosome-depleted regions (NDRs), which is predicted to affect local nucleosome positioning by limiting how nucleosomes can be packed between NDRs. Consistent with this prediction, serial deletions that altered the distance between heterochromatic NDRs revealed a striking oscillatory relationship between inter-NDR distance and defects in nucleosome positioning. Furthermore, conditions that caused poor nucleosome positioning also led to defects in both heterochromatin stability and the ability of cells to generate and inherit epigenetic transcriptional states. These findings strongly suggest that nucleosome positioning can contribute to formation and maintenance of functional heterochromatin and point to previously unappreciated roles of NDR positioning within heterochromatic domains.

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

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

Chromatin remodeler Ino80C acts independently of H2A.Z to evict promoter nucleosomes and stimulate transcription of highly expressed genes in yeast

The chromatin remodelers SWI/SNF and RSC function in evicting promoter nucleosomes at highly expressed yeast genes, particularly those activated by transcription factor Gcn4. Ino80 remodeling complex (Ino80C) can establish nucleosome-depleted regions (NDRs) in reconstituted chromatin, and was implicated in removing histone variant H2A.Z from the -1 and +1 nucleosomes flanking NDRs; however, Ino80C's function in transcriptional activation in vivo is not well understood. Analyzing the cohort of Gcn4-induced genes in ino80Δ mutants has uncovered a role for Ino80C on par with SWI/SNF in evicting promoter nucleosomes and transcriptional activation. Compared to SWI/SNF, Ino80C generally functions over a wider region, spanning the -1 and +1 nucleosomes, NDR and proximal genic nucleosomes, at genes highly dependent on its function. Defects in nucleosome eviction in ino80Δ cells are frequently accompanied by reduced promoter occupancies of TBP, and diminished transcription; and Ino80 is enriched at genes requiring its remodeler activity. Importantly, nuclear depletion of Ino80 impairs promoter nucleosome eviction even in a mutant lacking H2A.Z. Thus, Ino80C acts widely in the yeast genome together with RSC and SWI/SNF in evicting promoter nucleosomes and enhancing transcription, all in a manner at least partly independent of H2A.Z editing.

Genome wide distribution of G-quadruplexes and their impact on gene expression in malaria parasites.

Mechanisms of transcriptional control in malaria parasites are still not fully understood. The positioning patterns of G-quadruplex (G4) DNA motifs in the parasite’s AT-rich genome, especially within the var gene family which encodes virulence factors, and in the vicinity of recombination hotspots, points towards a possible regulatory role of G4 in gene expression and genome stability. Here, we carried out the most comprehensive genome-wide survey, to date, of G4s in the Plasmodium falciparum genome using G4Hunter, which identifies G4 forming sequences (G4FS) considering their G-richness and G-skewness. We show an enrichment of G4FS in nucleosome-depleted regions and in the first exon of var genes, a pattern that is conserved within the closely related Laverania Plasmodium parasites. Under G4-stabilizing conditions, i.e., following treatment with pyridostatin (a high affinity G4 ligand), we show that a bona fide G4 found in the non-coding strand of var promoters modulates reporter gene expression. Furthermore, transcriptional profiling of pyridostatin-treated parasites, shows large scale perturbations, with deregulation affecting for instance the ApiAP2 family of transcription factors and genes involved in ribosome biogenesis. Overall, our study highlights G4s as important DNA secondary structures with a role in Plasmodium gene expression regulation, sub-telomeric recombination and var gene biology.

Fine Chromatin-Driven Mechanism of Transcription Interference by Antisense Noncoding Transcription

Eukaryotic genomes are almost entirely transcribed by RNA polymerase II. Consequently, the transcription of long noncoding RNAs often overlaps with coding gene promoters, triggering potential gene repression through a poorly characterized mechanism of transcription interference. Here, we propose a comprehensive model of chromatin-based transcription interference in Saccharomyces cerevisiae (S.cerevisiae). By using a noncoding transcription-inducible strain, we analyze the relationship between antisense elongation and coding sense repression, nucleosome occupancy, and transcription-associated histone modifications using near-base pair resolution techniques. We show that antisense noncoding transcription leads to the deacetylation of a subpopulation of -1/+1 nucleosomes associated with increased H3K36me3. Reduced acetylation results in the decreased binding of the RSC chromatin remodeler at -1/+1 nucleosomes and subsequent sliding into the nucleosome-depleted region hindering pre-initiation complex association. Finally, we extend our model by showing that natural antisense noncoding transcription significantly represses ∼20% of S.cerevisiae genes through this chromatin-based transcription interference mechanism.

Cryo-EM structure of SWI/SNF complex bound to a nucleosome.

The chromatin remodeling complex SWI/SNF is highly conserved and plays critical roles in various cellular processes including transcription and DNA damage repair1,2. It hydrolyzes ATP to remodel chromatin structure by sliding and evicting histone octamers3-8, creating DNA regions that become accessible to other essential factors. However, our mechanistic understanding of the remodeling activity is largely hindered by the lack of a high-resolution structure of any complex from this family. Here we report the first structure of SWI/SNF from the yeast S. cerevisiae bound to a nucleosome at near atomic resolution determined by cryo-electron microscopy (cryo-EM). In the structure, the Arp module is sandwiched between the ATPase and the rest of the complex, with the Snf2 HSA domain connecting all modules. The body contains an assembly scaffold composed of conserved subunits Snf12 (SMARCD/BAF60), Snf5 (SMARCB1/BAF47/ INI1) and an asymmetric dimer of Swi3 (SMARCC/BAF155/170). Another conserved subunit Swi1 (ARID1/BAF250) resides in the core of SWI/SNF, acting as a molecular hub. We also observed interactions between Snf5 and the histones at the acidic patch, which could serve as an anchor during active DNA translocation. Our structure allows us to map and rationalize a subset of cancer-related mutations in the human SWI/SNF complex and propose a model of how SWI/SNF recognizes and remodels the +1 nucleosome to generate nucleosome-depleted regions during gene activation9.

The histone variant H2A.Z in yeast is almost exclusively incorporated into the +1 nucleosome in the direction of transcription.

Transcription start sites (TSS) in eukaryotes are characterized by a nucleosome-depleted region (NDR), which appears to be flanked upstream and downstream by strongly positioned nucleosomes incorporating the histone variant H2A.Z. H2A.Z associates with both active and repressed TSS and is important for priming genes for rapid transcriptional activation. However, the determinants of H2A.Z occupancy at specific nucleosomes and its relationship to transcription initiation remain unclear. To further elucidate the specificity of H2A.Z, we determined its genomic localization at single nucleosome resolution, as well as the localization of its chromatin remodelers Swr1 and Ino80. By analyzing H2A.Z occupancy in conjunction with RNA expression data that captures promoter-derived antisense initiation, we find that H2A.Z’s bimodal incorporation on either side of the NDR is not a general feature of TSS, but is specifically a marker for bidirectional transcription, such that the upstream flanking −1 H2A.Z-containing nucleosome is more appropriately considered as a +1 H2A.Z nucleosome for antisense transcription. The localization of H2A.Z almost exclusively at the +1 nucleosome suggests that a transcription-initiation dependent process could contribute to its specific incorporation.

CTCF-dependent chromatin boundaries formed by asymmetric nucleosome arrays with decreased linker length.

The CCCTC-binding factor (CTCF) organises the genome in 3D through DNA loops and in 1D by setting boundaries isolating different chromatin states, but these processes are not well understood. Here we investigate chromatin boundaries in mouse embryonic stem cells, defined by the regions with decreased Nucleosome Repeat Length (NRL) for ∼20 nucleosomes near CTCF sites, affecting up to 10% of the genome. We found that the nucleosome-depleted region (NDR) near CTCF is asymmetrically located >40 nucleotides 5′-upstream from the centre of CTCF motif. The strength of CTCF binding to DNA and the presence of cohesin is correlated with the decrease of NRL near CTCF, and anti-correlated with the level of asymmetry of the nucleosome array. Individual chromatin remodellers have different contributions, with Snf2h having the strongest effect on the NRL decrease near CTCF and Chd4 playing a major role in the symmetry breaking. Upon differentiation, a subset of preserved, common CTCF sites maintains asymmetric nucleosome pattern and small NRL. The sites which lost CTCF upon differentiation are characterized by nucleosome rearrangement 3′-downstream, with unchanged NDR 5′-upstream of CTCF motifs. Boundaries of topologically associated chromatin domains frequently contain several inward-oriented CTCF motifs whose effects, described above, add up synergistically.

Genome-Wide Analysis of the Nucleosome Landscape in Individuals with Coffin-Siris Syndrome

The switch/sucrose non-fermenting (SWI/SNF) complex is an ATP-dependent chromatin remodeller that regulates the spacing of nucleosomes and thereby controls gene expression. Heterozygous mutations in genes encoding subunits of the SWI/SNF complex have been reported in individuals with Coffin-Siris syndrome (CSS), with the majority of the mutations in ARID1B. CSS is a rare congenital disorder characterized by facial dysmorphisms, digital anomalies, and variable intellectual disability. We hypothesized that mutations in genes encoding subunits of the ubiquitously expressed SWI/SNF complex may lead to alterations of the nucleosome profiles in different cell types. We performed the first study on CSS-patient samples and investigated the nucleosome landscapes of cell-free DNA (cfDNA) isolated from blood plasma by whole-genome sequencing. In addition, we studied the nucleosome landscapes of CD14⁺ monocytes from CSS-affected individuals by nucleosome occupancy and methylome-sequencing (NOMe-seq) as well as their expression profiles. In cfDNA of CSS-affected individuals with heterozygous ARID1B mutations, we did not observe major changes in the nucleosome profile around transcription start sites. In CD14⁺ monocytes, we found few genomic regions with different nucleosome occupancy when compared to controls. RNA-seq analysis of CD14⁺ monocytes of these individuals detected only few differentially expressed genes, which were not in proximity to any of the identified differential nucleosome-depleted regions. In conclusion, we show that heterozygous mutations in the human SWI/SNF subunit ARID1B do not have a major impact on the nucleosome landscape or gene expression in blood cells. This might be due to functional redundancy, cell-type specificity, or alternative functions of ARID1B.

Differential nucleosome spacing in neurons and glia

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

Regulation of Gene Expression and Replication Initiation by Non-Coding Transcription: A Model Based on Reshaping Nucleosome-Depleted Regions: Influence of Pervasive Transcription on Chromatin Structure.

RNA polymerase II (RNAP II) non-coding transcription is now known to cover almost the entire eukaryotic genome, a phenomenon referred to as pervasive transcription. As a consequence, regions previously thought to be non-transcribed are subject to the passage of RNAP II and its associated proteins for histone modification. This is the case for the nucleosome-depleted regions (NDRs), which provide key sites of entry into the chromatin for proteins required for the initiation of coding gene transcription and DNA replication. In this review, recent data on the effects of pervasive transcription through NDRs are summarized and a model is proposed to explain how RNAP II-driven transcription is able to modify the nucleosomes flanking the NDRs, leading to nucleosome repositioning and NDR closure. Even though much of the mechanistic detail underlying these events remains to be elucidated, such a model provides a basis to explain how non-coding transcription through NDRs can regulate the initiation of coding gene expression and DNA replication.

Pervasive Transcription Represses Coding Gene Expression by Closing Their Nucleosome-Depleted Regions.

Pervasive Transcription Represses Coding Gene Expression by Closing Their Nucleosome-Depleted Regions.

A high-resolution 3D epigenomic map reveals insights into the creation of the prostate cancer transcriptome

To better understand the impact of chromatin structure on regulation of the prostate cancer transcriptome, we develop high-resolution chromatin interaction maps in normal and prostate cancer cells using in situ Hi-C. By combining the in situ Hi-C data with active and repressive histone marks, CTCF binding sites, nucleosome-depleted regions, and transcriptome profiling, we identify topologically associating domains (TADs) that change in size and epigenetic states between normal and prostate cancer cells. Moreover, we identify normal and prostate cancer-specific enhancer-promoter loops and involved transcription factors. For example, we show that FOXA1 is enriched in prostate cancer-specific enhancer-promoter loop anchors. We also find that the chromatin structure surrounding the androgen receptor (AR) locus is altered in the prostate cancer cells with many cancer-specific enhancer-promoter loops. This creation of 3D epigenomic maps enables a better understanding of prostate cancer biology and mechanisms of gene regulation.

Arabidopsis SWR1-associated protein methyl-CpG-binding domain 9 is required for histone H2A.Z deposition

Deposition of the histone variant H2A.Z by the SWI2/SNF2-Related 1 chromatin remodeling complex (SWR1-C) is important for gene regulation in eukaryotes, but the composition of the Arabidopsis SWR1-C has not been thoroughly characterized. Here, we aim to identify interacting partners of a conserved Arabidopsis SWR1 subunit ACTIN-RELATED PROTEIN 6 (ARP6). We isolate nine predicted components and identify additional interactors implicated in histone acetylation and chromatin biology. One of the interacting partners, methyl-CpG-binding domain 9 (MBD9), also strongly interacts with the Imitation SWItch (ISWI) chromatin remodeling complex. MBD9 is required for deposition of H2A.Z at a distinct subset of ARP6-dependent loci. MBD9 is preferentially bound to nucleosome-depleted regions at the 5’ ends of genes containing high levels of activating histone marks. These data suggest that MBD9 is a SWR1-C interacting protein required for H2A.Z deposition at a subset of actively transcribing genes.

NF-Y controls fidelity of transcription initiation at gene promoters through maintenance of the nucleosome-depleted region

Faithful transcription initiation is critical for accurate gene expression, yet the mechanisms underlying specific transcription start site (TSS) selection in mammals remain unclear. Here, we show that the histone-fold domain protein NF-Y, a ubiquitously expressed transcription factor, controls the fidelity of transcription initiation at gene promoters in mouse embryonic stem cells. We report that NF-Y maintains the region upstream of TSSs in a nucleosome-depleted state while simultaneously protecting this accessible region against aberrant and/or ectopic transcription initiation. We find that loss of NF-Y binding in mammalian cells disrupts the promoter chromatin landscape, leading to nucleosomal encroachment over the canonical TSS. Importantly, this chromatin rearrangement is accompanied by upstream relocation of the transcription pre-initiation complex and ectopic transcription initiation. Further, this phenomenon generates aberrant extended transcripts that undergo translation, disrupting gene expression profiles. These results suggest NF-Y is a central player in TSS selection in metazoans and highlight the deleterious consequences of inaccurate transcription initiation.

The NSL complex-mediated nucleosome landscape is required to maintain transcription fidelity and suppression of transcription noise.

Nucleosomal organization at gene promoters is critical for transcription, with a nucleosome-depleted region (NDR) at transcription start sites (TSSs) being required for transcription initiation. How NDRs and the precise positioning of the +1 nucleosomes are maintained on active genes remains unclear. Here, we report that the Drosophila nonspecific lethal (NSL) complex is necessary to maintain this stereotypical nucleosomal organization at promoters. Upon NSL1 depletion, nucleosomes invade the NDRs at TSSs of NSL-bound genes. NSL complex member NSL3 binds to TATA-less promoters in a sequence-dependent manner. The NSL complex interacts with the NURF chromatin remodeling complex and is necessary and sufficient to recruit NURF to target promoters. Not only is the NSL complex essential for transcription, but it is required for accurate TSS selection for genes with multiple TSSs. Furthermore, loss of the NSL complex leads to an increase in transcriptional noise. Thus, the NSL complex establishes a canonical nucleosomal organization that enables transcription and determines TSS fidelity.

Contrasting roles of the RSC and ISW1/CHD1 chromatin remodelers in RNA polymerase II elongation and termination.

Most yeast genes have a nucleosome-depleted region (NDR) at the promoter and an array of regularly spaced nucleosomes phased relative to the transcription start site. We have examined the interplay between RSC (a conserved essential SWI/SNF-type complex that determines NDR size) and the ISW1, CHD1, and ISW2 nucleosome spacing enzymes in chromatin organization and transcription, using isogenic strains lacking all combinations of these enzymes. The contributions of these remodelers to chromatin organization are largely combinatorial, distinct, and nonredundant, supporting a model in which the +1 nucleosome is positioned by RSC and then used as a reference nucleosome by the spacing enzymes. Defective chromatin organization correlates with altered RNA polymerase II (Pol II) distribution. RSC-depleted cells exhibit low levels of elongating Pol II and high levels of terminating Pol II, consistent with defects in both termination and initiation, suggesting that RSC facilitates both. Cells lacking both ISW1 and CHD1 show the opposite Pol II distribution, suggesting elongation and termination defects. These cells have extremely disrupted chromatin, with high levels of closely packed dinucleosomes involving the second (+2) nucleosome. We propose that ISW1 and CHD1 facilitate Pol II elongation by separating closely packed nucleosomes.