Transcriptional activation of yeast genes disrupts intragenic nucleosome phasing

Histone N-tails modulate sequence-specific positioning of nucleosomes

The precise mechanisms governing sequence-dependent positioning of nucleosomes on DNA remain unknown in detail. Existing algorithms, taking into account the sequence-dependent deformability of DNA and its interactions with the histone globular domains, predict rotational setting of only 65% of human nucleosomes mapped in vivo. To uncover novel factors responsible for the nucleosome positioning, we analyzed potential involvement of the histone N-tails in this process. To this aim, we reconstituted the H2A/H4 N-tailless nucleosomes on human BRCA1 DNA (~100 kb) and compared their positions and sequences with those of the wild-type nucleosomes. In the case of H2A tailless nucleosomes, the AT content of DNA sequences is changed locally at superhelical location (SHL) ±4, while maintaining the same rotational setting as their wild-type counterparts. Conversely, the H4 tailless nucleosomes display widespread changes of the AT content near SHL ±1 and SHL ±2, where the H4 N-tails interact with DNA. Furthermore, a substantial number of H4 tailless nucleosomes exhibit rotational setting opposite to that of the wild-type nucleosomes. Thus, our findings strongly suggest that the histone N-tails are operative in selection of nucleosome positions, which may have wide-ranging implications for epigenetic modulation of chromatin states.

Genome-wide maps of rare and atypical UV photoproducts reveal distinct patterns of damage formation and mutagenesis in yeast chromatin

Ultraviolet (UV) light induces different classes of mutagenic photoproducts in DNA, namely cyclobutane pyrimidine dimers (CPDs), 6-4 photoproducts (6-4PPs), and atypical thymine-adenine photoproducts (TA-PPs). CPD formation is modulated by nucleosomes and transcription factors (TFs), which has important ramifications for Ultraviolet (UV) mutagenesis. How chromatin affects the formation of 6-4PPs and TA-PPs is unclear. Here, we use UV damage endonuclease-sequencing (UVDE-seq) to map these UV photoproducts across the yeast genome. Our results indicate that nucleosomes, the fundamental building block of chromatin, have opposing effects on photoproduct formation. Nucleosomes induce CPDs and 6-4PPs at outward rotational settings in nucleosomal DNA but suppress TA-PPs at these settings. Our data also indicate that DNA binding by different classes of yeast TFs causes lesion-specific hotspots of 6-4PPs or TA-PPs. For example, DNA binding by the TF Rap1 generally suppresses CPD and 6-4PP formation but induces a TA-PP hotspot. Finally, we show that 6-4PP formation is strongly induced at the binding sites of TATA-binding protein (TBP), which is correlated with higher mutation rates in UV-exposed yeast. These results indicate that the formation of 6-4PPs and TA-PPs is modulated by chromatin differently than CPDs and that this may have important implications for UV mutagenesis.

Distinct functions of three chromatin remodelers in activator binding and preinitiation complex assembly

The nucleosome remodeling complexes (CRs) SWI/SNF, RSC, and Ino80C cooperate in evicting or repositioning nucleosomes to produce nucleosome depleted regions (NDRs) at the promoters of many yeast genes induced by amino acid starvation. We analyzed mutants depleted of the catalytic subunits of these CRs for binding of transcriptional activator Gcn4 and recruitment of TATA-binding protein (TBP) during preinitiation complex (PIC) assembly. RSC and Ino80 were found to enhance Gcn4 binding to both UAS elements in NDRs upstream of promoters and to unconventional binding sites within nucleosome-occupied coding sequences; and SWI/SNF contributes to UAS binding when RSC is depleted. All three CRs are actively recruited by Gcn4 to most UAS elements and appear to enhance Gcn4 binding by reducing nucleosome occupancies at the binding motifs, indicating a positive regulatory loop. SWI/SNF acts unexpectedly in WT cells to prevent excessive Gcn4 binding at many UAS elements, indicating a dual mode of action that is modulated by the presence of RSC. RSC and SWI/SNF collaborate to enhance TBP recruitment at Gcn4 target genes, together with Ino80C, in a manner associated with nucleosome eviction at the TBP binding sites. Cooperation among the CRs in TBP recruitment is also evident at the highly transcribed ribosomal protein genes, while RSC and Ino80C act more broadly than SWI/SNF at the majority of other constitutively expressed genes to stimulate this step in PIC assembly. Our findings indicate a complex interplay among the CRs in evicting promoter nucleosomes to regulate activator binding and stimulate PIC assembly.

