Is There a Histone Code for Cellular Quiescence?

Many of the cells in our bodies are quiescent, that is, temporarily not dividing. Under certain physiological conditions such as during tissue repair and maintenance, quiescent cells receive the appropriate stimulus and are induced to enter the cell cycle. The ability of cells to successfully transition into and out of a quiescent state is crucial for many biological processes including wound healing, stem cell maintenance, and immunological responses. Across species and tissues, transcriptional, epigenetic, and chromosomal changes associated with the transition between proliferation and quiescence have been analyzed, and some consistent changes associated with quiescence have been identified. Histone modifications have been shown to play a role in chromatin packing and accessibility, nucleosome mobility, gene expression, and chromosome arrangement. In this review, we critically evaluate the role of different histone marks in these processes during quiescence entry and exit. We consider different model systems for quiescence, each of the most frequently monitored candidate histone marks, and the role of their writers, erasers and readers. We highlight data that support these marks contributing to the changes observed with quiescence. We specifically ask whether there is a quiescence histone “code,” a mechanism whereby the language encoded by specific combinations of histone marks is read and relayed downstream to modulate cell state and function. We conclude by highlighting emerging technologies that can be applied to gain greater insight into the role of a histone code for quiescence.

cells can be distinguished from other types of non-dividing cells such as senescent or terminally differentiated cells by their temporary exit from the cell cycle and high likelihood of proliferating in response to a triggering stimulus (Sang and Coller, 2009;Cheung and Rando, 2013;Terzi et al., 2016). Quiescent cells must therefore preserve the ability to proliferate at a later time, and protect themselves from entering irreversible states (Coller et al., 2006;Sang and Coller, 2009;Sang et al., 2010;Bjornson et al., 2012). Cellular quiescence has been studied experimentally in multiple systems including yeast, cultured primary cells, and stem cells (Mitra et al., 2018a;Spain et al., 2018;Yang and Chi, 2018) ( Table 1). Some of the gene expression, signaling, and functional changes observed with quiescence are likely specific for a cell type, while others are shared. Transcriptional changes with quiescence have been analyzed using cDNA libraries (Schneider et al., 1988;Coppock et al., 1993), microarrays (Venezia et al., 2004;Coller et al., 2006;Suh et al., 2012;Johnson et al., 2017), next generation sequencing (van Velthoven et al., 2017;Mitra et al., 2018b;Srivastava et al., 2018), and single-cell RNA sequencing methods (Kalakonda et al., 2008;Coller, 2019a). These studies demonstrated widespread gene expression changes with quiescence, some of which are functionally important for the quiescent state (Suh et al., 2012;Johnson et al., 2017;Lee H.N. et al., 2018;Mitra et al., 2018b). These gene expression changes include downregulation of genes involved in cell cycle progression and upregulation of stress response genes (Lemons et al., 2010;Legesse-Miller et al., 2012;Valcourt et al., 2012;Coller, 2019b). Other gene expression changes allow the cell to re-organize metabolic pathways in quiescent cells to better match the availability of nutrients and the metabolic needs of the cell (Lemons et al., 2010;Valcourt et al., 2012;Coller, 2019b). Disruption of these cellular mechanisms can contribute to the occurrence and progression of pathologies related to aging, developmental defects, and cancer (Tumpel and Rudolph, 2019).
In addition to gene expression changes, quiescence is also associated with changes in the packaging of DNA into chromatin. Eukaryotic chromatin can take on two forms-a more condensed and transcriptionally silent form called heterochromatin and a less condensed and more transcriptionally active form called euchromatin (DesJarlais and Tummino, 2016). Within these states, the extent of compaction can vary, for instance, mitotic chromosomes are extremely condensed. Studies using imaging, flow cytometry, Hi-C, and other methods have shown that entry to a quiescent state involves changes in nuclear size, chromatin compaction and 3D genome architecture (Bridger et al., 2000;Evertts et al., 2013a;Guidi et al., 2015;Criscione et al., 2016;Swygert et al., 2019Swygert et al., , 2021. In yeast, the transition from exponential phase growth to stationary phase, a quiescent state achieved when yeast deplete their nutrients, is associated with downregulation of gene expression, a more condensed chromatin state (Martinez et al., 2004;Schafer et al., 2008), and more long-range chromosomal interactions (Swygert et al., 2019). In mammals, activation of quiescent lymphocytes is associated with an unpacking of condensed chromatin in a process that can be visualized with electron microscopy (Tokuyasu et al., 1968;Dardick et al., 1983;Setterfield et al., 1983;Grigoryev et al., 2004).
Using circular dichroism, Chiu and Baserga reported a likely change to a more open chromatin structure as quiescent fibroblasts re-enter the cell cycle (Chiu and Baserga, 1975). In contrast, in one study, bovine fibroblasts were reported to have a more relaxed chromatin state in G 0 (quiescent) compared with G 1 cells (Kallingappa et al., 2016).
In addition to changes in gene expression and chromatin compaction, quiescence is also associated with a change in the positioning of chromosomes within the nucleus. In yeast, hyperclustering of telomeres has been reported with quiescence (Guidi et al., 2015;Laporte et al., 2016). When serum was removed from the culture medium of human fibroblasts, chromosomes were repositioned within 15 min in a process that required ATP, actin polymerization, and myosin (Mehta et al., 2010). In another study in human dermal fibroblasts, gene-poor chromosome 18 was found near the edge of the nucleus and gene-rich chromosome 19 was found in the center of the nucleus in proliferating cells. In serum-starved, quiescent fibroblasts, chromosome 18 shifted to a more central location in the nucleus, and there was no longer a difference in the positioning of chromosomes 18 and 19 (Bridger et al., 2000). Taken together, these findings demonstrate changes in gene expression, chromatin compaction and chromosome positioning within the nucleus in quiescent cells.

HISTONE POST-TRANSLATIONAL MODIFICATIONS AS A POSSIBLE BIOLOGICAL CODE Nucleosome Structure and Histone Marks
Eukaryotic genomic DNA is tightly packed inside the nucleus. For mammalian chromosomes, this tight packing results in a 10,000-fold reduction in length (Kornberg and Lorch, 2020). The DNA in chromatin forms complexes with histone proteins that assemble the DNA strands into nucleosomes in a structure that resembles "beads on a string" with the nucleosomes (beads) representing the basic repeating unit of chromatin (Cutter and Hayes, 2015;Zhou et al., 2019;Ghoneim et al., 2021). Each core nucleosome consists of ∼147 base pairs (bp) of DNA in a left-handed super-helical conformation wrapped around an octamer of histone proteins (Zhou et al., 2019). The octamer consists of two copies each of histone proteins H2A, H2B, H3, and H4 with each of the two dimers of H2A-H2B interacting with either end of a (H3-H4) 2 tetramer (H4-H3:H3-H4). The core nucleosome is flanked by 10-70 bp of linker DNA and usually a linker histone (H1) (Cutter and Hayes, 2015). The disordered N-terminal tails of all four histone proteins as well as the C-terminal tail of H2A protrude out from the nucleosome core and are sites of diverse post translational modifications (PTMs) or marks such as lysine and arginine methylation, lysine acetylation, and serine and threonine phosphorylation (Bannister and Kouzarides, 2011;Greer and Shi, 2012;Cutter and Hayes, 2015). These histone tails modulate charge, hydrophobicity, and steric access to chromatin (Ghoneim et al., 2021). Histone PTMs are added and removed by enzymatic proteins referred to as "writers" and "erasers, " respectively (Soshnev et al., 2016;Hyun et al., 2017;Husmann and Gozani, 2019). Histone PTMs serve as recognition sites for proteins ("readers") that site-specifically bind to chromatin. In some cases, a single protein contains multiple domains and can act as both a reader for one type of PTM and a writer for another PTM (Smeenk and Mailand, 2016). The amino acid residues in the histone globular core can also be post-translationally modified and these core PTMs likely modulate interactions between histones and between histones and DNA (Tessarz and Kouzarides, 2014).

What Are the Properties or Functions of Histone Marks?
Histone H3K4me3, H3K36me3, and H3K79me3 and ubiquitylation of H2B marks are often associated with active transcription (Black et al., 2012;Hyun et al., 2017), whereas H3K9me3, H3K27me3, H2A ubiquitylation on lysine 119, and H4K20 methylation are indicators of a silenced chromatin state with reduced gene expression (Black et al., 2012;Hyun et al., 2017). These properties of histone marks are related to the way they interact with chromatin and chromatin binding proteins. Histone PTMs can be envisioned to function by at least two broad categories of mechanisms (Bannister and Kouzarides, 2011). The first involves direct structural effects on the biomechanical properties of DNA (Bannister and Kouzarides, 2011). In this role, histone PTMs can affect the accessibility of DNA, and thus the binding of transcription factors or other proteins that bind enhancers and affect transcription (Bannister and Kouzarides, 2011). Such effects can occur, for instance, when histone PTMs disrupt electrostatic interactions between histones and DNA. Nucleosome core particle has an overall charge of −150 electrons that is contributed by DNA (−294 electrons) and histones (+ 144 electrons) (Norton et al., 1989;Cortini, 2016;Prakash and Fournier, 2018). Some histone PTMs that are associated with a more open chromatin state and increased gene expression reduce the positive charges on histones thereby leading to less effective screening of the negative charges on DNA (Prakash and Fournier, 2018). Acetylation, in particular, impairs the affinity of histones for DNA by neutralizing the positive charges and disrupting the ionic interactions between histones and DNA. This results in increased histone mobility and a more open chromatin conformation (Allfrey et al., 1964;Cosgrove et al., 2004). An open chromatin conformation facilitates access to transcription factors and other chromatin binding proteins (Bannister and Kouzarides, 2011). While some modifications such as acetylation may be expected to alter the ionic charge and thus chromatin compaction, others, including methylation, may have more modest impacts on charge and chromatin structure (Bannister and Kouzarides, 2011).
The second way in which histone PTMs can exert a functional effect is by regulating the binding of different chromatin factors. As one example, proteins with PHD fingers and Tudor family of domains can bind lysine methylations (Maurer-Stroh et al., 2003;Champagne and Kutateladze, 2009;Bannister and Kouzarides, 2011). In some cases, multiple different domains can recognize the same lysine methylation (Bannister and Kouzarides, 2011). Another example of histone PTM recognition is the binding of dimeric Heterochromatin Protein 1 (HP1) to the H3K9me3 mark via the chromodomain. This is associated with repressive architecture and chromatin compaction (Bannister et al., 2001;Lachner et al., 2001;Bannister and Kouzarides, 2011).
Genome-wide studies of histone marks have revealed that combinations of histone marks can be used to classify chromatin into different states (Black et al., 2012). In Arabidopsis, four chromatin states were identified (Roudier et al., 2011); in Drosophila, 5-9 states have been reported Riddle et al., 2011); while in human cells, up to 51 chromatin states have been defined (Ernst and Kellis, 2010). In human cells, chromatin states include promoter-associated, transcription-associated, active intergenic, large-scale repressed, and repeat-associated states, each of which have distinct histone marks and biological roles (Ernst and Kellis, 2010). Different promoter states were defined by patterns of H3K4 methylation, H3K79 methylation, H4K20 methylation, and acetylation (Black et al., 2012). One particular chromatin state that involves a specific combination of histone marks is the bivalent mark (Bernhart et al., 2016). Bivalent marks are often found in the promoters of developmentally regulated genes (Bernstein et al., 2006;Lesch et al., 2013), and are defined by the simultaneous presence of activating marks such as H3K4me1 or H3K4me3, and the repressive chromatin mark H3K27me3 (Voigt et al., 2013). Genes with bivalent marks are repressed, but pre-loaded with RNA polymerase that is "poised" for rapid expression in response to a relevant trigger (Mikkelsen et al., 2007;Margaritis and Holstege, 2008;Gaertner et al., 2012). More recent studies have identified combinatorial marks that establish zones within the nucleus that can be identified by combinations of proteins and histone marks (Takei et al., 2021). These nuclear zones include nuclear speckles, active chromatin, heterochromatin zones and zones within the nucleolus (Takei et al., 2021). The active chromatin zone, for example, was characterized by histone H3K9ac, H3K27ac, H4K16ac, RNA polymerase II serine 5 phosphorylation, and SF3A66 (Takei et al., 2021).
Consistent with the concept of nuclear zones, histone marks may allow for the patterning of chromatin into regions of approximately 0.5-1 megabases with similar properties, termed topological domains as identified by the Hi-C technique (Prakash and Fournier, 2018). Histone modifications have been found to cluster at the genome scale as DNA tends to fold into domains in which the all of the DNA in that domain is labeled with similar histone marks (Dixon et al., 2012;Rao et al., 2014;Barbieri et al., 2017). For instance, H3K4me3-rich, H3K27me3rich and H3K9me3 rich regions have been found to segregate from each other, and to mark active genes, repressed genes and inactive chromatin, respectively (Prakash and Fournier, 2018). Thus, histone PTMs may be associated with multiple aspects of chromatin including the extent of local compaction, the extent of gene expression and the formation of chromatin domains.

