A Novel Pyroptosis-Related Signature for Predicting Prognosis and Indicating Immune Microenvironment Features in Osteosarcoma

Development and experimental validation of hypoxia-related gene signatures for osteosarcoma diagnosis and prognosis based on WGCNA and machine learning

Osteosarcoma (OS) is the most common primary malignant tumour of the bone with high mortality. Here, we comprehensively analysed the hypoxia signalling in OS and further constructed novel hypoxia-related gene signatures for OS prediction and prognosis. This study employed Gene Set Enrichment Analysis (GSEA), Weighted correlation network analysis (WGCNA) and Least absolute shrinkage and selection operator (LASSO) analyses to identify Stanniocalcin 2 (STC2) and Transmembrane Protein 45A (TMEM45A) as the diagnostic biomarkers, which further assessed by Receiver Operating Characteristic (ROC), decision curve analysis (DCA), and calibration curves in training and test dataset. Univariate and multivariate Cox regression analyses were used to construct the prognostic model. STC2 and metastasis were devised to forge the OS risk model. The nomogram, risk score, Kaplan Meier plot, ROC, DCA, and calibration curves results certified the excellent performance of the prognostic model. The expression level of STC2 and TMEM45A was validated in external datasets and cell lines. In immune cell infiltration analysis, cancer-associated fibroblasts (CAFs) were significantly higher in the low-risk group. And the immune infiltration of CAFs was negatively associated with the expression of STC2 (P < 0.05). Pan-cancer analysis revealed that the expression level of STC2 was significantly higher in Esophageal carcinoma (ESCA), Head and Neck squamous cell carcinoma (HNSC), Kidney renal clear cell carcinoma (KIRC), Lung squamous cell carcinoma (LUSC), and Stomach adenocarcinoma (STAD). Additionally, the higher expression of STC2 was associated with the poor outcome in those cancers. In summary, this study identified STC2 and TMEM45A as novel markers for the diagnosis and prognosis of osteosarcoma, and STC2 was shown to correlate with immune infiltration of CAFs negatively.

Exploring the Role of Neutrophil-Related Genes in Osteosarcoma via an Integrative Analysis of Single-Cell and Bulk Transcriptome

BACKGROUND

The involvement of neutrophil-related genes (NRGs) in patients with osteosarcoma (OS) has not been adequately explored. In this study, we aimed to examine the association between NRGs and the prognosis as well as the tumor microenvironment of OS.

METHODS

The OS data were obtained from the TARGET-OS and GEO database. Initially, we extracted NRGs by intersecting 538 NRGs from single-cell RNA sequencing (scRNA-seq) data between aneuploid and diploid groups, as well as 161 up-regulated differentially expressed genes (DEGs) from the TARGET-OS datasets. Subsequently, we conducted Least Absolute Shrinkage and Selection Operator (Lasso) analyses to identify the hub genes for constructing the NRG-score and NRG-signature. To assess the prognostic value of the NRG signatures in OS, we performed Kaplan-Meier analysis and generated time-dependent receiver operating characteristic (ROC) curves. Gene enrichment analysis (GSEA) and gene set variation analysis (GSVA) were utilized to ascertain the presence of tumor immune microenvironments (TIMEs) and immunomodulators (IMs). Additionally, the KEGG neutrophil signaling pathway was evaluated using ssGSEA. Subsequently, PCR and IHC were conducted to validate the expression of hub genes and transcription factors (TFs) in K7M2-induced OS mice.

RESULTS

FCER1G and C3AR1 have been identified as prognostic biomarkers for overall survival. The findings indicate a significantly improved prognosis for OS patients. The effectiveness and precision of the NRG signature in prognosticating OS patients were validated through survival ROC curves and an external validation dataset. The results clearly demonstrate that patients with elevated NRG scores exhibit decreased levels of immunomodulators, stromal score, immune score, ESTIMATE score, and infiltrating immune cell populations. Furthermore, our findings substantiate the potential role of SPI1 as a transcription factor in the regulation of the two central genes involved in osteosarcoma development. Moreover, our analysis unveiled a significant correlation and activation of the KEGG neutrophil signaling pathway with FCER1G and C3AR1. Notably, PCR and IHC demonstrated a significantly higher expression of C3AR1, FCER1G, and SPI1 in Balb/c mice induced with K7M2.

CONCLUSIONS

Our research emphasizes the significant contribution of neutrophils within the TIME of osteosarcoma. The newly developed NRG signature could serve as a good instrument for evaluating the prognosis and therapeutic approach for OS.

Identification of a novel monocyte/macrophage-related gene signature for predicting survival and immune response in acute myeloid leukemia

Acute myeloid leukemia (AML) is a heterogeneous hematological tumor with poor immunotherapy effect. This study was to develop a monocyte/macrophage-related prognostic risk score (MMrisk) and identify new therapeutic biomarkers for AML. We utilized differentially expressed genes (DEGs) in combination with single-cell RNA sequencing to identify monocyte/macrophage-related genes (MMGs). Eight genes were selected for the construction of a MMrisk model using univariate Cox regression analysis and LASSO regression analysis. We then validated the MMrisk on two GEO datasets. Lastly, we investigated the immunologic characteristics and advantages of immunotherapy and potential targeted drugs for MMrisk groups. Our study identified that the MMrisk is composed of eight MMGs, including HOPX, CSTB, MAP3K1, LGALS1, CFD, MXD1, CASP1 and BCL2A1. The low MMrisk group survived longer than high MMrisk group (P < 0.001). The high MMrisk group was positively correlated with B cells, plasma cells, CD4 memory cells, Mast cells, CAFs, monocytes, M2 macrophages, Endothelial, tumor mutation, and most immune checkpoints (PD1, Tim-3, CTLA4, LAG3). Furthermore, drug sensitivity analysis showed that AZD.2281, Axitinib, AUY922, ABT.888, and ATRA were effective in high-risk MM patients. Our research shows that MMrisk is a potential biomarker which is helpful to identify the molecular characteristics of AML immunology.

