Single-nucleosome Resolution Research Articles

Abstract 1648: RB represses cohesin-dependent loop formation and activates E2F-independent transcription

Abstract The retinoblastoma protein (RB) suppresses the activity of the E2F transcription factor family, controlling cellular proliferation. Recent studies, however, indicate that RB's role in chromatin organization, which is not yet fully understood, might be distinct from this E2F-dependent regulation. Employing chromosome conformation capture at a single nucleosome resolution, we identified RB as a cell cycle-regulated repressor of cohesin-dependent loop formation at topologically associating domain (TAD) boundaries. RB depletion increased the number and size of cohesin-dependent loops and strengthened topologically associating domains (TADs). This phenomenon was specific to the G1 phase and was not observed in the S phase, indicating that this novel RB function is regulated by the cell cycle. Mechanistically, RB showed extensive colocalization with cohesin in the human genome, and it impacted cohesin’s distribution on the chromatin. Active RB reduced cohesin from RB-bound TAD boundaries and decreased cohesin activity therein, as assessed by the increased K105/106 acetylation of the cohesin subunit SMC3. This led to reduced insulation in chromosome conformation capture assays. Importantly, by weakening the insulation activity of the adjacent insulators, RB non-canonically enhanced the expression of non-E2F target genes salient for cell adhesion and extracellular matrix organization. When RB was lost, cells showed a more rapid cellular detachment rate and an elevated migration rate, indicating that this novel RB function controls specific transcriptional programs rather than arbitrary genes to regulate the process of cell adhesion and migration. Overall, we conclude that RB has a central role in the interplay between cell cycle and chromatin organization, by repressing the cohesin-dependent loop formation at TAD boundaries. This RB function safeguards E2F-independent transcriptional programs driven by active enhancers and helps maintain cellular adhesion to the extracellular matrix. Citation Format: Hanjun Lee, Ioanna-Maria Gkotinakou, Badri Krishnan, Nicholas J. Dyson, Michael S. Lawrence, Ioannis Sanidas. RB represses cohesin-dependent loop formation and activates E2F-independent transcription [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 1648.

Hi-TrAC reveals division of labor of transcription factors in organizing chromatin loops

The three-dimensional genomic structure plays a critical role in gene expression, cellular differentiation, and pathological conditions. It is pivotal to elucidate fine-scale chromatin architectures, especially interactions of regulatory elements, to understand the temporospatial regulation of gene expression. In this study, we report Hi-TrAC as a proximity ligation-free, robust, and sensitive technique to profile genome-wide chromatin interactions at high-resolution among regulatory elements. Hi-TrAC detects chromatin looping among accessible regions at single nucleosome resolution. With almost half-million identified loops, we reveal a comprehensive interaction network of regulatory elements across the genome. After integrating chromatin binding profiles of transcription factors, we discover that cohesin complex and CTCF are responsible for organizing long-range chromatin loops, related to domain formation; whereas ZNF143 and HCFC1 are involved in structuring short-range chromatin loops between regulatory elements, which directly regulate gene expression. Thus, we introduce a methodology to identify a delicate and comprehensive network of cis-regulatory elements, revealing the complexity and a division of labor of transcription factors in organizing chromatin loops for genome organization and gene expression.

Cohesin and CTCF complexes mediate contacts in chromatin loops depending on nucleosome positions

The spatial organization of the eukaryotic genome plays an important role in regulating transcriptional activity. In the nucleus, chromatin forms loops that assemble into fundamental units called topologically associating domains that facilitate or inhibit long-range contacts. These loops are formed and held together by the ring-shaped cohesin protein complex, and this can involve binding of CCCTC-binding factor (CTCF). High-resolution conformation capture experiments provide the frequency at which two DNA fragments physically associate in three-dimensional space. However, technical limitations of this approach, such as low throughput, low resolution, or noise in contact maps, make data interpretation and identification of chromatin intraloop contacts, e.g., between distal regulatory elements and their target genes, challenging. Herein, an existing coarse-grained model of chromatin at single-nucleosome resolution was extended by integrating potentials describing CTCF and cohesin. We performed replica-exchange Monte Carlo simulations with regularly spaced nucleosomes and experimentally determined nucleosome positions in the presence of cohesin-CTCF, as well as depleted systems as controls. In fully extruded loops caused by the presence of cohesin and CTCF, the number of contacts within the formed loops was increased. The number and types of these contacts were impacted by the nucleosome distribution and loop size. Microloops were observed within cohesin-mediated loops due to thermal fluctuations without additional influence of other factors, and the number, size, and shape of microloops were determined by nucleosome distribution and loop size. Nucleosome positions directly affect the spatial structure and contact probability within a loop, with presumed consequences for transcriptional activity.

