Abstract
About 5% of patients with cutaneous squamous cell carcinoma (cSCC) have a poor prognosis which is associated with a loss of tumor differentiation, invasion and metastasis, all of which are linked to the process of epithelial-to-mesenchymal plasticity (EMP). Here, we showed that the EMP-associated transcription factor ZEB2 drives cSCC heterogeneity which resembles biphasic carcinosarcoma-like tumors. Single cell RNA sequencing revealed distinct subpopulations ranging from fully epithelial (E) to intermediate (EM) to fully mesenchymal (M), associated with the gradual loss of cell surface markers EPCAM, CDH1, ITGB4, and CD200. Mesenchymal features were associated with a higher metastatic capacity and anoikis resistance, yet this comes with a sensitivity toward TNF-induced cell death. Altogether we provide insights in cSCC heterogeneity and modes to target mesenchymal-metastasis inducing cells.