The Clinical Implications and Molecular Mechanism of CX3CL1 Expression in Urothelial Bladder Cancer

Background: CX3CL1 is a chemokine that may play important roles in cancer immune regulation. Its mechanism in bladder cancer (BCa) is poorly understood. The objective of the current study was to evaluate the association between CX3CL1 and BCa and the related biological mechanisms.

Methods: A total of 277 patients with BCa were enrolled in the present study. The association between CX3CL1 expression and disease outcome was evaluated. In vitro and in vivo experiments were performed using the TCCSUP cell line to investigate the function of CX3CL1 in BCa.

Results: Compared with low expression, high expression of CX3CL1 was significantly associated with poorer progression-free survival (hazard ratio [HR]=2.03, 95% confidence interval [95% CI]: 1.26-3.27, P=0.006), cancer-specific survival (HR=2.16, 95% CI: 1.59-2.93, P<0.001), and overall survival (HR=1.55, 95% CI: 1.08-2.24, P=0.039). Multivariable Cox regression analysis suggested that CX3CL1 was an independent prognostic factor for BCa outcomes. In vitro and in vivo experiments indicated that high expression of CX3CL1 was significantly associated with cell proliferation (P<0.001) and invasion (P<0.001). Gene expression profiling results showed that after CX3CL1 knockdown, CDH1 was significantly upregulated, while ETS1, RAF1, and EIF4E were significantly downregulated. Pathway enrichment analysis suggested that the ERK/MAPK signaling pathway was significantly inhibited (P<0.001).

Conclusions: CX3CL1 is an independent predictor of a poor prognosis in BCa and can promote the proliferation and invasion of BCa cells.

Introduction

With an estimated 573,000 new cases and 212,000 deaths in 2020, urinary bladder cancer (BCa) has become one of the most common malignancies worldwide (1). In the United States, BCa accounted for 7% of all malignancy cases and ranked fourth among all types of malignancies in 2020 (2). Chemical components from cigarettes and industrial products are well-known risk factors; however, the etiology of BCa is poorly understood.

Previous evidence has shown that the immune system might play important roles in the development and progression of BCa, involving immune surveillance and immune tolerance. In immunoregulatory pathways, chemokines serve as key factors (3). Tumor immunotherapy represented by PD1/PD-L1 has attracted great attention in recent years. The effect of chemokines on tumors has been proven to produce a two-way regulation pattern: first, chemokines may enhance antitumor immunity by promoting immune cell recruitment to the tumor site (4); second, autonomous secretion of chemokines and receptors from tumor cells may promote tumor growth by inhibiting antitumor immunity (5).

A unique class of chemokines-the CX3C family has been reported playing a critical role in tumor immune regulation, in which CX3CL1 is the only member (6). It contains a unique 29nm long pedicle structure (mucus domain) in N prime/C prime. This structure allows chemokines to extend from epithelial cells to contact flowing white blood cells, and is independent of selectins and integrins (7). In addition, it can also induce intracellular calcium flow and chemotaxis without activating integrins (7). The receptor of CX3CL1, CX3CR1, is a transmembrane domain receptor, which mediates cell adhesion to CX3CL1 without G protein activation (8). Therefore, the unique structure and function of CX3CL1 has attracted the attention of many researchers.

As an important molecule in the chemokine family, C-X3-C Motif Chemokine Ligand 1 (CX3CL1) was found to be significantly associated with the development and progression of multiple cancers (9–13). However, the roles of CX3CL1 and its receptors are still controversial in terms of the biological regulation of tumors. Studies have suggested that CX3CL1 exerts its tumor-suppressive effect by activating immune cells, including natural killer T cells, and inhibiting tumor cell aggregation and adhesion (14, 15). For example, in a study of breast cancer, the high expression of CX3CL1 might recruit immune cells, such as CD8+ T cells and NK cells, to exert anti-tumor activity (16). Another study showed that patients with high expression of CX3CL1 and CX3CL1 receptors would have a significantly better prognosis in liver cancer (15). However, CX3CL1 may also produce a cancer-promotive effect as an adhesion molecule (9, 13). It was reported that high expression of CX3CL1 can induce tumor cell proliferation and promote the transformation of the cell cycle to S phase in gastric cancer (17). Shulby et al. indicated that CX3CL1 involved in the invasion process of prostate cancer metastasis to bone marrow which could be significantly reduced by anti-CX3CL1 antibody (13). This two-way regulation pattern of CX3CL1 is still unclear in BCa.

The objective of the current study was to evaluate the association between CX3CL1 and BCa and to investigate the underlying biological mechanisms involving CX3CL1 and its receptors in BCa.

Materials and Methods

Fresh Tissue Sample

Seven fresh tumor tissue samples of bladder cancer and normal bladder mucosa tissue were collected by TURBT (transurethral resection of bladder tumor). In the end, the specimens due to severe burning and too small volume were eliminated. The qualified samples were divided into two groups (4 cases in the normal tissue group; 5 cases in the tumor tissue group) for Real-Time PCR experiment to determine the difference in CX3CL1 mRNA between the two groups.

Clinical Samples

Patients with BCa undergoing surgical treatment were enrolled in the present study between December 2011 and December 2016 at Huashan Hospital (Shanghai, China). The inclusion criteria were as follows: (a) sufficient clinical information and follow-up information were available, (b) pathological diagnosis and grading were independently confirmed by three pathologists as transitional cell carcinoma of the bladder, and (c) sufficient formalin-fixed paraffin-embedded (FFPE) specimens could be obtained for further investigation. Informed consent was obtained from all patients for institutional biobank, and the study was approved by the Institutional Review Board at Huashan Hospital. Pathological grading and staging were performed according to the WHO/ISUP (2004) and UICC-TNM (2002) guidelines (18, 19). Ta-T1 was defined as non-muscle-invasive BCa (NMIBC), and T2-T4 were defined as muscle-invasive BCa (MIBC).

