Emerging next-generation sequencing-based discoveries for targeted osteosarcoma therapy

Discovery of core genes and intercellular communication role in osteosarcoma

Osteosarcoma is a primary malignant bone tumor that affects children and young adults. Understanding the molecular mechanisms underlying osteosarcoma is critical to develop effective treatments. This study aimed to identify core genes and explore the role of intercellular communication in osteosarcoma. We used GSE87437 and GSE152048 dataset to conduct a weighted correlation network analysis (WGCNA) and identify co-expression modules. The enriched biological processes and cellular components of the genes in the steelblue module were analyzed. Next, we explored the expression, diagnostic value, correlation, and association with immune infiltrate of CCSER1 and LOC101929154. Finally, we utilized CIBERSORT algorithm to predict the infiltrated immune cells in osteosarcoma tissues. Our results identified 44 co-expression modules, and the steelblue module was mainly associated with axon development, axonogenesis, and innervation. CCSER1 and LOC101929154 were significantly upregulated in osteosarcoma tissues with poor response to preoperative chemotherapy. Moreover, the expressions of CCSER1 and LOC101929154 were positively correlated. The area under the receiver operating characteristic curve of CCSER1 and LOC101929154 was 0.800 and 0.773, respectively. The expression of CCSER1 was negatively correlated with follicular helper T cells and positively correlated with M0 macrophages, while LOC101929154 was negatively correlated with activated mast cells. Besides, CD4 memory-activated T cells were observed at lower levels in patients who responded well to chemotherapy. Our study identified core genes CCSER1 and LOC101929154 and provided insight into the intercellular communication profile in osteosarcoma. Our results suggested that targeting CCSER1, LOC101929154, and CD4 memory-activated T cells may be a promising strategy for the treatment of osteosarcoma.

Targeting the mitochondrial protein YME1L to inhibit osteosarcoma cell growth in vitro and in vivo

Exploring novel diagnostic and therapeutic biomarkers is extremely important for osteosarcoma. YME1 Like 1 ATPase (YME1L), locating in the mitochondrial inner membrane, is key in regulating mitochondrial plasticity and metabolic activity. Its expression and potential functions in osteosarcoma are studied in the present study. We show that YME1L mRNA and protein expression is significantly elevated in osteosarcoma tissues derived from different human patients. Moreover, its expression is upregulated in various primary and immortalized osteosarcoma cells. The Cancer Genome Atlas database results revealed that YME1L overexpression was correlated with poor overall survival and poor disease-specific survival in sarcoma patients. In primary and immortalized osteosarcoma cells, silencing of YME1L through lentiviral shRNA robustly inhibited cell viability, proliferation, and migration. Moreover, cell cycle arrest and apoptosis were detected in YME1L-silenced osteosarcoma cells. YME1L silencing impaired mitochondrial functions in osteosarcoma cells, causing mitochondrial depolarization, oxidative injury, lipid peroxidation and DNA damage as well as mitochondrial respiratory chain complex I activity inhibition and ATP depletion. Contrarily, forced YME1L overexpression exerted pro-cancerous activity and strengthened primary osteosarcoma cell proliferation and migration. YME1L is important for Akt-S6K activation in osteosarcoma cells. Phosphorylation of Akt and S6K was inhibited after YME1L silencing in primary osteosarcoma cells, but was strengthened with YME1L overexpression. Restoring Akt-mTOR activation by S473D constitutively active Akt1 mitigated YME1L shRNA-induced anti-osteosarcoma cell activity. Lastly, intratumoral injection of YME1L shRNA adeno-associated virus inhibited subcutaneous osteosarcoma xenograft growth in nude mice. YME1L depletion, mitochondrial dysfunction, oxidative injury, Akt-S6K inactivation, and apoptosis were detected in YME1L shRNA-treated osteosarcoma xenografts. Together, overexpressed YME1L promotes osteosarcoma cell growth, possibly by maintaining mitochondrial function and Akt-mTOR activation.

Advances in prognostic models for osteosarcoma risk

The risk prognosis model is a statistical model that uses a set of features to predict whether an individual will develop a specific disease or clinical outcome. It can be used in clinical practice to stratify disease severity and assess risk or prognosis. With the advancement of large-scale second-generation sequencing technology, along Prognosis models for osteosarcoma are increasingly being developed as large-scale second-generation sequencing technology advances and clinical and biological data becomes more abundant. This expansion greatly increases the number of prognostic models and candidate genes suitable for clinical use. This article will present the predictive effects and reliability of various prognosis models, serving as a reference for their evaluation and application.

ACSL4 accelerates osteosarcoma progression via modulating TGF-β/Smad2 signaling pathway

Osteosarcoma (OS) is a malignant bone sarcoma arising from mesenchymal stem cells. The biological role of Acyl-CoA synthetase long-chain family member 4 (ACSL4), recently identified as an oncogene in numerous tumor types, remains largely unclear in OS. In this study, we investigated the expression of ACSL4 in OS tissues using immunohistochemistry staining (IHC) staining of a human tissue microarray and in OS cells by qPCR assay. Our findings revealed a significant up-regulation of ACSL4 in both OS tissues and cells. To further understand its biological effects, we conducted a series of loss-of-function experiments using ACSL4-depleted MNNG/HOS and U-2OS cell lines, focusing on OS cell proliferation, migration, and apoptosis in vitro. Our results demonstrated that ACSL4 knockdown remarkably suppressed OS cell proliferation, arrested cells in the G2 phase, induced cell apoptosis, and inhibited cell migration. Additionally, a subcutaneous xenograft mice model was established to validate the in vivo impact of ACSL4, revealing ACSL4 silencing impaired tumor growth in the OS xenograft mice. Additionally, we discovered that ACSL4 could regulate the phosphorylation level of Smad2 through cooperative interactions, and treatment with a TGF-β inhibitor weakened the promoting effects of ACSL4 overexpression. In short, ACSL4 regulated OS progression by modulating TGF-β/Smad2 signaling pathway. These findings underscore ACSL4 as a promising therapeutic target for OS patients and contribute novel insights into the pathogenesis of OS.

Targeting POLRMT by a first-in-class inhibitor IMT1 inhibits osteosarcoma cell growth in vitro and in vivo

Osteosarcoma (OS) is a highly aggressive form of bone cancer that predominantly affects adolescents and young adults. In this study, we have undertaken an investigation into the potential anti-OS cell activity of IMT1 (inhibitor of mitochondrial transcription 1), a first-in-class inhibitor of RNA polymerase mitochondrial (POLRMT). IMT1 exhibited a profound inhibitory effect on cell survival, proliferation, cell cycle progression, and migration in primary and immortalized OS cells. Furthermore, this POLRMT inhibitor elicited apoptosis in the OS cells, without, however, inducing cytotoxicity in human osteoblasts or osteoblastic cells. IMT1 disrupted mitochondrial functions in OS cells, resulting in mitochondrial depolarization, oxidative injury, lipid peroxidation, and ATP reduction in OS cells. Silencing POLRMT using targeted shRNA closely mimicked the actions of IMT1 and exerted potent anti-OS cell activity. Importantly, IMT1's effectiveness was diminished in POLRMT-silenced OS cells. Subsequent investigations revealed that IMT1 suppressed the activation of the Akt-mammalian target of rapamycin (mTOR) cascade in OS cells. IMT1 treatment or POLRMT silencing in primary OS cells led to a significant reduction in Akt1-S6K-S6 phosphorylation. Conversely, it was enhanced upon POLRMT overexpression. The restoration of Akt-mTOR activation through the introduction of a constitutively active S473D mutant Akt1 (caAkt1) mitigated IMT1-induced cytotoxicity in OS cells. In vivo, oral administration of IMT1 robustly curtailed the growth of OS xenografts in nude mice. Furthermore, IMT1 suppressed POLRMT activity, impaired mitochondrial function, repressed Akt-mTOR activation, and induced apoptosis within xenograft tissues. Collectively, these findings underscore the potent growth-inhibitory effects attributed to IMT1 via targeted POLRMT inhibition. The utilization of this POLRMT inhibitor carries substantial therapeutic promise in the context of OS treatment.

Osteosarcoma Explorer: A Data Commons With Clinical, Genomic, Protein, and Tissue Imaging Data for Osteosarcoma Research

PURPOSE

Osteosarcoma research advancement requires enhanced data integration across different modalities and sources. Current osteosarcoma research, encompassing clinical, genomic, protein, and tissue imaging data, is hindered by the siloed landscape of data generation and storage.

MATERIALS AND METHODS

Clinical, molecular profiling, and tissue imaging data for 573 patients with pediatric osteosarcoma were collected from four public and institutional sources. A common data model incorporating standardized terminology was created to facilitate the transformation, integration, and load of source data into a relational database. On the basis of this database, a data commons accompanied by a user-friendly web portal was developed, enabling various data exploration and analytics functions.

RESULTS

The Osteosarcoma Explorer (OSE) was released to the public in 2021. Leveraging a comprehensive and harmonized data set on the backend, the OSE offers a wide range of functions, including Cohort Discovery, Patient Dashboard, Image Visualization, and Online Analysis. Since its initial release, the OSE has experienced an increasing utilization by the osteosarcoma research community and provided solid, continuous user support. To our knowledge, the OSE is the largest (N = 573) and most comprehensive research data commons for pediatric osteosarcoma, a rare disease. This project demonstrates an effective framework for data integration and data commons development that can be readily applied to other projects sharing similar goals.

CONCLUSION

The OSE offers an online exploration and analysis platform for integrated clinical, molecular profiling, and tissue imaging data of osteosarcoma. Its underlying data model, database, and web framework support continuous expansion onto new data modalities and sources.

Integration of single-cell sequencing and bulk expression data reveals chemokine signaling pathway in proliferating cells is associated with the survival outcome of osteosarcoma

BACKGROUND

Osteosarcoma, as the most common primary bone malignancy, is urgent to be well-studied on the biomarkers and therapeutic targets to improve the five-year survival rate. Transcriptomic analysis using single-cell RNA or bulk RNA sequencing has been developed to detect biomarkers in various cancer types.

METHODS AND RESULTS

We applied Scissor to combine single-cell RNA-seq data and bulk transcriptome data of osteosarcoma, providing cell-level information and sample phenotypes to identify the survival-associated cell subpopulations. By investigating the differences between the survival-associated cell subpopulations, we identified CCL21, CCL22, CCL24, CXCL11, CXCL12, CXCL13, GNAI2, and RAC2 in the proliferating cells that are significantly associated with osteosarcoma patient outcome. Then we assigned the risk score for each sample based on the cell proportion-normalized gene expression and validated it in the public dataset.

CONCLUSIONS

This study provides the clinical insight that chemokine signaling pathway genes (CCL21, CCL22, CCL24, CXCL11, CXCL12, CXCL13, GNAI2, and RAC2) in proliferating cells might be the potential biomarkers for treatment of osteosarcoma.

Cancer chemotherapy and beyond: Current status, drug candidates, associated risks and progress in targeted therapeutics

Cancer is an abnormal state of cells where they undergo uncontrolled proliferation and produce aggressive malignancies that causes millions of deaths every year. With the new understanding of the molecular mechanism(s) of disease progression, our knowledge about the disease is snowballing, leading to the evolution of many new therapeutic regimes and their successive trials. In the past few decades, various combinations of therapies have been proposed and are presently employed in the treatment of diverse cancers. Targeted drug therapy, immunotherapy, and personalized medicines are now largely being employed, which were not common a few years back. The field of cancer discoveries and therapeutics are evolving fast as cancer type-specific biomarkers are progressively being identified and several types of cancers are nowadays undergoing systematic therapies, extending patients' disease-free survival thereafter. Although growing evidence shows that a systematic and targeted approach could be the future of cancer medicine, chemotherapy remains a largely opted therapeutic option despite its known side effects on the patient's physical and psychological health. Chemotherapeutic agents/pharmaceuticals served a great purpose over the past few decades and have remained the frontline choice for advanced-stage malignancies where surgery and/or radiation therapy cannot be prescribed due to specific reasons. The present report succinctly reviews the existing and contemporary advancements in chemotherapy and assesses the status of the enrolled drugs/pharmaceuticals; it also comprehensively discusses the emerging role of specific/targeted therapeutic strategies that are presently being employed to achieve better clinical success/survival rate in cancer patients.

Comprehensive analysis of immune implication and prognostic value of DHX33 in sarcoma

DEAH-box helicase 33 (DHX33) is an RNA helicase that has been identified to promote the progression of a variety of cancers. However, the relationship between DHX33 and sarcoma remains unknown. RNA expression data with clinical information for the sarcoma project was collected from TCGA database. The association between the differential expression of DHX33 and the prognosis for sarcoma was assessed using survival analysis. CIBERSORT was used to evaluate the immune cell infiltration in sarcoma sample tissues. We then further investigated the association between DHX33 and tumor-infiltrating immune cells in sarcoma using the TIMER database. Finally, the immune/cancer-related signaling pathways involved in DHX33 were analyzed using gene set enrichment analysis. High DHX33 expression was discovered to be a poor prognostic indicator in TCGA-SARC. Immune subpopulations in the TCGA-SARC microenvironment are dramatically altered compared to normal tissues. The tumor immune estimation resource analysis revealed a strong correlation between the expression of DHX33 and the abundance of CD8+ T cells and dendritic cells. Changes in copy number also affected neutrophils, macrophages, and CD4+ T cells. According to gene set enrichment analysis, DHX33 may be involved in a number of cancer- and immune-related pathways, such as the JAK/STAT signaling pathway, P53 signaling pathway, chemokine signaling pathway, T cell receptor signaling pathway, complement and coagulation cascades, and cytokine-cytokine receptor interaction. Our study emphasized that DHX33 may be involved in the immune microenvironment of sarcoma and play an important role. As a result, it is possible that DHX33 might serve as an immunotherapeutic target for sarcoma.

Expanding individualized therapeutic options via genoproteomics

Clinical next-generation sequencing (NGS)² tests have enabled treatment recommendations for cancer patients with driver gene mutations. Targeted therapy options for patients without driver gene mutations are currently unavailable. Herein, we performed NGS and proteomics tests on 169 formalin-fixed paraffin-embedded (FFPE)³ samples of non-small cell lung cancers (NSCLC, 65),⁴ colorectal cancers (CRC, 61),⁵ thyroid carcinomas (THCA, 14),⁶ gastric cancers (GC, 2),⁷ gastrointestinal stromal tumors (GIST, 11),⁸ and malignant melanomas (MM, 6).⁹ Of the 169 samples, NGS detected 14 actionable mutated genes in 73 samples, providing treatment options for 43% of the patients. Proteomics identified 61 actionable clinical drug targets approved by the FDA or undergoing clinical trials in 122 samples, providing treatment options for 72% of the patients. In vivo experiments demonstrated that the Mitogen-Activated Protein Kinase (MEK)¹⁰ inhibitor induced the overexpression of MEK1 (Map2k1) to block lung tumor growth in mice. Therefore, protein overexpression is a potentially feasible indicator for guiding targeted therapies. Collectively, our analysis suggests that combining NGS and proteomics (genoproteomics) could expand the targeted treatment options to 85% of cancer patients.

A novel molecular classification method for osteosarcoma based on tumor cell differentiation trajectories

Subclassification of tumors based on molecular features may facilitate therapeutic choice and increase the response rate of cancer patients. However, the highly complex cell origin involved in osteosarcoma (OS) limits the utility of traditional bulk RNA sequencing for OS subclassification. Single-cell RNA sequencing (scRNA-seq) holds great promise for identifying cell heterogeneity. However, this technique has rarely been used in the study of tumor subclassification. By analyzing scRNA-seq data for six conventional OS and nine cancellous bone (CB) samples, we identified 29 clusters in OS and CB samples and discovered three differentiation trajectories from the cancer stem cell (CSC)-like subset, which allowed us to classify OS samples into three groups. The classification model was further examined using the TARGET dataset. Each subgroup of OS had different prognoses and possible drug sensitivities, and OS cells in the three differentiation branches showed distinct interactions with other clusters in the OS microenvironment. In addition, we verified the classification model through IHC staining in 138 OS samples, revealing a worse prognosis for Group B patients. Furthermore, we describe the novel transcriptional program of CSCs and highlight the activation of EZH2 in CSCs of OS. These findings provide a novel subclassification method based on scRNA-seq and shed new light on the molecular features of CSCs in OS and may serve as valuable references for precision treatment for and therapeutic development in OS.

Establishment of Rapid and Accurate Screening System for Molecular Target Therapy of Osteosarcoma

Introduction Comprehensive analyses using clinical sequences subcategorized osteosarcoma (OS) into several groups according to the activated signaling pathways. Mutually exclusive co-occurrences of gene amplification (PDGFRA/KIT/KDR, VEGFA/CCND3, and MDM2/CDK4) have been identified in approximately 40% of OS, representing candidate subsets for clinical evaluation of additional therapeutic options. Thus, it would be desirable to evaluate the specific gene amplification before starting therapy in patients with OS. Materials and Methods This is a retrospective study. We examined 13 cases of clinical OS samples using NanoString-based copy number variation (CNV) analysis. Decalcification and chemotherapeutic effects on this analysis were also assessed. Results First, the accuracy of this system was validated by showing that amplification/deletion data obtained from this system using various types of cancer cell lines almost perfectly matched to that from the Cancer Cell Line Encyclopedia (CCLE). We identified potentially actionable alterations in CDK4/MDM2 amplification in 10% of samples and potential additional therapeutic targets (PDGFRA/KIT/KDR and VEGFA/CCND3) in 20% of samples, which is consistent with the reported frequencies. Furthermore, this assay could identify these potential therapeutic targets regardless of the sample status (frozen vs formalin-fixed paraffin-embedded [FFPE] tissues). Conclusion We established a NanoString-based rapid and cost-effective method with a short turnaround time (TAT) to examine gene amplification status in OS. This CNV analysis using FFPE samples is recommended where the histological evaluation of viable tumor cells is possible, especially for tumors after chemotherapy with higher chemotherapeutic effects.

Superenhancers activate the autophagy-related genes Beclin1 and LC3B to drive metastasis and drug resistance in osteosarcoma

Metastasis and drug resistance are the leading causes of poor prognosis in patients with osteosarcoma. Identifying the relevant factors that drive metastasis and drug resistance is the key to improving the therapeutic outcome of osteosarcoma. Here, we reported that autophagy was highly activated in metastatic osteosarcoma. We found increased autophagolysosomes in metastatic osteosarcoma cell lines by using electron microscopy, Western blot, and immunofluorescence experiments. We further examined the expression of the autophagy-related genes Beclin1 and LC3B in 82 patients through immunohistochemistry and found that Beclin1 and LC3B were highly related to unfavorable prognosis of osteosarcoma. Knockdown of Beclin1 and LC3B reduced invasion, metastasis, and proliferation in metastatic osteosarcoma cells. In vitro and in vivo studies also demonstrated that inhibiting by 3-MA inhibited cell growth and metastasis. Moreover, we demonstrated that autophagy-related genes were activated by SEs and that the inhibition of SEs by JQ-1 decreased the metastasis of osteosarcoma. Overall, our findings highlighted the association of autophagy with osteosarcoma progression and shed new light on autophagy-targeting therapy for osteosarcoma.

ADCK1 is a potential therapeutic target of osteosarcoma

We here showed that ADCK1 (AarF domain-containing kinase 1), a mitochondrial protein, is upregulated in human osteosarcoma (OS) tissues and OS cells. In primary and established OS cells, ADCK1 shRNA or CRISPR/Cas9-induced ADCK1 knockout (KO) remarkably inhibited cell viability, proliferation and migration, and provoked apoptosis activation. Conversely, ectopic ADCK1 overexpression exerted pro-cancerous activity by promoting OS cell proliferation and migration. ADCK1 depletion disrupted mitochondrial functions in OS cells and induced mitochondrial membrane potential reduction, ATP depletion, reactive oxygen species production. Significantly, ADCK1 silencing augmented doxorubicin-induced apoptosis in primary OS cells. mTOR activation is important for ADCK1 expression in OS cells. The mTOR inhibitors, rapamycin and AZD2014, as well as mTOR shRNA, potently decreased ADCK1 expression in primary OS cells. In nude mice, the growth of subcutaneous pOS-1 xenografts was largely inhibited when bearing ADCK1 shRNA or ADCK1 KO construct. Moreover, ADCK1 KO largely inhibited pOS-1 xenograft in situ growth in proximal tibia of nude mice. ADCK1 depletion, apoptosis activation and ATP reduction were detected in pOS-1 xenografts bearing ADCK1 shRNA or ADCK1 KO construct. Together, the mitochondrial protein ADCK1 is required for OS cell growth and is a novel therapeutic target of OS.

Transcriptome Profiling of Different State Callus Induced from Immature Embryo in Maize

Embryogenic and regenerable tissue cultures are widely used in plant transformation. To dissect the molecular mechanism of embryogenesis, we used inbred line A188 as the material; the immature embryo of kernels (15 day after pollination, 15DAP) was isolated and cultured in inducing medium and subjected to RNA-Seq. The results revealed that 5,076 differentially expressed genes (DEGs) were involved in morphological and histological changes and endogenous indole-3-acetic acid (IAA) alteration. Functional analysis showed that the DEGs were related to metabolic pathways and biosynthesis of secondary metabolites. In particular, ARF16 and ARF8 genes of auxin response factors (ARF) were upregulated from EC to IDC and EC to IRC. Meanwhile, BBM2, SERK1, and SERK2 genes of the embryogenic pathway were upregulated, and WIP2 and ESR genes of the wound-inducible were upregulated from EC to IDC and EC to IRC. These changes can improve conversion efficiency from EC to IRC, which is important for elucidating the underlying molecular mechanisms of callus formation.

Propranolol blocks osteosarcoma cell cycle progression, inhibits angiogenesis and slows xenograft growth in combination with cisplatin-based chemotherapy

Osteosarcoma is still associated with limited response to standard-of-care therapy and alarmingly elevated mortality rates, especially in low- and middle-income countries. Despite multiple efforts to repurpose β-blocker propranolol in oncology, its potential application in osteosarcoma management remains largely unexplored. Considering the unsatisfied clinical needs of this aggressive disease, we evaluated the antitumoral activity of propranolol using different in vitro and in vivo osteosarcoma preclinical models, alone or in addition to chemotherapy. Propranolol significantly impaired cellular growth in β2-adrenergic receptor-expressing MG-63 and U-2OS cells, and was capable of blocking growth-stimulating effects triggered by catecholamines. siRNA-mediated ADRB2 knockdown in MG-63 cells was associated with decreased cell survival and a significant attenuation of PPN anti-osteosarcoma activity. Direct cytostatic effects of propranolol were independent of apoptosis induction and were associated with reduced mitosis, G0/G1 cell cycle arrest and a significant down-regulation of cell cycle regulator Cyclin D1. Moreover, colony formation, 3D spheroid growth, cell chemotaxis and capillary-like tube formation were drastically impaired after propranolol treatment. Interestingly, anti-migratory activity of β-blocker was associated with altered actin cytoskeleton dynamics. In vivo, propranolol treatment (10 mg/kg/day i.p.) reduced the early angiogenic response triggered by MG-63 cells in nude mice. Synergistic effects were observed in vitro after combining propranolol with chemotherapeutic agent cisplatin. Sustained administration of propranolol (10 mg/kg/day i.p., five days a week), alone and especially in addition to low-dose metronomic cisplatin (2 mg/kg/day i.p., three times a week), markedly reduced xenograft progression. After histological analysis, propranolol and cisplatin combination resulted in low tumor mitotic index and increased tumor necrosis. β-blockade using propranolol seems to be an achievable and cost-effective therapeutic approach to modulate osteosarcoma aggressiveness. Further translational studies of propranolol repurposing in osteosarcoma are warranted.

Colchicine as a novel drug for the treatment of osteosarcoma through drug repositioning based on an FDA drug library

Background

Colchicine is a traditional medication that is currently approved to treat gout and familial Mediterranean fever (FMF). However, colchicine has a wide range of anti-inflammatory activities, and several studies have indicated that it may be useful in a variety of other conditions, such as rheumatic disease, cardiac disease, and cancer. Osteosarcoma, the most common type of bone sarcoma, is derived from primitive bone-forming mesenchymal cells. In this study, we investigated whether colchicine could be used to treat osteosarcoma through the regulation of cell cycle signaling.

Methods

Two human osteosarcoma cell lines, U2OS and Saos-2, were used. A clonogenic assay was used to determine the antiproliferative effects of colchicine on osteosarcoma cells. Reactive oxygen species (ROS) production and apoptosis were measured by flow cytometry. Migration and invasion assays were performed to investigate the inhibitory effects of colchicine. The signaling pathways related to colchicine treatment were verified by GO biological process (GOBP) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses.

Results

Colchicine was selected as the lead compound based on the results of initial screening and cell viability assays conducted in Saos-2 and U2Os cells. Colchicine reduced the viability of Saos-2 and U2OS cells in a concentration-dependent manner. It also significantly inhibited colony-forming ability and induced ROS production and apoptosis. It also inhibited the migration and invasion of both Saos-2 and U2OS cells. GOBP and KEGG enrichment analyses indicated the involvement of microtubule-based processes and cancer-related pathways.

Conclusions

These findings suggest that colchicine has therapeutic potential in osteosarcoma.

m6A-Related lncRNAs Predict Overall Survival of Patients and Regulate the Tumor Immune Microenvironment in Osteosarcoma

Background

m6A-related lncRNAs have demonstrated great potential tumor diagnostic and therapeutic targets. The goal of this work was to find m6A-regulated lncRNAs in osteosarcoma patients.

