Genome-wide nucleosome and transcription factor responses to genetic perturbations reveal chromatin-mediated mechanisms of transcriptional regulation.

Epigenetic mechanisms contribute to gene regulation by altering chromatin accessibility through changes in transcription factor (TF) and nucleosome occupancy throughout the genome. Despite numerous studies focusing on changes in gene expression, the intricate chromatin-mediated regulatory code remains largely unexplored on a comprehensive scale. We address this by employing a factor-agnostic, reverse-genetics approach that uses MNase-seq to capture genome-wide TF and nucleosome occupancies in response to the individual deletion of 201 transcriptional regulators in Saccharomyces cerevisiae, thereby assaying nearly one million mutant-gene interactions. We develop a principled approach to identify and quantify chromatin changes genome-wide, observing differences in TF and nucleosome occupancy that recapitulate well-established pathways identified by gene expression data. We also discover distinct chromatin signatures associated with the up- and downregulation of genes, and use these signatures to reveal regulatory mechanisms previously unexplored in expression-based studies. Finally, we demonstrate that chromatin features are predictive of transcriptional activity and leverage these features to reconstruct chromatin-based transcriptional regulatory networks. Overall, these results illustrate the power of an approach combining genetic perturbation with high-resolution epigenomic profiling; the latter enables a close examination of the interplay between TFs and nucleosomes genome-wide, providing a deeper, more mechanistic understanding of the complex relationship between chromatin organization and transcription.

Spatiotemporal kinetics of CAF-1-dependent chromatin maturation ensures transcription fidelity during S-phase.

Proper maintenance of epigenetic information after replication is dependent on the rapid assembly and maturation of chromatin. Chromatin Assembly Complex 1 (CAF-1) is a conserved histone chaperone that deposits (H3-H4)2 tetramers as part of the replication-dependent chromatin assembly process. Loss of CAF-1 leads to a delay in chromatin maturation, albeit with minimal impact on steady-state chromatin structure. However, the mechanisms by which CAF-1 mediates the deposition of (H3-H4)2 tetramers and the phenotypic consequences of CAF-1-associated assembly defects are not well understood. We used nascent chromatin occupancy profiling to track the spatiotemporal kinetics of chromatin maturation in both wild-type (WT) and CAF-1 mutant yeast cells. Our results show that loss of CAF-1 leads to a heterogeneous rate of nucleosome assembly, with some nucleosomes maturing at near WT kinetics and others showing significantly slower maturation kinetics. The slow-to-mature nucleosomes are enriched in intergenic and poorly transcribed regions, suggesting that transcription-dependent assembly mechanisms can reset the slow-to-mature nucleosomes following replication. Nucleosomes with slow maturation kinetics are also associated with poly(dA:dT) sequences, which implies that CAF-1 deposits histones in a manner that counteracts resistance from the inflexible DNA sequence, promoting the formation of histone octamers as well as ordered nucleosome arrays. In addition, we show that the delay in chromatin maturation is accompanied by a transient and S-phase-specific loss of gene silencing and transcriptional regulation, revealing that the DNA replication program can directly shape the chromatin landscape and modulate gene expression through the process of chromatin maturation.

Spatiotemporal kinetics of CAF-1-dependent chromatin maturation ensures transcription fidelity during S-phase.

Proper maintenance of epigenetic information after replication is dependent on the rapid assembly and maturation of chromatin. Chromatin Assembly Complex 1 (CAF-1) is a conserved histone chaperone that deposits (H3-H4)2 tetramers as part of the replication-dependent chromatin assembly process. Loss of CAF-1 leads to a delay in chromatin maturation, albeit with minimal impact on steady-state chromatin structure. However, the mechanisms by which CAF-1 mediates the deposition of (H3-H4)2 tetramers and the phenotypic consequences of CAF-1-associated assembly defects are not well understood. We used nascent chromatin occupancy profiling to track the spatiotemporal kinetics of chromatin maturation in both wild-type (WT) and CAF-1 mutant yeast cells. Our results show that loss of CAF-1 leads to a heterogeneous rate of nucleosome assembly, with some nucleosomes maturing at near WT kinetics and others exhibiting significantly slower maturation kinetics. The slow-to-mature nucleosomes are enriched in intergenic and poorly transcribed regions, suggesting that transcription-dependent assembly mechanisms can reset the slow-to-mature nucleosomes following replication. Nucleosomes with slow maturation kinetics are also associated with poly(dA:dT) sequences, which implies that CAF-1 deposits histones in a manner that counteracts resistance from the inflexible DNA sequence, promoting the formation of histone octamers as well as ordered nucleosome arrays. In addition, we demonstrate that the delay in chromatin maturation is accompanied by a transient and S-phase specific loss of gene silencing and transcriptional regulation, revealing that the DNA replication program can directly shape the chromatin landscape and modulate gene expression through the process of chromatin maturation.

Cell-Cycle-Dependent Chromatin Dynamics at Replication Origins.

Origins of DNA replication are specified by the ordered recruitment of replication factors in a cell-cycle-dependent manner. The assembly of the pre-replicative complex in G1 and the pre-initiation complex prior to activation in S phase are well characterized; however, the interplay between the assembly of these complexes and the local chromatin environment is less well understood. To investigate the dynamic changes in chromatin organization at and surrounding replication origins, we used micrococcal nuclease (MNase) to generate genome-wide chromatin occupancy profiles of nucleosomes, transcription factors, and replication proteins through consecutive cell cycles in Saccharomyces cerevisiae. During each G1 phase of two consecutive cell cycles, we observed the downstream repositioning of the origin-proximal +1 nucleosome and an increase in protected DNA fragments spanning the ARS consensus sequence (ACS) indicative of pre-RC assembly. We also found that the strongest correlation between chromatin occupancy at the ACS and origin efficiency occurred in early S phase, consistent with the rate-limiting formation of the Cdc45-Mcm2-7-GINS (CMG) complex being a determinant of origin activity. Finally, we observed nucleosome disruption and disorganization emanating from replication origins and traveling with the elongating replication forks across the genome in S phase, likely reflecting the disassembly and assembly of chromatin ahead of and behind the replication fork, respectively. These results provide insights into cell-cycle-regulated chromatin dynamics and how they relate to the regulation of origin activity.

Disruption of origin chromatin structure by helicase activation in the absence of DNA replication.

Prior to initiation of DNA replication, the eukaryotic helicase, Mcm2-7, must be activated to unwind DNA at replication start sites in early S phase. To study helicase activation within origin chromatin, we constructed a conditional mutant of the polymerase α subunit Cdc17 (or Pol1) to prevent priming and block replication. Recovery of these cells at permissive conditions resulted in the generation of unreplicated gaps at origins, likely due to helicase activation prior to replication initiation. We used micrococcal nuclease (MNase)-based chromatin occupancy profiling under restrictive conditions to study chromatin dynamics associated with helicase activation. Helicase activation in the absence of DNA replication resulted in the disruption and disorganization of chromatin, which extends up to 1 kb from early, efficient replication origins. The CMG holohelicase complex also moves the same distance out from the origin, producing single-stranded DNA that activates the intra-S-phase checkpoint. Loss of the checkpoint did not regulate the progression and stalling of the CMG complex but rather resulted in the disruption of chromatin at both early and late origins. Finally, we found that the local sequence context regulates helicase progression in the absence of DNA replication, suggesting that the helicase is intrinsically less processive when uncoupled from replication.

Stochastic initiation of DNA replication across the human genome.

Wang et al. (2021) comprehensively map DNA replication initiation events across the human genome using single-molecule optical resolution mapping and find that initiation events are randomly distributed across broad initiation zones that are only utilized in a stochastic fashion across a population of cells.

RoboCOP: jointly computing chromatin occupancy profiles for numerous factors from chromatin accessibility data.

Chromatin is a tightly packaged structure of DNA and protein within the nucleus of a cell. The arrangement of different protein complexes along the DNA modulates and is modulated by gene expression. Measuring the binding locations and occupancy levels of different transcription factors (TFs) and nucleosomes is therefore crucial to understanding gene regulation. Antibody-based methods for assaying chromatin occupancy are capable of identifying the binding sites of specific DNA binding factors, but only one factor at a time. In contrast, epigenomic accessibility data like MNase-seq, DNase-seq, and ATAC-seq provide insight into the chromatin landscape of all factors bound along the genome, but with little insight into the identities of those factors. Here, we present RoboCOP, a multivariate state space model that integrates chromatin accessibility data with nucleotide sequence to jointly compute genome-wide probabilistic scores of nucleosome and TF occupancy, for hundreds of different factors. We apply RoboCOP to MNase-seq and ATAC-seq data to elucidate the protein-binding landscape of nucleosomes and 150 TFs across the yeast genome, and show that our model makes better predictions than existing methods. We also compute a chromatin occupancy profile of the yeast genome under cadmium stress, revealing chromatin dynamics associated with transcriptional regulation.

