Skip to main content

Decreased CX3CR1 messenger RNA expression is an independent molecular biomarker of early and late mortality in critically ill patients

Abstract

Background

Chemokine (C-X3-C motif) receptor 1 (CX3CR1) was identified as the most differentially expressed gene between survivors and non-survivors in two independent cohorts of septic shock patients and was proposed as a marker of sepsis-induced immunosuppression. Whether such a biomarker is associated with mortality in the heterogeneous group of critically ill patients is unknown. The primary objective of this study was to evaluate the association between CX3CR1 messenger RNA (mRNA) expression and mortality in intensive care unit (ICU) patients. The secondary objective was to evaluate similar endpoints in the subgroup of septic shock patients.

Methods

We performed a prospective, multicentre, non-interventional study in six ICUs of university hospitals in Lyon, France. Every consecutive adult patient with systemic inflammatory response syndrome and an expected length of stay in the ICU over 2 days was included. Whole-blood CX3CR1 mRNA expression was measured by quantitative real-time polymerase chain reaction at day 1 (D1) and D3 after inclusion.

Results

In ICU patients (n = 725), decreased CX3CR1 mRNA expression at D1 was associated with high D7 mortality (AUC 0.70, adjusted OR [aOR] 2.03, 95 % CI 1.19–3.46), while decreased expression at D3 was associated with increased D28 mortality (AUC 0.64, aOR 2.34, 95 % CI 1.45–3.77). In septic shock patients (n = 279), similar associations were observed between decreased D1 CX3CR1 mRNA expression and D7 mortality (AUC 0.69, aOR 2.76, 95 % CI 1.32–5.75) as well as decreased D3 expression and D28 mortality (AUC 0.72, aOR 3.98, 95 % CI 1.72–9.23). These associations were independent of lactacidaemia, Simplified Acute Physiology Score II, Sepsis-related Organ Failure Assessment score and Charlson comorbidity index.

Conclusions

This study represents the largest evaluation of such an mRNA marker in a heterogeneous cohort of severely injured patients. Our results show that decreased CX3CR1 mRNA expression is associated with increased mortality in ICU patients. This suggests a link between injury-induced immunosuppression and mortality in critically ill patients. In this context, the monitoring of such a host response molecular biomarker could prove very helpful for the identification of patients at high risk of death in the ICU.

Background

Intensive care units (ICUs) treat a heterogeneous group of patients with life-threatening injuries requiring close monitoring and support with high-technology medical equipment and multiple medications. In this context, the search for prognostic biomarkers that would help clinicians predict patients’ outcomes has been ongoing for years [1].

Several aspects of risk stratification biomarkers could be of interest in the ICU, especially in the context of financial constraints faced by the health care system and the strong need for optimization of limited medical resources. Indeed, besides use as benchmarking tools and in inter-ICU comparative studies, accurate outcome prediction based on biomarkers could be extremely useful for improving clinical decision-making. For example, biomarkers could help clinicians to dispense relevant treatments to patients who might benefit from them [2].

Transcriptomics has recently emerged as a novel tool to identify such biomarkers [38]. Most important, the availability of fully automated and user-friendly molecular platforms offers concrete solutions for routine clinical applications, even in ICU and emergency settings.

We previously identified chemokine (C-X3-C motif) receptor 1 (CX3CR1) as the most differentially expressed gene between survivors and non-survivors after septic shock [9, 10]. CX3CR1 is the sole receptor for fractalkine (CX3CL1) expressed on leucocyte subpopulations. Interestingly, we observed an association between decreased CX3CR1 expression on monocytes after septic shock and decreased functional responses of monocytes ex vivo [10]. This provides a link between decreased CX3CR1 and immune dysfunction after sepsis and suggests that the measure of such a host response messenger RNA (mRNA) marker may be useful in evaluation of injury-induced immune alterations. However, whether such a biomarker could be useful in the heterogeneous group of ICU patients is unknown. Therefore, we conducted a prospective, multicentre, non-interventional study with the primary objective of studying the association between CX3CR1 mRNA expression and risk of death in a large cohort of critically ill patients. The secondary objective was to evaluate the association between CX3CR1 mRNA expression and risk of death in the subgroup of septic shock patients.

Methods

Experimental design

This prospective, multicentre, non-interventional study was conducted in six ICUs in Lyon, France. It was approved by our institutional ethical review board (Comité d’Ethique des Centres d’investigation Clinique de l’Inter-Région Rhône-Alpes Auvergne – IRB 5044).

Follow-up and outcomes

The primary outcome was the association of day 1 (D1) and D3 CX3CR1 mRNA expression and risk of death during the hospital stay of ICU patients with systemic inflammatory response syndrome (SIRS). The secondary outcome was the association of D1 and D3 CX3CR1 mRNA expression and risk of death during the hospital stay in the subgroup of septic shock patients. In both cases, mortality was evaluated at 7 days (D7) and at 28 days (D28).

Since half of the deaths in our cohort occurred before D7 and the median ICU stay was 7 days, this time point (i.e., D7) was used to define early mortality among the 725 patients eligible at D1. Conversely, D28 was used to define late mortality after ICU admission among the 515 patients eligible at D3.

Patients