ANXA5: A Key Regulator of Immune Cell Infiltration in Hepatocellular Carcinoma

BACKGROUND Hepatocellular carcinoma (HCC) poses a significant threat to human life and is the most prevalent form of liver cancer. The intricate interplay between apoptosis, a common form of programmed cell death, and its role in immune regulation stands as a crucial mechanism influencing tumor metastasis. MATERIAL AND METHODS Utilizing HCC samples from the TCGA database and 61 anoikis-related genes (ARGs) sourced from GeneCards, we analyzed the relationship between ARGs and immune cell infiltration in HCC. Subsequently, we identified long non-coding RNAs (lncRNAs) associated with ARGs, using the least absolute shrinkage and selection operator (LASSO) regression analysis to construct a robust prognostic model. The predictive capabilities of the model were then validated through examination in a single-cell dataset. RESULTS Our constructed prognostic model, derived from lncRNAs linked to ARGs, comprised 11 significant lncRNAs: NRAV, MCM3AP-AS1, OTUD6B-AS1, AC026356.1, AC009133.1, DDX11-AS1, AC108463.2, MIR4435-2HG, WARS2-AS1, LINC01094, and HCG18. The risk score assigned to HCC samples demonstrated associations with immune indicators and the infiltration of immune cells. Further, we identified Annexin A5 (ANXA5) as the pivotal gene among ARGs, with it exerting a prominent role in regulating the lncRNA gene signature. Our validation in a single-cell database elucidated the involvement of ANXA5 in immune cell infiltration, specifically in the regulation of mononuclear cells. CONCLUSIONS This study delves into the intricate correlation between ARGs and immune cell infiltration in HCC, culminating in the development of a novel prognostic model reliant on 11 ARGs-associated lncRNAs. Furthermore, our findings highlight ANXA5 as a promising target for immune regulation in HCC, offering new perspectives for immune therapy in the context of HCC.

A programmed cell death-related gene signature to predict prognosis and therapeutic responses in liver hepatocellular carcinoma

BACKGROUND

Programmed cell death (PCD) functions critically in cancers and PCD-related genes are associated with tumor microenvironment (TME), prognosis and therapeutic responses of cancer patients. This study stratified hepatocellular carcinoma (HCC) patients and develop a prognostic model for predicting prognosis and therapeutic responses.

METHODS

Consensus clustering analysis was performed to subtype HCC patients in The Cancer Genome Atlas (TCGA) database. Differentially expressed genes (DEGs) among the subtypes were filtered and subjected to the least absolute shrinkage and selection operator (LASSO) regression analysis and univariate Cox regression analysis to filter prognostic genes. A PCD-related prognostic gene signature in TCGA was constructed and validated in ICGC-LIRI-JP and GSE14520 datasets. TME was analyzed using CIBERSORT, MCP-counter, TIMER and EPIC algorithms. Drug sensitivity was predicted by oncoPredict package. Spearman analysis was used to detect correlation.

RESULTS

Four molecular subtypes were categorized based on PCD-related genes. Subtype C1 showed the poorest prognosis, the most infiltration of Fibroblasts, dentritic cell (DC) and cancer-associated fibroblasts (CAFs), and the highest TIDE score. C4 had a better prognosis survival outcome, and lowest immune cell infiltration. The survival outcomes of C2 and C3 were intermediate. Next, a total of 69 co-DEGs were screened among the four subtypes and subsequently we identified five prognostic genes (MCM2, SPP1, S100A9, MSC and EPO) for developing the prognostic model. High-risk patients not only had unfavorable prognosis, higher clinical stage and grade, and more inflammatory pathway enrichment, but also possessed higher possibility of immune escape and were more sensitive to Cisplatin and 5. Fluorouracil. The robustness of the prognostic model was validated in external datasets.

CONCLUSION

This study provides new insights into clinical subtyping and the PCD-related prognostic signature may serve as a useful tool to predict prognosis and guide treatments for patients with HCC.

Gasdermins: a dual role in pyroptosis and tumor immunity

The gasdermin (GSDM) protein family plays a pivotal role in pyroptosis, a process critical to the body's immune response, particularly in combatting bacterial infections, impeding tumor invasion, and contributing to the pathogenesis of various inflammatory diseases. These proteins are adept at activating inflammasome signaling pathways, recruiting immune effector cells, creating an inflammatory immune microenvironment, and initiating pyroptosis. This article serves as an introduction to the GSDM protein-mediated pyroptosis signaling pathways, providing an overview of GSDMs' involvement in tumor immunity. Additionally, we explore the potential applications of GSDMs in both innovative and established antitumor strategies.

Novel pyroptosis-related lncRNAs and ceRNAs predict osteosarcoma prognosis and indicate immune microenvironment signatures

Objective

To study pyroptosis-related biomarkers that are associated with the prognosis and immune microenvironment characteristics of osteosarcoma (OS). The goal is to establish a foundation for the prognosis and treatment of OS.

Methods

We retrieved transcriptome and clinical data from The Cancer Genome Atlas (TCGA) database for 88 OS patients. Using this data, we constructed a prognostic model to identify pyroptosis-related genes (PRGs) associated with OS prognosis. To further explore the biological function of these PRGs, we performed enrichment analysis. To identify pyroptosis-related long non-coding RNAs (PRLncs) associated with the prognosis of OS, we performed co-expression analysis. Subsequently, a risk prognostic model was constructed using these PRLncs to generate a risk score, termed as PRLncs-score, thereby obtaining PRLncs associated with the prognosis of OS. The accuracy of the prognostic model was verified through survival analysis, risk curve, independent prognostic analysis, receiver operating characteristic (ROC) curve, difference analysis between high- and low-risk groups, and clinical correlation analysis. And to determine whether PRLncs-score is independent prognostic factor for OS. In addition, we further conducted external and internal validation for the risk prognosis model. Further analyses of immune cell infiltration and tumor microenvironment were performed. A pyroptosis-related competitive endogenous RNA (PRceRNA) network was constructed to obtain PRceRNAs associated with the prognosis of OS and performed gene set enrichment analysis (GSEA) on PRceRNA genes.

