The Immune Landscape of Osteosarcoma: Implications for Prognosis and Treatment Response

Proteomic and transcriptomic analyses identify apo-transcobalamin-II as a biomarker of overall survival in osteosarcoma

Background

The large-scale proteomic platform known as the SomaScan® assay is capable of simultaneously measuring thousands of proteins in patient specimens through next-generation aptamer-based multiplexed technology. While previous studies have utilized patient peripheral blood to suggest serum biomarkers of prognostic or diagnostic value in osteosarcoma (OSA), the most common primary pediatric bone cancer, they have ultimately been limited in the robustness of their analyses. We propose utilizing this aptamer-based technology to describe the systemic proteomic milieu in patients diagnosed with this disease.

Methods

To determine novel biomarkers associated with overall survival in OSA, we deployed the SomaLogic SomaScan® 7k assay to investigate the plasma proteomic profile of naive primary, recurrent, and metastatic OSA patients. Following identification of differentially expressed proteins (DEPs) between 2-year deceased and survivor cohorts, publicly available databases including Survival Genie, TIGER, and KM Plotter Immunotherapy, among others, were utilized to investigate the significance of our proteomic findings.

Results

Apo-transcobalamin-II (APO-TCN2) was identified as the most DEP between 2-year deceased and survivor cohorts (Log2 fold change = 6.8, P-value = 0.0017). Survival analysis using the Survival Genie web-based platform indicated that increased intratumoral TCN2 expression was associated with better overall survival in both OSA (TARGET-OS) and sarcoma (TCGA-SARC) datasets. Cell-cell communication analysis using the TIGER database suggested that TCN2+ Myeloid cells likely interact with marginal zone and immunoglobin-producing B lymphocytes expressing the TCN2 receptor (CD320) to promote their proliferation and survival in both non-small cell lung cancer and melanoma tumors. Analysis of publicly available OSA scRNA-sequencing datasets identified similar populations in naive primary tumors. Furthermore, circulating APO-TCN2 levels in OSA were then associated with a plasma proteomic profile likely necessary for robust B lymphocyte proliferation, infiltration, and formation of intratumoral tertiary lymphoid structures for improved anti-tumor immunity.

Conclusions

Overall, APO-TCN2, a circulatory protein previously described in various lymphoproliferative disorders, was associated with 2-year survival status in patients diagnosed with OSA. The relevance of this protein and apparent immunological function (anti-tumor B lymphocyte responses) was suggested using publicly available solid tumor RNA-sequencing datasets. Further studies characterizing the biological function of APO-TCN2 and its relevance in these diseases is warranted.

Global research trends and hotspots in metabolomics of osteosarcoma: a decade-spanning bibliometric and visualized analysis

Object

Osteosarcoma is a malignant tumor originating from the bones, commonly found in children and adolescents, especially in rapidly growing bone areas such as the knees and upper arms. In this study, we aim to delineate the evolution and convergence of research themes in osteosarcoma metabolomics over the past decade, identify major contributors, and forecast emerging trends that could direct future research efforts.

Method

The bibliometric method has been applied to systematically analyze the literature in the field of osteosarcoma metabolomics. The relevant literatures were collected from the Web of Science Core Collection, spanning from January 1, 2014, to December 31, 2023. Tools such as CiteSpace, Bibliometrix, and VOSviewer were used for the visual analysis of the collected literatures. The focused information includes institutions, journals, countries, authors, keywords, and citations.

Result

Various aspects in the field of osteosarcoma metabolism were analyzed. Shanghai Jiao Tong University has published the most papers in the past ten years, followed by Central South University and Zhejiang University. Among the sources, the international journal of molecular sciences publishes the most articles, and oncotarget is the journal with the highest H index. According to Bradford's law, there are 34 core journals identified. A total of 5501 authors participated in the creation of papers in this field. The distribution of authors follows Lotka`s Law, and 85.3% of authors have only one article. 46% of the corresponding authors are from China, but most of these corresponding authors are not good at international cooperation. China also has the largest number of publications, followed by the United States. It can be confirmed that China dominates this field. Among the keywords, "expression" is the keyword that has received the most attention in the past ten years. All keywords can be divided into 9 clusters. Based on the explosive words and hot topics each year, we speculate that future research will focus on the tumor microenvironment, molecular mechanisms and autophagy, targeted therapies and inhibitors.

Conclusion

In summary, this study comprehensively analyzed the current state of research in the field of osteosarcoma metabolism through bibliometric methods. The findings revealed the development trends and research hotspots in this field, which may provide valuable references for future research directions.

Spatial multiplexed immunofluorescence analysis reveals coordinated cellular networks associated with overall survival in metastatic osteosarcoma

Patients diagnosed with advanced osteosarcoma, often in the form of lung metastases, have abysmal five-year overall survival rates. The complexity of the osteosarcoma immune tumor microenvironment has been implicated in clinical trial failures of various immunotherapies. The purpose of this exploratory study was to spatially characterize the immune tumor microenvironment of metastatic osteosarcoma lung specimens. Knowledge of the coordinating cellular networks within these tissues could then lead to improved outcomes when utilizing immunotherapy for treatment of this disease. Importantly, various cell types, interactions, and cellular neighborhoods were associated with five-year survival status. Of note, increases in cellular interactions between T lymphocytes, positive for programmed cell death protein 1, and myeloid-derived suppressor cells were observed in the 5-year deceased cohort. Additionally, cellular neighborhood analysis identified an Immune-Cold Parenchyma cellular neighborhood, also associated with worse 5-year survival. Finally, the Osteosarcoma Spatial Score, which approximates effector immune activity in the immune tumor microenvironment through the spatial proximity of immune and tumor cells, was increased within 5-year survivors, suggesting improved effector signaling in this patient cohort. Ultimately, these data represent a robust spatial multiplexed immunofluorescence analysis of the metastatic osteosarcoma immune tumor microenvironment. Various communication networks, and their association with survival, were described. In the future, identification of these networks may suggest the use of specific, combinatory immunotherapeutic strategies for improved anti-tumor immune responses and outcomes in osteosarcoma.

Oral microbiota distinguishes patients with osteosarcoma from healthy controls

Objective

The human microbiota plays a key role in cancer diagnosis, pathogenesis, and treatment. However, osteosarcoma-associated oral microbiota alterations have not yet been unraveled. The aim of this study was to explore the characteristics of oral microbiota in osteosarcoma patients compared to healthy controls, and to identify potential microbiota as a diagnostic tool for osteosarcoma.

Methods

