Single Cell and Open Chromatin Analysis Reveals Molecular Origin of Epidermal Cells of the Skin

Molecular and Cellular Function of p63 in Skin Development and Genetic Diseases

The transcription factor p63 is a master regulator of multiple ectodermal derivatives. During epidermal commitment, p63 interacts with several chromatin remodeling complexes to transactivate epidermal-specific genes and repress transcription of simple epithelial and nonepithelial genes. In the postnatal epidermis, p63 is required to control the proliferative potential of progenitor cells, maintain epidermal integrity, and contribute to epidermal differentiation. Autosomal dominant sequence variant in p63 cause a spectrum of syndromic disorders that affect several tissues, including or derived from stratified epithelia. In this review, we describe the recent studies that have provided novel insights into disease pathogenesis and potential therapeutic targets.

ΔNp63 regulates Sfrp1 expression to direct salivary gland branching morphogenesis

Branching morphogenesis is a complex process shared by many organs including the lungs, kidney, prostate, as well as several exocrine organs including the salivary, mammary and lacrimal glands. This critical developmental program ensures the expansion of an organ's surface area thereby maximizing processes of cellular secretion or absorption. It is guided by reciprocal signaling from the epithelial and mesenchymal cells. While signaling pathways driving salivary gland branching morphogenesis have been relatively well-studied, our understanding of the underlying transcriptional regulatory mechanisms directing this program, is limited. Here, we performed in vivo and ex vivo studies of the embryonic mouse submandibular gland to determine the function of the transcription factor ΔNp63, in directing branching morphogenesis. Our studies show that loss of ΔNp63 results in alterations in the differentiation program of the ductal cells which is accompanied by a dramatic reduction in branching morphogenesis that is mediated by dysregulation of WNT signaling. We show that ΔNp63 modulates WNT signaling to promote branching morphogenesis by directly regulating Sfrp1 expression. Collectively, our findings have revealed a novel role for ΔNp63 in the regulation of this critical process and offers a better understanding of the transcriptional networks involved in branching morphogenesis.

ΔNp63 Regulates Homeostasis, Stemness, and Suppression of Inflammation in the Adult Epidermis

The p63 transcription factor is critical for epidermis formation in embryonic development, but its role in the adult epidermis is poorly understood. In this study, we show that acute genetic ablation of ΔNp63, the main p63 isoform, in adult epidermis disrupts keratinocyte proliferation and self-maintenance and, unexpectedly, triggers an inflammatory psoriasis-like condition. Mechanistically, single-cell RNA sequencing revealed the downregulation of cell cycle genes, upregulation of differentiation markers, and induction of several proinflammatory pathways in ΔNp63-ablated keratinocytes. Intriguingly, ΔNp63-ablated cells disappear by 3 weeks after ablation, at the expense of the remaining nonablated cells. This is not associated with active cell death and is likely due to reduced self-maintenance and enhanced differentiation. Indeed, in vivo wound healing, a physiological readout of the epidermal stem cell function, is severely impaired upon ΔNp63 ablation. We found that the Wnt signaling pathway (Wnt10A, Fzd6, Fzd10) and the activator protein 1 (JunB, Fos, FosB) factors are the likely ΔNp63 effectors responsible for keratinocyte proliferation/stemness and suppression of differentiation, respectively, whereas IL-1a, IL-18, IL-24, and IL-36γ are the likely negative effectors responsible for suppression of inflammation. These data establish ΔNp63 as a critical node that coordinates epidermal homeostasis, stemness, and suppression of inflammation, upstream of known regulatory pathways.

Molecular and spatial landmarks of early mouse skin development

A wealth of specialized cell populations within the skin facilitates its hair-producing, protective, sensory, and thermoregulatory functions. How the vast cell-type diversity and tissue architecture develops is largely unexplored. Here, with single-cell transcriptomics, spatial cell-type assignment, and cell-lineage tracing, we deconstruct early embryonic mouse skin during the key transitions from seemingly uniform developmental precursor states to a multilayered, multilineage epithelium, and complex dermal identity. We identify the spatiotemporal emergence of hair-follicle-inducing, muscle-supportive, and fascia-forming fibroblasts. We also demonstrate the formation of the panniculus carnosus muscle (PCM), sprouting blood vessels without pericyte coverage, and the earliest residence of mast and dendritic immune cells in skin. Finally, we identify an unexpected epithelial heterogeneity within the early single-layered epidermis and a signaling-rich periderm layer. Overall, this cellular and molecular blueprint of early skin development-which can be explored at https://kasperlab.org/tools-establishes histological landmarks and highlights unprecedented dynamic interactions among skin cells.

MOF-mediated histone H4 Lysine 16 acetylation governs mitochondrial and ciliary functions by controlling gene promoters

Histone H4 lysine 16 acetylation (H4K16ac), governed by the histone acetyltransferase MOF, orchestrates gene expression regulation and chromatin interaction. However, the roles of MOF and H4K16ac in controlling cellular function and regulating mammalian tissue development remain unclear. Here we show that conditional deletion of Mof in the skin, but not Kansl1, causes severe defects in the self-renewal of basal epithelial progenitors, epidermal differentiation, and hair follicle growth, resulting in barrier defects and perinatal lethality. MOF-regulated genes are highly enriched for essential functions in the mitochondria and cilia. Genetic deletion of Uqcrq, an essential subunit for the electron transport chain (ETC) Complex III, in the skin, recapitulates the defects in epidermal differentiation and hair follicle growth observed in MOF knockout mouse. Together, this study reveals the requirement of MOF-mediated epigenetic mechanism for regulating mitochondrial and ciliary gene expression and underscores the important function of the MOF/ETC axis for mammalian skin development.

Single cell analysis of transcriptome and open chromatin reveals the dynamics of hair follicle stem cell aging

