Clinicopathological, microenvironmental and genetic determinants of molecular subtypes in KEAP1/NRF2-mutant lung cancer

Further knowledge and developments in resistance mechanisms to immune checkpoint inhibitors

The past decade has witnessed a revolution in cancer treatment, shifting from conventional drugs (chemotherapies) towards targeted molecular therapies and immune-based therapies, in particular immune-checkpoint inhibitors (ICIs). These immunotherapies release the host's immune system against the tumor and have shown unprecedented durable remission for patients with cancers that were thought incurable, such as metastatic melanoma, metastatic renal cell carcinoma (RCC), microsatellite instability (MSI) high colorectal cancer and late stages of non-small cell lung cancer (NSCLC). However, about 80% of the patients fail to respond to these immunotherapies and are therefore left with other less effective and potentially toxic treatments. Identifying and understanding the mechanisms that enable cancerous cells to adapt to and eventually overcome therapy can help circumvent resistance and improve treatment. In this review, we describe the recent discoveries on the onco-immunological processes which govern the tumor microenvironment and their impact on the resistance to PD-1/PD-L1 checkpoint blockade.

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.

Glutathione peroxidase 2 is a metabolic driver of the tumor immune microenvironment and immune checkpoint inhibitor response

BACKGROUND

The existence of immunologically 'cold tumors' frequently found across a wide spectrum of tumor types represents a significant challenge for cancer immunotherapy. Cold tumors have poor baseline pan-leukocyte infiltration, including a low prevalence of cytotoxic lymphocytes, and not surprisingly respond unfavorably to immune checkpoint (IC) inhibitors. We hypothesized that cold tumors harbor a mechanism of immune escape upstream and independent of ICs that may be driven by tumor biology rather than differences in mutational neoantigen burden.

METHODS

Using a bioinformatic approach to analyze TCGA (The Cancer Genome Atlas) RNA sequencing data we identified genes upregulated in cold versus hot tumors across four different smoking-related cancers, including squamous carcinomas from the oral cavity (OCSCC) and lung (LUSC), and adenocarcinomas of the bladder (BLCA) and lung (LUAD). Biological significance of the gene most robustly associated with a cold tumor phenotype across all four tumor types, glutathione peroxidase 2 (GPX2), was further evaluated using a combination of in silico analyses and functional genomic experiments performed both in vitro and in in vivo with preclinical models of oral cancer.

RESULTS

Elevated RNA expression of five metabolic enzymes including GPX2, aldo-keto reductase family 1 members AKR1C1, AKR1C3, and cytochrome monoxygenases (CP4F11 and CYP4F3) co-occurred in cold tumors across all four smoking-related cancers. These genes have all been linked to negative regulation of arachidonic acid metabolism-a well-established inflammatory pathway-and are also known downstream targets of the redox sensitive Nrf2 transcription factor pathway. In OCSCC, LUSC, and LUAD, GPX2 expression was highly correlated with Nrf2 activation signatures, also elevated in cold tumors. In BLCA, however, GPX2 correlated more strongly than Nrf2 signatures with decreased infiltration of multiple leukocyte subtypes. GPX2 inversely correlated with expression of multiple pro- inflammatory cytokines/chemokines and NF-kB activation in cell lines and knockdown of GPX2 led to increased secretion of prostaglandin E2 (PGE2) and interleukin-6. Conversely, GPX2 overexpression led to reduced PGE2 production in a murine OCSCC model (MOC1). GPX2 overexpressing MOC1 tumors had a more suppressive tumor immune microenvironment and responded less favorably to anti-cytotoxic T-lymphocytes-associated protein 4 IC therapy in mice.

CONCLUSION

GPX2 overexpression represents a novel potentially targetable effector of immune escape in cold tumors.

Genomic landscape of chemical-induced lung tumors under Nrf2 different expression levels

