Transcriptional Programming of Normal and Inflamed Human Epidermis at Single-Cell Resolution

Skin in the game: a review of single-cell and spatial transcriptomics in dermatological research

Single-cell RNA sequencing (scRNA-seq) and spatial transcriptomics (ST) are two emerging research technologies that uniquely characterize gene expression microenvironments on a cellular or subcellular level. The skin, a clinically accessible tissue composed of diverse, essential cell populations, serves as an ideal target for these high-resolution investigative approaches. Using these tools, researchers are assembling a compendium of data and discoveries in healthy skin as well as a range of dermatologic pathophysiologies, including atopic dermatitis, psoriasis, and cutaneous malignancies. The ongoing advancement of single-cell approaches, coupled with anticipated decreases in cost with increased adoption, will reshape dermatologic research, profoundly influencing disease characterization, prognosis, and ultimately clinical practice.

The Immunology of Psoriasis-Current Concepts in Pathogenesis

Psoriasis is one of the most common inflammatory skin diseases with a chronic, relapsing-remitting course. The last decades of intense research uncovered a pathological network of interactions between immune cells and other types of cells in the pathogenesis of psoriasis. Emerging evidence indicates that dendritic cells, TH17 cells, and keratinocytes constitute a pathogenic triad in psoriasis. Dendritic cells produce TNF-α and IL-23 to promote T cell differentiation toward TH17 cells that produce key psoriatic cytokines IL-17, IFN-γ, and IL-22. Their activity results in skin inflammation and activation and hyperproliferation of keratinocytes. In addition, other cells and signaling pathways are implicated in the pathogenesis of psoriasis, including TH9 cells, TH22 cells, CD8+ cytotoxic cells, neutrophils, γδ T cells, and cytokines and chemokines secreted by them. New insights from high-throughput analysis of lesional skin identified novel signaling pathways and cell populations involved in the pathogenesis. These studies not only expanded our knowledge about the mechanisms of immune response and the pathogenesis of psoriasis but also resulted in a revolution in the clinical management of patients with psoriasis. Thus, understanding the mechanisms of immune response in psoriatic inflammation is crucial for further studies, the development of novel therapeutic strategies, and the clinical management of psoriasis patients. The aim of the review was to comprehensively present the dysregulation of immune response in psoriasis with an emphasis on recent findings. Here, we described the role of immune cells, including T cells, B cells, dendritic cells, neutrophils, monocytes, mast cells, and innate lymphoid cells (ILCs), as well as non-immune cells, including keratinocytes, fibroblasts, endothelial cells, and platelets in the initiation, development, and progression of psoriasis.

Dendrite intercalation between epidermal cells tunes nociceptor sensitivity to mechanical stimuli in Drosophila larvae

An animal's skin provides a first point of contact with the sensory environment, including noxious cues that elicit protective behavioral responses. Nociceptive somatosensory neurons densely innervate and intimately interact with epidermal cells to receive these cues, however the mechanisms by which epidermal interactions shape processing of noxious inputs is still poorly understood. Here, we identify a role for dendrite intercalation between epidermal cells in tuning sensitivity of Drosophila larvae to noxious mechanical stimuli. In wild-type larvae, dendrites of nociceptive class IV da neurons intercalate between epidermal cells at apodemes, which function as body wall muscle attachment sites, but not at other sites in the epidermis. From a genetic screen we identified miR-14 as a regulator of dendrite positioning in the epidermis: miR-14 is expressed broadly in the epidermis but not in apodemes, and miR-14 inactivation leads to excessive apical dendrite intercalation between epidermal cells. We found that miR-14 regulates expression and distribution of the epidermal Innexins ogre and Inx2 and that these epidermal gap junction proteins restrict epidermal dendrite intercalation. Finally, we found that altering the extent of epidermal dendrite intercalation had corresponding effects on nociception: increasing epidermal intercalation sensitized larvae to noxious mechanical inputs and increased mechanically evoked calcium responses in nociceptive neurons, whereas reducing epidermal dendrite intercalation had the opposite effects. Altogether, these studies identify epidermal dendrite intercalation as a mechanism for mechanical coupling of nociceptive neurons to the epidermis, with nociceptive sensitivity tuned by the extent of intercalation.

Short cell cycle duration is a phenotype of human epidermal stem cells

BACKGROUND

A traditional view is that stem cells (SCs) divide slowly. Meanwhile, both embryonic and pluripotent SCs display a shorter cell cycle duration (CCD) in comparison to more committed progenitors (CPs).

METHODS

We examined the in vitro proliferation and cycling behavior of somatic adult human cells using live cell imaging of passage zero keratinocytes and single-cell RNA sequencing.

RESULTS

We found two populations of keratinocytes: those with short CCD and protracted near exponential growth, and those with long CCD and terminal differentiation. Applying the ergodic principle, the comparative numbers of cycling cells in S phase in an enriched population of SCs confirmed a shorter CCD than CPs. Further, analysis of single-cell RNA sequencing of cycling adult human keratinocyte SCs and CPs indicated a shortening of both G1 and G2M phases in the SC.

CONCLUSIONS

Contrary to the pervasive paradigm, SCs progress through cell cycle more quickly than more differentiated dividing CPs. Thus, somatic human adult keratinocyte SCs may divide infrequently, but divide rapidly when they divide. Additionally, it was found that SC-like proliferation persisted in vitro.

Identification and validation of RNA-binding protein SLC3A2 regulates melanocyte ferroptosis in vitiligo by integrated analysis of single-cell and bulk RNA-sequencing

BACKGROUND

The pathogenesis of vitiligo remains unclear. The genes encoding vitiligo-related RNA-binding proteins (RBPs) and their underlying pathogenic mechanism have not been determined.

RESULTS

Single-cell transcriptome sequencing (scRNA-seq) data from the CNCB database was obtained to identify distinct cell types and subpopulations and the relative proportion changes in vitiligo and healthy samples. We identified 14 different cell types and 28 cell subpopulations. The proportion of each cell subpopulation significantly differed between the patients with vitiligo and healthy groups. Using RBP genes for unsupervised clustering, we obtained the specific RBP genes of different cell types in vitiligo and healthy groups. The RBP gene expression was highly heterogeneous; there were significant differences in some cell types, such as keratinocytes, Langerhans, and melanocytes, while there were no significant differences in other cells, such as T cells and fibroblasts, in the two groups. The melanocyte-specific RBP genes were enriched in the apoptosis and immune-related pathways in the patients with vitiligo. Combined with the bulk RNA-seq data of melanocytes, key RBP genes related to melanocytes were identified, including eight upregulated RBP genes (CDKN2A, HLA-A, RPL12, RPL29, RPL31, RPS19, RPS21, and RPS28) and one downregulated RBP gene (SLC3A2). Cell experiments were conducted to explore the role of the key RBP gene SLC3A2 in vitiligo. Cell experiments confirmed that melanocyte proliferation decreased, whereas apoptosis increased, after SLC3A2 knockdown. SLC3A2 knockdown in melanocytes also decreased the SOD activity and melanin content; increased the Fe2+, ROS, and MDA content; significantly increased the expression levels of TYR and COX2; and decreased the expression levels of glutathione and GPX4.

CONCLUSION

We identified the RBP genes of different cell subsets in patients with vitiligo and confirmed that downregulating SLC3A2 can promote ferroptosis in melanocytes. These findings provide new insights into the pathogenesis of vitiligo.

Translational frontiers: insight from lymphatics in skin regeneration

The remarkable regenerative ability of the skin, governed by complex molecular mechanisms, offers profound insights into the skin repair processes and the pathogenesis of various dermatological conditions. This understanding, derived from studies in human skin and various model systems, has not only deepened our knowledge of skin regeneration but also facilitated the development of skin substitutes in clinical practice. Recent research highlights the crucial role of lymphatic vessels in skin regeneration. Traditionally associated with fluid dynamics and immune modulation, these vessels are now recognized for interacting with skin stem cells and coordinating regeneration. This Mini Review provides an overview of recent advancements in basic and translational research related to skin regeneration, focusing on the dynamic interplay between lymphatic vessels and skin biology. Key highlights include the critical role of stem cell-lymphatic vessel crosstalk in orchestrating skin regeneration, emerging translational approaches, and their implications for skin diseases. Additionally, the review identifies research gaps and proposes potential future directions, underscoring the significance of this rapidly evolving research arena.