Aging is defined as the functional decline of tissues and organisms, leading to many human conditions, such as cancer, neurodegenerative diseases, and hair loss. Although stem cell exhaustion is widely recognized as a hallmark of aging, our understanding of cell state changes-specifically, the dynamics of the transcriptome and open chromatin landscape, and their relationship with aging-remains incomplete. Here we present a longitudinal, single-cell atlas of the transcriptome and open chromatin landscape for epithelia cells of the skin across various hair cycle stages and ages in mice. Our findings reveal fluctuating hair follicle stem cell (HF-SC) states, some of which are associated with the progression of the hair cycle during aging. Conversely, inner bulge niche cells display a more linear progression, seemingly less affected by the hair cycle. Further analysis of the open chromatin landscape, determined by single-cell Assay for Transposase-Accessible Chromatin (ATAC) sequencing, demonstrates that reduced open chromatin regions in HF-SCs are associated with differentiation, whereas gained open chromatin regions in HF-SCs are linked to the transcriptional control of quiescence. These findings enhance our understanding of the transcriptional dynamics in HF-SC aging and lay the molecular groundwork for investigating and potentially reversing the aging process in future experimental studies.

Wnt/β-catenin signaling controls mouse eyelid growth by mediating epithelial-mesenchymal interactions

PURPOSE

To investigate the role of Wnt/β-catenin signaling in mouse eyelid development.

METHODS

Wnt/β-catenin signaling was disrupted by deleting supraorbital mesenchymal β-catenin or epithelial Wls. p63 was removed to determine whether the expression of Wnts is affected. The eyelid morphology was examined at different stages. Proliferation, apoptosis, and expression of Wnt ligands and their target genes were analyzed via immunofluorescence staining, TUNEL assay, and in situ hybridization.

RESULTS

Deletion of β-catenin in supraorbital mesenchyme abolishes eyelid growth by causing decreased proliferation in supraorbital epithelium and underlying mesenchyme. Inhibition of Wnt secretion by deleting Wls in supraorbital epithelium results in failure of eyelid development, similar to the effects of deleting mesenchymal β-catenin. Knockout of p63 results in formation of hypoplastic eyelids and reduced expression of several Wnt ligands in eyelid epithelium.

CONCLUSIONS

Epithelial Wnt ligands activate mesenchymal Wnt/β-catenin signaling to control eyelid growth and their expression is partially regulated by p63.

MicroRNA-205 promotes hair regeneration by modulating mechanical properties of hair follicle stem cells

Stiffness and actomyosin contractility are intrinsic mechanical properties of animal cells required for the shaping of tissues. However, whether tissue stem cells (SCs) and progenitors located within SC niche have different mechanical properties that modulate their size and function remains unclear. Here, we show that hair follicle SCs in the bulge are stiff with high actomyosin contractility and resistant to size change, whereas hair germ (HG) progenitors are soft and periodically enlarge and contract during quiescence. During activation of hair follicle growth, HGs reduce contraction and more frequently enlarge, a process that is associated with weakening of the actomyosin network, nuclear YAP accumulation, and cell cycle reentry. Induction of miR-205, a novel regulator of the actomyosin cytoskeleton, reduces actomyosin contractility and activates hair regeneration in young and old mice. This study reveals the control of tissue SC size and activities by spatiotemporally compartmentalized mechanical properties and demonstrates the possibility to stimulate tissue regeneration by fine-tuning cell mechanics.

Live imaging reveals chromatin compaction transitions and dynamic transcriptional bursting during stem cell differentiation in vivo

Stem cell differentiation requires dramatic changes in gene expression and global remodeling of chromatin architecture. How and when chromatin remodels relative to the transcriptional, behavioral, and morphological changes during differentiation remain unclear, particularly in an intact tissue context. Here, we develop a quantitative pipeline which leverages fluorescently-tagged histones and longitudinal imaging to track large-scale chromatin compaction changes within individual cells in a live mouse. Applying this pipeline to epidermal stem cells, we reveal that cell-to-cell chromatin compaction heterogeneity within the stem cell compartment emerges independent of cell cycle status, and instead is reflective of differentiation status. Chromatin compaction state gradually transitions over days as differentiating cells exit the stem cell compartment. Moreover, establishing live imaging of Keratin-10 (K10) nascent RNA, which marks the onset of stem cell differentiation, we find that Keratin-10 transcription is highly dynamic and largely precedes the global chromatin compaction changes associated with differentiation. Together, these analyses reveal that stem cell differentiation involves dynamic transcriptional states and gradual chromatin rearrangement.

Differentiation of pluripotent stem cells for modeling human skin development and potential applications

The skin of mammals is a multilayered and multicellular tissue that forms an environmental barrier with key functions in protection, regulation, and sensation. While animal models have long served to study the basic functions of the skin in vivo, new insights are expected from in vitro models of human skin development. Human pluripotent stem cells (PSCs) have proven to be invaluable tools for studying human development in vitro. To understand the mechanisms regulating human skin homeostasis and injury repair at the molecular level, recent efforts aim to differentiate PSCs towards skin epidermal keratinocytes, dermal fibroblasts, and skin appendages such as hair follicles and sebaceous glands. Here, we present an overview of the literature describing strategies for human PSC differentiation towards the components of skin, with a particular focus on keratinocytes. We highlight fundamental advances in the field employing patient-derived human induced PSCs (iPSCs) and skin organoid generation. Importantly, PSCs allow researchers to model inherited skin diseases in the search for potential treatments. Skin differentiation from human PSCs holds the potential to clarify human skin biology.

ID1 and CEBPA coordinate epidermal progenitor cell differentiation