The oral microbiota was analyzed in osteosarcoma patients (n = 45) and matched healthy controls (n = 90) using 16S rRNA MiSeq sequencing technology.

Results

The microbial richness and diversity of the tongue coat were increased in osteosarcoma patients as estimated by the abundance-based coverage estimator indices, the Chao, and observed operational taxonomy units (OTUs). Principal component analysis delineated that the oral microbial community was significant differences between osteosarcoma patients and healthy controls. 14 genera including Rothia, Halomonas, Rhodococcus, and Granulicatella were remarkably reduced, whereas Alloprevotella, Prevotella, Selenomonas, and Campylobacter were enriched in osteosarcoma. Eventually, the optimal four OTUs were identified to construct a microbial classifier by the random forest model via a fivefold cross-validation, which achieved an area under the curve of 99.44% in the training group (30 osteosarcoma patients versus 60 healthy controls) and 87.33% in the test group (15 osteosarcoma patients versus 30 healthy controls), respectively. Notably, oral microbial markers validated strong diagnostic potential distinguishing osteosarcoma patients from healthy controls.

Conclusion

This study comprehensively characterizes the oral microbiota in osteosarcoma and reveals the potential efficacy of oral microbiota-targeted biomarkers as a noninvasive biological diagnostic tool for osteosarcoma.

Clinical and translational implications of immunotherapy in sarcomas

Immunotherapy has emerged as promising treatment in sarcomas, but the high variability in terms of histology, clinical behavior and response to treatments determines a particular challenge for its role in these neoplasms. Tumor immune microenvironment (TiME) of sarcomas reflects the heterogeneity of these tumors originating from mesenchymal cells and encompassing more than 100 histologies. Advances in the understanding of the complexity of TiME have led to an improvement of the immunotherapeutic responsiveness in sarcomas, that at first showed disappointing results. The proposed immune-classification of sarcomas based on the interaction between immune cell populations and tumor cells showed to have a prognostic and potential predictive role for immunotherapies. Several studies have explored the clinical impact of immune therapies in the management of these histotypes leading to controversial results. The presence of Tumor Infiltrating Lymphocytes (TIL) seems to correlate with an improvement in the survival of patients and with a higher responsiveness to immunotherapy. In this context, it is important to consider that also immune-related genes (IRGs) have been demonstrated to have a key role in tumorigenesis and in the building of tumor immune microenvironment. The IRGs landscape in soft tissue and bone sarcomas is characterized by the connection between several tumor-related genes that can assume a potential prognostic and predictive therapeutic role. In this paper, we reviewed the state of art of the principal immune strategies in the management of sarcomas including their clinical and translational relevance.

Fibroblast Subpopulations in Systemic Sclerosis: Functional Implications of Individual Subpopulations and Correlations with Clinical Features

Fibroblasts constitute a heterogeneous population of cells. In this study, we integrated single-cell RNA-sequencing and bulk RNA-sequencing data as well as clinical information to study the role of individual fibroblast populations in systemic sclerosis (SSc). SSc skin demonstrated an increased abundance of COMP+, COL11A1+, MYOC+, CCL19+, SFRP4/SFRP2+, and PRSS23/SFRP2+ fibroblasts signatures and decreased proportions of CXCL12+ and PI16+ fibroblast signatures in the Prospective Registry of Early Systemic Sclerosis and Genetics versus Environment in Scleroderma Outcome Study cohorts. Numerical differences were confirmed by multicolor immunofluorescence for selected fibroblast populations. COMP+, COL11A1+, SFRP4/SFRP2+, PRSS23/SFRP2+, and PI16+ fibroblasts were similarly altered between normal wound healing and patients with SSc. The proportions of profibrotic COMP+, COL11A1+, SFRP4/SFRP2+, and PRSS23/SFRP2+ and proinflammatory CCL19+ fibroblast signatures were positively correlated with clinical and histopathological parameters of skin fibrosis, whereas signatures of CXCL12+ and PI16+ fibroblasts were inversely correlated. Incorporating the proportions of COMP+, COL11A1+, SFRP4/SFRP2+, and PRSS23/SFRP2+ fibroblast signatures into machine learning models improved the classification of patients with SSc into those with progressive versus stable skin fibrosis. In summary, the profound imbalance of fibroblast subpopulations in SSc may drive the progression of skin fibrosis. Specific targeting of disease-relevant fibroblast populations may offer opportunities for the treatment of SSc and other fibrotic diseases.

The immune cell infiltration-associated molecular subtypes and gene signature predict prognosis for osteosarcoma patients

Host immune dysregulation involves in the initiation and development of osteosarcoma (OS). However, the exact role of immune cells in OS remains unknown. We aimed to distinguish the molecular subtypes and establish a prognostic model in OS patients based on immunocyte infiltration. The gene expression profile and corresponding clinical feature of OS patients were obtained from TARGET and GSE21257 datasets. MCP-counter and univariate Cox regression analyses were applied to identify immune cell infiltration-related molecular subgroups. Functional enrichment analysis and immunocyte infiltration analysis were performed between two subgroups. Furthermore, Cox regression and LASSO analyses were performed to establish the prognostic model for the prediction of prognosis and metastasis in OS patients. The subgroup with low infiltration of monocytic lineage (ML) was related to bad prognosis in OS patients. 435 DEGs were screened between the two subgroups. Functional enrichment analysis revealed these DEGs were involved in immune- and inflammation-related pathways. Three important genes (including TERT, CCDC26, and IL2RA) were identified to establish the prognostic model. The risk model had good prognostic performance for the prediction of metastasis and overall survival in OS patients. A novel stratification system was established based on ML-related signature. The risk model could predict the metastasis and prognosis in OS patients. Our findings offered a novel sight for the prognosis and development of OS.

Immunotherapy Innovations in the Fight against Osteosarcoma: Emerging Strategies and Promising Progress

Immunosuppressive elements within the tumor microenvironment are the primary drivers of tumorigenesis and malignant advancement. The presence, as well as the crosstalk between myeloid-derived suppressor cells (MDSCs), osteosarcoma-associated macrophages (OS-Ms), regulatory T cells (Tregs), and endothelial cells (ECs) with osteosarcoma cells cause the poor prognosis of OS. In addition, the consequent immunosuppressive factors favor the loss of treatment potential. Nanoparticles offer a means to dynamically and locally manipulate immuno-nanoparticles, which present a promising strategy for transforming OS-TME. Additionally, chimeric antigen receptor (CAR) technology is effective in combating OS. This review summarizes the essential mechanisms of immunosuppressive cells in the OS-TME and the current immune-associated strategies. The last part highlights the limitations of existing therapies and offers insights into future research directions.

