The outbreak of coronavirus disease 2019 (COVID-19), which is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has rapidly spread globally. On March 11, the World Health Organization (WHO) declared it as a public health emergency. 1
The clinical manifestations of COVID-19 are wide in spectrum, which range from asymptomatic or mildly symptomatic infections to severe pneumonia, acute respiratory distress syndrome (ARDS), and respiratory failure. Patients with severe COVID-19 have a substantial risk of prolonged critical illness and death. However, most previous studies of COVID-19 focused primarily on epidemiological and clinical characteristics; only limited evidence is available on risk factors for poor clinical outcomes. Reference Zhou, Yu and Du2-Reference Cummings, Baldwin and Abrams5
Here, we present details of the clinical manifestations, laboratory findings, imaging features, and clinical outcomes of patients with COVID-19 admitted to Huanggang Central Hospital, in Hubei. We aim to identify risk factors associated with in-hospital death and construct a clinical risk model to predict the fatal outcome of patients with COVID-19 upon admission.
Patients and Methods
Study Participants
Between January 23, 2020, and March 5, 2020, all patients with confirmed and probable COVID-19 in Huanggang Central Hospital were admitted to the hospital. Our study enrolled 150 consecutive inpatients with confirmed COVID-19 who had a definite outcome (discharged or died) during this period. All patients were diagnosed with COVID-19 pneumonia according to the WHO interim guidance. 6 Huanggang Central Hospital, a tertiary hospital with 2500 beds, was a designated hospital for patients with COVID-19. The patients were mainly from Huangzhou District and surrounding towns. The ethics committee of Huanggang Central Hospital approved this study and granted a waiver of informed consent from the study participants.
Data Collection
Data on epidemiological and demographic characteristics, underlying diseases, clinical manifestations, laboratory findings, chest computed tomography (CT) imaging, and outcomes of enrolled 150 patients with COVID-19 were obtained from electronic medical records. Test results at baseline rather than the worst value during hospitalization were used to predict clinical outcome. A team of experienced respiratory clinicians reviewed, abstracted, and cross-checked the data.
SARS-CoV-2 RT-PCR Test
Throat-swab samples were obtained from all patients on admission and tested for the presence of SARS-CoV-2 infection by using real-time reverse transcriptase-polymerase chain reaction assays (RT-PCR) according to a previously described protocol. Reference Huang, Wang and Li7 Throat-swab specimens were obtained for RT-PCR re-examination every 2 to 3 d. For patients with repeated PCR tests, the first date of the result was recorded if the patients had consecutive negative results, while the latest result and date were recorded for patients who had inconsistent results of the consecutive tests.
Discharge
The criteria for discharge were absence of fever for at least 7 d, basically normal blood oxygen saturation without supplemental oxygen, substantial improvement in chest CT scans, and 2 consecutive throat-swab samples negative for SARS-CoV-2 RT-PCR test taken at least 24 h apart.
Definitions
Fever was defined as axillary temperature of at least 37.3°C. Sepsis and septic shock were defined according to the 2016 Third International Consensus Definition for Sepsis and Septic Shock. Reference Shankar-Hari, Phillips and Levy8 Acute kidney injury (AKI) was diagnosed according to the Kidney Disease Improving Global Guidelines (KDIGO). Reference Khwaja9 ARDS was diagnosed according to the Berlin Definition. Reference Ranieri, Rubenfeld and Thompson10 Acute cardiac injury was diagnosed if the serum levels of cardiac biomarkers (eg, cardiac troponin I) were above the 99th percentile upper reference limit or if new abnormalities were shown in electrocardiography and echocardiography. Reference Huang, Wang and Li7 The illness severity of COVID-19 was defined according to the Chinese management guideline for COVID-19 (version 6.0). 11 Coagulopathy was defined as a 3-s extension of prothrombin time or a 5-s extension of activated partial thromboplastin time. Hypoproteinemia was defined as blood albumin of less than 25 g/L. Exposure history was defined as exposure to people with confirmed SARS-CoV-2 infection or residence or travel history in Wuhan in the last month, especially in the last 2 wk. The worldwide accepted pneumonia severity scoring systems such as CURB-65.12] and Pneumonia Severity Index (PSI) Reference Lim, van der Eerden and Laing13 ] were used to assess pneumonia severity. Critical illness evaluation systems, including Sequential Organ Failure Assessment (SOFA) Reference Vincent14 ] and Acute Physiology and Chronic Health Evaluation II (APACHE II), Reference Knaus, Draper and Wagner15 were used to assess disease severity.
Statistical Analysis
Categorical variables were presented as numbers (percentages) and analyzed using χ2 test or Fisher’s exact test. Continuous variables with skewed distribution were shown as median (interquartile ranges) and analyzed using Mann-Whitney U test. If the missingness of continuous variables did not exceed 30%, the median was used to fill in. The variables, including baseline characteristics, laboratory findings, were excluded those with more than 30% missing values and unconfirmed indicators (eg, exposure, which was self-reported). Logistic regression was used for univariable and multivariable analyses to explore the risk factors associated with in-hospital death. Variables with P value <0.1 entered into logistic multivariable analysis. Continuous variables were dichotomized. It was determined that whether or not with multicollinearity among independent variables. Forward-stepwise regression was used to select variables. Area under curve (AUC), sensitivity, and specificity were analyzed using the receiver operating characteristic curve (ROC). A nomogram was established based on the results of multivariable analysis. The C-index, decision curve, and clinical impact curve were used to verify the nomogram. A 2-sided P < 0.05 was considered statistically significant. All statistical analyses were performed using SPSS (version 25.0), R program (version 3.5.1), and Graphpad Prism (version 8.0).
Results
Demographics and Clinical Characteristics
A total of 150 patients were enrolled in this study, including 31 dead patients and 119 discharged patients. As shown in Table 1, the nonsurvivors were older than the survivors (73 [IQR 62-79] y vs 48 [IQR 37-57]) y, P < 0.001). There were more male nonsurvivors (61.3%), while the survivors were dominated by females (58.0%). We noted that more nonsurvivors had hypertension, cardiovascular diseases, and cerebrovascular diseases than the survivors. The most common symptoms on admission were fever and cough, followed by fatigue. Dyspnea was more common in the nonsurvivors. Approximately one-third of the patients had a definite history of exposure before illness onset, and the median time from exposure to illness onset was 5.5 (IQR 4.8-10) d. General patients (83.2%) accounted for the majority in the survivor group, while severe and critical patients accounted for 93.5% in the nonsurvivor group. In addition, the nonsurvivors had higher CURB-65, PSI, SOFA, and APACHE II scores (Figure 1).
Abbreviations: PSI, Pneumonia Severity Index; SOFA, Sequential Organ Failure Assessment; APACHE||, Acute Physiology and Chronic Health Evaluation.
* Data available for 30 patients, including 4 non-survivors and 26 survivors.
Laboratory Findings and Imaging Features
Laboratory findings on hospital admission are summarized in Table 2. White blood cell counts and neutrophil counts were elevated, while lymphocyte counts and platelet counts were decreased in the nonsurvivor group compared with those in the survivor group. The nonsurvivors showed higher levels of total bilirubin, blood ureanitrogen (BUN), creatinine, fasting blood glucose, lactate dehydrogenase (LDH), and D-dimer, but lower albumin levels than the survivors. C-Reactive protein and procalcitonin (PCT) were also different between the nonsurvivors and survivors.
* Data available for 141 patients, missing for 9 survivors.
† Data available for 111 patients, including 21 non-survivors and 90 survivors.
‡ Data available for 143 patients, including 25 non-survivors and 118 survivors.
§ Data available for 93 patients, including 15 non-survivors and 78 survivors.
During the period from admission to death, laboratory indicators showed a dynamic change in the nonsurvivors, especially in the levels of lymphocyte counts, LDH, BUN, and D-dimer (Figure 2). The serial biomarker results were based on a very small subcohort because most patients were not continuously monitored during hospitalization.
The most common imaging features were ground glass opacification (86.0%) and bilateral pulmonary infiltration (82.7%). Pleural effusion (22.6% vs 2.5%, P < 0.001) occurred more frequently in the nonsurvivors.
Clinical Outcomes
As shown in Table 3, of the 150 patients, sepsis was the most frequently observed complication (37.3%), followed by respiratory failure (26.7%). All of the nonsurvivors experienced sepsis. The common complications of the nonsurvivors included respiratory failure (93.5%), ARDS (83.9%), coagulopathy (77.4%), acidosis (64.5%), septic shock (61.3%), and AKI (61.3%). For the nonsurvivors, the median time from illness onset to respiratory failure was 9 (interquartile ration [IQR], 7-10) d; the median time from illness onset to ARDS was 10.5 (IQR 9-12) d; the median time from illness onset to septic shock was 18 (IQR, 14-23) d; the median time from illness onset to AKI was 16 (IQR, 10-22) d.
Abbreviations: ARDS, acute respiratory distress syndrome.
The median time from illness onset to discharge was 25 (IQR, 19-30) d, whereas the median time to death was 22 (IQR 13-25) d. The median hospital durations of stay were 17 (IQR, 12-23) d for the survivors and 13 (IQR, 6-17) d for the nonsurvivors, respectively. The median duration of viral shedding of the nonsurvivors were longer than that of the survivors (12 [IQR, 7.8-15] d vs 6 [IQR, 4-10]) d; P < 0.001), even the virus was continuously detectable until death in 17 nonsurvivors.
Risk Factors for In-hospital Death
On the univariable analysis, the risk factors associated with in-hospital death were age, hypertension, diabetes, and cardiovascular disease. The elevated respiratory rate, PSI score, SOFA score, APACHE II score, white blood cell counts, total bilirubin, BUN, creatinine, LDH, D-dimer, and the decreased white blood cell counts, lymphocyte counts, platelet counts, and albumin were also associated with in-hospital death (Table 4).
Abbreviations: PSI, Pneumonia Severity Index; SOFA, Sequential Organ Failure Assessment; APACHE||, Acute Physiology and Chronic Health Evaluation||.
* Per 1 unit increase.
The multivariable logistic analysis indicated that higher SOFA score (odds ratio [OR], 3.077; 95% confidence interval [CI]: 1.848-5.122; P < 0.001), history of diabetes (OR, 10.474; 95% CI: 1.554-70.617; P = 0.016), and LDH greater than 245 U/L (OR, 13.169; 95% CI: 2.934-59.105; P = 0.001) were risk factors for in-hospital death (Table 4). The predictive equation for in-hospital death: logit(P) = −5.594 + 1.124 × SOFA + 2.578 × LDH ([LDH ≤ 245] = 0, [LDH > 245] = 1) + 2.349 × diabetes (without diabetes = 0, with diabetes = 1), P = elogit(P) /1 + elogit(P).
Development and Validation of a Nomogram for Predicting In-hospital Death
The independent predictors were used to establish a nomogram, and the points of each variable are shown in Figure 3A. The ROC curve analysis revealed that the AUC of the nomogram was 0.970 (95% CI: 0.947-0.992), showing good accuracy in predicting the risk of in-hospital death (Figure 3B). The sensitivity was 90.3%, and the specificity was 90.8%. The favorable calibration curve indicated that the prediction by the nomogram was highly consistent with the actual observation (Figure 3C). Moreover, the clinical impact curve (Figure 3D) indicated that the nomogram had good net benefits for the identification of the fatal outcome of patients with COVID-19.
Fifty-two of 512 hospitalized patients with COVID-19 between March 5, 2020, and May 1, 2020, including all 12 dead patients and 40 discharged patients who were randomly selected according to hospitalization number, were enrolled to perform an external validation. The AUC of the nomogram was 0.923 (95% CI: 0.828-1.000), and the predictive accuracy was 0.942. The model was confirmed to be reliable.
Discussion
Most of patients with COVID-19 present with mild flu-like symptoms and recover quickly. However, many severe patients show rapid progression and develop multiple organ dysfunction, even death, indicating that recognition of risk factors are essential to identify those potentially needing critical care and management at an early stage. Several studies focused on risk factors associated with poor prognosis, including elder age, comorbidities, lymphopenia, D-dimer greater than 1 μg/L, elevated CRP, high LDH, high hypersensitive troponin I, high interleukin-6, hyperglycemia, and hypoproteinemia. However, these clinical and laboratory indicators are independent, Reference Zhou, Yu and Du2,Reference Cummings, Baldwin and Abrams5,Reference Li, Xu and Yu16-Reference Zhang, Yu and Tong18 which cannot comprehensively reflect the prognosis of COVID-19.
To our knowledge, no mature and reliable scoring system has been established to assess the risk of death in hospital with COVID-19. SOFA score is a low-cost, 2-min bedside clinical tool introduced to facilitate early recognition of sepsis and multi-organ dysfunction. SOFA score is a consistent and convincing tool that can be used across different patient populations and clinical settings. In the current study, we determined that sepsis occurred in 37.3% of patients with COVID-19 due to viral infection. We also identified SOFA score as an independent risk factor for death in adults who were hospitalized due to COVID-19.
Reports from the Chinese Center for Disease Control and Prevention from 44,672 confirmed cases of COVID-19 showed that patients with diabetes had a 3-fold higher overall case-fatality rate than those without diabetes (7.3% vs 2.3%, respectively). Reference Wu and McGoogan19 According to our results, we believed that increasing odds of death in hospital was associated with history of diabetes. The elevated LDH was confirmed as a predictive biomarker of death, consistent with those in previous reports. Reference Li, Xu and Yu16,Reference Cecconi, Piovani and Brunetta20 Furthermore, our study developed a nomogram model based on risk factors selected by multivariable logistic regression analysis to accurately predict fatal outcomes of patients with COVID-19. The nomogram model performed well in predicting in-hospital death, supported by external validation.
Few studies reported the long shedding of SARS-CoV-2 RNA, especially in severe patients. Reference Zhou, Yu and Du2,Reference Xiao, Tong and Zhang21 Among male patients, delayed admission to hospital after illness onset, and invasive mechanical ventilation were associated with prolonged SARS-CoV-2 RNA shedding. Reference Xu, Chen and Yuan22 In the current study, we found that the duration of SARS-CoV-2 viral shedding had median values of 6 d in the survivor group and 12 d in the nonsurvivor group; 17 patients were nucleic acid positive until death. Although most studies have used qualitative or quantitative PCR tests as a diagnostic marker for infectious SARS-CoV-2, caution is required when applying such data to assess the duration of viral shedding and infection potential because PCR does not distinguish between infectious virus and noninfectious nucleic acid.
COVID-19 mainly injures the respiratory system, and some patients rapidly progress to ARDS. We found that the incidence of ARDS was 17.3%, and the median time from illness onset to ARDS was 10.5 (IQR, 9-12) d, consistent with previous reports. Reference Zhou, Yu and Du2,Reference Wang, Hu and Hu23 At present, the mechanisms of COVID-19-related ARDS remain unclear. The entry of pathogenic SARS-CoV-2 in humans leads to the activation of inflammatory cells, specifically CD4 lymphocytes, which subsequently transform into T helper 1 (Th1) cells. Activated inflammatory cells (Th1 cells and macrophages) enter the pulmonary circulation and induce cytokines (ie, “cytokine storm”), which lead to rapid and wide-spread damage of the pulmonary epithelium and alveolar cells.
A recent study of the single-cell transcriptome analysis found that ACE2 genes were significantly expressed in podocytes and proximal convoluted tubules as potential hosts cells targeted by SARS-CoV-2; this work suggests that the kidney might be an important target organ for SARS-CoV-2. Reference Pan, Xu and Zhang24 A retrospective cohort study reported that mortality was higher in COVID-19 with AKI versus COVID-19 patients without AKI (60.5% vs 27.4%), and AKI was an independent predictor of mortality. Reference Kolhe, Fluck and Selby25 Our study revealed that AKI incidence rates of 14.0% in all hospitalized patients and 61.3% in the nonsurvivors. These data provide robust evidence to support that patients with COVID-19 should be closely monitored for the development of AKI and measures taken to prevent it.
This study has several limitations. First, our study was conducted in 1 hospital, thereby potentially limiting the generalizability to hospital settings, especially in terms of the demographic characteristics of the patient population and the external validation of the prediction models. Second, missing data on some variables, such as erythrocyte sedimentation rate (ESR), PCT, and cardiac troponin I, may cause bias in the estimation and reduce the representativeness of the samples. Finally, interpretation of our findings might be limited by the small sample size.
In conclusion, we found that high SOFA score, history of diabetes, and LDH greater than 245 U/L were risk factors for in-hospital death of adult patients with COVID-19. The nomogram proposed in our study objectively predicted the prognosis of patients with COVID-19.
Author Contributions
Yuanyuan Niu, Zan Zhan, and Jianfeng Li contributed equally to this work.
Conflicts of Interest
The authors declare no conflict of interests.