Single-cell analysis of psoriasis resolution reveals an inflammatory fibroblast state targeted by IL-23 blockade

Luc Francis, Daniel McCluskey, Clarisse Ganier, Treasa Jiang, Xinyi Du-Harpur, Jeyrroy Gabriel, Pawan Dhami, Yogesh Kamra, Sudha Visvanathan, Jonathan N. Barker, Catherine H. Smith, Francesca Capon, Satveer K. Mahil
doi: https://doi.org/10.1101/2023.09.11.23295356

Abstract

Biologics targeting the IL-23/IL-17 axis have transformed the treatment of psoriasis. However, the early mechanisms of action of these drugs remain poorly understood. Here, we performed longitudinal single-cell RNA-sequencing in affected individuals receiving IL-23 inhibitor therapy. By profiling skin at baseline, day 3 and day 14 of treatment, we demonstrated that IL-23 blockade causes marked gene expression shifts, with fibroblast and myeloid populations displaying the most extensive changes at day 3. We also identified a transient WNT5A+/IL24+ fibroblast state, which was only detectable in lesional skin. In-silico and in-vitro studies indicated that signals stemming from these WNT5A+/IL24+ fibroblasts upregulated multiple inflammatory genes in keratinocytes. Importantly, the abundance of WNT5A+/IL24+ fibroblasts was significantly reduced after treatment. This observation was validated in-silico, by deconvolution of multiple transcriptomic datasets, and experimentally, by RNA in-situ hybridization. These findings demonstrate that the evolution of inflammatory fibroblast states is a key feature of resolving psoriasis skin.

INTRODUCTION

Skin is an archetypal barrier tissue which protects from external insults, while orchestrating innate and adaptive responses to microbial agents 1. Our understanding of these functions has been transformed by advances in single-cell transcriptomics, which have uncovered key differences between the cellular ecosystems of healthy and inflamed skin 2, 3, 4, 5. Profiling of skin lesions, however, provides only a snapshot of end-stage disease, with limited insight into the sequence of events leading to the onset and resolution of pathology. In this context, a longitudinal analysis of resolving disease can illuminate early mechanisms of drug action, identifying the fundamental events that promote a return to homeostasis. While this potential is being realized in oncology (e.g., by studying evolving cell states during immunotherapy) 6, 7, inflammatory conditions remain relatively underexplored 8, 9.

Psoriasis is a chronic, immune-mediated skin disorder affecting up to 2% of individuals of European descent. It presents with red, disfiguring skin plaques and is associated with multi-morbidity, poor quality of life and increased cardiovascular mortality 10. Excessive activation of the IL-23/IL-17 axis is a key disease driver. IL-23 produced by dendritic cells and macrophages promotes the differentiation of T17 cells, leading to the release of IL-17. This cytokine, which acts in synergy with TNF, causes keratinocytes to proliferate and produce chemokines that attract further T17 cells to sites of inflammation. The ensuing feedback loop sustains chronic inflammation and epidermal hyperplasia, two hallmarks of psoriasis 11.

The IL-23/IL-17 immune axis has been successfully targeted by powerful biologics inhibiting IL-17 or IL-23. IL-23 antagonists have proven particularly efficacious, with rates of skin clearance (remission) exceeding 60% after 1 year 12. However, the early events leading to drug-induced remission are poorly understood. Skin biopsies have been profiled during trials of IL-17 and IL-23 antagonists 13, 14, but these studies have been limited by the poor resolution of bulk RNA-sequencing and the use of samples obtained weeks after the onset of clinical benefit.

Here, we sought to elucidate the early events underlying the effects of IL-23 antagonists, with a view to uncovering key processes regulating skin homeostasis. We performed longitudinal single-cell RNA-sequencing (scRNA-seq) of psoriasis skin, in individuals receiving the IL-23 inhibitor risankizumab. We generated a high-resolution atlas of early disease resolution, identifying a pro-inflammatory WNT5A+/IL24+ fibroblast state that is enriched in psoriasis lesions. This population, which shows a marked upregulation of TNF and IL-17 signaling, starts to decline in abundance 3 days after treatment initiation. Importantly, an early reduction in WNT5A+/IL24+ fibroblasts was also detectable in patients receiving other systemic and topical therapeutics, suggesting that the normalization of this cell state is critical to the resolution of psoriatic inflammation.

RESULTS

Single-cell RNA sequencing identifies key cell populations responding to IL-23 blockade

To investigate the early effects of IL-23 inhibition, we obtained serial, full-thickness skin biopsies from five individuals with psoriasis who had achieved skin clearance (remission) by 16 weeks of risankizumab treatment (Supplementary Table 1, Fig. 1, A and B). We sampled lesional skin at baseline (day 0, when we also obtained non-lesional biopsies) and at two early time points (day 3 and day 14 of treatment) before full resolution of inflammation was clinically apparent. Following scRNA-seq, we generated expression profiles from 164,553 viable cells (Supplementary Table 2). We analyzed this dataset with Seurat 15, identifying cell clusters corresponding to the main cell populations found in the epidermis (keratinocytes and melanocytes), dermis (fibroblasts, vascular and lymphatic endothelial cells, pericytes, epithelial cells) and in the psoriasis immune infiltrate (mast cells, myeloid cells and T cells) (Fig. 1, C and D; Supplementary Fig. 1).

Fig. 1:
Fig. 1: Identification of cell populations showing an early response to IL-23 blockade.

(A) Experimental design of the study, showing timepoints for the sampling of non-lesional (day 0) and lesional (day 0, day 3, day 14) skin biopsies. (B) Left: Percentage reduction of disease severity (measured by the psoriasis area and severity index, PASI) after 3 and 14 days of risankizumab treatment. Right: representative skin images showing clinical improvement during risankizumab treatment. (C) Unifold Manifold Approximation and Projection (UMAP) of 164,553 single cells (lesional and non-lesional skin), visualised as clusters (left) or skin cell types (right). (D) Dot plot showing the expression of marker genes used for the annotation of cell identities. (E) UMAP visualization of lesional skin showing the number of differentially expressed genes (DEG) observed in each cluster, after 3 (left) and 14 (right) days of treatment. (F) Heatmap depicting the number of ligand-receptor interactions between cell types. The cell types expressing the ligand are labelled on the y-axis, whereas those expressing the receptor are on the x-axis. KC, keratinocytes; FB, fibroblasts; PE, pericytes; LEC, lymphatic endothelial cells; VEC, vascular endothelial cells; EpC, epithelial cells; MEL, melanocytes; MAST, mast cells; MYEL, myeloid cells; TC, T cells.

Differential expression analysis (day 3 vs baseline) showed that the majority of day 3 changes were observed in myeloid cells (25 differential expressed genes, DEG) and fibroblasts (25 DEG), identifying these populations as early responders to IL-23 inhibition (Fig. 1E). Ligand-receptor analyses undertaken with CellChat 16 demonstrated a widespread reduction in the number of cell-cell interactions after 3 days of treatment (Fig. 1F). In keeping with the results of the differential expression analysis, the number of interactions originating from fibroblasts and myeloid cells was markedly reduced.

By day 14, myeloid cells and keratinocytes were displaying the strongest gene expression changes (121 and 115 DEG, respectively) (Fig. 1E). The ligand-receptor analyses were broadly consistent with these observations, showing that the number of interactions driven by myeloid cells was further reduced after 14 days of treatment (Fig. 1F).

Thus, myeloid cells and fibroblasts are the cell types showing the earliest response to IL-23 inhibition, with keratinocyte changes occurring shortly after. Importantly, the day 3 cellular shifts were observed before treatment effects became apparent. This mirrors previous observations of gene expression changes preceding clinical remission (weeks 2-4 of treatment), reported in psoriasis patients receiving other biologics17, 18.

IL-23 blockade affects cell proliferation and IL-17 signaling in keratinocytes