Exploration of anti‑osteosarcoma activity of asiatic acid based on network pharmacology and in vitro experiments

Osteosarcomas are malignant bone tumors that typically originate in the epiphyses of the long bones of the extremities in adolescents. Asiatic acid has been reported to possess anti‑inflammatory, neuroprotective, antidiabetic, antitumor and antimicrobial activities. The present study used a combination of network pharmacological prediction and in vitro experimental validation to explore the potential pharmacological mechanism of asiatic acid against osteosarcoma. A total of 78 potential asiatic acid targets in osteosarcoma were identified using databases. Kyoto Encyclopedia of Genes and Genomes analysis indicated that the PI3K/AKT and MAPK signaling pathways are essential in the treatment of osteosarcoma with asiatic acid. Molecular docking revealed binding of asiatic acid to EGFR, Caspase‑3, ESR1, HSP90AA1, IL‑6 and SRC proteins. asiatic acid inhibited proliferation through G2/M cell cycle arrest in osteosarcoma cells. In addition, asiatic acid induced mitochondria‑dependent apoptosis as demonstrated by increases in Bax and VDAC1 expression, and a decrease in Bcl‑2 protein expression. The increased autophagosomes, increased LC3‑II/I ratios and decreased p62 expression in the treatment group indicated that asiatic acid triggered autophagy. In addition, asiatic acid decreased the levels of phosphorylated (p‑)PI3K/PI3K and p‑AKT/AKT, increased reactive oxygen species (ROS) and upregulated the levels of p‑ERK1/2/ERK1/2, p‑p38/p38 and p‑JNK/JNK in osteosarcoma cells. These results demonstrated that asiatic acid inhibited osteosarcoma cells proliferation by inhibiting PI3K/AKT and activating ROS/MAPK signaling pathways, suggesting asiatic acid is a potential agent against osteosarcoma.

Unveiling the Protective Role of Melatonin in Osteosarcoma: Current Knowledge and Limitations

Melatonin, an endogenous neurohormone produced by the pineal gland, has received increased interest due to its potential anti-cancer properties. Apart from its well-known role in the sleep-wake cycle, extensive scientific evidence has shown its role in various physiological and pathological processes, such as inflammation. Additionally, melatonin has demonstrated promising potential as an anti-cancer agent as its function includes inhibition of tumorigenesis, induction of apoptosis, and regulation of anti-tumor immune response. Although a precise pathophysiological mechanism is yet to be established, several pathways related to the regulation of cell cycle progression, DNA repair mechanisms, and antioxidant activity have been implicated in the anti-neoplastic potential of melatonin. In the current manuscript, we focus on the potential anti-cancer properties of melatonin and its use in treating and managing pediatric osteosarcoma. This aggressive bone tumor primarily affects children and adolescents and is treated mainly by surgical and radio-oncological interventions, which has improved survival rates among affected individuals. Significant disadvantages to these interventions include disease recurrence, therapy-related toxicity, and severe/debilitating side effects that the patients have to endure, significantly affecting their quality of life. Melatonin has therapeutic effects when used for treating osteosarcoma, attributed to its ability to halt cancer cell proliferation and trigger apoptotic cell death, thereby enhancing chemotherapeutic efficacy. Furthermore, the antioxidative function of melatonin alleviates harmful side effects of chemotherapy-induced oxidative damage, aiding in decreasing therapeutic toxicities. The review concisely explains the many mechanisms by which melatonin targets osteosarcoma, as evidenced by significant results from several in vitro and animal models. Nevertheless, if further explored, human trials remain a challenge that could shed light and support its utility as an adjunctive therapeutic modality for treating osteosarcoma.

C-C Motif Chemokine Ligand 5 (CCL5): A Potential Biomarker and Immunotherapy Target for Osteosarcoma

BACKGROUND

Osteosarcoma (OS) is the most common primary malignant tumor of bone tissue, which has an insidious onset and is difficult to detect early, and few early diagnostic markers with high specificity and sensitivity. Therefore, this study aims to identify potential biomarkers that can help diagnose OS in its early stages and improve the prognosis of patients.

METHODS

The data sets of GSE12789, GSE28424, GSE33382 and GSE36001 were combined and normalized to identify Differentially Expressed Genes (DEGs). The data were analyzed by Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genome (KEGG) and Disease Ontology (DO). The hub gene was selected based on the common DEG that was obtained by applying two regression methods: the Least Absolute Shrinkage and Selection Operator (LASSO) and Support vVector Machine (SVM). Then the diagnostic value of the hub gene was evaluated in the GSE42572 data set. Finally, the correlation between immunocyte infiltration and key genes was analyzed by CIBERSORT.

RESULTS

The regression analysis results of LASSO and SVM are the following three DEGs: FK501 binding protein 51 (FKBP5), C-C motif chemokine ligand 5 (CCL5), complement component 1 Q subcomponent B chain (C1QB). We evaluated the diagnostic performance of three biomarkers (FKBP5, CCL5 and C1QB) for osteosarcoma using receiver operating characteristic (ROC) analysis. In the training group, the area under the curve (AUC) of FKBP5, CCL5 and C1QB was 0.907, 0.874 and 0.676, respectively. In the validation group, the AUC of FKBP5, CCL5 and C1QB was 0.618, 0.932 and 0.895, respectively. It is noteworthy that these genes were more expressed in tumor tissues than in normal tissues by various immune cell types, such as plasma cells, CD8+ T cells, T regulatory cells (Tregs), activated NK cells, activated dendritic cells and activated mast cells. These immune cell types are also associated with the expression levels of the three diagnostic genes that we identified.

CONCLUSION

We found that CCL5 can be considered an early diagnostic gene of osteosarcoma, and CCL5 interacts with immune cells to influence tumor occurrence and development. These findings have important implications for the early detection of osteosarcoma and the identification of novel therapeutic targets.

Bone-targeting exosome nanoparticles activate Keap1 / Nrf2 / GPX4 signaling pathway to induce ferroptosis in osteosarcoma cells

BACKGROUND

In recent years, the development of BMSCs-derived exosomes (EXO) for the treatment of osteosarcoma (OS) is a safe and promising modality for OS treatment, which can effectively deliver drugs to tumor cells in vivo. However, the differences in the drugs carried, and the binding of EXOs to other organs limit their therapeutic efficacy. Therefore, improving the OS-targeting ability of BMSCs EXOs and developing new drugs is crucial for the clinical application of targeted therapy for OS.