The Reign of Follistatin in Tumors and Their Microenvironment: Implications for Drug Resistance

Follistatin (FST) is a potent neutralizer of the transforming growth factor-β superfamily and is associated with normal cellular programs and various hallmarks of cancer, such as proliferation, migration, angiogenesis, and immune evasion. The aberrant expression of FST by solid tumors is a well-documented observation, yet how FST influences tumor progression and therapy response remains unclear. The recent surge in omics data has revealed new insights into the molecular foundation underpinning tumor heterogeneity and its microenvironment, offering novel precision medicine-based opportunities to combat cancer. In this review, we discuss these recent FST-centric studies, thereby offering an updated perspective on the protean role of FST isoforms in shaping the complex cellular ecosystem of tumors and in mediating drug resistance.

Significance of stress keratin expression in normal and diseased epithelia

A group of keratin intermediate filament genes, the type II KRT6A-C and type I KRT16 and KRT17, are deemed stress responsive as they are induced in keratinocytes of surface epithelia in response to environmental stressors, in skin disorders (e.g., psoriasis) and in carcinomas. Monitoring stress keratins is widely used to identify keratinocytes in an activated state. Here, we analyze single-cell transcriptomic data from healthy and diseased human skin to explore the properties of stress keratins. Relative to keratins occurring in healthy skin, stress-induced keratins are expressed at lower levels and show lesser type I-type II pairwise regulation. Stress keratins do not "replace" the keratins expressed during normal differentiation nor reflect cellular proliferation. Instead, stress keratins are consistently co-regulated with genes with roles in differentiation, inflammation, and/or activation of innate immunity at the single-cell level. These findings provide a roadmap toward explaining the broad diversity and contextual regulation of keratins.

Skin single-cell transcriptomics reveals a core of sebaceous gland-relevant genes shared by mice and humans

BACKGROUND

Single-cell RNA sequencing (scRNA-seq) has been widely applied to dissect cellular heterogeneity in normal and diseased skin. Sebaceous glands, essential skin components with established functions in maintaining skin integrity and emerging roles in systemic energy metabolism, have been largely neglected in scRNA-seq studies.

METHODS

Departing from mouse and human skin scRNA-seq datasets, we identified gene sets expressed especially in sebaceous glands with the open-source R-package oposSOM.

RESULTS

The identified gene sets included sebaceous gland-typical genes as Scd3, Mgst1, Cidea, Awat2 and KRT7. Surprisingly, however, there was not a single overlap among the 100 highest, exclusively in sebaceous glands expressed transcripts in mouse and human samples. Notably, both species share a common core of only 25 transcripts, including mitochondrial and peroxisomal genes involved in fatty acid, amino acid, and glucose processing, thus highlighting the intense metabolic rate of this gland.

CONCLUSIONS

This study highlights intrinsic differences in sebaceous lipid synthesis between mice and humans, and indicates an important role for peroxisomal processes in this context. Our data also provides attractive starting points for experimentally addressing novel candidates regulating sebaceous gland homeostasis.

Generating detailed intercellular communication patterns in psoriasis at the single-cell level using social networking, pattern recognition, and manifold learning methods to optimize treatment strategies

Psoriasis, a complex and recurrent chronic inflammatory skin disease involving various inflammatory cell types, requires effective cell communication to maintain the homeostatic balance of inflammation. However, patterns of communication at the single-cell level have not been systematically investigated. In this study, we employed social network analysis tools, pattern recognition, and manifold learning to compare molecular communication features between psoriasis cells and normal skin cells. Utilizing a process that facilitates the discovery of cell type-specific regulons, we analyzed internal regulatory networks among different cells in psoriasis. Advanced techniques for the quantitative detection of non-targeted proteins in pathological tissue sections were employed to demonstrate protein expression. Our findings revealed a synergistic interplay among the communication signals of immune cells in psoriasis. B-cells were activated, while Langerhans cells shifted into the primary signaling output mode to fulfill antigen presentation, mediating T-cell immunity. In contrast to normal skin cells, psoriasis cells shut down numerous signaling pathways, influencing the balance of skin cell renewal and differentiation. Additionally, we identified a significant number of active cell type-specific regulons of resident immune cells around the hair follicle. This study unveiled the molecular communication features of the hair follicle cell-psoriasis axis, showcasing its potential for therapeutic targeting at the single-cell level. By elucidating the pattern of immune cell communication in psoriasis and identifying new molecular features of the hair follicle cell-psoriasis axis, our findings present innovative strategies for drug targeting to enhance psoriasis treatment.

Interaction of human keratinocytes and nerve fiber terminals at the neuro-cutaneous unit

Traditionally, peripheral sensory neurons are assumed as the exclusive transducers of external stimuli. Current research moves epidermal keratinocytes into focus as sensors and transmitters of nociceptive and non-nociceptive sensations, tightly interacting with intraepidermal nerve fibers at the neuro-cutaneous unit. In animal models, epidermal cells establish close contacts and ensheath sensory neurites. However, ultrastructural morphological and mechanistic data examining the human keratinocyte-nerve fiber interface are sparse. We investigated this exact interface in human skin applying super-resolution array tomography, expansion microscopy, and structured illumination microscopy. We show keratinocyte ensheathment of afferents and adjacent connexin 43 contacts in native skin and have applied a pipeline based on expansion microscopy to quantify these parameter in skin sections of healthy participants versus patients with small fiber neuropathy. We further derived a fully human co-culture system, visualizing ensheathment and connexin 43 plaques in vitro. Unraveling human intraepidermal nerve fiber ensheathment and potential interaction sites advances research at the neuro-cutaneous unit. These findings are crucial on the way to decipher the mechanisms of cutaneous nociception.

Multiscale spatial mapping of cell populations across anatomical sites in healthy human skin and basal cell carcinoma

Our understanding of how human skin cells differ according to anatomical site and tumour formation is limited. To address this, we have created a multiscale spatial atlas of healthy skin and basal cell carcinoma (BCC), incorporating in vivo optical coherence tomography, single-cell RNA sequencing, spatial global transcriptional profiling, and in situ sequencing. Computational spatial deconvolution and projection revealed the localisation of distinct cell populations to specific tissue contexts. Although cell populations were conserved between healthy anatomical sites and in BCC, mesenchymal cell populations including fibroblasts and pericytes retained signatures of developmental origin. Spatial profiling and in silico lineage tracing support a hair follicle origin for BCC and demonstrate that cancer-associated fibroblasts are an expansion of a POSTN+ subpopulation associated with hair follicles in healthy skin. RGS5+ pericytes are also expanded in BCC suggesting a role in vascular remodelling. We propose that the identity of mesenchymal cell populations is regulated by signals emanating from adjacent structures and that these signals are repurposed to promote the expansion of skin cancer stroma. The resource we have created is publicly available in an interactive format for the research community.

Single-cell sequencing highlights heterogeneity and malignant progression in actinic keratosis and cutaneous squamous cell carcinoma

Cutaneous squamous cell carcinoma (cSCC) is the second most frequent of the keratinocyte-derived malignancies with actinic keratosis (AK) as a precancerous lesion. To comprehensively delineate the underlying mechanisms for the whole progression from normal skin to AK to invasive cSCC, we performed single-cell RNA sequencing (scRNA-seq) to acquire the transcriptomes of 138,982 cells from 13 samples of six patients including AK, squamous cell carcinoma in situ (SCCIS), cSCC, and their matched normal tissues, covering comprehensive clinical courses of cSCC. We identified diverse cell types, including important subtypes with different gene expression profiles and functions in major keratinocytes. In SCCIS, we discovered the malignant subtypes of basal cells with differential proliferative and migration potential. Differentially expressed genes (DEGs) analysis screened out multiple key driver genes including transcription factors along AK to cSCC progression. Immunohistochemistry (IHC)/immunofluorescence (IF) experiments and single-cell ATAC sequencing (scATAC-seq) data verified the expression changes of these genes. The functional experiments confirmed the important roles of these genes in regulating cell proliferation, apoptosis, migration, and invasion in cSCC tumor. Furthermore, we comprehensively described the tumor microenvironment (TME) landscape and potential keratinocyte-TME crosstalk in cSCC providing theoretical basis for immunotherapy. Together, our findings provide a valuable resource for deciphering the progression from AK to cSCC and identifying potential targets for anticancer treatment of cSCC.