Results

We obtained five PRGs (CHMP4C, BAK1, GSDMA, CASP1, and CASP6) that predicted OS prognosis and seven PRLncs (AC090559.1, AP003119.2, CARD8-AS1, AL390728.4, SATB2-AS1, AL133215.2, and AC009495.3) and one PRceRNA (CARD8-AS1-hsa-miR-21-5p-IL1B) that predicted OS prognosis and indicated characteristics of the OS immune microenvironment. The PRLncs-score, in combination with other clinical features, was established as an independent prognostic factor for OS patients. Subsequent scrutiny of the tumor microenvironment and immune infiltration indicated that patients with low-PRLncs-scores were associated with reduced metastatic risk, improved survival rates, heightened levels of immune cells and stroma, and increased immune activity compared to those with high-PRLncs-scores.

Conclusion

The study's findings offer insight into the prognosis of OS and its immune microenvironment, and hold promise for improving early diagnosis and immunotherapy.

Ferroptosis-related lncRNAs guiding osteosarcoma prognosis and immune microenvironment

OBJECTIVE

To investigate the ferroptosis-related long non-coding RNAs (FRLncs) implicated in influencing the prognostic and immune microenvironment in osteosarcoma (OS), and to establish a foundational framework for informing clinical decision making pertaining to OS management.

METHODS

Transcriptome data and clinical data pertaining to 86 cases of OS, the GSE19276, GSE16088 and GSE33382 datasets, and a list of ferroptosis-related genes (FRGs) were used to establish a risk prognostic model through comprehensive analysis. The identification of OS-related differentially expressed FRGs was achieved through an integrated analysis encompassing the aforementioned 86 OS transcriptome data and the GSE19276, GSE16088 and GSE33382 datasets. Concurrently, OS-related FRLncs were ascertained via co-expression analysis. To establish a risk prognostic model for OS, Univariate Cox regression analysis and Lasso Cox regression analysis were employed. Subsequently, a comprehensive evaluation was conducted, comprising risk curve analysis, survival analysis, receiver operating characteristic curve analysis and independent prognosis analysis. Model validation with distinct clinical subgroups was performed to assess the applicability of the risk prognostic model to diverse patient categories. Moreover, single sample gene set enrichment analysis (ssGSEA) was conducted to investigate variations in immune cell populations and immune functions within the context of the risk prognostic model. Furthermore, an analysis of immune checkpoint differentials yielded insights into immune checkpoint-related genes linked to OS prognosis. Finally, the risk prognosis model was verified by dividing the samples into train group and test group.

RESULTS

We identified a set of seven FRLncs that exhibit potential as prognostic markers and influence factors of the immune microenvironment in the context of OS. This ensemble encompasses three high-risk FRLncs, denoted as APTR, AC105914.2 and AL139246.5, alongside four low-risk FRLncs, designated as DSCR8, LOH12CR2, AC027307.2 and AC025048.2. Furthermore, our analysis revealed notable down-regulation in the high-risk group across four distinct immune cell types, namely neutrophils, natural killer cells, plasmacytoid dendritic cells and tumor-infiltrating lymphocytes. This down-regulation was also reflected in four key immune functions, antigen-presenting cell (APC)-co-stimulation, checkpoint, cytolytic activity and T cell co-inhibition. Additionally, we identified seven immune checkpoint-associated genes with significant implications for OS prognosis, including CD200R1, HAVCR2, LGALS9, CD27, LAIR1, LAG3 and TNFSF4.

CONCLUSION

The findings of this study have identified FRLncs capable of influencing OS prognosis and immune microenvironment, as well as immune checkpoint-related genes that are linked to OS prognosis. These discoveries establish a substantive foundation for further investigations into OS survival and offer valuable insights for informing clinical decision making in this context.

Large Scale Comparative Deconvolution Analysis of the Canine and Human Osteosarcoma Tumor Microenvironment Uncovers Conserved Clinically Relevant Subtypes

Osteosarcoma is a relatively rare but aggressive cancer of the bones with a shortage of effective biomarkers. Although less common in humans, Osteosarcomas are fairly common in adult pet dogs and have been shown to share many similarities with their human analogs. In this work, we analyze bulk transcriptomic data of 213 primary and 100 metastatic Osteosarcoma samples from 210 pet dogs enrolled in nation-wide clinical trials to uncover three Tumor Microenvironment (TME)-based subtypes: Immune Enriched (IE), Immune Enriched Dense Extra-Cellular Matrix-like (IE-ECM) and Immune Desert (ID) with distinct cell type compositions, oncogenic pathway activity and chromosomal instability. Furthermore, leveraging bulk transcriptomic data of canine primary tumors and their matched metastases from different sites, we characterize how the Osteosarcoma TME evolves from primary to metastatic disease in a standard of care clinical setting and assess its overall impact on clinical outcomes of canines. Most importantly, we find that TME-based subtypes of canine Osteosarcomas are conserved in humans and predictive of progression free survival outcomes of human patients, independently of known prognostic biomarkers such as presence of metastatic disease at diagnosis and percent necrosis following chemotherapy. In summary, these results demonstrate the power of using canines to model the human Osteosarcoma TME and discover novel biomarkers for clinical translation.

Construction of a T-cell exhaustion-related gene signature for predicting prognosis and immune response in hepatocellular carcinoma

BACKGROUND

Hepatocellular carcinoma (HCC) is a heterogeneous malignancy with a rising prevalence worldwide. Immunotherapy has been shown to improve treatment outcomes for HCC. We aimed to construct a T-cell exhaustion-related gene prognostic model (TEXPM) for HCC and to elucidate the immunologic characteristics and advantages of immunotherapy in T-cell exhaustion-Related Gene-defined HCC groups.

METHODS

