Necroptosis, a type of necrotic cell death independent of caspase regulation, is mainly mediated by receptor interacting serine/threonine kinase 1 (RIPK1), receptor interacting serine/threonine kinase 3 (RIPK3) and mixed lineage kinase domain-like (MLKL). Necroptosis plays an essential role in many tumors. However, the potential roles of necroptosis in tumor microenvironment (TME) of sarcoma (SARC) remain unknown. This study analyzed the expression, prognosis, genetic alterations of necroptosis genes in SARC. We identified two subtypes (cluster A and B) by performing unsupervised consensus clustering. Cluster A and B greatly differed in prognosis and immune infiltration, with cluster A showing more favorable prognosis, higher immune infiltration and higher expression levels of necroptosis genes than cluster B. Based on the differentially expressed genes (DEGs) between two clusters, a necroptosis scoring system was developed for predicting overall survival of SARC patients. Patients with high necroptosis score had worse survival status, with a decreased infiltration level of most immune cells. Our findings demonstrated the potential role of necroptosis in regulating tumor microenvironment and the prognostic value of necroptosis-related genes for SARC patients.

Autophagy is a catabolic pathway involved in the regulation of bone homeostasis. We explore clinical correlation of autophagy-related key molecules to establish risk signature for predicting the prognosis, tumor microenvironment (TME), and immunotherapy response of osteosarcoma. Single cell RNA sequencing data from GSE162454 dataset distinguished malignant cells from normal cells in osteosarcoma. Autophagy-related genes (ARGs) were extracted from the established risk signature of the Molecular Signatures Database of Gene Set Enrichment Analysis (GSEA) by univariate Cox and least absolute shrinkage and selection operator (LASSO) Cox regression analysis. Overall survival (OS), TME score, abundance of infiltrating immune cells, and response to immune-checkpoint blockade (ICB) treatment in patients with different risks were compared based on risk score. Nine ARGs were identified and risk signature was constructed. In all osteosarcoma datasets, the OS was significantly longer in the high-risk patients than low-risk onset. Risk signature significantly stratified clinical outcomes, including OS, metastatic status, and survival status. Risk signature was an independent variable for predicting osteosarcoma OS and showed high accuracy. A nomogram based on risk signature and metastases was developed. The calibration curve confirmed the consistency in 1-year, 3-year, and 5-year predicted OS and the actual OS. The risk score was related to 6 kinds of T cells and macrophages, myeloid-derived suppressor cell, natural killer cell, immune score, and stromal score in TME. The risk signature helped in predicting patients’ response to anti-PD1/anti-PD-L1 treatment. The ARGs-led risk signature has important value for survival prediction, risk stratification, tumor microenvironment, and immune response evaluation of osteosarcoma.

Aging is an inevitable process that biological changes accumulate with time and results in increased susceptibility to different tumors. But currently, aging-related genes (ARGs) in osteosarcoma were not clear. We investigated the potential prognostic role of ARGs and established an ARG-based prognostic signature for osteosarcoma. The transcriptome data and corresponding clinicopathological information of patients with osteosarcoma were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. Molecular subtypes were generated based on prognosis-related ARGs obtained from univariate Cox analysis. With ARGs, a risk signature was built by univariate, least absolute shrinkage and selection operator (LASSO), and multivariate Cox regression analyses. Differences in clinicopathological features, immune infiltration, immune checkpoints, responsiveness to immunotherapy and chemotherapy, and biological pathways were assessed according to molecular subtypes and the risk signature. Based on risk signature and clinicopathological variables, a nomogram was established and validated. Three molecular subtypes with distinct clinical outcomes were classified based on 36 prognostic ARGs for osteosarcoma. A nine-ARG-based signature in the TCGA cohort, including BMP8A, CORT, SLC17A9, VEGFA, GAL, SSX1, RASGRP2, SDC3, and EVI2B, has been created and developed and could well perform patient stratification into the high- and low-risk groups. There were significant differences in clinicopathological features, immune checkpoints and infiltration, responsiveness to immunotherapy and chemotherapy, cancer stem cell, and biological pathways among the molecular subtypes. The risk signature and metastatic status were identified as independent prognostic factors for osteosarcoma. A nomogram combining ARG-based risk signature and metastatic status was established, showing great prediction accuracy and clinical benefit for osteosarcoma OS. We characterized three ARG-based molecular subtypes with distinct characteristics and built an ARG-based risk signature for osteosarcoma prognosis, which could facilitate prognosis prediction and making personalized treatment in osteosarcoma.