Transcriptomic analysis of cutaneous squamous cell carcinoma reveals a multigene prognostic signature associated with metastasis.

BACKGROUND: Metastasis of cutaneous squamous cell carcinoma (cSCC) is uncommon. Current staging methods are reported to have sub-optimal performances in metastasis prediction. Accurate identification of patients with tumors at high risk of metastasis would have a significant impact on management. OBJECTIVE: To develop a robust and validated gene expression profile signature for predicting primary cSCC metastatic risk using an unbiased whole transcriptome discovery-driven approach. METHODS: Archival formalin-fixed paraffin-embedded primary cSCC with perilesional normal tissue from 237 immunocompetent patients (151 nonmetastasizing and 86 metastasizing) were collected retrospectively from four centers. TempO-seq was used to probe the whole transcriptome and machine learning algorithms were applied to derive predictive signatures, with a 3:1 split for training and testing datasets. RESULTS: A 20-gene prognostic model was developed and validated, with an accuracy of 86.0%, sensitivity of 85.7%, specificity of 86.1%, and positive predictive value of 78.3% in the testing set, providing more stable, accurate prediction than pathological staging systems. A linear predictor was also developed, significantly correlating with metastatic risk. LIMITATIONS: This was a retrospective 4-center study and larger prospective multicenter studies are now required. CONCLUSION: The 20-gene signature prediction is accurate, with the potential to be incorporated into clinical workflows for cSCC.

Validation of four cutaneous squamous cell carcinoma staging systems using nationwide data.

