No abstract

Osteosarcoma is an aggressive malignant bone sarcoma worldwide. A causal gene network with specific functions underlying both the development and progression of OS was still unclear. Here we firstly identified the differentially expressed genes (DEGs) between control and OS samples, and then defined the hub genes and top clusters in the protein–protein interaction (PPI) network of these DEGs. By focusing on the hub gene TYROBP in the top 1 cluster, a conserved TYROBP co-expression network was identified. Then the effect of the network on OS overall survival was analyzed. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses and Gene Set Enrichment Analysis (GSEA) were used to explore the functions of the network. XCell platform and ssGSEA algorithm were conducted to estimate the status of immune infiltration. ChEA3 platform, GSEA enrichment analysis, and Drug Pair Seeker (DPS) were used to predict the key transcription factor and its upstream signal. We identified the downregulated SPI1-TYROBP-FCER1G network in OS, which were significantly enriched in immune-related functions. We also defined a two-gene signature (SPI1/FCER1G) that can predict poorer OS overall survival and the attenuated immune infiltration when downregulated. The SPI1-TYROBP-FCER1G network were potentially initiated by transcription factor SPI1 and would lead to the upregulated CD86, MHC-II, CCL4/CXCL10/CX3CL1 and hence increased immune infiltrations. With this study, we could better explore the mechanism of OS oncogenesis and metastasis for developing new therapies.