Glycolysis-Related Gene Expression Profiling Screen for Prognostic Risk Signature of Pancreatic Ductal Adenocarcinoma

Abstract

Objective: Pancreatic ductal adenocarcinoma (PDAC) is highly lethal. Although progress has been made in the treatment of PDAC, its prognosis remains unsatisfactory. This study aimed to develop novel prognostic genes related to glycolysis in PDAC and to apply these genes to new risk stratification.Methods: In this study, based on the Cancer Genome Atlas (TCGA) PAAD cohort, the expression level of glycolysis-related gene at mRNA level in PAAD and its relationship with prognosis were analyzed. Non-negative matrix decomposition (NMF) clustering was used to cluster PDAC patients according to glycolytic genes. Prognostic glycolytic genes, screened by univariate Cox analysis and LASSO regression analysis were established to calculate risk scores. The differentially expressed genes (DEGs) in the high-risk group and the low-risk group were analyzed, and the signal pathway was further enriched to analyze the correlation between glycolysis genes. In addition, based on RNA-seq data, CIBERSORT was used to evaluate the infiltration degree of immune cells in PDAC samples, and ESTIMATE was used to calculate the immune score of the samples.Results: A total of 319 glycolysis-related genes were retrieved, and all PDAC samples were divided into two clusters by NMF cluster analysis. Survival analysis showed that PDAC patients in cluster 1 had shorter survival time and worse prognosis compared with cluster 2 samples (P < 0.001). A risk prediction model based on 11 glycolysis genes was constructed, according to which patients were divided into two groups, with significantly poorer prognosis in high-risk group than in low-risk group (P < 0.001). Both internal validation and external dataset validation demonstrate good predictive ability of the model (AUC = 0.805, P < 0.001; AUC = 0.763, P < 0.001). Gene aggregation analysis showed that DEGs highly expressed in high-risk group were mainly concentrated in the glycolysis level, immune status, and tumor cell proliferation, etc. In addition, the samples in high-risk group showed immunosuppressed status and infiltrated by relatively more macrophages and less CD8+T cell.Conclusions: These findings suggested that the gene signature based on glycolysis-related genes had potential diagnostic, therapeutic, and prognostic value for PDAC.