Skip to main content

A 4-gene-based hypoxia signature is associated with tumor immune microenvironment and predicts the prognosis of pancreatic cancer patients



Pancreatic cancer (PAC) is one of the most devastating cancer types with an extremely poor prognosis, characterized by a hypoxic microenvironment and resistance to most therapeutic drugs. Hypoxia has been found to be one of the factors contributing to chemoresistance in PAC, but also a major driver of the formation of the tumor immunosuppressive microenvironment. However, the method to identify the degree of hypoxia in the tumor microenvironment (TME) is incompletely understood.


The mRNA expression profiles and corresponding clinicopathological information of PAC patients were downloaded from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database, respectively. To further explore the effect of hypoxia on the prognosis of patients with PAC as well as the tumor immune microenvironment, we established a hypoxia risk model and divided it into high- and low-risk groups in line with the hypoxia risk score.


We established a hypoxia risk model according to four hypoxia-related genes, which could be used to demonstrate the immune microenvironment in PAC and predict prognosis. Moreover, the hypoxia risk score can act as an independent prognostic factor in PAC, and a higher hypoxia risk score was correlated with poorer prognosis in patients as well as the immunosuppressive microenvironment of the tumor.


In summary, we established and validated a hypoxia risk model that can be considered as an independent prognostic indicator and reflected the immune microenvironment of PAC, suggesting the feasibility of hypoxia-targeted therapy for PAC patients.


Pancreatic cancer (PAC) is one of the most aggressive and lethal malignancies. In data released by the American Cancer Society in 2019, PAC has a 5-year overall survival rate of only 9% and is one of the tumors with the lowest survival rate [1]. Despite the major advances in surgical techniques and adjuvant medical therapy, it is clear that the survival rate of PAC has not improved significantly. Meanwhile, more and more studies have shown that the high mortality and poor prognosis of PAC are closely correlated to changes in the TME [2]. Recently, bioinformatics-based therapy has become a promising tool in modern oncology. As such, it is pressing to apply to PAC and identify potential therapeutic targets.

Hypoxia is an inherent feature in the microenvironment of solid tumors, which can promote tumor survival and also lead to tumor proliferation and metabolism [3]. It is well-known that PAC possesses a complex TME that comprises myofibroblasts and immune cells [4]. Due to the strong association between hypoxia and fibrosis, PAC is also considered to be the most hypoxic tumor. The highly hypoxic environment produces a large extent of alterations in the tissue structure and cells of PAC in the TME [5]. For example, hypoxia is an important activator of pancreatic stellate cells (PSCs), which not only exacerbates the deterioration of the tumor hypoxic microenvironment [6] but also reduces the migration of natural killer (NK) cells and CD4+ and CD8+ T cells in the tumor stroma and accelerates the differentiation of myeloid-derived immunosuppressive cells to disrupt the balance of tumor immune microenvironment [7].

It has been found that hypoxia was served as a primary factor in the formation of tumor immunosuppressive microenvironment, which enhances tumor immune evasion by suppressing the anti-tumor immune responses [8]. Generally, hypoxia promotes the development of immunosuppressive cell populations (MDSC, Treg cells, M2-like macrophages, and immunosuppressive cytokines) in the microenvironment but also reduces the killing, survival, and migration of effector cells ( NK cells, CD4+ and CD8+ T cells), thereby impairing the regulation of anti-tumor immunity [9]. In recent years, many clinical trials have been performed on immunotherapy for tumors, but the efficacy of clinical treatment achieved on PAC are limited [10]. Although the specific mechanism of PAC immunotherapy resistance is still unclear, we believe that it is inseparable from the “barren soil” caused by hypoxia.

In the present study, we downloaded the mRNA profile data of PAC patients from the TCGA and GEO database and extracted genetic data related to hypoxia to construct a hypoxia risk model. We first conducted prognostic prediction for patients with PAC based on this risk model. On this basis, we further searched for the interrelation between hypoxia and tumor immune microenvironment. In the future, this approach will be crucial for researchers to develop new combination treatment strategies.


Data acquisition

