Next Article in Journal
Risk Factors for Perineal Wound Breakdown in Early Postpartum: A Retrospective Case–Control Study
Next Article in Special Issue
Kinetics of C-Reactive Protein and Procalcitonin in the Early Identification of ICU-Acquired Infections in Critically Ill COVID-19 Patients
Previous Article in Journal
Determination of Correlated Color Temperature in Ex Vivo Porcine Eyes during Intraocular Illumination
Previous Article in Special Issue
Mortality, Intensive Care Unit Admission, and Intubation among Hospitalized Patients with COVID-19: A One-Year Retrospective Study in Jordan
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of Distinct Clinical Phenotypes of Critically Ill COVID-19 Patients: Results from a Cohort Observational Study

by
José Pedro Cidade
1,2,*,
Vicente Cés de Souza Dantas
3,
Alessandra de Figueiredo Thompson
4,
Renata Carnevale Carneiro Chermont de Miranda
4,
Rafaela Mamfrim
4,
Henrique Caroli
4,
Gabriela Escudini
4,
Natalia Oliveira
4,
Taiza Castro
3 and
Pedro Póvoa
1,2,5,*
1
Intensive Care Unit 4, Department of Intensive Care São Francisco Xavier Hospital, CHLO, Lisbon, 1449-005 Lisbon, Portugal
2
Nova Medical School, Clinical Medicine, CHRC, New University of Lisbon, 1169-056 Lisbon, Portugal
3
Instituto D’Or de Pesquisa e Ensino, Rio de Janeiro 22281-100, Brazil
4
Hospital Copa D´Or, Rio de Janeiro 22031-011, Brazil
5
Center for Clinical Epidemiology, Research Unit of Clinical Epidemiology, OUH Odense University Hospital, 5000 Odense C, Denmark
*
Authors to whom correspondence should be addressed.
J. Clin. Med. 2023, 12(8), 3035; https://doi.org/10.3390/jcm12083035
Submission received: 9 March 2023 / Revised: 6 April 2023 / Accepted: 20 April 2023 / Published: 21 April 2023
(This article belongs to the Special Issue Intensive Care Management of COVID-19 Patients)

Abstract

:
Purpose: COVID-19 presents complex pathophysiology, and evidence collected points towards an intricate interaction between viral-dependent and individual immunological mechanisms. Identifying phenotypes through clinical and biological markers may provide a better understanding of the subjacent mechanisms and an early patient-tailored characterization of illness severity. Methods: A multicenter prospective cohort study was performed in 5 hospitals in Portugal and Brazil for one year between 2020–2021. All adult patients with an Intensive Care Unit admission with SARS-CoV-2 pneumonia were eligible. COVID-19 was diagnosed using clinical and radiologic criteria with a SARS-CoV-2 positive RT-PCR test. A two-step hierarchical cluster analysis was made using several class-defining variables. Results: 814 patients were included. The cluster analysis revealed a three-class model, allowing for the definition of three distinct COVID-19 phenotypes: 407 patients in phenotype A, 244 patients in phenotype B, and 163 patients in phenotype C. Patients included in phenotype A were significantly older, with higher baseline inflammatory biomarkers profile, and a significantly higher requirement of organ support and mortality rate. Phenotypes B and C demonstrated some overlapping clinical characteristics but different outcomes. Phenotype C patients presented a lower mortality rate, with consistently lower C-reactive protein, but higher procalcitonin and interleukin-6 serum levels, describing an immunological profile significantly different from phenotype B. Conclusions: Severe COVID-19 patients exhibit three different clinical phenotypes with distinct profiles and outcomes. Their identification could have an impact on patients’ care, justifying different therapy responses and inconsistencies identified across different randomized control trial results.

1. Introduction

The SARS-CoV-2 infection causes a wide spectrum of clinical manifestations, ranging from the complete absence of symptoms to a severe acute respiratory syndrome with a high fatality rate [1].
The natural history and pathophysiology of COVID-19 involves complex mechanisms, including viral infection and replication, immunological mechanisms, and a broad dysregulation of clotting, host defense, and endothelial dysfunction [2,3]. The evidence collected shows that many clinical and biochemical parameters are altered in COVID-19 patients, and this seems to be correlated with the severity of the disease and its prognosis [4].
On that regard, the characterization of clinical phenotypes within COVID-19 infected patients may represent a paramount clinical step, allowing us to increase our knowledge about the mechanisms of this disease and to comprehend inherent changes in patients’ diagnoses, prognoses, and treatment responses [5]. Azoulay et al. and Fish et al. have already provided insights into the biological plausibility of these phenotypes and their unique associations between several different biomarkers and patients’ prognoses [6,7]. Several clinical studies have further extended the characterization of these COVID-19 clinical phenotypes. Bruse et al. aimed to apply previously established clinical sepsis phenotypes to patients with COVID-19. The authors identified different phenotype distribution and biomarkers profiles when comparing to those previously established, confirming the need to conclusively establish clear clinical phenotypes to aid prognostication and prediction of treatment efficacy [8]. Siepel et al. attempted to confirm these clinical phenotypes using an unsupervised machine learning analysis. The authors conclude that the dynamic interchange of patients through phenotypes over time could explain the heterogeneity of results when aiming to describe these COVID-19 clinical phenotypes [9]. This evidence was in agreement with more recent studies that found similar clinical phenotypes across different countries [10,11,12,13,14,15,16].
The role of different biomarkers profiles in defining these COVID-19 phenotypes in critically ill patients has yet to be confirmed, but it could potentially aid to characterize these clinical phenotypes and to predict patients’ outcomes. In fact, biomarkers have been shown to be able to provide an assessment of several organ dysfunctions in COVID-19 patients, such as the risk of acute respiratory distress syndrome and disseminated intravascular coagulation [17,18]. Patients who present with lymphopenia, elevated D-dimer and troponin, hyperferritinemia, and increased lactic dehydrogenase usually have an unfavorable course of the disease [19]. Other non-specific biomarkers of COVID-19, such as the presence of neutrophilia, thrombocytopenia, hypoalbuminemia, elevation of liver enzymes and creatinine, as well as inflammatory markers, such as C-reactive protein (CRP) and interleukin 6 (IL-6), were also associated with a worse prognosis [20,21,22,23,24,25,26].
Blair et al. tried to stratify COVID-19 patients in different clusters using distinct blood inflammatory biomarkers. Their results pointed to three distinct inflammatory biomarker patterns that could stratify a heterogeneous population and that were associated with comorbid diseases and illness severity of those patients. This exploratory step gave important insights regarding the inclusion of biomarkers in the definition of clinical phenotypes, suggesting that these could be used for a personalized approach to the triage of care and therapeutic indications [27].
It has become progressively clearer that this biomarker phenotyping could lead the way to a more precise and personalized care for COVID-19 patients, highlighting the possibility of future treatments and interventions taking into account the disease’s particularities and patients’ predisposition and response [5,28,29]. Notably, the refinement of trials through prognostic enrichment is expected to increase the likelihood for beneficial effects of an intervention to emerge. Following these developments, this precision medicine paradigm may enable the tailoring of treatment to apparently similar patients with vastly different biochemical profiles and outcomes [30,31,32].
The main aim of our study is to assess and identify COVID-19 patients’ phenotypes using clinical and biomarker profiles, allowing the establishment of plausible correlations with clinical outcomes.

