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.