**Methods:** The authors performed single-cell RNA sequencing (scRNA-seq) on 138,982 cells from 13 samples of six patients, including three AK, one SCCIS, three cSCC (low-risk and high-risk), and six patient-matched normal skin samples. Single-cell ATAC sequencing (scATAC-seq) was also performed on one poorly differentiated cSCC sample. Cell types were identified using known markers, and copy number variation (CNV) analysis was used to define malignant cells. Differentially expressed genes (DEGs) were identified across stages, and key candidates were validated by immunohistochemistry (IHC) and immunofluorescence (IF) in independent cohorts (20 AK vs. normal pairs; 15 normal, 15 SCCIS, 60 cSCC). Functional experiments (siRNA knockdown, CCK8, apoptosis, migration, invasion assays) were performed in HaCaT, A431, SCL-I, and SCL-II cell lines. Cell-cell communication was analyzed using CellChat.
**Key Results:**
- Nine major cell types were identified in normal skin, with subtypes in basal (Basal1, Basal2), spinous (Spinous1, Spinous2), and follicular (Follicular1, Follicular2) cells. Basal1 expressed stemness markers (COL17A1, TP63, ITGB1), while Basal2 showed increased inflammatory genes (KRT16, S100A8, S100A9).
- From normal to AK, 549, 305, and 434 up-regulated DEGs were found in Basal1, Basal2, and Pro KCs, respectively. ALDH3A1 and IGFBP2 were significantly up-regulated in AK (IF: 85% and 75% of AK tissues positive, respectively; p<0.001) but not in SCCIS or cSCC. IGFBP2 overexpression promoted proliferation and invasion in HaCaT and A431 cells (p<0.05).
- In the P2 patient (AK and SCCIS), 21 monotonically up-regulated genes were identified in basal cells from normal to AK to SCCIS, including KLF6, FOSL1, JUND, SQSTM1, and RHOB. Monotonically down-regulated genes included SPRR2A, S100A14, and ENO1.
- In SCCIS, basal cell proportion increased significantly. CNV analysis identified a malignant basal subpopulation (Basal-SCCIS-tumor) with higher CNV scores, enriched for DNA damage response genes (PCNA, MCM7), heat shock proteins (HSPA1A/B, HSP90AA1, HSPA6), activated keratins (KRT6A/B/C, KRT16, KRT17, KRT19), and S100 family genes (S100A7, S100A8, S100A9). MAGEA4 was strongly positive in 73.33% of SCCIS and 76.67% of cSCC (p<0.001 vs. normal). ITGA6 was up-regulated in SCCIS (p<0.05) and cSCC (p<0.001). siRNA silencing of MAGEA4 in A431 cells reduced proliferation, migration, and invasion (p<0.01) and increased apoptosis (p<0.01). ITGA6 silencing similarly reduced proliferation, migration, and invasion in A431, SCL-I, and SCL-II cells (p<0.01) and increased apoptosis (p<0.01).
- In cSCC, CNV scores correlated with malignant grade (well-differentiated: low; moderately differentiated: moderate; poorly differentiated: high). LGALS1, IFITM3, and FTH1 were significantly overexpressed in cSCC vs. normal (IHC: p<0.001 for all; IFITM3 and FTH1 also higher in poorly differentiated vs. well-differentiated, p<0.05). siRNA silencing of LGALS1, IFITM3, FTH1, BST2, and SAT1 inhibited proliferation, increased apoptosis (p<0.01), and reduced migration (except FTH1) and invasion (p<0.01) in cSCC cell lines.
- TME analysis revealed 11 lymphocyte subpopulations, 8 DC/Mac subpopulations, and 6 fibroblast subpopulations. Cell-cell communication strength increased with malignancy (p<0.0001). MHC-II, laminin, and TNF signaling pathways were enhanced in poorly differentiated cSCC.
**Clinical Implications:** This study provides a comprehensive single-cell atlas of cSCC progression, identifying early malignant cells in SCCIS and key driver genes (ALDH3A1, IGFBP2, MAGEA4, ITGA6, LGALS1, IFITM3, FTH1, BST2, SAT1) that may serve as biomarkers or therapeutic targets. The findings highlight the importance of basal cell subpopulations and TME interactions in cSCC development, offering potential avenues for early detection and targeted therapy.