Nucleation and spreading of a heterochromatic domain in fission yeast

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.

Heterochromatin protein ERH represses alternative cell fates during early mammalian differentiation

During development, H3K9me3 heterochromatin is dynamically rearranged, silencing repeat elements and protein coding genes to restrict cell identity. Enhancer of Rudimentary Homolog (ERH) is an evolutionarily conserved protein originally characterized in fission yeast and recently shown to be required for H3K9me3 maintenance in human fibroblasts, but its function during development remains unknown. Here, we show that ERH is required for proper segregation of the inner cell mass and trophectoderm cell lineages during mouse development by repressing totipotent and alternative lineage programs. During human embryonic stem cell (hESC) differentiation into germ layer lineages, ERH is crucial for silencing naïve and pluripotency genes, transposable elements, and alternative lineage genes. Strikingly, ERH depletion in somatic cells reverts the H3K9me3 landscape to an hESC state and enables naïve and pluripotency gene and transposable element activation during iPSC reprogramming. Our findings reveal a role for ERH in initiation and maintenance of developmentally established gene repression.

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.

HP1-driven phase separation recapitulates the thermodynamics and kinetics of heterochromatin condensate formation

The spatial segregation of pericentromeric heterochromatin (PCH) into distinct, membrane-less nuclear compartments involves the binding of Heterochromatin Protein 1 (HP1) to H3K9me2/3-rich genomic regions. While HP1 exhibits liquid-liquid phase separation properties in vitro, its mechanistic impact on the structure and dynamics of PCH condensate formation in vivo remains largely unresolved. Here, using a minimal theoretical framework, we systematically investigate the mutual coupling between self-interacting HP1-like molecules and the chromatin polymer. We reveal that the specific affinity of HP1 for H3K9me2/3 loci facilitates coacervation in nucleo and promotes the formation of stable PCH condensates at HP1 levels far below the concentration required to observe phase separation in purified protein assays in vitro. These heterotypic HP1-chromatin interactions give rise to a strong dependence of the nucleoplasmic HP1 density on HP1-H3K9me2/3 stoichiometry, consistent with the thermodynamics of multicomponent phase separation. The dynamical cross talk between HP1 and the viscoelastic chromatin scaffold also leads to anomalously slow equilibration kinetics, which strongly depend on the genomic distribution of H3K9me2/3 domains and result in the coexistence of multiple long-lived, microphase-separated PCH compartments. The morphology of these complex coacervates is further found to be governed by the dynamic establishment of the underlying H3K9me2/3 landscape, which may drive their increasingly abnormal, aspherical shapes during cell development. These findings compare favorably to 4D microscopy measurements of HP1 condensate formation in live Drosophila embryos and suggest a general quantitative model of PCH formation based on the interplay between HP1-based phase separation and chromatin polymer mechanics.

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.

The molecular basis of heterochromatin assembly and epigenetic inheritance

Heterochromatin plays a fundamental role in gene regulation, genome integrity, and silencing of repetitive DNA elements. Histone modifications are essential for the establishment of heterochromatin domains, which is initiated by the recruitment of histone-modifying enzymes to nucleation sites. This leads to the deposition of histone H3 lysine-9 methylation (H3K9me), which provides the foundation for building high-concentration territories of heterochromatin proteins and the spread of heterochromatin across extended domains. Moreover, heterochromatin can be epigenetically inherited during cell division in a self-templating manner. This involves a "read-write" mechanism where pre-existing modified histones, such as tri-methylated H3K9 (H3K9me3), support chromatin association of the histone methyltransferase to promote further deposition of H3K9me. Recent studies suggest that a critical density of H3K9me3 and its associated factors is necessary for the propagation of heterochromatin domains across multiple generations. In this review, I discuss the key experiments that have highlighted the importance of modified histones for epigenetic inheritance.

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.

The Fission Yeast Mating-Type Switching Motto: "One-for-Two" and "Two-for-One"

Schizosaccharomyces pombe is an ascomycete fungus that divides by medial fission; it is thus commonly referred to as fission yeast, as opposed to the distantly related budding yeast Saccharomyces cerevisiae. The reproductive lifestyle of S. pombe relies on an efficient genetic sex determination system generating a 1:1 sex ratio and using alternating haploid/diploid phases in response to environmental conditions. In this review, we address how one haploid cell manages to generate two sister cells with opposite mating types, a prerequisite to conjugation and meiosis. This mating-type switching process depends on two highly efficient consecutive asymmetric cell divisions that rely on DNA replication, repair, and recombination as well as the structure and components of heterochromatin. We pay special attention to the intimate interplay between the genetic and epigenetic partners involved in this process to underscore the importance of basic research and its profound implication for a better understanding of chromatin biology.

