
Bladder cancer (BLCA) has a high rate of morbidity and mortality, and is considered as one of the most malignant tumors of the urinary system. Tumor cells interact with surrounding interstitial cells, playing a key role in carcinogenesis and progression, which is partly mediated by chemokines. CXC chemokines exert anti-tumor biological roles in the tumor microenvironment and affect patient prognosis. Nevertheless, their expression and prognostic values patients with BLCA remain unclear.
We used online tools, including Oncomine, UALCAN, GEPIA, GEO databases, cBioPortal, GeneMANIA, DAVID 6.8, Metascape, TRUST (version 2.0), LinkedOmics, TCGA, and TIMER2.0 to perform the relevant analysis.
The mRNA levels of C-X-C motif chemokine ligand (CXCL)1, CXCL5, CXCL6, CXCL7, CXCL9, CXCL10, CXCL11, CXCL13, CXCL16, and CXCL17 were increased significantly increased, and those of CXCL2, CXCL3, and CXCL12 were decreased significantly in BLCA tissues as assessed using the Oncomine, TCGA, and GEO databases. GEO showed that high levels of CXCL1, CXCL6, CXCL10, CXCL11, and CXCL13 mRNA expression are associated significantly with the poor overall survival (all p < 0.05), and similarly, those of CXCL2 and CXCL12 in the TCGA database (p < 0.05). The predominant signaling pathways involving the differentially expressed CXC chemokines are cell cycle, chemokine, and cytokine-cytokine receptor interaction. Moreover, transcription factors such as Sp1 transcription factor (SP1), nuclear factor kappa B subunit 1 (NFKB1), and RELA proto-oncogene, NF-KB subunit (RELA) were likely play critical roles in regulating CXC chemokine expression. LYN proto-oncogene, src family tyrosine kinase (LYN) and LCK proto-oncogene, src family tyrosine kinase (LCK) were identified as the key targets of these CXC chemokines. MicroRNAs miR200 and miR30 were identified as the main microRNAs that interact with several CXC chemokines through an miRNA-target network. The expression of these chemokines is closely associated with the infiltration of six categories of immune cells.
We explored the CXC chemokines superfamily-based biomarkers associated with BLCA prognosis using public databases, and provided possible chemokine targets for patients with BLCA.
Citation: Xiaoqi Sun, Qunxi Chen, Lihong Zhang, Jiewei Chen, Xinke Zhang. Exploration of prognostic biomarkers and therapeutic targets in the microenvironment of bladder cancer based on CXC chemokines[J]. Mathematical Biosciences and Engineering, 2021, 18(5): 6262-6287. doi: 10.3934/mbe.2021313
[1] | Caiyun Wu, Cong Ma, Jing Yuan, Pei Zhou . Exploration of potential therapeutic and prognostic value of CXC chemokines in cervical squamous cell carcinoma and endocervical adenocarcinoma based on bioinformatics analysis. Mathematical Biosciences and Engineering, 2021, 18(6): 8201-8222. doi: 10.3934/mbe.2021407 |
[2] | Guobing Lin, Cunming Zhang, Xuanyu Chen, Junwei Wang, Song Chen, Siyuan Tang, Tianqiang Yu . Identification of circulating miRNAs as novel prognostic biomarkers for bladder cancer. Mathematical Biosciences and Engineering, 2020, 17(1): 834-844. doi: 10.3934/mbe.2020044 |
[3] | Zekun Xin, Yang Li, Lingyin Meng, Lijun Dong, Jing Ren, Jianlong Men . Elevated expression of the MYB proto-oncogene like 2 (MYBL2)-encoding gene as a prognostic and predictive biomarker in human cancers. Mathematical Biosciences and Engineering, 2022, 19(2): 1825-1842. doi: 10.3934/mbe.2022085 |
[4] | Zhongwei Zhao, Xiaoxi Fan, Lili Yang, Jingjing Song, Shiji Fang, Jianfei Tu, Minjiang Chen, Liyun Zheng, Fazong Wu, Dengke Zhang, Xihui Ying, Jiansong Ji . The identification of a common different gene expression signature in patients with colorectal cancer. Mathematical Biosciences and Engineering, 2019, 16(4): 2942-2958. doi: 10.3934/mbe.2019145 |
[5] | Yonghui Wang, Jianhui Huang, Zhifeng Tian, Yuefen Zhou, Jie Yang . The role of CXC cytokines as biomarkers and potential targets in hepatocellular carcinoma. Mathematical Biosciences and Engineering, 2020, 17(2): 1381-1395. doi: 10.3934/mbe.2020070 |
[6] | Yuan Yang, Lingshan Zhou, Xi Gou, Guozhi Wu, Ya Zheng, Min Liu, Zhaofeng Chen, Yuping Wang, Rui Ji, Qinghong Guo, Yongning Zhou . Comprehensive analysis to identify DNA damage response-related lncRNA pairs as a prognostic and therapeutic biomarker in gastric cancer. Mathematical Biosciences and Engineering, 2022, 19(1): 595-611. doi: 10.3934/mbe.2022026 |
[7] | Eugene Kashdan, Svetlana Bunimovich-Mendrazitsky . Hybrid discrete-continuous model of invasive bladder cancer. Mathematical Biosciences and Engineering, 2013, 10(3): 729-742. doi: 10.3934/mbe.2013.10.729 |
[8] | Svetlana Bunimovich-Mendrazitsky, Yakov Goltser . Use of quasi-normal form to examine stability of tumor-free equilibrium in a mathematical model of bcg treatment of bladder cancer. Mathematical Biosciences and Engineering, 2011, 8(2): 529-547. doi: 10.3934/mbe.2011.8.529 |
[9] | K. E. Starkov, Svetlana Bunimovich-Mendrazitsky . Dynamical properties and tumor clearance conditions for a nine-dimensional model of bladder cancer immunotherapy. Mathematical Biosciences and Engineering, 2016, 13(5): 1059-1075. doi: 10.3934/mbe.2016030 |
[10] | Xin Lin, Xingyuan Li, Binqiang Ma, Lihua Hang . Identification of novel immunomodulators in lung squamous cell carcinoma based on transcriptomic data. Mathematical Biosciences and Engineering, 2022, 19(2): 1843-1860. doi: 10.3934/mbe.2022086 |
Bladder cancer (BLCA) has a high rate of morbidity and mortality, and is considered as one of the most malignant tumors of the urinary system. Tumor cells interact with surrounding interstitial cells, playing a key role in carcinogenesis and progression, which is partly mediated by chemokines. CXC chemokines exert anti-tumor biological roles in the tumor microenvironment and affect patient prognosis. Nevertheless, their expression and prognostic values patients with BLCA remain unclear.
We used online tools, including Oncomine, UALCAN, GEPIA, GEO databases, cBioPortal, GeneMANIA, DAVID 6.8, Metascape, TRUST (version 2.0), LinkedOmics, TCGA, and TIMER2.0 to perform the relevant analysis.
The mRNA levels of C-X-C motif chemokine ligand (CXCL)1, CXCL5, CXCL6, CXCL7, CXCL9, CXCL10, CXCL11, CXCL13, CXCL16, and CXCL17 were increased significantly increased, and those of CXCL2, CXCL3, and CXCL12 were decreased significantly in BLCA tissues as assessed using the Oncomine, TCGA, and GEO databases. GEO showed that high levels of CXCL1, CXCL6, CXCL10, CXCL11, and CXCL13 mRNA expression are associated significantly with the poor overall survival (all p < 0.05), and similarly, those of CXCL2 and CXCL12 in the TCGA database (p < 0.05). The predominant signaling pathways involving the differentially expressed CXC chemokines are cell cycle, chemokine, and cytokine-cytokine receptor interaction. Moreover, transcription factors such as Sp1 transcription factor (SP1), nuclear factor kappa B subunit 1 (NFKB1), and RELA proto-oncogene, NF-KB subunit (RELA) were likely play critical roles in regulating CXC chemokine expression. LYN proto-oncogene, src family tyrosine kinase (LYN) and LCK proto-oncogene, src family tyrosine kinase (LCK) were identified as the key targets of these CXC chemokines. MicroRNAs miR200 and miR30 were identified as the main microRNAs that interact with several CXC chemokines through an miRNA-target network. The expression of these chemokines is closely associated with the infiltration of six categories of immune cells.
We explored the CXC chemokines superfamily-based biomarkers associated with BLCA prognosis using public databases, and provided possible chemokine targets for patients with BLCA.
Bladder cancer (BLCA) is characterized by a high rate of morbidity and mortality, and is diagnosed in approximately 430,000 people per annum, of whom around 165,000 die of this cancer each year [1]. More than 90% of BLCA is considered to be urinary urothelial cancer according to its pathological histological subtype and smoking as a risk factor caused its formation [2]. Although the results of a prior clinical trial showed 5- and 10-year survival rates of approximately 50 and 36%, respectively, for muscle-invasive BLCA using neoadjuvant chemotherapy [3], other studies showed 5-year survival rates ranging from 15 to 20% and high recurrence rates [4,5]. Therefore, for patients with high-risk non-metastatic BLCA, the recommended therapeutic strategy is radical cystectomy with lymph node dissection because of the poor efficacy of drug treatment, including maintenance bacille Calmette-Guérin (BCG) [6,7]. To date, tumor stage and lymph node (LN) status have been identified as the most important prognostic factors after radical cystectomy [8], whereas in clinical practice, the value of other prognostic risk factors, including neutrophil-to-lymphocyte ratio, remains controversial [9,10]. Recently, molecular subtypes of BLCA derived from The Cancer Genome Atlas (TCGA) database attracted our attention. However, limited clinical value have restricted their application [11]. To date, although promising predictive biomarkers have been identified, they are not used routinely in clinical practice because of their finite values, and decision-making concerning therapy cannot be based only on these molecular markers, for example the tumor mutation burden, ribonucleic acid subtypes, and neutrophil-to-lymphocyte ratio [12]. Therefore, it is urgent to identify more therapeutic targets and prognostic biomarkers for use in patients with BLCA.
CXC chemokines represent the main component of approximately 50 chemokine family members, and their C-termini harbor a heparin-binding domain, which plays vital roles in the regulation of tumor-associated angiogenesis, tumor growth, and metastatic potential [13,14], in addition to mediating the migration of diversified leukocytes in non-tumor/tumor microenvironment (TME) [15,16,17]. The key role of CXC chemokines has been established as linking tumor development and metastasis; therefore, CXC chemokines and their receptors have been identified as therapeutic target in several cancers. Mounting experimental evidence shows that monoclonal antibodies against CXC chemokines or their receptors could inhibit tumor growth and/or metastasis, including breast cancer [18], non-Hodgkin lymphoma [19] brain tumors [20], pancreatic tumors, and bladder tumors [21]. Therefore, CXC chemokines are likely to become therapeutic targets or prognostic factors in many tumors including BLCA.
Previous studies have revealed the functional role of several CXC chemokines in BLCA [22,23,24,25,26]. However, the identification of optimal prognostic biomarkers and therapeutic targets originating from CXC chemokines in BLCA remains an urgent and unresolved problem. Therefore, with the aim of fully comprehending the function of CXC chemokines in BLCA, the present study investigated the expression of CXC chemokines and their potential as prognostic biomarkers and therapeutic targets using exhaustive bioinformatic analyses based on several established databases. In addition, we provided evidence that will allow clinicians to choose appropriate therapeutic targets and accurately predict clinical outcome in patients with BLCA.
Oncomine (www.oncomine.org) integrates RNA and DNA-sequencing data from sources such as GEO, TCGA, and the published literature [27]. Screening of mRNA data allowed us to investigate CXC chemokine expression in BLCA. We set the significance thresholds of a p-value less than 1E-4, two-fold change, and overexpression or underexpression gene rank in top the 10%. The differential expression of CXC chemokine mRNAs in BLCA was analyzed using Student's t test.
UALCAN (http://ualcan.path.uab.edu/analysis.html) is an effective online analysis and mining website of tumor data, mainly based on the TCGA database, and can be used to perform biomarker identification, expression profile analysis, and survival analysis [28]. In this study, CXC chemokine mRNA expression data was obtained using the "Expression Analysis" module of the "BLCA" dataset. The differential expression of CXC chemokine mRNAs in BLCA was analyzed using Student's t test.
GEPIA (http://gepia.cancer-pku.cn/index.html) is an analysis website of RNA sequencing data based on TCGA and GTEx databases, which was developed by Tang et al. [29]. We investigated the relationship between pathological stage and CXC chemokines with differential mRNA expression and their prognostic significance by screening Oncomine and UALCAN databases using the "Single Gene Analysis" module of GEPIA in BLCA.
GEO (Gene Expression Omnibus) databases collect and organize various expression chip data, and other chips, such as methylation chips, long noncoding RNA (lncRNA) chips, miRNA chips, chromosome number variation (CNV) chips, and even high-throughput sequencing data. CXC chemokine mRNA differential expression in BLCA was analyzed using Student's t test and the corresponding prognostic analysis was performed using the statistical software GraphPad Prism 8.0.1 (GraphPad Inc., La Jolla, CA, USA).
The data in cBioPortal (www.cbioportal.org) is derived from the International Cancer Genome Consortium (ICGC), TCGA, GEO, and other databases, which includes DNA methylation data, limited clinical data mRNA and microRNA expression data, non-synonymous mutations, protein level and phosphoprotein level (reverse-phase protein array (RPPA)) data, and DNA copy number data [30]. Based on TGCA data, cBioPortal provided genetic alterations and co-expression data for CXC chemokines. Bladder Urothelial Carcinoma (TCGA, Firehose Legacy) data from 413 cases were analyzed. Protein expression z scores (RPPA) and mRNA expression z scores (RNA Seq V2 (RNA-Seq by expectation-maximization (RSEM)) were entered using a z score threshold of ± 2.0.
The GeneMANIA plug-in of Cytoscape software is used to study protein-protein interactions [31].
STRING (https://string-db.org/) is an important software used to analyze protein-protein interactions [32]. To investigate the interactions of CXC chemokines, we constructed a network analysis of protein-protein interactions using STRING.
DAVID 6.8 (https://david.ncifcrf.gov/home.jsp) is used mainly for functional and pathway enrichment analysis for differentially expressed genes [33]. Differentially expressed CXC chemokines were analyzed using Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis and Gene Ontology (GO) enrichment analysis. The top 50 mutated genes of Bladder Urothelial Carcinoma (TCGA, Firehose Legacy; p < 0.05) were analyzed using DAVID 6.8. The "ggplot2" package in the R project visualized the data.
Metascape (http://metascape.org) is a powerful gene function annotation analysis tool, which can perform annotation enrichment for large quantities of genes or proteins, and construct protein-protein interaction networks [34]. Our study utilized the "Express Analysis" module to validate the enriched CXC chemokines and the top 50 mutated genes of Bladder Urothelial Carcinoma (TCGA, Firehose Legacy).
The TRRUST (version 2) database (https://www.grnpedia.org/trrust/) records transcription factor regulatory relationships, including the target genes corresponding to transcription factors, and the regulatory relationships among transcription factors [35].
TIMER2.0 (http://timer.cistrome.org/) is a tumor immunity related database. "Gene module" refers to analyzing the association between CXC chemokine expression and immune cell types. "Clinical module" assesses the relationship among clinical outcome and tumor infiltrating immune cells and CXC chemokines [36].
LinkedOmics (http://www.linkedomics.org/) was used to analyze tumor multi-omics data for 32 TCGA tumor types [37]. "LinkInterpreter module" was used to obtain biological insights into transcription factor, miRNA, and kinase targets of enriched CXC chemokines. A minimum number of genes of 3 and a simulation of 500 were performed in Gene Set Enrichment Analysis (GSEA) for the BLCA dataset.
Heatmap of the correlations of differentially expressed CXC chemokines and bubble plot of enrichment analysis of the 50 top mutated genes and differentially expressed CXC chemokines in bladder cancer was conducted using R version 3.6.3. Visualization of gene-gene interaction network of the 50 top mutated genes and the differentially expressed CXC chemokines was performed with Cytoscape software. CXC chemokine mRNA differential expression in BLCA was analyzed using Student's t test and the corresponding survival analysis was performed using the statistical software GraphPad Prism 8.0.1 (GraphPad Inc., La Jolla, CA, USA). The plots of Figures 1 and 2 were achieved from online database UALCAN (http://ualcan.path.uab.edu/analysis.html).
ONCOMINE database analysis allowed us to analyze the mRNA expression levels of 16 CXC chemokines (CXCL1 to CXCL14, CXCL16 and CXCL17) in BLCA and normal bladder mucosa tissues. Table 1 and Figure 3 show the results. ONCOMINE data demonstrated the mRNA expression levels of CXCL13, CXCL10, CXCL9, and CXCL6 increased significantly and those of CXCL2, CXCL3, and CXCL12 decreased significantly in BLCA tissues compared with those in normal bladder mucosa tissues. This agreed with the results of Sanchez-Carbayo et al., who reported significantly increased mRNA levels of CXCL13, CXCL10, CXCL9, and CXCL6 in infiltrating bladder urothelial carcinoma and obviously decreased mRNA levels of CXCL2, CXCL3, and CXCL12 in superficial bladder cancer [38]. Lee et al. observed lower CXCL12 expression was in superficial bladder cancer (fold change = -4.495 and p = 2.52 e-21) and infiltrating bladder urothelial carcinoma (fold change = -3.371 and p = 1.41 e-12) compared with that in normal bladder mucosa [39].
TLR | Type | Fold change | P-value | t-test | References |
CXCL2 | Superficial Bladder Cancer | -12.716 | 7.74E-17 | -11.618 | [38] |
CXCL3 | Superficial Bladder Cancer | -3.564 | 2.44E-08 | -6.122 | [38] |
CXCL6 | Infiltrating Bladder Urothelial Carcinoma | 2.272 | 1.34E-06 | 5.02 | [38] |
CXCL9 | Infiltrating Bladder Urothelial Carcinoma | 2.84 | 3.92E-09 | 6.197 | [38] |
CXCL10 | Infiltrating Bladder Urothelial Carcinoma | 3.039 | 6.79E-11 | 7.04 | [38] |
CXCL12 | Superficial Bladder Cancer | -4.495 | 2.52E-21 | -12.211 | [39] |
Infiltrating Bladder Urothelial Carcinoma | -3.371 | 1.41E-12 | -7.727 | [39] | |
Superficial Bladder Cancer | -2.531 | 1.92E-12 | -8.846 | [38] | |
CXCL13 | Infiltrating Bladder Urothelial Carcinoma | 6.369 | 1.07E-10 | 7.171 | [38] |
Next, we used UALCAN to evaluate the mRNA expression levels of CXC chemokines in BLCA and normal bladder mucosa tissues, which showed markedly increased mRNA expression of CXCL1 (p = 8.50e-03), CXCL5 (p = 1.12e-06), CXCL7 (p = 1.42e-02), CXCL9 (p = 4.06e-08), CXCL10 (p = 3.91e-09), CXCL11 (p = 6.29e-08), CXCL16 (p = 1.44e-03), and CXCL17 (p = 3.56e-02) in BLCA tissues, and significantly decreased mRNA levels of CXCL2 (p = 2.30e-02) and CXCL12 (p = 2.82e-04) (Figure 1). We also revealed that CXC17 had the highest mRNA expression level in BLCA tissues based on the TCGA database (Figure 2). To investigate the role of the CXC chemokines in BLCA more comprehensively, including tumorigenesis and clinical prognosis, we assessed all the differentially expressed CXC chemokines that were acquired by the Oncomine and TCGA databases including CXCL1, CXCL2, CXCL3, CXCL5, CXCL6, CXCL7, CXCL9, CXCL10, CXCL11, CXCL12, CXCL13, CXCL16, and CXCL17.
Next, we evaluated the relationship between the pathological stage of patients with BLCA and the differentially expressed CXC chemokines. We observed a significant association between CXCL2 (p = 0.0127) and CXCL12 (p = 8.68e-07) expression and pathological stage (Figure 4). CXCL2 and CXCL12 mRNA expression levels were higher in more advanced pathological stages, indicating that they might associated with BLCA tumorigenesis. Meanwhile, GEPIA was used to assess the impact of differentially expressed CXC chemokines on clinical outcome. Patients with BLCA with high CXCL2 (p = 0.17), CXCL16 (p = 0.14), and CXCL13 (p = 0.12) mRNA expression had a statistical tendency towards longer disease-free survival (Figure 5). In addition, high mRNA levels of CXCL1 (p = 0.17) and CXCL12 (p = 0.0075) had a statistical tendency towards shorter overall survival (Figure 6).
GEO database analysis demonstrated that CXCL3 mRNA expression (p = 0.034) was higher in BLCA tissues than in normal tissues. CXCL12 (p < 0.001), CXCL9 (p = 0.002), and CXCL10 (p = 0.015) mRNA expression levels were higher in normal tissues than in BLCA tissues (Figure 7). Survival analysis showed that high levels of CXCL1, CXCL6, CXCL10, CXCL11, and CXCL13 mRNA expression were associated significantly with the poor overall survival (all p < 0.05). High levels of CXCL3 mRNA had the statistic tendency towards predicting shorter overall survival (p = 0.053) (Figure 8).
We analyzed comprehensively the molecular characteristics of the differentially expressed CXC chemokines. First, we used TCGA datasets to analyze their genetic alterations, which showed the rates of genetic alteration of CXCL1, CXCL2, CXCL3, CXCL5, CXCL6, CXCL7, CXCL9, CXCL10, CXCL11, CXCL12, CXCL13, CXCL16, and CXCL17 were 6, 5, 4, 5, 5, 4, 5, 6, 5, 2.7, 3, 5, and 11%, respectively in the BLCA samples queried (Figure 9A). In the BLCA samples, the most common alterations were amplification and high mRNA expression. Differentially expressed CXC chemokine co-expression analysis identified strong correlations among CXCL9, CXCL10, and CXCL11 expression; moderate to strong correlations among CXCL1, CXCL2, CXCL3, CXCL5, and CXCL6; and weak to moderate correlations among CXCL7, CXCL12, CXCL13, CXCL16, and CXCL17 expression (Figure 9B). Secondly, PPI network analysis with STRING obtained 13 nodes and 71 edges in the analysis of the interactions among the differentially expressed CXC chemokines (Figure 9D). Enrichment pathways identified in STRING for these CXC chemokines included the TNF signaling pathway, the IL-17 signaling pathway, cytokine-cytokine receptor interactions, and the chemokine signaling pathway. The results of analysis at GeneMANIA revealed that their functions were related to chemokine receptor binding, cell chemotaxis, and chemokine activity (Figure 9E).
In addition, the top 50 mutated genes of 200 most mutated cases with BLCA were identified using cBioPortal (Figure 9C), representing mutations in which lead to the formation of a tumor-associated antigen that is associated with the TME and is closely associated with the CXC chemokines. Next, we analyzed the associations of these mutated genes and the CXC chemokines. These top 50 mutated genes have close interactions with the differentially expressed CXC chemokines in BLCA (Figure 9F).
The functions of differentially expressed CXC chemokines and the top 50 mutated genes in BLCA were analyzed by DAVID6.8 and Metascape. The top 10 GO enriched terms in the biological process (BP), cellular component (CC), and molecular function (MF) categories were analyzed using DAVID 6.8. Chemokine-mediated signaling pathway, positive regulation of leukocyte chemotaxis, response to lipopolysaccharide, regulation of cell proliferation, inflammatory response, immune response, positive regulation of cAMP metabolic process, and cell-cell signaling were related to the development and progression of BLCA (Figure 10A). The top 10 GO terms in the CC category included, Z disc, sarcolemma, extracellular space, sarcoplasmic reticulum membrane, sarcoplasmic reticulum, extracellular region, T-tubule, MLL3/4 complex, costamere, and histone methyltransferase complex (Figure 10B). Chemokine activity and CXCR chemokine receptor binding were the two most significantly enriched GO terms in the MF category (Figure 10C). The top 10 KEGG pathways were the chemokine signaling pathway, cytokine-cytokine receptor interaction, cell cycle, the TNF signaling pathway, and microRNAs in cancer, HIF-1 signaling pathway, pathways in cancer, viral carcinogenesis, central carbon metabolism in cancer, which were intimately associated with BLCA tumorigenesis (Figure 10D).
Metascape analysis showed that the functions of the differentially expressed CXC chemokines and top 50 mutated genes in BLCA were predominantly enriched in chemokine receptor bind chemokines, T cell chemotaxis, and microRNAs in cancer (Figure 11A, B). Next, we extracted protein-protein interaction (PPI) and MCODE components to explore the correlations between BLCA and CXC chemokines in detail; the list of identified genes is shown in Figure 11C, D. Three most significant MCODE components were obtained and showed that the biological functions were involved in chemokine receptor binding to chemokines, chemokine-mediated signaling pathways, chemokine signaling pathways, activation of homeobox (HOX) genes during differentiation, and histone modification (Figure 11E).
Next, we investigated the transcription factor and kinase targets of the differentially expressed chemokines by analyses in the TRRUST and LinkedOmics databases. TRRUST analysis of CXCL1, CXCL2, CXCL5, CXCL7, CXCL10, and CXCL12 identified three transcription factors: Sp1 transcription factor (SP1), nuclear factor kappa B subunit 1 (NFKB1), and RELA proto-oncogene, NF-KB subunit (RELA) that were possibly related CXC chemokine regulation (Table 2). Critically, NFKB1 and RELA mediated the expression of CXCL1, CXCL2, CXCL5, CXCL10, and CXCL12, and SP1 mediated the expression of CXCL1 and CXCL5. Next, LinkedOmics database analysis was used to identify the top two kinase targets of these CXC chemokines (Table 3). LYN proto-oncogene, src family tyrosine kinase (LYN) was identified as the kinase target of CXCL1, CXCL2, CXCL3, CXCL5, CXCL6, CXCL9, CXCL10, CXCL11, CXCL12, and CXCL13 in the kinase-target network, and LCK proto-oncogene, src family tyrosine kinase (LCK) was identified as the kinase target of CXCL2, CXCL3, CXCL5, CXCL6, CXCL9, CXCL10, CXCL11, CXCL13, and CXCL16. Inhibitor of nuclear factor kappa B kinase subunit beta (IKBKB), protein kinase CGMP-dependent 1 (PRKG1), and Janus kinase 2 (JAK2) were identified as key kinase targets of CXCL1, CXCL12, and CXCL16, respectively, in the kinase-target network. Moreover, protein kinase N1 (PKN1), adrenergic, beta, receptor kinase 1 (ADRBK1) and checkpoint kinase 1 (CHEK1), polo like kinase 1 (PLK1) were identified as the top two target kinase for CXCL7 and CXCL17, respectively. Finally, analysis at the LinkedOmics database allowed us to identify miRNAs that targeted the mRNAs encoding CXC chemokines (Table 4). The miRNA miR200A interacted mainly with the CXCL1 and CXCL5 miRNA-target network. Chemokine mRNAs CXCL1 and CXCL6 might be targeted by miR-30E-5P, miR-30B, miR-30D, miR-30C, and miR-30A-5P. CXCL3 expression might be modulated by miR-409-3P and miR-526B. CXCL5 and CXCL6 might be targeted by miR-141 and miR-374, respectively. CXCL7 might be targeted by miR-154, miR-487, miR-525, and miR-524s. Finally, CXCL17 might be targeted by miR-490.
Key TF | Description | Overlapped genes | P value | Q value | List of overlapped genes |
RELA | v-rel reticuloendotheliosis viral oncogene homolog A (avian) | 5 | 1.16E-06 | 1.79E-06 | CXCL10, CXCL12, CXCL2, CXCL5, CXCL1 |
NFKB1 | nuclear factor of kappa light polypeptide gene enhancer in B-cells 1 | 6 | 1.16E-06 | 1.79E-06 | CXCL12, CXCL10, CXCL5, CXCL1, CXCL2 |
SP1 | Sp1 transcription factor | 2 | 0.0406 | 0.0406 | CXCL5, CXCL1 |
CXC chemokines | Kinase target | Description | Leading Edge Number | P Value |
CXCL1 | Kinase_LYN | LYN proto-oncogene, Src family tyrosine kinase | 17 | 0 |
Kinase_IKBKB | inhibitor of nuclear factor kappa B kinase subunit beta | 9 | 0 | |
CXCL2 | Kinase_LCK | LCK proto-oncogene, Src family tyrosine kinase | 24 | 0 |
Kinase_LYN | LYN proto-oncogene, Src family tyrosine kinase | 21 | 0 | |
CXCL3 | Kinase_LYN | LYN proto-oncogene, Src family tyrosine kinase | 19 | 0 |
Kinase_LCK | LCK proto-oncogene, Src family tyrosine kinase | 25 | 0 | |
CXCL5 | Kinase_LCK | LCK proto-oncogene, Src family tyrosine kinase | 24 | 0 |
Kinase_LYN | LYN proto-oncogene, Src family tyrosine kinase | 18 | 0 | |
CXCL6 | Kinase_SYK | spleen associated tyrosine kinase | 16 | 0 |
Kinase_LYN | LYN proto-oncogene, Src family tyrosine kinase | 24 | 0 | |
CXCL7 | Kinase_PKN1 | protein kinase N1 | 5 | 0.009 |
Kinase_ADRBK1 | G protein-coupled receptor kinase 2 | 10 | 0 | |
CXCL9 | Kinase_LYN | LYN proto-oncogene, Src family tyrosine kinase | 18 | 0 |
Kinase_LCK | LCK proto-oncogene, Src family tyrosine kinase | 20 | 0 | |
CXCL10 | Kinase_LCK | LCK proto-oncogene, Src family tyrosine kinase | 20 | 0 |
Kinase_LYN | LYN proto-oncogene, Src family tyrosine kinase | 16 | 0 | |
CXCL11 | Kinase_LCK | LCK proto-oncogene, Src family tyrosine kinase | 20 | 0 |
Kinase_LYN | LYN proto-oncogene, Src family tyrosine kinase | 17 | 0 | |
CXCL12 | Kinase_LYN | LYN proto-oncogene, Src family tyrosine kinase | 20 | 0 |
Kinase_PRKG1 | protein kinase, cGMP-dependent, type I | 15 | 0.004 | |
CXCL13 | Kinase_LCK | LCK proto-oncogene, Src family tyrosine kinase | 22 | 0 |
Kinase_LYN | LYN proto-oncogene, Src family tyrosine kinase | 19 | 0 | |
CXCL16 | Kinase_LCK | LCK proto-oncogene, Src family tyrosine kinase | 22 | 0 |
Kinase_JAK2 | Janus kinase 2 | 7 | 0 | |
CXCL17 | Kinase_CHEK1 | checkpoint kinase 1 | 45 | 0 |
Kinase_PLK1 | polo like kinase 1 | 41 | 0 |
CXC chemokines | MIRNA targets | Leading Edge Number | P Value |
CXCL1 | GTAAGAT, MIR-200A | 20 | 0.009 |
TGTTTAC, MIR-30A-5P, MIR-30C, MIR-30D, MIR-30B, MIR-30E-5P | 133 | 0 | |
CXCL3 | AACATTC, MIR-409-3P | 47 | 0.009 |
CXCL5 | CTCAAGA, MIR-526B | 24 | 0.019 |
CAGTGTT, MIR-141, MIR-200A | 72 | 0.015 | |
CXCL6 | TATTATA, MIR-374 | 89 | 0 |
TGTTTAC, MIR-30A-5P, MIR-30C, MIR-30D, MIR-30B, MIR-30E-5P | 192 | 0 | |
CXCL7 | GTATGAT, MIR-154, MIR-487 | 21 | 0 |
GCGCCTT, MIR-525, MIR-524 | 7 | 0.044 | |
CXCL17 | CCAGGTT, MIR-490 | 10 | 0.025 |
CXC chemokines influence immune cell infiltration and inflammatory responses in the TME, thus possibly affecting the clinical therapy and outcome in patients with BLCA. Using the TIMER database, we investigated comprehensively the correlation between infiltration of various immune cells and the differentially expressed CXC chemokines. CXCL1 expression correlated positively with the infiltration of dendritic cells (Cor = 0.427, p = 1.43e-17, neutrophils (Cor = 0.444, p = 5.46e-19), CD4+ T cells (Cor = 0.112, p = 3.26e-02), and CD8+ T cells (Cor = 0.303, p = 3.39e-09); Figure 12A). Similarly, CXCL2 expression correlated positively with the infiltration of dendritic cells (Cor = 0.345, p = 1.22e-11, neutrophils (Cor = 0.355, p = 3.05e-12), CD4+ T cells (Cor = 0.128, p = 1.46e-02), and CD8+ T cells (Cor = 0.255, p = 7.61e-07; Figure 12B). CXCL3 expression correlated negatively with the infiltration of dendritic cells (Cor = 0.4, p = 1.93e-15, neutrophils (Cor = 0.424, p = 2.86e-17), CD4+ T cells (Cor = 0.108, p = 3.89e-02), CD8+ T cells (Cor = 0.302, p = 3.55e-09), and B cells (Cor = -0.103, p = 4.88e-02); Figure 12C). CXCL5 expression correlated positively with the infiltration of dendritic cells (Cor = 0.43, p = 7.83e-18), neutrophils (Cor = 0.385, p = 2.70e-14), CD4+ T cells (Cor = 0.116, p = 2.63e-02), CD8+ T cells (Cor = 0.288, p = 2.10e-08; Figure 12D). CXCL6 expression correlated positively with the infiltration of CD8+ T cells (Cor = 0.219, p = 2.42e-05), neutrophils (Cor = 0.339, p = 3.50e-11), and dendritic cells (Cor = 0.43, p = 7.83e-18; Figure 12E). A similar correlation tendency was found between CXCL7 expression and dendritic cells (Cor = 0.239, p = 3.74e-06), neutrophils (Cor = 0.174, p = 8.87e-04), and CD8+ T cells (Cor = 0.164, p = 1.60e-03; Figure 12F). CXCL9, CXCL10, and CXCL11 mRNA expression levels correlated positively with the infiltration of dendritic cells, neutrophils, CD4+ T cells, and CD8+ T cells (all p < 0.05; Figure 12G-I). CXCL12 expression correlated positively with the infiltration of dendritic cells (Cor = 0.156, p = 2.88e-03), neutrophils (Cor = 0.208, p = 6.27e-05), macrophages (Cor = 0.464, p = 7.25e-21), CD4+ T cells (Cor = 0.146, p = 5.21e-03), and CD8+ T cells (Cor = 0.134, p = 1.05e-02; Figure 12J). With the exception of macrophages, CXCL13 expression correlated positively associated with the other immune cell types consisted of dendritic cells, CD4+ T cells, neutrophils, CD8+ T cells, and B cells (p < 0.05; Figure 12K). CXCL16 expression correlated positively with the infiltration of CD8+ T cells, CD4+ T cells, macrophages, neutrophils, and dendritic cells (all p < 0.05; Figure 12L). CXCL17 mRNA expression correlated positively with B cells (Cor = 0.156, p = 2.83e-03), and negatively with dendritic cells (Cor = -0.167, p = 1.35e-03; Figure 12M). A Cox proportional hazard model was used identify factors associated with prognosis of patients with BLCA, which identified B cells (p = 0.032), macrophages (p = 0.001), and CXCL2 mRNA expression (p = 0.010). CXCL7 mRNA expression (p = 0.066) had significantly statistical tendency towards predicting the clinical outcome of BLCA patients (Table 5).
coef | HR | 95%CI_l | 95%CI_u | p value | sig | |
B_cell | -3.430 | 0.032 | 0.001 | 0.739 | 0.032 | * |
CD8_Tcell | 1.675 | 5.337 | 0.294 | 96.758 | 0.257 | |
CD4_Tcell | -0.733 | 0.480 | 0.008 | 29.113 | 0.726 | |
Macrophage | 4.651 | 104.684 | 6.900 | 1588.200 | 0.001 | ** |
Neutrophil | -3.390 | 0.034 | 0.000 | 17.688 | 0.289 | |
Dendritic | -0.073 | 0.929 | 0.171 | 5.054 | 0.933 | |
CXCL1 | 0.101 | 1.106 | 0.968 | 1.265 | 0.138 | |
CXCL2 | -0.279 | 0.756 | 0.612 | 0.934 | 0.010 | * |
CXCL3 | 0.192 | 1.212 | 0.929 | 1.581 | 0.157 | |
CXCL5 | 0.034 | 1.035 | 0.895 | 1.196 | 0.645 | |
CXCL6 | -0.083 | 0.921 | 0.801 | 1.059 | 0.247 | |
CXCL7 | 0.138 | 1.148 | 0.991 | 1.329 | 0.066 | |
CXCL9 | -0.028 | 0.973 | 0.808 | 1.171 | 0.770 | |
CXCL10 | 0.048 | 1.050 | 0.843 | 1.306 | 0.665 | |
CXCL11 | -0.060 | 0.942 | 0.762 | 1.166 | 0.584 | |
CXCL12 | 0.089 | 1.093 | 0.961 | 1.242 | 0.176 | |
CXCL13 | 0.040 | 1.041 | 0.919 | 1.178 | 0.531 | |
CXCL16 | -0.091 | 0.913 | 0.753 | 1.107 | 0.353 | |
CXCL17 | -0.010 | 0.990 | 0.926 | 1.057 | 0.756 |
Preclinical data indicate that BLCA is one of the most immunogenic tumors [40], which is assumed to make it sensitive to favorable targeted therapies, not including PD-1/PDL-1 checkpoint inhibitors. Previous studies have unveiled a tumor microenvironment (TME)-related role of CXC chemokines [41], and functions in the immunotherapy of cancer [42]. Chemokines are involved in the process of cancer biology, including tumor proliferation and metastasis, and targeting of chemokines and chemokine receptors have been evaluated in several preclinical studies and clinical trials, such as in breast cancer [43], pancreatic cancer [44], renal tumor [45], and prostate cancer [46]. Nonetheless, in the prognosis and biology of BLCA, the significance of CXC chemokines are poorly understood.
The present study characterized the mRNA expression patterns CXC chemokines and investigated their pathological relationship with BLCA stages. Oncomine and TCGA database analysis identified 13 CXC chemokines that were differentially expressed in BLCA tissues compared with that in normal tissues, including upregulated expression (CXCL1, CXCL5, CXCL6, CXCL7, CXCL9, CXCL10, CXCL11, CXCL13, CXCL16, and CXCL17) and downregulated expression (CXCL2, CXCL3, and CXCL12). In addition, the TCGA analysis results showed that high expression of CXCL2 and CXCL12 was associated with advanced pathological stage, and analysis of the GEO and TCGA databases revealed that high levels of CXCL1, CXCL6, CXCL10, CXCL11, CXCL12, and CXCL13 mRNA expression were associated significantly with poor overall survival, indicating these differentially expressed CXC chemokines exert important roles in the proliferation and development of BLCA. This agreed partially with the results of previous studies, which reported that CXCL5 promotes chemoresistance in BLCA cells [47] and CXCL12 is associated with the genesis and development of BLCA [23].
We next investigated the differentially expressed CXC chemokines for their molecular characteristics in BLCA. Their most frequent genetic alterations were high mRNA expression and amplification, following by deep deletion, missense mutation, and truncating mutation, suggesting that genetic alterations might play a vital role in tumorigenesis and progression of BLCA. In particular, high mRNA expression and amplification of CXC chemokine genes could have the predominant impact on the BLCA processes. Previous studies reported that overexpression of several chemokines, including CXCL1, CXCL2, CXCL3, CXCL4, CXCL5, CXCL6, CXCL7, and CXCL8, was involved in breast cancer and melanoma metastasis, resulting from enhanced gene transcription [48] and copy number amplification of their encoding genes [49]. Moreover, mutants of CXCL8 had lower potency and efficacy in chemotaxis assays using neutrophils [50]. Together, our results and those of previous studies suggest the critical role of CXC chemokines in the tumorigenesis and cancer progression.
We next used GO and KEGG analysis to provide clues to the biological functions of these CXC chemokines. The results showed the CXC chemokines and top 50 mutated genes in BLCA are functionally associated with the TNF pathway, cell cycle, cytokine-cytokine receptor interactions, and the chemokine signaling pathway. Tumor metastasis, immune evasion, angiogenesis, and proliferation are reported to be associated critically with chemokine signaling pathways [51,52]. The TNF signaling pathway responds to cellular stress and inflammatory signals to activate cytokine cascades, which further affect the immune microenvironment. In addition, TNF interacts with different receptors to mediate cellular apoptosis or growth [53,54]. Moreover, the TNF signaling pathway contributes to tumor growth and metastasis in various cancers and acts as a therapeutic target [55,56,57,58]. These pieces of evidence showed that CXC chemokines might be developed as effective treatment targets for patients with BLCA in the future.
We explored the transcription factor targets of the differentially expressed CXC chemokines, their kinase targets, and miRNAs that targeted them. We observed that SP1, NFKB1, and RELA might regulate the expression of CXC chemokines. RELA regulates breast cancer progression via the CXCL13/CXCR5 axis [59]. In addition, via the IRF1-CDK4 signaling axis, RELA-induced Interferon-γ inhibits the proliferation of breast cancer [60]. RELA also regulates oncogene-induced senescence to decrease the proliferation in murine Kras-driven pancreatic cancer via the CXCL1/CXCR2 axis [61]. Additionally, RELA acts as an important transcription factor in the activation and stability of regulatory T cells [62] and regulates Th17 differentiation by miR-30a in experimental autoimmune encephalomyelitis [63]. NFKB1 is involved widely in carcinogenesis, including driving tumor progression or acting as a tumor-suppressor by phosphorylation and ubiquitination pathways [64]. Polymorphisms of the NFKB1 promoter are associated with a higher risk of bladder cancer [65], and aberrant NFKB1 activity leads to impaired M1 polarization of macrophages by repressing the expression of CXCL10 and IL-12, which forms an immunosuppressive microenvironment to induce the immune escape of tumors [64]. Our study provided evidence that an association is showed between CXC chemokines and the NF-κB signaling pathway during BLCA tumorigenesis and development. Our results also suggested Src family tyrosine kinases (LCK and LYN), PKN1, PRKG1, and CHEK1 might be the targets of the differential CXC chemokines, which play important roles in cell growth, division, migration, and survival signaling pathways [66,67,68,69], and mediate the carcinogenesis of several tumors [70,71,72]. Therefore, these differentially expressed CXC chemokines might modulate BLCA development and progression dependent on the above-mentioned signaling pathways by regulating these kinases.
We found the miR-30 family is the main miRNA family targeting CXC chemokines such as CXCL1 and CXCL6 in BLCA, suggesting that this miRNA family plays crucial regulatory roles in BLCA carcinogenesis and progression by acting on mRNAs encoding CXC chemokines. This hypothesis was supported partially by previous reports that the miR-30 family is associated with the development of organs and the pathogenesis of diseases, including various cancers, such as breast cancer and colorectal cancer [73,74,75].
Chemokines function by mediating the migration of immune cells by interacting with their cognate ligands [41]. Immune cell infiltration, especially CD8+ T cells, could affect tumor progression, as well as predicting the sensitivity to immunotherapy and clinical prognosis [76,77]. The present study indicated that a distinct correlation is found between infiltrating immune cells (e.g., dendritic cells, neutrophils, macrophages, CD4+ T cells, CD8+ T cells, and B cells) and the differentially expressed CXC chemokines, suggesting that these chemokines acting as crucial components reflected the immune status in the TME and could provide evidence of drug targets in BLCA immunotherapy in the future.
We analyzed the relationship between the CXC chemokines and the 50 top mutated genes in patients with BLCA and explored the CXC chemokines superfamily-based biomarkers associated with BLCA prognosis by public databases. The CXC chemokines superfamily-based biomarkers gene signature has prognostic value in BLCA patients and thus provides and guides the treatment of BLCA.
Xinke Zhang and Jiewei Chen designed this study. Xiaoqi Sun and Qunxi Chen analyzed the data and wrote the manuscript. Lihong Zhang analyzed the data. Xinke Zhang revised the manuscript. All authors approved the final version for submission.
This study was supported by the Medical Science Funding of Guangdong province (A2020139).
The authors declare that they have no conflict of interest.
[1] |
S. Antoni, J. Ferlay, I. Soerjomataram, A. Znaor, A. Jemal, F. Bray, Bladder cancer incidence and mortality: a global overview and recent trends, Eur. Urol., 71 (2017), 96-108. doi: 10.1016/j.eururo.2016.06.010
![]() |
[2] | J. Ferlay, I. Soerjomataram, R. Dikshit, S. Eser, C. Mathers, M. Rebelo, et al., Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012, Int. J. Cancer, 136 (2015), E359-E386. |
[3] | International Collaboration of Trialists, International phase Ⅲ trial assessing neoadjuvant cisplatin, methotrexate, and vinblastine chemotherapy for muscle-invasive bladder cancer: long-term results of the BA06 30894 trial, J. Clin. Oncol., 29 (2011), 2171-2177. |
[4] | Bladder cancer: diagnosis and management of bladder cancer: (c) NICE (2015) Bladder cancer: diagnosis and management of bladder cancer, BJU Int., 120 (2017), 755-765. |
[5] | D. H. Aggen, C. G. Drake, Biomarkers for immunotherapy in bladder cancer: a moving target, J. Immunother. Cancer, 5 (2017), 94. |
[6] | M. Babjuk, M. Burger, E. M. Comperat, P. Gontero, A. H. Mostafid, J. Palou, et al., European association of urology guidelines on non-muscle-invasive bladder cancer (TaT1 and carcinoma in situ) - 2019 update, Eur. Urol., 76 (2019), 639-657. |
[7] | S. Cambier, R. J. Sylvester, L. Collette, P. Gontero, M. A. Brausi, G. van Andel, et al., EORTC nomograms and risk groups for predicting recurrence, progression, and disease-specific and overall survival in non-muscle-invasive stage Ta-T1 urothelial bladder cancer patients treated with 1-3 years of maintenance bacillus calmette-guerin, Eur. Urol., 69 (2016), 60-69. |
[8] | R. Dutta, A. Abdelhalim, J. W. Martin, S. L. Vernez, B. Faltas, Y. Lotan, et al., Effect of tumor location on survival in urinary bladder adenocarcinoma: a population-based analysis, Urol. Oncol., 34 (2016), 531.e1-531.e6. |
[9] | S. Wu, X. Zhao, Y. Wang, Z. Zhong, L. Zhang, J. Cao, et al., Pretreatment neutrophil-lymphocyte ratio as a predictor in bladder cancer and metastatic or unresectable urothelial carcinoma patients: a pooled analysis of comparative studies, Cell. Physiol. Biochem., 46 (2018), 1352-1364 |
[10] | E. Ojerholm, A. Smith, W. T. Hwang, B. C. Baumann, K. N. Tucker, S. P. Lerner, et al., Neutrophil-to-lymphocyte ratio as a bladder cancer biomarker: Assessing prognostic and predictive value in SWOG 8710, Cancer, 123 (2017), 794-801. |
[11] | A. Kamoun, A. de Reynies, Y. Allory, G. Sjodahl, A. G. Robertson, R. Seiler, et al., A consensus molecular classification of muscle-invasive bladder cancer, Eur. Urol., 77 (2020), 420-433. |
[12] | J. A. Witjes, H. M. Bruins, R. Cathomas, E. M. Comperat, N. C. Cowan, G. Gakis, et al., European association of urology guidelines on muscle-invasive and metastatic bladder cancer: summary of the 2020 guidelines, Eur. Urol., 79 (2021) 82-104, |
[13] | P. A. Gerber, A. Hippe, B. A. Buhren, A. Muller, B. Homey, Chemokines in tumor-associated angiogenesis, Biol. Chem., 390 (2009), 1213-1223. |
[14] |
E. C. Keeley, B. Mehrad, R. M. Strieter, CXC chemokines in cancer angiogenesis and metastases, Adv. Cancer Res., 106 (2010), 91-111. doi: 10.1016/S0065-230X(10)06003-3
![]() |
[15] |
A. D. Luster, Chemokines--chemotactic cytokines that mediate inflammation, N. Engl. J. Med., 338 (1998), 436-445. doi: 10.1056/NEJM199802123380706
![]() |
[16] |
S. Struyf, P. Proost, J. Van Damme, Regulation of the immune response by the interaction of chemokines and proteases, Adv. Immunol., 81 (2003), 1-44. doi: 10.1016/S0065-2776(03)81001-5
![]() |
[17] |
E. C. Keeley, B. Mehrad, R. M. Strieter, Chemokines as mediators of neovascularization, Arterioscler. Thromb. Vasc. Biol., 28 (2008), 1928-1936. doi: 10.1161/ATVBAHA.108.162925
![]() |
[18] | A. Muller, B. Homey, H. Soto, N. Ge, D. Catron, M. E. Buchanan, et al., Involvement of chemokine receptors in breast cancer metastasis, Nature, 410 (2001), 50-56. |
[19] | F. Bertolini, C. Dell'Agnola, P. Mancuso, C. Rabascio, A. Burlini, S. Monestiroli, et al., CXCR4 neutralization, a novel therapeutic approach for non-Hodgkin's lymphoma, Cancer Res., 62 (2002), 3106-3112. |
[20] | J. B. Rubin, A. L. Kung, R. S. Klein, J. A. Chan, Y. Sun, K. Schmidt, et al., A small-molecule antagonist of CXCR4 inhibits intracranial growth of primary brain tumors, Proc. Natl. Acad. Sci. U S A, 100 (2003), 13513-13518. |
[21] |
J. Vandercappellen, J. Van Damme, S. Struyf, The role of CXC chemokines and their receptors in cancer, Cancer Lett., 267 (2008), 226-244. doi: 10.1016/j.canlet.2008.04.050
![]() |
[22] | M. Miyake, S. Hori, Y. Morizawa, Y. Tatsumi, Y. Nakai, S. Anai, et al., CXCL1-mediated interaction of cancer cells with tumor-associated macrophages and cancer-associated fibroblasts promotes tumor progression in human bladder cancer, Neoplasia, 18 (2016), 636-646. |
[23] |
A. Nazari, H. Khorramdelazad, G. Hassanshahi, Biological/pathological functions of the CXCL12/CXCR4/CXCR7 axes in the pathogenesis of bladder cancer, Int. J. Clin. Oncol., 22 (2017), 991-1000. doi: 10.1007/s10147-017-1187-x
![]() |
[24] | H. Zhang, Y. L. Ye, M. X. Li, S. B. Ye, W. R. Huang, T. T. Cai, et al., CXCL2/MIF-CXCR2 signaling promotes the recruitment of myeloid-derived suppressor cells and is correlated with prognosis in bladder cancer, Oncogene, 36 (2017), 2095-2104. |
[25] | A. Nazari, Z. Ahmadi, G. Hassanshahi, M. Abbasifard, Z. Taghipour, S. K. Falahati-Pour, et al., Effective treatments for bladder cancer affecting CXCL9/CXCL10/CXCL11/CXCR3 axis: a review, Oman Med. J., 35 (2020), e103. |
[26] |
J. Zheng, X. Zhu, J. Zhang, CXCL5 knockdown expression inhibits human bladder cancer T24 cells proliferation and migration, Biochem. Biophys. Res. Commun., 446 (2014), 18-24. doi: 10.1016/j.bbrc.2014.01.172
![]() |
[27] | D. R. Rhodes, J. Yu, K. Shanker, N. Deshpande, R. Varambally, D. Ghosh, et al., ONCOMINE: a cancer microarray database and integrated data-mining platform, Neoplasia, 6 (2004), 1-6. |
[28] | D. S. Chandrashekar, B. Bashel, S. A. H. Balasubramanya, C. J. Creighton, I. Ponce-Rodriguez, B. Chakravarthi, et al., UALCAN: a portal for facilitating tumor subgroup gene expression and survival analyses, Neoplasia, 19 (2017), 649-658. |
[29] | Z. Tang, C. Li, B. Kang, G. Gao, C. Li, Z. Zhang, GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses, Nucleic. Acids. Res., 45 (2017), W98-W102. |
[30] | J. Gao, B. A. Aksoy, U. Dogrusoz, G. Dresdner, B. Gross, S. O. Sumer, et al., Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal, Sci. Signal., 6 (2013), pl1. |
[31] | D. Warde-Farley, S. L. Donaldson, O. Comes, K. Zuberi, R. Badrawi, P. Chao, et al., The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function, Nucleic Acids Res., 38 (2010), W214-220. |
[32] | D. Szklarczyk, A. L. Gable, D. Lyon, A. Junge, S. Wyder, J. Huerta-Cepas, et al., STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets, Nucleic Acids Res., 47 (2019), D607-D613. |
[33] |
W. Huang da, B. T. Sherman, R. A. Lempicki, Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources, Nat. Protoc., 4 (2009), 44-57. doi: 10.1038/nprot.2008.211
![]() |
[34] | Y. Zhou, B. Zhou, L. Pache, M. Chang, A. H. Khodabakhshi, O. Tanaseichuk, et al., Metascape provides a biologist-oriented resource for the analysis of systems-level datasets, Nat. Commun., 10 (2019), 1-10. |
[35] | H. Han, J. W. Cho, S. Lee, A. Yun, H. Kim, D. Bae, et al., TRRUST v2: an expanded reference database of human and mouse transcriptional regulatory interactions, Nucleic Acids Res., 46 (2018), D380-D386. |
[36] | T. Li, J. Fu, Z. Zeng, D. Cohen, J. Li, Q. Chen, et al., TIMER2.0 for analysis of tumor-infiltrating immune cells, Nucleic Acids Res., 48 (2020), W509-W514. |
[37] | S. V. Vasaikar, P. Straub, J. Wang, B. Zhang, LinkedOmics: analyzing multi-omics data within and across 32 cancer types, Nucleic Acids Res., 46 (2018), D956-D963. |
[38] |
M. Sanchez-Carbayo, N. D. Socci, J. Lozano, F. Saint, C. Cordon-Cardo, Defining molecular profiles of poor outcome in patients with invasive bladder cancer using oligonucleotide microarrays, J. Clin. Oncol., 24 (2006), 778-789. doi: 10.1200/JCO.2005.03.2375
![]() |
[39] | J. S. Lee, S. H. Leem, S. Y. Lee, S. C. Kim, E. S. Park, S. B. Kim, et al., Expression signature of E2F1 and its associated genes predict superficial to invasive progression of bladder tumors, J. Clin. Oncol., 28 (2010), 2660-2667. |
[40] | M. S. Lawrence, P. Stojanov, P. Polak, G. V. Kryukov, K. Cibulskis, A. Sivachenko, et al., Mutational heterogeneity in cancer and the search for new cancer-associated genes, Nature, 499 (2013), 214-218. |
[41] | K. H. Susek, M. Karvouni, E. Alici, A. Lundqvist, The role of CXC chemokine receptors 1-4 on immune cells in the tumor microenvironment, Front. Immunol., 9 (2018), 2159. |
[42] | V. Mollica Poeta, M. Massara, A. Capucetti, R. Bonecchi, chemokines and chemokine receptors: new targets for cancer immunotherapy, Front. Immunol., 10 (2019), 379. |
[43] | A. F. Schott, L. J. Goldstein, M. Cristofanilli, P. A. Ruffini, S. McCanna, J. M. Reuben, et al., Phase Ib pilot study to evaluate reparixin in combination with weekly paclitaxel in patients with HER-2-negative metastatic breast cancer, Clin. Cancer Res., 23 (2017), 5358-5365. |
[44] | T. M. Nywening, A. Wang-Gillam, D. E. Sanford, B. A. Belt, R. Z. Panni, B. M. Cusworth, et al., Targeting tumour-associated macrophages with CCR2 inhibition in combination with FOLFIRINOX in patients with borderline resectable and locally advanced pancreatic cancer: a single-centre, open-label, dose-finding, non-randomised, phase 1b trial, Lancet Oncol., 17 (2016), 651-662. |
[45] | C. Berlato, M. N. Khan, T. Schioppa, R. Thompson, E. Maniati, A. Montfort, et al., A CCR4 antagonist reverses the tumor-promoting microenvironment of renal cancer, J. Clin. Invest., 127 (2017), 801-813. |
[46] | D. Di Mitri, A. Toso, J. J. Chen, M. Sarti, S. Pinton, T. R. Jost, et al., Tumour-infiltrating Gr-1+ myeloid cells antagonize senescence in cancer, Nature, 515 (2014), 134-137. |
[47] |
C. Wang, A. Li, S. Yang, R. Qiao, X. Zhu, J. Zhang, CXCL5 promotes mitomycin C resistance in non-muscle invasive bladder cancer by activating EMT and NF-kappaB pathway, Biochem. Biophys. Res. Commun., 498 (2018), 862-868. doi: 10.1016/j.bbrc.2018.03.071
![]() |
[48] | I. Bieche, C. Chavey, C. Andrieu, M. Busson, S. Vacher, L. Le Corre, et al., CXC chemokines located in the 4q21 region are up-regulated in breast cancer, Endocr. Relat. Cancer, 14 (2007), 1039-1052. |
[49] | X. R. Yang, K. Brown, M. T. Landi, P. Ghiorzo, C. Badenas, M. Xu, et al., Duplication of CXC chemokine genes on chromosome 4q13 in a melanoma-prone family, Pigment Cell Melanoma Res., 25 (2012), 243-247. |
[50] | M. Lusti-Narasimhan, C. A. Power, B. Allet, S. Alouani, K. B. Bacon, J. J. Mermod, et al., Mutation of Leu25 and Val27 introduces CC chemokine activity into interleukin-8, J. Biol. Chem., 270 (1995), 2716-2721. |
[51] |
S. Y. Lim, A. E. Yuzhalin, A. N. Gordon-Weeks, R. J. Muschel, Targeting the CCL2-CCR2 signaling axis in cancer metastasis, Oncotarget, 7 (2016), 28697-28710. doi: 10.18632/oncotarget.7376
![]() |
[52] |
P. J. Sarvaiya, D. Guo, I. Ulasov, P. Gabikian, M. S. Lesniak, Chemokines in tumor progression and metastasis, Oncotarget, 4 (2013), 2171-2185. doi: 10.18632/oncotarget.1426
![]() |
[53] | F. C. Kischkel, D. A. Lawrence, A. Tinel, H. LeBlanc, A. Virmani, P. Schow, et al., Death receptor recruitment of endogenous caspase-10 and apoptosis initiation in the absence of caspase-8, J. Biol. Chem., 276 (2001), 46639-46646. |
[54] |
X. Chen, J. J. Subleski, H. Kopf, O. M. Howard, D. N. Mannel, J. J. Oppenheim, Cutting edge: expression of TNFR2 defines a maximally suppressive subset of mouse CD4+CD25+FoxP3+ T regulatory cells: applicability to tumor-infiltrating T regulatory cells, J. Immunol., 180 (2008), 6467-6471. doi: 10.4049/jimmunol.180.10.6467
![]() |
[55] | T. Tanaka, T. Imamura, M. Yoneda, A. Irie, H. Ogi, M. Nagata, et al., Enhancement of active MMP release and invasive activity of lymph node metastatic tongue cancer cells by elevated signaling via the TNF-alpha-TNFR1-NF-kappaB pathway and a possible involvement of angiopoietin-like 4 in lung metastasis, Int. J. Oncol., 49 (2016), 1377-1384. |
[56] | Y. Guo, Y. Feng, H. Liu, S. Luo, J. W. Clarke, P. G. Moorman, et al., Potentially functional genetic variants in the TNF/TNFR signaling pathway genes predict survival of patients with non-small cell lung cancer in the PLCO cancer screening trial, Mol. Carcinog., 58 (2019), 1094-1104. |
[57] | I. Zidi, S. Mestiri, A. Bartegi, N. B. Amor, TNF-alpha and its inhibitors in cancer, Med. Oncol., 27 (2010), 185-198. |
[58] |
X. Yuan, A. Gajan, Q. Chu, H. Xiong, K. Wu, G. S. Wu, Developing TRAIL/TRAIL death receptor-based cancer therapies, Cancer Metastasis Rev., 37 (2018), 733-748. doi: 10.1007/s10555-018-9728-y
![]() |
[59] |
S. Biswas, S. Roy Chowdhury, G. Mandal, S. Purohit, A. Gupta, A. Bhattacharyya, RelA driven co-expression of CXCL13 and CXCR5 is governed by a multifaceted transcriptional program regulating breast cancer progression, Biochim. Biophys. Acta Mol. Basis Dis., 1865 (2019), 502-511. doi: 10.1016/j.bbadis.2018.12.002
![]() |
[60] | B. S. Kochupurakkal, Z. C. Wang, T. Hua, A. C. Culhane, S. J. Rodig, K. Rajkovic-Molek, et al., RelA-induced interferon response negatively regulates proliferation, PLoS One, 10 (2015), e0140243. |
[61] | M. Lesina, S. M. Wormann, J. Morton, K. N. Diakopoulos, O. Korneeva, M. Wimmer, et al., RelA regulates CXCL1/CXCR2-dependent oncogene-induced senescence in murine Kras-driven pancreatic carcinogenesis, J. Clin. Invest., 126 (2016), 2919-2932. |
[62] | E. Ronin, M. Lubrano di Ricco, R. Vallion, J. Divoux, H. K. Kwon, S. Gregoire, et al., The NF-kappaB RelA transcription factor is critical for regulatory T cell activation and stability, Front. Immunol., 10 (2019), 2487. |
[63] | X. Qu, J. Han, Y. Zhang, X. Wang, H. Fan, F. Hua, et al., TLR4-RelA-miR-30a signal pathway regulates Th17 differentiation during experimental autoimmune encephalomyelitis development, J. Neuroinflammation, 16 (2019), 1-13. |
[64] | J. Concetti, C. L. Wilson, NFKB1 and cancer: friend or foe?, Cells, 7 (2018), 133. |
[65] | P. Li, J. Gu, X. Yang, H. Cai, J. Tao, X. Yang, et al., Functional promoter -94 ins/del ATTG polymorphism in NFKB1 gene is associated with bladder cancer risk in a Chinese population, PLoS One, 8 (2013), e71604. |
[66] | R. Roskoski Jr., Src protein-tyrosine kinase structure and regulation, Biochem. Biophys. Res. Commun., 324 (2004), 1155-1164. |
[67] | V. B. Venkadakrishnan, A. D. DePriest, S. Kumari, D. Senapati, S. Ben-Salem, Y. Su, et al., Protein Kinase N1 control of androgen-responsive serum response factor action provides rationale for novel prostate cancer treatment strategy, Oncogene, 38 (2019), 4496-4511. |
[68] | A. Franko, M. Kovarova, S. Feil, R. Feil, R. Wagner, M. Heni, et al., cGMP-dependent protein kinase I (cGKI) modulates human hepatic stellate cell activation, Metabolism, 88 (2018), 22-30. |
[69] | F. Schuler, S. Afreen, C. Manzl, G. Hacker, M. Erlacher, A. Villunger, Checkpoint kinase 1 is essential for fetal and adult hematopoiesis, EMBO Rep., 20 (2019), e47026. |
[70] |
D. D. Browning, I. K. Kwon, R. Wang, cGMP-dependent protein kinases as potential targets for colon cancer prevention and treatment, Future Med. Chem., 2 (2010), 65-80. doi: 10.4155/fmc.09.142
![]() |
[71] |
W. Lin, S. Chen, Checkpoint Kinase 1 Is Overexpressed during HPV16-Induced Cervical Carcinogenesis, Gynecol. Obstet. Invest., 83 (2018), 299-305. doi: 10.1159/000487943
![]() |
[72] | J. M. Lamar, Y. Xiao, E. Norton, Z. G. Jiang, G. M. Gerhard, S. Kooner, et al., SRC tyrosine kinase activates the YAP/TAZ axis and thereby drives tumor growth and metastasis, J. Biol. Chem., 294 (2019), 2302-2317. |
[73] | L. Mao, S. Liu, L. Hu, L. Jia, H. Wang, M. Guo, et al., miR-30 family: a promising regulator in development and disease, Biomed. Res. Int., 2018 (2018), 9623412. |
[74] | H. Yin, Y. Wang, Y. Wu, X. Zhang, X. Zhang, J. Liu, et al., EZH2-mediated Epigenetic Silencing of miR-29/miR-30 targets LOXL4 and contributes to tumorigenesis, metastasis, and immune microenvironment remodeling in breast cancer, Theranostics, 10 (2020), 8494-8512. |
[75] |
S. Jiang, D. Miao, M. Wang, J. Lv, Y. Wang, J. Tong, MiR-30-5p suppresses cell chemoresistance and stemness in colorectal cancer through USP22/Wnt/beta-catenin signaling axis, J. Cell. Mol. Med., 23 (2019), 630-640. doi: 10.1111/jcmm.13968
![]() |
[76] | B. Farhood, M. Najafi, K. Mortezaee, CD8(+) cytotoxic T lymphocytes in cancer immunotherapy: a review, J. Cell. Physiol., 234 (2019), 8509-8521. |
[77] | R. M. Bremnes, L. T. Busund, T. L. Kilvaer, S. Andersen, E. Richardsen, E. E. Paulsen, et al., The role of tumor-infiltrating lymphocytes in development, progression, and prognosis of non-small cell lung cancer, J. Thorac. Oncol., 11 (2016), 789-800. |
1. | Cláudia Martins-Lima, Ugo Chianese, Rosaria Benedetti, Lucia Altucci, Carmen Jerónimo, Margareta P. Correia, Tumor microenvironment and epithelial-mesenchymal transition in bladder cancer: Cytokines in the game?, 2023, 9, 2296-889X, 10.3389/fmolb.2022.1070383 | |
2. | Anqi Jin, Jin Zhou, Pengcheng Yu, Shichong Zhou, Cai Chang, High Expression of THBS1 Leads to a Poor Prognosis in Papillary Thyroid Cancer and Suppresses the Anti-Tumor Immune Microenvironment, 2022, 21, 1533-0346, 153303382210853, 10.1177/15330338221085360 | |
3. | Jianhui Li, Ye Nie, Weili Jia, Wenlong Wu, Wenjie Song, Yongxiang Li, Effect of Tertiary Lymphoid Structures on Prognosis of Patients with Hepatocellular Carcinoma and Preliminary Exploration of Its Formation Mechanism, 2022, 14, 2072-6694, 5157, 10.3390/cancers14205157 | |
4. | Lei Gu, Chunhui Jiang, Chunjie Xu, Ye Liu, Hong Zhou, Based on Molecular Subtypes, Immune Characteristics and Genomic Variation to Constructing and Verifying Multi-Gene Prognostic Characteristics of Colorectal Cancer, 2022, 10, 2296-634X, 10.3389/fcell.2022.828415 | |
5. | Yongli Situ, Xiaoyong Lu, Yongshi Cui, Qinying Xu, Li Deng, Hao Lin, Zheng Shao, Jv Chen, Aziz ur Rehman Aziz, Systematic Analysis of CXC Chemokine–Vascular Endothelial Growth Factor A Network in Colonic Adenocarcinoma from the Perspective of Angiogenesis, 2022, 2022, 2314-6141, 1, 10.1155/2022/5137301 | |
6. | Qianmiao Wu, Huaijun Tu, Jian Li, Multifaceted Roles of Chemokine C-X-C Motif Ligand 7 in Inflammatory Diseases and Cancer, 2022, 13, 1663-9812, 10.3389/fphar.2022.914730 | |
7. | Jan Korbecki, Mateusz Bosiacki, Iwona Szatkowska, Patrycja Kupnicka, Dariusz Chlubek, Irena Baranowska-Bosiacka, The Clinical Significance and Involvement in Molecular Cancer Processes of Chemokine CXCL1 in Selected Tumors, 2024, 25, 1422-0067, 4365, 10.3390/ijms25084365 | |
8. | Shi Fu, Zhiyong Tan, Hongjin Shi, Junhao Chen, Yawei Zhang, Chunming Guo, Wei Feng, Haole Xu, Jiansong Wang, Haifeng Wang, Development of a stemness-related prognostic index to provide therapeutic strategies for bladder cancer, 2024, 8, 2397-768X, 10.1038/s41698-024-00510-3 | |
9. | Zhenshen Bao, Minzhen Liao, Wanqi Dong, Yanhao Huo, Xianbin Li, Peng Xu, Wenbin Liu, A Pan-cancer Analysis Reveals the Tissue Specificity and Prognostic Impact of Angiogenesis-associated Genes in Human Cancers, 2023, 18, 15748936, 670, 10.2174/1574893618666230518163353 | |
10. | Jie Shen, Haidong Xia, Xiaojing Zhou, Lei Zhang, Qian Gao, Kan He, Dahai Liu, Bei Huang, Selenium enhances photodynamic therapy of C-phycocyanin against lung cancer via dual regulation of cytotoxicity and antioxidant activity, 2023, 55, 1672-9145, 1925, 10.3724/abbs.2023159 | |
11. | Guang‐hua Nie, Cheng‐yi Liu, Zhao Tian, Comprehensive analysis of prognostic and immunological role of basement membrane‐related genes in soft tissue sarcoma, 2024, 12, 2050-4527, 10.1002/iid3.70037 | |
12. | Tao Yin, Shuanzhu Mou, Haiyu Zhang, Ying Dong, Bing Yan, Weisheng Huang, Yuhan Liu, Hongbing Mei, CXCL10 could be a prognostic and immunological biomarker in bladder cancer, 2024, 15, 2730-6011, 10.1007/s12672-024-00982-6 | |
13. | Lakshmi Thangavelu, Ahsas Goyal, Muhammad Afzal, Ehssan Moglad, Sushama Rawat, Imran Kazmi, Sami I. Alzarea, Waleed Hassan Almalki, Richa Rani, Pusuluri Madhubabu, Pranchal Rajput, Pooja Bansal, Pyroptosis in lung cancer: The emerging role of non-coding RNAs, 2024, 263, 03440338, 155619, 10.1016/j.prp.2024.155619 | |
14. | Zaineb Hassouneh, Michelle E. Kim, Natalia Bowman, Manjeet Rao, Nu Zhang, Gang Huang, Robert S. Svatek, Neelam Mukherjee, The Laws of Attraction: Chemokines as Critical Mediators in Cancer Progression and Immunotherapy Response in Bladder Cancer, 2024, 16, 2072-6694, 3303, 10.3390/cancers16193303 | |
15. | Yuxuan Bao, Chang Tong, Xiangyang Xiong, CXCL3: A key player in tumor microenvironment and inflammatory diseases, 2024, 348, 00243205, 122691, 10.1016/j.lfs.2024.122691 | |
16. | Jiaqi Wang, Xinting Ouyang, Weijian Zhu, Qiang Yi, Jinghua Zhong, The Role of CXCL11 and its Receptors in Cancer: Prospective but Challenging Clinical Targets, 2024, 31, 1073-2748, 10.1177/10732748241241162 | |
17. | María García-Pola, Lucía Rodríguez-Fonseca, Claudia Llorente-Álvarez, Santiago Llorente-Pendás, History of Breast Cancer in Patients with Oral Lichen Planus: A Case–Control Study, 2024, 13, 2077-0383, 7208, 10.3390/jcm13237208 | |
18. | Shiqiang Dong, Lili Wang, Xinyu Liu, Dingkun Hou, Qing Liu, Ji Zheng, Haitao Wang, CELF2 inhibits bladder cancer progression by decreasing the stability of CXCL5, 2025, 00243205, 123585, 10.1016/j.lfs.2025.123585 | |
19. | Jan Korbecki, Mateusz Bosiacki, Maciej Pilarczyk, Marcin Kot, Piotr Defort, Ireneusz Walaszek, Dariusz Chlubek, Irena Baranowska-Bosiacka, The CXCL1-CXCR2 Axis as a Component of Therapy Resistance, a Source of Side Effects in Cancer Treatment, and a Therapeutic Target, 2025, 17, 2072-6694, 1674, 10.3390/cancers17101674 |
TLR | Type | Fold change | P-value | t-test | References |
CXCL2 | Superficial Bladder Cancer | -12.716 | 7.74E-17 | -11.618 | [38] |
CXCL3 | Superficial Bladder Cancer | -3.564 | 2.44E-08 | -6.122 | [38] |
CXCL6 | Infiltrating Bladder Urothelial Carcinoma | 2.272 | 1.34E-06 | 5.02 | [38] |
CXCL9 | Infiltrating Bladder Urothelial Carcinoma | 2.84 | 3.92E-09 | 6.197 | [38] |
CXCL10 | Infiltrating Bladder Urothelial Carcinoma | 3.039 | 6.79E-11 | 7.04 | [38] |
CXCL12 | Superficial Bladder Cancer | -4.495 | 2.52E-21 | -12.211 | [39] |
Infiltrating Bladder Urothelial Carcinoma | -3.371 | 1.41E-12 | -7.727 | [39] | |
Superficial Bladder Cancer | -2.531 | 1.92E-12 | -8.846 | [38] | |
CXCL13 | Infiltrating Bladder Urothelial Carcinoma | 6.369 | 1.07E-10 | 7.171 | [38] |
Key TF | Description | Overlapped genes | P value | Q value | List of overlapped genes |
RELA | v-rel reticuloendotheliosis viral oncogene homolog A (avian) | 5 | 1.16E-06 | 1.79E-06 | CXCL10, CXCL12, CXCL2, CXCL5, CXCL1 |
NFKB1 | nuclear factor of kappa light polypeptide gene enhancer in B-cells 1 | 6 | 1.16E-06 | 1.79E-06 | CXCL12, CXCL10, CXCL5, CXCL1, CXCL2 |
SP1 | Sp1 transcription factor | 2 | 0.0406 | 0.0406 | CXCL5, CXCL1 |
CXC chemokines | Kinase target | Description | Leading Edge Number | P Value |
CXCL1 | Kinase_LYN | LYN proto-oncogene, Src family tyrosine kinase | 17 | 0 |
Kinase_IKBKB | inhibitor of nuclear factor kappa B kinase subunit beta | 9 | 0 | |
CXCL2 | Kinase_LCK | LCK proto-oncogene, Src family tyrosine kinase | 24 | 0 |
Kinase_LYN | LYN proto-oncogene, Src family tyrosine kinase | 21 | 0 | |
CXCL3 | Kinase_LYN | LYN proto-oncogene, Src family tyrosine kinase | 19 | 0 |
Kinase_LCK | LCK proto-oncogene, Src family tyrosine kinase | 25 | 0 | |
CXCL5 | Kinase_LCK | LCK proto-oncogene, Src family tyrosine kinase | 24 | 0 |
Kinase_LYN | LYN proto-oncogene, Src family tyrosine kinase | 18 | 0 | |
CXCL6 | Kinase_SYK | spleen associated tyrosine kinase | 16 | 0 |
Kinase_LYN | LYN proto-oncogene, Src family tyrosine kinase | 24 | 0 | |
CXCL7 | Kinase_PKN1 | protein kinase N1 | 5 | 0.009 |
Kinase_ADRBK1 | G protein-coupled receptor kinase 2 | 10 | 0 | |
CXCL9 | Kinase_LYN | LYN proto-oncogene, Src family tyrosine kinase | 18 | 0 |
Kinase_LCK | LCK proto-oncogene, Src family tyrosine kinase | 20 | 0 | |
CXCL10 | Kinase_LCK | LCK proto-oncogene, Src family tyrosine kinase | 20 | 0 |
Kinase_LYN | LYN proto-oncogene, Src family tyrosine kinase | 16 | 0 | |
CXCL11 | Kinase_LCK | LCK proto-oncogene, Src family tyrosine kinase | 20 | 0 |
Kinase_LYN | LYN proto-oncogene, Src family tyrosine kinase | 17 | 0 | |
CXCL12 | Kinase_LYN | LYN proto-oncogene, Src family tyrosine kinase | 20 | 0 |
Kinase_PRKG1 | protein kinase, cGMP-dependent, type I | 15 | 0.004 | |
CXCL13 | Kinase_LCK | LCK proto-oncogene, Src family tyrosine kinase | 22 | 0 |
Kinase_LYN | LYN proto-oncogene, Src family tyrosine kinase | 19 | 0 | |
CXCL16 | Kinase_LCK | LCK proto-oncogene, Src family tyrosine kinase | 22 | 0 |
Kinase_JAK2 | Janus kinase 2 | 7 | 0 | |
CXCL17 | Kinase_CHEK1 | checkpoint kinase 1 | 45 | 0 |
Kinase_PLK1 | polo like kinase 1 | 41 | 0 |
CXC chemokines | MIRNA targets | Leading Edge Number | P Value |
CXCL1 | GTAAGAT, MIR-200A | 20 | 0.009 |
TGTTTAC, MIR-30A-5P, MIR-30C, MIR-30D, MIR-30B, MIR-30E-5P | 133 | 0 | |
CXCL3 | AACATTC, MIR-409-3P | 47 | 0.009 |
CXCL5 | CTCAAGA, MIR-526B | 24 | 0.019 |
CAGTGTT, MIR-141, MIR-200A | 72 | 0.015 | |
CXCL6 | TATTATA, MIR-374 | 89 | 0 |
TGTTTAC, MIR-30A-5P, MIR-30C, MIR-30D, MIR-30B, MIR-30E-5P | 192 | 0 | |
CXCL7 | GTATGAT, MIR-154, MIR-487 | 21 | 0 |
GCGCCTT, MIR-525, MIR-524 | 7 | 0.044 | |
CXCL17 | CCAGGTT, MIR-490 | 10 | 0.025 |
coef | HR | 95%CI_l | 95%CI_u | p value | sig | |
B_cell | -3.430 | 0.032 | 0.001 | 0.739 | 0.032 | * |
CD8_Tcell | 1.675 | 5.337 | 0.294 | 96.758 | 0.257 | |
CD4_Tcell | -0.733 | 0.480 | 0.008 | 29.113 | 0.726 | |
Macrophage | 4.651 | 104.684 | 6.900 | 1588.200 | 0.001 | ** |
Neutrophil | -3.390 | 0.034 | 0.000 | 17.688 | 0.289 | |
Dendritic | -0.073 | 0.929 | 0.171 | 5.054 | 0.933 | |
CXCL1 | 0.101 | 1.106 | 0.968 | 1.265 | 0.138 | |
CXCL2 | -0.279 | 0.756 | 0.612 | 0.934 | 0.010 | * |
CXCL3 | 0.192 | 1.212 | 0.929 | 1.581 | 0.157 | |
CXCL5 | 0.034 | 1.035 | 0.895 | 1.196 | 0.645 | |
CXCL6 | -0.083 | 0.921 | 0.801 | 1.059 | 0.247 | |
CXCL7 | 0.138 | 1.148 | 0.991 | 1.329 | 0.066 | |
CXCL9 | -0.028 | 0.973 | 0.808 | 1.171 | 0.770 | |
CXCL10 | 0.048 | 1.050 | 0.843 | 1.306 | 0.665 | |
CXCL11 | -0.060 | 0.942 | 0.762 | 1.166 | 0.584 | |
CXCL12 | 0.089 | 1.093 | 0.961 | 1.242 | 0.176 | |
CXCL13 | 0.040 | 1.041 | 0.919 | 1.178 | 0.531 | |
CXCL16 | -0.091 | 0.913 | 0.753 | 1.107 | 0.353 | |
CXCL17 | -0.010 | 0.990 | 0.926 | 1.057 | 0.756 |
TLR | Type | Fold change | P-value | t-test | References |
CXCL2 | Superficial Bladder Cancer | -12.716 | 7.74E-17 | -11.618 | [38] |
CXCL3 | Superficial Bladder Cancer | -3.564 | 2.44E-08 | -6.122 | [38] |
CXCL6 | Infiltrating Bladder Urothelial Carcinoma | 2.272 | 1.34E-06 | 5.02 | [38] |
CXCL9 | Infiltrating Bladder Urothelial Carcinoma | 2.84 | 3.92E-09 | 6.197 | [38] |
CXCL10 | Infiltrating Bladder Urothelial Carcinoma | 3.039 | 6.79E-11 | 7.04 | [38] |
CXCL12 | Superficial Bladder Cancer | -4.495 | 2.52E-21 | -12.211 | [39] |
Infiltrating Bladder Urothelial Carcinoma | -3.371 | 1.41E-12 | -7.727 | [39] | |
Superficial Bladder Cancer | -2.531 | 1.92E-12 | -8.846 | [38] | |
CXCL13 | Infiltrating Bladder Urothelial Carcinoma | 6.369 | 1.07E-10 | 7.171 | [38] |
Key TF | Description | Overlapped genes | P value | Q value | List of overlapped genes |
RELA | v-rel reticuloendotheliosis viral oncogene homolog A (avian) | 5 | 1.16E-06 | 1.79E-06 | CXCL10, CXCL12, CXCL2, CXCL5, CXCL1 |
NFKB1 | nuclear factor of kappa light polypeptide gene enhancer in B-cells 1 | 6 | 1.16E-06 | 1.79E-06 | CXCL12, CXCL10, CXCL5, CXCL1, CXCL2 |
SP1 | Sp1 transcription factor | 2 | 0.0406 | 0.0406 | CXCL5, CXCL1 |
CXC chemokines | Kinase target | Description | Leading Edge Number | P Value |
CXCL1 | Kinase_LYN | LYN proto-oncogene, Src family tyrosine kinase | 17 | 0 |
Kinase_IKBKB | inhibitor of nuclear factor kappa B kinase subunit beta | 9 | 0 | |
CXCL2 | Kinase_LCK | LCK proto-oncogene, Src family tyrosine kinase | 24 | 0 |
Kinase_LYN | LYN proto-oncogene, Src family tyrosine kinase | 21 | 0 | |
CXCL3 | Kinase_LYN | LYN proto-oncogene, Src family tyrosine kinase | 19 | 0 |
Kinase_LCK | LCK proto-oncogene, Src family tyrosine kinase | 25 | 0 | |
CXCL5 | Kinase_LCK | LCK proto-oncogene, Src family tyrosine kinase | 24 | 0 |
Kinase_LYN | LYN proto-oncogene, Src family tyrosine kinase | 18 | 0 | |
CXCL6 | Kinase_SYK | spleen associated tyrosine kinase | 16 | 0 |
Kinase_LYN | LYN proto-oncogene, Src family tyrosine kinase | 24 | 0 | |
CXCL7 | Kinase_PKN1 | protein kinase N1 | 5 | 0.009 |
Kinase_ADRBK1 | G protein-coupled receptor kinase 2 | 10 | 0 | |
CXCL9 | Kinase_LYN | LYN proto-oncogene, Src family tyrosine kinase | 18 | 0 |
Kinase_LCK | LCK proto-oncogene, Src family tyrosine kinase | 20 | 0 | |
CXCL10 | Kinase_LCK | LCK proto-oncogene, Src family tyrosine kinase | 20 | 0 |
Kinase_LYN | LYN proto-oncogene, Src family tyrosine kinase | 16 | 0 | |
CXCL11 | Kinase_LCK | LCK proto-oncogene, Src family tyrosine kinase | 20 | 0 |
Kinase_LYN | LYN proto-oncogene, Src family tyrosine kinase | 17 | 0 | |
CXCL12 | Kinase_LYN | LYN proto-oncogene, Src family tyrosine kinase | 20 | 0 |
Kinase_PRKG1 | protein kinase, cGMP-dependent, type I | 15 | 0.004 | |
CXCL13 | Kinase_LCK | LCK proto-oncogene, Src family tyrosine kinase | 22 | 0 |
Kinase_LYN | LYN proto-oncogene, Src family tyrosine kinase | 19 | 0 | |
CXCL16 | Kinase_LCK | LCK proto-oncogene, Src family tyrosine kinase | 22 | 0 |
Kinase_JAK2 | Janus kinase 2 | 7 | 0 | |
CXCL17 | Kinase_CHEK1 | checkpoint kinase 1 | 45 | 0 |
Kinase_PLK1 | polo like kinase 1 | 41 | 0 |
CXC chemokines | MIRNA targets | Leading Edge Number | P Value |
CXCL1 | GTAAGAT, MIR-200A | 20 | 0.009 |
TGTTTAC, MIR-30A-5P, MIR-30C, MIR-30D, MIR-30B, MIR-30E-5P | 133 | 0 | |
CXCL3 | AACATTC, MIR-409-3P | 47 | 0.009 |
CXCL5 | CTCAAGA, MIR-526B | 24 | 0.019 |
CAGTGTT, MIR-141, MIR-200A | 72 | 0.015 | |
CXCL6 | TATTATA, MIR-374 | 89 | 0 |
TGTTTAC, MIR-30A-5P, MIR-30C, MIR-30D, MIR-30B, MIR-30E-5P | 192 | 0 | |
CXCL7 | GTATGAT, MIR-154, MIR-487 | 21 | 0 |
GCGCCTT, MIR-525, MIR-524 | 7 | 0.044 | |
CXCL17 | CCAGGTT, MIR-490 | 10 | 0.025 |
coef | HR | 95%CI_l | 95%CI_u | p value | sig | |
B_cell | -3.430 | 0.032 | 0.001 | 0.739 | 0.032 | * |
CD8_Tcell | 1.675 | 5.337 | 0.294 | 96.758 | 0.257 | |
CD4_Tcell | -0.733 | 0.480 | 0.008 | 29.113 | 0.726 | |
Macrophage | 4.651 | 104.684 | 6.900 | 1588.200 | 0.001 | ** |
Neutrophil | -3.390 | 0.034 | 0.000 | 17.688 | 0.289 | |
Dendritic | -0.073 | 0.929 | 0.171 | 5.054 | 0.933 | |
CXCL1 | 0.101 | 1.106 | 0.968 | 1.265 | 0.138 | |
CXCL2 | -0.279 | 0.756 | 0.612 | 0.934 | 0.010 | * |
CXCL3 | 0.192 | 1.212 | 0.929 | 1.581 | 0.157 | |
CXCL5 | 0.034 | 1.035 | 0.895 | 1.196 | 0.645 | |
CXCL6 | -0.083 | 0.921 | 0.801 | 1.059 | 0.247 | |
CXCL7 | 0.138 | 1.148 | 0.991 | 1.329 | 0.066 | |
CXCL9 | -0.028 | 0.973 | 0.808 | 1.171 | 0.770 | |
CXCL10 | 0.048 | 1.050 | 0.843 | 1.306 | 0.665 | |
CXCL11 | -0.060 | 0.942 | 0.762 | 1.166 | 0.584 | |
CXCL12 | 0.089 | 1.093 | 0.961 | 1.242 | 0.176 | |
CXCL13 | 0.040 | 1.041 | 0.919 | 1.178 | 0.531 | |
CXCL16 | -0.091 | 0.913 | 0.753 | 1.107 | 0.353 | |
CXCL17 | -0.010 | 0.990 | 0.926 | 1.057 | 0.756 |