Association of TP53 and CDKN2A Mutation Profile with Tumor Mutation Burden in Head and Neck Cancer

Camrelizumab-based induction chemoimmunotherapy in locally advanced stage hypopharyngeal carcinoma: phase II clinical trial

This phase II trial aimed to determine the efficacy and safety of induction chemoimmunotherapy of camrelizumab plus modified TPF in locally advanced hypopharyngeal squamous cell carcinoma (LA HSCC) (NCT04156698). The primary endpoint was objective response rate (ORR), and secondary endpoints were 3-year overall survival (OS), progression-free survival (PFS), larynx preservation rate (LPR), and metastasis-free survival (MFS). Patients (cT3-4aN0-2M0), regardless of sex, received induction chemoimmunotherapy for three cycles: camrelizumab 200 mg d1, docetaxel 75 mg/m2 d1, cisplatin 25 mg/m2 d1-3, and capecitabine 800 mg/m2 bid d1-14, q21d. Patients were assigned to radioimmunotherapy if they had a complete or partial response, those with stable or progressive disease underwent surgery and adjuvant (chemo)radiotherapy. Camrelizumab was maintained post-radioimmunotherapy. Fifty-one patients were enrolled with a median follow-up duration of 23.7 months. After induction therapy, the ORR was 82.4% (42/51), meeting the prespecified endpoint. Grade 3/4 adverse events occurred in 26 patients, and no treatment-related death occurred. As three-year outcomes were immature, two-year OS, PFS and LPR were reported. As no distant metastatic event had occurred, MFS was not reported here. The two-year OS, PFS, and LPR rates were 83.0%, 77.1%, and 70.0%, respectively. The induction chemoimmunotherapy of camrelizumab plus TPF showed a high ORR rate with an acceptable safety profile in LA HSCC.

Clinical significance and prospective mechanism of increased CDKN2A expression in small cell lung cancer

BACKGROUND

Although it has been shown that cyclin dependent kinase inhibitor 2A (CDKN2A) plays a significant role in a number of malignancies, its clinicopathological value and function in small cell lung cancer (SCLC) is unclear and warrants additional research.

METHODS

The clinical significance of CDKN2A expression in SCLC was examined by multiple methods, including comprehensive integration of mRNA level by high throughput data, Kaplan-Meier survival analysis for prognostic value, and validation of its protein expression using in-house immunohistochemistry.

RESULTS

The expression of CDKN2A mRNA in 357 cases of SCLC was evidently higher than that in the control group (n = 525) combing the data from 20 research centers worldwide. The standardized mean difference (SMD) was 3.07, and the area under the curve (AUC) of summary receiver operating characteristic curve (sROC) was 0.97 for the overexpression of CDKN2A. ACC, COAD, KICH, KIRC, PCPG, PRAD, UCEC, UVM patients with higher CDKN2A expression had considerably worse overall survival rates than those with lower CDKN2A expression with the hazard ratio (HR) > 1.

CONCLUSION

CDKN2A upregulation extensively enhances the carcinogenesis and progression of SCLC.

Pan-cancer analysis of CDKN2A alterations identifies a subset of gastric cancer with a cold tumor immune microenvironment

BACKGROUND

Although CDKN2A alteration has been explored as a favorable factor for tumorigenesis in pan-cancers, the association between CDKN2A point mutation (MUT) and intragenic deletion (DEL) and response to immune checkpoint inhibitors (ICIs) is still disputed. This study aims to determine the associations of CDKN2A MUT and DEL with overall survival (OS) and response to immune checkpoint inhibitors treatment (ICIs) among pan-cancers and the clinical features of CDKN2A-altered gastric cancer.

METHODS

This study included 45,000 tumor patients that underwent tumor sequencing across 33 cancer types from four cohorts, the MSK-MetTropism, MSK-IMPACT, OrigiMed2020 and TCGA cohorts. Clinical outcomes and genomic factors associated with response to ICIs, including tumor mutational burden, copy number alteration, neoantigen load, microsatellite instability, tumor immune microenvironment and immune-related gene signatures, were collected in pan-cancer. Clinicopathologic features and outcomes were assessed in gastric cancer. Patients were grouped based on the presence of CDKN2A wild type (WT), CDKN2A MUT, CDKN2A DEL and CDKN2A other alteration (ALT).

RESULTS

Our research showed that CDKN2A-MUT patients had shorter survival times than CDKN2A-WT patients in the MSK MetTropism and TCGA cohorts, but longer OS in the MSK-IMPACT cohort with ICIs treatment, particularly in patients having metastatic disease. Similar results were observed among pan-cancer patients with CDKN2A DEL and other ALT. Notably, CDKN2A ALT frequency was positively related to tumor-specific objective response rates to ICIs in MSK MetTropism and OrigiMed 2020. Additionally, individuals with esophageal carcinoma or stomach adenocarcinoma who had CDKN2A MUT had poorer OS than patients from the MSK-IMPACT group, but not those with adenocarcinoma. We also found reduced levels of activated NK cells, T cells CD8 and M2 macrophages in tumor tissue from CDKN2A-MUT or DEL pan-cancer patients compared to CDKN2A-WT patients in TCGA cohort. Gastric cancer scRNA-seq data also showed that CDKN2A-ALT cancer contained less CD8 T cells but more exhausted T cells than CDKN2A-WT cancer. A crucial finding of the pathway analysis was the inhibition of three immune-related pathways in the CDKN2A ALT gastric cancer patients, including the interferon alpha response, inflammatory response, and interferon gamma response.

CONCLUSIONS

This study illustrates the CDKN2A MUT and DEL were associated with a poor outcome across cancers. CDKN2A ALT, on the other hand, have the potential to be used as a biomarker for choosing patients for ICI treatment, notably in esophageal carcinoma and stomach adenocarcinoma.

Identification and validation of a novel risk model based on cuproptosis‑associated m6A for head and neck squamous cell carcinoma

BACKGROUND

Head and neck squamous cell carcinoma (HNSCC) is a prevalent cancer with a poor survival rate due to anatomical limitations of the head and a lack of reliable biomarkers. Cuproptosis represents a novel cellular regulated death pathway, and N6-methyladenosine (m6A) is the most common internal RNA modification in mRNA. They are intricately connected to tumor formation, progression, and prognosis. This study aimed to construct a risk model for HNSCC using a set of mRNAs associated with m6A regulators and cuproptosis genes (mcrmRNA).

METHODS

