Mutation profiles in early-stage lung squamous cell carcinoma with clinical follow-up and correlation with markers of immune function

Molecular insights and clinical implications for the tumor suppressor role of SCFFBXW7 E3 ubiquitin ligase

FBXW7 is one of the most well-characterized F-box proteins, serving as substrate receptor subunit of SKP1-CUL1-F-box (SCF) E3 ligase complexes. SCFFBXW7 is responsible for the degradation of various oncogenic proteins such as cyclin E, c-MYC, c-JUN, NOTCH, and MCL1. Therefore, FBXW7 functions largely as a major tumor suppressor. In keeping with this notion, FBXW7 gene mutations or downregulations have been found and reported in many types of malignant tumors, such as endometrial, colorectal, lung, and breast cancers, which facilitate the proliferation, invasion, migration, and drug resistance of cancer cells. Therefore, it is critical to review newly identified FBXW7 regulation and tumor suppressor function under physiological and pathological conditions to develop effective strategies for the treatment of FBXW7-altered cancers. Since a growing body of evidence has revealed the tumor-suppressive activity and role of FBXW7, here, we updated FBXW7 upstream and downstream signaling including FBXW7 ubiquitin substrates, the multi-level FBXW7 regulatory mechanisms, and dysregulation of FBXW7 in cancer, and discussed promising cancer therapies targeting FBXW7 regulators and downstream effectors, to provide a comprehensive picture of FBXW7 and facilitate the study in this field.

Molecular and immune characterization of Chinese early-stage non-squamous non-small cell lung cancer: a multi-omics cohort study

Background

Albeit considered with superior survival, around 30% of the early-stage non-squamous non-small cell lung cancer (Ns-NSCLC) patients relapse within 5 years, suggesting unique biology. However, the biological characteristics of early-stage Ns-NSCLC, especially in the Chinese population, are still unclear.

Methods

Multi-omics interrogation of early-stage Ns-NSCLC (stage I-III), paired blood samples and normal lung tissues (n=76) by whole-exome sequencing (WES), RNA sequencing, and T-cell receptor (TCR) sequencing were conducted.

Results

An average of 128 exonic mutations were identified, and the most frequently mutant gene was EGFR (55%), followed by TP53 (37%) and TTN (26%). Mutations in MUC17, ABCA2, PDE4DIP, and MYO18B predicted significantly unfavorable disease-free survival (DFS). Moreover, cytobands amplifications in 8q24.3, 14q13.1, 14q11.2, and deletion in 3p21.1 were highlighted in recurrent cases. Higher incidence of human leukocyte antigen loss of heterozygosity (HLA-LOH), higher tumor mutational burden (TMB) and tumor neoantigen burden (TNB) were identified in ever-smokers than never-smokers. HLA-LOH also correlated with higher TMB, TNB, intratumoral heterogeneity (ITH), and whole chromosomal instability (wCIN) scores. Interestingly, higher ITH was an independent predictor of better DFS in early-stage Ns-NSCLC. Up-regulation of immune-related genes, including CRABP2, ULBP2, IL31RA, and IL1A, independently portended a dismal prognosis. Enhanced TCR diversity of peripheral blood mononuclear cells (PBMCs) predicted better prognosis, indicative of a noninvasive method for relapse surveillance. Eventually, seven machine-learning (ML) algorithms were employed to evaluate the predictive accuracy of clinical, genomic, transcriptomic, and TCR repertoire data on DFS, showing that clinical and RNA features combination in the random forest (RF) algorithm, with area under the curve (AUC) of 97.5% and 83.3% in the training and testing cohort, respectively, significantly outperformed other methods.

Conclusions

This study comprehensively profiled the genomic, transcriptomic, and TCR repertoire spectrums of Chinese early-stage Ns-NSCLC, shedding light on biological underpinnings and candidate biomarkers for prognosis development.

Molecular characterization of the evolution of premalignant lesions in the upper aerodigestive tract

Introduction

Early relapse and development of metastatic disease are some of the primary reasons for the poor prognosis of patients with head and neck squamous cell carcinoma (HNSCC). HNSCC is a heterogeneous disease which may develop in large premalignant fields of genetically altered cells. Yet knowing which individuals will progress and develop clinically significant cancers during their lifetimes remains one of the most important challenges of reducing HNSCC morbidity and mortality. To further elucidate the molecular mechanisms, we performed a focused analysis of the genome and immune microenvironment from multiple, matched normal squamous tissue, premalignant lesions, as well as primary and recurrent tumors from seven patients with p16-negative HNSCC.

Methods

We performed targeted panel Next Generation Sequencing (161 genes) to analyze somatic variants from sequentially collected, matched formalin-fixed paraffin-embedded tissue (normal, premalignant, HNSCC) from two patients. These samples plus samples from five additional patients were analyzed with the Nanostring PanCancer Immune Panel. In addition, we performed shallow whole genome sequencing (0.5x coverage on average) on samples from three of these patients. Patients were, apart from one case, primarily treated with curative-intent surgery, and received subsequent adjuvant treatment, if indicated.

Results

The most frequently mutated genes were TP53 and NOTCH1. Other mutated genes included NOTCH3 and CDKN2A, among others. A significant number of mutations were private to dysplasia and invasive carcinoma, respectively, however, almost 20% were shared between them. Increasing genomic instability was observed when comparing histologically normal squamous mucosa with higher levels of dysplasia. High-grade dysplasia showed similarly rearranged genomes as invasive carcinoma. Pathways related to interferon alpha and gamma response were upregulated even in moderate dysplastic lesions with increasing expression in higher grades of dysplasia and carcinoma. SPINK5, a known tumor suppressor gene in HNSCC, was already downregulated in low-grade dysplastic lesions, indicating an early deactivation in the evolution of the disease.

Conclusion

Genomic alterations as well as aberrant immune gene expression can be observed early in the evolution of tumors of the upper aerodigestive tract, highlighting the potential for targeting early mechanisms of disease progression.

A Novel Signature Based on m6A Regulator-Mediated Genes Along Glycolytic Pathway Predicts Prognosis and Immunotherapy Responses of Gastric Cancer Patients

N6-methyladenosine (m6A) modification is involved in many aspects of gastric cancer (GC). Moreover, m6A and glycolysis-related genes (GRGs) play important roles in immunotherapeutic and prognostic implication of GC. However, GRGs involved in m6A regulation have never been analyzed comprehensively in GC. Herein, the study aims to identify and validate a novel signature based on m6A-related GRGs in GC patients. Therefore, a m6A-related GRGs signature is established, which can predict the survival of patients with GC and remain an independent prognostic factor in multivariate analyses. Clinical significance of the model is well validated in internal cohort and independent validation cohort. In addition, the expression levels of risk model-related GRGs in clinical samples are validated. Consistent with the database results, all model genes are up-regulated in expression except DCN. After regrouping the patients based on this risk model, the study can effectively distinguish between them in respect to immune-cell infiltration microenvironment and immunotherapeutic response. Additionally, candidate drugs targeting risk model-related GRGs are confirmed. Finally, a nomogram combining risk scores and clinical parameters is created, and calibration plots show that the nomogram can accurately predict survival. This risk model can serve as a reliable assessment tool for predicting prognosis and immunotherapeutic responses in GC patients.

Prognostic analysis of lung squamous cell carcinoma patients with second primary malignancies: a SEER database study

Background

As lung squamous cell carcinoma (LUSC) patients are at increased risk of developing a second primary cancer, this complicates the patient's condition and thus makes prognostic assessment more difficult, posing a significant prognostic challenge for clinicians. Our goal was to assess the prognosis of LUSC patients with a second primary tumor, and provide insights into appropriate therapy and monitoring strategies.

Methods

Data was obtained for LUSC patients from the Surveillance, Epidemiology, and End Results (SEER) database. The LUSC patients were divided into three groups (LS-SPM, OT-LUSC and LUSC-only). Univariate and stratified analyses were performed for the baseline and clinical characteristics of the participants. Multiple regression and Kaplan-Meier survival analyses were also performed, followed by a final life table analysis.

Results

In our sample of 101,626 patients, the HR for OS in the LS-SPM group was 0.40 in univariate analysis. Kaplan-Meier survival curves showed that LS-SPM patients had considerably longer lifespans compared to the other groups. The LS-SPM patients had median and mean survival times of 64 months and 89.11 months. Unadjusted and adjusted multiple regression analyses showed that LS-SPM patients had a superior survival compared to LUSC-only and OT-LUSC groups.

Conclusion

LS-SPM patients have a good prognosis with aggressive therapy and immune monitoring. The present study offers novel insights into the pathophysiological causes and treatments for LS-SPM.

Immunological subtyping of salivary gland cancer identifies histological origin-specific tumor immune microenvironment

Gene expression analysis enhances proper cancer subtyping, a better understanding of the molecular characteristics of cancer, and strategies for precision medicine. However, salivary gland cancer (SGC) subtyping remains largely unexplored because of its rarity and diverse histopathological and immunological characteristics. This study aimed to determine whether the histological origin and immunological characteristics of SGC subtypes are intrinsic tumor immunity factors. We performed immune profiling of 94 RNA-seq of SGC tissues and found that the SGCs that originated from the excretory duct (ED), such as the salivary duct and mucoepidermoid carcinomas, exhibit higher immunity than those from the intercalated duct (ID), such as the adenoid cystic and myoepithelial carcinomas, based on the computationally predicted immune score (p < 0.001), immune cell enrichment in the tumor immune microenvironment (TIME) (p < 0.001), T-cell receptor diversity (p < 0.001), and expression of signal I (major histocompatibility complex, MHC, p < 0.001) and signal II (co-stimulatory, p < 0.001 and co-inhibitory, p < 0.001) genes. Further analysis revealed that tolerogenic dendritic cell-induced dysfunctional T-cell populations and T-cell exclusion in the TIME are the major immune evasive mechanisms of the ED-and ID-derived SGCs, respectively.