Histone deacetylation primes self-propagation of heterochromatin domains to promote epigenetic inheritance

Heterochromatin assembly, involving histone H3 lysine-9 methylation (H3K9me), is nucleated at specific genomic sites but can self-propagate across extended domains and, indeed, generations. Self-propagation requires Clr4/Suv39h methyltransferase recruitment by pre-existing H3K9 tri-methylation (H3K9me3) to perpetuate H3K9me deposition and is dramatically affected by chromatin context. However, the mechanism priming self-propagation of heterochromatin remains undefined. We show that robust chromatin association of fission yeast class II histone deacetylase Clr3 is necessary and sufficient to support heterochromatin propagation in different chromosomal contexts. Efficient targeting of Clr3, which suppresses histone turnover and maintains H3K9me3, enables self-propagation of an ectopic heterochromatin domain via the Clr4/Suv39h read-write mechanism requiring methylated histones. The deacetylase activity of Clr3 is necessary and, when inactivated, heterochromatin propagation can be recapitulated by removing two major histone acetyltransferases. Our results show that histone deacetylation, a conserved heterochromatin feature, preserves H3K9me3 that transmits epigenetic memory for stable propagation of silenced chromatin domains through multiple generations.

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.

Nucleosome sliding can influence the spreading of histone modifications

Nucleosomes are the fundamental building blocks of chromatin that not only help in the folding of chromatin, but also in carrying epigenetic information. It is known that nucleosome sliding is responsible for dynamically organizing chromatin structure and the resulting gene regulation. Since sliding can move two neighboring nucleosomes physically close or away, can it play a role in the spreading of histone modifications? We investigate this by simulating a stochastic model that couples nucleosome dynamics with the kinetics of histone modifications. We show that the sliding of nucleosomes can affect the modification pattern as well as the time it takes to modify a given region of chromatin. Exploring different nucleosome densities and modification kinetic parameters, we show that nucleosome sliding can be important for creating histone modification domains. Our model predicts that nucleosome density coupled with sliding dynamics can create an asymmetric histone modification profile around regulatory regions. We also compute the probability distribution of modified nucleosomes and relaxation kinetics of modifications. Our predictions are comparable with known experimental results.

Chromatin dynamics controls epigenetic domain formation

In multicellular organisms, nucleosomes carry epigenetic information that defines distinct patterns of gene expression, which are inherited over multiple generations. The enhanced capacity for information storage arises by nucleosome modifications, which are triggered by enzymes. Modified nucleosomes can transfer the mark to others that are in proximity by a positive-feedback (modification begets modification) mechanism. We created a generic polymer model, referred to as 3DSpreader, in which each bead, representing a nucleosome, stochastically switches between unmodified (U) and modified (M) states depending on the states of the neighbors. Modification begins at a specific nucleation site (NS) that is permanently in the M state, and could spread to other loci that is dictated by chromatin dynamics. Transfer of marks among the non-nucleation loci occurs stochastically as chromatin evolves in time. If the spreading rate is slower than the chromatin relaxation rate, which is biologically pertinent, then finite-sized domains form, driven by contacts between nucleosomes through a three-dimensional looping mechanism. Surprisingly, simulations based on the 3DSpreader model result in finite bounded domains that arise without the need for any boundary elements. Maintenance of spatially and temporally stable domains requires the presence of the NS, whose removal eliminates finite-sized modified domains. The theoretical predictions are in excellent agreement with experimental data for H3K9me3 spreading in mouse embryonic stem cells.

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.

DNA sequence-dependent formation of heterochromatin nanodomains

The mammalian epigenome contains thousands of heterochromatin nanodomains (HNDs) marked by di- and trimethylation of histone H3 at lysine 9 (H3K9me2/3), which have a typical size of 3–10 nucleosomes. However, what governs HND location and extension is only partly understood. Here, we address this issue by introducing the chromatin hierarchical lattice framework (ChromHL) that predicts chromatin state patterns with single-nucleotide resolution. ChromHL is applied to analyse four HND types in mouse embryonic stem cells that are defined by histone methylases SUV39H1/2 or GLP, transcription factor ADNP or chromatin remodeller ATRX. We find that HND patterns can be computed from PAX3/9, ADNP and LINE1 sequence motifs as nucleation sites and boundaries that are determined by DNA sequence (e.g. CTCF binding sites), cooperative interactions between nucleosomes as well as nucleosome-HP1 interactions. Thus, ChromHL rationalizes how patterns of H3K9me2/3 are established and changed via the activity of protein factors in processes like cell differentiation.

