Skip to main content

Comprehensive analysis of CXCR family members in lung adenocarcinoma with prognostic values

Abstract

Background

The expression profiles and molecular mechanisms of CXC chemokine receptors (CXCRs) in Lung adenocarcinoma (LUAD) have been extensively explored. However, the comprehensive prognostic values of CXCR members in LUAD have not yet been clearly identified.

Methods

Multiple available datasets, including Oncomine datasets, the cancer genome atlas (TCGA), HPA platform, GeneMANIA platform, DAVID platform and the tumor immune estimation resource (TIMER) were used to detect the expression of CXCRs in LUAD, as well as elucidate the significance and value of novel CXCRs-associated genes and signaling pathways in LUAD.

Results

The mRNA and/or protein expression of CXCR1, CXCR2, CXCR3, CXCR4, CXCR5 and CXCR6 displayed predominantly decreased in LUAD tissues as compared to normal tissues. On the contrary, compared with the normal tissues, the expression of CXCR7 was significantly increased in LUAD tissues. Subsequently, we constructed a network including CXCR family members and their 20 related genes, and the related GO functions assay showed that CXCRs connected with these genes participated in the process of LUAD through several signal pathways including Chemokine signaling pathway, Cytokine-cytokine receptor interaction and Neuroactive ligand-receptor interaction. TCGA and Timer platform revealed that the mRNA expression of CXCR family members was significantly related to individual cancer stages, cancer subtypes, patient’s gender and the immune infiltration level. Finally, survival analysis showed that low mRNA expression levels of CXCR2 (HR = 0.661, and Log-rank P = 1.90e−02), CXCR3 (HR = 0.674, and Log-rank P = 1.00e−02), CXCR4 (HR = 0.65, and Log-rank P = 5.01e−03), CXCR5 (HR = 0.608, and Log-rank P = 4.80e−03) and CXCR6 (HR = 0.622, and Log-rank P = 1.85e−03) were significantly associated with shorter overall survival (OS), whereas high CXCR7 mRNA expression (HR = 1.604, and Log-rank P = 4.27e−03) was extremely related with shorter OS in patients.

Conclusion

Our findings from public databases provided a unique insight into expression characteristics and prognostic values of CXCR members in LUAD, which would be benefit for the understanding of pathogenesis, diagnosis, prognosis prediction and targeted treatment in LUAD.

Peer Review reports

Background

Lung cancer is the most common type of human malignancy and the leading cause of cancer death [1]. Non-small cell lung cancer (NSCLC) is a common type of lung cancer, accounting for approximately 90% of all malignancies of the lung [2]. The two most common subtypes of non-small cell lung cancer are adenocarcinoma (LUAD) and squamous cell carcinoma (LUSC), which account for 50% and 40% of the total, respectively [3, 4]. Lung cancer consists of subpopulations of cells or clones that have distinct molecular characteristics of tumors and therefore genetically characterize the malignant behavior of most malignant tumors. Clonal mutations in patients with LUAD may increase the probability of postoperative recurrence, implying a higher propensity for early development of metastasis in lung adenocarcinoma and increasing cancer heterogeneity. [5] Similarly there are studies on the mechanisms of drug resistance in cancer showing that resistance to some drugs in lung cancer patients is objectively present before treatment. [6, 7]

Chemokines are a class of cytokines associated with cellular secretion and structure; their initial discovery and most important role is to make up the extracellular matrix of tumors and have the ability to directly influence cancer cell proliferation and metastasis. [8,9,10] Among them, CXC chemokines are a 7-transmembrane G protein-coupled receptor protein localized on the cell membrane with two cysteine residues near the N terminus. [11, 12] CXCR binds to its cognate ligand, changes the conformation of the receptor, and then activates the coupled G protein to start the corresponding signaling pathway to act in the cell. Most CXCRs are typical receptors coupled to G proteins, but CXCR7 is an atypical receptor coupled to a β receptor, also known as atypical chemokine receptor 3 (ACKR3). [13] In addition, some CXCRs do not bind only one ligand; for example, CXCR2 can bind to multiple ligands to exert different effects. [14] Following the action of the corresponding ligands that bind CXCRs, the receptors are usually degraded or restored to the plasma membrane by clathrin-mediated endocytosis. [15] CXCRs facilitate communication links between the intracellular and extracellular microenvironments, which in turn affect cellular behaviors such as cellular transport, proliferation and invasion. [16, 17] In addition, tumor angiogenesis and tumor immunity have emerged as one of the research areas related to CXCR in recent years.