Immunohistochemistry

Immunohistochemistry staining of FFPE samples was performed using an anti-CX3CL1 antibody (Abcam, Cambridge, UK; 1:100) and the avidin-biotin-complex method. The extent and intensity of tumor regions (without necrotic tissue) were further evaluated under a Nikon 80i microscope (Tokyo, Japan). Staining intensity was categorized as 0 (absent), 1 (weak), 2 (moderate), or 3 (strong). An H-score was then derived by multiplying the staining intensity and the percentage of cells stained. Given that there is still no established cutoff value for CX3CL1 expression, a cutoff value of 10% immuno-positive tumor cells was recommended by the 3 experienced pathologists (Figure S1). Disagreement as to whether >10% of cells were positively stained was resolved by consensus.

RNA Isolation and Quantitative PCR (qPCR)

Total RNA was extracted from tissues with TRIzol (Pufei, Shanghai, China) and retrotranscribed using PrimeScript RT Master Mix (Perfect Real Time; TaKaRa, Shiga, Japan) following the manufacturer’s recommendations. For qPCR, we used SYBR Premix Ex Taq (Tli RNaseH Plus; TaKaRa) and the 7500 Real-Time PCR System with Dell Tower (Applied Biosystems, Forster City, USA) system. The primers used to amplify CX3CL1, GAPDH, etc. were purchased from Ji Kai (Shanghai, China). The sequences of the primers for CX3CL1 were 5’-GTAGCTTTGCTCATCCACTATCA-3’ (upstream) and 5’-GACCACAGACTCGTCCATTCC-3’ (downstream). For each sample, the average value of the threshold cycle was normalized to the GAPDH level with the 2^-ΔΔCt method.

Cell Lines and Viral Infection

T24 (RRID: CVCL_0554), 5637 (RRID: CVCL_0126) and TCCSUP (RRID: CVCL_1738) bladder urothelial cancer cells were purchased from Cell Bank of Chinese Academy of Sciences (Shanghai, China) and cultured in DMEM (Gibco, Waltham, USA) supplemented with 10% fetal bovine serum (FBS; PAA).

Approximately 1×10⁶ cells were cultured for 12 h to allow attachment to plates. Cells were rinsed once with PBS, and 10 mL of serum-free medium was added to the plate. We used an adenovirus (Ji Kai, Shanghai, China) that overexpressed shRNA targeting CX3CL1 to infect the cells at a multiplicity of infection (MOI) of 10. Sixteen hours after infection, the medium was replaced with 10 mL of complete medium.

Cell Viability Assays and Transwell Assays

Cell viability was assessed at the indicated time points using an MTT kit (Gen-View, Calimesa, USA). MTT reagent (20 μL/well) was added to the medium, and the absorbance was read using a microplate reader at 490 nm for the test length and 570 nm for the reference. The percentage of viable cells was calculated based on the absorbance of the PBS control. Inserts in 24-well Transwell plates (Corning Costar, New York, USA) were coated with 500 μL/well Matrigel (BD Bioscience, San Jose, USA). Rehydration was performed with 500 μL/well serum-free medium for 2 h in a 37°C incubator (Sanyo, Japan). Samples were seeded into the inserts after removal of the rehydration medium. After 24 h of incubation, the cells that invaded through the membrane were stained with crystal violet (Beyotime, Shanghai, China) and observed at 200× magnification using a light microscope (Olympus, Tokyo, Japan).

RNA-Seq Expression Profiling and Data Analysis

CX3CL1 was knocked down in TCCSUP cells by lentiviral infection [virus strain: LV-CX3CL1-RNAi (52620-1)]. RNA was extracted using the TRIzol method used for qPCR. cDNA was obtained by reverse transcription (Promega M-MLV kit) via a real-time qPCR detection system (Type Roche-LightCycler 480). The knockdown efficiency of the CX3CL1 gene was also evaluated by qPCR.

GeneChip PrimeView detection: RNA samples were subjected to quality inspection analysis (Agilent 2100 Bioanalyzer), and the GeneChip 3’IVT Express Kit was used to prepare amplified RNA (aRNA). That is, cDNA was obtained through first-strand synthesis, a double-stranded DNA template was obtained through second-strand synthesis, and then biotinylated aRNA was obtained by in vitro inversion. The aRNA was purified, fragmented, hybridized, washed with a chip probe (GeneChip Hybridization Wash and Stain Kit), and finally scanned to obtain images and raw data.

Statistical analysis and verification: Quality control was performed on the original data, statistical testing was used to identify differential genes between paired samples, and scatter plots, volcano plots, cluster plots and other chip analysis result plots were drawn. Ingenuity Pathway Analysis (IPA) database software was used to perform disease and function enrichment analysis and pathway enrichment analysis of the differentially expressed gene list; the results were verified by Western blot analysis.

Xenografts and Animal Models

All experimental animals were obtained from Shanghai SLAC Laboratory Animal Co., Ltd. This work was carried out with 16 (8 male and 8 female) BALB/C nude mice aged 4-6 weeks and weighing 15 to 20 grams. All mice were kept under SPF living conditions. Mice were divided randomly and equally into 2 groups. The control group was named the negative control (NC) group and injected with normal TCCSUP cells, and the experimental group was named the KD group and injected with TCCSUP CX3CL1-knockdown cells. Each mouse was injected with 200 µl of cell suspension in the shoulder, which contained 1×10⁷ cells. We measured tumor volume every 5 days after injection, and calculated tumor volume with the following formula: Volume = 0.524 × L × W2 (L: long diameter; W: short diameter). The nude mice were sacrificed 35 days after injection, and tumor weight was measured. The subcutaneous tumors from each nude mouse were evaluated in pathological paraffin-embedded sections.

