Prediction of Prognosis and Immunotherapy Response with a Novel Natural Killer Cell Marker Genes Signature in Osteosarcoma

Background: Natural killer (NK) cells are characterized by their antitumor efficacy without previous sensitization, which have attracted attention in tumor immunotherapy. The heterogeneity of osteosarcoma (OS) has hindered therapeutic application of NK cell-based immunotherapy. The authors aimed to construct a novel NK cell-based signature to identify certain OS patients more responsive to immunotherapy. Materials and Methods: A total of eight publicly available datasets derived from patients with OS were enrolled in this study. Single-cell RNA sequencing data obtained from the Gene Expression Omnibus (GEO) database were analyzed to screen NK cell marker genes. Least Absolute Shrinkage and Selection Operator (LASSO) Cox regression analysis was used to construct an NK cell-based prognostic signature in the TARGET-OS dataset. The differences in immune cell infiltration, immune system-related metagenes, and immunotherapy response were evaluated among risk subgroups. Furthermore, this prognostic signature was experimentally validated by reverse transcription-quantitative real-time PCR (RT-qPCR). Results: With differentially expressed NK cell marker genes screened out, a five-gene NK cell-based prognostic signature was constructed. The prognostic predictive accuracy of the signature was validated through internal clinical subgroups and external GEO datasets. Low-risk OS patients contained higher abundances of infiltrated immune cells, especially CD8 T cells and naive CD4 T cells, indicating that T cell exhaustion states were present in the high-risk OS patients. As indicated from correlation analysis, immune system-related metagenes displayed a negative correlation with risk scores, suggesting the existence of immunosuppressive microenvironment in OS. In addition, based on responses to immune checkpoint inhibitor therapy in two immunotherapy datasets, the signature helped predict the response of OS patients to anti-programmed cell death protein 1 (PD-1) or anti-programmed cell death ligand 1 (PD-L1) therapy. RT-qPCR results demonstrated the roughly consistent relationship of these five gene expressions with predicting outcomes. Conclusions: The NK cell-based signature is likely to be available for the survival prediction and the evaluation of immunotherapy response of OS patients, which may shed light on subsequent immunotherapy choices for OS patients. In addition, the authors revealed a potential link between immunosuppressive microenvironment and OS.