Method

The Cancer Genome Atlas (TCGA) database was used to retrieve RNA sequencing and medical information from osteosarcoma sufferers. The Pearson's correlation test was used to identify the m6A-related lncRNAs. A risk model was built using univariate and multivariable Cox regression analysis. Kaplan–Meier survival analysis and receiver functional requirements were used to assess the risk model's performance (ROC). By using the CIBERSORT method, the associations between the relative risks and different immune cell infiltration were investigated. Lastly, the bioactivities of high-risk and low-risk subgroups were investigated using Gene Set Enrichment Analysis (GSEA).

Result

A total of 531 m6A-related lncRNAs were obtained from TCGA. Seven lncRNAs have demonstrated prognostic values. A total of 88 OS patients were separated into cluster 1, cluster 2, and cluster 3. The overall survival rate of OS patients in cluster 3 was more favorable than that of those in cluster 1 and cluster 2. The average Stromal score was much higher in cluster 1 than in cluster 2 and cluster 3 (P < 0.05). The expression levels of lncRNAs used in the construction of the risk prediction model in the high-risk group were generally lower than those in the low-risk group. Analysis of patient survival indicated that the survival of the low-risk group was higher than that of the high-risk group (P < 0.0001) and the area under the curve (AUC) of the ROC curve was 0.719. Using the CIBERSORT algorithm, the results revealed that Macrophages M0, Macrophages M2, and T cells CD4 memory resting accounted for a large proportion of immune cell infiltration. By GSEA analysis, our results implied that the high-risk group was mainly involved in unfolded protein response, DNA repair signaling, and epithelial-mesenchymal transition signaling pathway and glycolysis pathway; meanwhile, the low-risk group was mainly involved in estrogen response early and KRAS signaling pathway.

Conclusion

Our investigation showed that m6A-related lncRNAs remained tightly connected to the immunological microenvironment of osteosarcoma tumors, potentially influencing carcinogenesis and development. The immune microenvironment and immune-related biochemical pathways can be changed by regulating the transcription of M6A modulators or lncRNAs. In addition, we looked for risk-related signaling of m6A-related lncRNAs in osteosarcomas and built and validated the risk prediction system. The findings of our current analysis will facilitate the assessment of outcomes and the development of immunotherapies for sufferers of osteosarcomas.

Developing a 5-Gene Signature Related to Pyroptosis for Osteosarcoma Patients

Although the incidence of osteosarcoma (OS) is relatively low compared with other cancer types, the overall survival of metastatic OS was less than 30%. This study aimed to reveal the role of pyroptosis in osteosarcoma and develop a prognostic model related to pyroptosis. Weighted correlation network analysis (WGCNA) was applied to identify key gene modules related to pyroptosis. Univariate Cox regression analysis was used to screen prognostic genes related to pyroptosis. The least absolute shrinkage and selection operator (LASSO) and stepwise Akaike information criterion (stepAIC) were employed to optimize and construct a prognostic model. Five prognostic genes (COL13A1, TNFRSF1A, LILRA6, CTNNBIP1, and CD180) related to pyroptosis were identified. According to the 5-gene signature, OS samples were divided into high- and low-PPRS groups with differential prognosis. Immune-related pathways were more activated in the low-PPRS group. The 5-gene signature was effective and robust to predict OS prognosis. These five prognostic genes were involved in OS development and may serve as new targets for developing therapeutic drugs.

Emerging Applications of Deep Learning in Bone Tumors: Current Advances and Challenges

Deep learning is a subfield of state-of-the-art artificial intelligence (AI) technology, and multiple deep learning-based AI models have been applied to musculoskeletal diseases. Deep learning has shown the capability to assist clinical diagnosis and prognosis prediction in a spectrum of musculoskeletal disorders, including fracture detection, cartilage and spinal lesions identification, and osteoarthritis severity assessment. Meanwhile, deep learning has also been extensively explored in diverse tumors such as prostate, breast, and lung cancers. Recently, the application of deep learning emerges in bone tumors. A growing number of deep learning models have demonstrated good performance in detection, segmentation, classification, volume calculation, grading, and assessment of tumor necrosis rate in primary and metastatic bone tumors based on both radiological (such as X-ray, CT, MRI, SPECT) and pathological images, implicating a potential for diagnosis assistance and prognosis prediction of deep learning in bone tumors. In this review, we first summarized the workflows of deep learning methods in medical images and the current applications of deep learning-based AI for diagnosis and prognosis prediction in bone tumors. Moreover, the current challenges in the implementation of the deep learning method and future perspectives in this field were extensively discussed.

No Effect on Tumorigenesis in MG63 Cells Induced by Co-Cultured Mesenchymal Stem Cells

Osteosarcoma is a kind of bone tumor with an extremely high malignant degree and often occurs in adolescents. Mesenchymal stem cells are believed to play an important role in the microenvironment of osteosarcoma, but whether they promote or inhibit the development of osteosarcoma is controversial. In this study, the coexpression of mesenchymal stem cells (MSCs) with osteosarcoma cell MG63 was used to explore the effect of MSCs on MG63. We found that co-culture of MSCs with MG63 did not affect the proliferation, invasion, and migration of MG63 cells, nor did it significantly affect the epithelial- and glial-mesenchymal transformation of MG63 cells. Therefore, in this study, we obtained a new concept that MSCs neither promote nor inhibit the occurrence and development of osteosarcoma.

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.

Single-Cell Profiling of Tumor Microenvironment Heterogeneity in Osteosarcoma Identifies a Highly Invasive Subcluster for Predicting Prognosis

Osteosarcoma is the most common malignant bone tumor in adolescents, and metastasis is the key reason for treatment failure and poor prognosis. Once metastasis occurs, the 5-year survival rate is only approximately 20%, and assessing and predicting the risk of osteosarcoma metastasis are still difficult tasks. In this study, cellular communication between tumor cells and nontumor cells was identified through comprehensive analysis of osteosarcoma single-cell RNA sequencing (scRNA-seq) and bulk RNA-seq data, illustrating the complex regulatory network in the osteosarcoma microenvironment. In line with the heterogeneity of osteosarcoma, we found subpopulations of osteosarcoma cells that highly expressed COL6A1, COL6A3 and MIF and were closely associated with lung metastasis. Then, BCDEG, a reliable risk regression model that could accurately assess the metastasis risk and prognosis of patients, was established, providing a new strategy for the diagnosis and treatment of osteosarcoma.

Transcriptomic Analysis of Canine Osteosarcoma from a Precision Medicine Perspective Reveals Limitations of Differential Gene Expression Studies

Despite significant advances in cancer diagnosis and treatment, osteosarcoma (OSA), an aggressive primary bone tumor, has eluded attempts at improving patient survival for many decades. The difficulty in managing OSA lies in its extreme genetic complexity, drug resistance, and heterogeneity, making it improbable that a single-target treatment would be beneficial for the majority of affected individuals. Precision medicine seeks to fill this gap by addressing the intra- and inter-tumoral heterogeneity to improve patient outcome and survival. The characterization of differentially expressed genes (DEGs) unique to the tumor provides insight into the phenotype and can be useful for informing appropriate therapies as well as the development of novel treatments. Traditional DEG analysis combines patient data to derive statistically inferred genes that are dysregulated in the group; however, the results from this approach are not necessarily consistent across individual patients, thus contradicting the basis of precision medicine. Spontaneously occurring OSA in the dog shares remarkably similar clinical, histological, and molecular characteristics to the human disease and therefore serves as an excellent model. In this study, we use transcriptomic sequencing of RNA isolated from primary OSA tumor and patient-matched normal bone from seven dogs prior to chemotherapy to identify DEGs in the group. We then evaluate the universality of these changes in transcript levels across patients to identify DEGs at the individual level. These results can be useful for reframing our perspective of transcriptomic analysis from a precision medicine perspective by identifying variations in DEGs among individuals.

Novel Immunotherapies for Osteosarcoma

Osteosarcoma (OS) is the most common primary malignant bone sarcoma mainly affecting adolescents and young adults, which often progresses to pulmonary metastasis and leads to the death of OS patients. OS is characterized as a highly heterogeneous cancer type and the underlying pathologic mechanisms triggering tumor progress and metastasis are incompletely recognized. Surgery combined with neoadjuvant and postoperative chemotherapy has elevated 5-year survival to over 70% for patients with localized OS tumors, as opposed to only 20% of patients with recurrence and/or metastasis. Therefore, novel therapeutic strategies are needed to overcome the drawbacks of conventional treatments. Immunotherapy is gaining momentum for the treatment of OS with an increasing number of FDA-approved therapies for malignancies resistant to conventional therapies. Here, we review the OS tumor microenvironment and appraise the promising immunotherapies available in the management of OS.

