Prognostic value of CC-chemokine receptor seven expression in patients with metastatic renal cell carcinoma treated with tyrosine kinase inhibitor

Background CC-chemokine receptor seven (CCR7), a G-protein coupled receptor normally facilitating immune cells lymphatic homing, has recently been identified on several cancer cells in promoting invasion and lymphatic specific metastasis by mimicking normal leukocytes. As tyrosine kinase inhibitors for metastatic renal cell carcinoma (mRCC) mostly emphasized on vascular inhibition, whether the CCR7 expressing tumor cells with potential lymphatic invasion function could have an impact on mRCC patient’s drug response and survival, was unknown. Methods In this study, in a clinical aspect, we retrospectively investigated the prognostic and predictive impact of tumoral CCR7 expression in 110 mRCC patients treated with sunitinib and sorafenib, and its correlation with pre- or post-administration lymphatic involvement. Immunohistochemistry on tissue microarrays were conducted for CCR7 expression evaluation. Results Kaplan-Meier and univariate analyses suggested high tumoral CCR7 expression as an adverse prognosticator for mRCC patients’ overall survival (OS), which was further confirmed in the multivariate analyses (P = 0.002, P = 0.003 for bootstrap). This molecule could be combined with Heng’s risk model for better patient OS prediction. High tumoral CCR7 expression was also an independent dismal predictor for patients’ progression free survival (PFS) (P = 0.010, P = 0.013 for bootstrap), and correlated with poorer best drug response. Moreover, a possible correlation of CCR7 high expression and patients’ baseline and post-administration lymph node metastasis was found. Conclusions High tumoral CCR7 expression correlated with potential lymphatic involvement and poor prognosis of mRCC patients treated with tyrosine kinase inhibitors. Further external validations and basic researches were needed to confirm these results. Electronic supplementary material The online version of this article (doi:10.1186/s12885-017-3065-3) contains supplementary material, which is available to authorized users.


Background
For patients with metastatic renal cell carcinoma (mRCC), therapeutic options have expanded significantly these years, since vascular endothelial growth factor (VEGF)-targeted tyrosine kinase inhibitor (TKI) drugs such as sunitinib and pazopanib have been established as first-line therapy [1]. Several clinical based prognostic models, for example the Heng's risk criteria, have also achieved remarkable progress in mRCC patient survival prediction [2]. However, the objective response rates (ORRs) for most first-line TKI drugs were only around 30% [3], and the lack of validated molecular biomarkers impeded their personalized approach [4]. This was in contrast to many other tumor types, in which protein expression and mutation were used as basic accesses for drugs response and patient survival prediction [5,6].
CC-chemokine receptor 7 (CCR7), a G-protein coupled receptor (GPCR) mostly expressed on immune cells, was initially regarded as an important regulator facilitating leukocytes homing to the lymphatic structures, where its two ligands CC-chemokine ligand 19 (CCL19) and CCL21 are constitutively expressed [7]. However, in recent years, aberrant high CCR7 expression has also been identified on several tumor types, linking to a potential invasive phenotype [8]. It has been suggested that CCR7 positive tumor cells could mimic the normal lymphocyte homing function and interact with lymph vessels, leading to subsequent lymph node specific metastasis [9].
Lymphatic and hematogenous disseminations were two regular metastasis pathways for malignancy. For mRCC, although the most common metastatic site was the lung, possibly via a hematogenous approach, local or distant lymph node involvement at diagnosis was not rare [10]. Patient receiving TKIs could also develop new lymph node metastasis during the treatment, leading to a progressive disease (PD). Several theories of TKI drug resistance emphasized an increase of tumor cell invasiveness after drug administration [11]. These processes mostly accompanied with tumor cell migration and matrix metalloproteinase-9 (MMP-9) mediated extracellular matrix degradation [12], which was similar to the CCR7 mediated lymph vessel intravasation process [8]. As VEGF targeted therapies mostly focused on blood vessels inhibition, whether they could enhance the possibility of mRCC metastasis through other pathways, such as CCR7 mediated lymph vessel invasion and therapy resistance, was not known.
Thus, here through immunohistochemistry (IHC), we retrospectively evaluated the CCR7 expression in 110 primary tumor specimens of mRCC patients treated with sunitinib and sorafenib. The result suggested a positive correlation of CCR7 expression with patient baseline lymph node metastasis and TKI drugs response. CCR7 high expression could predict a poorer overall survival (OS) and progression free survival (PFS) for mRCC patients after TKIs.