Single-cell RNA sequencing data were used in conjunction with TCGA Differentially expressed genes (DEGs) to screen for T-cell exhaustion-Related Genes (TEXGs) for subsequent evaluation. Using univariate Cox regression analysis and LASSO regression analysis, five genes (FTL, GZMA, CD14, NPC2, and IER3) were subsequently selected for the construction of a TEXPM. Then, we evaluated the immunologic characteristics and advantages of immunotherapy in groups identified by TEXPM.

RESULTS

The TEXPM was formed with FTL, GZMA, CD14, NPC2, and IER3. The results of the training and validation team studies were consistent, with the low TEXPM group surviving longer than the high TEXPM group (P < 0.001). Multivariate Cox regression analysis demonstrated that TEXPM (HR: 2.347, 95%CI: 1.844-2.987; HR: 2.172, 95% CI: 1.689-2.793) was an independent prognostic variable for HCC patients. The low-TEXPM group was linked to active immunity, less aggressive phenotypes, strong infiltration of CD8+ T cells, CD4 + T cells, and M1 macrophages, and a better response to ICI treatment. A high TEXPM group, on the other hand, was associated with suppressive immunity, more aggressive phenotypes, a significant infiltration of B cells, M0 macrophages, and M2 macrophages, and a reduced response to ICI treatment. FTL is an independent prognostic variable in HCC patients and the knockdown of FTL can affect the biological behavior of hepatocellular carcinoma cells.

CONCLUSIONS

TEXPM is a promising prognostic biomarker connected to the immune system. Differentiating immunological and molecular features and predicting patient outcomes may be facilitated by TEXPM grouping. Furthermore, the expression of FTL was found to be an independent prognostic factor for HCC. Knockdown of FTL significantly inhibited proliferation, migration, and invasive activity in liver cancer cells.

A glycometabolic gene signature associating with immune infiltration and chemosensitivity and predicting the prognosis of patients with osteosarcoma

Background

Accumulating evidence has suggested that glycometabolism plays an important role in the pathogenesis of tumorigenesis. However, few studies have investigated the prognostic values of glycometabolic genes in patients with osteosarcoma (OS). This study aimed to recognize and establish a glycometabolic gene signature to forecast the prognosis, and provide therapeutic options for patients with OS.

Methods

Univariate and multivariate Cox regression, LASSO Cox regression, overall survival analysis, receiver operating characteristic curve, and nomogram were adopted to develop the glycometabolic gene signature, and further evaluate the prognostic values of this signature. Functional analyses including Gene Ontology (GO), kyoto encyclopedia of genes and genomes analyses (KEGG), gene set enrichment analysis, single-sample gene set enrichment analysis (ssGSEA), and competing endogenous RNA (ceRNA) network, were used to explore the molecular mechanisms of OS and the correlation between immune infiltration and gene signature. Moreover, these prognostic genes were further validated by immunohistochemical staining.

Results

A total of four genes including PRKACB, SEPHS2, GPX7, and PFKFB3 were identified for constructing a glycometabolic gene signature which had a favorable performance in predicting the prognosis of patients with OS. Univariate and multivariate Cox regression analyses revealed that the risk score was an independent prognostic factor. Functional analyses indicated that multiple immune associated biological processes and pathways were enriched in the low-risk group, while 26 immunocytes were down-regulated in the high-risk group. The patients in high-risk group showed elevated sensitivity to doxorubicin. Furthermore, these prognostic genes could directly or indirectly interact with other 50 genes. A ceRNA regulatory network based on these prognostic genes was also constructed. The results of immunohistochemical staining showed that SEPHS2, GPX7, and PFKFB3 were differentially expressed between OS tissues and adjacent normal tissues.

Conclusion

The preset study constructed and validated a novel glycometabolic gene signature which could predict the prognosis of patients with OS, identify the degree of immune infiltration in tumor microenvironment, and provide guidance for the selection of chemotherapeutic drugs. These findings may shed new light on the investigation of molecular mechanisms and comprehensive treatments for OS.

Anoikis patterns exhibit distinct prognostic and immune landscapes in Osteosarcoma

OBJECTIVES

Osteosarcoma is highly aggressive and prone to metastasis, with a poor prognosis. Increasing evidence identified anoikis has a critical effect in tumor metastasis and invasion. However, the prognostic value of anoikis-related genes (ANRGs) in osteosarcoma and their role in the immune landscape of osteosarcoma remain unclear.

METHODS

The RNA sequencing and clinical data of patients with osteosarcoma were extracted from the TARGET and GEO databases, and ANRGs were identified from the GeneCards database. Unsupervised clustering analysis was employed to identify anoikis-related patterns. The ESTIMATE, TIMER and ssGSEA algorithms were used to assess the immune microenvironment of different subtypes. A prognostic signature based on the identified ANRGs was constructed via univariate, LASSO and multivariate Cox regression analyses. KEGG, GO and GSEA were used for functional enrichment of genes associated with different risk subtypes. qPCR, WB and IHC were used to validate the expression of candidate genes.

RESULTS

Two anoikis-related patterns with distinct clinical features and immune statuses were identified based on prognosis-related ANRGs. Cluster 2 had more active immunogenicity and a better prognosis than Cluster 1. Subsequently, we developed and validated an anoikis prognostic signature demonstrating excellent predictive ability for the prognosis of osteosarcoma. Anoikis risk score was positively associated with osteosarcoma metastasis and was identified as an independent prognostic marker. Additionally, a nomogram was established to predict the 3- and 5-year survival probability of patients with osteosarcoma. Functional enrichment analysis revealed that immune dysregulation was correlated with poor prognosis. Besides, patients in the low-risk group had higher infiltration levels of immune cells and more active immune function than patients in the high-risk group. Drug sensitivity analysis revealed several chemotherapeutic agents for the treatment of different subtypes of osteosarcoma.

CONCLUSION

Our study demonstrated the role of ANRGs in osteosarcoma progression, providing insights into clinical decision making in osteosarcoma.

Identification and validation of an anoikis-related lncRNA signature to predict prognosis and immune landscape in osteosarcoma

