Establishment of heterochromatin in domain-size-dependent bursts

Theory of epigenetic switching due to stochastic histone mark loss during DNA replication

How much information does a cell inherit from its ancestors beyond its genetic sequence? What are the epigenetic mechanisms that allow this? Despite the rise in available epigenetic data, how such information is inherited through the cell cycle is still not fully understood. Often, epigenetic marks can display bistable behaviour and their bistable state is transmitted to daughter cells through the cell cycle, providing the cell with a form of memory. However, loss-of-memory events also take place, where a daughter cell switches epigenetic state (with respect to the mother cell). Here, we develop a framework to compute these epigenetic switching rates, for the case when they are driven by DNA replication, i.e. the frequency of loss-of-memory events due to replication. We consider the dynamics of histone modifications during the cell cycle deterministically, except at DNA replication, where nucleosomes are randomly distributed between the two daughter DNA strands, which is therefore implemented stochastically. This hybrid stochastic-deterministic approach enables an analytic derivation of the replication-driven switching rate. While retaining great simplicity, this framework can explain experimental switching rate data, establishing its biological importance as a framework to quantitatively study epigenetic inheritance.

The fork protection complex promotes parental histone recycling and epigenetic memory

The inheritance of parental histones across the replication fork is thought to mediate epigenetic memory. Here, we reveal that fission yeast Mrc1 (CLASPIN in humans) binds H3-H4 tetramers and operates as a central coordinator of symmetric parental histone inheritance. Mrc1 mutants in a key connector domain disrupted segregation of parental histones to the lagging strand comparable to Mcm2 histone-binding mutants. Both mutants showed clonal and asymmetric loss of H3K9me-mediated gene silencing. AlphaFold predicted co-chaperoning of H3-H4 tetramers by Mrc1 and Mcm2, with the Mrc1 connector domain bridging histone and Mcm2 binding. Biochemical and functional analysis validated this model and revealed a duality in Mrc1 function: disabling histone binding in the connector domain disrupted lagging-strand recycling while another histone-binding mutation impaired leading strand recycling. We propose that Mrc1 toggles histones between the lagging and leading strand recycling pathways, in part by intra-replisome co-chaperoning, to ensure epigenetic transmission to both daughter cells.

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.

Dissecting Mechanisms of Epigenetic Memory Through Computational Modeling

Understanding the mechanistic basis of epigenetic memory has proven to be a difficult task due to the underlying complexity of the systems involved in its establishment and maintenance. Here, we review the role of computational modeling in helping to unlock this complexity, allowing the dissection of intricate feedback dynamics. We focus on three forms of epigenetic memory encoded in gene regulatory networks, DNA methylation, and histone modifications and discuss the important advantages offered by plant systems in their dissection. We summarize the main modeling approaches involved and highlight the principal conceptual advances that the modeling has enabled through iterative cycles of predictive modeling and experiments. Lastly, we discuss remaining gaps in our understanding and how intertwined theory and experimental approaches might help in their resolution.

A zinc-finger protein Moc3 functions as a transcription activator to promote RNAi-dependent constitutive heterochromatin establishment in fission yeast

In fission yeast, Schizosaccharomyces pombe, constitutive heterochromatin defined by methylation of histone H3 lysine 9 (H3K9me) and its binding protein Swi6/HP1 localizes at the telomere, centromere, and mating-type loci. These loci contain DNA sequences called dg and dh, and the RNA interference (RNAi)-dependent system establishes and maintains heterochromatin at dg/dh. Bi-directional transcription at dg/dh induced by RNA polymerase II is critical in RNAi-dependent heterochromatin formation because the transcribed RNAs provide substrates for siRNA synthesis and a platform for assembling RNAi factors. However, a regulator of dg/dh transcription during the establishment of heterochromatin is not known. Here, we found that a zinc-finger protein Moc3 localizes dh and activates dh-forward transcription in its zinc-finger-dependent manner when heterochromatin structure or heterochromatin-dependent silencing is compromised. However, Moc3 does not localize at normal heterochromatin and does not activate the dh-forward transcription. Notably, the loss of Moc3 caused a retarded heterochromatin establishment, showing that Moc3-dependent dh-forward transcription is critical for RNAi-dependent heterochromatin establishment. Therefore, Moc3 is a transcriptional activator that induces RNAi to establish heterochromatin.

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.

Two-way feedback between chromatin compaction and histone modification state explains Saccharomyces cerevisiae heterochromatin bistability

