- Primary research
- Open Access
Analysis of therapeutic targets and prognostic biomarkers of CXC chemokines in cervical cancer microenvironment
Cancer Cell International volume 21, Article number: 399 (2021)
The tumor microenvironment (TME) has received an increasing amount of attention. CXC chemokines can regulate immune cell transport and tumor cell activity to exert anti-tumor immunity. However, studies on the expression and prognosis of CXC chemokines in cervical cancer (CC) are more limited.
The study investigated the role of CXC chemokines in TME of CC by using public databases. Moreover, quantitative real-time PCR (qRT-PCR) and immunohistochemistry (IHC) of CXC chemokines were performed to further verify.
The transcriptional levels of CXCL1/3/5/6/8/9/10/11/13/16/17 in CC tissues were significantly elevated while the transcriptional levels of CXCL12/14 were significantly reduced. We reached a consistent conclusion that the expression of CXCL9/10/11/13 was verified by quantitative real-time PCR and immunohistochemistry. Moreover, CC patients with low transcriptional levels of CXCL1/2/3/4/5/8 were significantly associated with longer overall survival (OS). The CCL family was related to CXC chemokines neighboring alteration. RELA, NFKB1, LCK and PAK2 were the key transcription factors and kinase targets of CXC chemokines, respectively. We also found there were significant correlations between the expression of CXCL9/10/11 and the infiltration of immune cells (CD8+ T cell, CD4+ T cell, neutrophils and dendritic cells).
In brief, we conducted a comprehensive analysis of CXC chemokines via clinical data and some online public databases. Our results may provide a new idea for the selection of immunotherapeutic targets and prognostic biomarkers for cervical cancer.
Cervical cancer accounts for approximately 12% of female cancers and is one of the major causes of death among women worldwide [1, 2]. Although the early detection and treatment of CC have improved, there are still about 604,127 new cases and 341,831 deaths in 2020, according to data from the World Health Organization . The actual situation shows that the majority of patients with cervical cancer are in the advanced stage of the disease, with limited access to appropriate treatment . As a result, the mortality rate is still high, and the median overall survival rate of advanced cervical cancer is only 16.8 months . Nowadays, an increasing number of studies have explored the kinase and immune checkpoint inhibitors in cancer . Moreover, epigenetics, the new regulation of specific genes, has also shown certain importance in the progression of cervical cancer . This provides a new idea for identifying the therapeutic targets and prognostic biomarkers of cervical cancer.
Chemokines are a type of secreted proteins with small molecular weights. Their role in autoimmune diseases, chronic inflammations and tumors has been continuously revealed , mainly mediating the migration of immune cells and the development of lymphoid tissues . In recent years, the role of chemokines in the tumor microenvironment has been continuously reported by researchers. As a major subfamily of the chemokines family, the function and mechanism of CXC chemokines have been discovered in tumors . CXC chemokines’ altered expression in cancers dictates immune cell recruitment, angiogenesis, tumorigenesis, cancer cell proliferation and metastasis [11, 12]. Previous studies have shown that there is a correlation between CXC chemokines, tumor microenvironment and tumor immunotherapy , which has been confirmed in some cancers [14,15,16]. This suggests that CXC chemokines may be potential therapeutic targets and prognostic biomarkers of cancer, by modulating tumor progression and immunotherapy efficacy. However, the function of the CXC chemokine family in CC has not been comprehensively described.
CXC chemokine is an important component of TME. Although previous studies have confirmed the expression and role of some members of the CXC chemokine family in cervical cancer, there is still a lack of comprehensive and systematic research. Therefore, it is worth exploring CXC chemokines as therapeutic targets and prognostic markers of cervical cancer. In this study, public databases were used to investigate the mRNA expression, prognosis, and related targets or kinase pathways of the CXC chemokine family in CC. Immunohistochemistry and qRT-PCR further verified the conclusion. Taken together, this study complements the function of CXC chemokines in cervical cancer, suggesting that certain CXC chemokines can be used as potential therapeutic targets and prognostic biomarkers for CC.
Materials and methods
In this study, clinical data and pathological specimens were collected retrospectively, and the patients’ informed consent was obtained before the pathological specimens were collected. Cervical cancer and para-carcinoma tissue samples for qRT-PCR and IHC analysis were taken from patients who underwent surgery from January 1, 2017 to December 31, 2018. We collected frozen tissues of 60 patients to analyze the mRNA levels of CXC chemokines by qRT-PCR and analyzed the protein expression levels of CXC chemokines by IHC staining in the paraffin tissues of 60 patients. This research work had been approved by the Academic Committee of The Third Clinical Medical College of Xinjiang Medical University (affiliated Tumor Hospital) and was carried out under the rules put forward in the Declaration of Helsinki. This study had the relevant informed consent exemption certificate issued by the Academic Committee. As for the public databases, neither ethics committee approval nor patient informed consent was needed for analyzing data.
Quantitative real-time PCR
Total RNA was isolated from tumor tissues and adjacent tissues using Trizol according to the manufacturer’s instructions. The extracted RNA was converted into cDNA with 5× primescript buffer, prime script RT enzyme mix I, oligo-dT primer and random 6 mers. The qRT-PCR was performed in the BioRad CFX96 Real-Time PCR Detection System machine in the presence of GAPDH, CXCL9, CXCL10, CXCL11, CXCL12 and CXCL13. We verified CXCL9/10/11/13 because these factors were expressed on the intersection of the three common databases including ONCOMINE, GEPIA and UALCAN. The detection of CXCL12 was because it was the only one that tended towards low expression in two databases. The transcription level of target genes was measured and normalized to GAPDH expression. The following primer sequences were used: GAPDH, 5′-GAAGGTGAAGGTCGGAGTC-3′ (forward) and 5′-GAAGATGGTGATGGGATTTC-3′ (reverse); CXCL9, 5′-TGAGAAAGGGTCGCTGTTCC-3′ (forward) and 5′-GGGCTTGGGGCAAATTGTTT-3′ (reverse); CXCL10, 5′-TGCCATTCTGATTTGCTGCC-3′ (forward) and 5′-TGCAGGTACAGCGTACAGTT-3′ (reverse); CXCL11, 5′-CCCTGGGGTAAAAGCAGTGA-3′ (forward) and 5′-TAAGCCTTGCTTGCTTCGAT-3′ (reverse); CXCL12, 5′-AGATGCCCATGCCGATTCTT-3′ (forward) and 5′-AGGGCACAGTTTGGAGTGTT-3′ (reverse); CXCL13, 5′-CGACATCTCTGCTTCTCATGC-3′ (forward) and 5′-ACTGAGCTCTCTTGGACACAT-3′ (reverse).
Formalin-fixed paraffin-embedded surgical specimens were used for immunohistochemical study. The sections were dried at 60 °C for 2 h, subsequently were dewaxed in xylene and graded alcohols, were hydrated and washed in phosphate-buffered saline. After antigen repair was treated in a microwave oven (15 min in citrate buffer, pH 6.0), the endogenous peroxidase was inhibited with 3% H2O2 for 30 min, then the sections were incubated with 10% normal goat serum for 40 min. Primary antibodies composed of rabbit anti-CXCL9 antibody (bs-2551R [Bioss], 1:100), rabbit anti-CXCL10 antibody (bs-1502R [Bioss], 1:150), rabbit anti-CXCL11 antibody (DF9917 [Affinity], 1:150) and rabbit anti-CXCL13 antibody (bs-2553R [Bioss], 1:100) were applied overnight in a moist room at 4°C. Then the tissues were incubated with secondary antibody (37 °C 50 min), stained with diaminobenzidine, and counterstained with hematoxylin. Positive staining was evaluated using computer-aided image analysis and Image J software. The average CXC chemokines infiltration was determined from three randomized fields by two independent pathologists who were unaware of the patients’ pathological and clinical status.
Transcription-related databases of CXC chemokines in patients of cervical cancer
We used the public databases ONCOMINE (http://www.oncomine.org) , GEPIA (http://gepia.cancer-pku.cn/index.html)  and UALCAN (http://ualcan.path.uab.edu)  that could provide cancer RNA sequence data and clinical data to analyze the differential expression of CXC chemokines in cervical cancer and adjacent cancer or normal tissues by using Student’s t-test. The cut-off of the p-value was 0.05. In ONCOMINE, the fold change was 2.0, and the gene rank was in the top 10%. In brief, we entered the target genes in the input box of the database, and then searched for them. In addition, we conducted a prognostic study of CXC chemokines in cervical cancer by Kaplan–Meier curve in the GEPIA database.
cBioPortal (http://www.cbioportal.org) is an online open-access website, which involves the exploration, visualization, and analysis of multidimensional cancer genomics data . Genetic alterations of CXC chemokines were obtained from cBioPortal based on The Cancer Genome Atlas (TCGA) database. In this study, 293 cervical squamous cell carcinoma samples were analyzed (TCGA, PanCancer Atlas). The z-score of mRNA expression (log RNA Seq V2 RSEM) was obtained using the threshold of ± 2.0.
CXC chemokines-related networks
We studied the gene-related networks and protein–protein interaction (PPI) of CXC chemokines by GeneMANIA (http://www.genemania.org) , a website which could provide information on the protein and genetic interactions, pathways and co-expression of submitted genes, and the STRING database (https://string-db.org/) , respectively. In our research, we entered the searched species and gene names in the input box of databases. Moreover, on the right side of the GeneMANIA website, we also set up bioinformatics methods such as co-expression, physical interaction, gene enrichment analysis, predictive interaction and pathway, etc.
Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis
DAVID 6.8 (https://david.ncifcrf.gov/home.jsp), a website that can clarify the biological functions of genes , was used to perform GO and KEGG pathway enrichment analysis of CXC chemokines and their neighboring 50 genes. Analysis of GO function mainly included biological process (BP), cell component (CC), molecular function (MF). Then, the “ggplot2” package of R software (4.0.2) was used for visualization. P-value < 0.05 was considered to be statistically different.
Target analysis of kinases and transcription factors
The TRRUST (https://www.grnpedia.org/trrust/) database can provide transcription factor (TF)–target regulatory relationship . The module of LinkInterpreter in the LinkedOmics (http://www.linkedomics.org/)  was used to obtain biological insights on the enrichment of kinase targets of CXC chemokines. In the LinkedOmics database, Gene Set Enrichment Analysis (GSEA) was investigated with a minimum number of genes (size) of 3 and a simulation of 500. The results were analyzed by the Spearman correlation test. The p-value cut-off was 0.05.
TIMER (https://cistrome.shinyapps.io/timer/) could provide a systematic assessment of the infiltration of different immune cells and their clinical effects . We conducted the module of correlation in the TIMER database to evaluate the interrelation between immune cell infiltration and CXC chemokines level (or CXCRs) by the purity-corrected partial Spearman method in the tumor microenvironment. Moreover, the correlations between cell markers of CD8+ T cells (CD8A and CD8B), natural killer (NK) cells (KLRK1, KIR2DL4, KIR3DL2, NCR1, and NCR3), T helper 1 (Th1) cells (TBX21, STAT1) and CXCL9–11/CXCR3 were explored.
The data of clinical patients were presented as mean ± SD, the Chi-square test and Wilcoxon rank-sum test was used to compare the data between the tumor group and the para-cancerous group (SPSS 26.0). The qRT-PCR and immunohistochemical data of 60 patients were analyzed using Student’s t-test and Mann–Whitney U test (Graphpad Prism 8.0). Differences with p values < 0.05 were considered statistically significant.
The clinicopathological characteristics of the patients
The clinicopathological features of 60 patients were studied. The results are shown in Table 1. There were no differences in ethnicity (p = 0.566), tumor size (p = 0.756), differentiation (p = 1.000), or FIGO stage (p = 0.378) between the tumor tissue group and the adjacent tissue group, except for first diagnostic age (p = 0.000).
Abnormal expression of CXC chemokines in CC patients
To investigate the transcription level of CXC chemokines between tumor and adjacent or normal tissues in CC, we performed an analysis using the ONCOMINE, GEPIA and UALCAN database. Sixteen CXC chemokines were retrieved using the ONCOMINE databases. The results are presented in Fig. 1 and Table 2. The transcriptional level of CXCL1/3/5/6/8/9/10/11/13/16 in cervical cancer tissues was significantly elevated, while the transcriptional level of CXCL12/14 was significantly lower than that of adjacent cancer tissue. In Zhai cervix statistics , CXCL1/3/5/6/8/13 were overexpressed in cervical tumor tissues versus tumor-adjacent tissues. In Scotto Cervix 2 Statistics , CXCL1/8/9/10/11/13 were elevated in CC tissues relative to the adjacent tissues of the tumor. Similarly, CXCL8/9/10/11 were overexpressed in CC tissues instead of para-carcinoma tissue in Biewenga Cervix Statistics . In the GEPIA database, the results indicated that the expressional levels of CXCL1/8/9/10/11/13/16/17 were increased in tumor tissues rather than normal tissues, while CXCL12 was reduced (Fig. 2a). In the UALCAN database, as expected, the transcriptional levels of CXCL6 (p = 1.18e−4), CXCL9 (p = 1.63e−12), CXCL10 (p = 7.77e−16), CXCL11 (p = 8.36e−13), CXCL13 (p = 2.81e−7) and CXCL17 (p = 1.74e−3) in cervical tissues were significantly elevated (Fig. 2b–g). Since the normal control of the study in the ONCOMINE database came from para-carcinoma tissues rather than normal tissues, it was rational to acquire diverse results from the three databases. In addition, the sample size was not sufficient to capture variability as there were only three normal cervical patients in the UALCAN database. This part of the data has yet to be further confirmed in clinical practice. Taken together, these data suggest that these CXC chemokines play a significant role in the tumorigenesis and progression of cervical cancer.
The mRNA expression of CXCL9/10/11/12/13 in cervical cancer tissues and adjacent cancer tissues was also examined by qRT-PCR (Fig. 3). We verified CXCL9/10/11/13 because these factors were expressed at the intersection of the three databases. CXCL12 was detected owing to it was the only one with a tendency of low expression in the two databases. CXCL9 (p = 0.027), CXCL10 (p < 0.001), CXCL11 (p = 0.002) and CXCL13 (p < 0.001) were significantly up-regulated in tumor tissues compared with tumor-adjacent tissues, while the expression of CXCL12 (p = 0.149) was down-regulated. We further performed IHC on CXCL9/10/11/13 to detect their protein expression in cervical cancer tissues and adjacent tissues. We found that the protein expression CXCL9 (p = 0.012), CXCL10 (p = 0.003), CXCL11 (p = 0.011), CXCL13 (p = 0.004) were higher in tumor tissues than that in adjacent tissues (Fig. 4). These were consistent with the findings in the databases and also further verified the transcriptional expression of CXCL9/10/11/13 in cervical cancer.
The prognostic value of CXC chemokines in cervical cancer
We conducted an analysis in the GEPIA database to explore the value of differentially expressed CXC chemokines in the progression of CC by evaluating the correlation between differentially expressed CXC chemokines and clinical outcomes. The low transcriptional levels of CXCL1 (p = 0.033), CXCL2 (p = 0.046), CXCL3 (p = 0.017), CXCL4 (p = 0.027), CXCL5 (p = 0.011) and CXCL8 (p = 1.5e−5) were significantly associated with longer OS of CC patients (Fig. 5).
Genetic alteration of CXC chemokine in cervical cancer
We analyzed the genetic alterations of differentially expressed CXC chemokines by using the TCGA dataset in the cBioPortal. We found that a high mutation rate of CXC chemokines (42%) was observed in CC patients. Among 293 cervical cancer patients, the genetic alteration was found in 124 cervical cancer patients. As a result, CXCL1, CXCL2, CXCL3, CXCL4, CXCL5, CXCL6, CXCL7, CXCL8, CXCL9, CXCL10, CXCL11, CXCL12, CXCL13, CXCL14, CXCL6 and CXCL17 were altered in 6%, 3%, 4%, 4%, 6%, 5%, 5%, 5%, 6%, 6%, 5%, 5%, 6%, 5%, 7% and 6% of the cervical cancer samples, respectively (Fig. 6a). In addition, we compared the overall survival, disease-free survival and progress-free survival between the CXC chemokine mutant group and the non-mutant group. Unfortunately, there were no difference in overall survival (p = 0.895), disease-free survival (p = 0.233) and progress-free survival (p = 0.175) between the CXC chemokine mutant group and the non-mutant group (Fig. 6b, c).
Analysis of adjacent gene network and interaction of CXC chemokines in cervical cancer
To further explore the potential interactions between differentially expressed CXC chemokines and adjacent genes, we performed a PPI network analysis in STRING (Fig. 7a). The results revealed that the function of the differentially expressed CXC chemokines was related to chemokine signaling pathways and inflammation response. Results on the GeneMANIA website also showed the differential expressed CXC chemokines’ function was primarily connected with cell chemotaxis, chemokine receptor binding, and chemokine activity (Fig. 7b).
Predicting the function and pathway of CXC chemokines in patients with cervical cancer
The function and pathway of the differentially expressed CXC chemokines and their neighboring genes were investigated using DAVID 6.8. As shown in Fig. 8a, GO term analysis showed that the differentially expressed CXC chemokines were mainly located in the extracellular space, extracellular region, cell, external side of plasma membrane and cell surface, where they participated in chemotaxis, inflammatory response, monocyte chemotaxis, cellular response to interleukin-1 and immune response. They acted as CCR chemokine receptor binding, CXCR chemokine receptor binding, heparin binding, chemokine receptor binding and chemoattractant activity. KEGG pathway analysis showed enrichment of chemokine signaling pathway, cytokine–cytokine receptor interaction, rheumatoid arthritis, Toll-like receptor (TLR) signaling pathway, TNF signaling pathway (Fig. 8b).
Analysis of kinase and transcription factor targets in patients with cervical cancer
The TRRUST and Linkedomics databases were used to explore the possible kinase and TF targets of differential CXC chemokines in the cervical cancer microenvironment. As shown in Table 3, RELA and NFKB1 were the key TFs for CXCL1/2/5/8/10/12 (p = 1.09e−07 and p = 1.14e−07, respectively), SP1 was the key TF for CXCL1/5/14 (p = 0.00683) in the TRRUST database. Moreover, the top two kinase targets of CXC chemokines were studied from the LinkedOmics database (Table 4). Src family tyrosine kinases (LYN and LCK) were the primary kinase target of CXCL10/11/16. The kinase target of CXCL5 was the PAK family kinases (PAK2 and PAK3).
CXC chemokines and immune cell infiltration
As a component of the tumor microenvironment, CXC chemokines were closely knitted with immune cells. Since the experimental validation was performed only on CXCL9/10/11/13, the current work primarily focused on the relationship between CXCL9/10/11 and their corresponding receptor CXCR3 in the TIMER database. The expression of CXCL9/10/11 and CXCR3 was negatively correlated with tumor purity, indicating that they were highly expressed in infiltrating immune cells rather than tumor cells in TME. The expression of CXCL9/10/11 was positively associated with the infiltration of CD8+ T cell, CD4+ T cell, neutrophils and dendritic cells (all p < 0.05, Additional file 1: Figure S1a–c). Moreover, there was a positive correlation between the expression of CXCR3 and the infiltration of B cell, CD8+ T cell, CD4+ T cell, macrophages, neutrophils and dendritic cells (all p < 0.05, Additional file 1: Figure S1d). It is universally accepted that CXCR3 is commonly expressed in CD8+ T cells, NK cells, and Th1 cells. We further evaluated the molecular markers of CXCL9/10/11, CXCR3 and CD8+ T cells, NK cells and Th1 cells. We found that the expression of CXCL9/10/11 and CXCR3 was positively correlated with these cell markers (Fig. 9).
Although there have been some improvements in the early detection and treatment of cervical cancer, it is still the leading cause of cancer-related deaths among females. CXC chemokines, as a component of TME, affect the occurrence and development of cancer. In our study, we comprehensively analyzed the prognostic value and biological function of CXC chemokines in CC using public databases. We found that the expression of CXCL1/3/5/6/8/9/10/11/13/16/17 was up-regulated in tumor tissues, while the expression of CXCL12/14 was down-regulated. The OS of CC patients with low transcription levels of CXCL1/2/3/4/5/8 was longer. In addition, we noticed that the CCL family was related to CXC chemokines neighboring alteration. The enrichment analysis of CXC chemokines was mainly connected with the chemokine signaling pathway, etc., RELA, NFKB1, LCK and PAK2 were the key TF and kinase targets of CXC chemokines, respectively. In short, these data indicate that differentially expressed CXC chemokines may play an important role in CC.
We have not reached a consistent conclusion on the expression of CXC chemokines in the three public databases. Considering that the normal control of this study in the ONCOMINE database was from cancer adjacent tissues rather than normal tissues, and there were only three normal cervical samples in the UALCAN database that were not sufficient to capture variability. Therefore, it was rational to obtain different results from the three databases. In addition, as for qRT-PCR analysis, the expression of CXCL9/10/11/12/13 was consistent with the data in the database. Unfortunately, there was no statistically significant difference in the expression of CXCL12 between tumor and para-cancerous tissues. One possible interpretation is that cervical expression of CXCL12 is heterogeneous, tending to be inhibited as the tumor grows .
The occurrence and development of cervical cancer is a complex and multi-faceted process. We investigated the molecular characteristics of CXC chemokines that were significantly differentially expressed in CC. The differentially expressed CXC chemokines in CC had frequent genetic alterations with a mutation rate of 42%, of which increased mRNA expression was the most alteration. The CCL family was associated with CXC chemokines neighboring alteration, including CCL1, CCL2, CCL8 and CCL21, etc. There were different but similar functions performed by these CCL factors in the tumor. It has been reported that CCL1, CCL2 and CCL21 affect tumor metastasis and progression through direct or indirect interactions with lymph nodes [31,32,33]. CCL8 was an independent prognostic factor for patients with CC, and increased secretion of CCL8 promoted tumor progression [34, 35]. These results revealed that CXC chemokines may inhibit or promote tumor growth by influencing the alterations of neighboring genes in the tumor microenvironment.
Then via GO enrichment analysis and KEGG pathway enrichment analysis, we found that the functions of these genes were primarily related to the chemokine signaling pathway, cytokine–cytokine receptor interaction, and the Toll-like receptor signaling pathway. Previous studies have shown that chemokine signaling pathways affect tumor proliferation, angiogenesis and metastasis [8, 36]. For example, the increase of CXCL1/2/8 level regulated by AKIP 1 played an important role in the angiogenesis and growth of cervical cancer . However, CXCL9/10/11, as anti-angiogenic chemokines, inhibited the angiogenesis of cervical cancer through its receptor CXCR3-B . CXCL13 and CXCR5 have also been revealed in some cancers. The hypermethylation of a single CpG dinucleotide in the promoter region of the CXCL13 gene promoted the migration of cervical cancer cells . Silencing of CXCL12 altered the homeostatic autocrine and paracrine CXCR4 signaling to endocrine communication in the tumor microenvironment, leading to tumor progression and metastases . In addition, TLRs are I intermembrane proteins, and TLR signaling pathways in tumor cells have been shown to influence cancer progression. Stimulating TLRs in tumor cells fostered chronic inflammation that drove cancer cell proliferation, migration and angiogenesis, and established a tumor microenvironment that impaired the immune system and facilitated tumor growth [39, 40]. These data indicated that the differentially expressed CXC chemokines in CC were potential therapeutic targets.
By analyzing the differential expression of CXC chemokines in cervical cancer, we identified three key transcription factors (RELA, NFKB1 and SP1). RELA phosphorylation resulted in cancer progression by regulating NF-κB signaling . Some studies have revealed an association between elevated RELA expression and poor survival in cancer . Moreover, NFKB1 was a tumor suppressor by inhibiting cell proliferation, colony formation and migration in cervical cancer, while the mutation could weaken the tumor-suppressing function of NFKB1 [43, 44]. Sp1 may lead to radioresistance by targeting CDK1 to inhibit G2/M phase block, thereby affecting the prognosis of patients . In addition, our data suggested that LCK, LYN, PAK2 and PAK3 may be kinase targets for differential expression of CXC chemokines. It was related to Src family tyrosine kinases and PAK family kinases. These kinases were involved in the occurrence and development of tumors by regulating the proliferation, migration, invasion and apoptosis of tumor cells [46,47,48]. All in all, differentially expressed CXC chemokines may affect tumor development and progression by regulating these transcription factors and kinases in cervical cancer.
In the TME of cervical cancer, CXC chemokines could control the migration and localization of immune cells, and immune cells could be a factor of immunotherapy and clinical outcome by affecting tumor development and progression [49, 50]. It was well known that CD8+ T cells and NK cells can be recruited via the CXCL9–11/CXCR3 axis, resulting in tumor-suppressing with a form of paracrine . In our study, CXCL9/10/11 and their ligand CXCR3 were positively related to CD8+ T cells and NK cells. The results suggested that CXCL9/10/11 in the tumor microenvironment may recruit CD8+ T cells and NK cells through the CXCL9–11/CXCR3 axis, thereby improving the survival of cervical cancer patients by inhibiting the development and progression of tumors. Drugs that augmented paracrine CXCL9/10/11 expression and deactivated CXCR3 expression on tumor cells, have shown anti-tumor activity in several tumor models. At present, the treatment of anti-CXCR 3 is promising. CXCR3 antagonist has been shown to inhibit the implantation and growth of tumor cells in vitro, and inhibit lung metastasis in a vivo model . However, it is still worth exploring in cervical cancer. In short, these results indicated that CXC9/10/11 could be regarded as anti-tumor target by attracting immune cells.
In our study, experiments and several online bioinformatic platforms have systematically analyzed the expression, mutations, and related pathway value of CXC chemokines in CC. We look forward to paying more attention to the therapeutic effect of CXC chemokines in cervical cancer in future research. Through a more profound and comprehensive exploration of the key transcription factors, kinase targets and immune cells, the mechanism of CXC chemokines in the tumor microenvironment can be better revealed. There were some limitations in our study. First of all, due to the limitation of the online database, there were only three normal cervical patients in the UALCAN database, which was insufficient to capture variability. Therefore, clinical studies are needed to verify the conclusions. Second, analysis on the transcriptional and translational level can certainly reflect some aspects of immune status, but not global changes. Third, we should collect more clinical specimen data for experimental verification. Since this study only verified the mRNA and protein expression of some CXC chemokines, more indicators are supposed to be further verified in the experiment.
In conclusion, we performed a comprehensive analysis of CXC chemokines via clinical data and some online tumor databases. Patients with low transcriptional levels of CXCL1/2/3/4/5/8 were significantly associated with better prognosis, suggesting a role in tumorigenesis and development. CXCL9–11 were highly correlated with CD8+ T cells, NK cells, and Th1 cells in the tumor microenvironment, indicating that they can play a certain role in the anti-tumor immune response. In addition, we found that RELA, NFKB1 and SP1 were the key transcription factors of CXC chemokines, and LCK, LYN, PAK2, and PAK3 were the kinase targets of CXC chemokines. These results suggested that the interaction between these factors and CXC chemokines was involved in the mechanism of progression in CC. We hope that the differentially expressed CXC chemokine data could provide ideas and insights on the selection of prognostic markers and immunotherapy for cervical cancer.
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article.
Quantitative real-time PCR
The Cancer Genome Atlas
Kyoto Encyclopedia of Genes and Genomes
Gene Set Enrichment Analysis
T helper 1
GBD 2017 Causes of Death Collaborators. Global, regional, and national age-sex-specific mortality for 282 causes of death in 195 countries and territories, 1980–2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet. 2018;392(10159):1736–88.
Naga ChP, Gurram L, Chopra S, Mahantshetty U. The management of locally advanced cervical cancer. Curr Opin Oncol. 2018;30(5):323–9.
Globocan, IAfRo. Cervix uteri fact sheet: Cervix uteri. 2020.
Small W, Bacon MA, Bajaj A, Chuang LT, Fisher BJ, Harkenrider MM, Jhingran A, Kitchener HC, Mileshkin LR, Viswanathan AN, et al. Cervical cancer: a global health crisis. Cancer. 2017;123(13):2404–12.
Orbegoso C, Murali K, Banerjee S. The current status of immunotherapy for cervical cancer. Rep Pract Oncol Radiother. 2018;23(6):580–8.
Menderes G, Black J, Schwab CL, Santin AD. Immunotherapy and targeted therapy for cervical cancer: an update. Expert Rev Anticancer Ther. 2016;16(1):83–98.
Zhu H, Zhu H, Tian M, Wang D, He J, Xu T. DNA methylation and hydroxymethylation in cervical cancer: diagnosis, prognosis and treatment. Front Genet. 2020;11(4):1–12.
Bikfalvi A, Billottet C. The CC and CXC chemokines: major regulators of tumor progression and the tumor microenvironment. Am J Physiol Cell Physiol. 2020;318(3):C542–54.
Nagarsheth N, Wicha MS, Zou W. Chemokines in the cancer microenvironment and their relevance in cancer immunotherapy. Nat Rev Immunol. 2017;17(9):559–72.
Zhu Q, Han X, Peng J, Qin H, Wang Y. The role of CXC chemokines and their receptors in the progression and treatment of tumors. J Mol Histol. 2012;43(6):699–713.
Santoni M, Bracarda S, Nabissi M, Massari F, Conti A, Bria E, Tortora G, Santoni G, Cascinu S. CXC and CC chemokines as angiogenic modulators in nonhaematological tumors. Biomed Res Int. 2014. https://doi.org/10.1155/2014/768758.
Lee NH, Nikfarjam M, He H. Functions of the CXC ligand family in the pancreatic tumor microenvironment. Pancreatology. 2018;18(7):705–16.
Poeta VM, Massara M, Capucetti A, Bonecchi R. Chemokines and chemokine receptors: new targets for cancer immunotherapy. Front Immunol. 2019;10(3):1664–3224.
Cabrero-De Las Heras S, Martínez-Balibrea E. CXC family of chemokines as prognostic or predictive biomarkers and possible drug targets in colorectal cancer. World J Gastroenterol. 2018;24(42):4738–49.
Zeng Q, Sun S, Li Y, Li X, Li Z, Liang H. Identification of therapeutic targets and prognostic biomarkers among CXC chemokines in the renal cell carcinoma microenvironment. Front Oncol. 2020;9:1555.
Chen E, Qin X, Peng K, Xu X, Li W, Cheng X, Tang C, Cui Y, Wang Z, Liu T. Identification of potential therapeutic targets among CXC chemokines in breast tumor microenvironment using integrative bioinformatics analysis. Cell Physiol Biochem. 2018;45(5):1731–46.
Rhodes DR, Yu J, Shanker K, Deshpande N, Varambally R, Ghosh D, Barrette T, Pandey A, Chinnaiyan AM. ONCOMINE: a cancer microarray database and integrated data-mining platform. Neoplasia. 2004;67(1):10163–72.
Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;45(W1):W98–102.
Chandrashekar DS, Bashel B, Balasubramanya SAH, Creighton CJ, Ponce-Rodriguez I, Chakravarthi BVSK, Varambally S. UALCAN: a portal for facilitating tumor subgroup gene expression and survival analyses. Neoplasia (United States). 2017;19(8):649–58.
Gao J, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO, Sun Y, Jacobsen A, Sinha R, Larsson E, et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal. 2013;6(269):Pl1.
Warde-Farley D, Donaldson SL, Comes O, Zuberi K, Badrawi R, Chao P, Franz M, Grouios C, Kazi F, Lopes CT, et al. The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res. 2010;38(SUPPL. 2):214–20.
Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, Simonovic M, Doncheva NT, Morris JH, Bork P, et al. STRING v11: protein–protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(D1):D607–13.
da Huang W, Sherman BT, Zheng X, Yang J, Imamichi T, Stephens R, Lempicki RA. Extracting biological meaning from large gene lists with DAVID. Curr Protoc Bioinform. 2009;27:1–13.
Han H, Cho JW, Lee S, Yun A, Kim H, Bae D, Yang S, Kim CY, Lee M, Kim E, et al. TRRUST v2: an expanded reference database of human and mouse transcriptional regulatory interactions. Nucleic Acids Res. 2018;46(D1):D380–6.
Vasaikar SV, Straub P, Wang J, Zhang B. LinkedOmics: analyzing multi-omics data within and across 32 cancer types. Nucleic Acids Res. 2018;46(D1):D956–63.
Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, Li B, Liu XS. TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. 2017. https://doi.org/10.1158/0008-5472.CAN-17-0307.
Zhai Y, Kuick R, Nan B, Ota I, Weiss SJ, Trimble CL, Fearon ER, Cho KR. Gene expression analysis of preinvasive and invasive cervical squamous cell carcinomas identifies HOXC10 as a key mediator of invasion. Cancer Res. 2007;77(21):e108–10.
Scotto L, Narayan G, Nandula SV, Arias-Pulido H, Subramaniyam S, Schneider A, Kaufmann AM, Wright JD, Pothuri B, Mansukhani M, et al. Identification of copy number gain and overexpressed genes on chromosome arm 20q by an integrative genomic approach in cervical cancer: potential role in progression. Genes Chromosomes Cancer. 2008;47(9):755–65.
Biewenga P, Buist MR, Moerland PD, Ver Loren van Themaat E, van Kampen AH, ten Kate FJ, Baas F. Gene expression in early stage cervical cancer. Gynecol Oncol. 2008;108(2008):520–6.
Yadav SS, Prasad SB, Prasad CB, Pandey LK, Pradhan S, Singh S, Narayan G. CXCL12 is a key regulator in tumor microenvironment of cervical cancer: an in vitro study. Clin Exp Metas. 2016;33(5):431–9.
Das S, Sarrou E, Podgrabinska S, Cassella M, Mungamuri SK, Feirt N, Gordon R, Nagi CS, Wang Y, Entenberg D, et al. Tumor cell entry into the lymph node is controlled by CCL1 chemokine expressed by lymph node lymphatic sinuses. J Exp Med. 2013;210(8):1509–28.
Kawai Y, Kaidoh M, Yokoyama Y, Sano K, Ohhashi T. Chemokine CCL2 facilitates ICAM-1-mediated interactions of cancer cells and lymphatic endothelial cells in sentinel lymph nodes. Cancer Sci. 2009;100(3):419–28.
Wang SQ, Zhang F, Peng LI. Relationship between expressions of chemokine CCR7, CCL21 and lymphatic metastasis in cervical carcinoma. Matern Child Health Care China. 2010.
Chen XJ, Deng YR, Wang ZC, Wei WF, Zhou CF, Zhang YM, Yan RM, Liang LJ, Zhong M, Liang L, et al. Hypoxia-induced ZEB1 promotes cervical cancer progression via CCL8-dependent tumour-associated macrophage recruitment. Cell Death Dis. 2019;10(7):1–11.
Budczies J, Denkert C, Győrffy B, Schirmacher P, Stenzinger A. Chromosome 9p copy number gains involving PD-L1 are associated with a specific proliferation and immune-modulating gene expression program active across major cancer types. BMC Med Genomics. 2017;10(1):74.
Roy I, Getschman AE, Volkman BF, Dwinell MB. Exploiting agonist biased signaling of chemokines to target cancer. Mol Carcinog. 2017;56(3):804–13.
Zhang W, Wu Q, Wang C, Yang L, Liu P, Ma C. AKIP1 promotes angiogenesis and tumor growth by upregulating CXC-chemokines in cervical cancer cells. Mol Cell Biochem. 2018;448(1–2):311–20.
Ma D, Fan SB, Hua N, Li GH, Chang Q, Liu X. Hypermethylation of single CpG dinucleotides at the promoter of CXCL13 gene promoting cell migration in cervical cancer. Curr Cancer Drug Targets. 2020;20(5):355–63.
McCall KD, Muccioli M, Benencia F. Toll-like receptors signaling in the tumor microenvironment. Adv Exp Med Biol. 2020:1223:81–97.
Tran TH, Tran TTP, Truong DH, Nguyen HT, Pham TT, Yong CS, Kim JO. Toll-like receptor-targeted particles: a paradigm to manipulate the tumor microenvironment for cancer immunotherapy. Acta Biomater. 2019;94:82–96.
Lu X, Yarbrough WG. Negative regulation of RelA phosphorylation: emerging players and their roles in cancer. Cytokine Growth Factor Rev. 2015;26(1):7–13.
Giles KM, Brown RA, Ganda C, Podgorny MJ, Candy PA, Wintle LC, Richardson KL, Kalinowski FC, Stuart LM, Epis MR, et al. microRNA-7-5p inhibits melanoma cell proliferation and metastasis by suppressing RelA/NF-κB. Oncotarget. 2016;7(22):31663–80.
Cartwright T, Perkins ND, Wilson CL. NFKB1: a suppressor of inflammation, ageing and cancer. FEBS J. 2016;283(10):1812–22.
Yang D, Zhang W, Liang J, Ma K, Chen P, Lu D, Hao W. Single cell whole genome sequencing reveals that NFKB1 mutation affects radiotherapy sensitivity in cervical cancer. Oncotarget. 2018;9(7):7332–40.
Deng YR, Chen XJ, Chen W, Wu LF, Jiang HP, Lin D, Wang LJ, Wang W, Guo SQ. Sp1 contributes to radioresistance of cervical cancer through targeting g2/m cell cycle checkpoint CDK1. Cancer Manag Res. 2019;11:5835–44.
Montero JC, Seoane S, Ocaña A, Pandiella A. Inhibition of Src family kinases and receptor tyrosine kinases by dasatinib: possible combinations in solid tumors. Clin Cancer Res. 2011;17(17):5546–52.
Bommhardt U, Schraven B, Simeoni L. Beyond TCR signaling: emerging functions of Lck in cancer and immunotherapy. Int J Mol Sci. 2019;20(14):1–18.
Chung EY, Mai Y, Shah UA, Wei Y, Ishida E, Kataoka K, Ren X, Pradhan K, Bartholdy B, Wei X, et al. PAK kinase inhibition has therapeutic activity in novel preclinical models of adult T-Cell leukemia/lymphoma. Clin Cancer Res. 2019;25(12):3589–601.
Griffith JW, Sokol CL, Luster AD. Chemokines and chemokine receptors: positioning cells for host defense and immunity. Annu Rev Immunol. 2014;32(1):659–702.
Bindea G, Mlecnik B, Tosolini M, Kirilovsky A, Waldner M, Obenauf AC, Angell H, Fredriksen T, Lafontaine L, Berger A, et al. Spatiotemporal dynamics of intratumoral immune cells reveal the immune landscape in human cancer. Immunity. 2013;39(4):782–95.
Tokunaga R, Zhang W, Naseem M, Puccini A, Berger MD, Soni S, McSkane M, Baba H, Lenz HJ. CXCL9, CXCL10, CXCL11/CXCR3 axis for immune activation—a target for novel cancer therapy. Cancer Treat Rev. 2018;63:40–7.
Pradelli E, Karimdjee-Soilihi B, Michiels JF, Ricci JE, Millet MA, Vandenbos F, Sullivan TJ, Collins TL, Johnson MG, Medina JC, et al. Antagonism of chemokine receptor CXCR3 inhibits osteosarcoma metastasis to lungs. Int J Cancer. 2009;125(11):2586–94.
Thanks to Xinjiang Medical University for helping the authors with the manuscript. And we sincerely thank all the authors and study participants for their support in this study.
This work was supported by Grants from Natural Science Foundation of Xinjiang Uygur Autonomous region (2021D01C402).
Ethics approval and consent to participate
The current research work received approval from the Academic Committee of The Third Clinical Medical College of Xinjiang Medical University (affiliated Tumor Hospital) and was carried out under the rules put forward in the Declaration of Helsinki. This study had the relevant exemption certificate of informed consent issued by the Academic Committee.
Consent for publication
All the listed authors have actively participated in the study and have read and approved the submitted manuscript.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Correlation between the expression of CXCL9–11, CXCR3 and immune cells. Immune cells include B cell, CD8+T cell, CD4+T cell, macrophage, neutrophil and dendritic cell. a CXCL9. b CXCL10. c CXCL11. d CXCR3.
About this article
Cite this article
Kong, W., Zhao, G., Chen, H. et al. Analysis of therapeutic targets and prognostic biomarkers of CXC chemokines in cervical cancer microenvironment. Cancer Cell Int 21, 399 (2021). https://doi.org/10.1186/s12935-021-02101-9
- Cervical cancer
- CXC chemokine
- Quantitative real-time PCR
- Tumor microenvironment