The transcription factor Nrf2 plays a crucial role in the anti-oxidative stress response, protection of DNA from injury, and DNA repair mechanisms. Nrf2 activity reduces cancer initiation, but how Nrf2 affects whole-genome alterations upon carcinogenic stimulus remains unexplored. Although recent genome-wide analysis using next-generation sequencing revealed landscapes of nucleotide mutations and copy number alterations in various human cancers, genomic changes in murine cancer models have not been thoroughly examined. We elucidated the relationship between Nrf2 expression levels and whole exon mutation patterns using an ethyl-carbamate (urethane)-induced lung carcinogenesis model employing Nrf2-deficient and Keap1-kd mice, the latter of which express high levels of Nrf2. Exome analysis demonstrated that single nucleotide and trinucleotide mutation patterns and the Kras mutational signature differed significantly and were dependent on the expression level of Nrf2. The Nrf2-deficient tumors exhibited fewer copy number alterations relative to the Nrf2-wt and Keap1-kd tumors. The observed trend in genomic alterations likely prevented the Nrf2-deficient tumors from progressing into malignancy. For the first time, we present whole-exome sequencing results for chemically-induced lung tumors in the Nrf2 gain or loss of function mouse models. Our results demonstrate that different Nrf2 expression levels lead to distinct gene mutation patterns that underly different oncogenic mechanisms in each tumor genotype.

Molecular Pathology of Lung Cancer

This overview of the molecular pathology of lung cancer includes a review of the most salient molecular alterations of the genome, transcriptome, and the epigenome. The insights provided by the growing use of next-generation sequencing (NGS) in lung cancer will be discussed, and interrelated concepts such as intertumor heterogeneity, intratumor heterogeneity, tumor mutational burden, and the advent of liquid biopsy will be explored. Moreover, this work describes how the evolving field of molecular pathology refines the understanding of different histologic phenotypes of non-small-cell lung cancer (NSCLC) and the underlying biology of small-cell lung cancer. This review will provide an appreciation for how ongoing scientific findings and technologic advances in molecular pathology are crucial for development of biomarkers, therapeutic agents, clinical trials, and ultimately improved patient care.

Epithelial-Mesenchymal Transition Induces GSDME Transcriptional Activation for Inflammatory Pyroptosis

GSDME is a newly recognized executor of cellular pyroptosis, and has been recently implicated in tumor growth and immunity. However, knowledge about the molecular regulators underlying GSDME abundance remains limited. Here, we performed integrative bioinformatics analyses and identified that epithelial-mesenchymal transition (EMT) gene signatures exhibited positive correlation with GSDME levels across human cancers. A causal role was supported by the observation that EMT dictated GSDME reversible upregulation in multiple experimental models. Mechanistically, transcriptional activation of GSDME was directly driven by core EMT-activating transcription factors ZEB1/2, which bound to the GSDME promoter region. Of functional importance, elevated GSDME in mesenchymally transdifferentiated derivatives underwent proteolytic cleavage upon antineoplastic drug exposure, leading to pyroptotic cell death and consequent cytokine release. Taken together, our findings pinpointed a key transcriptional machinery controlling GSDME expression and indicated potential therapeutic avenues to exploit GSDME-mediated inflammatory pyroptosis for the treatment of mesenchymal malignancies.

Thioredoxin reductase-1 levels are associated with NRF2 pathway activation and tumor recurrence in non-small cell lung cancer

Activating mutations in the KEAP1/NRF2 pathway characterize a subset of non-small cell lung cancer (NSCLC) associated with chemoresistance and poor prognosis. We herein evaluated the relationship between 64 oxidative stress-related genes and overall survival data from 35 lung cancer datasets. Thioredoxin reductase-1 (TXNRD1) stood out as the most significant predictor of poor outcome. In a cohort of NSCLC patients, high TXNRD1 protein levels correlated with shorter disease-free survival and distal metastasis-free survival post-surgery, including a subset of individuals treated with platinum-based adjuvant chemotherapy. Bioinformatics analysis revealed that NSCLC tumors harboring genetic alterations in the NRF2 pathway (KEAP1, NFE2L2 and CUL3 mutations, and NFE2L2 amplification) overexpress TXNRD1, while no association with EGFR, KRAS, TP53 and PIK3CA mutations was found. In addition, nuclear accumulation of NRF2 overlapped with upregulated TXNRD1 protein in NSCLC tumors. Functional cell assays and gene dependency analysis revealed that NRF2, but not TXNRD1, has a pivotal role in KEAP1 mutant cells' survival. KEAP1 mutants overexpress TXNRD1 and are less susceptible to the cytotoxic effects of the TXNRD1 inhibitor auranofin when compared to wild-type cell lines. Inhibition of NRF2 with siRNA or ML-385, and glutathione depletion with buthionine-sulfoximine, sensitized KEAP1 mutant A549 cells to auranofin. NRF2 knockdown and GSH depletion also augmented cisplatin cytotoxicity in A549 cells, whereas auranofin had no effect. In summary, these findings suggest that TXNRD1 is not a key determinant of malignant phenotypes in KEAP1 mutant cells, although this protein can be a surrogate marker of NRF2 pathway activation, predicting tumor recurrence and possibly other aggressive phenotypes associated with NRF2 hyperactivation in NSCLC.