2. Materials and Methods

2.1. Design and Setting

This was a prospective observational cohort study performed in five ICUs of Brazil and Portugal. Briefly, during the study period (March 2020 to June 2021), we evaluated every adult patient (≥18 year) who required ICU admission with PCR confirmed SARS-CoV-2 infection.

2.2. Definitions, Selection of Participants and Data Collection

The Ethics Committees of the Hospital Copa Star and Hospital Copa D’or in Rio de Janeiro, Brazil, and of the Portuguese Ethics Committee for Clinical Investigation in Lisbon, Portugal, approved the study (CAAE: 17079119.7.0000.5249; and REC: 2020_EO_02, respectively) and we received the waiver of the written informed consent.
Sample size was calculated considering COVID-19’s population prevalence from Portugal and Brazil at the time of the study’s proposal, using the highest limit of the 95% confidence level and considering a 5% marginal error. We collected demographic, clinical, and laboratory data using standardized case report forms and patients’ daily chart revisions, and included the Simplified Acute Physiology Score (SAPS) III [33] and the Sequential Organ Failure Assessment (SOFA) at admission [34].
Need of mechanical ventilation (MV), intensive care unit (ICU), hospital length of stay (LOS), and mortality rates from any cause were also assessed. All patients were followed up with until death or hospital discharge. One investigator conducted a revision of the data collected through a random check of 10% of all datasets in order to confirm internal validity and data consistency. After data collection, patients’ individual data with missing values per variable rate above 10% were excluded from the initial analysis. The remaining missing values were handled under the assumption of missing at random and no imputation methods were applied (missing values percentage negligibly small (1.257%)) [35].
CRP, procalcitonin (PCT), IL-6 and D-dimer levels, and leucocyte and lymphocyte counts were measured at ICU admission, considering their value as inflammatory biomarkers in COVID-19 patients [20,21,22,23,24,25,26]. Lymphocyte serum counts, regardless of the leucocyte serum counts analysis, were considered, due to the previously collected evidence of an independent association between disease severity and mortality in COVID-19-patients [17,18]. Biomarkers levels were determined with the Roche Cobas Integra 800 analyzer (Roche Diagnostic, Indianapolis, IN, USA).
The main outcome of interest was to correlate clinical and biomarker signatures at ICU admission to derive COVID-19 phenotypes and determine their correlation with clinical outcomes. Other measured outcomes were organ dysfunctions assessed by SOFA, ICU and hospital LOS, and mortality.

2.3. Data Processing and Statistical Analysis

Baseline clinical and serum data collected were considered for the Principal Component Analysis and clustering after data processing as previously described, confirming its completeness, consistency, and missing values.
Variables selected to be included for modeling were age, SAPS III score, SOFA score, CRP, D-dimer, PCT, IL-6 levels, and leucocyte and lymphocyte counts at ICU admission after standardization, in order to reduce magnitude effect.
We used agglomerative clustering on the selected variables using Ward’s method and applying squared Euclidean Distance as the similarity measure [36]. The classification was conducted without consideration of clinical outcomes. The aim of this statistical method was to find relatively homogeneous clusters of cases based on measured characteristics. The log-likelihood method was used to determine inter-subject distance and specific classification of participants. The model was produced using the Schwarz-Bayesian criterion, aiming to identify the most goodness of fit of the clustering model. We considered grouping patients into 2–6 phenotypes and determined the optimal number of three as the model with higher goodness of fit. The number of clusters was confirmed using an elbow plot, where the inflection point was used as a cutoff.
We then performed a Principal Component Analysis on the average covariance matrix to visualize the relationships among the three phenotypes and assess variable contributions.
Once the clusters were determined, we compared the association between classes and clinical outcomes (ICU and hospital LOS, mortality, MV need, vasopressor support, and renal replacement therapy) using logistic regression models. We used standard descriptive statistics and reported continuous variables as median [25–75% interquartile range]. Comparisons between groups were performed with two-tailed unpaired Student t-test, one-way ANOVA, Mann-Whitney U, or Kruskal-Wallis H tests for continuous variables according to data distribution. The Fisher exact test and Chi-square test were used to carry out comparisons between categorical variables as appropriate. All Gaussian distributed variables were expressed as mean (SD) and non-normally distributed variables as median (interquartile range [IQR]). Categorical variables were expressed as numbers and percentages.
We carried out all statistical analyses using the SPSS 23.0 software package (Chicago, IL, USA) and Prism 6.0 (Graphpad, Boston, MA, USA).

3. Results

3.1. Cluster Analysis

