Additive value of pre-operative and one-month post-operative lymphocyte count for death-risk stratification in patients with resectable pancreatic cancer: a multicentric study

Background Pancreatic adenocarcinoma (PDAC) incidence is increasing worldwide. Several studies have shown that lymphopenia was correlated with a poor prognosis but the potential interest to measure lymphopenia in the pre and post-operative setting as well as its added value among conventional prognostic factors was never investigated. Methods Data from two independent cohorts in whom patients underwent resection for pancreatic carcinoma were retrospectively recorded. We examined the association between perioperative findings, pre and post-operative lymphocyte counts and overall survival (OS) in univariate and multivariate analyses. Performance assessment and internal validation of the final model were evaluated with Harrell’s C-index, calibration plot and bootstrap sample procedures. Results Three hundred ninety patients were included in the analysis between 2000 and 2011. Pre and post-operative lymphocyte counts were independent prognostic factors associated with OS in multivariate analysis (p = 0.0128 and p = 0.0764, respectively). The addition of lymphocyte count variable to the conventional parameters identified in multivariate analysis (metastatic lymph node ratio, veinous emboli and adjuvant chemotherapy) significantly improved the model discrimination capacity (bootstrap mean difference = 0.04; 95 % CI, 0.01–0.06). The use of a threshold and combining the categorical (≥1000; <1000) information in pre and post lymphocyte counts permitted the identification of 4 subgroups of patients with different prognosis (p < 0.0001). Finally, the description of patients in long-term remission showed that only 3 of 65 (4.6 %) patients with post-operative lymphocyte count under 1000/mm3 were alive 4 years after surgery contrary to 54 of 236 (22.8 %) patients with a post-operative lymphocyte count above 1000/mm3. Conclusion Pre and post-operative lymphopenia are independent prognostic factors for OS and they have an additive value regarding conventional prognostic factors for death-risk stratification and to predict long-term survival. Lymphopenia should be included as stratification factors in future clinical trial assessing overall survival in pancreatic cancer patients. Electronic supplementary material The online version of this article (doi:10.1186/s12885-016-2860-6) contains supplementary material, which is available to authorized users.


Background
Pancreatic ductal adenocarcinoma incidence and mortality is increasing worldwide and is one of the leading causes of cancer-related death [1][2][3]. Surgical pancreatic resection is the only curative treatment. However, patients treated by surgery and adjuvant chemotherapy only achieved 10-20 % of long-term survival [4,5]. Then, the identification of prognostic factors correlated with the risk of early relapse is an important issue to better optimize neoadjuvant or adjuvant treatments.
Current recognized prognostic factors include mainly pathological parameters such as lymph node invasion, tumour grade and resection margin involvement [6]. However, these conventional clinical parameters as well as gene mutation status are insufficient to predict accurately deathrisk and the probability of long-term remission [7,8].
A growing body of evidence suggests that immunerelated biomarkers are correlated with survival in several cancers. A high infiltration of tumours by CD8+ T cell lymphocytes (TILs) was correlated to a good prognosis in colorectal and pancreatic cancers [9][10][11]. Conversely, cancers endowed with the ability to escape the immune system are expected to display a worse prognosis. Inflammation, recruitment of suppressive immune cells (regulatory T lymphocytes or myeloid derived suppressor cells) or induction of lymphocyte-apoptosis mediated by tumour cells, are potential mechanisms leading to immune escape [12]. Peripheral blood lymphocyte count to assess lymphopenia in cancer patients might be a convenient and clinically relevant option to identify cancers associated with an enhanced risk of tumour immune evasion and poor prognosis.
Several studies have shown that an elevated preoperative neutrophil-to-lymphocyte ratio and a decreased preoperative lymphocyte count allow the prediction of chemotherapy-related toxicities [13,14]. Moreover, preoperative lymphopenia might be a prognostic factor for several cancer patients [15][16][17]. In pancreatic cancer, immune suppression at baseline is also correlated with overall survival. However, despite the number of studies conducted so far, lymphocyte count parameter is not used in current practice, probably because its additive value regarding conventional prognostic markers is still not widely known [18][19][20][21][22]. Furthermore, recent advances in immunology evidenced that chemotherapy might promote an immunogenic cell death leading to an increased anti-tumour immunity [23]. The presence of a lymphopenia after surgery might impact the efficacy of adjuvant chemotherapy. Having observed in current clinical practice that most patients with pre-operative lymphopenia have a lymphocyte count above 1000/mm3 1-month post surgery, we hypothesized that the prognostic value of lymphopenia might be more accurate when assessed after pancreatic adenocarcinoma removal.
Consequently, we decided to conduct a study based on two independent cohorts including 390 patients treated by surgery for a localized pancreatic adenocarcinoma. The aims of this study were: i) to conduct a confirmatory study to validate the prognosis value of pre-operative lymphopenia, ii) to assess the impact of post-operative lymphopenia on overall survival, iii) to discriminate the additive value of pre and post-operative lymphocyte counts among conventional prognosis factors.