The regulatory circuits that coordinate epidermal differentiation during development are still not fully understood. Here, we report that the transcriptional regulator ID1 is enriched in mouse basal epidermal progenitor cells and find ID1 expression to be diminished upon differentiation. In utero silencing of Id1 impairs progenitor cell proliferation, leads to precocious delamination of targeted progenitor cells and enables differentiated keratinocytes to retain progenitor markers and characteristics. Transcriptional profiling suggests that ID1 acts by mediating adhesion to the basement membrane while inhibiting spinous layer differentiation. Co-immunoprecipitation reveals ID1 binding to transcriptional regulators of the class I bHLH family. We localize bHLH Tcf3, Tcf4 and Tcf12 to epidermal progenitor cells during epidermal stratification and establish TCF3 as a downstream effector of ID1-mediated epidermal proliferation. Finally, we identify crosstalk between CEBPA, a known mediator of epidermal differentiation, and Id1, and demonstrate that CEBPA antagonizes BMP-induced activation of Id1. Our work establishes ID1 as a key coordinator of epidermal development, acting to balance progenitor proliferation with differentiation and unveils how functional crosstalk between CEBPA and Id1 orchestrates epidermal lineage progression.

KLF5 governs sphingolipid metabolism and barrier function of the skin

Stem cells are fundamental units of tissue remodeling whose functions are dictated by lineage-specific transcription factors. Home to epidermal stem cells and their upward-stratifying progenies, skin relies on its secretory functions to form the outermost protective barrier, of which a transcriptional orchestrator has been elusive. KLF5 is a Krüppel-like transcription factor broadly involved in development and regeneration whose lineage specificity, if any, remains unclear. Here we report KLF5 specifically marks the epidermis, and its deletion leads to skin barrier dysfunction in vivo. Lipid envelopes and secretory lamellar bodies are defective in KLF5-deficient skin, accompanied by preferential loss of complex sphingolipids. KLF5 binds to and transcriptionally regulates genes encoding rate-limiting sphingolipid metabolism enzymes. Remarkably, skin barrier defects elicited by KLF5 ablation can be rescued by dietary interventions. Finally, we found that KLF5 is widely suppressed in human diseases with disrupted epidermal secretion, and its regulation of sphingolipid metabolism is conserved in human skin. Altogether, we established KLF5 as a disease-relevant transcription factor governing sphingolipid metabolism and barrier function in the skin, likely representing a long-sought secretory lineage-defining factor across tissue types.

Revisiting the embryogenesis of lip and palate development

Clefts of the lip and palate (CLP), the major causes of congenital facial malformation globally, result from failure of fusion of the facial processes during embryogenesis. With a prevalence of 1 in 500-2500 live births, CLP causes major morbidity throughout life as a result of problems with facial appearance, feeding, speaking, obstructive apnoea, hearing and social adjustment and requires complex, multi-disciplinary care at considerable cost to healthcare systems worldwide. Long-term outcomes for affected individuals include increased mortality compared with their unaffected siblings. The frequent occurrence and major healthcare burden imposed by CLP highlight the importance of dissecting the molecular mechanisms driving facial development. Identification of the genetic mutations underlying syndromic forms of CLP, where CLP occurs in association with non-cleft clinical features, allied to developmental studies using appropriate animal models is central to our understanding of the molecular events underlying development of the lip and palate and, ultimately, how these are disturbed in CLP.

A high-resolution route map reveals distinct stages of chondrocyte dedifferentiation for cartilage regeneration

Articular cartilage damage is a universal health problem. Despite recent progress, chondrocyte dedifferentiation has severely compromised the clinical outcomes of cell-based cartilage regeneration. Loss-of-function changes are frequently observed in chondrocyte expansion and other pathological conditions, but the characteristics and intermediate molecular mechanisms remain unclear. In this study, we demonstrate a time-lapse atlas of chondrocyte dedifferentiation to provide molecular details and informative biomarkers associated with clinical chondrocyte evaluation. We performed various assays, such as single-cell RNA sequencing (scRNA-seq), live-cell metabolic assays, and assays for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq), to develop a biphasic dedifferentiation model consisting of early and late dedifferentiation stages. Early-stage chondrocytes exhibited a glycolytic phenotype with increased expression of genes involved in metabolism and antioxidation, whereas late-stage chondrocytes exhibited ultrastructural changes involving mitochondrial damage and stress-associated chromatin remodeling. Using the chemical inhibitor BTB06584, we revealed that early and late dedifferentiated chondrocytes possessed distinct recovery potentials from functional phenotype loss. Notably, this two-stage transition was also validated in human chondrocytes. An image-based approach was established for clinical use to efficiently predict chondrocyte plasticity using stage-specific biomarkers. Overall, this study lays a foundation to improve the quality of chondrocytes in clinical use and provides deep insights into chondrocyte dedifferentiation.

AP-2α and AP-2β cooperatively function in the craniofacial surface ectoderm to regulate chromatin and gene expression dynamics during facial development

The facial surface ectoderm is essential for normal development of the underlying cranial neural crest cell populations, providing signals that direct appropriate growth, patterning, and morphogenesis. Despite the importance of the ectoderm as a signaling center, the molecular cues and genetic programs implemented within this tissue are understudied. Here, we show that removal of two members of the AP-2 transcription factor family, AP-2α and AP-2ß, within the early embryonic ectoderm of the mouse leads to major alterations in the craniofacial complex. Significantly, there are clefts in both the upper face and mandible, accompanied by fusion of the upper and lower jaws in the hinge region. Comparison of ATAC-seq and RNA-seq analyses between controls and mutants revealed significant changes in chromatin accessibility and gene expression centered on multiple AP-2 binding motifs associated with enhancer elements within these ectodermal lineages. In particular, loss of these AP-2 proteins affects both skin differentiation as well as multiple signaling pathways, most notably the WNT pathway. We also determined that the mutant clefting phenotypes that correlated with reduced WNT signaling could be rescued by Wnt1 ligand overexpression in the ectoderm. Collectively, these findings highlight a conserved ancestral function for AP-2 transcription factors in ectodermal development and signaling, and provide a framework from which to understand the gene regulatory network operating within this tissue that directs vertebrate craniofacial development.

