Construction of the prognostic enhancer RNA regulatory network in osteosarcoma

Osteosarcoma (OS) is a common malignant tumor in osteoarticular system, the 5-year overall survival of which is poor. Enhancer RNAs (eRNAs) have been implicated in the tumorigenesis of various cancer types, whereas their roles in OS tumorigenesis remains largely unclear.

Differentially expressed eRNAs (DEEs), transcription factors (DETFs), target genes (DETGs) were identified using limma (Linear Models for Microarray Analysis) package. Prognosis-related DEEs were accessed by univariate Cox regression analysis. A multivariate model was constructed to evaluate the prognosis of OS samples. Prognosis-related DEEs, DETFs, DETGs, immune cells, and hallmark gene sets were co-analyzed to construct an regulatory network. Specific inhibitors were also filtered by connectivity Map analysis. External validation and scRNA-seq analysis were performed to verify our key findings.

3,981 DETGs, 468 DEEs, 51 DETFs, and 27 differentially expressed hallmark gene sets were identified. A total of Multivariate risk predicting model based on 18 prognosis-related DEEs showed a high accuracy (area under curve (AUC) = 0.896). GW-8510 was the candidate inhibitor targeting prognosis-related DEEs (mean = 0.670, p < 0.001). Based on the OS tumorigenesis-related regulation network, we identified that CCAAT enhancer binding protein alpha (CEBPA, DETF) may regulate CD8A molecule (CD8A, DEE), thereby promoting the transcription of CD3E molecule (CD3E, DETG), which may affect allograft rejection based on CD8+ T cells.

We constructed an eRNA-based prognostic model for predicting the OS patients' prognosis and explored the potential regulation network for OS tumorigenesis by an integrated bioinformatics analysis, providing promising therapeutic targets for OS patients.