- 1Division of Pulmonary, Allergy, and Critical Care Medicine, Department of Medicine, NewYork-Presbyterian Hospital/Columbia University Irving Medical Center, New York, NY, United States
- 2Program for Hospital and Intensive Care Informatics, Department of Neurology, Columbia University Irving Medical Center, New York, NY, United States
- 3Department of Neurology, NewYork-Presbyterian Hospital/Columbia University Irving Medical Center, New York, NY, United States
- 4Department of Critical Care Medicine, University of Pittsburgh, Pittsburgh, PA, United States
- 5Department of Surgery, University of Pittsburgh, Pittsburgh, PA, United States
- 6Department of Biomedical Engineering, California State University Long Beach, Long Beach, CA, United States
Background: Characterization of coronavirus disease 2019 (COVID-19) endotypes may help explain variable clinical presentations and response to treatments. While risk factors for COVID-19 have been described, COVID-19 endotypes have not been elucidated.
Objectives: We sought to identify and describe COVID-19 endotypes of hospitalized patients.
Methods: Consensus clustering (using the ensemble method) of patient age and laboratory values during admission identified endotypes. We analyzed data from 528 patients with COVID-19 who were admitted to telemetry capable beds at Columbia University Irving Medical Center and discharged between March 12 to July 15, 2020.
Results: Four unique endotypes were identified and described by laboratory values, demographics, outcomes, and treatments. Endotypes 1 and 2 were comprised of low numbers of intubated patients (1 and 6%) and exhibited low mortality (1 and 6%), whereas endotypes 3 and 4 included high numbers of intubated patients (72 and 85%) with elevated mortality (21 and 43%). Endotypes 2 and 4 had the most comorbidities. Endotype 1 patients had low levels of inflammatory markers (ferritin, IL-6, CRP, LDH), low infectious markers (WBC, procalcitonin), and low degree of coagulopathy (PTT, PT), while endotype 4 had higher levels of those markers.
Conclusions: Four unique endotypes of hospitalized patients with COVID-19 were identified, which segregated patients based on inflammatory markers, infectious markers, evidence of end-organ dysfunction, comorbidities, and outcomes. High comorbidities did not associate with poor outcome endotypes. Further work is needed to validate these endotypes in other cohorts and to study endotype differences to treatment responses.
Introduction
Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-COV-2), has demonstrated a wide variety of clinical courses, including asymptomatic carriers (1), mild disease (2), brief hospitalizations (2), prolonged ICU courses (3, 4), and COVID-19 “long-haulers” with prolonged symptoms (5). The spectrum of disease seems broader than the spectrum caused by other respiratory viruses, such as non-SARS-COV-2 coronaviruses. The international scientific community is currently endeavoring to understand the biological constructs that influence the course of disease after COVID-19 infection. Improved understanding of the biological underpinnings of different COVID-19 courses could improve diagnosis, triage, management, and prognosis for patients.
Several patient characteristics are associated with more severe COVID-19 disease or worse outcomes, including older age (3, 6), male sex (3, 6), obesity (7), diabetes mellitus (DM) (8), cardiovascular disease (3, 6), chronic obstructive pulmonary disease (COPD) (9), and chronic kidney disease (CKD) (10). Knowledge of baseline characteristics (including demographics and/or initial laboratory values) can predict hospitalization and mortality (11). There may be a subset of patients with a hyperinflammatory response who are at increased risk of mortality (12, 13).
Understanding endotypes of disease can shed light on biological underpinnings of disease and identify those who are most susceptible. Endotypes are subtypes of a clinical condition which possess distinct functional or pathobiological mechanisms (with an implicit variable likelihood of response to therapies across endotypes). It is envisaged that patients with a specific endotype present themselves within phenotypic clusters of disease, and because of the mechanistic differentiation, show response to specific therapies. Endotypes consist of subsets of the disease itself, rather than biological constructs which may or may not progress to disease (14). This approach has been used to describe subgroups in asthma (15), sepsis (16–19), trauma (20), and acute respiratory distress syndrome (ARDS) (21).
In clinical practice, baseline comorbidities and/or initial lab values do not explain the full range of COVID-19 presentations that are seen. We hypothesize that COVID-19 endotypes identified based on observable characteristics of the entire hospitalization (age and a representation of laboratory values) will reveal unexpected clinical courses and outcomes that defy prediction using classic risk factors. This approach is in contrast to some initial reports of clustering COVID-19 patients including using initial laboratory values and clinical variables collected in the first 24 (22) and 72 h (23); clustering patients by demographics, comorbidity, and maximum laboratory value (24) and using principal component analysis (PCA) and k-means of 18 initial laboratory values resulting in six values used in final analysis (25). Additionally, clusters have been created from initial ICU clinical data for patients with COVID-19 ARDS (26) and from ICU patients using demographics, initial ICU labs, and other clinical variables (27). Finally, there have been descriptions of a hyperinflammatory phenotype identified by initial admission labs (28) or serial labs using cluster analysis of three laboratory values (29).
In this study, we sought to uncover endotypes of the hospitalized COVID-19 patient population using a robust clustering method (consensus clustering of ensemble classification) on patient age and laboratory values over the course of hospital admission. These endotypes were examined for insights into comorbidities, expected clinical courses, and outcomes including intubation, length of stay (LOS), and mortality.
Methods
Participants
Adults (18 years-old or older) admitted consecutively to a telemetry capable bed at NewYork-Presbyterian Hospital/Columbia University Irving Medical Center were included in the study if they had a positive SARS-COV-2 nasopharyngeal PCR test during their inpatient admission and were discharged between March 12, 2020 to July 15, 2020. Patients with multiple admissions with a positive SARS-COV-2 nasopharyngeal PCR test only had data included from the first admission. If a patient had a positive SARS-COV-2 test (any type) more than 21 days before the admission, the patient was excluded. Patients were identified prospectively for inclusion in the study cohort but had their laboratory information, outcomes, and past medical history retrospectively collected. The collection of clinical data was done before clustering, so the investigators were blinded to endotype at the time of data collection. This study was approved by the Columbia University Institutional Review Board.
Features Used for Clustering
The features that have been shown to be correlated to clinical course or outcomes of COVID-19 were considered. Laboratory values and age were used to identify endotypes (complete list available in Supplementary Material 1). Both the median and the IQR of all lab values for a patient during admission were used as features. Features missing more than 40% of patients were excluded from analysis.
Variables Used to Examine the Resulting Endotypes
Patient disposition was the primary outcome. Intubation status, length of intubation, length of stay, patient age, race, sex, comorbidities, and treatment with medications commonly used with COVID-19 patients were collected (complete list available in Supplementary Material 2).
Statistical Analyses
A schematic presentation of data collection and analysis can be seen in Figure 1. To discover endotypes, we relied on cluster analysis, which generally divides datasets into groups by minimizing the intra-group distance while maximizing the inter-group distance. Instead of using a single clustering algorithm, here we employed ensemble classification (30) by running multiple clustering algorithms (K-mean, Birch, Gaussian Mixture Model, and Agglomerative clustering) and integrating their results. Then, we applied consensus clustering (31) to the results of ensemble classification. Consensus clustering is a robust approach that relies on multiple iterations of the sampled dataset to derive more stable and meaningful clusters and has been widely used to identify biologically meaningful clusters. In our work, the consensus of the ensemble clustering was implemented with 50 bootstraps and 80% of the data.
Figure 1. Data collection and analysis schematic. Patients with positive SARS-COV-2 tests that were discharged between March 12, 2020 to July 15, 2020 were included in the study. Labs during hospitalization (median and IQR) and age were the features used for clustering and endotype discovery. Once endotypes were identified, they were analyzed for differences in demographics, outcomes, comorbidities, and treatments.
The stability of consensus matrices (when cluster number K changed from 2 to 10) were measured by obtaining their cumulative distribution function (CDF) as described by Monti et al. (31). Then for each K value, proportion of increase in area under the CDF (ΔK), Calinski Harabasz score (CH) (32) and Davies Bouldin score (DBS) (33) were calculated and compared to determine the optimal number of clusters. Finally, to visualize the underlying structure of the data, we generated the data dendrograms by applying hierarchical clustering on the consensus matrices. Pseudocode of our clustering approach is provided in Supplementary Material 3.
To compare the differences between endotypes, the Kruskal-Wallis test (34) and Dunn's multiple comparison test (35) were used for continuous variables, and chi-square tests were used for categorical variables. A significant p-value was defined as <0.05. The analysis was performed in MATLAB™ (The Math Works, Inc., Natick, MA) and Python (www.python.org) where we used Opensemble library (36) to perform the consensus clustering.
Results
Five hundred forty-four patients were identified prospectively for inclusion in the study. Sixteen patients were missing all laboratory data and therefore were excluded from analysis, leaving 528 patients in the final cohort. Baseline characteristics of the final cohort, their comorbidities and hospital characterizations are outlined in Table 1. In the study cohort, the median age was 66 (IQR 55-74), 209 (40%) were female, 103 (19.5%) were African American or Black, 1 (0.2%) was American Indian or Alaska Nation, 7 (1.3%) were Asian, 2 (0.4%) were Native Hawaiian or Other Pacific Islander, 179 (33.9%) were other combinations not described, 119 (22.5%) were White, and 117 (22.5%) declined to specify race. In the cohort, 223 (42.2%) were discharged home, 54 (10.2%) were discharged to rehab, 129 (24.4%) were discharged to a skilled nursing facility, 25 (4.7%) were discharged to hospice, and 97 (18.4%) died in the hospital. Comorbid CKD, ESRD, HTN, and DM were higher in endotype 2 and 4. Length of stay was a median of five days in endotype 1, nine days in endotype 2, 41 days in endotype 3 and 37 days in endotype 4. Percent of patients intubated was 1% in endotype 1, 6% in endotype 2, 72% in endotype 3, and 85% in endotype 4.
Endotype Descriptions
Features missing in more than 40% of patients were excluded from further analysis: blood pH, blood pCO2, blood pO2, β-d-Glucan, ionized calcium, and fibrinogen. After considering cluster quality and stability by examining CDF plot, measured ΔK, CH, DBS, and the underlying structure of the data using dendrograms (Supplementary Material 4), we opted for K = 4 which identified four endotypes.
Median values of the clustering features for each of the four endotypes are outlined in Table 2, Supplementary Material 5. All of the features were significantly different over the endotypes except for median bilirubin and age (p > 0.05). Characteristics of the endotypes are outlined in Table 3. Some comorbidities varied significantly across endotypes (i.e., CKD, ESRD, HTN, DM, COPD, heart failure with reduced ejection fraction [HFrEF], and obesity), while others (asthma, hyperlipidemia, HIV infection, history of stroke, heart failure with preserved ejection fraction, and heart failure with unknown EF) did not differ significantly. Treatments differed by endotype (p < 0.05) except for remdesivir and prednisone. Mortality and discharge from hospital rates also varied by endotype (Figure 2). Paired comparisons of characteristics are provided in Supplementary Material 6. A summary of the four endotypes is shown in Figure 3.
Figure 2. Cumulative survival and discharge out of hospital by endotype. For each endotype, the cumulative survival and cumulative discharge of surviving patients from hospital is displayed in days from admission.
Endotype 1 patients had a median age of 68 years, had the most women (46%), the lowest prevalence of mortality (1%), shortest hospital length-of-stay (median: 5 days), and fewest intubated patients (1%). This endotype had the lowest prevalence of HTN and DM and greatest prevalence of COPD. Endotype 1 patients had the lowest inflammatory markers (ferritin, IL-6, CRP, ESR, LDH), lowest infectious markers (WBC, procalcitonin), and lowest degree of coagulopathy (PT and PTT, but not significantly < endotype 2). Endotype 1 patients received the least of any endotype of the reviewed medications (except for enoxaparin) but overall had similar medication use as endotype 2 (except for hydroxychloroquine and methylprednisolone).
Endotype 2 patients had a median age of 67.5 years, included approximately the cohort average of women (42%), second-lowest mortality (6%), relatively short hospital length-of-stay (median: 9 days), and second-fewest intubations (6%). This patient subgroup had the most comorbidities (CKD, ESRD, HTN, DM, and HFrEF). Endotype 2 patients had similar inflammatory markers to endotype 3 (ferritin, CRP, ESR, and LDH but not IL-6 which was significantly lower than endotype 3), second lowest infectious markers (WBC and procalcitonin, although procalcitonin was not significantly < endotype 3), and second least degree of coagulopathy (PT and PTT, but not significantly more than endotype 1). Endotype 2 patients received less of the reviewed medications than endotypes 3 and 4 except for enoxaparin which was not significantly < endotype 4.
Endotype 3 patients had a median age of 66 years, included approximately the cohort average of women (42%), exhibited a mortality of 21%, had the longest hospital length-of-stay (median: 41 days), and had the second-highest prevalence of intubation (72%). Patients in this endotype had a relatively low number of comorbidities. Endotype 3 patients had similar inflammatory markers as endotype 2 (ferritin, CRP, ESR, and LDH, but not IL-6 which was significantly higher), second-highest infectious markers (WBC and procalcitonin, although procalcitonin was not significantly > endotype 2), and second-highest coagulopathy markers (PT and PTT, but PT was not significantly < endotype 4). Endotype 3 patients received reviewed medications at similar rates as patients in endotype 4 (except for enoxaparin, heparin, and hydrocortisone).
Endotype 4 patients had a median age of 64 years, included the fewest women (27%), greatest degree of mortality (43%), a fairly long hospital length-of-stay (median 37 days), and were the most intubated (85%). This endotype had moderate amounts of CKD and ESRD, higher amounts of HTN, and the most obesity. Endotype 4 patients had the highest inflammatory markers (ferritin, LDH were significantly higher than endotype 3 while IL-6 and CRP were similarly high as endotype 3), highest infectious markers (WBC, procalcitonin), and greatest degree of coagulopathy (PT and PTT, but PT was not significantly > endotype 3). The exception was ESR which was lower than endotypes 2 and 3. Endotype 4 patients received the most of the reviewed medications (except for enoxaparin and hydroxychloroquine). Of the medications, only hydrocortisone and heparin use were significantly more than in endotype 3.
Discussion
Our study has three main findings: first, four distinct groups of patients were identified though consensus clustering of ensemble classification using age and laboratory values over the entire hospitalization as features. The groups as a whole did not vary significantly by age or race but had differences in sex as well as comorbidities. We consider these patient subgroups to comprise endotypes (14) since the data used to segregate them include variables that are indicative of physiologic and inflammatory dysfunction. The endotypes were also treated with differing medications in the hospital. Endotype 1 and 2 exhibited low mortality and short length of stay. However, Endotype 2 had slightly worse outcomes and slightly higher inflammatory and organ damage markers. Endotypes 3 and 4 had more mortality and length of stay, with endotype 4 having a markedly high mortality at 43% and the highest levels markers of inflammation and end-organ dysfunction.
Second, we identified endotypes of COVID-19 patients with widely disparate outcomes that were not expected based on classic risk factors such as age, sex, and preexisting comorbidities (3, 6). We documented patients with lower-risk features who had worse courses than traditionally expected. Endotype 2 had the greatest number of comorbidities overall but a relatively low mortality. Focusing on comorbidities alone would have resulted in misclassification of endotype 2 patients. Along the same lines, endotype 3 had many fewer comorbidities than endotype 2, and yet endotype 3 had significantly worse outcomes. IL-6, d-dimer, and WBC are significantly higher in endotype 3 compared to endotype 2. Further examination of the different endotypes has potential to yield clinical and pathobiological insight into what is driving the vastly different clinical courses experienced by patients with COVID-19.
Third, consensus clustering of ensemble classification (37) supported the previously hypothesized existence of subgroups of COVID-19 manifestations. In part because elevated inflammatory markers such as C-reactive protein, ferritin, and IL-6 were associated with poor outcomes (38, 39), steroids were studied and proven effective at treating severe COVID-19 (5). Patients meeting a proposed criteria for COVID-19-associated hyperinflammatory syndrome (including fever; ferritin and d-dimer elevation; NLR elevation or anemia/thrombocytopenia; LDH or AST elevation; and IL-6, triglyceride, or CRP elevation) were shown recently to have higher risk of requiring mechanical ventilation and higher risk of mortality (13). The endotypes we identified that have higher levels of circulating inflammatory markers have worse outcomes than patient clusters with lower inflammatory markers. This appears to hold true even when patients are intubated, such as in endotypes 3 and 4 in which a higher number of patients were intubated, but where there were notably higher mortality and inflammatory markers in endotype 4. Endotype 4 patients also had notably higher procalcitonin levels, a potential indication that these patients with higher inflammatory markers may have experienced more (or more severe) bacterial infections.
Identification of endotypes has several potential useful functions. Endotypes may point to unique pathobiologic mechanisms of disease that warrant further investigation in each specific subset of patients. Different endotypes may respond differently to treatments and may explain the heterogeneity of disease course. Examining endotypes for differential response to treatments could identify subsets of patients where treatments are beneficial. If endotypes can be identified early in disease course, endotypes can offer prognostic and clinical management information. Future studies will need to validate these endotypes.
There are several limitations to our study. First, this is a single-center study that prospectively collected data from patients admitted to telemetry capable beds. We have not validated the endotypes in the setting of more recent SARS-COV-2 variants. However, in the setting of this fast-moving disease, validation of endotypes in the setting of the most recent variant will continue to be a challenge for any large COVID-19 cohort study. Second, there were some lab variables with a high amount of missing data. These variables were dropped which may have introduced some bias. Third, standard of care treatments for patients with COVID-19 changed over time. The treatments each endotype received may have been changing over time. Dosing data for medications was not available, therefore anticoagulation medications were not classified as prophylactic or therapeutic. Fourth, the admission criteria for patients with COVID-19 may have changed over time.
In conclusion, disease endotypes have the potential to describe a subset of patients that are undergoing shared biologic processes resulting in a similar phenotype of disease and may identify groups of patients with different clinical courses and responses to therapy. However, having certain high or low risk features does not guarantee association with a certain outcome; rather, patients with certain features appear to have one of multiple different clinical courses. In this cohort of patients hospitalized with COVID-19, we identified four unique endotypes of patients by using clustering of laboratory values throughout the hospitalization as well as patient age. The endotypes had differences in inflammatory markers, infectious markers, evidence of end-organ dysfunction, comorbidities, and outcomes. Further work is needed to validate these endotypes in other cohorts and study endotype differences to treatment response.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Materials, further inquiries can be directed to the corresponding author/s.
Ethics Statement
The studies involving human participants were reviewed and approved by Columbia University. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.
Author Contributions
BR, MM, KT, JC, MP, GC, YV, SA, and SP were involved in the conception and design of the study. BR, MM, KT, KD, MP, GC, YV, SA, and SP were involved in the acquisition and/or analysis. BR, MM, KT, MP, GC, YV, SA, and SP interpreted the data. All authors were involved in drafting the manuscript or critical appraisal and approved the final version of the manuscript.
Funding
The following grants supported this study: NIH R21 NS113055 (SP), AHA 20POST35210653 (MM), and NIH U01EB021960 (YV).
Conflict of Interest
YV is a co-founder of, and stakeholder in Immunetrics, Inc.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2021.770343/full#supplementary-material
References
1. Oran DP, Topol EJ. Prevalence of asymptomatic SARS-CoV-2 infection. Ann Intern Med. (2020) 173:362–7. doi: 10.7326/M20-3012
2. Gandhi RT, Lynch JB, del Rio C. Mild or moderate COVID-19 N Engl J Med. (2020) 383:1757–66. doi: 10.1056/NEJMcp2009249
3. Wu C, Chen X, Cai Y, Xia J, Zhou X, Xu S, et al. Risk factors associated with acute respiratory distress syndrome and death in patients with coronavirus disease 2019 pneumonia in Wuhan, China. JAMA Intern Med. (2020) 180:934–43. doi: 10.1001/jamainternmed.2020.0994
4. Cummings MJ, Baldwin MR, Abrams D, Jacobson SD, Meyer BJ, Balough EM, et al. Epidemiology, clinical course, and outcomes of critically ill adults with COVID-19 in New York City: a prospective cohort study. Lancet. (2020) 395:1763–70. doi: 10.1016/S0140-6736(20)31189-2
5. Carfì A, Bernabei R, Landi F, Group ftGAC-P-ACS. Persistent symptoms in patients after acute COVID-19. JAMA. (2020) 324:603–5. doi: 10.1001/jama.2020.12603
6. Chen T, Wu D, Chen H, Yan W, Yang D, Chen G, et al. Clinical characteristics of 113 deceased patients with coronavirus disease 2019: retrospective study. BMJ. (2020) 368:m1091. doi: 10.1136/bmj.m1091
7. Huang Y, Lu Y, Huang YM, Wang M, Ling W, Sui Y, et al. Obesity in patients with COVID-19: a systematic review and meta-analysis. Metabolism. (2020) 113:154378. doi: 10.1016/j.metabol.2020.154378
8. Apicella M, Campopiano MC, Mantuano M, Mazoni L, Coppelli A, Del Prato S. COVID-19 in people with diabetes: understanding the reasons for worse outcomes. Lancet Diabetes Endocrinol. (2020) 8:782–92. doi: 10.1016/S2213-8587(20)30238-2
9. Higham A, Mathioudakis A, Vestbo J, Singh D. COVID-19 and COPD: a narrative review of the basic science and clinical outcomes. Eur Respir Rev. (2020) 29:200199. doi: 10.1183/16000617.0199-2020
10. Bruchfeld A. The COVID-19 pandemic: consequences for nephrology. Nat Rev Nephrol. (2021) 17:81–2. doi: 10.1038/s41581-020-00381-4
11. Wynants L, van Calster B, Bonten MMJ, Collins GS, Debray TPA, de Vos M, et al. Prediction models for diagnosis and prognosis of COVID-19 infection: systematic review and critical appraisal. BMJ. (2020) 369:m1328. doi: 10.1101/2020.03.24.20041020
12. Henry BM, de Oliveira MHS, Benoit S, Plebani M, Lippi G. Hematologic, biochemical and immune biomarker abnormalities associated with severe illness and mortality in coronavirus disease 2019 (COVID-19): a meta-analysis. Clin Chem Lab Med. (2020) 58:1021–8. doi: 10.1515/cclm-2020-0369
13. Webb BJ, Peltan ID, Jensen P, Hoda D, Hunter B, Silver A, et al. Clinical criteria for COVID-19-associated hyperinflammatory syndrome: a cohort study. Lancet Rheumatol. (2020) 2:e754–e63. doi: 10.1016/S2665-9913(20)30343-X
14. Anderson GP. Endotyping asthma: new insights into key pathogenic mechanisms in a complex, heterogeneous disease. Lancet. (2008) 372:1107–19. doi: 10.1016/S0140-6736(08)61452-X
15. Kuruvilla ME, Lee FE-H, Lee GB. Understanding asthma phenotypes, endotypes, and mechanisms of disease. Clin Rev Allergy Immunol. (2019) 56:219–33. doi: 10.1007/s12016-018-8712-1
16. 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:2003–17. doi: 10.1001/jama.2019.5791
17. 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:816–26. doi: 10.1016/S2213-2600(17)30294-1
18. Sweeney TE, Azad TD, Donato M, Haynes WA, Perumal TM, Henao R, et al. Unsupervised analysis of transcriptomics in bacterial sepsis across multiple datasets reveals three robust clusters. Crit Care Med. (2018) 46:915–25. doi: 10.1097/CCM.0000000000003084
19. Bhavani SV, Carey KA, Gilbert ER, Afshar M, Verhoef PA, Churpek MM. Identifying novel sepsis subphenotypes using temperature trajectories. Am J Respir Crit Care Med. (2019) 200:327–35. doi: 10.1164/rccm.201806-1197OC
20. Schimunek L, Lindberg H, Cohen M, Namas RA, Mi Q, Yin J, et al. Computational derivation of core, dynamic human blunt trauma inflammatory endotypes. Front Immunol. (2020) 11:589304. doi: 10.3389/fimmu.2020.589304
21. Prescott HC, Calfee CS, Thompson BT, Angus DC, Liu VX. Toward smarter lumping and smarter splitting: rethinking strategies for sepsis and acute respiratory distress syndrome clinical trial design. Am J Respir Crit Care Med. (2016) 194:147–55. doi: 10.1164/rccm.201512-2544CP
22. Vasquez CR, Gupta S, Miano TA, Roche M, Hsu J, Yang W, et al. Identification of distinct clinical subphenotypes in critically ill patients with COVID-19. Chest. (2021) 160:929–43. doi: 10.1016/j.chest.2021.04.062
23. Lusczek ER, Ingraham NE, Karam BS, Proper J, Siegel L, Helgeson ES, et al. Characterizing COVID-19 clinical phenotypes and associated comorbidities and complication profiles. PLoS ONE. (2021) 16:e0248956. doi: 10.1371/journal.pone.0248956
24. da Silva JF, Hernandez-Romieu AC, Browning SD, Bruce BB, Natarajan P, Morris SB, 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. (2021) 8:ofaa596. doi: 10.1093/ofid/ofaa596
25. Ye W, Lu W, Tang Y, Chen G, Li X, Ji C, et al. Identification of COVID-19 clinical phenotypes by principal component analysis-based cluster analysis. Front Med. (2020) 7:570614. doi: 10.3389/fmed.2020.570614
26. Ranjeva S, Pinciroli R, Hodell E, Mueller A, Hardin CC, Thompson BT, et al. Identifying clinical and biochemical phenotypes in acute respiratory distress syndrome secondary to coronavirus disease-2019. EClinicalMedicine. (2021) 34:100829. doi: 10.1016/j.eclinm.2021.100829
27. Rodriguez A, Ruiz-Botella M, Martin-Loeches I, Jimenez Herrera M, Sole-Violan J, Gomez J, 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. doi: 10.1186/s13054-021-03487-8
28. Chen H, Xie J, Su N, Wang J, Sun Q, Li S, et al. Corticosteroid therapy is associated with improved outcome in critically ill patients with COVID-19 with hyperinflammatory phenotype. Chest. (2021) 159:1793–802. doi: 10.1016/j.chest.2020.11.050
29. Lee EE, Song KH, Hwang W, Ham SY, Jeong H, Kim JH, et al. Pattern of inflammatory immune response determines the clinical course and outcome of COVID-19: unbiased clustering analysis. Sci Rep. (2021) 11:8080. doi: 10.1038/s41598-021-87668-z
30. Rokach L. Ensemble-based classifiers. Artif Intell Rev. (2010) 33:1–39. doi: 10.1007/s10462-009-9124-7
31. Monti S, Tamayo P, Mesirov J, Golub T. Consensus clustering: a resampling-based method for class discovery and visualization of gene expression microarray data. Mach Learn. (2003) 52:91–118. doi: 10.1023/A:1023949509487
32. Caliński T, Harabasz J. A dendrite method for cluster analysis. Commun Stat. (1974) 3:1–27. doi: 10.1080/03610927408827101
33. Davies DL, Bouldin DW. A cluster separation measure. IEEE Trans Pattern Anal Mach Intell. (1979). PAMI-1:224–7. doi: 10.1109/TPAMI.1979.4766909
34. Kruskal WH, Wallis WA. Use of Ranks in One-Criterion Variance Analysis. J Am Stat Assoc. (1952) 47:583–621.
35. Dunn OJ. Multiple comparisons among means. J Am Stat Assoc. (1961) 56:52–64. doi: 10.1080/01621459.1961.10482090
36. Ronan T, Anastasio S, Qi Z, Tavares PHSV, Sloutsky R, Naegle KM. OpenEnsembles: a python resource for ensemble clustering. J Mach Learn Res. (2018) 19:1–6.
37. Merad M, Martin JC. Pathological inflammation in patients with COVID-19: a key role for monocytes and macrophages. Nat Rev Immunol. (2020) 20:355–62. doi: 10.1038/s41577-020-0331-4
38. Huang C, Wang Y, Li X, Ren L, Zhao J, Hu Y, et al. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet. (2020) 395:497–506. doi: 10.1016/S0140-6736(20)30183-5
Keywords: COVID-19, cluster analysis, endotype, phenotype, machine learning, treatment, survival
Citation: Ranard BL, Megjhani M, Terilli K, Doyle K, Claassen J, Pinsky MR, Clermont G, Vodovotz Y, Asgari S and Park S (2021) Identification of Endotypes of Hospitalized COVID-19 Patients. Front. Med. 8:770343. doi: 10.3389/fmed.2021.770343
Received: 03 September 2021; Accepted: 05 October 2021;
Published: 11 November 2021.
Edited by:
Zisis Kozlakidis, International Agency for Research on Cancer (IARC), FranceReviewed by:
Shant Ayanian, Mayo Clinic, United StatesWanlin Jiang, Massachusetts General Hospital and Harvard Medical School, United States
Sridhar Chilimuri, BronxCare Health System, United States
Copyright © 2021 Ranard, Megjhani, Terilli, Doyle, Claassen, Pinsky, Clermont, Vodovotz, Asgari and Park. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Soojin Park, sp3291@cumc.columbia.edu