Escape of hair follicle stem cells causes stem cell exhaustion during aging

Stem cell (SC) exhaustion is a hallmark of aging. However, the process of SC depletion during aging has not been observed in live animals, and the underlying mechanism contributing to tissue deterioration remains obscure. We find that, in aged mice, epithelial cells escape from the hair follicle (HF) SC compartment to the dermis, contributing to HF miniaturization. Single-cell RNA-seq and assay for transposase-accessible chromatin using sequencing (ATAC-seq) reveal reduced expression of cell adhesion and extracellular matrix genes in aged HF-SCs, many of which are regulated by Foxc1 and Nfatc1. Deletion of Foxc1 and Nfatc1 recapitulates HF miniaturization and causes hair loss. Live imaging captures individual epithelial cells migrating away from the SC compartment and HF disintegration. This study illuminates a hitherto unknown activity of epithelial cells escaping from their niche as a mechanism underlying SC reduction and tissue degeneration. Identification of homeless epithelial cells in aged tissues provides a new perspective for understanding aging-associated diseases.

Epigenetic regulation and signalling pathways in Merkel cell development

Merkel cells are specialized epithelial cells connected to afferent nerve endings responsible for light-touch sensations, formed at specific locations in touch-sensitive regions of the mammalian skin. Although Merkel cells are descendants of the epidermal lineage, little is known about the mechanisms responsible for the development of these unique mechanosensory cells. Recent studies have highlighted that the Polycomb group (PcG) of proteins play a significant role in spatiotemporal regulation of Merkel cell formation. In addition, several of the major signalling pathways involved in skin development have been shown to regulate Merkel cell development as well. Here, we summarize the current understandings of the role of developmental regulators in Merkel cell formation, including the interplay between the epigenetic machinery and key signalling pathways, and the lineage-specific transcription factors involved in the regulation of Merkel cell development.

Detection of differentially abundant cell subpopulations in scRNA-seq data

Comprehensive and accurate comparisons of transcriptomic distributions of cells from samples taken from two different biological states, such as healthy versus diseased individuals, are an emerging challenge in single-cell RNA sequencing (scRNA-seq) analysis. Current methods for detecting differentially abundant (DA) subpopulations between samples rely heavily on initial clustering of all cells in both samples. Often, this clustering step is inadequate since the DA subpopulations may not align with a clear cluster structure, and important differences between the two biological states can be missed. Here, we introduce DA-seq, a targeted approach for identifying DA subpopulations not restricted to clusters. DA-seq is a multiscale method that quantifies a local DA measure for each cell, which is computed from its k nearest neighboring cells across a range of k values. Based on this measure, DA-seq delineates contiguous significant DA subpopulations in the transcriptomic space. We apply DA-seq to several scRNA-seq datasets and highlight its improved ability to detect differences between distinct phenotypes in severe versus mildly ill COVID-19 patients, melanomas subjected to immune checkpoint therapy comparing responders to nonresponders, embryonic development at two time points, and young versus aging brain tissue. DA-seq enabled us to detect differences between these phenotypes. Importantly, we find that DA-seq not only recovers the DA cell types as discovered in the original studies but also reveals additional DA subpopulations that were not described before. Analysis of these subpopulations yields biological insights that would otherwise be undetected using conventional computational approaches.

ΔNp63 is a pioneer factor that binds inaccessible chromatin and elicits chromatin remodeling

BACKGROUND

ΔNp63 is a master transcriptional regulator playing critical roles in epidermal development and other cellular processes. Recent studies suggest that ΔNp63 functions as a pioneer factor that can target its binding sites within inaccessible chromatin and induce chromatin remodeling.

METHODS

In order to examine if ΔNp63 can bind to inaccessible chromatin and to determine if specific histone modifications are required for binding, we induced ΔNp63 expression in two p63-naïve cell lines. ΔNp63 binding was then examined by ChIP-seq and the chromatin at ΔNp63 targets sites was examined before and after binding. Further analysis with competitive nucleosome binding assays was used to determine how ΔNp63 directly interacts with nucleosomes.

RESULTS

Our results show that before ΔNp63 binding, targeted sites lack histone modifications, indicating ΔNp63's capability to bind at unmodified chromatin. Moreover, the majority of the sites that are bound by ectopic ΔNp63 expression exist in an inaccessible state. Once bound, ΔNp63 induces acetylation of the histone and the repositioning of nucleosomes at its binding sites. Further analysis with competitive nucleosome binding assays reveal that ΔNp63 can bind directly to nucleosome edges with significant binding inhibition occurring within 50 bp of the nucleosome dyad.

CONCLUSION

Overall, our results demonstrate that ΔNp63 is a pioneer factor that binds nucleosome edges at inaccessible and unmodified chromatin sites and induces histone acetylation and nucleosome repositioning.

Defining Epidermal Basal Cell States during Skin Homeostasis and Wound Healing Using Single-Cell Transcriptomics

Our knowledge of transcriptional heterogeneities in epithelial stem and progenitor cell compartments is limited. Epidermal basal cells sustain cutaneous tissue maintenance and drive wound healing. Previous studies have probed basal cell heterogeneity in stem and progenitor potential, but a comprehensive dissection of basal cell dynamics during differentiation is lacking. Using single-cell RNA sequencing coupled with RNAScope and fluorescence lifetime imaging, we identify three non-proliferative and one proliferative basal cell state in homeostatic skin that differ in metabolic preference and become spatially partitioned during wound re-epithelialization. Pseudotemporal trajectory and RNA velocity analyses predict a quasi-linear differentiation hierarchy where basal cells progress from Col17a1^Hi/Trp63^Hi state to early-response state, proliferate at the juncture of these two states, or become growth arrested before differentiating into spinous cells. Wound healing induces plasticity manifested by dynamic basal-spinous interconversions at multiple basal transcriptional states. Our study provides a systematic view of epidermal cellular dynamics, supporting a revised “hierarchical-lineage” model of homeostasis.