From December 2009 to June 2011, every consecutive patient aged ≥18 years with an expected length of stay in the ICU of more than 2 days was prospectively enrolled if the patient met the criteria for SIRS as described in the American College of Chest Physicians/Society of Critical Care Medicine (ACCP/SCCM) 1992 consensus statement [11] and if non-opposition to inclusion in the protocol was obtained.

Exclusion criteria comprised immunosuppression as defined by immunosuppressive treatment or corticoid treatment with dosage >10 mg/day or cumulative dose >700 mg equivalent prednisolone, aplasia as defined by number of circulating neutrophils <500 cells/mm3, primary innate immune deficiency, and extracorporeal circulation during the month before ICU admission.

No informed consent was needed, as this study was non-interventional and complementary blood samples in PAXgene® tubes (PreAnalytiX, Hombrechtikon, Switzerland) were obtained during patients’ routine blood sampling after completion of routine follow-up tests. Nevertheless, non-opposition to inclusion in the protocol was recorded from every patient or the patient’s next of kin.

Definitions

SIRS was defined as the presence of at least two of the following clinical criteria: temperature >38 °C or <36 °C, heart rate >90 beats/minute, respiratory rate >20 breaths per minute or PaCO2 < 32 mmHg (4.3 kPa), and leucocyte count >12,000/mm3 or <4000/mm3 [11].

Sepsis was defined as the presence of a proven (visible either clinically/surgically, radiologically or microbiologically) infection or a highly suspected infection at inclusion. Following the definitions proposed by Vincent et al., which partly rule out the former definitions of sepsis and severe sepsis described in the ACCP/SCCM 1992 consensus conference statement, sepsis was simply defined as an infection requiring ICU admission [12].

Shock was defined as persistent hypotension despite adequate fluid resuscitation requiring the use of epinephrine or norepinephrine at a dose >0.25 μg/kg/minute [13].

Antibiotic treatment was considered to be appropriate if all the pathogenic microorganisms were sensitive to at least one of the administered antimicrobial drugs as determined by an in vitro sensitivity pattern. In cases of culture-negative infection, treatment was considered appropriate if consistent with the local anti-infectious protocols.

Data collection at inclusion

The following demographic, clinical and biochemical data were collected at patient admission and/or inclusion in the protocol:

  1. 1.

    Demographic characteristics: age and sex

  2. 2.

    Clinical scores: Simplified Acute Physiology Score II (SAPS II) (range 0–163 [14], calculated at inclusion), Sepsis-related Organ Failure Assessment (SOFA) score (range 0–24 [15], calculated at inclusion and at D3) and Charlson comorbidity index (range 0–32 [16] for which age was not taken into account).

  3. 3.

    Lactic acid concentration (at D1 and at D3)

  4. 4.

    Need for respiratory (mechanical ventilation), haemodynamic (vasopressor treatment) or renal (renal replacement therapy) support

  5. 5.

    Presence of shock (defined as a cardiovascular SOFA score of 4)

  6. 6.

    Diagnostic category (planned or urgent, surgery or medical)

  7. 7.

    Presence of trauma

  8. 8.

    Presence of sepsis (either community- or hospital-acquired)

  9. 9.

    Characteristics of infection if relevant (site, type of detection, identified germs)

  10. 10.

    Presence of an anti-infectious treatment if relevant and its adequacy either to the antibiogram or to the local anti-infectious protocol

RNA extraction, reverse transcription and quantitative polymerase chain reaction

PAXgene® tubes were collected from patients within the first 12 h after inclusion in the protocol (D1) and at D3. Total RNA was extracted using the PAXgene® Blood RNA Kit (PreAnalytiX). Before RNA elution, the residual genomic DNA was digested using the RNase-Free DNase Set (QIAGEN, Hilden, Germany). Total RNA was reverse-transcribed into complementary DNA (cDNA) using the SuperScript VILO cDNA Synthesis Kit (Life Technologies, Carlsbad, CA, USA). The candidate gene was quantified using quantitative real-time polymerase chain reactions. Polymerase chain reactions (PCRs) were performed in a LightCycler instrument (Roche Diagnostics, Risch-Rotkreuz, Switzerland) using the standard TaqMan Fast Advanced Master Mix PCR kit according to the manufacturer’s instructions (Applied Biosystems, Foster City, CA, USA). Thermocycling was performed in a final volume of 20 μl containing 5 μM of required primers and 1 μM of required probe. PCR was performed with an initial denaturation step of 10 minutes at 95 °C, followed by 45 cycles of a touch-down PCR protocol (10 seconds at 95 °C, 29 seconds of annealing at 68 °C and a 1-second extension at 72 °C). The cDNA standards were prepared from purified PCR amplicons obtained with the corresponding primers (Additional file 1: Table S1). The second derivative maximum method was used with the LightCycler software to automatically determine the crossing point for individual samples, as previously described. Standard curves were generated by using quadruplicate cDNA standard. Relative standard curves describing the PCR efficiency of selected genes were created and used to perform efficiency-corrected quantification with the LightCycler Relative Quantification software (Roche Diagnostics). Gene expression normalization was performed using a selected housekeeping gene (hypoxanthine phosphoribosyltransferase 1 [HPRT1]), and results were expressed as calibrated normalized relative quantity (CNRQ) [17]. CX3CR1 expression and CX3CR1 CNRQ are equivalent terms.

