key: cord-0764597-es3xa9k0 authors: Bahl, Amit; Van Baalen, Morgan Nees; Ortiz, Laura; Chen, Nai-Wei; Todd, Courtney; Milad, Merit; Yang, Alex; Tang, Jonathan; Nygren, Madalyn; Qu, Lihua title: Early predictors of in-hospital mortality in patients with COVID-19 in a large American cohort date: 2020-09-24 journal: Intern Emerg Med DOI: 10.1007/s11739-020-02509-7 sha: aa839017fe18f4ac16b11562b99b53c17f898a94 doc_id: 764597 cord_uid: es3xa9k0 Coronavirus disease (COVID-19) has aggressively spread across the United States with numerous fatalities. Risk factors for mortality are poorly described. This was a multicentered cohort study identifying patient characteristics and diagnostic markers present on initial evaluation associated with mortality in hospitalized COVID-19 patients. Epidemiological, demographic, clinical, and laboratory characteristics of survivors and non-survivors were obtained from electronic medical records and a multivariable survival regression analysis was conducted to identify risk factors of in-hospital death. Of 1629 consecutive hospitalized adult patients with confirmed COVID-19 from March 1st thru March 31, 2020, 1461 patients were included in final analysis. 327 patients died during hospitalization and 1134 survived to discharge. Median age was 62 years (IQR 50.0, 74.0) with 56% of hospitalized patients under the age of 65. 47% were female and 63% identified as African American. Most patients (55%) had either no or one comorbidity. In multivariable analysis, older age, admission respiratory status including elevated respiratory rate and oxygen saturation ≤ 88%, and initial laboratory derangements of creatinine > 1.33 mg/dL, alanine aminotransferase > 40 U/L, procalcitonin > 0.5 ng/mL, and lactic acid ≥ 2 mmol/L increased risk of in-hospital death. This study is one of the largest analyses in an epicenter for the COVID-19 pandemic. Older age, low oxygen saturation and elevated respiratory rate on admission, and initial lab derangements including renal and hepatic dysfunction and elevated procalcitonin and lactic acid are risk factors for in-hospital death. These factors can help clinicians prognosticate and should be considered in management strategies. In December 2019, the novel coronavirus (severe acute respiratory syndrome coronavirus 2, SARS-CoV-2) was first detected in Wuhan, China. Since then, coronavirus disease 2019 has spread rapidly to nearly every continent around the globe [1] [2] [3] . After the first case of this illness was discovered in the United States in January 2020 in Seattle, WA, it spread to every state in under 2 months [4, 5] . The United States quickly became the world leader in number of cases and deaths. As of April 30, 2020, over 3.1 million cases and 227,000 fatalities were identified worldwide, of which over 60,000 fatalities were in the United States [6] . COVID-19 patients can present with a range of symptoms, varying from asymptomatic to severe cases with acute respiratory failure requiring hospitalization [7] [8] [9] . Presently, patients with more dire illness are apparently older and have complicating comorbidities such as diabetes, hypertension, and obesity [10, 11] . There is currently no cure for COVID-19. Hospitalization results in supportive treatment with oxygen support and a myriad of medications to prevent multiple organ dysfunction and death [12] . Despite the increasing number of cases, there is limited data regarding the clinical characteristics of patients admitted to the hospital, outside of China, specifically in a diverse American population. Furthermore, while it is evident that this disease is highly contagious and virulent with local community and person-to-person transmission [13, 14] , there is limited evidence that describes the predictors, particularly early in patients' hospital course, of poor outcomes and mortality for hospitalized patients. An improved understanding is required to guide decisions regarding admissions, hospital capacity and therapeutic and operational resources. As COVID-19 has spread, several regions in the United States have been disproportionately impacted. It is our objective to determine early risk factors for in-hospital death for COVID-19 patients in the largest health system in southeastern Michigan. This was an analysis of a multicentered cohort study assessing early predictors of mortality in hospitalized patients with COVID-19 comparing survivors to non-survivors. The study was conducted at Beaumont Health, an eight-hospital acute care regional health system caring for 2.2 million people across the communities within the Metro Detroit catchment area. The hospitals range from a large tertiary care academic center to intermediate-sized and smaller community hospitals. As the epidemic evolved, the health system converted one of the smaller hospitals to a complete COVID-19 center, including converting the emergency department into additional intensive care beds. As the surge in COVID-19 volume developed across the region, hospital systems collaborated to optimize the transfer process and accommodate capacity constraints. While many patients were transferred within Beaumont Health, many patients were also transferred to other hospitals. The study was approved by the Institutional Review Board at the home institution. Written informed consent requirement was waived due to the rapid emergence of COVID-19. Data were analyzed and interpreted by the authors. Data were obtained from the integrated electronic health record (EHR; Epic Systems, Verona, WI). Patients over 18 years of age who were admitted with COVID-19 from March 1st through March 31st 2020 were included. Data were collected until April 23, 2020. All patients had a laboratory confirmed case of COVID-19 as defined by a positive result on a reverse-transcriptase-polymerasechain-reaction (RT-PCR) test of nasopharyngeal swab. Exclusions consisted of patients who left the hospital against medical advice, transfers to external hospitals, or if hospital course was ongoing beyond April 23, 2020. Transfers within Beaumont Health were included as investigators had full access to medical records. Epidemiological, demographic, radiological, therapeutic, clinical, and outcomes data were extracted. The integrity of the data was verified by two attending emergency medicine physicians. Hospital admission was based on the clinical judgment of the treating emergency medicine provider. Laboratory and radiological testing was conducted at the discretion of the treating physicians. After initial COVID laboratory testing, patients were not routinely serially tested to evaluate for clearance of acute infection due to paucity of testing supplies. Discharge disposition post-hospitalization was based on patients' clinical condition. Patients were either discharged to home, skilled nursing facility or rehabilitation, hospice, or expired in the hospital. Admission data included demographics such as age, race, and gender. Clinical data included symptoms prompting presentation to the ED, time of onset (days) of symptoms prior to ED visit, comorbidities, body mass index (BMI), tobacco history, vital signs including lowest oxygen level on room air within first 24 h of admission and maximum temperature, laboratory findings, and admission floor type. Radiological data included chest x-ray and chest computed tomography (CT) impressions. Common laboratory analyses included complete blood count (CBC) with absolute lymphocyte count, metabolic chemistry panel (BMP/CMP), lactate dehydrogenase, lactic acid, procalcitonin, troponin, ferritin, and D-dimer. Laboratory and radiological data depicted the first test result occurring within the first 24 h of presentation to the ED. Hospital treatment data included use of investigational or adjunctive therapies, renal replacement therapy, extracorporeal membrane oxygenation, oxygen and ventilation therapy, and intensive care unit (ICU) admission. Fever was defined as an axillary temperature of ≥ 38.0 °C. Acute respiratory failure was defined as a partial pressure of oxygen (PaO 2 ) < 60 mm Hg on blood gas assessment. Pulse oximetry was used as a surrogate for PaO 2 , with ≤ 88% on room air meeting the criteria. Acute kidney injury was defined per KDIGO clinical practice guidelines [15] . Acute cardiac injury was defined as an elevated cardiac biomarker (troponin I) > 99% of the upper reference limit [16] . Septic shock was defined as a patient with diagnosed sepsis requiring vasopressor therapy to maintain mean arterial blood pressure > 65 mmHg. Patients with a quick Sequential Organ Failure Assessment score of two or more were considered to have sepsis [17] . Extubation failure was defined as need for reintubation within 72 h [18] . Descriptive analyses were used to summarize epidemiological characteristics and clinical findings stratified by in-hospital survival status. Continuous and categorical variables were expressed as medians (interquartile ranges; IQR) and frequencies (percentages), respectively. Kruskal-Wallis test (continuous variables) and Chi-squared or Fisher's exact test (categorical variables) were used to compare differences between non-survivors and survivors. To explore the association between risk factors and in-hospital mortality, discrete-time survival regression was employed for univariable and multivariable analyses. Missing data were imputed by the procedure of multiple imputation using PROC MI in SAS. To avoid collinearity and overfitting on multivariable model, the L1-penalized least absolute shrinkage and selection operator (LASSO) regression selection was employed to the analog with a logistic regression that fits the discrete-time survival, using a statistical package "glmnet" in R. Variables selected via LASSO regression and clinical representative variables in terms of the basis of relevant studies were applied to multivariable analysis. The effects of risk factors were combined from 20 imputed datasets, accounting for the additional variability introduced by the multiple imputation, through PROC MIANALYZE in SAS. The corresponding c-statistics and leave-one-out validation were further used to evaluate the performance of modeling in multivariable regression. To facilitate the assessment of risk on mortality from the multivariable model, a prediction risk score using categorical variables was built. The assigned score for each risk factor was determined according to estimates of regression coefficients multiplied by 10. Patients were classified into risk groups based on their risk score and the mortality rates of each risk group was calculated. The differences in survival among risk groups over study period were also assessed using Kaplan-Meier method with the log-rank test. A bootstrap cross-validation with 1000 bootstrap samples per imputed dataset was further employed to Hospice Discharge validate the performance of scoring predictivity by assessing the similarity of mortality rates on risk groups. All tests of statistical significance were two-sided with p value < 0.05 indicating a significant difference. Analyses were conducted using R-3.6.2 (R Foundation for Statistical Computing) and SAS v-9.4 (SAS Institute, Inc., Cary, NC). Respiratory failure was the most common patient outcome complicating 52.5% of admissions and present in 93.9% of non-survivors. High-flow oxygen therapy via nasal cannula was used in 28.7% of patients, more commonly employed in non-survivors (63.0%) compared to survivors (18.8%). Invasive mechanical ventilation was employed in 308 patients (21.1%). It was used in 67.0% of non-survivors compared to 7.9% of survivors (p < 0.001). Overall, 71.1% of patients did not survive to hospital discharge when invasive mechanical ventilation was utilized. Acute kidney injury was noted in 382 patients and developed in 71.3% of nonsurvivors compared to 13.1% of survivors (p < 0.001). Of these patients, 20.9% required hemodialysis. Acute cardiac injury was a less frequent complication occurring in 9.1% of cases but was a more common complication in non-survivors (26.0%). Septic shock requiring vasopressors occurred in 16.2% of patients, occurring in 59.9% of non-survivors compared to 3.6% of survivors (p < 0.001). Of the 237 patients classified as septic shock, only 30 patients (12.65%) had positive blood cultures. 25 patients had blood cultures that identified a bacterial source and six patients had blood cultures that grew Candida albicans, treated with an antifungal. 7 patients had more than one bacteria isolated, and one patient had bacteria and Candida isolated. 207 patients had The median length of illness (symptom onset) for survivors was 12.0 days compared to 15.0 days for non-survivors. The median hospital stay was 7 days. Intensive care admission was required in 25.6% of patients with a median ICU length of stay of 7.0 days ( Table 2) . The most commonly employed investigational treatment was hydroxychloroquine, used in 73.4% of patients with the average time from illness onset to therapy start of 7.0 days (IQR 5.0-9.0). Systemic glucocorticoids were utilized in 42.2% of patients with median time from illness onset to therapy start of 9.0 days (IQR 6.0-13.0). In univariable analysis, epidemiological and clinical characteristics listed on Table 1 were utilized to illustrate the unadjusted effects on the hazard of in-hospital patient death. In multivariable survival model 1, older age (HR 1.05, 95% CI 1.04 to 1.06), admission elevated respiratory rate (HR 1.06, 95% CI 1.04-1.09), pulse oximetry oxygen saturation less than or equal to 88% (HR 1.82, 95% CI 1.33-2.50) and initial laboratories: creatinine > 1.33 mg/dL (HR 1.41, 95% CI 1.08-1.85), ALT > 40 U/L (HR 1.36, 95% CI 1.05-1.77), procalcitonin > 0.5 ng/mL (HR 2.14, 95% CI 1.37-3.36), and lactic acid ≥ 2.0 mmol/L (HR 1.57, 95% CI 1.17 to 2.12) were predictors of in-hospital death. Similar results were found when using age groups as a risk factor in multivariable survival model 2 (Table 3) . To facilitate the assessment of risk on mortality, a multivariable survival model using categorical variables for predicting mortality is shown in Table 4 with a scoring system depicted in Table 5 . According to our model, a score of ≤ 12, 13-26, and ≥ 27 carries a mortality risk of 3.91%, 17.98%, and 51.50%, respectively, with the significant differences among risk-scoring groups (p < 0.001, Fig. 3 ). In this study, we described the characteristics of 1461 patients requiring hospitalization for laboratory-confirmed COVID-19 in Metro Detroit. We also delineated the average clinical course and identified risk factors for mortality that were present early in the hospitalization course. Overall, mortality exceeded 22% for hospitalized patients and is consistent with other published reports, ranging between 14 and 28% [14, [19] [20] [21] . Specifically older age, abnormal respiratory parameters on presentation, and initial laboratory results demonstrating renal insufficiency, liver injury, and elevated procalcitonin and lactic acid levels were risk factors for in-hospital death. These factors can help clinicians prognosticate so a scoring system was developed based on these variables among others to help predict inpatient mortality. It was found that a score of less than 12 was low risk of mortality, while a score ≥ 27 represented greater than 50% likelihood of death. This scoring system is based on the clinical features in our cohort and could help physicians estimate mortality and possibly aid in resource allocation, although has not been validated outside of this cohort. In this analysis, we identified a strong association of demographic and epidemiologic characteristics with hospitalization and mortality. Age was a significant predictor of death as median age of fatalities was 73 compared to 59 years for survivors. This is consistent with other reports citing that the geriatric population is at the highest risk in this pandemic [14] . Interestingly, though, while mortality predominantly occurs in the elderly population, younger adults remain at high risk for hospitalization as nearly 56% of patients in our study were between 18 and 64 years. Gender vulnerability of males has also been a commonly reported descriptive within investigations, with male gender comprising 73% of deaths and 82% of ICU admissions in Chinese and Italian publications, respectively [20, 21] . It is postulated that women may be partially protected due to genetic, biological, and hormonal factors [22] [23] [24] . We, however, did not find a male predominance to the COVID-19 illness within our large sample size. While males made ALT alanine aminotransferase, AST aspartate aminotransferase, HR hazard ratio, CI confidence interval a The model assessment indicated the predictive ability (c statistics) on in-hospital survival status was as high as 0.78 and the leave-one-out cross-validated results showed that there was no overfitting in the multivariable model up slightly more hospital admissions (53%) than females, there was no mortality difference between genders. Race was another relevant variable. In our study, we found that African Americans were over twice as likely to be hospitalized due to COVID-19 compared to Caucasians. In a quality analysis of demographics at Beaumont Health, African Americans accounted for only 22% of hospitalizations in 2019, yet they accounted for 63% of COVID-19 hospitalizations. The African American community had a higher prevalence of hospitalizations but not mortality, leading us to speculate that other factors such as inability to socially distance due to multi-generational households, the need to work outside the home, and other similar characteristics, put this population at greater risk to be exposed to the illness. Pre-existing conditions are described as significant predictors of COVID-19 disease outcomes in the literature [14, 19, 20, 25, 26] . Specifically, coronary artery disease and hypertension have been identified as risk factors for mortality. In our study, 34% of hospitalized patients had no comorbidities and 55% had one or no pre-existing medical conditions. This prevalence demonstrates the virulence of COVID-19 in a seemingly healthier population. Unlike other publications, our multivariable analysis did not provide supportive evidence for an association between specific comorbidities and mortality. Additionally, obesity has reportedly been associated with poor outcomes in COVID-19. In a large recent analysis based in New York, BMI greater than 30 was associated with critical illness and a threshold over 40 was implicated as a variable associated with hospitalization [14] . Our cohort had a median BMI of 31; however, while obesity is considered a pro-inflammatory state, in a multivariable analysis, this did not present a significant positive association with mortality. There are several clinical variables that increase the risk of in-hospital death. Laboratory analysis revealed that elevated lactic acid and procalcitonin levels were risk factors for death. These tests have been described in other investigations as predictors of mortality in COVID-19 patients likely due to relationship of these markers to severe sepsis and secondary bacterial infections [27] [28] [29] . Our study also identified renal insufficiency on initial laboratories as a prognosticator for severe illness and mortality. This has also been documented in Wuhan, where the presence of kidney disease was associated with in-hospital mortality. It has been postulated this is multifactorial but could be due to the novel coronavirus using angiotensin converting enzyme 2 (ACE2) receptor for renal cell entry [30, 31] . It is also thought that virus-induced cytokines could cause indirect effects on renal cells via shock or hypoxia [30, 32] . Not surprisingly, there was a strong association between hypoxia with both hospitalization and death. The hospitalized patient was hypoxic at baseline, with median admission pulse oximetry of 92%. Additionally, admission pulse oximetry ≤ 88% was predictive of demise. The use of high flow oxygenation and mechanical ventilation were common Table 4 Multivariable survival model using categorical variables for predicting mortality ALT alanine aminotransferase, HR hazard ratio, SE standard error a Estimates of regression coefficient (β) were pooled from survival models fitted on 20 imputed datasets. Points associated with the presence of a given level of a risk factor (variable) were assigned by multiplying β by 10 and rounding to the nearest integer. The reference level was assigned zero point. The theoretical ranges of total score were − 8 poor outcomes with a 71% case fatality rate in this cohort. In other smaller case series, reported mortality associated with invasive ventilation was also dismal with 97% mortality in one cohort of 32 patients and 100% in a 17-patient series [19, 20] . The utilization of invasive mechanical ventilation represents a complicated provider decision based on oxygen saturation, dyspnea, respiratory rate, chest x-ray, and other factors [33] . In our cohort this decision was also likely influenced by factors including guidelines supporting early intubation and cautioning the use of noninvasive mechanical ventilation strategies [34] . Noninvasive strategies such as bilevel positive airway pressure were used in less than 10% of cases due to the potential risk of viral aerosolization [35] . Other countries used non-invasive ventilation treatments more liberally and some research suggests a possible benefit [36] . In a Chinese study of 113 fatalities, only 6% of patients received invasive treatment compared to 37% receiving noninvasive ventilation [20] . It is unclear whether this strategy was employed purposefully or due to limited resources, and the impact of this choice on mortality is unclear. Given the high mortality rates associated with invasive mechanical ventilation, the overall Table 5 Mortality assessment using the predictor scoring system a The observed ranges of total score were − 4 and 58 on 20 imputed datasets. The median score was 19 and the 25th and 75th percentiles (IQR) were 12 and 26, respectively. The 10th and 90th percentiles were 7 and 33, respectively. Mortality rates and the corresponding estimates of hazard ratios were pooled from 20 imputed datasets b Using bootstrap cross-validation with 1000 bootstrap samples per imputed dataset, estimates of mortality rates from low risk to high risk were 3.92% (95% CI 2.03-6.05%), 17 . 3 Kaplan-Meier survival curve for mortality for three risk groups. Figure shows overall survival of study patients associated with three risk groups (score: ≤ 12, 13 to 26, ≥ 27) during study period. The estimated survival cures were pooled from 20 imputed datasets strategy for oxygenation and ventilation deserves additional investigation with full consideration of all modes of oxygen delivery, both invasive and noninvasive ventilation, delayed intubation with permissive hypoxia, positioning maneuvers, among other therapies. The study has several limitations. Patients are from a single geographic region, treated within a single health system and factors that are associated with hospitalization and oxygen requirements could differ in other regions, although the region is relatively diverse. While the health system services a range of urban to suburban populations, the rural community is relatively absent from this cohort. Second, ED provider discretion regarding need for admission created potential variability in admission characteristics of patients. Third, data from the EHR were occasionally incomplete. Finally, the case fatality rate likely underrepresents true mortality as transfers out of the health system often needed intensive care resources and ICU admission [19, 37] . Further, we quantified in-hospital mortality through the designated follow-up date of April 23rd. We recently queried clinical status for all patients and found that an additional 17 patients who were discharged home or to a skilled nursing facility had also expired. This study represents one of the largest cohort analyses in a region considered an epicenter for the COVID-19 pandemic. Older age, low admission saturation ≤ 88%, elevated respiratory rate, and initial laboratory derangements including renal and hepatic injury, and elevated procalcitonin and lactic acid levels were the most significant risk factors for in-hospital death. These factors, in addition to our scoring system, can help clinicians prognosticate and should be considered in management strategies. Conflict of interest This manuscript in part or in full has not been submitted or published anywhere. No authors have any relevant conflict of interest disclosures. The study was approved by the Institutional Review Board at the home institution. Informed consent Written informed consent requirement was waived due to the rapid emergence of COVID-19. China Novel Coronavirus investigating and research team. A novel coronavirus from patients with pneumonia in China COVID-19: diagnosis, management and prognosis: a new topical collection of internal and emergency medicine World Health Organization (2020) Novel coronavirus-China An interactive web-based dashboard to track COVID-19 in real time First case of 2019 novel coronavirus in the United States Coronavirus Disease (COVID-19)-statistics and research. Our World in Data 2020 Clinical features of patients infected with 2019 novel coronavirus in Wuhan Clinical characteristics of 138 hospitalized patients with 2019 novel coronavirus-infected pneumonia in Wuhan Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: a descriptive study COVID-19 with different severity: a multi-center study of clinical features A rapid advice guideline for the diagnosis and treatment of 2019 novel coronavirus (2019-nCoV) infected pneumonia (standard version) Early transmission dynamics in Wuhan, China, of novel coronavirus-infected pneumonia A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating person-to-person transmission: a study of a family cluster Factors associated with hospitalization and critical illness among 4,103 patients with COVID-19 disease in New York City KDIGO clinical practice guidelines for acute kidney injury Third universal definition of myocardial infarction Changing definitions of sepsis Cough peak flows and extubation outcomes Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study Clinical characteristics of 113 deceased patients with coronavirus disease 2019: a retrospective study Baseline characteristics and outcomes of 1591 patients infected with SARS-CoV-2 admitted to ICUs of the Lombardy region SARS-CoV-2 Cell Entry Depends on ACE2 and TMPRSS2 and Is Blocked by a Clinically Proven Protease Inhibitor Emerging markers in cardiovascular disease: where does angiotension-converting enzyme 2 fit in? Sex-based differences in susceptibility to SARS-CoV infection Preliminary estimates of the prevalence of selected underlying health conditions among patients with Coronavirus Disease 2019-United States Instituto SuperioreDiSanita Clinical characteristics and outcomes of 112 cardiovascular patients infected by 2019-nCoV Risk factors of critical & mortal COVID-19 cases: A systematic literature review and metaanalysis Serial Procalcitonin Predicts Mortality in Severe Sepsis Patients: Results From the Multicenter Procalcitonin MOnitoring SEpsis (MOSES) Study Kidney disease is associated with in-hospital death of patients with COVID-19 The SARS-CoV-2 receptor, ACE-2, is expressed on many different cell types: implications for ACE-inhibitor-and angiotensin II receptor blocker-based cardiovascular therapies Covid-19 and kidney failure in the acute care setting: our experience from Seattle Basing respiratory management of coronavirus on physiological principles Interim guidance: Clinical management of severe acute respiratory infection (SARI) when COVID-19 disease is suspected. Geneva: World Health Organization Why did outbreaks of severe acute respiratory syndrome occur in some hospital wards but not in others? Flowchart for non-invasive ventilation support in COVID-19 patients from a northern Italy Emergency Department Covid-19 in Critically Ill Patients in the Seattle Region-Case Series Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations