Linhua Zhu1†
Zhuoqun Lin2†Kai Wang2,3Jiaxin Gu2Xiaojing Chen2
Ruizhe Chen2Lingfang Wang2*
Xiaodong Cheng2*- 1Department of Obstetrics, Women’s Hospital, School of Medicine, Zhejiang University, Hangzhou, China
- 2Department of Gynecologic Oncology, Women’s Hospital, School of Medicine, Zhejiang University, Hangzhou, China
- 3Department of Obstetrics and Gynecology, Taizhou Hospital of Zhejiang Province Affiliated to Wenzhou Medical University, Linhai, China
Introduction: Ovarian cancer (OV) is a highly lethal gynecological malignancy with a poor prognosis. Lactate metabolism is crucial for tumor cell survival, proliferation, and immune evasion. Our study aims to investigate the role of lactate metabolism-related genes (LMRGs) in OV and their potential as biomarkers for prognosis, immune microenvironment, and immunotherapy response.
Methods: Ovarian samples were collected from the TCGA cohort. And 12 lactate-related pathways were identified from the MsigDB database. Differentially expressed genes within these pathways were designated as LMRGs, which undergo unsupervised clustering to identify distinct clusters based on LMRGs. Subsequently, we assessed survival outcomes, immune cell infiltration levels, Hallmaker pathway activation patterns, and chemotaxis among different subtypes. After conducting additional unsupervised clustering based on differentially expressed genes (DEGs), significant differences in the expression of LMRGs between the two clusters were observed. The differentially expressed genes were subjected to subsequent functional enrichment analysis. Furthermore, we construct a model incorporating LMRGs. Subsequently, the lactate score for each tumor sample was calculated based on this model, facilitating the classification of samples into high and low groups according to their respective lactate scores. Distinct groups examined disparities in survival prognosis, copy number variation (CNV), single nucleotide variation (SNV), and immune infiltration. The lactate score served as a quantitative measure of OV's lactate metabolism pattern and an independent prognostic factor.
Results: This study investigated the potential role of LMRGs in tumor microenvironment diversity and prognosis in OV, suggesting that LMRGs play a crucial role in OV progression and the tumor microenvironment, thus serving as novel indicators for prognosis, immune microenvironment status, and response to immunotherapy.
Introduction
Ovarian cancer (OV), the most fatal gynecological malignancy and the fifth leading cause of cancer-related mortality among women globally (1), predominantly presents at advanced stages with metastasis beyond the ovary, resulting in a dismal prognosis (2). Despite advancements in surgical techniques and chemotherapy options (3), the 5-year survival rate for ovarian cancer patients remains disappointingly low, hovering around 50% (4). Consequently, there is an urgent imperative to identify novel biomarkers and therapeutic targets specific to OV.
The occurrence and progression of tumors are intricately intertwined with the metabolic reprogramming exhibited by tumor cells, which autonomously regulate the fluxes of their diverse metabolic pathways to meet augmented bioenergetic and biosynthetic demands while mitigating oxidative stress crucial for cancer cell proliferation and survival (5, 6). Lactate is a metabolic by-product generated during glycolysis, the anaerobic conversion of glucose into energy (7). Tumor cells frequently exhibit enhanced glycolytic activity and lactate production, even in the presence of oxygen, which is commonly referred to as the Warburg effect (8). Additionally, The role of lactate as a metabolic feedback regulator and unique signaling molecule has been extensively investigated, shedding light on its involvement in various physiological and pathological processes, including the regulation of energy metabolism, immune response, memory formation, wound healing, and tumor development (7, 9). Moreover, lactate can modulate the tumor microenvironment (TME) by inducing angiogenesis, remodeling the extracellular matrix, and suppressing immune responses (10). Moreover, lactate can hinder the activity and functionality of immune cells such as natural killer cells, cytotoxic T cells, and dendritic cells through pH reduction and interference with signaling pathways (11). Furthermore, lactate can impact the expression and function of immune checkpoints like programmed cell death protein 1 (PD-1), which plays a crucial role in regulating immune response efficacy and immunotherapy effectiveness (12).
The TME encompasses not only tumor cells, signaling molecules, and extracellular matrix but also immune cells, fibroblasts, and bone marrow-derived inflammatory cells closely associated with tumor cells (13). Lactate is a crucial metabolite within the TME (9). As a result of the Warburg effect, tumor cells secrete a substantial amount of lactate into the extracellular microenvironment. Lactate and H+ released by tumor cells in the TME contribute to an acidic TME formation, modulate the metabolism of innate and adaptive immune cells, and hinder the activation and proliferation of CD8+ T cells, natural killer (NK) cells, as well as dendritic cells. Lactate can also affect the expression and function of immune checkpoints, such as PD-1, which are involved in the regulation of immune response and the efficacy of immunotherapy (12).
Given the pivotal role of lactate metabolism in ovarian cancer and its potential as a therapeutic target, our study aimed to investigate the expression and prognostic significance of lactate metabolism-related genes (LMRGs) in ovarian cancer. Additionally, we analyzed the correlation between LMRGs and immune microenvironment characteristics, as well as their impact on the effectiveness of immunotherapy. Our study provides novel insights into molecular and clinical characteristics of ovarian cancer subtypes based on lactate metabolism while suggesting potential biomarkers and therapeutic strategies.
Materials and methods
Acquisition of data and preprocessing
The TCGA OV RNAseq dataset, along with the phenotype and survival data of patients, was utilized in this study (14) (https://portal.gdc.cancer.gov/). Additionally, the Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE140082) GSE140082 queue provided OV RNAseq gene expression and patient survival information. Furthermore, the gencode annotation was obtained from the European Bioinformatics Institute (EBI) website (https://www.ebi.ac.uk/). 63 immune checkpoint genes were sourced from the literature (15). 58 chemokines and 12 lactate metabolism-related pathways were acquired from the Molecular Signatures Database (MsigDB) database (http://www.gsea-msigdb.org/gsea/msigdb). Table 1 summarizes detailed clinical data of OV patients in the TCGA database, while Table 2 presents comprehensive clinical information on OV patients in the GSE140082 database.
Table 1 Clinical characteristics of OV patients in TCGA database.
Table 2 Clinical characteristics of OV patients in GSE140082 database.
Identification of the LMRGs in OV
The MsigDB database (16) provided 12 lactate metabolism-related pathways, and the GSVA (17) and GSEABase packages(c5.all.v7.4.symbols.gmt) were used to calculate the ontology enrichment scores of each sample in the TCGA-OV database on these 12 pathways. T-test was employed to determine the enrichment differences between tumor and normal samples, resulting in the identification of 325 genes with significant differences. Among them, a total of 276 protein-coding genes were defined as LMRGs.
Unsupervised clustering of LMRGs
Based on the LMRGs, an unsupervised clustering method was employed to categorize OV patients in TCGA into geneCluster A and geneCluster B. The “ConsensusClusterPlus” package in R was utilized for performing a consensus clustering algorithm to determine the number of clusters and their stability (18). The optimal number of clusters was selected based on the consensus matrix, cumulative distribution function (CDF), and relative change in the area under the CDF curve.
Consensus clustering of DEGs
To gain a deeper understanding of the distinct pathways of the clusters, we performed an analysis using the “limma” package in R to compare the expression levels of differentially expressed genes (DEGs) within 2 lactate clusters (19). The significance threshold was set at | logFC | > 0.8 and p < 0.05. Subsequently, a Consensus Clustering method (clusterAlg= “hc”, distance= “spearman”) based on DEGs was employed to classify TCGA OV patients into various lactate gene clusters.
Functional enrichment analysis
The WebGestalt database (http://www.webgestalt.org/) was utilized for Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) (20)pathway enrichment analysis (21). The clusterProfiler (22) and org.Hs.eg.db packages were employed to conduct GO and KEGG pathway enrichment analysis on these differential genes.
Construction of lactatescore and evaluation of effectiveness
Univariate Cox regression analysis was conducted to identify differentially expressed genes between the two lactate clusters. This analysis revealed 9 genes significantly associated with prognosis (p<0.01). Subsequently, principal component analysis (PCA) was performed using the expression levels of these nine genes. The sum of the first two principal components for each sample (lactate_score=∑(PC1i+PC2i)) was utilized as the lactate score. The survminer package was utilized for the classification of patients into high and low-risk groups based on lactate score. Both survminer and survival packages were employed to generate survival curves depicting High lactate score and Low lactate score in TCGA-OV and GSE140082 datasets, respectively. Additionally, the timeROC package was utilized to plot the ROC curve of the TCGA-OV dataset over 1-3 years.
Comprehensive analysis
The RCircos package was utilized to plot the chromosomal locations of LMRGs. PCA analysis on normal and tumor samples in the TCGA-OV dataset, based on LMRGs, was performed using the pca3d package. The Cytoscape package was employed to visualize the association relationships among LMRGs, while autocorrelation assessment of LMRGs was conducted using the Hmisc package. Hallmark pathway enrichment scores for each sample in the TCGA-OV dataset were calculated using the GSVA and GSEABase packages (h.all.v7.4.symbols.gmt), along with gene signature-based (23) enrichment scores for each sample. Enrichment differences between lactate clusters were computed and compared using the ggplot2 and ggpubr packages.
Gene mutation analysis
The mutations of genes regarding copy number variation (CNV) and single-nucleotide variation (SNV) in OV patients were explored by the online database Gene Set Cancer Analysis (GISTIC)(https://cloud.genepattern.org/gp/pages/index.jsf?lsid=urn:lsid:broad.mit.edu:cancer.software.genepattern.module.analysis:00125:6.15.28).
Measurement of immunocyte infiltration
The CIBERSORT algorithm (24) was used to measure the infiltration score of 22 common immunocytes. The relationships between risk score and immunocyte infiltration were measured by the Spearman coefficients using the “Hmisc” R package. The R ggplot 和ggpubr package was used to measure the infiltration score of 21 common immunocytes.
Statical analysis
The R software (version 3.6.3) was utilized for data analysis. To assess the impact of identified LMRGs on patients’ prognosis, Kaplan-Meier plotter analysis was employed. Univariate and multivariate Cox regression analyses were conducted to evaluate the influence of various factors on patients’ prognosis. Statistical significance levels were defined as follows: ns denoting not significant; *P < 0.05; **P < 0.01; ***P < 0.001; ****P < 0.0001, with ns indicating no statistical significance.
Results
Identification of LMRGs
The MsigDB database was used to explore the lactate-related pathways, and a total of 12 lactate-related pathways were identified as follows: GOBP LACTATE METABOLIC PROCESS, GOBP LACTATE TRANSMEMBRANE TRANSPORT, GOMF LACTATE DEHYDROGENASE ACTIVITY, GOMF LACTATE TRANSMEMBRANE TRANSPORTER ACTIVITY HP ABSOLUTE BRAIN LACTATE LEVEL BY MRS, HP ABSOLUTE LACTATE DEHYDROGENASE LEVEL, HP CREATE CSF LACTATE, HP INCREASED LACTATE DEHYDROGENASE LEVEL, HP INCREASED SERUM LACTATE, HP LACTIC ACIDOSIS, HP LACTICACIDURIA, HP SEVERE LACTIC ACIDOSIS (Figure 1A). Their enrichment differences in tumor and normal samples from the OV-TCGA cohort were demonstrated in Figure 1B. Among them, the pathway of HP_SEVERE_LACTIC_ACDOSIS obtained the most significant differences between tumor and normal samples.
Figure 1 Identification of LMRGs. (A) Ontology enrichment score heat map of 12 lactate metabolism-related pathways in the TCGA-OV. (B) 10 lactate metabolism-related pathways exhibiting statistically significant differences (p<0.05). (C) Locations of the LMRGs on 23 chromosomes. (D) Boxplot of the expressions of the top 20 LMRGs in the TCGA-OV cohort. (E) Principal components analysis (PCA) of the normal samples and the tumor samples based on LMRGs. (F) Mutation frequency of the top 20 of 276 differentially expressed LMRGs in 432 patients with OV. (G) Copy number variations (CNVs) of the top 20 of 276 differentially expressed LMRGs. ****P < 0.0001.
A total of 325 genes involved in these lactate-related pathways were collected, out of which 276 were protein-coding genes. Their locations on the 24 human chromosomes are displayed in Figure 1C. The top 20 genes with significant differences are displayed in Figure 1D. Principal components analysis (PCA) revealed that the normal samples could be separated from the tumor samples based on these genes (Figure 1E). We then further measured their gene mutation characteristics in OV-TCGA samples (Figure 1F). Among them, TP53 was the most frequently mutated gene, with a mutation frequency of approximately 88%. Their CNV mutation characteristics were also measured (Figure 1G). The relationships between the regulation of genes related to lactate metabolism and their effects on patient prognosis are displayed in Figure 2A. The correlations between the top 20 genes related to lactate metabolism were displayed in Figure 2B, and their effects on immune cell infiltration were also measured (Figure 2C).
Figure 2 Biological characterizations of LMRGs. (A) Correlation and prognostic value of LMRGs in OV. The line connecting the LMRGs represents their correlation, with the line thickness indicating the strength of the correlation between LMRGs. Blue and red represent negative and positive correlations, respectively. (B) The autocorrelation of LMRGs is shown in a figure with the top 20 genes represented in both rows and columns. Red indicates a positive correlation, blue indicates a negative correlation. (C) The correlation between LMRGs and immune infiltrating cells, with red indicating positive correlation and blue indicating negative correlation. *P < 0.05; **P < 0.01; ***P < 0.001; ****P < 0.0001. ns, no significance.
Unsupervised clustering and subtype analysis based on LMRGs
Based on the identified LMRGs, we divided all OV-TCGA patients into two subtypes using unsupervised clustering (Figures 3A–C). The heatmap showed the expression of clusters A and B (Figure 3D), which indicated their vast differences in gene expression profile. Furthermore, Kaplain-Meier survival curve analysis revealed that patients of cluster A significantly obtained unfavorable prognosis compared to those of cluster B (P = 0.011, Figure 3E).
Figure 3 Unsupervised clustering based on LMRGs. (A) Cumulative distribution function (CDF) of consensus clustering at k = 2–4. (B) relative changes in area under the CDF curve. (C) Unsupervised clustering of the 204 differentially expressed LMRGs and optimal consensus matrices for k =2. (D) Differences in expression levels of LMRGs between lactate clusters A and B. (E) Kaplan-Meier analysis estimating the overall survival between lactate clusters A and B.
We then performed a GSVA score of 50 cancer hallmarks to reveal the difference between the two clusters. As shown in Figure 4A, there were obvious differences in certain cancer hallmarks between these two groups. For example, the immune-related pathways (e.g., interferon-alpha response, IL2 stat5 signaling, complement, inflammatory response) were active in cluster A, while the metabolism (e.g., glycolysis, bile acid metabolism, fatty acid metabolism) and cell-cycle related (e.g., G2M checkpoint, E2F targets) pathways were active in cluster B. Comparison of these two subtypes in 13 classic signaling pathways also revealed their vast differences in function (Figure 4B). A comparison of these two subtypes in immune infiltration revealed that the infiltering abundance of plasma cells, CD8+ T cells, Tregs, activated NK cells, and DC cells in cluster B were significantly higher than that in cluster A (Figure 4C). Correspondingly, the expression of certain chemokines was also differentially expressed in these two subtypes (Figure 4D).
Figure 4 Subtype analysis based on LMRGs. (A) GSVA and GSEA analysis focused on the differential enrichment of Hallmaker pathways between lactate clusters A and B. (B) Boxplot of the differential enrichment of 13 classic gene signatures between lactate clusters A and B. (C) Boxplot of the relative expression of 21 infiltrating immune cells between lactate clusters A and B. (D) Boxplot of the relative expression of 57 chemokines between lactate clusters A and B. *P < 0.05; **P < 0.01; ***P < 0.001; ****P < 0.0001. ns, no significance.
Identification of geneClusters in OV
We then identified 94 DEGs (Differentially expressed genes) between clusters A and B and constructed two geneClusters (geneCluster A and geneCluster B) for OV based on unsupervised clustering (Figure 5A). Further analysis identified 80 DEGs between geneCluster A and geneCluster B. Functional enrichment these 80 DEGs were mainly involved in the pathways of metabolism (oxidative phosphorylation, ATP synthesis coupled electron transport, ATP metabolic process, NADH dehydrogenase complex) regarding (Figures 5B–E). Kaplain-Meier survival curve analysis revealed that patients of geneCluster A significantly obtained better prognosis compared to that of geneCluster B (P = 0.0011, Figure 5F). Infiltrating abundance analysis of immune cells revealed that there were vast differences between geneCluster A and geneCluster B (Figure 5G). For example, CD8+T cells, and activated NK cells were more abundant in geneCluster A. While M2 macrophages were more abundant in geneCluste B. Based on the aforementioned findings, we propose that lactate plays a pivotal role in the characterization of immune infiltration.
Figure 5 Identification of geneClusters in OV. (A) Heatmap illustrating the differences in clinicopathologic features and expression levels of lactate-Phenotype genes between the lactate.Gene.cluster A and B. (B) The differential enrichment of GO BP between the lactate.Gene.cluster A and B. (C) The differential enrichment of GO CC between the lactate.Gene.cluster A and B. (D) The differential enrichment of MF between the lactate.Gene.cluster A and B. (E) The differential enrichment of KEGG pathways between the lactate.Gene.cluster A and B. (F) Kaplan-Meier analysis estimating the overall survival between lactate.Gene.cluster A and B. (G) Boxplot of the relative expression of 21 infiltrating immune cells between the lactate.Gene.cluster A and B. *P < 0.05; **P < 0.01; ***P < 0.001. ns, no significance.
Lactate score was an independent prognostic risk factor in OV
Based on the identified 94 DEGs between clusters A and B, nine pivotal genes: PSMA2, LRP1, TXNDC17, TMEM258, BLOC1S1, HINT2, SOGA1, COA6, and DPM3, which have a significant impact on patient prognosis, were carefully selected to calculate the patient lactate score. As shown in Figure 6A, patients with high lactate scores significantly obtained unfavorable prognosis in the OV-TCGA cohort (P = 5e-04) compared to those with low lactate score, which were further validated in the GSE140082 cohort (P = 0.019, Figure 6C). ROC curve displayed the lactate score had good diagnostic efficiency in predicant patients with 1/2/3-year survivals (Figure 6B). Multivariate COX regression analysis revealed that low lactate score was an independent prognostic factor for OV patients (HR = 0.67, 95%CI = 0.51-0.88, P = 0.004, Figure 6D). Further analysis of gene mutations revealed that patients in the high and low lactate score groups exhibited similar patterns of mutations. The TP53, TTN, and MUC16 genes were found to be the most frequently mutated genes in these groups (Figures 7A, B). Their CNV patterns were also measured (Figures 7C, D).
Figure 6 Construction of lactatescore. (A) Kaplan-Meier analysis in high versus low lactatescore patients. (B) ROC curves to predict sensitivity and specificity of the lactate score in predicting 1-, 2-, and 3-month survival. (C) Kaplan-Meier analysis of both groups in the GSE140082. (D) lactatescore, TCGA molecular classification, age, stage, and grade were analyzed by multivariate Cox regression analysis.
Figure 7 Genomic states of difference between high- and low-lactatescore. (A, B) Mutated genes in the TCGA-OV samples of low and high lactatescore subgroups, respectively (C, D) Amplify and delete in copy numbers in groups with low and high lactatescore subgroups.
Discussion
OV is a prevalent gynecological malignancy worldwide. From a cellular perspective, OV can be classified into epithelial and non-epithelial tumors, with epithelial ovarian cancer accounting for approximately 90% of cases (25, 26). The epithelial subtype further encompasses four distinct histological variants: serous, endometrioid, mucinous, and clear cell types (27, 28). Currently, available therapeutic approaches for OV are limited to radical surgery and chemotherapy regimens that primarily prolong disease-free intervals without significantly impacting overall patient survival rates (29, 30). Therefore, there is an urgent need to establish effective stratification mechanisms along with reliable prognostic indicators to facilitate informed selection among various treatment modalities.
In the tumor microenvironment, cancer cells undergo metabolic adaptations to sustain tumorigenicity and withstand the cytotoxic effects of chemotherapy and immune cell-mediated attacks (31, 32). Alterations in metabolic pathways not only facilitate the survival and proliferation of ovarian cancer cells but also confer abilities for metastasis, chemoresistance acquisition, maintenance of cancer stemness traits, and evasion from immunological surveillance against tumors (33, 34). Since the discovery of lactate in 1780, it has often been mistakenly perceived as a metabolic waste product under hypoxic conditions, leading to numerous detrimental effects and its association with hypoxia (35). Currently, our understanding of lactate has expanded beyond its role as a byproduct of glycolysis to encompass that Lactate is now acknowledged as a pivotal carbon source for cellular metabolism and a crucial signaling molecule in physiological, chronically inflamed, and neoplastic tissues (36). Studies have substantiated that lactate functions as a signaling molecule, playing an indispensable role in coordinating signal transmission among diverse cells, organs, and tissues (37). Lactate is an essential substrate for various physiological cellular functions and exerts regulatory influence on different aspects of energy metabolism and signal transduction (7, 38). Therefore, LMRGs are potential biomarkers and therapeutic targets for various types of cancer (39–44). In hepatocellular carcinoma (HCC), a novel prognostic signature composed of 6 key lactate metabolism-related genes (FKTN, PDSS1, PET117, PUS1, RARS1, and RNASEH1) was developed to predict the survival and tumor microenvironment of HCC patients (44). In endometrial cancer (EC), a risk signature based on 18 lactate metabolism-related genes was constructed to predict the clinical outcome and molecular characteristics of EC patients (40). The role of lactate metabolism-related genes in kidney renal clear cell carcinoma (KIRC) is explored by analyzing the expression patterns of 17 lactate metabolism-related genes in KIRC patient datasets (43). Three lactate metabolism-related genes (FBP1, HADH, and TYMP) were selected from the panel of 17 genes as they showed significant associations with prognosis in KIRC patients (43). These three genes were used to construct a lactate-related prognostic signature that could predict overall survival, tumor microenvironment status, and immune response in KIRC patients (43). In this study, we observed significant activation of metabolism-related pathways (e.g., glycolysis, bile acid metabolism, fatty acid metabolism) in cluster B, which exhibited a more favorable prognosis. Functional enrichment analysis of differentially expressed genes between gene clusters A and B revealed their predominant involvement in metabolic pathways (oxidative phosphorylation, ATP synthesis coupled electron transport, ATP metabolic process, NADH dehydrogenase complex). Kaplan-Meier survival curve analysis demonstrated a significant difference in the prognoses of gene clusters A and B. These findings suggest that metabolism-related pathways play a crucial role in tumor development.
The molecular subgroups of high-grade serous ovarian carcinoma (HGSOC) have been distinguished by research, differentiating between low- and high-OXPHOS phenotypes. High-OXPHOS HGSOCs rely on oxidative phosphorylation as their primary metabolic pathway, while low-OXPHOS HGSOCs exhibit a preference for glycolysis. Additionally, the high-OXPHOS cells demonstrate increased chemosensitivity and a more favorable prognosis (45). In fast-growing cancer cells, there is often a shift from less efficient ATP production through aerobic glycolysis to OXPHOS (46). However, the single-cell multi-omics sequencing analysis of HGSOC revealed a simultaneous enhancement in the activities of both the oxidative phosphorylation pathway and glycolysis pathway during tumor development. Interestingly, mitochondrial oxygen consumption was not completely shut down in tumor cells but rather maintained as a backup system for OXPHOS activity (5). Further analysis indicated that an increase in chromosome 8 copy number during HGSOC development resulted in elevated MYC gene dosage, thereby enhancing the activities of both oxidative phosphorylation and glycolysis pathways (47).
Tumor cells modulate lipid absorption and stimulate lipogenesis to augment intracellular lipid abundance (48), thereby meeting the heightened energy demands for cell proliferation (49). The high prevalence of omental metastasis in ovarian cancer suggests that a lipid-rich environment may facilitate ovarian growth. Further investigations have revealed the upregulation of fatty acid-binding protein 4 (FABP4) in ovarian cancer with momentum metastasis, and its depletion significantly impedes the dissemination of ovarian cancer (49) thus confirming the crucial role of lipid uptake in the progression and advancement of this disease (50).
The intricate TME exerts a profound impact on both the progression of tumors and their response to therapeutic interventions. Currently, numerous studies are dedicated to elucidating the regulatory effects of metabolites and metabolic processes on immune function (51). Among these, lactate emerges as the most significantly elevated metabolic substance in tumors (52). The generation of lactate is generally perceived as a mechanism employed by tumor cells to evade immune surveillance (53, 54). In this study, a comparison of immune infiltration between clusters A and B demonstrated that cluster B exhibited significantly higher levels of plasma cells, CD8+ T cells, regulatory T cells (Tregs), activated natural killer (NK) cells, and dendritic cells (DCs) compared to cluster A. For instance, tumor cells stimulate lactic acid production through LDHA, which disrupts the secretion of IFN-γ and other cytokines in tumor-infiltrating T cells and NK cells, thereby promoting tumor growth (55). Alternatively, it modulates RNA methyltransferase METTL3-mediated myeloid cell infiltration into tumors via protein lactation. N6-methyladenosine (m6A) modification of TIMs enhances their immunosuppressive function (54). However, within the TME, the impact of lactate on both tumors and immune cells can be confounded by acidic protons generated during glycolysis. To eliminate interference from protons dissociated from its acidic counterpart - lactic acid - sodium lactate was investigated for its effect on immunity. This study suggests that the influence of lactic acid on immunity is two-fold; lactate can enhance CD8+ T cell dryness (functionality), thus augmenting anti-tumor immune response (56).
Additionally, there were notable variations observed in the abundance of immune cell infiltration between geneCluster A and geneCluster B. It is worth mentioning that seven different cell types exhibited statistically significant variances (p<0.05). For instance, geneCluster A demonstrated a higher prevalence of CD8+T cells and activated NK cells, which were associated with a more favorable prognosis. Conversely, geneCluster B showed an increased presence of M2 macrophages. Additionally, a study elucidated the interplay between pituitary adenoma (PA) cells and the immune tumor microenvironment (TME), highlighting their role in promoting PA invasion through M2 polarization. Furthermore, lactate levels demonstrated a positive correlation with PA invasion (57).
In recent years, immunotherapy has effectively changed the treatment landscape of OV (58). Immune checkpoint inhibitors (ICIs) have shown clinical benefit for patients with OV. However, OV patient’s response to immunotherapy is limited (59). Thus, the evaluation of sensitive/resistant target treatment subpopulations based on stratification by tumor biomarkers may improve the predictiveness of response to immunotherapy. By analyzing the immunotherapy datasets, we proved that patients with high lactate score were resistant to ICI treatment in GSE126044 cohorts.
In summary, we initially defined genes that were differentially expressed in lactate-related pathways as LMRGs and subsequently stratified OV patients into two distinct clusters, exhibiting significant disparities in terms of survival prognosis, and immune infiltration, among other factors. Furthermore, we additionally developed the lactate score as an independent prognostic indicator for predicting the outcome of OV patients. These findings highlight the important clinical significance of LMRGs and provide a new idea for guiding the personalized treatment strategy of OV patients.
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 below: TCGA OV RNAseq dataset (https://portal.gdc.cancer.gov/), GEO (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE140082), EBI (https://www.ebi.ac.uk/), MsigDB (https://www.gsea-msigdb.org/gsea/msigdb).
Author contributions
XCheng: Conceptualization, Funding acquisition, Resources, Supervision, Writing – original draft, Writing – review & editing. LZ: Data curation, Methodology, Software, Writing – original draft. ZL: Formal analysis, Investigation, Project administration, Validation, Writing – original draft. KW: Formal analysis, Methodology, Software, Writing – original draft. JG: Data curation, Methodology, Writing – original draft. XChen: Methodology, Writing – original draft. RC: Formal analysis, Methodology, Writing – original draft. LW: Conceptualization, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the National Key R&D Program of China (2022YFC2704200,2022YFC2704203) and the National Natural Science Foundation of China, Grant No. 382273348.
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.
The reviewer SL declared a past co-authorship with the author XC to the handling editor.
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.
References
1. Siegel RL, Miller KD, Wagle NS, Jemal A. Cancer statistics, 2023. CA. Cancer J Clin. (2023) 73:17–48. doi: 10.3322/caac.21763
2. Narod S. Can advanced-stage ovarian cancer be cured? Nat Rev Clin Oncol. (2016) 13:255–61. doi: 10.1038/nrclinonc.2015.224
3. Jelovac D, Armstrong DK. Recent progress in the diagnosis and treatment of ovarian cancer. CA. Cancer J Clin. (2011) 61:183–203. doi: 10.3322/caac.20113
4. National Cancer Institute. Cancer stat facts: ovarian cancer. (2023). Available at: https://seer.cancer.gov/statfacts/html/ovary.html.
5. Martínez-Reyes I, Chandel NS. Cancer metabolism: looking forward. Nat Rev Cancer. (2021) 21:669–80. doi: 10.1038/s41568-021-00378-6
6. Kroemer G, Pouyssegur J. Tumor cell metabolism: cancer’s Achilles’ heel. Cancer Cell. (2008) 13:472–82. doi: 10.1016/j.ccr.2008.05.005
7. Li X, Yang Y, Zhang B, Lin X, Fu X, An Y, et al. Lactate metabolism in human health and disease. Signal Transduct Target Ther. (2022) 7(1):305. doi: 10.1038/s41392-022-01151-3
8. Warburg O. On the origin of cancer cells. Science. (1956) 123:309–14. doi: 10.1126/science.123.3191.309
9. Hirschhaeuser F, Sattler UGA, Mueller-Klieser W. Lactate: a metabolic key player in cancer. Cancer Res. (2011) 71:6921–5. doi: 10.1158/0008-5472.CAN-11-1457
10. Gottfried E, Kreutz M, Mackensen A. Tumor metabolism as modulator of immune response and tumor progression. Semin Cancer Biol. (2012) 22:335–41. doi: 10.1016/j.semcancer.2012.02.009
11. Nakagawa Y, Negishi Y, Shimizu M, Takahashi M, Ichikawa M, Takahashi H. Effects of extracellular pH and hypoxia on the function and development of antigen-specific cytotoxic T lymphocytes. Immunol Lett. (2015) 167:72–86. doi: 10.1016/j.imlet.2015.07.003
12. Kumagai S, Koyama S, Itahashi K, Tanegashima T, Lin YT, Togashi Y, et al. Lactic acid promotes PD-1 expression in regulatory T cells in highly glycolytic tumor microenvironments. Cancer Cell. (2022) 40:201–218.e9. doi: 10.1016/j.ccell.2022.01.001
13. Quail DF, Joyce JA. Microenvironmental regulation of tumor progression and metastasis. Nat Med. (2013) 19:1423–37. doi: 10.1038/nm.3394
14. Liu J, Lichtenberg T, Hoadley KA, Poisson LM, Lazar AJ, Cherniack AD, et al. An integrated TCGA pan-cancer clinical data resource to drive high-quality survival outcome analytics. Cell. (2018) 173(2):400–416.e11. doi: 10.1016/j.cell.2018.02.052
15. Hu F-F, Liu C-J, Liu L-L, Zhang Q, Guo A-Y. Expression profile of immune checkpoint genes and their roles in predicting immunotherapy response. Brief. Bioinform. (2021) 22:bbaa176. doi: 10.1093/bib/bbaa176
16. Liberzon A, Birger C, Thorvaldsdóttir H, Ghandi M, Mesirov JP, Tamayo P. “The Molecular Signatures Database (MSigDB) hallmark gene set collection,”. Cell Syst. (2015) 1:417–25. doi: 10.1016/j.cels.2015.12.004
17. Hänzelmann S, Castelo R, Guinney J. “GSVA: gene set variation analysis for microarray and RNA-seq data,”. BMC Bioinf. (2013) 14:7. doi: 10.1186/1471-2105-14-7
18. Wilkerson MD, Hayes DN. “ConsensusClusterPlus: a class discovery tool with confidence assessments and item tracking,”. Bioinforma. Oxf. Engl. (2010) 26:1572–3. doi: 10.1093/bioinformatics/btq170
19. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. “limma powers differential expression analyses for RNA-sequencing and microarray studies,”. Nucleic Acids Res. (2015) 43:e47. doi: 10.1093/nar/gkv007
20. Ogata H, Goto S, Sato K, Fujibuchi W, Bono H, Kanehisa M. “KEGG: kyoto encyclopedia of genes and genomes,”. Nucleic Acids Res. (1999) 27:29–34. doi: 10.1093/nar/27.1.29
21. Liao Y, Wang J, Jaehnig EJ, Shi Z, Zhang B. “WebGestalt 2019: gene set analysis toolkit with revamped UIs and APIs,”. Nucleic Acids Res. (2019) 47:W199–205. doi: 10.1093/nar/gkz401
22. Yu G, Wang L-G, Han Y, He Q-Y. “clusterProfiler: an R package for comparing biological themes among gene clusters,”. Omics J Integr Biol. (2012) 16:284–7. doi: 10.1089/omi.2011.0118
23. Mariathasan S, Turley SJ, Nickles D, Castiglioni A, Yuen K, Wang Y, et al. “TGFβ attenuates tumour response to PD-L1 blockade by contributing to exclusion of T cells,”. Nature. (2018) 554:544–8. doi: 10.1038/nature25501
24. Newman AM, Steen CB, Liu CL, Gentles AJ, Chaudhuri AA, Scherer F, et al. “Determining cell type abundance and expression from bulk tissues with digital cytometry,”. Nat Biotechnol. (2019) 37:773–82. doi: 10.1038/s41587-019-0114-2
25. Lu Z, Gao Y. “Screening differentially expressed genes between endometriosis and ovarian cancer to find new biomarkers for endometriosis,”. Ann Med. (2021) 53:1377–89. doi: 10.1080/07853890.2021.1966087
26. Xu H, Zou R, Li F, Liu J, Luan N, Wang S, et al. “MRPL15 is a novel prognostic biomarker and therapeutic target for epithelial ovarian cancer,”. Cancer Med. (2021) 10:3655–73. doi: 10.1002/cam4.3907
27. D’Alessandris N, Travaglino A, Santoro A, Arciuolo D, Scaglione G, Raffone A, et al. “TCGA molecular subgroups of endometrial carcinoma in ovarian endometrioid carcinoma: A quantitative systematic review,”. Gynecol. Oncol. (2021) 163:427–32. doi: 10.1016/j.ygyno.2021.08.011
28. Bao M, Zhang L, Hu Y. “Novel gene signatures for prognosis prediction in ovarian cancer,”. J Cell Mol Med. (2020) 24:9972–84. doi: 10.1111/jcmm.15601
30. Stewart C, Ralyea C, Lockwood S. “Ovarian cancer: an integrated review,”. Semin Oncol Nurs. (2019) 35:151–6. doi: 10.1016/j.soncn.2019.02.001
31. Zhang Y, Wang Y, Zhao G, Orsulic S, Matei D. “Metabolic dependencies and targets in ovarian cancer,”. Pharmacol Ther. (2023) 245:108413. doi: 10.1016/j.pharmthera.2023.108413
32. Finley LWS. “What is cancer metabolism?,”. Cell. (2023) 186:1670–88. doi: 10.1016/j.cell.2023.01.038
33. Vander Heiden MG, Cantley LC, Thompson CB. “Understanding the Warburg effect: the metabolic requirements of cell proliferation,”. Science. (2009) 324:1029–33. doi: 10.1126/science.1160809
34. Vander Heiden MG, DeBerardinis RJ. Understanding the intersections between metabolism and cancer biology. Cell. (2017) 168:657–69:2017. doi: 10.1016/j.cell.2016.12.039
35. Ferguson BS, Rogatzki MJ, Goodwin ML, Kane DA, Rightmire Z, Gladden LB. “Lactate metabolism: historical context, prior misinterpretations, and current understanding,”. Eur J Appl Physiol. (2018) 118:691–728. doi: 10.1007/s00421-017-3795-6
36. Zhang D, Li Y, Yang S, Wang M, Yao J, Zheng Y, et al. “Identification of a glycolysis-related gene signature for survival prediction of ovarian cancer patients,”. Cancer Med. (2021) 10:8222–37. doi: 10.1002/cam4.4317
37. Certo M, Llibre A, Lee W, Mauro C. “Understanding lactate sensing and signalling,”. Trends Endocrinol Metab TEM. (2022) 33:722–35. doi: 10.1016/j.tem.2022.07.004
38. Zheng M, Mullikin H, Hester A, Czogalla B, Heidegger H, Vilsmaier T, et al. “Development and validation of a novel 11-gene prognostic model for serous ovarian carcinomas based on lipid metabolism expression profile,”. Int J Mol Sci. (2020) 21:9169. doi: 10.3390/ijms21239169
39. Wang C, Lu T, Xu R, Luo S, Zhao J, Zhang L. “Multi-omics analysis to identify lung squamous carcinoma lactate metabolism-related subtypes and establish related index to predict prognosis and guide immunotherapy,”. Comput Struct Biotechnol J. (2022) 20:4756–70. doi: 10.1016/j.csbj.2022.08.067
40. Shi R, Li H, Wei S, Yu Z, Zhang J, Zhang Q, et al. “Lactate metabolism-related genes to predict the clinical outcome and molecular characteristics of endometrial cancer,”. BMC Cancer. (2023) 23:491. doi: 10.1186/s12885-023-10934-y
41. Zhao F, Wang Z, Li Z, Liu S, Li S. “Identifying a lactic acid metabolism-related gene signature contributes to predicting prognosis, immunotherapy efficacy, and tumor microenvironment of lung adenocarcinoma,”. Front Immunol. (2022) 13:980508. doi: 10.3389/fimmu.2022.980508
42. Yang L, Tan P, Sun H, Zeng Z, Pan Y. “Integrative dissection of novel lactate metabolism-related signature in the tumor immune microenvironment and prognostic prediction in breast cancer,”. Front Oncol. (2022) 12:874731. doi: 10.3389/fonc.2022.874731
43. Sun Z, Tao W, Guo X, Jing C, Zhang M, Wang Z, et al. “Construction of a lactate-related prognostic signature for predicting prognosis, tumor microenvironment, and immune response in kidney renal clear cell carcinoma,”. Front Immunol. (2022) 13:818984. doi: 10.3389/fimmu.2022.818984
44. Li Y, Mo H, Wu S, Liu X, Tu K. “A novel lactate metabolism-related gene signature for predicting clinical outcome and tumor microenvironment in hepatocellular carcinoma,”. Front Cell Dev Biol. (2021) 9:801959. doi: 10.3389/fcell.2021.801959
45. Gentric G, Kieffer Y, Mieulet V, Goundiam O, Bonneau C, Nemati F, et al. “PML-regulated mitochondrial metabolism enhances chemosensitivity in human ovarian cancers,”. Cell Metab. (2019) 29:156–173.e10. doi: 10.1016/j.cmet.2018.09.002
46. Cairns RA, Harris IS, Mak TW. “Regulation of cancer cell metabolism,”. Nat Rev Cancer. (2011) 11:85–95. doi: 10.1038/nrc2981
47. Wang Y, Xie H, Chang X, Hu W, Li M, Li Y, et al. “Single-cell dissection of the multiomic landscape of high-grade serous ovarian cancer,”. Cancer Res. (2022) 82:3903–16. doi: 10.1158/0008-5472.CAN-21-3819
48. Broadfield LA, Pane AA, Talebi A, Swinnen JV, Fendt S-M. “Lipid metabolism in cancer: New perspectives and emerging mechanisms,”. Dev Cell. (2021) 56:1363–93. doi: 10.1016/j.devcel.2021.04.013
49. Nieman KM, Kenny HA, Penicka CV, Ladanyi A, Buell-Gutbrod R, Zillhardt MR, et al. Adipocytes promote ovarian cancer metastasis and provide energy for rapid tumor growth. Nat Med. (2011) 17:1498–503. doi: 10.1038/nm.2492
50. Currie E, Schulze A, Zechner R, Walther TC, Farese RV. “Cellular fatty acid metabolism and cancer,”. Cell Metab. (2013) 18:153–61. doi: 10.1016/j.cmet.2013.05.017
51. Ye L, Jiang Y, Zhang M. “Crosstalk between glucose metabolism, lactate production and immune response modulation,”. Cytokine Growth Factor Rev. (2022) 68:81–92. doi: 10.1016/j.cytogfr.2022.11.001
52. Goveia J, Pircher A, Conradi LC, Kalucka J, Lagani V, Dewerchin M, et al. Meta-analysis of clinical metabolic profiling studies in cancer: challenges and opportunities. EMBO Mol Med. (2016) 8(10):1134–42. doi: 10.15252/emmm.201606798
53. Apostolova P, Pearce EL. “Lactic acid and lactate: revisiting the physiological roles in the tumor microenvironment,”. Trends Immunol. (2022) 43:969–77. doi: 10.1016/j.it.2022.10.005
54. Xiong J, He J, Zhu J, Pan J, Liao W, Ye H, et al. “Lactylation-driven METTL3-mediated RNA m6A modification promotes immunosuppression of tumor-infiltrating myeloid cells,”. Mol Cell. (2022) 82:1660–1677.e10. doi: 10.1016/j.molcel.2022.02.033
55. Brand A, Singer K, Koehl GE, Kolitzus M, Schoenhammer G, Thiel A, et al. “LDHA-associated lactic acid production blunts tumor immunosurveillance by T and NK cells,”. Cell Metab. (2016) 24:657–71. doi: 10.1016/j.cmet.2016.08.011
56. Feng Q, Liu Z, Yu X, Huang T, Chen J, Wang J, et al. “Lactate increases stemness of CD8 + T cells to augment anti-tumor immunity,”. Nat Commun. (2022) 13:4981. doi: 10.1038/s41467-022-32521-8
57. Zhang A, Xu Y, Xu H, Ren J, Meng T, Ni Y, et al. “Lactate-induced M2 polarization of tumor-associated macrophages promotes the invasion of pituitary adenoma by secreting CCL17,”. Theranostics. (2021) 11:3839–52. doi: 10.7150/thno.53749
58. Mei WJ, Mi M, Qian J, Xiao N, Yuan Y, Ding PR. Clinicopathological characteristics of high microsatellite instability/mismatch repair-deficient colorectal cancer: A narrative review. Front Immunol. (2022) 13:1019582. doi: 10.3389/fimmu.2022.1019582
Keywords: ovarian cancer, metabolism, lactate, immune microenvironment, prognostic signature
Citation: Zhu L, Lin Z, Wang K, Gu J, Chen X, Chen R, Wang L and Cheng X (2024) A lactate metabolism-related signature predicting patient prognosis and immune microenvironment in ovarian cancer. Front. Endocrinol. 15:1372413. doi: 10.3389/fendo.2024.1372413
Received: 18 January 2024; Accepted: 15 February 2024;
Published: 11 March 2024.
Edited by:
Fang Ren, First Affiliated Hospital of Zhengzhou University, ChinaReviewed by:
Yasen Maimaitiyiming, Xinjiang Medical University, ChinaShourong Wang, Shandong University, China
Shitong Lin, Huazhong University of Science and Technology, China
Copyright © 2024 Zhu, Lin, Wang, Gu, Chen, Chen, Wang and Cheng. 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: Xiaodong Cheng, chengxd@zju.edu.cn; Lingfang Wang, 11818387@zju.edu.cn
†These authors have contributed equally to this work