Clinical information and RNA-sequencing expression date of 186 PAC patients were collected from the TCGA ( as a training set. Subsequently, the corresponding information of 112 PAC patients from the GSE78229 and GSE57495 datasets were downloaded from the GEO database as a validation set ( The two chip datasets in GEO were combined for batch normalization, and the combined result was used for further analysis. Tumor immune gene set and hallmark gene set were downloaded from TIP ( and GSEA ( respectively. All P values were subjected to multiple testing, taking an adjusted P value < 0.05 as a threshold.

Construction of protein-protein interaction (PPI) network

String Database was used to construct the PPI network for downloaded hypoxia-related genes. R software language was utilized to analyze the correlation of hypoxia-related genes in the protein interaction relationship network, and then screened them as key core genes based on the number of interrelationships.

Constitution of a risk model

First, we performed univariable Cox regression analysis on the hypoxia core genes to obtain prognostic-related hypoxia genes. They were then analyzed by multivariable Cox regression to obtain the genes for building the model and their coefficients. The risk score of each patient was calculated according to the obtained genes, and then patients were divided into low and high hypoxia risk groups based on the median risk score. The risk score formula was constructed as:

$$ \mathrm{Risk}\ \mathrm{score}={\sum}_{i=1}^{\mathrm{n}}\left({Exp}_i\ast {Coe}_i\right) $$

where N = 4, Expi indicated the expression level for each hypoxia genes, and the Coei indicated the corresponding multivariable Cox regression coefficient.

Survival analysis

Overall survival (OS) analysis of PAC patients with low- and high- risk groups was performed by Kaplan-Meier method employing survminer and survival packages in R. Univariate and multivariate COX regression analysis was carried out to determine whether the risk score can be distinguished from other conventional clinical features as an independent prognostic factor for PAC patients. Bilateral P values less than 0.05 were considered statistically significant, and the hazard ratio (HR) was calculated for 95% confidence intervals. A ROC curve created by the survivalROC R package was used to evaluate the accuracy and reliability of risk model for predicting the patients’ OS.

To better assess the 1-, 3-, and 5-year survival probabilities of patients with PAC, we utilized the total independent prognostic factors to construct a nomogram.

Gene set enrichment analysis (GSEA)

GSEA was carried out in the enrichment of the MSigDB Collection (h.all.v7.0.cymbols.gmt) to determine the regulation of signaling pathways in the hypoxia gene set between low and high hypoxia risk groups. Gene set permutations were conducted 1000 times for each analysis. P < 0.05 was considered statistically significant.

Evaluation of immune cell type components

CIBERSORT is a tool for evaluating the proportion of various cell subtypes from a mixed cell samples through RNA-seq expression profiles and is the commonly cited method to estimate and analyze the immune cell infiltration [11]. We used CIBERSORT to evaluate the proportion of 22 immune cell subtypes in the high- and low-risk group, including CD8 T cells, resting memory CD4 T cells, and resting NK cells. The sum of the scores for the total immune cell types in a sample is equal to 1.


Characterization of hypoxia risk score to predict PAC prognosis

We downloaded the hypoxia-associated gene set from the GSEA (hallmark-hypoxia), which contains more than 200 genes involved in hypoxia regulation pathways. To further understand the relationship between these hypoxia-related genes, we utilized the string online data ( to conduct a network analysis of protein-protein interactions for these genes, and then extracted the most interconnected genes as core genes (Fig. 1a). Results showed that the 50 genes with the highest interaction levels were maintained, including GAPDH, VEGFA, IL6, EGFR, JUN, HK1, LDHA, ENO1, PGK1, ALDOA, indicating their crucial role in the hypoxia regulation process (Fig. 1b).

Fig. 1

Characterization of hypoxia risk score to predict prognosis of PAC. a Protein–protein interactions among 200 hypoxia-associated genes. b The 50 genes with the highest interaction degrees were extracted. c Sixteen hypoxia genes associated with patient’ OS were confirmed by univariate Cox regression. d Four hypoxia-related genes were chosen to establish a hypoxia risk model by multivariate Cox regression. e, f Spearman correlation analysis of four hypoxia genes in the TCGA and GEO databases

To build a hypoxia risk model that can predict the prognosis of PAC patients, we carried out univariate and multivariate Cox regression analysis on 50 hypoxia core genes in the TCGA training dataset (Tables 1 and 2). In the univariate Cox regression analysis, we confirmed 16 hypoxia genes associated with patient’ OS (Fig. 1c). In the following multivariate Cox regression analysis, the final 4 genes with a P value < 0.05 were chosen to establish the prognostic model, including PGK1, ENO3, LDHA, and PGM1 (Fig. 1d). The formula to calculate the risk score is defined as follows:

$$ \mathrm{Risk}\ \mathrm{score}=\left(-0.484\times \mathrm{PGK}1\right)+\left(-0.799\times \mathrm{ENO}3\right)+\left(1.026\times \mathrm{LDHA}\right)+\left(0.415\times \mathrm{PGM}1\right) $$
Table 1 A total of 16 hypoxia-related genes significantly associated with overall survival according to univariate Cox regression analysis. Hazard ratio (HR); low 95% confidence interval of hazard ratio (HR.95 L); high 95% confidence interval of hazard ratio (HR.95H)
Table 2 Details of the 4 hypoxia-related genes significantly associated with overall survival used to build the hypoxia risk model. Hazard ratio (HR); low 95% confidence interval of hazard ratio (HR.95 L); high 95% confidence interval of hazard ratio (HR.95H)

There was no significant correlation among the four hypoxia-related genes selected by us (Fig. 1e, f).

Prognostic efficacy of the hypoxia risk score in PAC patients

Studies have found that hypoxia can exacerbate the aggressiveness of tumors, so we further assessed the efficacy of hypoxia risk score on the prognosis of PAC. In the heatmap (Fig. 2a), we found that the expression of three hypoxia-associated genes in the high-risk group was significantly increased for both TCGA and GEO databases, indicating that PAC patients in the high-risk group may be more prone to form a hypoxic TME. Meanwhile, our research data found that hypoxia risk score was increased accompanying higher patient risk level (Fig. 2b), and the mortality rate of patients in the high-risk group was remarkably higher than that in the low-risk group (Fig. 2c, d). Besides, we used Kaplan-Meier analysis to estimate the prognostic value of hypoxia risk score in PAC patients. The results showed that high risk score was obviously associated with bad prognosis in the TCGA cohort, which was further validated from the GEO cohort (Fig. 2e).

Fig. 2

Prognostic efficacy of the hypoxia risk score in PAC. a The expression of four hypoxia-related genes in high and low hypoxia risk groups from the TCGA and GEO databases. b The relationship between hypoxia risk score and patient risk level. c Patient survival status distribution in the high and low hypoxia risk groups. The dot indicates the status of the patient ranked according to the increasing risk score. d The mortality rate of patients in high- and low- risk groups. e Kaplan-Meier overall survival analysis of PAC patients in high- and low-risk groups in the TCGA and GEO datasets. f Nomogram predicting the 1-, 3-, and 5-year survival probability of each individual patient

Moreover, we created a nomogram to further calculate the survival probability of each patient more accurately. Set the scoring criteria according to the regression coefficient of all independent factors, and then the score value of each independent factor was given to calculate the total scores in each individual. The 1-, 3-, and 5-year survival probabilities of each individual was achieved via the function conversion relationship of total scores (Fig. 2f).

Hypoxia risk score shows great feasibility for prognosis evaluation

We carried out the received operating characteristic (ROC) curve to further assess the predictive accuracy and reliability of hypoxia risk score for the survival rate. In the TCGA cohort, the area under ROC curve (AUC) was 0.701 at 1 year, 0.758 at 3 years, and 0.884 at 5 years, respectively, demonstrating the predictive power of our risk model (Fig. 3a). We further validated this result in an independent cohort using GEO (Fig. 3b).

Fig. 3

Prognostic value of the hypoxia risk score in PAC. a, b ROC curve showing the prognostic value of hypoxia risk scores in patients. cf Univariate and multivariate Cox analyses evaluating the independent prognostic value of the hypoxia risk score in terms of OS in PAC patients

We applied univariate and multivariate Cox regression analysis to examine whether the hypoxic risk score could be an independent prognostic factor for predicting survival in patients with PAC. First, univariate Cox regression analysis showed that tumor grade, N stage, age, and hypoxia risk score were associated with OS in TCGA cohort (Fig. 3c). Then, multivariate Cox regression analysis indicated that hypoxia risk score was independently correlated with unfavorable overall survival of PAC patients, which could act as an independent prognostic factor in PAC (Fig. 3d). Similarly, validation result of the GEO cohort also demonstrated that hypoxia risk score can be used as an independent prognostic factor (Fig. 3e, f).

GSEA identifies hypoxia signaling pathways

We carried out the GSEA to further study signaling pathways activated by hypoxia-related genes in the high-risk group. The results showed that more genes in the high-risk group of the TCGA dataset were significantly enriched in multiple pathways such as hypoxia, TGF-β signaling, epithelial-mesenchymal transition, and mTORC1 signaling (Fig. 4a). They were all associated with the processes including tumor metastasis, proliferation, and anti-apoptosis. This result was validated by GEO datasets, which completely echoed the result (Fig. 4b).

Fig. 4

GSEA enrichment between high and low hypoxia risk groups. a GSEA indicating that genes in the high hypoxia risk group were enriched for the characteristics of malignant tumors. b This result was validated by GEO datasets. Normalized enrichment score (NES) > 1 and nominal p value (NOM p-val) < 0.05 were considered significant gene sets

Analysis of immune cells in the high and low hypoxia risk group for PAC

Increasing studies indicated that hypoxia plays an integral role in tumor immune tolerance. It promotes the formation of the immunosuppressive microenvironment by means of abnormal immune cell activation, secretion of immunosuppressive factors, and downregulation of immune cells to protect tumor cells from immune cell attack.

The CIBERSORT method was utilized to analyze the infiltration of 22 subpopulations of immune cells between the low- and high-risk groups of PAC patients. The results obtained from 186 PAC patients in TCGA and 112 patients in GEO were summarized in Fig. 5a. It was demonstrated that the proportion of immunosuppressive cells (M2 macrophages), resting memory CD4+ T cells, and resting NK cells (Fig. 5b–d) was significantly higher in the high hypoxia risk group. Besides, the levels of CD8 T cells, plasma cells, and naive B cells (Fig. 5e–g) were significantly reduced in PAC patients with high-risk score.

Fig. 5

Analysis of immune cells in the high and low hypoxia risk group for PAC. a Relative proportion of immune infiltration in high and low hypoxia risk patients. bg Box plot visualizing the differentially infiltrated immune cells in the low- and high- risk patients.

Thus, it can be seen that targeting hypoxia research is very important for future immunotherapy for tumor patients.

High hypoxia risk tended to an immunosuppressive microenvironment

Tumor immunotherapy is a therapy that restores and enhances the host’s immune system to recognize and eradicate cancer. Immunotherapy involves a series of immune-related cell mediated tumor killing processes. Next, we examined the expression level of genes unfavorably regulating these processes in high- and low-risk groups. We then downloaded gene signatures on the Tracking Tumor Immunophenotype website ( [12].

The results showed that the expressions of VEGFA, HAVCR1, CXCL8, MICB, and ICAM1, which were positively associated with hypoxia risk score, were significantly higher in the high-risk group (Fig. 6a–g), and the analysis found that most of them were related to the regulation of HIF-1α. For the present study, most of the immunosuppressive cytokines and immune checkpoints were not found to present a significant increasing trend in the high hypoxia risk group. We thought this might indicate that they are not the main drivers of the hypoxia-induced immunosuppressive microenvironment in PAC.

Fig. 6

High hypoxia risk tended to an immunosuppressive microenvironment. a, b The heatmap displaying the expression of the gene set involved in the negative regulation of anti-tumor immunotherapy in the low and high hypoxia risk groups in the TCGA and GEO databases. cg The expression levels of VEGFA, CXCL8, ICAM1, MICB, and HAVCR1 in high and low hypoxia risk groups. Correlation between the expression of VEGFA, CXCL8, ICAM1, MICB, and HAVCRI and hypoxia risk score. hj The expression of genes that positively regulating T cells, dendritic cells, and MDSC cells in the low and high hypoxia risk groups.*P < 0.05, **P < 0.01, and ***P < 0.001

We further analyzed the expression of genes that positively regulating T cells, dendritic cells, and Myeloid-derived suppressor cells (MDSCs). As seen in Fig. 6h, i, the expression of regulatory genes of most anti-tumor immune effector cell was obviously downregulated in the high hypoxia risk group, including T cell regulatory genes CCL21, CXCR5, and CCR7 and dendritic cell regulatory genes CCL21, CCL4, and CCR7. Conversely, MDCS cell regulatory genes CXCL5 and CXCR2 were significantly upregulated in the high-risk group (Fig. 6j). This further confirms that the activity of effector cells in the immune microenvironment of PAC patients with high-risk score is reduced, while the immunosuppressive cell is increased.


Intratumoral hypoxia is a common characteristic of the TME, but also a prominent biological feature for PAC [3]. The hypoxic microenvironment of tumors is considered to be the main mechanism leading to tumor resistance to various treatments such as radiotherapy, chemotherapy, and immunotherapy [13, 14]. Immunotherapy refers to killing tumor cells by improving the antitumor immune responses of tumor patients, which is a promising tumor treatment method. Inhibition of immune checkpoint to enhance immune cell-mediated tumor killing has shown a promising result in many cancers (e.g., melanoma, colorectal cancer), but the effect in some solid tumors including PAC is disappointing [15, 16]. The existence of this phenomenon has been a subject of active debate. However, the potential effect of hypoxia as a key microenvironmental factor on the treatment of PAC is still being explored.

The risk model in our research was constructed by four hypoxia-related genes, and most of them were highly upregulated under hypoxic conditions. It was found that hypoxia-inducible factor 1/2α (HIF1/2α) can activate the expression of LDHA in PAC to achieve the purpose of promoting the proliferation and metastasis of PAC cells [17]. PGK1, a major enzyme in glycolysis, has been reported overexpressed in numerous malignancies. Studies have demonstrated that PGK1 is directly regulated by HIF-1α and acts as a promoter of metastasis in gastric cancer and colon cancer [18, 19]. Similarly, PGM1 was also shown to be upregulated under hypoxic environment [20]. In fact, hypoxia-associated signature for predicting the diagnosis, prognosis, and immune landscape has been used in a variety of tumors. For example, a hypoxia model developed by Zhang et al. could be regarded as a potential biomarker for diagnosis, prognosis, and recurrence of hepatocellular carcinoma [21]. And Mo et al. showed that the hypoxia-related gene signature could be considered as a prognostic factor and may guide the choice of immunotherapy in lung adenocarcinoma [22]. In this study, our hypoxic risk model composed of four genes can also be flexibly applied to PAC’s analysis.

Numerous studies have shown that hypoxia may induce tumor immune escape in a range of ways, including (A) reduce the activity of effector cells (NK cells, CD4+ and CD8+ T cells), (B) decrease the production and release of effector cytokines, (C) support the activity of immunosuppressive cells (Tregs, M2-like macrophage and MDSC), and (D) induce the expression of immunosuppressive cytokines.

The activation of NK cells and cytotoxic T lymphocytes (CTLs) is a crucial step in the tumor immune response. Accumulating evidence suggests that the low infiltration degree of NK cells, CD4+, and CD8+ T cells is a biomarker of poor prognosis and adverse clinical outcomes [23, 24]. However, it has been shown that hypoxia can suppress the activity of effector T cells and NK cells, which could result in decreased immune function. For example, hypoxia will enhance the process of glycolysis and glutamine decomposition, which will cause elevated levels of lactate production in the TME [25]. However, the acidic microenvironment not only inhibits the activity of T cells and the release of cytokines but also prevents NK cells from secreting TNF-α, IFN-y, perforin, and granzyme B [26, 27]. Our data illustrated that for patients in high-risk group, CD8+ T cells were decreased while resting NK cells and resting memory CD4+ T cells were increased, suggesting the formation of tumor immunosuppressive microenvironment.

In a recent study, three independent studies conducted by MD Anderson Cancer Center, INSERM, and Lund University all proposed that the appearance of B cells and tertiary lymphoid structures in tumor tissues is associated with better prognosis of patients undergoing immunotherapy [28,29,30]. The researchers speculate that B cells in the front line of anti-tumor may produce antibodies to effectively fight tumor cells and may also act by supporting CD8+ T cells [28]. Moreover, the study of Lee et al. [31] demonstrated that knockout of the HIF-1α gene in the PAC mouse model resulted in elevated expression levels of B cells chemokine, suggesting that hypoxia reduces the penetration of B cells in tumors, which is also consistent with our findings. In summary, although both anti-tumorigenic and pro-tumorigenic effects of B cells have been reported, it is still an indisputable fact that the result of hypoxia-induced B cell reduction in PAC patients may have a negative impact on the patient’s prognosis and immunotherapy.

TAMs originate from the response of macrophages in the blood to tumor signals and can be divided into M1 type and M2 type. Activation of the M2 phenotype has been proved to promote tumor angiogenesis and tumor cell metastasis through secretion of angiogenic factors (VEGF-A, IL-6, MMP) and immune suppressive factors (IL-10, TGF-β) [32]. Previous studies have shown that hypoxia develops a functionally immunosuppressive microenvironment by stimulating the differentiation of macrophages into M2 type [33]. In this study, we observed a significantly higher proportion of the M2-like phenotype in the high hypoxia risk group, suggesting that this hypoxia risk model has the ability to predict the immune microenvironment.

Furthermore, we analyzed the genes that were significantly increased in the high hypoxia risk group, including VEGFA, CXCL8, HAVCR1, MICB, and ICAM1, and found that they were all related to the regulation of HIF-1α [34,35,36,37]. Studies have found that VEGFA affects the response of effector T cells and the development of lymphocytes by promoting suppressive immune cell populations, which is closely related to hypoxia and immunosuppressive TME [38,39,40]. Meanwhile, Chao et al. [41] found that CXCL8 inhibits the function of CD8+ T cells via increasing the expression of PD-L1 on the surface of macrophages to participate in the formation of the immunosuppressive microenvironment of gastric cancer. However, the immune-related regulatory role of CXCL8 in PAC is unclear. In addition, ICAM-1 and HAVCR-1 are both overexpressed in various cancers and are related to the malignant potential of tumor cells [36, 42]. It is known that HIF-1 is an important regulator of the TME. Due to the specific hypoxic TME of pancreatic cancer, the expression of HIF-1α in PAC tissue is higher than that of other solid tumors [43, 44]. Therefore, HIF-1α seems to be implicated in the construction of an immunosuppressive microenvironment under hypoxia in PAC.

Recently, increasing studies observed that immune checkpoint plays a crucial part in tumor immune evasion and the formation of tumor immunosuppressive microenvironment [45]. Indeed, immune checkpoint blockade therapy has achieved remarkable results in some tumors but has little effect on PAC. Thus, it was proposed that PD1, PDL1, and CTLA4 are likely not the primary immune checkpoint molecules involved in immune suppression of PAC [5]. Consistent with this evidence, our data showed that associated immune checkpoint molecules were not significantly higher in the high hypoxia risk group. Taken together, we speculate that immune checkpoint may not be the major factor contributing to the immunosuppressive microenvironment in PAC under hypoxic conditions.

In this study, we successfully establish a 4-gene-based hypoxia risk model that effectively predicts the prognoses of patients with PAC. Here, our model served as an independent prognostic factor for PAC patients and describes how hypoxia status disrupts the immune microenvironment in PAC. However, it is undeniable that the interaction between hypoxia and the tumor’s immune microenvironment is complex, and more independent cohorts and functional experiments are needed to conduct a deeper discussion. We still believe this study will provide insights into the potential value of hypoxia-targeted therapies and help to design unconventional combinatorial methods to enhance the efficacy of PAC therapies.


By combining bioinformatics tools and related algorithms, we established and validated a hypoxia risk model to predict the prognosis of patients with pancreatic cancer and explored the changes in the tumor immune microenvironment under hypoxic conditions. This study found that patients with high hypoxia risk are associated with poor prognosis and the formation of an immunosuppressive microenvironment, which could bring new insights to tumor treatment.

Availability of data and materials

The data used to support the findings of this study are included within the article.



Pancreatic cancer


Tumor microenvironment


Overall survival


Received operating characteristic


Cytotoxic T lymphocytes


Gene set enrichment analysis

P :

P value


Hazard ratio


Myeloid-derived suppressor cells


Gene Expression Omnibus


The Cancer Genome Atlas


Hypoxia-inducible factor 1/2α


  1. 1.

    Ying H, Dey P, Yao W, Kimmelman AC, Draetta GF, Maitra A, et al. Genetics and biology of pancreatic ductal adenocarcinoma. Genes Dev. 2016;30(4):355–85.

  2. 2.

    Ren B, Cui M, Yang G, Wang H, Feng M, You L, et al. Tumor microenvironment participates in metastasis of pancreatic cancer. Mol Cancer. 2018;17(1):108.

  3. 3.

    Daniel SK, Sullivan KM, Labadie KP, Pillarisetty VG. Hypoxia as a barrier to immunotherapy in pancreatic adenocarcinoma. Clinical and Translational Medicine. 2019;8(1):10.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  4. 4.

    Shibuya KC, Goel VK, Xiong W, Sham JG, Pollack SM, Leahy AM, et al. Pancreatic ductal adenocarcinoma contains an effector and regulatory immune cell infiltrate that is altered by multimodal neoadjuvant treatment. Plos One. 2014;9(5):e96565.

  5. 5.

    Foucher ED, Ghigo C, Chouaib S, Galon J, Iovanna J, Olive D. Pancreatic ductal adenocarcinoma: a strong imbalance of good and bad immunological cops in the tumor microenvironment. Front Immunol. 2018;9:1044.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  6. 6.

    Erkan M, Kurtoglu M, Kleeff J. The role of hypoxia in pancreatic cancer: a potential therapeutic target? Expert Rev Gastroenterol Hepatol. 2016;10(3):301–16.

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    Mace TA, Ameen Z, Collins A, Wojcik S, Mair M, Young GS, et al. Pancreatic cancer-associated stellate cells promote differentiation of myeloid-derived suppressor cells in a STAT3-dependent manner. Cancer Res. 2013;73(10):3007–18.

  8. 8.

    Petrova V, Annicchiarico-Petruzzelli M, Melino G, Amelio I. The hypoxic tumour microenvironment. Oncogenesis. 2018;7(1):10.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  9. 9.

    Chouaib S, Noman MZ, Kosmatopoulos K, Curran MA. Hypoxic stress: obstacles and opportunities for innovative immunotherapy of cancer. Oncogene. 2017;36(4):439–45.

    CAS  Article  PubMed  Google Scholar 

  10. 10.

    Duffy MJ, Crown J. Biomarkers for predicting response to immunotherapy with immune checkpoint inhibitors in cancer patients. Clin Chem. 2019;65(10):1228–38.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, et al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. 2015;12(5):453–7.

  12. 12.

    Xu L, Deng C, Pang B, Zhang X, Liu W, Liao G, et al. TIP: a web server for resolving tumor immunophenotype profiling. Cancer Res. 2018;78(23):6575–80.

  13. 13.

    Harris AL. Hypoxia--a key regulatory factor in tumour growth. Nat Rev Cancer. 2002;2(1):38–47.

    CAS  Article  PubMed  Google Scholar 

  14. 14.

    Shannon AM, Bouchier-Hayes DJ, Condron CM, Toomey D. Tumour hypoxia, chemotherapeutic resistance and hypoxia-related therapies. Cancer Treat Rev. 2003;29(4):297–307.

    CAS  Article  PubMed  Google Scholar 

  15. 15.

    Ahn DH, Ramanathan RK, Bekaii-Saab T. Emerging therapies and future directions in targeting the tumor stroma and immune system in the treatment of pancreatic adenocarcinoma. Cancers (Basel). 2018;10(6).

  16. 16.

    McCormick KA, Coveler AL, Rossi GR, Vahanian NN, Link C, Chiorean EG. Pancreatic cancer: update on immunotherapies and algenpantucel-L. Hum Vaccin Immunother. 2016;12(3):563–75.

    Article  PubMed  Google Scholar 

  17. 17.

    Cui XG, Han ZT, He SH, Wu XD, Chen TR, Shao CH, et al. HIF1/2alpha mediates hypoxia-induced LDHA expression in human pancreatic cancer cells. Oncotarget. 2017;8(15):24840–52.

  18. 18.

    Ahmad SS, Glatzle J, Bajaeifer K, Buhler S, Lehmann T, Konigsrainer I, et al. Phosphoglycerate kinase 1 as a promoter of metastasis in colon cancer. International Journal of Oncology. 2013;43(2):586–90.

  19. 19.

    Zieker D, Konigsrainer I, Tritschler I, Loffler M, Beckert S, Traub F, et al. Phosphoglycerate kinase 1 a promoting enzyme for peritoneal dissemination in gastric cancer. International Journal of Cancer. 2010;126(6):1513–20.

  20. 20.

    Pelletier J, Bellot G, Gounon P, Lacas-Gervais S, Pouyssegur J, Mazure NM. Glycogen synthesis is induced in hypoxia by the hypoxia-inducible factor and promotes cancer cell survival. Front Oncol. 2012;2:18.

    Article  Google Scholar 

  21. 21.

    Zhang B, Tang B, Gao J, Li J, Kong L, Qin L. A hypoxia-related signature for clinically predicting diagnosis, prognosis and immune microenvironment of hepatocellular carcinoma patients. J Transl Med. 2020;18(1):342.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Mo Z, Yu L, Cao Z, Hu H, Luo S, Zhang S. Identification of a hypoxia-associated signature for lung adenocarcinoma. Front Genet. 2020;11:647.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  23. 23.

    Dehghani M, Sharifpour S, Amirghofran Z, Zare HR. Prognostic significance of T cell subsets in peripheral blood of B cell non-Hodgkin's lymphoma patients. Med Oncol. 2012;29(4):2364–71.

    Article  PubMed  Google Scholar 

  24. 24.

    Sanchez-Martinez D, Allende-Vega N, Orecchioni S, Talarico G, Cornillon A, Vo DN, et al. Expansion of allogeneic NK cells with efficient antibody-dependent cell cytotoxicity against multiple tumors. Theranostics. 2018;8(14):3856–69.

  25. 25.

    Romero-Garcia S, Moreno-Altamirano MM, Prado-Garcia H, Sanchez-Garcia FJ. Lactate contribution to the tumor microenvironment: mechanisms, effects on immune cells and therapeutic relevance. Front Immunol. 2016;7:52.

    Article  Google Scholar 

  26. 26.

    Calcinotto A, Filipazzi P, Grioni M, Iero M, De Milito A, Ricupito A, et al. Modulation of microenvironment acidity reverses anergy in human and murine tumor-infiltrating T lymphocytes. Cancer Res. 2012;72(11):2746–56.

  27. 27.

    Potzl J, Roser D, Bankel L, Homberg N, Geishauser A, Brenner CD, et al. Reversal of tumor acidosis by systemic buffering reactivates NK cells to express IFN-gamma and induces NK cell-dependent lymphoma control without other immunotherapies. International Journal of Cancer. 2017;140(9):2125–33.

  28. 28.

    Cabrita R, Lauss M, Sanna A, Donia M, Skaarup Larsen M, Mitra S, et al. Author correction: tertiary lymphoid structures improve immunotherapy and survival in melanoma. Nature. 2020;580(7801):E1.

  29. 29.

    Helmink BA, Reddy SM, Gao J, Zhang S, Basar R, Thakur R, et al. B cells and tertiary lymphoid structures promote immunotherapy response. Nature. 2020;577(7791):549–55.

  30. 30.

    Petitprez F, de Reynies A, Keung EZ, Chen TW, Sun CM, Calderaro J, et al. B cells are associated with survival and immunotherapy response in sarcoma. Nature. 2020;577(7791):556–60.

  31. 31.

    Lee KE, Spata M, Bayne LJ, Buza EL, Durham AC, Allman D, et al. Hif1a deletion reveals pro-neoplastic function of B cells in pancreatic neoplasia. Cancer Discov. 2016;6(3):256–69.

  32. 32.

    Van Woensel M, Mathivet T, Wauthoz N, Rosiere R, Garg AD, Agostinis P, et al. Sensitization of glioblastoma tumor micro-environment to chemo- and immunotherapy by Galectin-1 intranasal knock-down strategy. Sci Rep. 2017;7(1):1217.

  33. 33.

    Henze AT, Mazzone M. The impact of hypoxia on tumor-associated macrophages. J Clin Invest. 2016;126(10):3672–9.

    Article  PubMed  PubMed Central  Google Scholar 

  34. 34.

    Azoitei N, Becher A, Steinestel K, Rouhi A, Diepold K, Bobrovich S, et al. PKM2 promotes tumor angiogenesis by regulating HIF-1alpha through NF-kappaB activation. Cancer Research. 2017;77.

  35. 35.

    Li XP, Yang XY, Biskup E, Zhou J, Li HL, Wu YF, et al. Co-expression of CXCL8 and HIF-1alpha is associated with metastasis and poor prognosis in hepatocellular carcinoma. Oncotarget. 2015;6(26):22880–9.

  36. 36.

    Telford EJ, Jiang WG, Martin TA. HAVcR-1 involvement in cancer progression. Histol Histopathol. 2017;32(2):121–8.

    Article  PubMed  Google Scholar 

  37. 37.

    Chen Z, Zhao L, Zhao F, Yang G, Wang JJ. Tetrandrine suppresses lung cancer growth and induces apoptosis, potentially via the VEGF/HIF-1alpha/ICAM-1 signaling pathway. Oncol Lett. 2018;15(5):7433–7.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  38. 38.

    Colegio OR, Chu NQ, Szabo AL, Chu T, Rhebergen AM, Jairam V, et al. Functional polarization of tumour-associated macrophages by tumour-derived lactic acid. Nature. 2014;513(7519):559–63.

  39. 39.

    Palazon A, Aragones J, Morales-Kastresana A, de Landazuri MO, Melero I. Molecular pathways: hypoxia response in immune cells fighting or promoting cancer. Clin Cancer Res. 2012;18(5):1207–13.

    CAS  Article  PubMed  Google Scholar 

  40. 40.

    Terme M, Pernot S, Marcheteau E, Sandoval F, Benhamouda N, Colussi O, et al. VEGFA-VEGFR pathway blockade inhibits tumor-induced regulatory T-cell proliferation in colorectal cancer. Cancer Res. 2013;73(2):539–49.

  41. 41.

    Lin C, He HY, Liu H, Li RC, Chen YF, Qi YY, et al. Tumour-associated macrophages-derived CXCL8 determines immune evasion through autonomous PD-L1 expression in gastric cancer. Gut. 2019;68:1764.

  42. 42.

    Usami Y, Ishida K, Sato S, Kishino M, Kiryu M, Ogawa Y, et al. Intercellular adhesion molecule-1 (ICAM-1) expression correlates with oral cancer progression and induces macrophage/cancer cell adhesion. Int J Cancer. 2013;133(3):568–78.

  43. 43.

    Joshi S, Kumar S, Ponnusamy MP, Batra SK. Hypoxia-induced oxidative stress promotes MUC4 degradation via autophagy to enhance pancreatic cancer cells survival. Oncogene. 2016;35(45):5882–92.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  44. 44.

    Wang G, Wang JJ, Fu XL, Guang R, To ST. Advances in the targeting of HIF-1alpha and future therapeutic strategies for glioblastoma multiforme (Review). Oncol Rep. 2017;37(2):657–70.

    CAS  Article  PubMed  Google Scholar 

  45. 45.

    Pardoll DM. The blockade of immune checkpoints in cancer immunotherapy. Nat Rev Cancer. 2012;12(4):252–64.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

Download references


Thanks to the Fund (The Key Research and Development Plan Projects of Anhui Province) for supporting this research.


This work was supported by grants from the key research and development program of Anhui Province (201904a07020045).

Author information




Jianfeng Ding and Xiaobo He wrote the manuscript and methodology; Xiao Cheng and Guodong Cao processed the data analysis; Bo Chen organized and produced pictures; Sihan Chen and Maoming Xiong guided and revised the manuscript. All authors contributed equally to this manuscript. The author(s) read and approved the final manuscript.

Corresponding authors

Correspondence to Sihan Chen or Maoming Xiong.

Ethics declarations

Ethics approval and consent to participate

There were no cell, tissue, or animal studies. No ethical requirements are involved.

Consent for publication

Not applicable.

Competing interests

The authors declare that no competing or conflicting interests exist.

Additional information

Publisher’s Note

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

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Ding, J., He, X., Cheng, X. et al. A 4-gene-based hypoxia signature is associated with tumor immune microenvironment and predicts the prognosis of pancreatic cancer patients. World J Surg Onc 19, 123 (2021).

Download citation


  • Hypoxia
  • Hypoxia risk model
  • Pancreatic cancer
  • Tumor microenvironment
  • Immune response
  • Prognosis