Identification of molecular subtypes, risk signature, and immune landscape mediated by necroptosis-related genes in non-small cell lung cancer.

ABSTRACT

Background: Non-small cell lung cancer (NSCLC) is a highly heterogeneous malignancy with an extremely high mortality rate. Necroptosis is a programmed cell death mode mediated by three major mediators, RIPK1, RIPK3, and MLKL, and has been shown to play a role in various cancers. To date, the effect of necroptosis on NSCLC remains unclear. Methods: In The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases, we downloaded transcriptomes of lung adenocarcinoma (LUAD) patients and their corresponding clinicopathological parameters. We performed multi-omics analysis using consensus clustering based on the expression levels of 40 necroptosis-related genes. We constructed prognostic risk models and used the receiver operating characteristic (ROC) curves, nomograms, and survival analysis to evaluate prognostic models. Results: With the use of consensus clustering analysis, two distinct subtypes of necroptosis were identified based on different mRNA expression levels, and cluster B was found to have a better survival advantage. Correlation results showed that necroptosis was significantly linked with clinical features, overall survival (OS) rate, and immune infiltration. Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment analysis confirmed that these differential genes were valuable in various cellular and biological functions and were significantly enriched in various pathways such as the P53 signaling pathway and cell cycle. We further identified three genomic subtypes and found that gene cluster B patients had better prognostic value. Multivariate Cox analysis identified the 14 best prognostic genes for constructing prognostic risk models. The high-risk group was found to have a poor prognosis. The construction of nomograms and ROC curves showed stable validity in prognostic prediction. There were also significant differences in tumor immune microenvironment, tumor mutational burden (TMB), and drug sensitivity between the two risk groups. The results demonstrate that the 14 genes constructed in this prognostic risk model were used as tumor prognostic biomarkers to guide immunotherapy and chemotherapy. Finally, we used qRT-PCR to validate the genes involved in the signature. Conclusion: This study promotes our new understanding of necroptosis in the tumor microenvironment of NSCLC, mines prognostic biomarkers, and provides a potential value for guiding immunotherapy and chemotherapy.