
==== Front
Heliyon
Heliyon
Heliyon
2405-8440
Elsevier

S2405-8440(24)11750-1
10.1016/j.heliyon.2024.e35719
e35719
Research Article
A diagnostic signature developed based on the necroptosis-related genes and its association with immune infiltration in osteosarcoma
Bian Yiying
Shi Jixiang
Chen Ziyun
Fang Ji
Chen Weidong
Zou Yutong
Yao Hao
Tu Jian
Liao Yan
Xie Xianbiao xiexbiao@126.com
⁎
Shen Jingnan shenjn@mail.sysu.edu.cn
⁎⁎
Department of Musculoskeletal Oncology, The First Affiliated Hospital of Sun Yat-sen University, Guangzhou, 510080, China
⁎ Corresponding author. xiexbiao@126.com
⁎⁎ Corresponding author. shenjn@mail.sysu.edu.cn
10 8 2024
30 8 2024
10 8 2024
10 16 e3571927 5 2024
1 8 2024
2 8 2024
© 2024 The Authors
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).
Introduction

Osteosarcoma is a bone-derived malignancy that often leads to lung metastasis and death.

Material and methods

The RNA-seq data of TARGET-osteosarcoma were collected from TARGET database. GSE16088 and GSE12865 datasets of osteosarcoma x from Gene Expression Database (GEO) were donwloaded. ConsensusClusterPlus was used for molecular subtype classification. Univariate Cox and Lasso regression was employed to develop a risk model. To analyze the regulatory effects of model feature genes on the malignant phenotype of osteosarcoma cell lines, qRT-PCR, Transwell and wound healing assays were performed. The abundance of immune cell infiltration was assessed using MCP-Counter, Gene Set Enrichment Analysis (GSEA), and ESTIMATE. The Tumor Immune Dysfunction and Exclusion (TIDE) software was employed to evaluate immunotherapy and response to conventional chemotherapy drugs.

Results

Three clusters (C1, C2 and C3) were classified using 39 necroptosis score-associated genes. In general, C1 and C2 showed better prognosis outcome and lower death rate than C3. Specifically, C2 could benefit more from immunotherapy, while C3 was more sensitive to traditional medicines, and C1 had higher immune cell infiltration. Next, an 8-gene signature and a risk score model were developed, with a low risk score indicating better survival and immune cell infiltration. ROC analysis showed that 1-, 3-, and 5-year overall survival of osteosarcoma could be correctly predicted by the risk score model. Cellular experiments revealed that the model feature gene IFITM3 promoted the osteosarcoma cell migration and invasion. Furthermore, the overall survival of osteosarcoma patients from TARGET and validation datasets can be accurately evaluated using the nomogram model.

Conclusions

Our prognostic model developed using necroptosis genes could facilitate the prognostic prediction for patients suffering from osteosarcoma, offering potential osteosarcoma targets.

Keywords

Necroptosis
Osteosarcoma
Prognosis
Risk model
Tumor microenvironment (TME)
==== Body
pmc1 Introduction

Osteosarcoma is an aggressive bone tumor that frequently occurs to children and adolescents [1]. Lung metastasis is detected in around 15 %–20 % of patients with osteosarcoma at their first diagnosis [2,3]. According to histological findings, osteosarcoma can be divided into high grade, medium grade and low grade [4]. Neoadjuvant chemotherapy has provided more treatment choices than before, but the overall survival of patients with osteosarcoma has not been significantly improved [5]. Hence, finding new prognostic factors and therapeutic targets is necessary for managing osteosarcoma.

Necroptosis [6] is a non-caspase-dependent cell death characterized by both necrosis and apoptosis and is regulated by receptor-interacting protein (RIP). RIP1 fulfills a crucial function in the signal transduction of necrotizing apoptosis. Under abnormal conditions of the body, significantly increased production of caspase-8 will initiate the apoptotic signaling pathway [7], which can completely inhibit RIP1-mediated proteolytic activation of procaspase-8, block apoptosis signaling pathway, and then transform apoptosis into necrotizing apoptosis. Different stimuli are associated with necrotizing apoptosis, and TNF-α combined with Z-VAD-FMK (T/Z) is the classic model of inducing necroptosis [8,9]. In 2016, study reported that pancreatic ductal adenocarcinoma has high-expressed RIPK1 and RIPK3, which are the key components of necroptosis, but the suppression of RIPK3-mediated necrotizing apoptosis could delay the cancer progression [10]. Necroptosis is also present in the necrotic areas in mice with breast cancer, and the inflammatory response caused by necroptosis can cause lung metastasis [11]. Apart from these findings, numerous studies have been conducted on necroptosis in osteosarcoma [[12], [13], [14]].

In recent years, tumor diagnostic and prognostic biomarkers are increasingly discovered using public databases. Microarray technology is becoming increasingly useful in gene analysis. Gene expression profile database is a useful tool with significant clinical value in medical oncology [[15], [16], [17]]. Developing risk models for the adjuvant detection and management of multiple cancer types can improve the precision and effectiveness of treatment protocols [18,19]. Thus, in this study, TARGET-osteosarcoma, related osteosarcoma datasets of GSE12865 and GSE16088 were retrieved from the TARGET database and GEO database, respectively, and the potential key genes were identified to develop a prognostic risk model for the cancer. The patterns of immune infiltration were analyzed to provide a novel possibility for the pathogenesis and treatment strategy of osteosarcoma.