In total, 814 patients were included. The main demographics and clinical characteristics are depicted in Table 1. The cluster analysis revealed a three-class model as the most useful in the analyzed population, as presented in Table 2. The principal component analysis also defined three principal population-describing factors as appreciated through an elbow plot (Figure S1, Tables S1 and S2 of Supplementary Materials), supporting that three-class model after plotting for those factors (Figure 1).
This three-class model allocated 407 patients in phenotype A, 244 patients in phenotype B, and 163 patients in phenotype C, allowing for the definition of three distinct COVID-19 phenotypes (Table 1).

3.2. Phenotype’s Characterization and Clinical Outcomes

The main demographic, clinical, and analytical characteristics among the three phenotypes are presented in Table 1.
A clear distinction can be perceived in phenotype A versus phenotypes B and C. Phenotype A patients are significantly older, with a higher severity index (either by SOFA or SAPS III) and a higher requirement of organ support (either ventilatory or cardiovascular) (Table 1). Furthermore, these patients also showed significantly longer ICU and hospital LOS, and a higher in-hospital mortality rate.
Phenotypes B and C, on the other hand, demonstrate some overlapping of clinical and demographic characteristics but with markedly different outcomes. Although no difference is found between these two phenotypes regarding age, SOFA score, and SAPS III at ICU admission, and ventilatory, vasopressor, and renal supports, patients with phenotype C have higher ICU and hospital LOS (10 (3;10) days vs. 6 (5;11) days and 13 (2;14) days vs. 9 (1;14) days, respectively, p < 0.001) but a significantly lower in-hospital mortality rate (1.2% vs. 7%, respectively, p = 0.007), when compared to patients with phenotype B.

3.3. Phenotype’s Biomarker Profile

Higher baseline inflammatory markers were associated with a higher mortality rate (p = 0.001). Phenotype A patients presented with higher CRP, PCT, and IL-6 serum levels at ICU admission (p < 0.001). Furthermore, this phenotype also presented higher white blood cell counts, although no difference was registered in the lymphocyte count between groups.
Phenotypes B and C presented clearly different inflammatory profiles and mortality rates. Phenotype C patients presented higher PCT and IL-6 serum levels with lower CRP serum levels and mortality rate (1.2%). On the other hand, phenotype B patients presented lower PCT and IL-6 serum levels and higher CRP serum levels and a 4-fold higher mortality rate (7%).

4. Discussion

In our large cohort of severe COVID-19 patients, we were able to identify three distinct COVID-19 phenotypes. These phenotypes show significant differences in their clinical features, inflammatory profiles, and hospital outcomes. Our findings urge the necessity of comprehending COVID-19 as a disease that encompasses several patients with heterogenous outcomes. This may potentially prompt more patient-targeting therapeutic protocols and promote a medical practice based in more individualized and targeted approach to these patients.
Firstly, our cluster analysis identified a high-mortality rate phenotype (phenotype A) with an in-hospital mortality rate of 12%. This phenotype is clearly composed of older patients with high SOFA and SAPS III at ICU admission, a higher pro-inflammatory profile, and poor outcomes. This data reproduces previous evidence collected, describing a strong correlation between mortality rate and age, high CRP and PCT levels, and severity indexes at ICU admission in COVID-19 patients [37,38,39,40,41,42,43]. Clearly, this phenotype convenes the most severe COVID-19 patients, with significant differences to phenotypes B and C.
Phenotypes B and C, on the other hand, represent subpopulations of younger patients with lower severity scores and fewer organ support requirements. Albeit contrastingly different when comparing to phenotype A, these phenotypes B and C lose their apparent overlapping characteristics when their inflammatory biomarker profiles and mortality rates are considered. In fact, patients with phenotype C exhibit higher PCT and IL-6 serum levels with lower CRP serum levels and lower in-hospital mortality rates.
These results of COVID-19 patients’ phenotypes also demonstrate biological plausibility and are strikingly similar to those previously described by Azoulay et al. [6]. The authors reported the clinical and laboratory features of a cohort of 85 consecutive COVID-19 patients admitted to ICU. In this study, the hierarchical clustering also identified three clinical phenotypes of COVID-19 patients with a significant overlap with those identified in our study. Based on a bigger cohort of COVID-19 patients, our hierarchical clustering documents the same three phenotypes, with similar clinical profiles, inflammatory biomarker patterns, and hospital outcomes.
Our results are also consistently overlapping and complementary to those found in other COVID-19 phenotype-defining studies [10,11,12,13,14,15,16,27]. In a study by Lusczek et al., three clinical phenotypes were identified with similar clinical characteristics to those described in our study. A high-mortality rate phenotype composed of older patients with a higher number of comorbidities and higher in-hospital complications and CRP serum levels at admission was also described. Correspondingly, it also described two less-ill phenotypes of patients with lower mortality rates and lower inflammatory profiles [13]. The same results were also found in studies by Rodríguez et al. and Siepel et al., which deployed unsupervised clustering analysis of critically ill patients with COVID-19 in Spain and in the Netherlands, respectively [9,14]. The derived clinical phenotypes exhibit a marked similarity to our findings and significantly correlated with host-response patterns and ICU mortality. Across these studies, three phenotypes were identified with different high severity of illness and ICU mortality rates, and prominent dissimilarities between clinical phenotypes were associated to different inflammation profiles.
Although no statistical difference was found allowing for the definition of an inflammatory biomarker profile to segregate patient phenotypes B and C, this data is strikingly different from evidence previously collected [41,42]. Most of those results, suggesting a linear association between PCT, CRP, and IL-6 serum levels and patients’ mortality rate, stemmed from retrospective and observational analysis. Our results do not support a direct association between IL-6 and PCT serum levels and mortality rate in severe COVID-19 patients but suggest that CRP may be a more reliable surrogate of COVID-19 patient’s outcomes [44,45,46,47].
These findings could be somehow related to the heterogeneity of RCT results, namely of IL-6 receptor inhibitors [48]. Our and Azoulay’s results raised the possibility of a lack of effect or even harm of therapeutic options targeting this individual biomarker [7]. Pinpointing these isolated biomarkers as therapeutical targets may be insufficient to address all the inflammatory and mechanisms active in these patients. Furthermore, the prognostic value of this biomarker may be underlined in a more complex inflammatory profile, dependent on patient’s immunological variables, and subject to a different prevalence of these clinical phenotypes [49,50,51,52].
Our study conveys several strengths. It is composed of a large cohort of COVID-19 patients for a multinational and multicenter cluster analysis, circumventing some possible sources of bias. However, it is not without limitations. It is an observational study and, although it represents a large sample size of COVID-19 patients, it fails to document a clear segregation between the inflammatory profiles of the three clinical phenotypes that could reliably identify a subset of patients with higher mortality rates. A higher number of patients could also improve phenotype homogeneity. Furthermore, an external validation of the clusters is absent, undermining the potential of finding these phenotypes in other studied COVID-19 populations. Moreover, it convenes a study period of two variants and significant therapeutic recommendation changes, specifically, but not limited to, in the use of invasive and non-invasive mechanical ventilation, that could count for differences within the studied population. Furthermore, we recognize a significant absence of data concerning time from onset of symptoms and time until intubation and mechanical ventilation. Finally, COVID-19 treatment effects were not evaluated, and long-term follow-up of patients was not performed. Therefore, the impact of these phenotypes on their long-term mortality rate was not evaluated.

