Identification of Biomarkers and Study of Mechanisms Related to Metastatic of Osteosarcoma Based on Integrated Bioinformatics Analyses
Abstract BackgroundOsteosarcoma (OS) is a serious threat to public health. Because of high morbidity and fairly complicated pathogenesis. The study aim to identify candidate biomarkers and research the molecular mechanisms correlated of patients with metastatic OS. MethodsThe GSE21257 was downloaded from Gene Expression Omnibus(GEO) database, and the differentially expressed RNAs (DERs) were identified and functional enriched analysis by statistical soft-ware in R. Subsequently, the co-expression modules and its clinical characteristics of OS were identified by weighted gene co-expression network analysis (WGCNA) Following, the KEGG pathways directly related to metastatic OS was to researched by the Comparative Toxicogenomics Database 2019 update (CTD). Finally, the “survival” package in R was used to survival analysis and the DERs were verified using another independent profiling GSE14827. ResultsA total of 1,464 DERs were classified including 702 up-regulated and 762 down-regulated. In addition, a total of 1248 DERs were obtained by WGCNA analysis, the blue modules is the highest negative correlation (P=0) and the turquoise modules is highest positive correlation (P=3E-196) among all correlations with OS metastatic. The lncRNA-mRNA co-expression network including 4 lncRNAs and 507 mRNAs, and the cytokine-cytokine receptor interaction and JAK-STAT signaling pathway were found significantly correlation with metastatic. Finally, the increased expression levels of IFNGR1, lower DLEU1 and DLEU2 related to better prognosis. Which were significantly consistent in the another independent profiling GSE14827. ConclusionsA bioinformatics analysis related to the IFNGR1, DLEU1 and DLEU2 may as candidate biomarkers for metastatic OS.