The ability to predict epigenetic regulation is an important challenge in biology. Here the authors describe heterochromatin nanodomains (HNDs) and compare four different types of H3K9me2/3-marked HNDs in mouse embryonic stem cells. They further develop a computational framework to predict genome-wide HND maps from DNA sequence and protein concentrations, at single-nucleotide resolution.

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.

Spreading and epigenetic inheritance of heterochromatin require a critical density of histone H3 lysine 9 tri-methylation

In multicellular organisms, a single genome gives rise to a multitude of cell types by enforcing appropriate gene expression patterns. Epigenetic mechanisms involving modification of histones, including tri-methylation of histone H3 lysine 9 (H3K9me3), assemble and propagate repressive heterochromatin to prevent untimely gene expression. Dysregulation of epigenetic gene-silencing mechanisms is a hallmark of a variety of diseases including cancer. However, the requirements for epigenetic transmission of heterochromatin are not well understood. This study reveals the mechanism by which methylated histones provide an epigenetic template for heterochromatin propagation. We discover that a critical threshold of H3K9me3 is required for effective chromatin-association of the histone methyltransferase, which binds to and catalyzes additional H3K9me to propagate heterochromatin and enforce stable gene silencing.

Heterochromatin assembly requires methylation of histone H3 lysine 9 (H3K9me) and serves as a paradigm for understanding the importance of histone modifications in epigenetic genome control. Heterochromatin is nucleated at specific genomic sites and spreads across extended chromosomal domains to promote gene silencing. Moreover, heterochromatic structures can be epigenetically inherited in a self-templating manner, which is critical for stable gene repression. The spreading and inheritance of heterochromatin are believed to be dependent on preexisting H3K9 tri-methylation (H3K9me3), which is recognized by the histone methyltransferase Clr4/Suv39h via its chromodomain, to promote further deposition of H3K9me. However, the process involving the coupling of the “read” and “write” capabilities of histone methyltransferases is poorly understood. From an unbiased genetic screen, we characterize a dominant-negative mutation in histone H3 (H3^G13D) that impairs the propagation of endogenous and ectopic heterochromatin domains in the fission yeast genome. H3^G13D blocks methylation of H3K9 by the Clr4/Suv39h methyltransferase and acts in a dosage-dependent manner to interfere with the spreading and maintenance of heterochromatin. Our analyses show that the incorporation of unmethylatable histone H3^G13D into chromatin decreases H3K9me3 density and thereby compromises the read-write capability of Clr4/Suv39h. Consistently, enhancing the affinity of Clr4/Suv39h for methylated H3K9 is sufficient to overcome the defects in heterochromatin assembly caused by H3^G13D. Our work directly implicates methylated histones in the transmission of epigenetic memory and shows that a critical density threshold of H3K9me3 is required to promote epigenetic inheritance of heterochromatin through the read-write mechanism.

Establishment of heterochromatin in domain-size-dependent bursts

Methylation of histone H3K9 is a hallmark of epigenetic silencing in eukaryotes. Nucleosome modifications often rely on positive feedback where enzymes are recruited by modified nucleosomes. A combination of local and global feedbacks has been proposed to account for some dynamic properties of heterochromatin, but the range at which the global feedbacks operate and the exact mode of heterochromatin propagation are not known. We investigated these questions in fission yeast. Guided by mathematical modeling, we incrementally increased the size of the mating-type region and profiled heterochromatin establishment over time. We observed exponential decays in the proportion of cells with active reporters, with rates that decreased with domain size. Establishment periods varied from a few generations in wild type to >200 generations in the longest region examined, and highly correlated silencing of two reporters located outside the nucleation center was observed. On a chromatin level, this indicates that individual regions are silenced in sudden bursts. Mathematical modeling accounts for these bursts if heterochromatic nucleosomes facilitate a deacetylation or methylation reaction at long range, in a distance-independent manner. A likely effector of three-dimensional interactions is the evolutionarily conserved Swi6^HP1 H3K9me reader, indicating the bursting behavior might be a general mode of heterochromatin propagation.

Dri1 mediates heterochromatin assembly via RNAi and histone deacetylation