RNA-seq and clinical data of HNSCC patients from The Cancer Genome Atlas (TCGA) database were analyzed to develop a risk model through the least absolute shrinkage and selection operator (LASSO) analysis. Survival analysis and receiver operating characteristic (ROC) analysis were performed for the high- and low-risk groups. Additionally, the model was validated using the GSE41613 dataset from the Gene Expression Omnibus (GEO) database. GSEA and CIBERSORT were applied to investigate the immune microenvironment of HNSCC.

RESULTS

A risk model consisting of 32 mcrmRNA was developed using the LASSO analysis. The risk score of patients was confirmed to be an independent prognostic indicator by multivariate Cox analysis. The high-risk group exhibited a higher tumor mutation burden. Additionally, CIBERSORT analysis indicated varying levels of immune cell infiltration between the two groups. Significant disparities in drug sensitivity to common medications were also observed. Enrichment analysis further unveiled significant differences in metabolic pathways and RNA processing between the two groups.

CONCLUSIONS

Our risk model can predict outcomes for HNSCC patients and offers valuable insights for personalized therapeutic approaches.

Genomic Profiling of Rare Undifferentiated Sarcomatoid Subtypes of Pancreatic Carcinomas: In Search of Therapeutic Targets

PURPOSE

The highly aggressive undifferentiated sarcomatoid carcinoma (USC) subtype of pancreatic ductal adenocarcinoma (PDAC) remains poorly characterized because of its rarity. Previous case reports suggest that immune checkpoint inhibitors could be a promising treatment strategy, but the prevalence of established predictive biomarkers of response is largely unknown. The objective of this study was to leverage comprehensive genomic profiling of USC PDAC tumors to determine the prevalence of biomarkers associated with potential response to targeted therapies.

METHODS

USC tumors (n = 20) underwent central pathology review by a board-certified gastrointestinal pathologist to confirm the diagnosis. These samples were compared with non-USC PDAC tumors (N = 5,562). Retrospective analysis of DNA and RNA next-generation sequencing data was performed.

RESULTS

USC PDACs were more frequently PD-L1+ by immunohistochemistry than non-USC PDAC (63% v 16%, respectively, P < .001). Furthermore, USC PDAC had an increase in neutrophils (8.99% v 5.55%, P = .005) and dendritic cells (1.08% v 0.00%, q = 0.022) and an increased expression of PDCD1LG2 (4.6% v 1.3%, q = 0.001), PDCD1 (2.0% v 0.8%, q = 0.060), and HAVCR2 (45.9% v 21.7%, q = 0.107) than non-USC PDAC. Similar to non-USC PDAC, KRAS was the most commonly mutated gene (86% v 90%, respectively, P = 1).

CONCLUSION

To our knowledge, this work represents the largest molecular analysis of USC tumors to date and showed an increased expression of immune checkpoint genes in USC tumors. These findings provide evidence for further investigation into immune checkpoint inhibitors in USC tumors.

ARTS and small-molecule ARTS mimetics upregulate p53 levels by promoting the degradation of XIAP

Mutations resulting in decreased activity of p53 tumor suppressor protein promote tumorigenesis. P53 protein levels are tightly regulated through the Ubiquitin Proteasome System (UPS). Several E3 ligases were shown to regulate p53 stability, including MDM2. Here we report that the ubiquitin E3 ligase XIAP (X-linked Inhibitors of Apoptosis) is a direct ligase for p53 and describe a novel approach for modulating the levels of p53 by targeting the XIAP pathway. Using in vivo (live-cell) and in vitro (cell-free reconstituted system) ubiquitylation assays, we show that the XIAP-antagonist ARTS regulates the levels of p53 by promoting the degradation of XIAP. XIAP directly binds and ubiquitylates p53. In apoptotic cells, ARTS inhibits the ubiquitylation of p53 by antagonizing XIAP. XIAP knockout MEFs express higher p53 protein levels compared to wild-type MEFs. Computational screen for small molecules with high affinity to the ARTS-binding site within XIAP identified a small-molecule ARTS-mimetic, B3. This compound stimulates apoptosis in a wide range of cancer cells but not normal PBMC (Peripheral Blood Mononuclear Cells). Like ARTS, the B3 compound binds to XIAP and promotes its degradation via the UPS. B3 binding to XIAP stabilizes p53 by disrupting its interaction with XIAP. These results reveal a novel mechanism by which ARTS and p53 regulate each other through an amplification loop to promote apoptosis. Finally, these data suggest that targeting the ARTS binding pocket in XIAP can be used to increase p53 levels as a new strategy for developing anti-cancer therapeutics.

Identification of AURKA as a Biomarker Associated with Cuproptosis and Ferroptosis in HNSCC

Cuproptosis and ferroptosis represent copper- and iron-dependent forms of cell death, respectively, and both are known to play pivotal roles in head and neck squamous cell carcinoma (HNSCC). However, few studies have explored the prognostic signatures related to cuproptosis and ferroptosis in HNSCC. Our objective was to construct a prognostic model based on genes associated with cuproptosis and ferroptosis. We randomly assigned 502 HSNCC samples from The Cancer Genome Atlas (TCGA) into training and testing sets. Pearson correlation analysis was utilized to identify cuproptosis-associated ferroptosis genes in the training set. Cox proportional hazards (COX) regression and least absolute shrinkage operator (LASSO) were employed to construct the prognostic model. The performance of the prognostic model was internally validated using single-factor COX regression, multifactor COX regression, Kaplan-Meier analysis, principal component analysis (PCA), and receiver operating curve (ROC) analysis. Additionally, we obtained 97 samples from the Gene Expression Omnibus (GEO) database for external validation. The constructed model, based on 12 cuproptosis-associated ferroptosis genes, proved to be an independent predictor of HNSCC prognosis. Among these genes, the increased expression of aurora kinase A (AURKA) has been implicated in various cancers. To further investigate, we employed small interfering RNAs (siRNAs) to knock down AURKA expression and conducted functional experiments. The results demonstrated that AURKA knockdown significantly inhibited the proliferation and migration of HNSCC cells (Cal27 and CNE2). Therefore, AURKA may serve as a potential biomarker in HNSCC.

Translational research in head and neck cancer: Molecular and immunological updates

Head and neck squamous cell carcinoma (HNSCC) has a poor prognosis. Each year, approximately 880,000 patients are newly diagnosed with HNSCC worldwide, and 450,000 patients with HNSCC die. Risk factors for developing HNSCC have been identified, with cigarette smoking, alcohol consumption, and viral infections being the major factors. Owing to the prevalence of human papillomavirus infection, the number of HNSCC cases is increasing considerably. Surgery and chemoradiotherapy are the primary treatments for HNSCC. With advancements in tumor biology, patients are eligible for novel treatment modalities, namely targeted therapies, immunotherapy, and photoimmunotherapy. Because this area of research has rapidly progressed, clinicians should understand the basic biology of HNSCC to choose an appropriate therapy in the upcoming era of personalized medicine. This review summarized recent developments in tumor biology, focusing on epidemiology, genetic/epigenetic factors, the tumor microenvironment, microbiota, immunity, and photoimmunotherapy in HNSCC, as well as how these findings can be translated into clinical settings.

