Gene expression landscape of cutaneous squamous cell carcinoma progression

Background Cutaneous squamous cell carcinomas (cSCCs) are the second most common human cancer and have been characterized by RNA sequencing (RNA-Seq); however, the transferability of findings from individual studies may be limited by small sample sizes and diverse analysis protocols. Objectives To define the transcriptome landscape at different stages in the progression of normal skin to cSCC via a meta-analysis of publicly available RNA-Seq samples. Methods Whole-transcriptome data from 73 clinically normal skin samples, 46 actinic keratoses (AK) samples, 16 in situ SCC samples, 13 keratoacanthoma (KA) samples and 147 cSCC samples [including 30 samples from immunocompromised patients and 8 from individuals with recessive dystrophic epidermolysis bullosa (RDEB)] were uniformly processed to harmonize gene expression. Differential expression, fusion detection and cell-type deconvolution analyses were performed. Results Individual RNA-Seq studies of cSCC demonstrated study-specific clustering and varied widely in their differential gene expression detection. Following batch correction, we defined a consensus set of differentially expressed genes (DEGs), including those altered in the preinvasive stages of cSCC development, and used single-cell RNA-Seq data to demonstrate that DEGs are often - but not always - expressed by tumour-specific keratinocytes (TSKs). Analysis of the cellular composition of cSCC, KA and RDEB-cSCC identified an increase in differentiated keratinocytes in KA, while RDEB-cSCC contained the most TSKs. Compared with cSCC arising in immunocompetent individuals, cSCC samples from immunosuppressed patients demonstrated fewer memory B cells and CD8+ T cells. A comprehensive and unbiased search for fusion transcripts in cSCC and intermediate disease stages identified few candidates that recurred in >1% of all specimens, suggesting that most cSCC are not driven by oncogenic gene fusions. Finally, using Genotype-Tissue Expression (GTEx) data, we distilled a novel 300-gene signature of chronic sun exposure that affirms greater cumulative ultraviolet (UV) exposure in later stages of cSCC development. Conclusions Our results define the gene expression landscape of cSCC progression, characterize cell subpopulation heterogeneity in cSCC subtypes that contribute to their distinct clinical phenotypes, demonstrate that gene fusions are not a common cause of cSCC and identify UV-responsive genes associated with cSCC development.