Heterochromatin, a transcriptionally silenced chromatin domain, is important for genome stability and gene expression. Histone 3 lysine 9 methylation (H3K9me) and histone hypoacetylation are conserved epigenetic hallmarks of heterochromatin. In fission yeast, RNA interference (RNAi) plays a key role in H3K9 methylation and heterochromatin silencing. However, how RNAi machinery and histone deacetylases (HDACs) are coordinated to ensure proper heterochromatin assembly is still unclear. Previously, we showed that Dpb4, a conserved DNA polymerase epsilon subunit, plays a key role in the recruitment of HDACs to heterochromatin during S phase. Here, we identified a novel RNA-binding protein Dri1 that interacts with Dpb4. GFP-tagged Dri1 forms distinct foci mostly in the nucleus, showing a high degree of colocalization with Swi6/Heterochromatin Protein 1. Deletion of dri1+ leads to defects in silencing, H3K9me, and heterochromatic siRNA generation. We also showed that Dri1 physically associates with heterochromatic transcripts, and is required for the recruitment of the RNA-induced transcriptional silencing (RITS) complex via interacting with the complex. Furthermore, loss of Dri1 decreases the association of the Sir2 HDAC with heterochromatin. We further demonstrated that the C-terminus of Dri1 that includes an intrinsically disordered (IDR) region and three zinc fingers is crucial for its role in silencing. Together, our evidences suggest that Dri1 facilitates heterochromatin assembly via the RNAi pathway and HDAC.

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.

Epigenomic differences in the human and chimpanzee genomes are associated with structural variation

Structural variation (SV), including insertions and deletions (indels), is a primary mechanism of genome evolution. However, the mechanism by which SV contributes to epigenome evolution is poorly understood. In this study, we characterized the association between lineage-specific indels and epigenome differences between human and chimpanzee to investigate how SVs might have shaped the epigenetic landscape. By intersecting medium-to-large human-chimpanzee indels (20 bp-50 kb) with putative promoters and enhancers in cranial neural crest cells (CNCCs) and repressed regions in induced pluripotent cells (iPSCs), we found that 12% of indels overlap putative regulatory and repressed regions (RRRs), and 15% of these indels are associated with lineage-biased RRRs. Indel-associated putative enhancer and repressive regions are approximately 1.3 times and approximately three times as likely to be lineage-biased, respectively, as those not associated with indels. We found a twofold enrichment of medium-sized indels (20-50 bp) in CpG island (CGI)-containing promoters than expected by chance. Lastly, from human-specific transposable element insertions, we identified putative regulatory elements, including NR2F1-bound putative CNCC enhancers derived from SVAs and putative iPSC promoters derived from LTR5s. Our results show that different types of indels are associated with specific epigenomic diversity between human and chimpanzee.

SIR2 Expression Noise Can Generate Heterogeneity in Viability but Does Not Affect Cell-to-Cell Epigenetic Silencing of Subtelomeric URA3 in Yeast

Chromatin structure clearly modulates gene expression noise, but the reverse influence has never been investigated, namely how the cell-to-cell expression heterogeneity of chromatin modifiers may generate variable rates of epigenetic modification. Sir2 is a well-characterized histone deacetylase of the Sirtuin family. It strongly influences chromatin silencing, especially at telomeres, subtelomeres and rDNA. This ability to influence epigenetic landscapes makes it a good model to study the largely unexplored interplay between gene expression noise and other epigenetic processes leading to phenotypic diversification. Here, we addressed this question by investigating whether noise in the expression of SIR2 was associated with cell-to-cell heterogeneity in the frequency of epigenetic silencing at subtelomeres in Saccharomyces cerevisiae Using cell sorting to isolate subpopulations with various expression levels, we found that heterogeneity in the cellular concentration of Sir2 does not lead to heterogeneity in the epigenetic silencing of subtelomeric URA3 between these subpopulations. We also noticed that SIR2 expression noise can generate cell-to-cell variability in viability, with lower levels being associated with better viability. This work shows that SIR2 expression fluctuations are not sufficient to generate cell-to-cell heterogeneity in the epigenetic silencing of URA3 at subtelomeres in Saccharomyces cerevisiae but can strongly affect cellular viability.

Integrity of a heterochromatic domain ensured by its boundary elements

The fission yeast silent mating-type region provides an excellent system to ask how chromatic domains with opposite effects on gene expression coexist side by side along chromosomes and to investigate roles played by DNA elements and architectural proteins in the phenomenon. By showing that the IR-L and IR-R chromatin boundaries favor heterochromatin formation in the domain that separates them, dependent on each other and on binding sites for the architectural factor TFIIIC, our work brings to light an important function of these elements and supports the notion that similar types of interactions between boundaries might in other organisms as well stimulate heterochromatin formation in intervening chromosomal loops to actively shape gene expression landscapes.