ATP-dependent chromatin remodelers (CRs), including SWI/SNF and RSC in budding yeast, are thought to stimulate transcription by repositioning or evicting promoter nucleosomes, and we recently implicated the CR Ino80C in this process as well. The relative importance of these CRs in stimulating activator binding and recruitment of TATA-binding protein (TBP) to promoters is incompletely understood. Examining mutants depleted of the catalytic subunits of these CRs, we determined that RSC and Ino80C stimulate binding of transcription factor Gcn4 to nucleosome-depleted regions, or linkers between genic nucleosomes, at multiple target genes activated by Gcn4 in amino acid-starved cells, frequently via evicting nucleosomes from the Gcn4 binding motifs. At some genes, SWI/SNF functionally complements RSC, while opposing RSC at others to limit Gcn4 binding. The CRs in turn are recruited by Gcn4, consistent with a positive feedback loop that enhances Gcn4 binding. The three CRs also cooperate to enhance TBP recruitment, again involving nucleosome depletion, at both Gcn4 target and highly expressed ribosomal protein genes, whereas only RSC and Ino80C act broadly throughout the genome to enhance this key step in preinitiation complex assembly. Our findings illuminate functional cooperation among multiple CRs in regulating activator binding and promoter activation.

Dealing with transcription-blocking DNA damage: Repair mechanisms, RNA polymerase II processing and human disorders

Transcription-blocking DNA lesions (TBLs) in genomic DNA are triggered by a wide variety of DNA-damaging agents. Such lesions cause stalling of elongating RNA polymerase II (RNA Pol II) enzymes and fully block transcription when unresolved. The toxic impact of DNA damage on transcription progression is commonly referred to as transcription stress. In response to RNA Pol II stalling, cells activate and employ transcription-coupled repair (TCR) machineries to repair cytotoxic TBLs and resume transcription. Increasing evidence indicates that the modification and processing of stalled RNA Pol II is an integral component of the cellular response to and the repair of TBLs. If TCR pathways fail, the prolonged stalling of RNA Pol II will impede global replication and transcription as well as block the access of other DNA repair pathways that may act upon the TBL. Consequently, such prolonged stalling will trigger profound genome instability and devastating clinical features. In this review, we will discuss the mechanisms by which various types of TBLs are repaired by distinct TCR pathways and how RNA Pol II processing is regulated during these processes. We will also discuss the clinical consequences of transcription stress and genotype-phenotype correlations of related TCR-deficiency disorders.

Machine learning predicts nucleosome binding modes of transcription factors

Background

Most transcription factors (TFs) compete with nucleosomes to gain access to their cognate binding sites. Recent studies have identified several TF-nucleosome interaction modes including end binding (EB), oriented binding, periodic binding, dyad binding, groove binding, and gyre spanning. However, there are substantial experimental challenges in measuring nucleosome binding modes for thousands of TFs in different species.

Results

We present a computational prediction of the binding modes based on TF protein sequences. With a nested cross-validation procedure, our model outperforms several fine-tuned off-the-shelf machine learning (ML) methods in the multi-label classification task. Our binary classifier for the EB mode performs better than these ML methods with the area under precision-recall curve achieving 75%. The end preference of most TFs is consistent with low nucleosome occupancy around their binding site in GM12878 cells. The nucleosome occupancy data is used as an alternative dataset to confirm the superiority of our EB classifier.

Conclusions

We develop the first ML-based approach for efficient and comprehensive analysis of nucleosome binding modes of TFs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04093-9.

Temporal gene regulation by p53 is associated with the rotational setting of its binding sites in nucleosomes

The tumor suppressor protein p53 is a DNA-binding transcription factor (TF) that, once activated, coordinates the expression of thousands of target genes. Increased p53 binding to gene promoters occurs shortly after p53 activation. Intriguingly, gene transcription exhibits differential kinetics with some genes being induced early (early genes) and others being induced late (late genes). To understand pre-binding factors contributing to the temporal gene regulation by p53, we performed time-course RNA sequencing experiments in human colon cancer cell line HCT116 treated with fluorouracil to identify early and late genes. Published p53 ChIP fragments co-localized with the early or late genes were used to uncover p53 binding sites (BS). We demonstrate that the BS associated with early genes are clustered around gene starts with decreased nucleosome occupancy. DNA analysis shows that these BS are likely exposed on nucleosomal surface if wrapped into nucleosomes, thereby facilitating stable interactions with and fast induction by p53. By contrast, p53 BS associated with late genes are distributed uniformly across the genes with increased nucleosome occupancy. Predicted rotational settings of these BS show limited accessibility. We therefore propose a hypothetical model in which the BS are fully, partially or not accessible to p53 in the nucleosomal context. The partial accessibility of the BS allows subunits of a p53 tetramer to bind, but the resulting p53-DNA complex may not be stable enough to recruit cofactors, which leads to delayed induction. Our work highlights the importance of DNA conformations of p53 BS in gene expression dynamics.