Epigenetic MLH1 silencing concurs with mismatch repair deficiency in sporadic, naturally occurring colorectal cancer in rhesus macaques

BACKGROUND

Naturally occurring colorectal cancers (CRC) in rhesus macaques share many features with their human counterparts and are useful models for cancer immunotherapy; but mechanistic data are lacking regarding the comparative molecular pathogenesis of these cancers.

METHODS

We conducted state-of-the-art imaging including CT and PET, clinical assessments, and pathological review of 24 rhesus macaques with naturally occurring CRC. Additionally, we molecularly characterized these tumors utilizing immunohistochemistry (IHC), microsatellite instability assays, DNAseq, transcriptomics, and developed a DNA methylation-specific qPCR assay for MLH1, CACNA1G, CDKN2A, CRABP1, and NEUROG1, human markers for CpG island methylator phenotype (CIMP). We furthermore employed Monte-Carlo simulations to in-silico model alterations in DNA topology in transcription-factor binding site-rich promoter regions upon experimentally demonstrated DNA methylation.

RESULTS

Similar cancer histology, progression patterns, and co-morbidities could be observed in rhesus as reported for human CRC patients. IHC identified loss of MLH1 and PMS2 in all cases, with functional microsatellite instability. DNA sequencing revealed the close genetic relatedness to human CRCs, including a similar mutational signature, chromosomal instability, and functionally-relevant mutations affecting KRAS (G12D), TP53 (R175H, R273*), APC, AMER1, ALK, and ARID1A. Interestingly, MLH1 mutations were rarely identified on a somatic or germline level. Transcriptomics not only corroborated the similarities of rhesus and human CRCs, but also demonstrated the significant downregulation of MLH1 but not MSH2, MSH6, or PMS2 in rhesus CRCs. Methylation-specific qPCR suggested CIMP-positivity in 9/16 rhesus CRCs, but all 16/16 exhibited significant MLH1 promoter hypermethylation. DNA hypermethylation was modelled to affect DNA topology, particularly propeller twist and roll profiles. Modelling the DNA topology of a transcription factor binding motif (TFAP2A) in the MLH1 promoter that overlapped with a methylation-specific probe, we observed significant differences in DNA topology upon experimentally shown DNA methylation. This suggests a role of transcription factor binding interference in epigenetic silencing of MLH1 in rhesus CRCs.

CONCLUSIONS

These data indicate that epigenetic silencing suppresses MLH1 transcription, induces the loss of MLH1 protein, abrogates mismatch repair, and drives genomic instability in naturally occurring CRC in rhesus macaques. We consider this spontaneous, uninduced CRC in immunocompetent, treatment-naïve rhesus macaques to be a uniquely informative model for human CRC.

Radiomics Features on Enhanced Computed Tomography Predict FOXP3 Expression and Clinical Prognosis in Patients with Head and Neck Squamous Cell Carcinoma

Forkhead box P3 (FOXP3) has been identified as a novel molecular marker in various types of cancer. The present study assessed the expression of FOXP3 in patients with head and neck squamous cell carcinoma (HNSCC) and its potential as a clinical prognostic indicator, and developed a radiomics model based on enhanced computed tomography (CT) imaging. Data from 483 patients with HNSCC were downloaded from the Cancer Genome Atlas for FOXP3 prognostic analysis and enhanced CT images from 139 patients included in the Cancer Imaging Archives, which were subjected to the maximum relevance and minimum redundancy and recursive feature elimination algorithms for radiomics feature extraction and processing. Logistic regression was used to build a model for predicting FOXP3 expression. A prognostic scoring system for radiomics score (RS), FOXP3, and patient clinicopathological factors was established to predict patient survival. The area under the receiver operating characteristic (ROC) curve (AUC) and calibration curve and decision curve analysis (DCA) were used to evaluate model performance. Furthermore, the relationship between FOXP3 and the immune microenvironment, as well as the association between RS and immune checkpoint-related genes, was analyzed. Results of analysis revealed that patients with HNSCC and high FOXP3 mRNA expression exhibited better overall survival. Immune infiltration analysis revealed that FOXP3 had a positive correlation with CD4 + and CD8 + T cells and other immune cells. The 8 best radiomics features were selected to construct the radiomics model. In the FOXP3 expression prediction model, the AUC values were 0.707 and 0.702 for the training and validation sets, respectively. Additionally, the calibration curve and DCA demonstrated the positive diagnostic utility of the model. RS was correlated with immune checkpoint-related genes such as ICOS, CTLA4, and PDCD1. A predictive nomogram was established, the AUCs were 0.87, 0.787, and 0.801 at 12, 24, and 36 months, respectively, and DCA demonstrated the high clinical applicability of the nomogram. The enhanced CT radiomics model can predict expression of FOXP3 and prognosis in patients with HNSCC. As such, FOXP3 may be used as a novel prognostic marker to improve individualized clinical diagnosis and treatment decisions.

Developing and validating the model of tumor-infiltrating immune cell to predict survival in patients receiving radiation therapy for head and neck squamous cell carcinoma

Background

Radiotherapy (RT) is a mainstay of head and neck squamous cell carcinoma (HNSCC) treatment. Due to the influence of RT on tumor cells and immune/stromal cells in microenvironment, some studies suggest that immunologic landscape could shape treatment response. To better predict the survival based on genomic data, we developed a prognostic model using tumor-infiltrating immune cell (TIIC) signature to predict survival in patients undergoing RT for HNSCC.

Methods

Gene expression data and clinical information were downloaded from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. Data from HNSCC patients undergoing RT were extracted for analysis. TIICs prevalence in HNSCC patients was quantified by gene set variation analysis (GSVA) algorithm. TIICs and post-RT survival were analyzed using univariate Cox regression analysis and used to construct and validate a tumor-infiltrating cells score (TICS).

Results

Five of 26 immune cells were significantly associated with HNSCC prognosis in the training cohort (all P<0.05). Kaplan-Meier (KM) survival curves showed that patients in the high TICS group had better survival outcomes (log-rank test, P<0.05). Univariate analyses demonstrated that the TICS had independent prognostic predictive ability for RT outcomes (P<0.05). Patients with high TICS scores showed significantly higher expression of immune-related genes. Functional pathway analyses further showed that the TICS was significantly related to immune-related biological process. Stratified analyses supported integrating TICS and tumor mutation burden (TMB) into individualized treatment planning, as an adjunct to classification by clinical stage and human papillomavirus (HPV) infection.