Haensel et al. performed a comprehensive dissection of the cellular makeup of skin during homeostasis and wound healing and the molecular heterogeneity and cellular dynamics within its stem-cell-containing epidermal basal layer. Their work provides insights and stimulates further investigation into the mechanism of skin maintenance and repair.

Chromatin accessibility profiling methods

Chromatin accessibility, or the physical access to chromatinized DNA, is a widely studied characteristic of the eukaryotic genome. As active regulatory DNA elements are generally 'accessible', the genome-wide profiling of chromatin accessibility can be used to identify candidate regulatory genomic regions in a tissue or cell type. Multiple biochemical methods have been developed to profile chromatin accessibility, both in bulk and at the single-cell level. Depending on the method, enzymatic cleavage, transposition or DNA methyltransferases are used, followed by high-throughput sequencing, providing a view of genome-wide chromatin accessibility. In this Primer, we discuss these biochemical methods, as well as bioinformatics tools for analysing and interpreting the generated data, and insights into the key regulators underlying developmental, evolutionary and disease processes. We outline standards for data quality, reproducibility and deposition used by the genomics community. Although chromatin accessibility profiling is invaluable to study gene regulation, alone it provides only a partial view of this complex process. Orthogonal assays facilitate the interpretation of accessible regions with respect to enhancer-promoter proximity, functional transcription factor binding and regulatory function. We envision that technological improvements including single-molecule, multi-omics and spatial methods will bring further insight into the secrets of genome regulation.

Epidermal structure and differentiation

In this Primer, Moreci and Lechler follow the lifetime of an epidermal cell from its birth to its ultimate death, and detail how this journey is necessary for epidermal function.

Stem cell-intrinsic mechanisms regulating adult hair follicle homeostasis

Adult hair follicle stem cells (HFSCs) undergo dynamic and periodic molecular changes in their cellular states throughout the hair homeostatic cycle. These states are tightly regulated by cell-intrinsic mechanisms and by extrinsic signals from the microenvironment. HFSCs are essential not only for fuelling hair growth, but also for skin wound healing. Increasing evidence suggests an important role of HFSCs in organizing multiple skin components around the hair follicle, thus functioning as an organizing centre during adult skin homeostasis. Here, we focus on recent findings on cell-intrinsic mechanisms of HFSC homeostasis, which include transcription factors, histone modifications, DNA regulatory elements, non-coding RNAs, cell metabolism, cell polarity and post-transcriptional mRNA processing. Several transcription factors are now known to participate in well-known signalling pathways that control hair follicle homeostasis, as well as in super-enhancer activities to modulate HFSC and progenitor lineage progression. Interestingly, HFSCs have been shown to secrete molecules that are important in guiding the organization of several skin components around the hair follicle, including nerves, arrector pili muscle and vasculature. Finally, we discuss recent technological advances in the field such as single-cell RNA sequencing and live imaging, which revealed HFSC and progenitor heterogeneity and brought new light to understanding crosstalking between HFSCs and the microenvironment. The field is well on its way to generate a comprehensive map of molecular interactions that should serve as a solid theoretical platform for application in hair and skin disease and ageing.

Transcriptional and signalling regulation of skin epithelial stem cells in homeostasis, wounds and cancer

The epidermis and skin appendages are maintained by their resident epithelial stem cells, which undergo long-term self-renewal and multilineage differentiation. Upon injury, stem cells are activated to mediate re-epithelialization and restore tissue function. During this process, they often mount lineage plasticity and expand their fates in response to damage signals. Stem cell function is tightly controlled by transcription machineries and signalling transductions, many of which derail in degenerative, inflammatory and malignant dermatologic diseases. Here, by describing both well-characterized and newly emerged pathways, we discuss the transcriptional and signalling mechanisms governing skin epithelial homeostasis, wound repair and squamous cancer. Throughout, we highlight common themes underscoring epithelial stem cell plasticity and tissue-level crosstalk in the context of skin physiology and pathology.

Meeting Report: 68th Montagna Symposium on the Biology of Skin "Decoding Complex Skin Diseases: Integrating Genetics, Genomics, and Disease Biology"

The 68^th Montagna Symposium on the Biology of the Skin was held from 10 to 14 October 2019 at Salishan Lodge in Gleneden Beach, Oregon. The theme of the meeting was "Decoding Complex Skin Diseases: Integrating Genetics, Genomics, and Disease Biology." The meeting emphasized the integration of multiple themes and disciplines to better understand some of the most common skin diseases, ranging from psoriasis to alopecia areata to vitiligo to lupus erythematosus to atopic dermatitis and food allergy. Promising therapeutic strategies are emerging for all of these diseases, providing clues for ways to connect the bench to the bedside. A common thread was the success of GWASs, which have highlighted the importance of regulatory signals versus coding variation. These diseases also share an environmental component linked to immune system function. Hence, beyond GWASs, this meeting focused on gene regulatory mechanisms, the single-cell revolution, in vivo systems for dissection of disease pathogenesis, and the relationship between genetics and environment in the context of host defense. We concluded with a translational roundtable designed to explore how these interrelated fields can best be directed toward long-term disease control and, ultimately, a cure.

TP63 links chromatin remodeling and enhancer reprogramming to epidermal differentiation and squamous cell carcinoma development

