- Research
- Open access
- Published:
Inflammatory subphenotypes previously identified in ARDS are associated with mortality at intensive care unit discharge: a secondary analysis of a prospective observational study
Critical Care volume 28, Article number: 151 (2024)
Abstract
Background
Intensive care unit (ICU)-survivors have an increased risk of mortality after discharge compared to the general population. On ICU admission subphenotypes based on the plasma biomarker levels of interleukin-8, protein C and bicarbonate have been identified in patients admitted with acute respiratory distress syndrome (ARDS) that are prognostic of outcome and predictive of treatment response. We hypothesized that if these inflammatory subphenotypes previously identified among ARDS patients are assigned at ICU discharge in a more general critically ill population, they are associated with short- and long-term outcome.
Methods
A secondary analysis of a prospective observational cohort study conducted in two Dutch ICUs between 2011 and 2014 was performed. All patients discharged alive from the ICU were at ICU discharge adjudicated to the previously identified inflammatory subphenotypes applying a validated parsimonious model using variables measured median 10.6 h [IQR, 8.0–31.4] prior to ICU discharge. Subphenotype distribution at ICU discharge, clinical characteristics and outcomes were analyzed. As a sensitivity analysis, a latent class analysis (LCA) was executed for subphenotype identification based on plasma protein biomarkers at ICU discharge reflective of coagulation activation, endothelial cell activation and inflammation. Concordance between the subphenotyping strategies was studied.
Results
Of the 8332 patients included in the original cohort, 1483 ICU-survivors had plasma biomarkers available and could be assigned to the inflammatory subphenotypes. At ICU discharge 6% (n = 86) was assigned to the hyperinflammatory and 94% (n = 1397) to the hypoinflammatory subphenotype. Patients assigned to the hyperinflammatory subphenotype were discharged with signs of more severe organ dysfunction (SOFA scores 7 [IQR 5–9] vs. 4 [IQR 2–6], p < 0.001). Mortality was higher in patients assigned to the hyperinflammatory subphenotype (30-day mortality 21% vs. 11%, p = 0.005; one-year mortality 48% vs. 28%, p < 0.001). LCA deemed 2 subphenotypes most suitable. ICU-survivors from class 1 had significantly higher mortality compared to class 2. Patients belonging to the hyperinflammatory subphenotype were mainly in class 1.
Conclusions
Patients assigned to the hyperinflammatory subphenotype at ICU discharge showed significantly stronger anomalies in coagulation activation, endothelial cell activation and inflammation pathways implicated in the pathogenesis of critical disease and increased mortality until one-year follow up.
Background
Critically ill patients who survive the acute phase of their disease show persistent cognitive, physical and functional impairment, leading to lower quality of life, increased readmission rates and an enhanced risk of mortality [1]. Intensive care unit (ICU) survivors use more healthcare resources after hospital discharge compared to all hospitalized patients, and one-third of ICU patients discharged alive die within 5 years after discharge [2]. During and after ICU admission, especially for sepsis, profound and persistent alteration of the immune response in survivors has been observed [3] and this sustained dysregulation of the host immune response after hospital discharge is hypothesized to contribute to long-term sequelae [4]. Worsened physical, mental and neurocognitive status in ICU-survivors is called post-intensive care syndrome (PICS), and is not only related to the disease severity during ICU admission [1], calling for a more personalized approach.
Substantial progress has been made in the identification of subphenotypes of critical care syndromes, which have been prognostic of outcome and predictive of treatment response [5]. Subphenotypes have been identified for acute respiratory distress syndrome (ARDS) [6, 7], sepsis [8,9,10] and acute kidney injury (AKI) [11, 12], resulting in new clinical, biological and predictive insights. ARDS patients have been divided in a hyperinflammatory and a hypoinflammatory subphenotype using a latent class analysis (LCA) with clinical and plasma protein biomarker data [6]. These subphenotypes were identified within 36 h after the diagnosis of ARDS, and were found to respond differently to several treatments including administration of simvastatin, fluid management and positive and-expiratory pressure (PEEP) mechanical ventilation strategy in retrospective analyses [6, 13, 14]. These subphenotypes might not be disease specific; extending these ARDS subphenotypes to mechanically ventilated critically ill patients without ARDS [15] and in patients with sepsis [16] showed similar results, suggesting that these subphenotypes might be a universal predictor for outcome and treatment response in critically ill patients.
To date, only scarce data exist on patient stratification at ICU discharge. Whether these subphenotypes show association with short or long-term outcome is yet to be determined. We hypothesized that if these inflammatory subphenotypes previously identified among ARDS patients are assigned at ICU discharge, they are associated with worse clinical short and long-term outcomes and a more dysregulated host response. The present study seeks to identify subphenotypes in a more broadly defined population of critically ill ICU-survivors, consisting of sepsis, ARDS and non-infectious patients, at ICU discharge and investigates the association of these subphenotypes with one-year mortality.
Methods
Study design and patient selection
This study was conducted as a secondary analysis of the Molecular Diagnosis and Risk Stratification of Sepsis (MARS) cohort, a prospective observational cohort study in the mixed ICUs of 2 tertiary teaching hospitals (Amsterdam Universal Medical Centers, location AMC, in Amsterdam and University Medical Center Utrecht in Utrecht) [17]. The MARS study included all patients admitted to the ICU above 18 years of age between January 2011 and January 2014 with an expected length of stay longer than 24 h. Patients were included via an opt-out method approved by the medical ethical committees of the participating hospitals. Information on demographics, comorbidities, daily clinical, laboratory and outcome data were prospectively collected by trained and dedicated ICU researchers [18]. The diagnosis on admission was also collected. Patients were either classified as having sepsis or a noninfectious admission diagnosis. Sepsis was defined as the presence of an infection diagnosed within 24 h after ICU admission with a probable or definite likelihood [19] accompanied by a Sequential Organ Failure Assessment (SOFA) scores ≥ 2 according to the sepsis-3 criteria [20]. In patients with a noninfectious admission diagnosis (including those admitted with suspected infection but with a post hoc infection likelihood of none), the Acute Physiology and Chronic Health Evaluation IV (APACHE IV) primary admission diagnosis was used to classify patients [21]. In a retrospectively selected subset of the patients admitted between January 2011 and July 2013 plasma protein biomarkers were measured. This subset consisted of patients with an infection with a probable or definite likelihood, patients with ARDS or patients selected as non-infectious ICU controls [17, 22, 23] (infection likelihood criteria are described in [18], ARDS criteria in [24]). For more information of the selection of this subset see the Additional file 1: Supplementary methods. Patients were selected (1) who were discharged alive from the ICU, and (2) of whom plasma protein biomarkers were available and collected on the day of ICU discharge, or one or two day(s) before. ICU readmissions within the same hospital admission were excluded from the analyses. However, ICU readmissions were collected as outcome measurements reported in two ways; ICU readmissions within the same hospital admission and all ICU readmissions during the remaining study duration of MARS. Baseline characteristics, comorbidities, outcome data and severity indices such as APACHE IV [21] and SOFA scores [20] are reported. For information on the definitions of comorbidities and complications, see the Additional file 1: Supplementary methods.
Biomarker measurements
Host response biomarkers reflective of coagulation-, endothelial and inflammatory pathways were measured in ethylenediaminetetraacetic acid (EDTA) anticoagulated plasma, collected within 16 h after ICU admission, on day 2, on day 4 and on ICU discharge (-1 or -2 days). Biomarker panels consisted of angiopoietin-1, angiopoietin-2, D-dimer, interferon-γ, interleukin (IL)-6, IL-8, IL-10, fractalkine, matrix metalloproteinase (MMP)-8, protein C, soluble E-selectin and soluble intercellular adhesion molecule-1 (ICAM-1). C-reactive protein (CRP), creatinine lactate, platelets, prothrombin time (PT), bicarbonate and white blood cell count (WBC), amongst others, were measured for clinical purposes. For further information on sample handling and assays, see the Additional file 1: Supplementary methods.
Adjudication to subphenotypes
The probability of belonging to the hyperinflammatory or hypoinflammatory subphenotype [6] was estimated using a previously published and validated parsimonious 3-variable model using plasma levels of bicarbonate, IL-8 and Protein C at ICU discharge [25]. To assure a reliable subphenotype adjudication plasma bicarbonate results were aligned with the day of IL-8 and Protein C measurements. When the bicarbonate value was missing on the day of the IL-8 and protein C measurement, (1) the average of the bicarbonate value on the day before and the day after was used, or (2) if there was only one value available (on the day before or after), this value was carried forward or backward, or (3) when both values on the day before and after were missing, the patient was excluded. A probability of > 0.5 was used to assign patients to the hyperinflammatory subphenotype, and < 0.5 to the hypoinflammatory subphenotype.
Outcome measurements
The study aim was to analyze associations between the biological subphenotypes and the primary and secondary outcomes. The primary outcomes were short- and long-term mortality (30-day, 90-day mortality and one-year mortality) after ICU discharge. The secondary outcomes were ICU readmission after ICU discharge and biomarkers reflective of the coagulation-, endothelial and inflammatory pathways were compared between the different subphenotypes at ICU discharge.
Statistical analysis
Patient characteristics, plasma protein biomarker concentrations, and outcomes were compared using a t-test or one-way ANOVA for parametric data, a Mann–Whitney U test or Kruskal–Wallis test for nonparametric data, and with Chi-square test for categorical data, stratified by subphenotype. Subphenotype consistency was analyzed by adjudication of subphenotypes on ICU admission and readmission and assessing subphenotype overlap. The percentage of the patients in the subphenotypes was compared at different time points (ICU admission vs. ICU discharge, ICU discharge vs. ICU readmission), visualized with mosaic plots. Differences in individual plasma protein biomarker levels between the previously described subphenotypes were quantified and expressed as Hedges’ g, a commonly used effect size measure [26]. Survival was visualized using Kaplan–Meier curves and analyzed using a Cox proportional hazards model. Previously identified risk factors for long-term outcome after ICU survival were included in the Cox regression models [27]. These risks factors consisted of age, Charlson comorbidity index [28], chronic cardiovascular insufficiency, chronic heart failure, chronic kidney disease, malignancy, systolic blood pressure, temperature, platelets, white blood cell count and length of ICU stay at ICU discharge. As a sensitivity analysis to investigate the influence of admission diagnoses and complications, the adjusted Cox regression analysis was repeated with the following factors also included: complications upon ICU admission (ARDS, sepsis and AKI), admission diagnosis group (cardiovascular, gastrointestinal, metabolic, neurological, other, respiratory and trauma) and admission type (medical or surgical). A p value of < 0.05 was considered of statistical significance.
A sensitivity analysis was conducted using a LCA based on plasma protein biomarker concentrations in the coagulation, endothelial and inflammatory pathways (plasma protein biomarkers that were not part of the daily clinical laboratory measurements) at ICU discharge to identify classes of ICU survivors. The previously identified subphenotypes in ARDS patients were based on clinical and biological data, for this sensitivity analysis we aimed to analyze if similar patterns would emerge by only using non-daily clinical laboratory measurements. The process of model design and LCA followed the steps and considerations outlined by Sinha et al. [29]. First, missing plasma protein biomarkers used in the LCA were imputed using multivariate imputation by chained equations algorithm (see the Additional file 1: Supplementary methods) and transformed to resemble normally distributed data, which was verified using Shapiro–Wilk tests and histograms. If plasma protein biomarkers had a Spearman correlation coefficient > 0.6, one of the two was excluded from the LCA. Plasma biomarker levels were scaled by subtracting the mean in this cohort and dividing by the standard deviation. Five sequential models consisting of 1 to 5 classes were fit. As described [29] the best-fitting model was selected based on the Bayesian information criterion, the Vuong-Lo-Mendell-Rubin test [30] and entropy. Subphenotype characteristics were visualized using a profile plot, which displayed the mean standardized differences of subphenotype defining variables. To conduct LCA Mplus was used. Once the best model was selected, an individual patient’s class assignment was determined by the highest probability of class membership. The percentage of the patients in the subphenotypes was compared between the two clustering strategies at ICU discharge (adjudication of previously described subphenotypes vs. LCA), visualized with a mosaic plot.
Results
Population
Of the 8332 patients included in the original cohort, 1483 ICU survivors were selected for further subphenotyping (Fig. 1). The included patients had a median age of 62 [IQR 50–71] years and 61% was male (n = 907) (Table 1). On ICU admission, 292 (20%) patients had ARDS and 927 (63%) sepsis, of whom 257 (17%) had septic shock and 252 (17%) had both sepsis and ARDS (Table 1). The most common comorbidities were cardiovascular compromise (n = 440, 30%), hypertension (n = 438, 30%), immunocompromise (n = 295, 20%) and diabetes (n = 288, 19%) (Table 1). For the most common diagnoses on admission see Additional file 1: Supplementary Table 1. The APACHE IV score and the SOFA score upon ICU admission were 70 [IQR 54–90] and 6 [IQR 4–8] respectively.
Inflammatory subphenotypes previously identified among ARDS patients
In 959 (65%) of the included patients plasma biomarker levels were measured on the day of ICU discharge, in 419 (28%) the day before ICU discharge and in 105 (7%) patients 2 days before ICU discharge; this was approximately equally distributed in both subphenotypes (Additional file 1: Supplementary Table 2). For 307 (21%) this entailed that the plasma biomarker levels used for this analysis were taken upon ICU admission since no later sample was available. At ICU discharge the majority of patients (n = 1397, 94%) were adjudicated to the hypoinflammatory subphenotype whereas the hyperinflammatory subphenotype consisted of 86 patients (6%) (Fig. 1). Patients adjudicated to the hyperinflammatory subphenotype were discharged with more severe disease reflected by higher SOFA scores (7 [IQR 5–9] vs. 4 [IQR 2–6], p < 0.001). The number of ICU-acquired complications (shock, AKI, ARDS and sepsis) that had occurred prior to discharge was comparable between the two subphenotypes (Table 2). Patients within the hyperinflammatory subphenotype at discharge had been more severely ill upon admission (APACHE IV score 73 vs. 70, p = 0.044 and SOFA score 7 vs. 6, p = 0.027), were more often admitted with sepsis (79% vs. 62%, p = 0.002) and less often admitted with ARDS (7% vs. 21%, p = 0.004) (Table 1).
Subphenotype overlap at ICU discharge versus ICU admission
Of 1309 patients (88%) in our cohort, the corresponding subphenotypes could be adjudicated on admission (day 0 or 1). In this subset, of the patients with the hypoinflammatory subphenotype at ICU discharge (n = 1227), 1034 (84%) were also adjudicated to the hypoinflammatory subphenotype on ICU admission. In total, 193 patients (16%) changed from the hyperinflammatory at ICU admission to the hypoinflammatory subphenotype at ICU discharge. Of the patients with the hyperinflammatory subphenotype at ICU discharge (n = 83), 76 (92%) were also adjudicated to the hyperinflammatory subphenotype at ICU admission and only 6 (7%) came from the hypoinflammatory subphenotype at ICU admission (Fig. 2).
Outcome measurements
Mortality was significantly higher in ICU-survivors at ICU discharge adjudicated to the hyperinflammatory subphenotype compared to the hypoinflammatory subphenotype (30-day mortality 21% vs. 11%, p = 0.005; one-year mortality 48% vs. 28%, p < 0.001, Table 2 and Fig. 3). In a Cox proportional hazards model, adjusted for 11 variables [27], membership to the hyperinflammatory subphenotype at ICU discharge was independently associated with one-year mortality (adjusted hazard ratio (HR) = 2.02 (95% CI 1.37–2.96); p < 0.001). A similar result was found when, as a sensitivity analysis, the Cox proportional hazards model was expanded with complications upon ICU admission, admission diagnosis and admission type (HR = 2.10 (95% CI 1.41–3.13); p < 0.001). In this model age, Charlson comorbidity index, chronic cardiovascular insufficiency, malignancy, platelets at ICU discharge, the diagnosis on admission ‘other’ and the admission type ‘surgical’ were all significant associated with one-year mortality. For more information on the model see Additional file 1: Supplementary Table 3. In the hyperinflammatory subphenotype the ICU readmission rate during the remaining study duration of MARS between January 2011 and January 2014 (36% vs. 22%, p = 0.006) was significantly higher compared to the hypoinflammatory subphenotype (Table 2). In contrast, the length of ICU stay was shorter in the hyperinflammatory subphenotype (2 vs. 4 days, p < 0.001). ICU readmission within the same hospital admission was 24% in the hyperinflammatory vs. 16% in the hypoinflammatory subphenotype (p = 0.054). The differences in host response biomarkers between the inflammatory subphenotypes at ICU discharge can be found in the Additional file 1: Supplementary Results, Supplementary Table 4 and Supplementary Fig. 1.
Readmitted patients within the same hospital admission
In 158 (65%) of 243 patients who were readmitted to the ICU within the same hospital admission, plasma biomarkers were also measured on their corresponding ICU readmission. These 158 patients who were readmitted to the ICU had a SOFA score of 4 [IQR 3–7] at ICU discharge, had a median age of 63 [IQR 54–72] years and 65% (n = 103) was male. The median duration between the admission and readmission was 4 days [IQR 2–10].
In this subset, of the patients with the hypoinflammatory subphenotype at ICU discharge (n = 145), 122 (84%) were also adjudicated to the hypoinflammatory subphenotype at ICU readmission. 23 (16%) of the patients discharged with the hypoinflammatory subphenotype changed to the hyperinflammatory subphenotype at ICU readmission. Of the patients discharged with the hyperinflammatory subphenotype (n = 13), 8 (62%) were also adjudicated to the hyperinflammatory subphenotype at ICU readmission. 5 (38%) patients discharged as hyperinflammatory changed to the hypoinflammatory subphenotype at ICU readmission (Fig. 4).
Latent class analysis
The following thirteen plasma biomarkers measured at ICU discharge were considered as variables for the LCA: angiopoietin-1, angiopoietin-2, ratio angiopoietin-1/angiopoietin-2, D-dimer, E-Selectin, fractalkine, ICAM-1, INF-γ, IL-6, IL-8, IL-10, MMP-8 and protein C. Three variables, E-Selectin, ICAM-1 and fractalkine, had a low missing percentage (0.2%) before imputation, the other variables used in the LCA were complete. After excluding correlated variables (ratio angiopoietin-1/angiopoietin-2 and IL-8) and excluding angiopoietin-2, see the Additional file 1: Supplementary Results for the rationale, 10 variables remained (Fig. 5A; Additional file 1: Supplementary Fig. 2). Additional file 1: Supplementary Table 5 shows model-fitting statistics for LCA models consisting of 1 to 5 classes. The Vuong-Lo-Mendell-Rubin test showed a p < 0.001 for 2, 4 and 5 classes indicating that these class provided an improved fit. The entropy for all the number of classes was 0.75–0.79 and just did not reach the preferable > 0.8. The Bayesian information criterion (BIC) was the lowest for 5 classes, but with 3 or more classes the BIC did not lower significantly and so the increase of the model complexity of adding a third or more classes was considered unnecessary. Therefore, a 2-class latent model was deemed most suitable. The average latent class probability was calculated, which was 90.9% for the whole cohort; 87.8% for class 1 and 94.1% for class 2. The percentage of patients with a probability higher than 90% was 76.1%. In a sensitively analyses, in which the model was run with all five imputation datasets, outcomes and characteristics showed very similar results (data not shown), therefore an imputation dataset was chosen at random.
In the 2-class LCA model 349 (24%) patients were assigned to class 1 and 1134 (76%) to class 4 (Table 3). Class 1 was defined by a more dysregulated inflammatory profile, with the most distinct differences in IL-10 and IL-6 compared to class 2. Furthermore, class 1 was defined by higher endothelial markers, with the most distinct differences in ICAM-1 and E-Selectin. In class 2 patients were significantly older (63 years vs. 61 years, p = 0.010) and more often male (63% vs. 55%, p = 0.012, Table 3). Patients in class 1 were most severely ill as reflected by higher SOFA scores at ICU discharge (median 6 [IQR 4–8] vs. 3 [IQR 2–5], p < 0.001). More differences in medical history, diagnosis and disease severity on admission between patients in the four classes can be found in Additional file 1: Supplementary Table 6. Comparing biomarker levels at ICU-discharge, these results reflected the more dysregulated host response in patients in class 1 compared to patients in class 2, with among others, higher levels of CRP, IL-6, -8, -10, creatinine and lactate (Additional file 1: Supplementary Fig. 3). Mortality was significantly higher in ICU-survivors adjudicated to class 1 compared to class 2 (30-day mortality 17% vs. 10%, p = 0.001 and one-year mortality 41% vs. 25%, p < 0.001, Table 3 and Fig. 5B). In a Cox proportional hazards model, adjusted for 11 variables [29], membership to class 1 at ICU discharge was independently associated with one-year mortality (adjusted hazard ratio (HR) = 1.86 (95% CI 1.45–2.39); p < 0.001). In class 1 the ICU readmission rate was significantly higher (25% vs. 14%, p < 0.001) compared to class 2 (Table 3). In contrast, the length of ICU stay was shorter in class 1 (3 vs. 4 days, p < 0.001).
When assessing phenotype concordance; patients in the hyperinflammatory subphenotype almost completely belonged to class 1 (n = 79, 92%), only 7 (8%) belonged to class 2. Patients in the hypoinflammatory subphenotype most often belonged to class 2 (n = 1127, 81%), 270 (19%) belonged to class 1 (Fig. 5C), demonstrating a similar pattern in subphenotype adjudication with the largest differences being the distribution of the patients (6% in the hyperinflammatory subphenotype vs. 24% in class 1; 94% in the hyperinflammatory subphenotype vs. 76% in class 2).
Discussion
This study is the first to demonstrate that the inflammatory subphenotypes previously identified among ARDS patients at ICU discharge are associated with differences in short and long-term outcome. ICU survivors adjudicated to these inflammatory subphenotypes at ICU discharge had an increased mortality. Patients in this hyperinflammatory subphenotype and in the LCA class with the more dysregulated inflammatory profile (class 1) had more derailed markers of coagulation activation, endothelial cell activation and inflammation at ICU discharge.
These findings are in line with previous studies showing that persistent hyperinflammation at either ICU or hospital discharge in ICU survivors is associated with poor physical recovery [31] and higher mortality after hospital and ICU discharge in sepsis-survivors [32, 33]. Furthermore, patients in the hyperinflammatory subphenotype were discharged with significantly higher creatinine levels, which is in line with an earlier study demonstrating that elevated kidney biomarkers measured at ICU discharge in the general ICU population were associated with a poor 1-year outcome [34]. Persistent hyperinflammation at ICU discharge, as found in this study, or even after hospital discharge, plays a detrimental role in outcome [35]. Increased signs of inflammation, measured by CRP, is associated with poor physical recovery during the first 3 months post-ICU discharge [31]. Higher levels of CRP during rehabilitation after hospital admission was furthermore associated with institutionalization and mortality in geriatric patients [36]. Inflammation is a risk factor for physiologic alterations involved in persistent ICU-acquired weakness [37] and might play a role in the development of cognitive disorders by inducing neurotoxicity, endothelial injury and blood–brain barrier dysfunction [38]. The underlying mechanisms linking the acute critical illness to long lasting impairment are however not yet completely understood.
Long lasting impairment after acute critical illness is a frequently occurring health-problem. In a prospective multicenter cohort study conducted in 2345 patients, 50% was suffering from new physical, mental, and/or cognitive health problems 1 year after ICU admission; with 49% having physical impairment, 14% mental impairment, and 4% cognitive impairment [39]. Moreover, the general ICU population has an increased mortality risk after hospital discharge when compared to age and sex-matched population controls, with a mortality rate of 27.5% 3 years after hospital discharge versus 8.2% at 3 years for the general population [40]. There is limited beneficial effect of clinical studies investigating physical rehabilitation [35]. Interventions studied, such as home-based rehabilitation and follow-up clinics, did not yet improve outcome, which raises the question whether the target populations of the studied interventions were appropriate [35]. Using subphenotyping as an enrichment strategy could enhance the likelihood of detecting prevention and therapeutic strategies for specific groups of ICU survivors. Identification of ICU-survivors most at risk for adverse outcomes, could provide information on treatable traits that associate with readmission or death in this group of patients [41]. As an example of a targeted therapeutic strategy, in patients discharged alive from the ICU after experiencing AKI, angiotensin-converting enzyme inhibitor or receptor blocker use after discharge is associated with increased survival [42].
Notably, patients within the hyperinflammatory subphenotype as well as patients within the LCA class 1 with the more dysregulated inflammatory profile were discharged from the ICU significantly quicker, which could not be explained by patients who died during ICU stay since they were excluded from the analysis. This finding seems unrelated to clustering methods as in both techniques the patients with the most derailed host response biomarkers were discharged from the ICU more rapidly compared to the other group. A striking finding, since in other studies a longer ICU stay was associated with higher one-year mortality [27, 43]. The explanation for this shorter length of ICU-stay in patients who were discharged from the ICU with the hyperinflammatory subphenotype, could be related to their trajectory of recovery, i.e. that these patients were discharged to the ward prior to reaching the hypoinflammatory subphenotype. This could also be underlined by the finding of the more derailed host response at ICU discharge and the higher ICU readmission rated in the hyperinflammatory subphenotype. Patients admitted on the ward are less intensively monitored than patients admitted on the ICU, which could make them more susceptible for adverse outcomes. Investigating the value of subphenotyping for assessing the best moment for ICU discharge of individual patients has not been evaluated yet.
A strength of this study is the detailed data collected at ICU discharge including the availability of biomarkers reflective of coagulation activation, endothelial cell activation and inflammation, and the follow-up period of one-year. However, there are several limitations that should be taken into account. First, there might be selection bias in this cohort, as plasma biomarkers were only measured in patients with an infection likelihood of probable or definite, ARDS and matching non-infectious populations, therefore the incidence of specific disease etiologies like sepsis and ARDS might not be representative for the whole ICU population; the mortality might be higher in this cohort compared to a unselected ICU population. Second, the study was not externally validated in another cohort. Detailed information at ICU discharge is often lacking in other datasets, therefore external validation was not possible. Third, since this study has an observational design, any firm conclusions on casual relationships between subphenotypes and one-year mortality cannot be drawn. Last, the LCA model did not reach an entropy of 0.8, which is probably explained by the fact that the class separation at ICU discharge is not as clear as it would have been upon ICU admission. This could be expected, since patients assessed as ready for ICU discharge are probably a more homogenous group than compared to ICU admission.
Conclusion
In conclusion, this study shows an association between the hyperinflammatory subphenotype and LCA class 1 with the more dysregulated inflammatory profile upon ICU discharge and increased mortality during the subsequent year. This might be linked to the more derailed host response that was found in the hyperinflammatory subphenotype and the LCA class 1 at ICU discharge. The use of subphenotyping as an enrichment strategy could enhance the likelihood of detecting prevention and therapeutic strategies for the ICU survivors most at risk for a poor outcome.
Availability of data and materials
The data of our cohort (MARS) is available at the Gene Expression Omnibus public repository of NCBI under accession number GSE65682.
References
Rousseau AF, Prescott HC, Brett SJ, Weiss B, Azoulay E, Creteur J, et al. Long-term outcomes after critical illness: recent insights. Crit Care. 2021;25(1):108.
Hill AD, Fowler RA, Pinto R, Herridge MS, Cuthbertson BH, Scales DC. Long-term outcomes and healthcare utilization following critical illness–a population-based study. Crit Care. 2016;20:76.
Voiriot G, Oualha M, Pierre A, Salmon-Gandonniere C, Gaudet A, Jouan Y, et al. Chronic critical illness and post-intensive care syndrome: from pathophysiology to clinical challenges. Ann Intensive Care. 2022;12(1):58.
Prescott HC, Angus DC. Enhancing recovery from sepsis: a review. JAMA. 2018;319(1):62–75.
Reddy K, Sinha P, O’Kane CM, Gordon AC, Calfee CS, McAuley DF. Subphenotypes in critical care: translation into clinical practice. Lancet Respir Med. 2020;8(6):631–43.
Calfee CS, Delucchi K, Parsons PE, Thompson BT, Ware LB, Matthay MA, Network NA. Subphenotypes in acute respiratory distress syndrome: latent class analysis of data from two randomised controlled trials. Lancet Respir Med. 2014;2(8):611–20.
Bos LD, Schouten LR, van Vught LA, Wiewel MA, Ong DSY, Cremer O, et al. Identification and validation of distinct biological phenotypes in patients with acute respiratory distress syndrome by cluster analysis. Thorax. 2017;72(10):876–83.
Scicluna BP, van Vught LA, Zwinderman AH, Wiewel MA, Davenport EE, Burnham KL, et al. Classification of patients with sepsis according to blood genomic endotype: a prospective cohort study. Lancet Respir Med. 2017;5(10):816–26.
Davenport EE, Burnham KL, Radhakrishnan J, Humburg P, Hutton P, Mills TC, et al. Genomic landscape of the individual host response and outcomes in sepsis: a prospective cohort study. Lancet Respir Med. 2016;4(4):259–71.
Seymour CW, Kennedy JN, Wang S, Chang CH, Elliott CF, Xu Z, et al. Derivation, validation, and potential treatment implications of novel clinical phenotypes for sepsis. JAMA. 2019;321(20):2003–17.
Wiersema R, Jukarainen S, Vaara ST, Poukkanen M, Lakkisto P, Wong H, et al. Two subphenotypes of septic acute kidney injury are associated with different 90-day mortality and renal recovery. Crit Care. 2020;24(1):150.
Bhatraju PK, Zelnick LR, Herting J, Katz R, Mikacenic C, Kosamo S, et al. Identification of acute kidney injury subphenotypes with differing molecular signatures and responses to vasopressin therapy. Am J Respir Crit Care Med. 2019;199(7):863–72.
Famous KR, Delucchi K, Ware LB, Kangelaris KN, Liu KD, Thompson BT, et al. Acute respiratory distress syndrome subphenotypes respond differently to randomized fluid management strategy. Am J Respir Crit Care Med. 2017;195(3):331–8.
Calfee CS, Delucchi KL, Sinha P, Matthay MA, Hackett J, Shankar-Hari M, et al. Acute respiratory distress syndrome subphenotypes and differential response to simvastatin: secondary analysis of a randomised controlled trial. Lancet Respir Med. 2018;6(9):691–8.
Heijnen NFL, Hagens LA, Smit MR, Cremer OL, Ong DSY, van der Poll T, et al. Biological subphenotypes of acute respiratory distress syndrome show prognostic enrichment in mechanically ventilated patients without acute respiratory distress syndrome. Am J Respir Crit Care Med. 2021;203(12):1503–11.
Sinha P, Kerchberger VE, Willmore A, Chambers J, Zhuo H, Abbott J, et al. Identifying molecular phenotypes in sepsis: an analysis of two prospective observational cohorts and secondary analysis of two randomised controlled trials. Lancet Respir Med. 2023;11(11):965–74.
van Vught LA, Wiewel MA, Hoogendijk AJ, Frencken JF, Scicluna BP, Klein Klouwenberg PMC, et al. The host response in patients with sepsis developing intensive care unit-acquired secondary infections. Am J Respir Crit Care Med. 2017;196(4):458–70.
Klein Klouwenberg PM, Ong DS, Bos LD, de Beer FM, van Hooijdonk RT, Huson MA, et al. Interobserver agreement of centers for disease control and prevention criteria for classifying infections in critically ill patients. Crit Care Med. 2013;41(10):2373–8.
Klein Klouwenberg PM, Cremer OL, van Vught LA, Ong DS, Frencken JF, Schultz MJ, et al. Likelihood of infection in patients with presumed sepsis at the time of intensive care unit admission: a cohort study. Crit Care. 2015;19:319.
Singer M, Deutschman CS, Seymour CW, Shankar-Hari M, Annane D, Bauer M, et al. The third international consensus definitions for sepsis and septic shock (Sepsis-3). JAMA. 2016;315(8):801–10.
Zimmerman JE, Kramer AA, McNair DS, Malila FM. Acute physiology and chronic health evaluation (APACHE) IV: hospital mortality assessment for today’s critically ill patients. Crit Care Med. 2006;34(5):1297–310.
Scicluna BP, Klein Klouwenberg PM, van Vught LA, Wiewel MA, Ong DS, Zwinderman AH, et al. A molecular biomarker to diagnose community-acquired pneumonia on intensive care unit admission. Am J Respir Crit Care Med. 2015;192(7):826–35.
Scicluna BP, Wiewel MA, van Vught LA, Hoogendijk AJ, Klarenbeek AM, Franitza M, et al. Molecular Biomarker to assist in diagnosing abdominal sepsis upon ICU admission. Am J Respir Crit Care Med. 2018;197(8):1070–3.
Bernard GR, Artigas A, Brigham KL, Carlet J, Falke K, Hudson L, Lamy M, Legall JR, Morris A, Spragg R. The American-European Consensus Conference on ARDS Definitions mechanisms relevant outcomes and clinical trial coordination. Am J Respir Crit Care Med. 1994;149(3 Pt 1):818–24.
Sinha P, Delucchi KL, McAuley DF, O’Kane CM, Matthay MA, Calfee CS. Development and validation of parsimonious algorithms to classify acute respiratory distress syndrome phenotypes: a secondary analysis of randomised controlled trials. Lancet Respir Med. 2020;8(3):247–57.
Hedges LV. Distribution theory for Glass’s estimator of effect size and related estimators. J Educ Stat. 1981;6(2):107–28.
Gayat E, Cariou A, Deye N, Vieillard-Baron A, Jaber S, Damoisel C, et al. Determinants of long-term outcome in ICU survivors: results from the FROG-ICU study. Crit Care. 2018;22(1):8.
Charlson ME, Pompei P, Ales KL, MacKenzie CR. A new method of classifying prognostic comorbidity in longitudinal studies: development and validation. J Chron Dis. 1987;40(5):373–83.
Sinha P, Calfee CS, Delucchi KL. Practitioner’s guide to latent class analysis: methodological considerations and common pitfalls. Crit Care Med. 2021;49(1):e63–79.
Nylund KL, Asparouhov T, Muthén BO. Deciding on the number of classes in latent class analysis and growth mixture modeling: A Monte Carlo simulation study. Struct Equ Model. 2007;14(4):535–69.
Griffith DM, Lewis S, Rossi AG, Rennie J, Salisbury L, Merriweather JL, et al. Systemic inflammation after critical illness: relationship with physical recovery and exploration of potential mechanisms. Thorax. 2016;71(9):820–9.
Yende S, D’Angelo G, Kellum JA, Weissfeld L, Fine J, Welch RD, et al. Inflammatory markers at hospital discharge predict subsequent mortality after pneumonia and sepsis. Am J Respir Crit Care Med. 2008;177(11):1242–7.
Soussi S, Sharma D, Juni P, Lebovic G, Brochard L, Marshall JC, et al. Identifying clinical subtypes in sepsis-survivors with different one-year outcomes: a secondary latent class analysis of the FROG-ICU cohort. Crit Care. 2022;26(1):114.
Legrand M, Hollinger A, Vieillard-Baron A, Depret F, Cariou A, Deye N, et al. One-year prognosis of kidney injury at discharge from the ICU: A multicenter observational study. Crit Care Med. 2019;47(12):e953–61.
Elliott D, Puthucheary Z. Is rehabilitation post critical illness a new anti-inflammatory agent? Thorax. 2016;71(9):783–4.
Guan L, Tuttle CSL, Reijnierse EM, Lim WK, Maier AB. Unresolved inflammation during hospitalization is associated with post-discharge institutionalization and mortality in geriatric rehabilitation inpatients: the RESORT cohort. Exp Gerontol. 2021;156: 111597.
Langhans C, Weber-Carstens S, Schmidt F, Hamati J, Kny M, Zhu X, et al. Inflammation-induced acute phase response in skeletal muscle and critical illness myopathy. PLoS ONE. 2014;9(3): e92048.
Sasannejad C, Ely EW, Lahiri S. Long-term cognitive impairment after acute respiratory distress syndrome: a review of clinical impact and pathophysiological mechanisms. Crit Care. 2019;23(1):352.
Geense WW, Zegers M, Peters MAA, Ewalds E, Simons KS, Vermeulen H, et al. New physical, mental, and cognitive problems 1 year after ICU admission: a prospective multicenter study. Am J Respir Crit Care Med. 2021;203(12):1512–21.
Brinkman S, de Jonge E, Abu-Hanna A, Arbous MS, de Lange DW, de Keizer NF. Mortality after hospital discharge in ICU patients. Crit Care Med. 2013;41(5):1229–36.
Brown SM, Bose S, Banner-Goodspeed V, Beesley SJ, Dinglas VD, Hopkins RO, et al. Approaches to addressing post-intensive care syndrome among intensive care unit survivors: a narrative review. Ann Am Thorac Soc. 2019;16(8):947–56.
Gayat E, Hollinger A, Cariou A, Deye N, Vieillard-Baron A, Jaber S, et al. Impact of angiotensin-converting enzyme inhibitors or receptor blockers on post-ICU discharge outcome in patients with acute kidney injury. Intensive Care Med. 2018;44(5):598–605.
Moitra VK, Guerra C, Linde-Zwirble WT, Wunsch H. Relationship between ICU length of stay and long-term mortality for elderly ICU survivors. Crit Care Med. 2016;44(4):655–62.
Acknowledgements
MARS consortium: Amsterdam UMC location University of Amsterdam: Friso M. de Beer, Lieuwe D. J. Bos, Gerie J. Glas, Arie J. Hoogendijk, Roosmarijn T. M. van Hooijdonk, Janneke Horn, Mischa A. Huson, Laura R. A. Schouten, Marcus J. Schultz, Brendon P. Scicluna, Marleen Straat, Lonneke A. van Vught, Luuk Wieske, Maryse A. Wiewel, Esther Witteveen. University Medical Center Utrecht: Marc J.M. Bonten, Olaf M. Cremer, David S.Y. Ong, Jos F. Frencken, Peter M.C. Klein Klouwenberg, Maria E. Koster‐Brouwer, Kirsten van de Groep, Diana M. Verboom.
Funding
The MARS project was supported by the Center for Translational Molecular Medicine (http://www.ctmm.nl), project MARS (Grant 04I‐201). MAS was supported by the European Commission (Horizon 2020, ImmunoSep, Grant number 847422). LDJB received grants from Dutch lung foundation (Young investigator grant), Dutch lung foundation and Health Holland (Public–Private Partnership grant), Dutch lung foundation (Dirkje Postma Award), IMI COVID19 initiative, Amsterdam UMC fellowship, ZonMW COVID-19 Urgency grant, ERS Gold Metal for ARDS, Santhera and a VIDI Grant from ZonMW. LAvV was supported by a VENI Grant from ZonMW (Grant number 09150161910033). The funders had no role in the design and conduct of the study: collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.
Author information
Authors and Affiliations
Consortia
Contributions
OLC, TvdP and LAvV designed the study. OLC, TvdP and LAvV acquired the data. MAS drafted the frist manuscript. MAS and RBEvA performed the analyses. All authors provided intellectual input, interpreted the data and revised the manuscript. TvdP and LAvV supervised the project. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
This study was part of the Molecular Diagnosis and Risk Stratifcation of Sepsis (MARS) project (ClinicalTrials.gov identifer NCT0195033). The MARS project was a prospective observational study in two tertiary teaching hospitals in the Netherlands. Inclusion of patients was done using an opt-out method approved by the institutional review boards of both hospitals (IRB No. 10-056C).
Competing interests
LDJB received consultancy fees from Scailyte, Santhera and Janssen en Janssen, and was on the advisory board of Sobi, Exvastat, Pfizer and Astra Zeneca, all paid to the institution. The other authors have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1.
Supplementary methods. Supplementary methods on patient selection, comorbidities, complications, biomarker measurements and imputation. Supplementary results. Supplementary results on host response biomarkers and latent class analysis. Supplementary Table 1. Most common admission diagnoses per admission diagnosis group. Supplementary Table 2. Distribution of moment of sampling per inflammatory subphenotype. Supplementary Table 3. Results of the adjusted Cox proportional hazards model for associations with one-year mortality. Supplementary Table 4. Comparison of biomarkers in the complement-, endothelial and inflammatory pathways at ICU discharge between the hyperinflammatory and hypoinflammatory subphenotype. Supplementary Table 5. Model-fit statistics for different number of latent classes. Supplementary Table 6. Characteristics at ICU admission, cohort divided based on LCA. Supplementary Fig. 1. Biomarker differences at ICU discharge according to inflammatory subphenotypes. Supplementary Fig. 2. Correlation plot with variables intended to use in LCA. Supplementary Fig. 3. Biomarker plots stratified per LCA class.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.
About this article
Cite this article
Slim, M.A., van Amstel, R.B.E., Bos, L.D.J. et al. Inflammatory subphenotypes previously identified in ARDS are associated with mortality at intensive care unit discharge: a secondary analysis of a prospective observational study. Crit Care 28, 151 (2024). https://doi.org/10.1186/s13054-024-04929-9
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13054-024-04929-9