Statistical analysis

The descriptive analyses comparing surviving and non-surviving patients at D7 and D28 (all-cause mortality) after ICU admission were completed with the usual appropriate tests, such as the t test or Mann-Whitney U test for quantitative variables and Pearson’s χ2 test for qualitative variables. The Monte Carlo resampling method was used to estimate the p value of the χ2 test of independence in cases of an expected frequency <5. ICU survival curves were created using the Kaplan-Meier method. Every patient was followed for 28 days in the ICU and after ICU discharge. The results of the log-rank tests were associated with the survival representations. ROC curves were built for CX3CR1 expression levels at D1 and D3 according to D7 and D28 mortality, respectively. Best cut-off values (i.e., maximized sensitivity and specificity) for CX3CR1 mRNA expression prediction of mortality were identified according to the calculation of Youden indexes derived from ROC curve analyses. Univariate and multivariate analyses of mortality outcomes at D7 and D28 were studied with logistic regression. Parameters with a p value <0.1 in univariate analyses were kept in multivariate analyses. Parameters with >20 % missing values were not included in multivariate analyses. The 0.632+ bootstrap method was applied to adjust for potential overestimation of diagnostic performance and to estimate internal validity. For each model, 1000 resamplings with replacement were performed. The 95 % CI represents the 2.5th and 97.5th percentiles of estimates obtained in bootstrap resamples. All bootstrap estimations were performed using the ModelGood package in R. The level of significance was set at 5 %, and the results were described with a 95 % CI. The analyses were performed with IBM SPSS Statistics version 20 software (IBM, Armonk, NY, USA).

Results

Cohort description

Between December 2009 and June 2011, 852 ICU patients were screened for eligibility. After exclusion due to ICU discharge within the first 48 h, 749 patients were included. Because of technical issues, samples from 24 of these patients could not be processed at D1. Therefore, 725 patients were ultimately analysed (Fig. 1). Their clinical characteristics are described in Table 1. At D3, 515 patients were still alive and had samples available for analysis. The clinical characteristics of this subcohort of patients at admission were not different from those of the overall cohort.

Fig. 1
figure 1

Study flowchart. Between December 2009 and June 2011, a total of 852 critically ill patients were screened for eligibility. After exclusion for technical reasons or intensive care unit discharge within the first 48 h, 725 patients were ultimately included in the study. D day

Table 1 Clinical and demographic characteristics of critically ill patients in the study

In the total cohort (n = 749), 220 patients died within 28 days (D28 mortality 29 %) (Additional file 2: Figure S1). Half of these deaths occurred before D7 (D7 mortality 15 %) and one-third occurred before D3 (D3 mortality 10 %). The median ICU stay was 7 days.

In the D1–D7 period (n = 725), non-survivors were significantly older (p < 0.001) and presented with significantly higher SAPS II (p < 0.001) and SOFA (p < 0.001) scores and increased lactacidaemia at D1 (p < 0.001). Non-survivors were more frequently in shock at ICU admission (p < 0.001) but had infections less frequently than survivors (p = 0.023).

In the D3–D28 period in the subgroup of patients alive at D3 (n = 515), D28 mortality was 23 % (n = 117 non-survivors). Non-survivors were significantly older (p < 0.001) and had higher SAPS II (p < 0.001) and SOFA (p < 0.001) scores and increased lactacidaemia at D3 (p < 0.001).

Primary endpoint

At D1, CX3CR1 mRNA expression was significantly lower in non-survivors than in D7 survivors (p < 0.0001) (Fig. 2a).

Fig. 2
figure 2

Association between chemokine (C-X3-C motif) receptor 1 (CX3CR1) messenger RNA (mRNA) levels at day 1 (D1) or day 3 and early (day 7) or late (day 28) mortality in critically ill patients. CX3CR1 mRNA levels were measured by quantitative real-time polymerase chain reaction in whole-blood samples of a cohort of 725 adult intensive care unit (ICU) patients. The Mann-Whitney U test was used to compare CX3CR1 mRNA levels in survivors versus non-survivors. Kaplan-Meier survival curves were established after stratification based on thresholds defined in ROC curve analyses. A p value <0.05 was considered statistically significant. a CX3CR1 mRNA levels at D1 in ICU patients were significantly different between D7 survivors and non-survivors. b Log-rank test showed that patients with D1 CX3CR1 mRNA levels above the threshold had a significantly better survival than patients with CX3CR1 mRNA levels below the threshold. c CX3CR1 mRNA levels at D3 in ICU patients were significantly different between D28 survivors and non-survivors. d Log-rank test showed that patients with D3 CX3CR1 mRNA levels above the threshold had a significantly better survival than patients with CX3CR1 mRNA levels below the threshold

On the basis of ROC curve analysis, a threshold of 0.085 was defined as the best cut-off value for the prediction of D7 mortality (AUC 0.70, 95 % CI 0.65–0.76, sensitivity [Se] 55 %, specificity [Sp] 77 %, positive predictive value [PPV] 30 %, negative predictive value [NPV] 91 %, positive likelihood ratio [LR+] 2.42, negative likelihood ratio [LR−] 0.58). The internal validity of these results was confirmed using a bootstrapping technique (AUC 0.70, 95 % CI 0.63–0.77). Using this threshold, Kaplan-Meier survival curves were established (Fig. 2b). Log-rank test analysis showed that patients with D1 CX3CR1 mRNA expression over the threshold had better survival than patients with CX3CR1 expression below the cut-off (91 % and 70 %, respectively; p < 0.001).