In fission yeast, the inverted repeats IR-L and IR-R function as boundary elements at the edges of a 20-kb silent heterochromatic domain where nucleosomes are methylated at histone H3K9. Each repeat contains a series of B-box motifs physically associated with the architectural TFIIIC complex and with other factors including the replication regulator Sap1 and the Rix1 complex (RIXC). We demonstrate here the activity of these repeats in heterochromatin formation and maintenance. Deletion of the entire IR-R repeat or, to a lesser degree, deletion of just the B boxes impaired the de novo establishment of the heterochromatic domain. Nucleation proceeded normally at the RNA interference (RNAi)-dependent element cenH but subsequent propagation to the rest of the region occurred at reduced rates in the mutants. Once established, heterochromatin was unstable in the mutants. These defects resulted in bistable populations of cells occupying alternate “on” and “off” epigenetic states. Deleting IR-L in combination with IR-R synergistically tipped the balance toward the derepressed state, revealing a concerted action of the two boundaries at a distance. The nuclear rim protein Amo1 has been proposed to tether the mating-type region and its boundaries to the nuclear envelope, where Amo1 mutants displayed milder phenotypes than boundary mutants. Thus, the boundaries might facilitate heterochromatin propagation and maintenance in ways other than just through Amo1, perhaps by constraining a looped domain through pairing.

Chromatin state switching in a polymer model with mark-conformation coupling

We investigate the phase transition properties of the polymer-Potts model, a chain composed of monomers with magnetic degrees of freedom, with the motivation to study the conformation and mark switching dynamics of chromatin. By the mean-field approximation, we find that the phase transition between the swollen-disordered state and the compact-ordered state is discrete; it is first order as in the long-range Potts model, but with a significantly larger jump in magnetization (i.e., mark coherence) upon the ordering transition. The results imply how small changes in epigenetic writer concentrations can lead to a macroscopic switching of the chromatin state, suggesting a simple mechanism of discrete switching observed, for instance, in cell differentiation.

Epigenomics in 3D: importance of long-range spreading and specific interactions in epigenomic maintenance

Recent progresses of genome-wide chromatin conformation capture techniques have shown that the genome is segmented into hierarchically organized spatial compartments. However, whether this non-random 3D organization only reflects or indeed contributes—and how—to the regulation of genome function remain to be elucidated. The observation in many species that 3D domains correlate strongly with the 1D epigenomic information along the genome suggests a dynamic coupling between chromatin organization and epigenetic regulation. Here, we posit that chromosome folding may contribute to the maintenance of a robust epigenomic identity via the formation of spatial compartments like topologically-associating domains. Using a novel theoretical framework, the living chromatin model, we show that 3D compartmentalization leads to the spatial colocalization of epigenome regulators, thus increasing their local concentration and enhancing their ability to spread an epigenomic signal at long-range. Interestingly, we find that the presence of 1D insulator elements, like CTCF, may contribute greatly to the stable maintenance of adjacent antagonistic epigenomic domains. We discuss the generic implications of our findings in the light of various biological contexts from yeast to human. Our approach provides a modular framework to improve our understanding and to investigate in details the coupling between the structure and function of chromatin.

Storing memories: the distinct phases of Polycomb-mediated silencing of Arabidopsis FLC

Polycomb-mediated epigenetic silencing is central to correct growth and development in higher eukaryotes. The evolutionarily conserved Polycomb repressive complex 2 (PRC2) transcriptionally silences target genes through a mechanism requiring the histone modification H3K27me3. However, we still do not fully understand what defines Polycomb targets, how their expression state is switched from epigenetically ON to OFF and how silencing is subsequently maintained through many cell divisions. An excellent system in which to dissect the sequence of events underlying an epigenetic switch is the Arabidopsis FLC locus. Exposure to cold temperatures progressively induces a PRC2-dependent switch in an increasing proportion of cells, through a mechanism that is driven by the local chromatin environment. Temporally distinct phases of this silencing mechanism have been identified. First, the locus is transcriptionally silenced in a process involving cold-induced antisense transcripts; second, nucleation at the first exon/intron boundary of a Polycomb complex containing cold-induced accessory proteins induces a metastable epigenetically silenced state; third, a Polycomb complex with a distinct composition spreads across the locus in a process requiring DNA replication to deliver long-term epigenetic silencing. Detailed understanding from this system is likely to provide mechanistic insights important for epigenetic silencing in eukaryotes generally.

RNA-induced initiation of transcriptional silencing (RITS) complex structure and function

Heterochromatic regions of the genome are epigenetically regulated to maintain a heritable '"silent state"'. In fission yeast and other organisms, epigenetic silencing is guided by nascent transcripts, which are targeted by the RNA interference pathway. The key effector complex of the RNA interference pathway consists of small interfering RNA molecules (siRNAs) associated with Argonaute, assembled into the RNA-induced transcriptional silencing (RITS) complex. This review focuses on our current understanding of how RITS promotes heterochromatin formation, and in particular on the role of Argonaute-containing complexes in many other functions such as quelling, release of RNA polymerases, cellular quiescence and genome defense.

