SYSTEMATIC REVIEW article

Front. Oncol., 12 September 2022

Sec. Cancer Metabolism

Volume 12 - 2022 | https://doi.org/10.3389/fonc.2022.1004850

Prognostic value of glycolysis markers in pancreatic cancer: A systematic review and meta-analysis

  • 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

Abstract

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 (1419), 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 (1421, 2332), 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.

Figure 1

Table 1

Author yearCountryEthnicityGlycolysis markerSample sizeGender (M/F)Follow-up (months)PathologyDetection methodCut-off valueOutcomeHR [95% CI]NOS score
Yang 2016 (20)ChinaAsianGLUT15034/1617 (6–35)PCIHCScore≥9 (0-12)OS2.71 (0.93-7.91)8
Yang 2016 (20)ChinaAsianHK25034/1617 (6-35)PCIHCScore≥4 (0-12)OS0.94 (0.34-2.56)8
Takahashi 2020 (14)JapanAsianGLUT110159/42NAPDACIHCNAOS
DFS
OS: 3.68 (1.76-7.67)
DFS: 3.53 (1.99-6.26)
6
Boira 2020 (15)SpainCaucasianGLUT13923/1616 (9.7–39.2)PDACIHCPositive cells ≥ 80%OS
DFS
OS: 1.52 (0.57-4.1)
DFS: 2.23 (0.91-5.47)
7
Lu 2016 (16)ChinaAsianGLUT15329/24NAPCIHCScore≥2 (0-4)OS2.57 (0.46-14.32)7
Lyshchik 2007 (21)JapanAsianGLUT17427/47NAPCIHCindex ≥3OS1.52 (0.86-2.7)7
Lyshchik 2007 (21)JapanAsianHK27427/47NAPCIHCindex ≥ 3OS0.64 (0.34-1.18)7
Pizzi 2009 (17)ItalyCaucasianGLUT16030/30NAPDACIHCScore≥3 (1–4)OS2.81 (1.1-8.0)7
Chikamoto 2017 (18)JapanAsianGLUT113876/6227PCIHCScore≥11% (0-100%)OS
RFS
OS: 1.97 (1.140-3.488)
RFS: 1.89 (1.18-3.01)
8
Kitasato 2014 (19)JapanAsianGLUT14121/20NAPDACIHCIntensity score≥moderateOS1.86 (0.43-8.11)7
Baek 2014 (23)USACaucasianMCT4223121/102NAPDACIHCScore≥6 (0-9)OS2.264 (1.365-3.756)7
Ogawa 2015 (24)JapanAsianHK23621/15NAPDACIHCScore≥5 (0-10)OS
LRFS
DMFS
OS: 2.57 (0.89-8.39)
LRFS: 2.41 (0.69-8.46)
DMFS: 2.41 (0.74-7.86)
7
Ogawa 2015 (24)JapanAsianPKM23621/15NAPDACIHCScore≥5 (0-10)OS
LRFS
DMFS
OS: 2.16 (0.82-6.1)
LRFS: 2.22 (0.71-6.97)
DMFS: 2.84 (0.79-10.29)
7
Mohammad 2016 (25)UKCaucasianPKM27239/33NAPCIHCScore≥4 (1-9)OS1.55 (0.84-2.86)7
Calabretta 2016 (26)ItalyCaucasianPKM24222/20NAPDACIHCScore≥4 (0-5)RFS1.12 (1-4.4)7
Xu 2017 (27)ChinaAsianPKM26037/23NAPCIHC≥50%OS1.97 (1-3.89)7
Li 2016 (28)ChinaAsianPKM29057/33NAPDACIHCIntensity score≥moderateOS2.21 (0.93-5.23)6
Lockney 2015 (29)USAAfrican American,
Caucasian
PKM211563/52NAPDACIHCScore≥3 (0-12)OS0.57 (0.36-0.91)7
Hu 2020 (30)ChinaAsianPKM27751/26NAPDACIHCScore≥6 (0-12)OS2.117 (1.309-3.426)7
Wang 2019 (31)ChinaAsianENO157NANAPDACIHCScore≥6 (0-16)OS1.44 (0.50-4.18)7
Sun 2017 (32)ChinaAsianENO1100NANAPCIHCScore≥4 (0-9)OS2.469 (1.348-4.522)7

Characteristics of included studies.

M/F, male/female; HR, hazard ratio; CI, confidence interval; IHC, immunohistochemistry; OS, overall survival, DFS, disease free survival, RFS, recurrence-free survival; DMFS, distant metastasis- free survival; NOS, Newcastle–Ottawa Scale; NA, not available.

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

