Abstract:Abstract 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.

Ultradeep sequencing differentiates patterns of skin clonal mutations associated with sun-exposure status and skin cancer burden

Ultra-deep Sequencing Differentiates Patterns of Skin Clonal Mutations Associated with Sun-Exposure Status and Skin Cancer Burden

Large mitochondrial DNA deletions in ultraviolet B-induced cutaneous photodamage

Pervasive Induction of Regulatory Mutation Microclones in Sun-exposed Skin

Single-cell transcriptomic analysis reveals the critical molecular pattern of UV-induced cutaneous squamous cell carcinoma

UV-exposure, endogenous DNA damage, and DNA replication errors shape the spectra of genome changes in human skin

Spatial constraints govern competition of mutant clones in human epidermis

Gene expression landscape of cutaneous squamous cell carcinoma progression

Genome-wide comparison of human keratinocyte and squamous cell carcinoma responses to UVB irradiation: implications for skin and epithelial cancer

Driver gene combinations dictate cutaneous squamous cell carcinoma disease continuum progression

Single-cell sequencing highlights heterogeneity and malignant progression in actinic keratosis and cutaneous squamous cell carcinoma

Landscape dynamic network biomarker analysis reveals the tipping point of transcriptome reprogramming to prevent skin photodamage

High burden and pervasive positive selection of somatic mutations in normal human skin

Skin cancer risk is increased by somatic mutations detected noninvasively in healthy-appearing sun-exposed skin

Chronic UV radiation–induced RORγt+ IL-22–producing lymphoid cells are associated with mutant KC clonal expansion

Retrospective Single-Center Case Study of Clinical Variables and the Degree of Actinic Elastosis Associated with Rare Skin Cancers

Human Digital Twin: Automated Cell Type Distance Computation and 3D Atlas Construction in Multiplexed Skin Biopsies

Deciphering UV‐induced DNA Damage Responses to Prevent and Treat Skin Cancer

Photodynamic therapy reduces the burden of small ultraviolet-induced epidermal clones in human and mouse skin

Integrative analysis of RNA-sequencing and microarray for the identification of adverse effects of UVB exposure on human skin

The landscape of genetic alterations of UVB-induced skin tumors in DNA repair-deficient mice