Epigenetic switch reshapes epithelial progenitor cell signatures and drives inflammatory pathogenesis in hidradenitis suppurativa

Hidradenitis suppurativa (HS) is a complex inflammatory skin disease with undefined mechanistic underpinnings. Here, we investigated HS epithelial cells and demonstrated that HS basal progenitors modulate their lineage restriction and give rise to pathogenic keratinocyte clones, resulting in epidermal hyperproliferation and dysregulated inflammation in HS. When comparing to healthy epithelial stem/progenitor cells, in HS, we identified changes in gene signatures that revolve around the mitotic cell cycle, DNA damage response and repair, as well as cell-cell adhesion and chromatin remodeling. By reconstructing cell differentiation trajectory and CellChat modeling, we identified a keratinocyte population specific to HS. This population is marked by S100A7/8/9 and KRT6 family members, triggering IL1, IL10, and complement inflammatory cascades. These signals, along with HS-specific proinflammatory cytokines and chemokines, contribute to the recruitment of certain immune cells during the disease progression. Furthermore, we revealed a previously uncharacterized role of S100A8 in regulating the local chromatin environment of target loci in HS keratinocytes. Through the integration of genomic and epigenomic datasets, we identified genome-wide chromatin rewiring alongside the switch of transcription factors (TFs), which mediated HS transcriptional profiles. Importantly, we identified numerous clinically relevant inflammatory enhancers and their coordinated TFs in HS basal CD49fhigh cells. The disruption of the S100A enhancer using the CRISPR/Cas9-mediated approach or the pharmacological inhibition of the interferon regulatory transcription factor 3 (IRF3) efficiently reduced the production of HS-associated inflammatory regulators. Our study not only uncovers the plasticity of epidermal progenitor cells in HS but also elucidates the epigenetic mechanisms underlying HS pathogenesis.

Keratinocytes of the Upper Epidermis and Isthmus of Hair Follicles Express Hemoglobin mRNA and Protein

The epidermis, the keratinized stratified squamous epithelium surrounding the body surface, offers a valuable framework to investigate how terrestrial animals overcome environmental stresses. However, the mechanisms underlying epidermal barrier function remain nebulous. In this study, we examined genes highly expressed in the human and mouse upper epidermis, the outer frontier that induces various barrier-related genes. Transcriptome analysis revealed that the messenger RNA level of hemoglobin α (HBA), an oxygen carrier in erythroid cells, was enriched in the upper epidermis compared with that in the whole epidermis. Immunostaining analysis confirmed HBA protein expression in human and mouse keratinocytes (KCs) of the stratum spinosum and stratum granulosum. HBA was also expressed in hair follicle KCs in the isthmus region; its expression levels were more prominent than those in interfollicular KCs. HBA expression was not observed in noncutaneous keratinized stratified squamous epithelia of mice, for example, the vagina, esophagus, and forestomach. HBA expression was upregulated in human epidermal KC cultures after UV irradiation, a major cause of skin-specific oxidative stress. Furthermore, HBA knockdown increased UV-induced production of ROS in primary KCs. Our findings suggest that epidermal HBA expression is induced by oxidative stress and acts as an antioxidant, contributing to skin barrier function.

Single-cell RNA sequencing of human epidermis identifies Lunatic fringe as a novel regulator of the stem cell compartment

Single-cell RNA sequencing (scRNA-seq) of human skin provides a tool for validating observations from in vitro experimental models. By analyzing a published dataset of healthy adult epidermis, we confirm that the basal epidermal layer is heterogeneous, and three subpopulations of non-dividing cells can be distinguished. We show that Delta-like ligand 1 (DLL1) is expressed in a subset of basal cells previously identified as stem cells in cultured human keratinocytes and map the distribution of other Notch ligands and receptors to specific epidermal cell compartments. Although DLL1 is expressed at low levels, it is expressed in the same cell state as the Notch regulator, Lunatic -fringe (LFNG, O-fucosylpeptide 3-beta-N-acetylglucosaminyltransferase). Overexpression of LFNG amplifies the effects of DLL1 in cultured keratinocytes, increasing proliferation and colony-forming ability. We conclude that using scRNA-seq resources from healthy human skin not only validates previous experimental data but allows formulation of testable new hypotheses.

Architecture of skin inflammation in psoriasis revealed by spatial transcriptomics

Psoriasis is a chronic inflammatory skin disease, thought to be predominantly mediated by TH17 cells. Significance of other inflammatory pathways and the innate immune system is not well understood and the spatial heterogeneity of inflammation in the skin has largely been overlooked. Our aim was to create a comprehensive map of skin inflammation in psoriasis, exploring the tissue patterning of inflammation. In situ whole transcriptome sequencing (spatial sequencing) was performed on lesional psoriatic skin in four patients with moderate-to-severe disease to quantify all expressed genes within a tissue section. Transcriptional analysis revealed three major inflammatory niches in psoriasis skin, each with distinct cytokine circuits and chemokines: the hyperplastic epidermis, upper (papillary) dermis, and reticular dermis. Interestingly, key cytokines such as IL-23, IL-17 s, and TNFα were not notably present in the skin's transcriptomic signature. Unexpectedly, IL-32 showed strong expression in the dermis. Our findings underscore the complexity of psoriatic inflammation, highlighting its architectural heterogeneity and the roles of innate cytokines. Both IL-32 and IL-1 family cytokines appear to play critical roles in the dermal and epidermal inflammation, respectively, and may provide pharmacological targets to improve the control of the inflammatory process.

Single-Cell and Spatial Transcriptomic Analysis of Human Skin Delineates Intercellular Communication and Pathogenic Cells

Epidermal homeostasis is governed by a balance between keratinocyte proliferation and differentiation with contributions from cell-cell interactions, but conserved or divergent mechanisms governing this equilibrium across species and how an imbalance contributes to skin disease are largely undefined. To address these questions, human skin single-cell RNA sequencing and spatial transcriptomics data were integrated and compared with mouse skin data. Human skin cell-type annotation was improved using matched spatial transcriptomics data, highlighting the importance of spatial context in cell-type identity, and spatial transcriptomics refined cellular communication inference. In cross-species analyses, we identified a human spinous keratinocyte subpopulation that exhibited proliferative capacity and a heavy metal processing signature, which was absent in mouse and may account for species differences in epidermal thickness. This human subpopulation was expanded in psoriasis and zinc-deficiency dermatitis, attesting to disease relevance and suggesting a paradigm of subpopulation dysfunction as a hallmark of the disease. To assess additional potential subpopulation drivers of skin diseases, we performed cell-of-origin enrichment analysis within genodermatoses, nominating pathogenic cell subpopulations and their communication pathways, which highlighted multiple potential therapeutic targets. This integrated dataset is encompassed in a publicly available web resource to aid mechanistic and translational studies of normal and diseased skin.

Single-cell transcriptomic atlas-guided development of CAR-T cells for the treatment of acute myeloid leukemia

Chimeric antigen receptor T cells (CAR-T cells) have emerged as a powerful treatment option for individuals with B cell malignancies but have yet to achieve success in treating acute myeloid leukemia (AML) due to a lack of safe targets. Here we leveraged an atlas of publicly available RNA-sequencing data of over 500,000 single cells from 15 individuals with AML and tissue from 9 healthy individuals for prediction of target antigens that are expressed on malignant cells but lacking on healthy cells, including T cells. Aided by this high-resolution, single-cell expression approach, we computationally identify colony-stimulating factor 1 receptor and cluster of differentiation 86 as targets for CAR-T cell therapy in AML. Functional validation of these established CAR-T cells shows robust in vitro and in vivo efficacy in cell line- and human-derived AML models with minimal off-target toxicity toward relevant healthy human tissues. This provides a strong rationale for further clinical development.