Background

Anoikis is a specialized form of programmed apoptosis that occurs in two model epithelial cell lines and plays an important role in tumors. However, the prognostic value of anoikis-related lncRNA (ARLncs) in osteosarcoma (OS) has not been reported.

Methods

Based on GTEx and TARGET RNA sequencing data, we carried out a thorough bioinformatics analysis. The 27 anoikis-related genes were obtained from the Gene Set Enrichment Analysis (GSEA). Univariate Cox regression and least absolute shrinkage and selection operator (LASSO) analysis were successively used to screen for prognostic-related ARLncs. To create the prognostic signature of ARLncs, we performed multivariate Cox regression analysis. We calculated the risk score based on the risk coefficient, dividing OS patients into high- and low-risk subgroups. Additionally, the relationship between the OS immune microenvironment and risk prognostic models was investigated using function enrichment, including Gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), single-sample gene set enrichment analysis (ssGSEA), and GSEA analysis. Finally, the potential effective drugs in OS were found by immune checkpoint and drug sensitivity screening.

Results

A prognostic signature consisting of four ARLncs (AC079612.1, MEF2C-AS1, SNHG6, and TBX2-AS1) was constructed. To assess the regulation patterns of anoikis-related lncRNA genes, we created a risk score model. According to a survival analysis, high-risk patients have a poor prognosis as they progress. By using immune functional analysis, the lower-risk group demonstrated the opposite effects compared with the higher-risk group. GO and KEGG analysis showed that the ARLncs pathways and immune-related pathways were enriched. Immune checkpoints and drug sensitivity analysis might be used to determine the better effects of the higher group.

Conclusion

We identified a novel prognostic model based on a four-ARLncs signature that might serve as potential prognostic indicators that can be used to predict the prognosis of OS patients, and immunotherapy and drugs that may contribute to improving the overall survival of OS patients and advance our understanding of OS.

Development and validation of a cuproptosis-related lncRNA model correlated to the cancer-associated fibroblasts enable the prediction prognosis of patients with osteosarcoma

Background

Osteosarcoma is the most common primary pediatric and adolescent bone malignancy. An imbalance in copper homeostasis caused by copper ion accumulation could increase intracellular toxicity and regulate cancer cell growth. This study aimed to identify long non-coding RNAs (lncRNAs) associated with cuproptosis to predict prognosis and target drug use to improve patient survival.

Methods

RNA sequencing and relevant clinical information of ninety-three osteosarcoma patients were obtained from the TARGET database. We then identified thirteen prognostic cuproptosis-related lncRNAs(CRLncs) using coexpression and univariate Cox analyses. The prognostic risk model with three CRLncs was constructed using the least absolute shrinkage and selection operator (LASSO) and multivariate Cox regression analysis. Patients were divided into low-risk and high-risk subgroups using the median risk score. The tumor microenvironment (TME) and immune status of identified subgroups were analyzed using ESTIMATE, CIBERSORT, MCP-counter, xCELL, EPIC, and ssGSEA analyses. Functional analyses were conducted to elucidate the underlying mechanisms, including GO, KEGG, GSVA, and GSEA analyses. Also, the relationships between the model, tumor immunity, and drug sensitivity were explored. Lastly, the expression level of ZNF37BP, AL353759.1, and AC005034.5 was validated in vitro.

Results

We constructed a model containing three CRLncs (ZNF37BP, AL353759.1, and AC005034.5) and validated its excellent prognostic and predictive power. The AUC curves for 1-year, 3-year, and 5-year survival probabilities were 0.76, 0.84, and 0.89, respectively. Patients in the high-risk group had a shorter overall survival (OS) time than those in the low-risk. The stroma score and cancer-associated fibroblasts (CAFs) were significantly higher in the low-risk group. Immune cells such as T cells CD4 naive, T cells gamma delta, NK cells resting, dendritic cells resting, and mast cells activated were significantly upregulated in the high-risk group. Based on functional analyses, the PI3K-Akt pathway was identified as a critical metabolic pathway in osteosarcoma. Additionally, we obtained three potentially effective drugs for OS: erlotinib, MP470, and WH-4-023 targeting the PI3K-Akt pathway. The expression level of ZNF37BP was significantly elevated in OS cell lines than in normal osteoblast hFOB1.19 cells, and that of ATP7A, LIPT1, AL353759.1, and AC005034.5 were decreased considerably in OS cell lines.

Conclusion

Cuproptosis-related lncRNAs are correlated with the CAFs of osteosarcoma, and this could serve as a foundation for OS survival prediction and treatment.

Systematic Analysis of a Pyroptosis-Related Signature to Predict the Prognosis and Immune Microenvironment of Lower-Grade Glioma

Current treatments for lower-grade glioma (LGG) do not effectively improve life expectancy rates, and this is a major global health concern. Improving our knowledge of this disease will ultimately help to improve prevention, accurate prognosis, and treatment strategies. Pyroptosis is an inflammatory form of regulated cell death, which plays an important role in tumor progression and occurrence. There is still a lack of effective markers to evaluate the prognosis of LGG patients. We collected paraffin-embedded tissue samples and prognostic information from 85 patients with low-grade gliomas and fabricated them into a tissue microarray. Combining data from public databases, we explored the relationship between pyroptosis-related genes (PRGs) and the prognoses of patients with LGG and investigated their correlations with the tumor microenvironment (TME) by means of machine learning, single-cell, immunohistochemical, nomogram, GSEA, and Cox regression analyses. We developed a six-gene PRG-based prognostic model, and the results have identified CASP4 as an effective marker for LGG prognosis predictions. Furthermore, the effects on immune cell infiltration may also provide guidance for future immunotherapy strategies.

Identification of pyroptosis-related genes and long non-coding RNAs signatures in osteosarcoma