N6-Methyladenosine-Related lncRNAs Are Potential Prognostic Biomarkers and Correlated With Tumor Immune Microenvironment in Osteosarcoma

N6-methyladenosine (m6A) and long non-coding RNAs (lncRNAs) play vital roles in the prognostic value and immune microenvironment of malignant tumors. Here, we constructed a m6A-related lncRNA signature in osteosarcoma samples from TCGA dataset and analyzed the association of the signature with tumor immune microenvironment. m6A-related lncRNAs were identified by performing Pearson's correlation analysis and were used to construct a novel m6A-related lncRNA signature in osteosarcoma. Validation in testing and entire cohorts confirmed the satisfactory accuracy of the risk signature. Principal-component analysis verifies the grouping ability of the risk signature. Functional enrichment analyses connected immune with the risk signature based on the six m6A-related lncRNAs. When patients were separated into high- and low-risk group based on their risk scores, we found that patients in the high-risk group had lower stromal scores, immune scores, and ESTIMATE scores, while the tumor purity was higher in the high-risk group than that in the low-risk group. As for immune cell infiltration, the proportion of monocytes was significantly higher in the low-risk group than that in the high-risk group. Of the six lncRNAs, AC004812.2 was a protective factor in osteosarcoma and low expression of AC004812.2 predicted worse overall survival. Overexpression of AC004812.2 inhibited 143B cell proliferation and increased the expression levels of IGF2BP1 and YTHDF1. In all, our m6A-related lncRNA signature was a potential prognostic biomarker and correlated with tumor immune microenvironment and immune cell infiltration, and AC004812.2 might be an important regulator of m6A modification and a promising therapeutic target in osteosarcoma.

Development of Personalized Signature Based on the Immune Landscape to Predict the Prognosis of Osteosarcoma and the Response to Immunotherapy and Targeted Therapy

As a heterogeneous and aggressive disease, osteosarcoma (OS) faces great challenges to prognosis and individualized treatment. Hence, we explore the role of immune-related genes in predicting prognosis and responsiveness to immunotherapy and targeted therapies in patients with OS based on the immunological landscape of osteosarcoma. Based on the database of the Therapeutical Applicable Research to Generate Effective Treatments (TARGET), single-sample gene set enrichment analysis (ssGSEA) was used to obtain the enrichment scores of 29 immune characteristics. A series of bioinformatics methods were performed to construct the immune-related prognostic signature (IRPS). Gene set enrichment analysis and gene set variation analysis were used to explore the biological functions of IRPS. We also analyzed the relationship between IRPS and tumor microenvironment. Lastly, the reactivity of IRPS to immune checkpoint therapy and targeted drugs was explored. The ssGSEA algorithm was used to define two immune subtypes, namely Immunity_High and Immunity_Low. Immunity_High was associated with a good prognosis and was an independent prognostic factor of OS. The IRPS containing 7 genes was constructed by the least absolute shrinkage and selection operator Cox regression. The IRPS can divide patients into low- and high-risk patients. Compared with high-risk patients, low-risk patients had a better prognosis and were positively correlated with immune cell infiltration and immune function. Low-risk patients benefited more from immunotherapy, and the sensitivity of targeted drugs in high- and low-risk groups was determined. IRPS can be used to predict the prognosis of OS patients, and provide therapeutic responsiveness to immunotherapy and targeted therapy.

Profiles of immune cell infiltration and immune-related genes in the tumor microenvironment of osteosarcoma cancer

Backgrounds

Osteosarcomas are one of the most common primary malignant tumors of bone. It primarily occurs in children and adolescents, with the second highest incidence among people over 50 years old. Although there were immense improvements in the survival of patients with osteosarcoma in the past 30 years, targetable mutations and agents of osteosarcomas still have been generally not satisfactory. Therefore, it is of great importance to further explore the highly specialized immune environment of bone, genes related to macrophage infiltration and potential therapeutic biomarkers and targets.

Methods

The 11 expression data sets of OS tissues and the 11 data sets of adjacent non-tumorous tissues available in the GEO database GSE126209 were used to conduct immune infiltration analysis. Then, through WGCNA analysis, we acquired the co-expression modules related to Mast cells activated and performed the GO and KEGG enrichment analysis. Next, we did the survival prognosis analysis and plotted a survival curve. Finally, we analyzed the COX multivariate regression of gene expression on clinical parameters and drew forest maps for visualization by the forest plot package.

Results

OS disease-related immune cell populations, mainly Mast cells activated, have higher cell content (p = 0.006) than the normal group. Then, we identified co-expression modules related to Mast cells activated. In sum, a total of 822 genes from the top three strongest positive correlation module MEbrown4, MEdarkslateblue and MEnavajowhite2 and the strongest negative correlation module MEdarkturquoise. From that, we identified nine genes with different levels in immune cell infiltration related to osteosarcoma, eight of which including SORBS2, BAIAP2L2, ATAD2, CYGB, PAMR1, PSIP1, SNAPC3 and ZDHHC21 in their low abundance have higher disease-free survival probability than the group in their high abundances.

Conclusion

These results could assist clinicians to select targets for immunotherapies and individualize treatment strategies for patients with OS.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-09042-6.

Avenues of research in dietary interventions to target tumor metabolism in osteosarcoma

Osteosarcoma (OS) is the most frequent primary bone cancer, affecting mostly children and adolescents. Although much progress has been made throughout the years towards treating primary OS, the 5-year survival rate for metastatic OS has remained at only 20% for the last 30 years. Therefore, more efficient treatments are needed. Recent studies have shown that tumor metabolism displays a unique behavior, and plays important roles in tumor growth and metastasis, making it an attractive potential target for novel therapies. While normal cells typically fuel the oxidative phosphorylation (OXPHOS) pathway with the products of glycolysis, cancer cells acquire a plastic metabolism, uncoupling these two pathways. This allows them to obtain building blocks for proliferation from glycolytic intermediates and ATP from OXPHOS. One way to target the metabolism of cancer cells is through dietary interventions. However, while some diets have shown anticancer effects against certain tumor types in preclinical studies, as of yet none have been tested to treat OS. Here we review the features of tumor metabolism, in general and about OS, and propose avenues of research in dietary intervention, discussing strategies that could potentially be effective to target OS metabolism.

Molecular profiling of osteosarcoma in children and adolescents from different age groups using a next-generation sequencing panel

Osteosarcoma (OS) is a malignant bone tumor, with a peak of incidence in the second decade of life and possibly associated with the presence of germline mutations. Besides, clinicians have pointed to a second, rarer group of patients that develops OS before 10 years old. Here we access, through next-generation sequencing (NGS) strategy, the genetic alterations present in OS and blood samples from patients diagnosed before and during the second decade of life. A custom NGS panel, designed for the main alterations described in childhood and adolescence neoplasms, named Oncomine Childhood Cancer Research Assay (OCCRA©), was used. Of all 84 OS samples investigated, 42 (50%) presented some somatic variant, with TP53, MYC, CDK4, RB1 and PDGFRA genes harboring the most observed genetic variants. MYC CNVs were more frequent in tumors from patients diagnosed before 10 years old (X²₁= 5.18, p = 0.023). Additionally, patients diagnosed during the second decade of life presented a higher percentage of somatic and germline variants. Germline variants in TP53 and RB1 were found in 5 of the 11 (45.5%) patients analyzed. Clinical variables and tumor histopathological characteristics were also collected and correlated with our molecular findings.

Ultrathin 2D Inorganic Ancient Pigment Decorated 3D-Printing Scaffold Enables Photonic Hyperthermia of Osteosarcoma in NIR-II Biowindow and Concurrently Augments Bone Regeneration

Osteosarcoma (OS) is the primary malignant bone tumor. Despite therapeutic strategies including surgery, chemotherapy, and radiotherapy have been introduced into the war of fighting OS, the 5-year survival rate for patients still remains unchangeable for decades. Besides, the critical bone defects after surgery, drug-resistance and side effects also attenuate the therapeutic effects and predict poor prognosis. Recently, photothermal therapy (PTT) has attracted extensive attention featuring minimal invasiveness and high spatial-temporal precision characteristics. Herein, an ultrathin 2D inorganic ancient pigment Egyptian blue decorated 3D-printing scaffold (CaPCu) with profound PTT efficacy at the second near-infrared (NIR-II) biowindow against OS and enhanced osteogenesis performance is successfully constructed. Importantly, this work uncovers the underlying biological mechanisms that genes associated with cell death, proliferation, and bone development are regulated by CaPCu-scaffold-based therapy. This work not only elucidates the fascinating clinical translation prospects of CaPCu-scaffold-based PTT against OS in NIR-II biowindow, but also demonstrates the potential mechanisms and offers a novel strategy to develop the next-generation, multifunctional tissue-engineering biomaterials.

