Yucheng Fu, Guoyu He, Zhuochao Liu, Jun Wang, Zhusheng Zhang, Qiyuan Bao, Junxiang Wen, Zhijian Jin, Weibin Zhang Department of Orthopedics, Shanghai Key Laboratory for Prevention and Treatment of Bone and Joint Diseases, Shanghai Institute of Traumatology and Orthopedics, Ruijin Hospital, Shanghai Jiaotong University School of Medicine, Shanghai, 200025, People’s Republic of ChinaCorrespondence: Weibin ZhangDepartment of Orthopedics, Shanghai Key Laboratory for Prevention and Treatment of Bone and Joint Diseases, Shanghai Institute of Traumatology and Orthopedics, Ruijin Hospital, Shanghai Jiaotong University School of Medicine, 197 Ruijin Er Road, Shanghai, 200025, People’s Republic of ChinaTel +86-21-64370045, Ext 666983Fax +86-21-54660217Email zhangweibin10368@163.comBackground: Inflammatory response took part in the progression of tumor and was regarded as the hallmark of cancer. However, the prognostic relationship between osteosarcoma and inflammatory response-associated genes (IRGs) was unclear. This research aimed to explore the correlations between osteosarcoma prognosis and IRG signature.Methods: The inflammatory response-associated differentially expressed messenger RNAs (DEmRNAs) were screened out through Gene Expression Omnibus (GEO) and Molecular Signature Database (MSigDB) databases. Univariate and multivariate cox regression analyses were utilized to construct the IRG signature. The prognostic value of signature was investigated through Kaplan–Meier (KM) survival curve and nomogram. DEmRNAs among high and low inflammatory response-associated risks were identified and functional enrichment analyses were conducted. ESTIMATE, CIBERSORT and single-sample gene set enrichment analyses (ssGSEA) were implied to reveal the alterations in immune infiltration. All the above results were validated in Target database. The expression of IRGs was also validated in different cell lines by quantitative real-time PCR (qRT-PCR) and osteosarcoma patient samples by immunohistochemistry.Results: The IRG signature that consisted of two genes (MYC, CLEC5A) was established. In training and validation datasets, patients with lower risk scores survived longer and the IRG signature was confirmed as the independent prognostic factor in osteosarcoma. The nomogram was constructed and the calibration curves demonstrated the reliability of this model. Functional analysis of risk score-associated DEmRNAs indicated that immune-related pathways and functions were significantly enriched. ssGSEA revealed that 14 immune cells and 11 immune functions were significantly dysregulated. The qRT-PCR results indicated IRGs were significantly differently expressed in osteosarcoma and osteoblast cell lines. The immunohistochemistry analyses of patients’ samples revealed the same result.Conclusion: The novel osteosarcoma inflammatory response-associated prognostic signature was established and validated in this study. This model could serve as the biomarker and therapeutic target for osteosarcoma in the future.Keywords: osteosarcoma, inflammatory response, metastasis, prognosis, immune