Genomic insights into the mechanisms of FGFR1 dependency in squamous cell lung cancer

Although subsets of patients with lung squamous cell carcinoma (LSCC) benefit from immunotherapy, there are few effective molecularly targeted treatments for LSCC. Fibroblast growth factor receptor (FGFR) inhibitors provide a therapeutic option for patients with LSCC harboring FGFR aberrations, but their therapeutic efficacy has been limited to date. In this issue of the JCI, Malchers et al. identified tail-to-tail rearrangements, either within or near FGFR1, that are associated with FGFR1 dependency and sensitivity to FGFR inhibition in LSCC. These results may help improve the selection of patients with LSCC who are most likely to benefit from treatment with FGFR inhibitors.

Machine learning analysis of lung squamous cell carcinoma gene expression datasets reveals novel prognostic signatures

BACKGROUND

Lung squamous cell carcinoma (LUSC) patients are often diagnosed at an advanced stage and have poor prognoses. Thus, identifying novel biomarkers for the LUSC is of utmost importance.

METHODS

Multiple datasets from the NCBI-GEO repository were obtained and merged to construct the complete dataset. We also constructed a subset from this complete dataset with only known cancer driver genes. Further, machine learning classifiers were employed to obtain the best features from both datasets. Simultaneously, we perform differential gene expression analysis. Furthermore, survival and enrichment analyses were performed.

RESULTS

The kNN classifier performed comparatively better on the complete and driver datasets' top 40 and 50 gene features, respectively. Out of these 90 gene features, 35 were found to be differentially regulated. Lasso-penalized Cox regression further reduced the number of genes to eight. The median risk score of these eight genes significantly stratified the patients, and low-risk patients have significantly better overall survival. We validated the robust performance of these eight genes on the TCGA dataset. Pathway enrichment analysis identified that these genes are associated with cell cycle, cell proliferation, and migration.

CONCLUSION

This study demonstrates that an integrated approach involving machine learning and system biology may effectively identify novel biomarkers for LUSC.

Immune-active tumor-adjacent tissues are associated with favorable prognosis in stage I lung squamous cell carcinoma

The immunogenomic features of tumor-adjacent lungs (TALs) in stage I lung squamous cell carcinoma (LUSC) are not clear. Multiomics analyses of tumor tissues and paired TALs from 59 stage I LUSC patients were performed. Compared to tumors, TALs exhibited a better-preserved immune contexture indicated by upregulation of immune pathways, increased immune infiltration, and higher expression of immune effector molecules. Notably, TALs had no mutations in PTEN and KEAP1, a lower incidence of human leukocyte antigen (HLA) loss and higher expression of HLA class I genes, major histocompatibility complex (MHC) I chaperones, and interferon (IFN)-γ-associated genes. Digital spatial profiling validated the generally higher immune infiltration in TALs and revealed a higher level of immune heterogeneity in LUSC tumors. Importantly, patients with higher immune infiltration in TALs had significantly longer survival, while high immune heterogeneity was associated with inferior patient survival. Our work can be considered in the selection of patients for adjuvant therapy, especially immunotherapy.

To be, or not to be: the dilemma of immunotherapy for non-small cell lung cancer harboring various driver mutations

INTRODUCTION

Lung cancer is one of primary cancer type with high incidence and mortality, non-small cell lung cancer (NSCLC) is the most common type of lung cncer. For advanced lung cancer, traditional chemotherapy and targeted therapy become difficult to solve the dilemma of further progress. In recent years, with the clinical application of immunotherapy, the therapeutic strategy of lung cancer has changed dramatically. At present, immunotherapy has shown conspicuous efficacy in NSCLC patients with high expression of programmed death-ligand 1 (PD-L1) and high tumor mutational burden (TMB). The discovery of driver mutations brings delightful hope for targeted cancer therapy. However, it remains controversial whether immunotherapy can be used in NSCLC patients with these specific driver mutations.

METHOD

This article summarized the latest research progresses of immunotherapy in advanced NSCLC. We paid close attention to the relevance of various driver mutations and immunotherapy in NSCLC patients, and summarized the predictive effects of several driver mutations and immunotherapy.

RESULTS

The mutations of KRAS, KRAS+TP53, EPHA (especially EPHA5), ZFHX3, ZFHX3+TP53, NOTCH, BRAF and LRP1B+FAT3 have potential to be used as biomarkers to predict the positive effectiveness of immunotherapy. ZFHX3, ZFHX3+TP53, STKII/LKB1+KEAP1+SMARCA4+PBRM1 mutations in LUAD patients get more positive effect in immunotherapy. While the mutations of EGFR, KEAP1, STKII/LKB1+KRAS, EML4-ALK, MET exon 14 skipping mutation, PBRM1, STKII/LKB1+KEAP1+SMARCA4+PBRM1, ERBB2, PIK3CA and RET often indicate poor benefit from immunotherapy.

CONCLUSION

Many gene mutations have been shown to be associated with immunotherapy efficacy. Gene mutations should be combined with PD-L1, TMB, etc. to predict the effect of immunotherapy.

A cuproptosis-associated long non-coding RNA signature for the prognosis and immunotherapy of lung squamous cell carcinoma

Cuproptosis, a copper-induced mechanism of mitochondrial-related cell death, has been implicated as a breakthrough in the treatment of cancer and has become a new treatment strategy. Furthermore, long non-coding RNA (lncRNA) can change the biological activities of tumor cells. Worldwide, lung squamous cell carcinoma (LUSC) is among the most common annoying tumors. LncRNAs related to cuproptosis are not researched at LUSC. Our research intends to develop a signature on the basis of cuproptosis-associated lncRNAs, which can predict LUSC prognosis and investigate LUSC immunological features. The TCGA database was used to retrieve LUSC transcriptome, clinical, and gene mutation data. For statistical analysis, we utilized the R program. We created a signature consisting of three cuproptosis-related lncRNAs in this investigation (including AC002467.1, LINC01740, and LINC02345). Survival analyses and Receiver Operating Characteristic curves demonstrated that this signature possessed powerful predictive capability. The signature's ability to predict was confirmed by a Receiver Operating Characteristic curve and principal component analysis. Notably, the patient with a high-risk score and a high tumor mutation burden level had a lower survival time. Furthermore, the tumor immune dysfunction and exclusion analysis showed that these individuals with low-risk scores may benefit from immunotherapy. The signature constructed by three cuproptosis-associated lncRNAs may be prognostic markers of LUSC. It contributes to immunotherapy and offers LUSC's therapy a new treatment direction.

Integrated bioinformatics analyses identifying key transcriptomes correlated with prognosis and immune infiltrations in lung squamous cell carcinoma

Background

Lung Squamous Cell Carcinoma (LUSC) is a major subtype of lung malignancies and is associated with the cause of cancer-mediated mortality worldwide. However, identification of transcriptomic signatures associated with survival-prognosis and immunity of tumor remains lacking.

Method

The GSE2088, GSE6044, GSE19188, GSE21933, GSE33479, GSE33532, and GSE74706 were integrated for identifying differentially expressed genes (DEGs) with combined effect sizes. Also, the TCGA LUSC cohort was used for further analysis. A series of bioinformatics methods were utilized for conducting the whole study.

Results

The 831 genes (such as DSG3, PKP1, DSC3, TPX2, and UBE2C) were found upregulated and the 731 genes (such as ABCA8, SELENBP1, FAM107A, and CACNA2D2) were downregulated in the LUSC. The functional enrichment analysis identifies the upregulated KEGG pathways, including cell cycle, DNA replication, base excision repair, proteasome, mismatch repair, and cellular senescence. Also, the key hub genes (such as EGFR, HRAS, JUN, CDH1, BRCA1, CASP3, RHOA, HDAC1, HIF1A, and CCNA2) were identified along with the eight gene modules that were significantly related to the protein-protein interaction (PPI). The clinical analyses identified that the overexpression group of CDH3, PLAU, PKP3, STIL, CALU, LOXL2, POSTN, DPP3, GALNT2, LOX, and ITPA are substantially associated with a poor survival prognosis and the downregulated group of IL18R1 showed a similar trend. Moreover, our investigation demonstrated that the survival-associated genes were correlated with the stromal and immune scores in LUSC, indicating that the survival-associated genes regulate tumor immunity. The survival-associated genes were genetically altered in 27% of LUSC patients and showed excellent diagnostic efficiency. Finally, the consistent expression level of CDH3, PLAU, PKP3, STIL, CALU, LOXL2, POSTN, DPP3, GALNT2, and ITPA were found in the TCGA LUSC cohort.

Conclusions

The identification of key transcriptomic signatures can be elucidated by the crucial mechanism of LUSC carcinogenesis.

Phenotypic plasticity promotes lymph nodes metastasis and drug resistance in lung squamous cell carcinomas

Background

Phenotypic plasticity (PP) is a major promoter of tumor metastasis and drug resistance. Nevertheless, the molecular features and clinical significance of phenotypic plasticity in lung squamous cell carcinomas (LSCC) remained largely unexplored.

Methods