Theoretical analysis of Polycomb-Trithorax systems predicts that poised chromatin is bistable and not bivalent

Polycomb (PcG) and Trithorax (TrxG) group proteins give stable epigenetic memory of silent and active gene expression states, but also allow poised states in pluripotent cells. Here we systematically address the relationship between poised, active and silent chromatin, by integrating 73 publications on PcG/TrxG biochemistry into a mathematical model comprising 144 nucleosome modification states and 8 enzymatic reactions. Our model predicts that poised chromatin is bistable and not bivalent. Bivalent chromatin, containing opposing active and silent modifications, is present as an unstable background population in all system states, and different subtypes co-occur with active and silent chromatin. In contrast, bistability, in which the system switches frequently between stable active and silent states, occurs under a wide range of conditions at the transition between monostable active and silent system states. By proposing that bistability and not bivalency is associated with poised chromatin, this work has implications for understanding the molecular nature of pluripotency.

Polycomb and Trithorax group proteins regulate silent and active gene expression states, but also allow poised states in pluripotent cells. Here the authors present a mathematical model that integrates data on Polycomb/ Trithorax biochemistry into a single coherent framework which predicts that poised chromatin is not bivalent as previously proposed, but is bistable, meaning that the system switches frequently between stable active and silent states.

How can memories last for days, years, or a lifetime? Proposed mechanisms for maintaining synaptic potentiation and memory

With memory encoding reliant on persistent changes in the properties of synapses, a key question is how can memories be maintained from days to months or a lifetime given molecular turnover? It is likely that positive feedback loops are necessary to persistently maintain the strength of synapses that participate in encoding. Such feedback may occur within signal-transduction cascades and/or the regulation of translation, and it may occur within specific subcellular compartments or within neuronal networks. Not surprisingly, numerous positive feedback loops have been proposed. Some posited loops operate at the level of biochemical signal-transduction cascades, such as persistent activation of Ca²⁺/calmodulin kinase II (CaMKII) or protein kinase Mζ. Another level consists of feedback loops involving transcriptional, epigenetic and translational pathways, and autocrine actions of growth factors such as BDNF. Finally, at the neuronal network level, recurrent reactivation of cell assemblies encoding memories is likely to be essential for late maintenance of memory. These levels are not isolated, but linked by shared components of feedback loops. Here, we review characteristics of some commonly discussed feedback loops proposed to underlie the maintenance of memory and long-term synaptic plasticity, assess evidence for and against their necessity, and suggest experiments that could further delineate the dynamics of these feedback loops. We also discuss crosstalk between proposed loops, and ways in which such interaction can facilitate the rapidity and robustness of memory formation and storage.

Noncoding RNA-nucleated heterochromatin spreading is intrinsically labile and requires accessory elements for epigenetic stability

The heterochromatin spreading reaction is a central contributor to the formation of gene-repressive structures, which are re-established with high positional precision, or fidelity, following replication. How the spreading reaction contributes to this fidelity is not clear. To resolve the origins of stable inheritance of repression, we probed the intrinsic character of spreading events in fission yeast using a system that quantitatively describes the spreading reaction in live single cells. We show that spreading triggered by noncoding RNA-nucleated elements is stochastic, multimodal, and fluctuates dynamically across time. This lack of stability correlates with high histone turnover. At the mating type locus, this unstable behavior is restrained by an accessory cis-acting element REIII, which represses histone turnover. Further, REIII safeguards epigenetic memory against environmental perturbations. Our results suggest that the most prevalent type of spreading, driven by noncoding RNA-nucleators, is epigenetically unstable and requires collaboration with accessory elements to achieve high fidelity.

A Heterochromatin Domain Forms Gradually at a New Telomere and Is Dynamic at Stable Telomeres

Heterochromatin domains play important roles in chromosome biology, organismal development, and aging, including centromere function, mammalian female X chromosome inactivation, and senescence-associated heterochromatin foci. In the fission yeast Schizosaccharomyces pombe and metazoans, heterochromatin contains histone H3 that is dimethylated at lysine 9.