5. Conclusions

Severe COVID-19 patients exhibit three different phenotypes with distinct clinical and biochemical profiles and outcomes. Our data presents a clearly high-mortality group with a high level of organ requirements that should be promptly recognized in order to improve survivability and early organ support. It also suggests an association of CRP serum levels with patients’ outcomes. In contrast to our expectations, there was no evidence of any association between IL-6 serum levels and patients’ outcomes, namely in phenotype C, which presented intermediate levels of Il-6 serum levels and low patient mortality. Further validation of these phenotypes with external cohorts could allow the implementation of point-of-care biomarker determinations to guide patients’ management.
Our data also conveys that different patients’ profiles may translate to an early identification of patients with variable responses to certain COVID-19 directed therapies and may help understand heterogeneity of treatment effects. It highlights different prognoses in similar-appearing patients with COVID-19, reinforcing the potential use of biochemical phenotyping in the clinical characterization of these patients. Furthermore, it may provide impactful insights in future randomized clinical trial analyses and a concept for trial enrichment.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/jcm12083035/s1, Figure S1: Principal Component analysis using elbow method of the cluster dataset; Table S1: Principal Components analysis (PCA) of the cluster dataset; Table S2: Demographic, clinical and outcome variables of patients between Phenotypes B and C.

Author Contributions

Conceptualization: J.P.C., V.C.d.S.D. and P.P. Data collection: J.P.C., A.d.F.T., R.C.C.C.d.M., R.M., H.C., G.E., N.O. and T.C. Data analysis: J.P.C. and V.C.d.S.D.; Writing—original draft: J.P.C., V.C.d.S.D., A.d.F.T., R.C.C.C.d.M., R.M., H.C., G.E., N.O., T.C. and P.P. Writing (review and editing): J.P.C., V.C.d.S.D. and P.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The Ethics Committees of the Hospital Copa Star and Hospital Copa D’or in Rio de Janeiro, and of the Portuguese Ethics Committee for Clinical Investigation in Lisbon approved the study (CAAE: 17079119.7.0000.5249; and REC: 2020_EO_02, respectively) and we got the waiver of the written informed consent.

Informed Consent Statement

Patient consent was waived due to COVID-19 pandemic state, as stated Ethics Committees of the Hospital Copa Star and Hospital Copa D’or in Rio de Janeiro, and of the Portuguese Ethics Committee for Clinical Investigation in Lisbon.

Data Availability Statement

The datasets generated and/or analyzed during the current study are not publicly available due to privacy issues but are available from the corresponding author on reasonable request.

Conflicts of Interest

P.P. received fees for lecture from Gilead, Pfizer and Mundipharma, consulting from MSD and Sanofi and unrestricted research grant from Abionic. J.P.C., V.C.d.S.D., A.d.F.T., R.C.C.C.d.M., R.M., H.C., G.E., N.O. and T.C. declare no conflict of interest.