Topological polymorphism of nucleosome fibers and folding of chromatin

We discuss recent observations of polymorphic chromatin packaging at the oligonucleosomal level and compare them with computer simulations. Our computations reveal two topologically different families of two-start 30-nm fiber conformations distinguished by the linker length L; fibers with L ≈ 10n and L ≈ 10n+5 basepairs have DNA linking numbers per nucleosome of ΔLk ≈ -1.5 and -1.0, respectively (where n is a natural number). Although fibers with ΔLk ≈ -1.5 were observed earlier, the topoisomer with ΔLk ≈ -1.0 is novel. These predictions were confirmed experimentally for circular nucleosome arrays with precisely positioned nucleosomes. We suggest that topological polymorphism of chromatin may play a role in transcription, with the {10n+5} fibers producing transcriptionally competent chromatin structures. This hypothesis is consistent with available data for yeast and, partially, for fly. We show that both fiber topoisomers (with ΔLk ≈ -1.5 and -1.0) have to be taken into account to interpret experimental data obtained using new techniques: genome-wide Micro-C, Hi-CO, and RICC-seq, as well as self-association of nucleosome arrays in vitro. The relative stability of these topoisomers is likely to depend on epigenetic histone modifications modulating the strength of internucleosome interactions. Potentially, our findings may reflect a general tendency of functionally distinct parts of the genome to retain topologically different higher-order structures.

Global histone protein surface accessibility in yeast indicates a uniformly loosely packed genome with canonical nucleosomes

Background

The vast majority of methods available to characterize genome-wide chromatin structure exploit differences in DNA accessibility to nucleases or chemical crosslinking. We developed a novel method to gauge genome-wide accessibility of histone protein surfaces within nucleosomes by assessing reactivity of engineered cysteine residues with a thiol-specific reagent, biotin-maleimide (BM).

Results

Yeast nuclei were obtained from cells expressing the histone mutant H2B S116C, in which a cysteine resides near the center of the external flat protein surface of the nucleosome. BM modification revealed that nucleosomes are generally equivalently accessible throughout the S. cerevisiae genome, including heterochromatic regions, suggesting limited, higher-order chromatin structures in which this surface is obstructed by tight nucleosome packing. However, we find that nucleosomes within 500 bp of transcription start sites exhibit the greatest range of accessibility, which correlates with the density of chromatin remodelers. Interestingly, accessibility is not well correlated with RNA polymerase density and thus the level of gene expression. We also investigated the accessibility of cysteine mutations designed to detect exposure of histone surfaces internal to the nucleosome thought to be accessible in actively transcribed genes: H3 102, is at the H2A–H2B dimer/H3–H4 tetramer interface, and H3 A110C, resides at the H3–H3 interface. However, in contrast to the external surface site, we find that neither of these internal sites were found to be appreciably exposed.

Conclusions

Overall, our finding that nucleosomes surfaces within S. cerevisiae chromatin are equivalently accessible genome-wide is consistent with a globally uncompacted chromatin structure lacking substantial higher-order organization. However, we find modest differences in accessibility that correlate with chromatin remodelers but not transcription, suggesting chromatin poised for transcription is more accessible than actively transcribed or intergenic regions. In contrast, we find that two internal sites remain inaccessible, suggesting that such non-canonical nucleosome species generated during transcription are rapidly and efficiently converted to canonical nucleosome structure and thus not widely present in native chromatin.

Causes and consequences of RNA polymerase II stalling during transcript elongation

The journey of RNA polymerase II (Pol II) as it transcribes a gene is anything but a smooth ride. Transcript elongation is discontinuous and can be perturbed by intrinsic regulatory barriers, such as promoter-proximal pausing, nucleosomes, RNA secondary structures and the underlying DNA sequence. More substantial blocking of Pol II translocation can be caused by other physiological circumstances and extrinsic obstacles, including other transcribing polymerases, the replication machinery and several types of DNA damage, such as bulky lesions and DNA double-strand breaks. Although numerous different obstacles cause Pol II stalling or arrest, the cell somehow distinguishes between them and invokes different mechanisms to resolve each roadblock. Resolution of Pol II blocking can be as straightforward as temporary backtracking and transcription elongation factor S-II (TFIIS)-dependent RNA cleavage, or as drastic as premature transcription termination or degradation of polyubiquitylated Pol II and its associated nascent RNA. In this Review, we discuss the current knowledge of how these different Pol II stalling contexts are distinguished by the cell, how they overlap with each other, how they are resolved and how, when unresolved, they can cause genome instability.

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.

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

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

SWI/SNF and RSC cooperate to reposition and evict promoter nucleosomes at highly expressed genes in yeast

The nucleosome remodeling complex RSC functions throughout the yeast genome to set the positions of -1 and +1 nucleosomes and thereby determines the widths of nucleosome-depleted regions (NDRs). The related complex SWI/SNF participates in nucleosome remodeling/eviction and promoter activation at certain yeast genes, including those activated by transcription factor Gcn4, but did not appear to function broadly in establishing NDRs. By analyzing the large cohort of Gcn4-induced genes in mutants lacking the catalytic subunits of SWI/SNF or RSC, we uncovered cooperation between these remodelers in evicting nucleosomes from different locations in the promoter and repositioning the +1 nucleosome downstream to produce wider NDRs-highly depleted of nucleosomes-during transcriptional activation. SWI/SNF also functions on a par with RSC at the most highly transcribed constitutively expressed genes, suggesting general cooperation by these remodelers for maximal transcription. SWI/SNF and RSC occupancies are greatest at the most highly expressed genes, consistent with their cooperative functions in nucleosome remodeling and transcriptional activation. Thus, SWI/SNF acts comparably with RSC in forming wide nucleosome-free NDRs to achieve high-level transcription but only at the most highly expressed genes exhibiting the greatest SWI/SNF occupancies.

p53 binding sites in normal and cancer cells are characterized by distinct chromatin context

The tumor suppressor protein p53 interacts with DNA in a sequence-dependent manner. Thousands of p53 binding sites have been mapped genome-wide in normal and cancer cells. However, the way p53 selectively binds its cognate sites in different types of cells is not fully understood. Here, we performed a comprehensive analysis of 25 published p53 cistromes and identified 3,551 and 6,039 ‘high-confidence’ binding sites in normal and cancer cells, respectively. Our analysis revealed 2 distinct epigenetic features underlying p53-DNA interactions in vivo. First, p53 binding sites are associated with transcriptionally active histone marks (H3K4me3 and H3K36me3) in normal-cell chromatin, but with repressive histone marks (H3K27me3) in cancer-cell chromatin. Second, p53 binding sites in cancer cells are characterized by a lower level of DNA methylation than their counterparts in normal cells, probably related to global hypomethylation in cancers. Intriguingly, regardless of the cell type, p53 sites are highly enriched in the endogenous retroviral elements of the ERV1 family, highlighting the importance of this repeat family in shaping the transcriptional network of p53. Moreover, the p53 sites exhibit an unusual combination of chromatin patterns: high nucleosome occupancy and, at the same time, high sensitivity to DNase I. Our results suggest that p53 can access its target sites in a chromatin environment that is non-permissive to most DNA-binding transcription factors, which may allow p53 to act as a pioneer transcription factor in the context of chromatin.

Genome-wide maps of alkylation damage, repair, and mutagenesis in yeast reveal mechanisms of mutational heterogeneity

DNA base damage is an important contributor to genome instability, but how the formation and repair of these lesions is affected by the genomic landscape and contributes to mutagenesis is unknown. Here, we describe genome-wide maps of DNA base damage, repair, and mutagenesis at single nucleotide resolution in yeast treated with the alkylating agent methyl methanesulfonate (MMS). Analysis of these maps revealed that base excision repair (BER) of alkylation damage is significantly modulated by chromatin, with faster repair in nucleosome-depleted regions, and slower repair and higher mutation density within strongly positioned nucleosomes. Both the translational and rotational settings of lesions within nucleosomes significantly influence BER efficiency; moreover, this effect is asymmetric relative to the nucleosome dyad axis and is regulated by histone modifications. Our data also indicate that MMS-induced mutations at adenine nucleotides are significantly enriched on the nontranscribed strand (NTS) of yeast genes, particularly in BER-deficient strains, due to higher damage formation on the NTS and transcription-coupled repair of the transcribed strand (TS). These findings reveal the influence of chromatin on repair and mutagenesis of base lesions on a genome-wide scale and suggest a novel mechanism for transcription-associated mutation asymmetry, which is frequently observed in human cancers.

Understanding nucleosome dynamics and their links to gene expression and DNA replication

Advances in genomics technology have provided the means to probe myriad chromatin interactions at unprecedented spatial and temporal resolution. This has led to a profound understanding of nucleosome organization within the genome, revealing that nucleosomes are highly dynamic. Nucleosome dynamics are governed by a complex interplay of histone composition, histone post-translational modifications, nucleosome occupancy and positioning within chromatin, which are influenced by numerous regulatory factors, including general regulatory factors, chromatin remodellers, chaperones and polymerases. It is now known that these dynamics regulate diverse cellular processes ranging from gene transcription to DNA replication and repair.

The proto-chromatosome: A fundamental subunit of chromatin?

Eukaryotic DNA is packaged into regularly spaced nucleosomes, resembling beads on a string. Each bead contains ∼147 bp wrapped around a core histone octamer. Linker histone (H1) binds to the linker DNA to drive chromatin folding. Micrococcal nuclease (MNase) digestion studies reveal 2 mono-nucleosomal intermediates: the core particle (∼147 bp) and the chromatosome (∼160 bp; a core particle with additional DNA protected by H1). We have recently developed an improved method for mapping nucleosomes, using exonuclease III to remove residual linker (MNase-Exo-seq). (1) We discovered 2 new intermediate particles corresponding to core particles with ∼7 bp of linker protruding from one side (∼154 bp) or both sides (∼161 bp), which are formed in the absence of H1. We propose that these "proto-chromatosomes" are stabilized by core histone-DNA contacts in the linker, ∼7 bp from the nucleosome boundaries. These contacts may determine the topography of the H1 binding site.

Topological polymorphism of the two-start chromatin fiber

Specific details concerning the spatial organization of nucleosomes in 30 nm fibers remain unknown. To investigate this, we analyzed all stereochemically possible configurations of two-start nucleosome fibers with short DNA linkers L = 13-37 bp (nucleosome repeat length (NRL) = 160-184 bp). Four superhelical parameters-inclination of nucleosomes, twist, rise, and diameter-uniquely describe a regular symmetric fiber. The energy of a fiber is defined as the sum of four terms: elastic energy of the linker DNA, steric repulsion, electrostatics, and a phenomenological (H4 tail-acidic patch) interaction between two stacked nucleosomes. By optimizing the fiber energy with respect to the superhelical parameters, we found two types of topological transition in fibers (associated with the change in inclination angle): one caused by an abrupt 360° change in the linker DNA twisting (change in the DNA linking number, ΔLk = 1), and another caused by overcrossing of the linkers (ΔLk = 2). To the best of our knowledge, this topological polymorphism of the two-start fibers was not reported in the computations published earlier. Importantly, the optimal configurations of the fibers with linkers L = 10n and 10n + 5 bp are characterized by different values of the DNA linking number-that is, they are topologically different. Our results are consistent with experimental observations, such as the inclination 60° to 70° (the angle between the nucleosomal disks and the fiber axis), helical rise, diameter, and left-handedness of the fibers. In addition, we make several testable predictions, among them different degrees of DNA supercoiling in fibers with L = 10n and 10n + 5 bp, different flexibility of the two types of fibers, and a correlation between the local NRL and the level of transcription in different parts of the yeast genome.

The ISW1 and CHD1 ATP-dependent chromatin remodelers compete to set nucleosome spacing in vivo

Adenosine triphosphate-dependent chromatin remodeling machines play a central role in gene regulation by manipulating chromatin structure. Most genes have a nucleosome-depleted region at the promoter and an array of regularly spaced nucleosomes phased relative to the transcription start site. In vitro, the three known yeast nucleosome spacing enzymes (CHD1, ISW1 and ISW2) form arrays with different spacing. We used genome-wide nucleosome sequencing to determine whether these enzymes space nucleosomes differently in vivo We find that CHD1 and ISW1 compete to set the spacing on most genes, such that CHD1 dominates genes with shorter spacing and ISW1 dominates genes with longer spacing. In contrast, ISW2 plays a minor role, limited to transcriptionally inactive genes. Heavily transcribed genes show weak phasing and extreme spacing, either very short or very long, and are depleted of linker histone (H1). Genes with longer spacing are enriched in H1, which directs chromatin folding. We propose that CHD1 directs short spacing, resulting in eviction of H1 and chromatin unfolding, whereas ISW1 directs longer spacing, allowing H1 to bind and condense the chromatin. Thus, competition between the two remodelers to set the spacing on each gene may result in a highly dynamic chromatin structure.

Genome-wide cooperation by HAT Gcn5, remodeler SWI/SNF, and chaperone Ydj1 in promoter nucleosome eviction and transcriptional activation