Conclusions

The TICS model supports a personalized medicine approach to RT for HNSCC. Increased prevalence of TIIC within the tumor microenvironment (TME) confers a better prognosis for patients undergoing treatment for HNSCC.

Identification of a 5-Gene Cuproptosis Signature Predicting the Prognosis for Colon Adenocarcinoma Based on WGCNA

Background: Colorectal cancer is a highly aggressive malignant tumor that primarily affects the digestive system. It is frequently diagnosed at an advanced stage. Cuproptosis is a copper-dependent form cell death mechanism, distinct from all other known pathways underlying cell death, tumor progression, prognosis, and immune response. Although the role of cuproptosis in colorectal cancer has been investigated over time, there is still an urgent need to explore new methods and insights to understand its potential function. Methods: The Gene Expression Omnibus and The Cancer Genome Atlas gene expression data were systematically explored to investigate the role of cuproptosis in colon adenocarcinoma. The weighted gene coexpression network analysis was used to construct a gene coexpression network and identify the critical module and cuproptosis-related genes correlated with colon adenocarcinoma prognosis. A cuproptosis-related genes prognostic signature for colon adenocarcinoma was identified and validated. To validate the identified gene signature, quantitative reverse transcription-polymerase chain reaction was performed. Cell proliferation assays were analyzed by CCK8 and cell cycle detection. In addition, reactive oxygen species assay was also analyzed. Results: Five hub cuproptosis-related genes (Dihydrolipoamide S-acetyltransferase, Cyclin-dependent kinase inhibitor 2A, ATOX1, VEGFA, and ULK1) were screened and a prognostic risk model for predicting overall survival was established based on these genes. The model was successfully tested in the validation cohort and the GEPIA database. Colon adenocarcinoma patients were categorized into high-risk and low-risk groups based on risk scores. The study revealed that patients with higher risk scores were more likely to have a poor prognosis. Moreover, Dihydrolipoamide S-acetyltransferase was a tumor suppressor gene that can induce cell death and affected the redox reactions in the colon cancer cell line. Conclusions: These findings suggest that the newly identified 5-gene signature may serve as a more reliable prognostic factor than clinical factors such as age and stage of disease. These findings offer a theoretical foundation for further investigation into potential cuproptosis-related biomarkers for predicting colon adenocarcinoma prognosis in the future.

Predicting therapeutic responses in head and neck squamous cell carcinoma from TP53 mutation detected by cell-free DNA

Background

Head and neck squamous cell carcinoma (HNSCC) is an epithelial malignant tumor originating from the oral cavity, oropharynx, nasal cavity, sinuses, nasopharynx, hypopharynx, or larynx. Mutations in TP53 are the most common of all somatic genomic changes in HNSCC, and TP53 mutations are associated with the response to immunotherapy and chemotherapy. Tumor-derived circulating cell-free DNA (cfDNA) is a minimally invasive method to determining genetic alterations in cancer. This study aimed to explore the therapeutic responses of patients with HNSCC with TP53 mutation and the accuracy of cfDNA for detecting TP53 mutation.

Methods

Information on TP53 mutations, patient survival time, and clinical data in HNSCC were downloaded from The Cancer Genome Atlas database. The difference in immune infiltration between the TP53-mutant group and the wild-type group was compared. We applied the single-sample gene set enrichment analysis method on the transcriptome of HNSCC samples to assess the distribution of immune cell types between the two groups. The chemotherapy response was constructed using the R software package, "pRRophetic". Gene set enrichment analysis was performed based on the TP53 mutation. The next-generation sequencing was executed on cfDNA from nine patients with HNSCC to detect genetic alterations. Tumor biopsy (n=9) was sequenced using the same technique.

Results

TP53 was the most frequently mutated gene in HNSCC. The TP53 mutation was related to immune cells and the expression of immune-associated genes. The TP53 mutation group showed lower response to immunotherapy but high sensitivity to some chemotherapies compared with the wild-type group. TP53 was the most frequently mutated gene (6/9; 66.67%) in cfDNA. Only 27.27% of TP53 mutations in tumor tissue were detected outside of cfDNA.

Conclusions

TP53 mutation could be used as a specific predictor of treatment response in patients with HNSCC. Using cfDNA to detect the TP53 mutations in patients with HSNCC is a feasible method. The results suggested that the therapeutic response in patients could be predicted by detecting TP53 mutations in cfDNA, and large-scale and prospective studies are needed to validate this hypothesis.

A novel oxidative stress-related gene signature as an indicator of prognosis and immunotherapy responses in HNSCC

PURPOSE

To identify molecular subtypes of oxidative stress-related genes in head and neck squamous cell carcinoma (HNSCC) and to construct a scoring model of oxidative stress-related genes.

METHODS

R language based scRNA-seq and bulk RNA-seq analyses were used to identify molecular isoforms of oxidative stress-related genes in HNSCC. An oxidative stress-related gene scoring (OSRS) model was constructed, which were verified through online data and immunohistochemical staining of clinical samples.

RESULTS

Using TCGA-HNSCC datasets, nine predictive genes for overall patient survival, rarely reported in previous similar studies, were screened. AREG and CES1 were identified as prognostic risk factors. CSTA, FDCSP, JCHAIN, IFFO2, PGLYRP4, SPOCK2 and SPINK6 were identified as prognostic factors. Collectively, all genes formed a prognostic risk signature model for oxidative stress in HNSCC, which were validated in GSE41613, GSE103322 and PRJEB23709 datasets. Immunohistochemical staining of SPINK6 in nasopharyngeal cancer samples validated the gene panel. Subsequent analysis indicated that subgroups of the oxidative stress prognostic signature played important roles during cellular communication, the immune microenvironment, the differential activation of transcription factors, oxidative stress and immunotherapeutic responses.

CONCLUSIONS

The risk model might predict HNSCC prognosis and immunotherapeutic responses.

Has-miR-300-GADD45B promotes melanoma growth via cell cycle

