- Department of General Surgery, State Key Laboratory of Complex Severe and Rare Diseases, Peking Union Medical College Hospital, Chinese Academy of Medical Sciences, Peking Union Medical College, Beijing, China
Introduction: Previous studies have investigated the prognostic significance of glycolysis markers in pancreatic cancer; however, conclusions from these studies are still controversial.
Methods: PubMed, Embase, and Web of Science were systematically searched to investigate the prognostic role of glycolysis markers in pancreatic cancer up to May 2022. Pooled hazard ratios (HRs) with 95% confidence intervals (CIs) related to overall survival (OS), disease free survival (DFS), recurrence-free survival (RFS), and distant metastasis-free survival (DMFS) were calculated using the STATA 12.0 software.
Results: A total of 28 studies comprising 2010 patients were included in this meta-analysis. High expression of the five glycolysis markers was correlated with a poorer OS (HR = 1.72, 95% CI: 1.34-2.22), DFS (HR = 3.09, 95% CI: 1.91-5.01), RFS (HR = 1.73, 95% CI: 1.21-2.48) and DMFS (HR = 2.60, 95% CI: 1.09-6.20) in patients with pancreatic cancer. In subgroup analysis, it was shown that higher expression levels of the five glycolysis markers were related to a poorer OS in Asians (HR = 1.85, 95% CI: 1.46-2.35, P < 0.001) and Caucasians (HR = 1.97, 95% CI: 1.40-2.77, P < 0.001). Besides, analysis based on the expression levels of specific glycolysis markers demonstrated that higher expression levels of GLUT1 (HR = 2.11, 95% CI: 1.58-2.82, P < 0.001), MCT4 (HR = 2.26, 95% CI: 1.36-3.76, P = 0.002), and ENO1 (HR = 2.16, 95% CI: 1.28-3.66, P =0.004) were correlated with a poorer OS in patients with pancreatic cancer.
Conclusions: High expression of the five glycolysis markers are associated with poorer OS, DFS, RFS and DMFS in patients with pancreatic cancer, indicating that the glycolysis markers could be potential prognostic predictors and therapeutic targets in pancreatic cancer.
Introduction
Pancreatic cancer is the seventh leading cause of cancer-related death worldwide due to its poor prognosis with almost as many deaths (n=466,003) as cases (n=495,773) (1), which was estimated to become the second most common cause of cancer-related death by 2030 in the U.S (2). Pancreatic ductal adenocarcinoma (PDAC), the most common type of pancreatic cancer, is the deadliest malignancy with a 5-year survival rate of less than 8% (3). The worse clinical outcomes of PDAC are due to the lack of early typical symptoms and the highly aggressive biological characteristics, emphasizing the urgent need for searching promising prognostic markers for clinical practice.
Malignant proliferating cancer cells usually rewire energy metabolism manners to meet the requirements of rapid cell growth, division, invasion, and migration, which is an emerging hallmark of cancer known as “energy metabolism reprogramming” (4). It has been well known that tumor cells tend to increase their glucose uptake and lactate production even in the presence of ample oxygen, which is described as the “aerobic glycolysis” or “Warburg effect” (5). The enhancement of glycolysis has been demonstrated to play critical roles in pancreatic tumorigenesis and development (6, 7).
Previous studies have shown that multiple key glycolytic enzymes such as glucose transporter 1 (GLUT1), monocarboxylate transporter 4 (MCT4), hexokinase 2 (HK2), pyruvate kinase M2 (PKM2) and enolase 1 (ENO1) are usually overexpressed in tumors and play critical roles in glycolysis pathway (8). Transport of glucose across the plasma membrane is an important regulator of glucose metabolism and is mediated by glucose transporter (GLUT) family proteins. Increased rate of glucose uptake in malignant cells is associated with elevated expression of GLUTs, especially GLUT1 (9). Since enhanced glycolysis is associated with lactate production, tumor cells must eliminate excessive lactate out of the cell to prevent cellular acidification, which is achieved by the upregulation of MCT4, a proton-coupled lactate transporter (10). Besides cellular membrane transporters, several enzymes located in cytoplasmic fractions also contribute to glycolytic flux. As the most influential functional enzyme of glycolysis, overexpression of HK2 promotes glycolysis process by catalyzing phosphorylation of glucose to glucose-6-phosphate (11). PKM2 is the key rate-limiting enzyme involved in the final step of the glycolysis pathway and catalyzes the last irreversible step in glycolysis, the conversion of phosphoenolpyruvate to pyruvate, while ADP is phosphorylated to form ATP (12). Moreover, ENO1, also known as alpha-enolase, a critical glycolytic enzyme that dehydrates 2-phosphoglycerate to phosphoenolpyruvate, is present both in cell surface and cytoplasm (13). It has been well demonstrated that these glycolysis enzymes contribute to biological behaviors of cancer.
Although emerging evidence has focused on the correlation between various glycolysis markers and pancreatic cancer, the conclusions remain controversial. In terms of GLUT1, previous studies investigating the relationship between GLUT1 expression and clinical outcomes of patients with pancreatic cancer have yet yielded conflicting results. For example, some studies have demonstrated the independent prognostic role of GLUT1 in pancreatic cancer (14–19), while other investigators revealed that the expression of GLUT1 was not associated with clinicopathological characteristics and prognosis of pancreatic cancer (20, 21). Therefore, a comprehensive meta-analysis was performed to clarify the prognostic role of glycolysis markers GLUT1, MCT4, HK2, PKM2, and ENO1 in pancreatic cancer.
Materials and methods
Literature search
Eligible studies relevant to the association of glycolysis markers and prognosis of pancreatic cancer were systematically searched in databases, including PubMed, Embase, and Web of Science, dated until May 2022. The language was restricted to English. The retrieval strategy was listed as follow: (glucose transporter 1 OR GLUT1 OR monocarboxylate transporter 4 OR MCT4 OR hexokinase 2 OR HK2 OR pyruvate kinase M2 OR PKM2 OR Enolase 1 OR ENO1) AND (pancreatic OR pancreas) AND (cancer OR tumor OR tumor OR carcinoma OR adenocarcinoma OR neoplasia OR neoplasm) AND (prognosis OR prognostic OR prognoses OR survival OR outcome).
Inclusion and exclusion criteria
All included studies need to meet the following inclusion criteria (1): including patients with pathologically diagnosed pancreatic cancer; (2) glycolysis markers were measured using immunohistochemistry (IHC); (3) studies describing the correlation between the expression of glycolysis markers and survival outcome; (4) hazard ratio (HR) values and corresponding 95% confidence intervals (CIs) for overall survival (OS), disease free survival (DFS), recurrence-free survival (RFS) or distant metastasis-free survival (DMFS) were either reported or calculatable with the published data in studies. Exclusion criteria were as follows: (1) duplicated studies; (2) studies not based on humans subjects; (3) reviews, case reports, meta-analysis, and clinical trials; (4) studies without available data.
Data extraction
The following information was extracted from each study: first author’s name, publication year, country, ethnicity, types of glycolysis markers, sample size, gender, follow-up duration, cancer subtype, detection method, cut-off values of glycolysis markers and outcome. If HRs and 95% CIs were provided in both univariate and multivariate analyses, the latter was adopted. If a study only had Kaplan–Meier survival curves, survival data was extracted from the curves by Engauge Digitizer software.
Quality assessment
Newcastle-Ottawa Scale (NOS) was applied to comprehensively evaluate the quality of the included studies (22). NOS contained three domains: patient selection (0-4 points), comparability (0-2 points), and outcome (0-3 points). NOS scores ranged from 0-9 points. The studies with an NOS score of six or above were of high quality.
Statistical analysis
The association of glycolysis markers with OS, DFS, RFS and DMFS was evaluated by pooled HRs and 95% CIs. The pooled odds ratios (OR) and 95% CIs were used to investigate the association between the expression of glycolysis markers with the clinicopathologic characteristics of pancreatic cancer. Heterogeneity among various studies was evaluated using χ2 test and I2 statistic. The random-effects model was used when significant heterogeneity was detected (P < 0.1 or I2 > 50%); otherwise, the fixed-effects model was applied. Subgroup analysis was performed to explore the potential causes of heterogeneity. Publication bias was estimated using Begg’s funnel plot and Egger’s test. All statistics were analyzed using Stata software version 12.0 (STATA Corporation, College Station, TX, USA). P < 0.05 was considered statistically significant.
Results
Characteristics of included studies
The flow diagram for study selection is illustrated in Figure 1. The characteristics of included studies are listed in Table 1. The initial literature search retrieved a total of 1076 publications from three databases. Among them, 265 articles came from PubMed, 466 from Embase, and 345 from Web of Science. After removing duplicates and manually screening titles and abstracts, the full text of the remaining 32 studies were assessed for eligibility. Fourteen articles were subsequently excluded for the following reasons: 10 papers were excluded due to data insufficiency, two were bioinformatic analysis, one used irrelevant method, and one was about other tumor. Finally, a total of 18 studies published between the year 2007 to 2020 were included in this meta-analysis (14–21, 23–32), made up of 28 studies and 2010 patients. The quality of the included studies was assessed using NOS, with all scoring of six or above (Supplementary Table 1). Majority of these studies assessed GLUT1 (n=11), whereas the rest assessed PKM2 (n=9), HK2 (n=5), ENO1 (n=2) and MCT4 (n=1). The sample size varied from 36 to 223 patients. In all eligible studies, the expression of glycolysis markers was detected using IHC. HR values and 95% CIs were directly provided in 8 studies, while HRs and 95% CIs of other studies were extracted using Engauge Digitizer 4.1 software with reported survival curves. The prognostic value of specific glycolysis markers was investigated by evaluating the overall survival (OS) in 20 studies, disease-free survival (DFS) in 2 studies, recurrence-free survival (RFS) in 4 studies, and distant metastasis-free survival (DMFS) in 2 studies.
Sensitivity analysis indicated that no significant heterogeneity was present among the included studies (Figure 2). Then, the publication bias for OS was evaluated, where neither Begg’s test (Figure 3A) nor Egger’s test (Figure 3B) showed a significant publication bias (P=0.871 and P=0.245, respectively).
Figure 2 Sensitivity analysis between the expression levels of glycolysis markers and overall survival.
Figure 3 Evaluation of publication bias between the expression levels of glycolysis markers and overall survival. (A) Begg’s funnel plot, P=0.871; (B) Egger’s test, P=0.245.
Glycolysis markers and OS in pancreatic cancer
Next, the relationship between expression levels of the 5 glycolysis markers and OS in pancreatic cancer from 20 studies that included 1546 patients was evaluated. As described in Figure 4A, results indicated that a higher expression of glycolysis markers correlated with poor OS of pancreatic cancer (HR = 1.72, 95% CI: 1.34-2.22). Owing to a significant heterogeneity among studies (I2 = 56.3%, P=0.001), a random-effects model was applied. Then, subgroup analysis was applied to explore causes of heterogeneity. The regional subgroup analysis indicated that high expression levels of the five glycolysis markers correlated with poorer OS in Asians (HR = 1.85, 95% CI: 1.46-2.35, P < 0.001) and Caucasians (HR = 1.97, 95% CI: 1.40-2.77, P < 0.001) (Figure 4B). The decreased heterogeneity from 56.3% to 0%-27.8% following the subgroup analysis indicated that the differences in regions are the main cause of heterogeneity. In addition, it was shown that higher expression levels of GLUT1 (HR = 2.11, 95% CI: 1.58-2.82, P < 0.001), MCT4 (HR = 2.26, 95% CI: 1.36-3.76, P = 0.002), and ENO1 (HR = 2.16, 95% CI: 1.28-3.66, P =0.004) correlated with poorer OS in pancreatic cancer in the subgroup analysis based on the expression levels of specific glycolysis markers (Figure 4C). The various heterogeneity in different subgroups (GLUT1: I2 = 0%, HK2: I2 = 55.8%, PKM2: I2 = 75.1% and ENO1: I2 = 0%) suggested that different glycolysis markers may not be the main cause of heterogeneity.
Figure 4 Forest plots of pooled hazard ratio (HR) for the association between the expression levels of glycolysis markers and overall survival of pancreatic cancer patients. (A) The overall group; (B) Subgroup analysis for ethnicity; (C) Subgroup analysis for glycolysis markers types. A random-effects model was applied.
Glycolysis markers and DFS in pancreatic cancer
A total of 2 studies related to GLUT1 that included 140 patients were included in the pooled survival analysis of DFS. As shown in Figure 5A, results indicated that GLUT1 overexpression predicted a poor DFS of patients with pancreatic cancer (HR = 3.09, 95% CI: 1.91-5.01). A fixed-effects model was used since no significant heterogeneity exists between the 2 studies (I2 = 0%, P=0.398). However, subgroup analysis was not performed as the number of included studies was too small.
Figure 5 Forest plots of pooled hazard ratio (HR) for the association between the expression levels of glycolysis markers and (A) Disease free survival (B) Recurrence-free survival (C) Distant metastasis-free survival of pancreatic cancer patients. A fixed-effects model was used.
Glycolysis markers and RFS in pancreatic cancer
A total of 252 patients from 4 studies were evaluated to examine the correlation between expression levels of GLUT1, HK2 and PKM2 and RFS, and the meta-analysis of these studies showed that upregulation of the three glycolysis markers indicated a poor RFS of pancreatic cancer (HR = 1.73, 95% CI: 1.21-2.48) (Figure 5B). There was no significant heterogeneity between the studies (I2 = 0%, P=0.591), therefore, a fixed-effects model was used to calculate the pooled HRs. Likewise, subgroup analysis was not possible due to the limited number of studies that reported RFS.
Glycolysis markers and DMFS in pancreatic cancer
We then analyzed the data from 2 studies with 72 patients to determine the relationship between DMFS and the expression levels of HK2 and PKM2 in pancreatic cancer. We observed that significant correlation between expression levels of the two glycolysis markers and DMFS of patients with pancreatic cancer (HR = 2.60, 95% CI: 1.09-6.20) (Figure 5C). Moreover, there was no significant heterogeneity between the studies (I2 = 0%, P=0.854). Similarly, due to the limited number of included studies, subgroup analysis was not available.
Glycolysis markers and clinicopathologic characteristics of pancreatic cancer
The relationship between GLUT1 and PKM2 and clinicopathologic characteristics of pancreatic cancer was further analyzed. In terms of GLUT1, a total of 9 features from 8 studies including gender, age, tumor location, tumor differentiation, TNM stage, lymph node metastasis, perineural invasion, vascular invasion and resection margin were investigated in Table 2 and Supplementary Figure 1. However, the combined data indicated that the correlations between the expression of GLUT1 and the clinical characteristics were not significant. Moreover, we analyzed a total of 7 features from 9 studies that reported the association between the expression of PKM2 and the clinicopathologic features including gender, tumor location, tumor differentiation, TNM stage, lymph node metastasis, perineural invasion and resection margin (Table 3 and Supplementary Figure 2). Similarly, no statistical correlations were observed between PKM2 expression and the clinicopathologic characteristics of pancreatic cancer.
Discussion
Energy metabolism reprogramming including aberrant glucose, amino acid, lipid metabolism and other bioenergetic metabolism pathways confers tumor cells the ability to acquire necessary nutrients and biomass in order to survive and proliferate under a frequently severe and nutrient-limited microenvironment (33, 34). Among various metabolic manners, aerobic glycolysis which is the central pathway of dysregulation of glucose metabolism has been well demonstrated in many tumors. Notably, upregulation of key enzymes including GLUT1, HK2, phosphofructokinase (PFK), TP53-induced glycolysis and apoptosis regulator (TIGAR), PKM2, ENO1, lactate dehydrogenase A (LDHA) and MCT4 in glycolytic pathway contributes to the increased glycolytic flux and promotion of malignant behaviors of various cancers including pancreatic cancer (7, 35). In addition, therapeutic strategies for targeting glycolytic enzymes are emerging as promising anti-cancer candidates. 3-bromopyruvate (3-BrPA), an inhibitor of HK2, has been well proved to block the growth and progression of pancreatic cancer both in vitro and in vivo (36, 37). In addition to HK2, the other glycolytic rate-limiting enzyme PFK1 could be allosterically activated by fructose-2,6-bisphosphate (F2,6BP) mediated by well-established glucose metabolism regulators 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatases (PFKFBs) (38). PFK15, a PFKFB3 inhibitor, has been demonstrated to reduce glycolysis and plasma membrane calcium ATPases (PMCAs) activity, resulting in calcium overload and cell apoptosis in PDAC (39). Treating pancreatic cancer xenografts with FX11 which is a small-molecule inhibitor of LDHA suppressed tumor progression (40). Additionally, the combination of LDHA inhibitors and gemcitabine displayed synergistic cytotoxic activity against pancreatic cancer cells under hypoxic conditions (41). In terms of clinical applications, in a phase I clinical trial in patients with advanced solid tumors, including pancreatic cancer, the glycolytic inhibitor 2-deoxy-D-glucose (2DG) in combination with chemotherapeutic agent docetaxel produced tolerable adverse effects and feasible clinical benefit (NCT00096707) (42). Given the potential benefits for targeting glycolytic enzymes, more clinical trials of specific glycolysis inhibitors in patients with pancreatic cancer should be considered.
In the present meta-analysis, we analyzed data from 28 studies consisting of 2010 patients to determine the prognostic value of five glycolysis markers in pancreatic cancer. The results showed that high expression levels of glycolysis markers significantly correlated with worse OS, DFS, RFS and DMFS. The conclusion is consistent with the results of most studies included in this meta-analysis. In terms of subgroup analysis for OS, analysis by ethnicity showed that higher expression of glycolysis markers was correlated with shorter OS in both Asians and Caucasians. This implies that malignant glycolysis phenotypes of pancreatic cancer are common phenomenon in humans. Moreover, subgroup analysis of the expression levels of specific glycolysis markers highlighted that higher expression levels of GLUT1, MCT4, and ENO1 correlated with poorer OS in pancreatic cancer patients, while the expression levels of HK2 and PKM2 did not. These results indicated that the glycolysis pathways of pancreatic cancer may vary in diverse glycolysis markers and GLUT1, MCT4, and ENO1 may be potential therapeutic targets in treating pancreatic cancer. Due to the limited number of included studies, subgroup analysis for DFS, RFS and DMFS in pancreatic cancer patients was not performed in the present meta-analysis. Hence, more studies examining the correlation between expression levels of glycolysis markers and DFS, RFS and DMFS of pancreatic cancer patients need to be conducted. Moreover, we also analyzed the association between GLUT1 and PKM2 and clinicopathologic features including gender, age, tumor location, tumor differentiation, TNM stage, lymph node metastasis, perineural invasion, vascular invasion, and resection margin in pancreatic cancer. However, no significant correlations were observed. The reason may be attributed to limited number of studies and samples and the high levels of heterogeneity between various studies. Therefore, further studies are warranted to evaluate the role of glycolysis markers in clinicopathologic characteristics of pancreatic cancer patients in the future.
In recent years, 18F-fluorodeoxyglucose positron emission tomography (FDG-PET) based on the increased glucose metabolism of malignant cells has been being used as a valuable functional imaging technology in clinical practice for various malignancies. Altered expression levels of GLUT has been demonstrated to be associated with enhanced FDG uptake quantified by standardized uptake values (SUVs) in tumor cells (43). Interestingly, previous studies investigating the association between the expression of GLUT1 and SUVs in pancreatic cancer have yielded conflicting results. Chikamoto et al. suggested that a significant relationship existed between high maximum standardized uptake value (SUVmax) and the expression of GLUT1 (18), whereas other studies indicated that GLUT1 expression was not significantly correlated with SUVmax (14, 20). Similarly, a recent meta-analysis suggested that only a moderate association was identified between the expression of GLUT1 and SUVs derived from FDG-PET in multiple tumor types (44). Therefore, it can be deduced that the correlations between glucose hypermetabolism in vivo displayed by FDG-PET and enhanced glycolytic phenotype of tumor cells are more complex than the single expression of GLUT1. Other key proteins of the glucose metabolism such as other members of GLUT family and HK2 might contribute to the alterations in SUVs, which needs to be confirmed in further studies.
The pancreatic tumor microenvironment (TME) consisting of cancer cells, stromal cells, and extracellular components is considered to play a pivotal role in multiple malignant biological behaviors including carcinogenesis, proliferation, invasion, metastasis, angiogenesis, immunosuppression and chemoresistance. It should be noted that stromal cells could be induced by cancer cells to present a glycolytic phenotype and support mitochondrial oxidative phosphorylation (OXPHOS) in cancer cells to fuel cell proliferation and growth, which was described as “the reverse Warburg effect” (45). As a key component in pancreatic TME, pancreatic stellate cells (PSCs) could be induced a shift to glycolysis in Caveolin-1-ROS positive feedback signaling, presenting upregulated expression of glycolytic enzymes HK2, 6-phosphofructokinase (PFKP), PKM2 and GLUT1 with downregulated expression of OXPHOS enzymes (46). Cancer-associated fibroblasts (CAFs), which develops mostly from activated PSCs, showed elevated expression of lactate dehydrogenase (LDHA), PKM2 and enhanced glucose uptake capability and lactate production (47, 48). The enhanced glycolysis of both PSCs and CAFs contribute to the proliferative and invasive capacities of pancreatic cancer cells. Moreover, pancreatic cancer is characterized by extensive infiltration of immunosuppressive cells such as regulatory T cells (Tregs), tumor-associated macrophages (TAMs) and myeloid-derived suppressor cells (MDSCs). It has been well demonstrated that glycolytic by-product lactic acid not only blunts tumor immunosurveillance by T and NK cells but also enhances immunosuppressive phenotype of Tregs, TAMs and MDSCs (49–52). A recent study further indicated that a terminally differentiated subpopulation of tumor-associated neutrophils exhibited hyperactivated glycolytic activity in PDAC TME promotes pro-tumor and immunosuppression functions of neutrophils (53). Therefore, targeting enhanced glycolysis of multiple components in the TME, not just tumor cells, may contribute to the development of novel therapeutic strategies against PDAC.
Several limitations should be addressed for this meta-analysis. First, given the current lack of a universal cut-off value for various glycolysis markers, different cut-off values were applied in different included studies. Second, some of the HRs and 95% CIs were estimated by extracting data from the survival curves, which might produce statistical deviations. Third, many studies included did not provide follow-up durations and they may have contributed to the bias. Fourth, the sample size of some studies is small. Considering the limitations of the present study, additional high-quality, large-scale, long-term studies need to be conducted.
Conclusions
This meta-analysis demonstrated that high expression of the five glycolysis markers correlated with poorer OS, DFS, RFS and DMFS in pancreatic cancer patients, suggesting that glycolytic pathway enzymes are potential prognostic biomarkers and promising therapeutic targets in patients with pancreatic cancer. However, due to the limitations of this study, more high-quality and well-designed studies need to be performed in the future.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Author contributions
QX, CW, and RX conceived and designed the study. CW, RX, and JS drafted the manuscript, acquired the data, and performed the meta-analysis. YC, XY, and RR revised the manuscript for important intellectual contents. All authors read and approved the final manuscript.
Funding
This study was supported by the CAMS Innovation Fund for Medical Sciences (2021-1-I2M-002), National Natural Science Foundation of China (82102810 and 81970763) and fellowship of China Postdoctoral Science Foundation (2022T150067).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2022.1004850/full#supplementary-material
References
1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin (2021) 71(3):209–49. doi: 10.3322/caac.21660
2. Rahib L, Smith BD, Aizenberg R, Rosenzweig AB, Fleshman JM, Matrisian LM. Projecting cancer incidence and deaths to 2030: the unexpected burden of thyroid, liver, and pancreas cancers in the united states. Cancer Res (2014) 74(11):2913–21. doi: 10.1158/0008-5472.CAN-14-0155
3. Siegel RL, Miller KD, Fuchs HE, Jemal A. Cancer statistics, 2021. CA Cancer J Clin (2021) 71(1):7–33. doi: 10.3322/caac.21654
4. Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation. Cell (2011) 144(5):646–74. doi: 10.1016/j.cell.2011.02.013
5. Hsu PP, Sabatini DM. Cancer cell metabolism: Warburg and beyond. Cell (2008) 134(5):703–7. doi: 10.1016/j.cell.2008.08.021
6. Lunt SY, Vander Heiden MG. Aerobic glycolysis: meeting the metabolic requirements of cell proliferation. Annu Rev Cell Dev Biol (2011) 27:441–64. doi: 10.1146/annurev-cellbio-092910-154237
7. Yang J, Ren B, Yang G, Wang H, Chen G, You L, et al. The enhancement of glycolysis regulates pancreatic cancer metastasis. Cell Mol Life Sci (2020) 77(2):305–21. doi: 10.1007/s00018-019-03278-z
8. Altenberg B, Greulich KO. Genes of glycolysis are ubiquitously overexpressed in 24 cancer classes. Genomics (2004) 84(6):1014–20. doi: 10.1016/j.ygeno.2004.08.010
9. Macheda ML, Rogers S, Best JD. Molecular and cellular regulation of glucose transporter (GLUT) proteins in cancer. J Cell Physiol (2005) 202(3):654–62. doi: 10.1002/jcp.20166
10. Payen VL, Mina E, Van Hee VF, Porporato PE, Sonveaux P. Monocarboxylate transporters in cancer. Mol Metab (2020) 33:48–66. doi: 10.1016/j.molmet.2019.07.006
11. Garcia SN, Guedes RC, Marques MM. Unlocking the potential of HK2 in cancer metabolism and therapeutics. Curr Med Chem (2019) 26(41):7285–322. doi: 10.2174/0929867326666181213092652
12. Dayton TL, Jacks T, Vander Heiden MG. PKM2, cancer metabolism, and the road ahead. EMBO Rep (2016) 17(12):1721–30. doi: 10.15252/embr.201643300
13. Pancholi V. Multifunctional alpha-enolase: its role in diseases. Cell Mol Life Sci (2001) 58(7):902–20. doi: 10.1007/PL00000910
14. Takahashi M, Nojima H, Kuboki S, Horikoshi T, Yokota T, Yoshitomi H, et al. Comparing prognostic factors of glut-1 expression and maximum standardized uptake value by FDG-PET in patients with resectable pancreatic cancer. Pancreatology (2020) 20(6):1205–12. doi: 10.1016/j.pan.2020.07.407
15. Achalandabaso Boira M, Di Martino M, Gordillo C, Adrados M, Martin-Perez E. GLUT-1 as a predictor of worse prognosis in pancreatic adenocarcinoma: immunohistochemistry study showing the correlation between expression and survival. BMC Cancer (2020) 20(1):909. doi: 10.21203/rs.3.rs-18270/v1
16. Lu K, Yang J, Li DC, He SB, Zhu DM, Zhang LF, et al. Expression and clinical significance of glucose transporter-1 in pancreatic cancer. Oncol Lett (2016) 12(1):243–9. doi: 10.3892/ol.2016.4586
17. Pizzi S, Porzionato A, Pasquali C, Guidolin D, Sperti C, Fogar P, et al. Glucose transporter-1 expression and prognostic significance in pancreatic carcinogenesis. Histol Histopathol (2009) 24(2):175–85. doi: 10.14670/HH-24.175
18. Chikamoto A, Inoue R, Komohara Y, Sakamaki K, Hashimoto D, Shiraishi S, et al. Preoperative high maximum standardized uptake value in association with glucose transporter 1 predicts poor prognosis in pancreatic cancer. Ann Surg Oncol (2017) 24(7):2040–6. doi: 10.1245/s10434-017-5799-1
19. Kitasato Y, Yasunaga M, Okuda K, Kinoshita H, Tanaka H, Okabe Y, et al. Maximum standardized uptake value on 18F-fluoro-2-deoxy-glucose positron emission tomography/computed tomography and glucose transporter-1 expression correlates with survival in invasive ductal carcinoma of the pancreas. Pancreas (2014) 43(7):1060–5. doi: 10.1097/MPA.0000000000000185
20. Yang HJ, Xu WJ, Guan YH, Zhang HW, Ding WQ, Rong L, et al. Expression of glut-1 and HK-II in pancreatic cancer and their impact on prognosis and FDG accumulation. Transl Oncol (2016) 9(6):583–91. doi: 10.1016/j.tranon.2016.08.004
21. Lyshchik A, Higashi T, Hara T, Nakamoto Y, Fujimoto K, Doi R, et al. Expression of glucose transporter-1, hexokinase-II, proliferating cell nuclear antigen and survival of patients with pancreatic cancer. Cancer Invest (2007) 25(3):154–62. doi: 10.1080/07357900701208931
22. Stang A. Critical evaluation of the Newcastle-Ottawa scale for the assessment of the quality of nonrandomized studies in meta-analyses. Eur J Epidemiol (2010) 25(9):603–5. doi: 10.1007/s10654-010-9491-z
23. Baek G, Tse YF, Hu Z, Cox D, Buboltz N, McCue P, et al. MCT4 defines a glycolytic subtype of pancreatic cancer with poor prognosis and unique metabolic dependencies. Cell Rep (2014) 9(6):2233–49. doi: 10.1016/j.celrep.2014.11.025
24. Ogawa H, Nagano H, Konno M, Eguchi H, Koseki J, Kawamoto K, et al. The combination of the expression of hexokinase 2 and pyruvate kinase M2 is a prognostic marker in patients with pancreatic cancer. Mol Clin Oncol (2015) 3(3):563–71. doi: 10.3892/mco.2015.490
25. Mohammad GH, Olde Damink SW, Malago M, Dhar DK, Pereira SP. Pyruvate kinase M2 and lactate dehydrogenase a are overexpressed in pancreatic cancer and correlate with poor outcome. PloS One (2016) 11(3):e0151635. doi: 10.1371/journal.pone.0151635
26. Calabretta S, Bielli P, Passacantilli I, Pilozzi E, Fendrich V, Capurso G, et al. Modulation of PKM alternative splicing by PTBP1 promotes gemcitabine resistance in pancreatic cancer cells. Oncogene (2016) 35(16):2031–9. doi: 10.1038/onc.2015.270
27. Xu Y, Chu L, Yuan S, Yang Y, Yang Y, Xu B, et al. RGD-modified oncolytic adenovirus-harboring shPKM2 exhibits a potent cytotoxic effect in pancreatic cancer via autophagy inhibition and apoptosis promotion. Cell Death Dis (2017) 8(6):e2835. doi: 10.1038/cddis.2017.230
28. Li C, Zhao Z, Zhou Z, Liu R. PKM2 promotes cell survival and invasion under metabolic stress by enhancing warburg effect in pancreatic ductal adenocarcinoma. Dig Dis Sci (2016) 61(3):767–73. doi: 10.1007/s10620-015-3931-2
29. Lockney NA, Zhang M, Lu Y, Sopha SC, Washington MK, Merchant N, et al. Pyruvate kinase muscle isoenzyme 2 (PKM2) expression is associated with overall survival in pancreatic ductal adenocarcinoma. J Gastrointest Cancer (2015) 46(4):390–8. doi: 10.1007/s12029-015-9764-6
30. Hu H, Tu W, Chen Y, Zhu M, Jin H, Huang T, et al. The combination of PKM2 overexpression and M2 macrophages infiltration confers a poor prognosis for PDAC patients. J Cancer (2020) 11(8):2022–31. doi: 10.7150/jca.38981
31. Wang L, Bi R, Yin H, Liu H, Li L. ENO1 silencing impaires hypoxia-induced gemcitabine chemoresistance associated with redox modulation in pancreatic cancer cells. Am J Transl Res (2019) 11(7):4470–80.
32. Sun L, Guo C, Cao J, Burnett J, Yang Z, Ran Y, et al. Over-expression of alpha-enolase as a prognostic biomarker in patients with pancreatic cancer. Int J Med Sci (2017) 14(7):655–61. doi: 10.7150/ijms.18736
33. Pavlova NN, Thompson CB. The emerging hallmarks of cancer metabolism. Cell Metab (2016) 23(1):27–47. doi: 10.1016/j.cmet.2015.12.006
34. Sun L, Suo C, Li ST, Zhang H, Gao P. Metabolic reprogramming for cancer cells and their microenvironment: Beyond the warburg effect. Biochim Biophys Acta Rev Cancer (2018) 1870(1):51–66. doi: 10.1016/j.bbcan.2018.06.005
35. Baig MH, Adil M, Khan R, Dhadi S, Ahmad K, Rabbani G, et al. Enzyme targeting strategies for prevention and treatment of cancer: Implications for cancer therapy. Semin Cancer Biol (2019) 56:1–11. doi: 10.1016/j.semcancer.2017.12.003
36. Ota S, Geschwind JF, Buijs M, Wijlemans JW, Kwak BK, Ganapathy-Kanniappan S. Ultrasound-guided direct delivery of 3-bromopyruvate blocks tumor progression in an orthotopic mouse model of human pancreatic cancer. Target Oncol (2013) 8(2):145–51. doi: 10.1007/s11523-013-0273-x
37. Chapiro J, Sur S, Savic LJ, Ganapathy-Kanniappan S, Reyes J, Duran R, et al. Systemic delivery of microencapsulated 3-bromopyruvate for the therapy of pancreatic cancer. Clin Cancer Res (2014) 20(24):6406–17. doi: 10.1158/1078-0432.CCR-14-1271
38. Yi M, Ban Y, Tan Y, Xiong W, Li G, Xiang B. 6-Phosphofructo-2-kinase/fructose-2,6-biphosphatase 3 and 4: A pair of valves for fine-tuning of glucose metabolism in human cancer. Mol Metab (2019) 20:1–13. doi: 10.1016/j.molmet.2018.11.013
39. Richardson DA, Sritangos P, James AD, Sultan A, Bruce JIE. Metabolic regulation of calcium pumps in pancreatic cancer: role of phosphofructokinase-fructose-bisphosphatase-3 (PFKFB3). Cancer Metab (2020) 8:2. doi: 10.1186/s40170-020-0210-2
40. Le A, Cooper CR, Gouw AM, Dinavahi R, Maitra A, Deck LM, et al. Inhibition of lactate dehydrogenase a induces oxidative stress and inhibits tumor progression. Proc Natl Acad Sci U.S.A. (2010) 107(5):2037–42. doi: 10.1073/pnas.0914433107
41. Maftouh M, Avan A, Sciarrillo R, Granchi C, Leon LG, Rani R, et al. Synergistic interaction of novel lactate dehydrogenase inhibitors with gemcitabine against pancreatic cancer cells in hypoxia. Br J Cancer (2014) 110(1):172–82. doi: 10.1038/bjc.2013.681
42. Raez LE, Papadopoulos K, Ricart AD, Chiorean EG, Dipaola RS, Stein MN, et al. A phase I dose-escalation trial of 2-deoxy-D-glucose alone or combined with docetaxel in patients with advanced solid tumors. Cancer Chemother Pharmacol (2013) 71(2):523–30. doi: 10.1007/s00280-012-2045-1
43. Mochizuki T, Tsukamoto E, Kuge Y, Kanegae K, Zhao S, Hikosaka K, et al. FDG uptake and glucose transporter subtype expressions in experimental tumor and inflammation models. J Nucl Med (2001) 42(10):1551–5.
44. Meyer HJ, Wienke A, Surov A. Associations between GLUT expression and SUV values derived from FDG-PET in different tumors-a systematic review and meta analysis. PloS One (2019) 14(6):e0217781. doi: 10.1371/journal.pone.0217781
45. Pavlides S, Whitaker-Menezes D, Castello-Cros R, Flomenberg N, Witkiewicz AK, Frank PG, et al. The reverse warburg effect: aerobic glycolysis in cancer associated fibroblasts and the tumor stroma. Cell Cycle (2009) 8(23):3984–4001. doi: 10.4161/cc.8.23.10238
46. Shao S, Qin T, Qian W, Yue Y, Xiao Y, Li X, et al. Positive feedback in cav-1-ROS signalling in PSCs mediates metabolic coupling between PSCs and tumour cells. J Cell Mol Med (2020) 24(16):9397–408. doi: 10.1111/jcmm.15596
47. Shan T, Chen S, Chen X, Lin WR, Li W, Ma J, et al. Cancer-associated fibroblasts enhance pancreatic cancer cell invasion by remodeling the metabolic conversion mechanism. Oncol Rep (2017) 37(4):1971–9. doi: 10.3892/or.2017.5479
48. Chen S, Chen X, Shan T, Ma J, Lin W, Li W, et al. MiR-21-mediated metabolic alteration of cancer-associated fibroblasts and its effect on pancreatic cancer cell behavior. Int J Biol Sci (2018) 14(1):100–10. doi: 10.7150/ijbs.22555
49. Brand A, Singer K, Koehl GE, Kolitzus M, Schoenhammer G, Thiel A, et al. LDHA-associated lactic acid production blunts tumor immunosurveillance by T and NK cells. Cell Metab (2016) 24(5):657–71. doi: 10.1016/j.cmet.2016.08.011
50. Watson MJ, Vignali PDA, Mullett SJ, Overacre-Delgoffe AE, Peralta RM, Grebinoski S, et al. Metabolic support of tumour-infiltrating regulatory T cells by lactic acid. Nature (2021) 591(7851):645–51. doi: 10.1038/s41586-020-03045-2
51. Colegio OR, Chu NQ, Szabo AL, Chu T, Rhebergen AM, Jairam V, et al. Functional polarization of tumour-associated macrophages by tumour-derived lactic acid. Nature (2014) 513(7519):559–63. doi: 10.1038/nature13490
52. Yang X, Lu Y, Hang J, Zhang J, Zhang T, Huo Y, et al. Lactate-modulated immunosuppression of myeloid-derived suppressor cells contributes to the radioresistance of pancreatic cancer. Cancer Immunol Res (2020) 8(11):1440–51. doi: 10.1158/2326-6066.CIR-20-0111
Keywords: pancreatic cancer, glycolysis, prognosis, therapeutic targets, meta-analysis
Citation: Wang C, Xu R, Song J, Chen Y, Yin X, Ruze R and Xu Q (2022) Prognostic value of glycolysis markers in pancreatic cancer: A systematic review and meta-analysis. Front. Oncol. 12:1004850. doi: 10.3389/fonc.2022.1004850
Received: 27 July 2022; Accepted: 23 August 2022;
Published: 12 September 2022.
Edited by:
Wenkai Ni, Fudan University, ChinaReviewed by:
Mingbing Xiao, Affiliated Hospital of Nantong University, ChinaKetao Wang, Beijing Tongren Hospital, Capital Medical University, China
Copyright © 2022 Wang, Xu, Song, Chen, Yin, Ruze and Xu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Qiang Xu, xuqiang@pumch.cn
†These authors have contributed equally to this work