Chaperones, nucleosome remodeling complexes, and histone acetyltransferases have been implicated in nucleosome disassembly at promoters of particular yeast genes, but whether these cofactors function ubiquitously, as well as the impact of nucleosome eviction on transcription genome-wide, is poorly understood. We used chromatin immunoprecipitation of histone H3 and RNA polymerase II (Pol II) in mutants lacking single or multiple cofactors to address these issues for about 200 genes belonging to the Gcn4 transcriptome, of which about 70 exhibit marked reductions in H3 promoter occupancy on induction by amino acid starvation. Examining four target genes in a panel of mutants indicated that SWI/SNF, Gcn5, the Hsp70 cochaperone Ydj1, and chromatin-associated factor Yta7 are required downstream from Gcn4 binding, whereas Asf1/Rtt109, Nap1, RSC, and H2AZ are dispensable for robust H3 eviction in otherwise wild-type cells. Using ChIP-seq to interrogate all 70 exemplar genes in single, double, and triple mutants implicated Gcn5, Snf2, and Ydj1 in H3 eviction at most, but not all, Gcn4 target promoters, with Gcn5 generally playing the greatest role and Ydj1 the least. Remarkably, these three cofactors cooperate similarly in H3 eviction at virtually all yeast promoters. Defective H3 eviction in cofactor mutants was coupled with reduced Pol II occupancies for the Gcn4 transcriptome and the most highly expressed uninduced genes, but the relative Pol II levels at most genes were unaffected or even elevated. These findings indicate that nucleosome eviction is crucial for robust transcription of highly expressed genes but that other steps in gene activation are more rate-limiting for most other yeast genes.

Topological diversity of chromatin fibers: Interplay between nucleosome repeat length, DNA linking number and the level of transcription

The spatial organization of nucleosomes in 30-nm fibers remains unknown in detail. To tackle this problem, we analyzed all stereochemically possible configurations of two-start chromatin fibers with DNA linkers L = 10–70 bp (nucleosome repeat length NRL = 157–217 bp). In our model, the energy of a fiber is a sum of the elastic energy of the linker DNA, steric repulsion, electrostatics, and the H4 tail-acidic patch interaction between two stacked nucleosomes. We found two families of energetically feasible conformations of the fibers—one observed earlier, and the other novel. The fibers from the two families are characterized by different DNA linking numbers—that is, they are topologically different. Remarkably, the optimal geometry of a fiber and its topology depend on the linker length: the fibers with linkers L = 10n and 10n + 5 bp have DNA linking numbers per nucleosome ΔLk ≈ −1.5 and −1.0, respectively. In other words, the level of DNA supercoiling is directly related to the length of the inter-nucleosome linker in the chromatin fiber (and therefore, to NRL). We hypothesize that this topological polymorphism of chromatin fibers may play a role in the process of transcription, which is known to generate different levels of DNA supercoiling upstream and downstream from RNA polymerase. A genome-wide analysis of the NRL distribution in active and silent yeast genes yielded results consistent with this assumption.

Novel nucleosomal particles containing core histones and linker DNA but no histone H1

Eukaryotic chromosomal DNA is assembled into regularly spaced nucleosomes, which play a central role in gene regulation by determining accessibility of control regions. The nucleosome contains ∼147 bp of DNA wrapped ∼1.7 times around a central core histone octamer. The linker histone, H1, binds both to the nucleosome, sealing the DNA coils, and to the linker DNA between nucleosomes, directing chromatin folding. Micrococcal nuclease (MNase) digests the linker to yield the chromatosome, containing H1 and ∼160 bp, and then converts it to a core particle, containing ∼147 bp and no H1. Sequencing of nucleosomal DNA obtained after MNase digestion (MNase-seq) generates genome-wide nucleosome maps that are important for understanding gene regulation. We present an improved MNase-seq method involving simultaneous digestion with exonuclease III, which removes linker DNA. Remarkably, we discovered two novel intermediate particles containing 154 or 161 bp, corresponding to 7 bp protruding from one or both sides of the nucleosome core. These particles are detected in yeast lacking H1 and in H1-depleted mouse chromatin. They can be reconstituted in vitro using purified core histones and DNA. We propose that these ‘proto-chromatosomes’ are fundamental chromatin subunits, which include the H1 binding site and influence nucleosome spacing independently of H1.

Nucleosome adaptability conferred by sequence and structural variations in histone H2A-H2B dimers