In univariate and multivariate logistic regression analyses, low CX3CR1 expression (below threshold of 0.085) remained significantly associated with D7 mortality independently of confounding parameters (adjusted OR [aOR] 2.03, 95 % CI 1.19–3.46, p = 0.009) (Table 2). In these models, SAPS II and SOFA scores, presence of sepsis and lactacidaemia at D1 were also associated with D7 mortality. Of note, because the SOFA score was based on norepinephrine or epinephrine administration, a similar multivariate analysis including the presence of shock instead of the SOFA score was performed. In this analysis, neither shock nor Charlson comorbidity index was associated with D7 mortality (Additional file 3: Table S2A).

Table 2 Univariate and multivariate analyses of mortality according to CX3CR1 messenger RNA expression in critically ill patients

At D3, CX3CR1 mRNA expression was lower in patients who died within 28 days after admission (n = 117) versus survivors (n = 398) (p < 0.0001) (Fig. 2c). Based on ROC curve analysis, the best cut-off value for D3 CX3CR1 mRNA expression association with D28 mortality was 0.246 (AUC 0.64, 95 % CI 0.58–0.71, Se 73 %, Sp 56 %, PPV 85 %, NPV 38 %, LR+ 1.65, LR− 0.48). Again, the internal validity of these results was confirmed using a bootstrapping technique (AUC 0.64, 95 % CI 0.56–0.72). Kaplan-Meier survival curves were generated with this threshold. Log-rank test analysis showed that D28 survival was significantly higher in the group of patients with D3 CX3CR1 expression over threshold than in patients with CX3CR1 expression below the cut-off value (85 % versus 62 %, respectively; p < 0.001) (Fig. 2d). Finally, in multivariate analysis, low CX3CR1 expression at D3 was associated with increased risk of death at D28, with an aOR of 2.34 (95 % CI 1.45–3.77; p < 0.001) (Table 2B). This association was independent of SAPS II and SOFA scores as well as the Charlson comorbidity index. Of note, due to the high number of missing values for lactacidaemia at D3, this parameter was not included in the multivariate analysis. In multivariate analysis, replacement of SOFA score at D3 by presence of shock at D3 did not reveal any significant association between this last parameter and D28 mortality (Additional file 3: Table S2B).

Secondary endpoint: specificity of septic shock patients

The clinical characteristics of septic shock patients are presented in Table 3.

Table 3 Clinical and demographic characteristics of septic shock patients

At D1, CX3CR1 mRNA expression was lower in septic shock patients who died within the first week after ICU admission than in D7 survivors (p < 0.0001) (Fig. 3a).

Fig. 3
figure 3

Association between chemokine (C-X3-C motif) receptor 1 (CX3CR1) messenger RNA (mRNA) levels at day 1 (D1) or day 3 and early (day 7) or late (day 28) mortality in septic shock patients. CX3CR1 mRNA level was measured by quantitative real-time polymerase chain reaction in whole-blood samples in a cohort of 279 adult septic shock patients. The Mann-Whitney U test was used to compare CX3CR1 mRNA levels in survivors versus non-survivors. Kaplan-Meier survival curves were established after stratification on the basis of thresholds defined in ROC curve analyses. A p value <0.05 was considered statistically significant. a CX3CR1 mRNA levels at D1 in septic shock patients were significantly different between D7 survivors and non-survivors. b Log-rank test showed that patients with D1 CX3CR1 mRNA levels over the threshold had a significantly better survival than patients with CX3CR1 mRNA levels below threshold. c CX3CR1 mRNA levels at D3 in septic shock patients were significantly different between D28 survivors and non-survivors. d Log-rank test showed that patients with D3 CX3CR1 mRNA levels over the threshold had a significantly better survival than patients with CX3CR1 mRNA levels below the threshold. ICU intensive care unit

On the basis of ROC curve analysis, a threshold of 0.073 was defined as the best cut-off value for the prediction of D7 mortality (AUC 0.67, 95 % CI 0.59–0.76, Se 63 %, Sp 68 %, PPV 33 %, NPV 88 %, LR+ 1.94, LR− 0.55). Again, the internal validity of these results was confirmed using a bootstrapping technique (AUC 0.67, 95 % CI 0.59–0.76). Using this threshold, Kaplan-Meier survival curves were established (Fig. 3b). Log-rank test analysis showed that patients with CX3CR1 mRNA expression at D1 over the threshold had a better survival than patients with CX3CR1 expression below the cut-off value (88 % and 65 %, respectively; p < 0.001).

In univariate and multivariate logistic regression analyses, low CX3CR1 expression (below threshold of 0.073) remained associated with D7 mortality independently of confounding parameters (aOR 2.76, 95 % CI 1.32–5.75) (Table 4). In this model, SAPS II was also associated with D7 mortality, whereas the Charlson comorbidity index score and lactacidaemia did not show any significant association.