Osteosarcoma is a highly malignant tumor, with very high disability and fatality rates. However, the overall prognosis is not optimistic. Pyroptosis is a newly discovered cell death modality accompanied by inflammation, which is closely related to varieties of cancers. In this study, the RNA-seq data were downloaded from public databases, the differences in the expression of the pyroptosis-related genes (PRGs) were identified, and the six PRGs signature was established through the univariate and LASSO Cox analysis. The patients were grouped according to the PRGs signature, and the prognosis between the two groups was further compared. In addition, a ten pyroptosis-related lncRNAs (PRLs) prognostic signature was also constructed. Through functional analysis of the differentially expressed genes (DEGs), the immune-related pathways were found to be enriched. The Pearson correlation analysis showed a strong correlation between the pyroptosis-related biomarkers. Finally, we identified a promising biomarker, CHMP4C, which is highly expressed in osteosarcoma. Overexpression of CHMP4C promoted the proliferation, migration and invasion of the osteosarcoma cell. Our results thus provide new evidence for exploring prognostic biomarkers and therapeutic targets of osteosarcoma.

The emerging role of pyroptosis in pediatric cancers: from mechanism to therapy

Pediatric cancers are the driving cause of death for children and adolescents. Due to safety requirements and considerations, treatment strategies and drugs for pediatric cancers have been so far scarcely studied. It is well known that tumor cells tend to progressively evade cell death pathways, which is known as apoptosis resistance, one of the hallmarks of cancer, dominating tumor drug resistance. Recently, treatments targeting nonapoptotic cell death have drawn great attention. Pyroptosis, a newly specialized form of cell death, acts as a critical physiological regulator in inflammatory reaction, cell development, tissue homeostasis and stress response. The action in different forms of pyroptosis is of great significance in the therapy of pediatric cancers. Pyroptosis could be induced and consequently modulate tumorigenesis, progression, and metastasis if treated with local or systemic therapies. However, excessive or uncontrolled cell death might lead to tissue damage, acute inflammation, or even cytokine release syndrome, which facilitates tumor progression or recurrence. Herein, we aimed to describe the molecular mechanisms of pyroptosis, to highlight and discuss the challenges and opportunities for activating pyroptosis pathways through various oncologic therapies in multiple pediatric neoplasms, including osteosarcoma, neuroblastoma, leukemia, lymphoma, and brain tumors.

Insight into the interplay between mitochondria-regulated cell death and energetic metabolism in osteosarcoma

Osteosarcoma (OS) is a pediatric malignant bone tumor that predominantly affects adolescent and young adults. It has high risk for relapse and over the last four decades no improvement of prognosis was achieved. It is therefore crucial to identify new drug candidates for OS treatment to combat drug resistance, limit relapse, and stop metastatic spread. Two acquired hallmarks of cancer cells, mitochondria-related regulated cell death (RCD) and metabolism are intimately connected. Both have been shown to be dysregulated in OS, making them attractive targets for novel treatment. Promising OS treatment strategies focus on promoting RCD by targeting key molecular actors in metabolic reprogramming. The exact interplay in OS, however, has not been systematically analyzed. We therefore review these aspects by synthesizing current knowledge in apoptosis, ferroptosis, necroptosis, pyroptosis, and autophagy in OS. Additionally, we outline an overview of mitochondrial function and metabolic profiles in different preclinical OS models. Finally, we discuss the mechanism of action of two novel molecule combinations currently investigated in active clinical trials: metformin and the combination of ADI-PEG20, Docetaxel and Gemcitabine.

A ubiquitin-related gene signature for predicting prognosis and constructing molecular subtypes in osteosarcoma

Background: Ubiquitination is medicated by three classes of enzymes and has been proven to involve in multiple cancer biological processes. Moreover, dysregulation of ubiquitination has received a growing body of attention in osteosarcoma (OS) tumorigenesis and treatment. Therefore, our study aimed to identify a ubiquitin-related gene signature for predicting prognosis and immune landscape and constructing OS molecular subtypes.

Methods: Therapeutically Applicable Research to Generate Effective Treatments (TARGET) was regarded as the training set through univariate Cox regression, Lasso Cox regression, and multivariate Cox regression. The GSE21257 and GSE39055 served as the validation set to verify the predictive value of the signature. CIBERSORT was performed to show immune infiltration and the immune microenvironment. The NMF algorithm was used to construct OS molecular subtypes.

Results: In this study, we developed a ubiquitin-related gene signature including seven genes (UBE2L3, CORO6, DCAF8, DNAI1, FBXL5, UHRF2, and WDR53), and the gene signature had a good performance in predicting prognosis for OS patients (AUC values at 1/3/5 years were 0.957, 0.890, and 0.919). Multivariate Cox regression indicated that the risk score model and prognosis stage were also independent prognostic prediction factors. Moreover, analyses of immune cells and immune-related functions showed a significant difference in different risk score groups and the three clusters. The drug sensitivity suggested that IC50 of proteasome inhibitor (MG-132) showed a notable significance between the risk score groups (p < 0.05). Through the NMF algorithm, we obtained the three clusters, and cluster 3 showed better survival outcomes. The expression of ubiquitin-related genes (CORO6, UBE2L3, FBXL5, DNAI1, and DCAF8) showed an obvious significance in normal and osteosarcoma tissues.

Conclusion: We developed a novel ubiquitin-related gene signature which showed better predictive prognostic ability for OS and provided additional information on chemotherapy and immunotherapy. The OS molecular subtypes would also give a useful guide for individualized therapy.

A panel of necroptosis-related genes predicts the prognosis of pancreatic adenocarcinoma

•

The 5-NRGs signature can predict the prognosis of pancreatic adenocarcinoma.

•

The 5-NRGs signature reflects the immune landscape of pancreatic adenocarcinoma.

•

The 5-NRGs can be detected in exosomes of serum using RT-qPCR method.