Heterochromatin domains play important roles in chromosome biology, organismal development, and aging, including centromere function, mammalian female X chromosome inactivation, and senescence-associated heterochromatin foci. In the fission yeast Schizosaccharomyces pombe and metazoans, heterochromatin contains histone H3 that is dimethylated at lysine 9. While factors required for heterochromatin have been identified, the dynamics of heterochromatin formation are poorly understood. Telomeres convert adjacent chromatin into heterochromatin. To form a new heterochromatic region in S. pombe, an inducible DNA double-strand break (DSB) was engineered next to 48 bp of telomere repeats in euchromatin, which caused formation of a new telomere and the establishment and gradual spreading of a new heterochromatin domain. However, spreading was dynamic even after the telomere had reached its stable length, with reporter genes within the heterochromatin domain showing variegated expression. The system also revealed the presence of repeats located near the boundaries of euchromatin and heterochromatin that are oriented to allow the efficient healing of a euchromatic DSB to cap the chromosome end with a new telomere. Telomere formation in S. pombe therefore reveals novel aspects of heterochromatin dynamics and fail-safe mechanisms to repair subtelomeric breaks, with implications for similar processes in metazoan genomes.

Chromosomal barcoding as a tool for multiplexed phenotypic characterization of laboratory evolved lineages

Adaptive laboratory evolution is an important tool to evolve organisms to increased tolerance towards different physical and chemical stress. It is applied to study the evolution of antibiotic resistance as well as genetic mechanisms underlying improvements in production strains. Adaptive evolution experiments can be automated in a high-throughput fashion. However, the characterization of the resulting lineages can become a time consuming task, when the performance of each lineage is evaluated individually. Here, we present a novel method for the markerless insertion of randomized genetic barcodes into the genome of Escherichia coli using a novel dual-auxotrophic selection approach. The barcoded E. coli library allows multiplexed phenotyping of evolved strains in pooled competition experiments. We use the barcoded library in an adaptive evolution experiment; evolving resistance towards three common antibiotics. Comparing this multiplexed phenotyping with conventional susceptibility testing and growth-rate measurements we can show a significant positive correlation between the two approaches. Use of barcoded bacterial strain libraries for individual adaptive evolution experiments drastically reduces the workload of characterizing the resulting phenotypes and enables prioritization of lineages for in-depth characterization. In addition, barcoded clones open up new ways to profile community dynamics or to track lineages in vivo or situ.

Ten principles of heterochromatin formation and function

Heterochromatin is a key architectural feature of eukaryotic chromosomes, which endows particular genomic domains with specific functional properties. The capacity of heterochromatin to restrain the activity of mobile elements, isolate DNA repair in repetitive regions and ensure accurate chromosome segregation is crucial for maintaining genomic stability. Nucleosomes at heterochromatin regions display histone post-translational modifications that contribute to developmental regulation by restricting lineage-specific gene expression. The mechanisms of heterochromatin establishment and of heterochromatin maintenance are separable and involve the ability of sequence-specific factors bound to nascent transcripts to recruit chromatin-modifying enzymes. Heterochromatin can spread along the chromatin from nucleation sites. The propensity of heterochromatin to promote its own spreading and inheritance is counteracted by inhibitory factors. Because of its importance for chromosome function, heterochromatin has key roles in the pathogenesis of various human diseases. In this Review, we discuss conserved principles of heterochromatin formation and function using selected examples from studies of a range of eukaryotes, from yeast to human, with an emphasis on insights obtained from unicellular model organisms.

Shaping epigenetic memory via genomic bookmarking

Reconciling the stability of epigenetic patterns with the rapid turnover of histone modifications and their adaptability to external stimuli is an outstanding challenge. Here, we propose a new biophysical mechanism that can establish and maintain robust yet plastic epigenetic domains via genomic bookmarking (GBM). We model chromatin as a recolourable polymer whose segments bear non-permanent histone marks (or colours) which can be modified by 'writer' proteins. The three-dimensional chromatin organisation is mediated by protein bridges, or 'readers', such as Polycomb Repressive Complexes and Transcription Factors. The coupling between readers and writers drives spreading of biochemical marks and sustains the memory of local chromatin states across replication and mitosis. In contrast, GBM-targeted perturbations destabilise the epigenetic patterns. Strikingly, we demonstrate that GBM alone can explain the full distribution of Polycomb marks in a whole Drosophila chromosome. We finally suggest that our model provides a starting point for an understanding of the biophysics of cellular differentiation and reprogramming.

Inferring gene regulatory networks from single-cell data: a mechanistic approach

Background

The recent development of single-cell transcriptomics has enabled gene expression to be measured in individual cells instead of being population-averaged. Despite this considerable precision improvement, inferring regulatory networks remains challenging because stochasticity now proves to play a fundamental role in gene expression. In particular, mRNA synthesis is now acknowledged to occur in a highly bursty manner.

Results