RESULTS

In this study, we constructed a potential therapeutic nano platform by modifying BMSCs EXOs using the bone-targeting peptide SDSSD and encapsulated capreomycin (CAP) within a shell. These constructed nanoparticles (NPs) showed the ability of homologous targeting and bone-targeting exosomes (BT-EXO) significantly promotes cellular endocytosis in vitro and tumor accumulation in vivo. Furthermore, our results revealed that the constructed NPs induced ferroptosis in OS cells by prompting excessive accumulation of reactive oxygen species (ROS), Fe2+ aggregation, and lipid peroxidation and further identified the potential anticancer molecular mechanism of ferroptosis as transduced by the Keap1/Nrf2/GPX4 signaling pathway. Also, these constructed NP-directed ferroptosis showed significant inhibition of tumor growth in vivo with no significant side effects.

CONCLUSION

These results suggest that these constructed NPs have superior anticancer activity in mouse models of OS in vitro and in vivo, providing a new and promising strategy for combining ferroptosis-based chemotherapy with targeted therapy for OS.

A Disulfidptosis-Related Gene Signature Associated with Prognosis and Immune Cell Infiltration in Osteosarcoma

Osteosarcoma (OS) stands as a leading aggressive bone malignancy that primarily affects children and adolescents worldwide. A recently identified form of programmed cell death, termed Disulfidptosis, may have implications for cancer progression. Yet, its role in OS remains elusive. To elucidate this, we undertook a thorough examination of Disulfidptosis-related genes (DRGs) within OS. This involved parsing expression data, clinical attributes, and survival metrics from the TARGET and GEO databases. Our analysis unveiled a pronounced association between the expression of specific DRGs, particularly MYH9 and LRPPRC, and OS outcome. Subsequent to this, we crafted a risk model and a nomogram, both honed for precise prognostication of OS prognosis. Intriguingly, risks associated with DRGs strongly resonated with immune cell infiltration levels, myriad immune checkpoints, genes tethered to immunotherapy, and sensitivities to systematic treatments. To conclude, our study posits that DRGs, especially MYH9 and LRPPRC, hold potential as pivotal architects of the tumor immune milieu in OS. Moreover, they may offer predictive insights into treatment responses and serve as reliable prognostic markers for those diagnosed with OS.

Gossypetin Is a Novel Modulator of Inflammatory Cytokine Production and a Suppressor of Osteosarcoma Cell Growth

Osteosarcoma (OS) is a common childhood sarcoma, and its treatment is hindered by adverse effects, chemoresistance, and recurrence. Interleukin (IL)-6 production by tumors plays a significant role in inflammation, carcinogenesis, and metastasis. This study aimed to investigate the antiproliferative potential of luteolin derivatives in OS and to evaluate interleukin production. MG-63, Saos-2, HOS, and 143B human OS cell lines were incubated with luteolin and eight derivatives containing hydroxy, chlorine, or alkyl substitutions. The cell viability and growth were evaluated in the presence of these compounds. Apoptosis was also examined through the analysis of the Bax expression and caspase-3 activity. Finally, the gossypetin effects were measured regarding the production of proinflammatory cytokines interleukin (IL)-6, IL-1β, and IL-12p70. Our findings show that gossypetin was the most potent compound, with proliferation-suppressing activities that induced a series of critical events, including the inhibition of the cell viability and growth. Apoptosis was associated with enhanced caspase-3 activity and increased Bax expression, indicating the involvement of the intrinsic pathway of apoptosis. Moreover, pre-/co-treatment with gossypetin significantly reduced the autocrine production of proinflammatory cytokines. Further investigation is required; nevertheless, considering the link between inflammation, carcinogenesis, and metastasis in OS, our findings suggest that gossypetin exhibits anti-proliferative and anti-inflammatory properties that are potentially relevant in the clinical context.

HOXC13-driven TIMM13 overexpression promotes osteosarcoma cell growth

TIMM13 (translocase of inner mitochondrial membrane 13) located at the mitochondrial intermembrane space is vital for the integrity and function of mitochondria. We found that the mitochondrial protein TIMM13 is upregulated in human OS tissues and cells. In patient-derived primary OS cells and established cell lines, TIMM13 shRNA or knockout provoked mitochondrial dysfunction, causing mitochondrial depolarization, reactive oxygen species production, and oxidative injury, as well as lipid peroxidation, DNA damage, and ATP depletion. Moreover, TIMM13 depletion provoked OS cell apoptosis and inhibited cell proliferation and migration. Conversely, ectopic TIMM13 overexpression increased ATP contents, enhancing OS cell proliferation and migration. Moreover, we discovered that Akt-mTOR activation was inhibited with TIMM13 depletion in primary OS cells. Further studies revealed that HOXC13 (Homeobox C13)-dependent TIMM13 transcription was significantly increased in OS tissues and cells. Whereas TIMM13 transcription and expression were decreased following HOXC13 silencing in primary OS cells. In vivo, TIMM13 KO potently inhibited OS xenograft growth in the proximal tibia of nude mice. TIMM13 KO also induced Akt-mTOR inactivation, ATP depletion, oxidative injury, and apoptosis in the in situ OS tumors. Together, upregulation of the mitochondrial protein TIMM13 is important for OS cell growth, representing a novel and promising therapeutic target.

Rewiring innate and adaptive immunity with TLR9 agonist to treat osteosarcoma

BACKGROUND

Osteosarcoma (OS) is the most common primary bone tumor in children and adolescent. Surgery and multidrug chemotherapy are the standard of treatment achieving 60-70% of event-free survival for localized disease at diagnosis. However, for metastatic disease, the prognosis is dismal. Exploiting immune system activation in the setting of such unfavorable mesenchymal tumors represents a new therapeutic challenge.

METHODS

In immune competent OS mouse models bearing two contralateral lesions, we tested the efficacy of intralesional administration of a TLR9 agonist against the treated and not treated contralateral lesion evaluating abscopal effect. Multiparametric flow cytometry was used to evaluate changes of the tumor immune microenviroment. Experiments in immune-deficient mice allowed the investigation of the role of adaptive T cells in TLR9 agonist effects, while T cell receptor sequencing was used to assess the expansion of specific T cell clones.

RESULTS

