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

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.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.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.The findings of this study might help better understand AK pathogenesis.