Identification and validation of a six-gene signature associated with glycolysis to predict the prognosis of patients with cervical cancer
BMC Cancer volume 20, Article number: 1133 (2020)
Cervical cancer (CC) is one of the most common gynaecological cancers. The gene signature is believed to be reliable for predicting cancer patient survival. However, there is no relevant study on the relationship between the glycolysis-related gene (GRG) signature and overall survival (OS) of patients with CC.
We extracted the mRNA expression profiles of 306 tumour and 13 normal tissues from the University of California Santa Cruz (UCSC) Database. Then, we screened out differentially expressed glycolysis-related genes (DEGRGs) among these mRNAs. All patients were randomly divided into training cohort and validation cohort according to the ratio of 7: 3. Next, univariate and multivariate Cox regression analyses were carried out to select the GRG with predictive ability for the prognosis of the training cohort. Additionally, risk score model was constructed and validated it in the validation cohort.
Six mRNAs were obtained that were associated with patient survival. The filtered mRNAs were classified into the protective type (GOT1) and the risk type (HSPA5, ANGPTL4, PFKM, IER3 and PFKFB4). Additionally, by constructing the prognostic risk score model, we found that the OS of the high-risk group was notably poorer, which showed good predictive ability both in training cohort and validation cohort. And the six-gene signature is a prognostic indicator independent of clinicopathological features. Through the verification of PCR, the results showed that compared with the normal cervial tissuses, the expression level of six mRNAs were significantly higher in the CC tissue, which was consistent with our findings.
We constructed a glycolysis-related six-gene signature to predict the prognosis of patients with CC using bioinformatics methods. We provide a thorough comprehension of the effect of glycolysis in patients with CC and provide new targets and ideas for individualized treatment.
Cervical cancer (CC) is one of the most common gynaecological cancers, accounting for the fourth cause of cancer-related death in women . According to global cancer statistics, in 2018, there were nearly 570 thousand new cases of CC worldwide, with approximately 310 thousand deaths . In recent years, with the launch of the human papillomavirus (HPV) vaccination programme, the incidence of CC in developed countries has significantly decreased, but in developing countries, it is still on the rise , and the age of onset tends to be younger. In addition, a large proportion of patients with CC are found to be in an advanced stage, and at this time, treatment options are extremely limited, and side effects are more serious, with a 5-year survival rate of less than 20% [3,4,5]. Moreover, patients with the same clinical stage and pathological type tend to adopt the same treatment, but the prognosis of patients is different, which is mainly due to the genetic heterogeneity of patients. Therefore, it is necessary to identify effective biomarkers to predict the prognosis of patients with CC.
Aerobic glycolysis is a special mechanism of tumour cell metabolism, also known as the Warburg effect , which plays an important role in promoting the growth and metastasis of various tumours, including CC. Some studies have found that HPV protein can promote the development of cancer through the Warburg effect, and the Warburg effect may also contribute to the enhancement of virus replication ability in the early stage of HPV infection . In addition, total lesion glycolysis (TLG) is a measure of tumour metabolic activity, and some retrospective studies have found that TLG was significantly related to the recurrence-free survival (RFS) and overall survival (OS) of patients with CC [8,9,10]. Moreover, some glycolytic enzymes have also been proven to be related to the prognosis of CC. For example, hexokinase-2 (HK2), an enzyme that catalyses the conversion of glucose into glucose-6-phosphate, is overexpressed in a variety of cancers, has a promoting effect on the occurrence and development of CC, and is significantly related to the prognosis of patients . Lactate dehydrogenase A (LDHA) and phosphofructokinase P (PFKP) were found to be significantly correlated with progression-free survival (PFS) and OS in patients with CC, and the expression level of LDHA in recurrent tumours was significantly higher than that in nonrecurrent tumours . Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) is a classical glycolytic enzyme that has been reported to be significantly increased in CC . However, single gene biomarkers cannot achieve a good prediction effect, and some studies have suggested that gene signatures may be a better choice for predicting patient outcomes.
By mining public databases, many relevant studies have found that the glycolysis-related gene (GRG) signature is closely associated with the prognosis of patients with cancer, such as lung adenocarcinoma , liver cancer , pancreatic ductal carcinoma  and endometrial carcinoma . Nevertheless, there is no bioinformatics research on this in CC. Thus, in this study, we analysed the relationship between the GRG signature and CC through the University of California Santa Cruz (UCSC) Database, which helped us better assess the prognosis of patients and provided new insights for individualized treatment of patients with CC.
Acquisition of mRNA expression dataset
The workflow of the present study is displayed as Fig. 1. We extracted the mRNA expression profiles of 306 CC samples and 13 normal samples from the UCSC database (http://xena.ucsc.edu/) (Supplementary Table 1). For UCSC dataset, RNA-seq data (FPKM values) were normalized into log2 (FPKM+ 1). Then, patients with OS less than 30 days were excluded, and 273 CC patients were included. The clinical information of CC patients were collected which included age, grade, American Joint Committee on Cancer (AJCC) stage, T classification, N classification, M classification, and OS (Supplementary Table 2).
Identification and analysis of differentially expressed glycolysis-related genes (DEGRGs)
To obtain the cancer-related genes, we compared the expression levels of 16,208 mRNAs between CC and normal tissues using the limma package. A gene with false discovery rate (FDR) < 0.05 and |logFC| > 1 was defined as the DEG. Then, a list of GRGs (GO_GLYCOLYTIC_PROCESS, HALLMARK_GLYCOLYSIS, KEGG_GLYCOLYSIS_GLUCONEOGENESIS, REACTOME_GLYCOLYSIS, and BIOCARTA_GLYCOLYSIS_PATHWAY) was downloaded from the Molecular Signatures Database v5.1 (MSigDB). The Venn diagram was to identify the DEGRGs by combining GRGs and DEGs. Furthermore, to understand the potential function and related pathway of DEGRGs, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were carried out by the “clusterprofiler” package.
Construction and validation of the prognostic risk score signature
Previous researches indicated that GRGs have potential prognostic value for cancer patients, but the role in CC patients remains unclear. Therefore, the survival analyses of GRGs in CC patients were performed in the present study. First, according to the ratio of 7: 3, all CC patients were randomly divided into training cohort (193 patients) and validation cohort (80 patients) (Table 1). In our research, the prognostic signature was developed in the training cohort and validated in the validation cohort. To identify the prognostic value of DEGRGs in CC patients, we performed univariate Cox regression analysis to confirm OS-related DEGRGs in training cohort. Then, based on the OS-related DEGRGs, a stepwise model selection by the Akaike information criterion (AIC) was used to avoid overfitting, and significant genes were incorporated into the multivariate Cox analysis to construct the GRG-based prognostic signature. By linear combination of the expression values of the selected genes weighted with the regression coefficients of the multivariate Cox regression analysis, the prognostic risk score model was established as follows:
(βi is the coefficient of the gene; i in multivariate Cox analysis; Gi represents the expression value of gene i; n is the number of genes in the signature).
The time-dependent receiver operating characteristic (ROC) curve and area under the curve (AUC) value were used to evaluate the discrimination of the prognostic model. Additionally, the optimal cutoff value of risk score was determined by performing X-tile software, and patients were divided into low- and high-risk groups. Kaplan-Meier (K-M) survival curve and the log-rank test were adopted to compare the prognosis between two groups.
To further study the prognostic value of GRG signature in an independent cohort, the risk score of each patient in the validation cohort. Time-dependent ROC curve, AUC value, and K-M survival curve were used to show the prognostic ability of GRG signature in the validation cohort.
Development of a clinical-GRG nomogram for CC patients
As reported in previous studies, several clinical variables were confirmed as prognostic variables, such as AJCC, T classification, N classification and M classification. Therefore, to develop a comprehensive prognostic model for CC patients, we firstly studied the prognostic value of clinical data, including age, grade, ATCC stage and TNM classification in CC patients by univariate Cox analysis. Moreover, the OS-related clinical data and GRG signature were incorporated into the multivariate Cox analysis and the independent variables were used to develop a nomogram. Concordance-index(C-index) and calibration curve were used to evaluate the nomogram.
Ten cases of CC and 10 cases of normal cervical tissue were obtained from the Department of gynecology and obstetrics, the First Affiliated Hospital of Wenzhou Medical University. Informed consent was provided by all patients. This study was approved by the ethics committee of the First Affiliated Hospital of Wenzhou Medical University. All tissues were pathologically diagnosed and stored in − 80 °C for preservation.
Conducting quantitative real time (qRT)-PCR on tissues
Total RNA was refined from clinical specimens by TRIzol® reagent (Vazyme, Nanjing, China), then, following the manufacturer’s protocol, a PrimeScript RT reagent kit (Vazyme) was used to reverse transcription. qRT-PCR was performed to evaluate the DEGRGs expression in different specimens using SYBR-Green Premix (Vazyme) with specifc PCR primers shown in Supplementary Table 3. With GAPDH as the internal control, and the Ct method (2-ΔΔCt) was used to normalize the relative genes expression values.
Correlation analysis between immune cell infiltration and glycolysis
We adopted CIBERSORT , which is widely used to describe the immune cell composition of the gene expression profiles to explain the 22 immune cell subtypes. The global P value of each sample deconvolution was determined by CIBERSORT, and samples with CIBERSORT P < 0.05 was selected for further study. Thus, we included a total of 184 patients from 193 patients. To determine the relationship between immune cell infiltration and glycolysis, we divided 184 patients into high- and low-risk groups based on the prognostic risk score and compared the differences in immune cells between the two groups and displayed the results by heatmap and violin plot. Furthermore, we conducted univariate Cox regression, LASSO and multivariate Cox regression analyses to examine whether the risk score and immune cells can be used as independent factors to predict the prognosis of CC patients.
Preliminary screening of glycolysis-related mRNA
We obtained the mRNA expression profiles from the UCSC database, including 306 tumour samples and 13 normal samples. By comparing tumour and normal tissue samples, we screened 4772 differentially expressed mRNAs from 16,208 mRNAs. Then, Venn diagram software was applied to identify the DEGRGs in the 4772 differentially expressed mRNAs and 5 glycolytic gene sets, and the results showed that 96 DEGRGs were detected (Fig. 2).
To further analyse the relationship between these genes and glycolysis, we carried out GO analysis and KEGG pathway enrichment analysis on these genes. As a result, we found in GO analysis that the primary enriched terms were ADP metabolic process, pyruvate metabolic process and glycolytic process (Fig. 3a). For KEGG pathway enrichment analysis, glycolysis/gluconeogenesis, fructose and mannose metabolism and carbon metabolism were the most enriched. Based on these results, it has been shown that the genes we selected are indeed related to glycolysis (Fig. 3b).
Identification of DEGRGs associated with survival and visualization of the DEGRG status
We analysed 96 mRNAs selected above to screen out survival-related mRNAs in training cohort by univariate Cox regression analysis and obtained 19 mRNAs with P < 0.05 (Fig. 4a). Next, we performed a stepwise model selection by the AIC, and six mRNAs (HSPA5, ANGPTL4, PFKM, GOT1, IER3 and PFKFB4) were obtained. Finally, multivariate Cox analysis was carried out to further identify the relationship of selected mRNAs with CC patients’ prognosis and obtained the coefficients. The filtered mRNAs were classified into the protective type (GOT1), whose HR < 1 was related to better prognosis, and the risk type (HSPA5, ANGPTL4, PFKM, IER3 and PFKFB4), whose HR > 1 was related to poor prognosis (Table 2).
In addition, to examine the DEGRG status of the whole genome, we visualized the data by using Circos plots , which are shown in Fig. 4b. The outer layer includes chromosomes and 96 DEGRGs. The middle four layers, from the outside to the inside, are the average expression values of DEGRGs in normal tissues, the average expression values of DEGRGs in tumour tissues, the logFC of the difference analysis, and the 19 meaningful prognostic DEGRGs obtained by univariate Cox regression analysis in the training cohort. The inner layer is the protein-protein interaction (PPI) network of the DEGRGs. In addition, 0.4 was defined as the minimum required interaction score, and between 0.4 and 0.9 is displayed in blue, while greater than 0.9 is displayed in red.
Constructing and validating a six-mRNA signature to predict patient prognosis
Using the linear combination of the expression value of the selected genes and the regression coefficient of multivariate Cox regression analysis, the following predictive risk scoring model was established: risk score = − 0.5224 × expression of GOT1 + 0.6886 × expression of HSPA5 + 0.1495 × expression of ANGPTL4 + 0.4280 × expression of PFKM + 0.3741 × expression of PFKFB4 + 0.2287 × expression of IER3. Based on the prognosis risk score, 193 patients were divided into high- and low-risk groups by using optimal risk score cutoff identified by X-tile as the boundary value (Fig. 5a), and the respective survival status of 193 patients was presented (Fig. 5b). The K-M analysis showed that compared to the high-risk group, the OS in the low-risk group was dramatically better (p < 0.0001; Fig. 5c). The AUCs for 1-, 3- and 5-year OS were 0.791, 0.731 and 0.782, respectively (Fig. 5d), suggesting that the six-mRNA signature has excellent diagnostic significance for prognosis prediction. Additionally, we generated a heatmap to exhibit the expression profiles of the six mRNAs, from which we can see that, compared with the low-risk group, the expression level of risk-type mRNAs (HSPA5, ANGPTL4, PFKM, IER3 and PFKFB4) of the high-risk group was apparently higher, while the expression level of the protective-type mRNA GOT1 was opposite (Fig. 5e).
Next, we validated the predictive ability of the six-mRNA signature in the validation cohort. Eighty patients were divided into high- and low-risk groups by using the same optimal risk score cutoff which were used in the training cohort (Fig. 6a), and the respective survival status of 80 patients was presented (Fig. 6b). And the K-M analysis showed that the patients in the low-risk group had significantly better OS (P = 0.021) (Fig. 6c), which AUC values were 0.664, 0.635 and 0.661 for 1-, 3- and 5-year, respectively (Fig. 6d).
The risk score constructed by the six-mRNA signature is an independent prognostic indicator
To evaluate whether the six-mRNA signature was an independent prognostic indicator, we performed univariate and multivariate Cox regression analyses to compare the prognostic values of risk score with other clinical characteristics. We included a total of 193 patients, and among these patients, 124 (64.2%) were younger than 50 years old, 69 (35.8%) were older than 50 years old. Among 193 patients, 103 (53.4%) had grade 1–2 tumours, and 72 (37.3%) had grade 3–4 tumours, 139 (72.0%) had T1-T2 tumours, 17 (8.8%) had T3-T4 tumours. Furthermore, among these patients, 92 (47.6%) had no lymph node metastasis, 32 (16.5%) had lymph node metastasis, 69 (35.7%) had no distant metastasis, and 8 (4.1%) had distant metastasis. Additionally, among these patients, 150 (77.7%) had stage I-II disease, and 38 (19.7%) had stage III-IV disease. In all the above data, except for age and grade, AJCC stage (HR = 2.78, 95% CI 1.55–4.98, p < 0.001), T classification (HR = 5.15, 95% CI 2.55–10.42, p < 0.001), M classification (HR = 4.70, 95% CI 1.53–14.50, p = 0.007) and N classification (HR = 4.18, 95% CI 1.83–9.53 p < 0.001) showed significant differences in the univariate Cox regression analysis, which can be used as independent prognostic factors (Table 3). In the next multivariate Cox regression analysis, only the N classification showed significant differences (HR = 14.80, 95% CI 3.44–63.76, p < 0.001) (Table 3). However, regardless of univariate or multivariate Cox regression analysis, the risk score showed significant prognostic value (HR = 2.72, 95% CI 1.92–3.85, p < 0.001; HR = 2.57, 95% CI 1.47–4.49, p = 0.002) (Table 3). Based on these results above, we constructed a nomogram prediction model combined the risk score with N classification to predict CC patients’ OS (Fig. 7a), which C index is 0.83. The calibration plot showed that in the nomogram, the predicted values of OS at 1, 3 and 5 years for CC patients have a good correlation with the actual values (Fig. 7b-d).
Next, stratified analysis was performed to further validate whether the six-mRNA signature can be used as an independent prognostic factor, according to the grade, AJCC, T classification and age. First, we classified patients into the grade 1–2 dataset and grade 3–4 dataset, and the patients in the datasets were divided into high- and low-risk groups, respectively. As a result, we found that there were significant differences in OS between these two groups, and the OS of the low-risk group was obviously better (Fig. 8a). Next, we did the same thing for the other clinicopathological features, and we found that the OS of the high-risk group was significantly worse than that of the low-risk group regardless of the AJCC stage (stage I-II or stage III-IV) (Fig. 8b), T classification (T0 or T1) (Fig. 8c) and age (younger than 50 or older than 50) (Fig. 8d), further confirming the reliability of our analysis.
Validating the overexpression of the six mRNAs in CC tissues by qRT-PCR
We validated the expression level of six mRNAs in 10 CC tissues and 10 normal cervical tissues by using qRT-PCR. The results showed that compared with the normal cervial tissuses, the expression level of HSPA5, ANGPTL4, PFKM, GOT1, IER3 and PFKFB4 were significantly higher in the CC tissue (Fig. 9), making the bioinformatics analysis result much more reliable and precise.
The relationship between immune cell infiltration and glycolysis
We analysed a total of 184 patients’ immune cell composition from their gene expression profiles. Additionally, we divided these patients into two groups, and we can see that the expression levels of CD8 T cells, regulatory T cells (Tregs) and resting mast cells were significantly higher in the low-risk group compared to the high-risk group; however, in the high-risk group, the expression levels of neutrophils, M0 macrophages, activated mast cells and resting CD4 memory T cells were dramatically higher (Fig. 10).
In addition, we analyzed the relationship between immune cells and prognosis of CC patients. In the univariate Cox regression analysis, CD8 T cells, resting CD4 memory T cells, M0 macrophages, M2 macrophages, resting mast cells, activated mast cells and Neutrophils all showed significant differences with p < 0.05. Among these immune cells, the HR of CD8 T cells, M2 macrophages and resting mast cells were less than 1, which associated with better prognosis. In the next multivariate Cox regression analysis, only activated mast cells showed significant differences with p < 0.05. And both in univariate or multivariate Cox regression analysis, the risk score showed significant prognostic value (p < 0.001) (Table 4).
In recent years, increasing attention has been paid to the relationship between energy metabolism and tumours, among which the Warburg effect is a notable feature of the energy metabolism of tumour cells. Despite the low efficiency of glycolysis, tumour cells are still active in glycolysis, which can produce more energy and various metabolites in a short time so that tumour cells can benefit from glycolysis . More importantly, the rapid growth of tumours is often accompanied by hypoxia, and hypoxia inducible factor (HIF), as the key regulator of glycolysis, can solve the nutritional needs of tumour cells by inducing the expression of glucose and amino acid transporters (glucose transporter type 1 (GLUT1)) and L-type amino acid transporter 1 (LAT1) to further promote the progression of cancer [20, 21]. There have been many studies on the relationship between glycolysis and tumours, but the relationship between GRGs and the prognosis of patients with CC is still very limited. Because a single gene biomarker cannot provide a strong prediction effect, a more accurate and reliable gene signature was used to predict the clinical outcomes of patients. For instance, a glycolysis-related nine-gene signature was used to predict the survival of patients with endometrial cancer . A tumour microenvironment-related nine-gene signature was used to predict overall survival with ovarian cancer . Moreover, an autophagy-related five-gene signature was developed for prognosis prediction in patients with prostate cancer . In this study, we used the GRG signature for the first time to predict the survival of patients with CC and obtained a good prediction effect.
First, we screened out the differentially expressed glycolytic mRNAs from the CC dataset in the UCSC database and selected six DEGRGs with predictive ability for the prognosis of patients with CC in the training cohort through univariate and multivariate Cox regression analysis (P < 0.05). Subsequently, by constructing the prognostic risk score model, we found that there was a significant difference in OS between the high- and low-risk groups, and the OS of the high-risk group was evidently worse. Furthermore, the risk score model was validated in the validation cohort with good predictive ability. Finally, further verification was done through PCR, we found that the expression of these six genes in tumor tissue was indeed significantly higher than that in normal tissue.
From our analysis, we selected six DEGRGs (HSPA5, ANGPTL4, PFKM, GOT1, IER3 and PFKFB4). Heat shock protein A5 (HSPA5, also konwn as GRP78), a member of the HSP70 family, has been found to play an important role in cancer malignancy and anti-tumor therapy. And it has been reported that the high expression of HSPA5 can protect cancer cells from immune surveillance, and inhibition of its expression can promote cell apoptosis and inhibit tumor growth [24, 25]. Angiopoietin-like4 (ANGPTL4) is highly expressed in various tumors, which is closely related to tumor growth and metastasis, such as hepatocellular carcinoma , breast cancer , head and neck squamous cell carcinoma  and melanoma . And ANGPTL4 participates in the construction of GRG signature in lung adenocarcinoma to predict the survival and metastasis of patients . In addition, previous studies have found that ANGPTL4 is significantly associated with the susceptibility of CC, is a potential risk factor . Muscular phosphofructokinase (PFKM), a member of the phosphofructokinase (PFK) family, can promote the growth of muscle-infiltrating bladder cancer  and can be used as a new breast cancer gene , and its expression level can distinguish normal tissues from CC tissues . 6-Phosphofructo-2-kinase/fructose-2,6-biphosphatase 4 (PFKFB4), is a key kinase in Warburg pathway , has been found to be associated with a variety of cancers, including breast cancer , prostate cancer  and glioblastoma , promoting the progression and metastasis of cancer, and may become an effective molecular target of anti-tumor drugs. It was found that the expression of immediate early response gene 3 (IER3) was increased in advanced cancer [38, 39], but some studies have found that IER3 can also promote tumour cell apoptosis and has anti-tumour activity, such as lower expression in CC tissues , and increasing the expression of IER3 can enhance the sensitivity of CC cells to radiotherapy . In contrast, in our study, we found that IER3 exists as a risk factor. Moreover, glutamate oxaloacetate transaminase 1 (GOT1), is a gene that encodes cytoplasmic aspartate aminotransferase and a key aspartate-producing protein. Glutamine metabolism is essential for the proliferation of cancer cells in addition to glucose. Glutamine derived glutamic acid is used to generate nonessential amino acids by GOT1 in high proliferative cells, which plays an important role in cell proliferation [42, 43]. Moreover, HIFɑ can inhibit the proliferation of tumor cells by inhibiting the synthesis of aspartate ; however, in our study, we found that GOT1 was a protective factor, with higher expression in the low-risk group than in the high-risk group. Hence, further studies are needed to explore the relationship between these two genes (IER3 and GOT1) and the prognosis of patients with CC.
Interestingly, we also found some divergences in immune cell infiltration between the high- and low-risk groups, suggesting that glycolysis may be correlated with tumour immunity. In addition, several previous studies have elucidated the relationship between glycolysis and immunity. Excess glycolysis will lead to acidic tumor microenvironment, affect the infiltration of immune cells in varying degrees, and contribute to the survival of cancer cells . Cascone et al. found that in patients with melanoma and NSCLC, high glycolytic activity will be accompanied by poor immune cell infiltration, such as cytoxic T cells, T helper cells, memory T cells, macrophases, or natural killer (NK) cells reduction . Li et al. found that high glycolysis group had a higher immunescore and higher TIL percentage in breast cancer, however, the immune cells with anti-tumor effect were not enriched, and concluded that high glycolysis was associated with immunosuppression of tumor microenvironment . In our study, the expression levels of CD8 T cells and resting mast cells were significantly higher in the low-risk group. CD8 T cells, to our knowledge, play a pivotal role in the control of tumour cell growth, and their relationship with the prognosis of patients with CC has also been confirmed in a previous study [48, 49]. For mast cells, some studies have found that mast cell infiltration in CC indicates poor clinical prognosis [49,50,51,52], which was consistent with our findings that resting mast cells were highly expressed in the low-risk group and that activated mast cells were highly expressed in the high-risk group. These findings are consistent with previous studies and to some extent explain why patients in the low-risk group have better survival outcomes.
It is undeniable that this study does have some defects. On the one hand, we obtained clinical gene information for only 273 cases of CC, and the sample size was not large enough. On the other hand, it would be better to have an external validation cohort to verify the accuracy of the prediction model rather than the internal validation cohort. However, we verified these six DEGRGs between tumor and normal tissues, which to some extent compensated for this defect.
In summary, we first identified the relationship between the GRG signature and the prognosis of patients with CC using bioinformatics methods and found that patients in the high-risk group had significantly lower OS than those in the low-risk group. Furthermore, we found some relationships between the infiltration of immune cells and glycolysis.
In brief, we constructed a six-gene signature (HSPA5, ANGPTL4, PFKM, GOT1, IER3 and PFKFB4) to predict the prognosis of patients with CC, which was also a prognostic factor independent of clinicopathological features. These findings will provide us with new insights into the role of glycolysis in CC, guiding individualized treatment for patients with CC.
Total lesion glycolysis
University of California Santa Cruz
Differentially expressed glycolysis-related genes
Kyoto Encyclopedia of Genes and Genomes
American Joint Committee on Cancer
Receiver operating characteristic
Hypoxia inducible factor
Glucose transporter type 1
Epidermal growth factor receptors
Immediate early response gene 3
Glutamate oxaloacetate transaminase 1
Pancreatic ductal adenocarcinoma
Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.
Shrestha AD, Neupane D, Vedsted P, Kallestrup P. Cervical Cancer prevalence, incidence and mortality in low and middle income countries: a systematic review. Asian Pac J Cancer Prev. 2018;19(2):319–24.
Frenel JSLTC, O'Neil B, Ott PA, Piha-Paul SA, Gomez-Roca C, van Brummelen EMJ, Rugo HS, Thomas S, Saraf S. Safety and efficacy of Pembrolizumab in advanced, programmed death ligand 1-positive cervical Cancer_ results from the phase Ib KEYNOTE-028 trial. J Clin Oncol. 2017;35(36):4035–61.
Pfaendler KS, Tewari KS. Changing paradigms in the systemic treatment of advanced cervical cancer. Am J Obstet Gynecol. 2016;214(1):22–30.
Liontos M, Kyriazoglou A, Dimitriadis I, Dimopoulos MA, Bamias A. Systemic therapy in cervical cancer: 30 years in review. Crit Rev Oncol Hematol. 2019;137:9–17.
Vander Heiden MG, Cantley LC, Thompson CB. Understanding the Warburg effect: the metabolic requirements of cell proliferation. Science. 2009;324(5930):1029–33.
Martinez-Ramirez I, Carrillo-Garcia A, Contreras-Paredes A, Ortiz-Sanchez E, Cruz-Gregorio A, Lizano M. Regulation of cellular metabolism by high-risk human papillomaviruses. Int J Mol Sci. 2018;19(7):1839.
Leray H, Gabiache E, Courbon F, Brenot-Rossi I, Colineaux H, Lepage B, et al. FDG-PET/CT identifies predictors of survival in patients with locally advanced cervical carcinoma and Para-aortic lymph node involvement to increase treatment. J Nucl Med. 2020;61:1442.
Carpenter DJ, Jacobs CD, Wong TZ, Craciunescu O, Chino JP. Changes on midchemoradiation therapy fluorodeoxyglucose positron emission tomography for cervical cancer are associated with prognosis. Int J Radiation Oncol Biol Phys. 2019;105(2):356–66.
Leseur J, Roman-Jimenez G, Devillers A, Ospina-Arango JD, Williaume D, Castelli J, et al. Pre- and per-treatment 18F-FDG PET/CT parameters to predict recurrence and survival in cervical cancer. Radiother Oncol. 2016;120(3):512–8.
Liu C, Wang X, Zhang Y. The roles of HK2 on tumorigenesis of cervical Cancer. Technol Cancer Res Treat. 2019;18:1533033819871306.
Recondo G, Mezquita L, Bigot L, Galissant J, Frias RL, André F, et al. Preliminary results on mechanisms of resistance to ALK inhibitors: A prospective cohort from the MATCH-R study. Ann Oncol. 2018;29:vi21.
Kim JWKS, Han SM, Paik SY, Hur SY, Kim YW, Lee JM, Namkoong SE. Increased glyceraldehyde-3-phosphate dehydrogenase gene expression in human cervical cancers. Gynecol Oncol. 1998;71(2):266–9.
Zhang L, Zhang Z, Yu Z. Identification of a novel glycolysis-related gene signature for predicting metastasis and survival in patients with lung adenocarcinoma. J Transl Med. 2019;17(1):423.
Jiang LZL, Bi J, Guan Q, Qi A, Wei Q, He M, Wei M, Zhao L. Glycolysis gene expression profilings screen for prognostic risk signature of hepatocellular carcinoma. Aging (Albany NY). 2019;11(23):10861–82.
Tian G, Li G, Liu P, Wang Z, Li N. Glycolysis-based genes associated with the clinical outcome of pancreatic ductal adenocarcinoma identified by the Cancer genome atlas data analysis. DNA Cell Biol. 2020;39(3):417–27.
Wang ZH, Zhang YZ, Wang YS, Ma XX. Identification of novel cell glycolysis related gene signature predicting survival in patients with endometrial cancer. Cancer Cell Int. 2019;19:296.
Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, Hoang CD, Diehn M, Alizadeh AA. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. 2015;12(5):453–7.
Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA. Circos: an information aesthetic for comparative genomics. Genome Res. 2009;19(9):1639–45.
Maxwell PHDG, Gleadle JM, Nicholls LG, Harris AL, Stratford IJ, Hankinson O, Pugh CW, Ratcliffe PJ. Hypoxia-inducible factor-1 modulates gene expression in solid tumors and influences both angiogenesis and tumor growth. Proc Natl Acad Sci U S A. 1997;94(15):8104–9.
Parks SK, Chiche J, Pouyssegur J. Disrupting proton dynamics and energy metabolism for cancer therapy. Nat Rev Cancer. 2013;13(9):611–23.
Ding QDS, Wang R, Zhang K, Wang H, Zhou X, Wang J, Wong K, Long Y, Zhu S. A nine-gene signature related to tumor microenvironment predicts overall survival with ovarian cancer. Aging (Albany NY). 2020;12(6):4879–95.
Hu D, Jiang L, Luo S, Zhao X, Hu H, Zhao G, Tang W. Development of an autophagy-related gene expression signature for prognosis prediction in prostate cancer patients. J Transl Med. 2020;18(1):160.
AS L. The glucose-regulated proteins: stress induction and clinical applications. Trends Biochem Sci. 2001;26(8):504–10.
Jamora CDG, Lee AS. Inhibition of tumor progression by suppression of stress protein GRP78/BiP induction in fibrosarcoma B/C10ME. Proc Natl Acad Sci U S A. 1996;93(15):7690–4.
Li H, Ge C, Zhao F, Yan M, Hu C, Jia D, et al. Hypoxia-inducible factor 1 alpha-activated angiopoietin-like protein 4 contributes to tumor metastasis via vascular cell adhesion molecule-1/integrin beta1 signaling in human hepatocellular carcinoma. Hepatology. 2011;54(3):910–9.
Padua D, Zhang XH, Wang Q, Nadal C, Gerald WL, Gomis RR, Massague J. TGFbeta primes breast tumors for lung metastasis seeding through angiopoietin-like 4. Cell. 2008;133(1):66–77.
Chiang KH, Shieh JM, Shen CJ, Chang TW, Wu PT, Hsu JY, Tsai JP, Chang WC, Chen BK. Epidermal growth factor-induced COX-2 regulates metastasis of head and neck squamous cell carcinoma through upregulation of angiopoietin-like 4. Cancer Sci. 2020;111(6):2004–15.
Yang W, Khoury E, Guo Q, Prabhu SA, Emond A, Huang F, et al. MNK1 signaling induces an ANGPTL4-mediated gene signature to drive melanoma progression. Oncogene. 2020;39(18):3650–65.
Rahmani F, Hasanzadeh M, Hassanian SM, Khazaei M, Esmaily H, Asef-Agah SA, Naghipour A, AF G, Avan A. Association of a genetic variant in the angiopoietin-like protein 4 gene with cervical cancer. Pathol Res Pract. 2020;216(7):153011.
Sun CM, Xiong DB, Yan Y, Geng J, Liu M, Yao XD. Genetic alteration in phosphofructokinase family promotes growth of muscle-invasive bladder cancer. Int J Biol Markers. 2016;31(3):e286–93.
Ahsan H, Halpern J, Kibriya MG, Pierce BL, Tong L, Gamazon E, et al. A genome-wide association study of early-onset breast cancer identifies PFKM as a novel breast cancer gene and supports a common genetic spectrum for breast cancer at any age. Cancer Epidemiol Biomark Prev. 2014;23(4):658–69.
Mattarocci S, Abbruzzese C, Mileo AM, Carosi M, Pescarmona E, Vico C, et al. Identification of pivotal cellular factors involved in HPV-induced dysplastic and neoplastic cervical pathologies. J Cell Physiol. 2014;229(4):463–70.
Goncalves MD, Cantley LC. A glycolysis outsider steps into the Cancer spotlight. Cell Metab. 2018;28(1):3–4.
Dasgupta S, Rajapakshe K, Zhu B, Nikolai BC, Yi P, Putluri N, et al. Metabolic enzyme PFKFB4 activates transcriptional coactivator SRC-3 to drive breast cancer. Nature. 2018;556(7700):249–54.
Ros S, Santos CR, Moco S, Baenke F, Kelly G, Howell M, Zamboni N, Schulze A. Functional metabolic screen identifies 6-phosphofructo-2-kinase/fructose-2,6-biphosphatase 4 as an important regulator of prostate cancer cell survival. Cancer Discov. 2012;2(4):328–43.
Goidts V, Bageritz J, Puccio L, Nakata S, Zapatka M, Barbus S, et al. RNAi screening in glioma stem-like cells identifies PFKFB4 as a key molecule important for cancer cell survival. Oncogene. 2012;31(27):3235–43.
Arlt A, Schafer H. Role of the immediate early response 3 (IER3) gene in cellular stress response, inflammation and tumorigenesis. Eur J Cell Biol. 2011;90(6–7):545–52.
HM S. Gomori-positive astrocytes: biological properties and implications for neurologic and neuroendocrine disorders. Glia. 1991;4(4):365–77.
Jin H, Lee K, Kim YH, Oh HK, Maeng YI, Kim TH, Suh DS, Bae J. Scaffold protein FHL2 facilitates MDM2-mediated degradation of IER3 to regulate proliferation of cervical cancer cells. Oncogene. 2016;35(39):5106–18.
Jin H, Suh DS, Kim TH, Yeom JH, Lee K, Bae J. IER3 is a crucial mediator of TAp73beta-induced apoptosis in cervical cancer and confers etoposide sensitivity. Sci Rep. 2015;5:8367.
Abrego J, Gunda V, Vernucci E, Shukla SK, King RJ, Dasgupta A, et al. GOT1-mediated anaplerotic glutamine metabolism regulates chronic acidosis stress in pancreatic cancer cells. Cancer Lett. 2017;400:37–46.
Birsoy K, Wang T, Chen WW, Freinkman E, Abu-Remaileh M, Sabatini DM. An essential role of the mitochondrial Electron transport chain in cell proliferation is to enable aspartate synthesis. Cell. 2015;162(3):540–51.
Melendez-Rodriguez F, Urrutia AA, Lorendeau D, Rinaldi G, Roche O, Bogurcu-Seidel N, et al. HIF1alpha suppresses tumor cell proliferation through inhibition of aspartate biosynthesis. Cell Rep. 2019;26(9):2257–65 e2254.
Gill KS, Fernandes P, O'Donovan TR, McKenna SL, Doddakula KK, Power DG, Soden DM, Forde PF. Glycolysis inhibition as a cancer treatment and its role in an anti-tumour immune response. Biochim Biophys Acta. 2016;1866(1):87–105.
Cascone T, McKenzie JA, Mbofung RM, Punt S, Wang Z, Xu C, et al. Increased tumor glycolysis characterizes immune resistance to adoptive T cell therapy. Cell Metab. 2018;27(5):977–87 e974.
Li W, Xu M, Li Y, Huang Z, Zhou J, Zhao Q, et al. Comprehensive analysis of the association between tumor glycolysis and immune/inflammation function in breast cancer. J Transl Med. 2020;18(1):92.
de Vos van Steenwijk PJ, Ramwadhdoebe TH, Goedemans R, Doorduijn EM, van Ham JJ, Gorter A, et al. Tumor-infiltrating CD14-positive myeloid cells and CD8-positive T-cells prolong survival in patients with cervical carcinoma. Int J Cancer. 2013;133(12):2884–94.
Yang S, Wu Y, Deng Y, Zhou L, Yang P, Zheng Y, et al. Identification of a prognostic immune signature for cervical cancer to predict survival and response to immune checkpoint inhibitors. Oncoimmunology. 2019;8(12):e1659094.
Jiang LHY, Shen Q, Ding S, Jiang W, Zhang W, Zhu X. Role of mast cells in gynecological neoplasms. Front Biosci (Landmark Ed). 2013;18:773–81.
Ferrandina GRF, Legge F, Gessi M, Salutari V, Distefano MG, Lauriola L, Zannoni GF, Martinelli E, Scambia G. Prognostic role of the ratio between cyclooxygenase-2 in tumor and stroma compartments in cervical cancer. Clin Cancer Res. 2004;10(9):3117–23.
Ferrandina G, Lauriola L, Zannoni GF, Distefano MG, Legge F, Salutari V, et al. Expression of cyclooxygenase-2 (COX-2) in tumour and stroma compartments in cervical cancer: clinical implications. Br J Cancer. 2002;87(10):1145–52.
This study was supported by funds from the National Natural Science Foundation of China No. 81503293 (X. J. Y.). The funder had no role in study design, analysis, interpretation of data, the writing of the manuscript and the decision to submit the article for publication.
Ethics approval and consent to participate
Ethical approval for this study was obtained from the Ethics Committee of the First Affiliated Hospital of Wenzhou Medical University, Clinical research ethics review (2019) No.(118). And written informed consent was obtained from all patients.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Cai, L., Hu, C., Yu, S. et al. Identification and validation of a six-gene signature associated with glycolysis to predict the prognosis of patients with cervical cancer. BMC Cancer 20, 1133 (2020). https://doi.org/10.1186/s12885-020-07598-3