Abstract:Introduction: Osteosarcoma is a common type of bone cancer characterized by a poor prognosis due to its metastatic nature. The tumor microenvironment (TME) plays a critical role in tumor metastasis and therapy response. Therefore, our study aims to explore the metastatic mechanism of osteosarcoma, potentially opening new avenues for cancer treatment. Methods: In this study, we collected data from the GSE152048, GSE14359, and GSE49003 datasets. Differentially expressed genes (DEGs) were identified in osteosarcoma cases with primary and metastatic features using R software and the limma package. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed to investigate metastasis-related genes. A protein-protein interaction (PPI) network was established using the STRING database to further analyze these metastasis-associated genes. The abundances of different cell types with a mixed cell population were estimated using the CIBERSORT approach. The scRNA-seq data were analyzed by the Seurat package in R software, and intercellular communications were elucidated using the CellChat R package. Results: In this study, 92 DEGs related to metastasis were identified, including 41 upregulated and 51 downregulated genes in both the GSE14359 and GSE49003 datasets. Metastasis-associated pathways were identified, including those involving the cyclin-dependent protein kinase holoenzyme complex, transferase complex, transferring phosphorus-containing groups, SCF ubiquitin ligase complex, and the serine/threonine protein kinase complex. KEGG and PPI network analyses revealed 15 hub genes, including Skp2, KIF20A, CCNF, TROAP, PHB, CKS1B, MCM3, CCNA2, TRIP13, CENPM, Hsp90AB1, JUN, CKS2, TK1, and KIF4A. Skp2 has been known as an E3 ubiquitin ligase involved in osteosarcoma progression. The proportion of CD8+ T cells was found to be higher in metastatic osteosarcoma tissues, and high expression of PHB was associated with a favorable prognosis in osteosarcoma patients. Additionally, 23 cell clusters were classified into eight cell types, including chondrocytes, MSC, T cells, monocytes, tissue stem cells, neurons, endothelial cells, and macrophages. The 15 hub genes were expressed across various cell types, and interactions between different cell types were observed. Conclusion: Our study reveals the intricate communication between tumor microenvironment components and tumor metastasis in osteosarcoma.

Combined Analysis of ChIP-seq and Gene Microarray Datasets Identify the E2-mediated Genes in ERα-dependent Manner in Osteosarcoma.

Estrogen receptor β exhibited anti-tumor effects on osteosarcoma cells by regulating integrin, IAP, NF-kB/BCL-2 and PI3K/Akt signal pathway.

Estrogen Receptor Β Induces Autophagy of Osteosarcoma Through the Mtor Signaling Pathway.

Role of Estrogen Receptor-Mediated Anti-Tumor Effects in U2OS Cells.

RNA Sequencing Identifies Gene Expression Profile Changes Associated with Β-Estradiol Treatment in U2OS Osteosarcoma Cells.

Exosomes Containing Differential Expression of Microrna and Mrna in Osteosarcoma That Can Predict Response to Chemotherapy.

Identification and Analysis of Key Genes in Osteosarcoma Using Bioinformatics

Exploring the Molecular Mechanisms of Osteosarcoma by the Integrated Analysis of Mrnas and Mirna Microarrays.

Ott_a_195675 4153..4165

Identification of CD20, ECM, and ITGA as Biomarkers for Osteosarcoma by Integrating Transcriptome Analysis

Identification of Tumorigenicity-Associated Genes in Osteosarcoma Cell Lines Based on Bioinformatic Analysis and Experimental Validation.

Identification of Key Genes and Pathways in Ewing's Sarcoma Patients Associated with Metastasis and Poor Prognosis

Single-cell RNA sequencing reveals the communications between tumor microenvironment components and tumor metastasis in osteosarcoma

Predicting Novel Genes and Pathways Associated with Osteosarcoma by Using Bioinformatics Analysis.

Bioinformatics Analysis of Differentially Expressed Pathways Related to the Metastatic Characteristics of Osteosarcoma.

Identification of Potential Therapeutic Targets and Immune Cell Infiltration Characteristics in Osteosarcoma Using Bioinformatics Strategy.

Identification of Key Genes and Pathways in Osteosarcoma by Bioinformatics Analysis

Essential genes analysis reveals small ribosomal subunit protein eS28 may be a prognostic factor and potential vulnerability in osteosarcoma

Integrated Multiomic Analysis and High-Throughput Screening Reveal Potential Gene Targets and Synergetic Drug Combinations for Osteosarcoma Therapy

Construction of an ER stress-related prognostic signature for predicting prognosis and screening the effective anti-tumor drug in osteosarcoma

Comprehensive investigation of tumor immune microenvironment and prognostic biomarkers in osteosarcoma through integrated bulk and single-cell transcriptomic analysis