Integrated analysis of patients with KEAP1/NFE2L2/CUL3 mutations in lung adenocarcinomas

Objectives

To explore the clinical features, molecular characteristics, and immune landscape of lung adenocarcinoma patients with KEAP1/NFE2L2/CUL3 mutations.

Methods

The multi‐omics data from the GDC‐TCGA LUAD project of The Cancer Genome Atlas (TCGA) database were downloaded from the Xena browser. The estimate of the immune infiltration was implemented by using the GSVA analysis and CIBERSORT. The status of KEAP1/NFE2L2/CUL3 mutation in 50 LUAD samples of our department was detected by using Sanger sequencing, following the relative expression level of differentially expressed genes (DEGs), miRNAs (DEmiRNAs), and lncRNAs (DElncRNAs) was validated by IHC and real‐time quantitative polymerase chain reaction (RT‐qPCR).

Results

The Kaplan–Meier and multivariable Cox regression analyses demonstrated that KEAP1/NFE2L2/CUL3 mutations had independent prognostic value for OS and PFS in LUAD patients. The differential analysis detected 207 upregulated genes (like GSR/UGT1A6) and 447 downregulated genes (such as PIGR). GO, KEGG, and GSEA analyses demonstrated that DEGs were enriched in glutamate metabolism and the immune response. The constructed ceRNA network shows the linkage of differential lncRNAs and mRNAs. Three hundred and nine somatic mutations were detected, alterations in immune infiltration DNA methylations and stemness scores were also founded between the two groups. Eight mutated LUAD patients were detected by Sanger DNA sequencing in 50 surgical patients. GSR and UGT1A6 were validated to express higher in the Mut group, whereas the expression of PIGR was restrained. Furthermore, the IHC staining conducted on paraffin‐embedded tissue emphasizes the consistency of our result.

Conclusion

This research implemented the comprehensive analysis of KEAP1/NFE2L2/CUL3 somatic mutations in the LUAD patients. Compared with the wild type of LUAD patients, the Mut group shows a large difference in clinical features, RNA sequence, DNA methylation, and immune infiltrations, indicating complex mechanism oncogenesis and also reveals potential therapeutic targets.

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.

KEAP1/NFE2L2 Mutations of Liquid Biopsy as Prognostic Biomarkers in Patients With Advanced Non-Small Cell Lung Cancer: Results From Two Multicenter, Randomized Clinical Trials

PURPOSE

The KEAP1-NFE2L2 (Kelch-like ECH-associated protein 1 (KEAP1)-Nuclear factor (erythroid-derived 2)-like 2 (NFE2L2)) mutations are associated with resistance to chemotherapy or immunotherapy in non-small cell lung cancer (NSCLC). Conversely, it has been reported that NFE2L2 mutations potentiate improved clinical outcome with immunotherapy. However, therapeutic benefits for patients with KEAP1/NFE2L2 mutations remain unclear. The purpose of this study was to investigate the association between KEAP1/NFE2L2 and NSCLC prognosis, and to explore whether immunotherapy can improve prognosis in populations with KEAP1/NFE2L2 mutations.

EXPERIMENTAL DESIGN

The impact of KEAP1/NFE2L2 mutations on survival outcomes in NSCLC patients received immunotherapy and chemotherapy was verified in the randomized phase II/III POPLAR/OAK trials (blood-based sequencing, bNGS cohort, POPLAR (n = 211) and OAK (n = 642)). The Cancer Genome Atlas (TCGA) NSCLC cohort (n=998) and an in-house Chinese NSCLC cohort (n=733) was used For the analysis of immune-related markers.

RESULTS

Compared with KEAP1/NFE2L2 wild-type, patients with KEAP1/NFE2L2 mutations were significantly associated with poorer overall survival (OS, HR = 1.97, 95% CI 1.48-2.63, P < 0.001) on atezolizumab and docetaxel (HR = 1.66, 95% CI 1.28-2.16, P < 0.001). In KEAP1/NFE2L2 mutant group, there was no significant difference in median OS between atezolizumab and docetaxel (HR 0.74, 95% CI 0.53-1.03, P = 0.07). NFE2L2/KEAP1 mutations were significantly associated with higher TMB values and PD-L1 expression in the OAK/POPLAR and in-house Chinese NSCLC cohorts. GSEA revealed that KEAP1/NFE2L2mutant subgroup was associated with deficient infiltration of CD4+ T cells, NK T cells and natural Treg cells, and lower expression of DNA damage response genes in TCGA NSCLC cohort.

CONCLUSIONS

Our study revealed that patients with KEAP1/NFE2L2 mutations have a worse prognosis than wild-type patients, both on immunotherapy and chemotherapy. In addition, in patients with KEAP1/NFE2L2 mutations, immunotherapy did not significantly improve prognosis compared to chemotherapy.

Prognostic and clinicopathological significance of NRF2 expression in non-small cell lung cancer: A meta-analysis

Nuclear factor erythroid 2-related factor 2 (NRF2) functions as a transcription factor and regulates a wide array of antioxidant and stress-responsive genes. NRF2 has been widely implicated in different types of cancers, but only limited studies concerning the relationship between NRF2 expression and tumour invasion or prognosis in lung cancer. Therefore, we conducted a meta-analysis to determine the prognostic value of NRF2 in patients with non-small cell lung cancer (NSCLC). The relationship between NRF2 expression in NSCLC patients and clinicopathological features was also investigated. Overall survival (OS) and treatment response rate were evaluated using STATA software. Twenty eligible articles with 2530 lung cancer patients were included in this meta-analysis. The results revealed that high expression level of NRF2 was associated with pathologic distant metastasis (odds ratio (OR) = 2.64, 95% confidence interval (CI) 1.62-4.31; P < 0.001), lymph node metastasis (OR = 2.14, 95% CI: 1.53-3.00; P < 0.001), and tumour node metastasis (TNM) stage (OR = 1.95, 95% CI: 1.52-2.49, P < 0.001). High NRF2 expression was associated with low treatment response rate in platinum-based chemotherapy (HR = 0.11, 95% CI 0.02-0.51; P = 0.005). High expression level of NRF2 is predictive for poor overall survival rate (HR = 1.86, 95% CI 1.44-2.41, P < 0.001) and poor progression-free survival (PFS) (HR = 2.27, 95% CI 1.26-4.09, P = 0.006). Compared to patients with a low level of NRF2 expression, patients with high NRF2 expression levels were associated with worse OS and PFS when given the chemotherapy or EGFR-TKI. Together, our meta-analysis results suggest that NRF2 can act as a potential indicator of NSCLC tumour aggressiveness and help the prognosis and design of a better treatment strategy for NSCLC patients.

Identification, Characterization, and Mutational Analysis of a Probable KEAP1 Ortholog in Rice (Oryza sativa L.)

The Kelch-like ECH-associated protein 1 (KEAP1)-nuclear factor E2-related factor 2 (NRF2) module is a key component in the detoxification and antioxidant system in animals, which plays crucial roles in cell homeostasis and cytoprotection, and consequently in carcinogenesis and disease development. However, this system seems to have diverged throughout evolution across different organisms, and the question of whether a similar system exists in plants has thus far remained unresolved. In this study, a KEAP1 ortholog was identified in rice (Oryza sativa L., OsKEAP1) and its properties were characterized via in silico and laboratory studies. To reveal OsKEAP1's function, two knockdown mutants, oskeap1-1 and oskeap1-2, were generated by targeted mutagenesis in the 5' untranslated region (UTR) using the CRISPR-Cas9 system. In silico analysis showed that OsKEAP1 has a Kelch-repeat domain which is identical to those of animals and a plant-specific development and cell death (DCD) domain in place of the broad-complex, tramtrack, bric-a-brac (BTB) domain found in animals. Orthologs of OsKEAP1 are present across plant species and all have the DCD domain and the Kelch-repeat domain. OsKEAP1 was proven to be localized to both the cytoplasm and nucleus, in contrast to the exclusive cytoplasm localization of animal KEAP1. Single nucleotide insertions in the 5' UTR significantly reduced the transcription level of OsKEAP1 in the oskeap1-1 and oskeap1-2 mutants. The oskeap1 mutations greatly impaired plant growth and development, resulting in significant declines in a majority of agronomic and yield-related traits, i.e., plant height, panicle length, grain number per plant, and seed-set rate. The downregulation of OsKEAP1 increased the levels of H2O2, malondialdehyde, and proline while significantly decreasing the expression of two catalase genes in seedlings grown under normal and salt-stressed conditions. The changes in the above phenotypes are either positively or negatively correlated with the degree of OsKEAP1 downregulation. Altogether, we identified a probable KEAP1 ortholog in rice, revealed its unique subcellular localization, and demonstrated its important functions in vegetative and reproductive growth via regulation of the antioxidant response in plants.

T-cell exhaustion interrelates with immune cytolytic activity to shape the inflamed tumor microenvironment

Direct quantification of exhausted T cells in human cancer is lacking, and its predictive value for checkpoint-based treatment remains poorly investigated. We sought to systematically characterize the pan-cancer landscape and molecular hallmarks of T-cell dysfunction for the purpose of precision immunotherapy. Here, we defined a transcriptional signature for T-cell exhaustion through analyzing differential gene expression between PD-1-high and PD-1-negative CD8⁺ T lymphocytes from primary non-small cell lung cancer (NSCLC), followed by positive correlation tests with PDCD1 in TCGA lung carcinomas. A 78-gene signature for exhausted CD8⁺ T cells (GET) was identified and validated to reflect dysfunctional immune state spanning different species and disease models. We discovered that GET estimation significantly correlated with intratumoral immune cytolytic activity (CYT) and T-cell-inflamed gene expression profile (GEP) across 30 solid tumor types. Miscellaneous tumor-intrinsic and -extrinsic properties, in particular leukocyte proportions, genomic abnormalities, specific mutational signatures, and signaling pathways, were notably associated with GET levels. Furthermore, higher GET expression predicted an increased likelihood of clinical response to immune checkpoint inhibitors. These findings highlight the interrelation between T-cell exhaustion and immune cytolytic activity at the pan-cancer scale. The resulting inflamed tumor microenvironment may further crosstalk with other molecular and clinicopathological factors, which should be properly considered during immunotherapy biomarker development. © 2020 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Identification of novel Nrf2/Keap1 pathway mutations in pediatric acute lymphoblastic leukemia

Acute lymphoblastic leukemia (ALL) is a malignancy of lymphoid progenitor cells, characterized by a wide range of biological and clinical heterogeneity. Oxidative stress is a common problem observed in carcinogenesis and it is involved in developing treatment resistance. Nuclear Factor Erythroid-2-Like 2 (Nrf2) transcription factor is the main regulator of antioxidant responses. The levels of reactive oxygen species (ROS) are tightly controlled and regulated by Nrf2 and its suppressor protein Kelch-like ECH-associated protein 1 (Keap1). Recently, many studies have shown that most of the genes in the Nrf2/Keap1/nuclear factor kappa-B (NF-κB)/phosphotyrosine-independent ligand for the Lck SH2 domain Of 62 KDa (p62) pathway show abnormally high mutational variations in cancer. However, variations in the Nrf2/Keap1/NF-κB1/p62 pathway in pediatric ALL have not been thoroughly investigated, yet. Thirty children, who were diagnosed with pediatirc ALL were included in the study. The Nrf2/Keap1/NF-κB1/p62 pathway variants were analyzed by DNA sequencing analysis. The PolyPhen-2 program was used for identifying pathogenic mutations. Our study examined the molecular dynamics (MD) perspectives of the effect of A159T and E121K mutations on protein stability for the first time in the literature by using the GROMACS45 software package utilizing the OPSLAA force field. Of the detected 17 nucleotide changes, 6 were novel. The study predicted the potential pathological effect of two mutations p. A159T and p.E121K in the Keap1 gene. The MD perspectives revealed that the E121K mutant's observed structural behavior accounted for the key role of His-129 and E121K, where E121K exhibited much higher drift compared to His-129. For a future perspective, it would be meaningful to study the protein-small molecule interactions of the Keap1 protein to elaborate on the drug effects in patients carrying these mutations.