We propose to view the inference problem as a fitting procedure for a mechanistic gene network model that is inherently stochastic and takes not only protein, but also mRNA levels into account. We first explain how to build and simulate this network model based upon the coupling of genes that are described as piecewise-deterministic Markov processes. Our model is modular and can be used to implement various biochemical hypotheses including causal interactions between genes. However, a naive fitting procedure would be intractable. By performing a relevant approximation of the stationary distribution, we derive a tractable procedure that corresponds to a statistical hidden Markov model with interpretable parameters. This approximation turns out to be extremely close to the theoretical distribution in the case of a simple toggle-switch, and we show that it can indeed fit real single-cell data. As a first step toward inference, our approach was applied to a number of simple two-gene networks simulated in silico from the mechanistic model and satisfactorily recovered the original networks.

Conclusions

Our results demonstrate that functional interactions between genes can be inferred from the distribution of a mechanistic, dynamical stochastic model that is able to describe gene expression in individual cells. This approach seems promising in relation to the current explosion of single-cell expression data.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-017-0487-0) contains supplementary material, which is available to authorized users.

Noncoding RNAs in Polycomb and Trithorax Regulation: A Quantitative Perspective

The question of how noncoding RNAs are involved in Polycomb group (PcG) and Trithorax group (TrxG) regulation has been on an extraordinary journey over the last three decades. Favored models have risen and fallen, and healthy debates have swept back and forth. The field has recently reached a critical mass of compelling data that throws light on several previously unresolved issues. The time is ripe for a fruitful combination of these findings with two other long-running avenues of research, namely the biochemical properties of the PcG/TrxG system and the application of theoretical mathematical models toward an understanding of the system's regulatory properties. I propose that integrating our current knowledge of noncoding RNA into a quantitative biochemical and theoretical framework for PcG and TrxG regulation has the potential to reconcile several apparently conflicting models and identifies fascinating questions for future research.

Distinct phases of Polycomb silencing to hold epigenetic memory of cold in Arabidopsis

Gene silencing by Polycomb complexes is central to eukaryotic development. Cold-induced epigenetic repression of FLOWERING LOCUS C (FLC) in the plant Arabidopsis provides an opportunity to study initiation and maintenance of Polycomb silencing. Here, we show that a subset of Polycomb repressive complex 2 factors nucleate silencing in a small region within FLC, locally increasing H3K27me3 levels. This nucleation confers a silenced state that is metastably inherited, with memory held in the local chromatin. Metastable memory is then converted to stable epigenetic silencing through separate Polycomb factors, which spread across the locus after cold to enlarge the domain that contains H3K27me3. Polycomb silencing at FLC thus has mechanistically distinct phases, which involve specialization of distinct Polycomb components to deliver first metastable then long-term epigenetic silencing.

Models of life: epigenetics, diversity and cycles

This review emphasizes aspects of biology that can be understood through repeated applications of simple causal rules. The selected topics include perspectives on gene regulation, phage lambda development, epigenetics, microbial ecology, as well as model approaches to diversity and to punctuated equilibrium in evolution. Two outstanding features are repeatedly described. One is the minimal number of rules to sustain specific states of complex systems for a long time. The other is the collapse of such states and the subsequent dynamical cycle of situations that restitute the system to a potentially new metastable state.

Perspectives: using polymer modeling to understand the formation and function of nuclear compartments

Compartmentalization is a ubiquitous feature of cellular function. In the nucleus, early observations revealed a non-random spatial organization of the genome with a large-scale segregation between transcriptionally active—euchromatin—and silenced—heterochromatin—parts of the genome. Recent advances in genome-wide mapping and imaging techniques have strikingly improved the resolution at which nuclear genome folding can be analyzed and have revealed a multiscale spatial compartmentalization with increasing evidences that such compartment may indeed result from and participate to genome function. Understanding the underlying mechanisms of genome folding and in particular the link to gene regulation requires a cross-disciplinary approach that combines the new high-resolution techniques with computational modeling of chromatin and chromosomes. In this perspective article, we first present how the copolymer theoretical framework can account for the genome compartmentalization. We then suggest, in a second part, that compartments may act as a “nanoreactor,” increasing the robustness of either activation or repression by enhancing the local concentration of regulators. We conclude with the need to develop a new framework, namely the “living chromatin” model that will allow to explicitly investigate the coupling between spatial compartmentalization and gene regulation.

Coupling 1D modifications and 3D nuclear organization: data, models and function

Over the past decade, advances in molecular methods have strikingly improved the resolution at which nuclear genome folding can be analyzed. This revealed a wealth of conserved features organizing the one dimensional DNA molecule into tridimensional nuclear domains. In this review, we briefly summarize the main findings and highlight how models based on polymer physics shed light on the principles underlying the formation of these domains. Finally, we discuss the mechanistic similarities allowing self-organization of these structures and the functional importance of these in the maintenance of transcriptional programs.