Immunogenic cell death related mRNAs associated signature to predict immunotherapeutic response in osteosarcoma

Immunogenic cell death (ICD) is related to cancer prognosis, which has a synergic effect in combination with chemotherapy or immunotherapy. Yet, the relationship between ICD and osteosarcoma remained unclear.

Three osteosarcoma datasets including therapeutically applicable research to generate effective treatments (TARGET), GSE126209 and GSE21257 datasets were included. A protein-protein interaction network was constructed based on ICD-related genes. We performed unsupervised consensus clustering to classify molecular subtypes (clusters). Survival analysis, Estimation of stromal and immune cells in malignant tumour tissues using expression data (ESTIMATE), Cell-type identification by estimating relative subsets of RNA transcripts (CIBERSORT), and differential analysis were employed to characterize the molecular differences between different clusters. Univariate Cox regression analysis was conducted to confirm prognostic genes. Quantitative reverse transcription polymerase chain reaction (qRT-PCR) was used to demonstrate the aberrant expression of ICD-correlated signature genes in osteosarcoma. A series of cellular experiments, including cell counting kit-8 (CCK-8), transwell, and flow cytometry, were used to demonstrate the regulatory role of key genes in the ICD model on the malignant phenotype of osteosarcoma.

Three clusters (cluster1, 2, 3) were constructed and they showed distinct overall survival and immune infiltration. ICD-related genes were highly expressed in cluster1. Moreover, Cluster1 had the best prognosis, high immune score and high expression of human leukocyte antigen (HLA)-related genes. TLR4, LY96, IFNGR1, CD4, and CASP1 were identified as prognostic genes for establishing an ICD-related risk signature. According to the risk signature, two risk groups (high and low risks) showing differential prognosis and response to immunotherapy. The low risks group had a better prognosis but was not sensitive to immunotherapy. Molecular assays verified that prognostic genes were abnormally under-expressed in osteosarcoma. Cellular assays demonstrated that LY96, the most significantly down-regulated gene in osteosarcoma, inhibited the migration, invasion, and proliferation phenotypes of osteosarcoma cells and prolonged the cell cycle. Analysis of oxidative stress related pathway enrichment in tumor microenvironment was conducted by single-sample gene set enrichment analysis (ssGSEA).

This study demonstrated the prognostic significance of ICD-correlated genes in osteosarcoma patients. The five-gene risk signature facilitate prognostic evaluation and prediction of osteosarcoma patients' response to immunotherapy. The risk signature also offered a possibility for the exploit of novel ICD-related treatment.