An AhR-Ovol1-Id1 regulatory axis in keratinocytes promotes skin homeostasis against atopic dermatitis.

Skin is our outer permeability and immune defense barrier against myriad external assaults. Aryl hydrocarbon receptor (AhR) senses environmental factors and regulates barrier robustness and immune homeostasis. AhR agonist is in clinical trial for atopic dermatitis (AD) treatment, but the underlying mechanism of action remains ill-defined. Here we report OVOL1/Ovol1 as a conserved and direct transcriptional target of AhR in epidermal keratinocytes. We show that OVOL1/Ovol1 impacts AhR regulation of keratinocyte gene expression, and Ovol1 deletion in keratinocytes hampers AhR's barrier promotion function and worsens AD-like inflammation. Mechanistically, we identify Ovol1's direct downstream targets genome-wide, and provide in vivo evidence for Id1's critical role in barrier maintenance and disease suppression. Furthermore, our findings reveal an IL-1/dermal γδT cell axis exacerbating both type 2 and type 3 immune responses downstream of barrier perturbation in Ovol1 -deficient AD skin. Finally, we present data suggesting the clinical relevance of OVOL1 and ID1 function in human AD. Our study highlights a keratinocyte-intrinsic AhR-Ovol1-Id1 regulatory axis that promotes both epidermal and immune homeostasis against AD-like inflammation, implicating new therapeutic targets for AD.

Acquisition of Drug Resistance in Basal Cell Nevus Syndrome Tumors through Basal to Squamous Cell Carcinoma Transition.

Although basal cell carcinomas arise from ectopic Hedgehog pathway activation and can be treated with pathway inhibitors, sporadic basal cell carcinomas display high resistance rates, whereas tumors arising in patients with Gorlin syndrome with germline Patched (PTCH1) alterations are uniformly suppressed by inhibitor therapy. In rare cases, patients with Gorlin syndrome on long-term inhibitor therapy will develop individual resistant tumor clones that rapidly progress, but the basis of this resistance remains unstudied. In this study, we report a case of an SMO inhibitor-resistant tumor arising in a patient with Gorlin syndrome on suppressive SMO inhibitor for nearly a decade. Using a combination of multiomics and spatial transcriptomics, we define the tumor populations at the cellular and tissue level to conclude that Gorlin tumors can develop resistance to SMO inhibitors through the previously described basal to squamous cell carcinoma transition. Intriguingly, through spatial whole-exome genomic analysis, we nominate PCYT2, ETNK1, and the phosphatidylethanolamine biosynthetic pathway as genetic suppressors of basal to squamous cell carcinoma transition resistance. These observations provide a general framework for studying tumor evolution and provide important clinical insight into mechanisms of resistance to SMO inhibitors for not only Gorlin syndrome but also sporadic basal cell carcinomas.

Acquisition of drug resistance in basal cell nevus syndrome tumors through basal to squamous cell carcinoma transition.

While basal cell carcinomas (BCCs) arise from ectopic hedgehog pathway activation and can be treated with pathway inhibitors, sporadic BCCs display high resistance rates while tumors arising in Gorlin syndrome patients with germline Patched ( PTCH1 ) mutations are uniformly suppressed by inhibitor therapy. In rare cases, Gorlin syndrome patients on long-term inhibitor therapy will develop individual resistant tumor clones that rapidly progress, but the basis of this resistance remains unstudied. Here we report a case of an SMO i -resistant tumor arising in a Gorlin patient on suppressive SMO i for nearly a decade. Using a combination of multi-omics and spatial transcriptomics, we define the tumor populations at the cellular and tissue level to conclude that Gorlin tumors can develop resistance to SMO i through the previously described basal to squamous cell carcinoma transition (BST). Intriguingly, through spatial whole exome genomic analysis, we nominate PCYT2, ETNK1, and the phosphatidylethanolamine biosynthetic pathway as novel genetic suppressors of BST resistance. These observations provide a general framework for studying tumor evolution and provide important clinical insight into mechanisms of resistance to SMO i for not only Gorlin syndrome but sporadic BCCs as well.

Ovol1/2 loss-induced epidermal defects elicit skin immune activation and alter global metabolism.