BACKGROUND: Cutaneous squamous cell carcinoma (cSCC) is the second most common cancer worldwide with relatively low metastatic potential (2-5%). Developments in therapeutic options have highlighted the need to better identify high-risk patients who could benefit from closer surveillance, adjuvant therapies and baseline/follow-up imaging, while at the same time safely omitting low-risk patients from further follow-up. Controversy remains regarding the predictive performance of current cSCC staging systems and which methodology to adopt. OBJECTIVES: To validate the performance of four cSCC staging systems [American Joint Committee on Cancer 8th edition (AJCC8), Brigham and Women's Hospital (BWH), Tübingen and Salamanca T3 refinement] in predicting metastasis using a nationwide cohort. METHODS: A nested case-control study using data from the National Disease Registration Service, England, 2013-2015 was conducted. Metastatic cSCC cases were identified using an algorithm to identify all potential cases for manual review. These were 1 : 1 matched on follow-up time to nonmetastatic controls randomly selected from 2013. Staging systems were analysed for distinctiveness, homogeneity, monotonicity, specificity, positive predictive value (PPV), negative predictive value (NPV) and c-index. RESULTS: We included 887 metastatic cSCC cases and 887 nonmetastatic cSCC controls. The BWH system showed the highest specificity [92.8%, 95% confidence interval (CI) 90.8-94.3%, PPV (13.2%, 95% CI 10.6-16.2) and c-index (0.84, 95% CI 0.82-0.86). The AJCC8 showed superior NPV (99.2%, 95% CI 99.2-99.3), homogeneity and monotonicity compared with the BWH and Tübingen diameter and thickness classifications (P < 0.001). Salamanca refinement did not show any improvement in AJCC8 T3 cSCC staging. CONCLUSIONS: We validated four cSCC staging systems using the largest nationwide dataset of metastatic cSCC so far. Although the BWH system showed the highest overall discriminative ability, PPV was low for all staging systems, which shows the need for further improvement and refining of current cSCC staging systems.

Risk factors for metastatic cutaneous squamous cell carcinoma: Refinement and replication based on 2 nationwide nested case-control studies.

BACKGROUND: Risk factors for cutaneous squamous cell carcinoma (cSCC) metastasis have been investigated only in relatively small data sets. OBJECTIVE: To analyze and replicate risk factors for metastatic cSCC. METHODS: From English and Dutch nationwide cancer registry cohorts, metastatic cases were selected and 1:1 matched to controls. The variables were extracted from pathology reports from the National Disease Registration Service in England. In the Netherlands, histopathologic slides from the Dutch Pathology Registry were revised by a dermatopathologist. Model building was performed in the English data set using backward conditional logistic regression, whereas replication was performed using the Dutch data set. RESULTS: In addition to diameter and thickness, the following variables were significant risk factors for metastatic cSCC in the English data set (n = 1774): poor differentiation (odds ratio [OR], 4.56; 95% CI, 2.99-6.94), invasion in (OR, 1.69; 95% CI, 1.05-2.71)/beyond (OR, 4.43; 95% CI, 1.98-9.90) subcutaneous fat, male sex (OR, 2.59; 95% CI, 1.70-3.96), perineural/lymphovascular invasion (OR, 2.12; 95% CI, 1.21-3.71), and facial localization (OR, 1.57; 95% CI, 1.02-2.41). Diameter and thickness showed significant nonlinear relationships with metastasis. Similar ORs were observed in the Dutch data set (n = 434 cSCCs). LIMITATIONS: Retrospective use of pathology reports in the English data set. CONCLUSION: cSCC staging systems can be improved by including differentiation, clinical characteristics such as sex and tumor location, and nonlinear relationships for diameter and thickness.

Growth and Viability of Cutaneous Squamous Cell Carcinoma Cell Lines Display Different Sensitivities to Isoform-Specific Phosphoinositide 3-Kinase Inhibitors.

Cutaneous squamous cell carcinomas (cSCCs) account for about 20% of keratinocyte carcinomas, the most common cancer in the UK. Therapeutic options for cSCC patients who develop metastasis are limited and a better understanding of the biochemical pathways involved in cSCC development/progression is crucial to identify novel therapeutic targets. Evidence indicates that the phosphoinositide 3-kinases (PI3Ks)/Akt pathway plays an important role, in particular in advanced cSCC. Questions remain of whether all four PI3K isoforms able to activate Akt are involved and whether selective inhibition of specific isoform(s) might represent a more targeted strategy. Here we determined the sensitivity of four patient-derived cSCC cell lines to isoform-specific PI3K inhibitors to start investigating their potential therapeutic value in cSCC. Parallel experiments were performed in immortalized keratinocyte cell lines. We observed that pan PI3Ks inhibition reduced the growth/viability of all tested cell lines, confirming the crucial role of this pathway. Selective inhibition of the PI3K isoform p110α reduced growth/viability of keratinocytes and of two cSCC cell lines while affecting the other two only slightly. Importantly, p110α inhibition reduced Akt phosphorylation in all cSCC cell lines. These data indicate that growth and viability of the investigated cSCC cells display differential sensitivity to isoform-specific PI3K inhibitors.

A Unique Panel of Patient-Derived Cutaneous Squamous Cell Carcinoma Cell Lines Provides a Preclinical Pathway for Therapeutic Testing.

BACKGROUND: Cutaneous squamous cell carcinoma (cSCC) incidence continues to rise with increasing morbidity and mortality, with limited treatment options for advanced disease. Future improvements in targeted therapy will rely on advances in genomic/transcriptomic understanding and the use of model systems for basic research. We describe here the panel of 16 primary and metastatic cSCC cell lines developed and characterised over the past three decades in our laboratory in order to provide such a resource for future preclinical research and drug screening. METHODS: Primary keratinocytes were isolated from cSCC tumours and metastases, and cell lines were established. These were characterised using short tandem repeat (STR) profiling and genotyped by whole exome sequencing. Multiple in vitro assays were performed to document their morphology, growth characteristics, migration and invasion characteristics, and in vivo xenograft growth. RESULTS: STR profiles of the cSCC lines allow the confirmation of their unique identity. Phylogenetic trees derived from exome sequence analysis of the matched primary and metastatic lines provide insight into the genetic basis of disease progression. The results of in vivo and in vitro analyses allow researchers to select suitable cell lines for specific experimentation. CONCLUSIONS: There are few well-characterised cSCC lines available for widespread preclinical experimentation and drug screening. The described cSCC cell line panel provides a critical tool for in vitro and in vivo experimentation.

Mutations in AQP5, encoding a water-channel protein, cause autosomal-dominant diffuse nonepidermolytic palmoplantar keratoderma.

Autosomal-dominant diffuse nonepidermolytic palmoplantar keratoderma is characterized by the adoption of a white, spongy appearance of affected areas upon exposure to water. After exome sequencing, missense mutations were identified in AQP5, encoding water-channel protein aquaporin-5 (AQP5). Protein-structure analysis indicates that these AQP5 variants have the potential to elicit an effect on normal channel regulation. Immunofluorescence data reveal the presence of AQP5 at the plasma membrane in the stratum granulosum of both normal and affected palmar epidermis, indicating that the altered AQP5 proteins are trafficked in the normal manner. We demonstrate here a role for AQP5 in the palmoplantar epidermis and propose that the altered AQP5 proteins retain the ability to form open channels in the cell membrane and conduct water.

Type VII collagen regulates expression of OATP1B3, promotes front-to-rear polarity and increases structural organisation in 3D spheroid cultures of RDEB tumour keratinocytes.

Type VII collagen is the main component of anchoring fibrils, structures that are integral to basement membrane homeostasis in skin. Mutations in the gene encoding type VII collagen COL7A1 cause recessive dystrophic epidermolysis bullosa (RDEB) an inherited skin blistering condition complicated by frequent aggressive cutaneous squamous cell carcinoma (cSCC). OATP1B3, which is encoded by the gene SLCO1B3, is a member of the OATP (organic anion transporting polypeptide) superfamily responsible for transporting a wide range of endogenous and xenobiotic compounds. OATP1B3 expression is limited to the liver in healthy tissues, but is frequently detected in multiple cancer types and is reported to be associated with differing clinical outcome. The mechanism and functional significance of tumour-specific expression of OATP1B3 has yet to be determined. Here, we identify SLCO1B3 expression in tumour keratinocytes isolated from RDEB and UV-induced cSCC and demonstrate that SLCO1B3 expression and promoter activity are modulated by type VII collagen. We show that reduction of SLCO1B3 expression upon expression of full-length type VII collagen in RDEB cSCC coincides with acquisition of front-to-rear polarity and increased organisation of 3D spheroid cultures. In addition, we show that type VII collagen positively regulates the abundance of markers implicated in cellular polarity, namely ELMO2, PAR3, E-cadherin, B-catenin, ITGA6 and Ln332.

Novel CARD11 Mutations in Human Cutaneous Squamous Cell Carcinoma Lead to Aberrant NF-κB Regulation.

NF-κB signaling plays a crucial role in regulating proliferation and differentiation in the epidermis. Alterations in the NF-κB pathway can lead to skin pathologies with a significant burden to human health such as psoriasis and cutaneous squamous cell carcinoma (cSCC). Caspase recruitment domain (CARD)-containing scaffold proteins are key regulators of NF-κB signaling by providing a link between membrane receptors and NF-κB transcriptional subunits. Mutations in the CARD family member, CARD14, have been identified in patients with the inflammatory skin diseases psoriasis and pityriasis rubra pilaris. Here, we describe that the gene coding for another CARD scaffold protein, CARD11, is mutated in more than 38% of 111 cSCCs, and show that novel variants outside of the coiled-coil domain lead to constitutively activated NF-κB signaling. CARD11 protein expression was detectable in normal skin and increased in all cSCCs tested. CARD11 mRNA levels were comparable with CARD14 in normal skin and CARD11 mRNA was increased in cSCC. In addition, we identified CARD11 mutations in peritumoral and sun-exposed skin, suggesting that CARD11-mediated alterations in NF-κB signaling may be an early event in the development of cSCC.

The Promise of Genomics and the Development of Targeted Therapies for Cutaneous Squamous Cell Carcinoma.

Targeted therapies for melanoma and basal cell carcinoma have evolved from deciphering the molecular mechanisms involved in their tumorigenesis. Mutations in BRAF have led to clinical use of BRAF-inhibitors in advanced melanoma, and mutations in Hedgehog signaling to smoothened inhibitors in basal cell carcinoma. The development of tumor resistance to these treatments is leading to many new drug development initiatives and the exploration of multiple signalling pathways. Cutaneous squamous cell carcinoma continues to rise steeply in incidence with very limited therapeutic options for locally advanced or metastatic disease. New genetic technologies find significant levels of mutation in Notch gene family as well as other already recognized gene mutations, such as TP53. The mutational burden in cutaneous squamous cell carcinoma is massive, challenging the identification of driver genes and inhibiting translation from genomics to the clinic. Clinical experience with targeted therapies, such as epidermal growth factor receptor inhibitors, or immune modulatory drugs suggests that these agents may be of benefit to patients, while a more complete understanding of the mechanisms behind squamous cell carcinogenesis awaits further progress.

RHBDF2 mutations are associated with tylosis, a familial esophageal cancer syndrome.

Tylosis esophageal cancer (TOC) is an autosomal-dominant syndrome characterized by palmoplantar keratoderma, oral precursor lesions, and a high lifetime risk of esophageal cancer. We have previously localized the TOC locus to a small genomic interval within chromosomal region 17q25. Using a targeted capture array and next-generation sequencing, we have now identified missense mutations (c.557T>C [p.Ile186Thr] and c.566C>T [p.Pro189Leu] in RHBDF2, which encodes the inactive rhomboid protease RHBDF2 (also known as iRhom2), as the underlying cause of TOC. We show that the distribution of RHBDF2 in tylotic skin is altered in comparison with that in normal skin, and immortalized tylotic keratinocytes have decreased levels of total epidermal growth factor receptor (EGFR) and display an increased proliferative and migratory potential relative to normal cells, even when normal cells are stimulated with exogenous epidermal growth factor. It would thus appear that EGFR signaling is dysregulated in tylotic cells. Furthermore, we also show an altered localization of RHBDF2 in both tylotic and sporadic squamous esophageal tumors. The elucidation of a role of RHBDF2 in growth-factor signaling in esophageal cancer will help to determine whether targeting this pathway in chemotherapy for this and other squamous cell carcinomas will be effective.

Tyrosine dephosphorylation is required for Bak activation in apoptosis.

Activation of the cell-death mediator Bak commits a cell to mitochondrial apoptosis. The initial steps that govern Bak activation are poorly understood. To further clarify these pivotal events, we have investigated whether post-translational modifications of Bak impinge on its activation potential. In this study, we report that on apoptotic stimulation Bak undergoes dephosphorylation at tyrosine residue 108 (Y108), a critical event that is necessary but not sufficient for Bak activation, but is required both for early exposure of the occluded N-terminal domain and multimerisation. RNA interference (RNAi) screening identified non-receptor tyrosine phosphatases (PTPNs) required for Bak dephosphorylation and apoptotic induction through chemotherapeutic agents. Specifically, modulation of PTPN5 protein expression by siRNA and overexpression directly affected both Bak-Y108 phosphorylation and the initiation of Bak activation. We further show that MEK/ERK signalling directly affects Bak phosphorylation through inhibition of PTPN5 to promote cell survival. We propose a model of Bak activation in which the regulation of Bak dephosphorylation constitutes the initial step in the activation process, which reveals a previously unsuspected mechanism controlling the initiation of mitochondrial apoptosis.

ARNT controls the expression of epidermal differentiation genes through HDAC- and EGFR-dependent pathways.

Previously we showed that spatial and developmental modulation of ARNT (HIF1ß) expression in mouse epidermis is essential for maintenance of keratinocyte differentiation, proper formation of the barrier and normal desquamation. Here, using lentiviral suppression or induction of ARNT in TERT-immortalized (N-TERT) and HaCaT cells we assessed the nature and mechanisms of ARNT involvement in control of differentiation in human epidermal keratinocytes. ARNT depletion did not affect the levels of basal keratins K5 and K14, but significantly induced expression of several key differentiation markers (an effect abolished by EGF supplementation). Furthermore, ARNT deficiency resulted in the downregulation of amphiregulin (AREG) - the most highly expressed EGFR ligand in human keratinocytes - whereas upregulation of ARNT showed the opposite. In ARNT-deficient monolayer cultures and 3D epidermal equivalents, the downregulation of AREG was concurrent with a decline of EGFR and ERK1/2 phosphorylation. TSA, a potent suppressor of HDAC activity, abolished the effects of ARNT deficiency, implying a role for HDACs in ARNT-dependent modulation of the AREG-EGFR pathway and downstream epidermal genes. Total HDAC activity was significantly increased in ARNT-depleted cells and decreased with ARNT overexpression. ARNT-dependent shifts in HDAC activity were specifically attributed to significant changes in the levels of HDAC1, HDAC2 and HDAC3 proteins (but not mRNA) in both monolayer and 3D cultures. Collectively, our results suggest that ARNT controls AREG expression and the downstream EGFR-ERK pathway in keratinocytes, at least in part, by modulating HDAC activity. This novel regulatory pathway targeting advanced stages of epidermal differentiation might have important implications for skin pathology such as psoriasis, atopic dermatitis and cancer.

Loss-of-function mutations in Notch receptors in cutaneous and lung squamous cell carcinoma.

Squamous cell carcinomas (SCCs) are one of the most frequent forms of human malignancy, but, other than TP53 mutations, few causative somatic aberrations have been identified. We identified NOTCH1 or NOTCH2 mutations in ~75% of cutaneous SCCs and in a lesser fraction of lung SCCs, defining a spectrum for the most prevalent tumor suppressor specific to these epithelial malignancies. Notch receptors normally transduce signals in response to ligands on neighboring cells, regulating metazoan lineage selection and developmental patterning. Our findings therefore illustrate a central role for disruption of microenvironmental communication in cancer progression. NOTCH aberrations include frameshift and nonsense mutations leading to receptor truncations as well as point substitutions in key functional domains that abrogate signaling in cell-based assays. Oncogenic gain-of-function mutations in NOTCH1 commonly occur in human T-cell lymphoblastic leukemia/lymphoma and B-cell chronic lymphocytic leukemia. The bifunctional role of Notch in human cancer thus emphasizes the context dependency of signaling outcomes and suggests that targeted inhibition of the Notch pathway may induce squamous epithelial malignancies.

Mena regulates nesprin-2 to control actin-nuclear lamina associations, trans-nuclear membrane signalling and gene expression.

Interactions between cells and the extracellular matrix, mediated by integrin adhesion complexes, play key roles in fundamental cellular processes, including the sensing and transduction of mechanical cues. Here, we investigate systems-level changes in the integrin adhesome in patient-derived cutaneous squamous cell carcinoma cells and identify the actin regulatory protein Mena as a key node in the adhesion complex network. Mena is connected within a subnetwork of actin-binding proteins to the LINC complex component nesprin-2, with which it interacts and co-localises at the nuclear envelope. Moreover, Mena potentiates the interactions of nesprin-2 with the actin cytoskeleton and the nuclear lamina. CRISPR-mediated Mena depletion causes altered nuclear morphology, reduces tyrosine phosphorylation of the nuclear membrane protein emerin and downregulates expression of the immunomodulatory gene PTX3 via the recruitment of its enhancer to the nuclear periphery. We uncover an unexpected role for Mena at the nuclear membrane, where it controls nuclear architecture, chromatin repositioning and gene expression. Our findings identify an adhesion protein that regulates gene transcription via direct signalling across the nuclear envelope.

Driver gene combinations dictate cutaneous squamous cell carcinoma disease continuum progression.

The molecular basis of disease progression from UV-induced precancerous actinic keratosis (AK) to malignant invasive cutaneous squamous cell carcinoma (cSCC) and potentially lethal metastatic disease remains unclear. DNA sequencing studies have revealed a massive mutational burden but have yet to illuminate mechanisms of disease progression. Here we perform RNAseq transcriptomic profiling of 110 patient samples representing normal sun-exposed skin, AK, primary and metastatic cSCC and reveal a disease continuum from a differentiated to a progenitor-like state. This is accompanied by the orchestrated suppression of master regulators of epidermal differentiation, dynamic modulation of the epidermal differentiation complex, remodelling of the immune landscape and an increase in the preponderance of tumour specific keratinocytes. Comparative systems analysis of human cSCC coupled with the generation of genetically engineered murine models reveal that combinatorial sequential inactivation of the tumour suppressor genes Tgfbr2, Trp53, and Notch1 coupled with activation of Ras signalling progressively drives cSCC progression along a differentiated to progenitor axis. Taken together we provide a comprehensive map of the cSCC disease continuum and reveal potentially actionable events that promote and accompany disease progression.

Hypoxia-mediated control of HIF/ARNT machinery in epidermal keratinocytes.

Transcriptional activity of hypoxia-induced factor 1 (HIF1) - a heterodimer of HIF1α and ARNT (HIF1ß) - is essential for cellular adaptation to environmental stress and plays an important role in skin development, wound healing, tumorigenesis and barrier function. Using primary mouse and human epidermal keratinocytes at ambient or hypoxic (1% O(2)) conditions we studied effects of hypoxia upon HIF protein expression. Significant nuclear levels of ARNT and HIF1α along with high HIF1 activity in normoxic keratinocytes suggest an as yet uncharacterised oxygen-independent role for HIF pathway in the epidermis. Acute hypoxia results in an instant but transient increase of HIF1α protein accompanied by a gradual decrease in its mRNA, while ARNT expression remains unchanged. In prolonged (chronic) hypoxia both HIF1α and Arnt are downregulated along with decline of HIF1 activity. However, expression of classical HIF1 targets such as Selenbp1 and Vegfa remains high. Thus, keratinocytes respond to acute hypoxia with immediate block of HIF1α protein degradation and concomitant increase of HIF activity, while under chronic hypoxia pro-angiogenic signalling is maintained through HIF1-independent pathway(s). Decline of HIF1α during chronic exposure is controlled at both mRNA and protein levels, while Arnt is downregulated post-translationally. Distinct transcription levels of Hif1α and Hif3α splice variants under normoxia and their differential response to hypoxia suggest functional diversity of Hif-α isoforms and highlight the complexity of HIF machinery control in epidermal keratinocytes.

Metastatic cutaneous squamous cell carcinoma shows frequent deletion in the protein tyrosine phosphatase receptor Type D gene.

Cutaneous squamous cell carcinoma (cSCC) is the second most common form of nonmelanoma skin cancer (NMSC), and its incidence is increasing rapidly. Metastatic cSCC accounts for the majority of deaths associated with NMSC, but the genetic basis for cSCC progression remains poorly understood. A previous study identified small deletions (typically <1 Mb) in the protein tyrosine phosphatase receptor Type D (PTPRD) gene that segregated with more aggressive cSCC. To investigate the apparent association between deletion within PTPRD and cSCC metastasis, a series of 74 formalin-fixed paraffin-embedded tumors from 31 patients was analyzed using a custom Illumina 384 SNP microarray. Deletions were found in 37% of patients with metastatic cSCC and were strongly associated with metastatic tumors when compared to those that had not metastasized (p = 0.007). Subsequent mutation analysis revealed a higher mutation rate for PTPRD than has been reported in any other cancer type, with 37% of tumors harboring a somatic mutation. Conversely, bisulfite sequencing showed that methylation was not a mechanism of PTPRD disruption in cSCC. This is the first report to observe an association between deletion within PTPRD and metastatic disease and highlights the potential use of these deletions as a diagnostic biomarker for tumor progression. Combined with the high mutation rate observed in our study, PTPRD is one of the most commonly altered genes in cSCC and warrants further investigation to determine its significance for metastasis in other tumor types.

The reduction of water-soluble tetrazolium salt reagent on the plasma membrane of epidermal keratinocytes is oxygen dependent.

The water-soluble tetrazolium salt (WST-1) assay is frequently used to assess cell proliferation. However, our study showed that in normal and cancerous keratinocytes, this assay is more responsive to changes in oxygenation than to rates of cell growth. Stimulation of keratinocyte proliferation by low Ca(2+) and suppression of proliferation by nocodazole resulted in modest changes in WST-1 readings, whereas gradually reducing the level of oxygen in the cellular environment from ambient (21%) to near anoxic (0.1%) revealed a very strong negative correlation between cell oxygenation and WST-1 reagent reduction. In contrast, the very similar MTT [3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide] cell proliferation assay, which uses a different tetrazolium salt, showed no sensitivity to the level of oxygen. Unlike MTT, WST-1 reagent is reduced extracellularly through trans-plasma membrane transport (tPMET), thereby suggesting that tPMET is oxygen dependent. We propose that the WST-1 assay can be developed into a sensitive quantitative method to evaluate cell oxygenation in vitro and used to study the role of hypoxia and tPMET in homeostasis and disease (e.g., cancer). At the same time, WST-1 assay should be used cautiously to assess cell viability or proliferation because readings can be affected by certain extrinsic (low atmospheric oxygen or high density culture) or intrinsic (defects in oxygen-sensing pathways) factors.

The Genomic Landscape of Actinic Keratosis.

Actinic keratoses (AKs) are lesions of epidermal keratinocyte dysplasia and are precursors for invasive cutaneous squamous cell carcinoma (cSCC). Identifying the specific genomic alterations driving the progression from normal skin to skin with AK to skin with invasive cSCC is challenging because of the massive UVR-induced mutational burden characteristic at all stages of this progression. In this study, we report the largest AK whole-exome sequencing study to date and perform a mutational signature and candidate driver gene analysis on these lesions. We demonstrate in 37 AKs from both immunosuppressed and immunocompetent patients that there are significant similarities between AKs and cSCC in terms of mutational burden, copy number alterations, mutational signatures, and patterns of driver gene mutations. We identify 44 significantly mutated AK driver genes and confirm that these genes are similarly altered in cSCC. We identify azathioprine mutational signature in all AKs from patients exposed to the drug, providing further evidence for its role in keratinocyte carcinogenesis. cSCCs differ from AKs in having higher levels of intrasample heterogeneity. Alterations in signaling pathways also differ, with immune-related signaling and TGFß signaling significantly more mutated in cSCC. Integrating our findings with independent gene expression datasets confirms that dysregulated TGFß signaling may represent an important event in AK‒cSCC progression.