Stochastic models of nucleosome dynamics reveal regulatory rules of stimulus-induced epigenome remodeling.

The genomic positions of nucleosomes are a defining feature of the cell's epigenomic state, but signal-dependent transcription factors (SDTFs), upon activation, bind to specific genomic locations and modify nucleosome positioning. Here we leverage SDTFs as perturbation probes to learn about nucleosome dynamics in living cells. We develop Markov models of nucleosome dynamics and fit them to time course sequencing data of DNA accessibility. We find that (1) the dynamics of DNA unwrapping are significantly slower in cells than reported from cell-free experiments, (2) only models with cooperativity in wrapping and unwrapping fit the available data, (3) SDTF activity produces the highest eviction probability when its binding site is adjacent to but not on the nucleosome dyad, and (4) oscillatory SDTF activity results in high location variability. Our work uncovers the regulatory rules governing SDTF-induced nucleosome dynamics in live cells, which can predict chromatin accessibility alterations during inflammation at single-nucleosome resolution.

Prediction of histone post-translational modification patterns based on nascent transcription data.

The role of histone modifications in transcription remains incompletely understood. Here we examine the relationship between histone modifications and transcription using experimental perturbations combined with sensitive machine-learning tools. Transcription predicted the variation in active histone marks and complex chromatin states, like bivalent promoters, down to single-nucleosome resolution and at an accuracy that rivaled the correspondence between independent ChIP-seq experiments. Blocking transcription rapidly removed two punctate marks, H3K4me3 and H3K27ac, from chromatin indicating that transcription is required for active histone modifications. Transcription was also required for maintenance of H3K27me3 consistent with a role for RNA in recruiting PRC2. A subset of DNase-I hypersensitive sites were refractory to prediction, precluding models where transcription initiates pervasively at any open chromatin. Our results, in combination with past literature, support a model in which active histone modifications serve a supportive, rather than an essential regulatory, role in transcription.

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.

NUCLIZE for quantifying epigenome: generating histone modification data at single-nucleosome resolution using genuine nucleosome positions

BackgroundDefining histone modification at single-nucleosome resolution provides accurate epigenomic information in individual nucleosomes. However, most of histone modification data deposited in current databases, such as ENCODE and Roadmap, have low resolution with peaks of several kilo-base pairs (kb), which due to the technical defects of regular ChIP-Seq technology.ResultsTo generate histone modification data at single-nucleosome resolution, we developed a novel approach, NUCLIZE, using synergistic analyses of histone modification data from ChIP-Seq and high-resolution nucleosome mapping data from native MNase-Seq. With this approach, we generated quantitative epigenomics data of single and multivalent histone modification marks in each nucleosome. We found that the dominant trivalent histone mark (H3K4me3/H3K9ac/H3K27ac) and others showed defined and specific patterns near each TSS, indicating potential epigenetic codes regulating gene transcription.ConclusionsSingle-nucleosome histone modification data render epigenomic data become quantitative, which is essential for investigating dynamic changes of epigenetic regulation in the biological process or for functional epigenomics studies. Thus, NUCLIZE turns current epigenomic mapping studies into genuine functional epigenomics studies with quantitative epigenomic data.

Condensin-Dependent Chromatin Compaction Represses Transcription Globally during Quiescence

Quiescence is a stress-resistant state in which cells reversibly exit the cell cycle and suspend most processes. Quiescence is essential for stem cell maintenance, and its misregulation is implicated in tumor formation. One of the hallmarks of quiescent cells is highly condensed chromatin. Because condensed chromatin often correlates with transcriptional silencing, it has been hypothesized that chromatin compaction represses transcription during quiescence. However, the technology to test this model by determining chromatin structure within cells at gene resolution has not previously been available. Here, we use Micro-C XL to map chromatin contacts at single-nucleosome resolution genome-wide in quiescent Saccharomyces cerevisiae cells. We describe chromatin domains on the order of 10-60 kilobases that, only in quiescent cells, are formed by condensin-mediated loops. Condensin depletion prevents the compaction of chromatin within domains and leads to widespread transcriptional de-repression. Finally, we demonstrate that condensin-dependent chromatin compaction is conserved in quiescent human fibroblasts.

