A theoretical model of Polycomb/Trithorax action unites stable epigenetic memory and dynamic regulation

Chromatin enables precise and scalable gene regulation with factors of limited specificity

Biophysical constraints limit the specificity with which transcription factors (TFs) can target regulatory DNA. While individual nontarget binding events may be low affinity, the sheer number of such interactions could present a challenge for gene regulation by degrading its precision or possibly leading to an erroneous induction state. Chromatin can prevent nontarget binding by rendering DNA physically inaccessible to TFs, at the cost of energy-consuming remodeling orchestrated by pioneer factors (PFs). Under what conditions and by how much can chromatin reduce regulatory errors on a global scale? We use a theoretical approach to compare two scenarios for gene regulation: one that relies on TF binding to free DNA alone and one that uses a combination of TFs and chromatin-regulating PFs to achieve desired gene expression patterns. We find, first, that chromatin effectively silences groups of genes that should be simultaneously OFF, thereby allowing more accurate graded control of expression for the remaining ON genes. Second, chromatin buffers the deleterious consequences of nontarget binding as the number of OFF genes grows, permitting a substantial expansion in regulatory complexity. Third, chromatin-based regulation productively co-opts nontarget TF binding for ON genes in order to establish a "leaky" baseline expression level, which targeted activator or repressor binding subsequently up- or down-modulates. Thus, on a global scale, using chromatin simultaneously alleviates pressure for high specificity of regulatory interactions and enables an increase in genome size with minimal impact on global expression error.

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.

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.

Polycomb repression of Hox genes involves spatial feedback but not domain compaction or phase transition

Polycomb group proteins have a critical role in silencing transcription during development. It is commonly proposed that Polycomb-dependent changes in genome folding, which compact chromatin, contribute directly to repression by blocking the binding of activating complexes. Recently, it has also been argued that liquid-liquid demixing of Polycomb proteins facilitates this compaction and repression by phase-separating target genes into a membraneless compartment. To test these models, we used Optical Reconstruction of Chromatin Architecture to trace the Hoxa gene cluster, a canonical Polycomb target, in thousands of single cells. Across multiple cell types, we find that Polycomb-bound chromatin frequently explores decompact states and partial mixing with neighboring chromatin, while remaining uniformly repressed, challenging the repression-by-compaction or phase-separation models. Using polymer simulations, we show that these observed flexible ensembles can be explained by 'spatial feedback'-transient contacts that contribute to the propagation of the epigenetic state (epigenetic memory), without inducing a globular organization.

Application and prospects of single-cell and spatial omics technologies in woody plants

Over the past decade, high-throughput sequencing and high-resolution single-cell transcriptome sequencing technologies have undergone rapid development, leading to significant breakthroughs. Traditional molecular biology methods are limited in their ability to unravel cellular-level heterogeneity within woody plant tissues. Consequently, techniques such as single-cell transcriptomics, single-cell epigenetics, and spatial transcriptomics are rapidly gaining popularity in the study of woody plants. In this review, we provide a comprehensive overview of the development of these technologies, with a focus on their applications and the challenges they present in single-cell transcriptome research in woody plants. In particular, we delve into the similarities and differences among the results of current studies and analyze the reasons behind these differences. Furthermore, we put forth potential solutions to overcome the challenges encountered in single-cell transcriptome applications in woody plants. Finally, we discuss the application directions of these techniques to address key challenges in woody plant research in the future.

Tau protein modulates an epigenetic mechanism of cellular senescence in human SH-SY5Y neuroblastoma cells

Introduction: Progressive Tau deposition in neurofibrillary tangles and neuropil threads is the hallmark of tauopathies, a disorder group that includes Alzheimer's disease. Since Tau is a microtubule-associated protein, a prevalent concept to explain the pathogenesis of tauopathies is that abnormal Tau modification contributes to dissociation from microtubules, assembly into multimeric β-sheets, proteotoxicity, neuronal dysfunction and cell loss. Tau also localizes in the cell nucleus and evidence supports an emerging function of Tau in DNA stability and epigenetic modulation. Methods: To better characterize the possible role of Tau in regulation of chromatin compaction and subsequent gene expression, we performed a bioinformatics analysis of transcriptome data obtained from Tau-depleted human neuroblastoma cells. Results: Among the transcripts deregulated in a Tau-dependent manner, we found an enrichment of target genes for the polycomb repressive complex 2. We further describe decreased cellular amounts of the core components of the polycomb repressive complex 2 and lower histone 3 trimethylation in Tau deficient cells. Among the de-repressed polycomb repressive complex 2 target gene products, IGFBP3 protein was found to be linked to increased senescence induction in Tau-deficient cells. Discussion: Our findings propose a mechanism for Tau-dependent epigenetic modulation of cell senescence, a key event in pathologic aging.

