Abstract In the last decade, the interactions among histone modifications and DNA methylation and their effect on the DNA structure, i.e., chromatin state, have been identified as key mediators for the maintenance of cell identity, defined as epigenetic cell memory. In this paper, we determine how the positive feedback loops generated by the auto- and cross-catalysis among repressive modifications affect the temporal duration of the cell identity. To this end, we conduct a stochastic analysis of a recently published chromatin modification circuit considering two limiting behaviors: fast erasure rate of repressive histone modifications or fast erasure rate of DNA methylation. In order to perform this mathematical analysis, we first show that the deterministic model of the system is a singular singularly perturbed (SSP) system and use a model reduction approach for SSP systems to obtain a reduced one-dimensional model. We thus analytically evaluate the reduced system’s stationary probability distribution and the mean switching time between active and repressed chromatin states. We then add a computational study of the original reaction model to validate and extend the analytical findings. Our results show that the absence of DNA methylation reduces the bias of the system’s stationary probability distribution toward the repressed chromatin state and the temporal duration of this state’s memory. In the absence of repressive histone modifications, we also observe that the time needed to reactivate a repressed gene with an activating input is less stochastic, suggesting that repressive histone modifications specifically contribute to the highly variable latency of state reactivation.
more »
« less
Epigenetic cell memory: The gene’s inner chromatin modification circuit
Epigenetic cell memory allows distinct gene expression patterns to persist in different cell types despite a common genotype. Although different patterns can be maintained by the concerted action of transcription factors (TFs), it was proposed that long-term persistence hinges on chromatin state. Here, we study how the dynamics of chromatin state affect memory, and focus on a biologically motivated circuit motif, among histones and DNA modifications, that mediates the action of TFs on gene expression. Memory arises from time-scale separation among three circuit’s constituent processes: basal erasure, auto and cross-catalysis, and recruited erasure of modifications. When the two latter processes are sufficiently faster than the former, the circuit exhibits bistability and hysteresis, allowing active and repressed gene states to coexist and persist after TF stimulus removal. The duration of memory is stochastic with a mean value that increases as time-scale separation increases, but more so for the repressed state. This asymmetry stems from the cross-catalysis between repressive histone modifications and DNA methylation and is enhanced by the relatively slower decay rate of the latter. Nevertheless, TF-mediated positive autoregulation can rebalance this asymmetry and even confers robustness of active states to repressive stimuli. More generally, by wiring positively autoregulated chromatin modification circuits under time scale separation, long-term distinct gene expression patterns arise, which are also robust to failure in the regulatory links.
more »
« less
- NSF-PAR ID:
- 10338932
- Editor(s):
- Herrmann, Carl
- Date Published:
- Journal Name:
- PLOS Computational Biology
- Volume:
- 18
- Issue:
- 4
- ISSN:
- 1553-7358
- Page Range / eLocation ID:
- e1009961
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Epigenetic cell memory (ECM), the inheritance of gene expression patterns without changes in genetic sequence, is a critical property of multi-cellular organisms. Chromatin state, as dictated by histone covalent modifications, has recently appeared as a mediator of ECM. In this paper, we conduct a stochastic analysis of the histone modification circuit that controls chromatin state to determine key biological parameters that affect ECM. Specifically, we derive a one-dimensional Markov chain model of the circuit and analytically evaluate both the stationary probability distribution of chromatin state and the mean time to switch between active and repressed chromatin states. We then validate our analytical findings using stochastic simulations of the original higher dimensional circuit reaction model. Our analysis shows that as the speed of basal decay of histone modifications decreases compared to the speed of autocatalysis, the stationary probability distribution becomes bimodal and increasingly concentrated about the active and repressed chromatin states. Accordingly, the switching time between active and repressed chromatin states becomes larger. These results indicate that time scale separation among key constituent processes of the histone modification circuit controls ECM.more » « less
-
Epigenetic cell memory (ECM),the inheritance of gene expression patterns without changes in genetic sequence,is acritical property of multi-cellular organisms.Chromatin state, as dictated by histone covalent modifications, has recently appeared as a mediator of ECM. In this paper,we conduct a stochastic analysis of the histone modification circuit that controls chromatin state to determine key biological parameters that affect ECM. Specifically, we derive a one-dimensional Markov chain model of the circuit and analytically evaluate both the stationary probability distribution of chromatin state and the mean time to switch between active and repressed chromatin states.We then validate our analytical findings using stochastic simulations of the original higher dimen- sional circuit reaction model.Our analysis shows that as the speed of basal decay of histone modifications decreases compared to the speed of autocatalysis,the stationary probability distribution becomes bimodal and increasingly concentrated about the active and repressed chromatin states. Accordingly, the switching time between active and repressed chromatin states becomes larger.These results indicate that timescale separation among key constituent processes of the histone modification circuit controls ECM.more » « less
-
INTRODUCTION Neurons are by far the most diverse of all cell types in animals, to the extent that “cell types” in mammalian brains are still mostly heterogeneous groups, and there is no consensus definition of the term. The Drosophila optic lobes, with approximately 200 well-defined cell types, provides a tractable system with which to address the genetic basis of neuronal type diversity. We previously characterized the distinct developmental gene expression program of each of these types using single-cell RNA sequencing (scRNA-seq), with one-to-one correspondence to the known morphological types. RATIONALE The identity of fly neurons is determined by temporal and spatial patterning mechanisms in stem cell progenitors, but it remained unclear how these cell fate decisions are implemented and maintained in postmitotic neurons. It was proposed in Caenorhabditis elegans that unique combinations of terminal selector transcription factors (TFs) that are continuously expressed in each neuron control nearly all of its type-specific gene expression. This model implies that it should be possible to engineer predictable and complete switches of identity between different neurons just by modifying these sustained TFs. We aimed to test this prediction in the Drosophila visual system. RESULTS Here, we used our developmental scRNA-seq atlases to identify the potential terminal selector genes in all optic lobe neurons. We found unique combinations of, on average, 10 differentially expressed and stably maintained (across all stages of development) TFs in each neuron. Through genetic gain- and loss-of-function experiments in postmitotic neurons, we showed that modifications of these selector codes are sufficient to induce predictable switches of identity between various cell types. Combinations of terminal selectors jointly control both developmental (e.g., morphology) and functional (e.g., neurotransmitters and their receptors) features of neurons. The closely related Transmedullary 1 (Tm1), Tm2, Tm4, and Tm6 neurons (see the figure) share a similar code of terminal selectors, but can be distinguished from each other by three TFs that are continuously and specifically expressed in one of these cell types: Drgx in Tm1, Pdm3 in Tm2, and SoxN in Tm6. We showed that the removal of each of these selectors in these cell types reprograms them to the default Tm4 fate. We validated these conversions using both morphological features and molecular markers. In addition, we performed scRNA-seq to show that ectopic expression of pdm3 in Tm4 and Tm6 neurons converts them to neurons with transcriptomes that are nearly indistinguishable from that of wild-type Tm2 neurons. We also show that Drgx expression in Tm1 neurons is regulated by Klumpfuss, a TF expressed in stem cells that instructs this fate in progenitors, establishing a link between the regulatory programs that specify neuronal fates and those that implement them. We identified an intronic enhancer in the Drgx locus whose chromatin is specifically accessible in Tm1 neurons and in which Klu motifs are enriched. Genomic deletion of this region knocked down Drgx expression specifically in Tm1 neurons, leaving it intact in the other cell types that normally express it. We further validated this concept by demonstrating that ectopic expression of Vsx (visual system homeobox) genes in Mi15 neurons not only converts them morphologically to Dm2 neurons, but also leads to the loss of their aminergic identity. Our results suggest that selector combinations can be further sculpted by receptor tyrosine kinase signaling after neurogenesis, providing a potential mechanism for postmitotic plasticity of neuronal fates. Finally, we combined our transcriptomic datasets with previously generated chromatin accessibility datasets to understand the mechanisms that control brain wiring downstream of terminal selectors. We built predictive computational models of gene regulatory networks using the Inferelator framework. Experimental validations of these networks revealed how selectors interact with ecdysone-responsive TFs to activate a large and specific repertoire of cell surface proteins and other effectors in each neuron at the onset of synapse formation. We showed that these network models can be used to identify downstream effectors that mediate specific cellular decisions during circuit formation. For instance, reduced levels of cut expression in Tm2 neurons, because of its negative regulation by pdm3 , controls the synaptic layer targeting of their axons. Knockdown of cut in Tm1 neurons is sufficient to redirect their axons to the Tm2 layer in the lobula neuropil without affecting other morphological features. CONCLUSION Our results support a model in which neuronal type identity is primarily determined by a relatively simple code of continuously expressed terminal selector TFs in each cell type throughout development. Our results provide a unified framework of how specific fates are initiated and maintained in postmitotic neurons and open new avenues to understanding synaptic specificity through gene regulatory networks. The conservation of this regulatory logic in both C. elegans and Drosophila makes it likely that the terminal selector concept will also be useful in understanding and manipulating the neuronal diversity of mammalian brains. Terminal selectors enable predictive cell fate reprogramming. Tm1, Tm2, Tm4, and Tm6 neurons of the Drosophila visual system share a core set of TFs continuously expressed by each cell type (simplified). The default Tm4 fate is overridden by the expression of a single additional terminal selector to generate Tm1 ( Drgx ), Tm2 ( pdm3 ), or Tm6 ( SoxN ) fates.more » « less
-
more » « less
Polycomb Group (PcG) proteins are part of an epigenetic cell memory system that plays essential roles in multicellular development, stem cell biology, X chromosome inactivation, and cancer. In animals, plants, and many fungi, Polycomb Repressive Complex 2 (PRC2) catalyzes trimethylation of histone H3 lysine 27 (H3K27me3) to assemble transcriptionally repressed facultative heterochromatin. PRC2 is structurally and functionally conserved in the model fungus
Neurospora crassa , and recent work in this organism has generated insights into PRC2 control and function. To identify components of the facultative heterochromatin pathway, we performed a targeted screen ofNeurospora deletion strains lacking individual ATP-dependent chromatin remodeling enzymes. We found theNeurospora homolog of IMITATION SWITCH (ISW) is critical for normal transcriptional repression, nucleosome organization, and establishment of typical histone methylation patterns in facultative heterochromatin domains. We also found that stable interaction between PRC2 and chromatin depends on ISW. A functional ISW ATPase domain is required for gene repression and normal H3K27 methylation. ISW homologs interact with accessory proteins to form multiple complexes with distinct functions. Using proteomics and molecular approaches, we identified three distinctNeurospora ISW-containing complexes. A triple mutant lacking three ISW accessory factors and disrupting multiple ISW complexes led to widespread up-regulation of PRC2 target genes and altered H3K27 methylation patterns, similar to an ISW-deficient strain. Taken together, our data show that ISW is a key component of the facultative heterochromatin pathway inNeurospora , and that distinct ISW complexes perform an apparently overlapping role to regulate chromatin structure and gene repression at PRC2 target domains.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1371/journal.pcbi.1009961