Statistical Analysis

SPSS 21.0 (IBM, New York, USA) was used for statistical analyses. All data are presented as the mean ± standard deviation (SD). Student’s t-test and the Wilcoxon test were applied to compare gene expression and staining. Univariable and multivariable Cox/logistic regression analyses were used to evaluate the correlation between CX3CL1 expression and prognosis. Kaplan-Meier curve analyses and the log-rank test were also performed. A P value <0.05 was considered statistically significant.

Results

The expression level of CX3CL1 mRNA in bladder cancer tissues was significantly higher than that in normal tissues (Figure S2, P=0.0268). A total of 277 patients were enrolled in the present study, with an average follow-up time of 49 months (1 to 72 months). Among them, 50 (18.1%) were negative for CX3CL1 expression, and 227 (81.9%) were considered positive for CX3CL1 expression. The demographic and clinical characteristics of the patients are summarized in Table 1. The patients were divided into 2 groups according to their expression of CX3CL1. By univariable analysis, the expression status of CX3CL1 was significantly associated with tumor stage [odds ratio (OR)=2.55, 95% confidence interval (95%CI): 1.36-4.76, P=0.003], tumor size (OR=2.06, 95%CI: 1.10-3.86, P=0.025), tumor grade (OR=1.91, 95%CI: 1.02-3.55, P=0.043), disease recurrence (OR=2.05, 95%CI: 1.08-3.88, P=0.029), metastatic disease (OR=3.98, 95%CI: 1.95-8.13, P<0.001), Ki67 expression (OR=1.02, 95%CI: 1.01-1.04, P=0.002) and mortality (OR=4.90, 95%CI: 2.53-9.47, P<0.001, Table S1). Multivariable analysis showed that CX3CL1 expression remained significantly associated with tumor size, tumor stage, Ki67 expression, disease recurrence, and mortality (all P<0.05).

TABLE 1

Table 1 Clinical characteristics of entire cohort and subgroups by CX3CL1 expression.

CX3CL1 was an independent predictor of BCa recurrence in both a univariable analysis (OR=2.44, 95%CI: 1.24-2.83, P=0.001) and a multivariable analysis (OR=2.37, 95%CI: 1.18-2.75, P=0.006) after adjusting for multiple variables, as shown in Table 2. As shown in the cumulative hazard curve of disease recurrence (Figure 1A), the 6-year cumulative recurrence rate of the BCa CX3CL1-positive group was 62.9%, which was significantly higher than that of the BCa CX3CL1-negative group (31.8%, P=0.006). Similarly, the BCa-specific survival rate of CX3CL1-positive group was 23.2%, which was significantly lower than that of the negative expression group (41.1%, P<0.0001, Figure 1B). Patients with high expression of CX3CL1 had an ~2-fold increased risk of cancer-specific death based on univariable (HR=2.16, 95%CI: 1.59-2.93, P<0.001) and multivariable (HR=2.17, 95%CI: 1.08-2.39, P=0.006) analyses (Table 3). This indicated that CX3CL1 was also a significant and independent risk factor for disease-specific death. A total of 406 BCa samples were obtained from the TCGA database. Survival analysis suggested that the patients with high expression of CX3CL1 (FPKM≥1.0) were significantly associated with poor survival compared with the patients with low CX3CL1 expression (log-rank P value=0.039, Figure 1C). These results remained significant after adjusting for age at diagnosis and tumor stage via Cox regression (HR=1.76, 95%CI: 1.14-2.71, P=0.011)

TABLE 2

Table 2 Univariable and multivariable logistic regression analyses on tumor recurrence.

FIGURE 1

Figure 1 Survival curves. (A) Cumulative tumor recurrence rate: the blue line represents the cumulative recurrence rate of 277 patients with bladder cancer with positive CX3CL1 expression; the red line represents the CX3CL1 expression-negative group. (B) Cancer-specific survival: the blue line represents the cancer-specific survival of the bladder cancer CX3CL1-positive group; the red line represents the CX3CL1 expression-negative group. (C) Overall mortality in the TCGA database: the blue line represents the overall mortality associated with 406 bladder cancer samples obtained from the TCGA database with positive CX3CL1 expression; the red line represents the CX3CL1 expression-negative group. HR, hazard ratio; 95%CI, 95% confidence interval.

TABLE 3

Table 3 Univariable and multivariable-adjusted hazard ratios for cancer-specific mortality in Cox proportional hazards models.

To evaluate whether CX3CL1 expression is associated with disease aggressiveness with functional experiments, we evaluated the expression levels in 3 different types of cell lines, TCCSUP (grade IV cancer), 5637 (grade II cancer), and T24 cells (high-grade and invasive TCC). The mRNA expression of CX3CL1 in TCCSUP cells was significantly higher than that in 5637 and T24 cells (P<0.001, Figure 2A). No significant difference in CX3CL1 expression was observed between the 5637 and T24 cell lines (P=0.766). Thus, TCCSUP cells were used for further evaluation. CX3CL1 was then knocked down by lentiviral infection. The knockdown efficiency was then verified by qPCR (efficiency reached 56.5%, P=0.0032, Figure 2B). We then evaluated whether CX3CL1 expression is associated with cancer cell proliferation and invasiveness. MTT assay results showed that TCCSUP cell proliferation was significantly decreased in the KD group compared to the NC group (Figure 2C, P<0.001). Transwell assays suggested a significant decrease in cell invasion in the KD group compared with the NC group (Figures 2D, E, P<0.001).