Pancreatic adenocarcinoma (PAAD) has become one of the deadliest malignancies in the world. Since necroptosis plays a crucial role in regulating the immune system, it is necessary to develop novel prognostic biomarkers associated with necroptosis and explore its potential role in PAAD. The transcriptome RNA-seq data of PAAD were downloaded from the TCGA and GTEx databases. A prognostic signature was constructed by the least absolute shrinkage and selection operator (LASSO) Cox regression, and its prognostic value was evaluated by nomogram and validated in an independent GEO cohort. We identified a total of 24 differentially expressed NRGs in PAAD, and constructed a prognostic signature with 5 NRGs, which showed good performance in predicting the prognosis of PAAD patients. The ROC curves for 1-, 3-, and 5-year survival rate were 0.652, 0.778, and 0.817, respectively. This prognostic signature showed consistent prognosis prediction in an independent patient cohort. Furthermore, the correlations between 5-NRGs signature and TMB, MSI, histopathological classification, immune infiltration, immune types, and immunomodulators were all significant. Notably, the expression profiles of the five NRGs in exosomes of serum were consistent with their expression in tumor tissues. These data suggested that the 5-NRGs signature is a promising biomarker for predicting the prognosis of PAAD.

A novel signature to guide osteosarcoma prognosis and immune microenvironment: Cuproptosis-related lncRNA

Objective

Osteosarcoma (OS) is a common bone malignancy with poor prognosis. We aimed to investigate the relationship between cuproptosis-related lncRNAs (CRLncs) and the survival outcomes of patients with OS.

Methods

Transcriptome and clinical data of 86 patients with OS were downloaded from The Cancer Genome Atlas (TCGA). The GSE16088 dataset was downloaded from the Gene Expression Omnibus (GEO) database. The 10 cuproptosis-related genes (CRGs) were obtained from a recently published article on cuproptosis in Science. Combined analysis of OS transcriptome data and the GSE16088 dataset identified differentially expressed CRGs related to OS. Next, pathway enrichment analysis was performed. Co-expression analysis obtained CRLncs related to OS. Univariate COX regression analysis and least absolute shrinkage and selection operator (LASSO) regression analysis were used to construct the risk prognostic model of CRLncs. The samples were divided evenly into training and test groups to verify the accuracy of the model. Risk curve, survival, receiver operating characteristic (ROC) curve, and independent prognostic analyses were performed. Next, principal component analysis (PCA) and t-distributed stochastic neighbor embedding (t-SNE) analysis were performed. Single-sample gene set enrichment analysis (ssGSEA) was used to explore the correlation between the risk prognostic models and OS immune microenvironment. Drug sensitivity analysis identified drugs with potential efficacy in OS. Real-time quantitative PCR, Western blotting, and immunohistochemistry analyses verified the expression of CRGs in OS. Real-time quantitative PCR was used to verify the expression of CRLncs in OS.

Results

Six CRLncs that can guide OS prognosis and immune microenvironment were obtained, including three high-risk CRLncs (AL645608.6, AL591767.1, and UNC5B-AS1) and three low-risk CRLncs (CARD8-AS1, AC098487.1, and AC005041.3). Immune cells such as B cells, macrophages, T-helper type 2 (Th2) cells, regulatory T cells (Treg), and immune functions such as APC co-inhibition, checkpoint, and T-cell co-inhibition were significantly downregulated in high-risk groups. In addition, we obtained four drugs with potential efficacy for OS: AUY922, bortezomib, lenalidomide, and Z.LLNle.CHO. The expression of LIPT1, DLAT, and FDX1 at both mRNA and protein levels was significantly elevated in OS cell lines compared with normal osteoblast hFOB1.19. The mRNA expression level of AL591767.1 was decreased in OS, and that of AL645608.6, CARD8-AS1, AC005041.3, AC098487.1, and UNC5B-AS1 was upregulated in OS.

Conclusion

CRLncs that can guide OS prognosis and the immune microenvironment and drugs that may have a potential curative effect on OS obtained in this study provide a theoretical basis for OS survival research and clinical decision-making.

Expression of the Circadian Clock Gene ARNTL associated with DNA repair gene and prognosis of patient with osteosarcoma

PURPOSE

The study objects were to explore the correlation between the biological role of clock genes and clinical indicators in patients with osteosarcoma (OS).

METHODS

We acquired the clinical information and RNA sequencing data of OS samples from the TARGET database. The protein-protein interaction (PPI) network and expression correlation analysis of clock genes were performed. Then, the functional enrichment analysis of clock genes was analyzed. The survival analysis of clock genes in patients of OS was carried out by univariate cox regression, Kaplan-Meier (KM) curve and multivariate cox regression methods. Moreover, the spearmen correlation analysis was performed to explore the correlation between clock genes and DNA repair genes in patients with OS.

RESULTS

The PPI network and expression correlation analysis of clock genes indicated that the clock genes were highly correlated with each other. The survival analysis of clock genes found that clock gene ARNTL is a protective factor for the prognosis of patients with OS. We found that ARNTL was positively related to DNA repair genes and was involved in the biological process of DNA damage repair in patients with OS.

CONCLUSIONS

ARNTL may affect the prognosis and chemotherapy response of patients with OS by regulating DNA repair pathways.

Comprehensive Analysis of a Ferroptosis-Related lncRNA Signature for Predicting Prognosis and Immune Landscape in Osteosarcoma

Research on the implications of ferroptosis in tumors has increased rapidly in the last decades. There are evidences that ferroptosis is involved in several aspects of cancer biology, including tumor progression, metastasis, immunomodulation, and therapeutic response. Nonetheless, the interaction between ferroptosis-related lncRNAs (FRLs) and the osteosarcoma immune microenvironment is poorly understood. In this study, a risk model composed of FRLs was developed using univariate and LASSO Cox regression analyses. On the basis of this model, FRL scores were calculated to systematically explore the role of the model in predicting the prognosis and immune characteristics of osteosarcoma patients. Survival analysis showed that osteosarcoma samples with lower FRL-score had better overall survival. After predicting the abundance of immune cells in osteosarcoma microenvironment by single-sample gene-set enrichment analysis (ssGSEA) and ESTIMATE analysis, we found that the FRL-score could distinguish immune function, immune score, stromal score, tumor purity, and tumor infiltration of immune cells in different osteosarcoma patients. In addition, FRL-score was also associated with immune checkpoint gene expression and half-maximal inhibitory concentration of chemotherapeutic agents. Finally, we confirmed that knockdown of RPARP-AS1 suppressed the malignant activity of osteosarcoma cells in vitro experiments. In general, the FRL-based prognostic signature could promote our understanding of the immune microenvironment characteristics of osteosarcoma and guide more effective treatment regimens.

