Painters in chromatin: a unified quantitative framework to systematically characterize epigenome regulation and memory

Role of diffusion and reaction of the constituents in spreading of histone modification marks

Cells switch genes ON or OFF by altering the state of chromatin via histone modifications at specific regulatory locations along the chromatin polymer. These gene regulation processes are carried out by a network of reactions in which the histone marks spread to neighboring regions with the help of enzymes. In the literature, this spreading has been studied as a purely kinetic, non-diffusive process considering the interactions between neighboring nucleosomes. In this work, we go beyond this framework and study the spreading of modifications using a reaction-diffusion (RD) model accounting for the diffusion of the constituents. We quantitatively segregate the modification profiles generated from kinetic and RD models. The diffusion and degradation of enzymes set a natural length scale for limiting the domain size of modification spreading, and the resulting enzyme limitation is inherent in our model. We also demonstrate the emergence of confined modification domains without the explicit requirement of a nucleation site. We explore polymer compaction effects on spreading and show that single-cell domains may differ from averaged profiles. We find that the modification profiles from our model are comparable with existing H3K9me3 data of S. pombe.

Phase separation and inheritance of repressive chromatin domains

Polycomb-associated chromatin and pericentromeric heterochromatin form genomic domains important for the epigenetic regulation of gene expression. Both Polycomb complexes and heterochromatin factors rely on 'read and write' mechanisms, which, on their own, are not sufficient to explain the formation and the maintenance of these epigenetic domains. Microscopy has revealed that they form specific nuclear compartments separated from the rest of the genome. Recently, some subunits of these molecular machineries have been shown to undergo phase separation, both in vitro and in vivo, suggesting that phase separation might play important roles in the formation and the function of these two kinds of repressive chromatin. In this review, we will present the recent advances in the field of facultative and constitutive heterochromatin formation and maintenance through phase separation.

Dynamical phase transition in models that couple chromatin folding with histone modifications

Genomic regions can acquire heritable epigenetic states through unique histone modifications, which lead to stable gene expression patterns without altering the underlying DNA sequence. However, the relationship between chromatin conformational dynamics and epigenetic stability is poorly understood. In this paper, we propose kinetic models to investigate the dynamic fluctuations of histone modifications and the spatial interactions between nucleosomes. Our model explicitly incorporates the influence of chemical modifications on the structural stability of chromatin and the contribution of chromatin contacts to the cooperative nature of chemical reactions. Through stochastic simulations and analytical theory, we have discovered distinct steady-state outcomes in different kinetic regimes, resembling a dynamical phase transition. Importantly, we have validated that the emergence of this transition, which occurs on biologically relevant timescales, is robust against variations in model design and parameters. Our findings suggest that the viscoelastic properties of chromatin and the timescale at which it transitions from a gel-like to a liquidlike state significantly impact dynamic processes that occur along the one-dimensional DNA sequence.

Mathematical model for the role of multiple pericentromeric repeats on heterochromatin assembly

Although the length and constituting sequences for pericentromeric repeats are highly variable across eukaryotes, the presence of multiple pericentromeric repeats is one of the conserved features of the eukaryotic chromosomes. Pericentromeric heterochromatin is often misregulated in human diseases, with the expansion of pericentromeric repeats in human solid cancers. In this article, we have developed a mathematical model of the RNAi-dependent methylation of H3K9 in the pericentromeric region of fission yeast. Our model, which takes copy number as an explicit parameter, predicts that the pericentromere is silenced only if there are many copies of repeats. It becomes bistable or desilenced if the copy number of repeats is reduced. This suggests that the copy number of pericentromeric repeats alone can determine the fate of heterochromatin silencing in fission yeast. Through sensitivity analysis, we identified parameters that favor bistability and desilencing. Stochastic simulation shows that faster cell division and noise favor the desilenced state. These results show the unexpected role of pericentromeric repeat copy number in gene silencing and provide a quantitative basis for how the copy number allows or protects repetitive and unique parts of the genome from heterochromatin silencing, respectively.

Interplay between epigenome and 3D chromatin structure

Epigenetic mechanisms, primarily mediated through histone and DNA modifications, play a pivotal role in orchestrating the functional identity of a cell and its response to environmental cues. Similarly, the spatial arrangement of chromatin within the threedimensional (3D) nucleus has been recognized as a significant factor influencing genomic function. Investigating the relationship between epigenetic regulation and 3D chromatin structure has revealed correlation and causality between these processes, from the global alignment of average chromatin structure with chromatin marks to the nuanced correlations at smaller scales. This review aims to dissect the biological significance and the interplay between the epigenome and 3D chromatin structure, while also exploring the underlying molecular mechanisms. By synthesizing insights from both experimental and modeling perspectives, we seek to provide a comprehensive understanding of cellular functions. [BMB Reports 2023; 56(12): 633-644].

Uncoupling the distinct functions of HP1 proteins during heterochromatin establishment and maintenance

H3K9 methylation (H3K9me) marks transcriptionally silent genomic regions called heterochromatin. HP1 proteins are required to establish and maintain heterochromatin. HP1 proteins bind to H3K9me, recruit factors that promote heterochromatin formation, and oligomerize to form phase-separated condensates. We do not understand how these different HP1 properties are involved in establishing and maintaining transcriptional silencing. Here, we demonstrate that the S. pombe HP1 homolog, Swi6, can be completely bypassed to establish silencing at ectopic and endogenous loci when an H3K4 methyltransferase, Set1, and an H3K14 acetyltransferase, Mst2, are deleted. Deleting Set1 and Mst2 enhances Clr4 enzymatic activity, leading to higher H3K9me levels and spreading. In contrast, Swi6 and its capacity to oligomerize were indispensable during epigenetic maintenance. Our results demonstrate the role of HP1 proteins in regulating histone modification crosstalk during establishment and identify a genetically separable function in maintaining epigenetic memory.

Design principles of 3D epigenetic memory systems

Cells remember their identities, in part, by using epigenetic marks-chemical modifications placed along the genome. How can mark patterns remain stable over cell generations despite their constant erosion by replication and other processes? We developed a theoretical model that reveals that three-dimensional (3D) genome organization can stabilize epigenetic memory as long as (i) there is a large density difference between chromatin compartments, (ii) modifying "reader-writer" enzymes spread marks in three dimensions, and (iii) the enzymes are limited in abundance relative to their histone substrates. Analogous to an associative memory that encodes memory in neuronal connectivity, mark patterns are encoded in a 3D network of chromosomal contacts. Our model provides a unified account of diverse observations and reveals a key role of 3D genome organization in epigenetic memory.

Regulation of the heterochromatin spreading reaction by trans-acting factors

Heterochromatin is a gene-repressive protein-nucleic acid ultrastructure that is initially nucleated by DNA sequences. However, following nucleation, heterochromatin can then propagate along the chromatin template in a sequence-independent manner in a reaction termed spreading. At the heart of this process are enzymes that deposit chemical information on chromatin, which attracts the factors that execute chromatin compaction and transcriptional or co/post-transcriptional gene silencing. Given that these enzymes deposit guiding chemical information on chromatin they are commonly termed 'writers'. While the processes of nucleation and central actions of writers have been extensively studied and reviewed, less is understood about how the spreading process is regulated. We discuss how the chromatin substrate is prepared for heterochromatic spreading, and how trans-acting factors beyond writer enzymes regulate it. We examine mechanisms by which trans-acting factors in Suv39, PRC2, SETDB1 and SIR writer systems regulate spreading of the respective heterochromatic marks across chromatin. While these systems are in some cases evolutionarily and mechanistically quite distant, common mechanisms emerge which these trans-acting factors exploit to tune the spreading reaction.

Uncoupling the distinct functions of HP1 proteins during heterochromatin establishment and maintenance

H3K9 methylation (H3K9me) marks transcriptionally silent genomic regions called heterochromatin. HP1 proteins are required to establish and maintain heterochromatin. HP1 proteins bind to H3K9me, recruit factors that promote heterochromatin formation, and oligomerize to form phase-separated condensates. We do not understand how HP1 protein binding to heterochromatin establishes and maintains transcriptional silencing. Here, we demonstrate that the S.pombe HP1 homolog, Swi6, can be completely bypassed to establish silencing at ectopic and endogenous loci when an H3K4 methyltransferase, Set1 and an H3K14 acetyltransferase, Mst2 are deleted. Deleting Set1 and Mst2 enhances Clr4 enzymatic activity, leading to higher H3K9me levels and spreading. In contrast, Swi6 and its capacity to oligomerize were indispensable during epigenetic maintenance. Our results demonstrate the role of HP1 proteins in regulating histone modification crosstalk during establishment and identifies a genetically separable function in maintaining epigenetic memory.

4D epigenomics: deciphering the coupling between genome folding and epigenomic regulation with biophysical modeling

Recent experimental observations suggest a strong coupling between the 3D nuclear chromosome organization and epigenomics. However, the mechanistic and functional bases of such interplay remain elusive. In this review, we describe how biophysical modeling has been instrumental in characterizing how genome folding may impact the formation of epigenomic domains and, conversely, how epigenomic marks may affect chromosome conformation. Finally, we discuss how this mutual feedback loop between chromatin organization and epigenome regulation, via the formation of physicochemical nanoreactors, may represent a key functional role of 3D compartmentalization in the assembly and maintenance of stable - but yet plastic - epigenomic landscapes.

Dynamical modeling of the H3K27 epigenetic landscape in mouse embryonic stem cells

The Polycomb system via the methylation of the lysine 27 of histone H3 (H3K27) plays central roles in the silencing of many lineage-specific genes during development. Recent experimental evidence suggested that the recruitment of histone modifying enzymes like the Polycomb repressive complex 2 (PRC2) at specific sites and their spreading capacities from these sites are key to the establishment and maintenance of a proper epigenomic landscape around Polycomb-target genes. Here, to test whether such mechanisms, as a minimal set of qualitative rules, are quantitatively compatible with data, we developed a mathematical model that can predict the locus-specific distributions of H3K27 modifications based on previous biochemical knowledge. Within the biological context of mouse embryonic stem cells, our model showed quantitative agreement with experimental profiles of H3K27 acetylation and methylation around Polycomb-target genes in wild-type and mutants. In particular, we demonstrated the key role of the reader-writer module of PRC2 and of the competition between the binding of activating and repressing enzymes in shaping the H3K27 landscape around transcriptional start sites. The predicted dynamics of establishment and maintenance of the repressive trimethylated H3K27 state suggest a slow accumulation, in perfect agreement with experiments. Our approach represents a first step towards a quantitative description of PcG regulation in various cellular contexts and provides a generic framework to better characterize epigenetic regulation in normal or disease situations.