FIGURE 2

Figure 2 (A) The mRNA expression of CX3CL1 in TCCSUP, 5637 and T24 cells. The mRNA expression of CX3CL1 in TCCSUP cells was significantly higher than that in 5637 and T24 cells. (B) Difference in CX3CL1 mRNA levels between NC and KD TCCSUP cells. The CX3CL1 gene of TCCSUP cells in KD group was knocked down, and the knockdown rate was verified by qPCR. (C) Difference in TCCSUP cell proliferation between NC and KD groups determined by the MTT method. The results showed that cell proliferation in the KD group(Orange line) was significantly decreased. (D, E) Cell invasion of NC and KD cells determined by Transwell assays. The results showed that after CX3CL1 was knocked down, cell invasion in the KD group was significantly reduced. ***P < 0.001.

To further evaluate whether CX3CL1 expression is associated with BCa aggressiveness in vivo, we established xenografts in the shoulder of immunocompromised nude mice. Subcutaneous tumor formation in the nude mice was verified by pathologists, and the success rate was 100% (Figure 3). From the 10th day to the 35th day, the tumor volume in the KD group was significantly smaller than that in the NC group (Figure 3, days 10-35, P<0.05). The tumor weight in the KD group was significantly lower than that in the NC group (Figure 3, day 35, P=0.017).

FIGURE 3

Figure 3 Effects of CX3CL1 knockdown on (A) tumor volume and (B) weight changes in a nude mouse subcutaneous tumor model of bladder cancer. The results showed that compared with the control group, the volume and weight of subcutaneous tumors in the KD group decreased significantly. (C) Subcutaneous tumor in KD (up) and NC (down) group. (D) Hematoxylin & Eosin staining of the tumor cells. KD, knockdown; NC, negative control.

To investigate potential biological functions, we performed gene expression profiling using an mRNA microarray. Totals of 182 and 367 genes were found to be significantly upregulated and downregulated, respectively, in CX3CL1 KD TCCSUP cells. IPA-based disease and function analysis suggested that after CX3CL1 was inhibited, tumor cell proliferation (P=1.62×10^-11, z-score=-2.32), cell malignant transformation (P=6.47×10^-13, z-score=-2.289), tumor cell interphase (P=1.57×10^-12, z-score=-2.078) and differentiation of stem cells (P=4.02×10^-5, z-score=-2.465) were significantly inhibited, and tumor cell apoptosis (P=2.35×10^-7, z-score=2.236) was significantly increased (Figure 4A). CDH1, ETS1, RAF1, and EIF4E were genes with major effects in the enrichment analysis. In terms of classic pathway enrichment analysis based on IPA, the ERK/MAPK signaling pathway was significantly inhibited (P=2.19×10^-4, z-score=-1.604). ETS1, RAF1, and EIF4E were closely related to the ERK/MAPK pathway and significantly downregulated (Figure 4A and Table S2). A volcano plot and hierarchical clustering for significant difference analysis are shown in Figure 4B. Briefly, the tumor cell apoptosis (P=2.35×10^-7, z-score=2.236) pathway was significantly activated, while tumor cell proliferation (P=1.62×10^-11, z-score=-2.32), tumor cell interphase (P=1.57×10^-12, z-score=-2.078), cell viability of tumor cell lines (P=7.12×10^-14, z-score=-2.092), differentiation of stem cells (P=4.02×10^-5, z-score=-2.465), cell malignant transformation (P=6.47×10^-13, z-score=-2.289), etc. were significantly inhibited.

FIGURE 4

Figure 4 (A) Classic pathway enrichment analysis and disease and function analysis based on IPA by GeneChip PrimeView. Our experiments suggested that when CX3CL1 is inhibited, the ERK/MAPK signaling pathway is inhibited, the function of promoting tumor growth is inhibited, and tumor death is enhanced. (B) Chip volcano plot (left) and hierarchical clustering (right): in the volcano plot, the abscissa represents the multiple of the gene expression difference (transformed by log2), the ordinate represents the significant FDR of the difference (transformed by log10), the red dots were selected with |Fold Change|>1.5 and FDR<0.05 as the screening criteria for genes that were significantly differentially expressed, and the gray dots are other genes that were not significantly differentially expressed; in the hierarchical clustering data, red means that gene expression was relatively upregulated, green means that gene expression was relatively downregulated, black means that the degree of gene expression did not change significantly, and gray means that the signal intensity of the gene was not detected. (C) Relative expression of downstream genes and proteins in the ERK/MAPK signaling pathway between NC and KD cells (Western blot). NC, normal control group; KD, CX3CL1 gene knockdown group.

To further confirm whether the genes in the ERK/MAPK pathway are influenced by the expression of CX3CL1, the protein expression of the enriched genes including RAF1, EIF4E, ETS1, and CDH1 was evaluated by Western blotting. The results suggested that RAF1, EIF4E, and ETS1 protein expression was significantly downregulated by 67.2%, 52.0%, and 43.6%, respectively (all P<0.05), while CDH1 protein expression was significantly upregulated by 74.6% (P<0.05, Figure 4C).

Discussion

Tumor immunotherapy is considered as an important and hopeful curative therapy for cancers, such as PD-1/PD-L1. Immune surveillance and immune tolerance usually serve with a two-side effect throughout cancer development and progression. Scientists have found that chemokines family was the key factor in this process by affecting the immune system (3). Among different chemokines, CX3CL1 is considered as one of the most important one, and may enhance local anti-tumor immunity by polarizing immunocompetent cells around the tumor microenvironment, and therefore inhibit tumor growth (4); on the other hand, tumors spontaneously secrete chemokines and receptors, which can promote the growth of cancer cells (5). We found that high expression of CX3CL1 was significantly associated with higher risk of disease recurrence and cancer-specific death; through in vivo and in vitro functional study, we further confirmed that high expression of CX3CL1 was associated with tumor cell proliferation and invasion. Enrichment analysis and functional study based on expression array data suggested that the regulation of tumor growth by CX3CL1 might due to the enhance of the ERK/MAPK signaling pathway.