Linking the dynamics of chromatin occupancy and transcription with predictive models.

Though the sequence of the genome within each eukaryotic cell is essentially fixed, it exists within a complex and changing chromatin state. This state is determined, in part, by the dynamic binding of proteins to the DNA. These proteins-including histones, transcription factors (TFs), and polymerases-interact with one another, the genome, and other molecules to allow the chromatin to adopt one of exceedingly many possible configurations. Understanding how changing chromatin configurations associate with transcription remains a fundamental research problem. We sought to characterize at high spatiotemporal resolution the dynamic interplay between transcription and chromatin in response to cadmium stress. Whereas gene regulatory responses to environmental stress in yeast have been studied, how the chromatin state changes and how those changes connect to gene regulation remain unexplored. By combining MNase-seq and RNA-seq data, we found chromatin signatures of transcriptional activation and repression involving both nucleosomal and TF-sized DNA-binding factors. Using these signatures, we identified associations between chromatin dynamics and transcriptional regulation, not only for known cadmium response genes, but across the entire genome, including antisense transcripts. Those associations allowed us to develop generalizable models that predict dynamic transcriptional responses on the basis of dynamic chromatin signatures.

Local nucleosome dynamics and eviction following a double-strand break are reversible by NHEJ-mediated repair in the absence of DNA replication.

We interrogated at nucleotide resolution the spatiotemporal order of chromatin changes that occur immediately following a site-specific double-strand break (DSB) upstream of the PHO5 locus and its subsequent repair by nonhomologous end joining (NHEJ). We observed the immediate eviction of a nucleosome flanking the break and the repositioning of adjacent nucleosomes away from the break. These early chromatin events were independent of the end-processing Mre11-Rad50-Xrs2 (MRX) complex and preceded the MRX-dependent broad eviction of histones and DNA end-resectioning that extends up to ∼8 kb away from the break. We also examined the temporal dynamics of NHEJ-mediated repair in a G1-arrested population. Concomitant with DSB repair by NHEJ, we observed the redeposition and precise repositioning of nucleosomes at their originally occupied positions. This re-establishment of the prelesion chromatin landscape suggests that a DNA replication-independent mechanism exists to preserve epigenome organization following DSB repair.

Capturing the primordial Kras mutation initiating urethane carcinogenesis.

The environmental carcinogen urethane exhibits a profound specificity for pulmonary tumors driven by an oncogenic Q61L/R mutation in the gene Kras. Similarly, the frequency, isoform, position, and substitution of oncogenic RAS mutations are often unique to human cancers. To elucidate the principles underlying this RAS mutation tropism of urethane, we adapted an error-corrected, high-throughput sequencing approach to detect mutations in murine Ras genes at great sensitivity. This analysis not only captured the initiating Kras mutation days after urethane exposure, but revealed that the sequence specificity of urethane mutagenesis, coupled with transcription and isoform locus, to be major influences on the extreme tropism of this carcinogen.

RoboCOP: Multivariate State Space Model Integrating Epigenomic Accessibility Data to Elucidate Genome-Wide Chromatin Occupancy.

Chromatin is the tightly packaged structure of DNA and protein within the nucleus of a cell. The arrangement of different protein complexes along the DNA modulates and is modulated by gene expression. Measuring the binding locations and level of occupancy of different transcription factors (TFs) and nucleosomes is therefore crucial to understanding gene regulation. Antibody-based methods for assaying chromatin occupancy are capable of identifying the binding sites of specific DNA binding factors, but only one factor at a time. On the other hand, epigenomic accessibility data like ATAC-seq, DNase-seq, and MNase-seq provide insight into the chromatin landscape of all factors bound along the genome, but with minimal insight into the identities of those factors. Here, we present RoboCOP, a multivariate state space model that integrates chromatin information from epigenomic accessibility data with nucleotide sequence to compute genome-wide probabilistic scores of nucleosome and TF occupancy, for hundreds of different factors at once. RoboCOP can be applied to any epigenomic dataset that provides quantitative insight into chromatin accessibility in any organism, but here we apply it to MNase-seq data to elucidate the protein-binding landscape of nucleosomes and 150 TFs across the yeast genome. Using available protein-binding datasets from the literature, we show that our model more accurately predicts the binding of these factors genome-wide.

