Identification of Epigenetic-Dysregulated lncRNAs Signature in Osteosarcoma by Multi-Omics Data Analysis

Alterations of epigenetic modification patterns are potential markers of cancer. The current study characterized six histone modifications in osteosarcoma and identified epigenetically dysregulated long non-coding RNAs (epi-lncRNAs).

Multi-omics data were obtained from osteosarcoma cell line SJSA1 and a normal cell line. Differentially expressed lncRNAs (DElncRNAs) between osteosarcoma and normal skeletal muscle were analyzed using Limma. MACS2 was applied to identify the "peaks" modified by each histone in the cell. Promoters or enhancers of DElncRNA were overlapped with differential histone-modified regions (DHMR) to screen epi-lncRNAs. Univariate and multivariate Cox regression analysis were performed to detect the genes closely related to the prognosis of osteosarcoma and to construct risk models.

A total of 17 symbolic epi-lncRNA in osteosarcoma were screened, and 13 of them were differentially expressed between osteosarcoma and normal samples. Eight epi-lncRNAs were retained by Univariate Cox regression analysis. Four of these epi-lncRNAs were used to construct an epi-lncRNA signature. The risk score of each osteosarcoma sample in the high- or low-risk group was estimated according to the epi-lncRNA signature. The overall survival (OS) of the low-risk group was significantly better than that of the high-risk group. The area under the receiver operating characteristic (ROC) curve of the model was 0.79 and 0.82 for 1-, 3-, and 5-year OS, respectively.

Our results revealed the histone modification pattern in osteosarcoma and developed 4-epi-lncRNA signature to predict the prognosis of osteosarcoma, laying a foundation for the identification of highly specific epigenetic biomarkers for osteosarcoma.