Characterizing COVID-19 Clinical Phenotypes and Associated Comorbidities and Complication Profiles ================================================================================================== * Elizabeth R. Lusczek * Nicholas E. Ingraham * Basil S. Karam * Jennifer Proper * Lianne Siegel * Erika S. Helgeson * Sahar Lotfi-Emran * Emily J. Zolfaghari * Emma Jones * Michael G. Usher * Jeffrey G. Chipman * R. Adams Dudley * Bradley Benson * Genevieve B. Melton * Anthony Charles * Monica I. Lupei * Christopher J. Tignanelli ## 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). ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/14/2020.09.12.20193391/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2020/09/14/2020.09.12.20193391/F1) Figure 1: Score Plot: PC2 vs. PC1 The principal component scores for PC1 and PC2 are plotted. Each point represents a patient in the dataset. Colors represent the cluster (phenotype) that the patient was assigned to by consensus clustering. Ellipses around each cluster/phenotype specify 95% confidence intervals, assuming a bivariate normal distribution. Abbreviations: PC1 (principal component 1); PC2 (principal component 2) ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/14/2020.09.12.20193391/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2020/09/14/2020.09.12.20193391/F2) Figure 2: Contribution of Variables to PC1 The contributions of each of the 33 variables used in the clustering to principal component 1 are shown. The red line marks the expected average contribution of each variable if the contributions of the variables were uniform across the dataset. Variables contributing most to the observed pattern in PC1 are D-dimer and albumin. Abbreviations: PC1 (principal component 1); Abs\_Nphil_Ct (absolute neutrophil count); LDH (lactate dehydrogenase); CRP (C-reactive protein); WBC (white blood cell count); HCT (hematocrit); HGB (hemoglobin); Tbili (total bilirubin); RDW (red cell distribution width); AST (aspartate aminotransferase); Alk_phos (alkaline phosphatase); RR (respiratory rate); CA (calcium); TP (total protein); INR (internal normalized ratio of prothrombin time); CO2 (carbon dioxide); K (potassium); O2SAT (oxygen saturation); BMI (body mass index); PLT (platelet); PP (pulse pressure); Na (sodium); SBP (systolic blood pressure); Abs_mono_ct (absolute monocyte count); MCV (mean corpuscular volume) ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/14/2020.09.12.20193391/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2020/09/14/2020.09.12.20193391/F3) Figure 3: PCA Biplot: PC2 vs. PC1 The scores (points) and loadings (arrows) of PC1 and PC2 are plotted for each patient and variable in the model. 95% confidence ellipses for the scores are shown. The biplot facilitates interpretation of the scores and loadings, assigning context to the variables which prominently contribute to the phenotypes. Abbreviations: PC1 (principal component 1); PC2 (principal component 2); PCA (principal component analysis); Abs\_Nphil_Ct (absolute neutrophil count); LDH (lactate dehydrogenase); CRP (C-reactive protein); WBC (white blood cell count); HCT (hematocrit); HGB (hemoglobin); Tbili (total bilirubin); RDW (red cell distribution width); AST (aspartate aminotransferase); Alk_phos (alkaline phosphatase); RR (respiratory rate); CA (calcium); TP (total protein); INR (internal normalized ratio of prothrombin time); CO2 (carbon dioxide); K (potassium); O2SAT (oxygen saturation); BMI (body mass index); PLT (platelet); PP (pulse pressure); Na (sodium); SBP (systolic blood pressure); Abs_mono_ct (absolute monocyte count); MCV (mean corpuscular volume) ### 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**). View this table: [Table 1:](http://medrxiv.org/content/early/2020/09/14/2020.09.12.20193391/T1) Table 1: Baseline demographics, comorbidities, and clinical outcomes of hospitalized COVID-19 patients with clinical phenotypes I, II, and III * Categorical variables presented as count (%), continuous variables presented as median (interquartile range) unless otherwise specified. Continuous variables were evaluated with Kruskal-Wallis tests. Categorical variables were evaluated with chi-square tests or Fisher’s exact tests (counts < 5). Abbreviations: ADI, area deprivation index; BMI, body mass index; ECMO, Extracorporeal membrane oxygenation; ICU, Intensive Care Unit 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**). ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/14/2020.09.12.20193391/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2020/09/14/2020.09.12.20193391/F4) Figure 4: Relative Risk Ratio of Comorbidities to Clinical Phenotypes Relative Risk ratios of comorbidities of phenotypes I and II compared to the reference group phenotype III. ### 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. View this table: [Table 2:](http://medrxiv.org/content/early/2020/09/14/2020.09.12.20193391/T2) Table 2: Association of Clinical Phenotype with Clinical Outcomes Abbreviations: PH, proportional hazards; HR, hazard ratio; CI, confidence interval; OR, odds ratio; ICU, intensive care unit; IRR, incidence rate ratio; LOS, length of stay; LR, likelihood ratio Legend: Reference group for all models is Phenotype III. All models adjusted for sex, race/ethnicity, and Elixhauser Comorbidity Index. * LOS model only included patients that survived. ## 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 1. NIH NHLBI T32HL07741 (NEI) 2. 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. 3. 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) 1. NIH NHLBI T32HL07741 (NEI) 2. 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. 3. 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 * Received September 12, 2020. * Revision received September 12, 2020. * Accepted September 14, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission. ## References 1. 1.Wu Z, McGoogan JM. Characteristics of and Important Lessons From the Coronavirus Disease 2019 (COVID-19) Outbreak in China: Summary of a Report of 72314 Cases From the Chinese Center for Disease Control and Prevention. JAMA. 2020. 2. 2.Ingraham NE, Lotfi-Emran S, Thielen BK, et al. Immunomodulation in COVID-19. Lancet Respir Med. 2020;8(6):544–546. 3. 3.Ingraham NE, Tignanelli CJ. Fact Versus Science Fiction: Fighting Coronavirus Disease 2019 Requires the Wisdom to Know the Difference. Crit Care Explor. 2020;2(4):e0108. 4. 4.Group RC, Horby P, Lim WS, et al. Dexamethasone in Hospitalized Patients with Covid-19 - Preliminary Report. The New England journal of medicine. 2020. 5. 5.Wiersinga WJ, Rhodes A, Cheng AC, Peacock SJ, Prescott HC. Pathophysiology, Transmission, Diagnosis, and Treatment of Coronavirus Disease 2019 (COVID-19): A Review. JAMA. 2020. 6. 6.Li X, Ma X. Acute respiratory failure in COVID-19: is it “typical” ARDS? Crit Care. 2020;24(1): 198. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13054-020-02911-9&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F14%2F2020.09.12.20193391.atom) 7. 7.Diehl JL, Peron N, Chocron R, et al. Respiratory mechanics and gas exchanges in the early course of COVID-19 ARDS: a hypothesis-generating study. Ann Intensive Care. 2020;10(1):95. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13613-020-00716-1&link_type=DOI) 8. 8.Arentz M, Yim E, Klaff L, et al. Characteristics and outcomes of 21 critically ill patients with COVID-19 in Washington State. Jama. 2020;323(16):1612–1614. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2020.4326&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F14%2F2020.09.12.20193391.atom) 9. 9.Qin C, Zhou L, Hu Z, et al. Dysregulation of immune response in patients with COVID-19 in Wuhan, China. Clinical Infectious Diseases. 2020. 10. 10.Wu C, Chen X, Cai Y, et al. Risk factors associated with acute respiratory distress syndrome and death in patients with coronavirus disease 2019 pneumonia in Wuhan, China. JAMA internal medicine. 2020. 11. 11.Huang C, Wang Y, Li X, et al. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. The lancet. 2020;395(10223):497–506. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0146736(20)30183-5&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F14%2F2020.09.12.20193391.atom) 12. 12.Little RJA, Rubin DB. Statistical analysis with missing data. New York: Wiley; 1987. xiv, 278. 13. 13. Stef van Buuren, Karin Groothuis-Oudshoorn (2011). mice: Multivariate Imputation by Chained Equations in R. Journal of Statistical Software, 45(3), 1-67. URL [https://www.jstatsoft.org/v45/i03/](https://www.jstatsoft.org/v45/i03/). 14. 14. Derek Chiu and Aline Talhouk (2020). diceR: Diverse Cluster Ensemble in R. R package version 1.0.0. [https://CRAN.R-project.org/package=diceR](https://CRAN.R-project.org/package=diceR). 15. 15.Nassiri V, Lovik A, Molenberghs G, Verbeke G. On using multiple imputation for exploratory factor analysis of incomplete data. Behav Res Methods. 2018;50(2):501–517. 16. 16.Bramante C, Ingraham N, Murray T, et al. Observational Study of Metformin and Risk of Mortality in Patients Hospitalized with Covid-19. *medRxiv*. 2020:2020.2006.2019.20135095. 17. 17.Jin JM, Bai P, He W, et al. Gender Differences in Patients With COVID-19: Focus on Severity and Mortality. Front Public Health. 2020;8:152. 18. 18.Ingraham NE, Purcell LN, Karam BS, et al. Racial/Ethnic Disparities in Hospital Admissions from COVID-19 and Determining the Impact of Neighborhood Deprivation and Primary Language. *medRxiv*. 2020:2020.2009.2002.20185983. 19. 19.Elixhauser A, Steiner C, Harris DR, Coffey RM. Comorbidity measures for use with administrative data. Medical care. 1998:8-27. 20. 20.Gu Z, Gu L, Eils R, Schlesner M, Brors B. circlize implements and enhances circular visualization in R. Bioinformatics. 2014;30(19):2811–2812. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btu393&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24930139&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F14%2F2020.09.12.20193391.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000343082900017&link_type=ISI) 21. 21.R Core Team (2020). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. URL [https://www.R-project.org/](https://www.R-project.org/). 22. 22.Petrilli CM, Jones SA, Yang J, et al. Factors associated with hospital admission and critical illness among 5279 people with coronavirus disease 2019 in New York City: prospective cohort study. BMJ. 2020;369:m1966. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjE4OiIzNjkvbWF5MjJfMTUvbTE5NjYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wOS8xNC8yMDIwLjA5LjEyLjIwMTkzMzkxLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 23. 23.Zhou F, Yu T, Du R, et al. Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study. Lancet. 2020;395(10229):1054–1062. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(20)30566-3&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F14%2F2020.09.12.20193391.atom) 24. 24.Ye W, Chen G, Li X, et al. Dynamic changes of D-dimer and neutrophil-lymphocyte count ratio as prognostic biomarkers in COVID-19. Respir Res. 2020;21(1):169. 25. 25.Tian W, Jiang W, Yao J, et al. Predictors of mortality in hospitalized COVID-19 patients: A systematic review and meta-analysis. J Med Virol. 2020. 26. 26.Shi S, Qin M, Shen B, et al. Association of Cardiac Injury With Mortality in Hospitalized Patients With COVID-19 in Wuhan, China. JAMA Cardiol. 2020. 27. 27.Yang X, Yu Y, Xu J, et al. Clinical course and outcomes of critically ill patients with SARS-CoV-2 pneumonia in Wuhan, China: a single-centered, retrospective, observational study. Lancet Respir Med. 2020;8(5):475–481. 28. 28.Ingraham NE, Barakat AG, Reilkoff R, et al. Understanding the Renin-Angiotensin-Aldosterone-SARS-CoV-Axis: A Comprehensive Review. Eur Respir J. 2020:2000912. 29. 29.Tignanelli CJ, Ingraham NE, Sparks MA, et al. Antihypertensive drugs and risk of COVID-19? Lancet Respir Med. 2020;8(5):e30-e31. 30. 30.Kuo CP. L. Atkins JC, et al. COVID-19 severity is predicted by earlier evidence of accelerated aging. *Medrxiv*. 2020. 31. 31.Azoulay E, Zafrani L, Mirouse A, Lengline E, Darmon M, Chevret S. Clinical phenotypes of critically ill COVID-19 patients. Intensive Care Med. 2020;46(8):1651–1652. 32. 32.Yang M, Xie L, Liu X, Hao Q, Jiang J, Dong B. The gamma gap predicts 4-year all-cause mortality among nonagenarians and centenarians. Sci Rep. 2018;8(1): 1046.