Spatial transcriptomics landscape of lesions from non-communicable inflammatory skin diseases

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.

Dupilumab-associated head and neck dermatitis shows a pronounced type 22 immune signature mediated by oligoclonally expanded T cells

Dupilumab, an IL4R-blocking antibody, has shown clinical efficacy for atopic dermatitis (AD) treatment. In addition to conjunctivitis/blepharitis, the de novo appearance of head/neck dermatitis is now recognized as a distinct side effect, occurring in up to 10% of patients. Histopathological features distinct from AD suggest a drug effect, but exact underlying mechanisms remain unknown. We profiled punch biopsies from dupilumab-associated head and neck dermatitis (DAHND) by using single-cell RNA sequencing and compared data with untreated AD and healthy control skin. We show that dupilumab treatment was accompanied by normalization of IL-4/IL-13 downstream activity markers such as CCL13, CCL17, CCL18 and CCL26. By contrast, we found strong increases in type 22-associated markers (IL22, AHR) especially in oligoclonally expanded T cells, accompanied by enhanced keratinocyte activation and IL-22 receptor upregulation. Taken together, we demonstrate that dupilumab effectively dampens conventional type 2 inflammation in DAHND lesions, with concomitant hyperactivation of IL22-associated responses.

Spatial transcriptomics reveals altered lipid metabolism and inflammation-related gene expression of sebaceous glands in psoriasis and atopic dermatitis

Sebaceous glands drive acne, however, their role in other inflammatory skin diseases remains unclear. To shed light on their potential contribution to disease development, we investigated the spatial transcriptome of sebaceous glands in psoriasis and atopic dermatitis patients across lesional and non-lesional human skin samples. Both atopic dermatitis and psoriasis sebaceous glands expressed genes encoding key proteins for lipid metabolism and transport such as ALOX15B, APOC1, FABP7, FADS1/2, FASN, PPARG, and RARRES1. Also, inflammation-related SAA1 was identified as a common spatially variable gene. In atopic dermatitis, genes mainly related to lipid metabolism (e.g. ACAD8, FADS6, or EBP) as well as disease-specific genes, i.e., Th2 inflammation-related lipid-regulating HSD3B1 were differentially expressed. On the contrary, in psoriasis, more inflammation-related spatially variable genes (e.g. SERPINF1, FKBP5, IFIT1/3, DDX58) were identified. Other psoriasis-specific enriched pathways included lipid metabolism (e.g. ACOT4, S1PR3), keratinization (e.g. LCE5A, KRT5/7/16), neutrophil degranulation, and antimicrobial peptides (e.g. LTF, DEFB4A, S100A7-9). In conclusion, our results show that sebaceous glands contribute to skin homeostasis with a cell type-specific lipid metabolism, which is influenced by the inflammatory microenvironment. These findings further support that sebaceous glands are not bystanders in inflammatory skin diseases, but can actively and differentially modulate inflammation in a disease-specific manner.

A single-cell atlas of IL-23 inhibition in cutaneous psoriasis distinguishes clinical response

Psoriasis vulgaris and other chronic inflammatory diseases improve markedly with therapeutic blockade of interleukin-23 (IL-23) signaling, but the genetic mechanisms underlying clinical responses remain poorly understood. Using single-cell transcriptomics, we profiled immune cells isolated from lesional psoriatic skin before and during IL-23 blockade. In clinically responsive patients, a psoriatic transcriptional signature in skin-resident memory T cells was strongly attenuated. In contrast, poorly responsive patients were distinguished by persistent activation of IL-17-producing T (T17) cells, a mechanism distinct from alternative cytokine signaling or resistance isolated to epidermal keratinocytes. Even in IL-23 blockade-responsive patients, we detected a recurring set of recalcitrant, disease-specific transcriptional abnormalities. This irreversible immunological state may necessitate ongoing IL-23 inhibition. Spatial transcriptomic analyses also suggested that successful IL-23 blockade requires dampening of >90% of IL-17-induced response in lymphocyte-adjacent keratinocytes, an unexpectedly high threshold. Collectively, our data establish a patient-level paradigm for dissecting responses to immunomodulatory treatments.