Population
This retrospective analysis included data from two independent cohorts of patients with histologically confirmed pancreatic adenocarcinoma, who underwent a curative intent surgery. The first cohort included patients treated in a university hospital (Georges Pompidou European Hospital, henceforth referred as HEGP) and in a regional cancer institute (university hospital of Besançon, Belfort and Montbeliard general hospital) between January 1, 2000 and April 30, 2010. The second cohort included patients treated in the university hospital of Strasbourg (Hautepierre) between January 1, 2004 and December 31, 2011. We excluded patients with other histopathological type of cancers: cholangiocarcinoma, neuroendocrine tumor and patients for whom preoperative lymphocyte count parameter was not available. Data from the two cohorts were updated in june, 2013.

Data collection
The following data were collected at diagnosis for each patient: center and patient identification, age, sex, diabetes, level of albuminemia which is the presence of albumin in blood, pathological characteristics and prescription of adjuvant chemotherapy. Lymphocyte counts on routine blood tests were recorded the day before the surgery and 1 month after the surgery. Deaths were collected during the follow-up of the study for each patients. All data from the two cohorts were extracted from the paper files and the desktop folder for each patient. The software was ChimioProd and Dxcare, Millenium, Axigate and Dxcare for HEGP, Belfort and Montbeliard, Besançon and Strasbourg Hospitals respectively. Adjuvant chemotherapy data from Belfort, Montbeliard and Besançon hospitals were excerpted from a regional register (BPC software).

Statistical analysis
We provided the mean (SD) values and frequency (percentages) for the description of continuous and categorical variables, respectively. The means and the proportion were compared using Student's t test and the chi-squared test (or Fisher's exact test, if appropriate), respectively.
Due to the skewed lymphocyte count distribution we used for its description the median, and the interquartile range for the dispersion measurement, as recommended [24]. Wilcoxon rank-sum test was performed for lymphocyte count distribution comparison among the cohort set.
Overall Survival (OS) was calculated from the date of surgery to the date of death from any cause. Alive patients were censored at the last follow-up. OS was estimated using Kaplan Meier method and described using median or rate at specific time points with 95 % confidence intervals (CI). Follow-up was calculated using reverse Kaplan-Meier estimation.
The association of parameters at enrolment with OS was assessed using univariate Cox to produce the hazard ratio (HR) and 95 % confidence intervals firstly for preoperative parameters, tumoral, therapeutic and lymphocyte count factors. Separate multivariate cox-analysis in conventional factors block (pre-operative, tumoral and therapeutic) and lymphocyte count block were performed with stepwise backwardelimination with the inclusion of variable with p < 0.05 in univariate analysis. The factors identified in these analyses were thereafter included in a final multivariate model with stepwise backwardelimination. When used in continuous in the Cox modelisation, lymphocyte count variable had to be normalized by logarithmic transformation, considering its skewed distribution. The lymph node ratio was defined by the ratio of the number of involved lymph nodes reported to the total number of lymph nodes removed in the lymph node dissection. The threshold 0.2 was kept as proposed in the study of Yamamoto et al. [25]. Hazard proportionality was checked by plotting log-minus-log survival curves.
Accuracy of the final model was checked regarding two parameters: discrimination and calibration. The predictive value and the discrimination ability of the model were evaluated with Harrell's Concordance (C)-index. One thousand random samples of the population were used to derive 95 % CI for the Harrell's Concordance statistic. Calibration and goodness of fit of the model were assessed by using the extension of Hosmer-Lemeshow test and for survival analysis and p-value greater than 0.1 was considered as an indicator for acceptable agreement. Calibration was also assessed by visual examination of calibration plot. Internal validity of the model was assessed by a bootstrap sample procedure. Several approaches have been proposed to assess the performance in samples of the same population (internal validation). Sensitivity analyses were performed for univariate and multivariate Cox models with a stratified approach on the cohort set parameter that allowed considering the two cohort heterogeneities in the Cox modelisation.
The predictive value that lymphocytes counts variables added to a reference risk model (all parameters identified in the multivariate final model except lymphocytes counts variables) was evaluated with the use of C-statistic. This analysis was repeated 1000 times using bootstrap samples to derive 95 % CI for the difference in the C-statistics between the two models in order to finally, assess the improvement in discrimination of lymphocytes counts parameters among the other conventional parameters.
We also used net continuous reclassification improvement (NRI) and integrated discrimination improvement (IDI) to quantify the performance and the net benefit of the addition of lymphocyte count to the reference model for the prediction of 48 months death probability. Continuous NRI has several limitations but would give a consistent message and is therefore a descriptive marker. One should note, cNRI does not consider the magnitude of the change, but only the direction. This is done by the IDI. When significantly greater than 0, IDI and cNRI are in favour of a net benefit of the addition of the marker of interest to the reference model considered.
Then, we investigated the possibility to provide a simple implementation of lymphocyte count parameter in clinical practice, guided by the determination for a relevant cutoff in order to categorize patients at baseline and postoperative time. According to the study of Ray-Coquard and al the threshold chosen was 1000 cells/mm 3 [26]. In addition, in clinical current practice considerations, patients with a lymphocyte count lower than 1000 are commonly considered in a lymphopenic status.
Finally, considering the added value of lymphocyte count measurements for OS risk stratification among conventional factors and their independent association with OS, we investigated the interest for a combination of pre and post-operative lymphocyte count in clinical practice.
The analyses were conducted using SAS 9.2 (Statistical Analysis System, Cary, NC, USA) and R 3.0.2 (R foundation for Statistical Computing). All statistical tests were 2-sided, and probability values <0.05 were regarded as significant.

Population
Characteristics of the overall population and according to the two cohort sets are given in Table 1. A total of 390 patients with resectable pancreatic cancers were enrolled. There were 192 (49,2 %) patients in cohort 1 and 198 (50,8 %) in cohort 2. Pathological findings differed significantly between the two cohorts but the most frequent T stage was T3-T4 in both cohorts (85 and 94 % respectively in cohort 1 and 2; p < 0.01). However OS of the patients in the two cohorts were not significantly different (Additional file 1: Figure S1, HR = 0.870 95 % CI: 0.696-1.088 p = 0.2235).

Independent prognostic factors of OS
The association of pre-operative, tumoral and therapeutic post-operative factors, as well as lymphocyte counts with OS in univariate analysis is shown in Table 2. We identified 9 variables as prognostic factors for OS in the univariate analysis: age at surgery (p = 0.003), serum albumin (p = 0.009), lymph nodes ratio ≥ 0.  Table S1A). Similarly, a separate multivariate cox-analysis in lymphocyte count parameters identified pre and postoperative lymphocyte counts as factors independently associated with OS (Additional file 2: Table S1B; p = 0.02 and p = 0.0467 respectively). Factors identified in these two previous multivariate analyses were thereafter included in a final multivariate model presented in Table 3. The final multivariate model exhibited four parameters significantly independently associated with OS with a p value <0.05 and one parameter borderline probably due to a lack of power: lymph nodes ratio (p = 0.0001), venous emboli (p = 0.0114), adjuvant chemotherapy (p = 0.0014), pre and post-operative lymphocyte counts (p = 0.0128 and p = 0.0764 respectively). When considered as continuous pre and post-operative lymphocyte count variables, their non-parametric Spearman correlation coefficient is equal to 0.36446 (p-value < 0.0001). Then, there is a moderate correlation between the two parameters allowing their consideration in the final multivariate model development.