CX3CL1 is a large cytokine (373 amino acids) and the only known member of the CX3C chemokine family. It is also known as fractalkine (in humans) and neurotactin (in mice) (6, 20). It has been proved correlating with cell migration and interaction in multiple diseases. Shulby et al. presented the first evidence of the expression of CX3CR1 (the specific receptor of fractalkine) in human prostate cancer cells and the expression of fractalkine in human bone marrow endothelial cells. They also found that a neutralizing antibody against fractalkine might significantly reduce the adhesion of prostate cancer cells to human bone marrow endothelial cells, suggesting the important role of fractalkine in this process (13). By silencing CX3CR1, Jamieson et al. further defined the roles of this chemokine receptor in the adhesion of prostate cancer cells to bone marrow endothelial cells and prostate cancer cell migration towards media conditioned by different types of bone cells (21). Tardaguila et al. found that CX3CL1 expression was downregulated in HER2/neu tumors; however, paradoxically, adenovirus-mediated CX3CL1 expression in the tumor milieu enhanced mammary tumor numbers in a dose-dependent manner. Additionally, CX3CL1 triggered cell proliferation by inducing ErbB receptors through the proteolytic shedding of an ErbB ligand, and these findings support the conclusion that CX3CL1 acts as a positive modifier of breast cancer in concert with ErbB receptors (22). In another study on breast cancer, Tsang et al. reported that high CX3CL1 expression was detected in 33.3% of primary invasive cancers and that CX3CL1 expression correlated positively with increased tumor-infiltrating lymphocytes (TILs; P=0.005). Furthermore, adverse features in breast cancers, including lymph node involvement, α-β crystalline expression and high Ki67, were positively correlated with the expression of CX3CL1. Notably, in the same study, the authors showed that patients with high levels of CX3CL1 had poorer overall survival. The findings indicated an oncogenic role for CX3CL1, despite its previously suggested role in enhancing antitumor immunity, and highlighted the complicated functions of CX3CL1 in breast carcinogenesis (23). However, a study on the prognosis of patients with liver cancer showed that patients with high expression of CX3CL1 and its receptor had a better prognosis than those with low expression, with a significantly lower local recurrence rate and distant metastasis rate and longer disease-free survival and overall survival (15). These results indicate that CX3CL1 has complex functions in tumors. At present, there are few in-depth studies on the role of CX3CL1 in BCa, and the mechanism of action of the chemokine CX3CL1 is still unclear. We conducted a series of experimental studies to investigate the role of CX3CL1 in BCa.

Previous studies have shown that CX3CL1 may promote or suppress tumors in different tissues. Our research confirms the value and function of CX3CL1 in BCa. Compared with previous studies, our study identified possible signaling pathways and downstream genes of CX3CL1 in BCa through GeneChip PrimeView Human technology and explored the possible mechanisms in more depth. Through bioinformatic data analysis and statistics, the significant differences were further verified. We found that after CX3CL1 was inhibited, tumor cell proliferation, cell malignant transformation, tumor cell interphase, and differentiation of stem cells were significantly inhibited, and tumor cell apoptosis was significantly increased. These results suggest that CX3CL1 may play a role in promoting the proliferation of BCa cells. In this experiment, the ERK/MAPK pathway was significantly inhibited, and the downstream genes ETS1 (↓), RAF1 (↓), and EIF4E (↓), which are closely related to this pathway, were significantly downregulated. In addition, the downstream gene CDH1 (↑) was significantly upregulated. Further detection of the expression of the proteins encoded by the above genes was also confirmed. The RNA-Seq results, and protein expression results are consistent with a previous clinical prognostic analysis and in vivo and in vitro experimental results, which together indicate the tumor-promotive effect of CX3CL1 in BCa.

This study suggests that CX3CL1 might promote the cancer progression through the ERK/MAPK signaling pathway. Mitogen-activated protein kinase (MAPK) is a threonine/serine protein kinase that is widely present in cells. It is closely related to cell proliferation, apoptosis, tumor proliferation, invasion and metastasis. The MAPK pathway is one of the most studied signaling pathways in recent years (24, 25). Among MAPK family members, extracellular signal-regulated kinase (ERK) was the first discovered and is the most widely studied member. An abnormal ERK/MAPK pathway can lead to abnormal cell proliferation, which is related to tumorigenesis and can further affect the biological behavior of tumors (26). Current studies indicate that the ERK/MAPK signaling pathway plays an important role in tumor occurrence and development by promoting cancer cell proliferation, inhibiting tumor apoptosis, and inducing tumor neovascularization, invasion, and metastasis (27–30). Studies have also pointed out that the ERK/MAPK signaling pathway is also related to tumor resistance (31, 32). To date, many studies have found that the ERK/MAPK signaling pathway is related to many tumors, such as colon cancer, gastric cancer, liver cancer, melanoma, and BCa (33–38). According to previous research reports, among the four significantly regulated downstream genes found in this experiment, CDH1 was considered to be an anti-oncogene, and ETS1, RAF1 and EIF4E play roles in promoting cancer (39–48). Our research results are consistent with those in previous research reports.