Mapping Combinatorial Epigenetic Modifications at Single Nucleosome Resolution

Mapping Combinatorial Epigenetic Modifications at Single Nucleosome Resolution

Elucidating Combinatorial Chromatin States at Single-Nucleosome Resolution

SummaryChromatin immunoprecipitation followed by sequencing (ChIP-seq) has been instrumental to our current view of chromatin structure and function. It allows genome-wide mapping of histone marks, which demarcate biologically relevant domains. However, ChIP-seq is an ensemble measurement reporting the average occupancy of individual marks in a cell population. Consequently, our understanding of the combinatorial nature of chromatin states relies almost exclusively on correlation between the genomic distributions of individual marks. Here, we report the development of combinatorial-iChIP to determine the genome-wide co-occurrence of histone marks at single-nucleosome resolution. By comparing to a null model, we show that certain combinations of overlapping marks (H3K36me3 and H3K79me3) co-occur more frequently than would be expected by chance, while others (H3K4me3 and H3K36me3) do not, reflecting differences in the underlying chromatin pathways. We further use combinatorial-iChIP to illuminate aspects of the Set2-RPD3S pathway. This approach promises to improve our understanding of the combinatorial complexity of chromatin.

The complex pattern of epigenomic variation between natural yeast strains at single-nucleosome resolution.

BackgroundEpigenomic studies on humans and model species have revealed substantial inter-individual variation in histone modification profiles. However, the pattern of this variation has not been precisely characterized, particularly regarding which genomic features are enriched for variability and whether distinct histone marks co-vary synergistically. Yeast allows us to investigate intra-species variation at high resolution while avoiding other sources of variation, such as cell type or subtype.ResultsWe profiled histone marks H3K4me3, H3K9ac, H3K14ac, H4K12ac and H3K4me1 in three unrelated wild strains of Saccharomyces cerevisiae at single-nucleosome resolution and analyzed inter-strain differences statistically. All five marks varied significantly at specific loci, but to different extents. The number of nucleosomes varying for a given mark between two strains ranged from 20 to several thousands; +1 nucleosomes were significantly less subject to variation. Genes with highly evolvable or responsive expression showed higher variability; however, the variation pattern could not be explained by known transcriptional differences between the strains. Synergistic variation of distinct marks was not systematic, with surprising differences between functionally related H3K9ac and H3K14ac. Interestingly, H3K14ac differences that persisted through transient hyperacetylation were supported by H3K4me3 differences, suggesting stabilization via cross talk.ConclusionsQuantitative variation of histone marks among S. cerevisiae strains is abundant and complex. Its relation to functional characteristics is modular and seems modest, with partial association with gene expression divergences, differences between functionally related marks and partial co-variation between marks that may confer stability. Thus, the specific context of studies, such as which precise marks, individuals and genomic loci are investigated, is primordial in population epigenomics studies. The complexity found in this pilot survey in yeast suggests that high complexity can be anticipated among higher eukaryotes, including humans.Electronic supplementary materialThe online version of this article (doi:10.1186/s13072-015-0019-3) contains supplementary material, which is available to authorized users.

A high-throughput ChIP-Seq for large-scale chromatin studies.

We present a modified approach of chromatin immuno-precipitation followed by sequencing (ChIP-Seq), which relies on the direct ligation of molecular barcodes to chromatin fragments, thereby permitting experimental scale-up. With Bar-ChIP now enabling the concurrent profiling of multiple DNA–protein interactions, we report the simultaneous generation of 90 ChIP-Seq datasets without any robotic instrumentation. We demonstrate that application of Bar-ChIP to a panel of Saccharomyces cerevisiae chromatin-associated mutants provides a rapid and accurate genome-wide overview of their chromatin status. Additionally, we validate the utility of this technology to derive novel biological insights by identifying a role for the Rpd3S complex in maintaining H3K14 hypo-acetylation in gene bodies. We also report an association between the presence of intragenic H3K4 tri-methylation and the emergence of cryptic transcription in a Set2 mutant. Finally, we uncover a crosstalk between H3K14 acetylation and H3K4 methylation in this mutant. These results show that Bar-ChIP enables biological discovery through rapid chromatin profiling at single-nucleosome resolution for various conditions and protein modifications at once.

Natural Single-Nucleosome Epi-Polymorphisms in Yeast