Targeting Wnt/β-catenin signaling by microRNAs as a therapeutic approach in chemoresistant osteosarcoma

Osteosarcoma (OS) is an adolescent and young adult malignancy that mostly occurs in long bones. The treatment of OS is still a big challenge for clinicians due to increasing chemoresistance, and many efforts are being made today to find more beneficial treatments. In this regard, the use of microRNAs has shown a high capacity to develop promising therapies. By targeting cancer-involved signaling pathways, microRNAs reduce the cellular level of these protein pathways; thereby reducing the growth and invasion of tumors, and even leading cancer cells to apoptosis. One of these oncogenic pathways that play an important role in OS development and can be targeted by microRNAs is the Wnt/β-catenin signaling pathway. Hence, the first goal of this review article is to explain the cross-talk of microRNAs and the Wnt/β-catenin signaling in OS and then discussing recent findings of the use of microRNAs as a therapeutic approach in OS.

The survival and proliferation of osteosarcoma cells are dependent on the mitochondrial BIG3-PHB2 complex formation

Previous studies reported the critical role of the brefeldin A-inhibited guanine nucleotide exchange protein 3-prohibitin 2 (BIG3-PHB2) complex in modulating estrogen signaling activation in breast cancer cells, yet its pathophysiological roles in osteosarcoma (OS) cells remain elusive. Here, we report a novel function of BIG3-PHB2 in OS malignancy. BIG3-PHB2 complexes were localized mainly in mitochondria in OS cells, unlike in estrogen-dependent breast cancer cells. Depletion of endogenous BIG3 expression by small interfering RNA (siRNA) treatment led to significant inhibition of OS cell growth. Disruption of BIG3-PHB2 complex formation by treatment with specific peptide inhibitor also resulted in significant dose-dependent suppression of OS cell growth, migration, and invasion resulting from G2/M-phase arrest and in PARP cleavage, ultimately leading to PARP-1/apoptosis-inducing factor (AIF) pathway activation-dependent apoptosis in OS cells. Subsequent proteomic and bioinformatic pathway analyses revealed that disruption of the BIG3-PHB2 complex might lead to downregulation of inner mitochondrial membrane protein complex activity. Our findings indicate that the mitochondrial BIG3-PHB2 complex might regulate PARP-1/AIF pathway-dependent apoptosis during OS cell proliferation and progression and that disruption of this complex may be a promising therapeutic strategy for OS.

A Novel Six Metastasis-Related Prognostic Gene Signature for Patients With Osteosarcoma

Osteosarcoma is the most common malignant bone tumor, and although there has been significant progress in its management, metastases often herald incurable disease. Here we defined genes differentially expressed between primary and metastatic osteosarcoma as metastasis-related genes (MRGs) and used them to construct a novel six-MRG prognostic signature for overall survival of patients with osteosarcoma. Validation in internal and external datasets confirmed satisfactory accuracy and generalizability of the prognostic model, and a nomogram based on the signature and clinical variables was constructed to aid clinical decision-making. Of the six MRGs, FHIT is a well-documented tumor suppressor gene that is poorly defined in osteosarcoma. Consistent with tumor suppressor function, FHIT was downregulated in osteosarcoma cells and human osteosarcoma samples. FHIT overexpression inhibited osteosarcoma proliferation, migration, and invasion both in vitro and in vivo. Mechanistically, FHIT overexpression upregulate the epithelial marker E-cadherin while repressing the mesenchymal markers N-cadherin and vimentin. Our six-MRG signature represents a novel and clinically useful prognostic biomarker for patients with osteosarcoma, and FHIT might represent a therapeutic target by reversing epithelial to mesenchymal transition.

Osteosarcoma in One of Identical Twins: Three Cases Report and a Literature Review

Background

Osteosarcoma (OS) is the most common primary malignant bone tumor occurring mainly in children and young adults. OS is usually seen in sporadic cases, and it is an extremely rare phenomenon in blood relatives, particularly among identical twins.

Case Presentation

The present study reports three cases of OS occurring in only one of identical twins. The first case is a high‐grade OS in the left proximal tibia of a 16‐year‐old girl, treated with neo‐adjuvant chemotherapy, en bloc resection, and reconstruction with a modular knee tumor prosthesis. The second one is a high‐grade OS of the left proximal tibia of a 6‐year‐old girl. The patient was treated with neo‐adjuvant chemotherapy, en bloc resection, and reconstruction with inactived autograft. The third one is a conventional OS of the right proximal tibia of a 20‐year‐old woman. She was treated with neo‐adjuvant chemotherapy, en bloc resection, and reconstruction with a custom‐made prosthesis.

Conclusions

The occurrence of OS in one of identical twins is a relatively rare event but may present the best opportunity to understand the genetic mechanisms underlying the tumorigenesis and progression of this disease in humans. A longer follow‐up period and regular imaging evaluation are needed to confirm whether the identical twin of these patients will suffer OS in the future.

Clinicopathologic and genetic features of metaplastic breast cancer with osseous differentiation: a series of 6 cases

AIMS

Metaplastic breast cancer (MBC) comprises a heterogeneous group of tumors, of which MBC with osseous differentiation (MBC-OD) is extremely rare that only few cases have been reported. This study aimed to present the clinicopathologic and molecular features of MBC-OD.

METHODS

We collected 6 cases of MBC-OD from five different centers and described its clinicopathologic and molecular characteristics based on the next-generation sequencing. Another 11 cases from the literature were also reviewed to better characterize the tumor.

RESULTS

The tumor primarily showed an osteosarcoma-like appearance, which composed of high cellularity with spindle cells and osteoblast-like cells producing coarse lace-like neoplastic bone (4/6) or osteoid matrix (6/6). 55 somatic mutations including 39 missenses (70.9%), 9 frameshifts (16.4%), 3 splice sites (5.5%), 3 in-frame InDels (5.5%) and 1 nonsense (1.8%) were identified. TP53 was the most frequently mutated genes (5/6, 83.3%), followed by RB1 (3/6, 50.0%), BCOR (2/6, 33.3%), MED12 (2/6, 33.3%), PIK3CA (2/6, 33.3%), and TET2 (2/6, 33.3%). Genetic alterations suggesting therapies with clinical benefit, including mTOR inhibitors, tyrosine kinase inhibitors (TKI), and poly-ADP ribose polymerase inhibitor (PARPi), were observed in five cases. The median follow-up was 21 months (range, 3-80 months). Local recurrence was observed in two cases and three patients displayed distant metastasis. Two patients died of the disease at 3 months and 7 months, respectively.

CONCLUSIONS

Based on this series, MBC-OD is a highly aggressive breast tumor with osteosarcoma-like morphology and a high incidence of recurrent disease showing specific genetic profiles.

Efficacy and safety of apatinib in the treatment of osteosarcoma: a single-arm meta-analysis among Chinese patients

Background

Osteosarcoma is a relatively rare malignant tumor with a high incidence in young people. The development of tyrosine kinase inhibitors has brought the treatment of osteosarcoma into a new stage. Apatinib, a tyrosine kinase inhibitor specifically targeting VEGFR2, has been increasingly reported as a treatment for osteosarcoma with promising outcome parameters, but there has been no systematic analysis of the treatment of osteosarcoma by apatinib.

Methods

A single-arm meta-analysis was performed, and published literature from PubMed, Web of Science, Embase, Cochrane Library, CNKI and Wan Fang databases as of March 1, 2021 was systematically retrieved. Quality assessment is carried out in accordance with a 20 item checklist form prepared by the Institute of Health Economics (IHE). Double arcsine transformation is performed to stabilize the variance of the original ratio. When I² > 50%, the random effect model is used to calculate the pooled parameters; otherwise, the fixed effect model is used. We conducted subgroup analysis according to age and apatinib dose.

Results

This meta-analysis included 11 studies of 356 Chinese patients with osteosarcoma. The pooled objective remission rate (ORR) of advanced or metastatic osteosarcoma treated by oral apatinib in Chinese patients was 0.27(95%CI = 0.18–0.38). The pooled disease control rate (DCR) was 0.57 (95%CI = 0.42–0.72). The pooled median progression-free survival (mPFS) and median total survival (mOS) were 5.18 months (95%CI = 4.03–6.33) and 10.87 months (95% CI = 9.40–12.33), respectively. More than 70% of adverse reactions were mild, the most common adverse reaction was hand-foot syndrome (HFMD), with an incidence of 0.46 (95%CI = 0.35–0.58), the second was hypertension, with an incidence of 0.40 (95%CI = 0.29–0.51).

Conclusions

The efficacy of apatinib in the treatment of osteosarcoma is competitive with current evidence, and it is worth noting that its low cost can significantly improve patient compliance and increase therapeutic value.

The effect of miR-539 regulating TRIAP1 on the apoptosis, proliferation, migration and invasion of osteosarcoma cells

Objective

The purpose of this study is to explore the effect of miRNA-539 on osteosarcoma (OS) and the underlying mechanism, so as to find a new method for early diagnosis and treatment of osteosarcoma.

Method

miRNA-539 mimics was transfected into osteosarcoma cells 143b and MG-63 and upregulated the expression of miR-539. QT-PCR was used to detect transfection efficacy. CCK-8 method was used to detect proliferation of 143b and MG-63 osteosarcoma cells and flow cytometry was used to detect the apoptosis of osteosarcoma cells 143b and MG-63. Wound-healing test and Transwell test were used to detect the migration and invasion ability of osteosarcoma cells. TRIAP1 was found to be the potential target gene of miRNA-539 by online bioinformatics software and the expression level of TRIAP1 in osteosarcoma cells overexpressing miRNA-539 was detected by qT-PCR. Western blot was used to detect the level of expression of TRIAP1 and its downstream genes (p53, p21, apaf1 and caspase9) in osteosarcoma cells 143b and MG63 transfected with miR-539 mimics or miR-539 mimics-NC. A model of osteosarcoma subcutaneously transplanted in nude mice was constructed to observe the effect of miRNA-539 on the growth rate of osteosarcoma in vivo.

Results

After transfection of miRNA-539 mimics in osteosarcoma cells 143b and MG63, the proliferation level, migration ability, and invasion ability of the osteosarcoma cells were significantly lower than that in the control group, and the apoptosis level was significantly higher than that in the control group (P < 0.01). The dual luciferase reporter confirmed that TRIAP1 was the target of miR-539, and the expression level of TRIAP1 in 143b and MG63 transfected with miRNA-539 mimics was proved to be significantly lower than that in the control group (P < 0.01).The western blot showed the expression of genes targeted by TRIAP1 was upregulated when the expression of TRIAP1 was downregulated. In vivo, the osteosarcoma growth rate in the miRNA-539 mimics group was significantly slower than that in the control group (P < 0.01).

Conclusions

MiRNA-539 may inhibit the cell proliferation, migration and invasion of osteosarcoma cells and promote the apoptosis of osteosarcoma cells by targeting on TRIAP1.

The effects of MEX3A knockdown on proliferation, apoptosis and migration of osteosarcoma cells

Background

Osteosarcoma is an aggressive malignant tumor which has attracted worldwide attention. MEX3A may be associated with tumors while has not yet seen its coverage on osteosarcoma. Herein, this study was to investigate the correlation between MEX3A and the progression of osteosarcoma.

Methods

Firstly, we determined that expression of MEX3A was significantly higher in osteosarcoma tissues than that in marginal bone by immunohistochemical staining. Additionally, MEX3A expression was downregulated by the RNAi‐mediated knockdown. The functions of MEX3A knockdown on proliferation, apoptosis, cell cycle, migration was assessed by MTT assay, flow cytometry, wound-healing assay and Transwell assay, respectively. Knockdown of MEX3A resulted in suppressing cell proliferation, increasing cell apoptosis, inducing the G2 phase cell cycle arrest, and attenuating cellular migration. Furthermore, mouse xenograft model confirmed inhibitory effects of MEX3A knockdown on osteosarcoma formation.

Results

The preliminary exploration on the molecular mechanism of MEX3A in osteosarcoma cells showed that the induction of apoptosis needs the participation of a series of apoptosis- associated factors, such as upregulation of Caspase 3, Caspase 8 and HSP60, downregulation of HSP27 and XIAP.

Conclusions

In summary, these findings predicated that therapy directed at decreasing MEX3A expression is a potential osteosarcoma treatment.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-01882-3.

Dog breeds and body conformations with predisposition to osteosarcoma in the UK: a case-control study

Background

Osteosarcoma is an aggressive and painful bone neoplasm in dogs. Previous studies have reported epidemiological associations suggesting that large body mass, long bone length and the genetics of certain breeds including the Rottweiler are associated with elevated osteosarcoma risk. However, these studies were often limited by selection bias and confounding factors, and have rarely offered insights into breed-associated protection for osteosarcoma. The current study includes 1756 appendicular and axial osteosarcoma cases presenting to VPG Histology (Bristol, UK) compared against a control population of 905,211 dogs without osteosarcoma from primary care electronic patient records in the VetCompass™ dataset.

Methods and study design

Retrospective, case-control study. Multivariable logistic regression analysis explored associations between demographic risk factors (including breed, chondrodystrophy, age, sex/neuter status, skull-shape, and body mass) and osteosarcoma of all anatomical sites.

Results

We identified several breeds with increased and reduced odds of osteosarcoma. At highest risk were the Rottweiler and Great Dane, with > 10 times the odds of osteosarcoma compared with crossbreds, and the Rhodesian Ridgeback, which has not featured in previous lists of at-risk breeds for osteosarcoma, and had an odds ratio of 11.31 (95% confidence interval 7.37–17.35). Breeds at lowest risk of osteosarcoma (protected breeds) included the Bichon Frise, the French Bulldog and the Cavalier King Charles Spaniel, all with odd ratios of less than 0.30 compared with crossbreds. Body mass was strongly associated with osteosarcoma risk; dogs over 40 kg exhibited osteosarcoma odds of 45.44 (95% confidence interval 33.74–61.20) compared with dogs less than 10 kg. Chondrodystrophic breeds had an osteosarcoma odds ratio of 0.13 (95% confidence interval 0.11–0.16) compared with non-chondrodystrophic breeds.

Conclusions

This study provides evidence of strong breed-associated osteosarcoma risk and protection, suggesting a genetic basis for osteosarcoma pathogenesis. It highlights that breeds selected for long legs/large body mass are generally overrepresented amongst at-risk breeds, whilst those selected for short leg length/small body mass are generally protected. These findings could inform genetic studies to identify osteosarcoma risk alleles in canines and humans; as well as increasing awareness amongst veterinarians and owners, resulting in improved breeding practices and clinical management of osteosarcoma in dogs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40575-021-00100-7.

The histone acetyltransferase HBO1 functions as a novel oncogenic gene in osteosarcoma

HBO1 (KAT7 or MYST2) is a histone acetyltransferase that acetylates H3 and H4 histones.

Methods: HBO1 expression was tested in human OS tissues and cells. Genetic strategies, including shRNA, CRISPR/Cas9 and overexpression constructs, were applied to exogenously alter HBO1 expression in OS cells. The HBO1 inhibitor WM-3835 was utilized to block HBO1 activation.

Results:HBO1 mRNA and protein expression is significantly elevated in OS tissues and cells. In established (MG63/U2OS lines) and primary human OS cells, shRNA-mediated HBO1 silencing and CRISPR/Cas9-induced HBO1 knockout were able to potently inhibit cell viability, growth, proliferation, as well as cell migration and invasion. Significant increase of apoptosis was detected in HBO1-silenced/knockout OS cells. Conversely, ectopic HBO1 overexpression promoted OS cell proliferation and migration. We identified ZNF384 (zinc finger protein 384) as a potential transcription factor of HBO1. Increased binding between ZNF384 and HBO1 promoter was detected in OS cell and tissues, whereas ZNF384 silencing via shRNA downregulated HBO1 and produced significant anti-OS cell activity. In vivo, intratumoral injection of HBO1 shRNA lentivirus silenced HBO1 and inhibited OS xenograft growth in mice. Furthermore, growth of HBO1-knockout OS xenografts was significantly slower than the control xenografts. WM-3835, a novel and high-specific small molecule HBO1 inhibitor, was able to potently suppressed OS cell proliferation and migration, and led to apoptosis activation. Furthermore, intraperitoneal injection of a single dose of WM-3835 potently inhibited OS xenograft growth in SCID mice.

Conclusion: HBO1 overexpression promotes OS cell growth in vitro and in vivo.

Cyclin-dependent kinase 7 (CDK7) is an emerging prognostic biomarker and therapeutic target in osteosarcoma

Background:

Overexpression of cyclin-dependent kinase 7 (CDK7) is a well-known pathogenic feature of various malignancies and a sign of a more dismal prognosis. As relatively little is known about CDK7 in osteosarcoma, we elected to evaluate its expression, prognostic value, and function.

Methods:

We began by analyzing the publicly available data sets on CDK7 expression, including RNA sequencing data from the Therapeutically Applicable Research to Generate Effective Treatments on Osteosarcoma (TARGET-OS) and the Gene Expression database of Normal and Tumor tissues 2 (GENT2). The correlation between patient tissue CDK7 expression and their clinicopathological features and prognosis was assessed via immunohistochemical staining of a unique tissue microarray constructed from osteosarcoma specimens. Furthermore, we analyzed CDK7 expression in osteosarcoma cell lines and tissues by Western blot. CDK7-specific siRNA and a highly-selective CDK7 inhibitor, BS-181, were applied to determine the function of CDK7 on osteosarcoma cell growth and proliferation. In addition, the effect of CDK7 inhibition on clonogenicity was evaluated using a clonogenic assay, and a 3D cell culture model was used to mimic CDK7 effects in an in vivo environment.

Results:

Our results demonstrate that higher CDK7 expression significantly correlates with recurrence, metastasis, and shorter overall survival in osteosarcoma patients. Therapeutically, we show that CDK7 knockdown with siRNA or selective inhibition with BS-181 decreases proliferation and induces apoptosis of osteosarcoma cells.

Conclusion:

This study supports CDK7 overexpression as an independent predictor of poor prognosis and promising therapeutic target for osteosarcoma.

Precision Medicine in Osteosarcoma: MATCH Trial and Beyond

Osteosarcoma (OS) is a rare bone malignant tumour with a poor prognosis in the case of recurrence. So far, there is no agreement on the best systemic therapy for relapsed OS. The availability of next generation sequencing techniques has recently revolutionized clinical research. The sequencing of the tumour and its matched normal counterpart has the potential to reveal a wide landscape of genetic alterations with significant implications for clinical practice. The knowledge that the genomic profile of a patient’s tumour can be precisely mapped and matched to a targeted therapy in real time has improved the development of precision medicine trials (PMTs). PMTs aiming at determining the effectiveness of targeted therapies could be advantageous for patients with a tumour refractory to standard therapies. Development of PMTs for relapsed OS is largely encouraging and is in its initial phase. Assessing OS features, such as its rarity, its age distribution, the technical issues related to the bone tissue origin, and its complex genomic landscape, represents a real challenge for PMTs development. In this light, a multidisciplinary approach is required to fully exploit the potential of precision medicine for OS patients.

Identification of Immune-Related Genes MSR1 and TLR7 in Relation to Macrophage and Type-2 T-Helper Cells in Osteosarcoma Tumor Micro-Environments as Anti-metastasis Signatures

Metastasis of osteosarcoma (OS) is an essential factor affecting the prognosis and survival of patients. The tumor microenvironment, including tumor immune-infiltrating cells (TIIC), is closely related to tumor progression. The purpose of this study was to investigate the differences between metastatic and non-metastatic immune-infiltrating cells in OS and to identify key immune-related genes. The differences in immune infiltration in OS metastasis were calculated based on the ssGSEA algorithm of 28 immuno-infiltrating cells. Weighted gene co-expression network analysis (WGCNA) and intersection analysis were used to screen immune-related modules and hubgenes. Univariate/multivariate/Lasso Cox regressions were used for models construction and signatures screening. The receiver operating characteristic (ROC) and Kaplan-Meier (K-M) curves were constructed to observe the metastases of different groups. Both internal and external data were verified. We found that macrophages and Type-2 T-helper cells were significantly decreased in patients with OS metastases. The high-risk groups obtained from multivariate/Lasso Cox models constructed with 11 immune-related hubgenes almost all underwent distant metastases within 5 years. Interestingly and importantly, two genes, MSR1 and TLR7, appeared in various models and various hubgenes, which play an anti-metastasis role and may prolong overall survival in OS. Our study may help elucidate the impact of TIIC on OS metastasis outcomes and to identify biomarkers and therapeutic targets.

Single-cell RNA landscape of intratumoral heterogeneity and immunosuppressive microenvironment in advanced osteosarcoma

Osteosarcoma is the most frequent primary bone tumor with poor prognosis. Through RNA-sequencing of 100,987 individual cells from 7 primary, 2 recurrent, and 2 lung metastatic osteosarcoma lesions, 11 major cell clusters are identified based on unbiased clustering of gene expression profiles and canonical markers. The transcriptomic properties, regulators and dynamics of osteosarcoma malignant cells together with their tumor microenvironment particularly stromal and immune cells are characterized. The transdifferentiation of malignant osteoblastic cells from malignant chondroblastic cells is revealed by analyses of inferred copy-number variation and trajectory. A proinflammatory FABP4+ macrophages infiltration is noticed in lung metastatic osteosarcoma lesions. Lower osteoclasts infiltration is observed in chondroblastic, recurrent and lung metastatic osteosarcoma lesions compared to primary osteoblastic osteosarcoma lesions. Importantly, TIGIT blockade enhances the cytotoxicity effects of the primary CD3+ T cells with high proportion of TIGIT+ cells against osteosarcoma. These results present a single-cell atlas, explore intratumor heterogeneity, and provide potential therapeutic targets for osteosarcoma.

Recent advances in understanding osteosarcoma and emerging therapies

Osteosarcoma is the most common bone cancer in adolescents and young adults, but it is a rare cancer with no improvement in patient survival in the last four decades. The main problem of this bone tumor is its evolution toward lung metastatic disease, despite the current treatment strategy (chemotherapy and surgery). To further improve survival, there is a strong need for new therapies that control osteosarcoma cells with metastatic potential and their favoring tumor microenvironment (ME) from the diagnosis. However, the complexity and heterogeneity of those tumor cell genomic/epigenetic and biology, the diversity of tumor ME where it develops, the sparsity of appropriate preclinical models, and the heterogeneity of therapeutic trials have rendered the task difficult. No tumor- or ME-targeted drugs are routinely available in front-line treatment. This article presents up-to-date information from preclinical and clinical studies that were recently published or presented in recent meetings which we hope might help change the osteosarcoma treatment landscape and patient survival in the near future.

Expression and clinical implications of leucine-rich repeat containing 15 (LRRC15) in osteosarcoma

Leucine-rich repeat containing 15 (LRRC15) is a member of the leucine-rich repeat superfamily that is overexpressed in various cancers and associated with higher tumor grade and aggression. Despite its known tumorigenicity, its roles within osteosarcoma are unknown, prompting us to evaluate its expression and clinical significance within this rare yet aggressive cancer. Western blots showed differential expression of LRRC15 in the osteosarcoma cell lines MNNG/HOS, KHOS, 143B, MG63, Saos-2, and U2OS. We additionally validated this positive expression, as well as sublocalization to the cell membrane, with immunofluorescence. A tissue microarray constructed from 69 osteosarcoma patient tissues was immunohistochemically stained for LRRC15 expression, stratified, and used for clinicopathological analysis. Publicly available databases on LRRC15 expression, including RNA sequencing data from the Therapeutically Applicable Research to Generate Effective Treatments on Osteosarcoma (TARGET-OS) and the Gene Expression database of Normal and Tumor tissues 2 (GENT2) were also analyzed. We found 63 of the 69 (91.3%) patient tissues exhibited some degree of LRRC15 immunostaining, including no staining (6 of 69, 8.7%), 1+ staining (12 of 69, 17.4%), 2+ staining (25 of 69, 36.2%), and 3+ staining (26 of 69, 37.7%). The patients with osteosarcomas having elevated LRRC15 expression demonstrated comparatively increased metastasis, chemoresistance, and shorter 5-year survival rates. Our analysis of the TARGET-OS and GENT2 databases also showed increased LRRC15 gene expression in osteosarcoma. Taken together, our study supports LRRC15 as a prognostic biomarker and emerging therapeutic target in osteosarcoma.

p38γ overexpression promotes osteosarcoma cell progression

Osteosarcoma (OS) is the most common primary bone malignancy in the adolescent population. Recent studies demonstrate that p38 gamma (p38γ) phosphorylates retinoblastoma (Rb) to promote cyclin expression, cell-cycle entry and tumorigenesis. Studying the potential function of p38γ in human OS, we show that p38γ mRNA and protein expression are significantly elevated in OS tissues and OS cells, whereas its expression is relatively low in normal bone tissue and in human osteoblasts/osteoblastic cells. Knockdown of p38γ in established (U2OS) and primary human OS cells potently inhibited cell growth, proliferation, migration and invasion, while promoting cell apoptosis. Furthermore, CRISPR/Cas9-induced p38γ knockout inhibited human OS cell progression in vitro. Conversely, ectopic overexpression of p38γ in primary human OS cells augmented cell growth, proliferation and migration. Signaling studies show that retinoblastoma (Rb) phosphorylation and cyclin E1/cyclin A expression were decreased following p38γ shRNA knockdown and knockout, but increased after ectopic p38γ overexpression. Collectively, these results show that p38γ overexpression promotes human OS cell progression.