Phenotypic plasticity-related genes (PPRG) and clinical information of LSCC were downloaded from the cancer genome atlas (TCGA). The expression profiles of PPRG were compared between patients with and without lymph node metastasis. The prognostic signature was constructed, and survival analysis was performed based on phenotypic plasticity. Immunotherapy responses, chemotherapeutic drugs and targeted drug responses were investigated. In addition, the results were verified in an external cohort.

Results

Patients with and without lymph node metastasis exhibited significantly different genomic characteristics of phenotypic plasticity. Enrichment analysis showed that PP was strongly associated with cell responses and cell contraction. Survival analysis demonstrated that PPRG could serve as independent prognostic factor for overall survival. The phenotypic plasticity-related signature successfully divided patients into high- and low-PP score groups. Patients with low-PP scores were more sensitive to PD-L1, Cisplatin, Gefitinib, Obatoclax. Mesylate, Paclitaxel, Sorafenib and Vinorelbine (all p < 0.05). While, patients with low-PP scores were more sensitive to Axitinib and Camptothecin (all p < 0.05). Consistent with the results from TCGA, the external cohort validated the above findings.

Conclusions

Our study revealed that phenotypic plasticity may be involved in the lymph node metastasis in LSCC through regulating cell responses and cell contraction. Evaluation of phenotypic plasticity will assist clinicians in making treatment strategies.

Identification of novel biomarkers for lung squamous cell carcinoma

Lung squamous cell carcinoma (LUSC) is the second most common subtype of lung cancer, accounting for a majority of lung cancer-related deaths. Detection or diagnosis of cancer at an early stage is an unmet clinical need that is being actively explored. In this study, we aimed to identify potential biomarkers for LUSC, by screening expression status of all human genes against LUSC patient samples available with The Cancer Genome Atlas (TCGA). This led to the identification of several genes that are upregulated in LUSC. Further analysis revealed that many of these genes also show higher expression at the protein level not only in lung cancer but also in other cancers. Additionally, some of these genes show stage-dependent higher expression and are associated with statistically significant poor survival of LUSC patients. As per our results, more than 60 genes are overexpressed in LUSC at the level of mRNA and some at the protein level. Thus, we identified genes such as MCC1, MRPL47, CRYGS, HSP40, DNAJC19, GMPS and PARL as novel potential biomarkers for LUSC in this study. We believe that these genes hold great potential as LUSC biomarkers for early detection as the data are derived from patient samples.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03489-z.

Docetaxel Plus Nedaplatin or Carboplatin as Second-Line Chemotherapy for Advanced Lung Squamous Cell Carcinoma in Real-World Practice: A Single-Center Experience

Although single-drug chemotherapy regimens were used as second-line chemotherapy for advanced lung squamous cell carcinoma (LSCC) patients, there are still no standard guidelines for second-line chemotherapy. The purpose of this study was to compare the efficacy and safety of docetaxel combined with nedaplatin or carboplatin in the second-line treatment of advanced LSCC patients. One hundred and ninety-six LSCC patients receiving docetaxel plus nedaplatin (DN, n = 96) or carboplatin (DC, n = 100) were retrospectively collected until disease progression or unacceptable toxicity. The progression-free survival (PFS), overall survival (OS), objective response rate (ORR), disease control rate (DCR), and adverse events (AEs) were analyzed in the two groups. The ORR was 18.8% versus 16.0%, and the DCR was 39.6% versus 34.0% in DN group and DC group (P > .05 and P > .05), respectively. The PFS was 5.3 versus 3.8 months, and the OS was 8.5 and 6.7 months in DN group and DC group (P = .013 and P = .404), respectively. The rate of digestive reaction and hepatotoxicity was similar in DN and DC groups, whereas more patients in DC group than in DN group suffered from leucopenia (P < .05). Docetaxel combined with nedaplatin is an effective regimen for advanced LSCC patients. Compared with a similar regimen with carboplatin, the response rate was similar; however, nedaplatin regimen shows some superiority as regards survival and some treatment side effect.

APOBEC mutagenesis and selection for NFE2L2 contribute to the origin of lung squamous-cell carcinoma

Lung squamous-cell carcinoma originates as a consequence of oncogenic molecular variants arising from diverse mutagenic processes such as tobacco, defective homologous recombination, aging, and cytidine deamination by APOBEC proteins. Only some of the many variants generated by these processes actually contribute to tumorigenesis. Therefore, molecular investigation of mutagenic processes such as cytidine deamination by APOBEC should also determine whether the mutations produced by these processes contribute substantially to the growth and survival of cancer. Here, we determine the processes that gave rise to mutations of 681 lung squamous-cell carcinomas, and quantify the probability that each mutation was the product of each process. We then calculate the contribution of each mutation to increases in cellular proliferation and survival. We performed in vitro experiments to determine cytidine deamination activity of APOBEC3B against oligonucleotides corresponding with genomic sequences that give rise to variants of high cancer effect size. The largest APOBEC-related cancer effects are attributable to mutations in PIK3CA and NFE2L2. We demonstrate that APOBEC effectively deaminates NFE2L2 at the locations that confer high cancer effect. Overall, we demonstrate that APOBEC activity can lead to mutations in NFE2L2 that have large contributions to cancer cell growth and survival, and that NFE2L2 is an attractive potential target for therapeutic intervention.

Immunogenomic profiling of lung adenocarcinoma reveals poorly differentiated tumors are associated with an immunogenic tumor microenvironment

OBJECTIVES

Pathologists have routinely observed distinct histologic patterns of growth in early-stage lung adenocarcinoma (LUAD), which have been suggested to be associated with prognosis. Herein, we investigated the relationship between LUAD patterns of growth, as defined by the updated international association for the study of lung cancer (IASLC) grading criteria, and differences in the tumor immune microenvironment to identify predictors of response to immunotherapy.

METHODS

174 resected stage I-III LUAD tumors were classified by histologic pattern of growth (i.e. solid, micropapillary, acinar, papillary, and lepidic) and then grouped as well differentiated, moderately differentiated, and poorly differentiated. Comprehensive multiplatform analysis including whole exome sequencing, gene expression profiling, immunohistochemistry, CIBERSORT, and T-cell receptor sequencing was performed and groups were compared for differences in genomic drivers, immune cell infiltrate, clonality, and survival. Finally, multivariate analysis was performed adjusting for pathologic stage and smoking status.

RESULTS

Poorly differentiated tumors demonstrated a strong association with smoking relative to moderately differentiated or well differentiated tumors. However, unlike in prior reports, poorly differentiated tumors were not associated with a worse survival after curative-intent resection. Genomic analysis revealed that poorly differentiated tumors are associated with high tumor mutation burden but showed no association with oncogenic drivers. Immune analyses revealed that poorly differentiated tumors are associated with increased T-cell clonality, expression of PD-L1, and infiltration by cytotoxic CD8 T-cells, activated CD4 T-cells, and pro-inflammatory (M1) macrophages. Finally, multivariate analysis controlling for stage and smoking status confirmed independence of immune differences between IASLC grade groups.

CONCLUSIONS

Poorly differentiated tumors, as defined by the updated IASLC grading criteria, are associated with a distinct immunogenic tumor microenvironment that predicts for therapeutic response to immune agents, including checkpoint inhibitors, and should be included in the clinical trial design of immunotherapy studies in early-stage lung adenocarcinoma.

Co-dependencies in the tumor immune microenvironment

Activated oncogenes and disrupted tumor suppressor genes (TSGs) not only endow aspiring cancer cells with new biological capabilities but also influence the composition and function of host cells in the tumor microenvironment (TME). These non-cancer host cells can in turn provide cancer cells with growth support and protection from the anti-tumor immune response. In this ecosystem, geospatially heterogenous "subTME" adds to the complexity of the "global" TME which bestows tumors with increased tumorigenic ability and resistance to therapy. This review highlights how specific genetic alterations in cancer cells establish various symbiotic co-dependencies with surrounding host cells and details the cooperative role of the host cells in tumor biology. These essential interactions expand the repertoire of targets for the development of precision cancer treatments.

Recent Insight on Regulations of FBXW7 and Its Role in Immunotherapy

SCF^FBXW7 E3 ubiquitin ligase complex is a crucial enzyme of the ubiquitin proteasome system that participates in variant activities of cell process, and its component FBXW7 (F-box and WD repeat domain–containing 7) is responsible for recognizing and binding to substrates. The expression of FBXW7 is controlled by multiple pathways at different levels. FBXW7 facilitates the maturity and function maintenance of immune cells via functioning as a mediator of ubiquitination-dependent degradation of substrate proteins. FBXW7 deficiency or mutation results in the growth disturbance and dysfunction of immune cell, leads to the resistance against immunotherapy, and participates in multiple illnesses. It is likely that FBXW7 coordinating with its regulators and substrates could offer potential targets to improve the sensitivity and effects of immunotherapy. Here, we review the mechanisms of the regulation on FBXW7 and its tumor suppression role in immune filed among various diseases (mostly cancers) to explore novel immune targets and treatments.

Clinical significance of FBXW7 loss of function in human cancers