Do Histone Marks Create a Histone Code?
Biological codes that have been previously described include an input system that is translated into an output via adaptors (Prakash and Fournier, 2018). As one example, the genetic code translates sequences of nucleotide codons (input) into a sequence of amino acids (output) using the protein translation apparatus (adapter) (Prakash and Fournier, 2018). Histone PTMs have also been suggested to establish a biological code (Strahl and Allis, 2000;Allis and Jenuwein, 2016). According to the histone code hypothesis, the presence of specific histone marks, and in some cases, possible combinations of histone marks (inputs), provides information to reader proteins (adapters) that interpret the marks or combinations of marks to produce outputs such as gene activation or silencing, chromatin compaction, repair of DNA damage, cell division or differentiation (Strahl and Allis, 2000;Jenuwein and Allis, 2001;Turner, 2002;Smeenk and Mailand, 2016). Given that histone marks tend to be rapidly reestablished after cell division (Evertts et al., 2013a,b), information about a cell's state can be transmitted to descendant cells. Misreading of histone marks has been associated with cancer and developmental defects (Wang and Allis, 2009;Chi et al., 2010;Hyun et al., 2017).
One potential advantage of a histone code would be that combinations of histone marks could provide increased robustness to a system in which different inputs result in specific outcomes (Prakash and Fournier, 2018). Robustness can be achieved with cooperation and redundancy (Prakash and Fournier, 2018). A histone code has been hypothesized to provide a level of proofreading needed so that genes are not turned on or off inappropriately (Prakash and Fournier, 2018). If there are multiple independent histone marks that work in concert to achieve an outcome, then loss of one mark would have only a modest effect on the associated phenotypes (Prakash and Fournier, 2018). Further, comparing the use of histone marks in different species shows that histone modifications have evolutionarily conserved functions and play a similar functional role across eukaryotes (Ho et al., 2014;Prakash and Fournier, 2018).
Generating chromatin states with combinations of histone marks may reflect instances in which the presence of one histone PTM affects the recruitment of enzymes that create another PTM in the same or different histone tails resulting in reproducible output patterns (Figure 1). This can be achieved by multi-domain proteins that recognize the histone PTM through the reader domain and utilize a different domain for recruiting a histone writer. As one example, double-stranded DNA breaks serve as a signal to the ATM kinase, which leads to phosphorylation of the H2A variant H2A.X on its C terminal tail (Rogakou et al., 1998;Smeenk and Mailand, 2016). The presence of this mark, called γ-H2AX, creates binding sites for a reader for this mark, Mediator of DNA Damage Checkpoint Protein 1 (MDC1), a protein that recruits factors to the DNA damage site (Stucki et al., 2005), including E3 ligases that ubiquitinate histones (Huen et al., 2007;Kolas et al., 2007;Mailand et al., 2007;Doil et al., 2009). Histone ubiquitin modifications create recruitment platforms for DNA repair factors (Huen et al., 2007;Kolas et al., 2007;Mailand et al., 2007;Doil et al., 2009). Thus, histone modifications can transmit information about the site-specific presence of double strand breaks to affect an outcome, in this case, DNA repair.
Another example in which the presence of a histone mark affects the likelihood of other marks being added occurs during the mitosis phase of the cell cycle. Phosphorylation of Ser10 of histone H3 regulates transcription during interphase (Shimada et al., 2008) and chromosome condensation during mitosis (Wei et al., 1999). Histone H3S10 phosphorylation prevents phosphorylation of Thr6 and Thr11 on the same histone (Cosgrove, 2012;Liokatis et al., 2012). This hierarchy may ensure that phosphorylation of Ser10 during mitosis, which is required for chromosome condensation and separation, does not lead to subsequent formation of dually labeled histones with H3S10 phosphorylation and Thr6 or Thr11 phosphorylation during mitosis (Cosgrove, 2012;Liokatis et al., 2012). This example FIGURE 1 | Histone Code Hypothesis. A schematic highlighting the known interactions between chromatin modifying enzymes on the H3 and H4 tail as a potential hypothesis of the histone code in quiescence. Lysines that can be methylated or acetylated are gray, and are indicated to have methyl groups (red) or acetyl groups (green). H3K4 can be methylated by methyltransferase, Set1A, and demethylated by JARID1A, a demethylase. Methylation of lysine 4 prevents the methylation of H3K9 by inhibiting Suv39h2, the methyltransferase of lysine 9. Suv39h2 activity can also be prevented by acetylation of H3K14. Removal of H3K4 methylation by JARID1A enhances H3K9 methylation. On the H4 tail, MOF acetylates lysine 16 which inhibits activity of Suv4-20h2 in methylating H4K20. Removal of an acetyl group from H4K16 by SIRT1 enhances H4K20 methylation (Figure made in BioRender).
provides an instance in which histone marks antagonize each other, and combinations of histone marks are required to ensure a robust functional response.
Histone codes that lead to biologically significant outcomes have been suggested to play a role in neural plasticity (Farrelly and Maze, 2019), genome structure (Prakash and Fournier, 2018), and cancer (Godley and Le Beau, 2012). In neurons, serotonylation of histone H3 glutamine 5 works in conjunction with nearby H3K4me3 marks to regulate transcription Farrelly and Maze, 2019). The dual H3K4(me3)(Q5serotonin) mark was found to enhance binding of interacting proteins, including the transcription factor complex TFIID (Lauberth et al., 2013;Farrelly and Maze, 2019). In this way, the combination of histone marks results in increased transcription of specific nearby genes.
On the other hand, it is important to note that the model that specific combinations of histone tail PTMs lead to defined biological outcomes has been challenged (Henikoff and Shilatifard, 2011;Morgan and Shilatifard, 2020). One central argument is whether histone modifications are the cause of different transcriptional states or, instead, are formed as a consequence of transcription and other dynamic processes. It has been difficult to resolve this controversy because studying the direct role of a histone PTM is challenging. Traditional genetic tools such as knockdown or overexpression are not sufficient to differentiate the direct versus indirect effects of targeting a histone modifying enzyme that adds or removes a PTM because the enzyme may also act on non-histone substrates (Henikoff and Shilatifard, 2011;Cornett et al., 2019;Corvalan and Coller, 2021).
An alternative strategy to directly assess the role of a histone PTM is to mutate the amino acid residue that bears the PTM. While this strategy can be effective in some organisms such as yeast, it is not practical for higher eukaryotes due to the presence of multiple copies of the genes encoding the most frequently modified histones (Tripputi et al., 1986;Henikoff and Shilatifard, 2011;Soshnev et al., 2016;Corvalan and Coller, 2021).
Debate about the existence of a histone code has also centered on the nature of the histone code. The original paper describing a histone code suggested a code that has been described as hardwired and deterministic (Jenuwein and Allis, 2001), like the genetic code (Morgan and Shilatifard, 2020). With time, an alternative and more complex relationship between histone marks and functional outcomes has been described (Morgan and Shilatifard, 2020). In this potential representation, histone PTMs convey information in a context-dependent manner (Morgan and Shilatifard, 2020). A histone mark can have multiple potential outcomes, and the specific path would depend on multiple factors including the three-dimensional folding of the genome, the local chromatin environment, and the concentrations of the possible downstream effector molecules (Morgan and Shilatifard, 2020).
Is There a Histone Post-translational Modification Code for Quiescence? Does a histone code exist for quiescence? Are there specific patterns of histone tail PTMs that dictate or are associated with entry, exit, maintenance, or depth of a quiescent state? If so, do these histone PTMs modulate the physical properties of the DNA? Do these histone marks directly alter the chromatin accessibility of gene promoters and enhancers to induce molecular and phenotypic changes with quiescence? Alternatively, do these histone PTMs serve as binding sites for readers that recognize the PTMs and effect cellular changes during quiescence? If the histone marks serve as recognition sites, what are the most important effectors and what aspects of quiescence do they control? In this review, we address these questions and compare the findings from multiple experimental models of quiescence (Table 1). Advances in ChIP-seq technology (O'Geen et al., 2011) such as CUT&RUN (Hainer and Fazzio, 2019), CUT&Tag (Kaya-Okur et al., 2019), and HiChIP (Yan et al., 2014) have enabled fine resolution mapping of the genomic position of different histone tail marks ( Table 2). Mass spectrometry can be used to measure histone epigenetic mark abundance and dynamics in a multiplexed, parallel manner (Volker-Albert et al., 2018) (Table 2). Further, advances in imaging such as combining immunofluorescence Western blot Protein is isolated from sample, separated by weight and probed for on gel with specific antibodies with directly labeled histone modification-specific antibodies to monitor histone levels in single cells (Hayashi-Takanaka et al., 2020), sequential fluorescence in situ hybridization analysis (Takei et al., 2021), and improved imaging with Stochastic Optical Reconstruction Microscopy (STORM) , have permitted detailed analysis of the global and site-specific organization of histone marks (Takei et al., 2021). It is important to note that in this review we focus on the most intensively studied histone PTMs, but these are not the only possible candidates by which histones or a histone code could contribute to the changes observed with quiescence. For instance, there may be new histone marks that are specific for quiescence and have not yet been observed, despite mass spectrometry-based histone tail analysis (Evertts et al., 2013a). There are also linker histones and variants of core histones that affect nucleosome structure and function, and consequently chromatin architecture . Changes in the compositions of histones and histone linkers could potentially contribute to the functional attributes of quiescent cells, and these are not reviewed here. We conclude by identifying areas for future studies and methodologies that can be used to address existing gaps in our knowledge.

DIFFERENT QUIESCENCE MODEL SYSTEMS FOR STUDYING HISTONE MARKS Yeast Models of Quiescence
Multiple model systems have been used to study the molecular mechanisms of quiescence including yeast, mouse and human cells ( Table 1), all of which have different genomes, limiting our ability to make direct comparisons about histone modifications in specific genomic regions. Further, the signals that induce quiescence, and the quiescent state achieved in these model systems differs, which may contribute to differences in the levels of specific histone modifications (Valcourt et al., 2012;Coller, 2019a). Among these model systems, each has advantages and disadvantages, for instance, budding yeast in a haploid state contain only one copy each of the major core histone genes (Eriksson et al., 2012). In haploid yeast, it is possible to alter a single amino acid to test the importance of a specific histone PTM, thus making yeast a particularly attractive model system for such studies.
Budding yeast, such as the well-studied strain Saccharomyces cerevisiae (S. cerevisiae), participate in both symmetric mitotic cell divisions during the budding process, and meiotic cell divisions during yeast sporulation (Neiman, 2011). All microorganisms, including yeast, spend a majority of their life-cycle in a quiescent state due to a lack of resources in their natural environment (De Virgilio, 2012;Sagot and Laporte, 2019b). Diploid yeast cells can enter a quiescent state in response to nutrient depletion, stress, and even cell wall damage . Saccharomyces cerevisiae initiate quiescence following the exhaustion of nutrients, especially glucose, and have been widely used to study quiescence. Quiescent yeast share some similarities to the quiescent state of mammals (Gray et al., 2004;Dhawan and Laxman, 2015;Miles et al., 2021). When grown in the laboratory, yeast consume glucose present in their growth medium and when available nutrients have been depleted, the yeast undergo a change termed diauxic shift (Chu and Barnes, 2016) as their metabolic profiles transition from fermentation to respiration, resulting in a decreased growth rate. When no other carbon sources are readily available, the yeast enter stationary phase or quiescence (Galdieri et al., 2010). By fractionating the cells based on differing densities, the non-proliferating stationary phase yeast have been separated into a denser population that is long-lived, and a less dense subpopulation that has been termed "non-quiescent" (Allen et al., 2006). Differences in the accumulation of trehalose and lipids may contribute to the different densities of these populations (Sagot and Laporte, 2019b). Fractionation protocols that purify quiescent yeast have allowed for comparisons of histone modifications in quiescent and proliferative yeast samples (Mews et al., 2014). In budding yeast that initiate quiescence following nutrient exhaustion, there is a genomewide shift in gene expression, with transcriptional repression of genes involved in growth and proliferation including ribosomal genes (Werner-Washburne et al., 1996;Gray et al., 2004;Radonjic et al., 2005;McKnight et al., 2015).
Yeast can also form spores, which can serve as another model for quiescence. When diploid yeast cells undergo meiosis, the meiotic products can differentiate into dormant spores during the process of sporogenesis (Greig, 2009;Duina et al., 2014). The state achieved in dormant microbial spores shares some similarities to the quiescent state achieved by nutrient depletion in S. cerevisiae (Greig, 2009;Duina et al., 2014). For instance, spore formation, like quiescence, is reversible as spores germinate to form haploid cells when exposed to nutrients. Spores are distinguished from a quiescent state because quiescent cells maintain some metabolic capacity, maintain membrane potential and do not undergo a morphological differentiation (Rittershaus et al., 2013). Spore formation is characterized by a dramatic decrease in global transcription levels (Xu et al., 2012). As described below, both nutrient limitation and spore formation have been used as models to probe the role of histone PTMs in quiescence in S. cerevisiae.
Fission yeast like Schizosaccharomyces pombe (S. pombe) are also an excellent model for quiescence (Su et al., 1996). During meiosis in fission yeast, asymmetric division takes place in which inheritance of the mother cell's components give rise to four different, unique daughter cells (Higuchi-Sanabria et al., 2014). Removing nitrogen from S. pombe causes the yeast to mate with yeast of the opposite mating type followed by replication through meiosis (Freese et al., 1982). However, if there is only one mating type of yeast in the population, the fission yeast arrest in G1phase and enter quiescence (Nurse and Bissett, 1981;. These nitrogen-deprived fission yeast can remain viable for months. Quiescent fission yeast cells are metabolically active, engage stress-responsive signaling and are efficient in DNA damage repair (Su et al., 1996;Mochida and Yanagida, 2006;Ben Hassine and Arcangioli, 2009;Marguerat et al., 2012;Gangloff and Arcangioli, 2017). When fission yeast enter a state of quiescence as non-dividing spores, a gene regulatory program is activated that includes upregulation of genes needed to adapt to the quiescent state (Su et al., 1996;Sajiki et al., 2009;Takeda and Yanagida, 2010), the nucleus undergoes changes in chromatin compaction, and histone modifications are altered (Neiman, 2011). For all of these reasons, fission yeast represents a valuable model system for studying epigenetic changes with quiescence.

Fibroblast Models of Quiescence
In multicellular organisms, there are multiple different types of quiescent cells, such as quiescent fibroblasts, immune cells, and stem cells, that serve as model systems for the study of quiescence at the molecular level (Mitra et al., 2018a). Fibroblasts, which are normally quiescent in vivo, contribute to the physical form and biomechanics of tissue by secreting growth factors and extracellular matrix proteins. Fibroblasts are organizers of the wound healing process as they can proliferate and replenish dead cells at the wound site and secrete extracellular matrix proteins that contribute to the formation of granulation tissue and scars (Lynch and Watt, 2018). Fibroblasts isolated from different tissues such as skin or lung are relatively easy to culture, and quiescence can easily be achieved by serum starvation, contact inhibition, or loss of adhesion (Coller et al., 2006;Mitra et al., 2018a). When fibroblasts enter a quiescent state, there is a dramatic change in gene expression in which a large fraction of the genome is differentially regulated (Coller et al., 2006;Suh et al., 2012;Mitra et al., 2018b). This change in gene expression is accompanied by significant changes in the abundance and activity of microRNAs (Suh et al., 2012;Johnson et al., 2017), transcript decay rates (Johnson et al., 2017;Mitra et al., 2018b), splicing (Mitra et al., 2018b), and the use of proximal versus distal polyadenylation sites (Mitra et al., 2018b). Fibroblasts are a heterogeneous population of cells. They can be isolated from different locations within the skin including hair follicles, and locations within the dermal layer, such as the papillary and reticular dermis (Sorrell and Caplan, 2004), and they differ based on their location within tissue (Sorrell and Caplan, 2004;Lynch and Watt, 2018). Further, fibroblasts isolated from skin from different anatomical sites have distinct and characteristic transcriptional programs that include extracellular matrix synthesis, lipid metabolism and signaling pathways (Chang et al., 2002). Single-cell sequencing data has shed light on the heterogeneity of fibroblasts (Muhl et al., 2020). While the tissue of origin represents an important contributing factor to the differences among fibroblasts, single cell sequencing has also shown intra-organ heterogeneity (Muhl et al., 2020). Different fibroblast subpopulations have distinct characteristics and the contributions of each fibroblast population to physiology is being actively elucidated (Muhl et al., 2020).

Stem Cell Models of Quiescence
Adult stem cells are another widely used quiescence model. Stem cells have been used to study quiescence in the context of the tissue-specific niche in which they are located. Many types of stem cells are largely quiescent unless activated to proliferate and differentiate in order to maintain tissue homeostasis and tissue regeneration (Li and Bhatia, 2011;Cheung and Rando, 2013;Coller, 2019b;Urbán and Cheung, 2021). Dysregulation or loss of stem cell quiescence can result in depletion of a stem cell pool, which can impede tissue regeneration (Cheung and Rando, 2013). Cellular quiescence has been studied in different types of adult stem cells such as hair follicle stem cells (HFSCs) within the skin (Lien et al., 2011;Lee et al., 2016;Rodriguez and Nguyen, 2018), hematopoietic stem cells (HSCs) from bone marrow (Nakamura-Ishizu et al., 2014), neural stem cells (NSCs) in the brain (Basak et al., 2018), and skeletal muscle stem cells (MuSCs) (Fukada et al., 2007). In the skin, the hair follicles that anchor hair to the skin progress through a cycle. During the anagen phase, there is rapid cell proliferation and growth of a new hair follicle. In the catagen phase that follows, the hair stops growing and detaches from the blood supply. Finally, in the telogen phase or resting phase, a new hair grows beneath the existing hair. Quiescent, non-proliferative hair follicle stem cells, which can be identified with cell surface markers including CD34 and CD49 (Garza et al., 2011), reside within a portion of the hair follicle called the bulge during the hair follicle's resting stage, the telogen phase. During the transition from telogen to anagen, HFSCs are activated, exit the bulge and proliferate downward, creating a trail of rapidly proliferating cells (Nowak et al., 2008). These proliferating cells terminally differentiate to give rise to cells that form the new hair shaft and its channel. Hair follicle stem cells have been an important model system for understanding quiescence, including the epigenetics of quiescence. Another important model system for understanding quiescence is the hematopoietic stem cell (HSC) compartment. HSCs are the stem cells that give rise to other blood cells including both myeloid and lymphoid lineages. In adult animals, hematopoiesis occurs in the bone marrow and the stem cells are only a small fraction of all of the cells present. The HSCs with the greatest capacity for self-renewal in the mouse bone marrow are quiescent HSCs (Wilson et al., 2008;Tesio et al., 2015) which are long-term label retaining and are in a deeply quiescent state (Foudi et al., 2009;Tesio et al., 2015). They are reported to divide only 5 times per lifetime (Foudi et al., 2009;Tesio et al., 2015). In response to infection or chemotherapy, these cells enter the cell cycle and start to proliferate to replenish damaged or lost cells (Wilson et al., 2008;Tesio et al., 2015). HSCs can be identified and isolated from surrounding cells by combinations of cell surface markers, including the presence of CD34 and an absence of markers for specific cell lineages (Sieburg et al., 2006;Dykstra et al., 2007;Kent et al., 2007).
Neural stem cells (NSCs) are present in the developing brain where they generate neurons (Ma et al., 2009). In adult animals, specific regions in the brain contain NSCs that have the capacity to proliferate and generate new neurons, thereby allowing adults to learn and acquire new skills, for instance, the ability to smell new odors (Ma et al., 2009;Obernier et al., 2018;Kalamakis et al., 2019). Like quiescent HFSCs, NSCs in vivo are thought to be slowly dividing and can be identified based on their label retention (Ma et al., 2009). These NSCs are depleted as organisms age (Obernier et al., 2018), and recent studies in single cells have shown an increase in quiescent NSCs in older mice compared with younger mice (Kalamakis et al., 2019). Markers for NSCs include expression of glial fibrillary acidic protein and glycoprotein CD133, along with an absence of differentiated cell markers (Ma et al., 2009).
Finally, muscle stem cells (MuSCs), or muscle satellite cells, represent another valuable quiescence model as MuSCs are mostly quiescent in uninjured tissue . Muscle stem cells are crucial for the process of regenerating skeletal muscle (Boonsanay et al., 2016). Upon muscle injury, quiescent satellite cells are activated to re-enter the cell cycle and proliferate (Yin et al., 2013). The proliferating progeny of the muscle satellite cells can differentiate into myotubes that form muscle fibers (Dumont et al., 2015;Boonsanay et al., 2016). Other MuSCs return to quiescence and are available to assist in the repair of subsequent muscle injury events (Yin et al., 2013). In vivo analysis of MuSCs has revealed that MuSCs are primed for activation (van Velthoven et al., 2017). The quiescent state of MuSCs has been reported to be an "idling" state for stem cells because widespread, low-level transcription was observed and hypothesized to serve as a means to ensure that the transcription machinery is ready to respond when required (van Velthoven et al., 2017).
Stem cells in these different models are normally found in a quiescent state that is maintained by specific niches that are tightly regulated by intrinsic and extrinsic factors (So and Cheung, 2018). Within these complex niches and diverse tissue compartments, quiescent stem cells are identified by their low RNA content and lack of proliferative markers (Fukada et al., 2007). In many cases, quiescent stem cells can be isolated from the tissue of interest by monitoring the presence and levels of cell-surface markers with fluorescence-activated cell sorting (FACS), and label-incorporation assays that identify cells that haven't divided for an extended time (Cheung and Rando, 2013;Nakamura-Ishizu et al., 2014). While this general pattern of quiescent stem cells is reported in multiple tissues, it is important to note that not all stem cells are quiescent, and in some tissues, stem cells are proliferative (Barker et al., 2007).

H3 METHYLATION WITH QUIESCENCE
Histone H3 has more potential methylation sites than any other histone in the histone octamer. Lysine residues 4, 9, 27, and 36 within the tail region of H3 are the most frequently methylated amino acids (Hyun et al., 2017;Jambhekar et al., 2019). In addition, H3K79, a lysine located within the globular domain of the histone, rather than on the tail, can also be methylated (Farooq et al., 2016). Each of the H3 lysines can be mono-(me), di-(me2), or tri-methylated (me3) by histone methyltransferases (writers) (Husmann and Gozani, 2019), and the methylation marks can be removed by histone demethylases (erasers). Histone methyltransferases, with the exception of DOT1L, contain a Suppressor of variegation, Enhancer of zeste, and Trithorax (SET) domain that catalyzes methylation of the lysine ε-amino group (DesJarlais and Tummino, 2016). Methylations of different histone lysines have characteristic deposition patterns that support a possible role for them in establishing different types of chromatin either by modulating the accessibility of the DNA to proteins, or by serving as a binding site for readers that act as effectors. Below we discuss the evidence that each of these marks plays a role individually and in combination in quiescence. Many of these marks have been investigated individually, but not in combination, and thus in many instances, our understanding of how they contribute to a histone code is limited.

H3K4 Methylation
H3K4me1/me2/me3 are usually associated with gene activation (Jambhekar et al., 2019). H3K4me1, me2, and me3 methylation marks are enriched in enhancers, the 5' ends of genes, and promoters, respectively (Hyun et al., 2017;Jambhekar et al., 2019). Methylation of H3K4 is associated with gene activation and the presence and absence of H3K4 methylation marks provides insight into the genome-wide patterns of active and inactive genes, respectively. Changes in genome-wide H3K4 methylation patterns with quiescence provide insights into gene activation and repression with quiescence. Whether these changes are functionally important for the changes in gene expression with quiescence, or are only correlated with quiescence, and whether H3K4 marks act alone or in combination with other marks or effector proteins, is the subject of active investigation in multiple model systems.
Young and colleagues investigated these questions in budding yeast S. cerevisiae by inducing the yeast into quiescence by glucose depletion (diauxic shift) and analyzing the yeast at 7 and 14 days. Overall levels of H3K4me2 were similar between log (proliferating cells) and stationary phase cells containing a mixture of quiescent and non-quiescent cells ( Table 3) (Young et al., 2017). Levels of H3K4me3 decreased about 50% in quiescent and non-quiescent stationary phase cells compared with proliferating cells in these studies (Young et al., 2017). In another study focused on cell cycle entry from quiescence, S. cerevisiae were maintained in a nutrient depleted environment and then restimulated with complete medium (Mews et al., 2014). In this study, Mews and colleagues found that overall levels of H3K4me1, H3K4me2 and H3K4me3 were similar in log and a mixture of quiescent and non-quiescent stationary phase cells (Mews et al., 2014).
While H3K4me3 levels were similar or changed by 50% in the two studies comparing proliferating and stationary phase budding yeast, there was an observed difference in transcription rate between these two populations of cells. Levels of RNA Pol II CTD residues phospho-Ser5 and phospho-Ser2, which are indicators of transcriptional initiation and elongation, respectively, were high in day 3 post-diauxic cells, and subsequently decreased in day 5 post-diauxic cells (Young et al., 2017). In both proliferating and quiescent cells, H3K4me3 and RNA polymerase II were found at gene promoters, but the distribution of the H3K4me3 mark among promoters shifted with quiescence (Young et al., 2017). H3K4me3 was more abundant at the promoters of growth-associated genes in log phase yeast (Mews et al., 2014). Genes that retained the H3K4me3 mark and RNA Pol II at their promoters in quiescent yeast included genes responsible for stress response, protein catabolism, and energy production (Mews et al., 2014;Young et al., 2017). Thus, an association was observed between the presence of the H4K3me3 histone mark and activation of genes with quiescence, but further studies would be needed to determine whether this constitutes a "code" and if so, what the functional consequences of this mark are for quiescence.
Further studies have evaluated the functional importance of H3K4 methylation for quiescence in strains of S. cerevisiae. S. cerevisiae were genetically engineered so that lysine 4 of histone H3 was mutated to alanine. These mutant strains lost reproductive capacity over time, after being introduced into stationary phase, to a greater extent than wild-type yeast (Walter et al., 2014;Young et al., 2017). Further, in this mutant strain, the proportion of non-quiescent cells in the stationary phase increased (Walter et al., 2014;Young et al., 2017). These results indicate that methylation of H3K4 is required for the establishment or maintenance of a quiescent state initiated in response to nutrient depletion.
The effects of removal of H3K4 methylation marks (demethylation) have also been characterized in S. cerevisiae during spore production. In S. cerevisiae spores, which have very low transcription rates, there is an accumulation of the highly conserved H3K4 demethylase JARID1-family histone 2 (JHD2) during sporulation (Xu et al., 2012) (Figure 1). The spores of JHD2 mutant yeast strain, jhd2 , show a ∼2-fold increase in H3K4me3 and reduced levels of H3K4me1/me2 compared to wild-type spores, suggesting that absence of JHD2 leads to the conversion of H3K4me1/me2 to H3K4me3 (Xu et al., 2012). Studies with wild-type and jhd2 mutant strains indicated that JHD2 demethylases reduce intergenic transcription induced by H3K4me3 during spore differentiation, promoted protein coding gene transcription, and repressed nucleosome accumulation at transcription start sites (TSSs) of a large subset of ribosomal protein-coding genes. Mutants of jhd2 exhibited precocious differentiation and the spores formed were sensitive to stress. Since JHD2 needs alpha-ketoglutarate from the TCA cycle for enzymatic function, it is possible that this family of proteins can sense carbon metabolism activity, which in turn regulates transcription in response to nutrient availability (Xu et al., 2012). These findings, taken together, suggest that H3K4me3 may be a mark that can transform information about nutrient availability into a complex pattern of gene expression that has distinct effects on protein coding genes, ribosomal genes and non-coding RNAs.
The role of H3K4 methylation has been studied in the context of quiescent mammalian cells as well. Kallingappa and colleagues compared chromatin states of proliferating bovine adult ear skin fibroblasts in the G1 phase of the cell cycle with chromatin composition of serum-starved, quiescent (G 0 ) fibroblasts (Kallingappa et al., 2016). With fluorescence microscopy, quiescent nuclei were found to have a more relaxed or less compact chromatin state, and half the levels of H3K4me3 compared to G1 nuclei (Kallingappa et al., 2016). In mouse B cells, there was a dramatic shift in histone H3K4 methylation with quiescence, with much lower H3K4me2 levels in quiescent mouse B cells compared to cycling mouse B cells (Baxter et al., 2004). Reduction of H3K4me3 levels was also observed in quiescent mouse HSCs in late catagen stage in comparison to proliferating HSCs in early anagen stage (Lee et al., 2016;Kang et al., 2020). Based on ChIP-seq, the H3K4me3 signal generally decreased Decreased in sporulating yeast (Xu et al., 2012), mouse skin and HFSCs (Lee et al., 2016), and mouse B lymphocytes (Baxter et al., 2004) H3K4me3 Euchromatin Gene activation S. cerevisiae S. pombe Increase in adult mouse muscle stem cells (Liu et al., 2013) Decrease in quiescent yeast (Xu et al., 2012;Young et al., 2017), sporulating yeast (Xu et al., 2012), bovine fibroblasts (Kallingappa et al., 2016), and in HFSCs (Kang et al., 2020) H3K9me2 Heterochromatin Gene repression S. pombe increased in muscle quiescent stem cells (Cheedipudi et al., 2015) reduced in yeast (Oya et al., 2019) and bovine fibroblasts (Kallingappa et al., 2016) H3K9me3 Heterochromatin Gene repression S. pombe increased in MuSCs (Boonsanay et al., 2016) reduced in mouse skin and HFSCs (Lee et al., 2016), yeast (Oya et al., 2019), and bovine fibroblasts (Kallingappa et al., 2016) H3K9ac Euchromatin Gene activation S. cerevisiae S. pombe Global decrease in yeast (Mews et al., 2014) H3K14ac Gene activation S. cerevisiae S. pombe Decreased in muscle quiescent stem cells and yeast (Cheedipudi et al., 2015;Young et al., 2017) H3K27me3 Heterochromatin Gene repression Decrease in skin and hair stem cells for growth (Kang et al., 2020), in mouse skin and HFSCs (Lee et al., 2016), and in bovine fibroblasts (Kallingappa et al., 2016) H3K27ac Euchromatin Gene activation S. cerevisiae S. pombe Increase in neural stem cells (Martynoga et al., 2013) H3K36me1/me2/me3 Euchromatin Gene activation S. cerevisiae S. pombe me3 increase in mouse bone marrow cells  reduced in bovine fibroblasts (Meng et al., 2020) H3K79me1/me2 Heterochromatin S. cerevisiae S. pombe loss in yeast (Young et al., 2017) H4K5ac Euchromatin Gene activation S. cerevisiae S. pombe Global decrease in yeast (Mews et al., 2014) H4K16ac Euchromatin Gene activation S. cerevisiae S. pombe decreased in skeletal muscle stem cells and ESCs (Ryall et al., 2015;Kang et al., 2020) H4K20me3 Heterochromatin Gene repression S. pombe increased in primary human dermal fibroblasts (Evertts et al., 2013a), in MuSCs (Boonsanay et al., 2016) at transcription start sites and overall in quiescent compared with proliferating HSCs, although there was no clear correlation between H3K4me3 marks and gene expression changes with proliferation or quiescence (Lee et al., 2016). These findings suggest that H3K4 methylation may have other roles in addition to transcriptional regulation. Indeed, a recent article re-evaluating the role of histone-modifying enzymes argues that histone H3K4 methylation only has a minor role in transcriptional regulation (Rickels et al., 2017;Morgan and Shilatifard, 2020), but that H3K4 methylation has been associated with DNA recombination, repair, and replication (Daniel and Nussenzweig, 2012;Acquaviva et al., 2013;Kantidakis et al., 2016). Similar results were observed in MuSCs where neither the number nor the identify of genes marked by H3K4me3 at their transcription start site changed in activated compared with quiescent stem cells, and quiescence-specific genes retained their H3K4me3 mark at their transcription start site even upon activation (Liu et al., 2013). The authors conclude that H3K4me3 marks genes for transcriptional activation but its presence is not sufficient to determine whether a gene will be expressed (Liu et al., 2013).
Evidence for the importance of H3K4 demethylation as a regulator of quiescence is derived from studies of Retinoblastoma binding protein 2 (RBP2), which can demethylate all the methylation states of H3K4 in vivo (Klose et al., 2007). In vitro, RBP2 can processively remove the me3 and me2 marks on H3K4 to return the histone to a singly methylated state, but it cannot demethylate the me1 mark (Klose et al., 2007). Hematopoietic stem cells (HSCs) and myeloid progenitors isolated from Rbp2 knockout mice (Rbp2 −/− ) contain a higher proportion of cells exiting quiescence compared to WT. Further, in Rbp2 −/− cells, genes encoding cytokines were marked with higher levels of H3K4me3 and were expressed at higher levels, which could promote proliferation of the HSCs (Klose et al., 2007). In addition to the JmjC domain that is responsible for the demethylase activity, KDM5A/RBP2 also contains 2-3 plant homeodomain (PHD) domains (Klose et al., 2007). Binding of the PHD domain to unmodified H3 peptide activates the KDM5A/RBP2 catalytic activity and results in removal of methyl marks from histone H3K4me3 from a nearby nucleosome (Torres et al., 2015). By coupling the KDM5A/RBP2's ability to read unmodified H3K4 with demethylation of nearby H3K4me3, results in a positive feedback loop that allows the spreading of a chromatin state of demethylated histones (Torres et al., 2015).
The data that are currently available suggest that methylation of H3K4 may play a role in aspects of quiescence, including potentially transcription, but whether this mark is a determinant of gene expression or simply associated with an activated promoter is not yet clear. H3K4 methylation could also contribute to quiescence through its roles in replication, repair or recombination. Establishing whether H3K4 methylation affects quiescence entry or maintenance through direct effects on the biomechanical properties of the DNA or through readers, whether it is part of a histone code, and how this affects functional aspects of quiescence will require additional studies.

H3K9 Methylation
H3K9 methylation marks are primarily associated with gene repression in heterochromatin regions (Saksouk et al., 2015). H3K9me3 associates with regions of constitutive heterochromatin such as repeat regions of telomeres and centromeres (Saksouk et al., 2015). H3K9me3 is also deposited at some genomic regions in a tissue-specific manner and plays a role in cell identity (Ninova et al., 2019). Given the observation that in some systems, chromatin is more compact in quiescent than proliferating cells (Bridger et al., 2000;Evertts et al., 2013a;Guidi et al., 2015;Criscione et al., 2016;Swygert et al., 2019Swygert et al., , 2021, H3K9 methylation is of particular interest as a potential regulator of chromatin compaction with quiescence. In fission yeast, there is only one H3K9 methyltransferase, Clr4/Suv39H, which adds H3K9me1, me2, and me3 marks (Figure 1, Table 4). Inactivation of Clr4 resulted in viability similar to wild-type cells when nutrients were present, but reduced viability when quiescence was initiated following nitrogen starvation or glucose deprivation (Joh et al., 2016). This trend was also observed in yeast cells with a mutation that converted histone H3 lysine 9 to an alanine, suggesting that H3K9 methylation is important for survival during quiescence (Joh et al., 2016). Global ChIP-seq analysis showed that as fission yeast enter quiescence, the cells accumulate Clr4-dependent H3K9me2 and H3K9me3 marks at euchromatic genes whose transcriptional regulation has been shown to be important for establishing quiescence including genes involved in metabolism, ribosomal genes, cell cycle genes and stress response genes (Joh et al., 2016). In this study, a strong correlation was observed between genes with H3K9me2 marks and the set of genes repressed in G 0 (Joh et al., 2016). The enrichment of H3K9 methylation marks in euchromatic gene regions upon quiescence entry required the small RNAs (sRNAs) associated with RNA interference (RNAi) factor Argonaute (Ago1). Quiescent yeast had a distinct profile of these sRNAs (Joh et al., 2016). Before H3K9 is deposited in these euchromatic regions in quiescent yeast, Ago1-associated sRNAs are expressed from these regions (Joh et al., 2016), and these sRNAs may serve as guides for the deposition of H3K9 methylation marks (Joh et al., 2016). This may reflect a mechanism for regulating the expression of specific genes as yeast enter quiescence. In contrast to euchromatic regions, the levels of H3K9me2 in constitutive heterochromatic regions decline early during quiescence (8h and 24h after starvation) (Oya et al., 2019). These findings support the possibility of a combined histone-RNA code that controls gene expression and viability during quiescence. Surprisingly, the authors observed relatively little overlap between the H3K9me3-enriched genes in quiescent cells and the genes that are repressed by Clr4 in quiescence based on RNA-seq analysis of wild-type and clr4deleted yeast strains (Joh et al., 2016). This disconnect between the two genesets shows that additional studies will be needed to clearly determine whether H3K9 methylation plays a role in regulating transcription with quiescence, whether transcriptional regulation by H3K9 is crucial for viability in the quiescent state, and whether the effects of H3K9 in conjunction with short RNAs constitute part of a quiescence histone code. Similar to the findings in fission yeast, in hair follicles, quiescent HFSCs that were isolated from the late catagen stage of the hair follicle cycle contained considerably lower levels of the H3K9me3 mark compared to proliferating HFSCs isolated from the early anagen stage of the hair follicle cycle (Lee et al., 2016;Kang et al., 2020) (Table 3).
In contrast to other studies, Boonsanay et al. found no change in the levels of H3K9me3 between proliferating and quiescent mouse MuSCs (Boonsanay et al., 2016). This study did not measure the global distribution of H3K9 methyl marks, and thus the genomic regions where H3K9 methyl marks are found or how they are redistributed in proliferating and quiescent MuSCs is not known (Boonsanay et al., 2016).
Studies of fibroblasts have also revealed differences in the levels of H3K9 methyl marks with quiescence although the changes observed varies in fibroblasts from different sources. In adult ear skin fibroblasts, overall levels of H3K9me2 and me3 were roughly halved in quiescent cells compared to the same fibroblasts in the G1 phase of the cell cycle, while the levels of H3K9me1 were slightly elevated in quiescent relative to G1 cells (Kallingappa et al., 2016). In contrast, levels of H3K9me3 were modestly elevated in quiescent human dermal fibroblasts compared to proliferating fibroblasts (Evertts et al., 2013a), while levels of H3K9me2 and H3K9me1 were similar (Evertts et al., 2013a). In mouse B lymphocytes, the levels of H3K9me2 and me3 were lower in quiescent cells compared to activated cells (Baxter et al., 2004). Thus, in different quiescence model systems, changes in the levels of H3K9 methylation states are altered, but the specific changes reported have been different depending on the model system, and thus a consistent quiescence program of H3K9 methylation changes has not been observed.
In mammals, there are multiple methyltransferases that add methyl groups to histone H3 lysine 9 including the 17 members of the PRDM family of proteins (Steele-Perkins et al., 2001;Hyun et al., 2017). One of these members, PR domain-containing-2/Rb interacting zinc finger protein (PRDM2/RIZ) is expressed at high levels in quiescent mouse MuSCs in vivo (Cheedipudi et al., 2015) (Table 4). Knockdown and overexpression studies of PRDM2/RIZ indicated that in quiescent MuSCs, PRDM2/RIZ prevents lineage commitment and irreversible cell cycle arrest (Cheedipudi et al., 2015). Global analysis using ChIP coupled with DNA microarray (ChIP-Chip) showed that PRDM2 was associated with >4400 gene promoters in quiescent muscle cells that initiate quiescence in suspension culture, that is, loss of adhesion. Approximately 50% of these promoters were also marked with H3K9me2 (Cheedipudi et al., 2015). The PRDM2associated genes were enriched for differentiation, cell cycle, and developmental regulators (Cheedipudi et al., 2015). The levels of H3K9 methylation marks (me1, me2, and me3) did not change overall upon knockdown of PRDM2 (Cheedipudi et al., 2015). However, H3K9me2 levels were reduced at the MyoG promoter in G 0 cells, while PRDM2 knockdown cells showed reduced H3K9me2 at the same locus, suggesting that PRDM2 may add methyl groups to generate H3K9me2 at MyoG (Cheedipudi et al., 2015). Further, increased H3K14Ac was observed at the MyoG promoter upon PRDM2 knockdown, supporting a role for PRDM2 in regulating the expression of MyoG, a critical molecule for muscle differentiation (Cheedipudi et al., 2015). These findings support a model in which PRDM2 activation in G 0 ensures that two distinct possible outcomes-myogenesis and cell cycle progression-are poised for reactivation. The findings implicate H3K9 methylation, possibly in combination with other histone marks, in muscle stem cell quiescence and renewal (Cheedipudi et al., 2015).
Thus, while quiescence is generally associated with a reduction in transcription and number of active genes based on H3K4 methylation status as well as other markers, surprisingly in most, but not all, studies, there was not an increase in the H3K9me2 and H3K9me3 marks that might be expected. H3K9me3 is associated with constitutive heterochromatin (Saksouk et al., 2015) and it is possible that the increase in heterochromatin with quiescence reflects an increase in facultative heterochromatin. Nevertheless, even though H3K9me2/3 are not consistently observed to increase with quiescence, in yeast, there is evidence that these marks may be functionally important for the viability of quiescent cells as perturbations that reduce their levels reduce viability of quiescent cells with little effect on cells in full nutrient conditions. The data taken together support a possible role for H3K9 marks as contributors to changes that ensure the viability of quiescent cells, but these effects may be mediated through mechanisms other than transcriptional changes.

H3K27 Methylation
The H3K27me3 mark is found in multiple model organisms including Arabidopsis, Drosophila, worms, and the filamentous fungus Neurospora crassa, but not yeast (Jamieson et al., 2013). H3K27me3 is a reversible mark of facultative heterochromatin, chromatin that can become compact or open depending upon the circumstance and is not repetitive (Trojer and Reinberg, 2007). H3K27me3 decorates genes that are developmentally regulated and are switched on and off depending upon the stage of development (Jambhekar et al., 2019). The Enhancer of Zeste 1 (EZH1) and EZH2 histone lysine methyltransferases trimethylate H3K27, and are constituents of Polycomb Repressive Complex 2 (PRC2) (Margueron and Reinberg, 2011) (Figures 1, 2,  Table 4), while PRC1 stabilizes PRC2 binding to H3K27 and catalyzes monoubiquitination of histone H2A lysine 119 (Margueron and Reinberg, 2011). This ubiquitination mark represses transcription (Lavarone et al., 2019) and promotes chromatin compaction (Wiles and Selker, 2017;Lavarone et al., 2019). The Embryonic Ectoderm Development (EED) protein subunit of the PRC2 complex binds to H3K27me3 and this interaction increases the methyltransferase activity of the complex, resulting in a positive feedback loop that establishes zones of repressed chromatin (Margueron et al., 2009;DesJarlais and Tummino, 2016).
In some studies, H3K27me3 levels were reduced with quiescence. In bovine fibroblasts, which exhibit more open chromatin with quiescence, H3K27me3 levels were reduced by approximately a half when the fibroblasts entered quiescence in response to serum starvation (Kallingappa et al., 2016). Similarly, protein levels of the members of PRC2-EZH2, EED, and Suppressor of Zeste 12 (SUZ12)-were reduced by about half in G 0 nuclei (Kallingappa et al., 2016). Proteins in the PRC1 polycomb complex were also reduced in quiescent bovine fibroblasts, specifically Polyhomeotic Homolog (PHC1) and Ring Finger Protein 2 (RING2) (Kallingappa et al., 2016). H3K27me3 levels are also reduced in quiescent mouse B lymphocytes compared to activated and cycling B lymphocytes (Baxter et al., 2004).
In contrast, murine chondrocytes (cartilage cells) that entered into a quiescent state by the application of physiological hydrostatic pressure had lower levels of H3K9me3 and higher levels of H3K27me3 compared to control cells that were not subjected to hydraulic pressure (Maki et al., 2021). In primary human dermal fibroblasts, a mass spectrometry-based analysis of histone post-translational modifications revealed higher levels of H3K27me3 in contact-inhibited, quiescent fibroblasts than proliferating fibroblasts (Evertts et al., 2013a). Further dissection of the H3K27me3 signal revealed that levels were particularly high when H3K27me2 was found in combination with H3K36me2 or if H3K27me3 was found on the same histone tail as H3K36me1 and H3K36me2 (Evertts et al., 2013a). This pattern of histone modifications may reflect the fact that H3K36me3 can inhibit the ability of PRC2 to methylate H3K27 (Schmitges et al., 2011;Yuan et al., 2011).
H3K27me3 levels have been observed to be reduced in multiple quiescence model systems involving stem cells. Quiescent mouse HFSCs isolated in the late catagen (quiescent) stage experienced a global decrease in the levels of H3K27me3 along with a decrease in H3K4me3 and H3K9me3, as mentioned in the sections above, compared to cycling HFSCs in earlyanagen stage. This hypomethylation of quiescent HFSCs was confirmed using immunofluorescence, western blotting, and ChIP-seq methods (Lee et al., 2016). The levels of H3K27me3 were reduced in 64% of promoters with quiescence in HFSCs (Lee et al., 2016). Surprisingly, the changes in the levels of histone methylation marks, overall, did not correlate with changes in the levels of transcripts for the associated genes in quiescent versus proliferating hair follicle cells (Lee et al., 2016). There were exceptions: a larger than expected fraction of genes highly expressed at all hair cycle stages in the bulge had an increase in the levels of H3K27me3 and a set of genes defined as cell cycle regulators and tumor suppressors had almost no H3K27me3 in quiescent and proliferating HFSCs (Lee et al., 2016). To explore a possible role for decreased histone methylation levels during catagen, the bone morphogenetic protein (BMP) signal which normally maintains quiescence in vivo was inhibited, resulting in elevated levels of H3 K4/K9/K27 me3 in quiescent bulge HFSCs (Lee et al., 2016). These findings demonstrate that quiescent HFSCs in the bulge require active BMP signaling in order to maintain a hypomethylated H3 state (Lee et al., 2016). When keratinocytes, skin epithelial cells, were serum-starved in culture, transcript levels of several histone methyltransferases including EZH2 (forms H3K27me3), SUV39H1 and SUV39H2 (H3K9me3), decreased (Lee et al., 2016), while transcript levels of multiple histone demethylases increased (Lee et al., 2016). The latter include JMJD2a which catalyzes demethylation of histone H3 lysines 9 and 36 (Kim T.D. et al., 2012), UTX which demethylates H3K27me3 (Tang et al., 2017), and JARID1, which acts as a demethylase for H3K4me3 and H3K4me2 (Lee et al., 2016). After chemically inhibiting demethylases specific to the K4, K9, and K27 me3 marks with a cocktail applied to the mouse's skin, cells at the catagen stage failed to generate new hair follicles in the following hair cycle. Thus, the reduction in the H3K4/K9/K27 me3 levels observed in quiescent HFSCs was necessary for the ability of quiescent cells to re-enter the cell cycle (Lee et al., 2016). While these data support the importance of histone demethylation for the quiescent state of HFSCs, it remains unclear whether one, two or all of these marks was required for a functional state, whether these marks maintain chromatin structure or serve as binding sites for effectors, and whether the absence of the methylation marks affected the activity of other marks such as acetylation marks.
Three different studies have investigated the role of H3K27 methylation in quiescent mouse MuSCs (Liu et al., 2013;Cheedipudi et al., 2015;Boonsanay et al., 2016). Liu et al. found that H3K27me3 levels were low in quiescent MuSCs and dramatically increased in activated stem cells (Liu et al., 2013) ( Table 3). The transcription start sites of genes expressed at high levels in quiescent stem cells were marked with H3K4me3, but not H3K27me3 (Liu et al., 2013). The 2,019 genes that were marked by H3K27me3 at their transcription start site in quiescent stem cells displayed very low expression levels (Liu et al., 2013). Upon activation, there was a dramatic increase in H3K27me3 in the gene body and intergenic regions (Liu et al., 2013). Boonsanay and colleagues, in contrast, reported that the levels of H3K27me3 were not different between quiescent and proliferating MuSCs (Boonsanay et al., 2016).
Both Boonsanay et al. and Liu et al. discovered changes in the levels of writers and erasers of H3K27me3 as quiescent mouse MuSCs were activated. Liu and colleagues reported that the increase in H3K27me3 marks with activation was associated with higher transcript levels of Ezh2 and lower levels of the demethylase Jmjd3 (Liu et al., 2013). Boonsanay and colleagues also found higher levels of PRC2-Ezh2 in proliferating MuSCs and higher levels of PRC2-Ezh1 in the quiescent MuSCs (Boonsanay et al., 2016). These findings would be consistent with a study by Margueron and colleagues that showed that Ezh2 is more closely associated with proliferation, while Ezh1 is more abundant in non-proliferative adult tissues (Margueron et al., 2008). Margueron and colleagues discovered that PRC2-Ezh2 effectively catalyzes the formation of H3K27me2/3, while PRC2-Ezh1 directly represses transcription and compacts chromatin. Additional studies would be needed to determine whether the differences in expression of these methyltransferases results in changes in the genome-wide distribution of H3K27me2/3, whether there is a direct effect on transcription or chromosome compaction, or whether the shift in the relative abundance of these methyltransferases is functionally important for establishing, maintaining or reversing quiescence in MuSCs.
Cheedipudi and colleagues focused on H3K27me3 levels in MuSCs in the context of their co-existence with H3K4me3 and H3K9me2 and their co-regulation by the H3K9 methyltransferase PRDM2/RIZ. As described in the H3K9 methylation section above, PRDM2/RIZ, an H3K9 methyltransferase, is enriched in quiescent muscle cells in vitro, where it participates in stalling differentiation and cell cycle programs, while also maintaining genes involved in differentiation and proliferation poised for future activation (Cheedipudi et al., 2015). Cheedipudi and colleagues found that knocking down PRDM2 resulted in a reduction of H3K4me3 and H3K9me2, and higher levels of H3K27me3 at the cyclin A2 promoter in G 0 MuSCs (Cheedipudi et al., 2015). They conclude that PRDM2 may block the deposition of H3K27me3 silencing marks at the cyclin A2 promoter in G 0 , thereby preserving the gene's potential for reactivation (Cheedipudi et al., 2015). These studies highlight the combinatorial nature of the histone marks during the transition between proliferation and quiescence.
Thus, taken together, substantial changes in H3K27me3 have been observed in different models of quiescence. H3K27me3 may contribute to transcriptional repression in quiescent cells and may contribute to a poised quiescent state, characterized by an ability to re-enter the cell cycle. However, additional studies that probe the genomic location of H3K27me3 in proliferating and quiescent cells in different model system and their functional importance for quiescence will be needed to fully understand the specific role of H3K27 methylation individually or in combination with other marks in quiescence.

H3K36 Methylation
H3K36 methylation marks are usually deposited across the entire gene body (Hyun et al., 2017;Jambhekar et al., 2019), and are associated with transcriptional activation, dosage compensation, transcriptional repression, and DNA repair (Wagner and Carpenter, 2012). In the study by Young and colleagues described above, quiescent budding yeast S. cerevisiae had similar levels of H3K36me2 and H3K36me3 marks compared to proliferating cells (Young et al., 2017). However, the levels of SET Domain containing 2 (Set2) which deposits all three (mono, di, and tri) methylation marks in budding and fission yeast (Morris et al., 2005;Wagner and Carpenter, 2012), were reduced as the cells entered a quiescent state. These findings indicated that H3K36me3 marks were deposited prior to quiescence entry, concomitant with the activity of RNA Pol II.
While yeast has only one H3K36 methyltransferase, humans have eight, including SETD2, the human ortholog of yeast Set2, that generates the H3K36me3 mark in vivo, and enzymes of the Nuclear Receptor-binding Set Domain (NSD) family (NSD1/2/3) that deposit H3K36me1 and H3K36me2 marks (Wagner and Carpenter, 2012) (Figure 2 and Table 4). The trimethylation of H3K36 by SETD2 is highly efficient on an unmethylated H3K36 compared to a H3K36me2 substrate (Husmann and Gozani, 2019). In Setd2 conditional knockout mice with HSCspecific Setd2 inactivation, there is a reduction of the number of HSCs in a quiescent state . Setd2 knockout HSCs had a reduced G 0 fraction and increased G1 and S/G2/M phases of the cell cycle . In addition, knockout HSCs also exhibited higher levels of apoptosis, reduced stem cell identity, increased differentiation toward progenitors and reduced multiple-lineage terminal differentiation potential . Setd2 knockout mice had mild bone marrow fibrosis, increased erythroid progenitors, but a decreased population of other bone marrow progenitors. As a result, the HSCs from Setd2 conditional knockout mice were less able to repopulate the hematopoietic system upon transplantation . Setd2 knockout HSCs showed significant reduction of H3K36me3, increased levels of H3K36me1/me2, and increased Nsd1/2/3 at both transcript and protein levels . The levels of H3K4me3, H3K79me2, and phosphorylation of Ser2 residue of RNA pol II, all of which are associated with transcriptional elongation, increased in Setd2 knockout cells. Based on these studies, the authors proposed a model in which, in the absence of SETD2, NSD proteins promote the phosphorylation and elongation of RNA pol II on specific genes, leading to a loss of quiescence . The findings support an important role for SETD2 in maintaining the quiescent state of HSCs . It is unclear whether SETD2 plays this role through H3K36me3 or non-enzymatic functions such as by affecting cryptic transcription (Carvalho et al., 2013) or alternative splicing (Bhattacharya et al., 2021).

H3K79 Methylation
In contrast to the previously discussed histone marks that are located on the histone H3 tail, the H3K79 methylation mark is located in the globular domain of H3 (Farooq et al., 2016). H3K79me2 and H3K79me3 are mainly found in the bodies of active genes and are associated with transcription elongation (Mueller et al., 2007). H3K79me2/3 can also maintain enhancerpromoter interactions at a subset of enhancers (Godfrey et al., 2019). H3K79 methylation has also been implicated in telomere silencing (Singer et al., 1998), recombination, DNA repair and cell cycle progression (Nguyen and Zhang, 2011). Monomethylation and trimethylation of H3K79 has been associated with gene activation and gene repression, respectively, in some studies (Barski et al., 2007). Disruptor of Telomere Silencing-Dot1 in yeast and DOT1L in humans-is the only enzyme responsible for methylation of H3K79 in S. cerevisiae, Drosophila, and humans as knockout of DOT1 in these organisms results in a loss of H3K79 methylation (van Leeuwen et al., 2002;Lee S. et al., 2018) (Figure 2 and Table 4). Both Dot1 and DOTL1 can catalyze mono-, di-and trimethylation (Frederiks et al., 2008). Dot1 is the only known non-SET domain-containing methyltransferase (Farooq et al., 2016). S. cerevisiae Dot1 does not methylate free histones, only histones in chromatin, in contrast to other histone methyltransferases (Lacoste et al., 2002;Lee S. et al., 2018).
Yeast Dot1 also has histone chaperone activity and is particularly important for nucleosome dynamics and chromatin accessibility on transcribed regions of long genes (Lee S. et al., 2018).
After diauxic shift in S. cerevisiae, despite the general shut down of transcription, both the quiescent and non-quiescent yeast populations contained higher levels of the H3K79me3 mark than proliferating yeast cells (Young et al., 2017). However, the quiescent population had reduced H3K79me1 and H3K79me2 levels compared with the non-quiescent population (Young et al., 2017). Levels of Dot1 decreased by day 3 post-starvation and remained low throughout the time course that ended at day 7 (Young et al., 2017). The elevated levels of H3K79me3 and reduced levels of H3K79me1/2, while Dot1 was reduced could reflect that H3K79me3-containing nucleosomes were not turned over, or that demethylases were activated (Young et al., 2017). H3K79me3 was enriched in gene bodies in transcripts expressed specifically in growing cells or specifically in quiescent cells, with no redistribution to non-canonical locations in genes or intergenic regions (Young et al., 2017). Log and quiescent cells contained similar numbers of gene binding sites for H3K79me3 (Young et al., 2017). These marks were established soon after the shift to diauxic growth and then retained in quiescent cells even as transcription was reduced (Young et al., 2017). There was not a strong correlation between RNA polymerase II occupancy and H3K79me3 marks on specific genes in proliferating or quiescent cells (Young et al., 2017).
Mutant S. cerevisiae strains that are no longer able to methylate histone H3K4 or ubiquitinate histone H2B showed a shorter chronological lifespan, indicative of reduced ability to re-enter the cell cycle (Young et al., 2017). In contrast, a higher proportion of yeast entered quiescence upon glucose deprivation in H3K79 mutant, and the yeast with H3K79 mutations showed enhanced ability to re-enter the cell cycle (Young et al., 2017). The findings suggest that the presence of a lysine that can be methylated at H3K79 makes cells less able to re-enter the cell cycle after glucose deprivation, and surprisingly, makes the quiescent yeast less fit (Young et al., 2017).
The H3K79 mark has also been implicated in the quiescence of mouse HSCs. In a study described above, knockout of Setd2, a histone methyltransferase involved in the addition of the H3K36me3 mark, led to the loss of bone marrow reconstitution after transplantation, with the mouse HSCs exhibiting loss of quiescence, increased apoptosis, and inability to differentiate into multiple lineages . When studying the changes at the histone level, it was found that knockout of Setd2 led to increased expression of NSD1/2/3, which are also H3K36 methyltransferases that generate H3K36me1 and H3K36me2 marks . Loss-of-function in SetD2, combined with gain-of-function NSD1/2/3, led to an increase in the H3K79me2 mark through recruitment of DOTL1, the H3K79 methyltransferase . Together with increased recruitment of histone acetylase Brd4, which recruits the Super Elongation Complex (SEC), DOT1L enhanced RNA Polymerase II elongation and expression of target genes that promote apoptosis and differentiation of quiescent cells . This resulted in increased expression of genes including Myc, as inhibiting Brd4 or Dot1l reduced markers of transcription elongation and Myc expression . The findings support a model where multiple histone-modifying enzymes are co-regulated to affect the choice between quiescence and differentiation .
These studies, taken together, suggest H3K79 methylation may play a role in quiescence through mechanisms that involve transcription initiation or elongation, or alternative processes. In yeast, this mark reduces the ability of nutrient-starved cells to proliferate when provided nutrients. In mammalian HSCs, data support H3K79 methylation as part of a complex regulatory code that affects transcription elongation and the fate of quiescent HSCs.

H4 METHYLATION WITH QUIESCENCE
While the N-terminus of H3 tails exit the nucleosome near DNA entry and exit sites, the tails of histone H4 (residues 1-20) stick out from the nucleosome's face on either side (Ghoneim et al., 2021). DNA breathing dynamics are altered upon removal of either the H3 or H4 tail, suggesting that there is cross-talk between the H3 and H4 tails (Ghoneim et al., 2021). The region encompassing residues 16-23 of H4 contains a high density of arginine and lysine residues that constitutes a basic patch (Ghoneim et al., 2021). In crystal structures, this H4 basic patch interacts with DNA or the cluster of acidic residues called the H2A/H2B acidic patch on the same or nearby nucleosomes (Ghoneim et al., 2021). Histone H4 is mainly methylated at residue K20 in mammalian cells. The methylation of K20 requires the acidic patch, as H4K20 monomethylation is not detected in nucleosomes with a defective acidic patch (Ho et al., 2021). H4K20 methylation is conserved from fission yeast S. pombe to humans (Jorgensen et al., 2013). Histone H4K20me1 has been observed in budding yeast S. cerevisiae (Edwards et al., 2011).
A ChIP-seq study in murine ESCs found that H4K20me3 colocalizes with transcriptionally active marks like H3K4me3 and H3K36me3 (Xu and Kidder, 2018). This bivalent placement of H4K20me3 with activating histone marks suggest a potential combinatorial histone code for a "poised" state. These "poised" states have been suggested to allow for rapid activation of RNA pol II (Bernstein et al., 2006) during cell state transitions. Sims and colleagues report on another combinatorial code involving H4K20, but in conjunction with H3K9 (Sims et al., 2006). They found that H4K20me3 and H3K9me3 were both enriched in pericentric heterochromatin. H4K20me1 and H4K20me2 were found together with H3K9me1 and H3K9me2, respectively, but in different regions on chromosome arms (Sims et al., 2006). The authors further found that H4K20me1 and H3K9me1 were enriched in the same nucleosome (Sims et al., 2006). What role these bivalent marks play in quiescence is not known and requires additional exploration.
H4K20 methylation status is cell cycle dependent (Pesavento et al., 2008;Oda et al., 2009;Abbas et al., 2010;Centore et al., 2010;Adikes et al., 2020), at least in part because PR-Set-7, the enzyme that generates H4K20me1 (Fang et al., 2002;Nishioka et al., 2002b;Beck et al., 2012), is actively targeted for proteasome-mediated degradation in S phase (Rice et al., 2002;Julien and Herr, 2004;Jorgensen et al., 2007;Abbas et al., 2010;Tardat et al., 2010;Wu et al., 2010). Consistent with these previous reports, in primary human fibroblasts, unmodified H4K20 was most abundant in S phase, less abundant in G2/M cells, and present at low levels in G1 and cells that entered quiescence by contact inhibition for 14 days (Evertts et al., 2013a). H4K20me1 was present at low levels in S phase, increased in G2/M, and increased further in G1 (Evertts et al., 2013a). H4K20me3 represented a small fraction of all of the histones, and was strongly induced with quiescence compared with all other cell cycle phases (Evertts et al., 2013a) ( Table 3).
Boonsanay and colleagues investigated the levels of different histone marks in quiescent mouse MuSCs, proliferating MuSCs and differentiated myotubes (Boonsanay et al., 2016). Quiescent MuSCs contained a lower level of H4K20me1, but comparable levels of H4K20me2 to that present in proliferating MuSCs. The levels of H4K20me3 were in the order: differentiated myotubes > quiescent MuSCs > proliferating MuSCs (Boonsanay et al., 2016) (Table 3). In contrast, H3K9me3 and H3K27me3 and euchromatin marker H3K9ac did not change between proliferating and quiescent MuSCs (Boonsanay et al., 2016). Without ChIP-seq data, the distribution of histone H4K20 methylation marks, and the extent of overlap with other histone marks in these different cells is not known.
In mammals, H4K20me1 is catalyzed by PR-SET-7 (Nishioka et al., 2002a) and the methyltransferases KMT5B/Suv4-20h1 and KMT5C/Suv4-20h2, and possibly KMT3E/SMYD3, generate H4K20me2 and H4K20me3 (Rice et al., 2002;Schotta et al., 2004Schotta et al., , 2008 (Figures 1, 2, Table 4). Loss of both SUV4-20H enzymes leads to strongly elevated levels of H4K20me1 (Schotta et al., 2008). Electron micrographs of MuSC nuclei with inactivation of Suv4-20h1 revealed a reduction in condensed heterochromatin (Boonsanay et al., 2016). Loss of Suv4-20h1 lead to a reduced population of quiescent cells and increased population of differentiated cells expressing MyoD (Boonsanay et al., 2016). In Suv4-20h1-abrogated quiescent MuSCs, there was a reduction in H4K20me2 and nucleosome density, and an increase in H3K4me3, at the distal regulatory region (DRR) that controls MyoD expression (Boonsanay et al., 2016). This indicates that DRR is more accessible in Suv4-20h1 knockout quiescent cells and this may allow expression of differentiationrelated gene, MyoD. Thus, Suv4-20h1 reduces expression from the MyoD locus, resulting in the maintenance and preservation of stem cells in a quiescent state (Boonsanay et al., 2016). Further research will be required to determine whether this is a role for H4K20 methylation exclusively in MuSCs, or whether similar changes occur in other stem cells as well.
The H4K20 monomethyltransferase Pr-set-7 was discovered to regulate NSC quiescence in Drosophila (Huang et al., 2021). Targeted DNA adenine methyltransferase identification (TaDa) was used to determine genomic loci where Pr-set7 binds. This analysis revealed Pr-set-7 binds to the promoter and transcriptional start sites of Wnt pathway coactivator earthbound1 (Ebd1) and cyclin-dependent kinase 1 (Cdk1) (Huang et al., 2021). The mRNA levels of Ebd1 and Cdk1 were depleted in pr-set7 mutant brains (Huang et al., 2021). Thus, by increasing expression of Cdk1 and Ebd1 in NSCs, Pr-set7 modulates the cell cycle and Wnt signaling, and thereby promotes NSC reactivation from quiescence (Huang et al., 2021). In these studies, The presence or absence of H4K20me1 in Ebd1 and Cdk1 promoters was not determined in this study (Huang et al., 2021).
In primary human dermal fibroblasts, mass spectrometrybased analysis of histone tails revealed that the histone lysine that showed the strongest change in methylation levels with quiescence was H4K20 (Evertts et al., 2013a). H4K20me2, and especially H4K20me3, were highly induced in quiescent cells when compared with proliferating fibroblasts as well as the fibroblasts that were specifically in the G1 phase of the cell cycle and contained the same DNA content as quiescent fibroblasts (Evertts et al., 2013a). In contrast, H4K20me1 levels were reduced in quiescent compared with proliferating fibroblasts. Further knockdown of Suv4-20h1 and Suv4-20h2 methyltransferases that catalyze formation of di and tri-methylated H4K20 with both small hairpin RNAs and small interfering RNAs (siRNAs) resulted in less chromatin compaction, consistent with a potential role for H4K20 trimethylation in chromatin conformation (Evertts et al., 2013a). Knockdown of Suv4-20h2 specifically, with siRNAs, resulted in increased proliferation and more cells in S phase (Evertts et al., 2013a). These findings support a role for Suv4-20h2 in both the regulation of the quiescence-proliferation transition and chromatin compaction with quiescence. Whether there is a direct link between these two activities, and whether they are mediated through H4K20 methylation will require additional studies.
Altogether, these studies highlight the critical role of H4K20 methylation, and in particular, H4K20me3, in regulating the transition between quiescence and proliferation. The consistent upregulation of H3K20me3 in multiple different quiescent models indicates it might play a larger role in establishing the functional state of chromatin in quiescence and may regulate specific gene expression changes necessary for quiescence entry, exit, and maintenance. Given the available data, it is possible that the H4K20 methylation marks act in cooperation with other marks, for instance it has been detected as part of bivalent promoters in combination with H3K4me3 and H3K36me3 (Xu and Kidder, 2018), and as part of heterochromatin in combination with H3K9me3 (Schotta et al., 2004). Taken together, the data support the possibility that H4K20me3 represents one part of a combinatorial code that regulates quiescence.

HISTONE ACETYLATION WITH QUIESCENCE General Properties of Histone Acetylation
An array of conserved lysine residues is present in the N-terminal tails of the four histones forming the nucleosome and these positively charged lysines interact with the DNA and the negatively charged patch formed by H2A/H2B residues of the nearby nucleosome. The tail of histone H4 has the strongest effect on chromatin compaction, followed by the tail of histone H3, which is followed by the tails of histones H2A and H2B (Ghoneim et al., 2021). Histone acetylation weakens the interactions of tails with the DNA and the negative patch, thus making the chromatin more accessible to RNA Pol II (Park and Kim, 2020). In addition to affecting the structural state of the chromatin, acetylation can also affect the proteins bound to chromatin. Lysine acetylation marks on histone tails can be recognized by bromodomaincontaining reader proteins (such as chromatin remodelers) that are generally associated with transcriptional activation (Agalioti et al., 2002;Barnes et al., 2019). The removal of acetyl groups by histone deacetylases can modulate the chromatin state to make the chromatin more compact. Further, in some cases, an individual lysine residue can be acetylated, methylated or ubiquitinated, and this establishes a potential for competition among these different lysine modifications.
The N-terminal tail of histone H3 is mainly acetylated at residues K9, K14, K18, and K23, while the corresponding region of H4 is acetylated at residues K5, K8, K12, and K16. Acetylation on H4K16 is particularly linked to unpacking of the chromatin and transcriptional activation-mainly due to diminished interactions between the H4 tail and the H2A/H2B acidic patch (Shogren-Knaak et al., 2006). The periodic nature of the amino acid spacing between the acetylatable lysines in the H3 and H4 tails suggests a potential for cooperative behavior among these amino acids (Strahl and Allis, 2000). In particular, this spacing has been noted to be reminiscent of the 3.6 residues per turn of an α-helix, raising the possibility that these acetyl marks may act as part of a combinatorial code that provides information about which proteins should bind to chromatin at a specific genomic region (Strahl and Allis, 2000).
The combined presence of histone acetyl marks and other histone modifications have been shown to result in defined outcomes (Winter et al., 2008;Zippo et al., 2009). As one example of this combinatorial effect, Bromodomain PHD Finger Transcription Factor (BPTF), an ATP-dependent chromatin remodeling protein, has increased affinity for histones with two different marks. The presence of H3K4me3 and H4K16ac together results in a 2-fold stronger BPTF binding affinity than either H3K4me3 or H4K16ac alone (Ruthenburg et al., 2011;Rando, 2012). This effect was not limited to H4K16ac as peptide microarray studies showed an increase in binding of PHD-Bromodomain constructs to histones when H3K4me3 was combined with any of multiple different H3 acetyl states (Fuchs et al., 2011;Rando, 2012). As another example, some proteins bind preferentially when a single epigenetic mark is found by itself, in the absence of another mark (Rando, 2012). This is exemplified by a PHD-Bromo domain chromatin regulator TRIpartite-Motif containing 24 (TRIM24) as this protein binds to H3K23ac to activate estrogen-responsive genes, but this binding is inhibited if H3K4 is also methylated (Tsai et al., 2010;Rando, 2012). Taken together, these reports suggest that a histone acetylation-dependent combinatorial histone code may encode information through the presence of individual histone marks, combinations of marks, their spacing, and the specific readers that recognize the marks. In the sections below, we focus on acetylation of histones H3 and H4 as these have been studied the most in relation to cellular quiescence.

H3 Acetylation With Quiescence
In budding yeast S. cerevisiae, a global decrease in acetylation levels was observed after quiescence entry that parallels the repression of transcriptional activity in quiescence (McKnight et al., 2015). Upon entry into quiescence, the budding yeast cells have lower levels of histone H3 acetylation compared with log phase cells. Immunoblotting for histone H3K9ac, H3K23ac, and pan-acetyl histone H3 revealed a significant reduction in the levels of these marks in quiescent compared with proliferating yeast (McKnight et al., 2015). In a complementary set of studies, using mass spectrometry, Mews et al. showed that when quiescent yeast cells reentered the cell cycle upon nutrient replenishment, there was a burst of histone acetylation, specifically at H3K9 and H3K14, immediately upon cell cycle entry (Mews et al., 2014). In contrast, de novo histone methylation occurred at a later cell cycle stage after the acetylation burst. As compared to histone methylation, histone acetylation was more closely correlated with the transcriptional changes that took place soon after the cell cycle entry. Western blot analysis revealed that levels of acetylated histones H4K5ac, H4K8ac, and H3K9ac drop during quiescence and then robustly increase over a 240-min time course after nutrient refeeding. Using ChIP-seq, Mews and colleagues found that when yeast enter stationary phase, H3K9ac levels decreased at genes identified by microarrays as growth genes and increased at genes expressed when yeast initiated quiescence following nutrient exhaustion (Mews et al., 2014). After the nutrients were replenished, acetylation of H3K9 increased at growth genes and decreased at stress response genes (Mews et al., 2014). Histone acetylation rapidly responded to changes in metabolic state, while histone methylation levels were largely constant (Mews et al., 2014). Thus, S. cerevisiae undergo a dramatic increase in acetylation, and not methylation, when they re-enter the cell cycle from quiescence as nutrients are re-introduced (Mews et al., 2014;McKnight et al., 2015;Young et al., 2017).
The histone lysine deacetylase Rpd3 has been identified as an important regulator of quiescence entry and maintenance in S. cerevisiae in a study by McKnight and colleagues (McKnight et al., 2015). When the yeast entered quiescence, motifs associated with transcription factors inactivated by nutrient exhaustion were associated with more repressive chromatin structure including more nucleosomes positioned over the motif and lower levels of histone acetylation (McKnight et al., 2015). Motifs for transcriptional activators that function in stress conditions had more open chromatin structure in the quiescent cells (McKnight et al., 2015). A large fraction of transcription factors with binding motifs that exhibit changes in chromatin structure with quiescence function by recruiting Rdp3 lysine deacetylase (McKnight et al., 2015). Deleting Rpd3 did not affect growth in the log phase and Rdp3-deleted cells maintained high viability in log phase (McKnight et al., 2015). Rpd3-deficient cells arrest after glucose exhaustion similarly to wild-type cells. However, fewer quiescent cells are formed in an rpd3 mutant and the quiescent cells that are formed show reduced long-term survival (McKnight et al., 2015). Analysis of chromatin structure at gene promoters revealed that Rpd3 has a significant role in establishing the transcriptional profile of quiescent cells by affecting the density of histone H3 and the acetylation of histone H4 at the promoters of genes regulated with quiescence (McKnight et al., 2015). The findings, taken together, support histone lysine acetylation as an important regulator of gene expression and viability during quiescence in S. cerevisiae (McKnight et al., 2015).
Studies in mouse NIH 3T3 embryonic fibroblasts that were serum starved to induce quiescence and then restimulated with serum addition revealed changes in histone acetylation (Knosp et al., 1991). Using electrophoretic and fluorographic techniques, in 1991, Knosp et al. reported that addition of serum resulted in an increase in the acetylation rate of all core histones within 15 min (Knosp et al., 1991). The sharp increase in acetylation rate was followed by a gradual decline in acetylation rate that continued until 8 h post serum stimulation (Knosp et al., 1991). At 10-12 h after serum stimulation, there was a strong increase in the acetylation rate of all core histones, detected based on increased incorporation of tritium into histone H3 from labeled acetate, a methodology that doesn't distinguish between acetylation of different lysine residues. This increase in acetylation was followed by an increase in histone synthesis (Knosp et al., 1991). The pattern differed among histones: H3 > H2A ∼ H2BB > H4 (Knosp et al., 1991). In response to 24 h of serum withdrawal, there was a decrease in DNA synthesis, histone synthesis and histone acetylation (Knosp et al., 1991).
In a more recent study, proliferating and quiescent primary human fibroblasts, were incubated in 13 C-labeled acetate and the rate of histone lysine acetylation was monitored using mass spectrometry (Evertts et al., 2013b). The quiescent fibroblasts accumulated labeled acetylated histones more slowly than proliferating fibroblasts (Evertts et al., 2013b) and differential labeling rates were observed for the acetylation of H3K9, H3K14 and histone H4 (Evertts et al., 2013b). For histone H4, mass spectrometry was used to monitor the extent of acetylation of a peptide that contained lysines K5, K8, K12 and K16 was monitored. Levels of unmodified peptide, peptide with one acetyl group and peptides with two acetyl groups were determined (Evertts et al., 2013b). Incorporation of 13 C in acetyl groups in quiescent cells was approximately half of that for proliferating cells. Even though the rate of acetylation was faster in proliferating fibroblasts, steady state levels of histone acetylation were similar in proliferating and quiescent fibroblasts (Evertts et al., 2013a,b).

H4 Acetylation With Quiescence
In yeast, H4 acetylation correlates with transcriptional activation and plays an important role in chromatin decompaction (discussed in more detail in the later section). Using immunoblotting, the overall levels of H4K5ac and H4K8ac were found to decrease sharply during quiescence in budding yeast S. cerevisiae. This was followed by a rapid increase of these histone marks during cell cycle reentry when nutrients were reintroduced (Mews et al., 2014). Also using immunoblotting, a reduction in histone H4K12ac was observed in quiescent S. cerevisiae (McKnight et al., 2015). Similarly, by mass spectrometry and western blot analysis, two separate studies found that during stationary phase in S. cerevisiae, histone acetylation dropped dramatically at H4K5, H4K8, H4K12, and H4K16 compared to their levels in exponential growth phase (Sandmeier et al., 2002;Ngubo et al., 2011).
In mammals, reduced histone acetylation has been associated with a transition to quiescence in pluripotent cells (Khoa et al., 2020). Embryonic stem cells with pluripotent potential are characterized by high levels of histone acetylation, high chromatin accessibility and an active pluripotency transcriptional network. Mouse embryonic stem cells (ESCs) can exist in a proliferative, naïve ground state achieved by maintenance with a MAPK/ERK Kinase (MEK) inhibitor and a Glycogen Synthase Kinase 3 (GSK3) inhibitor (Boroviak et al., 2014). Deletion of histone acetyltransferase Males absent On the First (MOF), which acetylates histone H4 lysine 16 (Finley et al., 2018;Atlasi et al., 2019), resulted in quiescence in mouse ESCs (Figures 1, 2, Table 4). Among the > 200 histone posttranslational modifications that were monitored in these MOFdeleted cells, only H4K16ac was significantly decreased (Khoa et al., 2020). Comparing RNA-seq gene expression data with global H4K16ac distribution obtained by ChIP-seq revealed that many genes associated with metabolic processes were regulated in a MOF-dependent manner, especially β-oxidation of fatty acids (Khoa et al., 2020). The degradation of fatty acids provided metabolites for oxidative phosphorylation and energy production in the proliferative state (Khoa et al., 2020). Inhibiting fatty acid oxidation was also sufficient to induce the proliferative ESCs into a quiescent state (Khoa et al., 2020). Taken together, these studies highlight an important role for histone H4 acetylation in regulating the proliferation-quiescence transition.
While addition of H4 acetyl groups is associated with the transition from quiescence to proliferation, the removal of H4 acetyl groups by Sitruins, a class of histone deacetylases (HDACs), has also been associated with the proliferation-quiescence transition (Ryall et al., 2015) (Figure 1). When mouse MuSCs exit quiescence to begin proliferating, there is a decrease in NAD + levels and activity of NAD + -dependent sirtuin 1 (SIRTS), an HDAC (Michan and Sinclair, 2007). This reduction in SIRT1 activity results in increased H4K16 acetylation, activation of muscle gene transcription, premature myogenic differentiation, reduced myofiber size, impaired muscle regeneration, and derepression of muscle developmental genes (Ryall et al., 2015).
These studies taken together support a role for acetylation in modulating transcriptional activation at important genes in different quiescence models. Loss of acetylation marks at genes critical to cell cycle and cell growth define the functional state of the chromatin during quiescence. In yeast and MuSCs, the alterations in acetylation have been functionally linked to quiescence (McKnight et al., 2015;Ryall et al., 2015).

RNA-MEDIATED ADDITION OF HISTONE MARKS
The mechanisms that dictate where histone marks are deposited are not well-understood. There is an emerging literature that suggests that in some instances, long non-coding RNAs (lncRNAs) can help to direct the deposition of histone epigenetic marks (Batista and Chang, 2013;Hanly et al., 2018). LncRNAs can act as "address-codes" by directing the writers of histone marks to specific chromosomal locations (Batista and Chang, 2013;Hanly et al., 2018). The deposited marks can then affect chromatin conformation, nucleosome positioning and transcription at nearby genes (Faghihi and Wahlestedt, 2009;Guttman et al., 2011). The lncRNA Xist, for example, interacts with polycomb repressive complexes PRC1 and PRC2, which are responsible for the addition of the H3K27me3 mark, and recruiting these complexes is one part of the process by which Xist induces transcriptional silencing of the X chromosomes in female mammals (Kohlmaier et al., 2004;Bousard et al., 2019). Another lncRNA termed HOTAIR represses the HoxD gene locus by recruiting PRC2 for H3K27me3 addition (Rinn et al., 2007). LncRNAs have been confirmed to play critical roles in the establishment and maintenance of heterochromatic regions (Deng et al., 2009).
LncRNAs including H19 are emerging as critical regulators of quiescence and proliferation in HSCs (Venkatraman et al., 2013;Yildirim et al., 2013;Luo et al., 2015). The H19-Igf2 locus is an imprinted region that affects organism growth (DeChiara et al., 1991). The differentially methylated region (DMR) upstream of H19 is the imprinting control region that enforces the expression of H19 from the maternal allele only while Igf2 is expressed from the paternal allele (Bartolomei, 2009). Deletion of the maternal, but not the paternal, H19 DMR in adult mouse HSCs resulted in a activation of quiescent HSCs and reduced function of the HSCs (Venkatraman et al., 2013). When maternal H19 DMR was deleted, the Igf2-Igf1r pathway was activated, the FoxO3 transcription factor translocated from the nucleus to the cytoplasm, and quiescent HSCs were activated to proliferate, resulting in their exhaustion (Venkatraman et al., 2013).
PAPAS lncRNAs have been studied in the context of cell quiescence ( Figure 3A). PAPAS, or "promoter and pre-rRNA antisense, " is a heterogeneous group of lncRNAs that are generated when ribosomal DNA (rDNA) is transcribed in the antisense direction and localize to the pre-rRNA coding region and the rDNA promoter (Bierhoff et al., 2010). PAPAS interacts with Suv-20h2, the methyltransferase that generates H4K20me3, and mediates deposition of this mark at rDNA during quiescence (Bierhoff et al., 2010). PAPAS and H4K20me3 were found to be upregulated in quiescent breast cancer cell line MCF7 generated with an estrogen receptor antagonist (Bierhoff et al., 2014). PAPAS and H4K20me3 are also induced upon terminal differentiation of human colon cancer cells (Caco-2), and upon quiescence in C2C12 mouse myoblasts, and mouse fibroblast-like 3T3-L1 cells (Bierhoff et al., 2014). Knockdown of endogenous PAPAS with siRNAs in MEFs decreased H4K20me3 levels (Bierhoff et al., 2014). Gain-and loss-of-function experiments, RNA-protein binding assays, and chromatin accessibility analysis revealed that PAPAS guides Suv4-20h2 to nucleolar chromatin, reinforcing quiescence-dependent transcriptional repression of ribosomal RNAs through H4K20me3-dependent chromatin compaction (Bierhoff et al., 2014).
Intracisternal A particle (IAP)-specific lncRNAs were found to be involved in H4K20me3-mediated chromatin compaction at IAP retrotransposons (Bierhoff et al., 2014). At IAP elements, both H4K20me3 and Suv4-20h2 levels increased in serumstarved MEFs, and the IAP chromatin was more compact (Bierhoff et al., 2014). Knockdown of IAP lncRNA with locked nucleic acid technology resulted in upregulation of the IAP major 7 kb transcript, supporting a role for the lncRNA in causing repression of the IAP gene, while transfection of IAP lncRNA sequences that interact with Suv4-20h2 resulted in increased Suv4-20h2 and H4K20me3 at IAPs, with no change in H3K9me3 levels (Bierhoff et al., 2010). These studies suggest that lncRNAs may provide another dimensionality to the histone code by potentially providing an address code that directs histone marks to specific genomic positions.

HISTONE MARKS AND CHROMATIN CONFORMATION DURING QUIESCENCE
Multiple studies have reported changes in chromatin compaction with quiescence with most, but not all studies (Kallingappa et al., 2016), reporting that when cells enter quiescence, their chromatin becomes more compact (Chiu and Baserga, 1975;Setterfield et al., 1983;Evertts et al., 2013a;Swygert et al., 2019). In S. cerevisiae budding yeast, quiescence induction was found to induce global changes in chromosomal organization (Rutledge et al., 2015) using a global chromosome conformation capture highthroughput sequencing technology called Hi-C. By monitoring the frequency of interloci interactions, the study revealed an increase in long range cis interactions and a decrease in short range interactions in quiescent yeast cells (Rutledge et al., 2015). The study further found that inter-centromeric interactions decrease during quiescence, while inter-telomeric interactions increase in quiescence, indicating that quiescence maintenance is associated with substantial topological reorganization (Rutledge et al., 2015). Another study has also shown that when yeast enter the stationary phase, the telomeres hypercluster, that is, they colocalize in a single location (Laporte et al., 2016). Mutant yeast strains that lack linker histone Hho1, or condensin, or contain mutations in histone H4 lysine 16 are unable to form these FIGURE 3 | Histone Modifications and Metabolism | (A) During quiescence, the writer Suv420h2 adds methyl groups to H4K20 via mediation with the lncRNA PAPAS. PAPAS, Suv4-20h2, and H4K20 methylation levels all increase with quiescence (Q) compared to proliferation (P), as indicated by the arrows (Bierhoff et al., 2014). (B) Acetyl-CoA acts as an acetyl group donor for histone acetyltransferases, including MOF, which is responsible for H4K16 acetylation during proliferation (Khoa et al., 2020). (C) In fission yeast, synthesis of S-adenosylmethionine (SAM) from methionine with the help of S-Adenosylmethionine synthetase (Sam1) is implicated in histone methylation. SAM acts as a methyl donor for methyltransferases; thus, Sam1 is necessary for proper cell growth and proliferation, as well as quiescence maintenance (Hayashi et al., 2018)  telomere hyperclusters (Laporte et al., 2016). A contraction of the nucleolus  has also been observed in quiescent yeast cells as well.
A study in human diploid fibroblasts investigated changes in chromatin compaction with quiescence (Criscione et al., 2016). Using Formaldehyde Assisted Isolation of Regulatory Elements (FAIRE) to investigate global DNAse I sensitivity and chromatin accessibility, quiescent cells were found to be more resistant to DNAse I treatment, indicating more compact chromatin (Criscione et al., 2016). This study used Hi-C to show genes switching between A and B compartments with quiescence (Criscione et al., 2016). The A-type compartment has a more open chromatin structure and is enriched for activating marks such as H3K36me3, H3K79me2, H3K27ac, and H3K4me1 (Lieberman-Aiden et al., 2009;Rao et al., 2014). The B-type compartment, on the other hand, is characterized by more densely packed chromatin and correlates with repressive marks such as H3K27me3, H3K9me3, and H4K20me3 (Lieberman-Aiden et al., 2009;Rao et al., 2014). Genes associated with cell proliferation were enriched in the group of genes switching from A to B compartments as cells entered quiescence (Criscione et al., 2016).
HP1β is a dimeric protein that binds to the H3K9me3 mark in constitutive heterochromatin and can bridge two H3K9me3-containing nucleosomes (Machida et al., 2018). Both HP1β and H3K9me3 localized to constitutive heterochromatin regions in proliferating B lymphocytes, while in quiescent B lymphocytes they did not (Baxter et al., 2004). These findings suggest that the overall structure of DNA may be altered in quiescent B lymphocytes in a way that affects the accessibility of histone writers and readers, especially in constitutive heterochromatin.
One possible mechanism through which the histone code could affect chromatin state during quiescence involves the histone methyltransferase Suv4-20h2 that generates H4K20me3. As mentioned before, H4K20me3 is involved in heterochromatin formation and induces chromatin compaction (Evertts et al., 2013a;Hahn et al., 2013). Suv4-20h2 activity and the H4K20me3 mark are upregulated with quiescence (Evertts et al., 2013a;Bierhoff et al., 2014). A fluorescence in situ hybridization (FISH) analysis performed in primary human dermal fibroblasts found that contact-inhibited quiescent fibroblasts had more compact chromatin than their proliferating counterparts (Evertts et al., 2013a). Furthermore, knockdown of Suv20h2 resulted in decreased compaction (Evertts et al., 2013a). As described above, crystal structures of nucleosomes reconstituted with histones containing the H4K20me3 mark had alterations in higher order structure (Lu et al., 2008). Fluorescence Recovery After Photobleaching (FRAP) analysis revealed that Suv4-20h2 binds tightly to heterochromatin (Hahn et al., 2013). Other studies have shown that Suv4-20h2 associates with pericentric heterochromatin through multiple, independent interaction sites on its C-terminal domain that directly bind to multiple heterochromatin protein 1 (HP1) molecules (Schotta et al., 2004(Schotta et al., , 2008Watanabe et al., 2018). The HP1 protein family consists of three members: HP1α, HP1β, and HP1γ. HP1α localizes to heterochromatin, HP1β is found in both heterochromatic and euchromatic regions and HP1γ is associated with actively transcribed genes (Vakoc et al., 2005;Lomberk et al., 2006). Dimeric HP1α binds two H3K9me3 marks in adjacent nucleosomes and forms a bridge between them, thereby compacting the chromatin (Watanabe et al., 2018). Consequently, the combinatorial effects of the H3K9me3 and H4K20me3 marks together represent an opportunity for a quiescence combinatorial histone code that affects chromatin conformation and induces chromosomal compaction.
Chromatin-bound Suv4-20h2 also recruits the cohesin complex (Hahn et al., 2013), which is composed of rings of Smc1-Smc3 dimers connected by Scc1/Rad21 (Nasmyth, 2011). The cohesin complex can contribute to establishment of the chromatin loop domains at specific DNA regions, thereby inducing chromatin compaction (Maya- Miles et al., 2019). In Suv4-20h1/Suv4-20h2 double-knockout cells, cohesin was absent from heterochromatic regions in G0 phase, which showed that Suv4-20h enzymes are required for loading or maintaining cohesin at heterochromatin (Hahn et al., 2013). Reintroducing Suv4-20h2 or only the non-enzymatic clamp domain of Suv4-20h2 rescued the loss of heterochromatin-associated cohesin suggesting the effects of Suv4-20h2 may be mediated independent of its effects on H4K20me3 (Hahn et al., 2013). This ability of Suv4-20h2 to recruit cohesin and compact chromatin may be critical for chromatin compaction with quiescence, as Suv4-20h2-deficient MEFs synchronized in G 0 , contained virtually no cohesin in regions of heterochromatin (Hahn et al., 2013). Additional experiments will be required to understand the functional consequences of a lack of cohesin in heterochromatin in Suv4-20h2-deficient quiescent cells.

CROSS-TALK BETWEEN HISTONE MARKS AND METABOLISM
An important emerging theme in epigenetic regulation is the close association between histone marks and metabolism. In yeast, metabolic signals such as the presence or absence of glucose can determine whether cells proliferate or arrest (Laporte et al., 2011), and therefore, histone marks that signal the relative abundance of glucose can potentially transmit information about nutrient abundance in a locus-specific manner to affect gene regulation. The metabolite acetyl-CoA serves as the substrate for histone acetyltransferase enzymes that generate acetylated histones (Kaelin and McKnight, 2013) (Figure 3B). Acetyl CoA is formed in mitochondria when pyruvate generated by glycolysis is committed to the TCA cycle (Martinez-Reyes and Chandel, 2020). Citrate formed in the mitochondria can be exported and converted in the cytoplasm to acetyl CoA (Martinez-Reyes and Chandel, 2020). In proliferating yeast and mammalian cells, higher levels of glycolysis and increased export of acetyl CoA from the mitochondria have been observed in proliferating compared with quiescent cells (Frauwirth and Thompson, 2004;Vander Heiden et al., 2009;Lemons et al., 2010). Thus, the metabolic profiles of proliferating cells may facilitate the generation of histone acetylation, and the formation of more open chromatin (Frauwirth and Thompson, 2004;Lemons et al., 2010). The role of histone acetylation in S. pombe was investigated with a strain that exhibited a temperature-sensitive mutation in the catalytic region of phosphopantothenoylcysteine synthetase (designated Ppc1) (Nakamura et al., 2012), an enzyme in the biosynthetic pathway for acetyl-CoA. S. pombe with mutations that inactivate the ability to synthesize acetyl-CoA fail to acetylate histones and are unable to re-enter the cell cycle after initiating G 0 following nitrogen withdrawal (Nakamura et al., 2012).
In fibroblasts, the oncogene c-MYC has been identified as a transcription factor that can affect global chromatin remodeling (Morrish et al., 2010). In rat fibroblasts, MYC activity increased glucose metabolism and acetyl-CoA production (Morrish et al., 2010). The presence of MYC caused a forty percent increase in 13 C-labeled acetyl-CoA on H4-K16 during cell cycle entry with serum stimulation during the G 0 to S transition (Morrish et al., 2010). Metabolic tracing revealed that in MYC-stimulated cells, MYC increases acetylCoA levels (Morrish et al., 2010). Further, the GCN5 histone acetylase enzyme, an enzyme that is important for histone acetylation in response to nutrients in yeast (McMahon et al., 2000;Cai et al., 2011), is a target of MYC (Morrish et al., 2010).
In addition to a connection between metabolism and histone acetylation, there is also a connection between metabolism and histone methylation. S-adenosylmethionine (SAM) is a metabolite that acts as a donor of methyl groups, which are transferred by methyltransferase enzymes to histones (Kaelin and McKnight, 2013) (Figure 3C). In fission yeast, the enzyme responsible for synthesizing SAM, S-adenosylmethionine synthase 1 (Sam1), is required for proper cell growth, proliferation, and quiescence entry and exit (Hayashi et al., 2018). Loss of Sam1 results in reduced levels of H3K4me2 and H3K4me3, and a significant decrease in cell growth and defects in G2 cell cycle arrest (Hayashi et al., 2018). Yeast with mutations in Sam1 cannot survive in quiescence initiated following nitrogen starvation, and once released from G 0 , they are unable to increase in cell size and restart DNA replication in proliferative conditions (Hayashi et al., 2018).
Taken together, these studies exploring the crosstalk between metabolism and histone marks highlight the close association between nutrient uptake, proliferation and histone posttranslational modifications. How this crosstalk is established and maintained in quiescence is an important area for further exploration.

GAPS IN THE FIELD AND FUTURE STUDIES
The decision whether to proliferate or exit the proliferative cell cycle requires cells to integrate multiple types of information and execute a complex series of molecular changes that impact many of the cell's activities. Transitioning between a proliferative and quiescent state is associated with many alterations including changes in the expression and activity of specific genes, changes in the conformation of chromatin, and changes in the subnuclear localization of chromosomes. The molecular mechanisms that cells use to make these decisions are actively being investigated and there are likely levels of regulation that have yet to be uncovered. Here we explored the literature to find evidence that histone post-translational modifications serve as a code that interprets information about the activity of signaling pathways and transmits that information to enable the commitment to a proliferative or quiescent state and the functional changes required for this transition.
So, is there a histone code for cellular quiescence? As highlighted in this review, there are multiple lines of evidence that suggest that specific histone modifications alone and in combination are added or removed as cells transition between proliferation and quiescent in multiple model systems. Histone H3K4me3, for instance, is found at the promoters of genes that are activated when cells become quiescent. H3K36 methylation is important for establishment and maintenance of the quiescent state in HSCs as inactivation of the SETD2 methyltransferase results in depletion of quiescent HSCs (Wagner and Carpenter, 2012). In quiescent mouse MuSCs, a combination of H3K9 and H3K27 marks was found to regulate cyclin A2 expression and proliferation (Cheedipudi et al., 2015). As another example, H4K20me3 levels increase in quiescent cells and knockdown of Suv4-20h2 results in a more rapid cell cycle (Evertts et al., 2013a). Finally, changes in the levels of histone acetylation have been observed with the transition between proliferation and quiescence (Ryall et al., 2015), and the close association between nutrient availability and acetyl CoA levels suggests histone acetylation as a possible link between nutrient availability and proliferation (Cai et al., 2011). For each of the marks described, there is evidence that the mark plays a functional role in some aspect of quiescence: quiescence entry, quiescence exit, quiescence maintenance, or quiescence depth. In some cases, functional studies have shown a causative role for a specific reader, writer or eraser. However, in many studies, and for many of the histone PTMs described, it is not clear whether the role of the histone marks is correlative or causative. Also, the contribution of histone marks seems to be context-and cell type-dependent as indicated by our survey of different quiescence models and conditions. For many of the marks, the direction of change was not consistent in different model systems. For instance, H3K27me3 was induced in some quiescence models and repressed in others (Baxter et al., 2004;Evertts et al., 2013a;Kallingappa et al., 2016;Maki et al., 2021). Other histone marks appear to be regulated with quiescence in some species, but are absent in others, such as H3K27me3 which is absent from yeast (Jamieson et al., 2013). Furthermore, quiescence studies have usually focused on one or a few histone marks at a time. Additional studies will be needed to determine whether a well-defined system-dependent or system-independent combinatorial histone code exists for cellular quiescence and how this histone code varies among species and tissues.
The findings in this review show that certain histone marks, and even combinations of histone marks, are altered with quiescence and some of these changes are consistent among model systems. Future studies will be needed that systematically determine the genomewide deposition of histone marks and combinations of histone marks in proliferating and quiescent cells in multiple model systems. These investigations would cover both traditional and non-traditional histone marks in multiple quiescent states. Beyond just comparing proliferating and quiescent cells, these inquiries would look into histone marks in cells that initiate quiescence in response to different signals in the same cell type and cells that have been quiescent for different durations of time to achieve different depths of quiescence (Rodgers et al., 2014;Kwon et al., 2017). While many of the studies we review have not discovered a clear relationship between changes in individual histone marks and altered gene expression (Liu et al., 2013;Lee et al., 2016), a careful analysis of combinations of marks might determine if the changes in the combinatorial pattern of histone marks is a better indicator of gene expression changes with quiescence.
Many questions remain when studying histone marks in quiescence. What triggers the deposition of these marks? In particular, how do different histone writers and readers coordinate to establish marks at the appropriate time? In some cases, establishment and maintenance of histone marks associated with an open chromatin structure and active gene transcription may be achieved through positive feedback loops . These positive feedback loops can be formed when proteins that write histone marks also read the same mark . For instance, SETD1 complexes not only catalyze the formation of H3K4me3 marks, but may also recognize the same H3K4me3 mark, bind to it, and continue to generate additional H3K4me3 modifications (Shi et al., 2007;Murton et al., 2010;Zhang et al., 2015). Dot1, the H3K79 methyltransferase in yeast (Guan et al., 2013), recognizes modifications on the histone H2B tail and, in human, also binds phosphorylated forms of RNA polymerase II at the transcription start sites of actively transcribed genes (Kim S.K. et al., 2012). H3K4me3 can recruit histone acetyltransferases that add acetyl groups as well as deacetylases that remove acetyl groups , resulting in dynamic turnover of histone acetylation marks when H34me3, but not other marks such as H3K79me3 or H3K36me3, are present (Crump et al., 2011). For repressive chromatin, PRC2 not only generates H3K27me3 but also binds to it, resulting in a positive feedback loop in which local chromatin structure allows H3K27me3 to be deposited over chromatin regions to form domains Zhang et al., 2015). Crosstalk between H3K27me3 and monoubiquitinated H2A on lysine 119, H2AK119ul, has also been proposed as enzyme complexes that deposit each mark may recognize the other mark Cooper et al., 2014;Kalb et al., 2014), allowing for the reinforcement of heterochromatic regions . Negative feedback between histone marks also occurs as activating marks can inhibit the activity of enzymes that place repressive marks, and vice versa . For instance, activating marks H3K4me2/3 and H3K36me2/3 can inhibit the activity of PRC2 and prevent the deposition of H3K27me3 repressive marks (Schmitges et al., 2011). Experiments in which the impact of modulating a specific reader and writer on the levels of multiple marks in the context of quiescence models may shed light on this question. In particular, such experiments may shed light on how positive and negative feedback loops are interrupted at proliferation-associated genes and reestablished at quiescenceassociated genes.
The development and application of new technologies will facilitate future studies investigating how combinations of histone marks coordinate. One valuable approach will be the ability to visualize proliferation and quiescence decisions in organisms in real time. Toward this end, a recent paper describes the adaptation of a biosensor for CDK activity (Spencer et al., 2013) to monitor cell division in two model organisms, Caenorhabditis elegans and zebrafish (Adikes et al., 2020). CDK activity was higher at the end of a cell division in cases in which the cell went on to divide (Adikes et al., 2020). Such biosensors could be used in conjunction with visualization of histone marks to assess whether histone marks individually or in combination can predict whether a cell will proliferate.
We anticipate that CRISPR-Cas9 will prove to be a powerful methodology for understanding the impact of histone modifications. One important challenge in understanding the impact of different histone marks and combinations of histone marks has been developing specific systems to test their functional importance. Many studies to date have focused on investigating the role of specific readers and writers with knockdown and knockout approaches. Using CRISPR-Cas9, further studies will likely allow the inactivation of specific histone modifiers using protein degradation systems that allow the proteins to be degraded in proliferating or quiescent cells with defined timing thus permitting a more detailed dissection of their role (Wu et al., 2020). Functional dissection of histone readers and writers is sometimes complicated as they have non-histone targets as well (Cornett et al., 2019). In fact, nearly 3,000 human non-histone proteins have been reported to have a lysine that can be methylated (Hornbeck et al., 2015;Cornett et al., 2019). An alternative approach is to test for the functional consequences of modifying the histones themselves. To achieve this, CRISPR-Cas9 has been used in S. cerevisiae to generate yeast strains with different combinations of mutations at histone tail lysines for histone H3 and H4, allowing the investigators to assess the effects of loss of different combinations of histone marks (Fu et al., 2021). In Trypanosoma brucei, precise editing of genes in multicopy arrays was performed with CRISPR-Cas9, allowing for the replacement of histone H4K4 with H4R4 to mimic the constitutively non-acetylated state. The authors achieved 90% replacement of the 43 histone H4 copies to H4R4 (Vasquez et al., 2018).
CRISPR/Cas9 will also be valuable for its capacity to specifically target chromatin writers and readers to specific genomic regions. As an early example of this technology, "programmable chromatin kinase" dCas9-dMSK1 was generated by fusing nuclease-deficient CRISPR/Cas9 (or dCas9) to a histone H3 kinase . When this protein was targeted to specific promoters with guide RNAs, there was an increase in histone H3 serine 28 phosphorylation at the target genes' promoters and an increase in expression of the targeted gene. Such studies are likely to be valuable for defining causal connections between histone PTMs, their activity at specific genomic loci, and outcomes such as gene expression and proliferation.
We anticipate that the availability of new technologies that allow us to better visualize the relationships among histone marks and chromosomal organization on a cellby-cell basis will also benefit studies of the histone code in quiescence. A recent report built upon sequential FISH (seqFISH) and multiplexed FISH methods to target 3,660 loci in individual mouse ES cells (Takei et al., 2021). These studies revealed that nuclear zones were created by combinatorial chromatin patterns (Takei et al., 2021). Repressive histone marks H4K20me3, H3K9me3, and histone variant H2A1 were found together and colocalized with DAPI-rich regions (Takei et al., 2021). A second heterochromatic pattern of H4K20me2, H3K27me3 and H3K27me2 was also observed (Takei et al., 2021). Active histone marks, H3K9ac and H3K27ac, and RNA Polymerase II serine 5 phosphorylation localized to nuclear speckles and were excluded from heterochromatin and nuclear lamina (Takei et al., 2021). Similar studies comparing proliferating and quiescent cells could shed light on the role for a potential histone code in establishing a cell's proliferative fate. In particular, single cell analyses of proliferating and quiescent cells could reveal patterns that are consistent among cells versus those that are more variable from cell to cell.
We anticipate that with time the role of additional histone PTMs and histone variants and their roles in regulating chromatin structure and gene expression will become clearer. As one example, non-tail globular histone marks, such as H3K36 and H3K122 acetylation marks found in active gene promoters and a subset of enhancers, can contribute to the histone code and expand the possibilities for combinatorial histone PTMs that provide position-specific information to readers (Pradeepa et al., 2016). Non-enzymatic histone modifications, such as glycation, acylation and lipidation, generated by spontaneous chemical reactions also have the potential to alter chromatin structure and regulate genetic processes (Maksimovic and David, 2021). The non-enzymatic modifications, like the histone marks described in this review, can be "erased" by scavenger systems (Maksimovic and David, 2021). Improved mass spectrometry and chemoproteomics will likely provide important new insights into the role of these modifications in cellular decisions including the commitment to proliferation (Maksimovic and David, 2021). As another example, additional studies are likely to reveal that changes in the specific histone variants present at different positions in the chromatin cooperate with histone PTMs to alter chromatin state and reader proteins. Some histone variants such as histone H3.1 and histone H3.3 contain different residues, such as amino acid S31 in histone H3.3 and A31 in histone H3.1 that can alter the properties of the chromatin and its accessibility to the transcription apparatus (Armache et al., 2020). Linker histone H1 variants can also affect folding of nucleosome arrays and nucleosome compaction (Fyodorov et al., 2018).
Combinations of multiple distinct genome-wide, highthroughput analyses performed in models of proliferating and quiescent cells, will also be needed to allow us to dissect the role of histone modifications and other contributors to gene expression and functional changes with quiescence. Studies in which changes in the genome-wide localization of histone marks can be correlated with chromatin accessibility using ATAC-seq, DNA methylation, and Hi-C to assess A and B compartments and topologically associating domain boundaries have the potential to yield greater insight. Combining these datasets and analyzing them with deep learning algorithms, may allow scientists to predict which combinations of marks are associated with changes in chromatin accessibility and gene expression.

CONCLUSION
Based on the findings thus far, the specific histone marks discussed in this review, methylation of H3K4, H3K9, H3K27, H3K36, H3K79, H4K20, and acetylation of H3 and H4 have been discovered to be regulated with quiescence in different model systems. Readers, writers and erasers of these marks have been found to functionally contribute to the quiescence-proliferation transition and quiescence maintenance. The model systems employed for these studies include nutrient depletion and spore formation in yeast, cell culture models in which different antiproliferative signals are employed, and models in which quiescent stem cells are visualized in situ or isolated and characterized. One of the key limitations for the field is the fact that knockout or knockdown of readers, writers and erasers can impact not just the histone marks under study, but also the PTMs of other cellular proteins as well. Also, most studies have investigated the effect of a single modification in isolation rather than the impact of several modifications together and the relationships between marks. More research with emerging technologies will be needed to determine whether there is a quiescence histone code and if so, how changes in histones are used to create a complex and context-dependent grammar that incorporates not just levels of histone marks and their readers, but also the chromatin context (in 2D as well as in 3D). Future studies will likely address the causality of histone marks and phenotypic changes that serve as regulators of quiescence and landmarks of quiescence. These studies may define the functional consequences of these marks in terms of gene expression, chromatin conformation, and chromosomal positioning. These studies will also likely reveal the role of traditional and non-traditional histone PTMs and other changes to chromatin in the decision whether to proliferate, the mechanisms that cause cell cycle arrest, the maintenance of cells during quiescence, and the determination of quiescence depth. Altogether, these results indicate that the regulation of histone marks may help to maintain the delicate balance between quiescence, proliferation, differentiation, and cell death, with different model systems and cell types likely using both overlapping and distinct aspects of the information contained in these histone marks.

AUTHOR CONTRIBUTIONS
MM conceived the topic and design of the review. KB and KS prepared all the figures and tables. KB, KS, MM, and HC did the literature survey. MM, HC, KB, KS, and KA wrote the manuscript. MM and HC supervised the project. All authors were involved in proof-reading.