Identifications of genetic differences between metastatic and non-metastatic osteosarcoma samples based on bioinformatics analysis

Identification of potential crucial genes and key pathways in osteosarcoma

Background

The aim of this study is to identify the potential pathogenic and metastasis-related differentially expressed genes (DEGs) in osteosarcoma through bioinformatic analysis of Gene Expression Omnibus (GEO) database.

Results

Gene expression profiles of GSE14359, GSE16088, and GSE33383, in total 112 osteosarcoma tissue samples and 7 osteoblasts, were analyzed. Seventy-four normal-primary DEGs (NPDEGs) and 764 primary-metastatic DEGs (PMDEGs) were screened. VAMP8, A2M, HLA-DRA, SPARCL1, HLA-DQA1, APOC1 and AQP1 were identified continuously upregulating during the oncogenesis and metastasis of osteosarcoma. The enriched functions and pathways of NPDEGs include procession and presentation of antigens, activation of MHC class II receptors and phagocytosis. The enriched functions and pathways of PMDEGs include mitotic nuclear division, cell adhesion molecules (CAMs) and focal adhesion. With protein-protein interaction (PPI) network analyzed by Molecular Complex Detection (MCODE) plug-in of Cytoscape software, one hub NPDEG (HLA-DRA) and 7 hub PMDEGs (CDK1, CDK20, CCNB1, MTIF2, MRPS7, VEGFA and EGF) were eventually selected, and the most significant pathways in NPDEGs module and PMDEGs module were enriched in the procession and presentation of exogenous peptide antigen via MHC class II and the nuclear division, respectively.

Conclusions

By integrated bioinformatic analysis, numerous DEGs related to osteosarcoma were screened, and the hub DEGs identified in this study are possibly part of the potential biomarkers for osteosarcoma. However, further experimental studies are still necessary to elucidate the biological function and mechanism of these genes.

miR‑193b exhibits mutual interaction with MYC, and suppresses growth and metastasis of osteosarcoma

Emerging evidence has indicated that microRNAs (miRs) are involved in the malignant behavior of cancer. The present study explored the role of miR‑193b in the development and metastasis of osteosarcoma. Compared with F4 osteosarcoma cells, which have a relatively low metastatic potential, highly metastatic F5M2 cells exhibited a lower expression of miR‑193b. Furthermore, miR‑193b exerted negative effects on cell proliferation, colony formation, cell cycle progression, migration and invasion, and induced apoptosis. In vivo studies revealed negative influences of miR‑193b on tumorigenesis and metastasis. The tumor‑suppressive role of miR‑193b was achieved by targeting KRAS and stathmin 1 (STMN1). Notably, overexpression of KRAS and STMN1 attenuated the miR‑193b‑induced inhibition of malignant behaviors. There was a double‑negative regulatory loop between MYC and miR‑193b, with MYC inhibiting miR‑193b expression by directly binding to its promoter region and miR‑193b negatively influencing MYC expression indirectly through some unknown mechanism. Collectively, these findings indicated that miR‑193b may serve a tumor suppressive role in osteosarcoma by targeting KRAS and STMN1. The double‑negative regulatory loop between MYC and miR‑193b may contribute to the sustained upregulation of MYC, the downregulation of miR‑193b, and to the subsequently enhanced expression of KRAS and STMN1, which may eventually lead to the development and metastasis of osteosarcoma.