Table 4 Univariate and multivariate analyses of mortality according to CX3CR1 messenger RNA expression in septic shock patients

At D3, 192 septic shock patients were still alive. CX3CR1 mRNA expression was significantly lower in patients who died before D28 than in survivors (Fig. 3c). Based on ROC curve analysis, the best cut-off value for D3 CX3CR1 mRNA association with D28 mortality was 0.253 (AUC 0.72, 95 % CI 0.61–0.81, Se 59 %, Sp 83 %, PPV 92 %, NPV 38 %, LR+ 3.38, LR− 0.50). Again, the internal validity of these results was confirmed using a bootstrapping technique (AUC 0.71, 95 % CI 0.58–0.84). Kaplan-Meier survival curves were generated with this threshold. Log-rank test analysis showed that 28-day survival was higher in the group of septic shock patients with D3 CX3CR1 expression over the threshold than in patients with CX3CR1 expression below the cut-off value (89 % versus 60 %, respectively; p < 0.001) (Fig. 3d). In univariate and multivariate analyses, low CX3CR1 expression at D3 was associated with increased risk of death at D28, with an aOR of 3.98 (95 % CI 1.72–9.23; p = 0.001) (Table 4). This was independent of SAPS II and SOFA scores as well as the Charlson comorbidity index score.

Discussion

We demonstrate, for the first time to our knowledge, that decreased CX3CR1 mRNA expression is associated with early and late mortality in critically ill patients independently of initial severity, organ dysfunction, existing comorbidities and metabolic distress. This represents the largest study of such transcriptomic markers’ association with mortality in ICU patients. In addition, this study validates, in a large cohort of patients, our previous results showing that decreased CX3CR1 mRNA is an independent predictor of death after septic shock. This suggests that the monitoring of this host response biomarker in ICU patients may provide information complementary to the usual risk stratification parameters.

Researchers in numerous small observational studies have described the rapid modification of several immune parameters after injury and the association between the extent or duration of these alterations and increased risk of death in these patients [18]. Therefore, in the context of clinical trials already underway in which researchers are evaluating immunoadjuvant therapies in sepsis [19], host response biomarkers are being evaluated as risk stratification parameters in ICU patients. Considering the constantly evolving immune response after severe injury, such biomarkers could be used as a stratification tool or for treatment follow-up [20]. While some markers, such as decreased human leucocyte antigen (HLA)-DR expression on monocytes or lymphopenia, have been proposed as immunomonitoring tools [20], the monitoring of several different markers included in a panel is probably required to accurately assess patients’ immune status after injury [21]. In this context, considering the recent developments of standardized molecular biology techniques, the evaluation of transcriptomic markers presents the advantage of now being fully automated and therefore available with very low pre-analytical constraints [20]. This would most likely facilitate their use on a routine basis, as opposed to markers monitored with less standardized techniques. However, to date, the capacity of such mRNA marker to predict deleterious outcomes in ICU patients has never been evaluated in a large cohort of critically ill patients. In the present study, we evaluated the performance of CX3CR1 mRNA level in the prediction of mortality in a cohort of more than 720 ICU patients. This represents the largest clinical evaluation of a molecular host response marker in ICU patients.

CX3CR1 is a receptor expressed on leucocyte subpopulations. The CX3CR1 ligand (fractalkine) is the only member of the CX3C chemokine subfamily. CX3CR1 has a role in cell chemotaxis to sites of inflammation through interaction with fractalkine fixed on activated endothelial cells. In 2007, Auffray et al. showed that CX3CR1 is overexpressed on “patroller” monocytes and is required for their rapid tissue invasion at the infection site [22]. Interestingly, several studies investigated the expression of CX3CR1 after severe injuries. First, in a microarray study, we identified this gene as the most downregulated gene between survivors and non-survivors of septic shock [9]. This result was confirmed in a subsequent cohort of septic shock patients at both the mRNA and protein levels [10]. Interestingly, decreased CX3CR1 expression on monocytes from septic patients was associated with their decreased activation after ex vivo stimulation. This provides a link between decreased CX3CR1 and immune dysfunction after sepsis. In further studies in murine models of sepsis, investigators proposed a role of this molecule in bacterial killing and host defence against infection as well as in sepsis-induced organ dysfunction [2326]. This all suggests that measurement of such mRNA markers may help in evaluation of injury-induced immune alterations. Interestingly, besides sepsis, decreased CX3CR1 mRNA expression was also described in sterile inflammation, such as major vascular surgery [27]. Of note, in a recent study of more than 1100 patients, Hoogendijk et al. reported that plasma fractalkine level was positively correlated with severity in patients with sepsis admitted to the ICU [28]. Therefore, the tandem CX3CR1-fractalkine may represent a good prognostic marker in ICU patients. However, the performance of this marker as a risk stratification tool in total ICU patients has not been described yet. Therefore, in the present study, we evaluated the association between CX3CR1 mRNA expression and risk of death in a large cohort of critically ill patients. In addition, we confirmed the association between CX3CR1 mRNA expression and risk of death in the subgroup of septic shock patients.