To investigate the effects of IL-23 inhibition in further detail, we sought to define population subsets for each cluster. In the epidermis, melanocytes could not be further subdivided, but keratinocytes could be classified in five distinct sub-clusters (Fig. 2A and Supplementary Fig. 2A). These corresponded to basal (KRT15+), follicular (KRT17+), proliferating (MKI67+), spinous (KRT1+/SPRR2E-) and supra-spinous (KRT1+/SPRR2E+) populations (Fig. 2B). In keeping with the early clinical improvement of skin lesions (Fig. 1B), the abundance of proliferating keratinocytes was significantly reduced by day 14 (Fig. 2C). The largest number of DEG was observed in supra-spinous keratinocytes (n=228 at day 14 vs day 0) and spinous keratinocytes (n=193 at day 14 vs day 0) (Fig. 2D). The downregulated genes detected at day 14 in these populations included inflammatory molecules (e.g., IL36G, S100A7/A8/A9) and structural proteins (GJB2, KRT16/17) (Supplementary Fig. 2C), known to be over-expressed in the upper layers of the psoriasis epidermis5, 19. Conversely, the analysis of basal keratinocytes showed the upregulation of CCL27 and KRT15, two genes that are lowly expressed in lesional skin5, 19 (Supplementary Fig. 2C).

Fig. 2:
Fig. 2: Effects of IL-23 inhibition in keratinocytes.

(A) UMAP of 68,695 keratinocytes forming five distinct subclusters. (B) Dot plot showing the expression of marker genes used for the annotation of cell identities. (C) Abundance of proliferating cells (expressed as keratinocyte percentage) at different time points after risankizumab treatment. *P<0.05, **P<0.01 (repeated measures one-way ANOVA with Dunnett’s post-test). (D) UMAP visualization of keratinocytes from lesional skin, showing the number of differentially expressed genes (DEG) observed in each cluster, after 3 (left) and 14 (right) days of treatment. KC, keratinocytes; L, lesional skin; NL, non-lesional skin.

Pathway enrichment analysis of the DEG identified a significant over-representation of IL-17 related genes (“Role of IL-17A in psoriasis”; FDR <10-4 in supra-spinous keratinocytes, FDR<10-2 in spinous keratinocytes), confirming the suppressive effect of IL-23 inhibition on T17 activation. Likewise, an upstream regulator analysis demonstrated a significant enrichment for genes that are induced by IL-17 and its synergistic partner TNF (z< -2.50 for both in supra-spinous keratinocytes, FDR<10-5). Conversely, the impact of risankizumab on IFN-γ (another inflammatory cytokine induced downstream of IL-23) was limited, as the activation score (z= -1.97) did not reach statistical significance (z< -2.0). Thus, the main effect of IL-23 blockade in the epidermis is a marked downregulation of IL-17 signaling in spinous and supra-spinous keratinocytes.

IL-23 inhibition downregulates cytotoxic T cell activity and increases dendritic cells (DC)-3 signaling

Our initial clustering identified three immune cell types: mast cells, T cells and dendritic cells (DC). While mast cells could be subdivided into three sub-clusters (TPSAB1-/IL1RL1+, TPSAB1+/IL1RL1+, TPSAB1+/IL1RL1-) (Supplementary Fig. 3, A and B), the frequencies of these populations did not change with treatment (Supplementary Fig. 3C). The number of DEG observed in each cluster was also low (<15 at day 14), suggesting that IL-23 inhibition has minimal effects on mast-cell activity.

T cells formed six distinct populations, including T helper (CD40LG+), Treg (CTLA4+/FOXP3+), proliferating Treg (MKI67+/CTLA4+/FOXP3+), terminally differentiated effector memory (TEMRA) T cells (CD8+/CCL5+), cytotoxic T17 cells (CD8+/IL17A+) and activated PLCG2+ cells. A CD3-/XCL1+ cluster was also identified, corresponding to natural killer/innate lymphoid cells (NK/ILC) (Fig. 3, A and B, and Supplementary Fig. 3, D and E). Differential expression analysis identified TEMRA T cells as the population showing the strongest response to IL-23 blockade, with 28 DEG detected at day 3 and 30 at day 14 (Fig. 3C). These included several genes related to effector function (e.g., GZMA and GZMB at day 3; NKG7, KLRB1 and IL7R at day 14), suggesting a reduction in cytotoxic activity (Fig. 3D).

Fig. 3:
Fig. 3: Effects of IL-23 inhibition on lymphocytes and dendritic cells.

(A) UMAP of 6,263 lymphocytes forming seven distinct clusters. (B) Dot plot showing the expression of marker genes used for the annotation of T cell identities. (C) UMAP visualization of lymphocytes from lesional skin, showing the number of differentially expressed genes (DEG) observed in each cluster, after 3 (left) and 14 (right) days of treatment. (D) Downregulation of cytotoxicity markers in TEMRA T cells, following treatment with risankizumab. (E) UMAP of 18,544 myeloid cells, forming 11 distinct clusters. (F) Dot plot showing the expression of marker genes used for the annotation of myeloid cell identities. (G) UMAP visualization of myeloid cells from lesional skin, showing the number of differentially expressed genes (DEG) observed in each cluster, after 3 (left) and 14 (right) days of treatment. L, lesional skin.

The analysis of the myeloid compartment identified 10 LYZ+ sub-clusters. These broadly recapitulated the main mononuclear phagocyte populations described by Reynolds et al in psoriatic skin 5. In fact, we observed cDC1 (CLEC9A+) and proliferating cDC1 (CLEC9A+/MKI67+) cells, two cDC2 populations (CLEC10A+/CD1C+), two DC3 populations (IL1B+/IL23A+), mregDC (LAMP3+/BIRC3+), Langerhans cells (CLEC4A+/CD207+) and two macrophage populations (CD68+/CD163+). A LYZ-sub-cluster was also detected and identified as a plasmacytoid DC (pDC) population, based on IRF7 expression (Fig. 3, E and F, and Supplementary Fig. 3, F and G). Of note, the expression profile of myeloid subsets was in keeping with that described by Nakamizo et al in their single-cell dissection of psoriasis skin20. mregDC were the main source of IL15 and IL32, while the two DC3 subsets expressed IL23A, IL1B, CXCL8/IL8, IL10 and the glucose transporter SLC2A3. We were also able to reproduce the expansion of DC3, mregDC and cDC1 cells in lesional skin. Our analysis showed that the abundance of these cell subsets is not affected during the early stages of IL-23 blockade (Supplementary Fig. 3H). Conversely, day 14 vs day 0 differential expression analysis identified the two DC3 subsets as the cells showing the strongest response to risankizumab (180 DEG detected in DC3-2 and 155 in DC3-1) (Fig. 3G). DEG included IL23A, CXCL3 and IL1B, all of which were upregulated at day 14. Notably, the expression of TNFAIP3 (a negative regulator of IL-17 and NF-κB signaling) was also increased (Supplementary Fig. 3I), suggesting that the elevation in cytokine gene expression is a transient phenomenon, which could be caused by the modulation of negative feedback loops regulating IL-17/IL-23 signaling 21.

IL-23 blockade downregulates mediators of leukocyte diapedesis in pericytes and endothelial cells

While fibroblasts (described below) were the most abundant cell type in the dermis (Supplementary Fig 1B), we also observed epithelial cells, lymphatic endothelial cells, vascular endothelial cells and pericytes. Neither epithelial cells nor lymphatic endothelial cells could be further subdivided, while vascular endothelial cells formed five sub-clusters (Fig. 4A, Supplementary Fig. 4, A and B). These included arterioles (HEY1+/S100A4+), capillaries (PLVAP+/RGCC+; ACKR1+/FABP4+), post-capillary venules (ACKR1+/SELE+) and venules (ICAM1+/VWF+) (Fig. 4B). Although the number of DEG observed at day 14 was <40 in all populations, there was evidence for a significant downregulation of IL-17 signaling in lymphatic endothelial cells and ACKR1+/FABP4+ capillaries FDR<10-7 in both).

Fig. 4:
Fig. 4: Effects of IL-23 inhibition on vascular endothelial cells and pericytes.

