Background: The role of N6-methyladenosine (m6A)-related long non-coding RNAs (lncRNAs) in osteosarcoma (OS) has not been fully studied yet. We aimed to identify m6A-related lncRNAs that could act as prognostic biomarkers for OS. Methods: Pearson correlation was performed to identify m6A-related lncRNAs. Univariate and multivariate Cox regression analyses were performed to construct the risk model and assess whether the risk score was an independent prognostic factor for patients with OS. Gene Set Enrichment Analysis (GSEA) was performed to analyze the functions of genes in high-risk and low-risk groups. StarBase and Cytoscape were used to construct a competing endogenous RNA (ceRNA) network based on m6A-related prognostic lncRNA signature. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed to analyze the function of genes involved in the ceRNA network. Results: We extracted 122 common lncRNAs from TCGA and Gene Expression Omnibus (GEO) databases. Pearson correlation results revealed 59 significant m6A-related lncRNAs in The Cancer Genome Atlas (TCGA) database, from which 2 were screened to construct a risk signature in TCGA dataset, which was then validated in the GEO dataset. A corresponding risk score was calculated and shown to be an independent prognostic factor for patients with OS. Enrichment analysis indicated that cell proliferation-related biological processes were more common in the high-risk group, while immune-related biological processes were more common in the low-risk group. Moreover, we established a nomogram that had a good ability to predict the overall survival of patients with OS. Additionally, a ceRNA network based on small nucleolar RNA host gene 7 (SNHG7) and small nucleolar RNA host gene 12 (SNHG12) was constructed, with genes that were enriched in hepatocellular carcinoma, gastric cancer, and non-small-cell lung cancer pathways. Conclusion: Our study revealed the prognostic role of m6A-related lncRNAs in OS and identified SNHG7 and SNHG12 as potential biomarkers for predicting the prognosis of patients with OS. These findings have enriched our understanding of the role of m6A modification in the dysregulation of lncRNAs in OS.
Keywords: competing endogenous RNA; long non-coding RNA; m6A modification; osteosarcoma; prognosis.