TLR9 agonist strongly impaired the growth of locally-treated tumors and its therapeutic effect also extended to the contralateral, untreated lesion. Multiparametric flow cytometry showed conspicuous changes in the immune landscape of the OS immune microenvironment upon TLR9 engagement, involving a reduction in M2-like macrophages, paralleled by increased infiltration of dendritic cells and activated CD8 T cells in both lesions. Remarkably, CD8 T cells were needed for the induction of the abscopal effect, whereas they were not strictly necessary for halting the growth of the treated lesion. T cell receptor (TCR) sequencing of tumor infiltrating CD8 T cells showed the expansion of specific TCR clones in the treated tumors and, remarkably, their selected representation in the contralateral untreated lesions, providing the first evidence of the rewiring of tumor-associated T cell clonal architectures.

CONCLUSIONS

Overall these data indicate that the TLR9 agonist acts as an in situ anti-tumor vaccine, activating an innate immune response sufficient to suppress local tumor growth while inducing a systemic adaptive immunity with selective expansion of CD8 T cell clones, which are needed for the abscopal effect.

Signal Pathways and microRNAs in Osteosarcoma Growth and the Dual Role of Mesenchymal Stem Cells in Oncogenesis

The major challenges in Osteosarcoma (OS) therapy are its heterogeneity and drug resistance. The development of new therapeutic approaches to overcome the major growth mechanisms of OS is urgently needed. The search for specific molecular targets and promising innovative approaches in OS therapy, including drug delivery methods, is an urgent problem. Modern regenerative medicine focuses on harnessing the potential of mesenchymal stem cells (MSCs) because they have low immunogenicity. MSCs are important cells that have received considerable attention in cancer research. Currently, new cell-based methods for using MSCs in medicine are being actively investigated and tested, especially as carriers for chemotherapeutics, nanoparticles, and photosensitizers. However, despite the inexhaustible regenerative potential and known anticancer properties of MSCs, they may trigger the development and progression of bone tumors. A better understanding of the complex cellular and molecular mechanisms of OS pathogenesis is essential to identify novel molecular effectors involved in oncogenesis. The current review focuses on signaling pathways and miRNAs involved in the development of OS and describes the role of MSCs in oncogenesis and their potential for antitumor cell-based therapy.

Single-cell RNA-seq reveals intratumoral heterogeneity in osteosarcoma patients: A review

While primary bone malignancies make up just 0.2% of all cancers, osteosarcoma (OS) is the third most common cancer in adolescents. Due to its highly complex and heterogeneous tumor microenvironment (TME), OS has proven difficult to treat. There has been little to no improvement in therapy for this disease over the last 40 years. Even the recent success of immunotherapies in other blood-borne and solid malignancies has not translated to OS. With frequent recurrence and lung metastases continuing to pose a challenge in the clinic, recent advancements in molecular profiling, such as single-cell RNA sequencing (scRNA-seq), have proven useful in identifying novel biomarkers of OS tumors while providing new insight into this TME that could potentially lead to new therapeutic options. This review combines the analyses of over 150,000 cells from 18 lesions ranging from primary, recurrent, and metastatic OS lesions, revealing distinct cellular populations and gene signatures that exist between them. Here, we detail these previous findings and ultimately convey the intratumoral heterogeneity that exists within OS tumor specimens.

Identification of RALA as a Therapeutic Target and Prognostic Predictor of Osteosarcoma

Background

Osteosarcoma (OS) is the most common primary aggressive sarcoma of bone, with massive aberrant expression of oncogenes related to the development of OS. RALA, a kind of small Ras-like guanosine triphosphatases, has been identified as a potential therapeutic target in several types of tumor, but its role in OS remains largely unknown.

Methods

Abnormal expression of RALA was proven in the Cancer Genome Atlas (TCGA), Genotype-Tissue Expression (GTEx), Therapeutically Applicable Research to Generate Effective Treatments (TARGET), and RNA-sequence of samples and cell lines. The role of RALA in OS was analyzed in terms of DNA methylation, immune cell infiltration, and patient survival. The cancer-promoting effect of RALA was demonstrated in cell lines and xenograft osteosarcoma models. A prognostic scoring model incorporating RALA as an indicator was established with the clinical samples that we collected.

Results

The results showed that RALA was highly expressed in human OS tissues and cell lines. Survival analysis demonstrated that RALA was the sole independent risk factor for poor overall survival and disease-free survival in OS patients and impacted the proportion of infiltrating immune cells and DNA methylation in the OS tumor microenvironment. By gene-gene interaction analysis, we found that the expression of RALA was highly correlated to the expression of ABCE1. Similar to RALA, upregulated ABCE1 is correlated with poor survival outcome of OS patients. In addition, the functional experiment demonstrated that higher expression of RALA promoted the proliferation, migration, and invasion of OS cells. In vivo results were similar with the in vitro results. We examined m6a methylation-related genes and found that m6A methylation is responsible for the abnormal expression of RALA. Finally, the prognostic prediction model of RALA could be used to predict the long-term outcome of OS patients.

Conclusions

We identified RALA as an oncogene in OS, and RALA upregulation in a concerted manner with ABCE1 was significantly associated with worse outcomes of OS patients. Targeting RALA may prove to be a novel target for OS immunotherapy in future clinical practice.

Development of personalized classifier based on metastasis and the immune microenvironment to predict the prognosis of osteosarcoma patients

Background

Osteosarcoma is a common malignant bone tumor with a poor prognosis. The progression and metastasis of osteosarcoma are significantly influenced by the tumor microenvironment (TME). This study aimed to develop a personalized classifier based on metastasis and immune cells in the TME to achieve better prognostic prediction in osteosarcoma.

Methods

Firstly, osteosarcoma metastasis-related differentially expressed genes (DEGs) and infiltrating immune cells in the TME were analyzed using a series of bioinformatics methods. The metastasis-related gene signature (MRS) and TME score of osteosarcoma patients were then developed. On this basis, a personalized MRS-TME classifier was constructed and validated in other clinical cohorts and different subgroups. In addition, the relationship between the MRS-related genes and the immune microenvironment was also clarified. Finally, the signaling pathways and immune response genes in osteosarcoma patients among different MRS-TME subgroups were analyzed to explore the underlying molecular mechanism.

Results

We first identified the metastasis-related DEGs in osteosarcoma, which were primarily involved in the muscle system process, calcium ion homeostasis, cell chemotaxis, and leukocyte migration. A personalized MRS-TME classifier was then constructed by integrating the MRS (10 genes) and TME (six immune cells) scores. The MRS-TME classifier demonstrated a potent capacity of predicting the survival prognosis in diverse osteosarcoma cohorts as well as in the clinical feature subgroups. The MRS score was negatively associated with the TME score, and patients in the MRS^low/TME^high subgroup exhibited a better prognosis compared to all other subgroups. Significant differences existed between the cellular signaling pathways and immune response profiles among the different MRS-TME subgroups, especially in relation to the metabolism-related biological processes and the inflammatory response.

