Introduction
The pandemic of the novel corona virus disease 2019 (COVID-19), which originally emerged in Wuhan, China in December 2019 has spread around the world [Reference Huang1, Reference Chen2]. As of 5 June 2020, the WHO has reported a total of 6 535 354 COVID-19 cases and 387 155 deaths globally, with an average mortality of 5.92% and the person-to-person transmission is still continuing [3]. The clinical spectrum of COVID-19 appears to be wide, ranging from asymptomatic infection to mildly, severely, critically ill cases. Mild patients present only upper respiratory tract symptoms like cough and fever, however, respiratory failure, acute respiratory distress syndrome, heart failure, septic shock and even death can be observed in patients with critical conditions [Reference Zhou4]. Although most confirmed patients (81%) were classified as mild or moderate, 14% were severe and 5% were critical according to the largest investigation of 72 314 cases to date [Reference Wu5]. Accumulated evidences have indicated that older age, male, smoking, comorbidity, neutrophilia, coagulopathy, elevated D-dimer level and organ dysfunction were associated with increased risk of death from COVID-19 [Reference Wu5–Reference Wang10]. However, investigations of determinants of severity of COVID-19 are scarce. Early detecting cases with the potential deterioration of disease will contribute to optimising the use of limited health resources and allocating the proper care. Liang et al. developed a clinical risk score to predict the occurrence of critical COVID-19 based on severe or non-severe [Reference Liang11]. To our knowledge, no previous studies have been conducted to investigate the risk factors of severity of COVID-19 based on ordinal response, namely moderate, severe and critical illness. The estimation of risk factors of disease severity is therefore not very robust.
Here, we conducted a retrospective study based on COVID-19 patients from four designated hospitals in Shanghai, Hubei and Anhui provinces to describe the clinical features of COVID-19, and aimed to identify the predictors of multi-level response of severity from moderate, severe to critical illness.
Methods
Study design and participants
This multi-centre retrospective study encompassed COVID-19 patients classified as being moderately, severely and critically ill. The illness severity of COVID-19 was defined according to the Guideline on the Diagnosis and Treatment of COVID-19 by the National Health Commission (V.5) as described previously [Reference Feng12]. Patients were admitted to Shanghai Public Health Clinical Center, Wuhan Jinyintan Hospital and Tongji Hospital of Tongji Medical College HUST in Hubei province and Tongling Municipal People's Hospital in Anhui province from 1 January 2020 to 8 March 2020. All patients recruited in this study were laboratory-confirmed COVID-19. The study was approved by the Ethics Committees of these four hospitals, respectively. Written informed consent was waived owing to the need of rapid emergency response to this infectious disease.
Data collection
Medical records of COVID-19 patients were reviewed by the research team, and demographic, epidemiological, clinical, laboratory, treatment and outcome data were retrieved from electronic medical records using a standardised case report form. All data were cross-checked by two experienced doctors. To ascertain the medical histories or epidemiological data, which were unavailable from electronic medical records, the patients or their close relatives were interviewed by researchers. Data from the medical records were adopted if there was a discrepancy between the subjective description and the medical records.
Laboratory procedures
Method of laboratory confirmation of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been described elsewhere [Reference Huang1]. Simply, the Chinese Center for Disease Control and Prevention (CDC) and local CDC were in charge of detecting SARS-CoV-2 in throat-swab specimens from the upper respiratory tract by real-time reverse transcription polymerase chain reaction assay (RT-PCR). The criteria of discharge included absence of fever for at least 3 days, remission of respiratory symptoms, complete improvement in bilateral lungs in chest CT, together with negative for 2 times in throat-swab samples for SARS-CoV-2 RNA at least 24 h apart.
Initial clinical laboratory examinations involved complete blood count, serum biochemical tests (including liver and kidney functions, creatine kinase, lactate dehydrogenase (LDH) and electrolytes), myocardial enzymes, D-dimer and procalcitonin (PCT). Frequency of examinations was under the discretion of treating physicians. Chest computed tomographic (CT) scans were carried out for all COVID-19 patients. Two radiologists were invited to interpret chest CT scans independently and were blinded to the severity of the patient. When disagreement arose, a third radiologist was consulted to reach a final decision.
Statistical analysis
Continuous variables were presented as median with interquartile range (IQR) and the analysis of variance (ANOVA) or Kruskal−Wallis H test were used to compare the difference among three groups as appropriate. Categorical variables were expressed as frequency with percentages, and were analysed by Pearson's χ 2 test or Fisher's exact test. Bonferroni's correction was used for pairwise comparison. All patients were divided into moderate, severe and critical illness groups. Potential predictive variables included the following case characteristics on admission: demographic and epidemiological features, comorbidity, clinical signs and symptoms, laboratory findings and chest imaging results. To explore the risk factors associated with illness severity of COVID-19, namely moderately, severely and critically ill, which means the response variable was ordinally scaled, a cumulative logit model was used to investigate the effect of predictors of COVID-19 severity. Imputation for missing variables of some patients at hospital admission was considered if missing values were less than 20%, and imputation based on the expectation−maximisation algorithm method was used to replace missing values. Before ordinal logistic regression model was fitted, continuous variables of laboratory findings were transformed into categorical variables according to their reference values. The univariate and multivariate cumulative logit models were fitted with moderate illness as the reference level. Potential predictors of severity were investigated using univariate ordinal logistic regression firstly. We further conducted a backward stepwise multivariate ordinal logistic regression analysis excluding variables which were not significant in univariate cumulative logit model. Since missing rate of 34.6% occurred in the lung imaging results and over 40% existed in urine protein and urine glucose, these variables were excluded from multivariate ordinal logistic model. The overall survival (OS) was estimated using the method of Kaplan−Meier and the log-rank test was applied to compare the survival difference among different severity illness groups. The hazard ratio with 95% confidence interval (CI) was estimated with Cox proportional hazard model. A two-sided α of less than 0.05 was considered statistically significant. All statistical analyses were conducted using SAS software (V. 9.4) (SAS Institute Inc., USA).
Results
Demographics, laboratory findings and clinical course
As of 28 April 2020, data from 598 COVID-19 cases admitted to these four hospitals, including 400 (66.89%) moderate cases, 85 (14.21%) severe cases and 113 (18.90%) critical cases, had been collected to be incorporated into this study, of whom 79 cases had died during hospitalisation, with an average mortality of 13.21%, and 457 cases had recovered and been discharged. The remaining 62 cases were still in hospitals. The median age of the 598 patients was 57 years (IQR 42–66), ranging from 11 to 89 years, and 58.03% patients were male (Table 1). At least one comorbidity was present in 51.54% of patients, with hypertension being the most frequent comorbidity (33.90%), followed by diabetes (13.18%) and cardiovascular disease (6.51%). Few cases had a current (7.53%) or former (2.57%) smoking habit. The most common symptoms on admission were fever (81.22%) and dry cough (33.63%), followed by sputum production (30.77%) and shortness of breath (26.65%). Overall, the median interval between illness onset and confirmed diagnosis was 4 days (IQR 2–7), whereas the median interval between illness onset and admission was 7 days (IQR 3–11).
Clinical characteristics and comorbidities of 598 COVID-19 patients
COPD, chronic obstructive pulmonary disease.
Data are median (IQR) or n/ total (%). P value denotes the comparison among moderate, severe and critical illness group. * and †Signify P < 0.05 for pos-hoc comparison. *Refers to comparison between the critical group and the severe group or the moderate group. †Refers to comparison between the severe group and the moderate group.
The substantial differences of laboratory findings on admission among three groups of patients were observed (Table 2). The results indicated that significantly higher proportion of patients showing abnormal alanine aminotransferase (ALT) and aspartate aminotransferase (AST), and significantly higher levels of C-reactive protein, neutrophil count, serum potassium, cardiac troponin I (cTnI), myohaemoglobin, PCT, brain natriuretic peptide (BNP) and D-dimer were observed in the critically ill group and the severely ill group than in the moderately ill group, whereas, the levels of haemoglobin and serum albumin were significantly lower in critically ill group and severely ill group as compared with moderately ill group (P < 0.05). The critical group showed a significantly higher level of platelet count, fibrinogen and serum calcium as compared with severe group and moderate group. The severe group showed a significantly lower level of lymphocyte count than the moderate group. The normal ranges of laboratory indicators are shown in Supplementary Table S1. Abnormalities on chest radiographs on admission were seen in most patients (Table 2). Overall, typical findings on chest CT images were ground-glass opacity (94.37%), followed by pleural thickening (47.06%) and consolidation (35.74%). Chest CT scans showed significantly higher percentage of bilateral lungs involvement in critical group (92.45%) and severe group (95.24%) as compared with moderate group (80.07%). The median lung lobes involved in critical group (5, IQR 5–5) and severe group (5, IQR 2.75–5) were greater than those in moderate group (4, IQR 2–5).
Laboratory and chest CT findings on admission of 598 COVID-19 patients
ALT, alanine aminotransferase; AST, aspartate aminotransferase. LDH, lactate dehydrogenase; CK-MB, creatine kinase isoenzymes; cTnI, cardiac troponin I; PCT, procalcitonin; ESR, erythrocyte sedimentation rate; BNP, brain natriuretic pepti.
Data are median (IQR) or n/ total (%). P value denotes the comparison among moderate, severe and critical illness group. * and †Signify P < 0.05 for post-hoc comparison. *Refers to comparison between the critical group and the severe group or the moderate group. †Refers to comparison between the severe group and the moderate group.
Totally, 320 (61.07%) patients were given antivirals within 2 days after admission including lopinavir/ritonavir, arbidol, darunavir and chloroquine. In all, 369 (70.83%) patients received antibiotics and 117 (23.08%) received corticosteroids. More patients received corticosteroids in critical group and severe group as compared with the moderate group (P < 0.05). The proportions of patients accepting high-flow nasal cannula oxygen therapy and non-invasive mechanical ventilation, respectively, in critical and severe groups were significantly higher than in mode rate group (P < 0.05). Compared with the moderate group, the critical group and the severe group had a significantly lower rate of discharge and a higher mortality rate (P < 0.05), as shown in Table 3. The comparison of demographic and baseline characteristics, symptoms, laboratory parameters, lung image features, treatment and prognosis among moderately, severely and critically ill patients were shown in Tables 1–3.
Treatment and prognosis of 598 COVID-19 patients
Data are median (IQR) or n/ total (%). P value denotes the comparison among moderate, severe and critical illness group. * and †Signify P < 0.05 for post-hoc comparison. *Refers to comparison between the critical group and the severe group or the moderate group. †Refers to comparison between the severe group and the moderate group.
Determinants of illness severity
Fifty-three variables on admission were successively included in the univariate ordinal logistic regression, and 35 variables were found to be associated with illness severity, including age, gender, hypertension, diabetes, interval between illness onset and diagnosis, interval between illness onset and admission, pharyngodynia, shortness of breath, early administration of antiviral, C-reactive protein, white blood cell (WBC) count, neutrophil count, lymphocyte count, haemoglobin, platelet count, ALT, AST, albumin, blood urea nitrogen, creatinine, potassium, LDH, creatine kinase, myohaemoglobin, troponin I (cTnI), PCT, erythrocyte sedimentation rate (ESR), BNP, fibrinogen, D-dimer, bilateral lungs involved, consolidation, linear opacity, pleural thickening and lung lobes involved (Table 4).
Results of univariate ordinal logistic model using three levels of severity as response
COPD, chronic obstructive pulmonary disease; CRP, C-reactive protein; WBC, white blood cell; ALT, alanine aminotransferase; AST, aspartate aminotransferase; LDH, lactate dehydrogenase; CK-MB, creatine kinase isoenzymes; cTnI, cardiac troponin I; PCT, procalcitonin; ESR, erythrocyte sedimentation rate; BNP, brain natriuretic peptid.
Except five lung image variables, 30 significant predictors of severity in univariable analysis were included in a multivariable stepwise cumulative logit model, and seven variables retained in the final model which were statistically significant independent determinants of COVID-19 illness severity (Table 5). The results of multivariate model revealed that the risks of having more severe illness were 1.586 (95% CI: 0.824–3.053) and 3.419 (95% CI: 1.596–7.323) times higher among patients belonging to the age group 40–69 and 70+ years, respectively, when compared with patients of less than 40 years. Patients with hypertension had 3.372 (95% CI: 2.185–5.202) times greater risk of having worse severity of illness compared with patients without hypertension. The risk of having worse severity of illness was found higher for patients with ALT>50 μ/l (OR = 3.304; 95% CI: 2.107–5.180) when compared with those having ALT⩽50 μ/l. The risk of having worse severity of illness was found significantly higher for patients with higher cTnI (>0.04 ng/ml) than those with normal cTnI (⩽0.04 ng/ml) with OR being 7.464 (95% CI: 4.292–12.980). COVID-19 patients with myohaemoglobin>48.8 ng/ml at admission had a 2.214 (95% CI: 1.42–3.453)-fold greater risk of having worse severity of illness when comparison was made with patients having normal myohaemoglobin level. Table 5 also shows that interval between illness onset and diagnosis and interval between illness onset and admission were independent significant predictors of illness severity with OR being 1.056 (95% CI: 1.012–1.101) and 1.048 (95% CI: 1.009–1.087), respectively.
Results of multiple ordinal logistic model using three levels of severity as response
ALT, alanine aminotransferase; cTnI, cardiac troponin I.
Post-hoc comparison of survival among patient groups with different risk factors
In the critical illness group, 68 of 113 patients (60.18%) died, as compared with five of 85 (5.88%) in the severe illness group and six out of 400 (1.50%) in the moderate illness group. The median OS was 29 days in critical group, as compared with not attainable in severe group and moderate group. The 30-day OS rates were 97.7% (95% CI: 95.5–100%), 95.3% (95% CI: 90.2–100%) and 46.4% (95% CI: 37.4–57.6%) in moderate, severe and critical groups, respectively (P < 0.001), as shown in Figure 1. Patients of critical illness suffered from inferior survival, as compared with patients in severe group (HR = 14.309, 95% CI: 5.585–36.659) and those in moderate group (HR = 41.021, 95% CI: 17.588–95.678), representing 14.309 times greater risk of death and 41.021 times greater risk of death when compared with severe and moderate groups, respectively. As shown in Supplementary Table S2, COVID-19 patients of older age (HR = 9.823 for ⩾70 vs. <40; HR = 3.361 for 40–69 vs. <40), comorbidity of hypertension (HR = 3.161), abnormal ALT (HR = 1.657), abnormal cTnI (HR = 2.513) or abnormal myohaemoglobin (HR = 2.671) suffered from inferior survival. Also, the Kaplan−Meier survival curves (Supplementary Figs S1–S5) demonstrated that the differences of OS with respect to stratification by these risk factors were all statistically significant (P < 0.05).
Kaplan−Meier estimate of OS of COVID-19 patients according to severity of illness.
Discussion
COVID-19 has presented an unprecedented challenge to the health-care system across the world. With the increasingly scarce health resources, mortality is the most important issue when dealing with epidemics. Early identifying potential of severe and critical patients becomes the priority in minimising the mortality and contributes to allocation of limited critical care. Previously, a risk forecasting model to predict the occurrence of critical illness among hospitalised COVID-19 patients in China has been reported by Liang et al. [Reference Liang11]. The response variable in Liang's model had just two levels of severe and non-severe. Actually, the mortalities of COVID-19 patients with different severity are variant. This study demonstrated the mortality of 60.18% in critical cases, followed by 5.88% in severe cases and 1.50% in moderate cases. Kaplan−Meier analysis also showed that patients of critical illness suffered from inferior survival, as compared with patients of severe illness and moderate illness, indicating that illness severity is related to the prognosis. In our opinion, more precise classification of illness severity into three levels suits the clinical spectrum of COVID-19 and the predicting model based on multi-level ordinal response variables can be more practical. The current multi-centre retrospective study identified the risk factors of illness severity among COVID-19 inpatients from four hospitals of Hubei, Anhui and Shanghai. According to the Guideline on the Diagnosis and Treatment of COVID-19 by the National Health Commission, inpatients of COVID-19 in this study were divided into moderately ill, severely ill and critically ill groups, and the ordinal logistic regression model was fitted to identify the predicators of severity of illness. Based on current data, older age, comorbidity of hypertension, elevated levels of ALT, elevated cTnI, elevated myohaemoglobin, together with prolonged interval between illness onset and diagnosis and interval between illness onset and admission were independent determinants of severity of COVID-19 and represented higher odds of worse severity of illness.
Evidence is gradually accumulating with regard to the risk factors associated with severity of COVID-19. As described previously, older age has been reported as an important independent predictor of severity in COVID-19 patients [Reference Wu5, Reference Wu6, Reference Liang11], which is proved by this study. Comorbidity of hypertension was found to be associated with an increased risk of death. A meta-analysis including 46 248 patients with confirmed COVID-19 indicated that those with the most severe illness were more likely to have hypertension with OR of 2.36 (95% CI: 1.46–3.83). Hypertension was reported to increase the OR for death by 3.05 (95% CI: 1.57–5.92) in patients with COVID-19 [Reference Zhou4]. Similarly, patients with hypertension had 3.372 (95% CI: 2.185–5.202) times greater risk of developing worse severity of illness in comparison with those without hypertension in our study. The relationship between hypertension and COVID-19 may relate to the role of angiotensin converting enzyme (ACE2) [Reference Kreutz13]. As a key element in the renin–angiotensin–aldosterone system (RAAS), ACE2 is critically involved in the pathophysiology of hypertension. Studies demonstrated that inhibition of the RAAS with ACE inhibitors (ACEIs) or angiotensin II receptor blockers (ARBs) may result in a compensatory increase in tissue levels of ACE2 and poorer clinical course and prognosis, leading to suggestions that these drugs may be detrimental to COVID-19 patients [Reference Vaduganathan14].
It has been reported that COVID-19 had significant impact on the liver function [Reference Zhang15, Reference Lei16]. A meta-analysis including 20 retrospective studies with 3428 COVID-19 patients revealed that higher serum levels of AST (mean difference = 8.84 U/l, 95% CI: 5.97–11.71, P < 0.001) and ALT (mean difference = 7.35 U/l, 95% CI: 4.77–9.93, P < 0.001) and lower serum levels of albumin (mean difference = −4.24 g/l, 95% CI: −6.20 to −2.28, P < 0.001) were associated with a significant increase in the severity of COVID-19 [Reference Parohan17]. Our univariate logistic regression results also showed patients with elevated AST, ALT and decreased albumin had 3.410-, 3.783- and 5.734-fold greater risk of worse severity, which was consistent with this meta-analysis. Particularly, our multivariate ordinal regression demonstrated that ALT is an independent predictor of severity of COVID-19 patients. Since elevated liver injury indicators are strongly associated with the severity risk and subsequent death risk, the liver function should be monitored during hospitalisation. Acute myocardial injury is the most commonly described cardiovascular complication in COVID-19 [Reference Bansal18]. The overall incidence of acute myocardial injury has been variable but roughly 8–12% of COVID-19 patients are found to develop significant elevation of cTnI [Reference Lippi19]. The patients admitted to ICU or having severe/fatal illness have several-fold higher likelihood of cTnI elevation. Several recent studies indicated that higher concentration of cTnI and myohaemoglobin were associated with the severity and case fatality rate of COVID-19 [Reference Chen20–Reference Gaze22]. Chen et al. reported that elevated cTnI (OR = 26.909, 95%CI: 4.086–177.226, P = 0.001) were the independent risk factors of critical disease status [Reference Chen20]. Han et al., found there were statistically significant differences in the level and positive rate of cTnI and myohaemoglobin among the mild, severe and critical COVID-19 case groups [Reference Han23]. Our results are in agreement with the previous studies that the elevated myocardial injury markers such as cTnI and myohaemoglobin are independent determinants of illness severity in COVID-19 patients representing negative clinical course and potentially life-threatening prognosis. In particular, cTnI is the strongest predictor of worse severity with OR being 7.464 (95% CI: 4.292–12.980). Direct myocardial injury due to viral myocarditis or the effect of systemic inflammation appears to be the most common mechanisms of acute cardiac injury. Our study also found that interval between illness onset and diagnosis and interval between illness onset and admission were independently associated with illness severity. The risk of developing worse illness severity increased by 1.056-fold and 1.048-fold for each day delay of interval between illness onset and diagnosis and interval between illness onset and admission, respectively. Chen et al. [Reference Chen24] found the median time from symptom onset to admission was 10 days among deceased COVID-19 patients, one day longer than that of discharged patients. Considering the limited health resources and high case volume in some regions currently, adjusted tactics and strategies should be taken to maximise the availability of and accessibility to medical service to shorten the diagnosis delay or admission delay.
Our study had several limitations. First, this is a retrospective study design, which could be subject to recall bias and selection bias. Second, not all laboratory parameters were tested in all patients, including LDH and D-dimer. Although imputation technique was used to replace missing values, their role might be underestimated in predicting illness severity. Third, due to massive loss of chest CT results, the predicting role of chest CT abnormalities could not be evaluated in this study. Its role of predicting critical illness of COVID-19 has been demonstrated by other studies. Last but not least, generalisability of our findings might be limited by the sample size, and the results need to be validated based on a much larger patient population.
In this study, we identified older age, presence of hypertension, elevated ALT, cTnI and myohaemoglobin, prolonged interval between illness onset and diagnosis and admission as the independent determinants to predict the risk of developing more severe illness among COVID-19 patients. Given the ongoing global pandemic of COVID-19, this study will contribute to early identifications of patients with high risk of developing critical illness and optimising the arrangement of health resources.
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1017/S0950268820001533.
Acknowledgements
We acknowledge all health-care workers involved in the diagnosis and treatment of patients in four sites and we thank Dr Qiqi Cao and Dr Weixia Li for assistance in the analysis and interpretation of the lung image data.
Author contributions
Dr Jian Li had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis. Study design: Jian Li; acquisition of data: Kandi Xu, Dexiang Yang, Yun Ling, Kui Liu, Tao Bai; analysis and interpretation of data: Jian LI, Kandi Xu, Min Zhou, Zenghui Chen; funding acquisition: Min Zhou; drafting of the manuscript: Jian Li; statistical analysis: Jian Li, Kandi Xu; critical revision of the manuscript: all authors.
Financial support
This study was funded by grants from the National R&D Program of China, Shanghai Key Discipline for Respiratory Diseases, Shanghai Shenkang Hospital Development Center Clinical Science and Technology Innovation Project and Zhejiang University Special Scientific Research Fund for COVID-19 Prevention and Control.
Conflict of interests
We declare no competing interests.
Data availability statements
The data that support the findings of this study are available from Clinical Research Center, Ruijin Hospital, Shanghai Jiao Tong University School of Medicine. Restrictions apply to the availability of these data, which were used under license for this study.
Open access