References

  1. Guan, W.J.; Ni, Z.Y.; Hu, Y.; Liang, W.H.; Ou, C.Q.; He, J.X.; Liu, L.; Shan, H.; Lei, C.L.; Hui, D.S.C.; et al. Clinical Characteristics of Coronavirus Disease 2019 in China. N. Engl. J. Med. 2020, 382, 1708–1720. [Google Scholar] [CrossRef] [PubMed]
  2. Lucas, C.; Wong, P.; Klein, J.; Castro, T.B.R.; Silva, J.; Sundaram, M.; Ellingson, M.K.; Mao, T.; Oh, J.E.; Israelow, B.; et al. Longitudinal analyses reveal immunological misfiring in severe COVID-19. Nature 2020, 584, 463–469. [Google Scholar] [CrossRef] [PubMed]
  3. Overmyer, K.A.; Shishkova, E.; Miller, I.J.; Balnis, J.; Bernstein, M.N.; Peters-Clarke, T.M.; Meyer, J.G.; Quan, Q.; Muehlbauer, L.K.; Trujillo, E.A.; et al. Large-Scale Multi-omic Analysis of COVID-19 Severity. Cell Syst. 2021, 12, 23–40.e27. [Google Scholar] [CrossRef] [PubMed]
  4. Del Valle, D.M.; Kim-Schulze, S.; Huang, H.H.; Beckmann, N.D.; Nirenberg, S.; Wang, B.; Lavin, Y.; Swartz, T.H.; Madduri, D.; Stock, A.; et al. An inflammatory cytokine signature predicts COVID-19 severity and survival. Nat. Med. 2020, 26, 1636–1643. [Google Scholar] [CrossRef]
  5. Arish, M.; Naz, F. Personalized therapy: Can it tame the COVID-19 monster? Per. Med. 2021, 18, 583–593. [Google Scholar] [CrossRef]
  6. Azoulay, E.; Zafrani, L.; Mirouse, A.; Lengliné, E.; Darmon, M.; Chevret, S. Clinical phenotypes of critically ill COVID-19 patients. Intensive Care Med. 2020, 46, 1651–1652. [Google Scholar] [CrossRef]
  7. Fish, M.; Rynne, J.; Jennings, A.; Lam, C.; Lamikanra, A.A.; Ratcliff, J.; Cellone-Trevelin, S.; Timms, E.; Jiriha, J.; Tosi, I.; et al. Coronavirus disease 2019 subphenotypes and differential treatment response to convalescent plasma in critically ill adults: Secondary analyses of a randomized clinical trial. Intensive Care Med. 2022, 48, 1525–1538. [Google Scholar] [CrossRef]
  8. Bruse, N.; Kooistra, E.J.; Jansen, A.; van Amstel, R.B.E.; de Keizer, N.F.; Kennedy, J.N.; Seymour, C.; van Vught, L.A.; Pickkers, P.; Kox, M. Clinical sepsis phenotypes in critically ill COVID-19 patients. Crit. Care 2022, 26, 244. [Google Scholar] [CrossRef]
  9. Siepel, S.; Dam, T.A.; Fleuren, L.M.; Girbes, A.R.J.; Hoogendoorn, M.; Thoral, P.J.; Elbers, P.W.G.; Bennis, F.C. Evolution of Clinical Phenotypes of COVID-19 Patients During Intensive Care Treatment: An Unsupervised Machine Learning Analysis. J. Intensive Care Med. 2023. [Google Scholar] [CrossRef]
  10. Portacci, A.; Carpagnano, G.E.; Tummolo, M.G.; Santomasi, C.; Palma, L.; Fasano, D.; Resta, E.; Lozupone, M.; Solfrizzi, V.; Panza, F.; et al. COVID-19 clinical phenotypes and short-term outcomes: Differences between the first and the second wave of pandemic in Italy. Expert Rev. Respir. Med. 2021, 15, 1619–1625. [Google Scholar] [CrossRef]
  11. Da Silva, J.F.; Hernandez-Romieu, A.C.; Browning, S.D.; Bruce, B.B.; Natarajan, P.; Morris, S.B.; Gold, J.A.W.; Neblett Fanfair, R.; Rogers-Brown, J.; Rossow, J.; et al. COVID-19 Clinical Phenotypes: Presentation and Temporal Progression of Disease in a Cohort of Hospitalized Adults in Georgia, United States. Open Forum Infect. Dis. 2020, 8. [Google Scholar] [CrossRef] [PubMed]
  12. Gattinoni, L.; Chiumello, D.; Caironi, P.; Busana, M.; Romitti, F.; Brazzi, L.; Camporota, L. COVID-19 pneumonia: Different respiratory treatments for different phenotypes? Intensive Care Med. 2020, 46, 1099–1102. [Google Scholar] [CrossRef] [PubMed]
  13. Lusczek, E.R.; Ingraham, N.E.; Karam, B.S.; Proper, J.; Siegel, L.; Helgeson, E.S.; Lotfi-Emran, S.; Zolfaghari, E.J.; Jones, E.; Usher, M.G.; et al. Characterizing COVID-19 clinical phenotypes and associated comorbidities and complication profiles. PLoS ONE 2021, 16, e0248956. [Google Scholar] [CrossRef] [PubMed]
  14. Rodríguez, A.; Ruiz-Botella, M.; Martín-Loeches, I.; Jimenez Herrera, M.; Solé-Violan, J.; Gómez, J.; Bodí, M.; Trefler, S.; Papiol, E.; Díaz, E.; et al. Deploying unsupervised clustering analysis to derive clinical phenotypes and risk factors associated with mortality risk in 2022 critically ill patients with COVID-19 in Spain. Crit. Care 2021, 25, 63. [Google Scholar] [CrossRef] [PubMed]
  15. Su, L.; Zhang, Z.; Zheng, F.; Pan, P.; Hong, N.; Liu, C.; He, J.; Zhu, W.; Long, Y.; Liu, D. Five novel clinical phenotypes for critically ill patients with mechanical ventilation in intensive care units: A retrospective and multi database study. Respir. Res. 2020, 21, 325. [Google Scholar] [CrossRef]
  16. de Prost, N.; Audureau, E.; Heming, N.; Gault, E.; Pham, T.; Chaghouri, A.; de Montmollin, N.; Voiriot, G.; Morand-Joubert, L.; Joseph, A.; et al. Clinical phenotypes and outcomes associated with SARS-CoV-2 variant Omicron in critically ill French patients with COVID-19. Nat. Commun. 2022, 13, 6025. [Google Scholar] [CrossRef]
  17. Lippi, G.; Plebani, M. The critical role of laboratory medicine during coronavirus disease 2019 (COVID-19) and other viral outbreaks. Clin. Chem. Lab. Med. 2020, 58, 1063–1069. [Google Scholar] [CrossRef]
  18. Van der Zee, P.; Rietdijk, W.; Somhorst, P.; Endeman, H.; Gommers, D. A systematic review of biomarkers multivariately associated with acute respiratory distress syndrome development and mortality. Crit. Care 2020, 24, 243. [Google Scholar] [CrossRef]
  19. Letelier, P.; Encina, N.; Morales, P.; Riffo, A.; Silva, H.; Riquelme, I.; Guzmán, N. Role of biochemical markers in the monitoring of COVID-19 patients. J. Med. Biochem. 2021, 40, 115–128. [Google Scholar] [CrossRef]
  20. Leisman, D.E.; Ronner, L.; Pinotti, R.; Taylor, M.D.; Sinha, P.; Calfee, C.S.; Hirayama, A.V.; Mastroiani, F.; Turtle, C.J.; Harhay, M.O.; et al. Cytokine elevation in severe and critical COVID-19: A rapid systematic review, meta-analysis, and comparison with other inflammatory syndromes. Lancet Respir. Med. 2020, 8, 1233–1244. [Google Scholar] [CrossRef]
  21. Fajgenbaum, D.C.; June, C.H. Cytokine Storm. N. Engl. J. Med. 2020, 383, 2255–2273. [Google Scholar] [CrossRef] [PubMed]
  22. McElvaney, O.J.; McEvoy, N.L.; McElvaney, O.F.; Carroll, T.P.; Murphy, M.P.; Dunlea, D.M.; Ní Choileáin, O.; Clarke, J.; O’Connor, E.; Hogan, G.; et al. Characterization of the Inflammatory Response to Severe COVID-19 Illness. Am. J. Respir. Crit. Care Med. 2020, 202, 812–821. [Google Scholar] [CrossRef] [PubMed]
  23. Su, Y.; Chen, D.; Yuan, D.; Lausted, C.; Choi, J.; Dai, C.L.; Voillet, V.; Duvvuri, V.R.; Scherler, K.; Troisch, P.; et al. Multi-Omics Resolves a Sharp Disease-State Shift between Mild and Moderate COVID-19. Cell 2020, 183, 1479–1495.e1420. [Google Scholar] [CrossRef]
  24. Wang, J.; Jiang, M.; Chen, X.; Montaner, L.J. Cytokine storm and leukocyte changes in mild versus severe SARS-CoV-2 infection: Review of 3939 COVID-19 patients in China and emerging pathogenesis and therapy concepts. J. Leukoc. Biol. 2020, 108, 17–41. [Google Scholar] [CrossRef] [PubMed]
  25. Li, L.; Li, J.; Gao, M.; Fan, H.; Wang, Y.; Xu, X.; Chen, C.; Liu, J.; Kim, J.; Aliyari, R.; et al. Interleukin-8 as a Biomarker for Disease Prognosis of Coronavirus Disease-2019 Patients. Front. Immunol. 2020, 11, 602395. [Google Scholar] [CrossRef] [PubMed]
  26. Ma, A.; Zhang, L.; Ye, X.; Chen, J.; Yu, J.; Zhuang, L.; Weng, C.; Petersen, F.; Wang, Z.; Yu, X. High Levels of Circulating IL-8 and Soluble IL-2R Are Associated With Prolonged Illness in Patients With Severe COVID-19. Front. Immunol. 2021, 12, 626235. [Google Scholar] [CrossRef]
  27. Blair, P.W.; Brandsma, J.; Chenoweth, J.; Richard, S.A.; Epsi, N.J.; Mehta, R.; Striegel, D.; Clemens, E.G.; Alharthi, S.; Lindholm, D.A.; et al. Distinct blood inflammatory biomarker clusters stratify host phenotypes during the middle phase of COVID-19. Sci. Rep. 2022, 12, 22471. [Google Scholar] [CrossRef]
  28. Bakkerus, L.; Pickkers, P. Personalized medicine in COVID-19. Intensive Care Med. 2022, 48, 1607–1610. [Google Scholar] [CrossRef]
  29. Dopazo, J.; Maya-Miles, D.; García, F.; Lorusso, N.; Calleja, M.; Pareja, M.J.; López-Miranda, J.; Rodríguez-Baño, J.; Padillo, J.; Túnez, I.; et al. Implementing Personalized Medicine in COVID-19 in Andalusia: An Opportunity to Transform the Healthcare System. J. Pers. Med. 2021, 11, 475. [Google Scholar] [CrossRef]
  30. Teodori, L.; Osimani, B.; Isidoro, C.; Ramakrishna, S. Mass versus personalized medicine against COVID-19 in the “system sciences” era. Cytom. A 2022, 101, 995–999. [Google Scholar] [CrossRef]
  31. Zhou, A.; Sabatello, M.; Eyal, G.; Lee, S.S.-J.; Rowe, J.W.; Stiles, D.F.; Swanson, A.; Appelbaum, P.S. Is precision medicine relevant in the age of COVID-19? Genet. Med. 2021, 23, 999–1000. [Google Scholar] [CrossRef] [PubMed]
  32. DeMerle, K.; Angus, D.C.; Seymour, C.W. Precision Medicine for COVID-19: Phenotype Anarchy or Promise Realized? JAMA 2021, 325, 2041–2042. [Google Scholar] [CrossRef] [PubMed]
  33. Le Gall, J.R.; Lemeshow, S.; Saulnier, F. A new Simplified Acute Physiology Score (SAPS II) based on a European/North American multicenter study. JAMA 1993, 270, 2957–2963. [Google Scholar] [CrossRef] [PubMed]
  34. Vincent, J.L.; Moreno, R.; Takala, J.; Willatts, S.; De Mendonça, A.; Bruining, H.; Reinhart, C.K.; Suter, P.M.; Thijs, L.G. The SOFA (Sepsis-related Organ Failure Assessment) score to describe organ dysfunction/failure. On behalf of the Working Group on Sepsis-Related Problems of the European Society of Intensive Care Medicine. Intensive Care Med. 1996, 22, 707–710. [Google Scholar] [CrossRef]
  35. O’Kelly, M.; Ratitch, B. Clinical Trials with Missing Data: A Guide for Practitioners; John Wiley & Sons: Hoboken, NJ, USA, 2014. [Google Scholar]
  36. Hagen, D.R.; Tidor, B. Efficient Bayesian estimates for discrimination among topologically different systems biology models. Mol. Biosyst. 2015, 11, 574–584. [Google Scholar] [CrossRef]
  37. Henkens, M.; Raafs, A.G.; Verdonschot, J.A.J.; Linschoten, M.; van Smeden, M.; Wang, P.; van der Hooft, B.H.M.; Tieleman, R.; Janssen, M.L.F.; Ter Bekke, R.M.A.; et al. Age is the main determinant of COVID-19 related in-hospital mortality with minimal impact of pre-existing comorbidities, a retrospective cohort study. BMC Geriatr. 2022, 22, 184. [Google Scholar] [CrossRef]
  38. Bonanad, C.; García-Blas, S.; Tarazona-Santabalbina, F.; Sanchis, J.; Bertomeu-González, V.; Fácila, L.; Ariza, A.; Núñez, J.; Cordero, A. The Effect of Age on Mortality in Patients With COVID-19: A Meta-Analysis With 611,583 Subjects. J. Am. Med. Dir. Assoc. 2020, 21, 915–918. [Google Scholar] [CrossRef]
  39. Lentner, J.; Adams, T.; Knutson, V.; Zeien, S.; Abbas, H.; Moosavi, R.; Manuel, C.; Wallace, T.; Harmon, A.; Waters, R.; et al. C-reactive protein levels associated with COVID-19 outcomes in the United States. J. Osteopath. Med. 2021, 121, 869–873. [Google Scholar] [CrossRef]
  40. Malik, P.; Patel, U.; Mehta, D.; Patel, N.; Kelkar, R.; Akrmah, M.; Gabrilove, J.L.; Sacks, H. Biomarkers and outcomes of COVID-19 hospitalisations: Systematic review and meta-analysis. BMJ Evid. Based Med. 2021, 26, 107–108. [Google Scholar] [CrossRef]
  41. Loomba, R.S.; Villarreal, E.G.; Farias, J.S.; Aggarwal, G.; Aggarwal, S.; Flores, S. Serum biomarkers for prediction of mortality in patients with COVID-19. Ann. Clin. Biochem. 2022, 59, 15–22. [Google Scholar] [CrossRef] [PubMed]
  42. Fayed, M.; Patel, N.; Angappan, S.; Nowak, K.; Vasconcelos Torres, F.; Penning, D.H.; Chhina, A.K. Sequential Organ Failure Assessment (SOFA) Score and Mortality Prediction in Patients With Severe Respiratory Distress Secondary to COVID-19. Cureus 2022, 14, e26911. [Google Scholar] [CrossRef] [PubMed]
  43. Beigmohammadi, M.T.; Amoozadeh, L.; Rezaei Motlagh, F.; Rahimi, M.; Maghsoudloo, M.; Jafarnejad, B.; Eslami, B.; Salehi, M.R.; Zendehdel, K. Mortality Predictive Value of APACHE II and SOFA Scores in COVID-19 Patients in the Intensive Care Unit. Can. Respir. J. 2022, 2022, 5129314. [Google Scholar] [CrossRef] [PubMed]
  44. Tong-Minh, K.; van der Does, Y.; Engelen, S.; de Jong, E.; Ramakers, C.; Gommers, D.; van Gorp, E.; Endeman, H. High procalcitonin levels associated with increased intensive care unit admission and mortality in patients with a COVID-19 infection in the emergency department. BMC Infect. Dis. 2022, 22, 165. [Google Scholar] [CrossRef] [PubMed]
  45. Twe, C.W.; Khoo, D.K.Y.; Law, K.B.; Ahmad Nordin, N.S.B.; Sathasivan, S.; Lim, K.C.; Atikah, S.K.; Syed Badaruddin, S.; Chidambaram, S.K. The role of procalcitonin in predicting risk of mechanical ventilation and mortality among moderate to severe COVID-19 patients. BMC Infect. Dis. 2022, 22, 378. [Google Scholar] [CrossRef]
  46. Smilowitz, N.R.; Kunichoff, D.; Garshick, M.; Shah, B.; Pillinger, M.; Hochman, J.S.; Berger, J.S. C-reactive protein and clinical outcomes in patients with COVID-19. Eur. Heart J. 2021, 42, 2270–2279. [Google Scholar] [CrossRef] [PubMed]
  47. Ali, N. Elevated level of C-reactive protein may be an early marker to predict risk for severity of COVID-19. J. Med. Virol. 2020, 92, 2409–2411. [Google Scholar] [CrossRef]
  48. Gomez, C.A.; Kalil, A.C. Immunomodulators in Mechanically Ventilated Patients With COVID-19: Lessons Learned From Underpowered Trials. Crit. Care Med. 2022, 50, 1821–1828. [Google Scholar] [CrossRef]
  49. Avni, T.; Leibovici, L.; Cohen, I.; Atamna, A.; Guz, D.; Paul, M.; Gafter-Gvili, A.; Yahav, D. Tocilizumab in the treatment of COVID-19-a meta-analysis. QJM 2021, 114, 577–586. [Google Scholar] [CrossRef]
  50. Salama, C.; Han, J.; Yau, L.; Reiss, W.G.; Kramer, B.; Neidhart, J.D.; Criner, G.J.; Kaplan-Lewis, E.; Baden, R.; Pandit, L.; et al. Tocilizumab in Patients Hospitalized with COVID-19 Pneumonia. N. Engl. J. Med. 2021, 384, 20–30. [Google Scholar] [CrossRef]
  51. Boppana, T.K.; Mittal, S.; Madan, K.; Mohan, A.; Hadda, V.; Guleria, R. Tocilizumab for COVID-19: A systematic review and meta-analysis of randomized controlled trials. Monaldi Arch. Chest Dis. 2022, 92. [Google Scholar] [CrossRef]
  52. Mahroum, N.; Watad, A.; Bridgewood, C.; Mansour, M.; Nasr, A.; Hussein, A.; Khamisy-Farah, R.; Farah, R.; Gendelman, O.; Lidar, M.; et al. Systematic Review and Meta-Analysis of Tocilizumab Therapy versus Standard of Care in over 15,000 COVID-19 Pneumonia Patients during the First Eight Months of the Pandemic. Int. J. Env. Res. Public Health 2021, 18, 9149. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Factor map of the Principal Component Analysis showing the distribution of each patient in each cluster.