Skin epidermis constitutes the outer permeability barrier that protects the body from dehydration, heat loss, and myriad external assaults. Mechanisms that maintain barrier integrity in constantly challenged adult skin and how epidermal dysregulation shapes the local immune microenvironment and whole-body metabolism remain poorly understood. Here, we demonstrate that inducible and simultaneous ablation of transcription factor-encoding Ovol1 and Ovol2 in adult epidermis results in barrier dysregulation through impacting epithelial-mesenchymal plasticity and inflammatory gene expression. We find that aberrant skin immune activation then ensues, featuring Langerhans cell mobilization and T cell responses, and leading to elevated levels of secreted inflammatory factors in circulation. Finally, we identify failure to gain body weight and accumulate body fat as long-term consequences of epidermal-specific Ovol1/2 loss and show that these global metabolic changes along with the skin barrier/immune defects are partially rescued by immunosuppressant dexamethasone. Collectively, our study reveals key regulators of adult barrier maintenance and suggests a causal connection between epidermal dysregulation and whole-body metabolism that is in part mediated through aberrant immune activation.

Skin basal cell carcinomas assemble a pro-tumorigenic spatially organized and self-propagating Trem2+ myeloid niche.

Cancer immunotherapies have revolutionized treatment but have shown limited success as single-agent therapies highlighting the need to understand the origin, assembly, and dynamics of heterogeneous tumor immune niches. Here, we use single-cell and imaging-based spatial analysis to elucidate three microenvironmental neighborhoods surrounding the heterogeneous basal cell carcinoma tumor epithelia. Within the highly proliferative neighborhood, we find that TREM2+ skin cancer-associated macrophages (SCAMs) support the proliferation of a distinct tumor epithelial population through an immunosuppression-independent manner via oncostatin-M/JAK-STAT3 signaling. SCAMs represent a unique tumor-specific TREM2+ population defined by VCAM1 surface expression that is not found in normal homeostatic skin or during wound healing. Furthermore, SCAMs actively proliferate and self-propagate through multiple serial tumor passages, indicating long-term potential. The tumor rapidly drives SCAM differentiation, with intratumoral injections sufficient to instruct naive bone marrow-derived monocytes to polarize within days. This work provides mechanistic insights into direct tumor-immune niche dynamics independent of immunosuppression, providing the basis for potential combination tumor therapies.

LY6D marks pre-existing resistant basosquamous tumor subpopulations.

Improved response to canonical therapies requires a mechanistic understanding of dynamic tumor heterogeneity by identifying discrete cellular populations with enhanced cellular plasticity. We have previously demonstrated distinct resistance mechanisms in skin basal cell carcinomas, but a comprehensive understanding of the cellular states and markers associated with these populations remains poorly understood. Here we identify a pre-existing resistant cellular population in naive basal cell carcinoma tumors marked by the surface marker LY6D. LY6D+ tumor cells are spatially localized and possess basal cell carcinoma and squamous cell carcinoma-like features. Using computational tools, organoids, and spatial tools, we show that LY6D+ basosquamous cells represent a persister population lying on a central node along the skin lineage-associated spectrum of epithelial states with local environmental and applied therapies determining the kinetics of accumulation. Surprisingly, LY6D+ basosquamous populations exist in many epithelial tumors, such as pancreatic adenocarcinomas, which have poor outcomes. Overall, our results identify the resistant LY6D+ basosquamous population as an important clinical target and suggest strategies for future therapeutic approaches to target them.

Wound healing in aged skin exhibits systems-level alterations in cellular composition and cell-cell communication.

Delayed and often impaired wound healing in the elderly presents major medical and socioeconomic challenges. A comprehensive understanding of the cellular/molecular changes that shape complex cell-cell communications in aged skin wounds is lacking. Here, we use single-cell RNA sequencing to define the epithelial, fibroblast, immune cell types, and encompassing heterogeneities in young and aged skin during homeostasis and identify major changes in cell compositions, kinetics, and molecular profiles during wound healing. Our comparative study uncovers a more pronounced inflammatory phenotype in aged skin wounds, featuring neutrophil persistence and higher abundance of an inflammatory/glycolytic Arg1Hi macrophage subset that is more likely to signal to fibroblasts via interleukin (IL)-1 than in young counterparts. We predict systems-level differences in the number, strength, route, and signaling mediators of putative cell-cell communications in young and aged skin wounds. Our study exposes numerous cellular/molecular targets for functional interrogation and provides a hypothesis-generating resource for future wound healing studies.

