Abstract
Cutaneous squamous cell carcinomas (cSCCs) arise from keratinocytes in the skin, but the molecular changes driving this transformation remain unclear. To better understand this process, we perform multi-omic profiling of keratinocytes, actinic keratoses, and cSCCs. Single-cell mutational analyses reveal that most keratinocytes have remarkably low mutation burdens; however, keratinocytes with TP53 or NOTCH1 mutations exhibit substantially higher burdens. These findings suggest that keratinocytes can withstand high dosages of cumulative ultraviolet radiation, but certain pathogenic mutations break these adaptive mechanisms, inducing a mutator phenotype. Mutational profiling of cSCCs adjacent to actinic keratoses reveals TERT promoter and CDKN2A mutations emerge in actinic keratoses, whereas additional mutations that inactivate ARID2 and activate the mitogen-activated protein kinase pathway delineate the transition to cSCC. Surprisingly, actinic keratoses are often not related to their neighboring cSCC. Spatial analyses reveal gene expression heterogeneity, including checkpoint molecule enrichment at invasive fronts, highlighting tumor and immune cell interactions.