Squamous cell carcinoma (SCC) is an aggressive malignancy that can originate from various organs. TP63 is a master regulator that plays an essential role in epidermal differentiation. It is also a lineage-dependent oncogene in SCC. ΔNp63α is the prominent isoform of TP63 expressed in epidermal cells and SCC, and overexpression promotes SCC development through a variety of mechanisms. Recently, ΔNp63α was highlighted to act as an epidermal-specific pioneer factor that binds closed chromatin and enhances chromatin accessibility at epidermal enhancers. ΔNp63α coordinates chromatin-remodeling enzymes to orchestrate the tissue-specific enhancer landscape and three-dimensional high-order architecture of chromatin. Moreover, ΔNp63α establishes squamous-like enhancer landscapes to drive oncogenic target expression during SCC development. Importantly, ΔNp63α acts as an upstream regulator of super enhancers to activate a number of oncogenic transcripts linked to poor prognosis in SCC. Mechanistically, ΔNp63α activates genes transcription through physically interacting with a number of epigenetic modulators to establish enhancers and enhance chromatin accessibility. In contrast, ΔNp63α also represses gene transcription via interacting with repressive epigenetic regulators. ΔNp63α expression is regulated at multiple levels, including transcriptional, post-transcriptional, and post-translational levels. In this review, we summarize recent advances of p63 in epigenomic and transcriptional control, as well as the mechanistic regulation of p63.

Murine interfollicular epidermal differentiation is gradualistic with GRHL3 controlling progression from stem to transition cell states

The interfollicular epidermis (IFE) forms a water-tight barrier that is often disrupted in inflammatory skin diseases. During homeostasis, the IFE is replenished by stem cells in the basal layer that differentiate as they migrate toward the skin surface. Conventionally, IFE differentiation is thought to be stepwise as reflected in sharp boundaries between its basal, spinous, granular and cornified layers. The transcription factor GRHL3 regulates IFE differentiation by transcriptionally activating terminal differentiation genes. Here we use single cell RNA-seq to show that murine IFE differentiation is best described as a single step gradualistic process with a large number of transition cells between the basal and spinous layer. RNA-velocity analysis identifies a commitment point that separates the plastic basal and transition cell state from unidirectionally differentiating cells. We also show that in addition to promoting IFE terminal differentiation, GRHL3 is essential for suppressing epidermal stem cell expansion and the emergence of an abnormal stem cell state by suppressing Wnt signaling in stem cells.

p63 and Its Target Follistatin Maintain Salivary Gland Stem/Progenitor Cell Function through TGF-β/Activin Signaling

Multipotent ΔNp63-positive cells maintain all epithelial cell lineages of the embryonic and adult salivary gland (SG). However, the molecular mechanisms by which ΔNp63 regulates stem/progenitor (SP) cell populations in the SG remains elusive. To understand the role of ΔNp63 in directing cell fate choices in this gland, we have generated ΔNp63-deleted adult mice and primary salivary cell cultures to probe alterations in SP cell differentiation and function. In parallel, we have leveraged RNA-seq and ChIP-seq-based characterization of the ΔNp63-driven cistrome and scRNA-seq analysis to molecularly interrogate altered SG cellular identities and differentiation states dependent on ΔNp63. Our studies reveal that ablation of ΔNp63 results in a loss of the SP cell population and skewed differentiation that is mediated by Follistatin-dependent dysregulated TGF-β/Activin signaling. These findings offer new revelations into the SP cell gene regulatory networks that are likely to be relevant for normal or diseased SG states.

An updated classification of hair follicle morphogenesis

Hair follicle (HF) formation in developing embryonic skin requires stepwise signalling between the epithelial epidermis and mesenchymal dermis, and their specialized derivatives, the placode/germ/peg and dermal condensate/papilla, respectively. Classically, distinct stages of HF morphogenesis have been defined, in the mouse model, based on (a) changes in cell morphology and aggregation; (b) expression of few known molecular markers; (c) the extent of follicle downgrowth; and (d) the presence of differentiating cell types. Refined genetic strategies and recent emerging technologies, such as live imaging and transcriptome analyses of isolated cell populations or single cells, have enabled a closer dissection of the signalling requirements at different stages of HF formation, particularly early on. They have also led to the discovery of precursor cells for placode, dermal condensate and future bulge stem cells that, combined with molecular insights into their fate specification and subsequent formation, serve as novel landmarks for early HF morphogenetic events and studies of the signalling networks mediating these processes. In this review, we integrate the emergence of HF precursor cell states and novel molecular markers of fate and formation to update the widely used 20-year-old seminal classification guide of HF morphogenetic stages by Paus et al. We then temporally describe the latest insights into the early cellular and molecular events and signalling requirements for HF morphogenesis in relation to one another in a holistic manner.

The KEAP1/NRF2 Signaling Pathway in Keratinization

Keratinization is a tissue adaptation, but aberrant keratinization is associated with skin disorders such as ichthyoses, atopic dermatitis, psoriasis, and acne. The disease phenotype stems from the interaction between genes and the environment; therefore, an understanding of the adaptation machinery may lead to a new appreciation of pathomechanisms. The KEAP1/NRF2 signaling pathway mediates the environmental responses of squamous epithelial tissue. The unpredicted outcome of the Keap1-null mutation in mice allowed us to revisit the basic principle of the biological process of keratinization: sulfur metabolism establishes unparalleled cytoprotection in the body wall of terrestrial mammals. We summarize the recent understanding of the KEAP1/NRF2 signaling pathway, which is a thiol-based sensor-effector apparatus, with particular focuses on epidermal differentiation in the context of the gene-environment interaction, the structure/function principles involved in KEAP1/NRF2 signaling, lessons from mouse models, and their pathological implications. This synthesis may provide insights into keratinization, which provides physical insulation and constitutes an essential innate integumentary defense system.

Single-cell transcriptional profiling identifies a cluster of potential metastasis-associated UBE2C+ cells in immature ovarian teratoma

