Single-cell Sequencing Highlights Heterogeneity and Malignant Progression in Actinic Keratosis and Cutaneous Squamous Cell Carcinoma
Summary:
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-seq (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 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 (DEG) analysis screened out multiple key driver genes including transcription factors (TFs) 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.
Overall Design:
single-cell RNA-seq profiles of 13 samples from 6 patients presenting for surgical resection using the 10X Genomics Chromium platform. All these samples included 3 AK samples, 1 SCCIS tumor samples, 3 cSCC tumor samples (low-risk and high-risk) and patient-matched 6 normal skin samples without any treatment
Cell Clustering
Color UMAP by
Switch between fine-grained and coarse annotations
Loading UMAP…
Enrichment Analysis
Enrichment mode
Select enrichment collection
Top pathways
Ranked by absolute NES
Differentially Expressed Genes
p-value threshold0.05
Maximum adjusted p-value
Log fold change1.0
Minimum log₂ fold change
Cell type group
Filter by cell type
Gene
logFC
p-value
Score
Group
CellPhoneDB Cell-Cell Communication Analysis
Dynamic Analysis
Select Cell Types
Choose 2 or more cell types to analyze ligand-receptor interactions