Response to oncogenic factors like UV, GADD45 family in skin participates in scavenging ROS, DNA repair and cell cycle control. Because of this, the previous study of the chronic UVB injury model has found that hsa-miR-300 can conduct intercellular transport by exosomes and target regulation of GADD45B. Whether the hsa-miR-300-GADD45B still regulates tumor development by cell cycle pathway is unclear. Through transcriptomic analysis of primary (n=39) and metastatic (n=102) melanoma, it was confirmed that in metastatic samples, some of the 97 down-regulated genes participate in maintaining skin homeostasis while 42 up-regulated genes were enriched in cancer-related functions. Furthermore, CDKN1A, CDKN2A, CXCR4 and RAD51 in the melanoma pathway, were also differentially expressed between normal skin and melanoma. CDKN1A and CDKN2A were also found to be involved in TP53-dependent cell cycle regulation. In conclusion, it was speculated that CDKN1A, CDKN2A, TP53, GADD45B and hsa-miR-300 may have regulatory relationships. It was demonstrated that there is a bidirectional regulation between hsa-miR-300 and TP53. In addition, miR-300 can regulate CDKN1A by GADD45B/TP53 and promote melanoma growth by accelerating the cell cycle transition from G1/S to G2 phase.

Comprehensive Genome profile testing in head and neck cancer

OBJECTIVE

Head and neck cancer (HNC) is a tumor occurring in various primary sites with limited chemotherapy options for its treatment. Recently, comprehensive genomic profiling (CGP) testing has become clinically widespread. In this study, we examined the utility of CGP in diagnosing and treating HNC.

METHODS

This study included 29 patients with HNC who underwent CGP testing at the Gifu University Hospital between December 2019 and April 2022. We analyzed the types of gene mutations and tumor mutational burden (TMB) based on the CGP results. Squamous cell carcinoma accounted for 55.2%, and other cancers accounted for 44.8%. And we investigated the correlation of prognosis with gene mutations and TMB.

RESULTS

Gene mutations were detected in TP53(48.3%), CDKN2A (27.6%), CDKN2B (17.2%), NOTCH1 (17.2%), PIK3CA (17.2%), ARID1A (13.8%), and NF1 (13.8%). TP53, CDKN2A and CDKN2B mutations significantly decreased survival rate in HNC. Five cases (17.2%) were TMB-high and 82.8% were TMB-low. In SCC cases treated with immune checkpoint inhibitors, TMB-high had better Overall survival than TMB-low. And all patients with TMB-high were oropharyngeal cancer.

CONCLUSION

Although there were no cases in which effective treatment was actually performed based on the results of CGP, many gene mutations have been detected and several gene mutations correlated with prognosis. Furthermore, TMB can be used as a biomarker to predict the therapeutic effects of immune checkpoint inhibitors in cases of SCC.

Next-generation sequencing identifies CDKN2A alterations as prognostic biomarkers in recurrent or metastatic head and neck squamous cell carcinoma predominantly receiving immune checkpoint inhibitors

Background

This study aimed to identify potential biomarkers in patients with recurrent or metastatic head and neck squamous cell carcinoma (HNSCC) and further probe the prognostic implications of CDKN2A mutations, particularly within a subset receiving immunotherapy.

Methods

In this retrospective single-center study, we evaluated the next-generation sequencing (NGS) data from Foundation Medicine (FM) for patients with recurrent or metastatic HNSCC between January 1, 2019, and December 31, 2021. Patients were stratified based on CDKN2A loss-of-function (LOF) versus wild-type (WT) categorizations, with a focused subgroup analysis on those administered immunotherapy.

Results

The study encompassed 77 patients, of which 62 had undergone immunotherapy. The median duration of follow-up was 22.6 months. For the CDKN2A LOF group, the median overall survival (OS) was 16.5 months, contrasted with 30.0 months in the CDKN2A WT group (P=0.014). Notably, female gender (hazard ratio [HR]=4.526, 95% confidence interval [CI]: 1.934-10.180, P=0.0003) and CDKN2A LOF (HR=2.311, 95% CI: 1.156-4.748, P=0.019) emerged as independent risk factors for mortality in patients with recurrent or metastatic HNSCC. Within the immunotherapy subset, the median OS was 11.7 months for the CDKN2A LOF group, and 22.5 months for the CDKN2A WT group (P=0.017). Further, the female gender (HR=4.022, 95% CI: 1.417-10.710, P=0.006), CDKN2A LOF (HR=4.389, 95% CI: 1.782-11.460, P=0.002), and a combined positive score below 1 (HR=17.20, 95% CI: 4.134-79.550, P<0.0001) were identified as significant predictors of mortality among patients with recurrent or metastatic HNSCC receiving immunotherapy.

Conclusion

Alterations manifesting as LOF in the CDKN2A gene stand as robust indicators of unfavorable survival outcomes in HNSCC patients, including the subset that underwent immunotherapy.

Different Mutational Landscapes in Human Papillomavirus-Induced and Human Papillomavirus-Independent Invasive Penile Squamous Cell Cancers

Penile squamous cell carcinomas (SCC) are rare cancers that arise after transforming human papillomavirus (HPV) infections or independent of HPV in the background of chronic dermatoses. Limited knowledge about genetic alterations driving penile carcinogenesis comes from studies of mainly small cohorts of typically mixed etiology. In this comparative genetic study of HPV-induced and HPV-independent invasive penile SCC of 156 patients from a single institution in a low-incidence country, hotspots of 50 cancer-relevant genes were analyzed with targeted next-generation sequencing. Seventy-nine of 156 SCC were classified as HPV induced, and 77 of 156 SCC arose independent of HPV. Only 28 (35%) of 79 HPV-induced penile SCC, but 69 (90%) of 77 HPV-independent SCC carried somatic gene mutations. PIK3CA, FGFR3, and FBXW7 mutations occurred in both groups in similar numbers as seen in other human cancers. In contrast, mutations in TP53 (44/77; 57%), CDKN2A (35/77; 45%), and HRAS (13/77; 17%) genes occurred with one exception of a HIV positive patient exclusively in HPV-independent SCC with a frequent co-occurrence of TP53 and CDKN2A mutations (28/77; 42%). Mutations in multiple genes occurred in 9 (11%) of 79 HPV-induced SCC versus 47 (62%) of 77 HPV-independent SCC (χ2; P < .001). More than one mutation per gene (multi hits) was characteristic for HPV-independent SCC in 14 (18%) of 77 compared with only 3 (4%) of 79 HPV-induced SCC (χ2; P < .001). The total number of mutations in HPV-induced penile SCC (47 mutations) was significantly lower than that in HPV-independent SCC (143 mutations; Welsh test; P < .001). The presence of somatic driver gene mutations did not correlate with the age of patients, histology, or tumor stage of the primary SCC in either etiologic group, suggesting that acquisition of driver gene mutations is an early event after invasion. This large cohort analysis identified characteristic differences in mutational landscapes for the 2 etiologies. While genetic mutations in tumor suppressor genes drive HPV-independent penile carcinogenesis, oncogenic action of E6 and E7 substitute for mutations in HPV-induced SCC. A subgroup of patients with advanced SCC may be candidates for targeted therapy and clinical trials, although the majority of advanced penile SCC remain a therapeutic challenge.