(A) UMAP of 16,420 vascular endothelial cells forming five clusters. (B) Dot plot showing the expression of marker genes used for the annotation of vascular endothelial cell identities. (C) Risankizumab-induced downregulation of proinflammatory molecules and collagen genes, in PLVAP+/RGCC+ capillaries. (D) UMAP of 13,507 pericytes forming two clusters. (E) Dot plot showing the expression of marker genes used for the annotation of pericyte cell identities. (F) UMAP visualization of vascular endothelial cells from lesional skin, showing the number of differentially expressed genes (DEG) observed in each cluster, after 3 (left) and 14 (right) days of treatment. (G) Genes promoting leukocyte extravasation are downregulated in pericyte 1, following risankizumab treatment. L, lesional skin; NL, non-lesional skin.

Downregulation of genes promoting leukocyte chemotaxis (S100A7, S100A8, S100A9) and extravasation (COL4A1, COL4A2) was also observed in PLVAP+/RGCC+ capillaries (Fig. 4C).

Pericytes formed two clusters (pericyte 1 and 2), distinguishable based on PDGFRB and MYL9 expression (Fig. 4, D and E, and Supplementary Fig. 4, C and D). Of note, pericytes 1 expressed IL17A, similarly to the PDGRFB+ pericyte-like fibroblasts described by Gao et al. Pericyte 1 cells also showed the largest number of DEG (n=63, day 14 vs day 0) (Fig. 4F), including several genes implicated in leukocyte chemotaxis (CXCL8, S100A7, S100A8, S100A9) and diapedesis (VIM, ACTA2, SPARC) (Fig. 4G).

Thus, IL-23 blockade causes an early downregulation of the intermediate filaments, extracellular matrix proteins and chemokines that facilitate the migration of leukocytes towards sites of skin inflammation.

IL-23 inhibition reduces the abundance of a pro-inflammatory, WNT5A+/IL24+ fibroblast population

Fibroblasts formed 11 distinct sub-clusters (Fig. 5A and Supplementary Fig. 5, A and B). These included CCN5+/PI16+, COL11A1+/DPEP1+ and ANGPTL7+ populations, as well as three COMP+ sub-clusters, three APOE+ populations and a WNT5A+/IL24+ sub-cluster (Fig. 5B). While CCN5+/PI16+ cells showed a strong response to treatment at day 3 (143 DEG), the number of differentially expressed genes at day 14 was lower (<60) in all subclusters (Supplementary Fig. 5C). In keeping with this observation, we were able to reproduce the previously reported upregulation of CCL19 in psoriasis fibroblasts5. However, we noted that CCL19 expression was not significantly reduced by IL-23 blockade, suggesting that the modulation of this chemokine is not essential for early-stage remission (Supplementary Fig. 5D).

Fig. 5:
Fig. 5: Effects of IL-23 inhibition on fibroblasts.

(A) UMAP of 31,765 fibroblasts forming eleven clusters. (B) Dot plot showing the expression of marker genes used for the annotation of cell identities. (C) Plot showing the decline in WNT5A+/IL24+ cell abundance, following treatment with risankizumab. Every line represents a patient. *P<0.05; **P<0.01 (repeated measures ANOVA, with Dunnett’s post-test). (D) Real time PCR analysis showing that key markers for WNT5A+/IL24+ cells are upregulated in primary human fibroblasts stimulated with IL-17A and TNF. N=5; *P<0.05; **P<0.01, ***P<0.001 (Friedman test with Dunn’s post-test). (E) Real time PCR analysis of inflammatory markers in human primary keratinocytes cultured with unconditioned medium (no FB grown), supernatants from IL-17A/TNF stimulated fibroblasts (Stim), control medium (supernatant from unstimulated fibroblasts, Unstim) or control medium supplemented with IL-24. N=2-4; *P<0.05; **P<0.01 (Kruskal Wallis test with Dunn’s post-test). (F) Plot showing the SCENIC regulons detected in WNT5A+/IL24+ fibroblasts. The three regulons with the highest specificity scores are indicated by labels, with their size (number of regulated genes) in brackets. (G) Boxplots showing fibroblast clusters ordered by pseudotime. (H) Dot plot showing the expression of the top regulons in the various fibroblast clusters. L, lesional skin; NL, non-lesional skin.

Conversely, we found that the treatment had a marked impact on cell abundance, causing a significant decrease in the frequency of WNT5A+/IL24+ cells. This phenomenon mirrors the absence of WNT5A+/IL24+ cells from non-lesional skin (Fig. 5C), identifying the sub-cluster as a lesional cell state. Accordingly, an analysis of the 237 genes that were differentially expressed in WNT5A+/IL24+ cells compared to all other fibroblast clusters demonstrated a marked enrichment for the TNF and IL-17 signaling pathways (FDR<10-10 and 10-4, respectively). This was reflected by the expression of TNFRSF1A, TNFRSF1B, IL17RA and IL17RC (but not IL23R) in WNT5A+/IL24+ fibroblasts (Supplementary Fig. 6A).

To further validate these findings, we stimulated primary human fibroblasts with IL-17A and TNF. We found that the treatment caused a significant upregulation of key marker genes that are preferentially expressed in the WNT5A+/IL24+ population (WNT5A, IL24, CXCL8/IL8 and TBX3) (Fig. 5D). Thus, IL-17 synergizes with TNF in shaping the expression profiles of WNT5A+/IL24+ cells.

In keeping with the specificity of WNT5A+/IL24+ cells, IL-24 signaling was only detectable in day 0 lesional skin, where CellChat ligand-receptor analysis uncovered an interaction between the cytokine (originating from WNT5A+/IL24+ fibroblasts) and its receptor on spinous keratinocytes (Supplementary Fig. 5E). The effects of risankizumab on IL-24 signaling were independently validated with NicheNet 22, which showed that reduced IL24 levels at day 14 correlated with the differential expression of IL-24 target genes in spinous keratinocytes (Supplementary Fig. 5F). Of note, the DEG included several transcription factors regulating keratinocyte proliferation and differentiation (ID1, CEBPB, STAT1), as well as inflammatory mediators (S100A8, S100A9, SERPINB3) and structural proteins (KRT15). These observations suggest that IL-24 produced by WNT5A+/IL24+ fibroblasts promote the activation of spinous keratinocytes.

To further investigate these findings, we re-analyzed spatial transcriptomic profiles generated in lesional and non-lesional psoriasis skin. By examining WNT5A and COMP as markers for our population of interest (IL24 expression was undetectable), we confirmed that it was located in the upper dermis and expanded in lesional skin (Supplementary Fig. 6, B and C). We then identified the genes that were preferentially expressed in WNT5A+/COMP+ spots. We found that these included several markers of spinous and supra-spinous keratinocytes (CALML5, KRT1, KRT10, SLURP1, SPRR2E), as well as pro-inflammatory mediators detected in the in-silico ligand-receptor analysis (S100A8, S100A9, SERPINB3) (Supplementary Fig. 6D). These findings indicate that WNT5A+/IL24+ fibroblasts are located in an inflammatory microenvironment, in close proximity to spinous/supra-spinous keratinocytes.

To recapitulate the ligand-receptor interactions inferred in-silico, we cultured IL-17A/TNF stimulated primary human fibroblasts as an in-vitro model for the WNT5A+/IL24+ state. We then transferred supernatants from these cells to human primary keratinocytes. While the supernatants did not affect the expression of genes related to keratinocyte proliferation (MKI67, STAT1) or differentiation (ID1, CEBPB) (Supplementary Fig. 6E), they caused the upregulation of key inflammatory mediators uncovered by the in-silico ligand-receptor analysis (S100A8, S100A9 and SERPINB3) (Fig. 5E). Of note, these changes in gene expression were not fully recapitulated by IL-24 supplementation of control medium (supernatant from unstimulated primary fibroblasts) (Fig. 5E), indicating a likely synergy with other cytokines.

To further characterize WNT5A+/IL24+ fibroblasts, we carried out an analysis of gene regulatory networks with SCENIC 23. This demonstrated an enrichment for FOXQ1, TBX3 and PRDM1 regulons in WNT5A+/IL24+ cells (Fig. 5F). Further investigation of fibroblast trajectories using Monocle 3 24 identified the three COMP+ subsets as the populations mapping closest to WNT5A+/IL24+ fibroblasts (Fig. 5G). This is in line with the progressively decreasing expression of WNT5A (Fig. 5B) and FOXQ1-regulated genes in COMP+3, COMP+2 and COMP+1 cells (Fig. 5H) and their absence from all other clusters. Given the role of FOXQ1 in activating Wnt signaling 25, this suggests that the drug-induced reduction in WNT5A+/IL24+ cell numbers reflects a progressive downregulation of Wnt-related genes, leading to the re-establishment of a physiological COMP+ phenotype. Accordingly, we found that IL-23 inhibitor treatment caused a significant increase in the abundance of COMP+ fibroblasts (Supplementary Fig. 5G).