Identification of a Solute Carrier Family-Based Signature for Predicting Overall Survival in Osteosarcoma

Given the important role of SLC family in essential physiological processes including nutrient uptake, ion transport, and waste removal, and that their dysregulation was found in distinct forms of cancer, here we identified a novel gene signature of SLC family for patient risk stratification in osteosarcoma. Gene expression data and relevant clinical materials of osteosarcoma samples were retrieved from The Cancer Genome Atlas (TCGA) database. Prognosis-related SLC genes were identified by performing univariate Cox regression analysis and were utilized to construct a four-SLC gene signature in osteosarcoma. It allowed patients to be classified into high- and low-risk groups, and Kaplan-Meier survival analysis in the training, testing, entire, and external GSE21257 cohorts suggested that the overall survival of patients in high-risk group was consistently worse than that in low-risk group, suggesting the promising accuracy and generalizability of the SLC-based signature in predicting the prognosis of patients with osteosarcoma. Moreover, univariate and multivariate Cox regression analyses indicated that the derived risk score was the only independent prognostic factor for osteosarcoma patients in TCGA and GSE21257 cohorts. Besides, a prognostic nomogram comprising the derived risk score and clinical features including gender and age was developed for clinical decision-making. Functional enrichment analyses of the differentially expressed genes between high- and low-risk group revealed that immune-related biological processes and pathways were significantly enriched. Estimation of tumor immune microenvironment using ESTIMATE algorithm revealed that patients with lower risk score had higher stromal, immune, and ESTIMATE score, and lower tumor purity. ssGSEA analyses indicated that the scores of various immune subpopulations including CD8+ T cells, DCs, and TIL were lower in high-risk group than these in low-risk group in both cohorts. As for the related immune functions, the scores of APC co-inhibition, CCR, check-point, T cell co-stimulation, and Type II IFN response were lower in high-risk group than these in low-risk group in both cohorts. In all, we identified a novel prognostic signature based on four SLC family genes that accurately predicted overall survival in osteosarcoma patients. Furthermore, the signature is linked to differences in immunological status and immune cell infiltrations in the tumor microenvironment.

Integrated analyses of an RNA binding protein-based signature related to tumor immune microenvironment and candidate drugs in osteosarcoma

OBJECTIVE

Osteosarcoma is the most frequent primary bone malignancy, associated with frequent recurrence and lung metastasis. RNA-binding proteins (RBPs) are pivotal in regulating several aspects of cancer biology. Nonetheless, interaction between RBPs and the osteosarcoma immune microenvironment is poorly understood. We investigated whether RBPs can predict prognosis and immunotherapy response in osteosarcoma patients.

METHODS

We constructed an RBP-related prognostic signature (RRPS) by univariate coupled with multivariate analyses and verified the independent prognostic efficacy of the signature. Single-sample Gene Set Enrichment Analysis (ssGSEA) along with ESTIMATE analysis were carried out to investigate the variations in immune characteristics between subgroups with various RRPS-scores. Furthermore, we investigatedpossible small molecule drugs using the connectivity map database and validated the expression of hub RBPs by qRT-PCR.

RESULTS

The RRPS, consisting of seven hub RBPs, was an independent prognostic factor compared to traditional clinical features. The RRPS could distinguish immune functions, immune score, stromal score, tumor purity and tumor infiltration by immune cells in different osteosarcoma subjects. Additionally, patients with high RRPS-scores had lower expression of immune checkpoint genes than patients with low RRPS-scores. We finally identified six small molecule drugs that may improve prognosis in osteosarcoma patients and substantiated notable differences in the contents of these RBPs.

CONCLUSION

We evaluated the prognostic value and clinical application of an RBPs-based prognostic signature and identified promising biomarkers to predict immune cell infiltration and immunotherapy response in osteosarcoma.

Comprehensive Analysis of Pyroptosis-Related Genes and Tumor Microenvironment Infiltration Characterization in Papillary Renal Cell Carcinoma

Background: Immunotherapy has emerged as an important technique for treating a variety of cancers. The dynamic interplay between tumor cells and invading lymphocytes in the tumor microenvironment is responsible for the good response to immunotherapy (TME). Pyroptosis, or inflammation-induced cell death, is closely linked to a number of cancers. However, in papillary renal cell carcinoma (KIRP), the association between pyroptosis and clinical prognosis, immune cell infiltration, and immunotherapy impact remains unknown.

Methods: We carefully investigated the link between pyroptosis and tumor growth, prognosis, and immune cell infiltration by evaluating 52 pyroptosis-related genes. The PRG score was utilized to measure a single tumor patient’s pyroptosis pattern. After that, we looked at how well these values predicted prognoses and therapy responses in KIRP.

Results: We discovered that PRG differences between subgroups were linked to clinical and pathological aspects, prognosis, and TME in two separate genetic subtypes. After that, a PRG score for estimating overall survival (OS) was developed, and its predictive potential in KIRP patients was confirmed. As a result, we developed a very precise nomogram to improve the PRG score’s clinical usefulness. A low PRG score, which is determined by mutation load and immune activation, suggests a good chance of survival. Furthermore, the PRG score was linked to chemotherapeutic drug sensitivity in a substantial way.

Conclusions: The possible functions of PRGs in the TME, clinical and pathological characteristics, and prognosis were established in our thorough investigation of PRGs in KIRP. These results might help us better understand PRGs in KIRP and offer a new avenue for prognostic evaluation and the development of more effective immunotherapy treatments.