TP53 or CDKN2A/B covariation in ALK/RET/ROS1-rearranged NSCLC is associated with a high TMB, tumor immunosuppressive microenvironment and poor prognosis

INTRODUCTION

ALK-rearranged lung adenocarcinomas with TP53 mutations have more unstable genomic features, poorer ALK-TKI efficacy and a worse prognosis than ALK-rearranged lung adenocarcinomas with wild-type TP53. Here, we examine the gene variations that co-occur with ALK/RET/ROS1 rearrangements in NSCLC and the corresponding tumor immune microenvironment, as well as their association with prognosis.

METHODS

A total of 155 patients with ALK/RET/ROS1 fusions were included retrospectively. Tumor genome mutation analysis was performed by next-generation sequencing. PD-L1 expression and tumor-infiltrating lymphocytes were assessed by multiplex immunohistochemistry. The correlations among gene covariation, the tumor immune microenvironment, and clinicopathological characteristics were analyzed.

RESULTS

Among the 155 patients, concomitant TP53 mutation appeared most frequently (31%), followed by CDKN2A/B copy number loss (15%). The ALK/RET/ROS1 fusion and TP53 or CDKN2A/B covariation group had more males and patients with stage IV disease (p < 0.001, p = 0.0066). Patients with TP53 or CDKN2A/B co-occurrence had higher tumor mutation burdens and more neoantigens (p < 0.001, p = 0.0032). PD-L1 expression was higher in the tumor areas of the TP53 or CDKN2A/B co-occurring group (p = 0.00038). However, the levels of CD8+, CD8+PD1-, and CD8+PD-L1- TILs were lower in the tumor areas of this group (p = 0.043, p = 0.029, p = 0.025). In the TCGA NSCLC cohorts, the top 2 mutated genes were CDKN2A/B (24%) and TP53 (16%). The TP53 or CDKN2A/B co-occurring group had higher tumor mutation burdens and shorter OS (p < 0.001, p < 0.001).

CONCLUSIONS

Patients with co-occurring TP53/CDKN2A/B variations and ALK/RET/ROS1 rearrangements are associated with high TMB, more neoantigens, an immunosuppressive microenvironment and a worse prognosis.

Mutant p53 in head and neck squamous cell carcinoma: Molecular mechanism of gain‑of‑function and targeting therapy (Review)

Head and neck squamous cell carcinoma (HNSCC) is one of the most widespread malignancies worldwide. p53, as a transcription factor, can play its role in tumor suppression by activating the expression of numerous target genes. However, p53 is one of the most commonly mutated genes, which frequently harbors missense mutations. These missense mutations are nucleotide substitutions that result in the substitution of an amino acid in the DNA binding domain. Most p53 mutations in HNSCC are missense mutations and the mutation rate of p53 reaches 65‑85%. p53 mutation not only inhibits the tumor suppressive function of p53 but also provides novel functions to facilitate tumor recurrence, called gain‑of‑function (GOF). The present study focused on the prevalence and clinical relevance of p53 mutations in HNSCC, and further described how mutant p53 accumulates. Moreover, mutant p53 in HNSCC can interact with proteins, RNA, and exosomes to exert effects on proliferation, migration, invasion, immunosuppression, and metabolism. Finally, several treatment strategies have been proposed to abolish the tumor‑promoting function of mutant p53; these strategies include reactivation of mutant p53 into wild‑type p53, induction of mutant p53 degradation, enhancement of the synthetic lethality of mutant p53, and treatment with immunotherapy. Due to the high frequency of p53 mutations in HNSCC, a further understanding of the mechanism of mutant p53 may provide potential applications for targeted therapy in patients with HNSCC.

Low-grade papillary Schneiderian carcinoma with TP53 mutation: a case report and review of the literature

BACKGROUND

Low-grade papillary Schneiderian carcinoma (LGPSC) is a relatively new entity of the sinonasal tract and is characterized by a bland morphology simulating sinonasal papilloma, invasive growth pattern with pushing borders, and aggressive clinical behavior with multiple recurrences and metastatic potential. Recently, DEK::AFF2 fusions were identified in LGPSC. However, some LPGSCs lack DEK::AFF2 fusion, and the molecular features of these tumors have not been clarified.

CASE PRESENTATION

A 69-year-old man presented with a discharge of pus from his left cheek. Computed tomography revealed a mass involving the left maxillary sinus, ethmoid sinus, and nasal cavity with the destruction of the orbital wall. The biopsy specimens showed that the tumor had a predominantly exophytic, papillary growth and did not have an apparent stromal invasion. The tumor was composed of multilayered epithelium that showed bland morphology with a round to polygonal shape, abundant eosinophilic cytoplasm, and uniform nuclei. Dense neutrophilic infiltrates were focally present. Immunohistochemically, CK5/6 was strongly and diffusely positive, and p16 was negative. p63 was mainly positive in the basal layer, and EMA was predominantly expressed in the outermost cell layer. DNA-based targeted sequencing showed TP53 R175H mutation, whereas neither EGFR nor KRAS mutation was identified. Reverse transcription polymerase chain reaction and fluorescence in situ hybridization revealed no DEK::AFF2 fusion.

CONCLUSIONS

We describe the first case of TP53-mutant LGPSC and review the literature. LGPSC is a genetically heterogeneous entity, and the recognition of this rare entity and comprehensive assessment of clinicopathological and molecular findings are crucial for the correct pathological diagnosis and clinical management.

Identification of p53-target genes in human papillomavirus-associated head and neck cancer by integrative bioinformatics analysis

Introduction

Head and neck cancer (HNC) is a highly prevalent and heterogeneous malignancy. Although extensive efforts have been made to advance its treatment, the prognosis remained poor with increased mortality. Human papillomaviruses (HPV) have been associated with high risk in HNC. TP53, a tumor suppressor, is the most frequently altered gene in HNC, therefore, investigating its target genes for the identification of novel biomarkers or therapeutic targets in HPV-related HNC progression is highly recommended.

Methods

Transcriptomic profiles from three independent gene expression omnibus (GEO) datasets, including 44 HPV+ and 70 HPV- HNC patients, were subjected to integrative statistical and Bioinformatics analyses. For the top-selected marker, further in-silico validation in TCGA and GTEx databases and experimental validation in 65 (51 HPV- and 14 HPV+) subjects with histologically confirmed head and neck squamous cell carcinoma (HNSCC) have been performed.

Results

A total of 498 differentially expressed genes (DEGs) were identified including 291 up-regulated genes and 207 down-regulated genes in HPV+ compared to HPV- HNSCC patients. Functional annotations and gene set enrichment analysis (GSEA) showed that the up-regulated genes were significantly involved in p53-related pathways. The integrative analysis between the Hub-genes identified in the complex protein-protein network and the top frequent genes resulting from GSEA showed an intriguing correlation with five biomarkers which are EZH2, MDM2, PCNA, STAT5A and TYMS. Importantly, the MDM2 gene showed the highest gene expression difference between HPV+ and HPV- HNSCC (Average log2FC = 1.89). Further in-silico validation in a large HNSCC cohort from TCGA and GTEx databases confirmed the over-expression of MDM2 in HPV+ compared to HPV- HNSCC patients (p = 2.39E-05). IHC scoring showed that MDM2 protein expression was significantly higher in HPV+ compared to HPV- HNSCC patients (p = 0.031).

Discussion

Our findings showed evidence that over-expression of MDM2, proto-oncogene, may affect the occurrence and proliferation of HPV-associated HNSCC by disturbing the p53-target genes and consequently the p53-related pathways.

Circulating Tumor DNA in Head and Neck Squamous Cell Carcinoma

Tumors shed cell-free DNA (cfDNA) into the plasma. “Liquid biopsies” are a diagnostic test to analyze cfDNA in order to detect minimal residual cancer, profile the genomic tumor landscape, and monitor cancers non-invasively over time. This technique may be useful in patients with head and neck squamous cell carcinoma (HNSCC) due to genetic tumor heterogeneity and limitations in imaging sensitivity. However, there are technical challenges that need to be overcome for the widespread use of liquid biopsy in the clinical management of these patients. In this review, we discuss our current understanding of HNSCC genetics and the role of cfDNA genomic analyses as an emerging precision diagnostic tool.

Identification and Biological Validation of a Chemokine/Chemokine Receptor-Based Risk Model for Predicting Immunotherapeutic Response and Prognosis in Head and Neck Squamous Cell Carcinoma

Over 80% of head and neck squamous cell carcinoma (HNSCC) patients failed to respond to immunotherapy, which can likely be attributed to the tumor microenvironment (TME) remolding mediated by chemokines/chemokine receptors (C/CR). This study aimed to establish a C/CR-based risk model for better immunotherapeutic responses and prognosis. After assessing the characteristic patterns of the C/CR cluster from the TCGA-HNSCC cohort, a six-gene C/CR-based risk model was developed to stratify patients by LASSO Cox analysis. The screened genes were multidimensionally validated by RT-qPCR, scRNA-seq, and protein data. A total of 30.4% of patients in the low-risk group had better responses to anti-PD-L1 immunotherapy. A Kaplan-Meier analysis showed that patients in the low-risk group had longer overall survival. A time-dependent receiver operating characteristic curve and Cox analyses indicated that risk score served as an independent predictive indicator. The robustness of the immunotherapy response and prognosis prediction was also validated in independent external datasets. Additionally, the TME landscape revealed that the low-risk group was immune activated. Furthermore, the cell communication analysis on the scRNA-seq dataset revealed that cancer-associated fibroblasts were the main communicators within the C/CR ligand-receptor network of TME. Collectively, The C/CR-based risk model simultaneously predicted immunotherapeutic response and prognosis, potentially optimizing personalized therapeutic strategies of HNSCC.

Identification of a cuproptosis-related lncRNA signature to predict the prognosis and immune landscape of head and neck squamous cell carcinoma

Background

Cuproptosis is considered a novel copper-induced cell death model regulated by targeting lipoylated TCA cycle proteins. In this study, we established a novel signature based on cuproptosis-related lncRNAs (crlncRNAs) to predict the prognosis and immune landscape of head and neck squamous cell carcinoma.

Methods

RNA-seq matrix, somatic mutation files, and clinical data were obtained from The Cancer Genome Atlas database. After dividing patients into two sets, a crlncRNA signature was established based on survival related crlncRNAs, which were selected by the univariate Cox analysis and least absolute shrinkage and selection operator Cox regression. To evaluate the model, Kaplan-Meier survival analysis and time-dependent receiver operating characteristic (ROC) were utilized, and a nomogram was established for survival prediction. Immune landscape analysis, drug sensitivity, cluster analysis, tumor mutation burden (TMB) and ceRNA network analysis were conducted subsequently.

Results

A crlncRNA related prognosis signature was finally constructed with 12 crlncRNAs. The areas under the ROC curves (AUCs) were 0.719, 0.705 and 0.693 respectively for 1, 3, and 5-year's overall survival (OS). Patients in the low-risk group behaved a better prognosis, lower TMB, higher immune function activity and scores. In addition, patients from cluster 2 were more sensitive to chemotherapy and immunotherapy.

Conclusion

In this study, we constructed a novel crlncRNA risk model to predict the survival of HNSCC patients. This reliable and acceptable prognostic signature may guide and promote the progress of novel treatment strategies for HNSCC patients.

John M, Fuentes-Bayne H, Patel SP, Mardiros A, Xu H, Ng EW, Go WY, Wong DJ, Sunwoo JB, Welch JS. Genomic Immune Evasion: Diagnostic and Therapeutic Opportunities in Head and Neck Squamous Cell Carcinomas

Head and neck squamous cell cancers (HNSCCs) represent a diverse group of tumors emerging within different mucosal surfaces of the oral cavity, nasopharynx, oropharynx, larynx, and hypopharynx. HNSCCs share common clinical risk factors and genomic features, including smoking, alcohol, age, male sex, aneuploidy, and TP53 mutations. Viral initiating and contributing events are increasingly recognized in HNSCCs. While both Epstein-Barr Virus (EBV) and human papilloma virus (HPV) are observed, EBV is more frequently associated with nasopharyngeal cancers whereas HPV is associated with oropharyngeal cancers. HNSCCs are associated with high tumor mutational burden and loss of tumor suppressor gene function, especially in TP53 and X-linked genes. Multiple lines of evidence suggest that HNSCCs are subject to immunologic surveillance and immune-induced evolutionary pressure that correlate with negative clinical outcomes. This review will discuss genomic mechanisms related to immune-mediated pressures and propose prognostic and therapeutic implications of detectable immune escape mechanisms that drive tumorigenesis and disease progression.

Identification of cuproptosis-related subtypes, cuproptosis-related gene prognostic index in hepatocellular carcinoma

