Introduction {#S1}
============

The cell nucleus is a complex organelle, in which genetic material (DNA) is compacted up to several thousand fold and organized into chromosomes in a form that allows the genome to be replicated, repaired and transcribed in a cell-type specific manner ([@R66]; [@R67]). In three-dimensional space, the interphase nucleus represents a highly compartmentalized organelle that consists of the nuclear membrane, individual chromosomes occupying the distinct territories, as well as a number of nuclear bodies (nucleoli, Cajal bodies, promyelocytic leukaemia (PML) bodies, nuclear speckles, Polycomb bodies, etc.) located in the inter-chromosomal compartments ([@R70]; [@R83]) (**[Fig. 1 A](#F1){ref-type="fig"}**). In contrast to the cytoplasmic organelles, distinct nuclear compartments are not separated by the membranes and are formed by the dynamically-changing protein-protein and/or nucleic acid-protein complexes ([@R39]).

In the mammalian epidermis, progenitor cells of the basal layer proliferate and differentiate, forming multi-layered epithelium, as well as give rise to a variety of skin appendages (hair follicles, nails, glands, etc.) ([@R14]; [@R58]; [@R79]). Keratinocyte differentiation in the epidermis and skin appendages is driven by coordinated changes in the expression of the lineage-specific genes including those that constitute several genomic loci (Keratin type I and II loci, Epidermal Differentiation Complex (EDC), keratin-associated protein locus, etc.), and encode the essential structural components of the epidermal barrier or appendage-specific cellular products ([@R9]; [@R127]). Gene expression in these loci is regulated in a cell type-specific and differentiation stage-specific manner by a number of signaling pathways and transcription factors ([@R14]). However, mechanisms that regulate the balance of expression between different genes in these loci and how cell type-specific and differentiation stage-specific patterns of gene expression are achieved at the levels of whole loci remain unclear.

During transition of keratinocytes from basal epidermal layer to granular layer, the nucleus undergoes programmed transformation from highly active status, associated with execution of the genetic program of epidermal barrier formation, to fully inactive condition and, through a process of DNA degradation, becomes a part of the keratinized cells of the cornified epidermal layer ([@R53]). During last two decades, a tremendous progress has been achieved in understanding the genetic basis of skin disorders and identification of a large number of mutations underlying the development of many pathological skin conditions associated with alterations of the epidermal differentiation and integrity, hair loss and pigmentation \[for review, see ([@R131]; [@R137])\]. However, despite the fact that the genome-wide association studies and comparative genome analyses brought invaluable insights into the physiological role of genetic information, the question how genome of keratinocytes is organized beyond its linear sequence and which epigenetic mechanisms control gene expression programs in healthy and diseased epidermis remain largely unknown.

Currently, epigenetics is defined as a discipline that, in broad terms, investigates the heritable changes in the genome functions occurring without alterations to the DNA sequence ([@R10]). Epigenetic inheritance is often based on the inheritance of relatively stable chromatin states in the distinct genomic regions ([@R10]). Despite the fact that not all chromatin structural states are epigenetically transmitted during cell cycle ([@R117]), such mechanisms play an important role in the control of cellular functions in living organisms and are currently considered as a driving force of the development, phenotypic plasticity and evolutionary adaptation ([@R47]).

During last few years, it was shown that epigenetic mechanisms are involved in the control of epidermal development and terminal keratinocyte differentiation. In particular, an involvement of the DNA methytransferase DNMT1, histone demethylase JMJD3 and histone methyltransferase Setd8, histone deacetylases 1/2, ATP-dependent chromatin modifying enzymes (Brg1, Mi-2β), the polycomb group proteins (Bmi1, Cbx4, Ezh2, Jarid2) and genome organizer Satb1 in establishing tissue-specific differentiation programs in the epidermis was demonstrated ([@R38]; [@R40]; [@R44]; [@R45]; [@R52]; [@R71]; [@R73]; [@R84]; [@R91]; [@R93]; [@R129]). However, many aspects of the epigenetic control of gene expression programs in keratinocytes including an integration of the distinct levels of regulation in the context of three-dimensional nuclear organization and its dynamics associated with terminal keratinocyte differentiation remain to be elucidated.

Here, we overview data on epigenetic mechanisms that control tissue-specific gene expression programs in keratinocytes. We focused our attention on normal epidermal keratinocytes, whereas epigenetic mechanisms underlying the response of keratinocytes to ionizing radiation, as well as changes of keratinocyte epigenetic status in pathological skin conditions are outside of the scope of this review. Also, the roles of non-coding and microRNAs in the control of keratinocyte differentiation are not discussed here, because of several recent reviews published elsewhere ([@R15]; [@R142]; [@R150]). However, we propose the directions on how to further explore the role of epigenetic mechanisms and their cross-talk with other regulatory systems in the control of keratinocyte growth and differentiation in normal and diseased skin.

Distinct levels of chromatin organization and the control of gene expression programs in keratinocytes {#S2}
======================================================================================================

Eukaryotic genome is organized into distinct chromosomes that are formed by a nucleo-protein complex called chromatin. Chromatin consists of the nucleosomes, in which 147 base-pairs of genomic DNA form a loop around the histone octamer core represented by two molecules of each core histone (H2A, H2B, H3 and H4) ([@R3]). Nucleosomes are connected by non-histone associated linker DNA of approximately 10-80 bp length resulting in a formation of the "beads on the string" structures, which are folded into 30 nm fiber ([@R49]). Histone H1 binds to DNA as it exits the nucleosome, interacts with the linker DNA and plays an essential role for the condensation of nucleosomes and 30 nm fibers into higher order loop structures. These loop structures are less rigid compared to the nucleosomes and are folded via still poorly understood mechanisms into larger loops to accommodate the genome within the 3D nuclear space ([@R64]) (**[Fig. 2 A](#F2){ref-type="fig"}**).

Epigenetic mechanisms of the control of gene expression include several levels of regulation which are based on: Distinct patterns of DNA methylation and recently discovered DNA hydroxymethylation \[reviewed in ([@R50])\].Variability in distribution of the core histone variants and histone covalent modifications in the nucleosomes at the specific genome regions \[reviewed in ([@R24])\];Changes in the nucleosome positioning relatively to specific DNA sequences (translational nucleosome positioning) and ATP-dependent modulation of the DNA-histone interactions within nucleosomes (ATP dependent chromatin remodeling) \[reviewed in ([@R29])\];Changes in higher-order chromatin folding and genome organization within 3D nuclear space (higher-order chromatin remodeling) \[reviewed in ([@R33]; [@R43])\].

Execution of cell-type specific differentiation program in multi-potent progenitor cells requires a high degree of coordination between the distinct levels of chromatin organization ([@R72]; [@R83]; [@R125]). Below, we review the mechanisms that control chromatin remodeling at distinct levels of its organization with special emphasis on execution of keratinocyte-specific gene expression programs in developing and postnatal mammalian skin (data on the genetically-modified in vivo models illustrating the roles of chromatin remodeling factors in the control of skin development are shown in **Table 1**).

DNA methylation and hydroxymethylation {#S3}
======================================

In plants and vertebrates, DNA methylation is a post-replication modification resulting in addition of the methyl group at C5 position of cytosine to form 5-methylthytocyne (5mC). In vertebrates, this modification occurs predominantly at CpG dinucleotides ([@R50]). Approximately 60-70% of CpG sites are methylated in mammalian genome, except the relatively short regions with high CpG density called CpG islands mainly located at the promoters and first exons of the majority of housekeeping genes ([@R119]) (**[Fig. 2 B](#F2){ref-type="fig"}**). CpG sites are often differentially methylated in a cell type-specific fashion in gene regulatory regions beyond proximal promoters, particularly in recently discovered CpG island shores, located at approximately 2 kb from the transcription start sites or even further ([@R50]).

DNA methylation is usually associated with gene repression and is catalyzed by DNA methyltransferases (DNMTs) ([@R50]). In mammals, DNMT1 is the major methyltransferase that maintain the methylation status of the genome via methylation of cytosine within hemi-methylated CpG sites after DNA replication ([@R50]). DNMT3a and DNMT3b are the major vertebrate de novo DNA methylases, while DNMT3l, which lacks the catalytic domain, binds to DNMT3A/b and facilitates their chromatin targeting ([@R62]). DNA methylation patterns are significantly changed during distinct developmental stages (in germ cells and pre-implantation embryos), as well as during adult somatic cell differentiation, thus representing important epigenetic mechanisms in establishing cell-type specific programs driving tissue development and homeostasis ([@R50]; [@R119]).

In general, DNA methylation represses gene expression either through inhibiting transcription factor binding to DNA, or via interaction with methyl-DNA binding proteins that target repressive chromatin remodeling complexes to the methylated genome regions ([@R50]; [@R77]; [@R119]). However, DNA methylation is also required for induction of the subset of C/EBPα target genes by creating new binding sites from half-CRE sequences for this transcription factor ([@R120]), demonstrating that in some cases, DNA methylation contributes to gene activation.

An essential role of the DNMT1 in maintenance of the epidermal progenitor cells and the epidermal tissue renewal was recently demonstrated in human 3D culture and after grafting of DNMT1-deficient keratinocytes onto *nude* mice ([@R75]; [@R129]). The lack of DNTM1 in epidermal progenitor cells results in the severe defects in cell proliferation and loss of the tissue self-renewal capacity ([@R129]). Importantly, analysis of changes in gene expression program after DNMT1 ablation in epidermal progenitor cells showed abnormal induction of the genes associated with cell cycle arrest (cycling-dependent kinase inhibitors p15^INK4B^ and p16^INK4A^), while co-expression of Cdk4 and Cdk6 partially rescued the effect of DNMT1 inactivation on cell proliferation ([@R129]). Expressions of DNMT1, DNMT3A/B in keratinocytes are strongly elevated after UV irradiation, which contributes to gene silencing during UV-induced photocarcinogenesis ([@R109]). These results clearly demonstrate the importance of maintenance of the DNA methylation patterns in preserving epidermal progenitor cell identity and the remodeling of this pattern during terminal differentiation and UV response.

However, some gene promoters remain active being methylated, and in undifferentiated primary mouse keratinocytes the promoters of about 20% actively transcribed genes are methylated ([@R120]). Methylation of these promoters is important for binding of the transcription factor C/EBPα and subsequent activation of gene expression, while inhibition of DNA methylation by 5-aza-cytidine results in failure to induce the subset of these genes during calcium-induced differentiation in culture ([@R120]). 5-aza-cytidine treatment also differentially affects gene expression within the EDC locus in normal human keratinocytes: while expressions of SPRR1/2 and involucrin were markedly increased, the expression of S100A2 was decreased compared to controls ([@R42]). These data suggests that DNA methylation is capable of regulating both gene repression and activation and that the basal status of promoter methylation is important for the control of expression of individual genes in keratinocytes.

Interestingly, a novel DNA modification - the hydroxylation of the 5mC in DNA into 5-hydroxymethylcytocine (5hymC) regulated by TET enzymes, has recently been demonstrated both in vitro and in vivo ([@R144]). DNA hydroxymethylation leads to conversion of 5mC into non-modified C and result in changes of the DNA methylation patterns ([@R50]; [@R144]). It has been proposed that hydroxylation of the 5mC occurs at the regulatory regions of those genes that are programmed to be activated during development and this mechanism plays a role in orchestrating the balance between pluripotent and lineage committed states in embryonic stem cells ([@R146]). However, the significance of this novel DNA modification in the control of lineage-specific gene expression programs in keratinocyte progenitor cells remains unclear.

Covalent histone modifications {#S4}
==============================

Histone molecules are subject to numerous post-translational covalent modifications (PTM), particular in their N-terminal tails. Such modifications include acetylation, methylation, phosphorylation, ADP-ribosylation, ubiquitination, sumolyation, argenine deamination and proline isomerysation ([@R1]). Histone PTMs have profound effect on the chromatin structure and are actively involved in regulation of gene transcription via at least three major mechanisms that modulate: a) DNA-histone contacts within the core nucleosome particles; b) Inter-nucleosomal contacts and higher-order chromatin folding; c) Changes in interaction of the chromatin with non-histone regulatory proteins ([@R1]; [@R24]; [@R136]).

Modification of histones is catalyzed by specific groups of enzymes, including lysine acetyl-transferases, lysine and argenine methyl-transferases, kinases etc. Some modificatons (for example lysine acetylation, methylation and ubiquitination) are reversible. The removal of the acetyl and/or methyl groups from the histone molecules is also catalyzed by specific enzymes (including lysine or histone deacetylases, lysine or histone demethylases, etc.) ([@R1]). These enzymes do not usually have absolute substrate specificity (for example, CBP/P300 acetylates several lysine residues in all core histones) and could also have non-histone substrates ([@R1]). We review here the data on two most studied histone modifications (methylation and acetylation), whereas the impact of other modifications in the control of gene expression remain to be defined.

***Histone methylation*** is dynamically regulated by histone methyltransferases and demethylases ([@R143]). Analysis of the genome-wide patterns of histone PTMs, RNA polymerase II occupation and gene expression revealed the complex correlation between histone methylation and gene transcription status ([@R8]). These studies showed a strong correlation between the histone methylation patterns and the transcription activity in the distinct genomic regions (**[Fig. 2 B](#F2){ref-type="fig"}**). Lysine methylation in the histone 3 molecule could be associated with either transcription activation or repression depending on the extent (mono-, di-, or tri-methylation) and position of this modification in the histone tails ([@R8]).

Actively transcribed genes are characterized by trimethylation of the histone 3 at lysine 4 position (H3K4me3) that marks gene transcription start sites (**[Fig. 2 B](#F2){ref-type="fig"}**). H3K4me3 also enriched in the genes with the paused RNA polymerase II and some inactive genes without the polymerase II bound to the promoter, while H3K4Me2/1 occupies broader area around the transcription start site and 5' region of active genes ([@R143]). In addition, H3K79me2/3 is enriched through the 5' portion of actively transcribed genes, and H3K36me3 is associated with 3' parts of such genes ([@R143]).

In contrast, chromatin of the inactive genes is enriched in H3K9me2/3, H3K27me2/3 and H4K20me3 modifications ([@R143]) (**[Fig. 2 B](#F2){ref-type="fig"}**). The role of H3K27 trimethylation in human epidermal progenitor cells has been revealed by inactivation of H3K27me3 demethylase JMJD3 leading to block of the progenitor cell differentiation ([@R130]). The genes normally induced during epidermal differentiation (such as *KRT1, KRT10, IVL* and *S100A8)* were repressed in epidermal progenitor cells lacking JMJD3 protein and their induction during differentiation were associated with H3K27me3 de-methylation in the promoter regions ([@R130]).

Histone H4 mono-methylation at lysine 20 is regulated by Setd8 histone methytransferase, which ablation in the epidermis results in impairment of both cell proliferation and differentiation ([@R38]). Interestingly, Setd8 is a transcriptional target of c-Myc and an essential mediator of Myc-induced epidermal differentiation. However, Setd8 also serve as upstream regulator of p63 expression, which is lost in Setd8-null epidermis ([@R38]). Setd8 inhibits apoptosis in the skin and its activity is essential for epidermal stem cell survival, proliferation and differentiation, at least, in part, via regulation of p63 expression ([@R38]).

The analysis of global histone modifications reveals that in addition to gene promoters enriched in either active (H3K4me3) or repressive (H3K27me3) marks, some gene promoters are simultaneously enriched by both types of histone modifications and called "bivalent domains" ([@R34]) (**[Fig. 2 B](#F2){ref-type="fig"}**). Bivalent promoters are frequently seen in embryonic stem cells and depict the specific "paused" state for a number of genes. Since H3K4me3 is often associated with transcription initiation, while H3K27me3 is associated with gene repression ([@R143]), differentiation of stem cells towards distinct cell lineages is accompanied by either activation of selected "paused" genes and loss of H3K27me3 mark, or by complete repression of other genes and loss of the H3K4me3 active mark ([@R4]).

Although bivalent chromatin domains were identified in hematopoetic and neuronal stem cells and some T cells ([@R54]), recent studies identified only a small number of the bivalent chromatin domains in the hair follicle quiescent and activated stem cells compared to embryonic stem cells ([@R91]). These data suggest that bivalency is rather required for maintenance of broader differentiation potentials in pluripotent stem cells and is no longer required in adult tissue-specific stem cells with more restricted differentiation options ([@R91]).

In adult epithelial stem cells in the skin, actively transcribed genes that maintain stemness show both H3K4me3 and H3K79me2 histone marks, while repressed non-epidermal genes and genes activated during cell differentiation are enriched by H3K27me3 ([@R91]). Stem cell transition towards hair follicle-specific differentiation is accompanied by loss of H3K27me3 and appearance of H3K4me3 and H3K79me2 in the chromatin of those genes (Lef1, Bmp4, Wnt5a, Msx1, etc.) that become active, whereas the chromatin of key stemness genes (Cd34, Sox9, Nfatc1) that become repressed in transient-amplifying cells showed loss of H3K4me3/H3K79me2 and appearance of H3K27me3 ([@R91]). However, chromatin of the cell cycle-associated genes that are repressed in stem cells do not contain H3K27me3 and show instead H3K4me3, while their activation in transient amplifying cells is accompanied, in addition to H3K4me3, by the appearance of H3K79me2 ([@R91]). These data suggest that more complex regulatory mechanisms are involved in the control of repression of cell cycle-associated genes in stem cells, which may include the activity of anti-proliferative signaling pathways, lack of pro-proliferative signaling molecules, etc.

Histone acetylation {#S5}
-------------------

Acetylation state of lysine residues of histones affects gene expression through influence on chromatin conformation and is controlled by two counteracting enzyme families, histone acetyl transferases (HATs) and histone deacetylases (HDACs) ([@R134]). Histone hyper-acetylation has been correlated with transcription activation of numerous genes ([@R4]), and several transcription co-activators including CBP/P300 and PCAF/CGN5 are, in fact, histone acetyl-transferases ([@R134]). On the other hand, histone hypo-acetylation is usually correlated with gene repression and some transcription co-repressors possess histone deacetylase activity ([@R1]) (**[Fig. 2 B](#F2){ref-type="fig"}**).

The promoter regions of actively transcribed genes are enriched in histone H3 acetylated at lysine 9 (H3K9ac), H3K18ac and H2B12ac, whereas H4K12ac and H4K16ac are broadly distributed through the whole body of the active genes ([@R143]). Quiescent stem cells in the interfollicular epidermis and the hair follicle bulge are hypo-acetylated at histone H4, whereas c-Myc-induced exit from stem cell niche correlated with increased histone H4 acetylation ([@R56]). The significance of global histone acetylation levels for keratinocyte-specific gene expression programs was also demonstrated *in vitro* using histone deacetylase inhibitors trichostatin A or sodium butyrate, which both induced growth arrest in primary human keratinocytes via inhibition of expression of Cdk1 and stimulation of expression of differentiation markers involucrin and transglutaminase 1 ([@R42]; [@R124]). Consistently with these data, trichostatin A treatment of cultured epidermal explants resulted in inhibition of cell proliferation and stimulation of terminal differentiation ([@R100]). However, the effects of trichostatin A on gene expression in the EDC locus were gene-specific and also resulted in the inhibition of the expression of profillagrin and loricrin genes, suggesting that more complex mechanisms are involved in the control of expression of each individual gene in the EDC locus ([@R100]; [@R121]).

Histone deacetylases (HDACs) are part of a multiprotein family, in which each member has its specialized functions ([@R21]). HDACs play pivotal roles in coordinating the interaction of signaling pathways with chromatin remodeling and transcription factor function to modulate gene expression ([@R21]). Repressive action of HDAC1/3 complex is also achieved via a specific histone H2A monoubiquitination at lysine 119 by 2A-HUB ubiquitin ligase, which blocks RNA PolII release at the early stage of elongation ([@R151]). Conditional deletion of HDAC1/2 in the mouse epidermis demonstrated their crucial role in the control of the gene expression program in the epidermal progenitor cells and epidermal differentiation ([@R86]). The HDAC1/2 -/- epidermis failed to differentiate and remained single-layered, similarly to p63-null mice ([@R86]). The analysis of gene expression program in HDAC1/2 deficient mouse epidermis indicated that expression of p63 gene and its positive gene targets were not affected, whereas selected genes normally repressed by dNp63 including p21 and p16INK4A were up-regulated in the absence of the HDAC1/2 activity ([@R86]). Correspondingly, upregulation of HDAC1/2/3 in ARNT-deficient human epidermal keratinocytes resulted in suppression of AREG/EGFR pathway and concomitant enhancement of terminal differentiation ([@R121]).

Polycomb group (PcG) proteins {#S6}
=============================

PcG proteins were originally identified in Drosophila as negative regulators of Hox gene expression essential for the fly body patterning (reviewed in ([@R132]). In mammals, PcG proteins form two families of the polycomb repressive complexes (PRCs) called PRC1 and PRC2, which compact chromatin and promote transcriptional silencing ([@R40]; [@R132]). These complexes are crucial for the development of many organs including skin, as well as for maintenance and regulation of stem cell activity ([@R132]). PRC1 complex contains Ring, Bmi1, Cbx and Phi core subunits and might include other proteins, while PRC2 consist of H3K27 methyltransferases Ezh1 or Ezh2, Suz12, one of the Eed variants, Rbap48 and might also contain additional subunits ([@R6]; [@R132]).

Polycomb complexes exhibit their repressive effects on transcription via at least three mechanisms: 1) The component of PRC2 complex histone methylase Ezh2 promotes tri-methylation of H3K27 followed by targeting the Cbx proteins as a part of PRC1 complex to H3K27me3 and subsequent ubiquitinilation of H2A at lysine 119 catalyzed by the PRC1 component Ring1B ([@R88]). This compact the chromatin and inhibits transcription through preventing binding of transcription machinery to gene promoters; 2) PRC complexes also inhibit transcription elongation via binding to initiating RNA polymerase II that result in a production of the short, but not full length RNA, indicating that PRC complex inhibits their transcription at post-initiation steps ([@R46]); 3) Polycomb complexes control higher-order chromatin remodelling and promote formation of the repressive hubs where repressed genes might be grouping ([@R6]).

It was recently shown that a component of the PRC1 complex Cbx4 plays an important role in the maintenance of the quiescent state in human epidermal progenitor cells ([@R93]). Cbx4 protects them from senescence through PRC-dependent mechanisms by repressing *INK4a* locus, as well as controls their activation and differentiation through PRC-independent mechanisms ([@R93]). Other component of the PRC1 complex Bmi1 is also implicated in the control of keratinocyte senescence and inhibits the expression of p16^INK4A^ and TGF-beta signaling ([@R31]; [@R76]). Bmi-1 is expressed in the basal and suprabasal layers of human epidermis, stimulates keratinocyte proliferation, promotes their survival and protects from apoptosis ([@R40]; [@R87]).

Recent studies also provided important insights in the role of PRC2 complexes in skin development. Conditional inactivation of *Ezh2* in mouse epidermis led to premature epidermal differentiation during morphogenesis ([@R45]). Ezh2 controls proliferative potentials of basal progenitors by repressing the INK4A-INK4B locus and prevents their premature differentiation by inhibiting recruitment of AP1 transcriptional activator to the promoters of terminal differentiation-associated genes including those located in the EDC ([@R45]).

*Ezh2* expression is restricted to basal keratinocytes in embryonic skin, and *Ezh1* is expressed in supra-basal keratinocytes ([@R44]). Whereas single inactivation of *Ezh1* or *Ezh2* in keratinocytes does not lead to any defects in postnatal skin development, their simultaneous inactivation results in the arrest of hair follicle morphogenesis and subsequent degeneration due to defective proliferation and increased apoptosis ([@R44]). These effects were accompanied by complete lack of H3K27me3 in epithelial stem cells and resulted in up-regulation of expression of many non-skin genes, as well as of *INK4a/INK4b* and *Mdm2* genes ([@R44]). Surprisingly, *Ezh1/2* knockout mice did not show any epidermal defects, thus suggesting that loss of their activities in epidermal keratinocytes might be compensated by other PRC2 components ([@R44]).

Repressive effects of Ezh1/2 on transcription are potentiated by a member of Jumonji family of chromatin-remodelling factors Jarid2, which also serves as a component of PRC2 complex ([@R84]). Ablation of Jarid2 in mouse epidermis inhibits proliferation and promotes differentiation of progenitor cells postnatally, without affecting epidermal development. Jarid2 deficiency in keratinocytes reduces H3K27 trimethylation and results in delayed hair follicle cycling as a consequence of decreased proliferation of HF stem cells and their progeny ([@R104]). Thus, similarly to Ezh1/2, Jarid2 is required for the maintenance of cell proliferation and inhibition of differentiation in epidermal stem and progenitor cells ([@R104]).

ATP-dependent chromatin remodeling {#S7}
==================================

ATP-dependent chromatin remodeling complexes alter DNA-histone interactions by using an energy of the ATP hydrolysis ([@R29]). Their action could lead to changes of the nucleosome translational positioning and/or nucleosome conformation and could also facilitate eviction and deposition of nucleosomes ([@R16]). ATP-dependent chromatin remodeling activity resides in multi-subunit complexes, which always contain core ATPase. The additional subunits could modulate the activity of the core ATPase and facilitate targeting of the complexes to the specific genomic regions. Based on their presence, ATP-dependent chromatin remodeling complexes are divided into several (SWI/SNF, ISWI and CHD) groups ([@R29]).

ATPases of the SWI/SNF group contain bromo-domain which binds to acetylated histone tails. In mammals, SWI/SNF group ATPases are represented by two major subunits: Brm (Smarca 2 or Snf2α) and Brg1 (Smarca4 or Snf2β), which form complexes with other subunits called BAFs ([@R29]). The genetics studies demonstrated a crucial role of Brg1 in the development of different organs including skin, while Brm appears to be mostly dispensable ([@R35]; [@R71]; [@R107]). In the developing skin, Brg1 deficiency does not alter keratinocyte proliferation or early differentiation, but lead to the defects in terminal differentiation and formation of non-functional epidermal barrier ([@R71]). qRT-PCR analysis demonstrate decrease of expression of transglutaminase 1/3 genes, as well as of corneodesmosin gene in Brg1-deficient epidermis ([@R71]). Interestingly, inactivation of Brm gene in addition to Brg1 increased severity of epidermal defects in comparison to inactivation of Brg1 gene alone, thus suggesting about a partial redundancy of their functions in keratinocytes ([@R71]).

ATP-dependent chromatin remodelers of the CHD group contain 2 chromo-domains in their N-terminal part ([@R98]). These domains could directly bind to DNA, RNA and methylated histone H3. CHD group proteins are divided into several subgroups based on presence of additional domains together with the chromo-domain. One of such subfamilies includes Chd3 (or Mi-2α) and Chd4 (also known as Mi-2β) proteins that harbor paired PHD (plant homeo domain) domains in their N-terminus involving in chromatin remodeling and binding to methylated histones. Chd3 and Chd4 protein are found in multi-subunit complexes that posses both histone deacetylase and ATP-dependent chromatin remodeling activities. These complexes are called NURD (nucleosome remodeling and histone deacetylase) and contain several proteins usually associated with transcriptional repression: HDAC1 and HDAC2, Rba48 and Rba46, MTA1, MTA2 and MTA3, and MBD3 ([@R36]).

Keratinocyte-specific inactivation of *Chd4* demonstrated that its role in the control of epidermal development is stage-dependent: ablation of *Chd4* during early steps of skin development (E10.5) results in the defects in the basal epidermal layer formation, while loss of *Chd4* during later stages (E13.5) did not cause any significant epidermal defects, but, instead, led to impaired induction and development of the hair follicles ([@R73]). The expression of the hair follicle placode markers, including *Edar,* β*-catenin, Lef1, Shh, Patched1 and Bmp2,* were also consistently reduced in the *Chd4* mutant skin ([@R73]). However, the targets of *Chd4* and other ATP-dependent chromatin remodellers in keratinocytes remain to be identified.

Higher-order chromatin remodeling and 3D genome organization {#S8}
============================================================

The term "higher order chromatin folding or structure" is used to describe spatial organization of chromatin beyond "beads on the string" repeat \[reviewed in ([@R51])\]. Numerous studies demonstrated that chromatin organization in 3D nuclear space is not random. Light microscopic analysis of the nucleus performed over hundred years ago revealed that some chromatin regions are stained with basic dyes brighter and called heterochromatin, while the others are stained scarcer and represent euchromatin. Heterochromatin is more condensed, gene-poor, replicates late during cell cycle, generally associates with low transcription activity and could be sub-divided into constitutive and facultative, while euchromatin is less condensed, gene-rich, replicates relatively early and more transcriptionally active ([@R37]). Constitutive heterochromatin incorporates centromeric, telomeric and other silenced repetitive regions of higher eukaryotic genomes, while facultative heterochromatin contains specific genes and gene loci that could be repressed in some cells and activated in the others in different contexts ([@R37]; [@R95]).

Recent progress in analyses of 3D genome organization has been driven by using two approaches: i) confocal microscopy after fluorescent in situ hybridization (FISH) or gene loci labeling with transgenic fluorescent chimeric proteins containing specific DNA binding domains ([@R72]; [@R80]), and ii) chromatin conformation capture (3C and its variations) technologies based on the restriction of cross-linked chromatin followed by ligation at high dilution to facilitate the formation of intra-molecular but not inter-molecular products ([@R111]; [@R123]). Both approaches revealed that in interphase nucleus chromosomes occupy separate areas (chromosome territories) and genes within distinct chromosomes are non-randomly positioned relatively to each other and other nuclear sub-organelles (reviewed in ([@R33]; [@R111]; [@R123])).

Chromosomes and chromosomal territories {#S9}
---------------------------------------

Chromosomes are the largest units of genome organization occupying distinct territories in the interphase nucleus ([@R32]; [@R33]) (**[Fig. 1](#F1){ref-type="fig"}**). Each chromosome contains a centromere, pericentromeric chromatin enriched in α-satellite repetitive sequences, chromosome arms containing the gene-rich and gene-poor domains enriched in the GC- and AT-sequences and visualized as the light and dark bands by Gimsa staining, respectively, as well as the telomeres ([@R59]).

Although highly variable, the relative positioning of chromosomes in the nucleus is not random, and there is evidence suggested that positioning of the chromosomes could be cell-type specific and also dependent on their size ([@R33]). In cultured human keratinocytes, positions of the chromosomes 18 and 19 are changed during calcium-induced differentiation ([@R97]). However, data obtained from mouse epidermis show that position of the chromosome 3 in keratinocytes of the basal and suprabasal layers is more peripheral compared to the chromosomes 11 and 15 and does not show changes during embryonic and post-natal development (**[Fig. 1 B](#F1){ref-type="fig"}**, our unpublished observations).

Chromosome positioning in post-mitotic nuclei is controlled either through tethering of specific chromatin regions to the nuclear lamina via specialized lamina-associated domains or chromatin-associated proteins (BAF, heterochromatin protein 1, etc.), or through association of the chromosomes bearing the nucleolar-organizing region domains with nucleoli \[reviewed in ([@R72]; [@R101]; [@R106])\]. Several indications suggest that distinct chromosomes may form clusters in the nuclei of differentiated cells in a tissue-specific manner, thus explaining an increased frequency of translocations between the distinct chromosomal parts in the corresponding tumors ([@R18]; [@R74]; [@R113]; [@R122]). However, it is unclear whether non-random distribution of chromosomes may have any functional significance for execution of tissue-specific gene expression programs and whether genes from neighboring chromosomes may share common regulatory mechanisms required for their transcription.

Recently developed approach with using three-dimensional ultra-high resolution light microscopy revealed that each chromosome territory resembles a sponge-like structure and consists of the chromatin domains permeated by interchromatin channels connected with a network of larger channels and lacunas separating distinct chromosomes and harboring a number of nuclear bodies ([@R99]). The network of interchromatin channels starts at nuclear pores and expands throughout the nuclear space, while chromatin domains in each territory are separated from the interchromatin channels by a 100-200-nm layer of decondensed chromatin, called the perichromatin and enriched by nascent DNA and RNA, RNA polymerase II (RNA Pol II), as well as by histone modifications specific for transcriptionally active chromatin ([@R99]).

These observations provide a solid background for the chromosome territory-interchromatin compartment (CT-IC) model ([@R33]) and suggest a large degree of flexibility in the positioning of distinct chromatin domains inside of each chromosome territory. These data also help in better understanding the results published previously and demonstrating the intermingling between the chromatin domains of neighboring chromosome territories ([@R17]), as well as the phenomenon of gene looping outside the chromosome territory and formation of associations in 3D-nuclear space between the genes from different chromosomes ([@R125]). However, functional significance of such associations and mechanisms of their formation remain mostly unknown.

Transcription factories, nuclear speckles and transcriptionally-active chromatin domains {#S10}
----------------------------------------------------------------------------------------

It was shown that actively transcribed genes tend to be associated closely to each other and form active domains within chromosomal territories, whereas inactive genes are also clustered to form silenced chromosomal domains ([@R90]). Moreover, many indications suggest that gene activation or silencing are associated with changes of their positioning in the nucleus relatively to neighboring genes, distinct chromosomal domains and/or nuclear bodies ([@R41]; [@R55]; [@R106]; [@R125]). The dynamic changes in higher-order chromatin structure associated with cell differentiation are quite consistent with data obtained with super-resolution confocal microscopy that revealed a network of channels and lacunas inside of both euchromatin and heterochromatin, thus implicating a large degree of flexibility in the positioning of genes and chromosomal domains relatively to each other ([@R99]) (**[Fig. 1 A](#F1){ref-type="fig"}**).

Many indications suggest that actively transcribed genes are associated with transcription factories visualized in the nucleus by immunostaining to phosphorylated at Serine-5 and/or Serine-2 RNA Polymerase II isoforms that depict the initiating or elongating forms of enzyme, respectively ([@R30]; [@R126]). It was proposed that several actively transcribed genes might share single transcription factory, which is formed through the aggregation of RNA Pol II-containing pre-initiation complexes assembled next to each other in the nuclear space ([@R118]). However, super-resolution confocal microscopy revealed that single PolII sites are also frequently seen in the nucleus, indicating that transcription might also occur outside of such factories ([@R99]).

There is evidence that distinct transcription factories might contain unique combinations of transcription factors and/or other molecules facilitating transcription of co-regulated genes from the same or different chromosomes ([@R43]). Data obtained from erythroid cells demonstrate that tissue-specific Hbα/β globin genes are co-localized with hundreds of other transcribed genes from different chromosomes in the transcription factories enriched by Klf1 transcription factor that promote erythroid differentiation ([@R126]). This supports an idea that distinct transcription factories might serve as hubs for inter- and intra-chromosomal interactions between many genes to achieve an efficient and coordinated transcriptional control required for proper execution of cell type-specific differentiation programs ([@R25]).

During development, many lineage-specific genes change their positioning and move away from the nuclear periphery towards nuclear interior \[reviewed in ([@R41])\]. Several examples (selected Hox genes in embryonic stem cells, Shh gene in the developing limb bud, etc.) demonstrate that some genes and gene loci form loops and become localized outside of the corresponding chromosomal territories after activation or in highly active state ([@R2]; [@R11]; [@R51]). In human primary keratinocytes, EDC locus also shows looping outside of the chromosomal territory 1 ([@R145]). However, many active genes are also found inside or on the border of the chromosome territories, suggesting that repositioning might not be essential for maintenance of their expression ([@R81]; [@R94]).

Developmentally-regulated gene repositioning might be explained by the necessity for selected genes to reach a permissive nuclear environment required for maintenance or co-regulation of their expression ([@R83]). Nuclear speckles localized in the interchromatin compartments and enriched by the components of the RNA splicing machinery might provide such permissive environment, and, indeed, many highly-transcribed genes including the EDC locus in keratinocytes show a tendency to cluster around speckles ([@R19]; [@R69]; [@R133]; [@R135]) (**[Fig. 1 A](#F1){ref-type="fig"}**). Recent data demonstrate that relocation of the proliferation-associated genes from the repressive Polycomb bodies to nuclear speckles depends on methylation/demethylation of the PRC1 protein Cbx4: methylated form of Cbx4 binds to the non-coding RNA TUG1, while demethylated Cbx4 binds to non-coding RNA MALAT1/NEAT2, located in the Polycomb bodies or speckles, respectively ([@R149]). Repositioning of genes to speckles results in their activation through MALAT1/NEAT2-mediated interactions with multiple co-activators of transcription and/or splicing factors ([@R149]). Frequent associations between EDC and nuclear speckles are also seen in the epidermal keratinocytes (**[Fig. 1 B](#F1){ref-type="fig"}**). However, functional significance of keratinocyte-specific gene loci co-localization with nuclear speckles remains to be determined.

Long-range chromatin interactions and establishing specific chromatin conformations in tissue-specific genomic loci {#S11}
-------------------------------------------------------------------------------------------------------------------

Increased evidence of data suggest that inter- and intra-chromosomal long-range associations between distinct genomic regions are essential for maintenance of cell type-specific transcriptional profiles and form "transcriptional interactome", which is quite dynamic in its nature ([@R125]). In a number of tissue-specific gene loci (beta-globin locus, T~H~2-cytokine locus, etc.), long-range interactions help to establish specific loop structures, in which distant locus control elements are brought in close proximity to gene promoters to regulate transcription, whereas alterations in spatial organization of such loop structures significantly impact gene expression programs in differentiating cells ([@R41]). In differentiating keratinocytes, calcium stimulation result in increased physical proximity of the enhancer and the promoter regions of the peptidylarginine deaminase 3 gene that control metabolism of the filaggrin ([@R26]).

Repositioning of the genes and long-range interactions between genomic loci are energy-dependent processes and are regulated, at least in part, by actin-myosin nuclear motor complexes ([@R28]; [@R69]; [@R102]). Formation of the chromatin loop-like structures inside of the chromosomes (in cis) and between different chromosomes (in trans) is also controlled by several proteins that function as topological genome organizers, such as CTCF, cohesin and SATB1 ([@R23]; [@R60]; [@R105]; [@R112]). CTCF is DNA-binding zinc-finger protein implicated in the control of many important processes in the nucleus ([@R112]). CTCF binding sites are spread all over the genome, thus providing the possibility to form thousands of inter-chromosomal loops of different size, as well as to bring genes located on different chromosomes into close proximity ([@R112]; [@R152]). Recent studies revealed that CTCF recruits the mitosis regulator cohesin, which is essential for proper chromosome segregation, to specific chromosomal locations in the interphase nucleus ([@R63]). In T-cell receptor locus, CTCF-cohesin complex is required for long-range promoter-enhancer interactions, maintenance of transcription and demarcation of the locus from interspersed elements and neighbouring housekeeping genes ([@R128]). Thus, in addition to its global role in mitosis regulation, cohesin controls chromosome conformation and gene expression ([@R105]).

AT-rich binding protein Satb1 also controls establishing specific three-dimensional conformations in tissue-specific gene loci ([@R22]; [@R23]). Satb1 targets chromatin remodelling enzymes and transcription factors to specific genomic regions, establishes region-specific epigenomic modification status, and plays a fundamental role in the execution of tissue-specific gene expression programs ([@R23]). In epidermal keratinocytes, Satb1 regulates conformation of the central domain of the EDC locus containing genes activated during terminal keratinocyte differentiation ([@R52]). Ablation of Satb1 results in alterations of 3D-chromatin structure and elongation of this domain, accompanied by marked alterations of gene expression and epidermal morphology ([@R52]). This suggests that three-dimensional organization of multi-gene tissue-specific loci is indeed important for maintenance of the proper balance of gene expression and outcomes of the transcription programs (**[Fig. 2 C](#F2){ref-type="fig"}**).

Heterochromatin, transcriptionally-repressive domains and gene silencing {#S12}
------------------------------------------------------------------------

In almost all mammalian somatic cells, heterochromatin is associated with nuclear membrane, nucleoli, and centromeric repeat clusters, whereas small regions of facultative heterochromatin are also incorporated into euchromatin ([@R37]). Key features of heterochromatin are its associations with repressive epigenetic marks including DNA methylation and H3K9 tri-methylation, as well as with Heterochromatin Protein 1 (HP1) that interacts with a number of other proteins involved in the control of chromatin remodeling ([@R82]; [@R116]).

Currently, three distinct types of transcriptionally inactive heterochromatin domains are described: 1) Domains associated with centromeres and marked by DNA methylation, H3K9me3, H3K20me3 and α/β/γ HP1 ιsoforms; 2) Chromatin associated with nuclear lamina via Lamina-associated domains that enriched by H3K9me2 and HP1α/β; 3) Chromatin marked by H3K27me3 and Polycomb group proteins ([@R103]; [@R116]). Although all three types of heterochromatin domains contain actively transcribed genes, transcription of the majority of genes associated with these domains is repressed either constitutively (at the centromere regions) or temporarily in the developmentally-regulated or context-dependent manners (lamina-associated and Polycomb-repressed genes) ([@R103]).

Heterochromatin Protein 1 plays a critical role in organization and formation of heterochromatin ([@R82]). In mammals, HP1 consists of three isoforms (α, β, γ), encoded by three different chromobox genes (Cbx5, Cbx1 and Cbx3, respectively), which are quite distinct in their functions: HP1α/β bind H3K9me3 and target HP1 to heterochromatin, whereas HP1γ via histone chaperone FACT binds to elongating form of PolII and interacts with euchromatin ([@R82]). HP1α/β interactions with H3K9me3 is important for condensing pericentric chromatin regions to maintain transcriptional repression, as well as for binding Lamin B receptor and maintenance of lamina-associated heterochromatin ([@R13]). In contrast to HP1α/β, HP1γ is concentrated in 3'-end of active genes and might positively regulate transcription ([@R82]).

Polycomb proteins also play crucial role in maintenance of the chromatin repressed states via interaction with Polycomb response elements and gene promoters and formation of topologically complex structures (repressive chromatin hubs) that promotes gene silencing ([@R85]). Repressive chromatin loops might be organized in 3D nuclear space into Polycomb bodies that are considered as chromatin silencing factories and contain a number of transcriptional repressors including non-coding RNA TUG1 ([@R7]; [@R68]; [@R96]; [@R149]). However, it is unclear whether such factories exist in keratinocytes and how they contribute to terminal differentiation-associated gene silencing in the epidermis remain to be clarified.

Cross-talk between the epigenetic and transcription factor-mediated regulatory mechanisms in the control of keratinocyte differentiation {#S13}
========================================================================================================================================

The program of epidermal differentiation and barrier formation is tightly controlled by a number of transcription regulators including the p63, c-Myc, AP-2, Klf4, Arnt, PPAR-alpha, m-Ovo, etc. ([@R14]; [@R61]; [@R79]; [@R138]). The p63 transcription factor serves as a master regulator of epidermal development, and p63 knockout (-/-) mice fail to form stratified epithelium and to express a number of epidermis-specific genes \[reviewed in ([@R79]; [@R139])\]. In keratinocytes, p63 regulates a large number of genes that encode distinct adhesion/signalling molecules, transcription factors, cell cycle-associated proteins, as well as tissue-specific proteins, such as keratins, involucrin, and loricrin \[reviewed in ([@R140])\].

Recent data revealed that in the epidermis of p63-null mice the expression of a large number of genes involved in the control of nuclear/chromatin assembly and remodeling is significantly changed compared to wild-type controls ([@R52]). These genes include several key regulators of the chromatin structure and remodeling: histone methytransferases Ehmt1/2, components of the Polycomb complexes Cbx2/4/6/7 and Jarid2, ATP-dependent chromatin remodellers Mi-2α/β, Brg1, Chd1 and Brm, as well as genome organizer Satb1 ([@R52]). Whereas involvement of p63 in regulation of the majority of these genes remain to be clarified, it is shown that p63 directly controls the expression of Satb1 and promotes establishing of specific higher-order chromatin structure in central domain of the EDC required for maintenance of the proper balance of gene expression during terminal keratinocyte differentiation ([@R52]). These data demonstrate that Satb1 is an important, hitherto unrecognized component of the global regulatory network of p63, and suggest a novel fundamental mechanism on how master regulators of tissue morphogenesis establish tissue-specific patterns of higher-order chromatin structure and remodeling to control the fate of multi-potent progenitor cells.

p63 expression in keratinocytes is regulated by histone methyltransferase Setd8, which, in turn, is a target of c-Myc transcription factor and mediates its effects on epidermal differentiation ([@R38]). c-Myc is an important regulator of epidermal homeostasis and promotes keratinocyte proliferation and differentiation into epidermal and sebaceous gland lineages ([@R92]; [@R110]; [@R141]). In addition to direct involvement in the control of expression of epidermis-specific genes, c-Myc also exhibits a potent activity as a modulator of local chromatin structure: it regulates expression of histone acetyltransferase Gcn5 and promotes histone acetylation ([@R78]). Furthermore, c-Myc regulates expression of RNA methyltransferase Misu, which modulate c-Myc effects on cell proliferation and differentiation ([@R38]; [@R57]).

Transcription factors also directly interact with chromatin remodeling factors and target them to gene promoters: c-Myc targets histone acetyltransferases Gcn5 to chromatin and promotes global euchromatisation ([@R78]), whereas NFATc1 transcription factor forms complex with ATP-dependent chromatin remodeler Brg1 and stimulates formation of DNAse1 hypersensitive sites and recruitment of other transcription factors to target genes ([@R114]; [@R147]). However, transcription factors might also operate as the repressors of transcription by interacting with HDACs and targeting repressive complexes to gene promoters: for instance, transcriptional repressor NFX1 interacts with mSin3A/HDAC and recruits the repression complex to the hTERT promoter in keratinocytes, while p53 inhibits expression of its target genes during DNA damage response via interaction and targeting Sin3B/HDAC complex to the promoters ([@R5]; [@R148]).

Thus, chromatin remodeling factors operate as integral part of the genetic programs governed by transcription factors and modulate their effects on local or global chromatin structure (**[Fig. 3](#F3){ref-type="fig"}**). Moreover, transcription factors might directly interact with chromatin remodelers to target them to the specific genomic sites and regulate gene expression. However, chromatin remodeling factors might, in turn, regulate expression of the distinct transcription factors via formation of the active or repressive local chromatin structure at their promoter regions.

Conclusions and future directions {#S14}
=================================

In summary, epigenetic regulators modulate both local and higher-order chromatin structure in epidermal keratinocytes and exhibit their effects on cell proliferation and differentiation via at least three mechanisms: 1) Controlling the expression of cell cycle-associated genes; 2) Regulating expression of the terminal differentiation-associated genes; 3) Modulating gene expression programs governed by distinct transcription factors, such as p63, c-Myc, NF-kappaB, etc.

Group of factors that predominantly operate as repressive chromatin regulators (DNA methyltransferase DNMT1, histone deacetylases HDAC1/2, Polycomp components Cbx4, Bmi1, Ezh1/2) stimulate proliferation of progenitor cells via repression of the genes encoding cell-cycle inhibitors p16^INK4A^, p15^INK4B^, p19^ARF^ and p21. Some of the repressive chromatin regulators (DNMT1, Ezh1/2, Jarid2) also inhibit premature activation of terminal differentiation-associated genes, whereas HDAC2 inhibits expression of deltaNp63 and interfere with deltaNp63-regulated gene expression program in keratinocytes.

Vice versa, a group of chromatin remodelers that generally *support active transcription* (histone demethylase Jmjd3, ATP-dependent chromatin remodeler Brg1, genome organizer Satb1) promote terminal keratinocyte differentiation and also exhibit effects on cell proliferation. It is likely that expression of terminal differentiation-associated genes including those that constitute EDC in epidermal keratinocytes requires tightly controlled coordination of action of several epigenetic regulators operating at the levels of both local (DNMT1, Ezh1/2, Bmi1, Jarid2, Jmjd3, Brg1) and higher-order (Satb1) chromatin structures (**[Fig. 2 C](#F2){ref-type="fig"}; [Fig. 3](#F3){ref-type="fig"}**).

Epigenetic modifications play an important role in the development of many disorders, such as cancer, immunodeficiency, and autoimmune diseases ([@R12]; [@R48]). Epigenetic status of skin epithelial and non-epithelial cells is altered in various types of skin cancers and psoriasis ([@R20]; [@R27]; [@R89]; [@R108]). Comprehensive knowledge of epigenetic mechanisms and their alterations in health and disease, has become a growing priority in biomedical research ([@R115]). However, additional efforts are required to fully understand how epigenetic mechanisms contribute to interpret different environmental insults applied to the skin, as well as how these mechanisms are altered in skin disorders.

In particular, at least the following important issues need to be resolved: Cross-talk between distinct signaling pathways that control keratinocyte proliferation/differentiation/migration and epigenetic regulatory mechanisms need to be investigated;Mechanisms that control specificity of distinct classes of epigenetic regulators in targeting different genes/gene loci in the genome of keratinocytes remain unclear;Roles of the non-coding and microRNAs and their interaction with epigenetic regulators in modulating gene expression programs need to be clarified;Topological organization of the keratinocyte genome, links between the higher-order and local chromatin organization and the role of each nuclear sub-organelle in keratinocyte differentiation process, as well as in the response to various environmental stimuli need to be carefully studied;Roles of distinct classes of epigenetic regulators in the pathological skin conditions associated with alterations of keratinocyte differentiation, their expansion or loss, as well as during DNA damage response remain to be further defined.

Recent advances of pharmacogenomics resulted in the development of a number of molecules that are capable of modulating distinct epigenetic regulatory mechanisms, whereas some of them are already under clinical trials \[reviewed in ([@R65])\]. Future research in this direction will help to bridge the gap between our current knowledge of basic epigenetic mechanisms and potential applications of epigenetic modulators, which hopefully will serve as new generation of drugs for treatment of skin disorders, protection of skin against environmental stressors and aging.

Work on this review was supported from the grants from Medical Research Council (UK) and NIH to V.A.B. and via K01 Mentored Scientist Award to A.A.S. Critical comments on the manuscript from Dr. Barbara A. Gilchrest and Dr. Chris Lewis are gratefully acknowledged.

The authors state no conflict of interest.

DNMT

:   DNA methyltransferase

EDC

:   Epidermal Differentiation Complex

HDAC

:   Histone deacetylase

Pol II

:   RNA Polymerase II

![Schematic structure of interphase nucleus\
**A** - Nucleus is surrounded by nuclear envelope that consists of outer and inner membranes, nuclear pore complexes and lamina located beneath the inner nuclear membrane. Chromosomes occupy distinct territories, in which distinct chromatin domains are permeated by interchromatin channels connected with a network of larger channels and lacunas separating distinct chromosomes and harboring a number of nuclear bodies including speckles (inset). Nucleus also contains a number of nuclear bodies including nucleolus, PML bodies, speckles and Polycomb bodies.\
**B**: Images of nuclei of murine basal epidermal keratinocytes showing chromosome territory 11 (left image, large arrow) with keratin 15 gene (small arrow), as well as Epidermal Differentiation Complex locus (right image, red color) contacting SC-35-positive speckle (inset, green color, arrow).](nihms-373606-f0001){#F1}

![Distinct levels of chromatin organization, characteristics of the active and inactive chromatin domains and epigenetic control of gene expression in the Epidermal Differentiation Complex locus\
**A** -- Levels of chromatin folding and size of the corresponding chromatin domains of the interphase chromosome. The lowest unit of chromatin organization is nucleosome, in which DNA form 1 and ¾ loops around the histone octamer core. Nucleosomes are connected by linker DNA and form the "beads on the string" structures, which are folded into 30 nm fiber followed by further folding into higher-order loop structures (modified from [@R49]).\
**B** -- Distribution of different epigenetic marks among the distinct regions of actively transcribed, paused and inactive genes. Actively transcribed genes are characterized by high levels of Polymerase II occupancy, low levels of DNA methylation, high levels of H3K4me3 and H3K9acetyl/H3K18acetyl at the transcription start site, as well as by high levels of H3K79me2/3 and H3K36me3 at the gene body region. Genes with paused transcription show high levels of Polymerase II occupancy, low levels of DNA methylation, high levels of H3K4me3 and H3K27me2/3 at the transcription start site. Inactive genes show lack of Polymerase II, high levels of DNA methylation and H3K9me2/3 at the transcription start site (modified from [@R143]).\
**C** -- Epigenetic control of expression of terminal differentiation-associated genes in the Epidermal Differentiation Complex locus occupying about 3 Mb in mouse chromosome 3. Scheme of the locus shows chromosomal localization of the distinct genes and gene families. Genes activated during terminal keratinocyte differentiation constitute a central domain of the locus. Active transcription in the central domain is maintained by coordinated involvement of the regulators of local chromatin structure (ATP-dependent chromatin remodelers Brg1 and Mi-2β, histone demethylase JMJD3) and genome organizers Satb1 that support formation of specific higher-order chromatin folding in this region. Inhibition of transcription is associated with presence of DNA methyltransferase 1, HDAC1/2 and components of the Polycomb complexes Cbx4, Ezh1/2 and Bmi1, which promote formation of the repressive chromatin structure, as well as with lack of specific 3D-organization of the locus.](nihms-373606-f0002){#F2}

![Model illustrating cross-talk between transcription factors and epigenetic regulators in the control of expression of tissue-specific genes in keratinocytes\
In addition to direct regulation of tissue-specific genes, p63 and c-Myc transcription factors promote formation of the permissive chromatin structure in tissue-specific gene loci. p63 via regulation of genome organizer Satb1 support formation of specific higher-order organization of the tissue-specific EDC locus, whereas c-Myc via regulation of histone acetyltrasferase Gcn5 and histone methyltransferase Setd8 promotes formation of the transcription-friendly local chromatin structure. In turn, Setd8 is capable of positively influence p63 expression (based on data published ([@R38]; [@R52]; [@R78]; [@R110]).](nihms-373606-f0003){#F3}

[^1]: This review is dedicated to the 90th Anniversary of the Corresponding Member of Russian Academy of Natural Sciences Prof. Dina S. Gordon, whose wisdom, knowledge and encouragement supported V.A.B., A.A.S. and A.N.M. in early stages of their scientific careers.
