Combining Single-Cell and Transcriptomic Data Revealed the Prognostic Significance of Glycolysis in Pancreatic Cancer

Background: Pancreatic cancer (PC), the most common fatal solid malignancy, has a very dismal prognosis. Clinical computerized tomography (CT) and pathological TNM staging are no longer sufficient for determining a patient’s prognosis. Although numerous studies have suggested that glycolysis is important in the onset and progression of cancer, there are few publications on its impact on PC. Methods: To begin, the single-sample gene set enrichment analysis (ssGSEA) approach was used to quantify the glycolysis pathway enrichment fraction in PC patients and establish its prognostic significance. The genes most related to the glycolytic pathway were then identified using weighted gene co-expression network analysis (WGCNA). The glycolysis-associated prognostic signature in PC patients was then constructed using univariate Cox regression and lasso regression methods, which were validated in numerous external validation cohorts. Furthermore, we investigated the activation of the glycolysis pathway in PC cell subtypes at the single-cell level, performed a quasi-time series analysis on the activated cell subtypes and then detected gene changes in the signature during cell development. Finally, we constructed a decision tree and a nomogram that could divide the patients into different risk subtypes, according to the signature score and their different clinical characteristics and assessed the prognosis of PC patients. Results: Glycolysis plays a risky role in PC patients. Our glycolysis-related signature could effectively discriminate the high-risk and low-risk patients in both the trained cohort and the independent externally validated cohort. The survival analysis and multivariate Cox analysis indicated this gene signature to be an independent prognostic factor in PC. The prognostic ROC curve analysis suggested a high accuracy of this gene signature in predicting the patient prognosis in PC. The single-cell analysis suggested that the glycolytic pathway may be more activated in epithelial cells and that the genes in the signature were also mainly expressed in epithelial cells. The decision tree analysis could effectively identify patients in different risk subgroups, and the nomograms clearly show the prognostic assessment of PC patients. Conclusion: Our study developed a glycolysis-related signature, which contributes to the risk subtype assessment of patients with PC and to the individualized management of patients in the clinical setting.