Involvement of transcribed lncRNA uc.291 in hyperproliferative skin disorders

The uc.291 transcript controls keratinocytes differentiation by physical interaction with ACTL6A and subsequent induction of transcription of the genes belonging to the epidermal differentiation complex (EDC). Uc.291 is also implicated in the dedifferentiation phenotype seen in poorly differentiated cutaneous squamous cell carcinomas. Here, we would like to investigate the contribution of uc.291 to the unbalanced differentiation state of keratinocytes observed in hyperproliferative skin disorders, e. g., psoriasis. Psoriasis is a multifactorial inflammatory disease, caused by alteration of keratinocytes homeostasis. The imbalanced differentiation state, triggered by the infiltration of immune cells, represents one of the events responsible for this pathology. In the present work, we explore the role of uc.291 and its interactor ACTL6A in psoriasis skin, using quantitative real-time PCR (RT-qPCR), immunohistochemistry and bioinformatic analysis of publicly available datasets. Our data suggest that the expression of the uc.291 and of EDC genes loricrin and filaggrin (LOR, FLG) is reduced in lesional skin compared to nonlesional skin of psoriatic patients; conversely, the mRNA and protein level of ACTL6A are up-regulated. Furthermore, we provide evidence that the expression of uc.291, FLG and LOR is reduced, while ACTL6A mRNA is up-regulated, in an in vitro psoriasis-like model obtained by treating differentiated keratinocytes with interleukin 22 (IL-22). Furthermore, analysis of a publicly available dataset of human epidermal keratinocytes treated with IL-22 (GSE7216) confirmed our in vitro results. Taken together, our data reveal a novel role of uc.291 and its functional axis with ACTL6A in psoriasis disorder and a proof of concept that biological inhibition of this molecular axis could have a potential pharmacological effect against psoriasis and, in general, in skin diseases with a suppressed differentiation programme.

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.

The Use of Single-Cell RNA-Sequencing and Spatial Transcriptomics in Understanding the Pathogenesis and Treatment of Skin Diseases

The development of multiomic profiling tools has rapidly expanded in recent years, along with their use in profiling skin tissues in various contexts, including dermatologic diseases. Among these tools, single-cell RNA-sequencing (scRNA-seq) and spatial transcriptomics (ST) have emerged as widely adopted and powerful assays for elucidating key cellular components and their spatial arrangement within skin disease. In this paper, we review the recent biological insights gained from the use of scRNA-seq and ST and the advantages of combining both for profiling skin diseases, including aberrant wound healing, inflammatory skin diseases, and cancer. We discuss the role of scRNA-seq and ST in improving skin disease treatments and moving toward the goal of achieving precision medicine in dermatology, whereby patients can be optimally matched to treatments that maximize therapeutic response.

Differentiation of IL-26+ TH17 intermediates into IL-17A producers via epithelial crosstalk in psoriasis

Interleukin (IL)-26 is a TH17 cytokine with known antimicrobial and pro-inflammatory functions. However, the precise role of IL-26 in the context of pathogenic TH17 responses is unknown. Here we identify a population of blood TH17 intermediates that produce high levels of IL-26 and differentiate into IL-17A-producing TH17 cells upon TGF-β1 exposure. By combining single cell RNA sequencing, TCR sequencing and spatial transcriptomics we show that this process occurs in psoriatic skin. In fact, IL-26+ TH17 intermediates infiltrating psoriatic skin induce TGF-β1 expression in basal keratinocytes and thereby promote their own differentiation into IL-17A-producing cells. Thus, our study identifies IL-26-producing cells as an early differentiation stage of TH17 cells that infiltrates psoriatic skin and controls its own maturation into IL17A-producing TH17 cells, via epithelial crosstalk involving paracrine production of TGF-β1.

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.

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.