Figure 3

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

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

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.

Table 2

FeaturesNo. of studiesTest for associationTest for heterogeneity
OR95% CIPChi2I2 (%)P
Gender (Male vs. Female)51.030.72-1.480.8752.4100.661
Age (<60 years vs. >60 years)21.030.49-2.180.941000.947
Tumor location (Head vs. Body and tail of pancreas)31.170.73-1.860.5190.1500.927
Tumor differentiation (poor vs. well/moderate)51.200.82-1.750.3591.1500.887
TNM stage (III + IV vs. I + II)21.470.66-3.260.3410.0700.799
Lymph node metastasis (Present vs. Absent)41.130.76-1.670.5431.4500.695
Perineural invasion (Present vs. Absent)31.030.53-1.980.9390.5100.775
Vascular invasion (Present vs. Absent)30.870.49-1.540.6330.5800.749
Resection margin (R1 vs. R0)21.060.52-2.170.8790.0400.841

The relationship between GLUT1 and clinicopathologic characteristics.

Table 3

FeaturesNo. of studiesTest for associationTest for heterogeneity
OR95% CIPChi2I2 (%)P
Gender (Male vs. Female)61.010.72-1.410.9751.2100.944
Tumor location (Head and neck vs. Body and tail of pancreas)30.970.61-1.530.8961.5900.452
Tumor differentiation (poor vs. well/moderate)61.230.87-1.740.2474.8200.438
TNM stage (II +III + IV vs. I)41.620.97-2.700.0634.1227.30.248
Lymph node metastasis (Present vs. Absent)51.340.81-2.200.2520.6300.960
Perineural invasion (Present vs. Absent)41.030.68-1.560.8892.0100.570
Resection margin (R1 vs. R0)20.760.42-1.380.3680.1200.729

The relationship between PKM2 and clinicopathologic characteristics.

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 (4952). 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.

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).

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.