To the best of our knowledge, this was the first study investigating the association between CX3CL1 expression and BCa. The results showed that CX3CL1 overexpression was independently and significantly associated with poor prognosis (increased risk of recurrence and cancer-specific death rate). In addition, we were able to confirm our results in TCGA dataset. This finding may be applied to clinical practice as a valuable prognostic biomarker via additional translational research approaches.

Our study has carried out a series of elaboration on the role of CX3CL1 in BCa. However, several limitations should be noted in the present study. First, the sample size was relatively small. However, we were able to validate our results in TCGA database. The results should be investigated in a large-scale study in the future. Second, the subcutaneous xenograft animal model could not fully simulate the tumor microenvironment. Although orthotopic bladder cancer model would be a better option, the current method was the most common used model in BCa research. Last but not least, only one cell line was chosen in the present study. However, the main objective of this study was to evaluate the association of CX3CL1 expression and BCa prognosis.

Conclusion

CX3CL1 was a significant and independent predictor for bladder cancer prognosis. It could promote the proliferation and invasion of bladder cancer cells via inhibiting the ERK/MAPK signaling pathway.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The studies involving human participants were reviewed and approved by The Institutional Review Board at Huashan Hospital. The written informed consent was obtained from all patients for institutional biobank. The animal study was reviewed and approved by The Institutional Review Board at Huashan Hospital.

Author Contributions

KX, RN, and NZ conceived and designed the study. GJ, HW, YW, WD, and QZ acquired the data. GJ, HW, and DH analyzed and interpreted the data. GJ, HW, and DH drafted the manuscript. KX, RN, NZ, and QD contributed to the critical revision of the manuscript. KX and RN supervised the study. All authors contributed to the article and approved the submitted version.

Funding

This research was supported by the National Natural Science Foundation of China (Grant No. 81472379 and 81402339).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The reviewer HW declared a shared affiliation, with no collaboration, with several of the authors HW, YW, WD, QZ, QD, KX to the handling editor at the time of review.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We thank all of the subjects included in this study.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2021.752860/full#supplementary-material

References

1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin (2021) 71(3):209–49. doi: 10.3322/caac.21660

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Siegel RL, Miller KD, Jemal A. Cancer Statistics, 2020. CA Cancer J Clin (2020) 70(1):7–30. doi: 10.3322/caac.21590

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Pardoll D. Does the Immune System See Tumors as Foreign or Self? Annu Rev Immunol (2003) 21(1):807–39. doi: 10.1146/annurev.immunol.21.120601.141135

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Salcedo R, Ponce ML, Young HA, Wasserman K, Ward JM, Kleinman HK, et al. Human Endothelial Cells Express CCR2 and Respond to MCP-1: Direct Role of MCP-1 in Angiogenesis and Tumor Progression. Blood (2000) 96(1):34–40. doi: 10.1182/blood.V96.1.34

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Schneider GP, Salcedo R, Welniak LA, Howard OZ, Murphy WJ. The Diverse Role of Chemokines in Tumor Progression: Prospects for Intervention. Int J Mol Med (2001) 8(3):235–44. doi: 10.3892/ijmm.8.3.235

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Bazan JF, Bacon KB, Hardiman G, Wang W, Soo K, Rossi D, et al. A New Class of Membrane-Bound Chemokine With a CX 3 C Motif. Nature (1997) 385(6617):640–4. doi: 10.1038/385640a0

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Fong AM, Erickson HP, Zachariah JP, Poon S, Schamberg NJ, Imai T, et al. Ultrastructure and Function of the Fractalkine Mucin Domain in CX(3)C Chemokine Domain Presentation. J Biol Chem (2000) 275(6):3781–6. doi: 10.1074/jbc.275.6.3781

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Haskell CA, Cleary MD, Charo IF. Molecular Uncoupling of Fractalkine-Mediated Cell Adhesion and Signal Transduction. Rapid Flow Arrest of CX3CR1-Expressing Cells is Independent of G-Protein Activation. J Biol Chem (1999) 274(15):10053–8. doi: 10.1074/jbc.274.15.10053

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Andre F, Cabioglu N, Assi H, Sabourin JC, Delaloge S, Sahin A, et al. Expression of Chemokine Receptors Predicts the Site of Metastatic Relapse in Patients With Axillary Node Positive Primary Breast Cancer. Ann Oncol (2006) 17(6):945–51. doi: 10.1093/annonc/mdl053

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Celesti G, Di Caro G, Bianchi P, Grizzi F, Marchesi F, Basso G, et al. Early Expression of the Fractalkine Receptor CX3CR1 in Pancreatic Carcinogenesis. Br J Cancer (2013) 109(9):2424–33. doi: 10.1038/bjc.2013.565

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Kim M, Rooper L, Xie J, Kajdacsy-Balla AA, Barbolina MV. Fractalkine Receptor CX3CR1 is Expressed in Epithelial Ovarian Carcinoma Cells and Required for Motility and Adhesion to Peritoneal Mesothelial Cells. Mol Cancer Res (2012) 10(1):11–24. doi: 10.1158/1541-7786.MCR-11-0256

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Yao X, Qi L, Chen X, Du J, Zhang Z, Liu S. Expression of CX3CR1 Associates With Cellular Migration, Metastasis, and Prognosis in Human Clear Cell Renal Cell Carcinoma. Urol Oncol (2014) 32(2):162–70. doi: 10.1016/j.urolonc.2012.12.006

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Shulby SA, Dolloff NG, Stearns ME, Meucci O, Fatatis A. CX3CR1-Fractalkine Expression Regulates Cellular Mechanisms Involved in Adhesion, Migration, and Survival of Human Prostate Cancer Cells. Cancer Res (2004) 64(14):4693–8. doi: 10.1158/0008-5472.CAN-03-3437

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Ohta M, Tanaka F, Yamaguchi H, Sadanaga N, Inoue H, Mori M. The High Expression of Fractalkine Results in a Better Prognosis for Colorectal Cancer Patients. Int J Oncol (2005) 26(1):41–7. doi: 10.3892/ijo.26.1.41

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Matsubara T, Ono T, Yamanoi A, Tachibana M, Nagasue N. Fractalkine-CX3CR1 Axis Regulates Tumor Cell Cycle and Deteriorates Prognosis After Radical Resection for Hepatocellular Carcinoma. J Surg Oncol (2007) 95(3):241–9. doi: 10.1002/jso.20642

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Park MH, Lee JS, Yoon JH. High Expression of CX3CL1 by Tumor Cells Correlates With a Good Prognosis and Increased Tumor-Infiltrating CD8+ T Cells, Natural Killer Cells, and Dendritic Cells in Breast Carcinoma. J Surg Oncol (2012) 106(4):386–92. doi: 10.1002/jso.23095

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Tang J, Chen Y, Cui R, Li D, Xiao L, Lin P, et al. Upregulation of Fractalkine Contributes to the Proliferative Response of Prostate Cancer Cells to Hypoxia via Promoting the G1/S Phase Transition. Mol Med Rep (2015) 12(6):7907–14. doi: 10.3892/mmr.2015.4438

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Eble JN, Sauter G, Epstein JI, Sesterhenn IA. Pathology and Genetics of Tumours of the Urinary System and Male Genital Organs: World Health Organization Classification of Tumours. Int Agency Res Cancer (2004) 7:89–154.

Google Scholar

19. Edge SB, Compton CC. The American Joint Committee on Cancer: The 7th Edition of the AJCC Cancer Staging Manual and the Future of TNM. Ann Surg Oncol (2010) 17(6):1471–4. doi: 10.1245/s10434-010-0985-4

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Pan Y, Lloyd C, Zhou H, Dolich S, Deeds J, Gonzalo J-A, et al. Neurotactin, a Membrane-Anchored Chemokine Upregulated in Brain Inflammation. Nature (1997) 387(6633):611–7. doi: 10.1038/42491

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Jamieson WL, Shimizu S, D'Ambrosio JA, Meucci O, Fatatis A. CX3CR1 is Expressed by Prostate Epithelial Cells and Androgens Regulate the Levels of CX3CL1/fractalkine in the Bone Marrow: Potential Role in Prostate Cancer Bone Tropism. Cancer Res (2008) 68(6):1715–22. doi: 10.1158/0008-5472.CAN-07-1315

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Tardáguila M, Mira E, García-Cabezas MA, Feijoo AM, Quintela-Fandino M, Azcoitia I, et al. CX3CL1 Promotes Breast Cancer via Transactivation of the EGF Pathway. Cancer Res (2013) 73(14):4461–73. doi: 10.1158/0008-5472.CAN-12-3828

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Tsang JYS, Ni Y-B, Chan S-K, Shao M-M, Kwok Y-K, Chan K-W, et al. CX3CL1 Expression is Associated With Poor Outcome in Breast Cancer Patients. Breast Cancer Res Treat (2013) 140(3):495–504. doi: 10.1007/s10549-013-2653-4

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Shen H, Xu W, Luo W, Zhou L, Yong W, Chen F, et al. Upregulation of Mdr1 Gene is Related to Activation of the MAPK/ERK Signal Transduction Pathway and YB-1 Nuclear Translocation in B-Cell Lymphoma. Exp Hematol (2011) 39(5):558–69. doi: 10.1016/j.exphem.2011.01.013

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Parakh S, Murphy C, Lau D, Cebon JS, Andrews MC. Response to MAPK Pathway Inhibitors in BRAF V600M-Mutated Metastatic Melanoma. J Clin Pharm Ther (2014) 40(1):121–3. doi: 10.1111/jcpt.12229

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Aebersold DM, Burri P, Beer KT, Laissue J, Djonov V, Greiner RH, et al. Expression of Hypoxia-Inducible Factor-1α: A Novel Predictive and Prognostic Parameter in the Radiotherapy of Oropharyngeal Cancer. Cancer Res (2001) 61(7):2911–6.

PubMed Abstract | Google Scholar

27. Wang YF, Jiang CC, Kiejda KA, Gillespie S, Zhang XD, Hersey P. Apoptosis Induction in Human Melanoma Cells by Inhibition of MEK is Caspase-Independent and Mediated by the Bcl-2 Family Members PUMA, Bim, and Mcl-1. Clin Cancer Res (2007) 13(16):4934–42. doi: 10.1158/1078-0432.CCR-07-0665

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Chen Y-J, Tsai Y-M, Kuo C-D, Ku K-L, Shie H-S, Liao H-F. Norcantharidin is a Small-Molecule Synthetic Compound With Anti-Angiogenesis Effect. Life Sci (2009) 85(17-18):642–51. doi: 10.1016/j.lfs.2009.09.003

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Yue P, Gao Z-H, Xue X, Cui S-X, Zhao C-R, Yuan YI, et al. Des-γ-Carboxyl Prothrombin Induces Matrix Metalloproteinase Activity in Hepatocellular Carcinoma Cells by Involving the ERK1/2 MAPK Signalling Pathway. Eur J Cancer (2011) 47(7):1115–24. doi: 10.1016/j.ejca.2011.01.017

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Corona G, Deiana M, Incani A, Vauzour D, Dessì MA, Spencer JPE. Hydroxytyrosol Inhibits the Proliferation of Human Colon Adenocarcinoma Cells Through Inhibition of ERK1/2 and Cyclin D1. Mol Nutr Food Res (2009) 53(7):897–903. doi: 10.1002/mnfr.200800269

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Ochi N, Takigawa N, Harada D, Yasugi M, Ichihara E, Hotta K, et al. Src Mediates ERK Reactivation in Gefitinib Resistance in non-Small Cell Lung Cancer. Exp Cell Res (2014) 322(1):168–77. doi: 10.1016/j.yexcr.2014.01.007

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Fan D-P, Zhang Y-M, Hu X-C, Li J-J, Zhang W. Activation of AKT/ERK Confers non-Small Cell Lung Cancer Cells Resistance to Vinorelbine. Int J Clin Exp Pathol (2014) 7(1):134.