Dual functions of TET1 in germ layer lineage bifurcation distinguished by genomic context and dependence on 5-methylcytosine oxidation

Gastrulation begins when the epiblast forms the primitive streak or becomes definitive ectoderm. During this lineage bifurcation, the DNA dioxygenase TET1 has bipartite functions in transcriptional activation and repression, but the mechanisms remain unclear. By converting mouse embryonic stem cells (ESCs) into neuroprogenitors, we defined how Tet1-/- cells switch from neuroectoderm fate to form mesoderm and endoderm. We identified the Wnt repressor Tcf7l1 as a TET1 target that suppresses Wnt/β-catenin and Nodal signalling. ESCs expressing catalytic dead TET1 retain neural potential but activate Nodal and subsequently Wnt/β-catenin pathways to generate also mesoderm and endoderm. At CpG-poor distal enhancers, TET1 maintains accessible chromatin at neuroectodermal loci independently of DNA demethylation. At CpG-rich promoters, DNA demethylation by TET1 affects the expression of bivalent genes. In ESCs, a non-catalytic TET1 cooperation with Polycomb represses primitive streak genes; post-lineage priming, the interaction becomes antagonistic at neuronal genes, when TET1's catalytic activity is further involved by repressing Wnt signalling. The convergence of repressive DNA and histone methylation does not inhibit neural induction in Tet1-deficient cells, but some DNA hypermethylated loci persist at genes with brain-specific functions. Our results reveal versatile switching of non-catalytic and catalytic TET1 activities based on genomic context, lineage and developmental stage.

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.

Regulation, functions and transmission of bivalent chromatin during mammalian development

Cells differentiate and progress through development guided by a dynamic chromatin landscape that mediates gene expression programmes. During development, mammalian cells display a paradoxical chromatin state: histone modifications associated with gene activation (trimethylated histone H3 Lys4 (H3K4me3)) and with gene repression (trimethylated H3 Lys27 (H3K27me3)) co-occur at promoters of developmental genes. This bivalent chromatin modification state is thought to poise important regulatory genes for expression or repression during cell-lineage specification. In this Review, we discuss recent work that has expanded our understanding of the molecular basis of bivalent chromatin and its contributions to mammalian development. We describe the factors that establish bivalency, especially histone-lysine N-methyltransferase 2B (KMT2B) and Polycomb repressive complex 2 (PRC2), and consider evidence indicating that PRC1 shapes bivalency and may contribute to its transmission between generations. We posit that bivalency is a key feature of germline and embryonic stem cells, as well as other types of stem and progenitor cells. Finally, we discuss the relevance of bivalent chromtin to human development and cancer, and outline avenues of future research.

Epigenetic factor competition reshapes the EMT landscape

The emergence of and transitions between distinct phenotypes in isogenic cells can be attributed to the intricate interplay of epigenetic marks, external signals, and gene-regulatory elements. These elements include chromatin remodelers, histone modifiers, transcription factors, and regulatory RNAs. Mathematical models known as gene-regulatory networks (GRNs) are an increasingly important tool to unravel the workings of such complex networks. In such models, epigenetic factors are usually proposed to act on the chromatin regions directly involved in the expression of relevant genes. However, it has been well-established that these factors operate globally and compete with each other for targets genome-wide. Therefore, a perturbation of the activity of a regulator can redistribute epigenetic marks across the genome and modulate the levels of competing regulators. In this paper, we propose a conceptual and mathematical modeling framework that incorporates both local and global competition effects between antagonistic epigenetic regulators, in addition to local transcription factors, and show the counterintuitive consequences of such interactions. We apply our approach to recent experimental findings on the epithelial-mesenchymal transition (EMT). We show that it can explain the puzzling experimental data, as well as provide verifiable predictions.

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.

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

In eukaryotes, many stable and heritable phenotypes arise from the same DNA sequence, owing to epigenetic regulatory mechanisms relying on the molecular cooperativity of 'reader-writer' enzymes. In this work, we focus on the fundamental, generic mechanisms behind the epigenome memory encoded by post-translational modifications of histone tails. Based on experimental knowledge, we introduce a unified modeling framework, the painter model, describing the mechanistic interplay between sequence-specific recruitment of chromatin regulators, chromatin-state-specific reader-writer processes and long-range spreading mechanisms. A systematic analysis of the model building blocks highlights the crucial impact of tridimensional chromatin organization and state-specific recruitment of enzymes on the stability of epigenomic domains and on gene expression. In particular, we show that enhanced 3D compaction of the genome and enzyme limitation facilitate the formation of ultra-stable, confined chromatin domains. The model also captures how chromatin state dynamics impact the intrinsic transcriptional properties of the region, slower kinetics leading to noisier expression. We finally apply our framework to analyze experimental data, from the propagation of γH2AX around DNA breaks in human cells to the maintenance of heterochromatin in fission yeast, illustrating how the painter model can be used to extract quantitative information on epigenomic molecular processes.

