Whole-exome sequencing and immune profiling of early-stage lung adenocarcinoma with fully annotated clinical follow-up

Validation of a Novel Cuproptosis-Related Prognostic Gene Marker and Differential Expression Associated with Lung Adenocarcinoma

BACKGROUND

Cuproptosis induction is seen as a promising alternative for immunotherapies and targeted therapies in breast cancer. The objective of this research was to examine the prognostic and biological importance of cuproptosis-related genes (CRGs) in lung adenocarcinoma (LUAD).

METHODS

The following methods were used: GSE10072 dataset and TCGA database analysis, differential expression analysis of CRGs, and biological function (BP) and signaling pathway enrichment analysis, prognostic analysis and clinical analysis of CRGs, construction of the prognostic signature and RNA modified genes and miRNA analysis of CRGs in LUAD, immunoinfiltration analysis and immunohistochemical staining of DβH, UBE2D3, SOD1, UBE2D1 and LOXL2.

RESULTS

AOC1, ATOX1, CCL8, CCS, COX11, CP, LOXL2, MAP2K2, PDK1, SCO2, SOD1, UBE2D1, UBE2D3 and VEGFA showed significantly higher expression, while ATP7B, DβH, PDE3B, SLC31A2, UBE2D2, UBE2D4 and ULK2 showed lower expression in LUAD tissues than normal tissues. We also found that ATP7B (4%), AOC1 (3%) PDE3B (2%), DβH (2%), CP (1%), ULK2 (1%), PDK1 (1%), LOXL2 (1%) and UBE2D3 (1%) showed higher mutation frequencies. The univariate Cox analysis was used to identify CRGs that have prognostic value. It identified 21 genes that showed significant prognostic value, containing DβH, UBE2D3, SOD1, UBE2D1 and LOXL2. Patients with DβH up-expression have a longer survival time and patients with UBE2D3, SOD1, UBE2D1 and LOXL2 down-expression also have a longer survival time. hsa-miR-29c-3p, hsa-miR-29a-3p, hsa-miR-181c-5p, hsa-miR-1245a, etc., play an important role in the miRNA regulatory network, and in LUAD, miR-29a, miR-29c and miR-181c high expression survival was longer, and miR-1245a low expression survival was longer. We also performed an analysis to examine the relationships between DβH, LOXL2, SOD1, UBE2D1 and UBE2D3 and immune infiltration in LUAD, including B cells, CD8+ T cells, CD4+ T cells, macrophages, neutrophils and DCs.

CONCLUSION

DβH, UBE2D3, SOD1, UBE2D1, and LOXL2 are potential candidates implicated in LUAD and can be further explored for their application as diagnostic, prognostic, and therapeutic biomarkers for LUAD.

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.

Cuproptosis predicts the risk and clinical outcomes of lung adenocarcinoma

Copper is an essential microelement for the body and a necessary coregulator for enzymatic reactions, yet an unbalanced copper level promotes reactive oxidation and cytotoxicity, which ultimately induces cell death. Several small molecules targeting copper-induced cell death have been investigated, yet few showed promising therapeutic effects in clinical trials. In March 2022, Science first introduced the concept and mechanisms of cuproptosis, suggesting that copper-induced cell death targets the tricarboxylic acid (TCA) cycle via protein lipoylation. Does this novel form of cell death take part in tumorigenesis or tumor progression? Is cuproptosis related to clinical outcomes of diseases? Is there a cuproptosis-related panel for clinical practice in cancer treatment? Herein, based on 942 samples of lung adenocarcinoma (LUAD), we analyzed on gene set level the existence and predictive value of cuproptosis in disease diagnosis and treatment. We screened out and identified the “cupLA” panel which indicates the risk of LUAD occurrence, clinicopathological features of LUAD patients, and could guide clinicians to refine LUAD subtypes and make treatment choices.

Clinical Significance and Immune Landscape of Recurrence-Associated Ferroptosis Signature in Early-Stage Lung Adenocarcinoma

BACKGROUND

The prevalence of patients newly diagnosed with early-stage lung adenocarcinoma (LUAD) is growing alongside significant advances in screening approaches. This study aimed to construct ferroptosis-related gene score (FRGscore) for predicting recurrence, explore immune-molecular characteristics, and determine the benefit of immunotherapy in distinct ferroptosis-based patterns and FRGscore-defined subgroups.

METHODS

A total of 1,085 early-stage LUAD patients from four independent cohorts were included. Consensus clustering analysis was performed using 217 co-expressed FRGs to explore different ferroptosis-mediated patterns. An FRG scoring system was established to predict relapse, quantify ferroptosis-mediated patterns, and evaluate the response to immunotherapy in individual patients based on Lasso-penalized and stepwise Cox regression analyses. Immune landscape involving multiple parameters was further evaluated, stratified by cluster subtypes and FRGscore subgroups.

RESULTS

Two ferroptosis-mediated patterns were identified and verified, which were characterized by significantly distinct prognosis and immune profiles. Analyses of immune characteristics showed that identified ferroptosis patterns were characterized as immune-inflamed phenotype and immune-exhausted phenotype. The FRG scoring model based on 11 FRG-derived signatures panel classified patients into the FRGscore-high and FRGscore-low subgroups. Significantly longer recurrence-free survival (RFS) and overall survival (OS) were observed in the FRGscore-low subgroup. FRGscore-low patients were characterized by higher tumor mutational burden (TMB), immunoscore, immunophenoscore, and PD-L1 expression level and were associated with lower Tumor Immune Dysfunction and Exclusion (TIDE) score, whereas the opposite was observed in FRGscore-high patients. Immune-active pathways were remarkably enriched in the FRGscore-low subgroup. This scoring model remained highly predictive of prognosis across different clinical, molecular, and immune subgroups. Further analysis indicated that FRGscore-low patients exhibited higher response to anti-PD-1/PD-L1 immunotherapy and better clinical benefits based on two independent immunotherapy cohorts.

CONCLUSION

The proposed FRGscore could highly distinguish the recurrence patterns and molecular and immune characteristics and could predict immunotherapy prognosis, potentially representing a powerful prognostic tool for further optimization of individuated treatment and management strategies in early-stage LUAD.

Immune Microenvironment Signatures as Biomarkers to Predict Early Recurrence of Stage Ia-b Lung Cancer

Introduction

Approximately 30% of patients diagnosed with stage Ia-b NSCLC die of recurrent disease after surgery. This study aimed to identify immune-related biomarkers that might predict tumor recurrence in stage Ia-b NSCLC within 40 months after curative resection.

Methods

Gene expression data of stage Ia-b NSCLC samples was retrieved from the TCGA database, the GEO databases, and the Second Xiangya hospital (XXEYY) database. 22 types of tumors infiltrating immune cells and the expression of immune-associated genes were investigated using CIBERSORT, immunohistochemical staining, and GSEA analyses in a total of 450 patients (80 in the training cohort and 370 in the validation cohorts). Recurrence-related immune features were selected based on the LASSO Cox regression model.

Results

High density of Tregs, Macrophages M0 and M1 cell could be observed in recurrence group while the memory B cell was more frequently enriched in controls, yet Tregs alone was significantly associated with tumor early recurrence in TCGA cohort, XYEYY cohort and GSE37745 dataset. A handful of immune-related genes were identified in the recurrence group. Based on Lasso regression analysis, the expressions of five immune-related genes, RLTPR, SLFN13, MIR4500HG, HYDIN and TPRG1 were closely correlated with tumor early recurrence. In the training cohort (TCGA), the combination of these five genes has sensitivity and specificity of 85% and 85%, with AUC of 0.91 (95% CI 0.84-0.98) for lung cancer early recurrence prediction, whereas in validation cohorts, the sensitivity and specificity using this panel was 61-89% and 54-82%, with AUC of 0.62-0.84.

Conclusion

Our study demonstrated that the immune microenvironment signatures were closely related to tumor early recurrence. Compared to tumor-infiltrating lymphocytes, the expression of five immune-related genes could be robust biomarkers to predict early recurrence of stage Ia-b NSCLC after curative resection.

The Resistance Mechanisms of Lung Cancer Immunotherapy

Immunotherapy has revolutionized lung cancer treatment in the past decade. By reactivating the host’s immune system, immunotherapy significantly prolongs survival in some advanced lung cancer patients. However, resistance to immunotherapy is frequent, which manifests as a lack of initial response or clinical benefit to therapy (primary resistance) or tumor progression after the initial period of response (acquired resistance). Overcoming immunotherapy resistance is challenging owing to the complex and dynamic interplay among malignant cells and the defense system. This review aims to discuss the mechanisms that drive immunotherapy resistance and the innovative strategies implemented to overcome it in lung cancer.

Insights Into Lung Cancer Immune-Based Biology, Prevention, and Treatment

