Comprehensive Analysis of the Relationship Between Metabolic Reprogramming and Immune Function in Prostate Cancer

Abstract

Purpose

Prostate cancer is the most common malignant urinary tumor among men. Treatments are currently unsatisfactory for advanced prostate cancer. Cancer biology remains the basis for developing new antitumor drugs. Therefore, it is crucial to study the metabolic reprogramming, immune microenvironment, and immune evasion of tumors. This study aimed to clarify the relationship between tumor glycolysis and immune function in prostate cancer.

Materials and Methods

We downloaded the gene expression matrix and clinical data of prostate cancer from The Cancer Genome Atlas. We studied the expression profiles and prognostic significance of glycolysis-related genes and used CIBERSORT to identify the proportion of tumor-infiltrating immune cells. Through differential gene expression analysis, gene ontology analysis, Kyoto Encyclopedia of Genes and Genomes analysis, gene set enrichment analysis, and correlation analysis, we further explored the relationship between glycolytic activity and immune function. We also performed immunohistochemistry, Western blot and RT-qPCR experiments using human prostate cancer tissue and cell lines to verify the expression of some glycolytic genes, macrophage infiltration and polarization.

Results

Among glycolysis-related genes, the expression of SLC16A3 in prostate cancer tissues was lower than that in normal tissues, but its high expression was associated with poor prognosis. In the high SLC16A3 expression group, several glycolysis-related genes also showed high expression, which was confirmed by immunohistochemistry experiments and Western blot. In high-glycolysis group, the expression of immune-related genes and the interleukin-17 (IL-17) signaling pathway were upregulated. CD8⁺ T cells, regulatory T cells, macrophages, and other immune cells were highly enriched. Among them, M2 macrophage infiltration was associated with poor prognosis.

Conclusion

The enhanced glycolytic activity of prostate cancer may contribute to the formation of a pro-tumor immune microenvironment. The IL-17 signaling pathway may play an important mediating role in the interaction between tumor glycolysis and immune function.

Keywords:

• glycolysis
• tumor-infiltrating immune cells
• IL-17 signaling pathway
• TCGA

Introduction

Prostate cancer is a malignant tumor of the urinary system in men with high incidence. In 2018, there were about 1.3 million new cases of prostate cancer worldwide, resulting in 359,000 deaths and making it the fifth leading cause of cancer death.¹ Prostate cancer is the most common malignant tumor and the second leading cause of death in the United States, accounting for 174,650 estimated new cases and 31,620 estimated deaths in 2019.² In China, with the increasing use of ultrasound-guided prostate biopsy, the incidence of prostate cancer has increased sharply,³ with an estimated 72,000 new cases and 31,000 deaths in 2015.⁴ Compared with other malignant tumors, the progression of prostate cancer is relatively slow, resulting in a 5-year survival rate of > 80%.^5,6 However, although hormone deprivation therapy can inhibit progression of advanced prostate cancer, this effect cannot be sustained for a long time. Since then, the serum level of prostate-specific antibody rises rapidly again, and the tumor progresses to the stage of castration-resistant prostate cancer (CRPC), which eventually leads to death.

The metabolism reprogramming is an important sign of cancer. Most cancer cells exhibit an increase in glucose consumption and glycolysis rate and the formation of the final product lactate, even in aerobic conditions. This phenomenon is called “Warburg effect”,^7,8 which provides cancer cells with multiple survival advantages and regulates the metabolism and function of neighboring cells in the tumor microenvironment (TME). However, the phenomenon called the “reverse Warburg effect” also exists. Some cancer cells may exhibit oxidative phenotype, while stromal cells can become glycolytic. The two kinds of cells work together to boost tumor progression.⁹ Regardless of the metabolism mode, lactate is an important mediator of metabolic interaction between cancer cells and TME.¹⁰ Early prostate cancer tends to enhance oxidative phosphorylation and fat production but limit glycolysis. However, with tumor development, prostate cancer cells begin to exhibit the Warburg effect. As tumor aerobic glycolysis promotes malignant transformation and tumor progression, it has a potential significance for prognosis judgement and cancer treatment.^11–14

As an important part of the TME, tumor-infiltrating immune cells (TIICs) can change the immune microenvironment of tumors. Studies have shown that TIICs can interact with tumor cells and influence carcinogenesis and tumor progression, which makes targeting TIICs a promising method of antitumor therapy.^15,16 The interaction between prostate cancer and TIIC has been proved by experiments. Nardone et al found that in prostate cancer patients who underwent salvage radiotherapy after biochemical recurrence, the expression of FoxP3, CCR7, or PD-1 by tumor-infiltrating T lymphocytes was associated with progression-free survival and overall survival.¹⁷ Erlandsson et al reported that a large number of M2 macrophages in the prostate cancer microenvironment increased the probability of prostate cancer mortality.¹⁸

