Yuyang Gu1†
Tongpeng Xu
Zhenghua Fei- 1Department of Oncology, The First Affiliated Hospital of Wenzhou Medical University, Wenzhou, China
- 2Yancheng Third People’s Hospital, The Sixth Affiliated Hospital of Nantong University, Yancheng, China
- 3Department of Oncology, The First Affiliated Hospital of Nanjing Medical University, Nanjing, China
Background: Gastric cancer (GC) is a leading cause of cancer-related deaths worldwide, accounting for high rates of morbidity and mortality in the population. The tumor microenvironment (TME), which plays a crucial role in GC progression, may serve as an optimal prognostic predictor of GC. In this study, we identified CXC motif chemokine receptor 4 (CXCR4) as a TME-related gene among thousands of differentially expressed genes (DEGs). We showed that CXCR4 can be used to predict the effect of immunotherapy in patients with GC.
Methods: GC samples obtained from The Cancer Genome Atlas (TCGA) were analyzed for the presence of stroma (stromal score), the infiltration of immune cells (immune score) in tumor tissues, and the tumor purity (estimate score) using the ESTIMATE (Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data) algorithm. DEGs were sorted based on differences in the values of the three scores. Furthermore, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed to determine the biological processes and pathways enriched in these DEGs. The correlations of scores with clinicopathological features and overall survival (OS) of patients with GC were assessed by the Kaplan–Meier survival and Cox regression analyses. Through subsequent protein–protein interaction (PPI) network and univariate Cox regression analyses, CXCR4 was identified as a TME-related gene. Gene Set Enrichment Analysis (GSEA) was performed to assess the role of CXCR4 in the TME of GC. The CIBERSORT algorithm was used to further explore the correlation between tumor-infiltrating immune cells (TIICs) and CXCR4. Finally, the TISIDB database was used to predict the efficacy of immunotherapy in patients with GC.
Results: We extracted 1231 TME-related DEGs and by an overlapping screening of PPI network and univariate Cox regression, CXCR4 was identified as a biomarker of TME, which deeply engaged in immune-related biological processes of gastric cancer and have close association with several immunocompetent cells.
Conclusion: CXCR4 may be a useful biomarker of prognosis and an indicator of the TME in GC.
Introduction
Gastric cancer (GC) is one of the most common malignant tumors of the alimentary system, with growing incidences worldwide. Globally, GC is the sixth most frequently diagnosed cancer and the third leading cause of cancer-related deaths (Sitarz et al., 2018). More than 1,000,000 new cases are diagnosed, and approximately 783,000 deaths occur annually (Bray et al., 2018; Global Burden of Disease Cancer Collaboration, Fitzmaurice et al., 2019). Although surgical resection remains the primary curative treatment for GC, the incidence of postoperative tumor recurrence is high. In particular, the 5-year survival rate of patients with stage II, III, and IV GC is approximately 31, 13, and 3%, respectively (Akhondi-Meybodi et al., 2017). Despite remarkable progress in treatment modalities in recent years, the mortality rate of patients with GC remains high (Rawla and Barsouk, 2019). Therefore, it is necessary and urgent to develop novel strategies for the early diagnosis and prognostic prediction to reduce the high mortality and recurrence rates of patients with GC.
Previous investigations have demonstrated that the characteristics of the tumor microenvironment (TME) are closely associated with the progression and prognosis of GC (Wang et al., 2019). The significance of the TME in cancer initiation and progression has drawn increasing attention in recent years. Research has shown that the TME is an active promoter of cancer progression, as opposed to its previous designation as a silent bystander during cancer (Bussard et al., 2016). Emerging evidence has indicated that the TME, which is mainly composed of the extracellular matrix, stromal cells, blood vessels, and lymphatic networks, plays a key role in tumor development and metastasis (Hanahan and Coussens, 2012; Junttila and de Sauvage, 2013; Quail and Joyce, 2013). The type and proportion of stromal cells are related to the physiological state of the TME (Alkasalias et al., 2018). Moreover, tumor-infiltrating immune cells (TIICs), such as CD8+ T cells, regulatory T cells (Tregs), and tumor-associated macrophages (TAMs), positively affect the clinical outcome of patients with various malignancies, including melanoma, lung cancer, breast cancer, and GC (Adams et al., 2013; Massi et al., 2015; Bremnes et al., 2016; Jiang et al., 2017). The dynamic interplay between stromal cells and immune cells in the TME involves several cellular events and physiological processes (Lee et al., 2014). Further investigations of various components and pathways of GC in the TME may facilitate targeted therapy.
Recently, ESTIMATE (Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data), a novel algorithm, has been developed to calculate stromal and immune scores, which are used to assess the extent of stromal and immune cells infiltrating into tumor tissues. The ESTIMATE algorithm helps present a better picture of the numbers of stromal and immune cells in the TME (Yoshihara et al., 2013). Thus, based on the scores calculated by ESTIMATE, the clinical outcomes of patients with GC may be predicted (Liu et al., 2018).
In this article, we collected the gene expression profiles of patients with GC from The Cancer Genome Atlas (TCGA) and used the ESTIMATE algorithm to calculate immune and stromal scores of the TME in GC. Moreover, we investigated the correlation between the risk scores obtained from differentially expressed genes (DEGs) and the clinicopathological characteristics of patients with GC. Furthermore, we constructed a protein–protein interaction (PPI) network and conducted a functional enrichment analysis of the identified DEGs to explore their potential correlations with TIICs.
Materials and Methods
Microarray Data Collection and Processing
From the TCGA dataset, transcriptome and relevant clinical data of 373 patients with GC (343 tumor samples vs. 30 normal samples) were collected, and the ESTIMATE algorithm was employed to evaluate the composition of the TME. The results were represented as three scores, namely immune score, stromal score, and estimate score, corresponding to the proportion of immune cells, stromal cells, and both immune and stromal cells, respectively. All genes of tumor samples were ranked by their expression levels, and DEGs were screened out using the “limma” package in R. DEGs were identified based on the following criteria: an absolute value of log2 fold change (| log2FC|) > 1 and false discovery rate (FDR) < 0.05. Furthermore, the “VennDiagram” package was used to screen for genes with similar expression levels in both stromal and immune cells. The “pheatmap” package was used to produce heatmaps of TME-related DEGs.
Functional Enrichment Analysis of TME-Related DEGs
tumor microenvironment-related DEGs were performed Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) analysis, which revealed the function of DEGs in the biologic process, molecular function, and showing the pathway enrichment result. The “ggplot2,” “enrichplot,” and “clusterProfiler” packages in R were used to perform GO and KEGG analyses. Statistical significance was set at P < 0.05 and q < 0.05.
Correlation Between Scores With Clinicopathological Characteristics and Survival
The clinicopathological characteristics of each sample were evaluated by the Wilcoxon rank-sum and Kruskal–Wallis rank-sum test, clarifying the correlation between scores and the clinical stage. Samples were divided into high- and low-score groups by compared to media value and executed survival analysis. R packages “survival” and “survminer” were applied and P < 0.05 was identified as significant difference.
PPI Network and Cox Regression Analysis
Next, to explore the relationship among DEGs, the STRING platform1 was used to establish a PPI network, and nodes were employed to reconstruct the network with the confidence of interactive relationship greater than 0.95. The Cytoscape software was used to identify the top 30 hub genes. Univariate Cox regression was performed using the “survival” package in R to select DEGs associated with the prognosis of GC. The top 50 genes ranked according to log-rank test P-values in univariate Cox analysis are shown in the plot. Finally, based on the results of the intersection analysis of the PPI network and Cox regression analysis, only the CXC motif chemokine receptor 4 (CXCR4) gene was found to meet all the above-mentioned metrics.
Correlation Between CXCR4 Expression and Clinicopathological Characteristics
Kaplan–Meier survival analysis was performed to illustrate the differences in the overall survival (OS) between the GC groups with low and high expression. Next, correlation analysis was performed between clinical characteristics and CXCR4 expression levels, which were contrasted by univariate analysis. Statistical analysis was performed using SPSS 22.0, and statistical significance was set at P < 0.05.
Further Analysis of the Relationship Between CXCR4 and Tumor Immunoreaction
To explore the role of CXCR4 in the TME of GC, Gene Set Enrichment Analysis (GSEA) was performed to verify the results of KEGG pathway enrichment analysis using GSEA version 3.0 (Broad Institute, Cambridge, MA, United States). Differences were considered significant if NOM P-value < 0.05 and FDR < 0.25. In addition, to determine the relative abundance of TIICs in GC samples, the extent of infiltration was estimated using the CIBERSORT algorithm. Samples with P < 0.05 were identified to have significantly different immune cell infiltration between the two groups. Furthermore, correlation analyses between the expression of CXCR4 and immune cell infiltration in the TME were performed. Additionally, the correlation of the expression levels of CXCR4 with those of immune checkpoint molecules in GC was identified by the TISIDB web portal2. Statistical significance was set at P < 0.05.
Results
Identification and Functional Analysis of DEGs
ImmuneScore, StromalScore, and ESTIMATEScore were dissected by Kaplan–Meier survival analysis. The high- and low-score samples were analyzed and compared to determine the differences in gene expression patterns in immune and stromal components. A total of 2143 DEGs were obtained based on the immune score, out of which 1553 genes were upregulated and 580 genes were downregulated (Figures 1A,C,D). Similarly, 2454 DEGs were acquired based on the stromal score, out of which 2152 genes were upregulated and 302 genes were downregulated (Figures 1B,C,D). Furthermore, the Venn plot identified 1051 upregulated genes and 180 downregulated genes in both the immune and stromal components. These 1231 DEGs were identified as TME-related DEGs. The results of GO analysis demonstrated that these DEGs were mostly engaged in immune-related functions, such as the regulation of lymphocyte activation and lymphocyte-mediated immunity (Figures 2A,C). Moreover, the results of KEGG analysis demonstrated the involvement of DEGs in certain immune-related functions, including cytokine-cytokine receptor interaction and chemokine signaling pathway (Figures 2B,D).
Figure 1. Heatmaps and Venn plots of differentially expressed genes (DEGs). (A) A heatmap of immune-related DEGs between the high- and low CXCR4 expression groups. (B) A heatmap of stromal-related DEGs between the high- and low CXCR4 expression groups. (C) A Venn diagram of commonly upregulated DEGs in the stromal and immune components. (D) A Venn diagram of commonly downregulated DEGs in the stromal and immune components.
Figure 2. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of differentially expressed genes (DEGs) (A–D).
Correlation Between Scores and Clinicopathological Features of Patients With GC
To ascertain the correlation between stromal and immune scores with the clinicopathological characteristics of patients with GC, commensurable clinical features of patients with GC acquired from TCGA were analyzed. The immune score was found to be significantly and positively correlated with the T stage of the tumor (P = 0.00086) (Figure 3D) and the cancer grade (P = 0.016) (Figure 3M). Moreover, the stromal scores were significantly correlated with the T stage of the tumor (P < 0.001) (Figure 3E) and the cancer stage (P = 0.02) (Figure 3B). Furthermore, the estimate score was closely associated with the T and N stages of the tumor (P = 0.024, P < 0.001, and P = 0.036, respectively) (Figures 3C,F,L). However, all three scores were not significantly associated with the M stage of the tumor (P = 0.49, P = 0.61, and P = 0.04, respectively) (Figures 3G–I) while immune score was not statistically significantly correlated with N stage (P = 0.068) (Figure 3J) or the stage of tumor (P = 0.13) (Figure 3A). Neither stromal score or estimate score were associated with gender (P = 0.49, P = 0.37, respectively) (Figures 3N,O). Besides, there was no significant correlation between stromal score and N stage (P = 0.067) (Figure 3K). The correlation of immune, stromal, and estimate scores with patient survival was analyzed by the Kaplan–Meier survival method, and the correlation of each score with the survival rate was assessed. As shown in Figure 4C, the amount of the stromal constituent was negatively correlated with the OS of patients with GC (P = 0.005). However, the immune and estimate scores had no significant correlation with the OS (Figure 4A, P = 0.233; Figure 4B, P = 0.476). These results suggest that the stromal components in the TME are significantly associated with the prognosis of patients with GC.
Figure 3. Correlation of clinicopathological characteristics with immune, stromal, and estimate scores (A–O).
Figure 4. Correlation of each estimate, immune, and stromal score with the survival outcome of patients with gastric cancer (A–C).
Intersection Analysis of Univariate Cox Regression and PPI Network
Next, we thoroughly investigated the interactions among DEGs. We constructed a PPI network using the STRING database in Cytoscape. The correlation between each DEG and the top 30 genes ranked by the nodes was displayed in Figures 5A,B. The results of univariate Cox regression analysis showed that 50 genes were associated with the prognosis of GC (Figure 5C). The results of intersection analysis between the 30 hub genes and 50 prognostic DEGs revealed CXCR4 as the only overlapping gene (Figure 5D).
Figure 5. Protein–protein interaction (PPI) network and univariate Cox regression analysis. (A) Construction of the PPI network of 1231 differentially expressed genes (DEGs). (B) The top 30 genes ranked by the number of adjacent nodes of the PPI network. (C) Results of univariate Cox regression analysis with selected DEGs, with the top 50 genes displayed in the forest plot. (D) A Venn diagram showing CXCR4 as the most commonly expressed DEG in gastric cancer by combination analysis of the top 30 genes in the PPI network and the top 50 prognostic genes from the results of Cox regression analysis.
Association Between CXCR4 Expression With Clinicopathological Factors and Disease Progression
In the present study, GC samples were divided into high CXCR4 and low CXCR4 groups using the median expression value of CXCR4 as the threshold value. Survival analysis showed that GC patients with low CXCR4 expression had a longer survival time than those with high CXCR4 expression (Figure 6C, P = 0.010). Results of the Wilcoxon rank-sum test demonstrated that the level of CXCR4 in tumor tissues was significantly higher than that in healthy tissues in both paired or unpaired samples (Figure 6B, P = 0.045; Figure 6A, P = 0.008). In addition, the expression of CXCR4 was significantly different between different age groups (Figure 6I, P = 0.042). Additionally, the level of CXCR4 was inextricably linked to the T stage (Figure 6D, P = 2.4e-06) and the stage of cancer (Figures 6G, P = 0.0036). However, there were no significant differences between the expression level of CXCR4 and the M stage (Figures 6F, P = 0.089), N stage (Figures 6E, P = 0.089), grade (Figure 6H, P = 0.17), or sex (Figures 6J, P = 0.15).
Figure 6. The relationship between CXCR4 expression with survival rate, clinical and pathological characteristics of patients with gastric cancer. (A) unpaired samples, (B) paired samples, (C) survival rate, (D) T classification, (E) N classification, (F) M classification, (G) stage, (H) grade, (I) age, and (J) sex.
Role of CXCR4 in the TME of Gastric Cancer
As shown in Figure 7, CXCR4 was mainly engaged in immune-related activities, such as an intestinal immune network for IgA production, JAK-STAT signaling pathway, natural killer cell-mediated cytotoxicity, and Toll-like receptor signaling pathway. Using the CIBERSORT algorithm, we identified the infiltrating profiles of 22 different types of immune cells in tumor tissues (Figures 8A,B). A total of six types of TIICs were found to be strongly correlated with CXCR4 expression in the TME of GC cells (Figure 9). Immune cells such as memory B cells, resting dendritic cells, CD8+ T cells, monocytes, and Tregs were positively related with the expression of CXCR4, while decreased activation of mast cells was negatively correlated with CXCR4 expression. Furthermore, we characterized the interactions between CXCR4 with 22 immune control genes. As shown in Figure 10, the expression level of CXCR4 was positively correlated with that of 20 immune checkpoint molecules, including CD274, CTLA-4, and LAG3, among others. Therefore, these results indicate that CXCR4 plays an important role in the immune evasion of GC cells.
Figure 7. Gene set enrichment analysis of CXCR4.
Figure 8. Profiles of tumor-infiltrating immune cells (TIICs) in gastric cancer (GC) samples and their correlation analysis. (A) Barplot showing the proportion of 21 different types of TIICs in GC tumor samples. (B) A heatmap showing the correlation between 21 different types of TIICs.
Figure 9. Relation between the expression of CXCR4 and proportion of tumor-infiltrating immune cells (TIICs). (A) A violin plot showing the differences in the proportions of 21 different types of immune cells in GC tumor samples with low or high CXCR4 expression. (B) A Venn plot displaying the differences and correlation between two types of TIICs associated with CXCR4 expression. (C–H) A scatter plot showing the correlation of the proportions of TIICs with CXCR4 expression.
Figure 10. Correlation analysis of the expression of CXCR4 expression with that of immune checkpoint genes. (A,B) Correlation analysis of the expression levels of CXCR4 with those of 22 common immune checkpoint genes in gastric cancer. (C–X) The relationship between the expression of CXCR4 with that of IDO1, PVRL2, ADORA2A, BTLA, CD96, CD160, CD244, PDCD1, CD274, CTLA4, HAVCR2, LAG3, LGALS9, VTCN1, CSF1R, IL10RB, KDR, IL10, TGFBR1, TIGIT, PDCD1LG2, and TGFB1.
Discussion
Great advances in whole-genome sequencing have facilitated the development of molecular classification systems and treatment strategies for cancer. In the present study, we identified CXCR4 as a TME-related gene associated with survival and TMN-stage classification in GC samples gathered from TCGA database.
CXC motif chemokine receptor 4, the receptor for chemokine CXCL12/SDF-1 and a member of the G protein-coupled receptor superfamily, is overexpressed in various types of solid cancers, including non-small-cell lung cancer (NSCLC), breast cancer, colorectal cancer, and GC. Previous studies have demonstrated that CXCR4 accelerates the metastasis, invasion, growth, and therapeutic resistance of cancer (Zhao et al., 2010; Otsuka et al., 2011; Ying et al., 2012; Mukherjee and Zhao, 2013; Xu et al., 2018). Moreover, CXCR4 affects the migration of GC cells via the ERK/Akt signaling pathway (Cheng et al., 2017). Cheng et al. (2020) found that the positive crosstalk between CXCR4 and EGFR promotes GC metastasis via the NF-kB pathway. Another study revealed that CXCR4 activates the NF-kB pathway and upregulates the expression of serine proteinase inhibitor clade B member 3, thereby facilitating the migration and invasion of GC cells (Gong et al., 2020). CXCR4 also plays a critical role in tumor angiogenesis in GC by activating the JAK2/STAT3 (Zhang et al., 2017). Furthermore, it promotes the proliferation and invasion processes via the Wnt/β-catenin pathway (Lin et al., 2017). Thus, inhibition of CXCR4 can disrupt multiple processes that facilitate the growth and spread of GC tumors. Therefore, given its multiple functions, CXCR4 may prove to be a promising target for immunotherapy.
Tumor progression is determined not only by cancer cells but also by the TME, which is the internal environment of malignant tumor progression. The TME can reduce the resistance of cancer cells to chemotherapy and immunotherapy (Russi et al., 2019). Epithelial-mesenchymal transition (EMT), a major modulator of tumor metastasis, may be involved in the interaction between tumor cells and the TME (Zheng et al., 2015; Suarez-Carmona et al., 2017; Hsieh et al., 2018). Compelling evidence indicates that CXCR4 regulates tumor EMT together with the c-MET signaling pathway (Quail and Joyce, 2013). The stromal component, another key component of the TME, primarily consists of cancer-associated fibroblasts (CAFs), which drive the growth, metastasis, and malignancy of cancer cells (Cheng et al., 2018). One study showed that CXCL12 secreted from CAFs promotes GC cell invasion by enhancing the clustering of integrin β1 in GC cells (Daisuke et al., 2015). In addition, CXCL12 mainly stems from the stromal part, which directly stimulates the proliferation and migration of CXCR4-expressing cells. Thus, the specific interaction between stromal cells and tumor cells might be one of the causes of drug resistance. Fortunately, AMD300, a CXCR4 inhibitor approved by the United States Food and Drug Administration (FDA), has already been shown to disrupt tumor-stromal interactions, sensitizing cancer cells to docetaxel-based chemotherapy in prostate cancer (Domanska et al., 2012). Currently, AMD300 is the most frequently used drug targeting the CXCL12-CXCR4 axis in clinical trials for solid gastrointestinal tumors, and therefore, CXCR4 might be a promising candidate target for GC immunotherapy.
In this study, the CIBERSORT algorithm was used to analyze the proportion of TIICs. The results showed that Tregs, dendritic cells, eosinophils, and CD8+ T cells were significantly positively correlated with CXCR4 expression in patients with GC, which is concurrent with our hypothesis that CXCR4 may be the hub gene of the TME in GC. Previous studies have revealed that CD8+ T cells are associated with poor prognosis in GC (Thompson et al., 2017). Similarly, Tregs are known to play an immunosuppressive role in the TME. Further reports have indicated that Tregs derived from patients with cancer usually express diverse chemokine receptors, which contributes to their migration into tumors in response to the signals stemming from the TME (Yan et al., 2011). However, some studies argued that Treg infiltration predicts favorable outcomes for patients with GC. For example, Li et al. (2016) reported that lower FOXP3+ and GARP+ Treg levels after neoadjuvant chemotherapy are associated with good outcomes in progressive GC. Thus, Treg infiltration may play a subtle yet vital role in GC progression. Nevertheless, more subset and related molecule regulation mechanisms of Tregs should to investigated thoroughly to better evaluate the prognosis of patients with GC (Liu et al., 2019).
Moreover, we have attempted to characterized the correlation between the expression levels of CXCR4 with that of immune checkpoint molecules. The results revealed moderate positive correlations between the expression of CXCR4 with that of PD-L1 (CD274) or CTLA4 in GC, which, in some ways, may be exploited for improving the efficacy of immunotherapy. Previous evidence indicated that suppression of CXCR4 promotes anti-PD-1/PD-L1 efficacy by reshaping the TME in hepatocellular carcinoma (Chen et al., 2015). A similar conclusion was observed in osteosarcoma, where Feig et al. (2013) showed that the SDF-1/CXCR4 axis facilitates the accumulation of myeloid-derived suppressor cells in the TME to abate the response to anti-PD-1 therapy. Another study discovered that in patients with colorectal cancer liver metastases, the expression of CXCR4 CXCR7, TLR2/TLR4, and PD-1/PD-L1 was significantly increased in metastatic liver tissues compared to unaffected liver tissues (D’Alterio et al., 2016). Interestingly, our study found that CXCR4 was closely related to the Toll-like receptor signaling pathway (Figure 10). However, more experiments should be conducted to elaborate the interaction of CXCR4, TLRs, and PD-1/PD-L1 in GC. Furthermore, Pep R, a novel CXCR4 antagonist, has been found to enhance the efficiency of anti-PD-1 in various models (D’Alterio et al., 2019). X4-136, another CXCR4 inhibitor, could serve as a monotherapy or combined with immune checkpoint inhibitors in the treatment of melanoma and renal cell carcinoma (Saxena et al., 2020). In addition, a few studies have demonstrated that the combined blockade of CXCL12-CXCR4 and PD-1- PD-L1 pathways could provide survival benefits by regulating the TME of various solid tumors (Feig et al., 2013; Wu et al., 2019; Zeng et al., 2019; Seo et al., 2019) shedding light on CXCR4/PD-1-targeting combination therapy in GC. Given these advances, although the correlation between the expression of CXCR4 and PD-1/PD-L1 was moderate, we theorized that it has far-reaching clinical implications and relevance, which needs further experimental verification.
This study has some limitations that need to be acknowledged. First, as the clinical data were mainly acquired from TCGA database, result biases were unavoidable. Additionally, we did not conduct experimental research to inspect the function of CXCR4 in GC. The combined application of CXCR4 blocker and PD-1 inhibitor may prolong the survival time of GC. However, more evidence is needed to prove the mechanism of combined immunotherapy.
In summary, we employed the ESTIMATE algorithm to identify genes that were associated with the TME in GC samples gathered from TCGA database. Consequently, CXCR4 was discovered as a promising prognostic target for patients with GC. Nevertheless, more experimental research is warranted to explore the underlying molecular mechanisms and potential clinical value of CXCR4 for the early diagnosis of tumor micrometastasis.
Conclusion
Overall, the ESTIMATE algorithm was used to calculate the immune, stromal, and estimate scores of GC samples acquired from TCGA. Stromal and immune cells that infiltrated into the TME were closely related to tumor growth. We identified a few TME-related DEGs, out of which CXCR4, was significantly associated with the regulation of the immune-active status in the TME. Therefore, CXCR4 might be a latent biomarker in GC, which determines the efficacy of cancer immunotherapy. The results of the present study may provide new insights into the development of effective therapeutic strategies targeted against GC.
Data Availability Statement
Publicly available datasets were analyzed in this study. This data can be found here: the datasets analyzed for this study can be found in the TCGA (https://cancergenome.nih.gov/) and ImmPort (https://www.immport.org).
Author Contributions
TX and ZF conceived and devised the study and gave aid in writing the manuscript. YG and WG performed the data analyses and contributed to the writing of the manuscript. ZC and RX reviewed the original manuscript. All authors reviewed the manuscript.
Funding
This study was supported by grant from the Natural Science Foundation of Zhejiang Province (No. LY21H280012).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Footnotes
References
Adams, S., Gray, R. J., Demaria, S., Goldstein, L. J., and Badve, S. S. (2013). Abstract S1-07: prognostic value of tumor-infiltrating lymphocytes (TILs) in two phase III randomized adjuvant breast cancer trials: ECOG 2197 and ECOG 1199. Cancer Res. 73, S1–S07.
Akhondi-Meybodi, M., Ghane, M., Akhondi-Meybodi, S., and Dashti, G. (2017). Five-year survival rate for gastric cancer in Yazd Province, Central Iran, from 2001 to 2008. Middle East J. Dig. Dis. 9, 39–48. doi: 10.15171/mejdd.2016.50
Alkasalias, T., Moyano-Galceran, L., Arsenian-Henriksson, M., and Lehti, K. (2018). Fibroblasts in the tumor microenvironment: shield or spear. Int. J. Mol. Sci. 19:1532. doi: 10.3390/ijms19051532
Bray, F., Ferlay, J., Soerjomataram, I., Siegel, R. L., Torre, L. A., and Jemal, A. (2018). Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 68, 394–424. doi: 10.3322/caac.21492
Bremnes, R. M., Busund, L. T., Kilvær, T. L., Andersen, S., Richardsen, E., Paulsen, E. E., et al. (2016). The role of tumor-infiltrating lymphocytes in development, progression, and prognosis of non-small cell lung cancer. J. Thorac. Oncol. 11, 789–800. doi: 10.1016/j.jtho.2016.01.015
Bussard, K. M., Mutkus, L., Stumpf, K., Gomez-Manzano, C., and Marini, F. C. (2016). Tumor-associated stromal cells as key contributors to the tumor microenvironment. Breast Cancer Res. 18:84. doi: 10.1186/s13058-016-0740-2
Chen, Y., Ramjiawan, R. R., Reiberger, T., Ng, M. R., Hato, T., Huang, Y., et al. (2015). CXCR4 inhibition in tumor microenvironment facilitates anti-programmed death receptor-1 immunotherapy in sorafenib-treated hepatocellular carcinoma in mice. Hepatology 61, 1591–1602. doi: 10.1002/hep.27665
Cheng, Y., Che, X., Zhang, S., Guo, T., He, X., Liu, Y., et al. (2020). Positive cross-talk between CXC chemokine receptor 4 (CXCR4) and epidermal growth factor receptor (EGFR) promotes gastric cancer metastasis via the nuclear factor kappa B (NF-kB)-dependent pathway. Med. Sci. Monit. 26:e925019. doi: 10.12659/MSM.925019
Cheng, Y., Qu, J., Che, X., Xu, L., Song, N., Ma, Y., et al. (2017). CXCL12/SDF-1α induces migration via SRC-mediated CXCR4-EGFR cross-talk in gastric cancer cells. Oncol. Lett. 14, 2103–2110. doi: 10.3892/ol.2017.6389
Cheng, Y., Song, Y., Qu, J., Che, X., Song, N., Fan, Y., et al. (2018). The chemokine receptor CXCR4 and c-MET cooperatively promote epithelial-mesenchymal transition in gastric cancer cells. Transl. Oncol. 11, 487–497. doi: 10.1016/j.tranon.2018.02.002
Daisuke, I., Ishimoto, T., Miyake, K., Sugihara, H., Eto, K., Sawayama, H., et al. (2015). CXCL12/CXCR4 activation by cancer-associated fibroblasts promotes integrin β1 clustering and invasiveness in gastric cancer. Int. J. Cancer 138, 1207–19.
D’Alterio, C., Buoncervello, M., Ieranó, C., Napolitano, M., Portella, L., Rea, G., et al. (2019). Targeting CXCR4 potentiates anti-PD-1 efficacy modifying the tumor microenvironment and inhibiting neoplastic PD-1. J. Exp. Clin. Cancer Res. 38:432. doi: 10.1186/s13046-019-1420-8
D’Alterio, C., Nasti, G., Polimeno, M., Ottaiano, A., Conson, M., Circelli, L., et al. (2016). CXCR4-CXCL12-CXCR7, TLR2-TLR4, and PD-1/PD-L1 in colorectal cancer liver metastases from neoadjuvant-treated patients. Oncoimmunology 5:e1254313. doi: 10.1080/2162402X.2016.1254313
Domanska, U. M., Timmer-Bosscha, H., Nagengast, W. B., Oude Munnink, T. H., Kruizinga, R. C., Ananias, H. J., et al. (2012). CXCR4 inhibition with AMD3100 sensitizes prostate cancer to docetaxel chemotherapy. Neoplasia 14, 709–718. doi: 10.1593/neo.12324
Feig, C., Jones, J. O., Kraman, M., Wells, R. J., Deonarine, A., Chan, D. S., et al. (2013). Targeting CXCL12 from FAP-expressing carcinoma-associated fibroblasts synergizes with anti-PD-L1 immunotherapy in pancreatic cancer. Proc. Natl. Acad. Sci. U.S.A. 110, 20212–20217. doi: 10.1073/pnas.1320318110
Global Burden of Disease Cancer Collaboration, Fitzmaurice, C., Abate, D., Abbasi, N., Abbastabar, H., Abd-Allah, F., et al. (2019). Global, regional, and national cancer incidence, mortality, years of life lost, years lived with disability, and disability-adjusted life-years for 29 cancer groups, 1990 to 2017: a systematic analysis for the global burden of disease study. JAMA Oncol. 5, 1749–1768. doi: 10.1001/jamaoncol.2019.2996
Gong, J., Song, Y., Xu, L., Che, X., Hou, K., Guo, T., et al. (2020). Upregulation of serine proteinase inhibitor clade B member 3 (SERPINB3) expression by stromal cell-derived factor (SDF-1)/CXCR4/nuclear factor kappa B (NF-κB) promotes migration and invasion of gastric cancer cells. Med. Sci. Monit. 26:e927411. doi: 10.12659/MSM.927411
Hanahan, D., and Coussens, L. M. (2012). Accessories to the crime: functions of cells recruited to the tumor microenvironment. Cancer Cell 21, 309–322. doi: 10.1016/j.ccr.2012.02.022
Hsieh, C. H., Tai, S. K., and Yang, M. H. (2018). Snail-overexpressing cancer cells promote M2-like polarization of tumor-associated macrophages by delivering MiR-21-abundant exosomes. Neoplasia 20, 775–788. doi: 10.1016/j.neo.2018.06.004
Jiang, W., Liu, K., Guo, Q., Cheng, J., Shen, L., Cao, Y., et al. (2017). Tumor-infiltrating immune cells and prognosis in gastric cancer: a systematic review and meta-analysis. Oncotarget 8, 62312–62329. doi: 10.18632/oncotarget.17602
Junttila, M. R., and de Sauvage, F. J. (2013). Influence of tumour micro-environment heterogeneity on therapeutic response. Nature 501, 346–354. doi: 10.1038/nature12626
Lee, K., Hwang, H., and Nam, K. T. (2014). Immune response and the tumor microenvironment: how they communicate to regulate gastric cancer. Gut Liver 8, 131–139. doi: 10.5009/gnl.2014.8.2.131
Li, K., Chen, F., and Xie, H. (2016). Decreased FOXP3+ and GARP+ Tregs to neoadjuvant chemotherapy associated with favorable prognosis in advanced gastric cancer. Onco Targets Ther. 9, 3525–3533. doi: 10.2147/OTT.S101884
Lin, X. L., Xu, Q., Tang, L., Sun, L., Han, T., Wang, L. W., et al. (2017). Regorafenib inhibited gastric cancer cells growth and invasion via CXCR4 activated Wnt pathway. PLoS One 12:e0177335. doi: 10.1371/journal.pone.0177335
Liu, W., Ye, H., Liu, Y., Xu, C., Zhong, Y., Tian, T., et al. (2018). Transcriptome-derived stromal and immune scores infer clinical outcomes of patients with cancer. Oncol. Lett. 15, 4351–4357.
Liu, X., Zhang, Z., and Zhao, G. (2019). Recent advances in the study of regulatory T cells in gastric cancer. Int. Immunopharmacol. 73, 560–567. doi: 10.1016/j.intimp.2019.05.009
Massi, D., Brusa, D., Merelli, B., Falcone, C., Xue, G., Carobbio, A., et al. (2015). The status of PD-L1 and tumor-infiltrating immune cells predict resistance and poor prognosis in BRAFi-treated melanoma patients harboring mutant BRAFV600. Ann. Oncol. 26, 1980–1987. doi: 10.1093/annonc/mdv255
Mukherjee, D., and Zhao, J. (2013). The Role of chemokine receptor CXCR4 in breast cancer metastasis. Am. J. Cancer Res. 3, 46–57.
Otsuka, S., Klimowicz, A. C., Kopciuk, K., Petrillo, S. K., Konno, M., Hao, D., et al. (2011). CXCR4 overexpression is associated with poor outcome in females diagnosed with stage IV non-small cell lung cancer. J. Thorac. Oncol. 6, 1169–1178. doi: 10.1097/JTO.0b013e3182199a99
Quail, D. F., and Joyce, J. A. (2013). Microenvironmental regulation of tumor progression and metastasis. Nat. Med. 19, 1423–1437. doi: 10.1038/nm.3394
Rawla, P., and Barsouk, A. (2019). Epidemiology of gastric cancer: global trends, risk factors and prevention. Prz. Gastroenterol. 14, 26–38. doi: 10.5114/pg.2018.80001
Russi, S., Verma, H. K., Laurino, S., Mazzone, P., Storto, G., Nardelli, A., et al. (2019). Adapting and surviving: intra and extra-cellular remodeling in drug-resistant gastric cancer cells. Int. J. Mol. Sci. 20:3736. doi: 10.3390/ijms20153736
Saxena, R., Wang, Y., and Mier, J. W. (2020). CXCR4 inhibition modulates the tumor microenvironment and retards the growth of B16-OVA melanoma and Renca tumors. Melanoma Res. 30, 14–25. doi: 10.1097/CMR.0000000000000639
Seo, Y. D., Jiang, X., Sullivan, K. M., Jalikis, F. G., Smythe, K. S., Abbasi, A., et al. (2019). Mobilization of CD8 + T cells via CXCR4 blockade facilitates PD-1 checkpoint therapy in human pancreatic cancer. Clin. Cancer Res. 25, 3934–3945. doi: 10.1158/1078-0432.CCR-19-0081
Sitarz, R., Skierucha, M., Mielko, J., Offerhaus, G., Maciejewski, R., and Polkowski, W. P. (2018). Gastric cancer: epidemiology, prevention, classification, and treatment. Cancer Manag. Res. 10, 239–248. doi: 10.2147/CMAR.S149619
Suarez-Carmona, M., Lesage, J., Cataldo, D., and Gilles, C. (2017). EMT and inflammation: inseparable actors of cancer progression. Mol. Oncol. 11, 805–823. doi: 10.1002/1878-0261.12095
Thompson, E. D., Zahurak, M., Murphy, A., Cornish, T., Cuka, N., Abdelfatah, E., et al. (2017). Patterns of PD-L1 expression and CD8 T cell infiltration in gastric adenocarcinomas and associated immune stroma. Gut 66, 794–801. doi: 10.1136/gutjnl-2015-310839
Wang, H., Wu, X., and Chen, Y. (2019). Stromal-immune score-based gene signature: a prognosis stratification tool in gastric cancer. Front. Oncol. 9:1212. doi: 10.3389/fonc.2019.01212
Wu, A., Maxwell, R., Xia, Y., Cardarelli, P., Oyasu, M., Belcaid, Z., et al. (2019). Combination anti-CXCR4 and anti-PD-1 immunotherapy provides survival benefit in glioblastoma through immune cell modulation of tumor microenvironment. J. Neurooncol. 143, 241–249. doi: 10.1007/s11060-019-03172-5
Xu, C., Zheng, L., Li, D., Chen, G., Gu, J., Chen, J., et al. (2018). CXCR4 overexpression is correlated with poor prognosis in colorectal cancer. Life Sci. 208, 333–340. doi: 10.1016/j.lfs.2018.04.050
Yan, M., Jene, N., Byrne, D., Millar, E. K., O’Toole, S. A., McNeil, C. M., et al. (2011). Recruitment of regulatory T cells is correlated with hypoxia-induced CXCR4 expression, and is associated with poor prognosis in basal-like breast cancers. Breast Cancer Res. 13:R47. doi: 10.1186/bcr2869
Ying, J., Xu, Q., Zhang, G., Liu, B., and Zhu, L. (2012). The expression of CXCL12 and CXCR4 in gastric cancer and their correlation to lymph node metastasis. Med. Oncol. 29, 1716–1722. doi: 10.1007/s12032-011-9990-0
Yoshihara, K., Shahmoradgoli, M., Martínez, E., Vegesna, R., Kim, H., Torres-Garcia, W., et al. (2013). Inferring tumour purity and stromal and immune cell admixture from expression data. Nat. Commun. 4:2612. doi: 10.1038/ncomms3612
Zeng, Y., Li, B., Liang, Y., Reeves, P. M., Qu, X., Ran, C., et al. (2019). Dual blockade of CXCL12-CXCR4 and PD-1-PD-L1 pathways prolongs survival of ovarian tumor-bearing mice by prevention of immunosuppression in the tumor microenvironment. FASEB J. 33, 6596–6608. doi: 10.1096/fj.201802067RR
Zhang, Q., Xu, F., Shi, Y., Chen, Y. W., Wang, H. P., Yu, X., et al. (2017). C-X-C motif chemokine receptor 4 promotes tumor angiogenesis in gastric cancer via activation of JAK2/STAT3. Cell Biol. Int. 41, 854–862. doi: 10.1002/cbin.10794
Zhao, C., Lu, X., Bu, X., Zhang, N., and Wang, W. (2010). Involvement of tumor necrosis factor-alpha in the upregulation of CXCR4 expression in gastric cancer induced by Helicobacter pylori. BMC Cancer 10:419. doi: 10.1186/1471-2407-10-419
Keywords: CXCR4, tumor microenvironment, gastric cancer, ESTIMATE algorithm, CIBERSORT algorithm
Citation: Gu Y, Gu W, Xie R, Chen Z, Xu T and Fei Z (2021) Role of CXCR4 as a Prognostic Biomarker Associated With the Tumor Immune Microenvironment in Gastric Cancer. Front. Cell Dev. Biol. 9:654504. doi: 10.3389/fcell.2021.654504
Received: 16 January 2021; Accepted: 19 August 2021;
Published: 08 September 2021.
Edited by:
Liu Yang, Zhejiang Provincial People’s Hospital, ChinaReviewed by:
Raphael Carmo Valente, Federal University of Rio de Janeiro, BrazilXiaodong Liang, Zhejiang Provincial People’s Hospital, China
Copyright © 2021 Gu, Gu, Xie, Chen, Xu and Fei. 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: Tongpeng Xu, tongpeng_xu_njmu@163.com; Zhenghua Fei, 940077455@qq.com