De Novo Polycomb Recruitment and Repressive Domain Formation

Every cell of an organism shares the same genome; even so, each cellular lineage owns a different transcriptome and proteome. The Polycomb group proteins (PcG) are essential regulators of gene repression patterning during development and homeostasis. However, it is unknown how the repressive complexes, PRC1 and PRC2, identify their targets and elicit new Polycomb domains during cell differentiation. Classical recruitment models consider the pre-existence of repressive histone marks; still, de novo target binding overcomes the absence of both H3K27me3 and H2AK119ub. The CpG islands (CGIs), non-core proteins, and RNA molecules are involved in Polycomb recruitment. Nonetheless, it is unclear how de novo targets are identified depending on the physiological context and developmental stage and which are the leading players stabilizing Polycomb complexes at domain nucleation sites. Here, we examine the features of de novo sites and the accessory elements bridging its recruitment and discuss the first steps of Polycomb domain formation and transcriptional regulation, comprehended by the experimental reconstruction of the repressive domains through time-resolved genomic analyses in mammals.

The Role of Polycomb Proteins in Cell Lineage Commitment and Embryonic Development

Embryonic development is a highly intricate and complex process. Different regulatory mechanisms cooperatively dictate the fate of cells as they progress from pluripotent stem cells to terminally differentiated cell types in tissues. A crucial regulator of these processes is the Polycomb Repressive Complex 2 (PRC2). By catalyzing the mono-, di-, and tri-methylation of lysine residues on histone H3 tails (H3K27me3), PRC2 compacts chromatin by cooperating with Polycomb Repressive Complex 1 (PRC1) and represses transcription of target genes. Proteomic and biochemical studies have revealed two variant complexes of PRC2, namely PRC2.1 which consists of the core proteins (EZH2, SUZ12, EED, and RBBP4/7) interacting with one of the Polycomb-like proteins (MTF2, PHF1, PHF19), and EPOP or PALI1/2, and PRC2.2 which contains JARID2 and AEBP2 proteins. MTF2 and JARID2 have been discovered to have crucial roles in directing and recruiting PRC2 to target genes for repression in embryonic stem cells (ESCs). Following these findings, recent work in the field has begun to explore the roles of different PRC2 variant complexes during different stages of embryonic development, by examining molecular phenotypes of PRC2 mutants in both in vitro (2D and 3D differentiation) and in vivo (knock-out mice) assays, analyzed with modern single-cell omics and biochemical assays. In this review, we discuss the latest findings that uncovered the roles of different PRC2 proteins during cell-fate and lineage specification and extrapolate these findings to define a developmental roadmap for different flavors of PRC2 regulation during mammalian embryonic development.

Chromatin Hubs: A biological and computational outlook

This review discusses our current understanding of chromatin biology and bioinformatics under the unifying concept of “chromatin hubs.” The first part reviews the biology of chromatin hubs, including chromatin–chromatin interaction hubs, chromatin hubs at the nuclear periphery, hubs around macromolecules such as RNA polymerase or lncRNAs, and hubs around nuclear bodies such as the nucleolus or nuclear speckles. The second part reviews existing computational methods, including enhancer–promoter interaction prediction, network analysis, chromatin domain callers, transcription factory predictors, and multi-way interaction analysis. We introduce an integrated model that makes sense of the existing evidence. Understanding chromatin hubs may allow us (i) to explain long-unsolved biological questions such as interaction specificity and redundancy of mechanisms, (ii) to develop more realistic kinetic and functional predictions, and (iii) to explain the etiology of genomic disease.

Investigating Histone Modification Dynamics by Mechanistic Computational Modeling