Patient selection
We initially screened a total of 138 mRCC patients treated with sunitinib or sorafenib between March 2005 and June 2014 at the Department of Urology, Zhongshan Hospital, Fudan University. The inclusion criteria were: pathologically proven RCC patient with metastatic lesion, treated with sunitinib or sorafenib at first without further second-line treatment, had enough Formalin Fixed Paraffin Embedded (FFPE) specimens, and had detailed laboratorial, imaging and survival information. Patients who had former malignant history, received metastasectomy or those with tumor necrosis area >80%, unavailable data were excluded. At last, 28 patients were excluded and 110 patients were selected for the study, in which three were excluded from PFS analysis for incomplete drug response information. This study was approved by the Clinical Research Ethics Committee of Zhongshan Hospital, Fudan University (Shanghai, China) (B2015-030).
Patients' OS was defined as the time from therapy initiation to the time of death, or was censored at the last follow-up. PFS was calculated from the time of therapy initiation to the time of progression, according to the RECIST 1.1 criteria [13], or was censored at the last follow-up. All data were collected retrospectively from medical records and electronic databases using uniform database templates, and the last follow-up time was December 2015. According to the 2014 EAU guidelines [14] and 2012 ISUP consensus [15], two urologic pathologists (Yuan J. and Jun H.) independently reviewed all the H & E slides of patient samples and confirmed the RCC diagnosis and Fuhrman grade classification. Initial stage at diagnosis was reclassified based on the 2010 AJCC TNM classification [10]. Heng's risk model was applied as previously reported [2].

Tissue microarray and immunohistochemistry
Two representative tumor cores 3 mm in diameter from each sample were selected for tissue micro array (TMA) construction. Anti-CCR7 antibody (ab32527, Abcam, diluted 1/1000) and Dako EnVision Detection System were applied in the immunohistochemistry procedure [16]. Through western blot in RCC cell lines, the specificity of antibody was confirmed. Negative control was performed without applying primary antibody. Olympus CDD camera, Nikon eclipse Ti-s microscope (×200 magnification) and NIS-Elements F3.2 software were used to record the staining results. Using Image-Pro Plus version 6.0 software (Media Cybernetics Inc., USA), an integrated optical density (IOD) score could be calculated for each scan. Two urologists unaware of the patients' clinical data evaluated these slides. Each person took three independent shots with the strongest staining for each core, and the IOD mean of each patient's two cores (six scans) were calculated. Kappa value was calculated for evaluating inter-observer agreement.

Statistical analysis
Univariate analysis was carried out to explore the prognostic and predictive value of continuous CCR7 IOD score. The smooth estimates of hazard ratio (HR) of IOD on patient survival were displayed using R software, "phenoTest" package [17]. For clinical usage, we dichotomized the IOD into high/low expression through minimum p value method using x-tile software [18], and because the p values obtained might be overestimated, they were corrected using the formula proposed by Altman and colleagues [19]. The smooth HR curves after dichotomizing were illustrated through R software, "smoothHR" package [20]. After this, χ2 test, Fisher's exact method and Cochran-Mantel-Haenszel χ2 test were applied for assessing the relationship between CCR7 expression and patients' clinicopathological parameters. Kaplan-Meier, univariate and subsequent multivariate analysis were performed, in which 1000 bootstrap resamples was used for reducing overfitting bias. Finally, time dependent receiver operating characteristic (ROC) analysis was done to analyze the adding prognostic value of CCR7 expression to the Heng's risk model. GraphPad Prism 6 (GraphPad Software Inc., USA), SPSS 21.0 (SPSS Inc., USA), X-tile 3.6.1 (Robert L Camp, USA) and R software 3.1.2 (R Foundation for Statistical Computing, Austria) were used in these procedures. A two-sided P-value < 0.05 was regarded as statistically significant.

CCR7 staining and cut off point choosing
CCR7 expression in the RCC sample was variable, and mostly on the membrane and cytoplasm of tumor cells (Additional file 1: Figure S1A and B). Its expression in peritumoral tissue was relatively low (Additional file 1: Figure S1C). Inter-observer agreement was acceptable according to the kappa value 0.745, thus the CCR7 IOD means from the two observers were again averaged as the final IOD, and the distributions were 73-584 for range; 238 ± 97 for mean and SD; 223 (173-287) for median and IQR.
In order to select an appropriate cut off value for clinical usage, first we performed a univariate analysis using the continuous CCR7 IOD score (Additional file 2: Table S1). Result suggested CCR7 as a significant adverse prognostic marker for mRCC patients' OS and PFS (P < 0.001 for  Figure  S1D and F). Then, through minimum p value method using log rank test, IOD = 215 was chosen as the cut off point (P < 0.001 for both), and the p values were still significant after being corrected (OS, P < 0.001, PFS, P = 0.001). The smooth HR curve also displayed significant and stable prognostic differences after dichotomizing (Additional file 1: Figure S1E and G).