Cuproptosis is a novel form of cell death, correlated with the tricarboxylic acid (TCA) cycle. However, the metabolic features and the benefit of immune checkpoint inhibitor (ICI) therapy based on cuproptosis have not yet been elucidated in Hepatocellular carcinoma (HCC). First, we identified and validated three cuproptosis subtypes based on 10 cuproptosis-related genes (CRGs) in HCC patients. We explored the correlation between three cuproptosis subtypes and metabolism-related pathways. Besides, a comprehensive immune analysis of three cuproptosis subtypes was performed. Then, we calculated the cuproptosis-related gene prognostic index (CRGPI) score for predicting prognosis and validated its predictive capability by Decision curve analysis (DCA). We as well explored the benefit of ICI therapy of different CRGPI subgroups in two anti-PD1/PD-L1 therapy cohorts (IMvigor210 cohort and GSE176307). Finally, we performed the ridge regression algorithm to calculate the IC50 value for drug sensitivity and Gene set enrichment analysis (GSEA) analysis to explore the potential mechanism. We found that cluster A presented a higher expression of FDX1 and was correlated with metabolism, glycolysis, and TCA cycle pathways, compared with the other two clusters. HCC patients with high CRGPI scores had a worse OS probability, and we further found that the CRGPI-high group had high expression of PD1/PDL1, TMB, and better response (PR/CR) to immunotherapy in the IMvigor210 cohort and GSE176307. These findings highlight the importance of CRGPI serving as a potential biomarker for both prognostic and immunotherapy for HCC patients. Generally, our results provide novel insights about cuproptosis into immune therapeutic strategies.

Immunization Combined with Ferroptosis Related Genes to Construct a New Prognostic Model for Head and Neck Squamous Cell Carcinoma

Simple Summary

Immunity combined with ferroptosis is being considered as a new tumor treatment modality, and its regulation in head and neck squamous cell carcinoma is still unknown. The purpose of this study was to look into the potential molecular biological roles of immune ferroptosis genes in head and neck squamous cell carcinoma. The 12-IFRM signatures were successfully constructed and classified into high- and low-risk groups using the TCGA database and related data resources. In patients with head and neck squamous cell carcinoma, feature-based risk scores were more predictive of survival than traditional clinicopathological features. Furthermore, the expression of CD8⁺T cells and macrophage M0 differed significantly between the two groups. The expression of TNFSF9 and CD44 in the high-risk groups was significantly increased compared with the low-risk groups. Next, we found a higher proportion of high-risk mutations than in the low-risk group. In addition, the high-risk group was more sensitive to some chemotherapy drugs. Finally, we performed correlation analysis on the model genes. In this paper, the 12-IFRM signatures was developed with promising application prospects for predicting the clinical outcomes and treatment outcomes in head and neck squamous cell carcinoma.

Abstract

Ferroptosis is a new type of programmed cell death that plays a pivotal role in a variety of tumors. Moreover, immunity is closely related to ferroptosis. However, immune-ferroptosis-related mRNAs (IFRMs) are still not fully understood in the regulation of head and neck squamous cell carcinoma (HNSC). The purpose of this paper was to investigate the IFRMs prediction of HNSC and its possible molecular biological role. RNA-Seq and related clinical data were mined from the TCGA database, ImmPort database, GeneCards database, FerrDb database, and previous data. In R software, the “DESeq2” package was used to analyze the differential expression of IFRMs. We used univariate Cox analysis to judge the prognosis of the IFRMs. Using the least absolute shrinkage and selection operator (LASSO) and Cox regression, a prediction model for 12 IFRMs was established. In this study, the Kaplan–Meier survival curve and receiver operating characteristic (ROC) curve analysis were used to evaluate the prediction results. Moreover, factors such as immune landscape, somatic mutations, and drug susceptibility are also discussed. We successfully constructed the signature of 12-IFRMs. The two risk groups were classified according to the risk score obtained by this signature. Compared with conventional clinicopathological features, the characteristic-based risk score was more predictive of survival in patients with HNSC. Furthermore, the expression of CD8⁺T cells and macrophage M0 differed significantly between the two groups. Moreover, the expression of TNFSF9 and CD44 in high-risk groups was significantly increased compared with the low-risk groups. Then, we found a higher proportion of high-risk mutations than in the low-risk group. Next, the high-risk group was more sensitive to chemotherapy drugs such as bosutinib, docetaxel, erlotinib, gefitinib, imatinib, lapatinib, and sorafenib. Finally, an in-depth analysis of the association and potential value of the 12 genes was performed. In summary, the 12-IFRM signatures established in this paper had good application prospects and could be effectively used to predict the clinical outcome and treatment response of head and neck squamous cell carcinoma.

Whole genome analysis reveals the genomic complexity in metastatic cutaneous squamous cell carcinoma

Metastatic cutaneous squamous cell carcinoma (CSCC) is a highly morbid disease requiring radical surgery and adjuvant therapy, which is associated with a poor prognosis. Yet, compared to other advanced malignancies, relatively little is known of the genomic landscape of metastatic CSCC. We have previously reported the mutational signatures and mutational patterns of CCCTC-binding factor (CTCF) regions in metastatic CSCC. However, many other genomic components (indel signatures, non-coding drivers, and structural variants) of metastatic CSCC have not been reported. To this end, we performed whole genome sequencing on lymph node metastases and blood DNA from 25 CSCC patients with regional metastases of the head and neck. We designed a multifaceted computational analysis at the whole genome level to provide a more comprehensive perspective of the genomic landscape of metastatic CSCC. In the non-coding genome, 3′ untranslated region (3′UTR) regions of EVC (48% of specimens), PPP1R1A (48% of specimens), and ABCA4 (20% of specimens) along with the tumor-suppressing long non-coding RNA (lncRNA) LINC01003 (64% of specimens) were significantly functionally altered (Q-value < 0.05) and represent potential non-coding biomarkers of CSCC. Recurrent copy number loss in the tumor suppressor gene PTPRD was observed. Gene amplification was much less frequent, and few genes were recurrently amplified. Single nucleotide variants driver analyses from three tools confirmed TP53 and CDKN2A as recurrently mutated genes but also identified C9 as a potential novel driver in this disease. Furthermore, indel signature analysis highlighted the dominance of ID signature 13 (ID13) followed by ID8 and ID9. ID9 has previously been shown to have no association with skin melanoma, unlike ID13 and ID8, suggesting a novel pattern of indel variation in metastatic CSCC. The enrichment analysis of various genetically altered candidates shows enrichment of “TGF-beta regulation of extracellular matrix” and “cell cycle G1 to S check points.” These enriched terms are associated with genetic instability, cell proliferation, and migration as mechanisms of genomic drivers of metastatic CSCC.