To dissect the disease heterogeneity and identify the underlying cellular and molecular events related to metastasis of immature ovarian teratoma in children, single-cell RNA sequencing was performed for a 2-year-old patient with liver metastases from immature ovarian teratoma. A total of 5976 cells were obtained for further analysis, with a median unique molecular identifier count of 6011 per cell and a median number of 1741 genes detected per cell. Fourteen clusters were recognized, with the main lineages comprising epithelial cells, macrophages, fibroblasts, glial cells, and dendritic cells. Ten subclusters of epithelial cells were further defined, originating from the urinary tract, esophagus, bronchus, lung, skin, and gastrointestinal tract. An undefined UBE2C + population in an active state of proliferation was also identified and its biological processes were related to meiosis and maturation of oocytes. Pseudotime analysis revealed different distributions of epithelial cells in the development trajectory. In conclusion, a cluster of UBE2C + epithelial cells in an active state of proliferation was identified in an immature ovarian teratoma in a child, and may contribute to metastasis by regulating epithelial-mesenchymal transition. These findings help toward understanding the origin of the malignant behaviors, offer a potential biomarker for early determination of the tumor nature, and provide new ideas for the therapy of immature ovarian teratoma in children.

Applications of single-cell sequencing for the field of otolaryngology: A contemporary review

OBJECTIVES

Single-cell RNA sequencing (scRNA-Seq) is a new technique used to interrogate the transcriptome of individual cells within native tissues that have already resulted in key discoveries in auditory basic science research. Rapid advances in scRNA-Seq make it likely that it will soon be translated into clinical medicine. The goal of this review is to inspire the use of scRNA-Seq in otolaryngology by giving examples of how it can be applied to patient samples and how this information can be used clinically.

METHODS

Studies were selected based on the scientific quality and relevance to scRNA-Seq. In addition to mouse auditory system (inner ear including hair cells and supporting cells, spiral ganglion neurons, and inner ear organoids), recent studies using human primary cell samples are discussed. We also perform our own analysis on publicly available, published scRNA-Seq data from oral head and neck squamous cell carcinoma (HNSCC) samples to serve as an example of a clinically relevant application of scRNA-Seq.

RESULTS

Studies focusing on patient tissues show that scRNA-Seq reveals tissue heterogeneity and rare-cell types responsible for disease pathogenesis. The heterogeneity detected by scRNA-Seq can result in both the identification of known or novel disease biomarkers and drug targets. Our analysis of HNSCC data gives an example for how otolaryngologists can use scRNA-Seq for clinical use.

CONCLUSIONS

Although there are limitations to the translation of scRNA-Seq to the clinic, we show that its use in otolaryngology can give physicians insight into the tissue heterogeneity within their patient's diseased tissue giving them information on disease pathogenesis, novel disease biomarkers or druggable targets, and aid in selecting patient-specific drug cocktails.

Clinical implications of chromatin accessibility in human cancers

Assay for transposase-accessible chromatin using sequencing (ATAC-seq) has not yet been widely used in cancer research. Clinical implications of chromatin accessibility assessed by ATAC-seq profiling in human cancers especially in a large patient cohort is largely unknown. In this study, we analyzed ATAC-seq data in 404 cancer patients from the Cancer Genome Atlas, representing the largest cancer patient cohort with ATAC-seq data, and correlated chromatin accessibility with patient demographics, tumor histology, molecular subtypes, and survival. Our results showed that chromatin accessibility varies from chromosome to chromosome, and is different in different genomic regions along the same chromosome. Chromatin accessibility especially on the X chromosome is strongly dependent on patient sex, but not much on patient age or tumor stage. Striking difference in chromatin accessibility is observed between lung adenocarcinoma and lung squamous cell carcinoma, the two most common histological subgroups in lung cancer. Furthermore, chromatin accessibility was different between basal and non-basal breast cancer. Finally, we identified prognostic peaks in the promoter regions that were significantly correlated with survival. In particular, we identified six peaks in the ESR1 gene promoter region in the ATAC-seq profiling and found that the peak about 247 bp away from the transcription start site was significantly associated with better survival. In conclusion, our study provides an alternative mechanism underlying tumor prognosis.

Unraveling cancer lineage drivers in squamous cell carcinomas

Cancer hijacks embryonic development and adult wound repair mechanisms to fuel malignancy. Cancer frequently originates from de-regulated adult stem cells or progenitors, which are otherwise essential units for postnatal tissue remodeling and repair. Cancer genomics studies have revealed convergence of multiple cancers across organ sites, including squamous cell carcinomas (SCCs), a common group of cancers arising from the head and neck, esophagus, lung, cervix and skin. In this review, we summarize our current knowledge on the molecular drivers of SCCs, including these five major organ sites. We especially focus our discussion on lineage dependent driver genes and pathways, in the context of squamous development and stratification. We then use skin as a model to discuss the notion of field cancerization during SCC carcinogenesis, and cancer as a wound that never heals. Finally, we turn to the idea of context dependency widely observed in cancer driver genes, and outline literature support and possible explanations for their lineage specific functions. Through these discussions, we aim to provide an up-to-date summary of molecular mechanisms driving tumor plasticity in squamous cancers. Such basic knowledge will be helpful to inform the clinics for better stratifying cancer patients, revealing novel drug targets and providing effective treatment options.

Single-cell RNA-seq identifies a reversible mesodermal activation in abnormally specified epithelia of p63 EEC syndrome

Mutations in transcription factor p63 are associated with developmental disorders that manifest defects in stratified epithelia including the epidermis. The underlying cellular and molecular mechanism is however not yet understood. We established an epidermal commitment model using human induced pluripotent stem cells (iPSCs) and characterized differentiation defects of iPSCs derived from ectrodactyly, ectodermal dysplasia, and cleft lip/palate (EEC) syndrome patients carrying p63 mutations. Transcriptome analyses revealed stepwise cell fate transitions during epidermal commitment: Specification from multipotent simple epithelium to basal stratified epithelia and ultimately to the mature epidermal fate. Differentiation defects of EEC iPSCs caused by p63 mutations occurred during the specification switch from the simple epithelium to the basal-stratified epithelial fate. Single-cell transcriptome and pseudotime analyses of cell states identified mesodermal activation that was associated with the deviated commitment route of EEC iPSCs. Integrated analyses of differentially regulated genes and p63-dependent dynamic genomic enhancers during epidermal commitment suggest that p63 directly controls epidermal gene activation at the specification switch and has an indirect effect on mesodermal gene repression. Importantly, inhibitors of mesodermal induction enhanced epidermal commitment of EEC iPSCs. Our findings demonstrate that p63 is required for specification of stratified epithelia, and that epidermal commitment defects caused by p63 mutations can be reversed by repressing mesodermal induction. This study provides insights into disease mechanisms underlying stratified epithelial defects caused by p63 mutations and suggests potential therapeutic strategies for the disease.