Patient baseline characteristics and its association with dichotomous CCR7 expression
The 110 patients' baseline clinical characteristics were shown in Table 1. All patients were East Asian. The medium age at TKIs initiation was 59 years old (range 14-78). All patients have received radical, partial or cytoreductive nephrectomy. Results showed that the dichotomous CCR7 expression associated with TNM stage at initial diagnosis (P = 0.004), lymph node involvement (P = 0.031) and marginally with histology (P = 0.086) and number of metastatic sites (P = 0.070). It also correlated with patients' best response of TKIs (P < 0.001).

Impact of baseline characteristics, including dichotomous CCR7 expression, on OS in mRCC patients receiving TKIs
Within this cohort, 64.5% (71/110) patients died during the follow-up and the median OS was 23.5 months. Kaplan-Meier analysis revealed that patients with high CCR7 expression had a significantly poorer OS (P < 0.001 after correction) (Fig. 1a). Univariate analysis confirmed this significance (P < 0.001) (Additional file 2: Table S1). Multivariate Cox analysis was further performed and suggested high CCR7 expression as an independent adverse prognostic factor for mRCC patients' OS prediction (HR 2.256, 95% CI 1.336-3.809, P = 0.002; P = 0.003 after 1000 bootstrap), together with tumor histology and Heng's risk group (Table 2). Lung and lymph node involvement were excluded from the model for being potential confounding factors for number of metastatic sites. Stratified analysis were further performed, and we found that CCR7 expression could discriminate most patient groups' overall survival except those in the nonclear cell type or Heng's risk favorable/poor groups (Additional file 3: Table S2). But after incorporating the high/low CCR7 expression directly into Heng's risk to form a new model, the OS between different groups displayed vigorously discriminative consequences (P < 0.001) (Fig. 1b). Furthermore, ROC analysis was carried out at the time of 12 and 24-month follow-up, and the new model showed better prognostic power than using Heng's risk model alone in both the ccRCC and all patient groups (Fig. 2). Incorporating CCR7 IOD score as a continuous variable also displayed similar results (Additional file 4: Figure S2).

Impact of baseline characteristics, including dichotomous CCR7 expression, on PFS in mRCC patients receiving TKIs
During the follow-up period, 85.0% (91/107) patients have developed disease progression. The median PFS was 9.8 months. Patients' best response and its correlation with CCR7 were shown in Table 1. Figure 3a revealed that RCC in the PD group displayed a significantly higher CCR7 expression compared to partial response (PR) and stable disease (SD) groups. Kaplan-Meier analysis suggested an adverse predictive effect of high CCR7 expression in patients receiving TKIs (P = 0.001 after correction) (Fig. 3b) and was also confirmed in a multivariate model (HR 1.835, 95% CI 1.156-2.912, P = 0.010; P = 0.013 after 1000 bootstrap) ( Table 2). After incorporating CCR7 into the Heng's model, patients in the new model  Bootstrapping with 1000 resamples were used displayed significant PFS divergence between different groups (P < 0.001) (Fig. 3c). Since the Heng's risk criteria was initially designed for OS prediction, further ROC comparison was not performed.

