Abstract
Osteosarcoma is a malignant bone cancer that originates from the bone with the strongest invasiveness. Tumor formation strongly correlates with immune cell infiltration into the tumor immune microenvironment (TIME). Therefore, we aimed to identify TIME-related biomarkers as potential prognostic markers of osteosarcoma. The mRNA and long noncoding RNA (lncRNA) transcriptome data of 88 patients with osteosarcoma and the expression profile of GSE99671 were downloaded from The Cancer Genome Atlas and Gene Expression Omnibus, respectively. Immune infiltration scores and types were evaluated using ESTIMATE and CIBERSORT. A linear model was established to identify the differentially expressed genes (DEGs) and lncRNAs (DElncRNAs). Functional enrichment analysis of DEGs was conducted by Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, gene set enrichment analysis, and gene set variation analysis. DElncRNAs were analyzed using a weighted gene co-expression network. Least absolute shrinkage and selection operator regression was applied to screen for prognostic markers. Patient survival was predicted by the risk score and analyzed by receiver operating characteristic curve. Clinical features affecting patient survival were assessed. Immune infiltration positively correlated with osteosarcoma patient survival. Different immune cell infiltrates in patients with osteosarcma may serve as prognostic indicators and targets for immunotherapy. In total, 1125 DEGs, 80 DElncRNAs, and 11 pairs of co-expressed lncRNA-mRNAs were identified. DEGs in the three modules were associated with immune infiltration into the TIME. Four DElncRNAs, namely AC015819.1, AC015911.3, AL365361.1, and USP30-AS1, showed good prognostic ability for osteosarcoma and were positively correlated with the immune score. Tumor metastasis and risk scores alone were good prognostic indicators, and a combination of the two variables can better predict the prognosis of osteosarcoma. We identified four lncRNAs, AC015819.1, AC015911.3, AL365361.1, and USP30-AS1, as potential biomarkers for osteosarcoma prognosis.