FBXW7 (F-Box and WD Repeat Domain Containing 7) (also referred to as FBW7 or hCDC4) is a component of the Skp1-Cdc53 / Cullin-F-box-protein complex (SCF/β-TrCP). As a member of the F-box protein family, FBXW7 serves a role in phosphorylation-dependent ubiquitination and proteasome degradation of oncoproteins that play critical role(s) in oncogenesis. FBXW7 affects many regulatory functions involved in cell survival, cell proliferation, tumor invasion, DNA damage repair, genomic instability and telomere biology. This thorough review of current literature details how FBXW7 expression and functions are regulated through multiple mechanisms and how that ultimately drives tumorigenesis in a wide array of cell types. The clinical significance of FBXW7 is highlighted by the fact that FBXW7 is frequently inactivated in human lung, colon, and hematopoietic cancers. The loss of FBXW7 can serve as an independent prognostic marker and is significantly correlated with the resistance of tumor cells to chemotherapeutic agents and poorer disease outcomes. Recent evidence shows that genetic mutation of FBXW7 differentially affects the degradation of specific cellular targets resulting in a distinct and specific pattern of activation/inactivation of cell signaling pathways. The clinical significance of FBXW7 mutations in the context of tumor development, progression, and resistance to therapies as well as opportunities for targeted therapies is discussed.

Real-world survival analysis by tumor mutational burden in non-small cell lung cancer: a multisite U.S. study

BACKGROUND

Tumor mutational burden (TMB) is a potential biomarker to predict tumor response to immuno-oncology agents in patients with metastatic non-small cell lung cancer (NSCLC).

MATERIALS AND METHODS

A multi-site cohort study evaluated patients diagnosed with stage IV NSCLC between 2012 and 2019 who had received comprehensive genomic profiling (CGP) and any NSCLC-related treatment at 9 U.S. cancer centers. Baseline characteristics and clinical outcomes were compared between patients with TMB <10 and TMB ≥10.

RESULTS

Among the 667 patients with CGP results, most patients received CGP from Foundation Medicine (64%) or Caris (20%). Patients with TMB ≥10 (vs. TMB <10) were associated with a positive smoking history. TMB was associated with ALK (p = 0.01), EGFR (p < 0.01), and TP53 (p < 0.05) alterations. TMB >10 showed a significant association towards longer overall survival (OS) (HR: 0.43, 95% CI: 0.21-0.88, p = 0.02) and progression-free survival (PFS) (HR: 0.43, 95% CI: 0.21-0.85, p = 0.02) in patients treated with first-line immunotherapy and tested by Foundation Medicine or Caris at treatment initiation.

CONCLUSIONS

TMB levels greater than or equal to 10 mut/Mb, when tested by Foundation Medicine or Caris at treatment initiation, were significantly associated with improved OS and PFS among patients treated with first-line immunotherapy-containing regimens. Additional prospective research is warranted to validate this biomarker along with PD-L1 expression.

Deciphering the Immune-Tumor Interplay During Early-Stage Lung Cancer Development via Single-Cell Technology

Lung cancer is the leading cause of cancer-related death worldwide. Cancer immunotherapy has shown great success in treating advanced-stage lung cancer but has yet been used to treat early-stage lung cancer, mostly due to lack of understanding of the tumor immune microenvironment in early-stage lung cancer. The immune system could both constrain and promote tumorigenesis in a process termed immune editing that can be divided into three phases, namely, elimination, equilibrium, and escape. Current understanding of the immune response toward tumor is mainly on the "escape" phase when the tumor is clinically detectable. The detailed mechanism by which tumor progenitor lesions was modulated by the immune system during early stage of lung cancer development remains elusive. The advent of single-cell sequencing technology enables tumor immunologists to address those fundamental questions. In this perspective, we will summarize our current understanding and big gaps about the immune response during early lung tumorigenesis. We will then present the state of the art of single-cell technology and then envision how single-cell technology could be used to address those questions. Advances in the understanding of the immune response and its dynamics during malignant transformation of pre-malignant lesion will shed light on how malignant cells interact with the immune system and evolve under immune selection. Such knowledge could then contribute to the development of precision and early intervention strategies toward lung malignancy.

Construction and Validation of Prognostic Regulation Network Based on RNA-Binding Protein Genes in Lung Squamous Cell Carcinoma

Lung squamous cell carcinoma (LUSC) is a common histologic subtype of non-small cell lung cancer with a poor prognosis. RNA-binding proteins (RBPs) are key modulators in the posttranscriptional regulation and RBP alterations are commonly found in various cancer types. However, its roles in predicting the tumorigenesis and prognosis have not been identified in LUSC. To identify the roles of RBPs in the tumorigenesis and prognosis of LUSC, the RNA sequencing data of patients with LUSC were retrieved from The Cancer Genome Atlas (TCGA) databases. The differential expressed genes (DEGs) were evaluated and identified. The intersection of manually curated RBPs and tumorigenesis-related DEGs was filtered to the univariate Cox regression analysis. The intersection genes with prognostic value were defined as prognostic RNA-binding protein genes (PRBPGs). Based on them, the predicted model was constructed. Its accuracy was tested by the area under the curve (AUC) of the receiver operator characteristic curve and the risk score. In addition, to explore the key regulatory network, the relationship among PRBPGs, target RNA, and absolute quantification of 50 hallmarks of cancer was also identified by Pearson correlation analysis. A total of 311 genes were filtered as the intersection of 1542 manually curated RBPs and tumorigenesis-related DEGs and the results revealed 17 PRBPGs. Based on them, we constructed the predict model with a relatively high accuracy (AUC: 0.739). The Kaplan-Meier survival curve showed the significant prognostic value of risk score (p < 0.001). Moreover, we uncovered the regulatory networks of PHF5A-TOMM22-oxidative phosphorylation, TLR3-CTSO inflammation-related pathway, SECISBP2L-targeted RNA (ADGRF5, TGFBR2, CD302, AC096921.2, AHCYL2, RPS6KA2, SLC34A2, and SFTPB) angiogenesis, and SECISBP2L-AKAP13 signaling (DNA repair, MTORC1 signaling, and MYC targets). The regulation mechanisms and cellular location of key PRBPGs were validated by assay for targeting accessible chromatin with high-throughput sequencing and single-cell RNA sequencing. Our study identifies PRBPGs as reliable indexes in predicting the tumorigenesis and prognosis of patients with LUSC and provides a well-applied model for predicting the overall survival for patients with LUSC. Besides, we also identified the regulatory network among PRBPGs, target RNA, and cancer gene sets in mediating the LUSC tumorigenesis.

Association of PTPRT mutations with immune checkpoint inhibitors response and outcome in melanoma and non-small cell lung cancer

Purpose

Protein tyrosine phosphatase receptor type T (PTPRT), which is a well‐known phosphatase and mutates frequently in melanoma and non‐small cell lung cancer (NSCLC). Our research aims to elucidate its mutation association with immune checkpoint inhibitors (ICI) efficacy.

Methods

We integrated whole‐exome sequencing (WES)‐based somatic mutation profiles and clinical characteristics of 631 melanoma samples received ICI agents from eight studies and 109 NSCLC samples from two studies. For validation, 321 melanoma and 350 NSCLC immunotherapy samples with targeted next‐generation sequencing (NGS) were employed. Besides, an independent NSCLC cohort contained 240 samples was also collected for further corroboration. Distinct immune infiltration was evaluated according to the PTPRT mutational status.

Results

In the WES melanoma cohort, patients with PTPRT mutations harbored a significantly elevated ICI response rate (40.5% vs. 28.6%, p = 0.036) and a prolonged survival outcome (35.3 vs. 24.9 months, p = 0.006). In the WES NSCLC cohort, the favorable response and immunotherapy survival were also observed in PTPRT‐mutated patients (p = 0.036 and 0.019, respectively). For the validation cohorts, the associations of PTRPT mutations with better prognoses were identified in melanoma, NSCLC, and pan‐cancer patients with targeted‐NGS (all p < 0.05). Moreover, immunology analyses showed the higher mutation burden, increased lymphocyte infiltration, decreased‐ activated‐stroma, and immune response pathways were detected in patients with PTPRT mutations.

Conclusion

Our investigation indicates that PTPRT mutations may be considered as a potential indicator for assessing ICI efficacy in melanoma and NSCLC, even across multiple cancers. Further prospective validation cohorts are warranted.

Mutational signature analysis in non-small cell lung cancer patients with a high tumor mutational burden

BACKGROUND

Lung cancer is the leading cause of cancer death worldwide. With the growing number of targeted therapies and the introduction of immuno-oncology (IO), personalized medicine has become standard of care in patients with metastatic disease. The development of predictive and prognostic biomarkers is of great importance. Mutational signatures harbor potential clinical value as predictors of therapy response in cancer. Here we set out to investigate particular mutational processes by assessing mutational signatures and associations with clinical features, tumor mutational burden (TMB) and targetable mutations.

METHODS

In this retrospective study, we studied tumor DNA from patients with non-small cell lung cancer (NSCLC) irrespective of stage. The samples were sequenced using a 2 megabase (Mb) gene panel. On each sample TMB was determined and defined as the total number of single nucleotide mutations per Mb (mut/Mb) including non-synonymous mutations. Mutational signature profiling was performed on tumor samples in which at least 30 somatic single base substitutions (SBS) were detected.

RESULTS

In total 195 samples were sequenced. Median total TMB was 10.3 mut/Mb (range 0-109.3). Mutational signatures were evaluated in 76 tumor samples (39%; median TMB 15.2 mut/Mb). SBS signature 4 (SBS4), associated with tobacco smoking, was prominently present in 25 of 76 samples (33%). SBS2 and/or SBS13, both associated with activity of the AID/APOBEC family of cytidine deaminases, were observed in 11 of 76 samples (14%). SBS4 was significantly more present in early stages (I and II) versus advanced stages (III and IV; P = .005).

CONCLUSION