We chose to evaluate both early (i.e., before D7) and late (i.e., before D28) mortality, hypothesizing that parameters associated with mortality after ICU admission may evolve over time. Indeed, we observed that the clinical parameters associated with mortality were different when we considered early versus late deaths. This was the case both in the total cohort of ICU patients and in the subgroup of septic shock patients. For example, in the multivariate analyses, Charlson comorbidity index, illustrating patients’ comorbidities, was not associated with early mortality but rather with late deaths. Similarly, SOFA score at D1 or the presence of shock at D1 was not associated with early death, while the SOFA score at D3 was significantly associated with D28 mortality. This suggests that the later course of ICU patients with SIRS may be impacted principally by the patient’s medical background (such as comorbidities) and the persistence of organ dysfunction. Conversely, early deaths may be related mostly to the severity of the initial insult. Interestingly, host immune response assessed by CX3CR1 mRNA expression constantly remained associated with both early and late outcomes. This may illustrate that alterations of the host immune response after severe injury may play a role in mortality in the ICU and therefore supports the use of immunoadjuvant therapies in sepsis.

Different CX3CR1 thresholds were selected for prediction of D7 versus D28 mortality. Indeed, cut-off values for the predictive capacity of CX3CR1 at D1 were very low compared with cut-offs at D3. This may illustrate the evolution of the host immune response over time after severe aggression. Severe injury strongly impacts the host immune response by inducing a “genomic storm” [4, 6]. However, in patients who survive this initial aggression, their immune status should return to normal. Failure to return to normal may lead to late mortality [18, 19]. Interestingly, these cut-off values were very similar between ICU patients and the subcohort of septic shock patients. This suggests that underlying pathophysiological mechanisms associated with clinical severity are shared by sterile inflammation and sepsis [4, 12, 29, 30]. In our results, the association between lower CX3CR1 expression and fatal outcome was independent of the presence of either an infection or a shock. Interestingly, non-infectious SIRS seems to be as deleterious as septic shock [31, 32]. In the subgroup of non-infected patients with SIRS, we also observed a significant association between decreased CX3CR1 mRNA level at D1 and D7 mortality (AUC 0.75, 95 % CI 0.66–0.84). In univariate and multivariate logistic regression analyses, low CX3CR1 expression (below threshold of 0.12) in non-infected patients remained significantly associated with D7 mortality independently of confounding parameters (aOR 3.63, 95 % CI 1.43–9.20, p = 0.007). However, in this group of patients, no association between D3 CX3CR1 mRNA level and D28 mortality was observed.

This study has some limitations. In particular, we did not include measurement of other markers of injury-induced immunosuppression, such as monocytic HLA-DR. Therefore, we could not compare the performance of our transcriptomic marker versus such a “gold standard” parameter in the present cohort. In addition, as no specific immune function tests were performed, we cannot definitively prove that decreased CX3CR1 mRNA is associated with immune dysfunction in ICU patients. However, considering the interesting results obtained, these two aspects need to be confirmed in a dedicated clinical study.

Conclusions

We observed that decreased CX3CR1 mRNA expression is a very early phenomenon in critically ill patients and is associated with mortality independently of the presence of shock or sepsis. The recent progress in molecular biology techniques and the development of standardized, fully automated molecular biology platforms should ensure the availability of measurement of this marker around the clock for monitoring of host response in ICU patients. We suggest that the measurement of this marker should be included in a panel of host response markers in forthcoming clinical trials in severely injured patients to identify the group of patients with a high risk of death after initial injury.

Abbreviations

ACCP/SCCM, American College of Chest Physicians/Society of Critical Care Medicine; aOR, adjusted odds ratio; cDNA, complementary DNA; CNRQ, calibrated normalized relative quantity; CX3CR1, chemokine (C-X3-C motif) receptor 1; D, day; HLA-DR, human leucocyte antigen; ICU, intensive care unit; LR, likelihood ratio; mRNA, messenger RNA; NPV, negative predictive value; PCR, polymerase chain reaction; PPV, positive predictive value; qRT-PCR, quantitative real-time polymerase chain reaction; SAPS II, Simplified Acute Physiology Score II; Se, sensitivity; SIRS, systemic inflammatory response syndrome; SOFA, Sepsis-related Organ Failure Assessment; Sp, specificity