Compact chromatin is closely linked with gene silencing in part by sterically masking access to promoters, inhibiting transcription factor binding and preventing polymerase from efficiently transcribing a gene. However, a broader hypothesis suggests that chromatin compaction can be both a cause and a consequence of the locus histone modification state, with a tight bidirectional interaction underpinning bistable transcriptional states. To rigorously test this hypothesis, we developed a mathematical model for the dynamics of the HMR locus in Saccharomyces cerevisiae, that incorporates activating histone modifications, silencing proteins, and a dynamic, acetylation-dependent, three-dimensional locus size. Chromatin compaction enhances silencer protein binding, which in turn feeds back to remove activating histone modifications, leading to further compaction. The bistable output of the model was in good agreement with prior quantitative data, including switching rates from expressed to silent states (and vice versa), and protein binding/histone modification levels within the locus. We then tested the model by predicting changes in switching rates as the genetic length of the locus was increased, which were then experimentally verified. Such bidirectional feedback between chromatin compaction and the histone modification state may be a widespread and important regulatory mechanism given the hallmarks of many heterochromatic regions: physical chromatin compaction and dimerizing (or multivalent) silencing proteins.

Safeguarding the epigenome through the cell cycle: a multitasking game

Sustaining cell identity and function across cell division is germane to human development, healthspan, and cancer avoidance. This relies significantly on propagation of chromatin organization between cell generations, as chromatin presents a barrier to cell fate and cell state conversions. Inheritance of chromatin states across the many cell divisions required for development and tissue homeostasis represents a major challenge, especially because chromatin is disrupted to allow passage of the DNA replication fork to synthesize the two daughter strands. This process also leads to a twofold dilution of epigenetic information in histones, which needs to be accurately restored for faithful propagation of chromatin states across cell divisions. Recent research has identified distinct multilayered mechanisms acting to propagate epigenetic information to daughter strands. Here, we summarize key principles of how epigenetic information in parental histones is transferred across DNA replication and how new histones robustly acquire the same information postreplication, representing a core component of epigenetic cell memory.

Diffusion controls local versus dispersed inheritance of histones during replication and shapes epigenomic architecture

The dynamics of inheritance of histones and their associated modifications across cell divisions can have major consequences on maintenance of the cellular epigenomic state. Recent experiments contradict the long-held notion that histone inheritance during replication is always local, suggesting that active and repressed regions of the genome exhibit fundamentally different histone dynamics independent of transcription-coupled turnover. Here we develop a stochastic model of histone dynamics at the replication fork and demonstrate that differential diffusivity of histones in active versus repressed chromatin is sufficient to quantitatively explain these recent experiments. Further, we use the model to predict patterns in histone mark similarity between pairs of genomic loci that should be developed as a result of diffusion, but cannot originate from either PRC2 mediated mark spreading or transcriptional processes. Interestingly, using a combination of CHIP-seq, replication timing and Hi-C datasets we demonstrate that all the computationally predicted patterns are consistently observed for both active and repressive histone marks in two different cell lines. While direct evidence for histone diffusion remains controversial, our results suggest that dislodged histones in euchromatin and facultative heterochromatin may exhibit some level of diffusion within "Diffusion-Accessible-Domains" (DADs), leading to redistribution of epigenetic marks within and across chromosomes. Preservation of the epigenomic state across cell divisions therefore might be achieved not by passing on strict positional information of histone marks, but by maintaining the marks in somewhat larger DADs of the genome.

Two-way feedback between chromatin compaction and histone modification state explains S. cerevisiae heterochromatin bistability

Compact chromatin is closely linked with gene silencing in part by sterically masking access to promoters, inhibiting transcription factor binding and preventing polymerase from efficiently transcribing a gene. Here, we propose a broader view: chromatin compaction can be both a cause and a consequence of the histone modification state, and this tight bidirectional interaction can underpin bistable transcriptional states. To test this theory, we developed a mathematical model for the dynamics of the HMR locus in S. cerevisiae, that incorporates activating histone modifications, silencing proteins and a dynamic, acetylation-dependent, three-dimensional locus size. Chromatin compaction enhances silencer protein binding, which in turn feeds back to remove activating histone modifications, leading to further compaction. The bistable output of the model was in good agreement with prior quantitative data, including switching rates from expressed to silent states, and vice versa, and protein binding levels within the locus. We then tested the model by predicting changes in switching rates as the genetic length of the locus was increased, which were then experimentally verified. This bidirectional feedback between chromatin compaction and the histone modification state may be an important regulatory mechanism at many loci.

A dual, catalytic role for the fission yeast Ccr4-Not complex in gene silencing and heterochromatin spreading

Heterochromatic gene silencing relies on combinatorial control by specific histone modifications, the occurrence of transcription, and/or RNA degradation. Once nucleated, heterochromatin propagates within defined chromosomal regions and is maintained throughout cell divisions to warrant proper genome expression and integrity. In the fission yeast Schizosaccharomyces pombe, the Ccr4-Not complex partakes in gene silencing, but its relative contribution to distinct heterochromatin domains and its role in nucleation versus spreading have remained elusive. Here, we unveil major functions for Ccr4-Not in silencing and heterochromatin spreading at the mating type locus and subtelomeres. Mutations of the catalytic subunits Caf1 or Mot2, involved in RNA deadenylation and protein ubiquitinylation, respectively, result in impaired propagation of H3K9me3 and massive accumulation of nucleation-distal heterochromatic transcripts. Both silencing and spreading defects are suppressed upon disruption of the heterochromatin antagonizing factor Epe1. Overall, our results position the Ccr4-Not complex as a critical, dual regulator of heterochromatic gene silencing and spreading.