The maintenance of transcriptional states regulated by histone modifications and controlled switching between these states are fundamental concepts in our understanding of nucleosome-mediated epigenetic memory. Any approach relying on genome-wide bioinformatic analyses alone offers limited scope for dissecting the molecular mechanisms involved in maintenance and switching. Mechanistic mathematical models-describing the dynamics of histone modifications at individual genomic loci-offer an alternative way to investigate these mechanisms. These models, in conjunction with quantitative experimental data-ChIP data, quantification of mRNA levels, and single-cell fluorescence tracking in clonal lineages-can generate predictions that drive more targeted experiments, allowing us to understand mechanisms that would be challenging to unravel by a purely experimental approach. In this chapter, we describe a generic stochastic modeling framework that can be used to capture histone modification dynamics and associated molecular processes-including transcription and read-write feedback by chromatin modifying complexes-at individual genomic loci. Using a specific example-transcriptional silencing by Polycomb-mediated H3K27 methylation-we demonstrate how to construct and simulate a stochastic histone modification model. We provide a step-by-step guide to programming simulations for such a model and discuss how to analyze the simulation output.

The molecular principles of gene regulation by Polycomb repressive complexes

Precise control of gene expression is fundamental to cell function and development. Although ultimately gene expression relies on DNA-binding transcription factors to guide the activity of the transcription machinery to genes, it has also become clear that chromatin and histone post-translational modification have fundamental roles in gene regulation. Polycomb repressive complexes represent a paradigm of chromatin-based gene regulation in animals. The Polycomb repressive system comprises two central protein complexes, Polycomb repressive complex 1 (PRC1) and PRC2, which are essential for normal gene regulation and development. Our early understanding of Polycomb function relied on studies in simple model organisms, but more recently it has become apparent that this system has expanded and diverged in mammals. Detailed studies are now uncovering the molecular mechanisms that enable mammalian PRC1 and PRC2 to identify their target sites in the genome, communicate through feedback mechanisms to create Polycomb chromatin domains and control transcription to regulate gene expression. In this Review, we discuss and contextualize the emerging principles that define how this fascinating chromatin-based system regulates gene expression in mammals.

The Flowering Season-Meter at FLOWERING LOCUS C Across Life Histories in Crucifers

Many plant species overwinter before they flower. Transition to flowering is aligned to the seasonal transition as a response to the prolonged cold in winter by a process called vernalization. Multiple well-documented vernalization properties in crucifer species with diverse life histories are derived from environmental regulation of a central inhibitor of the flowering gene, Flowering Locus C (FLC). Episode(s) of flowering are prevented during high FLC expression and enabled during low FLC expression. FLC repression outlasts the winter to coincide with spring; this heterochronic aspect is termed "winter memory." In the annual Arabidopsis thaliana, winter memory has long been associated with the highly conserved histone modifiers Polycomb and Trithorax, which have antagonistic roles in transcription. However, there are experimental limitations in determining how dynamic, heterogenous histone modifications within the FLC locus generate the final transcriptional output. Recent theoretical considerations on cell-to-cell variability in gene expression and histone modifications generating bistable states brought support to the hypothesis of chromatin-encoded memory, as with other experimental systems in eukaryotes. Furthermore, these advances unify multiple properties of vernalization, not only the winter memory. Similarly, in the perennial Arabidopsis halleri ssp. gemmifera, recent integration of molecular with mathematical and ecological approaches unifies FLC chromatin features with the all-year-round memory of seasonal temperature. We develop the concept of FLC season-meter to combine existing information from the contrasting annual/perennial and experimental/theoretical sectors into a transitional framework. We highlight simplicity, high conservation, and discrete differences across extreme life histories in crucifers.

Using computational modelling to reveal mechanisms of epigenetic Polycomb control

The Polycomb system is essential for stable gene silencing in many organisms. This regulation is achieved in part through addition of the histone modifications H3K27me2/me3 by Polycomb Repressive Complex 2 (PRC2). These modifications are believed to be the causative epigenetic memory elements of PRC2-mediated silencing. As these marks are stored locally in the chromatin, PRC2-based memory is a cis-acting system. A key feature of stable epigenetic memory in cis is PRC2-mediated, self-reinforcing feedback from K27-methylated histones onto nearby histones in a read-write paradigm. However, it was not clear under what conditions such feedback can lead to stable memory, able, for example, to survive the perturbation of histone dilution at DNA replication. In this context, computational modelling has allowed a rigorous exploration of possible underlying memory mechanisms and has also greatly accelerated our understanding of switching between active and silenced states. Specifically, modelling has predicted that switching and memory at Polycomb loci is digital, with a locus being either active or inactive, rather than possessing intermediate, smoothly varying levels of activation. Here, we review recent advances in models of Polycomb control, focusing on models of epigenetic switching through nucleation and spreading of H3K27me2/me3. We also examine models that incorporate transcriptional feedback antagonism and those including bivalent chromatin states. With more quantitative experimental data on histone modification kinetics, as well as single-cell resolution data on transcription and protein levels for PRC2 targets, we anticipate an expanded need for modelling to help dissect increasingly interconnected and complex memory mechanisms.