In a large proportion of NSCLC patients tissue panel sequencing with a 2 Mb panel can be used to determine the mutational signatures. In general, mutational signature SBS4 was more often found in early versus advanced stages of NSCLC. Further studies are needed to determine the clinical utility of mutational signature analyses.

CALCR knockdown inhibits the development and progression of non-small-cell lung cancer

G-protein-coupled receptors (GPCRs) have been reported to participate in the occurrence and development of a variety of human cancers. CALCR is one of the hundreds of GPCRs, but its expression level and functional importance have never been investigated in non-small-cell lung cancer (NSCLC). In the present study, the protein expression level of CALCR was detected by immunohistochemical staining and western blot analysis. The Celigo cell counting assay was used to assess cell proliferation. Both the wound-healing assay and the transwell assay were performed to evaluate cell migration. Flow cytometric analysis was utilized to detect cell apoptosis and cell cycle. A mouse xenograft model was constructed to conduct the in vivo experiments. The results indicated that the CALCR expression was abundantly up-regulated in NSCLC and positively related to tumor infiltrate. Besides, CALCR knockdown could significantly suppress cell proliferation, migration, enhance apoptosis and arrest cell cycle. The in vivo study verified the inhibitory effects of CALCR knockdown on NSCLC tumorigenesis. The abovementioned results provided a reference for the treatment of NSCLC, that was, CALCR knockdown might be a considerable therapeutic strategy.

TP53 and CDKN2A mutations in patients with early-stage lung squamous cell carcinoma: an analysis of the correlations and prognostic outcomes

Background

Lung squamous cell carcinoma (LUSC) is characterized by frequent mutations of tumor protein p53 (TP53) and cyclin dependent kinase inhibitor 2A (CDKN2A). However, to date, the impact of TP53/CDKN2A status on the clinical outcome of patients with early-stage LUSC is unclear.

Methods

Tissue samples from 16 early-stage, surgically resected LUSCs were analyzed by next-generation sequencing (NGS). Information regarding TP53 and CDKN2A alterations and patient survival time was downloaded from The Cancer Genome Atlas (TCGA) database. The associations between TP53 and CDKN2A status and tumor characteristics, outcomes including overall survival (OS) and disease-free survival (DFS), and mutation counts were investigated.

Results

TP53 and CDKN2A exhibited a high frequency of somatic mutations in early-stage LUSC in our center. Data for 1,176 samples were collected from TCGA. CDKN2A mutation status was associated with TP53 mutation status (P=0.040). TP53 mutation was a favorable prognostic factor for early-stage LUSC. The OS times of patients with wild-type and mutated TP53 were 28.94 and 60.48 months, respectively (P=0.002). In contrast, CDKN2A mutations were significantly associated with a shorter survival time in early-stage LUSC. The OS times for wild-type and mutated CDKN2A patients were 62.81 and 37.55 months, respectively (P=0.026). Patients with TP53 mutations had higher total mutation counts compared to patients with wild-type TP53. Furthermore, OS was significantly shorter in patients with a low mutation count compared to patients with a median or high mutation count.

Conclusions

Early-stage LUSC patients with TP53 mutations had a longer OS, while those with CDKN2A mutations had a shorter OS. Furthermore, patients with TP53 mutation/CDKN2A wild-type status had a longer OS. CDKN2A mutation is a vital indicator for prognostic assessment according to TP53 status. The prolonged survival of patients with TP53 mutations may be due to their high mutation counts. Larger datasets are required to validate these observations.

Multiomics Differences in Lung Squamous Cell Carcinoma Patients with High Radiosensitivity Index Compared with Those with Low Radiosensitivity Index

Objectives

Radiosensitivity Index (RSI) can predict intrinsic radiotherapy sensitivity. We analyzed multiomics characteristics in lung squamous cell carcinoma between high and low RSI groups, which may help understand the underlying molecular mechanism of radiosensitivity and guide optional treatment for patients in the future.

Methods

The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) data were used to download clinical data, mRNA, microRNA, and lncRNA expression. Differential analyses, including mRNA, miRNA, lncRNA, and G.O. and KEGG, and GSVA analyses, were performed with R. Gene set enrichment analysis was done by GSEA. miRNA-differentially expressed gene network and ceRNA network were analyzed and graphed by the Cytoscape software.

Results

In TCGA data, 542 patients were obtained, including 171 in the low RSI group (LRSI) and 371 in the high RSI group (HRSI). In RNAseq, 558 significantly differentially expressed genes (DEGs) were obtained. KRT6A was the most significantly upregulated gene and IDO1 was the most significantly downregulated gene. In miRNAseq, miR-1269a was the most significantly upregulated. In lncRNAseq, LINC01871 was the most upregulated. A 66-pair interaction between differentially expressed genes and miRNAs and an 11-pair interaction between differential lncRNAs and miRNAs consisted of a ceRNA network, of which miR-184 and miR-490-3p were located in the center. In the GEO data, there were 40 DEGs. A total of 17 genes were founded in both databases, such as ADAM23, AHNAK2, BST2, COL11A1, CXCL13, FBN2, IFI27, IFI44L, MAGEA6, and PTGR1. GSVA analysis revealed 31 significant pathways. GSEA found 87 gene sets enriched in HRSI and 91 gene sets in LRSI. G.O. and KEGG of RNA expression levels revealed that these genes were most enriched in T cell activation and cytokine-cytokine receptor interaction.

Conclusions

Patients with lung squamous cell carcinoma have different multiomics characteristics between two groups. These differences may have an essential significance with radiotherapy effect.

Predictive Value of the TP53/PIK3CA/ATM Mutation Classifier for Patients With Bladder Cancer Responding to Immune Checkpoint Inhibitor Therapy

Background

Only a proportion of patients with bladder cancer may benefit from durable response to immune checkpoint inhibitor (ICI) therapy. More precise indicators of response to immunotherapy are warranted. Our study aimed to construct a more precise classifier for predicting the benefit of immune checkpoint inhibitor therapy.

Methods

This multi-cohort study examined the top 20 frequently mutated genes in five cohorts of patients with bladder cancer and developed the TP53/PIK3CA/ATM mutation classifier based on the MSKCC ICI cohort. The classifier was then validated in a validation set consisting of IMvigor210 cohort and Broad/Dana-Farber cohort. The molecular profile and immune infiltration characteristics in each subgroup as defined by this classifier were explored.

Results

Among all 881 patients with bladder cancer, the mutation frequency of TP53, PIK3CA, and ATM ranked in the top 20 mutated genes. The TP53/PIK3CA/ATM mutation classifier was constructed based on the Memorial Sloan Kettering Cancer Center (MSKCC) ICI cohort and only showed predictive value for patients with bladder cancer who received ICI therapy (median overall survival: low-risk group, not reached; moderate-risk group, 13.0 months; high-risk group, 8.0 months; P<0.0001). Similar results were found in subgroups of MSKCC ICI cohort defined by tumor mutation burden. Multivariate Cox analysis revealed that the risk group defined by the classifier served as an independent prognostic factor for overall survival in patients with bladder cancer. Efficacy of the classifier was verified in a validation set consisting of IMvigor210 cohort and Broad/Dana-Farber cohort. Lower expression of PD-1/PD-L1 and less tumor immune infiltration were observed in the high-risk group than the other two groups of the TCGA cohort and the IMvigor210 cohort.

Conclusion

Our study constructed a TP53/PIK3CA/ATM mutation classifier to predict the benefit of immune checkpoint inhibitor therapy for patients with bladder cancer. This classifier can potentially complement the tumor mutation burden and guide clinical ICI treatment decisions according to distinct risk levels.

Phase 1, pharmacogenomic, dose-expansion study of pegargiminase plus pemetrexed and cisplatin in patients with ASS1-deficient non-squamous non-small cell lung cancer

Introduction

We evaluated the arginine‐depleting enzyme pegargiminase (ADI‐PEG20; ADI) with pemetrexed (Pem) and cisplatin (Cis) (ADIPemCis) in ASS1‐deficient non‐squamous non‐small cell lung cancer (NSCLC) via a phase 1 dose‐expansion trial with exploratory biomarker analysis.

Methods

Sixty‐seven chemonaïve patients with advanced non‐squamous NSCLC were screened, enrolling 21 ASS1‐deficient subjects from March 2015 to July 2017 onto weekly pegargiminase (36 mg/m²) with Pem (500 mg/m²) and Cis (75 mg/m²), every 3 weeks (four cycles maximum), with maintenance Pem or pegargiminase. Safety, pharmacodynamics, immunogenicity, and efficacy were determined; molecular biomarkers were annotated by next‐generation sequencing and PD‐L1 immunohistochemistry.

Results

ADIPemCis was well‐tolerated. Plasma arginine and citrulline were differentially modulated; pegargiminase antibodies plateaued by week 10. The disease control rate was 85.7% (n = 18/21; 95% CI 63.7%–97%), with a partial response rate of 47.6% (n = 10/21; 95% CI 25.7%–70.2%). The median progression‐free and overall survivals were 4.2 (95% CI 2.9–4.8) and 7.2 (95% CI 5.1–18.4) months, respectively. Two PD‐L1‐expressing (≥1%) patients are alive following subsequent pembrolizumab immunotherapy (9.5%). Tumoral ASS1 deficiency enriched for p53 (64.7%) mutations, and numerically worse median overall survival as compared to ASS1‐proficient disease (10.2 months; n = 29). There was no apparent increase in KRAS mutations (35.3%) and PD‐L1 (<1%) expression (55.6%). Re‐expression of tumoral ASS1 was detected in one patient at progression (n = 1/3).

Conclusions