Final multivariate model performance assessment
Accuracy of the model was checked regarding two parameters: discrimination and calibration, which measure the ability to separate patients with different prognosis and to provide unbiased survival predictions in groups of similar patients, respectively. Our final multivariable Cox model exhibited good calibration as shown in the calibration plot (Additional file 3: Figure S2) and acceptable discrimination (C-statistic 0.64; 95 % CI: 0.60-0.69).
With the replicated datasets (n = 1000) derived from the bootstrap sample procedure, uncertainties around hazard ratio estimates can be measured. Bootstrapping results for the internal validation reflect the robustness of the final model as presented in Table 3. A sensitivity analysis was performed to validate the robustness of our final model with a stratified approach. By forcing prognostic factors not involved in the multivariate analysis (T, N and age) results remained similar reflecting the robustness of our final model (Additional file 4: Table S2).
Additional value of pre and post-operative lymphocyte count parameters for OS prediction The inclusion of pre and post-operative lymphocyte count parameters in the reference model (including only conventional parameters) was found to significantly improve the discriminative ability of the model, because the C statistic increased significantly from 0.60-0.64 (bootstrap mean difference = 0.04; 95 % CI, 0.01-0.06). These results show that the addition of lymphocyte count parameter to clinical conventional parameters improved the stratification of patients at risk for death and then the model discrimination capacity. Similarly, the addition of lymphocyte count block to the conventional parameter block adequately reclassified at 48 months patients at lower risk for death and those at higher risk, as shown by a continuous net reclassification improvement of 0.3355 (95 % CI, 0.0719-0.5991; p = 0.01261; Fig. 1) and the integrated discrimination improvement was 0.03 (95 % CI, 0.01-0.06; p = 0.00339).
Medians (minimal-maximal) of pre and post-operative lymphocyte counts in our population were 1320 (150-8350) and 1410 (200-32000), respectively. Thus, the total lymphocyte count was categorized using a threshold of 1000 cells/mm 3 . Among the 390 patients involved in the final analysis, 110 (28 %) had lymphocyte count below 1000/mm3 at baseline and exhibited different prognostic profiles for OS (p = 0.0028). Post-operative lymphocyte count parameter was available for 301 patients. 65 (22 %) of them had post-operative lymphopenia. A pejorative correlation with OS was also evidenced (Table 2, Fig. 2a and b, p < 0.0001).
By combining the categorical (>1000; ≤1000) information in pre and post-operative lymphocyte counts we identified 4 groups of patients with different prognostic profiles. Within patient's pre and post-operative data available, 219 had a baseline lymphocyte count above 1000/mm 3 . 37 (17 %) of them were classified as lymphopenic 1 month after surgery and had poor prognosis (HR = 2.2; 95 % CI: 1.476-3.317 p < 0.0001). In addition, among the 83 patients of this cohort who displayed pre-operative lymphopenia, 55 (66 %) normalized their lymphocyte count after surgery and had a better prognosis (HR = 1.5; 95 % CI: 1.058; 2.088 p < 0.0001). Patients with the worst outcomes following surgery were those with pre and post-operative lymphopenia (HR = 2.340; 95 % CI: 1.524-3.593; p < 0.0001). Finally, we observed that the existence of post-operative lymphopenia categorizes patients into one group with poor prognosis (HR more than 2) whatever the pre-operative lymphocyte count. (Table 2, Fig. 2c and d).

Long-term survivor patient's description
Among the 390 patients enrolled, 67 (17 %) and 28 (7 %) were alive at 48 and 72 months respectively (Table 4). Lymph node ratio, venous emboli and adjuvant chemotherapy were identified as predictive factors for longterm survival. Among the 67 patients alive at 48 months, pre and post-operative lymphocyte counts were available for 67 and 57 patients respectively. 56 (84 %) and 54 (95 %) patients alive at 48 months had pre and post-operative lymphocyte counts above 1000/mm3 respectively. 82 (27,2 %) patients of the overall population displayed baseline lymphopenia. Interestingly, patients who recovered an absolute lymphocyte count above 1000/mm 3 1 month after surgery (n = 54) had a similar probability to be alive at 6 years compared to the 182 patients who were never classified as lymphopenic (6 year survival rate of 10 % and 9.75 % respectively). Conversely, the 6year survival rates of patients with lymphopenia before/ after surgery or who became lymphopenic post-surgery were respectively 0 and 2.7 %. Altogether, the analysis of patients in long-term remission showed that 54 (22.8 %) and 23-236 (9.75 %) with a 1-month lymphocyte count above 1000/mm 3 were alive at 4 and 6 years after surgery, respectively. By contrast, the probability of survival at 4 and 6 year was only 4.6 % (n = 3) and 1.5 % (n = 1) among patients with post-operative lymphopenia compared to 11 (n = 9) and 6 % (n = 5) for patients with pre-operative lymphopenia (Table 4).

Discussion
The prognostic value of lymphopenia has been previously identified in different malignancies including pancreatic cancers [15][16][17][18][19]27]. However, these results did not sustain the use of lymphopenia as a stratification criterion for clinical trials or to predict overall survival. Our study contributes to better determine the added value of lymphopenia in localized pancreatic cancer.  Neutrophil-Lymphocyte ratio (NLR) was previously proposed as an independent prognostic factor for overall survival both in localized [18][19][20] and in metastatic ductal pancreatic adenocarcinoma [28]. Nevertheless, NLR has some limitations. It includes two potentially independent biological factors. While neutrophils are related to inflammation, lymphocytes are directly involved in immune regulation. Moreover, NLR cut off differs from one study to another. In addition, other biological parameters related to inflammation, such as C-reactive protein, were shown to be significantly correlated to a poor clinical outcome [29].
We recorded 9 studies addressing the potential prognostic impact of pre-operative lymphopenia in localized pancreatic cancers (Additional file 5: Table S3). Only three of them identified NLR as an independent prognostic factor in multivariate analysis [18,19,30]. Our study confirms with statistical robustness that pre-operative lymphocyte count is an independent prognostic factor in pancreatic cancer on a larger scale (Table 3 and Fig. 2a; HR of 0,64; 95 % CI 0.450-0.909; p = 0.0013). The median of lymphocyte count in our cohort is in line with those observed previously in the studies of Garcea et al., and Stotz et al. [18,31].
Having observed that 66 % of the patients in the pre-operative lymphopenia group had a total lymphocyte count above 1000/mm 3 1 month after surgery, the prognostic value of post-operative lymphopenia was also investigated and demonstrated. (Table 3 and  Romano et al., have shown that post-operative immunodepression was significantly higher in pancreatic cancers than in colorectal and gastric cancers. Interestingly, recovery of normal post-operative lymphocyte count was longer in pancreatic cancers [32]. Only one small-scale study reported a negative relation between post-operative day 1 lymphopenia and overall survival for 111 patients with pancreatic cancer (p = 0.0029) [33]. However, most of the patients recovered from their decreased lymphocyte count several days following surgery and we postulated that lymphopenia monitored 1-month after the surgery might be more relevant to explore its influence on long term survival.
Of note, there are some limitations in our study. First, there are some differences between the two cohorts. Despite these differences the survival prognosis of the patients in the two cohorts was not significantly different (Additional file 1: Figure S1, Log-rank p value = 0.2236). From a statistical point of view, the assessment of model performance measures such as discrimination, calibration and internal validation strengthen the present investigation.
Addition of lymphocyte count variable to the conventional parameter block, in multivariate analysis, significantly improved the model discrimination capacity because the C statistic increased significantly from 0.60 to 0.64 (bootstrap mean difference = 0.04; 95 % CI, 0.01-0.06) demonstrating the additive value of lymphopenia to other conventional parameters. The use of a threshold offered better discrimination than the use of lymphocyte count because it allows death-risk stratification. In addition, combining the categorical (≤1000; >1000) information in pre and post-operative lymphocyte counts permitted the identification of several subgroups of patients with different prognoses. Patients with the worse prognosis were those with pre and post-operative lymphopenia (HR = 2.340; 95 % CI: 1.524-3.593; p < 0.0001). Patients who displayed lymphopenia prior to surgery and recovered absolute lymphocyte count above 1000/mm 3 1 month following pancreatic cancer removal had a better prognostic than patients without correction of lymphopenia (p < 0.0001). Consequently, the measure of lymphopenia seems more discriminant in the post-operative setting and has improved additive value for death risk stratification. Finally, the analysis of long-term survival patients showed that 23 of the 24 patients alive 6 years after the surgery had 1-month lymphocyte count above 1000/mm 3 .
These results suggest that lymphopenia is one of the most important prognostic factors to predict long term overall survival. The impact of lymphopenia on longterm survival was also reported in metastatic patients. Indeed, a recent study including patients treated with nabpaclitaxel and gemcitabine or with gemcitabine alone reported that the number of patients alive at least 24 months after treatment initiation was increased if NLR was below 5 [28]. Such results support the inclusion of lymphopenia as a risk stratification criterion in clinical trials and in models to determine the probability of overall survival.
Prospective immunological monitoring of those patients is needed to better explain the precise mechanisms involved in lymphocyte homeostasis in pancreatic cancer patients. The role of the immune system was pointed out by studies investigating the influence of TIL on pancreatic cancer prognosis. The frequency of CD8 + T lymphocytes was correlated to favourable clinical outcomes and prolonged survival [34][35][36]. The polarization of CD4+ T lymphocytes is another possible relevant immunological parameter correlated to patients' survival in several cancers [37]. GATA3/Tbet ratio and TH2 infiltrates were proposed as an independent prognostic factor for pancreatic cancer patients treated by surgery [38]. In this study, a predominant TH2 infiltrate was observed in peritumoral stroma suggesting a skewing of local immune responses toward TH2 polarization [38,39]. Moreover, regulatory T cell infiltration in pancreatic cancer tissue increases during disease progression and was evidenced as a prognostic factor in resected pancreatic cancers [40]. On the other hand, lymphopenia might reflect a global metabolism alteration such as malnutrition [41]. Albuminemia was monitored in 191 out of the 390 (49 %) patients included in our cohort. We observed an influence of hypoalbuminemia on OS in univariate analysis (HR = 0.97, p = 0.009). Albuminemia was initially droped out of multivariate analysis due to it's high rate of missing data. When performing the multivariate analysis with albuminemia (n = 152) similar results were obtained for veinous emboli and chemotherapy parameter (p = 0.0035 and p = 0.0010, respectively), even if the low number of data prompted further investigations to conclude. However, the direct role of albuminemia as an independent prognostic factor remained unclarified, as this parameter is mostly not significant in multivariate analysis [6]. The hypothesis of a direct influence of the tumour on lymphocyte homeostasis is plausible and supported by the capacity of the tumour to secret immunosuppressive cytokines like IL-10 and TGFβ [37]. Another possible hypothesis might be the production of lymphocyteapoptosis inducers such as Fas-Ligand, by pancreatic ductal adenocarcinoma [42].
One-month post-operative lymphopenia has independent and additive values for death risk stratification in localized pancreatic adenocarcinoma. The clinical significance of lymphopenia after surgery is highlighted by its negative correlation with the probability of longterm survival. The development of strategies based on immunonutrition [43], recombinant IL-7 [44], to expand CD4 T cells and the preliminary results of novel immunotherapies [45], offer new therapeutic endpoints to be assessed in pancreatic cancer patients.