CCR7 expression and its correlation with lymph node involvement
In Fig. 4a, the CCR7 IOD score of mRCC patients with different baseline metastatic sites were plotted, and revealed a potential higher expression of CCR7 in patients with baseline lymph node metastasis, in accordance with the χ2 test in Table 1, though the Kruskal-Wallis test did not meet statistical significance (P = 0.083). For exploring the possible impact of CCR7 on lymphatic invasion during the drug treatment period, we found that four patients within this cohort have developed disease progression due to new lymph node lesions development, and all their tumor samples displayed CCR7 high expression (Fig. 4b).
Discussion CCR7 was naturally a homeostatic chemokine receptor expressed on various subtypes of immune cells encompassing T cells, B cells, natural killer cells and dendritic cells, enabling them to circulate through the CCL-19/21 positive lymphatic highways [7]. Its expression on cancer cells was first recognized on hematogenous malignancies, in which the connection between high CCR7 expression and lymphoid organ involvement was discovered [21]. Subsequently, numerous laboratory studies confirmed this pro-invasion, mostly pro-lymph vessel metastatic function in various cancer types, including breast cancer, melanoma, non-small cell lung cancer, prostate cancer, headand-neck cancer and gastrointestinal cancer [9]. Here, our results suggested that CCR7 also expressed on several RCC specimens and associated with various patient baseline characteristics (Table 1, Additional file 1: Figure S1). The clinical prognostic value of CCR7 has been studied in several other cancer types. In a meta-analysis including 1697 gastric cancer patients, high CCR7 correlated with a worse 5-year overall survival rate [22]. Studies in melanoma and colorectal cancer also displayed similar results [23,24]. In this study, we have found that high CCR7 staining intensity could be used as an adverse prognosticator for mRCC patients' overall survival. Moreover, a new model integrating CCR7 expression into Heng's risk criteria performed better than using Heng's risk alone, and the adding prognostic value mostly came from the intermediate patient groups (Additional file 3: Table S2, Fig. 2a and b). This indicated that a substantial number of patients might switch between risk groups with a consequence for choice of treatment strategy after CCR7 status was considered. ROC analysis for the ccRCC patient group were also carried out, considering the poor performance of CCR7 in non-ccRCC patients ( Fig. 2c and d).
For PFS analysis, to our knowledge, this study was the first to report an association between tumoral CCR7 expression and PFS in mRCC patients following TKIs. As the most widely used systemic therapy at present, TKI drugs only reached ORRs for about 30% for mRCC patients, and the theories of drug resistance was complicated [25]. It is increasingly evident that in some tumors, in which angiogenesis is thwarted genetically or pharmacologically, cancer cells could adapt by migrating more aggressively into normal tissue, based on several pre-existing invasion programmes such as epithelial-mesenchymal transition (EMT) and MMP2/MMP9 secretion [26,27], or by switching on several distinctive programmes which were currently unknown [28]. Since the CCR7 mediated lymphatic specific migration and metastasis was also based on several above mentioned programmes [8,29], and as a molecule which could give tumor itself survival signal besides [30], we hypothesized that CCR7 positive RCC cells might have a potential to migrate into the adjacent lymphatic tissue for survival after TKI drugs administration, subsequently leading to resistance and disease progression. As a result, our study did identify an independent significant adverse predictive effect of high CCR7 expression on mRCC patients' PFS (Fig. 3b) and its association with poorer best drug response (Table 1, Fig. 3a). Patients' baseline lymph node involvement status was also correlated ( Table 1, Fig. 4a) and four patients who developed PD due to new lymph node metastasis all represented high CCR7 staining for their tumor tissues. All these results above indicated an impact of tumoral CCR7 expression on patients' PFS and lymphatic involvement status, making this molecule a potential predictor for mRCC patients receiving TKIs.
The major limitations of this pilot study were its retrospective design and relatively small sample size, with patients from a single center and same ethnic background. Although central patient data review and bootstrap validation were performed for minimizing inter-observer and over-fitting bias, further external prospective validations were required, and related basic researches were Fig. 3 Impact of tumoral CCR7 expression on patients' best drug response for tyrosine kinase inhibitors and PFS. a Patients' best drug response according to tumoral CCR7 expression; b PFS according to tumoral CCR7 expression; c Heng's risk model expanded with tumoral CCR7 expression needed for verifying our hypothesis. Samples from the metastatic site after TKI therapies might be excellent candidates for in-depth investigation. Besides, although we took two cores and six scans from each patient's sample, intratumoral heterogeneity still might confound the results. Patients in this study received sunitinib or sorafenib as first-line therapy because other agents were not available in China at that time, and other VEGFbased TKI drugs such as pazopanib, axitinib and cabozantinib should also be taken into further consideration.

Conclusions
Our study indicated that high tumoral CCR7 expression correlated with potential lymphatic involvement and poor prognosis in mRCC patients treated with TKIs. This biomarker could also be combined with the Heng's risk model for better risk stratification.

Additional files
Additional file 1: Figure S1.  Availability of data and materials The basic patient information and IHC staining results have been shown in Table 1 and Additional file 1: Figure S1. And we are sorry that the detailed dataset supporting the conclusions could not be provided at present, because the cohort information was currently under an update. We are now adding more mRCC patients (diagnosed from 2014 to 2016) into the cohort and refreshing the follow up information (till 2016). Once this is done, we would like to combine CCR7 with other biomarkers for better mRCC patient survival prediction, and at that time those updated information would be available to the readers.
Authors' contributions YX carried out and conducted experiments, performed statistical analysis and drafted the manuscript. LL participated in the collection of patient materials and drafting of the manuscript. YX carried out laboratory work and data analysis. JW and QB performed laboratory work and participated in the correction of words in the manuscript. WX and YQ participated in the study design and collection of related articles. JX took charge of the study design and revising manuscript critically for important intellectual content. JG conceived of the study, and led the data analysis and oversaw the drafting of the manuscript. All authors read and approved the final manuscript.