Single-cell RNA-seq reveals actinic keratosis-specific keratinocyte subgroups and their crosstalk with secretory-papillary fibroblasts

BACKGROUND AND AIM

Actinic keratosis (AK) represents an intraepidermal malignant neoplasm with the proliferation of atypical keratinocytes. AK lesions are regarded as early in situ squamous cell carcinomas (SCCs) having the potential to progress into invasive SCC (iSCC) and metastasize, causing death. This study aimed to investigate the heterogeneity of keratinocytes and how this heterogeneity promoted AK development and progression.

METHODS

We employed single-cell RNA sequencing (scRNA-seq) to examine the heterogeneity of keratinocytes and dermal fibroblast clusters in AKs and adjacent normal skins. Cell clustering, pseudotime trajectory construction, gene ontology enrichment analysis, transcription factor network analysis, and cell-cell communication were used to investigate the heterogeneity of keratinocytes in AK. The cellular identity and function were verified by immunohistochemical and immunofluorescence staining.

RESULTS

Using scRNA-seq, we revealed 13 keratinocyte subgroups (clusters 0-12) in AK tissues and characterized 2 AK-specific clusters. Cluster 9 displayed high levels of IL1R2 and WFDC2, and cluster 11 showed high levels of FADS2 and FASN. The percentages of cells in these two clusters significantly increased in AK compared with normal tissues. The existence and spatial localization of AK-specific IL1R2+WFDC2+ cluster were verified by immunohistochemical and immunofluorescence staining. Functional studies indicated that the genes identified in the IL1R2+WFDC2+ cluster were crucial for epithelial cell proliferation, migration, and angiogenesis. Further immunofluorescent staining revealed the interactions between AK-specific keratinocytes and secretory-papillary fibroblasts mainly through ANGPTL4-ITGA5 signalling pathway rarely seen in normal tissues.

CONCLUSION

The findings of this study might help better understand AK pathogenesis.

Multi-modal skin atlas identifies a multicellular immune-stromal community associated with altered cornification and specific T cell expansion in atopic dermatitis

In healthy skin, a cutaneous immune system maintains the balance between tolerance towards innocuous environmental antigens and immune responses against pathological agents. In atopic dermatitis (AD), barrier and immune dysfunction result in chronic tissue inflammation. Our understanding of the skin tissue ecosystem in AD remains incomplete with regard to the hallmarks of pathological barrier formation, and cellular state and clonal composition of disease-promoting cells. Here, we generated a multi-modal cell census of 310,691 cells spanning 86 cell subsets from whole skin tissue of 19 adult individuals, including non-lesional and lesional skin from 11 AD patients, and integrated it with 396,321 cells from four studies into a comprehensive human skin cell atlas in health and disease. Reconstruction of human keratinocyte differentiation from basal to cornified layers revealed a disrupted cornification trajectory in AD. This disrupted epithelial differentiation was associated with signals from a unique immune and stromal multicellular community comprised of MMP12 + dendritic cells (DCs), mature migratory DCs, cycling ILCs, NK cells, inflammatory CCL19 + IL4I1 + fibroblasts, and clonally expanded IL13 + IL22 + IL26 + T cells with overlapping type 2 and type 17 characteristics. Cell subsets within this immune and stromal multicellular community were connected by multiple inter-cellular positive feedback loops predicted to impact community assembly and maintenance. AD GWAS gene expression was enriched both in disrupted cornified keratinocytes and in cell subsets from the lesional immune and stromal multicellular community including IL13 + IL22 + IL26 + T cells and ILCs, suggesting that epithelial or immune dysfunction in the context of the observed cellular communication network can initiate and then converge towards AD. Our work highlights specific, disease-associated cell subsets and interactions as potential targets in progression and resolution of chronic inflammation.

IL-31-dependent neurogenic inflammation restrains cutaneous type 2 immune response in allergic dermatitis

Despite robust literature associating IL-31 with pruritic inflammatory skin diseases, its influence on cutaneous inflammation and the interplay between inflammatory and neurosensory pathways remain unmapped. Here, we examined the consequences of disrupting Il31 and its receptor Il31ra in a mouse model of house dust mite (HDM)-induced allergic dermatitis. Il31-deficient mice displayed a deficit in HDM dermatitis-associated scratching, consistent with its well-established role as a pruritogen. In contrast, Il31 deficiency increased the number and proportion of cutaneous type 2 cytokine-producing CD4+ T cells and serum IgE in response to HDM. Furthermore, Il4ra+ monocytes and macrophages capable of fueling a feedforward type 2 inflammatory loop were selectively enriched in Il31ra-deficient HDM dermatitis skin. Thus, IL-31 is not strictly a proinflammatory cytokine but rather an immunoregulatory factor that limits the magnitude of type 2 inflammatory responses in skin. Our data support a model wherein IL-31 activation of IL31RA+ pruritoceptors triggers release of calcitonin gene-related protein (CGRP), which can mediate neurogenic inflammation, inhibit CD4+ T cell proliferation, and reduce T cell production of the type 2 cytokine IL-13. Together, these results illustrate a previously unrecognized neuroimmune pathway that constrains type 2 tissue inflammation in the setting of chronic cutaneous allergen exposure and may explain paradoxical dermatitis flares in atopic patients treated with anti-IL31RA therapy.

Keratin 17- and PKCα-dependent transient amplification of neutrophil influx after repeated stress to the skin

Neutrophils contribute to the pathogenesis of chronic inflammatory skin diseases. Little is known about the source and identity of the signals mediating their recruitment in inflamed skin. We used the phorbol ester TPA and UVB, alone or in combination, to induce sterile inflammation in mouse skin and assess whether keratinocyte-derived signals impact neutrophil recruitment. A single TPA treatment results in a neutrophil influx in the dermis that peaks at 12h and resolves within 24h. A second TPA treatment or a UVB challenge, when applied at 24h but not 48h later, accelerates, amplifies, and prolongs neutrophil infiltration. This transient amplification response (TAR) is mediated by local signals in inflamed skin, can be recapitulated in ex vivo culture, and involves the K17-dependent sustainment of protein kinase Cα (PKCα) activity and release of neutrophil chemoattractants by stressed keratinocytes. We show that K17 binds RACK1, a scaffold essential for PKCα activity. Finally, analyses of RNAseq data reveal the presence of a transcriptomic signature consistent with TAR and PKCα activation in chronic inflammatory skin diseases. These findings uncover a novel, transient, and keratin-dependent mechanism that amplifies neutrophil recruitment to the skin under stress, with direct implications for inflammatory skin disorders.

The mycosis fungoides cutaneous microenvironment shapes dysfunctional cell trafficking, antitumor immunity, matrix interactions, and angiogenesis

Malignant T lymphocyte proliferation in mycosis fungoides (MF) is largely restricted to the skin, implying that malignant cells are dependent on their specific cutaneous tumor microenvironment (TME), including interactions with non-malignant immune and stromal cells, cytokines, and other immunomodulatory factors. To explore these interactions, we performed a comprehensive transcriptome analysis of the TME in advanced-stage MF skin tumors by single-cell RNA sequencing. Our analysis identified cell-type compositions, cellular functions, and cell-to-cell interactions in the MF TME that were distinct from those from healthy skin and benign dermatoses. While patterns of gene expression were common among patient samples, high transcriptional diversity was also observed in immune and stromal cells, with dynamic interactions and crosstalk between these cells and malignant T lymphocytes. This heterogeneity mapped to processes such as cell trafficking, matrix interactions, angiogenesis, immune functions, and metabolism that affect cancer cell growth, migration, and invasion, as well as antitumor immunity. By comprehensively characterizing the transcriptomes of immune and stromal cells within the cutaneous microenvironment of individual MF tumors, we have identified patterns of dysfunction common to all tumors that represent a resource for identifying candidates with therapeutic potential as well as patient-specific heterogeneity that has important implications for personalized disease management.

Characterization of 3D organotypic epithelial tissues reveals tonsil-specific differences in tonic interferon signaling