References

  1. Power GS, Harrison DA. Why try to predict ICU outcomes? Curr Opin Crit Care. 2014;20(5):544–9.

    Article  PubMed  Google Scholar 

  2. Kahn JM. Predicting outcome in critical care: past, present and future. Curr Opin Crit Care. 2014;20(5):542–3.

    Article  PubMed  Google Scholar 

  3. Marini JJ, Vincent JL, Annane D. Critical care evidence—new directions. JAMA. 2015;313(9):893–4.

    Article  CAS  PubMed  Google Scholar 

  4. Xiao W, Mindrinos MN, Seok J, Cuschieri J, Cuenca AG, Gao H, et al. A genomic storm in critically injured humans. J Exp Med. 2011;208(13):2581–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Wong HR, Cvijanovich NZ, Allen GL, Thomas NJ, Freishtat RJ, Anas N, et al. Validation of a gene expression-based subclassification strategy for pediatric septic shock. Crit Care Med. 2011;39(11):2511–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Cazalis MA, Lepape A, Venet F, Frager F, Mougin B, Vallin H, et al. Early and dynamic changes in gene expression in septic shock patients: a genome-wide approach. Intensive Care Med Exp. 2014;2(1):20.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Cazalis MA, Friggeri A, Cavé L, Demaret J, Barbalat V, Cerrato E, et al. Decreased HLA-DR antigen-associated invariant chain (CD74) mRNA expression predicts mortality after septic shock. Crit Care. 2013;17(6):R287.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Maslove DM, Wong HR. Gene expression profiling in sepsis: timing, tissue, and translational considerations. Trends Mol Med. 2014;20(4):204–13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Pachot A, Lepape A, Vey S, Bienvenu J, Mougin B, Monneret G. Systemic transcriptional analysis in survivor and non-survivor septic shock patients: a preliminary study. Immunol Lett. 2006;106(1):63–71.

    Article  CAS  PubMed  Google Scholar 

  10. Pachot A, Cazalis MA, Venet F, Turrel F, Faudot C, Voirin N, et al. Decreased expression of the fractalkine receptor CX3CR1 on circulating monocytes as new feature of sepsis-induced immunosuppression. J Immunol. 2008;180(9):6421–9.

    Article  CAS  PubMed  Google Scholar 

  11. Bone RC, Balk RA, Cerra FB, Dellinger RP, Fein AM, Knaus WA, et al. Definitions for sepsis and organ failure and guidelines for the use of innovative therapies in sepsis. Chest. 1992;101(6):1644–55.

    Article  CAS  PubMed  Google Scholar 

  12. Vincent JL, Opal SM, Marshall JC, Tracey KJ. Sepsis definitions: time for change. Lancet. 2013;381(9868):774–5.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Annane D, Sebille V, Charpentier C, Bollaert PE, Francois B, Korach JM, et al. Effect of treatment with low doses of hydrocortisone and fludrocortisone on mortality in patients with septic shock. JAMA. 2002;288(7):862–71.

    Article  CAS  PubMed  Google Scholar 

  14. Le Gall JR, Lemeshow S, Saulnier F. A new Simplified Acute Physiology Score (SAPS II) based on a European/North American multicenter study. JAMA. 1993;270(24):2957–63.

    Article  PubMed  Google Scholar 

  15. Vincent JL, de Mendonca A, Cantraine F, Moreno R, Takala J, Suter PM, et al. Use of the SOFA score to assess the incidence of organ dysfunction/failure in intensive care units: results of a multicenter, prospective study. Crit Care Med. 1998;26(11):1793–800.

    Article  CAS  PubMed  Google Scholar 

  16. Charlson ME, Pompei P, Ales KL, MacKenzie CR. A new method of classifying prognostic comorbidity in longitudinal studies: development and validation. J Chronic Dis. 1987;40(5):373–83.

    Article  CAS  PubMed  Google Scholar 

  17. Hellemans J, Mortier G, De Paepe A, Speleman F, Vandesompele J. qBase relative quantification framework and software for management and automated analysis of real-time quantitative PCR data. Genome Biol. 2007;8(2):R19.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Hotchkiss RS, Monneret G, Payen D. Sepsis-induced immunosuppression: from cellular dysfunctions to immunotherapy. Nat Rev Immunol. 2013;13(12):862–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Hotchkiss RS, Monneret G, Payen D. Immunosuppression in sepsis: a novel understanding of the disorder and a new therapeutic approach. Lancet Infect Dis. 2013;13(3):260–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Venet F, Lukaszewicz AC, Payen D, Hotchkiss R, Monneret G. Monitoring the immune response in sepsis: a rational approach to administration of immunoadjuvant therapies. Curr Opin Immunol. 2013;25(4):477–83.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. van der Poll T. Future of sepsis therapies. Crit Care. 2016;20(1):106.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Auffray C, Fogg D, Garfa M, Elain G, Join-Lambert O, Kayal S, et al. Monitoring of blood vessels and tissues by a population of monocytes with patrolling behavior. Science. 2007;317(5838):666–70.

    Article  CAS  PubMed  Google Scholar 

  23. Ishida Y, Hayashi T, Goto T, Kimura A, Akimoto S, Mukaida N, et al. Essential involvement of CX3CR1-mediated signals in the bactericidal host defense during septic peritonitis. J Immunol. 2008;181(6):4208–18.

    Article  CAS  PubMed  Google Scholar 

  24. Raspe C, Hocherl K, Rath S, Sauvant C, Bucher M. NF-kB-mediated inverse regulation of fractalkine and CX3CR1 during CLP-induced sepsis. Cytokine. 2013;61(1):97–103.

    Article  CAS  PubMed  Google Scholar 

  25. Chousterman BG, Boissonnas A, Poupel L. Baudesson de Chanville C, Adam J, Tabibzadeh N, et al. Ly6Chigh monocytes protect against kidney damage during sepsis via a CX3CR1-dependent adhesion mechanism. J Am Soc Nephrol. 2016;27(3):792–803.

    Article  PubMed  Google Scholar 

  26. An JL, Ishida Y, Kimura A, Tsokos M, Kondo T. Immunohistochemical detection of CCR2 and CX3CR1 in sepsis-induced lung injury. Forensic Sci Int. 2009;192(1-3):e21–5.

    Article  CAS  PubMed  Google Scholar 

  27. Feezor RJ, Baker HV, Xiao W, Lee WA, Huber TS, Mindrinos M, et al. Genomic and proteomic determinants of outcome in patients undergoing thoracoabdominal aortic aneurysm repair. J Immunol. 2004;172(11):7103–9.

    Article  CAS  PubMed  Google Scholar 

  28. Hoogendijk AJ, Wiewel MA, van Vught LA, Scicluna BP, Belkasim-Bohoudi H, Horn J, et al. Plasma fractalkine is a sustained marker of disease severity and outcome in sepsis patients. Crit Care. 2015;19:412.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Johnson SB, Lissauer M, Bochicchio GV, Moore R, Cross AS, Scalea TM. Gene expression profiles differentiate between sterile SIRS and early sepsis. Ann Surg. 2007;245(4):611–21.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Matzinger P. The danger model: a renewed sense of self. Science. 2002;296(5566):301–5.

    Article  CAS  PubMed  Google Scholar 

  31. Sprung CL, Sakr Y, Vincent JL, Le Gall JR, Reinhart K, Ranieri VM, et al. An evaluation of systemic inflammatory response syndrome signs in the Sepsis Occurrence in Acutely ill Patients (SOAP) study. Intensive Care Med. 2006;32(3):421–7.

    Article  PubMed  Google Scholar 

  32. Dulhunty JM, Lipman J, Finfer S. Does severe non-infectious SIRS differ from severe sepsis? Results from a multi-centre Australian and New Zealand intensive care unit study. Intensive Care Med. 2008;34(9):1654–61.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