The highly acidic microenvironment caused by tumor glycolysis may affect the infiltration of immune cells, ultimately severely damaging the function of the immune system itself and promoting the growth, dissemination, and immune escape of primary tumors.^19,20 In in vitro experiments, lactate could damage the maturation and differentiation of dendritic cells (DCs) and induce a suppressor phenotype.^21,22 In vitro experiments also proved that an acidic pH can inhibit T cell responses, including the elimination of IFN-γ and TNF-α secretion.²³ However, based on clinical data in prostate cancer research, there are few comprehensive analyses of the relationship between tumor glycolysis, TIICs, and clinical characteristics. This study, therefore, aimed to investigate the status of glycolysis and immunity in prostate cancer using transcriptome and clinical data from The Cancer Genome Atlas (TCGA) project, to identify the relationship between tumor glycolysis and TIICs and reveal the potential role of TIICs in the prognosis of prostate cancer.

Materials and Methods

Data Preparation for Preprocessing Transcriptome and Clinical Data

Overall gene expression profiles and clinical data of patients were downloaded from the TCGA database (https://portal.gdc.cancer.gov/). We implemented a multi-array average algorithm and preprocessed the raw data. We then used the voom method to convert mRNA sequencing data to a microarray result.²⁴ Furthermore, we used “limma” package of R to standardize the mRNA sequencing results. Subsequently, further analysis was performed with 15,514 gene expression datasets obtained from 52 normal tissue samples and 498 tumor tissues. Clinicopathological data, such as tumor stage, Gleason grade, survival status, and survival time, were obtained from the TCGA database.

Differentially Expressed Genes and Enriched Pathways in Prostate Cancer

The tumor samples were divided into low-glycolysis and high-glycolysis groups based on the median expression of SLC16A3. The differential gene expression analysis was obtained by using limma R package. Heatmap and volcano map were used to display differentially expressed genes (DEGs). The threshold of DEGs was |log2FC|≥1, with the adjusted P-value was <0.05. The exception was that when obtaining DEGs related to glycolysis between normal and tumor samples, the threshold of DEGs was |log2FC|≥0.6. DAVID (https://david.ncifcrf.gov/) was used for gene ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis.²⁵ Gene set enrichment analysis (GSEA) based on WebGestalt (http://www.webgestalt.org/),²⁶ with a false discovery rate (FDR) of <0.05 as the cutoff criterion, was performed to determine the specific signal pathway between the above two groups.

Enumeration of TIICs

CIBERSORT is an analytical tool that can use gene expression characteristics of 547 genes to estimate the abundance of cell types in mixed cell populations.²⁷ It describes the characteristics of each immune cell subtype and uses a deconvolution algorithm to accurately quantify the different immune cell components. In this study, the LM22 signature file including 22 immune cell subtypes was used to illustrate the immune infiltration of each sample. After analyzing the differential gene expression profile on the CIBERSORT web (http://cibersort.stanford.edu), we obtained the distribution of 22 TIIC subtypes and the correlation coefficients, P-values, and root mean square errors (RMSEs).

Patients and Specimens

Human prostate cancer tissues were obtained from 60 patients of the Shenzhen University General Hospital between January 2018 and December 2020. All research procedures were approved by the Medical Ethical Committee of Shenzhen University Medical School according to the Declaration of Helsinki. The informed consents were obtained from all patients before tissue specimen collection.

Cell Culture

Human prostate cancer cell line PC3 and RAW246.7 macrophages were purchased from BeNa Culture Collection (Beijing, China). PC3 cells and RAW246.7 macrophages were cultured in PRMI 1640 medium (Invitrogen, Gaithersburg, MD, USA) containing 10% fetal bovine serum (Invitrogen, USA) under the condition of 37°C and 5% CO₂. IL-4 (10 ng/mL) was used to stimulate RAW246.7 macrophages to obtain anti-inflammatory macrophages (M2),²⁸ while control macrophages (M0) were cultured in medium alone. Pro-inflammatory macrophages (M1) were induced by treating RAW246.7 macrophages with 50 ng/mL interferon (IFN)-γ and 10 ng/mL lipopolysaccharide (LPS).²⁹

PC3 Cell Transfection of siRNA Adenovirus

PC3 cell lines were transfected with adenovirus vector packed with siRNA-MCT4 (RiboBio Corporation, Guangzhou, China), and polycoagulant amine (RiboBio Corporation, Guangzhou, China) was used as virus transfection efficiency enhancer. The specific sequence of interfering RNA fragment in siRNA-MCT4: 5ʹ-CCUACUCCGUCUACCUCUUTT-3ʹ (sense) and 5ʹ-AAGAGGUAGACGGAGUAGGTT-3ʹ (antisense); negative control (siRNA-NC) specific sequence: 5ʹ-UUCUCCGAACGUGUCACGUTT-3ʹ (sense) and 5ʹ-ACGUGACACGUUCGGAGAATT-3ʹ (antisense). The PC3 cell lines were laid on a 12-well plate and incubated overnight with 1mL culture medium in each well. After 24 hours of incubation, the culture medium was removed and the mixture of 1mL polycoagulant with culture medium was added to each well, and the concentration of polycondensed amine was 5 ug/mL. Adenoviral particles were melted on the ice and the virus transfection began after complete melting. Virus particles with 1×10⁷ viruses in each well were added to the culture medium to transfect cells. Then the 12-well plate was centrifuged at a low speed for 1.5 hours and incubated at 37°C for 24 hours. The culture medium was then removed and 1mL culture medium without polyamine was added to each well, and the cells were incubated overnight. After 4 days of culture, the cells were collected and cryopreserved.

PC3 and RAW264.7 Cocultured

Tumor-associated macrophages (TAMs) were obtained by incubating RAW246.7 macrophages with PC3 cells or si-MCT4-PC3 cells (PC3-TAM or si-MCT4-PC3-TAM).³⁰ We developed a co-culture model by seeding RAW246.7 macrophages (4×10⁴/well) in the lower chambers of the transwell plate with PC3 cells (1×10⁴/well) in the upper chambers. After co-culturing for 84 hours, RAW246.7 macrophages were collected and mRNA analysis was performed by real-time quantitative PCR.

Immunohistochemistry

Human prostate cancer tissues were cut into 3-mm sections. The paraffin sections were baked at 56°C for 2 h for dewaxing, boiled in citrate buffer to recover the antigen, and blocked with 3% hydrogen peroxide. Antibodies against monocarboxylate transporter 4 (MCT4) (dilution 1:200; ProteinTech Group, Chicago, IL, USA), PFKFB3 (dilution 1:100; Abcam, UK), ENO2 (dilution 1:300; Servicebio, Wuhan, China), and M2 macrophage marker CD163 (dilution 1:100; Servicebio, Wuhan, China) were used to detect protein expression at 4°C overnight, which was followed by incubation with HRP-conjugated secondary antibodies for 1 h at 37°C. The 3ʹ-diaminobenzidine kit (DAB; Dako Corporation, Carpinteria, CA, USA) was used for immune complex visualization. ImageJ software was used to measure the immunohistochemical staining intensity and the proportion of positive cells. Immunohistochemistry scores were independently evaluated by two doctors based on the method described in the article in a blinded manner.³¹ Scores of 1+ and 2+ were considered as low expression, while scores of 3+, 4+, and 5+ were regarded as high expression.²⁵ The number of M2 macrophages was calculated for the entire tissue specimen under low magnification (×10).

Western Blotting

The cells were divided into siRNA interference group (si-MCT4 group), negative control group (si-NC group) and blank control group (NC group). The cells were harvested and lysed in RIPA Lysis and Extraction Buffer (Invitrogen, USA). The protein concentration of the lysate was determined by Pierce BCA protein assay kit (Thermo, Waltham, MA, USA). The lysate was electrophoretic on 8% SDS polyacrylamide gel and transferred to PVDF membrane (Millipore, Boston, MA, USA). The membrane was probed with the primary antibodies overnight at 4°C. The primary antibodies included MCT4 (dilution 1:1000; ProteinTech Group, Chicago, IL, USA), PFKFB3 (dilution 1:1000; Abcam, UK), ENO2 (dilution 1:1000; Servicebio, Wuhan, China) and β-actin (dilution 1:1000; CST, Boston, MA, USA). Then the membrane was incubated with the horseradish peroxidase-conjugated secondary antibodies (Zhongshan Jinqiao, Beijing, China) and scanned by Odyssey imaging system. The relative quantitative analysis was completed by lmageJ software, and the experiment was repeated three times.

Real-Time Quantitative PCR (RT-qPCR)

RNA was extracted from PC3 cells or macrophages by using Trizol reagent (Invitrogen, USA), and cDNA was generated by using RT-PCR kit (TaKaRa, Kyoto, Japan). Real-time quantitative PCR was performed by using the SYBR Green QuantiTect RT-PCR kit (Roche, South San Francisco, CA), and each of the samples was analyzed three times. The sequences of the primers used for RT-qPCR were shown as follows:

MCT4, 5ʹ-TTTTGCTGCTGGGCAACTTCTTCTG-3ʹ (sense) and 5ʹ-TCACGTTGTCTCGAAGCATGGGTTT-3ʹ (antisense);

CD206, 5ʹ-GGCGAGCATCAAGAGTAAAGA-3ʹ (sense) and 5ʹ-CATAGGTCAGTCCCAACCAAA-3ʹ (antisense);

Arg1, 5ʹ-TTTTTCCAGCAGACCAGCTT-3ʹ (sense) and 5ʹ-CATGAGCTCCAAGCCAAAGT-3ʹ (antisense);

iNOS, 5ʹ-GAAACGCTTCACTTCCAATG-3ʹ (sense) and 5ʹ-AATCCACAACTCGCTCCAA-3ʹ (antisense);

IL-6, 5ʹ-AACGATGATGCACTTGCAGA-3ʹ (sense) and 5ʹ-TGGTACTCCAGAAGACCAGAGG-3ʹ (antisense);

β-actin, 5ʹ-CTAAGGCCAACCGTGAAAAG-3ʹ (sense) and 5ʹ-GGTACGACCAGAGGCATACA-3ʹ (antisense).

Statistical Analysis

All analyses were performed using R version 3.6.3, and all statistical tests were two-sided. Differences in continuous data between the groups were compared using analysis of variance. Categorical data were analyzed using the chi-square test or Fisher’s exact test. The prognostic significance of genes was analyzed by using Cox proportional hazards regression model. Spearman correlation analysis was used to assess the potential correlations. The Kaplan-Meier method was used to estimate the disease-free survival rate, and the Log rank test to evaluate the differences between the groups. P<0.05 was considered statistically significant.

Results

Expression and Prognostic Significance of Glycolysis-Related Genes

We selected 45 glycolysis-related genes and compared the expression profiles of these genes between normal tissues and tumor tissues. Results were showed in the heatmap and volcano map (Figure 1A and B). Compared with expression in normal prostate tissue, the expression of SLC2A9, ENO2, SLC2A4, SLC16A7, ALDOC, PFKFB3, LDHB, and SLC16A3 in tumor tissue decreased, while that of HK2, SLC16A8, and FBP1 increased.

Figure 1 Expression and prognostic significance of glycolysis-related genes in prostate cancer. (A) Heatmap analysis of expression levels of 45 glycolysis-related genes. (B) Volcano plot for differentially expressed glycolysis-related genes. (C and D) Univariate and multivariate Cox regression of glycolysis-related genes. (E) Kaplan-Meier analysis of disease-free survival comparing low and high SLC2A4 expression. (F) Kaplan-Meier analysis of disease-free survival comparing low and high SLC16A3 expression.

To find out the candidate glycolysis-related enzymes that may play an important role in the progression of prostate cancer, we investigated the relationship between clinical manifestations of prostate cancer and the expression of glycolysis-related genes. Univariate and multivariate Cox proportional hazards regression analyses were used to evaluate the gene-associated hazard ratio of disease-free survival. The results showed that SLC2A4 and SLC16A3 were independent risk factors for prostate cancer progression (Figure 1C and D). According to the SLC2A4 or SLC16A3 expression, the 550 prostate cancer patients from TCGA were divided into two groups. Survival analysis showed that the disease-free survival time was significantly shortened in the group with high SLC16A3 expression, which clarified the prognostic value of SLC16A3 (Figure 1E). In contrast, the group with high SLC2A4 expression predicted better disease-free survival (Figure 1F). We also found that a high SLC16A3 expression was positively correlated with a higher Gleason grade and worse disease-free survival (Table 1).

Table 1 Clinical Characteristics of 489 Samples Base on SLC16A3 Expression

Characteristics	SLC16A3
	Low Expression	High Expression	P-value
Age	60.30±7.01	61.66±6.60	0.028
T stage			0.320
T2	102	83
T3	136	152
T4	5	5
Tx	2	4
N stage			0.109
N0	173	167
N1	31	46
Nx	41	31
M stage			0.343
M0	223	225
M1	0	2
Mx	22	17
Gleason grade group			0.005
1	28	17
2	84	60
3	48	52
4	33	30
5	52	85
Overall survival			1
Alive	243	242
Dead	2	2
Disease-free survival			0.003
Disease-free	212	186
Recurred/progressed	33	58

SLC16A3 Expression and Glycolytic Activity

SLC16A3 belongs to the solute carrier 16 family of genes and encodes MCT4. MCT4 shows high affinity for lactate, but low affinity for pyruvate, which leads to the conversion of pyruvate to lactate before export, promotes the outflow of lactate from glycolysis cells, and maintains high glycolysis flux in tumor cells.³² In view of the important role of SLC16A3 in glycolysis, we chose SLC16A3 for further studies. In order to further clarify whether the SLC16A3 expression can represent the level of glycolysis, we compared the expression of other glycolysis-related genes between the low SLC16A3 and high SLC16A3 expression group. The results demonstrated that PFKFB3, ALDOA, ALDOC, ENO2, ENO3, LDHA, LDHB, PCK2, PFKL, PFKP and PGM1 were significantly upregulated in the high SLC16A3 expression group (Figure 2). All these genes are key genes in the glycolysis pathway. All above results suggested that SLC16A3 expression can reflect glycolytic activity, so the prostate cancer patients from TCGA can be divided into the low-glycolysis and high-glycolysis group according to the median expression of SLC16A3, respectively.

Figure 2 The expression of SLC16A3 can represent the glycolytic activity of prostate cancer to some extent. Box plots showed the expression of glycolysis-related genes ALDOA, ALDOC, ENO2, ENO3, LDHA, LDHB, PCK2, PFKFB3, PFKL, PFKP, and PGM1 in the low and high SLC16A3 expression groups.

Correlation Between Glycolytic Activity and Immune Function

To study the biological differences between the low- and high-glycolysis groups, we analyzed the profiles of 260 DEGs between the two groups, including 239 upregulated genes and 21 downregulated genes (Figure 3A). Many cytokines and inflammatory mediators in the high-glycolysis group were dysregulated and participated in the regulation of immune cell function (Figure 3B).

Figure 3 Correlation between glycolytic activity and immune function in prostate cancer. (A) Volcano plot for DEGs between the low- and high-glycolysis groups. (B) Heatmap analysis of dysregulated inflammation-related genes between the low- and high-glycolysis groups. (C and D) The top 20 most significant enrich Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) terms with 260 DEGs analyses.

The enrichment analyses of GO and KEGG were carried out to reveal the gene functions and biological pathways. The functions related to immunity included GO:0006958 complement activation, classical pathway, GO:0006955 immune response, GO:0006956 complement activation, GO:0050871 positive regulation of B cell activation, GO:0006911 phagocytosis, engulfment, GO:0038096 Fc-gamma receptor signaling pathway involved in phagocytosis, GO:0006910 phagocytosis, recognition, GO:0050776 regulation of immune response, GO:0050853 B cell receptor signaling pathway, GO:0045087 innate immune response, GO:0006898 receptor-mediated endocytosis, GO:0006954 inflammatory response, GO:0006935 chemotaxis, GO:0042742 defense response to bacterium, GO:0038095 Fc-epsilon receptor signaling pathway, etc. (Figure 3C). The enriched signaling pathways were the hsa04621 NOD-like receptor signaling pathway, hsa04060 cytokine-cytokine receptor interaction, and hsa04668 TNF signaling pathway (Figure 3D).

Upregulation of Immune-Related Pathways

In view of the above findings, a GSEA study was performed to investigate the change of immune-related pathways in the high-glycolysis group. FDR of <0.05 was considered statistically significant. We found that two immune-related pathways were significantly upregulated, including the hsa04060 cytokine-cytokine receptor interaction and hsa04657 interleukin-17 signaling pathway (Figure 4A–C). Taken together, these findings indicated that high glycolytic activity was related to enhanced immune activity in prostate cancer.

Figure 4 Upregulation of immune-related pathways in the high-glycolysis group. (A) The downregulated and upregulated pathways were obtained through gene set enrichment analysis (GSEA). (B and C) The GSEA charts showed the interleukin-17 (IL-17) signaling pathway and cytokine-cytokine receptor interaction were significantly upregulated.

Relationship Between Glycolytic Activity and Immune Cell Infiltration

To verify the relationship between glycolytic activity and immune function, we compared the degree of enrichment of immune cell types in the low- and high-glycolysis groups. In the high-glycolysis group, the degree of enrichment of CD8⁺ T cells, Tregs, monocytes, M1 macrophages, M2 macrophages, resting DCs, and activated mast cells was significantly higher. In contrast, the degree of enrichment of naïve B cells, plasma cells, and eosinophils in the high-glycolysis group was lower than that in the low-glycolysis group (Figure 5A). Additionally, SLC16A3 expression was positively correlated with the proportion of resting DCs (r=0.147, P=0.001), M1 macrophages (r=0.164, P<0.001), M2 macrophages (r=0.132, P=0.003), activated mast cells (r=0.116, P=0.009), monocytes (r=0.211, P<0.001), CD8⁺ T cells (r=0.243, P<0.001), and Tregs (r=0.284, P<0.001), but negatively correlated with that of naïve B cells (r=−0.205, P<0.001), plasma cells (r=−0.15, P<0.001), and eosinophils (r=−0.146, P=0.001) (Figure 5B). To distinguish tumor-promoting from antitumor immune cells, survival analysis was carried out to investigate the relationship between the survival time and the enrichment fraction of immune cells. The results showed that M2 macrophages were highly enriched in the high-glycolysis group, which predicted poor disease-free survival in prostate cancer (Figure 5C).

Figure 5 Relationship between immune cell infiltration and glycolytic activity. (A) Violin plot showed the degree of enrichment of tumor-infiltrating immune cell subtypes between the low- and high-glycolysis groups. (B) Spearman correlation test was used to analyze the correlation between SLC16A3 expression and immune cell enrichment score. The Spearman correlation coefficients between SLC16A3 expression and naïve B cells, resting dendritic cells, eosinophils, M1 macrophages, M2 macrophages, activated mast cells, monocytes, plasma cells, CD8⁺ T cells, and regulatory T cells (Tregs) were remarkable. (C) Kaplan-Meier analysis showed that the group with highly enriched M2 macrophages had poorer disease-free survival rates.

Immunohistochemistry Results

Considering that the high SLC16A3 expression was positively correlated with a higher Gleason grade, we evaluated the protein expression of MCT4 in human prostate cancer tissues. MCT4 expression was higher in tissues with a Gleason grade of 4–5 than in other tissues (Figure 6A). According to the expression score of MCT4, we defined the high- and low-expression groups in prostate cancer tissues as the high-glycolysis group and low-glycolysis group, respectively. Meanwhile, ENO2 protein expression was higher in the high-glycolysis group than in the low-glycolysis group (Figure 6B). However, there was no difference in PFKFB3 expression between both groups (Figure 6C). Immunohistochemical analysis of CD163 showed that M2 macrophages were mainly located in the stroma near the tumor cells. The number of M2 macrophages in prostate cancer was significantly higher in the high-glycolysis group than in the low-glycolysis group (Figure 6D).

Figure 6 Immunohistochemistry results of MCT4, ENO2, PFKFB3, and CD163. (A) Representative MCT4 immunohistochemical staining patterns in prostate cancer tissues and distribution of MCT4 immunohistochemical staining score in prostate cancer tissues according to Gleason grade. (B) Representative ENO2 immunohistochemical staining patterns in prostate cancer tissues and distribution of ENO2 immunohistochemical staining score in prostate cancer tissues according to glycolysis level. (C) Representative PFKFB3 immunohistochemical staining patterns in prostate cancer tissues. (D) Representative CD163 immunohistochemical staining patterns in prostate cancer tissues with red circles demonstrating the selection of CD163⁺ cells (M2 macrophages) under low magnification. Box plots showed the number of M2 macrophages in the low- and high-glycolysis groups in the entire prostate cancer tissue specimen.

Western Blot and RT-qPCR Results

Furthermore, we conducted verification experiments on PC3 cell lines. We successfully constructed the si-NC and si-MCT4 PC3 cells, which were confirmed by RT-qPCR, Western blot and fluorescence transfection experiments (Figure 7A, Supplementary Figure S1). The results (Figure 7A) showed that the gene and protein expression level of MCT4 in PC3 cells in si-MCT4 group was both significantly lower than that in NC group and si-NC group (all P<0.001). The results suggested that the expression of MCT4 in PC3 cells decreased significantly after si-MCT4 transfection, and the PC3 cells silencing MCT4 were successfully established. The results of Western blot also showed that protein expression levels of ENO2 and PFKFB3 in PC3 cells in si-MCT4 group were significantly lower than those in si-NC group and NC group (P<0.001 and P<0.01, respectively). The results suggested that the inhibition of MCT4 could reduce the glycolysis level of PC3 cells to some extent.

Figure 7 The relative expressions of MCT4, ENO2 and PFKFB3 in PC3 cells and the relative gene expressions of Arg1, CD206, IL-6 and iNOS in macrophages. (A) RT-qPCR analysis of MCT4 gene expression of PC3 cells in NC group, si-NC group and si-MCT4 group. Western blot analysis of MCT4, ENO2 and PFKFB3 protein expressions of PC3 cells in NC group, si-NC group and si-MCT4 group. (B) Bar plots showed the gene expressions of Arg1, CD206, IL-6 and iNOS measured by RT-qPCR of macrophages in M0, M1, M2, PC3-TAM and si-MCT4-PC3-TAM groups.The data are presented as the mean ± SD values, and the error bars represent data from triplicate biological experiments. ns, P≥0.05; **P<0.01; ***P<0.001.

In order to determine whether PC3 cells and si-MCT4-PC3 cells could change the phenotype of co-cultured macrophages, we constructed M1 and M2 macrophages for further experiments. As seen in Figure 7B, the gene expressions of anti-inflammatory markers Arg1 and CD206 in macrophages induced by IL-4 were significantly increased (both P<0.001), while the gene expressions of pro-inflammatory markers IL-6 and iNOS were significantly increased in macrophages induced by IFN-γ/LPS compared with M0 macrophages (P<0.001 and P<0.01, respectively). Moreover, we cultured PC3 cells and RAW246.7 macrophages together in a transwell plate (Supplementary Figure S2). After macrophages were co-cultured with PC3 cells, the gene expressions of Arg1 and CD206 in macrophages increased significantly compared with M0 macrophages (both P<0.001), but not IL-6 or iNOS, which indicated that PC3 cells promoted M2 polarization of macrophages. Moreover, when we silenced the expression of MCT4 in PC3 cells, the gene expressions of CD206 and Arg1 in si-MCT4-PC3-TAM were significantly downregulated compared with PC3-TAM (both P<0.001). These results suggested that the expression of MCT4 in PC3 cells played an important role in inducing M2 phenotype of macrophages.

Discussion

Tumorigenesis depends on the reprogramming of cell metabolism, which is characterized by a markedly increased consumption of glucose to generate energy by glycolysis, even under conditions of sufficient oxygen. The metabolic state of tumor cells also affects other cells in the TME, such as tumor-related fibroblasts and immune cells. Metabolites derived from tumor cells can affect the activation and function of immune cells, which ultimately promote the survival, proliferation and resistance to anti-tumor therapy of tumor cells.³³ For example, in prostate cancer, circular RNA could promote tumor progression and radioresistance by upregulating glycolysis.³⁴ Therefore, it will be a promising treatment to reverse the metabolic state of tumor with drugs. Studies have shown that sulforaphane can prevent the occurrence of prostate cancer by inhibiting glycolysis and fatty acid metabolism.^35,36 However, the search for key targets in the metabolic pathway may contribute more to precision treatment of prostate cancer. In this study, we demonstrated that SLC16A3, a key enzyme in the glycolysis pathway, was positively correlated with a higher Gleason grade and worse disease-free survival. The immunohistochemistry results also revealed that increased expression of MCT4, a plasma membrane transporter protein encoded by SLC16A3, correlated with a higher Gleason grade in prostate cancer tissues. The upregulation of glycolysis increases the levels of intracellular lactate and H⁺ ions, thus reducing the intracellular pH of tumor cells. To maintain the dynamic balance of pH and avoid the inhibition of glycolysis by negative feedback, lactate and H⁺ ions are excreted from the cytoplasm of tumor cells through monocarboxylate transporters, which leads to extracellular environment acidification.³⁷ The efflux of cellular lactate/H⁺ ions is thought to be mediated mainly by MCT4.³⁸ It has been demonstrated that MCT4 is upregulated in tumors and is associated with poor prognosis in several types of cancer,^39–42 including prostate cancer.^43–46 MCT4 promotes the invasion of prostate cancer cells by regulating the molecules related to the VEGF/CD147/MMP9 signaling pathway.⁴⁶ Knockdown of MCT4 inhibits glycolysis and tissue invasion/migration of prostate cancer cells.⁴⁴ MCT4-mediated efflux disturbance usually leads to decreased levels of pyruvic and/or lactate released by glycolysis cancer cells and increased levels of intracellular pyruvate and/or lactate, resulting in cytoplasmic acidification and inhibition of glycolysis.^32,47 The above evidence indicates that SLC16A3 or MCT4 plays an important role in glycolysis. This study also confirmed that the expressions of some glycolysis-related genes were significantly higher in the high SLC16A3 expression group by immunohistochemistry and Western blot. Therefore, in this study, SLC16A3 expression was considered to represent glycolytic activity.

Many studies have been conducted on tumor glycolysis, but it remains unclear how glycolysis in prostate cancer affects its immune regulation. In this study, tumor tissues were divided into high-glycolysis and low-glycolysis groups according to the expression of SLC16A3. The gene expression profiles of the two groups were compared, and functional enrichment analyses were performed. The results demonstrated that in the high-glycolysis group, the expression of many cytokines and inflammatory mediators was upregulated. Further GSEA analysis revealed that the most abundant immune-related pathway in the high-glycolysis group was the IL-17 signaling pathway. IL-17 is a powerful pro-inflammatory cytokine, which has been proven to be closely related to the formation, growth, and metastasis of many types of malignant tumors.^48–51 Chronic IL-17 activity results in the production of a population of pathogenic myeloid cells.⁵² These myeloid cell populations play a key role in promoting the formation of pro-tumor immune microenvironment. In a series of murine cancer models, induction of IL-17-dependent pathogenic myeloid cells was associated with tumor progression, including breast cancer,⁵³ lung cancer,⁵⁰ colon cancer,⁵⁴ and liver cancer.⁵⁵ However, the specific mechanism by which glycolytic tumors affect the IL-17 signaling pathway remains unclear. Some studies have demonstrated that lactate from tumor cells would prompt CD4⁺ T cells to produce large amounts of IL-17, which in turn activates the IL-17 signaling pathway.^56–58 Interestingly, the KEGG analysis also proved that the TNF signaling pathway was enriched in the high-glycolysis group. Previous study has found that IL-17 and TNF-α cooperatively promoted tumor glycolysis by regulating the expression of the glucose transporter SLC2A1 and hexokinase-2.⁵⁹ IL-17 may stabilize the target gene mRNA induced by TNF-α.⁶⁰ For the relationship between tumor glycolysis and IL-17 signaling pathways, further experimental studies are needed.

Many immune cell subgroups are essential for antitumor immunotherapy in the TME. Depending on their function and phenotype, immune cells can promote or hinder antitumor immune responses. Thus far, therapies based on immune checkpoints have not yet achieved ideal results in prostate cancer.⁶¹ Therefore, it is particularly necessary to explore immune cell infiltration in the TME of prostate cancer and its relationship with glycolysis. In this study, we compared the enrichment of immune cells between the low- and high-glycolysis groups. According to the results, the degree of enrichment of CD8⁺ T cells, Tregs, monocytes, M1 macrophages, M2 macrophages, resting DCs, and activated mast cells was significantly higher, while that of naïve B cells, plasma cells, and eosinophils was significantly lower in the high-glycolysis group than in the low-glycolysis group. We found that the increased infiltration of immune cells in the high-glycolysis group did not lead to a good prognosis. It seemed that monocytes, macrophages, DCs, mast cells, and T lymphocytes established an immunosuppressive TME in prostate cancer. The existence of T lymphocytes in the TME is related to the prognosis and clinical significance of many types of cancer. Therefore, the density of CD8⁺ and CD3⁺ T cell infiltration at the tumor core and the invasive margin served as a basis for establishing a scoring system, known as the “Immunoscore”, to evaluate the tumor immune status.⁶² Due to their cytotoxic function, CD8⁺ T cells are considered a key component of antitumor immunity. In many solid tumors, highly activated CD8⁺ T cell infiltration often indicates better prognosis.^63,64 However, the correlation between CD8⁺ T cell infiltration and prognosis is controversial in prostate cancer. A study reported that the high infiltration of CD8⁺ T cell predicted shorter biochemical failure-free survival in prostate cancer patients.⁶⁵ Another study also found that a high proportion of CD8⁺ T cells was significantly associated with poor prognosis after examining the density of CD8⁺ T cells in 285 cases of prostate cancer.⁶⁶ An analysis of gene expression data from 1567 prostate cancer samples with long-term clinical results also demonstrated that high tumor CD8⁺ T cell infiltration was associated with poor distant metastasis-free survival.⁶⁷ In our study, we found that prostate cancer tissues in the high-glycolysis group with worse disease-free survival also had higher CD8⁺ T cell infiltration, which is consistent with the results of previous studies. Evidence has shown that the TME can inhibit CD8⁺ T cell immunity, resulting in exhaustion and dysfunction of CD8⁺ T cells.⁶⁸ Persistent antigen stimulation is the key factor leading to the exhaustion of CD8⁺ T cells. One of the main characteristics of exhaustion is the high level of co-expression of a variety of inhibitory receptors, including PD-1, CTLA-4 and others.⁶⁹ The exhaustion of CD8⁺ T cells finally suppresses its control of tumor. Therefore, high glycolysis of prostate cancer may indicate the exhaustion of CD8⁺ T cells infiltration, but the regulatory mechanism between them remains to be further elucidated.

In the TME, tumor cells with a highly glycolytic phenotype are likely to participate in metabolic competition with immune cells. The high level of lactate secreted by highly glycolytic tumor cells hinders the output of lactate by immune cells, thereby interfering with their metabolism and function.⁷⁰ Moreover, some immune cells may take up lactate, which impairs glycolytic function.⁷¹ On this basis, a number of studies have shown that tumor-derived lactate inhibits the proliferation, motility, and cytotoxic activity of immune cells, which ultimately inhibits immune response to tumor cells.^56,72,73 Our results showed that the enrichment of M2 macrophages in the high-glycolysis group was significantly associated with poor disease-free survival in prostate cancer. Immunohistochemical experiment further confirmed that there were large amounts of M2 macrophage infiltration in prostate cancer tissue of high-glycolysis group. Moreover, we found that PC3 cells could promote M2 polarization of co-cultured macrophages and once MCT4 was knocked down, this function of PC3 cells was significantly weakened. This suggested that high glycolytic prostate cancer cells played an important role in promoting the phenotype induction of macrophages in TME. Previous studies have shown that lactate can induce phenotypic changes in macrophages, turning them into a tumor-promoting phenotype (M2-like). Lactate produced by bladder cancer cells induced M2 polarization of TAMs, while blockade of lactate flux between cancer cells and macrophages could inhibit the acquisition of this phenotype.⁷⁴ Similarly, in prostate cancer, the acidic TME enhances a series of functions associated with the M2 polarization of TAMs.²⁹ Although no further experiments were conducted to verify the effect of MCT4 expression in macrophages on M2 polarization in this article, an interesting study found that M2 macrophage polarization mainly depended on active oxidative phosphorylation rather than glycolysis.⁷⁵ We hypothesize that highly glycolytic prostate cancer cells export lactate to the TME through MCT4, which in turn inhibits glycolysis of macrophages. Macrophages then acquire energy through oxidative phosphorylation, which eventually induces M2 polarization. Even though highly glycolytic tumor cells are prone to cause acidity and immunosuppression in the TME,¹⁰ further in vivo and in vitro experiments are still required to verify the mechanism of interaction between tumor glycolysis and TIICs. Our results based on public databases and bioinformatics should be interpreted with caution.

In summary, tumor glycolysis is closely related to the infiltration of immune cells in the TME. Since highly glycolytic tumors have higher immune cell infiltration, targeting tumor glycolysis may help to promote the efficacy of tumor immunotherapy by changing the function of immune cells and enhancing their antitumor immune response.⁷⁶

Furthermore, in view of the importance of MCT4 in glycolysis and tumorigenesis, MCT4 can serve as a potential therapeutic target for prostate cancer. Currently, there is no effective MCT4 inhibitor for treating malignant tumors in clinical applications. However, through the use of siRNA and antisense oligonucleotides for RNA interference experiments, it has been confirmed that inhibiting the expression of MCT4 can reduce tumor growth.^41,44,77 Therefore, the development of selective targeted inhibitors of MCT4 may become a new therapeutic strategy for malignant tumor.

Conclusion

This study revealed that the enhancement of glycolytic activity in prostate cancer was related to increased immune cell infiltration. The IL-17 signaling pathway may play an important regulatory role in the interaction between immune cell infiltration and tumor glycolysis. Infiltration of tumor-promoting immune cells can suppress antitumor immunity in prostate cancer with high level of glycolysis. Therefore, drugs targeting glycolysis may have the effect of inhibiting tumor growth. As an important enzyme in the glycolysis pathway, MCT4 is responsible for exporting lactate from tumor cells, thereby creating an acidic TME that may influence immune cell infiltration and function. Therefore, MCT4 plays a vital role in the development of prostate cancer, which made it a potential therapeutic target. However, since most of the results in this article were based on public database through bioinformatics, the interpretation of the results requires further experiments to verify.

Data Sharing Statement

All data reported in the manuscript are from the public database, and the website can be found in Materials and methods. There are no additional unpublished data from the study.

Acknowledgments

Thanks are due to Professor Fu for assistance with the experiments.

Disclosure

The authors report no conflicts of interest in this work.

Additional information

Funding

The study was sponsored by the Science and Technology Innovation Commission Foundation of Shenzhen (Grant No. JCYJ20190808141013454 and JCYJ20180305124827261).