ADIPemCis was safe and highly active in patients with ASS1‐deficient non‐squamous NSCLC, however, survival was poor overall. ASS1 loss was co‐associated with p53 mutations. Therapies incorporating pegargiminase merit further evaluation in ASS1‐deficient and treatment‐refractory NSCLC.

PTPRT Could Be a Treatment Predictive and Prognostic Biomarker for Breast Cancer

The role of PTPRT in breast cancer was not comprehensively explored and well analyzed. Our study comprehensively searched available databases to analyze the clinical role of PTPRT in breast cancer. We found PTPRT was an antioncogene and could be used to distinguish different stages, age groups, molecular types, and grades for breast cancer. PTPRT might be primary resistance biomarkers for taxane, anthracycline, and ixabepilone but not be acquired resistance biomarkers. Higher PTPRT expression levels were associated with longer overall survival and recurrence-free survival. PTPRT was negatively associated with Ki67 and CDK4/6 but positively associated with BCL-2. PTPRT might be associated with cell cycle and microtubule, and tumor infiltration in B cell and macrophage cell. PTPRT could predict chemotherapy effectiveness and prognosis for breast cancer patients. PTPRT might inhibit tumor growth via disrupting the microtubule dynamics and cell cycle in breast cancer.

Identification and Validation of PIK3CA as a Marker Associated with Prognosis and Immune Infiltration in Renal Clear Cell Carcinoma

Background

Kidney renal clear cell carcinoma (KIRC) is the most prevalent renal malignancy. The therapeutic strategies for advanced KIRC are very few, with only sunitinib being widely approved. Mutations in the PIK3CA gene can affect tumor cell proliferation, metastasis, and patients' survival.

Methods

Bioinformatics analysis was performed to explore the expression and clinical significance of PIK3CA in KIRC. Moreover, qRT-PCR was conducted to verify the result.

Results

Subgroup analyses of KIRC tissue based on gender, tumor grade, and cancer stage indicated downregulation of PIK3CA mRNA expression. The KIRC patients with high PIK3CA expression indicated a better overall survival, progression-free survival, and disease-free survival. A predictive nomogram was constructed and demonstrated that the calibration plots for the 3-year and 5-year OS rates were predicted relatively well compared with an ideal model in the TCGA KIRC cohort. The validation study revealed that downregulation of PIK3CA in KIRC tissues and low PIK3CA expression had a poor overall survival with an AUC of 0.775 in the ROC curve. Moreover, Cox regression analysis revealed that PIK3CA expression and clinical stage were independent factors affecting the prognosis of KIRC patients. PIK3CA expression was found to be significantly associated with the abundance of immune cells and immune biomarker sets. PIK3CA and associated genes were found to be mainly associated with immune response and the JAK-STAT signaling pathway.

Conclusion

We identified PIK3CA as a potential biomarker for prognosis correlated with immune infiltrates in KIRC. Further studies should focus on the functions of PIK3CA in KIRC carcinogenesis.

Role of NRF2 in Lung Cancer

The gene expression program induced by NRF2 transcription factor plays a critical role in cell defense responses against a broad variety of cellular stresses, most importantly oxidative stress. NRF2 stability is fine-tuned regulated by KEAP1, which drives its degradation in the absence of oxidative stress. In the context of cancer, NRF2 cytoprotective functions were initially linked to anti-oncogenic properties. However, in the last few decades, growing evidence indicates that NRF2 acts as a tumor driver, inducing metastasis and resistance to chemotherapy. Constitutive activation of NRF2 has been found to be frequent in several tumors, including some lung cancer sub-types and it has been associated to the maintenance of a malignant cell phenotype. This apparently contradictory effect of the NRF2/KEAP1 signaling pathway in cancer (cell protection against cancer versus pro-tumoral properties) has generated a great controversy about its functions in this disease. In this review, we will describe the molecular mechanism regulating this signaling pathway in physiological conditions and summarize the most important findings related to the role of NRF2/KEAP1 in lung cancer. The focus will be placed on NRF2 activation mechanisms, the implication of those in lung cancer progression and current therapeutic strategies directed at blocking NRF2 action.

Using bioinformatics approaches to investigate driver genes and identify BCL7A as a prognostic gene in colorectal cancer

Colorectal cancer (CRC) results from the uncontrolled growth of cells in the colon, rectum, or appendix. The 5-year relative survival rate for patients with CRC is 65% and is correlated with the stage at diagnosis (being 91% for stage I at diagnosis versus 12% for stage IV). This study aimed to identify CRC driver genes to assist in the design of a cancer panel to detect gene mutations during clinical early-stage screening and identify genes for use in prognostic assessments and the evaluation of appropriate treatment options. First, we utilized bioinformatics approaches to analyze 354 paired sequencing profiles from The Cancer Genome Atlas (TCGA) to identify CRC driver genes and analyzed the sequencing profiles of 38 patients with >5 years of follow-up data to search for prognostic genes. The results revealed eight driver genes and ten prognostic genes. Next, the presence of the identified gene mutations was verified using tissue and blood samples from Taiwanese CRC patients. The results showed that the set identified gene mutations provide high coverage for driver gene screening, and APC, TP53, PIK3CA, and FAT4 could be detected in blood as ctDNA test targets. We further found that BCL7A gene mutation was correlated with prognosis in CRC (log-rank p-value = 0.02), and that mutations of BCL7A could be identified in ctDNA samples. These findings may be of value in clinical early cancer detection, disease monitoring, drug development, and treatment efforts in the future.

Identification of SLC38A7 as a Prognostic Marker and Potential Therapeutic Target of Lung Squamous Cell Carcinoma

BACKGROUND

No effective molecular targeted therapy has been established for SCC. We conducted a comprehensive study of SCC patients using RNA-sequencing and TCGA dataset to clarify the driver oncogene of SCC.

METHOD

Forty-six samples of 23 patients were totally analyzed with RNA-sequencing. We then searched for candidate-oncogenes of SCC using the TCGA database. To identify candidate oncogenes, we used the following 2 criteria: (1) the genes of interest were overexpressed in tumor tissues of SCC patients in comparison to normal tissues; and (2) using an integrated mRNA expression and DNA copy number profiling analysis using the TCGA dataset, the DNA copy number of the genes was positively correlated with the mRNA expression.

RESULT

We identified 188 candidate-oncogenes. Among those, the high expression of SLC38A7 was a strong prognostic marker that was significantly associated with a poor prognosis in terms of both overall survival (OS) and recurrence-free survival in the TCGA dataset (P < 0.05). Additionally, 202 resected SCC specimens were also subjected to an immunohistochemical analysis. Patients with the high expression of SLC38A7 (alternative name is sodium-coupled amino acid transporters 7) protein showed significantly shorter OS in comparison to those with the low expression of SLC38A7 protein [median OS 3.9 years (95% confidence interval, 2.4-6.4 years) vs 2.2 years (95% confidence interval, 1.9-4.1 years); log rank test: P = 0.0021].

CONCLUSION

SLC38A7, which is the primary lysosomal glutamine transporter required for the extracellular protein-dependent growth of cancer cells, was identified as a candidate therapeutic target of SCC.

EGFR-Mutated Squamous Cell Lung Cancer and Its Association With Outcomes

Background

The therapeutic efficacy of epidermal growth factor receptor tyrosine kinase inhibitors (EGFR-TKIs) in advanced EGFR-mutant lung squamous cell carcinoma (SCC) patients remains uncertain. Furthermore, the factors underlying the responsiveness have not been fully investigated. We therefore investigated the link between genomic profiles and EGFR-TKI efficacy.

Material and Methods

We consecutively enrolled stage IV, EGFR-mutant, and EGFR-TKI–treated patients with SCC. Patients with EGFR wild-type lung SCC and EGFR-mutant lung adenocarcinoma were consecutively enrolled as controls, and next-generation sequencing (NGS) was performed.

Results

In total, 28 EGFR-mutant lung SCC, 41 EGFR-mutant lung adenocarcinoma, and 40 EGFR wild-type lung SCC patients were included. Among the patients with EGFR mutations, shorter progression-free survival (PFS) was observed in SCC compared to adenocarcinoma (4.6 vs. 11.0 months, P<0.001). Comparison of the genomic profiles revealed that EGFR-mutant SCC patients had similar mutation characteristics to EGFR-mutant adenocarcinoma patients, but differed from those with EGFR wild-type SCC. Further exploration of EGFR-mutant SCC revealed that mutations in CREBBP (P = 0.005), ZNF217 (P = 0.016), and the Wnt (P = 0.027) pathway were negatively associated with PFS. Mutations in GRM8 (P = 0.025) were associated with improved PFS.

Conclusions

EGFR-mutant lung SCC has a worse prognosis than EGFR-mutant adenocarcinoma. Mutations in other genes, such as CREBBP, ZNF217, GRM8, or Wnt that had implications on PFS raise the possibility of understanding mechanisms of resistance to EGFR-TKI in lung SCC, which will aid identification of potential beneficial subgroups of patients with EGFR-mutant SCCs receiving EGFR-TKIs.

Female Gender Predicts Augmented Immune Infiltration in Lung Adenocarcinoma

INTRODUCTION

Immune infiltration in lung adenocarcinomas (LUADs) has been associated with response to immune checkpoint inhibitors. Clinical features underlying differential responses of patients with LUADs to immunotherapy are not well understood. Here, we analyzed the association between LUAD immune infiltration and clinicopathologic variables.

MATERIALS AND METHODS