MIP-REA study group

Dr. Guillaume Marcotte and Dr. Christian Guillaume (Hospices Civils de Lyon, Department of Anaesthesiology and Critical Care Medicine, Groupement Hospitalier Edouard Herriot, University Claude Bernard Lyon 1, Lyon, France); Dr. Romain Hernu, Sylvie De La Salle, Dr. Marie Simon, and Dr. Thomas Baudry (Hospices Civils de Lyon, Medical Intensive Care Unit, Groupement Hospitalier Edouard Herriot, Lyon, France); Elisabeth Cerrato, Emmanuelle Gallet-Gorius, Audrey Larue, and Christine Alberti-Segui (Hospices Civils de Lyon - bioMérieux Joint Research Unit, Groupement Hospitalier Edouard Herriot, Lyon, France); Nathalie Panel and Marion Provent (Hospices Civils de Lyon, Intensive Care Unit, Centre Hospitalier Lyon Sud, Pierre Bénite, France); Prof. Jean-Paul Viale (Hospices Civils de Lyon, Intensive Care Unit, Hôpital de la Croix Rousse, Lyon, France); and Anne Portier (Hospices Civils de Lyon, Immunology Laboratory, Groupement Hospitalier Edouard Herriot, Lyon, France).

Funding

This work was supported by funds from bioMérieux and Hospices Civils de Lyon. The other authors declare that they have no competing interests. This project is part of Advanced Diagnostics for New Therapeutic Approaches (ADNA), a program dedicated to personalized medicine that is coordinated by Institut Mérieux and supported by the French public agency BPI France.

Authors’ contributions

All authors were involved in the analysis and interpretation of data as well as drafting the manuscript or revising it critically for important intellectual content. AF, MAC, AP, AL, GM, JT and FV made substantial contributions to study conception and design and acquisition of data. SM, SB, DMB, AF, MAC, AP, AL, GM, JT and FV participated actively in statistical analysis and interpretation of data. All authors read and approved the final manuscript.

Competing interests

MAC, AP, SB and JT are employees of bioMérieux. bioMérieux shares with Hospices Civils de Lyon the intellectual property of two patents related to the use of CX3CR1 as a monitoring tool in ICU patients. MAC, AP, GM and AL are coinventors on these patents.

Ethical approval and consent to participate

This prospective, multicentre, non-interventional study was approved by our institutional ethical review board (Comité d’Ethique des Centres d’investigation Clinique de l’Inter-Région Rhône-Alpes Auvergne – IRB 5044). No informed consent was needed, as this study was non-interventional and complementary analysis was done on residual blood samples after completion of routine follow-up tests. Nevertheless, non-opposition to inclusion in the protocol was recorded for every patient or the patient’s next of kin.

Author information

Authors and Affiliations

Authors

Consortia

Corresponding author

Correspondence to Fabienne Venet.

Additional information

Arnaud Friggeri and Marie-Angélique Cazalis as Co-first authors.

Additional files

Additional file 1: Table S1.

Primer designs. (DOC 29 kb)

Additional file 2: Figure S1.

Survival curve of the total cohort of intensive care unit patients. (PPT 83 kb)

Additional file 3: Table S2.

Univariate and multivariate analyses of mortality according to CX3CR1 mRNA expression in critically ill patients. A- D7 mortality and CX3CR1 mRNA expression at D1. B- D28 mortality and CX3CR1 mRNA expression at D3. (DOC 53 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Friggeri, A., Cazalis, MA., Pachot, A. et al. Decreased CX3CR1 messenger RNA expression is an independent molecular biomarker of early and late mortality in critically ill patients. Crit Care 20, 204 (2016). https://doi.org/10.1186/s13054-016-1362-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13054-016-1362-x

Keywords