Abstract
Background There is limited understanding of heterogeneity in outcomes across hospitalized patients with coronavirus disease 2019 (COVID-19). Identification of distinct clinical phenotypes may facilitate tailored therapy and improve outcomes.
Objective Identify specific clinical phenotypes across COVID-19 patients and compare admission characteristics and outcomes.
Design, Settings, and Participants Retrospective analysis of 1,022 COVID-19 patient admissions from 14 Midwest U.S. hospitals between March 7, 2020 and August 25, 2020.
Methods Ensemble clustering was performed on a set of 33 vitals and labs variables collected within 72 hours of admission. K-means based consensus clustering was used to identify three clinical phenotypes. Principal component analysis was performed on the average covariance matrix of all imputed datasets to visualize clustering and variable relationships. Multinomial regression models were fit to further compare patient comorbidities across phenotype classification. Multivariable models were fit to estimate the association between phenotype and in-hospital complications and clinical outcomes.
Main outcomes and measures Phenotype classification (I, II, III), patient characteristics associated with phenotype assignment, in-hospital complications, and clinical outcomes including ICU admission, need for mechanical ventilation, hospital length of stay, and mortality.
Results The database included 1,022 patients requiring hospital admission with COVID-19 (median age, 62.1 [IQR: 45.9-75.8] years; 481 [48.6%] male, 412 [40.3%] required ICU admission, 437 [46.7%] were white). Three clinical phenotypes were identified (I, II, III); 236 [23.1%] patients had phenotype I, 613 [60%] patients had phenotype II, and 173 [16.9%] patients had phenotype III.
When grouping comorbidities by organ system, patients with respiratory comorbidities were most commonly characterized by phenotype III (p=0.002), while patients with hematologic (p<0.001), renal (p<0.001), and cardiac (p<0.001) comorbidities were most commonly characterized by phenotype I. The adjusted odds of respiratory (p<0.001), renal (p<0.001), and metabolic (p<0.001) complications were highest for patients with phenotype I, followed by phenotype II. Patients with phenotype I had a far greater odds of hepatic (p<0.001) and hematological (p=0.02) complications than the other two phenotypes. Phenotypes I and II were associated with 7.30-fold (HR: 7.30, 95% CI: (3.11-17.17), p<0.001) and 2.57-fold (HR: 2.57, 95% CI: (1.10-6.00), p=0.03) increases in the hazard of death, respectively, when compared to phenotype III.
Conclusion In this retrospective analysis of patients with COVID-19, three clinical phenotypes were identified. Future research is urgently needed to determine the utility of these phenotypes in clinical practice and trial design.
Introduction
The coronavirus disease 2019 (COVID-19), a disease caused by the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), has infected over 18 million and led to over 700,000 deaths since first appearing in late 2019.1 Researchers are rapidly attempting to understand the natural history of and immune response to COVID-19.2 Despite intense research since the arrival of this novel coronavirus3, only one pharmaco-therapeutic agent, dexamethasone, has been associated with reduced mortality in at-risk individuals.4 COVID-19 results in a constellation of symptoms, laboratory derangement, immune dysregulation, and clinical complications.5
Emergency department presentation varies widely, suggesting various clinical phenotypes exist and, importantly, it is likely these various phenotypes respond differently to treatment. To illustrate, two early phenotypes of respiratory failure likely exist in COVID-19. A classic ARDS phenotype exists with poorly compliant lungs and poor gas exchange; however, a phenotype with normal lung compliance also exists in COVID-19 and is hypothesized to be driven by shunting secondary to pulmonary microthrombi.6,7 An intricate, multidimensional view is required to adequately understand the disease and account for the variation in clinical outcomes. Furthermore, patients could benefit from phenotype-specific medical care, which may differ from established standards of care
Despite this need, few studies have characterized COVID-19 clinical phenotypes and evaluated their association with complications and clinical outcomes. The aim of this study was to characterize clinical phenotypes in COVID-19 according to disease-system factors using electronic health record (EHR) data pooled from 14 U.S. Midwest hospitals between March 7, 2020 and August 25, 2020.
Methods
Data Collection
The data source for this study included EHR reports from 14 U.S. Midwest hospitals and 60 primary care clinics. Patient and hospital-level data were available for 7,538 patients with PCR-confirmed COVID-19. Of these, 1,022 required hospital admission and were included in this analysis. The database included all comorbidities reported since March 29, 1997 for each patient and prior to their COVID-19 diagnosis. The database also included home medications, laboratory values, clinic visits, social history, and patient demographics (age, gender, race/ethnicity, language spoken, zip code, socioeconomic status indicators). For each COVID-19 hospitalization the database included all laboratory values, vitals, orders, medications, complications, length of stay, and hospital disposition. State death certificate data was linked with the database to enable capture of out-of-hospital death. Additionally, the database allowed linkage across the 14 hospitals, facilitating the tracking of transfers.
This study was approved by the University of Minnesota institutional review board (STUDY00001489).
Participants
Patient-level data were obtained from the COVID-19 database from March 7, 2020 to August 25, 2020. The inclusion criterion was as follows: PCR-positive COVID-19 test requiring inpatient hospital admission to one of the 14 hospitals providing data. No hospitalized patients were excluded in this analysis to maximize generalizability. Follow-up data were available for a minimum of two weeks following admission for all patients.
Clinical Variables for Phenotyping
We selected 33 variables for clustering based on their association with COVID-19 mortality, known COVID-19 pathophysiology, and presence in the database (no more than 50% missingness).8-11 The following variables were included: age, body mass index (BMI), heart rate, respiratory rate, oxygen saturation, pulse pressure, systolic blood pressure, total protein, red cell distribution width, mean corpuscular volume, alkaline phosphatase, calcium, anion gap, bicarbonate, hematocrit, aspartate aminotransferase, glucose, absolute monocyte count, absolute neutrophil count, absolute lymphocyte count, white blood cell count, platelet, albumin, bilirubin, international normalized ratio (INR), lactate dehydrogenase, potassium, sodium, D-dimer, hemoglobin, C-reactive protein (CRP), creatinine, and gamma gap. For each variable we selected the first recorded value within the first 72 hours of the emergency department (ED) presentation that ultimately resulted in their hospitalization.
Comorbidities
We selected 68 comorbidities documented for each patient from March 29, 1997 preceding their COVID-19 hospital admission in their electronic health record (Supplemental Table 3). All comorbidities were identified based on ICD-9, ICD-10, or problem list documentation within the electronic health record. An indicator variable was created for each comorbidity to denote the presence of the selected ICD-9, ICD-10, or problem list documentation at any time in the medical record. To facilitate analysis, comorbidities were grouped by organ system into the following categories: cardiac, respiratory, hematologic, metabolic, renal, hepatic, autoimmune, cancer, and cerebrovascular disease. Complications and Clinical Outcomes We selected 30 in-hospital complications measured during each patient’s hospital stay for COVID-19 categorized into the following systems: cardiovascular, respiratory, hematologic, renal, hepatic, metabolic, and infectious (Supplemental Table 4). If applicable, complications could span multiple organ system variables. For example, ventilator associated pneumonia was included in both infectious and respiratory complications. Additional clinical outcomes included hospital length of stay (LOS), need for intensive care unit (ICU) admission, need for mechanical ventilation, and mortality. Mortality was defined as any in-hospital or out-of-hospital death based on death certificate data. All complications and outcomes were followed for a minimum of 2 weeks following hospital admission.
Statistical Analysis
The overall rate of missingness of the 33 variables used for phenotyping, which included the first vitals and labs recorded for each inpatient within 72 hours of admission, was 19% (range 0% - 50%). We imputed missing values using multivariate imputations by chained equations implemented with the mice package (v.3.10.0).12,13 Data were log-transformed before imputing missing values with predictive mean matching. A total of 40 imputed datasets were generated. The diceR package (v.1.0.0)14 was used to perform k-means-based consensus clustering on each imputed dataset using 80% subsamples and 1,000 iterations. We considered grouping patients into 2-7 phenotypes and determined the optimal number was 3 by evaluating the consensus cumulative distribution function (CDF) plot, the delta area plot, and the consensus matrix heatmap. These figures were generated using the consensus clustering results for each imputed dataset, and all figures were qualitatively similar across datasets. For visualization purposes, these images are provided for a randomly selected dataset in Supplemental Figures 1-4. The final assignment of each patient into one of the three phenotypes was determined by majority voting across the 40 consensus clustering results. Principal component analysis (PCA) was performed on the average covariance matrix to visualize the relationships among the three phenotypes and assess variable contributions.15
Continuous variables were summarized using the median and interquartile range (IQR) and compared across phenotypes using a Kruskal-Wallis test. Categorical characteristics and outcomes were summarized using counts and proportions and compared across phenotypes using a Pearson’s chi-squared test or Fisher’s exact test. Multinomial regression models were fit to further compare patient comorbidities across phenotype classification.
We next evaluated the relationship between phenotype and subsequent outcomes using both unadjusted and adjusted models. The adjusted models included sex,16,17 race and ethnicity (white, Black, Asian, Hispanic, other, not reported),18 and Elixhauser Comorbidity Index,19 since these are known risk factors for the outcomes of interest and were not included in the clustering analysis. The associations between phenotype and complications, ICU admission and need for mechanical ventilation, were estimated using logistic regression models. Mortality was compared across phenotypes using Cox proportional hazard models and patients were censored at the last date of data collection, August 25, 2020. Hospital length of stay was compared across phenotypes using negative binomial regression models. The primary negative binomial model included individuals who died during hospitalization for whom length of stay was defined as the number of days until death. We performed a sensitivity analysis to assess the impact of mortality as a competing risk by refitting the length of stay model after removing the 127 patients who died. Two-sided p-values < 0.05 were considered statistically significant. P-values were not adjusted for multiple comparisons. Visualizations of comorbidities, complications, and outcomes by clinical phenotype were performed using the circlize package for R.20 Comorbidities and complications were grouped into separate organ systems and the prevalence of each complication/comorbidity type was calculated as a percentage for each phenotype. All analyses were conducted using R version 3.6.321 and Stata version 16.1 (StataCorp).
Result
The database included 1,022 patients requiring hospital admission with COVID-19. Among these patients, the median age was 62.1 [IQR: 45.9, 75.8] years; 481 [48.6%] male, 412 [40.3%] required ICU admission). Additionally, 437 [46.7%] were white, 188 [20.1%] were Black, 159 [17.0%] were Asian, 103 [11.0%] were Hispanic, 20 [2.1%] reported other race, and 28 [2.9%] did not report. Three clinical phenotypes were identified (I, II, III); 236 [23.1%] patients had phenotype I, 613 [60%] patients had phenotype II, and 173 [16.9%] patients had phenotype III.
Variable Contributions to Clustering
The first two principal components (PCs) from PCA were used to visualize the relationship between phenotypes. PC1 and PC2 captured approximately 11% and 9% of the variance in the clustering variables, respectively. Thirteen components were needed to explain 70% of the variance (Supplemental Figure 5). While phenotypes II and III overlay substantially, phenotype I is more clearly defined in the right-hand side of the score plot of the first two principal components (Figure 1). Notably, this figure shows that distinctions between phenotypes are primarily driven by variation in PC1 as opposed to PC2. The variable contributions to PC1 (Figure 2) demonstrate that the largest contributors to the variation in PC1 are from LDH, absolute neutrophil count, and D-dimer. These variables therefore prominently contribute to separating the three phenotypes as shown in the biplot (Figure 3). Univariate tests showed that D-dimer and neutrophil count are highest in phenotype I and albumin is highest in phenotype III. Other variables influential to phenotype clustering are lactate dehydrogenase (highest in I), C-reactive protein (highest in I), white cell count (highest in I), red cell distribution width (highest in I), bilirubin (highest in I), AST (highest in I), hematocrit (highest in III), and hemoglobin (highest in III).
Phenotype Characteristics
Differences across phenotypes with respect to patient demographics, admission vitals and labs, complications, comorbidities, and clinical outcomes are presented in Table 1. Patients with phenotype I were older than patients in phenotypes II and III (67.2 [52.9, 79.0] years vs. 60.9 [45.9, 75.4] and 58.6 [34.8, 71.3] years respectively, p < 0.001). Patients with phenotype III were more often female than patients with phenotype I or II (57.6% vs. 41.6% and 53.4%, respectively, p = 0.004). Patients with phenotype I were more less likely white (38.8% vs. 45.6% vs. 60.7%, respectively, p = 0.002) and more likely to be non-English speaking (47.9% vs. 39.2% vs. 23.7%, respectively, p <0.001). There were no statistically significant differences in BMI or socioeconomic status, as measured using the area deprivation index, between phenotypes (Table 1). Patients that presented with phenotype III had a more frequent history of smoking, alcohol abuse, neutropenia. Patients that presented with phenotype II had a less frequent history of hepatic disease than phenotypes I or III (Table 1).
When grouping comorbidities by organ system, cardiac (p <0.001), respiratory (p =0.002), hematologic (p <0.001), and renal (p <0.001) comorbidities were found to be significantly associated with phenotype. Cancer, hepatic, autoimmune, cerebrovascular, and metabolic comorbidities were not significantly associated with phenotype (Table 1, Supplemental Figure 6). Based on the estimated relative risk ratios, patients with renal (RRR 2.35; 95% CI 1.5-3.67; p <0.001), hematologic (RRR 2.64; 95% CI 1.75-3.98; p <0.001), and cardiac comorbidities (RRR 2.65; 95% CI: 1.68-4.17; p <0.001) were more likely to have phenotype I vs. III (Figure 4). Patients with respiratory comorbidities were 0.47 (95% CI: 0.31-0.72; p <0.001) times as likely to have phenotype I vs. III and 0.74 (95% CI: 0.52-1.04 p =0.09) times as likely to have phenotype II vs. III (Figure 4).
Association between Phenotype and Clinical Outcomes
Clinical phenotypes I and II were associated with increased odds of respiratory (I: OR: 2.98, 95% CI 1.58 - 5.59; II: OR: 2.32, 95% CI: 1.29-4.17; p<0.001), renal (I: OR: 7.04, 95% CI 3.11-15.9; II: OR: 2.57, 95% CI: 1.15-5.74; p <0.001), and metabolic (I: OR: 4.85, 95% CI: 2.78-8.45; II: OR: 2.57, 95% CI: 1.52-4.34; p <0.001) complications, compared to phenotype III after adjusting for sex, race, and Elixhauser Comorbidity Index (Supplemental Table 2). There was a trend towards increased odds of hematologic complications among patients with phenotype I (I: OR: 2.11, 95% CI: 0.99-4.48, p =0.05) compared to III. Phenotype was associated with hepatic complications (p <0.001); however, while phenotype I was associated with a 8.35-fold (OR: 8.35, 95% CI: 1.93-36.11, p < 0.001) increase in the odds of hepatic complication, phenotype II did not differ significantly from phenotype III (OR: 0.56, 95% CI: 0.10-3.09, p = 0.51). This is not surprising since only 4 individuals in phenotype II and 2 in phenotype III experienced hepatic complications during hospitalization (Table 1). Phenotype was also significantly associated with the rate of infectious complications (p <0.001) for phenotype 1 (OR 2.57, 95% CI 1.57 - 4.21; <0.001) but not did not reach statistical significance for phenotype 2 (OR 1.51, 95% CI 0.96 - 2.38; p = 0.07) (Supplemental Table 2 and Supplemental Figure 7).
Clinical phenotypes differed in odds of ICU admission (p <0.001) and mechanical ventilation (p <0.001), hospital LOS (p <0.001), and risk of mortality (<0.001) on adjusted analysis which accounted for sex, race, and Elixhauser Comorbidity Index (Table 2, Supplemental Figure 8). Controlling for these risk factors and compared to phenotype III, phenotypes I and II were associated with 7.88-fold (OR: 7.88, 95% CI: 4.65-13.37) and 2.32-fold (OR: 2.32, 95% CI: 1.46-3.68) increases in the odds of ICU admission, respectively. Phenotypes I and II were associated with 25.59-fold (OR: 25.59, 95% CI: 7.69,85.17) and 7.45-fold (OR: 7.45, 95% CI: 2.27-24.43) increases in the odds of requiring mechanical ventilation. Phenotypes I and II were associated with 1.74-fold (IRR: 1.74, 95% CI: 1.45-2.10, p<0.001) and 1.22-fold (IRR: 1.22, 95% CI: 1.05-1.43, p = 0.01) increases in hospital LOS. Phenotype I was associated with a 7.30-fold (HR: 7.30, 95% CI: 3.11-17.17, p <0.001) increase in risk of mortality, and Phenotype II had a 2.57-fold (HR: 2.57, 95% CI: 1.10-6.00, p=0.03) increase in the hazard of death compared to Phenotype 3. We performed a sensitivity analysis to assess the impact of mortality as a competing risk by fitting the LOS model before and after removing the 127 patients who died. The estimated effect sizes were similar between these two models (data not shown). Table 2 includes the LOS model with only survivors.
Discussion
This is one of the first studies to report on clinical phenotypes associated with COVID-19. We identified three clinical phenotypes for patients with COVID-19 on hospital presentation. Most patients presented with phenotype II, which is associated with a moderate course and an approximately 10% mortality. A subset of patients presented with the more severe phenotype I, which is associated with a staggering 27% mortality. Patients with cardiac, hematologic, and renal comorbidities were most likely to be characterized by phenotype I. Surprisingly, respiratory comorbidities appeared less related to phenotypes I or II and were most associated with phenotype III, which had the most indolent course. Despite this indolent course, patients with phenotype III had the highest rate of readmission which is likely in part due to the high survival rate. This also suggests patients with pre-existing respiratory comorbidities, while not at highest risk for mortality, may be at highest risk for long term sequalae following COVID-19. Patients that presented with phenotype I were most associated with the development of respiratory, hematologic, renal, metabolic, hepatic, and infectious complications. Surprisingly, cardiovascular complications did not significantly differ between phenotypes.
Elucidating patient risk factors and severe COVID-19 disease markers may allow early treatment implementation that may improve the patient’s outcome. Multiple studies have documented COVID-19 risk factors; however, most have done so from a homogenous lens. For example, a prospective cohort study from New York City identified that the most considerable risks for hospital admission were age, male sex, heart failure, chronic kidney disease, and high BMI.22 A large observational study conducted in the UK reported that increasing age, male gender, comorbidities such as cardiac disease, chronic lung disease, chronic kidney disease, and obesity were associated with higher mortality in COVID-19 positive patients admitted to the hospital.14 A study from China found that increased odds of in-hospital death due to COVID-19 were associated with older age, higher SOFA score and D-dimers > 1.0 μg/mL on admission.23 Another retrospective study reported that patients with severe COVID-19 disease and diabetes had increased leucocytes, neutrophils count, and increased C-reactive protein (CRP), D-dimers, fibrinogen levels.24 A systematic review and meta-analysis found that the biomarkers associated with increased mortality include higher CRP, higher D-dimers, increased creatinine, and lower albumin levels.25 However it is well known that patients do not have a singular natural history of disease. Multiple studies including this study found that only half of patients suffer a primarily respiratory disease.26,27 Patients suffer a constellation of cardiovascular, hematologic, renal, or hepatic progression of disease following COVID-19. It is likely patient baseline risk factors related to the virus,28 home medications,16,29 genetic predisposition30, race/ethnicity,18 and other factors predispose patients to one of the various clinical manifestations and natural history of COVID-19.
Treatment of hospitalized patients should be tailored based on the clinical courses most likely for a patient given their a priori risk. For example, phenotypes with a higher risk of thrombotic events, may benefit from more aggressive anticoagulation. Phenotypes more prone to infectious complications, may benefit from more targeted immunomodulation instead of broad and systemic steroid therapy. A key first step to evaluate these treatment decisions is to characterize and describe clinical phenotypes requiring hospitalization. In this analysis we identified three clinical phenotypes for patients that required hospitalization for COVID-19. Few studies to date have attempted to elucidate clinical phenotypes. One study attempted to characterize clinical phenotypes at ICU admission using a dataset of 85 critically ill patients.31 Similar to our analysis, they identified three distinct clinical phenotypes. Their low mortality cluster which they called cluster 1 was very similar to our phenotype III with a predominance of females, lower mortality rate, lower D-dimer and CRP levels. Similarly, their high mortality cluster was predominantly male, with elevated inflammation markers on ICU presentation. In this study, we not only characterized three clinical phenotypes, but extended findings outside of the ICU by characterizing the association of comorbidities with clinical phenotype and the association of clinical phenotypes with in-hospital complication and clinical outcomes.
Phenotype I can be termed the “Adverse phenotype” and was associated with the worst clinical outcomes. LDH, Absolute Neutrophil Count, D-dimer, AST, and CRP were most influential in phenotype I determination. The strong association of RDW with phenotype I was interesting. RDW was strongly associated with genetic age which is hypothesized to be a risk factor in Covid-19.30 As people age, variability in red blood cell volumes increases. Similarly, Gamma Gap, a marker of immunoglobulin levels, was elevated in all three phenotypes (median >3.5).32 However, patients with clinical phenotype I were noted to have the largest increase in Gamma Gap. In this scenario elevated Gamma Gap was likely an indicator of systemic inflammation and has been associated in other inflammatory disease processes with prognosis. Other groups have previously reported on the importance of the Absolute Neutrophil to Absolute Lymphocyte count, here we noted that ANC/ALC was lowest for phenotype III and highest for phenotype I, in line with previous reports. Patients with cardiac, hematologic and renal comorbidities were most prone to develop phenotype I. Phenotype I was associated with numerous complications (hematologic, hepatic, metabolic, renal, respiratory, and infectious) when compared to other phenotypes. It is interesting to note despite a higher rate of baseline cardiac comorbidities phenotype I was not associated with increased cardiac complications.
Phenotype III was associated with the best clinical outcomes and can be termed the “Favorable Phenotype”. Surprisingly, patients with phenotype III had a very high rate of respiratory comorbidities and the best clinical outcomes. What is most surprising is despite the lowest complication rate and mortality, this phenotype was associated with a greater than 10% rate of hospital readmission. It is possible that patients pre-existing respiratory comorbidities predisposed them to longer term sequelae which may have resulted in this readmission rate, although additional studies are needed to better elucidate these findings, specifically controlling for differences in survival. Patients with respiratory comorbidities such as asthma and COPD routinely use medications which may be protective in SARS-CoV-2 pathogenesis which may explain this protective effect. For example, our group has previously identified reduced mortality in COVID-19 for patients with asthma treated with beta2-agonists.16 Patients with phenotype III were more likely to use inhaled steroids, nasal fluticasone, albuterol, and antihistamines.
Ultimately, a deeper investigation into clinical phenotypes and associated genomic, transcriptomic, and proteomic is needed. The ability to classify patients into clinical phenotypes can facilitate the linkage of exome data to better understand SARS-CoV-2 pathogenesis and natural history. Understanding the COVID-19 severity, the biomarkers, and the risk factors is paramount during the COVID-19 pandemic.
Limitations
Our study has several limitations, including that this is a retrospective study and therefore results may be biased or subject to residual confounding. Second, patients were followed for variable lengths of time. Patients that were admitted in March 2020 thus had approximately 5 months of follow-up whereas patients admitted in late August had limited time. We accounted for this by conducting a Cox proportional hazard analysis when analyzing in- and out-of hospital mortality. Additionally, when the data were pulled, only 54 patients (5%) remained hospitalized. While most patients developed complications within their first 2 weeks of hospital admission, it is possible that they may still develop clinical complications which is not reflected in this analysis.
Conclusion
In this retrospective analysis of patients with COVID-19, three clinical phenotypes were identified. Future research is urgently needed to determine the utility of these phenotypes in clinical practice and trial design.
Data Availability
Data are not available.
Conflicts of Interest and Source of Funding
NIH NHLBI T32HL07741 (NEI)
This research was supported by the Agency for Healthcare Research and Quality (AHRQ) and Patient-Centered Outcomes Research Institute (PCORI), grant K12HS026379 (CJT) and the National Institutes of Health’s National Center for Advancing Translational Sciences, grant KL2TR002492 and UL1TR002494.
NIH NHLBI T32HL129956 (JP, LS)
Author Contribution
Concept and design: All authors
Acquisition, analysis, or interpretation of data: Lusczek, Proper, Siegel, Helgeson, Usher, Tignanelli
Drafting of the manuscript: All authors
Critical revision of the manuscript for important intellectual content: All authors
Conflicts of Interest and Funding Source(s)
NIH NHLBI T32HL07741 (NEI)
This research was supported by the Agency for Healthcare Research and Quality (AHRQ) and Patient-Centered Outcomes Research Institute (PCORI), grant K12HS026379 (CjT) and the National Institutes of Health’s National Center for Advancing Translational Sciences, grant KL2TR002492 and UL1TR002494.
NIH NHLBI T32HL129956 (JP, LS)
Reprints will not be available from the authors. All authors significantly contributed to developing, writing, and revising this manuscript.
Supplemental Figure 1: Consensus Cumulative Distribution Functions
Cumulative distribution functions (CDF) for a randomly selected imputed dataset are shown. A range of phenotypes (2-7) were considered, and the optimal choice of phenotypes is 3.
Supplemental Figure 2: Delta Area
The relative change in delta area under the cumulative distribution function is shown for the range of phenotypes (k=2-7) for a randomly selected imputed dataset. The optimal choice of phenotypes is 3.
Abbreviations: CDF (cumulative distribution function)
Supplemental Figure 3: Consensus matrix with 3 clusters
A consensus matrix heatmap is shown for a randomly selected imputed dataset clustered into 3 phenotypes. The heatmap allows visualization of consensus cluster assignments to evaluate cluster stability. Darker shades of green indicate higher stability.
Supplemental Figure 4: Consensus matrix with 4 clusters
A consensus matrix heatmap is shown for a randomly selected imputed dataset clustered into 4 phenotypes. The heatmap allows visualization of consensus cluster assignments to evaluate cluster stability. Darker shades of green indicate higher stability. The choice of 4 clusters shows less stability than 3 clusters (see Supplemental Figure 3).
Supplemental Figure 5: Cumulative Proportion of Variance Explained
The proportion of variance explained by each principal component is summed over all principal components. For example, PC1 and PC2 cumulatively explain 20% of the variation in the dataset.
Abbreviations: PC1 (principal component 1); PC2 (principal component 2)
Supplemental Figure 6: Comorbidities by Phenotype
Chord diagram illustrates the prevalence of comorbidities (% observed) for the three clinical phenotypes.
Supplemental Figure 7: Complications by Phenotype
Chord diagram illustrates the prevalence of complications (% observed) for the three clinical phenotypes.
Supplemental Figure 8: Clinical Outcomes by Phenotype
Chord diagram illustrates the prevalence of clinical outcomes (% observed) for the three clinical phenotypes.
Abbreviations: ICU (intensive care unit); Vent (mechanical ventilation); Readmit (readmission to hospital or ICU); ECMO (extracorporeal membrane oxygenation).
Supplemental Table 1: Home medications and hospital day 5 laboratory values of hospitalized COVID-19 patients with clinical phenotypes I, II, and III.
* Categorical variables presented as count (%), continuous variables presented as median (interquartile range).
† Loop Diuretics include furosemide, torsemide, budesonide.
Abbreviations: ALC, Absolute Lymphocyte Count; ANC, Absolute Neutrophil Count; ARB, Angiotensin receptor blocker; AST, Aspartate transaminase; CRP, C-Reactive Protein; ESR, Erythrocyte Sedimentation Rate; EKG, Electrocardiogram; IL, Interleukin; LDH, Lactate Dehydrogenase; Pro-BNP, proB-type Natriuretic Peptide; PLT, Platelets; PTT, Partial Thromboplastin Time; TNF-alpha, Tumor Necrosis Factor-alpha; WBC, White Blood Cell;
Supplemental Table 2: Association of Clinical Phenotype with In-Hospital Complications
Abbreviations: OR, odds ratio; CI, confidence interval; LR, likelihood ratio Legend: Reference group for all models is Phenotype III. All models adjusted for sex, race/ethnicity, and Elixhauser Comorbidity Index.
Supplemental Table 3: Comorbidities based on ICD-10 codes
Supplementary Table 3: List of ICD 10 codes that were used to classify diagnosis.
Abbreviations: T1DM: Type 1 diabetes mellitus, T2DM Type 2 diabetes mellitus, HFpEF: heart failure with preserved ejection fraction, HF: heart failure, CAD: coronary artery disease, NAFLD: non-alcoholic fatty liver disease, NASH: non-alcoholic steatohepatitis, VTE: venous thromboembolism, HIT: heparin induced thrombocytopenia, DIC: disseminated intravascular coagulation, ITP: idiopathic thrombocytopenia, HLH: hemophagocytic lymphohistiocytosis, MI: myocardial infarction, COPD: chronic obstructive lung disease, ILD: interstitial lung disease, AICD: Automatic Implantable Cardioverter Defibrillator, VAD: ventricular assist device, CKD: chronic kidney disease, ESRD: end stage renal disease, Afib: atrial fibrillation, HIV: human immunodeficiency virus, Flu: influenza, IBD: inflammatory bowel disease.
Supplemental Table 4: Classification of complications by Group
Supplementary Table 4: Complications were grouped into different organ/system-based classifications.
Acknowledgements
The authors thank Eric Murray and Fairview IT for collection of data
Footnotes
NEI ingra107{at}umn.edu
TS schac008{at}umn.edu
JGC chipm001{at}umn.edu
BB benso040{at}umn.edu
CTJ ctignane{at}umn.edu
MIL lupei001{at}umn.edu
GMM gmelton{at}umn.edu
DW wack0012{at}umn.edu
RR rreilkof{at}umn.edu
SE saharle{at}umn.edu
BT thie0149{at}umn.edu
ERL lusc0006{at}umn.edu
JP prope012{at}umn.edu
AC anthchar{at}med.unc.edu
LS siege245{at}umn.edu
MHS syed0021{at}umn.edu
JR rhei0005{at}umn.edu
MGU mgusher{at}umn.edu
SI ikram001{at}umn.edu
MA abass004{at}umn.edu