Statements

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.

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.

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

    SungHFerlayJSiegelRLLaversanneMSoerjomataramIJemalAet 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

    RahibLSmithBDAizenbergRRosenzweigABFleshmanJMMatrisianLM. 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

    SiegelRLMillerKDFuchsHEJemalA. Cancer statistics, 2021. CA Cancer J Clin (2021) 71(1):733. doi: 10.3322/caac.21654

  • 4

    HanahanDWeinbergRA. Hallmarks of cancer: the next generation. Cell (2011) 144(5):646–74. doi: 10.1016/j.cell.2011.02.013

  • 5

    HsuPPSabatiniDM. Cancer cell metabolism: Warburg and beyond. Cell (2008) 134(5):703–7. doi: 10.1016/j.cell.2008.08.021

  • 6

    LuntSYVander HeidenMG. 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

    YangJRenBYangGWangHChenGYouLet 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

    AltenbergBGreulichKO. 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

    MachedaMLRogersSBestJD. 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

    PayenVLMinaEVan HeeVFPorporatoPESonveauxP. Monocarboxylate transporters in cancer. Mol Metab (2020) 33:4866. doi: 10.1016/j.molmet.2019.07.006

  • 11

    GarciaSNGuedesRCMarquesMM. Unlocking the potential of HK2 in cancer metabolism and therapeutics. Curr Med Chem (2019) 26(41):7285–322. doi: 10.2174/0929867326666181213092652

  • 12

    DaytonTLJacksTVander HeidenMG. PKM2, cancer metabolism, and the road ahead. EMBO Rep (2016) 17(12):1721–30. doi: 10.15252/embr.201643300

  • 13

    PancholiV. Multifunctional alpha-enolase: its role in diseases. Cell Mol Life Sci (2001) 58(7):902–20. doi: 10.1007/PL00000910

  • 14

    TakahashiMNojimaHKubokiSHorikoshiTYokotaTYoshitomiHet 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 BoiraMDi MartinoMGordilloCAdradosMMartin-PerezE. 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

    LuKYangJLiDCHeSBZhuDMZhangLFet 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

    PizziSPorzionatoAPasqualiCGuidolinDSpertiCFogarPet 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

    ChikamotoAInoueRKomoharaYSakamakiKHashimotoDShiraishiSet 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

    KitasatoYYasunagaMOkudaKKinoshitaHTanakaHOkabeYet 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

    YangHJXuWJGuanYHZhangHWDingWQRongLet 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

    LyshchikAHigashiTHaraTNakamotoYFujimotoKDoiRet 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

    StangA. 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

    BaekGTseYFHuZCoxDBuboltzNMcCuePet 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

    OgawaHNaganoHKonnoMEguchiHKosekiJKawamotoKet 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

    MohammadGHOlde DaminkSWMalagoMDharDKPereiraSP. 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

    CalabrettaSBielliPPassacantilliIPilozziEFendrichVCapursoGet 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

    XuYChuLYuanSYangYYangYXuBet 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

    LiCZhaoZZhouZLiuR. 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

    LockneyNAZhangMLuYSophaSCWashingtonMKMerchantNet 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

    HuHTuWChenYZhuMJinHHuangTet 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

    WangLBiRYinHLiuHLiL. 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

    SunLGuoCCaoJBurnettJYangZRanYet 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

    PavlovaNNThompsonCB. The emerging hallmarks of cancer metabolism. Cell Metab (2016) 23(1):2747. doi: 10.1016/j.cmet.2015.12.006

  • 34

    SunLSuoCLiSTZhangHGaoP. Metabolic reprogramming for cancer cells and their microenvironment: Beyond the warburg effect. Biochim Biophys Acta Rev Cancer (2018) 1870(1):5166. doi: 10.1016/j.bbcan.2018.06.005

  • 35

    BaigMHAdilMKhanRDhadiSAhmadKRabbaniGet al. Enzyme targeting strategies for prevention and treatment of cancer: Implications for cancer therapy. Semin Cancer Biol (2019) 56:111. doi: 10.1016/j.semcancer.2017.12.003

  • 36

    OtaSGeschwindJFBuijsMWijlemansJWKwakBKGanapathy-KanniappanS. 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

    ChapiroJSurSSavicLJGanapathy-KanniappanSReyesJDuranRet 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

    YiMBanYTanYXiongWLiGXiangB. 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:113. doi: 10.1016/j.molmet.2018.11.013

  • 39

    RichardsonDASritangosPJamesADSultanABruceJIE. 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

    LeACooperCRGouwAMDinavahiRMaitraADeckLMet 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

    MaftouhMAvanASciarrilloRGranchiCLeonLGRaniRet 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

    RaezLEPapadopoulosKRicartADChioreanEGDipaolaRSSteinMNet 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

    MochizukiTTsukamotoEKugeYKanegaeKZhaoSHikosakaKet al. FDG uptake and glucose transporter subtype expressions in experimental tumor and inflammation models. J Nucl Med (2001) 42(10):1551–5.

  • 44

    MeyerHJWienkeASurovA. 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

    PavlidesSWhitaker-MenezesDCastello-CrosRFlomenbergNWitkiewiczAKFrankPGet al. The reverse warburg effect: aerobic glycolysis in cancer associated fibroblasts and the tumor stroma. Cell Cycle (2009) 8(23):39844001. doi: 10.4161/cc.8.23.10238

  • 46

    ShaoSQinTQianWYueYXiaoYLiXet 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

    ShanTChenSChenXLinWRLiWMaJet 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

    ChenSChenXShanTMaJLinWLiWet 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

    BrandASingerKKoehlGEKolitzusMSchoenhammerGThielAet 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

    WatsonMJVignaliPDAMullettSJOveracre-DelgoffeAEPeraltaRMGrebinoskiSet 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

    ColegioORChuNQSzaboALChuTRhebergenAMJairamVet al. Functional polarization of tumour-associated macrophages by tumour-derived lactic acid. Nature (2014) 513(7519):559–63. doi: 10.1038/nature13490

  • 52

    YangXLuYHangJZhangJZhangTHuoYet 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

  • 53

    WangLLiuYDaiYTangXYinTWangCet al. Single-cell RNA-seq analysis reveals BHLHE40-driven pro-tumour neutrophils with hyperactivated glycolysis in pancreatic tumour microenvironment. Gut (2022). doi: 10.1136/gutjnl-2021-326070

Summary

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

Volume

12 - 2022

Edited by

Wenkai Ni, Fudan University, China

Reviewed by

Mingbing Xiao, Affiliated Hospital of Nantong University, China; Ketao Wang, Beijing Tongren Hospital, Capital Medical University, China

Updates

Copyright

*Correspondence: Qiang Xu,

†These authors have contributed equally to this work

This article was submitted to Cancer Metabolism, a section of the journal Frontiers in Oncology

Disclaimer

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics