Clinical phenotyping uncovers heterogeneous associations between corticosteroid treatment and survival in critically ill COVID-19 patients
Niklas Bruse, Ana Motos, Rombout Van Amstel, Eckart De Bie, Emma J. Kooistra, Aron Jansen, Dirk van Lier, Jason Kennedy, Daniel Schwarzkopf, Daniel Thomas‐Rüddel, Jesús F. Bermejo-Martín, Ferrán Barbé, Nicolette F. de Keizer, Michael Bauer, Johannes G. van der Hoeven, Antoní Torres, Christopher W. Seymour, Lonneke van Vught, Peter Pickkers, Matthijs Kox
Abstract
PURPOSE: Disease heterogeneity in coronavirus disease 2019 (COVID-19) may render the current one-size-fits-all treatment approach suboptimal. We aimed to identify and immunologically characterize clinical phenotypes among critically ill COVID-19 patients, and to assess heterogeneity of corticosteroid treatment effect. METHODS: We applied consensus k-means clustering on 21 clinical parameters obtained within 24 h after admission to the intensive care unit (ICU) from 13,279 COVID-19 patients admitted to 82 Dutch ICUs from February 2020 to February 2022. Derived phenotypes were reproduced in 6225 COVID-19 ICU patients from Spain (February 2020 to December 2021). Longitudinal immunological characterization was performed in three COVID-19 ICU cohorts from the Netherlands and Germany, and associations between corticosteroid treatment and survival were assessed across phenotypes. RESULTS: ratio, and a 90-day in-hospital mortality rate of 18%. COVIDICU2 patients (37%) had the lowest BMI, were older and had higher APACHE scores and mortality rate (24%) than COVIDICU1. Patients with COVIDICU3 (20%) were the eldest with the most comorbidities, the highest APACHE scores, acute kidney injury and metabolic dysregulations, and the highest mortality rate (47%). These patients also displayed the most pronounced inflammatory response. Corticosteroid therapy started at day 5 [2-9] after ICU admission and administered for 5 [3-7] days was associated with an increased risk for 90-day mortality in patients with the COVIDICU1 and COVIDICU2 phenotypes (hazard ratio [HR] 1.59 [1.09-2.31], p = 0.015 and HR 1.79 [1.42-2.26], p < 0.001, respectively), but not in patients with the COVIDICU3 phenotype (HR 1.08 [0.76-1.54], p = 0.654). CONCLUSION: Our multinational study identified three distinct clinical COVID-19 phenotypes, each exhibiting marked differences in demographic, clinical, and immunological features, and in the response to late and short-term corticosteroid treatment.