Nascent chromatin occupancy profiling reveals locus- and factor-specific chromatin maturation dynamics behind the DNA replication fork.

Proper regulation and maintenance of the epigenome is necessary to preserve genome function. However, in every cell division, the epigenetic state is disassembled and then reassembled in the wake of the DNA replication fork. Chromatin restoration on nascent DNA is a complex and regulated process that includes nucleosome assembly and remodeling, deposition of histone variants, and the re-establishment of transcription factor binding. To study the genome-wide dynamics of chromatin restoration behind the DNA replication fork, we developed nascent chromatin occupancy profiles (NCOPs) to comprehensively profile nascent and mature chromatin at nucleotide resolution. Although nascent chromatin is inherently less organized than mature chromatin, we identified locus-specific differences in the kinetics of chromatin maturation that were predicted by the epigenetic landscape, including the histone variant H2AZ, which marked loci with rapid maturation kinetics. The chromatin maturation at origins of DNA replication was dependent on whether the origin underwent initiation or was passively replicated from distal-originating replication forks, suggesting distinct chromatin assembly mechanisms surrounding activated and disassembled prereplicative complexes. Finally, we identified sites that were only occupied transiently by DNA-binding factors following passage of the replication fork, which may provide a mechanism for perturbations of the DNA replication program to shape the regulatory landscape of the genome.

Sir2 suppresses transcription-mediated displacement of Mcm2-7 replicative helicases at the ribosomal DNA repeats.

Repetitive DNA sequences within eukaryotic heterochromatin are poorly transcribed and replicate late in S-phase. In Saccharomyces cerevisiae, the histone deacetylase Sir2 is required for both transcriptional silencing and late replication at the repetitive ribosomal DNA arrays (rDNA). Despite the widespread association between transcription and replication timing, it remains unclear how transcription might impinge on replication, or vice versa. Here we show that, when silencing of an RNA polymerase II (RNA Pol II)-transcribed non-coding RNA at the rDNA is disrupted by SIR2 deletion, RNA polymerase pushes and thereby relocalizes replicative Mcm2-7 helicases away from their loading sites to an adjacent region with low nucleosome occupancy, and this relocalization is associated with increased rDNA origin efficiency. Our results suggest a model in which two of the major defining features of heterochromatin, transcriptional silencing and late replication, are mechanistically linked through suppression of polymerase-mediated displacement of replication initiation complexes.

Chromatin conformation and transcriptional activity are permissive regulators of DNA replication initiation in Drosophila.

Chromatin structure has emerged as a key contributor to spatial and temporal control over the initiation of DNA replication. However, despite genome-wide correlations between early replication of gene-rich, accessible euchromatin and late replication of gene-poor, inaccessible heterochromatin, a causal relationship between chromatin structure and replication initiation remains elusive. Here, we combined histone gene engineering and whole-genome sequencing in Drosophila to determine how perturbing chromatin structure affects replication initiation. We found that most pericentric heterochromatin remains late replicating in H3K9R mutants, even though H3K9R pericentric heterochromatin is depleted of HP1a, more accessible, and transcriptionally active. These data indicate that HP1a loss, increased chromatin accessibility, and elevated transcription do not result in early replication of heterochromatin. Nevertheless, a small amount of pericentric heterochromatin with increased accessibility replicates earlier in H3K9R mutants. Transcription is de-repressed in these regions of advanced replication but not in those regions of the H3K9R mutant genome that replicate later, suggesting that transcriptional repression may contribute to late replication. We also explored relationships among chromatin, transcription, and replication in euchromatin by analyzing H4K16R mutants. In Drosophila, the X Chromosome gene expression is up-regulated twofold and replicates earlier in XY males than it does in XX females. We found that H4K16R mutation prevents normal male development and abrogates hyperexpression and earlier replication of the male X, consistent with previously established genome-wide correlations between transcription and early replication. In contrast, H4K16R females are viable and fertile, indicating that H4K16 modification is dispensable for genome replication and gene expression.

Nucleosomes influence multiple steps during replication initiation.

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

Temporal association of ORCA/LRWD1 to late-firing origins during G1 dictates heterochromatin replication and organization.