Nucleosome variability is essential for their functions in compacting the chromatin structure and regulation of transcription, replication and cell reprogramming. The DNA molecule in nucleosomes is wrapped around an octamer composed of four types of core histones (H3, H4, H2A, H2B). Nucleosomes represent dynamic entities and may change their conformation, stability and binding properties by employing different sets of histone variants or by becoming post-translationally modified. There are many variants of histones H2A and H2B. Specific H2A and H2B variants may preferentially associate with each other resulting in different combinations of variants and leading to the increased combinatorial complexity of nucleosomes. In addition, the H2A-H2B dimer can be recognized and substituted by chaperones/remodelers as a distinct unit, can assemble independently and is stable during nucleosome unwinding. In this review we discuss how sequence and structural variations in H2A-H2B dimers may provide necessary complexity and confer the nucleosome functional variability.

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.

Heavy transcription of yeast genes correlates with differential loss of histone H2B relative to H4 and queued RNA polymerases

Eukaryotic chromatin is composed of nucleosomes, which contain nearly two coils of DNA wrapped around a central histone octamer. The octamer contains an H3-H4 tetramer and two H2A-H2B dimers. Gene activation is associated with chromatin disruption: a wider nucleosome-depleted region (NDR) at the promoter and reduced nucleosome occupancy over the coding region. Here, we examine the nature of disrupted chromatin after induction, using MNase-seq to map nucleosomes and subnucleosomes, and a refined high-resolution ChIP-seq method to map H4, H2B and RNA polymerase II (Pol II) genome-wide. Over coding regions, induced genes show a differential loss of H2B relative to H4, which correlates with Pol II density and the appearance of subnucleosomes. After induction, Pol II is surprisingly low at the promoter, but accumulates on the gene and downstream of the termination site, implying that dissociation is very slow. Thus, induction-dependent chromatin disruption reflects both eviction of H2A-H2B dimers and the presence of queued Pol II elongation complexes. We propose that slow Pol II dissociation after transcription is a major factor in chromatin disruption and that it may be of critical importance in gene regulation.

Prediction of nucleosome rotational positioning in yeast and human genomes based on sequence-dependent DNA anisotropy

BACKGROUND

An organism's DNA sequence is one of the key factors guiding the positioning of nucleosomes within a cell's nucleus. Sequence-dependent bending anisotropy dictates how DNA is wrapped around a histone octamer. One of the best established sequence patterns consistent with this anisotropy is the periodic occurrence of AT-containing dinucleotides (WW) and GC-containing dinucleotides (SS) in the nucleosomal locations where DNA is bent in the minor and major grooves, respectively. Although this simple pattern has been observed in nucleosomes across eukaryotic genomes, its use for prediction of nucleosome positioning was not systematically tested.

RESULTS

We present a simple computational model, termed the W/S scheme, implementing this pattern, without using any training data. This model accurately predicts the rotational positioning of nucleosomes both in vitro and in vivo, in yeast and human genomes. About 65 - 75% of the experimentally observed nucleosome positions are predicted with the precision of one to two base pairs. The program is freely available at http://people.rit.edu/fxcsbi/WS_scheme/. We also introduce a simple and efficient way to compare the performance of different models predicting the rotational positioning of nucleosomes.

CONCLUSIONS

This paper presents the W/S scheme to achieve accurate prediction of rotational positioning of nucleosomes, solely based on the sequence-dependent anisotropic bending of nucleosomal DNA. This method successfully captures DNA features critical for the rotational positioning of nucleosomes, and can be further improved by incorporating additional terms related to the translational positioning of nucleosomes in a species-specific manner.

RSC-dependent constructive and destructive interference between opposing arrays of phased nucleosomes in yeast

RSC and SWI/SNF are related ATP-dependent chromatin remodeling machines that move nucleosomes, regulating access to DNA. We addressed their roles in nucleosome phasing relative to transcription start sites in yeast. SWI/SNF has no effect on phasing at the global level. In contrast, RSC depletion results in global nucleosome repositioning: Both upstream and downstream nucleosomal arrays shift toward the nucleosome-depleted region (NDR), with no change in spacing, resulting in a narrower and partly filled NDR. The global picture of RSC-depleted chromatin represents the average of a range of chromatin structures, with most genes showing a shift of the +1 or the -1 nucleosome into the NDR. Using RSC ChIP data reported by others, we show that RSC occupancy is highest on the coding regions of heavily transcribed genes, though not at their NDRs. We propose that RSC has a role in restoring chromatin structure after transcription. Analysis of gene pairs in different orientations demonstrates that phasing patterns reflect competition between phasing signals emanating from neighboring NDRs. These signals may be in phase, resulting in constructive interference and a regular array, or out of phase, resulting in destructive interference and fuzzy positioning. We propose a modified barrier model, in which a stable complex located at the NDR acts as a bidirectional phasing barrier. In RSC-depleted cells, this barrier has a smaller footprint, resulting in narrower NDRs. Thus, RSC plays a critical role in organizing yeast chromatin.