c-FOS drives reversible basal to squamous cell carcinoma transition.

While squamous transdifferentiation within subpopulations of adenocarcinomas represents an important drug resistance problem, its underlying mechanism remains poorly understood. Here, using surface markers of resistant basal cell carcinomas (BCCs) and patient single-cell and bulk transcriptomic data, we uncover the dynamic roadmap of basal to squamous cell carcinoma transition (BST). Experimentally induced BST identifies activator protein 1 (AP-1) family members in regulating tumor plasticity, and we show that c-FOS plays a central role in BST by regulating the accessibility of distinct AP-1 regulatory elements. Remarkably, despite prominent changes in cell morphology and BST marker expression, we show using inducible model systems that c-FOS-mediated BST demonstrates reversibility. Blocking EGFR pathway activation after c-FOS induction partially reverts BST in vitro and prevents BST features in both mouse models and human tumors. Thus, by identifying the molecular basis of BST, our work reveals a therapeutic opportunity targeting plasticity as a mechanism of tumor resistance.

OVOL1 Regulates Psoriasis-Like Skin Inflammation and Epidermal Hyperplasia.

Psoriasis is a common inflammatory skin disease characterized by aberrant inflammation and epidermal hyperplasia. Molecular mechanisms that regulate psoriasis-like skin inflammation remain to be fully understood. Here, we show that the expression of Ovol1 (encoding ovo-like 1 transcription factor) is upregulated in psoriatic skin, and its deletion results in aggravated psoriasis-like skin symptoms following stimulation with imiquimod. Using bulk and single-cell RNA sequencing, we identify molecular changes in the epidermal, fibroblast, and immune cells of Ovol1-deficient skin that reflect an altered course of epidermal differentiation and enhanced inflammatory responses. Furthermore, we provide evidence for excessive full-length IL-1α signaling in the microenvironment of imiquimod-treated Ovol1-deficient skin that functionally contributes to immune cell infiltration and epidermal hyperplasia. Collectively, our study uncovers a protective role for OVOL1 in curtailing psoriasis-like inflammation and the associated skin pathology.

AP-1 and TGFß cooperativity drives non-canonical Hedgehog signaling in resistant basal cell carcinoma.

Tumor heterogeneity and lack of knowledge about resistant cell states remain a barrier to targeted cancer therapies. Basal cell carcinomas (BCCs) depend on Hedgehog (Hh)/Gli signaling, but can develop mechanisms of Smoothened (SMO) inhibitor resistance. We previously identified a nuclear myocardin-related transcription factor (nMRTF) resistance pathway that amplifies noncanonical Gli1 activity, but characteristics and drivers of the nMRTF cell state remain unknown. Here, we use single cell RNA-sequencing of patient tumors to identify three prognostic surface markers (LYPD3, TACSTD2, and LY6D) which correlate with nMRTF and resistance to SMO inhibitors. The nMRTF cell state resembles transit-amplifying cells of the hair follicle matrix, with AP-1 and TGFß cooperativity driving nMRTF activation. JNK/AP-1 signaling commissions chromatin accessibility and Smad3 DNA binding leading to a transcriptional program of RhoGEFs that facilitate nMRTF activity. Importantly, small molecule AP-1 inhibitors selectively target LYPD3+/TACSTD2+/LY6D+ nMRTF human BCCs ex vivo, opening an avenue for improving combinatorial therapies.

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 Col17a1Hi/Trp63Hi 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.

Ex Vivo Imaging and Genetic Manipulation of Mouse Hair Follicle Bulge Stem Cells.

Stem cells that reside in the bulge of adult mouse hair follicles are a leading model of tissue stem cell research. Ex vivo culturing, molecular and cell biological characterizations, as well as genetic manipulation of fluorescence-activated cell sorting-isolated bulge stem cells offer a useful experimental pipeline to complement in vivo studies. Here we describe detailed methods for culturing, immunostaining, live cell imaging, and adenoviral infection of bulge stem cells for downstream applications such as in vitro clonal and in vivo patch assays.

An Ovol2-Zeb1 transcriptional circuit regulates epithelial directional migration and proliferation.

Directional migration is inherently important for epithelial tissue regeneration and repair, but how it is precisely controlled and coordinated with cell proliferation is unclear. Here, we report that Ovol2, a transcriptional repressor that inhibits epithelial-to-mesenchymal transition (EMT), plays a crucial role in adult skin epithelial regeneration and repair. Ovol2-deficient mice show compromised wound healing characterized by aberrant epidermal cell migration and proliferation, as well as delayed anagen progression characterized by defects in hair follicle matrix cell proliferation and subsequent differentiation. Epidermal keratinocytes and bulge hair follicle stem cells (Bu-HFSCs) lacking Ovol2 fail to expand in culture and display molecular alterations consistent with enhanced EMT and reduced proliferation. Live imaging of wound explants and Bu-HFSCs reveals increased migration speed but reduced directionality, and post-mitotic cell cycle arrest. Remarkably, simultaneous deletion of Zeb1 encoding an EMT-promoting factor restores directional migration to Ovol2-deficient Bu-HFSCs. Taken together, our findings highlight the important function of an Ovol2-Zeb1 EMT-regulatory circuit in controlling the directional migration of epithelial stem and progenitor cells to facilitate adult skin epithelial regeneration and repair.

Intermediate cell states in epithelial-to-mesenchymal transition.

The transition of epithelial cells into a mesenchymal state (epithelial-to-mesenchymal transition or EMT) is a highly dynamic process implicated in various biological processes. During EMT, cells do not necessarily exist in 'pure' epithelial or mesenchymal states. There are cells with mixed (or hybrid) features of the two, which are termed as the intermediate cell states (ICSs). While the exact functions of ICS remain elusive, together with EMT it appears to play important roles in embryogenesis, tissue development, and pathological processes such as cancer metastasis. Recent single cell experiments and advanced mathematical modeling have improved our capability in identifying ICS and provided a better understanding of ICS in development and disease. Here, we review the recent findings related to the ICS in/or EMT and highlight the challenges in the identification and functional characterization of ICS.

Multiscale modeling of layer formation in epidermis.

The mammalian skin epidermis is a stratified epithelium composed of multiple layers of epithelial cells that exist in appropriate sizes and proportions, and with distinct boundaries separating each other. How the epidermis develops from a single layer of committed precursor cells to form a complex multilayered structure of multiple cell types remains elusive. Here, we construct stochastic, three-dimensional, and multiscale models consisting of a lineage of multiple cell types to study the control of epidermal development. Symmetric and asymmetric cell divisions, stochastic cell fate transitions within the lineage, extracellular morphogens, cell-to-cell adhesion forces, and cell signaling are included in model. A GPU algorithm was developed and implemented to accelerate the simulations. These simulations show that a balance between cell proliferation and differentiation during lineage progression is crucial for the development and maintenance of the epidermal tissue. We also find that selective intercellular adhesion is critical to sharpening the boundary between layers and to the formation of a highly ordered structure. The long-range action of a morphogen provides additional feedback regulations, enhancing the robustness of overall layer formation. Our model is built upon previous experimental findings revealing the role of Ovol transcription factors in regulating epidermal development. Direct comparisons of experimental and simulation perturbations show remarkable consistency. Taken together, our results highlight the major determinants of a well-stratified epidermis: balanced proliferation and differentiation, and a combination of both short- (symmetric/asymmetric division and selective cell adhesion) and long-range (morphogen) regulations. These underlying principles have broad implications for other developmental or regenerative processes leading to the formation of multilayered tissue structures, as well as for pathological processes such as epidermal wound healing.

Epithelial-to-mesenchymal transition in cutaneous wound healing: Where we are and where we are heading.

Cutaneous wound healing occurs in distinct yet overlapping steps with the end goal of reforming a stratified epithelium to restore epidermal barrier function. A key component of this process is re-epithelialization, which involves the proliferation and migration of epidermal keratinocytes surrounding the wound. This spatiotemporally controlled process resembles aspects of the epithelial-to-mesenchymal transition (EMT) process and is thus proposed to involve a partial EMT. Here, we review current literature on the cellular and molecular changes that occur during, and the known or potential regulatory factors of cutaneous wound re-epithelialization and EMT to highlight their similarities and differences. We also discuss possible future directions toward a better understanding of the underlying regulatory mechanisms with implications for developing new therapeutics to improve wound repair in humans. Developmental Dynamics 247:473-480, 2018. © 2017 Wiley Periodicals, Inc.