Epigenomes commonly refer to the sequence of presence/absence of specific epigenetic marks along eukaryotic chromatin. Complete histone-borne epigenomes have now been described at single-nucleosome resolution from various organisms, tissues, developmental stages, or diseases, yet their intra-species natural variation has never been investigated. We describe here that the epigenomic sequence of histone H3 acetylation at Lysine 14 (H3K14ac) differs greatly between two unrelated strains of the yeast Saccharomyces cerevisiae. Using single-nucleosome chromatin immunoprecipitation and mapping, we interrogated 58,694 nucleosomes and found that 5,442 of them differed in their level of H3K14 acetylation, at a false discovery rate (FDR) of 0.0001. These Single Nucleosome Epi-Polymorphisms (SNEPs) were enriched at regulatory sites and conserved non-coding DNA sequences. Surprisingly, higher acetylation in one strain did not imply higher expression of the relevant gene. However, SNEPs were enriched in genes of high transcriptional variability and one SNEP was associated with the strength of gene activation upon stimulation. Our observations suggest a high level of inter-individual epigenomic variation in natural populations, with essential questions on the origin of this diversity and its relevance to gene x environment interactions.

Cell Cycle– and Chaperone-Mediated Regulation of H3K56ac Incorporation in Yeast

Acetylation of histone H3 lysine 56 is a covalent modification best known as a mark of newly replicated chromatin, but it has also been linked to replication-independent histone replacement. Here, we measured H3K56ac levels at single-nucleosome resolution in asynchronously growing yeast cultures, as well as in yeast proceeding synchronously through the cell cycle. We developed a quantitative model of H3K56ac kinetics, which shows that H3K56ac is largely explained by the genomic replication timing and the turnover rate of each nucleosome, suggesting that cell cycle profiles of H3K56ac should reveal most first-time nucleosome incorporation events. However, since the deacetylases Hst3/4 prevent use of H3K56ac as a marker for histone deposition during M phase, we also directly measured M phase histone replacement rates. We report a global decrease in turnover rates during M phase and a further specific decrease in turnover at several early origins of replication, which switch from rapidly replaced in G1 phase to stably bound during M phase. Finally, by measuring H3 replacement in yeast deleted for the H3K56 acetyltransferase Rtt109 and its two co-chaperones Asf1 and Vps75, we find evidence that Rtt109 and Asf1 preferentially enhance histone replacement at rapidly replaced nucleosomes, whereas Vps75 appears to inhibit histone turnover at those loci. These results provide a broad perspective on histone replacement/incorporation throughout the cell cycle and suggest that H3K56 acetylation provides a positive-feedback loop by which replacement of a nucleosome enhances subsequent replacement at the same location.

Dynamics of Replication-Independent Histone Turnover in Budding Yeast

Chromatin plays roles in processes governed by different time scales. To assay the dynamic behavior of chromatin in living cells, we used genomic tiling arrays to measure histone H3 turnover in G1-arrested Saccharomyces cerevisiae at single-nucleosome resolution over 4% of the genome, and at lower (approximately 265 base pair) resolution over the entire genome. We find that nucleosomes at promoters are replaced more rapidly than at coding regions and that replacement rates over coding regions correlate with polymerase density. In addition, rapid histone turnover is found at known chromatin boundary elements. These results suggest that rapid histone turnover serves to functionally separate chromatin domains and prevent spread of histone states.

Single-Nucleosome Mapping of Histone Modifications in S. cerevisiae

Covalent modification of histone proteins plays a role in virtually every process on eukaryotic DNA, from transcription to DNA repair. Many different residues can be covalently modified, and it has been suggested that these modifications occur in a great number of independent, meaningful combinations. Published low-resolution microarray studies on the combinatorial complexity of histone modification patterns suffer from confounding effects caused by the averaging of modification levels over multiple nucleosomes. To overcome this problem, we used a high-resolution tiled microarray with single-nucleosome resolution to investigate the occurrence of combinations of 12 histone modifications on thousands of nucleosomes in actively growing S. cerevisiae. We found that histone modifications do not occur independently; there are roughly two groups of co-occurring modifications. One group of lysine acetylations shows a sharply defined domain of two hypo-acetylated nucleosomes, adjacent to the transcriptional start site, whose occurrence does not correlate with transcription levels. The other group consists of modifications occurring in gradients through the coding regions of genes in a pattern associated with transcription. We found no evidence for a deterministic code of many discrete states, but instead we saw blended, continuous patterns that distinguish nucleosomes at one location (e.g., promoter nucleosomes) from those at another location (e.g., over the 3′ ends of coding regions). These results are consistent with the idea of a simple, redundant histone code, in which multiple modifications share the same role.