Pioneer and repressive functions of p63 during zebrafish embryonic ectoderm specification

The transcription factor p63 is a master regulator of ectoderm development. Although previous studies show that p63 triggers epidermal differentiation in vitro, the roles of p63 in developing embryos remain poorly understood. Here, we use zebrafish embryos to analyze in vivo how p63 regulates gene expression during development. We generate tp63-knock-out mutants that recapitulate human phenotypes and show down-regulated epidermal gene expression. Following p63-binding dynamics, we find two distinct functions clearly separated in space and time. During early development, p63 binds enhancers associated to neural genes, limiting Sox3 binding and reducing neural gene expression. Indeed, we show that p63 and Sox3 are co-expressed in the neural plate border. On the other hand, p63 acts as a pioneer factor by binding non-accessible chromatin at epidermal enhancers, promoting their opening and epidermal gene expression in later developmental stages. Therefore, our results suggest that p63 regulates cell fate decisions during vertebrate ectoderm specification.

p63 cooperates with CTCF to modulate chromatin architecture in skin keratinocytes

The transcription factor p63 regulates epidermal genes and the enhancer landscape in skin keratinocytes. Its molecular function in controlling the chromatin structure is, however, not yet completely understood. Here, we integrated multi-omics profiles, including the transcriptome, transcription factor DNA-binding and chromatin accessibility, in skin keratinocytes isolated from EEC syndrome patients carrying p63 mutations, to examine the role of p63 in shaping the chromatin architecture. We found decreased chromatin accessibility in p63- and CTCF-bound open chromatin regions that potentially contributed to gene deregulation in mutant keratinocytes. Cooperation of p63 and CTCF seemed to assist chromatin interactions between p63-bound enhancers and gene promoters in skin keratinocytes. Our study suggests an intriguing model where cell type-specific transcription factors such as p63 cooperate with the genome organizer CTCF in the three-dimensional chromatin space to regulate the transcription program important for the proper cell identity.

Epigenetic control in skin development, homeostasis and injury repair

Cell-type- and cell-state-specific patterns of covalent modifications on DNA and histone tails form global epigenetic profiles that enable spatiotemporal regulation of gene expression. These epigenetic profiles arise from coordinated activities of transcription factors and epigenetic modifiers, which result in cell-type-specific outputs in response to dynamic environmental conditions and signalling pathways. Recent mouse genetic and functional studies have highlighted the physiological significance of global DNA and histone epigenetic modifications in skin. Importantly, specific epigenetic profiles are emerging for adult skin stem cells that are associated with their cell fate plasticity and proper activity in tissue regeneration. We can now begin to draw a more comprehensive picture of how epigenetic modifiers orchestrate their cell-intrinsic role with microenvironmental cues for proper skin development, homeostasis and wound repair. The field is ripe to begin to implement these findings from the laboratory into skin therapies.

De Novo PITX1 Expression Controls Bi-Stable Transcriptional Circuits to Govern Self-Renewal and Differentiation in Squamous Cell Carcinoma

Basal tumor propagating cells (TPCs) control squamous cell carcinoma (SCC) growth by self-renewing and differentiating into supra-basal SCC cells, which lack proliferative potential. While transcription factors such as SOX2 and KLF4 can drive these behaviors, their molecular roles and regulatory interactions with each other have remained elusive. Here, we show that PITX1 is specifically expressed in TPCs, where it co-localizes with SOX2 and TRP63 and determines cell fate in mouse and human SCC. Combining gene targeting with chromatin immunoprecipitation sequencing (ChIP-seq) and transcriptomic analyses reveals that PITX1 cooperates with SOX2 and TRP63 to sustain an SCC-specific transcriptional feed-forward circuit that maintains TPC-renewal, while inhibiting KLF4 expression and preventing KLF4-dependent differentiation. Conversely, KLF4 represses PITX1, SOX2, and TRP63 expression to prevent TPC expansion. This bi-stable, multi-input network reveals a molecular framework that explains self-renewal, aberrant differentiation, and SCC growth in mice and humans, providing clues for developing differentiation-inducing therapeutic strategies.

Genetic and Epigenetic Alterations Underlie Oligodendroglia Susceptibility and White Matter Etiology in Psychiatric Disorders

Numerous genetic risk loci are found to associate with major neuropsychiatric disorders represented by schizophrenia. The pathogenic roles of genetic risk loci in psychiatric diseases are further complicated by the association with cell lineage- and/or developmental stage-specific epigenetic alterations. Besides aberrant assembly and malfunction of neuronal circuitry, an increasing volume of discoveries clearly demonstrate impairment of oligodendroglia and disruption of white matter integrity in psychiatric diseases. Nonetheless, whether and how genetic risk factors and epigenetic dysregulations for neuronal susceptibility may affect oligodendroglia is largely unknown. In this mini-review, we will discuss emerging evidence regarding the functional interplay between genetic risk loci and epigenetic factors, which may underlie compromised oligodendroglia and myelin development in neuropsychiatric disorders. Transcriptional and epigenetic factors are the major aspects affected in oligodendroglia. Moreover, multiple disease susceptibility genes are connected by epigenetically modulated transcriptional and post-transcriptional mechanisms. Oligodendroglia specific complex molecular orchestra may explain how distinct risk factors lead to the common clinical expression of white matter pathology of neuropsychiatric disorders.