A WNT5A+/IL24+ fibroblast state is also detectable in prurigo nodularis

To understand whether WNT5A+/IL24+ fibroblasts are present in other inflammatory skin diseases beyond psoriasis, we sought to detect these cells in publicly available scRNA-seq datasets. We therefore retrieved single-cell profiles generated in healthy skin, prurigo nodularis and atopic dermatitis26, and integrated them with our day 0 data (Supplementary Fig. 7, A-C). We found that WNT5A+/IL24+ signature genes were absent from healthy skin, weakly expressed in atopic dermatitis, but readily detectable in prurigo nodularis lesions (Supplementary Fig. 7, D). Interestingly, a closer inspection of single-cell profiles identified some differences between prurigo nodularis and psoriasis lesions, with WNT5A+/IL24+ fibroblasts expressing TBX3 and CXCL8/IL8 only in the latter (Supplementary Fig. 7E). Taken together, our observations identify WNT5A+/IL24+ fibroblasts as an IL-17/TNF dependent, inflammatory cell state, that in psoriasis, is rapidly normalized by IL-23 blockade.

A reduction in WNT5A+/IL24+ cell numbers is an early event in the response to psoriasis therapeutics

To validate the pathogenic involvement of WNT5A+/IL24+ fibroblasts, we sought to detect these cells in other psoriasis datasets generated in the first two weeks of treatment. We retrieved three transcriptomic datasets from public repositories and used a deconvolution algorithm (CIBERSORTx 27) to estimate the frequency of WNT5A+/IL24+ fibroblasts in patient skin. We then compared the inferred cell abundance in different sample groups.

The first dataset we examined was generated in individuals with psoriasis who were receiving guselkumab, a biologic therapy that blocks the same IL-23 subunit targeted by risankizumab 28. The analysis of these samples confirmed that WNT5A+/IL24+ fibroblasts are only detectable in lesional psoriasis skin, where their abundance begins to decrease after one week of treatment (Fig. 6A). The second dataset was obtained in a placebo-controlled clinical trial of the IL-17 inhibitor ixekizumab 29. Again, our analysis showed that the decline of WNT5A+/IL24+ fibroblasts occurred early (2 weeks) after treatment. Importantly, this reduction in cell numbers was not observed among the individuals receiving placebo (Fig. 6B). Finally, we examined RNA-seq data generated in individuals receiving a topical glucocorticoid (halometasone monohydrate 0.05% cream) for the treatment of psoriasis 30. We observed that the abundance of WNT5A+/IL24+ fibroblasts in lesional skin was significantly reduced after two weeks of effective treatment. Of note, this trend was reversed upon treatment withdrawal and clinical relapse (Fig. 6C).

Fig. 6:
Fig. 6: The abundance of WNT5A+/IL24+ fibroblasts declines rapidly following psoriasis treatment.

(A-C) Deconvolution of transcriptome datasets shows that the frequency of WNT5A+/IL24+ fibroblasts is significantly reduced by treatment with guselkumab (A) (N=8-22 per group), ixekizumab (B) (N=6-7 per group) or topical glucocorticoid (halomethasone monohydrate 0.05%) (N=7 per group), with the latter effect going into reverse upon drug withdrawal and clinical relapse (C). (D) The percentage of WNT5A+ fibroblasts (PDGFRA+ cells) detected by RNA in-situ hybridization decreases following treatment with risankizumab, N=3. (E) Representative confocal microscopy image showing the results of RNA in-situ hybridization (10x magnification, scale bar 100 µm). Arrows indicate cells triple-positive for PDGFRA (green), COMP (red) and WNT5A (magenta). (F) Diagram showing the early effects of therapeutic IL-23 inhibition (schematic created using BioRender.com). Left: in lesional skin, IL-23 drives the differentiation of T17 cells. These produce IL-17, leading to the activation of keratinocytes and the evolution of COMP+ fibroblasts into a WNT5A+/IL24+ state. The ensuing IL-24 production further reinforces keratinocyte activation. Right: following IL-23 blockade, T17 cell differentiation and IL-17 production are inhibited, while WNT5A+/IL24+ fibroblasts revert to their COMP+ state. As a result of diminished IL-17 and IL-24 production, keratinocyte homeostasis is re-established. Ixe, ixekizumab; L, lesional skin; NL, non-lesional skin; ns, not significant; *P<0.05; **P<0.01; ***P<0.0001 (mixed effects ANOVA with Dunnett’s post-test (A), paired t-test (B) or repeated measures ANOVA with Dunnet’s post-test (C-D)).

To complement these in-silico observations, we carried out RNA in-situ hybridization in the skin of three newly ascertained patients, who had been successfully treated with risankizumab. We sampled lesional psoriasis skin at baseline (day 0, when non-lesional biopsies were also obtained) and at day 14 of treatment. The analysis of these samples confirmed that our population of interest was present in the upper dermis of lesional psoriasis skin, where its abundance was markedly reduced by IL-23 blockade (Fig. 6, D and E).

Taken together, these observations identify the reduction of WNT5A+/IL24+ fibroblasts as an early event mediating the resolution of skin inflammation in psoriasis, following systemic or topical treatment (Fig. 6F).

DISCUSSION

We present a comprehensive single-cell atlas of resolving psoriasis, which captures the dynamic cell states mediating a return to skin homeostasis. We specifically investigated the effects of risankizumab, as this highly effective drug targets the central IL-23/T17 pathogenic axis 12.

We readily detected a marked downregulation of IL-17 signaling in keratinocytes, fibroblasts and vascular endothelial cells sampled after 14 days of risankizumab treatment. At this time point, the impact of IL-23 blockade on IFN-γ activity was modest, with inhibition scores only trending towards statistical significance in keratinocytes. These results complement those of a recent study where the systemic effects of IL-23 inhibition were examined by phenotyping individuals with IL23R deficiency. Interestingly, the analysis of this rare condition uncovered prominent defects of IFN-γ responses and only a limited impairment of IL-17 activity31. Thus, the consequences of IL-23 inhibition may be organ specific. This warrants the study of IL-23 antagonists in other tissues targeted by these drugs (e.g., the intestinal mucosa 32 and the synovium 33).

A key result from our study is the identification of a WNT5A+/IL24+ inflammatory fibroblast state, which is expanded in lesional psoriasis skin and rapidly shrinks after treatment initiation. Interestingly, other cell types that are altered in psoriasis (e.g., DC3, mregDC and cDC1 cells) were similarly expanded in lesional skin, but did not decrease in abundance within 2 weeks of treatment.

Our interrogation of external datasets also indicates that the decline in WNT5A+/IL24+ cell numbers is a drug-agnostic phenomenon, suggesting that the normalization of this cell state might contribute to early, treatment induced remission.

Inflammatory fibroblasts are typified by their ability to produce and respond to cytokines. As single-cell transcriptomics have illuminated the functional diversity of stromal cells, these inflammatory subsets have been implicated in the pathogenesis of several immune mediated conditions 3, 34, 35. For example, an integrated analysis of four immune-mediated disorders (rheumatoid arthritis, inflammatory bowel disease, interstitial lung disease, and Sjogren’s syndrome) identified two shared pro-inflammatory states (CXCL10+/CCL19+ and SPARC+/COL3A1+ fibroblasts) that were also detectable in atopic dermatitis 36.

