Abstract

The alteration of glycometabolism is a characteristic of cancer cells. Long non-coding RNAs (lncRNAs) have been documented to occupy a considerable position in glycometabolism regulation. This research aims to construct an effective prediction model for the prognosis of bladder cancer (BC) based on glycometabolism-associated lncRNAs (glyco-lncRNAs). Pearson correlation analysis was applied to get glyco-lncRNAs, and then, univariate cox regression analysis was employed to further filtrate survival time-associated glyco-lncRNAs. Multivariate cox regression analysis was utilized to construct the prediction model to divide bladder cancer (BC) patients into high- and low-risk groups. The overall survival (OS) rates of these two groups were analyzed using the Kaplan–Meier method. Next, gene set enrichment analysis and Cibersortx were used to explore the enrichment and the difference in immune cell infiltration, respectively. pRRophetic algorithm was applied to explore the relation between chemotherapy sensitivity and the prediction model. Furthermore, reverse transcriptase quantitative polymerase chain reaction was adopted to detect the lncRNAs constituting the prediction signature in tissues and urine exosomal samples of BC patients. A powerful model including 6 glyco-lncRNAs was proposed, capable of suggesting a risk score for each BC patient to predict prognosis. Patients with high-risk scores demonstrated a shorter survival time both in the training cohort and testing cohort, and the risk score could predict the prognosis without depending on the traditional clinical traits. The area under the receiver operating characteristic curve of the risk score was higher than that of other clinical traits (0.755 > 0.640, 0.485, 0.644, or 0.568). The high- and low-risk groups demonstrated very distinct immune cells infiltration conditions and gene set enriched terms. Besides, the