Lung cancer is the number one cause of cancer-related deaths. The malignancy is characterized by dismal prognosis and poor clinical outcome mostly due to advanced-stage at diagnosis, thereby inflicting a heavy burden on public health worldwide. Recent breakthroughs in immunotherapy have greatly benefited a subset of lung cancer patients, and more importantly, they are undauntedly bringing forth a paradigm shift in the drugs approved for cancer treatment, by introducing "tumor-type agnostic therapies". Yet, and to fulfill immunotherapy's potential of personalized cancer treatment, demarcating the immune and genomic landscape of cancers at their earliest possible stages will be crucial to identify ideal targets for early treatment and to predict how a particular patient will fare with immunotherapy. Recent genomic surveys of premalignant lung cancer have shed light on early alterations in the evolution of lung cancer. More recently, the advent of immunogenomic technologies has provided prodigious opportunities to study the multidimensional landscape of lung tumors as well as their microenvironment at the molecular, genomic, and cellular resolution. In this review, we will summarize the current state of immune-based therapies for cancer, with a focus on lung malignancy, and highlight learning outcomes from clinical and preclinical studies investigating the naïve immune biology of lung cancer. The review also collates immunogenomic-based evidence from seminal reports which collectively warrant future investigations of premalignancy, the tumor-adjacent normal-appearing lung tissue, pulmonary inflammatory conditions such as chronic obstructive pulmonary disease, as well as systemic microbiome imbalance. Such future directions enable novel insights into the evolution of lung cancers and, thus, can provide a low-hanging fruit of targets for early immune-based treatment of this fatal malignancy.

Effects of hub genes on the clinicopathological and prognostic features of lung adenocarcinoma

Lung adenocarcinoma (LUAD) is a common malignancy; however, the majority of its underlying molecular mechanisms remain unknown. In the present study, weighted gene co-expression network analysis was applied to construct gene co-expression networks for the GSE19804 dataset, in order to screen hub genes associated with the pathogenesis of LUAD. In addition, with the aid of the Database for Annotation, Visualization and Integrated Discovery, Gene Ontology, and Kyoto Encyclopedia of Genes and Genomes, pathway enrichment analyses were performed on the genes in the selected module. Using the GSE40791 dataset and The Cancer Genome Atlas database, the hub genes were identified. It was discovered that the turquoise module was the most significant module associated with the tumor stage of LUAD. After performing functional enrichment analyses, it was indicated that the turquoise module was mainly enriched in signal transduction. Additionally, at the transcriptional and translational level, nine hub genes were identified and validated: Carbonic anhydrase 4 (CA4), platelet and endothelial cell adhesion molecule 1 (PECAM1), DnaJ member B4 (DNAJB4), advanced glycosylation end-product specific receptor (AGER), GTPase, IMAP family member 6 (GIMAP6), chromosome 10 open reading frame 54 (C10orf54), dedicator of cytokinesis 4 (DOCK4), Golgi membrane protein 1 (GOLM1) and platelet activating factor acetylhydrolase 1b catalytic subunit 3 (PAFAH1B3). CA4, PECAM1, DNAJB4, AGER, GIMAP6, C10orf54 and DOCK4 were expressed at lower levels in the tumor samples, whereas GOLM1 and PAFAH1B3 were highly expressed in tumor samples. In addition, all hub genes were associated with prognosis. In conclusion, one module and nine genes were recognized to be associated with the tumor stage of LUAD. These findings may enhance the understanding of the progression and prognosis of LUAD.

Genome-Wide and Phenotypic Evaluation of Stem Cell Progenitors Derived From Gprc5a-Deficient Murine Lung Adenocarcinoma With Somatic Kras Mutations

Lung adenocarcinomas (LUADs) with somatic mutations in the KRAS oncogene comprise the most common molecular subtype of lung cancer in smokers and present with overall dismal prognosis and resistance to most therapies. Our group recently demonstrated that tobacco carcinogen-exposed mice with knockout of the airway lineage G-protein coupled receptor, Gprc5a, develop LUADs with somatic mutations in Kras. Earlier work has suggested that cancer stem cells (CSCs) play crucial roles in clonal evolution of tumors and in therapy resistance. To date, our understanding of CSCs in LUADs with somatic Kras mutations remains lagging. Here we derived CSCs (as spheres in 3D cultures) with self-renewal properties from a murine Kras-mutant LUAD cell line we previously established from a tobacco carcinogen-exposed Gprc5a^−/− mouse. Using syngeneic Gprc5a^−/− models, we found that these CSCs, compared to their parental isoforms, exhibited increased tumorigenic potential in vivo, particularly in female animals. Using whole-transcriptome sequencing coupled with pathways analysis and confirmatory PCR, we identified gene features (n = 2,600) differentially expressed in the CSCs compared to parental cells and that were enriched with functional modules associated with an augmented malignant phenotype including stemness, tumor-promoting inflammation and anti-oxidant responses. Further, based on in silico predicted activation of GSK3β in CSCs, we found that tideglusib, an irreversible inhibitor of the kinase, exhibited marked anti-growth effects in the cultured CSCs. Our study underscores molecular cues in the pathogenesis of Kras-mutant LUAD and presents new models to study the evolution, and thus high-potential targets, of this aggressive malignancy.

Genetic polymorphisms of histone methyltransferase SETD2 predicts prognosis and chemotherapy response in Chinese acute myeloid leukemia patients

BACKGROUND

SETD2, the single mediator of trimethylation of histone 3 at position lysine 36, has been reported associated with initiation progression and chemotherapy resistance in acute myeloid leukemia (AML). Whether polymorphisms of SETD2 affect prognosis and chemotherapy response of AML remains elusive.

METHODS

Three tag single-nucleotide polymorphisms (tagSNPs) of SETD2 were genotyped in 579 AML patients by using Sequenom Massarray system. Association of the SNPs with complete remission (CR) rate after Ara-C based induction therapy, overall survival (OS) and relapse-free survival (RFS) were analyzed.

RESULT

Survival analysis indicated that SETD2 rs76208147 TT genotype was significantly associated with poor prognosis of AML (TT vs. CC + CT hazard ratio: HR = 1.838, 95% confidence interval (CI) 1.005-3.360, p = 0.048). After adjusting for the known prognostic factors including risk stratification, age, allo-SCT, WBC count and LDH count, rs76208147 TT genotype was still associated with OS in the multivariate analysis (TT vs. CC + CT HR = 1.923, 95% CI 1.007-3.675, p = 0.048). In addition, after adjusting by other clinical features, patients with rs4082155  allele G carries showed higher rate of complete remission which indicated by CR rate (AG + GG vs. AA odd ratio (OR) = 0.544, 95% CI 0.338-0.876, p = 0.012).

CONCLUSIONS

SETD2 genetic polymorphism is associated with AML prognosis and chemotherapy outcome, suggesting the possibility for development in AML diagnostics and therapeutics towards SETD2.

Targeting the cell signaling pathway Keap1-Nrf2 as a therapeutic strategy for adenocarcinomas of the lung

INTRODUCTION

Kelch-like ECH associated protein 1/Nuclear factor erythroid 2-like factor 2 (Keap1-Nrf2) signaling plays a pivotal role in response to oxidative stress in lung cancer. Mutations in KEAP1/NFE2L2 genes always cause persistent Nrf2 activation in lung cancer cells that confer therapeutic resistance and aggressive tumorigenic activity, dictating either poor prognosis or short duration of response to chemotherapy in clinical observations. Areas covered: We provide a review of the mechanisms underlying the regulation of Keap1-Nrf2 at different stages, including genetic mutations, epigenetic modifications, translational/post-translational alterations, and protein-protein interactions. Based on the current knowledge, we discuss the possibilities of intervening Keap1-Nrf2 in lung adenocarcinoma as a therapeutic target. Expert opinion: It is prevalently conceived that Keap1-Nrf2 signaling plays different roles at diverse stages of cancer. Although various Nrf2 or Keap1 inhibitors have been reported during the last decades, none of these inhibitors are currently under clinical studies or in clinical applications, suggesting that sole inhibition of Nrf2 might not be sufficient to suppress tumor growth. On the basis of current studies, we suggest that the rational combination of Nrf2 suppression with chemical agents which cause enhanced oxidative imbalance or abnormal metabolism would be promising in the treatment of lung adenocarcinoma.

Development of Kras mutant lung adenocarcinoma in mice with knockout of the airway lineage-specific gene Gprc5a

Despite the urgency for prevention and treatment of lung adenocarcinoma (LUAD), we still do not know drivers in pathogenesis of the disease. Earlier work revealed that mice with knockout of the G-protein coupled receptor Gprc5a develop late onset lung tumors including LUADs. Here, we sought to further probe the impact of Gprc5a expression on LUAD pathogenesis. We first surveyed GPRC5A expression in human tissues and found that GPRC5A was markedly elevated in human normal lung relative to other normal tissues and was consistently downregulated in LUADs. In sharp contrast to wild-type littermates, Gprc5a-/- mice treated chronically with the nicotine-specific carcinogen NNK developed LUADs by 6 months following NNK exposure. Immunofluorescence analysis revealed that the LUADs exhibited abundant expression of surfactant protein C and lacked the clara cell marker Ccsp, suggesting that these LUADs originated from alveolar type II cells. Next, we sought to survey genome-wide alterations in the pathogenesis of Gprc5a-/- LUADs. Using whole exome sequencing, we found that carcinogen-induced LUADs exhibited markedly higher somatic mutation burdens relative to spontaneous tumors. All LUADs were found to harbor somatic mutations in the Kras oncogene (p. G12D or p. Q61R). In contrast to spontaneous lesions, carcinogen-induced Gprc5a-/- LUADs exhibited mutations (variants and copy number gains) in additional drivers (Atm, Kmt2d, Nf1, Trp53, Met, Ezh2). Our study underscores genomic alterations that represent early events in the development of Kras mutant LUAD following Gprc5a loss and tobacco carcinogen exposure and that may constitute targets for prevention and early treatment of this disease.