Nucleosomes shape DNA polymorphism and divergence

An estimated 80% of genomic DNA in eukaryotes is packaged as nucleosomes, which, together with the remaining interstitial linker regions, generate higher order chromatin structures [1]. Nucleosome sequences isolated from diverse organisms exhibit ∼10 bp periodic variations in AA, TT and GC dinucleotide frequencies. These sequence elements generate intrinsically curved DNA and help establish the histone-DNA interface. We investigated an important unanswered question concerning the interplay between chromatin organization and genome evolution: do the DNA sequence preferences inherent to the highly conserved histone core exert detectable natural selection on genomic divergence and polymorphism? To address this hypothesis, we isolated nucleosomal DNA sequences from Drosophila melanogaster embryos and examined the underlying genomic variation within and between species. We found that divergence along the D. melanogaster lineage is periodic across nucleosome regions with base changes following preferred nucleotides, providing new evidence for systematic evolutionary forces in the generation and maintenance of nucleosome-associated dinucleotide periodicities. Further, Single Nucleotide Polymorphism (SNP) frequency spectra show striking periodicities across nucleosomal regions, paralleling divergence patterns. Preferred alleles occur at higher frequencies in natural populations, consistent with a central role for natural selection. These patterns are stronger for nucleosomes in introns than in intergenic regions, suggesting selection is stronger in transcribed regions where nucleosomes undergo more displacement, remodeling and functional modification. In addition, we observe a large-scale (∼180 bp) periodic enrichment of AA/TT dinucleotides associated with nucleosome occupancy, while GC dinucleotide frequency peaks in linker regions. Divergence and polymorphism data also support a role for natural selection in the generation and maintenance of these super-nucleosomal patterns. Our results demonstrate that nucleosome-associated sequence periodicities are under selective pressure, implying that structural interactions between nucleosomes and DNA sequence shape sequence evolution, particularly in introns.

In eukaryotic cells, the majority of DNA is packaged in nucleosomes comprised of ∼147 bp of DNA wound tightly around the highly conserved histone octamer. Nucleosomal DNA from diverse organisms shows an anti-correlated ∼10 bp periodicity of AT-rich and GC-rich dinucleotides. These sequence features influence DNA bending and shape, facilitating structural interactions. We asked whether natural selection mediated through the periodic sequence preferences of nucleosomes shapes the evolution of non-protein-coding regions of D. melanogaster by examining the inter- and intra-species genomic variation relative to these fundamental chromatin building blocks. The sequence changes across nucleosome-bound regions on the melanogaster lineage mirror the observed nucleosome dinucleotide periodicities. Importantly, we show that the frequencies of polymorphisms in natural populations vary across these regions, paralleling divergence, with higher frequencies of preferred alleles. These patterns are most evident for intronic regions and indicate that non-protein coding regions are evolving toward sequences that facilitate the canonical association with the histone core. This result is consistent with the hypothesis that interactions between DNA and the core have systematic impacts on function that are subject to natural selection and are not solely due to mutational bias. These ubiquitous interactions with the histone core partially account for the evolutionary constraint observed in unannotated genomic regions, and may drive broad changes in base composition.

Rotational positioning of nucleosomes facilitates selective binding of p53 to response elements associated with cell cycle arrest

The tumor suppressor protein p53 exhibits high affinity to the response elements regulating cell cycle arrest genes (CCA-sites), but relatively low affinity to the sites associated with apoptosis (Apo-sites). This in vivo tendency cannot be explained solely by the p53-DNA binding constants measured in vitro. Since p53 can bind nucleosomal DNA, we sought to understand if the two groups of p53 sites differ in their accessibility when embedded in nucleosomes. To this aim, we analyzed the sequence-dependent bending anisotropy of human genomic DNA containing p53 sites. For the 20 CCA-sites, we calculated rotational positioning patterns predicting that most of the sites are exposed on the nucleosomal surface. This is consistent with experimentally observed positioning of human nucleosomes. Remarkably, the sequence-dependent DNA anisotropy of both the p53 sites and flanking DNA work in concert producing strong positioning signals. By contrast, both the predicted and observed rotational settings of the 38 Apo-sites in nucleosomes suggest that many of these sites are buried inside, thus preventing immediate p53 recognition and delaying gene induction. The distinct chromatin organization of the CCA response elements appears to be one of the key factors facilitating p53-DNA binding and subsequent activation of genes associated with cell cycle arrest.