Spatial transcriptomics in human skin research

Spatial transcriptomics is a revolutionary technique that enables researchers to characterise tissue architecture and localisation of gene expression. A plethora of technologies that map gene expression are currently being developed, aiming to facilitate spatially resolved, high-dimensional assessment of gene transcription in the context of human skin research. Knowing which gene is expressed by which cell and in which location within skin, facilitates understanding of skin function and dysfunction in both health and disease. In this review, we summarise the available spatial transcriptomic methods and we describe their application to a broad spectrum of dermatological diseases.

Spatial cancer systems biology resolves heterotypic interactions and identifies disruption of spatial hierarchy as a pathological driver event

Spatially annotated single-cell datasets provide unprecedented opportunities to dissect cell-cell communication in development and disease. Heterotypic signaling includes interactions between different cell types and is well established in tissue development and spatial organization. Epithelial organization requires several different programs that are tightly regulated. Planar cell polarity (PCP) is the organization of epithelial cells along the planar axis, orthogonal to the apical-basal axis. Here, we investigate PCP factors and explore the implications of developmental regulators as malignant drivers. Utilizing cancer systems biology analysis, we derive a gene expression network for WNT-ligands (WNT) and their cognate frizzled (FZD) receptors in skin cutaneous melanoma. The profiles supported by unsupervised clustering of multiple-sequence alignments identify ligand-independent signaling and implications for metastatic progression based on the underpinning developmental spatial program. Omics studies and spatial biology connect developmental programs with oncological events and explain key spatial features of metastatic aggressiveness. Dysregulation of prominent PCP factors such as specific representatives of the WNT and FZD families in malignant melanoma recapitulates the development program of normal melanocytes but in an uncontrolled and disorganized fashion.

Spatial cancer systems biology resolves heterotypic interactions and identifies disruption of spatial hierarchy as a pathological driver event

Spatially annotated single-cell datasets provide unprecedented opportunities to dissect cell-cell communication in development and disease. Heterotypic signaling includes interactions between different cell types and is well established in tissue development and spatial organization. Epithelial organization requires several different programs that are tightly regulated. Planar cell polarity is the organization of epithelial cells along the planar axis orthogonal to the apical-basal axis. In this study, we investigate planar cell polarity factors and explore the implications of developmental regulators as malignant drivers. Utilizing cancer systems biology analysis, we derive gene expression network for WNT-ligands (WNT) and their cognate frizzled (FZD) receptors in skin cutaneous melanoma. The profiles supported by unsupervised clustering of multiple-sequence alignments identify ligand-independent signaling and implications for metastatic progression based on the underpinning developmental spatial program. Omics studies and spatial biology connect developmental programs with oncological events and explain key spatial features of metastatic aggressiveness. Dysregulation of prominent planar cell polarity factors such specific representative of the WNT and FZD families in malignant melanoma recapitulates the development program of normal melanocytes but in an uncontrolled and disorganized fashion.

Single-Cell Genomics for Investigating Pathogenesis of Inflammatory Diseases

Recent technical advances have enabled unbiased transcriptomic and epigenetic analysis of each cell, known as "single-cell analysis". Single-cell analysis has a variety of technical approaches to investigate the state of each cell, including mRNA levels (transcriptome), the immune repertoire (immune repertoire analysis), cell surface proteins (surface proteome analysis), chromatin accessibility (epigenome), and accordance with genome variants (eQTLs; expression quantitative trait loci). As an effective tool for investigating robust immune responses in coronavirus disease 2019 (COVID-19), many researchers performed single-cell analysis to capture the diverse, unbiased immune cell activation and differentiation. Despite challenges elucidating the complicated immune microenvironments of chronic inflammatory diseases using existing experimental methods, it is now possible to capture the simultaneous immune features of different cell types across inflamed tissues using various single-cell tools. In this review, we introduce patient-based and experimental mouse model research utilizing single-cell analyses in the field of chronic inflammatory diseases, as well as multi-organ atlas targeting immune cells.