Figure 1. Factor map of the Principal Component Analysis showing the distribution of each patient in each cluster.
Jcm 12 03035 g001
Table 1. Demographic, clinical and outcome variables of patients according to the COVID-19 phenotype *.
Table 1. Demographic, clinical and outcome variables of patients according to the COVID-19 phenotype *.
PHENOTYPE A (n = 407; 60%)PHENOTYPE B
(n = 244; 24%)
PHENOTYPE C (n = 163; 16%)p ¥
Age, years, (median (IQR))81 (65–97)62 (47–79)63 (47–80)<0.001
Gender, males, (n, %)221 (54.3%)131 (53.7%)78 (47.9%)0.08
Previous Medical Comorbidities
Chronic Obstructive Pulmonary Disease (n, %)37 (9.1%)18 (7.4%)14 (8.9%)0.528
Asthma (n, %)18 (4.4%)8 (3.3%)5 (3.1%)0.347
Chronic Kidney Disease (n, %)89 (21.9%)29 (11.9%)23 (14.1%)0.125
Obesity (n, %)64 (15.7%)36 (14.8%)27 (16.7%)0.214
Diabetes Mellitus (n, %)129 (31.7%)82 (33.6%)47 (28.8%)0.08
Ischemic Cardiopathy (n, %)147 (36.1%)64 (26.2%)39 (23.9%)0.04
SOFA at admission (median (IQR))10 (5; 13)3 (2; 5)1 (0; 3)<0.001
SAPS III at admission (mean ± SD)78 ± 1050 ± 747 ± 12<0.001
Mechanical Ventilation (n, %)173 (42.5%)11 (4.6%)11 (6.9%)<0.001
Vasopressor Support (n, %)112 (27.5%)38 (15.6%)38 (14.4%)<0.001
Renal replacement therapy (n, %)58 (14.2%)19 (7.8%)19 (11.6%)0.152
Laboratory results
C reactive Protein at admission, mg/dL (median (IQR))32.3 (24.8; 81.5)20.0 (10.3; 40.6)17.20 (4.0; 23.7)0.013
Max registered C-Reactive protein, mg/dL (mean ± SD)32.3 ± 11.025.3 ± 10.418.6 ± 12.5<0.001
Procalcitonin at admission, ng/mL (median (IQR))3.30 (0.55; 3.35)0.17 (0.05; 0.23)0.22 (0.12; 0.23)<0.001
Max registered Procalcitonin, ng/mL (median (IQR))9.73 (0.86; 13.54)0.34 (0.06; 0.74)1.30 (0.70; 1.40)<0.001
D-dimer level at admission, ng/mL (median (IQR))1165 (587; 1663)610 (97; 753)202 (119; 262)0.003
Max D-dimer registered, ng/mL (median (IQR)2778 (875; 3822)655 (48; 1305)303 (78; 307)0.018
Minimum Leucocyte count registered, ×109 (mean ± SD)11.0 ± 7.095.0 ± 2.035.3 ± 2.2<0.001
Minimum Lymphocyte count registered, ×109 (median (IQR))0.52 (0.32; 0.62)0.36 (0.12; 0.39)0.57 (0.68; 1.09)0.146
IL-6 serum levels, mg/mL (median (IQR))57.9 (5.7; 61.0)35.4 (6.6; 42.7)41.0 (16.0; 49.0)0.01
Remdesivir, (n, %)359 (88.2%)243 (99.6%)147 (90.2%)0.167
Corticosteroid therapy
(n, %)
176 (43.2%)137 (97.9%)147 (90.0%)0.001
Ventilator-free days, days, (median (IQR))23 (20; 24)12 (4; 20)10 (4; 17)0.001
ICU length of stay, days, (median (IQR))14 (11; 15)6 (5; 11)10 (3; 10)<0.001
Hospital Length of stay, days, (mean ± sd)18 (7; 20)9 (1; 14)13 (2; 14)0.001
* IQR denotes Interquartile range and SD denotes standard deviation. ¥ p-values were determined using Chi-square test, ANOVA test and Kruskal-Wallis H test.
Table 2. Cluster class fitting analysis.
Table 2. Cluster class fitting analysis.
Number of Individuals per Class/Phenotype
Number of ClassesBIC ¥Entropy *N1N2N3N4p-Value **
26034.40.7778727 0.227
32094.10.87244163407 0.007
48056.70.52114475760.183
¥ BIC stands for Bayesian information criterion. * Entropy is an index of how well the classes are separated. It ranges from zero to one and values around 0.8 and up are generally considered a sign of a useful model. ** By likelihood ratio test, testing whether the number of classes improved the model fitness compared to a model using one fewer class.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Cidade, J.P.; de Souza Dantas, V.C.; de Figueiredo Thompson, A.; de Miranda, R.C.C.C.; Mamfrim, R.; Caroli, H.; Escudini, G.; Oliveira, N.; Castro, T.; Póvoa, P. Identification of Distinct Clinical Phenotypes of Critically Ill COVID-19 Patients: Results from a Cohort Observational Study. J. Clin. Med. 2023, 12, 3035. https://doi.org/10.3390/jcm12083035

AMA Style

Cidade JP, de Souza Dantas VC, de Figueiredo Thompson A, de Miranda RCCC, Mamfrim R, Caroli H, Escudini G, Oliveira N, Castro T, Póvoa P. Identification of Distinct Clinical Phenotypes of Critically Ill COVID-19 Patients: Results from a Cohort Observational Study. Journal of Clinical Medicine. 2023; 12(8):3035. https://doi.org/10.3390/jcm12083035

Chicago/Turabian Style

Cidade, José Pedro, Vicente Cés de Souza Dantas, Alessandra de Figueiredo Thompson, Renata Carnevale Carneiro Chermont de Miranda, Rafaela Mamfrim, Henrique Caroli, Gabriela Escudini, Natalia Oliveira, Taiza Castro, and Pedro Póvoa. 2023. "Identification of Distinct Clinical Phenotypes of Critically Ill COVID-19 Patients: Results from a Cohort Observational Study" Journal of Clinical Medicine 12, no. 8: 3035. https://doi.org/10.3390/jcm12083035

APA Style

Cidade, J. P., de Souza Dantas, V. C., de Figueiredo Thompson, A., de Miranda, R. C. C. C., Mamfrim, R., Caroli, H., Escudini, G., Oliveira, N., Castro, T., & Póvoa, P. (2023). Identification of Distinct Clinical Phenotypes of Critically Ill COVID-19 Patients: Results from a Cohort Observational Study. Journal of Clinical Medicine, 12(8), 3035. https://doi.org/10.3390/jcm12083035

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop