Abstract
Natural killer (NK) cells, the effectors of the innate immune system, have a remarkable influence on cancer prognosis and immunotherapy. In this study, a total of 1,816 samples from nine independent cohorts in public datasets were enrolled. We first conducted a comprehensive analysis of single-cell RNA-sequencing data of lung adenocarcinoma (LUAD) from the Gene Expression Omnibus (GEO) database and determined 189 NK cell marker genes. Subsequently, we developed a seven-gene prognostic signature based on NK cell marker genes in the TCGA LUAD cohort, which stratified patients into high-risk and low-risk groups. The predictive power of the signature was well verified in different clinical subgroups and GEO cohorts. With a multivariate analysis, the signature was identified as an independent prognostic factor. Low-risk patients had higher immune cell infiltration states, especially CD8(+) T cells and follicular helper T cells. There existed a negative association between inflammatory activities and risk score, and the richness and diversity of the T-cell receptor (TCR) repertoire was higher in the low-risk groups. Importantly, analysis of an independent immunotherapy cohort (IMvigor210) revealed that low-risk patients had better immunotherapy responses and prognosis than high-risk patients. Collectively, our study developed a novel signature based on NK cell marker genes, which had a potent capability to predict the prognosis and immunotherapy response of LUAD patients.