Here we observed both cell states in our samples (corresponding to APOE+/CCL19+ and COL11A1+ fibroblasts) and detected an upregulation of CCL19 in lesional fibroblasts. However, neither inflammatory population showed increased abundance in lesional skin or declined in frequency following IL-23 inhibitor treatment (Supplementary Fig. 5B). A WNT5A+/POSTN+ fibroblast subset has recently been identified in the inflammatory skin disease prurigo nodularis 37. In affected individuals, these cells were shifted towards a cancer-associated phenotype and showed a prominent activation of type I interferon responses. Neither of these features were observed in psoriasis WNT5A+/IL24+ cells, despite these expressing POSTN. Instead, our analysis indicated that a WNT5A+/IL24+ inflammatory state can be accompanied by the upregulation of disease-specific genes such as CXCL8/IL8 in psoriasis. As additional datasets enabling robust sub-clustering of fibroblasts become available, it will be possible to determine whether this phenomenon also applies to other IL-17/TNF mediated skin disorders. Likewise, the analysis of future datasets of treatment responses will establish if WNT5A+/IL24+ fibroblasts undergo dynamic, treatment-dependent shifts in other inflammatory skin diseases beyond psoriasis.

Both WNT5A and IL24 are overexpressed in psoriasis skin lesions 38, 39. Furthermore, the upregulation of IL-24 has specifically been detected in human dermal fibroblasts exposed to IL-17 cytokines 39, which is in keeping with our observation of enhanced IL-17 signaling in WNT5A+/IL24+ cells. Notably, prior research has demonstrated that both Wnt5a and IL-24 promote keratinocyte proliferation and cytokine release 38, 39. Accordingly, our ligand-receptor analyses showed that WNT5A+/IL24+ fibroblasts communicate with spinous keratinocytes in lesional skin. We also found that IL-23 inhibition abolishes this interaction, causing the downregulation of genes mediating keratinocyte activation. Our functional data supports an inflammatory crosstalk between WNT5A+/IL24+ fibroblasts and keratinocytes, although further experiments will be required to identify the ligands mediating cell-cell communication and fully characterize their pathogenic potential. The crosstalk was further substantiated by our in-situ hybridization experiments and re-analysis of spatial transcriptomic data, which mapped WNT5A/IL24+ fibroblasts to the upper dermis, in proximity to keratinocytes.

Our work also sheds new light on the relationships between fibroblast populations, suggesting that the WNT5A+/IL24+ state originates from COMP+ cells. As dermal fibroblast plasticity has been implicated in wound healing 40, fibrosis 41 and tumor responses 42, these results add to emerging evidence that the evolution of dynamic fibroblast states is also important for the onset and resolution of chronic inflammation.

Our study has some limitations. We were unable to obtain viable neutrophils and only captured a relatively small number of T cells (∼6,000 in total), which prevented us from resolving rarer subsets (e.g., MAIT cells) and cell states. Similar issues have been reported in other scRNA-seq studies 43 and likely reflect cell damage caused by the skin dissociation protocol. Given our stringent inclusion criteria, our cohort only included male individuals of European descent without comorbidities. Thus, further investigations will be required to establish whether our findings can be generalized to other demographics.

Since the aim of this study was to investigate the early effects of IL-23 inhibition, our patient sample collection deliberately targeted the first two weeks of treatment. Nonetheless, our analysis of publicly available transcriptome data indicates that the abundance of WNT5A+/IL24+ cells rapidly increases upon treatment withdrawal and clinical relapse. If this finding is confirmed in future cohorts, the frequency of WNT5A+/IL24+ cells could be monitored in the increasing number of patients receiving biologics, with the potential to inform personalized, drug tapering strategies in psoriasis remission.

METHODS

Participants

Written informed consent was obtained from all participants, in line with approval from the London - Westminster Research Ethics Committee (REC ref 11/H0802/7). All participants were biologic-naïve, adult males of European descent, with a dermatologist confirmed diagnosis of severe psoriasis (PASI > 10) (Supplementary Table 1). To minimize confounders, we also applied these additional inclusion criteria: no comorbidities, no systemic immune-modifying therapy for at least 4 weeks prior to risankizumab treatment, no topical steroids to the sampled area for at least 5 days prior to sampling. Six to eight mm full thickness punch biopsies were obtained from participants at baseline (day 0; lesional and non-lesional skin) and after 3 and 14 days of treatment (lesional skin only). All biopsies were taken from site-matched areas on the lower back/buttocks.

Whole skin dissociation, scRNA-seq library preparation and scRNA-sequencing

Skin biopsies were processed immediately after collection and incubated in 5U/ml dispase (Stem Cell Technologies) at 4°C overnight, to separate the epidermis from the dermis. The dermis was minced and then processed using the Whole Skin Dissociation Kit (Miltenyi Biotec) and a gentleMACS™ Tissue Dissociator (Miltenyi Biotec). The epidermis was minced and digested with a 50:50 trypsin-EDTA: Versene solution (Gibco) at 37°C for 15 minutes. Dissociated cells were then filtered through a 100μm cell strainer. Following incubation with Fc block (Biolegend), TotalSeq™ Hashtag (Biolegend) antibodies were added to delineate lesional dermis, lesional epidermis, non-lesional dermis and non-lesional epidermis cell suspensions. Following sorting of viable cells on a FACSAria™ flow cytometer (BD), dermal and epidermal suspensions were mixed at a 1:1 ratio. Libraries were prepared with the Chromium Single Cell 3’ kit (10X Genomics) and sequenced using a NextSeq 2000 instrument (Illumina) to generate 150 base pair paired end reads.

scRNA-seq data processing

FASTQ files were aligned against the human GRCh38 reference genome using Cell Ranger v4.0.0 (10X Genomics). Cells expressing <300 or >4000 unique genes, or with greater than 20% of genes of mitochondrial origin, were removed. The scDblFinder package (v3.1.6) 44 was used to examine each sample and cells identified as doublets were also excluded. To avoid batch effects, data generated in different NextSeq runs was integrated using the reciprocal principal component analysis workflow provided by the Seurat package (v4.1.0) 15. The data were normalized and then scaled by regressing out the following variables: number of unique molecular identifiers per cell, percentage of mitochondrial genes, day of single-cell dissociation, patient identifier, and difference between G2M and S phase score.

Cell clustering and cell type annotation

The batch-corrected coordinate space was used for linear dimensional reduction with Seurat. A K-Nearest Neighbour graph was constructed with the FindNeighbours function and unsupervised cell clustering was implemented using the FindClusters function at a resolution of 0.4). Subclustering was performed using the same approach, except the resolution was optimised separately for each population, using Clustree 45.

Genes that were differentially expressed (fold change >1.5; Bonferroni adjusted p-value <0.05) in a cluster compared to all the others were identified using the FindAllMarkers function. Cell identities were then annotated by cross referencing these cluster-specific genes with published signature genes for the relevant cell type. The same approach was applied to subclusters, except small subclusters expressing multiple lineage markers were considered as doublets and removed.

Differential expression analysis

Genes that were differentially expressed between day 0 and day 3 or day 0 and day 14 were identified with the FindMarkers function. Genes detected in <5% of cells from either sample group were removed. A negative binomial regression was applied, including patient id as a latent variable.

Identification of enriched pathways, upstream regulators and gene regulatory networks

The Ingenuity Pathway Analysis package (Qiagen) was used to identify canonical pathways and upstream regulators that were enriched (FDR<0.05) among the DEG identified with the FindAllMarkers and FindMarkers functions. DEG networks centred around IL-17, IL-17R and IL-22 were visualised with the igraph (v1.2.6) R Package 46.

The gene regulatory networks that were active in lesional fibroblast subclusters at day 0 were identified using the workflow implemented by SCENIC (v1.3.0) 23. Briefly, co-expression modules between transcription factor and target genes (regulons) were derived using GRNBoost (a scalable alternative to GENIE3), then RcisTarget was used to refine the regulons by inferring direct targets of the transcription factors. Finally, regulon activity scores were calculated for each cell, using AUCell.

Cell-cell interaction inference

Receptor-ligand interaction analysis was performed using CellChat (v1.4.0) 16. The main analysis was based on major cell type annotations. A separate run was performed for each sample group and the number of significant interactions was calculated for each cell type pair. For the analysis of IL-24 signalling, keratinocytes, fibroblasts, T cells, and myeloid cells were divided into subclusters. Ligand-receptor interactions originating from WNT5A+/IL24+ fibroblasts were then queried at day 0 and day 14. The interaction between IL-24 and its receptor was also confirmed by using NicheNet (v1.1.1)22 to analyse day 0 and day 14 lesional skin data. WNT5A+/IL24+ cells were set as source and the various keratinocyte subclusters were in turn set as recipient cells. Results were generated using the ‘nichenet_seruatobj_aggregate’ function.