Although a series of studies have elucidated the significant prognostic roles of some CXCR members in LUAD, the whole picture of the prognostic values of CXCR members remain inequitably characterized in LUAD. In this study, the clinical significance of CXCR family members in LUAD was analyzed in terms of mRNA expression level, protein expression level, immunity, interaction network, pathway analysis, clinical and prognostic aspects by using multiple databases and platforms based on the study of Xiaojuan Li et al. [18] in order to provide a general overview and expansion of the value of CXCRs in LUAD.

Methods

Differential study of CXCRs at the transcriptional level

Oncomine (https://www.oncomine.org/resource/login.html) is an integrated online gene microarray database and data mining platform that provides peer review, powerful analysis methods, and a robust set of analytical capabilities to calculate gene expression levels. [19] During the analysis of differential mRNA expression in LUAD tissues and their corresponding normal lung gland tissues, we selected data for CXCRs by the following criteria: P-value < 0.05, fold change = 1.5, gene rank = 10%, the Benjamini and Hochberg method was used to correct the P value.

The LUAD data in the Cancer Genome Atlas (TCGA) database (https://portal.gdc.cancer.gov/ repository) were analyzed using the UALCAN (http://ualcan.path.uab.edu/) platform. The database is a more comprehensive information platform containing gene expression databases and corresponding clinical information data. [20] UALCAN is an open access database developed on multiple platforms such as TCGA. [21] In this study, the differences in CXCRs mRNA expression levels expressed in normal samples and lung adenocarcinoma tissues were analyzed by comparing two data sets, UALCAN platform and oncomine.

Differentially study CXCRs at the protein level

HPA (https://www.proteinatlas.org/) is a platform containing Protein expression level data for a wide range of common cancers. [22] Immunohistochemistry (IHC) staining data of all CXCRs were obtained from the HPA database and analyzed. In the HPA database, protein expression scoring was evaluated by taking among staining, intensity and quantity, including not detected (negative, none), low (weak, < 25%), medium (moderate, 25–75%), and high (strong, > 75%).

Construction of protein interaction networks

GeneMANIA (http: //www.genemania.org/) is an interactive, visual online protein interaction prediction tool. [23] Given a list of gene queries, GeneMANIA uses available genomics and proteomics data to search for functionally similar genes to predict interacting genes for the target gene. In order to obtain the genes interacting with CXCRs and to construct a protein interaction network, seven CXCRs were studied as a whole in this study at GeneMANIA.

GO enrichment analysis and kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment analysis

DAVID (https://david.ncifcrf.gov/) is a functional enrichment analysis web tool with continuously updated and effectively reduce data redundancy. [24] In this study the biological processes, cellular components, molecular functions and Kyoto genes and genomic encyclopedia pathways involved in CXCR and the 20 genes associated with it were obtained by enrichment using DAVID, and then the results were visualized by creating enrichment bubble maps using R software.

Correlation analysis of CXCRs and immune infiltration

TIMER (https://cistrome.shinyapps.io/timer/) is an online tool for systematical analyses of immune infiltration of various cancers. [25] Then the correlation between the six immune cells and CXCR was analyzed by this tool. Their required correlations and P-values were automatically calculated and displayed in the graphs. And the Benjamini and Hochberg method was used to correct the P value.

Clinicopathological analysis associated with CXCRs

In addition, downloaded TCGA data were used to analyze the association between mRNA expression of CXCRs in LUAD tissues and their clinicopathological parameters (e.g., TNM stage, cell subtype, and gender of the patient). The results could be got by selecting R (4.02) integrated into the TCGA database.

Survival analysis

The TCGA-LUAD gene expression data and clinical data were read separately using R (v4.0.2) software, and the patients were divided into high and low expression groups using the median expression of the samples as the cut-off point. The "survival" package was applied to analyze the survival rate of the high and low gene expression groups. P < 0.05 was considered statically significant.

Results

MRNA expression of CXCR family members in LUAD patients

This study first used the Oncomine database to compare CXCRs mRNA expression levels from multiple cancers, and then to identify differences in CXCRs expression in LUAD patients. By comparing the expression of CXCRs in multiple cancer species, Oncomine data visualization results showed differences in their mRNA expression levels in LUAD patients. (Fig. 1) The differential expression of CXCRs in different subtypes of lung cancer was then obtained by pooling data from various sources in Oncomine as follows. The mRNA expression levels of CXCR1 (P = 4.71E−7, P = 1.19E−6), CXCR2 (P = 4.20E−16, P = 1.76E−11), CXCR4 (P = 1.63E−15) and CXCR7 (P = 1.31E−12, P = 1.70E−8) were significantly lower in LUAD than in adjacent tissues. (Table 1).

Fig. 1
figure 1

Transcriptional expressions of different CXCR family members in different types of cancers. The data were compared by the t-test and cut-off P-value and fold change were as following: P-value < 0.05, fold change = 1.5, gene rank = 10%. (Red indicates over-expression, blue indicates down-expression)

Table1 CXCR1, CXCR2, CXCR4, CXCR5 and CXCR7 displayed significant regulated expression in lung cancer compared to normal from Oncomine

The results derived from Oncomine were then validated by the UALCAN platform. The results showed that mRNA expression levels of CXCR1 and CXCR2 were significantly lower and CXCR3 was significantly upregulated in LUAD compared to adjacent tissues. (Fig. 2).

Fig. 2
figure 2

MRNA expressions of CXCR family members in patients with LUAD and normal lung tissues. The mRNA expressions of different CXCR family members were significantly regulated in patients with LUAD from the TCGA database (*P < 0.05, **P < 0.01, ***P < 0.001.)

Protein expression of CXCR family members in LUAD patients

We then found that the protein expression of CXCR2, CXCR3 and CXCR7 was unchanged in LUAD compared to normal samples. The LUAD tissues displayed low (3/3) CXCR1 staining and the CXCR1 staining in lung adenocarcinoma cells was usually not detected (11/11). While the CXCR5 staining in lung alveolar cells was usually not detected (3/6) or low (3/6), the LUAD tissues displayed not detected (10/22), low (6/22), medium (4/22) or high (2/22) CXCR1 staining. (Fig. 3).

Fig. 3
figure 3

Protein expressions of CXCR family members in patients with NSCLC and normal lung tissues. Compared to normal samples protein expressions of CXCR1, CXCR3 and CXCR5 were significantly down-regulated and the protein expressions of CXCR7 was up-regulated

Function enrichment of CXCR family members in LUAD

In this study, core proteins interacting with CXCR family members were identified through the GeneMANIA platform including NTSR2, P2RY10, CCR2, CCR4, CXCL13, CCR7, CCR8, CXCL5, ADRA1A, CD6, CCR1, PPBP, CCR3, CXCL3, PF4, ACKR2, CXCL2, FPR2, GPR183 and CCRL2. CXCL2, FPR2, GPR183, and CCRL2. and an interaction network was constructed using CXCR family members and their associated core genes to determine their relationships. For example, NTSR2 and CXCR7 share the same protein structural domain. CXCR1 and NTSR2 have a co-expression relationship. (Fig. 4A).

Fig. 4
figure 4

Function enrichment of CXCRs family members. A Network of CXCRs and their 20 related genes was analyzed. B Cellular component; C Biological processes; D Molecular functions; E KEGG pathway analysis

This study then performed GO function and pathway analysis of the above gene set via DAVID to predict the subsequent complex biological functions and signaling pathways. The biological processes such as G-protein coupled receptor signaling pathway, chemokine-mediated signaling pathway, inflammatory response, chemotaxis, immune response and positive regulation of cytosolic calcium ion concentration were remarkably regulated by the CXCRs in LUAD. (Fig. 4C) The cellular components including plasma membrane, integral component of membrane andintegral component of plasma membrane. (Fig. 4B) Besides, CXCRs also prominently affected the molecular functions, such as G-protein coupled receptor activity, chemokine receptor activity, C–C chemokine receptor activity, chemokine activity and C–X–C chemokine receptor activity were significantly associated with the CXCRs alterations. (Fig. 4D).

In KEGG analysis, we found that chemokine signaling pathways, cytokine–cytokine receptor interactions and neuroactive ligand—receptor interactions may have relevance to CXCR function in LUAD. (Fig. 4E).

Correlation of mRNA expression levels of CXCRs with immune infiltration levels in patients with LUAD

By examining the level of immune infiltration in the TIMER database, this study found that the mRNA expression levels of all CXCR family members correlated significantly with tumor purity and neutrophils in LUAD. The mRNA expression of all CXCR family members except CXCR7 correlated with B cells, CD8 + T cells, macrophages and dendritic cells at significant levels. The mRNA expression levels of all CXCR family members except CXCR1 correlated statistically significant levels with CD4 + T cells. (Table 2) In summary, CXCR7 was not significantly associated with the level of immune cell infiltration, whereas the other CXCRs were significantly and positively associated with the level of immune cells.

Table2 The mRNA expression of CXCR family members was significantly related to the immune infiltration level in LUAD

Correlation of mRNA expression levels of CXCRs with clinicopathological features in patients with LUAD

This study then analyzed the relationship between mRNA expression levels of CXCR family members and clinicopathological parameters of LUAD patients, including individual cancer stage, cancer subtype and patient gender, using TCGA samples.

The results of the analysis of TNM stage in LUAD patients showed that the mRNA expression levels of CXCR4, CXCR5 and CXCR6 gradually increased with the progression of TNM stage. Conversely, CXCR7 showed a gradual decrease in its expression level with the progression of TNM stage, while the difference between CXCR1, CXCR2 and CXCR3 in the progression of TNM stage was not significant. (Fig. 5) The results of analyzing the mRNA expression levels of CXCR family members in relation to the gender of LUAD patients showed that CXCR3 and CXCR5 had statistically significantly lower mRNA expression levels in males than in females in LUAD, while the differences between CXCR1, CXCR2, CXCR4, CXCR6 and CXCR7 were not significant. (Fig. 6).

Fig. 5
figure 5

Association of mRNA expression of CXCRs family members with cancer stages of LUAD patients. The mRNA expression of CXCRs family members in normal people or in LUAD patients in different stages (*P < 0.05; **P < 0.01; ***P < 0.001)

Fig. 6
figure 6

Association of mRNA expression of CXCRs family members with gender of LUAD patients. The mRNA expression of CXCR family members in normal people or in LUAD patients in male or female (*P < 0.05; **P < 0.01; ***P < 0.001)

The relationship between mRNA expressions of CXCR family members with cell subtypes of LUAD patients was analyzed by R. In LUAD, the expression of CXCR/3/4/5/6 in magnoid always lower than both bronchioid and squamoid. On the contrary, while CXCR7 in magnoid higher than both bronchioid and squamoid. (Fig. 7).

Fig. 7
figure 7

Association of mRNA expression of CXCR family members with cancer subtype of LUAD patients. The mRNA expression of CXCR family members in normal people or in LUAD patients in different subtype (*P < 0.05; **P < 0.01; ***P < 0.001)

Prognostic value of mRNA expression levels of family members in patients with LUAD

Kaplan–Meier analysis of the relationship between CXCRs and prognostic status of LUAD patients showed the following results. CXCR1/2/3///4/5/6 may play a role as beneficial factors in the prognosis of LUAD patients with HR values less than 1 and their P values of 7.70e−02, 1.90e−02, 1.00e−02, 5.01e−03, 4.80e−03, 1.85e−03, respectively. Conversely CXCR7 may play a role as a risk factor in the prognosis of LUAD patients with an HR of 1.604 and its P value of 4.27e−03. (Fig. 8) All of the above results suggest that the mRNA expression levels of CXCR family members are significantly associated with the prognosis of LUAD patients, with CXCR7 acting as a drug target for LUAD patients.

Fig. 8
figure 8

Prognostic value of mRNA expression of CXCR family members in LUAD patients about OS. OS comparing the high and low expression of CXCRs family members in LUAD patients

Discussion

CXCRs is a superfamily of seven-transmembrane G protein-coupled receptors and play critical roles in immune surveillance, inflammation, tissue development and homeostasis. Increasing evidence suggested that various CXCRs were structural expression in various malignancies and be involved in the initiation, progression, and outcome of various human tumors.

In our study, we found that the low expression of CXCR1/2/3/4/5/6 in LUAD were associated with poor prognosis through TCGA dataset and multiple public databases, while low expression of CXCR7 was correlated with favorable survival in LUAD patients. We also studied the association between tumor immune infiltrating cells and CXCRs expression, and the results revealed that the expression levels of the seven CXCRs were significantly correlated with CD8 + T cell, CD4 + T cell, neutrophil, and dendritic cell, which were related to tumor progression, metastasis, or prognosis. [26]

The results of TCGA analysis showed that CXCR1 and CXCR2 were lowly expressed in LUAD, and patients with the cancer subtype of magnoid tended to express high level of CXCR1 mRNA. In addition to CXCL8, CXCR1 and CXCR2 were also activated by other CXC-chemokines, including CXCL1, − 2, − 3, − 5, − 6 and − 7. It is widely accepted that CXCR1 and CXCR2 are essential for the initiation and growth of human tumors, thus serve as the novel therapeutic targets for many solid tumors, including lung cancer, breast cancer, prostate cancer, ovarian cancer, liver cancer and melanoma. [27,28,29,30] However, our findings are inconsistent with these studies, the results from our enrichment and immune infiltration assays speculated that CXCR1 and CXCR2 may inhibit the progression of LUAD through their abilities to increase the level of immune activation. Given the current studies on CXCR1 or CXCR2 in LUAD remained sparse, CXCR1 and CXCR2 are worthy of further research and exploration.

We also found that the mRNA level of CXCR3 was higher, however, its protein level was down-regulated in LUAD. In addition, the mRNA expression of CXCR3 in female patients was significantly higher than that in male patients. Some evidence suggests that CXCR3 secreted by immune cells can inhibit the development of gastric cancer through paracrine pathway. [31] On the contrary, it was demonstrated that CXCR3 promotes the proliferation, migration and vascular invasion of cancer cells, such as breast cancer cells, gastric cancer cells. [32, 33] Thus, in different stages of tumor differentiation and development, we propose that CXCR3 plays different roles in the immune response of cancer cells.

In the current study, we also demonstrated that the expression levels of CXCR4/5 mRNA and protein were predominantly down-regulated in LUAD as compared to normal tissues, whereas there is no obvious difference in the expression of CXCR6 between the two groups. Furthermore, low expression of CXCR4/5 was significantly correlated with shorter survival, CXCR4/5 expression levels were also positively correlated with immune cell infiltration,and the expression levels of CXCR4 and CXCR5 were gradually decreased with the increase of TNM staging in patients with LUAD by TCGA dataset. Coupled with this, it was found that the overall survival ability of patients with lung cancer was significantly improved following the increase of CXCR4 expression in the tumor stroma. [34] A previous study also demonstrated that several leukocyte including recirculating B cells, small subsets of CD4 + and CD8 + T cells have the ability to express CXCR5, which in turn controls migration, entry and exit of these leukocytes to lymph nodes through interaction with its CXCL13 ligand. [35] In addition, invariant NKT (iNKT) cells is known to trigger potent antitumor responses in vivo due to their homeostasis and activation, in which CXCR6 and its ligand CXCL16 have been shown to play critical roles. [36] However, others suggested that CXCR4 and CXCR5 might play a significant role in the occurrence and development of breast cancer and prostate cancer. [37,38,39]

Our data also showed that the protein level of CXCR7 in LUAD was up-regulated compared with normal tissues, the high expression of CXCR7 was associated with poor OS, and the expression levels of CXCR7 increased gradually with the progress of TNM staging. Several studies have demonstrated that CXCR7 may be involved in tumor-associated signaling pathways, including PLC/MAPK, ERK1/2, STAT3 and AKT pathways, which have been revealed to play a prominent role in tumor cell adhesion, invasion and metastasis. [40,41,42]

To explore the potential mechanism of CXCR in LUAD, functional enrichment analyses were performed, according to the results of network, our results demonstrated that CXCRs and their related genes may be involved in the immune process, angiogenesis, and tumor initiation and progression via a variety of signaling pathways (e.g., Chemokine signaling pathway) and biological processes (e.g., immune response). According to the previous literature and our research results, we speculated that CXCRs, expressed by tumor cells and immune cells in the microenvironment, have the ability to regulate the growth, invasion and metastasis of LUAD in a variety of ways, in which the exact molecular mechanisms underlying these need to be further explored.

In summary, our results provided supportive evidence for undestanding the complexity of CXCR1-7 and their related biological functions, which may provide a valuable insight for the development of CXCRs-mediated targeted therapy. However, there were some limitations to our study. First, all the data was based on the online databases, further in vivo and in vitro studies are required to verify these findings. Second, the underlying mechanisms of CXCRs in LUAD is still unknown. In addition, the current study was only a retrospective study, future detailed prospective studies need to be further explored.

Conclusions

Our findings from public databases provided a unique insight into expression characteristics and prognostic values of CXCR members in LUAD, which would be benefit for the understanding of pathogenesis, diagnosis, prognosis prediction and targeted treatment in LUAD.

Availability of data and materials

All data generated or analyzed in this study are contained in the multiple public databases mentioned in the article, including Oncomine (https://www.oncomine.org/resource/login.html), Human Protein Atlas (http://www.proteinatlas.org/), UALCAN (http://ualcan.path.uab.edu/), DAVID (https://david.ncifcrf.gov/), GeneMANIA (http://www.genemania.org/), TIMER (https://cistrome.shinyapps.io/timer/) and The Cancer Genome Atlas (TCGA) database (https://portal.gdc.cancer.gov/).

Abbreviations

CXCR:

CXC chemokine receptors

LUAD:

Lung adenocarcinoma

TCGA:

The cancer genome atlas

TIMER:

The tumor immune estimation resource

OS:

Overall survival

NSCLC:

Non-small cell lung cancer

LUSC:

Lung squamous cell carcinoma

GO:

Gene ontology

KEGG:

Kyoto encyclopedia of genes and genomes

TNM:

Tumor node metastasis

CXCL:

CXC chemokine ligands

NKT:

Natural killer T cells

PLC:

Phospholipase C

MAPK:

Mitogen-activated protein kinase

ERK:

Extracellular signal-regulated kinase

STAT3:

Signal transducer and activator of transcription 3

AKT:

Protein kinase B

References

  1. Alberg AJ, Brock MV, Samet JM. Epidemiology of lung cancer: looking to the future. J Clin Oncol. 2005;23(14):3175–85.

    Article  Google Scholar 

  2. Molina JR, Yang P, Cassivi SD, Schild SE, Adjei AA. Non-small cell lung cancer: epidemiology, risk factors, treatment, and survivorship. Mayo Clin Proc. 2008;83(5):584–94.

    Article  Google Scholar 

  3. Travis WD, Brambilla E, Noguchi M, et al. International association for the study of lung cancer/American thoracic society/European respiratory society: international multidisciplinary classification of lung adenocarcinoma: executive summary. Proc Am Thorac Soc. 2011;8(5):381–5.

    Article  Google Scholar 

  4. Chen Z, Fillmore CM, Hammerman PS, Kim CF, Wong KK. Non-small-cell lung cancers: a heterogeneous set of diseases [published correction appears in Nat Rev Cancer. 2015 Apr;15(4):247]. Nat Rev Cancer. 2014;14(8):535–546.

  5. Zhang J, Fujimoto J, Zhang J, et al. Intratumor heterogeneity in localized lung adenocarcinomas delineated by multiregion sequencing. Science. 2014;346(6206):256–9.

    Article  CAS  Google Scholar 

  6. Engelman JA, Zejnullahu K, Mitsudomi T, et al. MET amplification leads to gefitinib resistance in lung cancer by activating ERBB3 signaling. Science. 2007;316(5827):1039–43.

    Article  CAS  Google Scholar 

  7. Turke AB, Zejnullahu K, Wu YL, et al. Preexistence and clonal selection of MET amplification in EGFR mutant NSCLC. Cancer Cell. 2010;17(1):77–88.

    Article  CAS  Google Scholar 

  8. Griffith JW, Sokol CL, Luster AD. Chemokines and chemokine receptors: positioning cells for host defense and immunity. Annu Rev Immunol. 2014;32:659–702.

    Article  CAS  Google Scholar 

  9. Caronni N, Savino B, Recordati C, Villa A, Locati M, Bonecchi R. Cancer and chemokines. Methods Mol Biol. 2016;1393:87–96.

    Article  Google Scholar 

  10. Chow MT, Luster AD. Chemokines in cancer. Cancer Immunol Res. 2014;2(12):1125–31.

    Article  CAS  Google Scholar 

  11. Moser B, Wolf M, Walz A, Loetscher P. Chemokines: multiple levels of leukocyte migration control. Trends Immunol. 2004;25(2):75–84.

    Article  CAS  Google Scholar 

  12. Rajagopal S, Rajagopal K, Lefkowitz RJ. Teaching old receptors new tricks: biasing seven-transmembrane receptors. Nat Rev Drug Discov. 2010;9(5):373–86.

    Article  CAS  Google Scholar 

  13. Bachelerie F, Graham GJ, Locati M, et al. New nomenclature for atypical chemokine receptors. Nat Immunol. 2014;15(3):207–8.

    Article  CAS  Google Scholar 

  14. Varney ML, Singh S, Li A, Mayer-Ezell R, Bond R, Singh RK. Small molecule antagonists for CXCR2 and CXCR1 inhibit human colon cancer liver metastases. Cancer Lett. 2011;300(2):180–8.

    Article  CAS  Google Scholar 

  15. Marchese A. Endocytic trafficking of chemokine receptors. Curr Opin Cell Biol. 2014;27:72–7.

    Article  CAS  Google Scholar 

  16. Lee HJ, Song IC, Yun HJ, Jo DY, Kim S. CXC chemokines and chemokine receptors in gastric cancer: from basic findings towards therapeutic targeting. World J Gastroenterol. 2014;20(7):1681–93.

    Article  CAS  Google Scholar 

  17. Ben-Baruch A. The multifaceted roles of chemokines in malignancy. Cancer Metastasis Rev. 2006;25(3):357–71.

    Article  CAS  Google Scholar 

  18. Li X, Li H, Yang C, Liu L, Deng S, Li M. Comprehensive Analysis of ATP6V1s Family Members in Renal Clear Cell Carcinoma With Prognostic Values. Front Oncol. 2020;10:567970.

    Article  Google Scholar 

  19. Rhodes DR, Kalyana-Sundaram S, Mahavisno V, et al. Oncomine 3.0: genes, pathways, and networks in a collection of 18,000 cancer gene expression profiles. Neoplasia. 2007;9(2):166–80.

    Article  CAS  Google Scholar 

  20. Tomczak K, Czerwińska P, Wiznerowicz M. The cancer genome atlas (TCGA): an immeasurable source of knowledge. Contemp Oncol. 2015;19(1A):A68–77.

    Google Scholar 

  21. Chandrashekar DS, Bashel B, Balasubramanya SAH, et al. UALCAN: a portal for facilitating tumor subgroup gene expression and survival analyses. Neoplasia. 2017;19(8):649–58.

    Article  CAS  Google Scholar 

  22. Thul PJ, Ã…kesson L, Wiking M, et al. A subcellular map of the human proteome. Science. 2017;356(6340):eaal3321.

    Article  Google Scholar 

  23. Warde-Farley D, Donaldson SL, Comes O, et al. The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res. 2010;38(Web Server issue):W214–20.

    Article  CAS  Google Scholar 

  24. da Huang W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.

    Article  CAS  Google Scholar 

  25. Li T, Fan J, Wang B, et al. TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. 2017;77(21):e108–10.

    Article  CAS  Google Scholar 

  26. Renner K, Singer K, Koehl GE, et al. Metabolic hallmarks of tumor and immune cells in the tumor microenvironment. Front Immunol. 2017;8:248.

    Article  Google Scholar 

  27. Müller A, Homey B, Soto H, et al. Involvement of chemokine receptors in breast cancer metastasis. Nature. 2001;410(6824):50–6.

    Article  Google Scholar 

  28. Sun YX, Wang J, Shelburne CE, et al. Expression of CXCR4 and CXCL12 (SDF-1) in human prostate cancers (PCa) in vivo. J Cell Biochem. 2003;89(3):462–73.

    Article  CAS  Google Scholar 

  29. Scotton CJ, Wilson JL, Milliken D, Stamp G, Balkwill FR. Epithelial cancer cell migration: A role for chemokine receptors? Cancer Res. 2001;61(13):4961–5.

    CAS  PubMed  Google Scholar 

  30. Li A, Varney ML, Singh RK. Expression of interleukin 8 and its receptors in human colon carcinoma cells with different metastatic potentials. Clin Cancer Res. 2001;7(10):3298–304.

    CAS  PubMed  Google Scholar 

  31. Fangfang C, Jingping Y, et al. Chemokine receptor CXCR3 correlates with decreased M2 MACROPHAGE INFILTRATION AND FAVORABLE PROGNOSIS in gastric cancer. BioMed Res Int. 2019;2019:6832867.

    Google Scholar 

  32. Zhou H, Wu J, Wang T, et al. CXCL10/CXCR3 axis promotes the invasion of gastric cancer via PI3K/AKT pathway-dependent MMPs production. Biomed Pharmacother. 2016;82:479–88.

    Article  CAS  Google Scholar 

  33. Bronger H, Karge A, Dreyer T, Zech D, Kraeft S, Avril S, Kiechle M, Schmitt M. Induction of cathepsin B by the CXCR3 chemokines CXCL9 and CXCL10 in human breast cancer cells. Oncol Let. 2017;13(6):4224–30. https://doi.org/10.3892/ol.2017.5994.

    Article  CAS  Google Scholar 

  34. Chen IX, Chauhan VP, Posada J, et al. Blocking CXCR4 alleviates desmoplasia, increases T-lymphocyte infiltration, and improves immunotherapy in metastatic breast cancer. In: Proceedings of the national academy of sciences of the United States of America, 2019.

  35. Haibi C, Sharma PK, Singh R, et al. PI3Kp110-, Src-, FAK-dependent and DOCK2-independent migration and invasion of CXCL13-stimulated prostate cancer cells. Mol Cancer. 2010;9(1):85.

    Article  Google Scholar 

  36. Kapur N, Mir H, Sonpavde GP, et al. Prostate cancer cells hyper-activate CXCR6 signaling by cleaving CXCL16 to overcome effect of docetaxel. Cancer Lett. 2019;454:1–13.

    Article  CAS  Google Scholar 

  37. Burger JA, Kipps TJ. CXCR4: a key receptor in the crosstalk between tumor cells and their microenvironment. Blood. 2006;107(5):1761–7.

    Article  CAS  Google Scholar 

  38. Bürkle A, Niedermeier M, Schmitt-Gräff A, Wierda WG, Keating MJ, Burger JA. Overexpression of the CXCR5 chemokine receptor, and its ligand, CXCL13 in B-cell chronic lymphocytic leukemia. Blood. 2007;110(9):3316–25.

    Article  Google Scholar 

  39. Cullen R, Germanov E, Shimaoka T, et al. Enhanced tumor metastasis in response to blockade of the chemokine Receptor CXCR6 Is overcome by NKT Cell activation. J Immunol. 2009;183(9):5807.

    Article  CAS  Google Scholar 

  40. Wani N, Nasser MW, Ahirwar DK, et al. C-X-C motif chemokine 12/C-X-C chemokine receptor type 7 signaling regulates breast cancer growth and metastasis by modulating the tumor microenvironment. Breast Cancer Res. 2014;16(3):R54.

    Article  Google Scholar 

  41. Chen Q, Zhang M, Li Y, et al. CXCR7 mediates neural progenitor cells migration to CXCL12 independent of CXCR4. Stem Cells. 2015;33(8):2574–85.

    Article  CAS  Google Scholar 

  42. Yu Y, Li H, Xue B, et al. SDF-1/CXCR7 axis enhances ovarian cancer cell invasion by MMP-9 expression through p38 MAPK pathway. DNA Cell Biol. 2014;33(8):543–9.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

We are grateful to the TCGA database for providing reliable data.

Funding

The present study was supported by the National Natural Science Foundation of China (Grant No. 81272854), Key Projects of Natural Science Foundation of Heilongjiang Province (Grant No. ZD2019H008), Excellent Innovation Team Construction Project of Basic Scientific Research Business Fee of Provincial Colleges and Universities in Heilongjiang Province (Grant No. 2019-KYYWF-1334), Young innovative talents training project of regular undergraduate colleges and universities in Heilongjiang Province (Grant No. UNPYSCT-2020054), Jiamusi University doctoral Special Research Fund launch project (JMSUBZ2019-01), Shenyang Science and Technology Plan Project (No.19-112-4-092).

Author information

Authors and Affiliations

Authors

Contributions

LTH was involved in writing the paper and processing the data, WQW and YL were involved in supervising the paper, and WJD was involved in collecting the literature and processing some of the data. ZSC, LS, CBZ, SZL, JRW, QSR, XYD, WDL, XBL, and YTD were involved in processing some of the data. All authors have read and approved the final manuscript.

Corresponding authors

Correspondence to Yue Li or Wei-Qun Wang.

Ethics declarations

Ethics approval and consent to participate

All methods were carried out in accordance with relevant guidelines and regulations.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no conflicts of interest.

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 http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) 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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hu, LT., Deng, WJ., Chu, ZS. et al. Comprehensive analysis of CXCR family members in lung adenocarcinoma with prognostic values. BMC Pulm Med 22, 259 (2022). https://doi.org/10.1186/s12890-022-02051-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12890-022-02051-6

Keywords