2 Materials and methods

2.1 Raw data

We downloaded the RNA sequencing (RNA-seq) data of “TARGET-Osteosarcoma” from the TARGET (https://ocg.cancer.gov/programs/target) database. A sum of 79 osteosarcoma samples with 19,533 genes were retained by removing samples without survival time, survival status or clinical follow-up information while retaining mRNAs encoding proteins (Supplementary Table 1). Additionally, the dataset GSE21257 (45 osteosarcoma samples) and GSE39058 (36 osteosarcoma samples) were downloaded from GEO database by GPL10295 and GPL14951 [20].The data from the TARGET-osteosarcoma dataset were used as a training cohort, while those from the GSE21257 and GSE39058 datasets were used as a validation cohort. In total, 74 necroptosis-correlated genes were collected following a previous study [21].

2.2 Molecular subtypes

Firstly, necroptosis score was calculated for TARGET-osteosarcoma using ssGSEA. Genes (cor >0.5 and P < 0.05) related to the necroptosis score were selected by Hmisc package. Univariable Cox analysis was performed to screen the prognostic genes from TARGET-osteosarcoma dataset using the Coxph function of R package survival, and the liminal value was defined when p < 0.05. Using these genes, the samples in the TARGET-osteosarcoma were subjected to molecular subtyping using the R package Consensus Cluster Plus 1.52.0 [22]. Applying Pam arithmetic and “spearman” distance, 500 bootstraps were conducted with each bootstrap containing at least 80 % of the samples from the TARGET-osteosarcoma dataset. The Cluster number k was between the range of 2 and 10, and according to cumulative distribution function (CDF) and AUC, the optimum k was determined.

2.3 Analysis of immune infiltration

R software ESTIMATE [23] was employed to calculate overall immunocyte infiltration (Immune Score), stroma level (Stromal Score), and the combination (ESTIMATE Score) for the samples in the TARGET-osteosarcoma cohort using Wilcox.test analysis.

2.4 SsGSEA

The GSEA was used to estimate a total of 28 subpopulations of TILs including the major innate immunity-correlated cell types including natural killer T (NKT) cells, macrophages, monocytes, etc.

2.5 TIDE

TIDE [24,25] algorithm (http://tide.dfci.harvard.edu) was run to evaluate the myeloid suppressor cells (MDSC), dysfunction of tumor infiltration cytotoxic T lymphocytes (CTL) (Dysfunction), IFNG, Exclusion of CTL by immunosuppressive factors (Exclusion), and M2 subtypes of M2 (TAM), which are the three types of cells that restrict T cell invasion into tumors.

2.6 Drug sensitivity analysis

Sensitivities of Erlotinib, MG-132, Z-LLNle-CHO, Dasatinib, CGP-60474 and WH-4-023 to IC50 were predicted using pRRophetic [26].

2.7 Gene set enrichment analysis (GSEA)

GSEA was conducted on the basis of the H.all.v7.5.1.Entrez.GMT gene set in the MSigDB database [27], with FDR<0.05 indicating significant enrichment.

2.8 Development of a prognosis model for osteosarcoma

To filter DEGs and prognostic genes from the molecular subtypes classified, the limma package and univariable Cox analysis were employed. Lasso regression was conducted using the glmnet package [28] to reduce the range of genes. The following formula was calculated to assess the prognosis of osteosarcoma:RiskScore=∑k=0nβi×Expi

where βi and Expi refer to the Cox regression coefficient and expression of the i gene, respectively. The median RiskScore value was used to divide low- and high-risk groups of the samples in the training and independent validation datasets. The prognostic prediction was further assessed using KM survival curve and ROC.

2.9 MCP-counter and CIBERSORT algorithm

The 8 immune populations and 2 stromal populations were quantified using MCP-counter for all the samples [29]. The CIRERSORT algorithm was run to assess the landscape of immune infiltration of 22 immune cells in different risk groups [30].

2.10 Nomogram development

The model was further evaluated by developing a nomogram with other clinicopathological features of patients with osteosarcoma (e.g. age, gender, metastasis) using RMS software R package. The decision curve was plotted to test the nomogram prediction.

2.11 Cell culture and transfection

The osteosarcoma cell lines hFOB1.19 (BNCC338626) and Saos-2 (BNCC338485) were purchased from BNCC (Beijing). Dulbecco's Modified Eagle Medium (Gibco, 11,965–092) added with 10 % fetal bovine serum (Gibco, 26,140–095) and 1 % antibiotics (Gibco, 15,070–063) were used for cell culture at 37 °C in 5 % CO2. Next, cell transfection with negative control (NC) and IFITM3 siRNA (Sagon Co., China) was conducted with the use of Lipofectamine 2000 (Invitrogen, USA).

2.12 QRT-PCR

TRIzol reagent was used for total RNA extraction (Invitrogen, Carlsbad, California, USA) from cells. RNA was reverse-transcribed into cDNA using the Qiagen One-Step RT-PCR kit (Qiagen Gmbh, Germany) and subjected to qRT-PCR experiments. Amplification was performed using SYBR Green (Bimake, Houston, Texas, USA) on the ABI 7500 system (Thermo Fisher Scientific, USA). The primers for qRT-PCR are shown in Supplementary Table 2.

2.13 Wound-healing assay

Osteosarcoma Saos-2 cellss were seeded into 6-well plate to grow until they covered the entire bottom. A scratch was made vertically using a 200 μL pipette tip. After washing the cells twice with phosphate-buffered saline, the cells were photographed under an inverted microscope at 0 h and 48 h after the scratching. The wound healing rate was calculated as [(Width at 0 h - Width at 48 h)/Width at 0 h] × 100 %. The experiment was performed in triplicate.

2.14 Transwell analysis

Saos-2 cells were seeded into the upper chamber of the Transwell apparatus (Corning, USA) containing serum-free medium. Subsequently, the lower chamber was supplemented with PRMI-1640 medium with 10 % FBS. After 1-day incubation at 37 °C, cells remaining on the upper chamber were removed, whereas the rest cells were fixed by 4 % paraformaldehyde and subsequently dyed by crystal violet for 30 min. A microscope was employed to observe the cells and count the cell number.

3 Results

3.1 Three molecular subtypes were defined using the necroptosis genes

There were a total of 116 genes associated with necroptosis score in TARGET-osteosarcoma dataset. Next, 39 prognostic genes (34 protective genes and 5 risk genes) for osteosarcoma were screened by univariate cox survival analysis. Using the 39 genes, ConsensusClusterPlus was employed to group patients in TARGET-osteosarcoma dataset. As shown in the relative change in area under CDF curve (Fig. 1A and B), C1, C2 and C3 were clustered when K = 3 (Fig. 1C). Specifically, the prognosis were the worst in C3 in TARGET-osteosarcoma (Fig. 1D) and GSE39058 dataset (Fig. 1E), whereas the samples in C1 and C2 had better survival. Risk genes had high expressions in C3, while protective genes were high-expressed in C1 and C2 (Fig. 1F).Fig. 1 Classification of 3 molecular subtypes. A: Cumulative distribution function. B: Delta area of Cumulative distribution function. C: when k = 3, 3 molecular subtypes were identified. D: KM survival curve of 3 molecular subtypes in TARGET-osteosarcoma. E: KM survival curve of 3 molecular subtypes in GSE39058 dataset. F: Heatmap of necroptosis genes expressions in 3 molecular subtypes.

Fig. 1

The distribution of clinical characteristics in C1, C2 and C3 was shown in Fig. 2, which demonstrated that patients’ status in TARGET-osteosarcoma cohort was noticeably different.Fig. 2 Clinical features distribution in 3 molecular subtypes.

Fig. 2

3.2 Immune microenvironment and immunotherapy analysis for the three molecular subtypes

Immune cell infiltration in the TARGET-osteosarcoma cohort was assessed based on the expressions of immune cells. Compared to the other two subtypes, C1 had the highest StromalScore, ImmuneScore and ESTIMATEScore (Fig. 3A) and the highest infiltration of 22 out of 28 immune cells (Fig. 3B). We also observed that most of the immune checkpoint genes were high-expressed in C1 than in C2 and C3 (Fig. 3C). The results of TIDE showed higher Exclusion, TIDE, and MDSC in C2 group and higher IFNG and Dysfunctions in C1. In addition, C3 had more tumor-associated macrophage M2 subtypes than in C1 and C2 (Fig. 3D). Moreover, the IC50 of Erlotinib, MG-132, Z-LLNle-CHO, Dasatinib, CGP-60474 and WH-4-023 was higher in C3, suggesting a lower sensitivity of C3 patients to these drugs (Fig. 3E).Fig. 3 Immune microenvironment and immunotherapy analysis. A: Differences of the 3 molecular subtypes in ESTIMATEScore, StromalScore, and ImmuneScore. B: Differences of 28 kinds of immunes score in the 3 clusters. C: Differences of the expression of immune checkpoint genes in the 3 clusters. D: TIDE analysis among C1, C2, C3. E: The estimated IC50 for drug in TARGET-osteosarcoma was shown in the box plots. (*P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001).

Fig. 3

3.3 Pathways analysis for the three molecular subtypes

The GSEA showed that compared with genes C3, genes in C1 in the TARGET cohort were significantly enriched to 21 pathways, including epithelial mesenchymal transition (EMT), IFN-α response, interferon-γ response, etc. (Fig. 4A). Subsequently, we also compared the differential pathways between the three subtypes (including C1vsC2, C1vsC3, C2vsC3). Among the differential pathways between C2 and C3, the enrichment scores in EMT and UV response DN pathways were the highest (Fig. 4B). As shown in Fig. 4C, the immune-related pathways (including interferon response and IL6-JAK-STAT3 signaling pathway) in patients in C1 subtype were activated but cell cycle-related pathways were inhibited.Fig. 4 Pathways enrichment analysis. A: GSEA analysis of C1vsC3 in TARGET-osteosarcoma. B: GSEA analysis of C1vsC3, C1vs C2, C2vs C3 in TARGET-osteosarcoma. C: Radar plot of activated pathways in C1vsC2, C1vsC3, and C2vsC3 in the TARGET-osteosarcoma cohort.

Fig. 4

3.4 Establishment of a prognosis risk model

A total of 114 prognosis genes (48 risk genes and 66 protective genes) were selected out of 983 DEGs from C1, C2, and C3 (Fig. 5A). The trajectory and confidence interval of lambda in Lasso analysis was presented (Fig. 5B and C). When lambda = 0.1547, a risk model was developed using 8 genes: RiskScore=(−0.516*IFITM3)+0.068*ACTA2+(−0.108*GBP1)+(−0.355*APBB1IP)+(−0.26*GJA5)+0.116*CGREF1+(−0.309*CDK6)+0.26*TAC4.Fig. 5 Identification of hub genes. A: 114 genes associated with prognosis were identified. B: As lambda changes, the trajectory of the independent variable was shown. C: Confidence interval under lambda. D: The distribution of Lasso coefficients of genes related to in the necroptosis-based signature.

Fig. 5

3.5 Validation of the prognosis risk model

The risk score for all the patients in TARGET-osteosarcoma dataset was calculated with the formula (Fig. 6A). Subsequently, osteosarcoma patients were categorized into two risk groups (high and low risks) by the risk score based on the threshold of “0". In the TARGET-osteosarcoma cohort, the AUC for 1-year, 3-year and 5-year survival was 0.8, 0.9 and 0.88, respectively (Fig. 6B). The KM survival curve showed a better OS in low-risk group in TARGET-osteosarcoma dataset (Fig. 6C). In GSE39058 cohort, the AUC for 1-, 3- and 5-year survival was 0.8, 0.86 and 0.9, respectively (Fig. 6D), with the low-risk samples having better survival (Fig. 6E). In GSE21257 cohort, the AUC for 1-, 3- and 5-year survival was 0.73, 0.7 and 0.73, respectively (Fig. 6F), with the low-risk samples having longer survival (Fig. 6G).Fig. 6 Validation of the risk score. A: The risk score, survival status, and gene expression in TARGET-osteosarcoma dataset. B: ROC analysis of the risk score in TARGET-osteosarcoma. C: KM survival curve between low group and high group in TARGET-osteosarcoma. D: ROC analysis of the risk score in GSE39058 dataset. E: KM survival curve between the two groups in GSE39058 dataset. F: ROC analysis of the risk score in GSE21257 dataset. G: KM survival curve between the two groups in GSE21257 dataset.

Fig. 6

3.6 Effects of the model feature genes on the malignant phenotype of osteosarcoma cell lines

Cellular experiments were used to analyze the influences of the model genes on osteosarcoma progression of and to detect key genes involved in the malignant phenotype of osteosarcoma. Molecular detection results showed that IFITM3, APBB1IP, CDK6, GBP1, TAC4, and GJA5 expressions were relatively upregulated, and GJA5, CGREF1 mRNA levels were downregulated in osteosarcoma cell lines, with IFITM3 having the highest expression level (Fig. 7A–H). Wound healing and Transwell experiment demonstrated that silencing IFITM3 gene suppressed the migration (Fig. 7I–J) and invasion (Fig. 7K–L) of osteosarcoma cells, respectively. Overall, IFITM3 is a gene significantly affecting the prognostic outcomes of patients suffering from osteosarcoma and had a tangible regulatory effect on osteosarcoma cell lines, which verified the reliability of the prognostic model.Fig. 7 Effects of model feature genes on osteosarcoma cell. A–H: qRT-PCR results of model feature genes IFITM3, APBB1IP, CDK6, GBP1, TAC4, and GJA5. I–J: Scratch healing experiments and the quantified results of osteosarcoma cell line of si-IFITM3. K–L: Transwell experiment and the quantified results of osteosarcoma cell line of si-IFITM3.

Fig. 7

3.7 Distribution of the risk score in clinical characteristics

The risk score for samples in different clinical characteristics was computed, including metaststic, gender, age, cluster and status. The samples in metaststic, C3 and Dead had a higher risk score (Fig. 8A), while those in gender, age, metaststic had low risk score and better survival (Fig. 8B).Fig. 8 Clinical features analysis for both the high and low groups. A: The differences in the risk score between different clinicopathological groups in the TARGET-osteosarcoma cohort. B: KM curves for both the high-risk and low-risk groups among different clinicopathological groups in the TARGET-osteosarcoma cohort.

Fig. 8

3.8 Analysis on immune microenvironment between the two risk score groups

Immune microenvironment differences between the two groups were observed, with the group of low risk score having higher ESTIMATEScore, StromalScore, and ImmuneScore (Fig. 9A). MCP-Counter analysis demonstrated that the scores of 6 out of 10 types of immune cells in the low group were higher (Fig. 9B). The ssGSEA analysis showed that the scores of 15 of 28 types of immune cell scores were higher in the low group (Fig. 9C). Additionally, we calculated the infiltration of 22 immune cells in different risk groups using the CIBERSORT and obtained similar results to the MCP-Counter algorithm. When compared to the high-risk group, low-risk patients had higher CD8+T cells infiltration (Supplementary Fig. 1). Further study revealed that most immune cells from 28 immune cells showed negative correlation with the risk score (Fig. 9D).Fig. 9 Immune microenvironment analysis. A: The differences of the high group and low group in ESTIMATEScore, StromalScore, ImmuneScore. B: The differences of the two groups in 10 kinds of immunes score. C: The differences of the two groups in 28 kinds of immunes score. D: The correlation between 28 immune cells scores and the risk score. (*P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001).

Fig. 9

3.9 Pathways analysis for the two risk score groups

The ssGSEA analysis demonstrated a negative correlation between the immune pathways and risk score (Fig. 10A). GSEA analysis revealed that 12 pathways and 30 pathways were activated in the group with a high risk score in TARGET-osteosarcoma and GSE39058 dataset, respectively, and that most pathways were associated with cell cycle (Fig. 10B).Fig. 10 Pathways enrichment analysis. A: Results of correlation analysis between the risk score and KEGG pathways. NPRS showed a correlation higher than 0.4. B: By comparing high and low group, normalized enrichment scores of Hallmark pathways were shown in the heatmap.

Fig. 10

3.10 Clinical pathology features and the risk score to predict osteosarcoma survival

The decision tree was produced based on gender, age, metastatic and the risk score in the TARGET-osteosarcoma cohort. Five different risk subgroups (Fig. 11A) with significant differences in OS (Fig. 11B) were classified. Patients in C1, C2 and C3 were in low risk score group, while those in C4 and C5 were in high risk score group (Fig. 11C). The 5 subgroups had significance differences in survival status (Fig. 11D). Both metastatic and the risk score were validated as independent prognosis factors by univariate and multivariate Cox survival analysis (Fig. 11E and F). We also constructed a nomogram together with the risk score and metastatic to predict 1-, 3-, 5-year OS of osteosarcoma patients (Fig. 11G). The risk score had higher AUC (Fig. 10H) and the reliability of the nomogram was proven by the calibration curve (Fig. 11I). DCA further validated the risk score as an indicator for effectively predicting osteosarcoma prognosis in comparison to other clinical variables (Fig. 11J.).Fig. 11 Construction of nomogram. A: To optimize risk stratification, a survival decision tree was established according to the full-scale annotations of patients including age, metastatic, NPRS. B: The overall survival outcomes of the 5 risk subgroups were significantly different. C: The distribution of low and high group samples in the 5 risk subgroups. D: The status distribution in 5 risk subgroups. E: Univariate Cox analysis. F: Multivariate Cox analysis. G: The development of a nomogram. H: The survival prediction of the nomogram was the most accurate in comparison to other clinicopathological features. I: 1-, 3- and 5- years calibration curves for the nomogram. J: Decision curve of the nomogram.

Fig. 11

4 Discussion

Previous studies have confirmed crucial role of necroptosis in cancer development and metastasis [31,32]. In this work, 39 necroptosis-related genes in TARGET-osteosarcoma were screened to classify C1, C2, and C3 as three molecular subtypes of osteosarcoma. Patients in C1 and C2 had enhanced immune status and better survival but limited immunotherapy benefit. The prognostic outcomes of patients with osteosarcoma could be correctly evaluated by the 8 necroptosis-correlated prognostic signature genes. The prognostic outcomes were better for patients who had a low risk score, and the nomogram combining both the risk score and risk factors could produce an accurate prognostic prediction.

A total of 8 DEGs (ACTA2, APBB1IP, CDK6, CGREF1, GBP1, GJA5, IFITM3, and TAC4) were contained in the present model. The ROC curve and calibration curve showed that the model had an accurate prediction for osteosarcoma patients in the training and validation cohorts. Among all the 8 genes related to osteosarcoma prognosis, 4 genes (ACTA2, CDK6, IFITM3) were either involved in necroptosis-related pathogenesis of osteosarcoma or were important indicators of overall survival or relapse-free survival [[33], [34], [35]]. As an interferon-induced transmembrane protein, IFITM3 is high-expressed in the tumor region of gastric cancer. Depletion of IFITM3 can cause HGF-triggered inhibition of cell growth and migration by suppressing AKT/c-MYC signaling in gastric cancer [36]. Gan et al. found that knockdown of IFITM3 could impair the growth of the oral squamous cell carcinoma cell lines by inhibiting the proliferation of the cancer cells and inducing cell cycle arrest, senescence and apoptosis [37]. Consistently, we observed that silencing IFITM3 expression greatly affected the osteosarcoma cell migration and invasion. In addition, Tang et al. observed that patients with high-expressed ACTA2 have a noticeably more favorable prognosis than those with low-expressed ACTA2, suggesting that high-expressed ACTA2 in osteosarcoma patients with lung metastasis indicates a better prognosis [38]. Recent discoveries showed that CDK6 has a crucial function in the advancement of various human cancers, and that a downregulation in its expression is associated with a negative prognostic outcome. Based on an analysis of the GEPIA and STARBASE databases, Zhao et al. found that CDK6 is an oncogene in pancreatic cancer, and that CDK6 is significantly associated with PD-L1, PD-L2, and HAVCR2 (three immune checkpoints), immune cells infiltration, and immune biomarkers [39]. However, previous study was less concerned with the remaining 5 model genes in the osteosarcoma prognosis or their roles as new potential biomarkers for the cancer. Therefore, these key gene targets should be analyzed more extensively.

Past study has analyzed the roles of TME, in particular immune microenvironment [40,41]. Immune infiltration analysis demonstrated a negative relationship between the immune cell score and the risk score. Compared to the high-risk patients, those with a low risk had more immune cells and greater activation of immune pathways. Studies found that tumors could induce the generation of type II NKT cells and secrete IL-13 in turn, which causes MDSC aggregation in the tumor microenvironment, activates STAT6 signaling pathway and suppresses the function of CD8+ T cell [42,43]. Based on scRNA-seq transcriptomic data and immunohistochemical staining analysis, Xie et al. confirmed that EVI2B is present on CD8+ T cells derived from osteosarcoma patients. EVI2B enhances the expressions of granzyme A and K in CD8+ T cells, resulting in a robust cytotoxic effect against tumor cells [44]. NK cells within the TME of osteosarcoma are innate immune cells that display cytokine-secreting and cytotoxic properties along with cytotoxic and helper T cells [[45], [46], [47]]. Both CD8+ T cells and NK cells possess the ability to eliminate cancer cells through their cytotoxic activities, yet the synergistic relationship between them is intricate. MHC-I expression downregulated by specific cancer cells escape the detection of CD8+ T cells while enhancing NK cell activation by eliminating a significant inhibitory cue [48]. At present, the application of T cells and NK cells in tumor immunotherapy has attracted extensive research attention. Apart from hematological tumors, memory-like NK cells and modified NK cells have great potential in treating breast cancer, colorectal cancer, ovarian cancer, non-small cell lung cancer, liver cancer [49]. This suggested that the tumor microenvironment of low-risk patients may be more conducive to immune cell survival and function, which was strongly associated with a better prognosis for patients.

In conclusion, differentially expressed necroptosis-related genes were identified in this study to develop a prediction model for osteosarcoma prognosis, and the association between immune activities and the risk score was systematically assessed. However, the mechanisms of necroptosis in tumor immunology should be further investigated by future studies.

Ethical approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Ethical guidelines

Not applicable.

Availability of data and materials

The datasets generated and/or analyzed during the current study are available in the [GSE16088] repository, [https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE16088]; in [GSE12865] repository, [https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE12865]; in [GSE21257] repository, [https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE21257]; in [GSE39058] repository, [https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE39058].

Funding

10.13039/501100001809 National Natural Science Foundation of China (No.81772864 , No.81772861 ) supported this project.

CRediT authorship contribution statement

Yiying Bian: Writing – review & editing, Writing – original draft, Visualization, Validation, Software, Resources, Investigation, Data curation, Conceptualization. Jixiang Shi: Validation, Software, Resources, Methodology, Investigation. Ziyun Chen: Supervision, Software, Methodology, Investigation, Formal analysis, Data curation. Ji Fang: Visualization, Supervision, Resources, Formal analysis. Weidong Chen: Visualization, Supervision, Software, Project administration, Methodology. Yutong Zou: Visualization, Resources, Investigation, Formal analysis, Conceptualization. Hao Yao: Visualization, Software, Project administration, Methodology, Investigation. Jian Tu: Visualization, Validation, Supervision, Formal analysis, Data curation, Conceptualization. Yan Liao: Visualization, Supervision, Resources, Project administration, Investigation. Xianbiao Xie: Writing – original draft, Visualization, Validation, Software, Funding acquisition, Formal analysis. Jingnan Shen: Writing – review & editing, Visualization, Supervision, Project administration, Funding acquisition, Formal analysis.

Declaration of competing interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Abbreviations

(DEGs) differentially expressed genes

(LASSO) least absolute shrinkage and selection operator

(ROC) receiver operating characteristic

(GEO) Gene Expression Omnibus

(GSEA) gene set enrichment analysis

(TNF-α) tumor necrosis factor-α

(CDF) cumulative distribution function

(MDSC) myeloid suppressor cells

(TAM.M2) M2 subtypes of tumor-associated macrophages

(MCP-Counter) Microenvironment cell populations-counter

(RIP) receptor-interacting protein

(TIDE) Tumor Immune Dysfunction and Exclusion

(ESTIMATE) Estimation of STromal and Immune cells in MAlignant Tumors using Expression data

(AUC) area under the receiver operating characteristic curve

(KM) Kaplan-Meier

(CTL) cytotoxic T lymphocytes

Appendix A Supplementary data

The following are the Supplementary data to this article:Supplementary Fig. 1

Based on the CIBERSORT algorithm to assess differences in the levels of infiltration of 22 immune cells in patients from different osteosarcoma risk groups in the TARGET database.

Supplementary Fig. 1

Supplementary Table 1

Sample information based on osteosarcoma patients obtained from the TARGET database.

Supplementary Table 1

Supplementary Table 2

Target genes primers.

Supplementary Table 2

Acknowledgements

Not applicable.

Appendix A Supplementary data to this article can be found online at https://doi.org/10.1016/j.heliyon.2024.e35719.
==== Refs
References

1 Xie B. Li Y. Zhao R. Xu Y. Wu Y. Wang J. Identification of key genes and miRNAs in osteosarcoma patients with chemoresistance by bioinformatics analysis BioMed Res. Int. 2018 2018 4761064
2 Damron T.A. Ward W.G. Stewart A. Osteosarcoma, chondrosarcoma, and ewing's sarcoma: national cancer data base report Clin. Orthop. Relat. Res. 459 2007 40 47 17414166
3 Miller B.J. Cram P. Lynch C.F. Buckwalter J.A. Risk factors for metastatic disease at presentation with osteosarcoma: an analysis of the SEER database J. Bone Jt. Surg. Am. Vol. 95 13 2013 e89
4 Lindsey B.A. Markel J.E. Kleinerman E.S. Osteosarcoma overview Rheumatology and therapy 4 1 2017 25 43 27933467
5 Zhang Y. Yang J. Zhao N. Wang C. Kamar S. Zhou Y. Progress in the chemotherapeutic treatment of osteosarcoma Oncol. Lett. 16 5 2018 6228 6237 30405759
6 Degterev A. Huang Z. Boyce M. Li Y. Jagtap P. Mizushima N. Chemical inhibitor of nonapoptotic cell death with therapeutic potential for ischemic brain injury Nat. Chem. Biol. 1 2 2005 112 119 16408008
7 Wu Y.T. Tan H.L. Huang Q. Sun X.J. Zhu X. Shen H.M. zVAD-induced necroptosis in L929 cells depends on autocrine production of TNFα mediated by the PKC-MAPKs-AP-1 pathway Cell Death Differ. 18 1 2011 26 37 20539307
8 Sawai H. Ogiso H. Okazaki T. Differential changes in sphingolipids between TNF-induced necroptosis and apoptosis in U937 cells and necroptosis-resistant sublines Leuk. Res. 39 9 2015 964 970 26189109
9 Cao K. Tait S.W.G. Parkin inhibits necroptosis to prevent cancer Nat. Cell Biol. 21 8 2019 915 916 31358970
10 Park S. Hatanpaa K.J. Xie Y. Mickey B.E. Madden C.J. Raisanen J.M. The receptor interacting protein 1 inhibits p53 induction through NF-kappaB activation and confers a worse prognosis in glioblastoma Cancer Res. 69 7 2009 2809 2816 19339267
11 Wang Q. Chen W. Xu X. Li B. He W. Padilla M.T. RIP1 potentiates BPDE-induced transformation in human bronchial epithelial cells through catalase-mediated suppression of excessive reactive oxygen species Carcinogenesis 34 9 2013 2119 2128 23633517
12 Xiao H. Jensen P.E. Chen X. Elimination of osteosarcoma by necroptosis with graphene oxide-associated anti-HER2 antibodies Int. J. Mol. Sci. 20 18 2019
13 Li J. Yang Z. Li Y. Xia J. Li D. Li H. Cell apoptosis, autophagy and necroptosis in osteosarcoma treatment Oncotarget 7 28 2016 44763 44778 27007056
14 Eskandari A. Flamme M. Xiao Z. Suntharalingam K. The bulk osteosarcoma and osteosarcoma stem cell activity of a necroptosis-inducing nickel(II)-Phenanthroline complex Chembiochem : a European journal of chemical biology 21 19 2020 2854 2860 32415808
15 He P. Zhang Z. Liao W. Xu D. Fu M. Kang Y. Screening of gene signatures for rheumatoid arthritis and osteoarthritis based on bioinformatics analysis Mol. Med. Rep. 14 2 2016 1587 1593 27356888
16 Li Z. Wang Q. Chen G. Li X. Yang Q. Du Z. Integration of gene expression profile data to screen and verify hub genes involved in osteoarthritis BioMed Res. Int. 2018 2018 9482726
17 Lu Q.Y. Han Q.H. Li X. Li Z.C. Pan Y.T. Liu L. Analysis of differentially expressed genes between rheumatoid arthritis and osteoarthritis based on the gene co-expression network Mol. Med. Rep. 10 1 2014 119 124 24788818
18 Zou J. Chen L. Xu H. Unveiling and validation of a disulfidptosis determined prognostic model for osteosarcoma: new insights from prognosis to immunotherapy and chemotherapy Oncologie 25 4 2023 417 433
19 Chen Y. Huang L. Wei Z. Liu X. Chen L. Wang B. Development and validation of a nomogram model to predict the prognosis of intrahepatic cholangiocarcinoma Oncologie 24 2 2022 329 340
20 Yang J. Zhang J. Na S. Wang Z. Li H. Su Y. Integration of single-cell RNA sequencing and bulk RNA sequencing to reveal an immunogenic cell death-related 5-gene panel as a prognostic model for osteosarcoma Front. Immunol. 13 2022 994034
21 Xin S. Mao J. Duan C. Wang J. Lu Y. Yang J. Identification and quantification of necroptosis landscape on therapy and prognosis in kidney renal clear cell carcinoma Front. Genet. 13 2022 832046
22 Wilkerson M.D. Hayes D.N. ConsensusClusterPlus: a class discovery tool with confidence assessments and item tracking Bioinformatics 26 12 2010 1572 1573 20427518
23 Yang P. Chen W. Xu H. Yang J. Jiang J. Jiang Y. Correlation of CCL8 expression with immune cell infiltration of skin cutaneous melanoma: potential as a prognostic indicator and therapeutic pathway Cancer Cell Int. 21 1 2021 635 34844613
24 Fu J. Li K. Zhang W. Wan C. Zhang J. Jiang P. Large-scale public data reuse to model immunotherapy response and resistance Genome Med. 12 1 2020 21 32102694
25 Jiang P. Gu S. Pan D. Fu J. Sahu A. Hu X. Signatures of T cell dysfunction and exclusion predict cancer immunotherapy response Nat. Med. 24 10 2018 1550 1558 30127393
26 Geeleher P. Cox N. Huang R.S. pRRophetic: an R package for prediction of clinical chemotherapeutic response from tumor gene expression levels PLoS One 9 9 2014 e107468
27 Liberzon A. Birger C. Thorvaldsdóttir H. Ghandi M. Mesirov J.P. Tamayo P. The Molecular Signatures Database (MSigDB) hallmark gene set collection Cell systems 1 6 2015 417 425 26771021
28 Engebretsen S. Bohlin J. Statistical predictions with glmnet Clin. Epigenet. 11 1 2019 123
29 Liu F. Qian J. Ma C. MPscore: a novel predictive and prognostic scoring for progressive meningioma Cancers 13 5 2021
30 Chen B. Khodadoust M.S. Liu C.L. Newman A.M. Alizadeh A.A. Profiling tumor infiltrating immune cells with CIBERSORT Methods Mol. Biol. 1711 2018 243 259 29344893
31 Barbosa L.A. Fiuza P.P. Borges L.J. Rolim F.A. Andrade M.B. Luz N.F. RIPK1-RIPK3-MLKL-Associated necroptosis drives leishmania infantum killing in neutrophils Front. Immunol. 9 2018 1818 30154785
32 Jouan-Lanhouet S. Riquet F. Duprez L. Vanden Berghe T. Takahashi N. Vandenabeele P. Necroptosis, in vivo detection in experimental disease models Semin. Cell Dev. Biol. 35 2014 2 13 25160988
33 Yao F. Zhu Z.F. Wen J. Zhang F.Y. Zhang Z. Zhu L.Q. PODN is a prognostic biomarker and correlated with immune infiltrates in osteosarcoma Cancer Cell Int. 21 1 2021 381 34273970
34 Chen Y. Tang G. Qian H. Chen J. Cheng B. Zhou C. LncRNA LOC100129620 promotes osteosarcoma progression through regulating CDK6 expression, tumor angiogenesis, and macrophage polarization Aging 13 10 2021 14258 14276 34015762
35 Salas S. Jézéquel P. Campion L. Deville J.L. Chibon F. Bartoli C. Molecular characterization of the response to chemotherapy in conventional osteosarcomas: predictive value of HSD17B10 and IFITM2 Int. J. Cancer 125 4 2009 851 860 19449377
36 Chu P.Y. Huang W.C. Tung S.L. Tsai C.Y. Chen C.J. Liu Y.C. IFITM3 promotes malignant progression, cancer stemness and chemoresistance of gastric cancer by targeting MET/AKT/FOXO3/c-MYC axis Cell Biosci. 12 1 2022 124 35941699
37 Gan C.P. Sam K.K. Yee P.S. Zainal N.S. Lee B.K.B. Abdul Rahman Z.A. IFITM3 knockdown reduces the expression of CCND1 and CDK4 and suppresses the growth of oral squamous cell carcinoma cells Cell. Oncol. 42 4 2019 477 490
38 Tang L. Hu H. Zhou Y. Huang Y. Wang Y. Zhang Y. Expression and clinical significance of ACTA2 in osteosarcoma tissue Oncologie 24 4 2022 913 9125
39 Zhao Y.X. Xu B.W. Wang F.Q. Jiang F.Y. Xu J.W. Yu D.X. nc-RNA-mediated high expression of CDK6 correlates with poor prognosis and immune infiltration in pancreatic cancer Cancer Med. 12 4 2023 5110 5123 36457244
40 Xu D. Wang Y. Zhou K. Wu J. Zhang Z. Zhang J. Identification of an extracellular vesicle-related gene signature in the prediction of pancreatic cancer clinical prognosis Biosci. Rep. 40 12 2020
41 Li X.Y. Zhang L.Y. Li X.Y. Yang X.T. Su L.X. A pyroptosis-related gene signature for predicting survival in glioblastoma Front. Oncol. 11 2021 697198
42 Terabe M. Matsui S. Noben-Trauth N. Chen H. Watson C. Donaldson D.D. NKT cell-mediated repression of tumor immunosurveillance by IL-13 and the IL-4R-STAT6 pathway Nat. Immunol. 1 6 2000 515 520 11101874
43 Terabe M. Swann J. Ambrosino E. Sinha P. Takaku S. Hayakawa Y. A nonclassical non-Valpha14Jalpha18 CD1d-restricted (type II) NKT cell is sufficient for down-regulation of tumor immunosurveillance J. Exp. Med. 202 12 2005 1627 1633 16365146
44 Xie T. Feng W. He M. Zhan X. Liao S. He J. Analysis of scRNA-seq and bulk RNA-seq demonstrates the effects of EVI2B or CD361 on CD8(+) T cells in osteosarcoma Exp. Biol. Med. (Maywood, NJ, U. S.) 248 2 2023 130 145
45 Mason N.J. Comparative immunology and immunotherapy of canine osteosarcoma Adv. Exp. Med. Biol. 1258 2020 199 221 32767244
46 Cao H. Xue Y. Wang F. Li G. Zhen Y. Guo J. Identification of prognostic molecular subtypes and model based on CD8+ T cells for lung adenocarcinoma Biocell 48 3 2024 473 490
47 Li G. Li F. Wei N. Jia Q. Designing a risk prognosis model based on natural killer cell-linked genes to accurately evaluate the prognosis of gastric cancer Biocell 47 9 2023 2081 20899
48 Myers J.A. Miller J.S. Exploring the NK cell platform for cancer immunotherapy Nat. Rev. Clin. Oncol. 18 2 2021 85 100 32934330
49 Li Y. Yin J. Li T. Huang S. Yan H. Leavenworth J. NK cell-based cancer immunotherapy: from basic biology to clinical application Sci. China Life Sci. 58 12 2015 1233 1245 26588912