PubMed Abstract | Google Scholar

33. McCubrey JA, Steelman LS, Abrams SL, Chappell WH, Russo S, Ove R, et al. Emerging MEK Inhibitors. Expert Opin emerging Drugs (2010) 15(2):203–23. doi: 10.1517/14728210903282760

CrossRef Full Text | Google Scholar

34. Simeonova PP, Wang S, Hulderman T, Luster MI. C-Src-Dependent Activation of the Epidermal Growth Factor Receptor and Mitogen-Activated Protein Kinase Pathway by Arsenic Role in Carcinogenesis. J Biol Chem (2002) 277(4):2945–50. doi: 10.1074/jbc.M109136200

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Drobna Z, Jaspers I, Thomas DJ, Stýblo M. Differential Activation of AP-1 in Human Bladder Epithelial Cells by Inorganic and Methylated Arsenicals. FASEB J (2003) 17(1):67–9. doi: 10.1096/fj.02-0287fje

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Fujimori Y, Inokuchi M, Takagi Y, Kato K, Kojima K, Sugihara K. Prognostic Value of RKIP and P-ERK in Gastric Cancer. J Exp Clin Cancer Res (2012) 31(1):30. doi: 10.1186/1756-9966-31-30

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Chen C-C, Chen C-Y, Wang S-H, Yeh C-T, Su S-C, Ueng S-H, et al. Melatonin Sensitizes Hepatocellular Carcinoma Cells to Chemotherapy Through Long non-Coding RNA RAD51-AS1-Mediated Suppression of DNA Repair. Cancers (Basel) (2018) 10(9):320. doi: 10.3390/cancers10090320

CrossRef Full Text | Google Scholar

38. Gupta A, Hossain MM, Miller N, Kerin M, Callagy G, Gupta S. NCOA3 Coactivator is a Transcriptional Target of XBP1 and Regulates PERK–Eif2α–ATF4 Signalling in Breast Cancer. Oncogene (2016) 35(45):5860–71. doi: 10.1038/onc.2016.121

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Giroldi LA, Schalken JA. Decreased Expression of the Intercellular Adhesion Molecule E-Cadherin in Prostate Cancer: Biological Significance and Clinical Implications. Cancer Metastasis Rev (1993) 12(1):29–37. doi: 10.1007/BF00689788

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Giroldi LA, Bringuier PP, Schalken JA. Defective E-Cadherin Function in Urological Cancers: Clinical Implications and Molecular Mechanisms. Invasion Metastasis (1994) 14(1-6):71–81.

PubMed Abstract | Google Scholar

41. Nigam AK, Savage FJ, Boulos PB, Stamp GW, Liu D, Pignatelli M. Loss of Cell-Cell and Cell-Matrix Adhesion Molecules in Colorectal Cancer. Br J Cancer (1993) 68(3):507–14. doi: 10.1038/bjc.1993.377

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Chen M, Cooper JA. Ser-3 is Important for Regulating Mos Interaction With and Stimulation of Mitogen-Activated Protein Kinase Kinase. Mol Cell Biol (1995) 15(9):4727–34. doi: 10.1128/MCB.15.9.4727

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Dent P, Reardon DB, Morrison DK, Sturgill TW. Regulation of Raf-1 and Raf-1 Mutants by Ras-Dependent and Ras-Independent Mechanisms In Vitro. Mol Cell Biol (1995) 15(8):4125–35. doi: 10.1128/MCB.15.8.4125

PubMed Abstract | CrossRef Full Text | Google Scholar

44. von der Haar T, Gross JD, Wagner G, McCarthy JEG. The mRNA Cap-Binding Protein Eif4e in Post-Transcriptional Gene Expression. Nat Struct Mol Biol (2004) 11(6):503–11. doi: 10.1038/nsmb779

PubMed Abstract | CrossRef Full Text | Google Scholar

45. De Benedetti A, Graff JR. eIF-4E Expression and its Role in Malignancies and Metastases. Oncogene (2004) 23(18):3189–99. doi: 10.1038/sj.onc.1207545

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Oettgen P. The Role of Ets Factors in Tumor Angiogenesis. J Oncol (2010) 2010:767384. doi: 10.1155/2010/767384

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Ahn KS, Bae E, Han D, Choi HY. 211 Increasing Expression Level of ETS-1 via HIF-1a Induces the Angiogenesis and Invasiveness Human Bladder Cancer: It may be One of Useful Molecular Markers to Determine the Progression of Human Bladder Cancer. Eur Urol Suppl (2004) 3(2):55. doi: 10.1016/S1569-9056(04)90212-6

CrossRef Full Text | Google Scholar

48. Oikawa M, Abe M, Kurosawa H, Hida W, Shirato K, Sato Y. Hypoxia Induces Transcription Factor ETS-1 via the Activity of Hypoxia-Inducible Factor-1. Biochem Biophys Res Commun (2001) 289(1):39–43. doi: 10.1006/bbrc.2001.5927

PubMed Abstract | CrossRef Full Text | Google Scholar