Construction and validation of a glycolysis-related lncRNA signature for prognosis prediction in Stomach Adenocarcinoma

Background: Glycolysis is closely related to the occurrence and progression of gastric cancer (GC). Currently, there is no systematic study on using the glycolysis-related long non-coding RNA (lncRNA) as a model for predicting the survival time in patients with GC. Therefore, it was essential to develop a signature for predicting the survival based on glycolysis-related lncRNA in patients with GC. Materials and methods: LncRNA expression profiles, containing 375 stomach adenocarcinoma (STAD) samples, were obtained from The Cancer Genome Atlas (TCGA) database. The co-expression network of lncRNA and glycolysis-related genes was used to identify the glycolysis-related lncRNAs. The Kaplan-Meier survival analysis and univariate Cox regression analysis were used to detect the glycolysis-related lncRNA with prognostic significance. Then, Bayesian Lasso-logistic and multivariate Cox regression analyses were performed to screen the glycolysis-related lncRNA with independent prognostic significance and to develop the risk model. Patients were assigned into the low- and high-risk cohorts according to their risk scores. A nomogram model was constructed based on clinical information and risk scores. Gene Set Enrichment Analysis (GSEA) was performed to visualize the functional and pathway enrichment analyses of the glycolysis-related lncRNA. Finally, the robustness of the results obtained was verified in an internal validation data set. Results: Seven glycolysis-related lncRNAs (AL353804.1, AC010719.1, TNFRSF10A-AS1, AC005586.1, AL355574.1, AC009948.1, and AL161785.1) were obtained to construct a risk model for prognosis prediction in the STAD patients using Lasso regression and multivariate Cox regression analyses. The risk score was identified as an independent prognostic factor for the patients with STAD [HR = 1.315, 95% CI (1.056–1.130), p < 0.001] via multivariate Cox regression analysis. Receiver operating characteristic (ROC) curves were drawn and the area under curve (AUC) values of 1-, 3-, and 5-year overall survival (OS) were calculated to be 0.691, 0.717, and 0.723 respectively. Similar results were obtained in the validation data set. In addition, seven glycolysis-related lncRNAs were significantly enriched in the classical tumor processes and pathways including cell adhesion, positive regulation of vascular endothelial growth factor, leukocyte transendothelial migration, and JAK_STAT signaling pathway. Conclusion: The prognostic prediction model constructed using seven glycolysis-related lncRNA could be used to predict the prognosis in patients with STAD, which might help clinicians in the clinical treatment for STAD.