Glycolysis-Related LINC02432/Hsa-miR-98–5p/HK2 Axis Inhibits Ferroptosis and Predicts Immune Infiltration, Tumor Mutation Burden, and Drug Sensitivity in Pancreatic Adenocarcinoma

Abstract

Pancreatic adenocarcinoma (PAAD) is a malignant cancer with high incidence and mortality. Glycometabolic rearrangements (aerobic glycolysis) is a hallmark of PAAD and contributes to tumorigenesis and progression through numerous mechanisms. This study aimed to identify a novel glycolysis-related lncRNA-miRNA-mRNA ceRNA signature in PAAD and explore its potential molecular function. We first calculated the glycolysis score for each PAAD patient by the ssGSEA algorithm and found that patients with higher hallmark glycolysis scores had poorer prognosis. Subsequently, we obtained a novel glycolysis-related LINC02432/hsa-miR-98–5p/HK2 axis from the TCGA and GEO databases using comprehensive bioinformatics analysis and developed a nomogram to predict overall survival. Furthermore, functional characterization analysis revealed that LINC02432/hsa-miR-98–5p/HK2 axis risk score was negatively correlated with ferroptosis. The tumor immune infiltration analysis suggested positive correlations between ceRNA risk score and infiltrated M0 macrophage levels in PAAD. Correlation analysis found that ceRNA risk scores were positively correlated with four chemokines (CXCL3, CXCL5, CXCL8 and CCL20) and one immune checkpoint gene (SIGLEC15). Meanwhile, tumor mutation burden (TMB), an indicator for predicting response to immunotherapy, was positively correlated with ceRNA risk score. Finally, the drug sensitivity analysis showed that the high-risk score patients might be more sensitive to EGFR, MEK and ERK inhibitors than low-risk score patients. In conclusion, our study suggested that LINC02432/hsa-miR-98–5p/HK2 axis may serve as a novel diagnostic, prognostic, and therapeutic target in PAAD treatment.