Abstract
Osteosarcoma is a highly malignant bone cancer which primarily occurs in children and young adults. Increasing evidence indicates that long noncoding RNAs (lncRNAs), which function as competing endogenous RNAs (ceRNAs) that sponge microRNAs (miRNAs) and messenger RNAs (mRNAs), play a pivotal role in the pathogenesis and progression of cancers. The regulatory mechanisms of lncRNA-mediated ceRNAs in osteosarcoma have not been fully elucidated. In this study, we identified differentially expressed lncRNAs, miRNAs, and mRNAs in osteosarcoma based on RNA microarray profiles in the Gene Expression Omnibus (GEO) database. A ceRNA network was constructed utilizing bioinformatic tools. Kaplan-Meier survival analysis showed that lncR-C3orf35 and HMGB1 were associated with poor prognosis of osteosarcoma patients. Furthermore, results of Gene Set Enrichment Analysis (GSEA) suggested that lncR-C3orf35 may be involved in cellular invasion, the Toll-like receptor signaling pathway, and immune cell infiltration in the tumor microenvironment. Further analysis showed that patients with osteosarcoma metastasis expressed higher levels of lncR-C3orf35 and HMGB1 compared to metastasis-free patients. Moreover, the metastasis-free survival rate of the high lncR-C3orf35/HMGB1 expression group was significantly lower than that of the low expression group. The ESTIMATE algorithm was used to calculate the immune score and stromal scores for each sample. High lncR-C3orf35 and HMGB1 levels were correlated with low immune scores. ImmuCellAI analysis revealed that a low proportion of macrophage infiltration was associated with high lncR-C3orf35 and HMGB1 expression. The differential expression of lncR-C3orf35, miR-142-3p, and HMGB1 was further verified by quantitative real-time PCR. This study indicates that lncR-C3orf35 could be considered as a novel potential biomarker and therapeutic target of osteosarcoma, which may contribute to a better understanding of ceRNA regulatory mechanisms.