Conclusions

The MRS-TME classifier might be a beneficial tool to aid in the prognostic evaluation and risk stratification of osteosarcoma patients.

Characterization of osteosarcoma subtypes mediated by macrophage-related genes and creation and validation of a risk score system to quantitatively assess the prognosis of osteosarcoma and reflect the tumor microenvironment

Background

Macrophages are the main immune components in the microenvironment of osteosarcoma. The treatment strategy centered on macrophages has become a hot topic to improve cancer treatment. However, the research on the role of macrophages in the treatment of osteosarcoma is still in its infancy.

Methods

The data of osteosarcoma samples were downloaded from the Therapeutically Applicable Research to Generate Effective Treatments (TARGET) and GSE21257 datasets, and the macrophage enrichment fraction of osteosarcoma samples in TARGET was calculated by single-sample gene set enrichment analysis (ssGSEA) method to screen macrophage-related genes for consensus clustering. Differential expression analysis, univariable Cox, and least absolute shrinkage and selection operator (LASSO) regression were conducted to select reliable predictors and create a risk score system. The GSE21257 dataset was used as a verification set to verify the accuracy of risk score system.

Results

We identified 2 osteosarcoma clusters mediated by 22 macrophage score-related genes, namely cluster 1 (C1) and cluster 2 (C2). Compared with C2, C1 had a significant advantage in prognosis, and the degree of immune cell infiltration in tumor microenvironment (TME) was significantly higher, the expression of immune checkpoint molecules was significantly enhanced, and the Tumor Immune Dysfunction and Exclusion (TIDE) score was also significantly down-regulated. A robust risk score system was presented and validated, which demonstrated accuracy and independence in assessing the risk of death of osteosarcoma. The risk score system could also monitor TME infiltration in osteosarcoma samples and showed a close relationship with osteosarcoma biology, including metastasis and immunity.

Conclusions

We identified 2 types of clusters mediated by macrophage-related genes and helped to analyze the cluster suitable for immunotherapy. A new prognostic risk score system was created to quantitatively evaluate the prognosis and TME of osteosarcoma, and to provide a new entry point for the design of personalized treatment.

A novel immune cell signature for predicting osteosarcoma prognosis and guiding therapy

Dysregulation of immune cell infiltration in the tumor microenvironment contributes to the progression of osteosarcoma (OS). In the present study, we explored genes related to immune cell infiltration and constructed a risk model to predict the prognosis of and guide therapeutic strategies for OS. The gene expression profile of OS was obtained from TARGET and Gene Expression Omnibus, which were set as the discovery and verification cohorts. CIBERSORT and Kaplan survival analyses were used to analyze the effects of immune cells on the overall survival rates of OS in the discovery cohort. Differentially expressed gene (DEG) analysis and protein–protein interaction (PPI) networks were used to analyze genes associated with immune cell infiltration. Cox regression analysis was used to select key genes to construct a risk model that classified OS tissues into high- and low-risk groups. The prognostic value of the risk model for survival and metastasis was analyzed by Kaplan–Meier survival analyses, receiver operating characteristic curves, and immunohistochemical experiments. Immunological characteristics and response effects of immune checkpoint blockade (ICB) therapy in OS tissues were analyzed using the ESTIMATE and Tumor Immune Dysfunction and Exclusion algorithms, while sensitivity for both targeted and chemotherapy drugs was analyzed using the OncoPredict algorithm. It was demonstrated that the high infiltration of resting dendritic cells in OS tissues was associated with poor prognosis. A total of 225 DEGs were found between the high- and low-infiltration groups of OS tissues, while 94 genes interacted with others. Through COX analyses, among these 94 genes, four genes (including AOC3, CDK6, COL22A1, and RNASE6) were used to construct a risk model. This risk model showed a remarkable prognostic value for survival rates and metastasis in both the discovery and verification cohorts. Even though a high microsatellite instability score was observed in the high-risk group, the ICB response in the high-risk group was poor. Furthermore, using OncoPredict, we found that the high-risk group OS tissues were resistant to seven drugs and sensitive to 25 drugs. Therefore, our study indicates that the resting dendritic cell signature constructed by AOC3, CDK6, COL22A1, and RNASE6 may contribute to predicting osteosarcoma prognosis and thus therapy guidance.

Characterization of the Tumor Microenvironment in Osteosarcoma Identifies Prognostic- and Immunotherapy-Relevant Gene Signatures

The osteosarcoma (OS) microenvironment is composed of tumor cells, immune cells, and stromal tissue and is emerging as a pivotal player in OS development and progression. Thus, microenvironment-targeted strategies are urgently needed to improve OS treatment outcomes. Using principal component analysis (PCA), we systematically examined the tumor microenvironment (TME) and immune cell infiltration of 88 OS cases and constructed a TME scoring system based on the TMEscore high and TMEscore low phenotypes. Our analysis revealed that TMEscore high correlates with longer survival in OS patients, elevated immune cell infiltration, increased immune checkpoints, and increased sensitivity to chemotherapy. TMEscore low strongly correlated with immune exclusion. These observations were externally validated using a GEO dataset (GSE21257) from 53 OS patients. Our laboratory data also proved our findings. This finding enhances our understanding of the immunological landscape in OS and may uncover novel targeted therapeutic strategies.

Osteoclast-like stromal giant cells in breast cancer likely belong to the spectrum of immunosuppressive tumor-associated macrophages

Background: Breast cancer with osteoclast-like stromal giant cells (OSGC) is an exceedingly rare morphological pattern of invasive breast carcinoma. The tumor immune microenvironment (TIME) of these tumors is populated by OSGC, which resemble osteoclasts and show a histiocytic-like immunophenotype. Their role in breast cancer is unknown. The osteoclast maturation in the bone is regulated by the expression of cytokines that are also present in the TIME of tumors and in breast cancer tumor-associated macrophages (TAMs). TAMs-mediated anti-tumor immune pathways are regulated by miRNAs akin to osteoclast homeostasis. Here, we sought to characterize the different cellular compartments of breast cancers with OSGC and investigate the similarities of OSGC with tumor and TIME in terms of morphology, protein, and miRNA expression, specifically emphasizing on monocytic signatures.

Methods and Results: Six breast cancers with OSGC were included. Tumor-infiltrating lymphocytes (TILs) and TAMs were separately quantified. The different cellular populations (i.e., normal epithelium, cancer cells, and OSGC) were isolated from tissue sections by laser-assisted microdissection. After RNA purification, 752 miRNAs were analyzed using a TaqMan Advanced miRNA Low-Density Array for all samples. Differentially expressed miRNAs were identified by computing the fold change (log2Ratio) using the Kolmogorov-Smirnov test and p values were corrected for multiple comparisons using the false discovery rate (FDR) approach. As a similarity analysis among samples, we used the Pearson test. The association between pairs of variables was investigated using Fisher exact test. Classical and non-classical monocyte miRNA signatures were finally applied. All OSGC displayed CD68 expression, TILs (range, 45–85%) and high TAMs (range, 35–75%). Regarding the global miRNAs profile, OSGC was more similar to cancer cells than to non-neoplastic ones. Shared deregulation of miR-143-3p, miR-195-5p, miR-181a-5p, and miR-181b-5p was observed between OSGC and cancer cells. The monocyte-associated miR-29a-3p and miR-21-3p were dysregulated in OSGCs compared with non-neoplastic or breast cancer tissues.

Conclusion: Breast cancers with OSGC have an activated TIME. Shared epigenetic events occur during the ontogenesis of breast cancer cells and OSGC but the innumophenotype and miRNA profiles of the different cellular compartmens suggest that OSGC likely belong to the spectrum of M2 TAMs.

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.

The Role of NR4A1 in the Pathophysiology of Osteosarcoma: A Comprehensive Bioinformatics Analysis of the Single-Cell RNA Sequencing Dataset

Osteosarcoma is a kind of aggressive human malignancy, and the prognosis of the patients with osteosarcoma remains low. Studies have demonstrated that the tumor microenvironment plays a key role in regulating osteosarcoma progression. Recent studies have also shown that scRNA-seq plays an essential role in understanding the tumor heterogeneity and distinct subpopulations of tumors. In order to further understand the scRNA-seq data of osteosarcoma tissues, the present study further analyzed the scRNA-seq dataset (GSE152048) and explored the potential role of nuclear receptor-related genes in the pathophysiology of osteosarcoma. In our analysis, we identified 11 cell types in all the osteosarcoma tissues and nuclear receptors (NRs) were distributed in all types of cells. Further stratification analysis showed that NRs were mainly detected in “TIL” and “Osteoblastic” of the metastasis osteosarcoma, in “TIL”, “Myoblast”, “Endothelial”, and “Myeloid” of the primary osteosarcoma, and in “Chondroblastic”, “Osteoblast”, and “Pericyte” of the recurrent osteosarcoma. The NRs were also differentially expressed in different cell types among the metastasis, primary, and recurrent osteosarcoma. Furthermore, several NRs such as NR4A2, NR4A1, and NR3C1 have been found to be differentially expressed in most types of DEGs among metastasis, primary, and recurrent osteosarcoma. A high expression of NR4A1 in the osteosarcoma tissues was significantly correlated with a shorter 5-year overall survival of patients with osteosarcoma. On the other hand, there was no significant association between NR4A2 expression and the 5-year overall survival of patients with osteosarcoma. The expression of NR4A1 was significantly higher in the metastasis osteosarcoma tissues than in the primary osteosarcoma tissues as validated from GSE32981 and GSE154540. The expression of NR4A1 was significantly higher in osteosarcoma tissues from patients with poor chemosensitivity than that from patients with good chemosensitivity as validated from GSE154540. Further analysis of the scRNA-seq data revealed that the percentage of osteoblasts with a high NR4A1 expression was higher in the recurrent osteosarcoma tissues than that with a low NR4A1 expression. In conclusion, the present study may suggest that NR4A1 may be an important prognostic biomarker for osteosarcoma progression. However, further validation studies should be performed to confirm our findings.

Paediatric Strategy Forum for medicinal product development of multi-targeted kinase inhibitors in bone sarcomas: ACCELERATE in collaboration with the European Medicines Agency with participation of the Food and Drug Administration

The eighth Paediatric Strategy Forum focused on multi-targeted kinase inhibitors (mTKIs) in osteosarcoma and Ewing sarcoma. The development of curative, innovative products in these tumours is a high priority and addresses unmet needs in children, adolescents and adults. Despite clinical and investigational use of mTKIs, efficacy in patients with bone tumours has not been definitively demonstrated. Randomised studies, currently being planned or in progress, in front-line and relapse settings will inform the further development of this class of product. It is crucial that these are rapidly initiated to generate robust data to support international collaborative efforts. The experience to date has generally indicated that the safety profile of mTKIs as monotherapy, and in combination with chemotherapy or other targeted therapy, is consistent with that of adults and that toxicity is manageable. Increasing understanding of relevant predictive biomarkers and tumour biology is absolutely critical to further develop this class of products. Biospecimen samples for correlative studies and biomarker development should be shared, and a joint academic-industry consortium created. This would result in an integrated collection of serial tumour tissues and a systematic retrospective and prospective analyses of these samples to ensure robust assessment of biologic effect of mTKIs. To support access for children to benefit from these novel therapies, clinical trials should be designed with sufficient scientific rationale to support regulatory and payer requirements. To achieve this, early dialogue between academia, industry, regulators, and patient advocates is essential. Evaluating feasibility of combination strategies and then undertaking a randomised trial in the same protocol accelerates drug development. Where possible, clinical trials and development should include children, adolescents, and adults less than 40 years. To respond to emerging science, in approximately 12 months, a multi-stakeholder group will meet and review available data to determine future directions and priorities.

Polydatin Incorporated in Polycaprolactone Nanofibers Improves Osteogenic Differentiation

Polycaprolactone nanofibers are used as scaffolds in the field of tissue engineering for tissue regeneration or drug delivery. Polycaprolactone (PCL) is a biodegradable hydrophobic polyester used to obtain implantable nanostructures, which are clinically applicable due to their biological safety. Polydatin (PD), a glycosidic precursor of resveratrol, is known for its antioxidant, antitumor, antiosteoporotic, and bone regeneration activities. We aimed to use the osteogenic capacity of polydatin to create a biomimetic innovative and patented scaffold consisting of PCL-PD for bone tissue engineering. Both osteosarcoma cells (Saos-2) and mesenchymal stem cells (MSCs) were used to test the in vitro cytocompatibility of the PD-PCL scaffold. Reverse-phase (RP) HPLC was used to evaluate the timing release of PD from the PCL-PD nanofibers and the MTT assay, scanning electron microscopy, and alkaline phosphatase (ALP) activity were used to evaluate the proliferation, adhesion, and cellular differentiation in both osteosarcoma and human mesenchymal stem cells (MSCs) seeded on PD-PCL nanofibers. The proliferation of osteosarcoma cells (Saos-2) on the PD-PCL scaffold decreased when compared to cells grown on PLC nanofibers, whereas the proliferation of MSCs was comparable in both PCL and PD-PCL nanofibers. Noteworthy, after 14 days, the ALP activity was higher in both Saos-2 cells and MSCs cultivated on PD-PCL than on empty scaffolds. Moreover, the same cells showed a spindle-shaped morphology after 14 days when grown on PD-PCL as shown by SEM. In conclusion, we provide evidence that nanofibers appropriately coated with PD support the adhesion and promote the osteogenic differentiation of both human osteosarcoma cells and MSCs.

