Skip to main content
  • Primary research
  • Open access
  • Published:

Identification of novel cell glycolysis related gene signature predicting survival in patients with endometrial cancer

Abstract

Background

Endometrial cancer (EC) is one of the three major gynecological malignancies. Numerous biomarkers that may be associated with survival and prognosis have been identified through database mining in previous studies. However, the predictive ability of single-gene biomarkers is not sufficiently specific. Genetic signatures may be an improved option for prediction. This study aimed to explore data from The Cancer Genome Atlas (TCGA) to identify a new genetic signature for predicting the prognosis of EC.

Methods

mRNA expression profiling was performed in a group of patients with EC (n = 548) from TCGA. Gene set enrichment analysis was performed to identify gene sets that were significantly different between EC tissues and normal tissues. Cox proportional hazards regression models were used to identify genes significantly associated with overall survival. Quantitative real-time-PCR was used to verify the reliability of the expression of selected mRNAs. Subsequent multivariate Cox regression analysis was used to establish a prognostic risk parameter formula. Kaplan–Meier survival estimates and the log‐rank test were used to validate the significance of risk parameters for prognosis prediction.

Result

Nine genes associated with glycolysis (CLDN9, B4GALT1, GMPPB, B4GALT4, AK4, CHST6, PC, GPC1, and SRD5A3) were found to be significantly related to overall survival. The results of mRNA expression analysis by PCR were consistent with those of bioinformatics analysis. Based on the nine-gene signature, the 548 patients with EC were divided into high/low-risk subgroups. The prognostic ability of the nine-gene signature was not affected by other factors.

Conclusion

A nine-gene signature associated with cellular glycolysis for predicting the survival of patients with EC was developed. The findings provide insight into the mechanisms of cellular glycolysis and identification of patients with poor prognosis in EC.

Background

Endometrial cancer (EC) is one of the three major gynecological malignancies and fifth most common cancer among women (4.8% of female cancer cases) in the United States [1]. It is expected that 61,880 new cases will be diagnosed in 2019 (7% of all female cancer cases) and 12,160 deaths (4% of all female cancer deaths). In the past ten years, with the irregular use of hormones and changes in people’s living environment and lifestyle, the prevalence and mortality of endometrial cancer in China and abroad have been increasing annually [2]. Although most patients are diagnosed early, approximately 28% of patients are diagnosed with advanced disease. However, patients with the same degree of progression can show different prognoses and treatment responses. Therefore, effective EC biomarkers must be discovered for assessing prognosis and identifying potentially patients at a high risk of EC.

Numerous biomarkers for EC have been identified, such as the SIX1 and HER-2 genes [3, 4]. With advancements in high-throughput sequencing, researchers have established various patient genome databases to enable a more systematic understanding of genomic changes. Through database mining, we identified thousands of biomarkers that may be associated with the prognosis of patients with tumors [5, 6]. However, the predictive ability of single-gene biomarkers remains insufficient. Studies have shown that the evaluation of genetic traits, which involve multiple genes, may improve prognosis prediction [7, 8]. Multigenic prognostic features from primary tumor biopsy can guide more specific treatment strategies. Recent studies have explored the effects of multiple-gene signature on EC for assessing prognosis and identifying potentially patients at a high risk of EC [9, 10].

In this study, genes were selected by performing gene set enrichment analysis (GSEA). To identify biomarkers, differential analysis typically involves comparison of the expression differences between groups and focuses on genes whose expression levels are significantly regulated. However, this method can easily exclude genes that do not show obvious expression differences but may provide important biological information or exhibit biological significance. As an emerging computational method, GSEA does not require a clear differential gene threshold or extensive experience to test the overall expression of several genes. It reveals general trends in the data. Therefore, this approach improves the statistical analysis between biological expression and biological significance [11].

Accordingly, in the present study, we aimed to explore data from The Cancer Genome Atlas (TCGA) to identify a new genetic signature for predicting the prognosis of EC. To this end, we used mRNA expression data from TCGA to map the marker genome of 548 patients with EC. We identified 119 mRNAs significantly related to glycolysis and developed a nine-gene risk profile for effectively predicting patient outcomes. Interestingly, the risk factors associated with glycolysis can be used assess prognosis of high-risk patients independently. A novel cell glycolysis-related gene signature was identified and validated.