DNA replication requires the recruitment of a pre-replication complex facilitated by Origin Recognition Complex (ORC) onto the chromatin during G1 phase of the cell cycle. The ORC-associated protein (ORCA/LRWD1) stabilizes ORC on chromatin. Here, we evaluated the genome-wide distribution of ORCA using ChIP-seq during specific time points of G1. ORCA binding sites on the G1 chromatin are dynamic and temporally regulated. ORCA association to specific genomic sites decreases as the cells progressed towards S-phase. The majority of the ORCA-bound sites represent replication origins that also associate with the repressive chromatin marks H3K9me3 and methylated-CpGs, consistent with ORCA-bound origins initiating DNA replication late in S-phase. Further, ORCA directly associates with the repressive marks and interacts with the enzymes that catalyze these marks. Regions that associate with both ORCA and H3K9me3, exhibit diminished H3K9 methylation in ORCA-depleted cells, suggesting a role for ORCA in recruiting the H3K9me3 mark at certain genomic loci. Similarly, DNA methylation is altered at ORCA-occupied sites in cells lacking ORCA. Furthermore, repressive chromatin marks influence ORCA's binding on chromatin. We propose that ORCA coordinates with the histone and DNA methylation machinery to establish a repressive chromatin environment at a subset of origins, which primes them for late replication.

DNA replication origins-where do we begin?

For more than three decades, investigators have sought to identify the precise locations where DNA replication initiates in mammalian genomes. The development of molecular and biochemical approaches to identify start sites of DNA replication (origins) based on the presence of defining and characteristic replication intermediates at specific loci led to the identification of only a handful of mammalian replication origins. The limited number of identified origins prevented a comprehensive and exhaustive search for conserved genomic features that were capable of specifying origins of DNA replication. More recently, the adaptation of origin-mapping assays to genome-wide approaches has led to the identification of tens of thousands of replication origins throughout mammalian genomes, providing an unprecedented opportunity to identify both genetic and epigenetic features that define and regulate their distribution and utilization. Here we summarize recent advances in our understanding of how primary sequence, chromatin environment, and nuclear architecture contribute to the dynamic selection and activation of replication origins across diverse cell types and developmental stages.

Noncoding Transcription Is a Driving Force for Nucleosome Instability in spt16 Mutant Cells.

FACT (facilitates chromatin transcription) consists of two essential subunits, Spt16 and Pob3, and functions as a histone chaperone. Mutation of spt16 results in a global loss of nucleosomes as well as aberrant transcription. Here, we show that the majority of nucleosome changes upon Spt16 depletion are alterations in nucleosome fuzziness and position shift. Most nucleosomal changes are suppressed by the inhibition of RNA polymerase II (Pol II) activity. Surprisingly, a small subgroup of nucleosome changes is resistant to transcriptional inhibition. Notably, Spt16 and distinct histone modifications are enriched at this subgroup of nucleosomes. We also report 1,037 Spt16-suppressed noncoding transcripts (SNTs) and found that the SNT start sites are enriched with the subgroup of nucleosomes resistant to Pol II inhibition. Finally, the nucleosomes at genes overlapping SNTs are more susceptible to changes upon Spt16 depletion than those without SNTs. Taken together, our results support a model in which Spt16 has a role in maintaining local nucleosome stability to inhibit initiation of SNT transcription, which once initiated drives additional nucleosome loss upon Spt16 depletion.

Methylation of histone H4 lysine 20 by PR-Set7 ensures the integrity of late replicating sequence domains in Drosophila.

The methylation state of lysine 20 on histone H4 (H4K20) has been linked to chromatin compaction, transcription, DNA repair and DNA replication. Monomethylation of H4K20 (H4K20me1) is mediated by the cell cycle-regulated histone methyltransferase PR-Set7. PR-Set7 depletion in mammalian cells results in defective S phase progression and the accumulation of DNA damage, which has been partially attributed to defects in origin selection and activation. However, these studies were limited to only a handful of mammalian origins, and it remains unclear how PR-Set7 and H4K20 methylation impact the replication program on a genomic scale. We employed genetic, cytological, and genomic approaches to better understand the role of PR-Set7 and H4K20 methylation in regulating DNA replication and genome stability in Drosophila cells. We find that deregulation of H4K20 methylation had no impact on origin activation throughout the genome. Instead, depletion of PR-Set7 and loss of H4K20me1 results in the accumulation of DNA damage and an ATR-dependent cell cycle arrest. Coincident with the ATR-dependent cell cycle arrest, we find increased DNA damage that is specifically limited to late replicating regions of the Drosophila genome, suggesting that PR-Set7-mediated monomethylation of H4K20 is critical for maintaining the genomic integrity of late replicating domains.