Baicalin induces apoptosis and autophagy in human osteosarcoma cells by increasing ROS to inhibit PI3K/Akt/mTOR, ERK1/2 and β-catenin signaling pathways

•

Baicalin causes apoptosis and autophagy through accumulating ROS to suppress PI3K/Akt/mTOR, ERK1/2 and β-catenin pathways in OS cells.

•

Baicalin-induced autophagosome further triggers apoptosis.

•

Baicalin-induced ROS and Ca²⁺ interactions induce apoptosis.

•

Baicalin molecule targets PI3Kγ, inhibiting downstream effectors AKT and mTOR.

Baicalin, a flavonoid derivative, exerts antitumor activity in a variety of neoplasms. However, whether baicalin exerts antitumor effects on osteosarcoma cells remains to be elucidated. In this study, treatment with baicalin reduced the proliferation and invasive potential of osteosarcoma cells and reduced the mitochondrial membrane potential, which eventually caused mitochondrial apoptosis. In addition, baicalin increased intercellular Ca²⁺ and ROS concentrations. Baicalin-induced apoptosis was confirmed by enhanced Bax, cleaved caspase-3, and cleaved PARP levels and decreased Bcl-2 levels. The increase in LC3-II and p62 suggested that baicalin induced autophagosome formation but ultimately inhibited downstream autophagy. Moreover, apoptosis induced by baicalin was attenuated by the addition of 3-MA. Furthermore, we found that baicalin inhibited the PI3K/Akt/mTOR, ERK1/2 and β-catenin signaling pathways. Chelation of free Ca²⁺ by BAPTA-AM also inhibited both apoptosis induction and ROS concentration changes. Finally, NAC pretreatment reversed baicalin treatment outcomes, including the increase in Ca²⁺ concentration, induction of apoptosis and autophagy, and inhibition of the pathways. Molecular docking results indicated that baicalin might interact with the structural domain of PI3Kγ. Thus, baicalin may be considered a potential candidate for osteosarcoma treatment.

ATG16L1 is a Potential Prognostic Biomarker and Immune Signature for Osteosarcoma: A Study Based on Bulk RNA and Single-Cell RNA-Sequencing

Background

Osteosarcoma is a common solid malignancy of the bone in children and adolescents, and its metastasis and recurrence are the principal causes of poor treatment outcomes.

Methods

Autophagy-related genes were used to cluster osteosarcoma patients by consensus clustering analysis using the GSE21257 database. Differentially expressed genes (DEGs) were identified by limma package. Multiple-gene risk signature was constructed using least absolute shrinkage and selection operator (LASSO) analysis and Cox regression analyses. Quantitative reverse transcription polymerase chain reaction (qRT-PCR) was used to determine gene expression levels. Then, single-cell RNA-sequencing dataset GSE152048 were used to identify the correlation between the DEGs and effector molecules expressed in specific tumor-infiltrating immune cells.

Results

Two clusters were identified in the consensus clustering analysis, which were confirmed by principal component analysis. Limma analysis revealed that 15 genes were related, and 9 genes were screened using protein-protein interaction network and LASSO regression analysis. Cox regression analyses identified 5 genes. Combined with survival analysis, only the autophagy related 16 like 1 gene (ATG16L1) was significant. The results of qRT-PCR showed low expression levels of ATG16L1 in tumor cells group. Immune infiltration analysis revealed significantly lower expression of CD8⁺ T cells in the high ATG16L1 gene expression group. ScRNA-seq revealed that in the ATG16L1⁺CD8⁺ T cell group, the expression of GZMB was lower, whereas the expression of ITGA1 was higher. These results showed that ATG16L1 is an immune-related gene, which is associated with poor prognosis in patients with osteosarcoma.

Conclusion

ATG16L1 is a potential prognostic biomarker and immune signature and may be a therapeutic target for osteosarcoma.

Dissecting the effect of sphingolipid metabolism gene in progression and microenvironment of osteosarcoma to develop a prognostic signature

Sphingolipid metabolism (SM) fuels tumorigenesis and the malignant progression of osteosarcoma (OS), which leads to an unfavorable prognosis. Elucidating the molecular mechanisms underlying SM in osteosarcoma and developing a SM-based prognostic signature could be beneficial in the clinical setting. This study included 88 frozen OS samples to recognize the vital SM-relevant genes in the development of OS utilizing univariate Cox regression. The Least Absolute Shrinkage and Selection Operator (LASSO) regression analysis was conducted on the SM- relevant genes to minimize the risk of overfitting. The prognostic signature was generate utilizing the multivariable Cox regression analysis and was verified in the validation cohort. Moreover, cellular and molecular mechanisms associated with SM have an unfavorable prognosis for OS patients and have been widely studied. Resultantly, an SM-based prognostic risk model was established according to critical prognostic genes (CBS, GLB1, and HACD1), which had an excellent ability to predict the prognosis of OS patients (AUC for the train cohort was 0.887 and AUC for validation cohort was 0.737). The high-risk OS patients identified based on this prognostic signature had significantly poor immune microenvironment, indicated by significantly low immune score (mean=216.290 ± 662.463), reduced infiltrations of 25 immune cells, including NK cells (LogFC= -0.3597), CD8+T cells ((LogFC=-0.2346), Cytolytic activity ((LogFC=-0.1998), etc. The immunosuppressive microenvironment could be due to dysregulated SM of glycolipids. Further, a nomogram was constructed by integrating the SM-based prognostic signature and clinical paraments to facilitate clinical application. The nomogram could accurately predict the prognosis of OS invalids. Collectively, this study clarified the function of SM in the development of OS and helped develop a tool for risk stratification based on SM-related genes with application in clinical settings. The results of our study will aid in identifying high-risk patients and provide individualized treatments.