Cell culture

Primary dermal fibroblasts and primary keratinocytes were isolated from healthy skin discarded after plastic surgery. Fibroblasts were grown in Dulbecco Modified Eagle Medium (DMEM) supplemented with 10% fetal bovine serum and 1% penicillin-streptomycin (all reagents from Gibco). For the analysis of WNT5A+/IL24+ marker genes, fibroblasts were stimulated for 24 hours with 5 ng/mL TNF (R&D Systems), 50 ng/mL IL-17A (R&D Systems), 5 ng/mL TNF + 50 ng/mL IL-17A, or vehicle. Cells were then harvested for RNA isolation.

For fibroblast-keratinocyte crosstalk experiments, fibroblasts were first stimulated for 24 hours with 5 ng/mL TNF + 50 ng/mL IL-17A, or vehicle. Cells were then washed with phosphate buffered saline (Gibco) and incubated in complete DMEM medium for a further 24h. Meanwhile keratinocytes were grown in EpiLife Keratinocyte Medium supplemented with 60μM calcium, Supplement 7 and 1% penicillin-streptomycin (all reagents from Gibco). After an overnight incubation in Supplement 7-free medium, keratinocytes were incubated for 24h in complete Epilife Medium diluted 1:3 in either complete DMEM, culture medium from TNF/IL-17A stimulated fibroblasts, culture medium from unstimulated fibroblasts, or culture medium from unstimulated fibroblasts supplemented with IL-24 (2 ng/ml or 10 ng/ml, R&D Systems). Following RNA isolation from harvested keratinocytes, gene expression was assessed by real-time PCR using primers reported in Supplementary Table 3. Transcript levels were normalized to expression of GAPDH (fibroblasts) or PKG1 (keratinocytes).

Pseudotime analysis

Seurat fibroblast data from day 0 lesional skin were converted and imported into Monocle 3 (v1.3.1) 24 using the cell_data_set function from SeuratWrappers. Cells were separated into clusters by Leiden community detection (cluster_cells function). A cell trajectory was then fitted with the learn_graph function. Root nodes were selected within the WNT5A+/IL24A+ population and the order_cells function was applied.

Re-analysis of publicly available datasets

For the analysis of spatial transcriptomic profiles, the dataset generated by Castillo et al47 (GSE202011) was downloaded as an RDS file from https://zenodo.org/records/7562864. The data were analyzed using Seurat to visualize gene expression. Spots were considered positive for the expression of genes of interest, if the relevant UMI count was ≥1. FindMarkers was then used to identify genes that were upregulated in WNT5A+/COMP+ spots compared to all other spots.

For the analysis of WNT5A+/IL24+ cells in healthy, prurigo nodularis and atopic dermatitis skin, the scRNA-seq datasets generated by Alkon et al48 and Rindler et al49 were retrieved from the Gene Expression Omnibus, using identifier GSE222840. The datasets were then integrated with our day 0 lesional and non-lesional scRNA-seq data. Briefly, all datasets were normalized using the SCTransform function then merged into a single Seurat object. Following principal component analysis of the normalized data, batch correction was performed using the RunHarmony function from the Harmony package50 (v1.0.1), with sample ID and dataset variables. The main cell clusters were identified by DEG annotation after applying FindNeighbors and FindClusters. The fibroblast cluster was extracted, split by sample id and re-integrated by repeating the steps above. Finally, the fibroblasts were sub-clustered at a resolution of 0.3 and the Seurat AddModuleScore function was applied using the expression of 20 signature genes for WNT5A+/IL24+ cells (the top DEG observed in WNT5A+/IL24+ cells vs all other fibroblasts in baseline psoriasis skin).

For the deconvolution analysis, the scRNA-seq dataset produced in this study was down-sampled to 5000 cells and used to generate a reference signature matrix for CIBERSORTx 27. The transcriptomic datasets generated by Sofen et al., 28 Krueger et al., 29 and Cai et al., 30 were retrieved from GEO using their respective identifiers (GSE51440, GSE31652 and GSE114729). The reference matrix was then used to infer the relative abundance of WNT5A+/IL24+ fibroblasts.

RNA fluorescence in situ hybridization

16μm sections were cut from skin biopsies frozen in optimal cutting temperature compound (VWR International). RNA fluorescence in situ hybridization of Hs-PDGFRA, Hs-COMP-C2 and Hs-WNT5A-C4 probes (ACD Biologics) was performed using RNAscope® Multiplex Fluorescent Reagent Kit v2 (ACD Biologics) and 4-plex Ancillary kit (ACD Biologics). The manufacturer’s instructions were followed with minor modifications. Fixation in 4% paraformaldehyde/PBS (Santa Cruz Biotechnology) was extended to 60 minutes at 4°C, and TSA VIVID™ Fluorophore 520, 570 and 650 (ACD Biologics) were diluted at 1:1000. At least two representative images were captured for each section, using an Eclipse Ti-2 inverted microscope (Nikon).

Statistical analysis

All statistical tests were implemented in GraphPad Prism (version 9.3.0 for Windows). For the analysis of scRNA-seq and RNA in situ hybridisation results, the relative abundance of cell types in the four sample groups (non-lesional, lesional day 0, lesional day 3, lesional day 14) was compared using repeated measures analysis of variance (ANOVA), with a Dunnett’s post-test. For the analysis of cell fractions derived by deconvolution, the relative abundance of WNT5A+/IL24+ fibroblasts was compared across treatment groups, using a repeated measures ANOVA with a Dunnet’s post-test or a paired t-test, as appropriate. The real-time PCR data was analyzed using a Kruskal Wallis test or Friedman test, as appropriate. Both were followed by Dunn’s multiple comparison test. Details of all statistical tests are reported in the relevant figure legends. Plots are shown with mean and standard deviation error bars.

Data availability

The scRNA-seq data associated with this manuscript have been deposited in the NCBI Gene Expression Omnibus (series id: GSE228421). The rest of the data needed to evaluate the conclusions of the study are present in the paper or the Supplementary Materials. Source data are also provided with this paper.

Author contributions

F.C. and S.K.M conceived the study, designed the experiments with input from LF, and obtained the necessary funding. Further funding was provided by S.V; S.K.M. supervised the patient recruitment, which was also facilitated by C.H.S and J.N.B; Y.K. prepared the single-cell libraries under P.D.’s supervision. L.F. carried out the remaining experimental work with help from T.J. and input from S.K.M., C.G., J.G. and X.D.H; L.F. and D.Mc. carried out the computational analyses under F.C.’s supervision. F.C. and S.K.M. drafted the manuscript, which was then reviewed by LF with input from all co-authors.

Competing interests

SV is a Boehringer-Ingelheim employee. JNB has attended advisory boards and/or spoken at sponsored symposia and/or received research funding from: AbbVie, Almirall, Amgen, Boehringer-Ingelheim, Bristol Myers Squibb, Celgene, Janssen, Leo, Lilly, Novartis, Samsung, Sun

Pharma. CHS reports departmental research funding as investigator in EU-IMI consortia involving multiple industry partners (see biomap-imi.eu and hippocrates-imi.eu for details). FC has received grant support and consultancy fees from Boehringer Ingelheim. SKM reports departmental income from Abbvie, Almirall, Eli Lilly, Leo, Novartis, Sanofi, UCB, outside the submitted work.

Acknowledgments

We are grateful to all the patients who made this study possible. We also wish to acknowledge the support of King’s College London Advanced Cytometry Platform (Isabel Correa, Richard Ellis, Leanne Farnan, Anna Rose) and Genomics Core (Ulrich Kadolsky, Shichina Kannambath, Michelle Kleeman, Athul Menon, Rosamond Nuamah, Heli Vaikkinen). We thank Paola Di Meglio for her input on the manuscript, Hannah Toghill, Yong-Xin Sieh and Imogen Brooks for technical assistance.

This research was supported by the National Institute for Health and Care Research (NIHR) Biomedical Research Centre based at Guy’s and St Thomas’ NHS Foundation Trust and King’s College London (guysbrc-2012-1). The views expressed are those of the author(s) and not necessarily those of the NHS, the NIHR or the Department of Health and Social Care. This work was funded by the Psoriasis Association (PhD studentship ST2/21 RE18411 supporting LF and Grant BSTOP50/5 to CHS) and the Wellcome Trust (grant 096540/Z/11/Z). DMc and TJ were supported by the Medical Research Council (grants MR/R015643/1 and MR/W006820/1) and King’s College London as a member of the MRC Doctoral Training Partnership in Biomedical Sciences. SKM is funded by a NIHR Advanced Fellowship (NIHR302258). CHS is supported by a NIHR Senior Investigator Award.

Footnotes

  • Text updated to better relate findings to previous psoriasis single-cell papers; Additional experiments performed to functionally validate the IL24+/WNT5A+ fibroblast state; Section added on analysis of IL24+/WNT5A+ fibroblast state in Atopic Dermatitis, Prurigo Nodularis and Healthy Skin; Section added on analysis of IL24+/WNT5A+ fibroblast state in a public spatial transcriptomic dataset; Figure 3, 5 and 6 revised; Supplemental files updated; Author list updated.

References

  1. 1.↵
    Pasparakis M, Haase I, Nestle FO. Mechanisms regulating skin immunity and inflammation. Nat Rev Immunol 2014, 14(5): 289–301.
    OpenUrlCrossRefPubMed
  2. 2.↵
    Billi AC, Ma F, Plazyo O, Gharaee-Kermani M, Wasikowski R, Hile GA, et al. Nonlesional lupus skin contributes to inflammatory education of myeloid cells and primes for cutaneous inflammation. Science translational medicine 2022, 14(642): eabn2263.
    OpenUrl
  3. 3.↵
    Ko KI, Merlet JJ, DerGarabedian BP, Zhen H, Suzuki-Horiuchi Y, Hedberg ML, et al. NF-kappaB perturbation reveals unique immunomodulatory functions in Prx1(+) fibroblasts that promote development of atopic dermatitis. Science translational medicine 2022, 14(630): eabj0324.
    OpenUrlCrossRef
  4. 4.↵
    Liu Y, Wang H, Taylor M, Cook C, Martinez-Berdeja A, North JP, et al. Classification of human chronic inflammatory skin disease based on single-cell immune profiling. Sci Immunol 2022, 7(70): eabl9165.
    OpenUrl
  5. 5.↵
    Reynolds G, Vegh P, Fletcher J, Poyner EFM, Stephenson E, Goh I, et al. Developmental cell programs are co-opted in inflammatory skin disease. Science 2021, 371(6527).
  6. 6.↵
    Bachireddy P, Azizi E, Burdziak C, Nguyen VN, Ennis CS, Maurer K, et al. Mapping the evolution of T cell states during response and resistance to adoptive cellular therapy. Cell reports 2021, 37(6): 109992.
    OpenUrl
  7. 7.↵
    Yost KE, Satpathy AT, Wells DK, Qi Y, Wang C, Kageyama R, et al. Clonal replacement of tumor-specific T cells following PD-1 blockade. Nat Med 2019, 25(8): 1251–1259.
    OpenUrlCrossRefPubMed
  8. 8.↵
    Bangert C, Rindler K, Krausgruber T, Alkon N, Thaler FM, Kurz H, et al. Persistence of mature dendritic cells, T(H)2A, and Tc2 cells characterize clinically resolved atopic dermatitis under IL-4Ralpha blockade. Sci Immunol 2021, 6(55).
  9. 9.↵
    Damsky W, Wang A, Kim DJ, Young BD, Singh K, Murphy MJ, et al. Inhibition of type 1 immunity with tofacitinib is associated with marked improvement in longstanding sarcoidosis. Nature communications 2022, 13(1): 3140.
    OpenUrl
  10. 10.↵
    Burden AD, Kirby B. Psoriasis and related disorders. In: Griffiths CEM, Barker JN, Bleiker T, Chalmers RJ, Creamer D (eds). Rook’s Textbook of Dermatology. Wiley-Blackwell: Chichester, 2016.
  11. 11.↵
    Puig L, Costanzo A, Munoz-Elias EJ, Jazra M, Wegner S, Paul CF, et al. The biological basis of disease recurrence in psoriasis: a historical perspective and current models. Br J Dermatol 2022, 186(5): 773–781.
    OpenUrl
  12. 12.↵
    Ghoreschi K, Balato A, Enerback C, Sabat R. Therapeutics targeting the IL-23 and IL-17 pathway in psoriasis. Lancet 2021, 397(10275): 754–766.
    OpenUrl
  13. 13.↵
    Lu J, Wang Y, Li Y, Zhong X, Gong Y, Ding Y, et al. Based on Gene Expression Analysis: Low-Density Neutrophil Expression Is a Characteristic of the Fast Responders Treated With Guselkumab for Psoriasis. Frontiers in immunology 2022, 13: 865875.
    OpenUrl
  14. 14.↵
    Visvanathan S, Baum P, Vinisko R, Schmid R, Flack M, Lalovic B, et al. Psoriatic skin molecular and histopathologic profiles after treatment with risankizumab versus ustekinumab. The Journal of allergy and clinical immunology 2019, 143(6): 2158–2169.
    OpenUrl
  15. 15.↵
    Stuart T, Butler A, Hoffman P, Hafemeister C, Papalexi E, Mauck WM, 3rd, et al. Comprehensive Integration of Single-Cell Data. Cell 2019, 177(7): 1888–1902 e1821.
    OpenUrl
  16. 16.↵
    Jin S, Guerrero-Juarez CF, Zhang L, Chang I, Ramos R, Kuan CH, et al. Inference and analysis of cell-cell communication using CellChat. Nature communications 2021, 12(1): 1088.
    OpenUrl
  17. 17.↵
    Correa da Rosa J, Kim J, Tian S, Tomalin LE, Krueger JG, Suarez-Farinas M. Shrinking the Psoriasis Assessment Gap: Early Gene-Expression Profiling Accurately Predicts Response to Long-Term Treatment. J Invest Dermatol 2017, 137(2): 305–312.
    OpenUrl
  18. 18.↵
    Liu J, Chang HW, Grewal R, Cummins DD, Bui A, Beck KM, et al. Transcriptomic Profiling of Plaque Psoriasis and Cutaneous T-Cell Subsets during Treatment with Secukinumab. JID Innov 2022, 2(3): 100094.
    OpenUrl
  19. 19.↵
    Kim J, Lee J, Kim HJ, Kameyama N, Nazarian R, Der E, et al. Single-cell transcriptomics applied to emigrating cells from psoriasis elucidate pathogenic versus regulatory immune cell subsets. The Journal of allergy and clinical immunology 2021, 148(5): 1281–1292.
    OpenUrl
  20. 20.↵
    Nakamizo S, Dutertre CA, Khalilnezhad A, Zhang XM, Lim S, Lum J, et al. Single-cell analysis of human skin identifies CD14+ type 3 dendritic cells co-producing IL1B and IL23A in psoriasis. The Journal of experimental medicine 2021, 218(9).
  21. 21.↵
    Liu T, Li S, Ying S, Tang S, Ding Y, Li Y, et al. The IL-23/IL-17 Pathway in Inflammatory Skin Diseases: From Bench to Bedside. Frontiers in immunology 2020, 11: 594735.
    OpenUrl
  22. 22.↵
    Browaeys R, Saelens W, Saeys Y. NicheNet: modeling intercellular communication by linking ligands to target genes. Nature methods 2020, 17(2): 159–162.
    OpenUrl
  23. 23.↵
    Aibar S, Gonzalez-Blas CB, Moerman T, Huynh-Thu VA, Imrichova H, Hulselmans G, et al. SCENIC: single-cell regulatory network inference and clustering. Nature methods 2017, 14(11): 1083–1086.
    OpenUrl
  24. 24.↵
    Cao J, Spielmann M, Qiu X, Huang X, Ibrahim DM, Hill AJ, et al. The single-cell transcriptional landscape of mammalian organogenesis. Nature 2019, 566(7745): 496-502.
    OpenUrlCrossRefPubMed
  25. 25.↵
    Pizzolato G, Moparthi L, Soderholm S, Cantu C, Koch S. The oncogenic transcription factor FOXQ1 is a differential regulator of Wnt target genes. Journal of cell science 2022, 135(19).
  26. 26.↵
    Arnhold J. The Dual Role of Myeloperoxidase in Immune Response. Int J Mol Sci 2020, 21(21).
  27. 27.↵
    Newman AM, Steen CB, Liu CL, Gentles AJ, Chaudhuri AA, Scherer F, et al. Determining cell type abundance and expression from bulk tissues with digital cytometry. Nature biotechnology 2019, 37(7): 773–782.
    OpenUrlPubMed
  28. 28.↵
    Sofen H, Smith S, Matheson RT, Leonardi CL, Calderon C, Brodmerkel C, et al. Guselkumab (an IL-23-specific mAb) demonstrates clinical and molecular response in patients with moderate-to-severe psoriasis. The Journal of allergy and clinical immunology 2014, 133(4): 1032–1040.
    OpenUrlCrossRef
  29. 29.↵
    Krueger JG, Fretzin S, Suarez-Farinas M, Haslett PA, Phipps KM, Cameron GS, et al. IL-17A is essential for cell activation and inflammatory gene circuits in subjects with psoriasis. The Journal of allergy and clinical immunology 2012, 130(1): 145–154 e149.
    OpenUrlCrossRefPubMedWeb of Science
  30. 30.↵
    Cai Y, Xue F, Quan C, Qu M, Liu N, Zhang Y, et al. A Critical Role of the IL-1beta-IL-1R Signaling Pathway in Skin Inflammation and Psoriasis Pathogenesis. J Invest Dermatol 2019, 139(1): 146–156.
    OpenUrlCrossRef
  31. 31.↵
    Philippot Q, Ogishi M, Bohlen J, Puchan J, Arias AA, Nguyen T, et al. Human IL-23 is essential for IFN-gamma-dependent immunity to mycobacteria. Sci Immunol 2023, 8(80): eabq5204.
    OpenUrl
  32. 32.↵
    Sandborn WJ, D’Haens GR, Reinisch W, Panes J, Chan D, Gonzalez S, et al. Guselkumab for the Treatment of Crohn’s Disease: Induction Results From the Phase 2 GALAXI-1 Study. Gastroenterology 2022, 162(6): 1650–1664 e1658.
    OpenUrlPubMed
  33. 33.↵
    Huang X, Shentu H, He Y, Lai H, Xu C, Chen M, et al. Efficacy and safety of IL-23 inhibitors in the treatment of psoriatic arthritis: a meta-analysis based on randomized controlled trials. Immunol Res 2023.
  34. 34.↵
    Croft AP, Campos J, Jansen K, Turner JD, Marshall J, Attar M, et al. Distinct fibroblast subsets drive inflammation and damage in arthritis. Nature 2019, 570(7760): 246-251.
    OpenUrlCrossRefPubMed
  35. 35.↵
    Smillie CS, Biton M, Ordovas-Montanes J, Sullivan KM, Burgin G, Graham DB, et al. Intra- and Inter-cellular Rewiring of the Human Colon during Ulcerative Colitis. Cell 2019, 178(3): 714–730 e722.
    OpenUrlCrossRefPubMed
  36. 36.↵
    Korsunsky I, Wei K, Pohin M, Kim EY, Barone F, Major T, et al. Cross-tissue, single-cell stromal atlas identifies shared pathological fibroblast phenotypes in four chronic inflammatory diseases. Med 2022, 3(7): 481-518 e414.
    OpenUrl
  37. 37.↵
    Patel JR, Joel MZ, Lee KK, Kambala A, Cornman H, Oladipo O, et al. Single-cell RNA sequencing reveals dysregulated fibroblast subclusters in prurigo nodularis. bioRxiv 2023.
  38. 38.↵
    Tian F, Mauro TM, Li Z. The pathological role of Wnt5a in psoriasis and psoriatic arthritis. J Cell Mol Med 2019, 23(9): 5876–5883.
    OpenUrl
  39. 39.↵
    Xu X, Prens E, Florencia E, Leenen P, Boon L, Asmawidjaja P, et al. Interleukin-17A Drives IL-19 and IL-24 Expression in Skin Stromal Cells Regulating Keratinocyte Proliferation. Frontiers in immunology 2021, 12: 719562.
    OpenUrl
  40. 40.↵
    Abbasi S, Sinha S, Labit E, Rosin NL, Yoon G, Rahmani W, et al. Distinct Regulatory Programs Control the Latent Regenerative Potential of Dermal Fibroblasts during Wound Healing. Cell Stem Cell 2020, 27(3): 396–412 e396.
    OpenUrlCrossRef
  41. 41.↵
    Shaw TJ, Rognoni E. Dissecting Fibroblast Heterogeneity in Health and Fibrotic Disease. Current rheumatology reports 2020, 22(8): 33.
    OpenUrlCrossRef
  42. 42.↵
    Luo H, Xia X, Huang LB, An H, Cao M, Kim GD, et al. Pan-cancer single-cell analysis reveals the heterogeneity and plasticity of cancer-associated fibroblasts in the tumor microenvironment. Nature communications 2022, 13(1): 6619.
    OpenUrl
  43. 43.↵
    He H, Suryawanshi H, Morozov P, Gay-Mimbrera J, Del Duca E, Kim HJ, et al. Single-cell transcriptome analysis of human skin identifies novel fibroblast subpopulation and enrichment of immune subsets in atopic dermatitis. The Journal of allergy and clinical immunology 2020, 145(6): 1615–1628.
    OpenUrl
  44. 44.↵
    Germain PL, Lun A, Garcia Meixide C, Macnair W, Robinson MD. Doublet identification in single-cell sequencing data using scDblFinder. F1000Res 2021, 10: 979.
    OpenUrl
  45. 45.↵
    Zappia L, Oshlack A. Clustering trees: a visualization for evaluating clusterings at multiple resolutions. GigaScience 2018, 7(7).
  46. 46.↵
    Csardi G, T. N. The igraph software package for complex network research. Inter Journal 2006, Complex Systems: 1695.
  47. 47.↵
    Castillo RL, Sidhu I, Dolgalev I, Chu T, Prystupa A, Subudhi I, et al. Spatial transcriptomics stratifies psoriatic disease severity by emergent cellular ecosystems. Sci Immunol 2023, 8(84): eabq7991.
    OpenUrl
  48. 48.↵
    Alkon N, Assen FP, Arnoldner T, Bauer WM, Medjimorec MA, Shaw LE, et al. Single-cell RNA sequencing defines disease-specific differences between chronic nodular prurigo and atopic dermatitis. The Journal of allergy and clinical immunology 2023, 152(2): 420–435.
    OpenUrl
  49. 49.↵
    Rindler K, Jonak C, Alkon N, Thaler FM, Kurz H, Shaw LE, et al. Single-cell RNA sequencing reveals markers of disease progression in primary cutaneous T-cell lymphoma. Molecular cancer 2021, 20(1): 124.
    OpenUrl
  50. 50.↵
    Korsunsky I, Millard N, Fan J, Slowikowski K, Zhang F, Wei K, et al. Fast, sensitive and accurate integration of single-cell data with Harmony. Nature methods 2019, 16(12): 1289–1296.
    OpenUrl
Back to top
Posted November 20, 2023.
Download PDF

Supplementary Material

Data/Code
Email
Single-cell analysis of psoriasis resolution reveals an inflammatory fibroblast state targeted by IL-23 blockade
Luc Francis, Daniel McCluskey, Clarisse Ganier, Treasa Jiang, Xinyi Du-Harpur, Jeyrroy Gabriel, Pawan Dhami, Yogesh Kamra, Sudha Visvanathan, Jonathan N. Barker, Catherine H. Smith, Francesca Capon, Satveer K. Mahil
medRxiv 2023.09.11.23295356; doi: https://doi.org/10.1101/2023.09.11.23295356
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area