Jun Yang
Yuening Zhang2†
Haibin Yu- 1Phase Ⅰ Clinical Trial Ward, Beijing Youan Hospital, Capital Medical University, Beijing, China
- 2Beijing Youan Hospital, Capital Medical University, Beijing, China
Aim: To develop a glycolysis-related gene signature that correlated with the characteristics of the tumor immune microenvironment and had good predictive power for overall survival (OS) in hepatocellular carcinoma (HCC).
Methods: Gene expression profiles, RNA sequencing data, clinical characteristics and survival information for 407 patients with HCC and 58 healthy controls were downloaded from the TCGA database. GSEA 4.1.0 software was used to evaluate the glycolysis-related pathways enriched in HCC compared to normal liver tissue. Univariate Cox, Least Absolute Shrinkage, Selection Operator, and two-step multivariate Cox analyses were used to construct a glycolysis-related gene signature for prognostic prediction. The glycolysis-related gene signature was combined with clinical characteristics to generate a nomogram. Tumor-infiltrating immune cell profiles and PD-L1 protein expression in HCC tissues were investigated.
Results: The gene expression profiles of HCC tissues were enriched in glycolysis-related pathways. A glycolysis-related gene signature was used to categorize patients as high-risk or low-risk, where high-risk patients had significantly worse OS. Receiver operating characteristic curves confirmed the predictive capability of the glycolysis-related gene signature for OS (AUC >0.80). There was a significant difference in M0 macrophage (p = 0.017), dendritic cell (p = 0.043), B cell (p = 0.0018), CD4 T cell (p = 0.003), Treg (p = 0.01) and mast cell (p = 0.02) content and PD-L1 protein expression (p = 0.019) between HCC tissues in patients in the high-risk and low-risk groups.
Conclusion: We established a glycolysis-related gene signature for OS in HCC that was predictive in training and test TCGA cohorts and correlated with the characteristics of the HCC tumor immune microenvironment. The glycolysis-related gene signature may guide clinical decision-making concerning patient selection for immunotherapy in HCC.
Introduction
Hepatocellular carcinoma (HCC) is the most prevalent form of liver cancer. HCC is susceptible to metastasis and a frequent cause of cancer-related mortality (Siegel et al., 2020). Treatments for HCC include partial liver resection, interventional surgery, liver transplantation, and systemic treatment (Hua et al., 2018). These have benefit in early-stage patients, but the treatment of patients with advanced HCC is challenging (Forner et al., 2018).
Cancer progression is promoted by aerobic glycolysis, which is persistently activated in tumor cells (Ganapathy-Kanniappan, 2018). The liver is essential for regulating glucose metabolism (Kitamura et al., 2011), and HCC cells are characterized by increased glycolysis and lactic acid production (Liberti and Locasale, 2016). There are several reports describing a glycolysis-related gene-signature that can predict prognosis in patients with HCC (Zhou et al., 2020).
Novel immuno-oncology drugs are changing the outcomes of patients with HCC (Cheng et al., 2019; Anwanwan et al., 2020). Immune cells act as tumor suppressors or promoters and play a pivotal role in the tumor microenvironment (Lei et al., 2020). The characteristics of the tumor microenvironment greatly impact the response to immunotherapy (Zhang et al., 2019). There remains a critical unmet need to identify immune-related indicators that can predict response to immunotherapy and prognosis in patients with HCC.
To fill this evidence gap, we studied the correlation of clinical characteristics with glycolysis-related gene expression, explored multiomics differences in clinical characteristics and immune infiltration status, and constructed a prognostic signature related to glycolysis and immunity that has effective predictive value in patients with HCC.
Materials and Methods
Data Source
Gene expression profiles, RNA sequencing data, clinical characteristics and survival information for 407 patients with HCC and 58 healthy controls were extracted from The Cancer Genome Atlas (TCGA) database. Patients were randomly divided into training and test cohorts. Gene expression profiles were merged and normalized using the “sva” package (Leek et al., 2019) in R (v4.1.0).
Gene Set Enrichment Analysis of Glycolysis-Related Pathways
The Gene Set Enrichment Analysis (GSEA) website provided data on five glycolysis-related pathways. Glycolysis-related pathways enriched in HCC compared to normal liver tissue were identified using GSEA 4.1.0 software.
Establishment of a Prognostic Glycolysis-Related Gene Signature
The performance of glycolysis-related genes for predicting survival in patients with HCC was assessed with Least Absolute Shrinkage, Selection Operator (LASSO) and multivariate Cox regression analyses using the “glmnet” package in R.
A prognostic glycolysis-related gene signature was established and a signature-based risk score was calculated for each patient, as: score = esum (mRNA level X multivariate Cox regression coefficient ratio of each mRNA). Using the median risk score as the cutoff value, patients were classified as high-risk or low-risk. Overall survival (OS) and gene expression profiles for high-risk and low-risk patients were determined using the “survival” and “pheatmap” packages in R. Kaplan–Meier survival and receiver operating characteristic (ROC) curve analyses were used to estimate the sensitivity and specificity of the glycolysis-related gene signature.
Establishment of a Nomogram
A nomogram with independent predictors of prognosis was constructed. Goodness-of-fit of the nomgram was assessed with the area under the receiver operating characteristic curve and the Hosmer-Lemeshow test.
Validation of the Prognostic Glycolysis-Related Gene Signature
The cutoff value for high-risk vs. low-risk patients and differences in OS and clinicopathological factors, such as age, sex and tumor-node-metastasis (TNM) stage, between the high-risk and low-risk groups were validated using Kaplan–Meier survival and ROC curve analyses. Associations between the glycolysis-related gene signature and clinicopathological parameters were assessed using the χ2 or Wilcoxon signed-rank test. The potential for the glycolysis-related gene signature to represent an independent predictor of prognosis was explored using univariate and multivariate Cox regression analyses. These operations used the survival, glmnet, pbapply, survival ROC, survminer, pHeatmap, and ggupbr packages in R.
Tumor-Infiltrating Immune Cell Estimation
The association of the glycolysis-related genes with the HCC tumor microenvironment was explored by constructing a violin plot using the Fluidigm Singular Analysis Toolset 3.5.2 R package. This showed the contribution of immune cell infiltration in HCC tissues in patients in the high-risk and low-risk groups. Tumor-infiltrating immune cell (TIIC) profiles were assessed with scatter plots and Spearman correlations using the Tumor-Immune Estimation Resource (Li et al., 2017). TIICs in HCC tissues in patients in the high-risk and low-risk groups were compared using the limma and ggpubr packages in R.
PD-L1 Protein Expression in the HPA Database
PD-L1 protein expression was investigated in healthy liver tissues and liver cancer tissues from the HPA dataset (http://www.proteinatlas.org/) by immune histochemical staining.
Results
Study Design and Summary of Patients’ Information
The study process is shown in Figure 1. A total of 465 (n = 407 HCC tissues; n = 58 healthy liver tissues) gene expression profiles were obtained from the TCGA database (Table 1). An mRNA volcano map and an mRNA heatmap are shown in Figure 2.
FIGURE 1. Study process.
TABLE 1. Baseline data of all HCC samples.
FIGURE 2. Gene expression profiles of patients with HCC (A) Heatmap; (B) mRNA volcano map.
GSEA of Glycolysis-Related Pathways
HCC tissues were enriched in glycolysis-related genes (Figure 3). The top 17 Kyoto Encyclopedia of Genes and Genomes pathways are shown in Figure 3A.
FIGURE 3. Glycolysis-related gene signature (A) GSEA of 17 pathways; (B–F) KEGG pathways.
Construction of a Prognostic Glycolysis-Related Gene Signature
The expression profiles of the glycolysis-related genes were incorporated into a prognostic indicator using LASSO Cox regression analysis (Figures 4D,E). A glycolysis-related gene signature was constructed. All patients in the TCGA dataset were classified as high-risk (n = 202) or low-risk (n = 201) based on their signature-based risk scores. OS was significantly longer in low-risk patients compared to high-risk patients (Figure 4A). Univariate and multivariate analyses indicated that the glycolysis-related gene signature was an independent predictor of OS and several clinicopathological factors, including age, sex and tumor stage, in HCC (Figures 4B,C).
FIGURE 4. Prognostic glycolysis-related gene signature (A) Kaplan–Meier survival analysis of the low-risk and high-risk groups; (B) Univariate Cox regression analysis of the risk score and clinical characteristics; (C) Multivariate Cox regression analysis of the risk score and clinical characteristics; (D) LASSO regression; (E) Optimal lambda values.
Validation of the Prognostic Glycolysis-Related Gene Signature
In the training cohort, OS was significantly longer in low-risk patients compared to high-risk patients (p < 0.001; Figure 5A). Calibration curves showed good consistency between the prediction based on the glycolysis-related gene signature and the actual observation (Figure 5B). The predictive ability of the glycolysis-related gene signature, determined by the area under the ROC curve (AUC = 0.828; Figure 5C), was good. The distribution of risk scores, patient survival status, and a heatmap of glycolysis-related gene expression are shown in Figure 5F.
FIGURE 5. Training set (A) Kaplan-Meier curve; (B) Calibration curve; (C) Time-dependent ROC curve analysis; (D) Heat map of glycolysis-related gene expression; (E) Distribution of risk scores; (F) Survival time.
The test cohort was used to validate the prognostic glycolysis-related gene signature. OS was significantly longer in low-risk patients compared to high-risk patients (p = 0.001; Figure 6A). ROC curves, the distribution of risk scores, patient survival status, and a heatmap of glycolysis-related gene expression are shown in Figures 6B–F.
FIGURE 6. Test set (A) Kaplan-Meier curve; (B) Calibration curve; (C) Time-dependent ROC curve analysis; (D) Heat map of glycolysis-related gene expression; (E) Distribution of risk scores; (F) Survival time.
Stratification of the Prognostic Glycolysis-Related Gene Signature
The glycolysis-related gene signature was significantly associated with tumor stage (p = 0.0022), the T category of TNM staging (p = 0.0025) and tumor grade (p = 0.00012), but not the M category of TNM staging (p = 0.38), N category of TNM staging (p = 0.69), sex (p = 0.51), or age (p = 0.79) (Figures 7A–I).
FIGURE 7. Stratification of the prognostic glycolysis-related gene signature. Stage of HCC: TMN, high:high risk, low:low risk.
Nomogram Building and Validation
Multivariate Cox regression analysis identified the risk score, age, sex, and tumor stage as independent prognostic factors in the TCGA cohort. A novel nomogram incorporating these prognostic factors was constructed. The risk score had the widest range of variables when calculating the total nomogram score, indicating that it played a prominent role in determining OS (Figure 8A). Calibration curves showed good consistency between the nomogram prediction and the actual observation (Figure 8C). The predictive ability of the nomogram, determined by the area under the ROC curve (AUC = 0.80; Figure 8D), was good. Decision curve analysis indicated that the nomogram had the best net benefit (Figure 8B).
FIGURE 8. Nomogram (A) Nomogram incorporating prognostic factors (B) Decision curve analysis; (C) Calibration curves showed good consistency between the nomogram prediction and the actual observation; (D) Time-dependent ROC curve analysis.
Tumor Infiltration Immune Cell Profiling
CIBERSORT showed a significant difference in M0 macrophage (p = 0.017), dendritic cell (p = 0.043), B cell (p = 0.0018), CD4 T cell (p = 0.003), Treg (p = 0.01) and mast cell (p = 0.02) content in HCC tissues from patients in the high-risk group compared to the low-risk group (Figures 9A–G).
FIGURE 9. Tumor infiltration immune cell profiling.
Verification of PD-L1 Protein Expression in the HPA Database
Immunohistochemical staining showed a significant difference in PD-L1 expression level between normal liver tissue and liver cancer tissue (Figure 10B; p = 0.001) from the HPA dataset (http://www.proteinatlas.org/) and HCC tissues from patients in the high-risk group compared to the low-risk group in the TCGA cohort (Figure 10C; p = 0.019).
FIGURE 10. PD-L1 protein expression (A) Risk model; (B) Normal tissue and tumor tissue; (C) Tumor tissues from patients in the low-risk and high-risk groups; (D,E) Tumor tissues derived from the HPA database.
Discussion
HCC has a complex pathogenesis and involves dysregulation of the cell cycle and energy metabolism (Hsu et al., 2015). mRNAs are involved in tumorigenesis and progression by regulating the rate-limiting enzymes involved in energy metabolism. Glycolysis is an inefficient form of energy generation, but it produces reducing equivalents (Terabe et al., 2019) and biosynthetic substrates required for tumor cell proliferation (Liang et al., 2019). The objective of this study was to develop a novel prognostic glycolysis-related gene signature that correlated with the characteristics of the tumor immune microenvironment in HCC. Gene expression profiles and clinical information were obtained from the TCGA database. Glycolysis-related genes were identified using univariate Cox, LASSO, and two-step multivariate Cox analyses. Differentially expressed glycolysis-related genes were combined to construct a glycolysis-related gene signature with value for predicting prognosis. The glycolysis-related gene-signature was combined with clinical characteristics to generate a nomogram with a high prediction performance. The glycolysis-related gene signature correlated with TIIC profiles and PD-L1 protein expression in HCC tissues. Previous reports have described glycolysis-related gene-signatures that can predict prognosis in patients with HCC; however, our glycolysis-related gene-signature has greater capability for predicting survival (Zhou et al., 2020).
Glycolysis-related genes may represent biomarkers of survival and prognosis in HCC, likely because altered energy metabolism and epigenetic processes are key components of HCC pathophysiology (Chen et al., 2018). Specifically, the Warburg effect, characterized by high glucose uptake and increased lactate production, is a phenotype commonly exhibited by HCC cells (Anderson et al., 2018; Pascale et al., 2020). The Warburg effect results from mitochondrial dysfunction (Riera Leal et al., 2020), tumor adaptation (Ždralević et al., 2018), microenvironmental changes (Sun et al., 2018), oncogenes (Banks, 2013) and abnormalities in related signaling pathways. In the present study, GSEA showed HCC tissues from patients at high risk for poor prognosis were enriched in glycolysis-related genes. Glycolysis-related pathways have been associated with recurrence in liver cancer (Zhang et al., 2018).
A previous study reported that long coding RNAs were involved in HCC immune evasion (Peng et al., 2020). Consistent with this, findings from the present study showed a significant difference in the immune cell profile of HCC tumors in the high-risk and low-risk groups, risk characteristics could be used as immune indicators for HCC, immune features were activated when the risk score increased, and the proportion of Tregs increased with an increase in the risk score. These data indicate Tregs play a role in the HCC tumor microenvironment, and patients with low-risk scores may respond to immunotherapy.
TIICs have been identified as a key component of the tumor microenvironment and response to immunotherapy. An immune score is usually used to guide immunotherapy and assess prognosis. Immune checkpoint inhibitors have changed the therapeutic landscape in cancer (Salik et al., 2020). PD-L1 blockade has provided new insights into the clinical management of HCC (Ng et al., 2020). mRNAs play a key role in immune regulation, through antigen presentation and immune cell infiltration (Denaro et al., 2019). In the present study, there was a significant difference in PD-L1 expression level in HCC tissues from patients in the high-risk group compared to the low-risk group, suggesting patients with high-risk scores may be candidates for immunotherapy targeting the PD-1/PD-L1 pathway.
This study was limited. First, it was a retrospective analysis using information extracted from a single database. we are going to examine the conclusion with our own patient cohort in future. Second, in Figure 7B, although p value <0.05, the Pearson R value are mostly lower than 0.3. According to the data, it is hard to drawn any conclusion about the relationship between glycolysis-related pathway and immune infiltration. We are going to study the question in other data in future. Third, we are going to examine the impacts of genes involved in glycolysis-related signature on tumor progression and immunotherapy efficacy by establishing specific gene knockout hepatoma-bearing mice in future.
Future research should be large scale, prospective, and involve multiple centers.
Conclusion
We established a glycolysis-related gene signature for OS in HCC that correlated with the characteristics of the HCC tumor immune microenvironment and was predictive in training and test TCGA cohorts. The glycolysis-related gene signature may guide clinical decision-making concerning patient selection for immunotherapy in HCC.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
Conception and design of the study (ZH), collection of the data (HY, XH), analysis of the data (JD), and writing of the paper (JY, YZ).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
The author thanks the platform for help provided by the TCGA database (https://portal.gdc.cancer.gov/).
References
Anderson, N. M., Mucka, P., Kern, J. G., and Feng, H. (2018). The Emerging Role and Targetability of the TCA Cycle in Cancer Metabolism. Protein Cell 9 (2), 216–237. doi:10.1007/s13238–017–0451–110.1007/s13238-017-0451-1
Anwanwan, D., Singh, S. K., Singh, S., Saikam, V., and Singh, R. (2020). Challenges in Liver Cancer and Possible Treatment Approaches. Biochim. Biophys. Acta (Bba) - Rev. Cancer 1873 (1), 188314. doi:10.1016/j.bbcan.2019.188314
Banks, R. E. (2013). Oncogene-induced Cellular Senescence Elicits an Anti-warburg Effect. Proteomics 13 (17), 2542–2543. doi:10.1002/pmic.201300335
Chen, R., Zhu, S., Fan, X.-G., Wang, H., Lotze, M. T., Zeh, H. J., et al. (2018). High Mobility Group Protein B1 Controls Liver Cancer Initiation through Yes-Associated Protein -dependent Aerobic Glycolysis. Hepatology 67 (5), 1823–1841. doi:10.1002/hep.29663
Cheng, H., Sun, G., Chen, H., Li, Y., Han, Z., Li, Y., et al. (2019). Trends in the Treatment of Advanced Hepatocellular Carcinoma: Immune Checkpoint Blockade Immunotherapy and Related Combination Therapies. Am. J. Cancer Res. 9 (8), 1536–1545.
Denaro, N., Merlano, M. C., and Lo Nigro, C. (2019). Long Noncoding RNA S as Regulators of Cancer Immunity. Mol. Oncol. 13 (1), 61–73. doi:10.1002/1878–0261.12413
Forner, A., Reig, M., and Bruix, J. (2018). Hepatocellular Carcinoma. The Lancet 391, 1301–1314. doi:10.1016/s0140-6736(18)30010-2
Ganapathy-Kanniappan, S. (2018). Molecular Intricacies of Aerobic Glycolysis in Cancer: Current Insights into the Classic Metabolic Phenotype. Crit. Rev. Biochem. Mol. Biol. 53 (6), 667–682. doi:10.1080/10409238.2018.1556578
Hsu, C.-C., Wu, L.-C., Hsia, C.-Y., Yin, P.-H., Chi, C.-W., Yeh, T.-S., et al. (2015). Energy Metabolism Determines the Sensitivity of Human Hepatocellular Carcinoma Cells to Mitochondrial Inhibitors and Biguanide Drugs. Oncol. Rep. 34 (3), 1620–1628. doi:10.3892/or.2015.4092
Hua, S., Liu, C., Liu, L., and Wu, D. (2018). miR-142–3p Inhibits Aerobic Glycolysis and Cell Proliferation in Hepatocellular Carcinoma by Targeting LDHA. Biochem. Biophys. Rescommun 496, 947–954. doi:10.1016/j.bbrc.2018.01.112
Kitamura, K., Hatano, E., Higashi, T., Narita, M., Seo, S., Nakamoto, Y., et al. (2011). Proliferative Activity in Hepatocellular Carcinoma Is Closely Correlated with Glucose Metabolism but Not Angiogenesis. J. Hepatol. 55, 846–857. doi:10.1016/j.jhep.2011.01.038
Leek, J. T., Johnson, W. E., Parker, H. S., Fertig, E. J., Jaffe, A. E., Story, J. D., et al. (2019). Sva: Surrogate Variable Analysis. R. Package Version
Lei, X., Lei, Y., Li, J.-K., Du, W.-X., Li, R.-G., Yang, J., et al. (2020). Immune Cells within the Tumor Microenvironment: Biological Functions and Roles in Cancer Immunotherapy. Cancer Lett. 470, 126–133. doi:10.1016/j.canlet.2019.11.009
Li, T., Fan, J., Wang, B., Traugh, N., Chen, Q., Liu, J. S., et al. (2017). TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells. Cancer Res. 77, e108–e110. doi:doi:10.1158/0008–5472.CAN-17–0307
Liang, W., Zhang, Y., Song, L., and Li, Z. (2019). 2,3′4,4′,5-Pentachlorobiphenyl Induces Hepatocellular Carcinoma Cell Proliferation through Pyruvate Kinase M2-dependent Glycolysis. Toxicol. Lett. 313, 108–119. doi:10.1016/j.toxlet.2019.06.006
Liberti, M. V., and Locasale, J. W. (2016). The Warburg Effect: How Does it Benefit Cancer Cells? Trends Biochem. Sci. 41, 211–218. doi:10.1016/j.tibs.2015.12.001
Ng, H. H. M., Lee, R. Y., Goh, S., Tay, I. S. Y., Lim, X., Lee, B., et al. (2020). Immunohistochemical Scoring of CD38 in the Tumor Microenvironment Predicts Responsiveness to Anti-PD-1/pd-L1 Immunotherapy in Hepatocellular Carcinoma. J. Immunother. Cancer 8 (2). doi:10.1136/jitc-2020–000987
Pascale, R. M., Calvisi, D. F., Simile, M. M., Feo, C. F., and Feo, F. (2020). The Warburg Effect 97 Years after its Discovery. Cancers 12 (10), 2819. doi:10.3390/cancers12102819
Peng, L., Chen, Y., Ou, Q., Wang, X., and Tang, N. (2020). LncRNA MIAT Correlates with Immune Infiltrates and Drug Reactions in Hepatocellular Carcinoma. Int. Immunopharmacology 89, 107071. doi:10.1016/j.intimp.2020.107071
Riera Leal, A., Ortiz-Lazareno, P. C., Jave-Suárez, L. F., Ramírez De Arellano, A., Aguilar-Lemarroy, A., Ortiz-García, Y. M., et al. (2020). 17β-estradiol-induced M-itochondrial D-ysfunction and Warburg E-ffect in C-ervical C-ancer C-ells A-llow C-ell S-urvival under M-etabolic S-tress. Int. J. Oncol 56 (1), 33–46. doi:doi:10.3892/ijo.2019.4912
Salik, B., Smyth, M. J., and Nakamura, K. (2020). Targeting Immune Checkpoints in Hematological Malignancies. J. Hematol. Oncol. 13 (1), 111. doi:10.1186/s13045–020–00947–6
Siegel, R. L., Miller, K. D., and Jemal, A. (2020). Cancer Statistics, 2020. CA A. Cancer J. Clin. 70, 7–30. doi:10.3322/caac.21590
Sun, L., Suo, C., Li, S.-t., Zhang, H., and Gao, P. (2018). Metabolic Reprogramming for Cancer Cells and Their Microenvironment: beyond the Warburg Effect. Biochim. Biophys. Acta (Bba) - Rev. Cancer 1870 (1), 51–66. doi:doi:10.1016/j.bbcan.2018.06.005
Terabe, K., Ohashi, Y., Tsuchiya, S., Ishizuka, S., Knudson, C. B., and Knudson, W. (2019). Chondroprotective Effects of 4-methylumbelliferone and Hyaluronan Synthase-2 Overexpression Involve Changes in Chondrocyte Energy Metabolism. J. Biol. Chem. 294 (47), 17799–17817. doi:10.1074/jbc.ra119.009556
Ždralević, M., Brand, A., Di Ianni, L., Dettmer, K., Reinders, J., Singer, K., et al. (2018). Double Genetic Disruption of Lactate Dehydrogenases A and B Is Required to Ablate the "Warburg Effect" Restricting Tumor Growth to Oxidative Metabolism. J. Biol. Chem. 293 (41), 15947–15961. doi:10.1074/jbc.RA118.004180
Zhang, X.-P., Chen, Q., Liu, Q., Wang, Y., Wang, F., Zhao, Z. M., et al. (2018). Development and Validation of GlycolysisRelated Prognostic Score for Prediction of Prognosis and Chemosensitivity. Biochem. Biophys. Rescommun 25 (4), 5615–5627. doi:10.1111/jcmm.16573
Zhang, Q., He, Y., Luo, N., Patel, S. J., Han, Y., Gao, R., et al. (2019). Landscape and Dynamics of Single Immune Cells in Hepatocellular Carcinoma. Cell 179 (4), 829–845. doi:10.1016/j.cell.2019.10.003
Keywords: immune-related glycolysis gene, hepatocellular carcinoma, prognosis, immune microenvironment, PD-L1
Citation: Yang J, Zhang Y, Duan J, Huang X, Yu H and Hu Z (2022) A Glycolysis-Related Gene Signature Correlates With the Characteristics of the Tumor Immune Microenvironment and Predicts Prognosis in Patients With Hepatocellular Carcinoma. Front. Mol. Biosci. 9:834976. doi: 10.3389/fmolb.2022.834976
Received: 14 December 2021; Accepted: 16 March 2022;
Published: 25 April 2022.
Edited by:
Amaresh Kumar Sahoo, Indian Institute of Information Technology, Allahabad, IndiaReviewed by:
Yuan Wei, Sun Yat-sen University, ChinaTara Chand Yadav, Indian Institute of Technology Roorkee, India
Copyright © 2022 Yang, Zhang, Duan, Huang, Yu and Hu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zhongjie Hu, huzhongjie@ccmu.edu.cn
†These authors have contributed equally to this work