Three-dimensional (3D) culturing techniques can recapitulate the stratified nature of multicellular epithelial tissues. Organotypic 3D epithelial tissue culture methods have several applications, including the study of tissue development and function, drug discovery and toxicity testing, host-pathogen interactions, and the development of tissue-engineered constructs for use in regenerative medicine. We grew 3D organotypic epithelial tissues from foreskin, cervix, and tonsil-derived primary cells and characterized the transcriptome of these in vitro tissue equivalents. Using the same 3D culturing method, all three tissues yielded stratified squamous epithelium, validated histologically using basal and superficial epithelial cell markers. The goal of this study was to use RNA-seq to compare gene expression patterns in these three types of epithelial tissues to gain a better understanding of the molecular mechanisms underlying their function and identify potential therapeutic targets for various diseases. Functional profiling by over-representation and gene set enrichment analysis revealed tissue-specific differences: i.e., cutaneous homeostasis and lipid metabolism in foreskin, extracellular matrix remodeling in cervix, and baseline innate immune differences in tonsil. Specifically, tonsillar epithelia may play an active role in shaping the immune microenvironment of the tonsil balancing inflammation and immune responses in the face of constant exposure to microbial insults. Overall, these data serve as a resource, with gene sets made available for the research community to explore, and as a foundation for understanding the epithelial heterogeneity and how it may impact their in vitro use. An online resource is available to investigate these data (https://viz.datascience.arizona.edu/3DEpiEx/).

Single-cell transcriptional profiling reveals immunomodulatory properties of stromal and epithelial cells in periodontal immune milieu with diabetes in rats

Periodontitis is the sixth major complication of diabetes. Gingiva, as an important component of periodontal tissues, serves as the first defense barrier against infectious stimuli. However, relatively little is known about cellular heterogeneity and cell-specific changes in gingiva in response to diabetes-associated periodontitis. To characterize molecular changes linking diabetes with periodontitis, we profiled single-cell transcriptome analyses of a total of 45,259 cells from rat gingiva with periodontitis under normoglycemic and diabetic condition. The single-cell profiling revealed that stromal and epithelial cells of gingiva contained inflammation-related subclusters enriched in functions of immune cell recruitment. Compared to normoglycemic condition, diabetes led to a reduction in epithelial basal cells, fibroblasts and smooth muscle cells in gingiva with periodontitis. Analysis of differentially expressed genes indicated that stromal and epithelial populations were reprogrammed towards pro-inflammatory phenotypes promoting immune cell recruitment in diabetes-related periodontitis. In aspect of immune cells, diabetes prominently enhanced neutrophil and M1 macrophage infiltration in periodontitis lesions. Cell-cell communications revealed enhanced crosstalk between stromal/epithelial cells and immune cells mediating by chemokine/chemokine receptor interplay in diabetes-associated periodontitis. Our findings deconvolved cellular heterogeneity of rat gingiva associated with periodontitis and diabetes, uncovered altered immune milieu caused by the disease, and revealed immunomodulatory functions of stromal and epithelial cells in gingival immune niche. The present study improves the understanding of the link between the diabetes and periodontitis and helps in formulating precise therapeutic strategies for diabetes-enhanced periodontitis.

Leprosy-specific subsets of macrophages and Schwann cells identified by single-cell RNA-sequencing

In Mycobacterium leprae (M. leprae)-infection, inflammatory cells' subsets and dynamics as well as the interactions with Schwann cells have remained elusive. We investigated individual cells in M. leprae-inoculated nude mice by single-cell RNA-sequencing (scRNA-seq). For macrophages, we dissected two M1-like subsets and five M2-like subsets, where lipid-associated signatures were pervasive in both M1-like and M2-like subsets. There were four macrophage trajectories showing: (i) pro-inflammatory (M1), (ii) lipid metabolism-related (M2), (iii) anti-inflammatory (M2), and (iv) interferon-stimulated gene-related (M2) fates. They displayed early divergence without ever rejoining along the paths, suggesting simultaneous or continuous stimuli for macrophage activation in leprosy. The scRNA-seq predicted Schwann cell-macrophage interactions (Notch1-Jag1, Plxnb1-Sema4d interactions). An immature Schwann cell subset showing Tfap2a expression was identified, indicating Schwann cell dedifferentiation in leprosy tissues. Expressions of Notch1, Jag1, Plxnb1, Sema4d, and Tfap2a were validated in mouse or human leprosy tissues by immunohistochemistry. We identified both pro-inflammatory and inflammation-resolution signatures, where lipid-associated signatures were pervasive to the macrophages, representing leprosy-specific macrophage states for prolonged and repeated episodes of inflammation and resolution. Our study identified refined molecular states and interactions of macrophages and Schwann cells, suggesting novel insights into the pathogenesis of unhealed inflammation with neuropathy and potential therapeutic targets for leprosy.

Integrative single-cell transcriptomic investigation unveils long non-coding RNAs associated with localized cellular inflammation in psoriasis

Psoriasis is a complex, chronic autoimmune disorder predominantly affecting the skin. Accumulating evidence underscores the critical role of localized cellular inflammation in the development and persistence of psoriatic skin lesions, involving cell types such as keratinocytes, mesenchymal cells, and Schwann cells. However, the underlying mechanisms remain largely unexplored. Long non-coding RNAs (lncRNAs), known to regulate gene expression across various cellular processes, have been particularly implicated in immune regulation. We utilized our neural-network learning pipeline to integrate 106,675 cells from healthy human skin and 79,887 cells from psoriatic human skin. This formed the most extensive cell transcriptomic atlas of human psoriatic skin to date. The robustness of our reclassified cell-types, representing full-layer zonation in human skin, was affirmed through neural-network learning-based cross-validation. We then developed a publicly available website to present this integrated dataset. We carried out analysis for differentially expressed lncRNAs, co-regulated gene patterns, and GO-bioprocess enrichment, enabling us to pinpoint lncRNAs that modulate localized cellular inflammation in psoriasis at the single-cell level. Subsequent experimental validation with skin cell lines and primary cells from psoriatic skin confirmed these lncRNAs' functional role in localized cellular inflammation. Our study provides a comprehensive cell transcriptomic atlas of full-layer human skin in both healthy and psoriatic conditions, unveiling a new regulatory mechanism that governs localized cellular inflammation in psoriasis and highlights the therapeutic potential of lncRNAs in this disease's management.

Distinct mechanisms for sebaceous gland self-renewal and regeneration provide durability in response to injury

Sebaceous glands (SGs) release oils that protect our skin, but how these glands respond to injury has not been previously examined. Here, we report that SGs are largely self-renewed by dedicated stem cell pools during homeostasis. Using targeted single-cell RNA sequencing, we uncovered both direct and indirect paths by which resident SG progenitors ordinarily differentiate into sebocytes, including transit through a Krt5+PPARγ+ transitional basal cell state. Upon skin injury, however, SG progenitors depart their niche, reepithelialize the wound, and are replaced by hair-follicle-derived stem cells. Furthermore, following targeted genetic ablation of >99% of SGs from dorsal skin, these glands unexpectedly regenerate within weeks. This regenerative process is mediated by alternative stem cells originating from the hair follicle bulge, is dependent upon FGFR2 signaling, and can be accelerated by inducing hair growth. Altogether, our studies demonstrate that stem cell plasticity promotes SG durability following injury.

Dendrite intercalation between epidermal cells tunes nociceptor sensitivity to mechanical stimuli in Drosophila larvae

An animal's skin provides a first point of contact with the sensory environment, including noxious cues that elicit protective behavioral responses. Nociceptive somatosensory neurons densely innervate and intimately interact with epidermal cells to receive these cues, however the mechanisms by which epidermal interactions shape processing of noxious inputs is still poorly understood. Here, we identify a role for dendrite intercalation between epidermal cells in tuning sensitivity of Drosophila larvae to noxious mechanical stimuli. In wild-type larvae, dendrites of nociceptive class IV da neurons intercalate between epidermal cells at apodemes, which function as body wall muscle attachment sites, but not at other sites in the epidermis. From a genetic screen we identified miR-14 as a regulator of dendrite positioning in the epidermis: miR-14 is expressed broadly in the epidermis but not in apodemes, and miR-14 inactivation leads to excessive apical dendrite intercalation between epidermal cells. We found that miR-14 regulates expression and distribution of the epidermal Innexins ogre and Inx2 and that these epidermal gap junction proteins restrict epidermal dendrite intercalation. Finally, we found that altering the extent of epidermal dendrite intercalation had corresponding effects on nociception: increasing epidermal intercalation sensitized larvae to noxious mechanical inputs and increased mechanically evoked calcium responses in nociceptive neurons, whereas reducing epidermal dendrite intercalation had the opposite effects. Altogether, these studies identify epidermal dendrite intercalation as a mechanism for mechanical coupling of nociceptive neurons to the epidermis, with nociceptive sensitivity tuned by the extent of intercalation.

Peripheral mechanisms of peripheral neuropathic pain

Peripheral neuropathic pain (PNP), neuropathic pain that arises from a damage or disease affecting the peripheral nervous system, is associated with an extremely large disease burden, and there is an increasing and urgent need for new therapies for treating this disorder. In this review we have highlighted therapeutic targets that may be translated into disease modifying therapies for PNP associated with peripheral neuropathy. We have also discussed how genetic studies and novel technologies, such as optogenetics, chemogenetics and single-cell RNA-sequencing, have been increasingly successful in revealing novel mechanisms underlying PNP. Additionally, consideration of the role of non-neuronal cells and communication between the skin and sensory afferents is presented to highlight the potential use of drug treatment that could be applied topically, bypassing drug side effects. We conclude by discussing the current difficulties to the development of effective new therapies and, most importantly, how we might improve the translation of targets for peripheral neuropathic pain identified from studies in animal models to the clinic.

Single-cell analysis of human prepuce reveals dynamic changes in gene regulation and cellular communications

BACKGROUND

The cellular and molecular dynamics of human prepuce are crucial for understanding its biological and physiological functions, as well as the prevention of related genital diseases. However, the cellular compositions and heterogeneity of human prepuce at single-cell resolution are still largely unknown. Here we systematically dissected the prepuce of children and adults based on the single-cell RNA-seq data of 90,770 qualified cells.

RESULTS

We identified 15 prepuce cell subtypes, including fibroblast, smooth muscle cells, T/natural killer cells, macrophages, vascular endothelial cells, and dendritic cells. The proportions of these cell types varied among different individuals as well as between children and adults. Moreover, we detected cell-type-specific gene regulatory networks (GRNs), which could contribute to the unique functions of related cell types. The GRNs were also highly dynamic between the prepuce cells of children and adults. Our cell-cell communication network analysis among different cell types revealed a set of child-specific (e.g., CD96, EPO, IFN-1, and WNT signaling pathways) and adult-specific (e.g., BMP10, NEGR, ncWNT, and NPR1 signaling pathways) signaling pathways. The variations of GRNs and cellular communications could be closely associated with prepuce development in children and prepuce maintenance in adults.

CONCLUSIONS

Collectively, we systematically analyzed the cellular variations and molecular changes of the human prepuce at single-cell resolution. Our results gained insights into the heterogeneity of prepuce cells and shed light on the underlying molecular mechanisms of prepuce development and maintenance.

A Roadmap for a Consensus Human Skin Cell Atlas and Single-Cell Data Standardization

Single-cell technologies have become essential to driving discovery in both basic and translational investigative dermatology. Despite the multitude of available datasets, a central reference atlas of normal human skin, which can serve as a reference resource for skin cell types, cell states, and their molecular signatures, is still lacking. For any such atlas to receive broad acceptance, participation by many investigators during atlas construction is an essential prerequisite. As part of the Human Cell Atlas project, we have assembled a Skin Biological Network to build a consensus Human Skin Cell Atlas and outline a roadmap toward that goal. We define the drivers of skin diversity to be considered when selecting sequencing datasets for the atlas and list practical hurdles during skin sampling that can result in data gaps and impede comprehensive representation and technical considerations for tissue processing and computational analysis, the accounting for which should minimize biases in cell type enrichments and exclusions and decrease batch effects. By outlining our goals for Atlas 1.0, we discuss how it will uncover new aspects of skin biology.

Hedgehog receptors exert immune-surveillance roles in the epidermis across species

Hedgehog signaling plays pivotal roles in the development and homeostasis of epithelial barrier tissues. However, whether and how Hedgehog signaling directly regulates innate immunity in epithelial cells remains unknown. By utilizing C. elegans epidermis as the model, we found that several Hedgehog receptors are involved in cell-autonomous regulation of the innate immune response in the epidermis. Particularly, loss of the Patched family receptor induces aberrant up-regulation of epidermal antimicrobial peptides in a STAT-dependent manner. External or internal insult to the epidermis triggers rapid rearrangement of Patched distribution along the plasma membrane, indicating that the Hedgehog (Hh) receptor is likely involved in recognition and defense against epidermal damage. Loss of PTCH1 function in primary human keratinocytes and intact mouse skin also results in STAT-dependent immune activation. These findings reveal an evolutionally conserved immune-surveillance function of Hedgehog receptors and an insult-sensing and response strategy of epithelial tissues.

The dual role of cannabidiol on monocyte-derived dendritic cell differentiation and maturation

Introduction

Extracts and compounds isolated from hemp (Cannabis sativa) are increasingly gaining popularity in the treatment of a number of diseases, with topical formulations for dermatological conditions leading the way. Phytocannabinoids such as ( )-cannabidiol, ( )-cannabinol and ( )-Δ9-tetrahydrocannabivarin (CBD, CBN, and THCV, respectively), are present in variable amounts in the plant, and have been shown to have mostly anti-inflammatory effects both in vitro and in vivo, albeit dominantly in murine models. The role of phytocannabinoids in regulating responses of dendritic cells (DCs) remains unclear.

Methods

Our research aimed to investigate the effects of CBD, CBN, and THCV on human DCs differentiated from monocytes (moDCs). moDCs were treated with up to 10 μM of each phytocannabinoid, and their effects on viability, differentiation, and maturation were assessed both alone, and in conjunction with TLR agonists. The effects of CBD on cytokine production, T cell activation and polarization as well as the transcriptome of moDCs was also determined.

Results

Phytocannabinoids did not influence the viability of moDCs up to 10 μM, and only CBD had effects on maturational markers of moDCs, and neither compound influenced LPS-induced activation at 10 μM. Since only CBD had measurable effects on moDCs, in our subsequent experiments we tested the effect only of that pCB. On moDCs differentiated in the presence of CBD subsequent activation by LPS induced a markedly different, much more tolerogenic response. CBD-treated moDCs also produced significantly more interleukin (IL)-6, TNFα and, importantly, IL-10 in response to LPS, which shows a shift toward anti-inflammatory signaling, as well as a more robust secretory response in general. To rule out the possibility that these effects of CBD are specific to TLR4 signaling, we determined the effect of CBD on TLR7/8-induced maturation as well, and saw similar, although less marked responses. CBD-treated moDCs were also less efficient at activating naïve T cells after LPS stimulation, further supporting the tolerogenic effect of this phytocannabinoid on moDCs. Reactome pathway analysis showed an inflammatory response to LPS in moDCs, and to a lesser extent to CBD as well. In contrast CBD-treated moDCs responded to LPS with a shift towards a more tolerogenic phenotype, as IL-10 signaling was the most prominently induced pathway in this group.

Discussion

Our results show that CBD achieves an anti-inflammatory effect on adaptive immune responses only in the presence of an activating stimuli on moDCs by reprogramming cells during long-term treatment, and not through acute, short-term effects.

Non-IgE-reactive allergen peptides deteriorate the skin barrier in house dust mite-sensitized atopic dermatitis patients

Atopic dermatitis (AD) is a chronic inflammatory skin disease characterized by type 2 cytokine-driven skin inflammation and epithelial barrier dysfunction. The latter is believed to allow the increased penetration of chemicals, toxins, and allergens into the skin. House dust mite allergens, particularly Der p 2, are important triggers in sensitized individuals with AD; the precise actions of these allergens in epithelial biology remain, however, incompletely understood. In this study, we compared the effects of the protein allergen Der p 2 and a mix of non-IgE-reactive Der p 2 peptides on skin cells using patch tests in AD patients and healthy participants. We then analyzed mRNA expression profiles of keratinocytes by single-cell RNA-sequencing. We report that existing barrier deficiencies in the non-lesional skin of AD patients allow deep penetration of Der p 2 and its peptides, leading to local microinflammation. Der p 2 protein specifically upregulated genes involved in the innate immune system, stress, and danger signals in suprabasal KC. Der p 2 peptides further downregulated skin barrier genes, in particular the expression of genes involved in cell-matrix and cell-cell adhesion. Peptides also induced genes involved in hyperproliferation and caused disturbances in keratinocyte differentiation. Furthermore, inflammasome-relevant genes and IL18 were overexpressed, while KRT1 was downregulated. Our data suggest that Der p 2 peptides contribute to AD initiation and exacerbation by augmenting hallmark features of AD, such as skin inflammation, barrier disruption, and hyperplasia of keratinocytes.

Blank Spots in the Map of Human Skin: The Challenge for Xenotransplantation

Most of the knowledge about human skin homeostasis, development, wound healing, and diseases has been accumulated from human skin biopsy analysis by transferring from animal models and using different culture systems. Human-to-mouse xenografting is one of the fundamental approaches that allows the skin to be studied in vivo and evaluate the ongoing physiological processes in real time. Humanized animals permit the actual techniques for tracing cell fate, clonal analysis, genetic modifications, and drug discovery that could never be employed in humans. This review recapitulates the novel facts about mouse skin self-renewing, regeneration, and pathology, raises issues regarding the gaps in our understanding of the same options in human skin, and postulates the challenges for human skin xenografting.

Updated skin transcriptomic atlas depicted by reciprocal contribution of single-nucleus RNA sequencing and single-cell RNA sequencing

BACKGROUND

Single-cell RNA sequencing (scRNA-seq) has advanced our understanding of skin biology, but its utility is restricted by the requirement of fresh samples, inadequate dissociation-induced cell loss or death, and activation during tissue digestion. Single-nucleus RNA sequencing (snRNA-seq) can use frozen, hard-to-dissociate materials, which might be a promising method to circumvent the limitations of scRNA-seq for the skin tissue.

OBJECTIVE

To profile skin cells using snRNA-seq in parallel with scRNA-seq.

METHODS

We performed snRNA-seq in parallel with scRNA-seq for the bisected skin sample of one person and integrated previously published scRNA-seq data for analysis. We comparatively analyzed the differences in cell proportions and gene expression between the two methods. The differentiation trajectories of keratinocytes and fibroblasts were analyzed by Slingshot analysis.

RESULTS

snRNA-seq was less susceptible to contamination from mitochondrial and ribosomal RNA, and exhibited a greater capacity to detect transcription factors. snRNA-seq identified more spatially and functionally relevant keratinocyte clusters that constitute cell trajectories with expected differentiation dynamics. Novel markers, e.g., LYPD3, EMP2, and CSTB, were revealed for different differentiation stages of keratinocytes, and NFIB and GRHL1 were identified as transcription factors involving in the proliferation and functional differentiation of keratinocytes. Fibroblasts were found in a state of activation in scRNA-seq. And scRNA-seq detected a greater number of immune cells.

CONCLUSIONS

We generated an updated atlas of the skin transcriptome based on the reciprocal contribution of scRNA-seq and snRNA-seq.

Single-cell RNA sequencing defines disease-specific differences between chronic nodular prurigo and atopic dermatitis

BACKGROUND

Chronic nodular prurigo (CNPG) is an inflammatory skin disease that is maintained by a chronic itch-scratch cycle likely rooted in neuroimmunological dysregulation. This condition may be associated with atopy in some patients, and there are now promising therapeutic results from blocking type 2 cytokines such as IL-4, IL-13, and IL-31.

OBJECTIVES

This study aimed to improve the understanding of pathomechanisms underlying CNPG as well as molecular relationships between CNPG and atopic dermatitis (AD).

METHODS

We profiled skin lesions from patients with CNPG in comparison with AD and healthy control individuals using single-cell RNA sequencing combined with T-cell receptor sequencing.

RESULTS

We found type 2 immune skewing in both CNPG and AD, as evidenced by CD4+ helper T cells expressing IL13. However, only AD harbored an additional, oligoclonally expanded CD8A+IL9R+IL13+ cytotoxic T-cell population, and immune activation pathways were highly upregulated in AD, but less so in CNPG. Conversely, CNPG showed signatures of extracellular matrix organization, collagen synthesis, and fibrosis, including a unique population of CXCL14-IL24+ secretory papillary fibroblasts. Besides known itch mediators such as IL31 and oncostatin M, we also detected increased levels of neuromedin B in fibroblasts of CNPG lesions compared with AD and HC, with neuromedin B receptors detectable on some nerve endings.

CONCLUSIONS

These data show that CNPG does not harbor the strong disease-specific immune activation pathways that are typically found in AD but is rather characterized by upregulated stromal remodeling mechanisms that might have a direct impact on itch fibers.

Assessment of Treatment-Relevant Immune Biomarkers in Psoriasis and Atopic Dermatitis: Toward Personalized Medicine in Dermatology

Immunologically targeted therapies have revolutionized the treatment of inflammatory dermatoses, including atopic dermatitis and psoriasis. Although immunologic biomarkers hold great promise for personalized classification of skin disease and tailored therapy selection, there are no approved or widely used approaches for this in dermatology. This review summarizes the translational immunologic approaches to measuring treatment-relevant biomarkers in inflammatory skin conditions. Tape strip profiling, microneedle-based biomarker patches, molecular profiling from epidermal curettage, RNA in situ hybridization tissue staining, and single-cell RNA sequencing have been described. We discuss the advantages and limitations of each and open questions for the future of personalized medicine in inflammatory skin disease.

Microenvironmental and cell intrinsic factors governing human cDC2 differentiation and monocyte reprogramming

cDC2s occur abundantly in peripheral tissues and arise from circulating blood cDC2s. However, the factors governing cDC2 differentiation in tissues, especially under inflammatory conditions, remained poorly defined. We here found that psoriatic cDC2s express the efferocytosis receptor Axl and exhibit a bone morphogenetic protein (BMP) and p38MAPK signaling signature. BMP7, strongly expressed within the lesional psoriatic epidermis, cooperates with canonical TGF-β1 signaling for inducing Axl+cDC2s from blood cDC2s in vitro. Moreover, downstream induced p38MAPK promotes Axl+cDC2s at the expense of Axl+CD207+ Langerhans cell differentiation from blood cDC2s. BMP7 supplementation allowed to model cDC2 generation and their further differentiation into LCs from CD34+ hematopoietic progenitor cells in defined serum-free medium. Additionally, p38MAPK promoted the generation of another cDC2 subset lacking Axl but expressing the non-classical NFkB transcription factor RelB in vitro. Such RelB+cDC2s occurred predominantly at dermal sites in the inflamed skin. Finally, we found that cDC2s can be induced to acquire high levels of the monocyte lineage identity factor kruppel-like-factor-4 (KLF4) along with monocyte-derived DC and macrophage phenotypic characteristics in vitro. In conclusion, inflammatory and psoriatic epidermal signals instruct blood cDC2s to acquire phenotypic characteristics of several tissue-resident cell subsets.

Cux1+ proliferative basal cells promote epidermal hyperplasia in chronic dry skin disease identified by single-cell RNA transcriptomics

Pathological dry skin is a disturbing and intractable healthcare burden, characterized by epithelial hyperplasia and severe itch. Atopic dermatitis (AD) and psoriasis models with complications of dry skin have been studied using single-cell RNA sequencing (scRNA-seq). However, scRNA-seq analysis of the dry skin mouse model (acetone/ether/water (AEW)-treated model) is still lacking. Here, we used scRNA-seq and in situ hybridization to identify a novel proliferative basal cell (PBC) state that exclusively expresses transcription factor CUT-like homeobox 1 (Cux1). Further in vitro study demonstrated that Cux1 is vital for keratinocyte proliferation by regulating a series of cyclin-dependent kinases (CDKs) and cyclins. Clinically, Cux1+ PBCs were increased in patients with psoriasis, suggesting that Cux1+ PBCs play an important part in epidermal hyperplasia. This study presents a systematic knowledge of the transcriptomic changes in a chronic dry skin mouse model, as well as a potential therapeutic target against dry skin-related dermatoses.

Single cell and spatial sequencing define processes by which keratinocytes and fibroblasts amplify inflammatory responses in psoriasis

The immunopathogenesis of psoriasis, a common chronic inflammatory disease of the skin, is incompletely understood. Here we demonstrate, using a combination of single cell and spatial RNA sequencing, IL-36 dependent amplification of IL-17A and TNF inflammatory responses in the absence of neutrophil proteases, which primarily occur within the supraspinous layer of the psoriatic epidermis. We further show that a subset of SFRP2+ fibroblasts in psoriasis contribute to amplification of the immune network through transition to a pro-inflammatory state. The SFRP2+ fibroblast communication network involves production of CCL13, CCL19 and CXCL12, connected by ligand-receptor interactions to other spatially proximate cell types: CCR2+ myeloid cells, CCR7+ LAMP3+ dendritic cells, and CXCR4 expressed on both CD8+ Tc17 cells and keratinocytes, respectively. The SFRP2+ fibroblasts also express cathepsin S, further amplifying inflammatory responses by activating IL-36G in keratinocytes. These data provide an in-depth view of psoriasis pathogenesis, which expands our understanding of the critical cellular participants to include inflammatory fibroblasts and their cellular interactions.

Visualizing scRNA-Seq Data at Population Scale with GloScope

Increasingly scRNA-Seq studies explore the heterogeneity of cell populations across different samples and its effect on an organism's phenotype. However, relatively few bioinformatic methods have been developed which adequately address the variation between samples for such population-level analyses. We propose a framework for representing the entire single-cell profile of a sample, which we call its GloScope representation. We implement GloScope on scRNA-Seq datasets from study designs ranging from 12 to over 300 samples. These examples demonstrate how GloScope allows researchers to perform essential bioinformatic tasks at the sample-level, in particular visualization and quality control assessment.

Single-cell transcriptomics of human-skin-equivalent organoids

Several methods for generating human-skin-equivalent (HSE) organoid cultures are in use to study skin biology; however, few studies thoroughly characterize these systems. To fill this gap, we use single-cell transcriptomics to compare in vitro HSEs, xenograft HSEs, and in vivo epidermis. By combining differential gene expression, pseudotime analyses, and spatial localization, we reconstruct HSE keratinocyte differentiation trajectories that recapitulate known in vivo epidermal differentiation pathways and show that HSEs contain major in vivo cellular states. However, HSEs also develop unique keratinocyte states, an expanded basal stem cell program, and disrupted terminal differentiation. Cell-cell communication modeling shows aberrant epithelial-to-mesenchymal transition (EMT)-associated signaling pathways that alter upon epidermal growth factor (EGF) supplementation. Last, xenograft HSEs at early time points post transplantation significantly rescue many in vitro deficits while undergoing a hypoxic response that drives an alternative differentiation lineage. This study highlights the strengths and limitations of organoid cultures and identifies areas for potential innovation.

Decoding the Human Epidermal Complexity at Single-Cell Resolution

The epidermis is one of the largest tissues in the human body, serving as a protective barrier. The basal layer of the epidermis, which consists of epithelial stem cells and transient amplifying progenitors, represents its proliferative compartment. As keratinocytes migrate from the basal layer to the skin surface, they exit the cell cycle and initiate terminal differentiation, ultimately generating the suprabasal epidermal layers. A deeper understanding of the molecular mechanisms and pathways driving keratinocytes' organization and regeneration is essential for successful therapeutic approaches. Single-cell techniques are valuable tools for studying molecular heterogeneity. The high-resolution characterization obtained with these technologies has identified disease-specific drivers and new therapeutic targets, further promoting the advancement of personalized therapies. This review summarizes the latest findings on the transcriptomic and epigenetic profiling of human epidermal cells, analyzed from human biopsy or after in vitro cultivation, focusing on physiological, wound healing, and inflammatory skin conditions.

Distinct mechanisms for sebaceous gland self-renewal and regeneration provide durability in response to injury

Sebaceous glands (SGs) release oils that protect our skin, but how these glands respond to injury has not been previously examined. Here, we report that SGs are largely self-renewed by dedicated stem cell pools during homeostasis. Using targeted single cell RNA-sequencing, we uncovered both direct and indirect paths by which these resident SG progenitors ordinarily differentiate into sebocytes, including transit through a PPARγ+Krt5+ transitional cell state. Upon skin injury, however, SG progenitors depart their niche, reepithelialize the wound, and are replaced by hair follicle-derived stem cells. Furthermore, following targeted genetic ablation of >99% of SGs from dorsal skin, these glands unexpectedly regenerate within weeks. This regenerative process is mediated by alternative stem cells originating from the hair follicle bulge, is dependent upon FGFR signaling, and can be accelerated by inducing hair growth. Altogether, our studies demonstrate that stem cell plasticity promotes SG durability following injury.

Generation of Schwann cell derived melanocytes from hPSCs identifies pro-metastatic factors in melanoma

The neural crest (NC) is highly multipotent and generates diverse lineages in the developing embryo. However, spatiotemporally distinct NC populations display differences in fate potential, such as increased gliogenic and parasympathetic potential from later migrating, nerve-associated Schwann cell precursors (SCPs). Interestingly, while melanogenic potential is shared by both early migrating NC and SCPs, differences in melanocyte identity resulting from differentiation through these temporally distinct progenitors have not been determined. Here, we leverage a human pluripotent stem cell (hPSC) model of NC temporal patterning to comprehensively characterize human NC heterogeneity, fate bias, and lineage development. We captured the transition of NC differentiation between temporally and transcriptionally distinct melanogenic progenitors and identified modules of candidate transcription factor and signaling activity associated with this transition. For the first time, we established a protocol for the directed differentiation of melanocytes from hPSCs through a SCP intermediate, termed trajectory 2 (T2) melanocytes. Leveraging an existing protocol for differentiating early NC-derived melanocytes, termed trajectory 1 (T1), we performed the first comprehensive comparison of transcriptional and functional differences between these distinct melanocyte populations, revealing differences in pigmentation and unique expression of transcription factors, ligands, receptors and surface markers. We found a significant link between the T2 melanocyte transcriptional signature and decreased survival in melanoma patients in the cancer genome atlas (TCGA). We performed an in vivo CRISPRi screen of T1 and T2 melanocyte signature genes in a human melanoma cell line and discovered several T2-specific markers that promote lung metastasis in mice. We further demonstrated that one of these factors, SNRPB, regulates the splicing of transcripts involved in metastasis relevant functions such as migration, cell adhesion and proliferation. Overall, this study identifies distinct developmental trajectories as a source of diversity in melanocytes and implicates the unique molecular signature of SCP-derived melanocytes in metastatic melanoma.

Transcriptomic Heterogeneity of Skin Across Different Anatomic Sites

Multiomic studies, including RNA sequencing, single-cell RNA sequencing, and epigenomics, can provide insight into the connection between anatomically heterogeneous gene expression profile of the skin and dermatoses-predisposed sites, in which RNA sequencing is essential. Therefore, in this study, 159 skin samples collected mainly from discarded normal skin tissue during surgical treatment for benign skin tumors were used for RNA sequencing. On the basis of cluster analysis, the skin was divided into four regions, with each region showing specific physiological characteristics through differentially expressed gene analysis. The results showed that the head and neck region, perineum, and palmoplantar area were closely associated with lipid metabolism, hormone metabolism, blood circulation, and related neural regulation, respectively. Transcription factor enrichment indicated that different regions were associated with the development of adjacent tissues. Specifically, the head and neck region, trunk and extremities, perineum, and palmoplantar area were associated with the central nervous, axial, urogenital, and vascular systems, respectively. The results were imported into an open website (https://dermvis.github.io/) for retrieval. Our transcriptomic data elucidated that human skin exhibits transcriptomic heterogeneity reflecting physiological and developmental variation at different anatomic sites and provided guidance for further studies on skin development and dermatoses predisposed sites.