Methods

Clinical information and mRNA expression data set of patients

We extracted clinical data and the mRNA expression profiles of patients with endometrial cancer from TCGA (https://cancergenome.nih.gov/) [12]. The study included clinical information from 548 patients and enrolled matching age, stage, grade, radiation therapy, residual tumor, histological type, diabetes, new tumor events, and hypertension (Table 1).

Table 1 Clinical pathological parameters of patients with Endometrioid cancer in this study

Gene set enrichment analysis

We performed GSEA (http://www.broadinstitute.org/gsea/index.jsp) to determine if the identified gene sets were significantly different between the EC and normal groups. Next, we analyzed the expression levels of 24,991 mRNAs in EC samples and in adjacent noncancerous tissues. Finally, we determined functions for subsequent analysis by using normalized p values (p < 0.05).

Data processing and risk-parameter calculation

Log2 transformation was used to normalize each mRNA from among the expression profiles. Univariate Cox regression analysis was used to identify genes associated with overall survival (OS), which were then subjected to multivariable Cox regression to confirm the genes related to prognosis and obtain the coefficients. The selected mRNAs were then divided into the risky (hazard ratio, HR > 1) type and protective (0 < HR < 1) type. By linearly combining the expression values of filtered genes weighted by their coefficients, we constructed a risk-parameter formula as follows: Risk parameter = ∑ (βn × expression of gene n). Using the median risk parameter as a cut-off, the 548 patients were divided into high‐risk and low‐risk subgroups.

Specimens and patients of quantitative real-time (qRT)-PCR

A total of 20 EC tissues and 20 normal endometrial tissues were obtained from patients at the Department of Gynecology and Obstetrics, Shengjing Hospital of China Medical University, China. Normal tissues were obtained from patients who underwent hysterectomy for endometrial-irrelevant diseases. All patients provided informed consent, and this study was approved by the Ethics Committee of Shengjing Hospital of China Medical University. Histological diagnosis and grade were assessed by experienced pathologists in accordance with the FIGO 2009. No patient was administered systemic treatment preoperatively.

RNA extraction and qRT-PCR

Total RNA was extracted from tissues using TRIzol reagent (Vazyme, Nanjing, China). PrimeScript RT-polymerase (Vazyme) was used to reverse-transcribe cDNAs corresponding to the mRNAs of interest. qRT-PCR was performed using SYBR-Green Premix (Vazyme) with specific PCR primers (Sangon Biotech Co., Ltd, Shanghai, China). Glyceraldehyde-3-phosphate dehydrogenase was used as an internal control. The 2−ΔΔCt method was used to calculate fold-changes. Primer sequences are listed in Additional file 1: Table S1.

Statistical analysis

We used Kaplan–Meier survival curves and the log‐rank method to estimate the significance of the risk parameter. We performed multivariate Cox analysis and data stratification analysis to test whether the risk parameter was independent of the clinical features, including age, grade, stage, new event, residual tumor, and neoplasm cancer status, which were used as covariates. A p < 0.05 was considered as statistically significant. Statistical analyses were performed using GraphPad Prism7 software (GraphPad, Inc., La Jolla, CA, USA) and SPSS 20.0 software (SPSS, Inc., Chicago, IL, USA).

Results

Initial screening of genes by GSEA

We obtained the clinical features of 548 patients with EC, along with an expression data set for 24,991 mRNAs from TCGA database. We performed GSEA to determine if the identified gene sets were significantly different between EC tissues and endometrium tissues. We validated 26 gene sets that were upregulated in EC. Ten gene sets, G2 M checkpoints, MYC targets V1, glycolysis, MYC targets V2, MTORC1 signaling, oxidative phosphorylation, DNA repair, unfolded protein response, E2F targets, and UV response, were significantly enriched (Table 2, Fig. 1). We then filtered the top‐ranking function, glycolysis (p = 0.000), among 119 genes for subsequent analysis.

Table 2 Gene sets enriched in Endometrial cancer (548 samples)
Fig. 1
figure 1

Enrichment plots of nine gene sets which had significant difference between noncancerous tissues and EC tissues by performing GSEA

Identification of survival-associating glycolysis‐related mRNAs

First, we performed univariate Cox regression analysis of 119 genes for preliminary screening and obtained 21 genes (p < 0.05). Next, multivariate Cox regression analysis was performed to further examine the association between the 21 mRNA expression profiles and patient survival, and the stepwise elimination method was used to identify the most significant mRNAs combinations. Nine mRNAs (CLDN9, B4GALT1, GMPPB, B4GALT4, AK4, CHST6, PC, GPC1, and SRD5A3) were verified, as shown in Table 3, and six of the nine genes (CLDN9, B4GALT1, GMPPB, AK4, PC, and SRD5A3) were validated as independent prognostic markers of EC. The filtered mRNAs were divided into the risky type (CLDN9, AK4, PC, GPC1, and SRD5A3), with HR > 1 associated with poorer survival and the protective type (B4GALT1, GMPPB, B4GALT4, and CHST6), with HR < 1 associated with better survival (Table 3).

Table 3 The detailed information of nine prognostic mRNAs significantly associated with overall survival in patients with endometrial cancer

We then assessed the alterations in nine filtered genes by analyzing 548 EC samples in the cBioPortal database (http://cbioportal.org) [13]. The results showed that the queried genes were altered in 98 (18.3%) of the sequenced cases. The PC gene included ten amplification samples, 2 deep deletion samples, 7 missense mutations samples, and one sample with an in-frame mutation. The CLDN9 gene was altered in 2% of cases, showing various changes. The CHST6 gene was altered in 2.8% of cases, and the B4GALT1 and B4GALT4 genes were altered in 2.4% and 1.8% of cases, respectively (Fig. 2a).

Fig. 2
figure 2

Identification of mRNAs related to patients’ survival. a Selected genes’ alteration in 548 clinical samples. b Selected genes’ specific alteration in different pathological types of EC. c Different expression of nine selected genes

The specific alterations in the selected genes were significant in specific cancer types. In EC, 3.88% of changes were mutations, 7.17% were amplifications, 2.67% were deep deletions, and 0.16% were multiple alterations. In endometrioid adenocarcinoma, the only alteration was amplification (12.76% of patients). In mixed serous and endometrioid cases, the most eminent alteration was mutation (Fig. 2b).

Comparison of the expression differences of 9 genes between adjacent normal tissues and EC tissues was also performed. We found that the expression levels of the nine genes were significantly up- or downregulated in EC tissues (Fig. 2c).

Validation of TCGA expression results using qRT-PCR

We examined the expression of 9 mRNAs by qRT-PCR in 20 EC tissues and 20 normal endometrial tissues. We applied the unpaired t test to assess the differences between the two groups. The results showed that CLDN9, AK4, PC, GPC1, and SRD5A3 were upregulated in EC tissues compared to in normal endometrium tissues, whereas B4GALT1, GMPPB, B4GALT4, and CHST6 were downregulated in tumor tissues (Fig. 3). The mRNA expression results of qRT-PCR validation in 20 patients with EC were consistent with the bioinformatics results, which showed that the bioinformatics analysis was precise and gave significant results.

Fig. 3
figure 3

Expression of nine mRNAs in endometrial cancer tissues and normal tissues. a CLDN9, b B4GALT1, c GMPPB, d B4GALT4, e AK4, f CHST6, g PC, h GPC1, i SRD5A3

Construction of a nine‐mRNA signature to predict patient outcomes

By linearly combining the expression values of selected genes weighted by their coefficients derived from multivariate Cox regression analysis, we established the following prognostic risk-parameter formula. Risk parameter = 0.1059 × expression of CLDN9 + 0.3181 × expression of AK4 + 0.328 × expression of PC + 0.2056 × expression of GPC1 + 0.2345 × expression of SRD5A3 − 0.2504 × expression of B4GALT1 − 0.3041 × expression of B4GALT4 − 0.4346 × expression of GMPPB − 0.1191 × expression of CHST6. We calculated parameters and assigned one risk parameter to each patient. We then ranked the patients in ascending order by the parameter and divided the patients into high-risk and low-risk subgroups using the median (Fig. 4a). The survival time of each patient is shown in Fig. 4b. Patients in the high‐risk parameter group showed poorer survival, whereas patients with the low‐risk parameter had lower mortality rates. Additionally, a heatmap displayed the expression profiles of nine mRNAs (Fig. 4c). Compared to the low-risk group, the expression level of risky-type mRNA (CLDN9, AK4, PC, GPC1, and SRD5A3) was higher in the high-risk group. In contrast, the expression level of protective-type mRNA (CLDN9, AK4, PC, GPC1, and SRD5A3) in the high-risk group was lower than that in the low-risk group. The expression levels of the nine genes in the high and low risk groups are shown in Additional file 2: Figure S1.

Fig. 4
figure 4

The nine‐mRNA signature associated with risk parameter predicts OS in patients with endometrial cancer. a mRNA risk parameter distribution in each patient. b Survival days of EC patients in ascending order of risk parameters. c A heatmap of nine genes’ expression profile. d GO analysis and KEGG analysis of nine differentially expressed mRNAs

We performed Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses of differentially expressed mRNAs, including KEGG pathway enrichment analysis and GO functional annotation analysis of the 9 differentially expressed genes in the high- and low-risk groups (Fig. 4d). KEGG pathway enrichment analysis showed that the differentially expressed genes were involved in multiple signaling pathways and enriched in neuroactive ligand-receptor interaction, protein digestion and absorption, regulation of lipolysis in adipocytes, and another 12 pathways. The two groups of differentially expressed genes were annotated with GO functions and analyzed for their significant functions. The results revealed 15 GO items in the significantly enriched molecular function category, including ion channel activity, substance-specific channel activity, and passive transmembrane transporter activity, among which passive transmembrane transporter activity was the most significantly enriched GO term.

Risk parameter derived from nine-mRNA signature is an independent prognostic indicator

We compared the prognostic value of risk parameters with clinical pathology parameters by univariate and multivariate analyses (Table 1). Samples with well-established clinical data were selected. The median age of the 546 patients with EC was 66 years. Among the 548 patients, 63 (11.5%) had a new event during the follow‐up, 459 (75%) had endometrioid adenocarcinoma, and 31 (5.7%) suffered from hypertension. Among the 398 patients, 22 (5.5%) had residual tumors, and 99 (18.4%) had grade 1 tumors. Among 537 patients, 99(18.4%)had grade 1 tumors, 122 (22.7%) had grade 2 tumors, and the remaining 316 (58.6%) had grade 3 tumors. Furthermore, among the 548 patients with EC, 341 (62.2%) patients had stage I disease, and the remaining 207 (37.8) patients had stage II–IV disease. From the data set above, we identified the risk parameter, stage, grade, neoplasm cancer status, and new tumor event as independent prognostic indicators, as these factors showed significant differences in both univariate and multivariate analyses (Table 4). Notably, the risk parameters showed significant prognostic values with p < 0.05 (HR = 1.783).

Table 4 Univariable and multivariable analyses for each clinical feature

Verification of nine‐mRNA signature for prognosis prediction by K–M survival estimates

K–M survival estimates and the log‐rank test revealed that patients in the high-risk group had a poor prognosis (Fig. 5a). Univariate Cox regression analysis of OS revealed several clinicopathological parameters which were predictive of EC survival, including age, grade, stage, history of cancer, residual tumor, tumor recurrence (new events), and histological type. We then used the Kaplan–Meier survival estimates to validate the above conclusions, which showed consistent results, Patients older than 66 years, with tumor grade G2–3, disease stage II–IV, histological type of serous adenocarcinoma or mixed serous/endometrioid, and suffering from tumor recurrence (new events) and residual tumors were associated with poor prognosis (Fig. 5b, c). These results further confirmed the reliability of the analysis.

Fig. 5
figure 5

Kaplan–Meier survival analysis for EC patients in TCGA data set. a K–M survival curve for EC patients with high/low risk. b Clinical features including age, grade, stage, neoplasm cancer status, residual tumor and new event predict patients survival. c histological type predict patients survival

After further data mining and stratified analysis, the survival curves were not affected by stage (stage I or stage II–IV), and the nine-mRNA signature was validated as a reliable prognostic indicator for patients with EC. Patients in the high-risk group showed a poor prognosis (Fig. 6b). Similarly, despite the tumor grade, neoplasm cancer status, and tumor recurrence, the risk parameter based on the nine‐mRNA signature could be used to predict the prognosis of patients with EC (Fig. 6c–e). However, when we stratified patients with EC into two subgroups according to age (> 66 or ≤ 66 years), the risk parameter could no longer be independently used as an prognostic indicator for the subgroup of age < 66 years (Fig. 6a), indicating that the risk parameter is affected by the age of patients with EC; this point requires further exploration.

Fig. 6
figure 6

Kaplan–Meier curves for prognostic value of risk parameter signature for the patients divided by each clinical feature. a Age, b stage, c grade, d Person neoplasm cancer status, e new event

Discussion

Recent studies showed that clinicopathological features such as age and metastatic diagnosis are not sufficient to precisely predict the outcome of patients with cancer. Thus, an increasing number of mRNAs have been identified as biomarkers of tumor progression or prognosis, and the clinical significance of the biomarkers has been evaluated [14]. For example, Nadaraja et al. [15] confirmed that low expression of ARAP1 is an independent prognostic biomarker of shorter progression-free survival in older patients with ovarian high-grade serous adenocarcinoma being administered first-line platinum-based antineoplastic therapy. Similarly, multivariate Cox proportional hazards regression model analysis was used to verified that patients with cervical cancer who had high tumor protein p73 expression had better outcomes, and thus this protein was considered as a prognostic indicator in patients with cervical cancer [16]. However, these biomarkers were still not sufficient for independently predicting patient prognosis. Particularly, single gene expression levels can be affected by multiple factors, preventing these markers from being used as reliable and independent prognosis indicators. Thus, a statistical model comprised of genetic markers for multiple related genes, combined with the predictive effect of each constituent gene, was used to improve prediction. The model is significantly more accurate than using single biomarkers in assessing the prognosis of patients with tumors [17, 18], leading to widespread use of the model.

The rapid development of high-throughput genetic sequencing technology has established a foundation for large biological data research [19]. Large amounts of genomic data were extracted from individual specimens to identify new diagnostic, prognostic, or pharmacological biomarkers [20]. In recent studies, a new prognostic signature was constructed by using microarray and RNA-sequencing data for gene expression levels or mutations. A Cox proportional hazards regression model was used for identification and verification [21, 22]. In the current study, we identified 10 functions showing significant differences in GSEA. As described above, rather than wide-range exploration, we selected the top‐ranking function to filter genes related to patient survival prediction. Univariate and multivariate Cox regression analyses were performed to determine the prognostic value of the combination of nine genes for patients with EC. This selected risk profile may be a more targeted and powerful prognostic assessment for predicting positive clinical outcomes and may be a more effective classification tool for patients with EC compared to other known prognostic assessment markers.

In this study, bioinformatics methods were used to explore the characteristics of mRNA risk factors and their clinical significance, and a new method for mining of potential prognostic markers was explored. This study y complements the previous understanding of EC and provides a foundation for future EC research. We used the EC dataset in TCGA to collect glycolysis-related genes and compare data from normal and EC tissues. Kaplan–Meier survival estimates revealed that patients with low-risk parameters had a better prognosis. The detection and calculation of risk parameter in EC patients have important clinical implications. However, because of the lack of patient metastasis and recurrence information in TCGA database, we could only use OS to assess patient prognosis, which is one limitation of our research. Additionally, in stratified analysis, the risk parameter could predict the prognosis of patients with EC in all subgroups except for the subgroup of age < 66 years. The reason for this difference is unclear requires further examination.

In addition, the nine-gene signature and same analysis method in liver cancer and colon cancer were used to obtain and verify the corresponding risk parameter (Additional file 3: Figure S2; Additional file 4: Figure S3). The results showed that the risk parameter based on the nine genes is not an independent prognostic indicator for liver cancer and colon cancer, confirming that the nine-gene signature is particularly important in EC.

Tumors are characterized by uncontrolled cell proliferation, which not only eliminates control of the cell cycle but also promotes cellular energy metabolism and finally leads to tumor cell growth and differentiation. Cellular energy is mainly derived from sugar metabolism, and most energy is supplied by ATP. In the 1920s, the German biologist Otto Warburg discovered abnormalities in energy metabolism in hepatoma cells. Although oxygen is present, tumor cells mainly rely on glycolysis for metabolism and consume large amounts of glucose accompanied by lactic acid production. This phenomenon of abnormal glucose metabolism was named as aerobic glycolysis or the Warburg effect [23]. Studies have shown that tumor cells can precisely regulate ATP synthesis by regulating substrate uptake and enzymes related to glycolysis, enabling them to adapt to the nutrient microenvironment, meet the energy and nutrient requirements for malignant proliferation, rapidly proliferate. Moreover, cancer metabolic reprogramming, which is closely associated with the Warburg effect, plays an important role in maintaining the interaction between oxygen-sensing transcription factors and the nutrient-sensing signal pathway [24]. This indicates that aerobic glycolysis uses a complicated mechanism of action. Tumor cell proliferation proceeds at a pace exceeding cellular energy supply, and thus excessive consumption of oxygen and nutrients by the cells can cause the tumor microenvironment to be hypoxic, low in sugar, and acidic, which is more pronounced in solid tumors [25]. Although not all tumors exhibit the Warburg effect, cellular energy abnormalities are widely recognized as one of the characteristics of tumor cells. After more than 90 years of continuous exploration and research, the Warburg effect has been found to occur in many malignancies, such as lung cancer, breast cancer, colon cancer, and gastric cancer. Recent studies showed that aerobic glycolysis plays an important role in EC occurrence and development. Metabolic profiling of EC cells revealed higher rates of glycolysis and lower glucose oxidation, and tumor cells may rely on GLUT6-mediated glucose transport and glycolytic–lipogenic metabolism for survival [26]. Highly differentiated EC showed significantly lower GLUT1 and GLUT3 expression than poorly differentiated tumors [27]. Several studies have predicted the survival of patients with EC using genes associated with cellular glycolysis. For example, high mobility group protein 1 suppression effectively inhibits the development and progression of EC [28]. The expression of lactate dehydrogenase 5 in EC is an independent prognostic indicator strongly associated with poor prognosis [29]. However, glycolysis-related gene markers for predicting EC prognosis have not been established. Using bioinformatics methods, we determined the genetic characteristics associated with cellular aerobic glycolysis (CLDN9, B4GALT1, GMPPB, B4GALT4, AK4, CHST6, PC, GPC1, and SRD5A3) and demonstrated their prognostic value in EC.

Conclusion

We developed a nine-gene risk profile associated with cellular glycolysis which predicts the prognosis of patients with EC, with a higher risk parameter indicating poorer prognosis. The signature can be a used as a classification tool in clinical practice. These findings provide insight into the mechanisms of cellular glycolysis and identification of patients with poor prognosis in EC.

Availability of data and materials

The datasets generated and analysed during the current study are available in the TCGA (http://cancergenome.nih.gov/abouttcga) and cBioPortal (http://www.cbioportal.org).

Abbreviations

EC:

endometrial cancer

TCGA:

The Cancer Genome Atlas

GSEA:

gene set enrichment analysis

OS:

overall survival

MSigDB:

Molecular Signatures Database

FIGO:

International Federation of Gynecology and Obstetrics

SIX1:

oncofetal protein sine oculis-related homeobox 1

HER-2:

human epidermal growth factor receptor-2

ARAP1:

Ankyrin repeat and PH domain 1

GLUT1:

glucose transporter type 1

GLUT3:

glucose transporter type 3

References

  1. Morice P, Leary A, Creutzberg C, Abu-Rustum N, Darai E. Endometrial cancer. Lancet. 2016;387(10023):1094–108.

    Article  Google Scholar 

  2. Torre LA, Islami F, Siegel RL, Ward EM, Jemal A. Global cancer in women: burden and trends. Cancer Epidemiol Biomarkers Prev. 2017;26(4):444–57.

    Article  Google Scholar 

  3. Suen AA, Jefferson WN, Wood CE, Padilla-Banks E, Bae-Jump VL, Williams CJ. SIX1 oncoprotein as a biomarker in a model of hormonal carcinogenesis and in human endometrial cancer. Mol Cancer Res. 2016;14(9):849–58.

    Article  CAS  Google Scholar 

  4. Shulman LP. HER-2 is an independent prognostic factor in endometrial cancer: association with outcome in a large cohort of surgically staged patients. Yearb Obstetric Gynecol Women’s Health. 2007;2007:404–5.

    Article  Google Scholar 

  5. Liu S, Miao C, Liu J, Wang CC, Lu XJ. Four differentially methylated gene pairs to predict the prognosis for early stage hepatocellular carcinoma patients. J Cell Physiol. 2018;233(9):6583–90.

    Article  CAS  Google Scholar 

  6. Qixing M, Gaochao D, Wenjie X, Anpeng W, Bing C, Weidong M, Lin X, Feng J. Microarray analyses reveal genes related to progression and prognosis of esophageal squamous cell carcinoma. Oncotarget. 2017;8(45):78838–50.

    Article  Google Scholar 

  7. Cheng W, Ren X, Zhang C, Cai J, Liu Y, Han S, Wu A. Bioinformatic profiling identifies an immune-related risk signature for glioblastoma. Neurology. 2016;86(24):2226–34.

    Article  CAS  Google Scholar 

  8. Chen YL, Ge GJ, Qi C, Wang H, Wang HL, Li LY, Li GH, Xia LQ. A five-gene signature may predict sunitinib sensitivity and serve as prognostic biomarkers for renal cell carcinoma. J Cell Physiol. 2018;233(10):6649–60.

    Article  CAS  Google Scholar 

  9. Liu J, Zhou S, Li S, Jiang Y, Wan Y, Ma X, Cheng W. Eleven genes associated with progression and prognosis of endometrial cancer (EC) identified by comprehensive bioinformatics analysis. Cancer Cell Int. 2019;19:136.

    Article  Google Scholar 

  10. Wang Y, Ren F, Chen P, Liu S, Song Z, Ma X. Identification of a six-gene signature with prognostic value for patients with endometrial carcinoma. Cancer Med. 2018;7(11):5632–42.

    Article  CAS  Google Scholar 

  11. Thomas MA, Yang L, Carter BJ, Klaper RD. Gene set enrichment analysis of microarray data from Pimephales promelas (Rafinesque), a non-mammalian model organism. BMC Genomics. 2011;12:66.

    Article  CAS  Google Scholar 

  12. The Cancer Genome Atlas. https://tcga-data.nci.nih.gov/tcga/. Accessed 14 Feb 2019.

  13. cBioPortal for Cancer Genomics. http://cbioportal.org. Accessed 19 Mar 2019.

  14. Guo W, Chen X, Zhu L, Wang Q. A six-mRNA signature model for the prognosis of head and neck squamous cell carcinoma. Oncotarget. 2017;8(55):94528–38.

    PubMed  PubMed Central  Google Scholar 

  15. Nadaraja S, Schledermann D, Herrstedt J, Ostrup O, Ditzel HJ, Academy of Geriatric Cancer R. ARAP1 is an independent prognostic biomarker in older women with ovarian high-grade serous adenocarcinoma receiving first-line platinum-based antineoplastic therapy. Acta Oncol. 2019. https://doi.org/10.1080/0284186X.2019.1657941.

    Article  PubMed  Google Scholar 

  16. Ye H, Guo X. TP73 is a credible biomarker for predicting clinical progression and prognosis in cervical cancer patients. Biosci Rep. 2019. https://doi.org/10.1042/BSR20190095.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Bao ZS, Li MY, Wang JY, Zhang CB, Wang HJ, Yan W, Liu YW, Zhang W, Chen L, Jiang T. Prognostic value of a nine-gene signature in glioma patients based on mRNA expression profiling. CNS Neurosci Ther. 2014;20(2):112–8.

    Article  CAS  Google Scholar 

  18. Cheng W, Ren X, Cai J, Zhang C, Li M, Wang K, Liu Y, Han S, Wu A. A five-miRNA signature with prognostic and predictive value for MGMT promoter-methylated glioblastoma patients. Oncotarget. 2015;6(30):29285–95.

    PubMed  PubMed Central  Google Scholar 

  19. Peng PL, Zhou XY, Yi GD, Chen PF, Wang F, Dong WG. Identification of a novel gene pairs signature in the prognosis of gastric cancer. Cancer Med. 2018;7(2):344–50.

    Article  CAS  Google Scholar 

  20. Subramanian J, Simon R. What should physicians look for in evaluating prognostic gene-expression signatures? Nat Rev Clin Oncol. 2010;7(6):327–34.

    Article  Google Scholar 

  21. Shen S, Bai J, Wei Y, Wang G, Li Q, Zhang R, Duan W, Yang S, Du M, Zhao Y, et al. A seven-gene prognostic signature for rapid determination of head and neck squamous cell carcinoma survival. Oncol Rep. 2017;38(6):3403–11.

    CAS  PubMed  PubMed Central  Google Scholar 

  22. Zhao Y, Varn FS, Cai G, Xiao F, Amos CI, Cheng C. A P53-deficiency gene signature predicts recurrence risk of patients with early-stage lung adenocarcinoma. Cancer Epidemiol Biomarkers Prev. 2018;27(1):86–95.

    Article  CAS  Google Scholar 

  23. Koppenol WH, Bounds PL, Dang CV. Otto Warburg’s contributions to current concepts of cancer metabolism. Nat Rev Cancer. 2011;11(5):325–37.

    Article  CAS  Google Scholar 

  24. Yoshida GJ. Metabolic reprogramming: the emerging concept and associated therapeutic strategies. J Exp Clin Cancer Res. 2015;34:111.

    Article  Google Scholar 

  25. Ganapathy-Kanniappan S, Geschwind JF. Tumor glycolysis as a target for cancer therapy: progress and prospects. Mol Cancer. 2013;12:152.

    Article  Google Scholar 

  26. Byrne FL, Poon IK, Modesitt SC, Tomsig JL, Chow JD, Healy ME, Baker WD, Atkins KA, Lancaster JM, Marchion DC, et al. Metabolic vulnerabilities in endometrial cancer. Cancer Res. 2014;74(20):5832–45.

    Article  CAS  Google Scholar 

  27. Han J, Zhang L, Guo H, Wysham WZ, Roque DR, Willson AK, Sheng X, Zhou C, Bae-Jump VL. Glucose promotes cell proliferation, glucose uptake and invasion in endometrial cancer cells via AMPK/mTOR/S6 and MAPK signaling. Gynecol Oncol. 2015;138(3):668–75.

    Article  CAS  Google Scholar 

  28. Luan X, Ma C, Wang P, Lou F. HMGB1 is negatively correlated with the development of endometrial carcinoma and prevents cancer cell invasion and metastasis by inhibiting the process of epithelial-to-mesenchymal transition. Onco Targets Ther. 2017;10:1389–402.

    Article  CAS  Google Scholar 

  29. Giatromanolaki A, Sivridis E, Gatter KC, Turley H, Harris AL, Koukourakis MI. Tumour, angiogenesis research G: lactate dehydrogenase 5 (LDH-5) expression in endometrial cancer relates to the activated VEGF/VEGFR2(KDR) pathway and prognosis. Gynecol Oncol. 2006;103(3):912–8.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by the National Natural Science Foundation of China (No. 81872123); University innovation team of Liaoning Province; Special Professor of Liaoning Province; “Major Special Construction Plan” for Discipline Construction of China Medical University in 2018 (No. 3110118029); Outstanding Scientific Fund of Shengjing Hospital (No. 201601).

Author information

Authors and Affiliations

Authors

Contributions

ZHW participated in the design of the work, methodology, data interpretation, and analysis for the work; carried out the statistical analyses; and drafted the manuscript. YZZ carried out the statistical analyses. YSW participated in the methodology, data interpretation, and analysis for the work. XXM designed the study; participated in data interpretation, analysis for the work, and methodology. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Xiao-Xin Ma.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Ethics Committee of Shengjing Hospital of China Medical University.

Consent for publication

All listed authors have actively participated in the study and have read and approved the submitted manuscript.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1.

Primer sequence.

Additional file 2: Figure S1.

Expression of nine genes in high and low risk groups (* represents for p < 0.01, ** represent for p < 0.001, **** represent for p < 0.00001).

Additional file 3: Figure S2.

(a) The Kaplan–Meier curve for patients divided into high risk and low risk in liver cancer. 45(b) Different expression of nine selected genes in liver cancer.

Additional file 4: Figure S3.

a) The Kaplan–Meier curve for patients divided into high risk and low risk in colon cancer. (b) Different expression of nine selected genes in colon cancer.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wang, ZH., Zhang, YZ., Wang, YS. et al. Identification of novel cell glycolysis related gene signature predicting survival in patients with endometrial cancer. Cancer Cell Int 19, 296 (2019). https://doi.org/10.1186/s12935-019-1001-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12935-019-1001-0

Keywords