Epigenetic switching with asymmetric bridging interactions

Gene expression states are often stably sustained in cis despite massively disruptive events like DNA replication. This is achieved by on-going enzymatic activity that maintains parts of the DNA in either heterochromatic (packed) or euchromatic (free) states, each of which is stabilized by both positive feedback and bridging interactions between individual nucleosomes. In contrast to condensed matter, however, the dynamics is not only governed by equilibrium binding interactions but is also mediated by enzymes that recognize and act on specific amino acid tails of the nucleosomes. The mechanical result is that some nucleosomes can bind to one another and form tightly packed polymer configurations, whereas others remain unbound and form free, noncompact polymer configurations. Here, we study the consequences of such an asymmetric interaction pattern on the dynamics of epigenetic switching. We develop a 3D polymer model and show that traits associated with epigenetic switching, such as bistability and epigenetic memory, are permitted by such a model. We find, however, that the experimentally observed burst-like nature of some epigenetic switches is difficult to reproduce by this biologically motivated interaction. Instead, the behavior seen in experiments can be explained by introducing partial confinement, which particularly affects the euchromatic regions of the chromosome.

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.

A simulation model of heterochromatin formation at submolecular detail

Heterochromatin is a physical state of the chromatin fiber that maintains gene repression during cell development. Although evidence exists on molecular mechanisms involved in heterochromatin formation, a detailed structural mechanism of heterochromatin formation needs a better understanding. We made use of a simple Monte Carlo simulation model with explicit representation of key molecular events to observe molecular self-organization leading to heterochromatin formation. Our simulations provide a structural interpretation of several important traits of the heterochromatinization process. In particular, this study provides a depiction of how small amounts of HP1 are able to induce a highly condensed chromatin state through HP1 dimerization and bridging of sequence-remote nucleosomes. It also elucidates structural roots of a yet poorly understood phenomenon of a nondeterministic nature of heterochromatin formation and subsequent gene repression. Experimental chromatin in vivo assay provides an unbiased estimate of time scale of repressive response to a heterochromatin-triggering event.

The transcription factor Atf1 lowers the transition barrier for nucleosome-mediated establishment of heterochromatin

Transcription factors can exert opposite effects depending on the chromosomal context. The fission yeast transcription factor Atf1 both activates numerous genes in response to stresses and mediates heterochromatic gene silencing in the mating-type region. Investigating this context dependency, we report here that the establishment of silent heterochromatin in the mating-type region occurs at a reduced rate in the absence of Atf1 binding. Quantitative modeling accounts for the observed establishment profiles by a combinatorial recruitment of histone-modifying enzymes: locally by Atf1 at two binding sites and over the whole region by dynamically appearing heterochromatic nucleosomes, a source of which is the RNAi-dependent cenH element. In the absence of Atf1 binding, the synergy is lost, resulting in a slow rate of heterochromatin formation. The system shows how DNA-binding proteins can influence local nucleosome states and thereby potentiate long-range positive feedback on histone-modification reactions to enable heterochromatin formation over large regions in a context-dependent manner.

A Liquid State Perspective on Dynamics of Chromatin Compartments

The interior of the eukaryotic cell nucleus has a crowded and heterogeneous environment packed with chromatin polymers, regulatory proteins, and RNA molecules. Chromatin polymer, assisted by epigenetic modifications, protein and RNA binders, forms multi-scale compartments which help regulate genes in response to cellular signals. Furthermore, chromatin compartments are dynamic and tend to evolve in size and composition in ways that are not fully understood. The latest super-resolution imaging experiments have revealed a much more dynamic and stochastic nature of chromatin compartments than was appreciated before. An emerging mechanism explaining chromatin compartmentalization dynamics is the phase separation of protein and nucleic acids into membraneless liquid condensates. Consequently, concepts and ideas from soft matter and polymer systems have been rapidly entering the lexicon of cell biology. In this respect, the role of computational models is crucial for establishing a rigorous and quantitative foundation for the new concepts and disentangling the complex interplay of forces that contribute to the emergent patterns of chromatin dynamics and organization. Several multi-scale models have emerged to address various aspects of chromatin dynamics, ranging from equilibrium polymer simulations, hybrid non-equilibrium simulations coupling protein binding and chromatin folding, and mesoscopic field-theoretic models. Here, we review these emerging theoretical paradigms and computational models with a particular focus on chromatin’s phase separation and liquid-like properties as a basis for nuclear organization and dynamics.