Intratumoral CD3, CD8, and CD68 cell densities (tumor-associated immune cells [TAICs]) were immunohistochemically assessed in 146 surgically resected LUADs. LUADs were classified into 2 groups, low and high TAICs, based on the median values of cell densities for CD3, CD8, and CD68. Somatic mutation burden and driver gene mutation status were analyzed in a subset of the cases (n = 92). We statistically analyzed the association between the TAIC groups and various clinicopathologic and molecular variables by using the χ²/Fisher and Wilcoxon sum tests and multivariable logistic regression models.

RESULTS

Patient gender, tumor size, and STK11 mutations were significantly associated with TAIC levels in LUAD. Female patients exhibited significantly elevated TAIC levels (P = .005) compared with male patients. Tumor size was inversely associated with TAIC levels (P = .012). STK11 mutated tumors were associated with lower TAICs (P = .008). Higher TAICs were consistently observed in female patients with LUADs after adjusting for stage, tumor size, and age. Multivariable regression models confirmed female gender as an independent variable associated with TAIC levels in LUAD (P = .0141).

CONCLUSION

Immune infiltration in LUADs was significantly higher in female patients, warranting further exploration into the association between this clinical variable and immunotherapeutic response in LUAD.

Molecular characterization of lung squamous cell carcinoma tumors reveals therapeutically relevant alterations

INTRODUCTION

Unlike lung adenocarcinoma patients, there is no FDA-approved targeted-therapy likely to benefit lung squamous cell carcinoma patients.

MATERIALS AND METHODS

We performed survival analyses of lung squamous cell carcinoma patients harboring therapeutically relevant alterations identified by whole exome sequencing and mass spectrometry-based validation across 430 lung squamous tumors.

RESULTS

We report a mean of 11.6 mutations/Mb with a characteristic smoking signature along with mutations in TP53 (65%), CDKN2A (20%), NFE2L2 (20%), FAT1 (15%), KMT2C (15%), LRP1B (15%), FGFR1 (14%), PTEN (10%) and PREX2 (5%) among lung squamous cell carcinoma patients of Indian descent. In addition, therapeutically relevant EGFR mutations occur in 5.8% patients, significantly higher than as reported among Caucasians. In overall, our data suggests 13.5% lung squamous patients harboring druggable mutations have lower median overall survival, and 19% patients with a mutation in at least one gene, known to be associated with cancer, result in significantly shorter median overall survival compared to those without mutations.

CONCLUSIONS

We present the first comprehensive landscape of genetic alterations underlying Indian lung squamous cell carcinoma patients and identify EGFR, PIK3CA, KRAS and FGFR1 as potentially important therapeutic and prognostic target.

A retrospective study of clinicopathologic and molecular features of inoperable early-stage non-small cell lung cancer treated with stereotactic ablative radiotherapy

BACKGROUND/PURPOSE

Stereotactic ablative radiotherapy (SABR) is the treatment of choice for medically inoperable, early-stage non-small cell lung cancer (ES-NSCLC). The influence of oncogenic driver alterations and comorbidities are not well known. Here we present treatment outcomes based on clinicopathologic features and molecular profiles.

METHODS

We retrospectively analyzed patients treated with SABR for inoperable ES-NSCLC. Molecular features of oncogenic driver alterations included EGFR, ALK, and ROS1. Comorbidities were assessed using the age-adjusted Charlson Comorbidity Index (ACCI). Survival was calculated using the Kaplan-Meier method. The Cox regression model was performed for univariate and multivariate analyses of prognostic factors. Competing risk analysis was used to evaluate the cumulative incidence of disease progression.

RESULTS

From 2008 to 2020, 100 patients (median age: 82 years) were enrolled. The majority of patients were male (64%), ever-smokers (60%), and had adenocarcinoma (65%). With a median follow-up of 21.5 months, the median overall survival (OS) and real-world progression-free survival were 37.7 and 25.1 months, respectively. The competing-risk-adjusted 3-year cumulative incidences of local, regional, and disseminated failure were 8.2%, 14.5%, and 31.2%, respectively. An ACCI ≥7 was independently associated with inferior OS (hazard ratio [HR] 2.45, p = 0.03). Tumor size ≥4 cm (HR 4.16, p < 0.001) was the most important independent prognostic factor predicting real-world progression. EGFR mutation status had no impact on the outcomes.

CONCLUSION

SABR provides excellent local control in ES-NSCLC, although disseminated failures remains a major concern. ACCI is the best indicator for OS, while tumor sizes ≥4 cm predicts poor disease control.

A novel immune classification reveals distinct immune escape mechanism and genomic alterations: implications for immunotherapy in hepatocellular carcinoma

Background

The tumor immunological microenvironment (TIME) has a prominent impact on prognosis and immunotherapy. However, the heterogeneous TIME and the mechanisms by which TIME affects immunotherapy have not been elucidated in hepatocellular carcinoma (HCC).

Methods

A total of 2195 eligible HCC patients from TCGA and GEO database were collected. We comprehensively explored the different heterogeneous TIME phenotypes and its clinical significance. The potential immune escape mechanisms and what genomic alterations may drive the formation of different phenotypes were further investigated.

Results

We identified three phenotypes in HCC: TIME-1, the “immune-deficiency” phenotype, with immune cell depletion and proliferation; TIME-2, the “immune-suppressed” phenotype, with enrichment of immunosuppressive cells; TIME-3, the “immune-activated phenotype”, with abundant leukocytes infiltration and immune activation. The prognosis and sensitivity to both sorafenib and immunotherapy differed among the three phenotypes. We also underlined the potential immune escape mechanisms: lack of leukocytes and defective tumor antigen presentation capacity in TIME-1, increased immunosuppressive cells in TIME-2, and rich in immunoinhibitory molecules in TIME-3. The different phenotypes also demonstrated specific genomic events: TIME-1 characterized by TP53, CDKN2A, CTNNB1, AXIN1 and FOXD4 alterations; TIME-2 characterized by significant alteration patterns in the PI3K pathway; TIME-3 characterized by ARID1A mutation. Besides, the TIME index (TI) was proposed to quantify TIME infiltration pattern, and it was a superior prognostic and immunotherapy predictor. A pipeline was developed to classify single patient into one of these three subtypes and calculated the TI.

Conclusions

We identified three TIME phenotypes with different clinical outcomes, immune escape mechanisms and genomic alterations in HCC, which could present strategies for improving the efficacy of immunotherapy. TI as a novel prognostic and immunotherapeutic signature that could guide personalized immunotherapy and clinical management of HCC.

The underlying molecular mechanism and drugs for treatment in adrenal cortical carcinoma

Purpose: The study aimed to predict and explore the possible clinical value and mechanism of genetic markers in adrenal cortical carcinoma using a bioinformatics analysis method. Methods: The RNA-seqs and miRNAs data were downloaded from TCGA database to identify the differentially expressed genes and differentially expressed miRNAs. The hub-genes were screened by building protein-protein interaction sub-networks with 12 topological analysis methods. We conducted the receiver operating characteristic curve to elevate the diagnostic value of hub-genes in distinguishing the death and alive groups. The survival analysis of hub-genes and key miRNAs were conducted using Kaplan-Meier curves. Furthermore, most significant small molecules were identified as therapeutic candidates for adrenal cortical carcinoma by the CMap analysis. Results: Compared to survival group, we found 475 up-regulated genes and 354 genes and the key pathways leading to the death of different ACC individual patients. Then we used 12 topological analysis methods to found the most possible 22 hub-genes. Among these hub-genes, nine hub-genes (C3, CXCL5, CX3CR1, GRM8, HCAR2, HTR1B, SUCNR1, PTGER3 and SSTR1) could be used to distinguish the death and survival groups for patients. We also revealed that mRNA expressions of 12 genes (C3, CXCL8, CX3CR1, GNAT3, GNGT1, GRM8, HCAR2, HTR1B, HTR1D, PTGER3, SSTR1 and SUCNR1) and four key miRNAs (hsa-mir-330, hsa-mir-489, hsa-mir-508 and hsa-mir-513b) were related to survival. Three most small molecules were identified (H-9, AZ-628 and phensuximide) as potential therapeutic drugs for adrenal cortical carcinoma. Conclusion: The hub-genes expression was significant useful in adrenal cortical carcinoma, provide new diagnostic, prognosis and therapeutic approaches for adrenal cortical carcinoma. Furthermore, we also explore the possible miRNAs involved in regulation of hub-genes.

m6A regulator-based methylation modification patterns characterized by distinct tumor microenvironment immune profiles in colon cancer

Recent studies have highlighted the biological significance of RNA N⁶-methyladenosine (m⁶A) modification in tumorigenicity and progression. However, it remains unclear whether m⁶A modifications also have potential roles in immune regulation and tumor microenvironment (TME) formation. Methods: In this study, we curated 23 m⁶A regulators and performed consensus molecular subtyping with NMF algorithm to determine m⁶A modification patterns and the m⁶A-related gene signature in colon cancer (CC). The ssGSEA and CIBERSORT algorithms were employed to quantify the relative infiltration levels of various immune cell subsets. An PCA algorithm based m⁶Sig scoring scheme was used to evaluate the m⁶A modification patterns of individual tumors with an immune response. Results: Three distinct m6A modification patterns were identified among 1307 CC samples, which were also associated with different clinical outcomes and biological pathways. The TME characterization revealed that the identified m⁶A patterns were highly consistent with three known immune profiles: immune-inflamed, immune-excluded, and immune-desert, respectively. Based on the m⁶Sig score, which was extracted from the m⁶A-related signature genes, CC patients can be divided into high and low score subgroups. Patients with lower m⁶Sig score was characterized by prolonged survival time and enhanced immune infiltration. Further analysis indicated that lower m⁶Sig score also correlated with greater tumor mutation loads, PD-L1 expression, and higher mutation rates in SMGs (e.g., PIK3CA and SMAD4). In addition, patients with lower m⁶Sig scores showed a better immune responses and durable clinical benefits in three independent immunotherapy cohorts. Conclusions: This study highlights that m⁶A modification is significantly associated with TME diversity and complexity. Quantitatively evaluating the m⁶A modification patterns of individual tumors will strengthen our understanding of TME characteristics and promote more effective immunotherapy strategies.

The NRF2, Thioredoxin, and Glutathione System in Tumorigenesis and Anticancer Therapies

Cancer remains an elusive, highly complex disease and a global burden. Constant change by acquired mutations and metabolic reprogramming contribute to the high inter- and intratumor heterogeneity of malignant cells, their selective growth advantage, and their resistance to anticancer therapies. In the modern era of integrative biomedicine, realizing that a personalized approach could benefit therapy treatments and patients' prognosis, we should focus on cancer-driving advantageous modifications. Namely, reactive oxygen species (ROS), known to act as regulators of cellular metabolism and growth, exhibit both negative and positive activities, as do antioxidants with potential anticancer effects. Such complexity of oxidative homeostasis is sometimes overseen in the case of studies evaluating the effects of potential anticancer antioxidants. While cancer cells often produce more ROS due to their increased growth-favoring demands, numerous conventional anticancer therapies exploit this feature to ensure selective cancer cell death triggered by excessive ROS levels, also causing serious side effects. The activation of the cellular NRF2 (nuclear factor erythroid 2 like 2) pathway and induction of cytoprotective genes accompanies an increase in ROS levels. A plethora of specific targets, including those involved in thioredoxin (TRX) and glutathione (GSH) systems, are activated by NRF2. In this paper, we briefly review preclinical research findings on the interrelated roles of the NRF2 pathway and TRX and GSH systems, with focus given to clinical findings and their relevance in carcinogenesis and anticancer treatments.

Functional network analysis reveals potential repurposing of β-blocker atenolol for pancreatic cancer therapy

BACKGROUND

The survival rate of patients with pancreatic cancer is low; therefore, continuous discovery and development of novel pancreatic cancer drugs are required. Functional network analysis is an integrated bioinformatics approach based on gene, target, and disease networks interaction, and it is extensively used in drug discovery and development.

OBJECTIVE

This study aimed to identify if atenolol, a selective adrenergic inhibitor, can be repurposed for the treatment of pancreatic cancer using functional network analysis.

METHODS

Direct target proteins (DTPs) and indirect target proteins (ITPs) were obtained from STITCH and STRING databases, respectively. Atenolol-mediated proteins (AMPs) were collected from DTPs and ITPs and further analyzed for gene ontology, KEGG pathway enrichment, genetic alterations, overall survival, and molecular docking.

RESULTS

We obtained 176 AMPs that consisted of 10 DTPs and 166 ITPs. Among the AMPs involved in the pancreatic cancer pathways, several AMPs such as MAPK1, RELA, MAPK8, STAT1, and STAT3 were identified. Genetic alterations in seven AMPs were identified in 0.9%-16% of patients. Patients with high mRNA levels of MAPK1, RELA, STAT3, GNB1, and MMP9 had significantly worse overall survival rates compared with patients with low expression. Molecular docking studies showed that RELA and MMP9 are potential target candidates of atenolol in the treatment of patients with pancreatic cancer.

CONCLUSION

In conclusion, atenolol can potentially be repurposed to target pancreatic cancer cells by modulating MMP9 and NF-κB signaling. The results of this study need to be further validated in vitro and in vivo.

Advanced Squamous Cell Carcinoma of the Lung: Current Treatment Approaches and the Role of Afatinib

Options for the treatment of squamous cell lung carcinoma expanded in recent years with the introduction of the immune checkpoint inhibitors into routine clinical practice in both the first- and second-line settings but are still limited. As a result, pembrolizumab, given either alone or in combination with platinum-based chemotherapy, is now a standard first-line treatment for squamous cell lung cancer. However, few options exist once patients have progressed on immune checkpoint inhibitors and chemotherapy. In this setting, the irreversible ErbB family blocker, afatinib, has a potential role as second or subsequent therapy for some patients. The Phase III LUX-Lung 8 study demonstrated that afatinib significantly prolonged progression-free and overall survival compared with erlotinib in patients with squamous cell lung carcinoma. Notably, retrospective, ad-hoc biomarker analyses of a subset of patients from LUX-Lung 8 suggested that patients with ErbB family mutations derived particular benefit from afatinib, especially those with ErbB2 (HER2) mutations. Afatinib has a manageable and predictable safety profile, and adverse events can be managed with the use of a tolerability-guided dose modification protocol. Until more data are available, afatinib could be considered as a potential second-line treatment option for patients who have progressed on combined pembrolizumab and platinum-based chemotherapy and are ineligible for more established second-line options, or as a third-line option in patients who have received first-line immunotherapy, and second-line chemotherapy or chemotherapy and antiangiogenesis therapy. However, further data are required to support the use of afatinib following immunotherapy. Given that treatment options are limited in both of these settings, investigating an agent with an entirely new mechanism of action is warranted. If available, molecular analysis to identify ErbB family mutations or the use of proteomic profiling could help to further isolate patients who are likely to derive the most benefit from afatinib.

Clinical features, treatment, and survival outcome of primary pulmonary NUT midline carcinoma

Objective

NUT midline carcinoma (NMC), a rare type of squamous cell carcinoma, is genetically characterised by NUT midline carcinoma family member 1 (NUTM1) gene rearrangement. NMC can arise from the lungs; however, there is no standard for the management of primary pulmonary NMC. This study aimed to confirm the clinical features and report the treatments, especially with immune checkpoint inhibitors (ICIs), and outcomes of patients with primary pulmonary NMC.

Methods

A retrospective review of patients with primary pulmonary NMC was performed in the First Affiliated Hospital of Guangzhou Medical University between January 2015 and December 2018. Clinical manifestations as well as radiographic and pathological findings were recorded. Whole-exome sequencing (WES), a predictor for ICI response, was used to determine the tumour mutational burden (TMB). Treatments, especially by immune checkpoint blockade, and patient survival were analysed.

Results

Seven patients with primary pulmonary mass (four men and three women) with a mean age of 42 years (range, 23–74) who were diagnosed with NMC according to NUT immunohistochemistry staining were included for analysis. One patient had a rare fusion of CHRM5-NUTM1 by tumour sequencing. A wide range of TMB (1.75–73.81 mutations/Mbp) was observed. The initial treatments included chemotherapy (5/7, 71.4%), surgery (1/7, 14.3%), and radiotherapy (1/7, 14.3%). Five patients (5/7, 71.4%) received ICIs (programmed cell death protein 1 [PD1]/programmed cell death ligand 1 [PDL1] monoclonal antibody) as second- or higher-line treatments. The median overall survival (OS) was 4.1 months (range, 1.5–26.7 months).

Conclusions

Patients with primary pulmonary NMC have a poor prognosis and chemotherapy is often preferred. Checkpoint immunotherapy is a good option as the second- or higher-line treatment. TMB seems to be not associated with OS.

Whole genome sequencing analysis identifies recurrent structural alterations in esophageal squamous cell carcinoma

Esophageal squamous cell carcinoma (ESCC) is the predominant type of esophageal cancer in the Asian region, including Japan. A previous study reported mutational landscape of Japanese ESCCs by using exome sequencing. However, somatic structural alterations were yet to be explored. To provide a comprehensive mutational landscape, we performed whole genome sequencing (WGS) analysis of biopsy specimens from 20 ESCC patients in a Japanese population. WGS analysis identified non-silent coding mutations of TP53, ZNF750 and FAT1 in ESCC. We detected six mutational signatures in ESCC, one of which showed significant association with smoking status. Recurrent structural variations, many of which were chromosomal deletions, affected genes such as LRP1B, TTC28, CSMD1, PDE4D, SDK1 and WWOX in 25%–30% of tumors. Somatic copy number amplifications at 11q13.3 (CCND1), 3q26.33 (TP63/SOX2), and 8p11.23 (FGFR1) and deletions at 9p21.3 (CDKN2A) were identified. Overall, these multi-dimensional view of genomic alterations improve the understanding of the ESCC development at molecular level and provides future prognosis and therapeutic implications for ESCC in Japan.

Clinical and genomic characterization of neutral tumor evolution in Head and Neck Squamous Cell Carcinoma

Recent studies suggest that a significant proportion of cancers undergo neutral tumor evolution. We applied neutral evolution model in HNSCC patients from The Cancer Genome Atlas (TCGA). To ensure the accuracy of classification results, a sample with the purity of tumor <0.7 was excluded. A tumor sample was considered to evolve neutrally if R² ≥ 0.98. We found that about 16% of HNSCC patients undergo neutral tumor evolution. We showed that neutral evolution HNSCC patients have better prognosis and higher activities of immune response pathways, and the numbers of co-occurring mutation events and significantly positive selection mutations are significantly less than non-neutral evolution HNSCC patients. In conclusion, we described a comprehensive clinical and genomic characteristics of neutral tumor evolution in Head and Neck Squamous Cell Carcinoma (HNSCC), and provided evidence that the evolution history of HNSCC has both clinical and biological implications.