Introduction
The major clinical effects of coronavirus disease 2019 (COVID-19) infection are on the respiratory system although other systems can be affected [Reference Chang, Yan and Wang1, Reference Li, Bai and Hashikawa2]. COVID-19 may result in acute respiratory failure requiring mechanical ventilation and even leading to death [Reference Chen3–Reference MacLaren, Fisher and Brodie6]. A study in Italy reported rates of respiratory failure as high as 29%–40% [Reference Bartoletti7]. The 28-day mortality could occur in 26%–30% of patients with COVID-19 who had respiratory failure necessitating invasive mechanical ventilation (IMV) [Reference Ferguson8]. Since the COVID-19 pandemic is still evolving, the true mortality has not been defined, but the crude mortality ratio (the number of reported deaths divided by the number of reported cases) has been estimated to be 3%–4%, which appears to be higher than that for influenza [9]. However, the infection mortality rate (the number of reported deaths divided by the number of infections) is lower than the crude mortality ratio; and the mortality rate varies among different regions, demographic and socioeconomic characteristics, levels of healthcare access and quality, intervention methods and qualities of reported deaths and cases [9–Reference Sun11]. Identifying risk factors of respiratory failure in patients with COVID-19 could help clinicians recognise patients at high risk of respiratory failure and hence take active treatment for them to prevent further worse outcomes and reduce emergency intubation or cardiopulmonary resuscitation to protect medical staff from related infections.
For patients with COVID-19, risk factors of severe COVID-19, admission to intensive care unit (ICU) and death have been reported in many studies [Reference Zhao12–Reference Noor and Islam15], but studies on risk factors of respiratory failure remain scarce. To the best of our knowledge, only two studies are available in the literature identifying the predictors of respiratory failure and both studies were conducted in Italy [Reference Bartoletti7, Reference Rottoli16]. However, demographic characteristics may vary across different countries, and the levels and risk factors of respiratory failure may also differ [Reference Richardson17]. Therefore, information on the respiratory failure in other settings than Italy will be of great academic and clinical value.
Based on a large multicentre retrospective cohort with rich information on demographic, epidemiologic, clinical and laboratory features as well as CT imaging features, this study aims to assess the level and identify potential risk factors of respiratory failure in patients with COVID-19.
Methods
Study design and participants
The study design and participants have been described in our previous study [Reference Luo18]. Here we briefly describe it. This is a retrospective multicentre cohort study. Patient inclusion criterion was as of 15th March 2020, all patients diagnosed with COVID-19 at 24 hospitals designated to treat COVID-19 in Jiangsu province, China according to the diagnostic criteria of the ‘Diagnosis and treatment protocol for novel coronavirus pneumonia (trial version 7)’ released by National Health Commission and National Administration of Traditional Chinese Medicine of China [19]. The diagnosis of COVID-19 was based on epidemiological history, clinical and CT manifestations and laboratory confirmation (real-time reverse transcriptase-polymerase chain reaction assays [RT-PCR] to detect aetiological agent severe acute respiratory syndrome coronavirus-2, SARS-CoV-2, which caused COVID-19) [Reference Sun11, 19]. The exclusion criterion was patients with no available medical records. The criteria for discharge were: the patient's body temperature remained normal for more than 3 days, the symptoms were relieved (if there were symptoms), and the results of two consecutive RT-PCR assays were negative (throat swab samples, at least 1 day apart).
Data collection and definition of variables
The primary outcome was the occurrence of acute respiratory failure during 14-day follow-up (days 1, 2, 3, 4, 5, 6, 7 and 14 from admission). The respiratory failure was defined as oxygen saturation (SpO2) <93% and/or partial pressure of oxygen in arterial blood (PaO2) <60 mmHg on room air and/or requirement of high-flow nasal cannula oxygen therapy (HFNC), non-invasive or invasive mechanical ventilation. Type 1 (hypoxaemic) respiratory failure refers to the hypoxaemia (PaO2 <60 mmHg (8 kPa)) with the normal (normocapnia) or low (hypocapnia) partial pressure of carbon dioxide (PaCO2) in arterial blood. Type 2 (hypercapnic) respiratory failure refers to the hypoxaemia with the hypercapnia (PaCO2 >45 mm Hg (6 kPa)).
Demographic features analysed in this study included sex and age (year). Epidemiologic features included exposure types (imported case or local case) and types of disease onset (single onset or clustered onset). Clinical features included initial symptoms (fever, cough, sputum, shortness of breath, anorexia and diarrhoea), medical history (hypertension, coronary heart disease, diabetes, stroke, smoking and drinking) and vital signs (temperature (°C), heart rate (HR, beats/minute), systolic blood pressure (SBP, mmHg), diastolic blood pressure (DBP, mmHg), respiratory rate (breaths/min) and SpO2 (%)). Laboratory features included blood test parameters (white blood cell count (WBC, 109/l), neutrophil count (109/l), lymphocyte count (109/l), haemoglobin (g/l) and platelet count (109/l)), organ function parameters (albumin (g/l) and creatinine (μmol/l)), inflammatory factors (C-reactive protein (mg/l)) and coagulation function parameters (activated partial thromboplastin time (APTT, s), fibrinogen (g/l) and D-dimer (mg/l)). Radiologic features included lesion distribution (outer third of lung involved, middle third of lung involved or inner third of lung involved), lesion density (below 20% consolidation, 20−80% consolidation, or above 80% consolidation), lesion border (well-defined border, moderately defined border or ill-defined border), quadrant score and pulmonary opacity score. Data of features listed above were collected at admission. Data of treatments (supportive treatments and medication) were obtained within the whole study period.
Imported cases were defined as those who had been to the pandemic centre of China (Wuhan city), or had contact with people or patients with COVID-19 who had been to Wuhan; and other cases were classified as ‘local cases’. A clustered onset is defined as the occurrence of two or more confirmed COVID-19 cases in the same cluster/group within 14 days, such as family, community, hospital, workplace or public place. A clustered onset may occur from interpersonal transmission via close contact with or joint contact with a confirmed COVID-19 case. Other cases not meeting the criteria for a clustered onset were classified as ‘single onset’.
Radiologic features were evaluated visually, by two radiologists with more than 5 years’ working experience. The radiologists were blinded to patients’ other characteristics and would reach agreements on different assessments of radiologic features. Chest CT axial sections were divided into four quadrants (left, right, anterior and posterior) by drawing horizontal and vertical lines through the centre of the chest. The quadrant score was defined as the number of quadrants with pulmonary opacities extending from the proximal end to the distal end of the chest, ranging from 0 to 4; and pulmonary opacity score was defined as the percentage of pulmonary opacity area in the area of bilateral lungs, rounded to the nearest 5%.
Statistical analysis
Normally distributed continuous variables and skewed distributed continuous variables of patients were reported as mean (standard deviation, s.d.) and median (interquartile range, IQR) by group (patients with and without respiratory failure) and compared using Student's t test or Mann−Whitney U test, respectively. Categorical variables were summarised using frequency and percentage and compared by χ 2 or Fisher's exact test.
Logistic regression models were used to identify risk factors of having respiratory failure. Univariate logistic regression models were fitted first to evaluate associations between each variable measured at admission and respiratory failure on the complete cases (without missing value). Variables that were statistically significant in the univariate analysis were then included in the multivariate logistic regression model. In the multivariate analysis, missing covariates were imputed with multiple imputation using a Markov chain Monte Carlo simulation method with 10 iterations. A sensitivity analysis was performed on complete cases. Odds ratio (OR) of having respiratory failure for each variable was calculated along with 95% confidence interval (CI). The two-tailed P < 0.05 was considered as statistically significant. The analyses were performed using SAS 9.4 (SAS Institute).
Results
Of the 721 suspected cases with possible COVID-19 in Jiangsu province from 10th January to 15th March 2020, 90 cases were excluded since RT-PCR tests showed negative results and six cases were excluded due to no available medical records. Finally, 625 cases (52.6% male; median age 46 years old (IQR 32–57)) in 24 hospitals were included for analysis, mainly from the Second Hospital of Nanjing, Suzhou Infectious Disease Hospital and Huai'an No. 4 People's Hospital (113 (18.1%), 86 (13.8%) and 79 (12.6%), respectively). The remaining cases were disproportionately from the other hospitals. The hospital to which a patient was admitted was mainly determined by geographic location. Of the 625 patients, 56 (9%) had respiratory failure, mainly type I (hypoxaemic). At the study end point (15th March 2020), no patients died and all patients were discharged from hospitals.
At admission to hospital, compared with patients without respiratory failure, those with respiratory failure were significantly older (mean age 59.70 vs. 42.93 years, P < 0.0001); were more likely to be single onset (62.5% vs. 48.3%, P = 0.0498); were more likely to have symptoms including fever (82.1% vs. 64.3%, P = 0.0074), cough (71.4% vs. 53.4%, P = 0.0110), sputum (42.9% vs 25.0%, P = 0.0064) and shortness of breath (12.5% vs. 2.3%, P = 0.0010); were more likely to have prior histories of hypertension (32.1% vs. 13.9%, P = 0.0014), coronary heart disease (7.1% vs. 1.6%, P = 0.0224) and diabetes (17.9% vs. 5.3%, P = 0.0015); had higher mean temperature (37.25 vs. 37.03 °C, P = 0.0347), greater mean heart rate (90.71 vs. 86.82 beats/min, P = 0.0387), greater mean respiratory rate (21.13 vs. 18.88 breaths/min, P < 0.0001) and lower mean SpO2 (95.27% vs. 97.92%, P < 0.0001) (Table 1).
s.d., standard deviation; HR, heart rate; SBP, systolic blood pressure; DBP, diastolic blood; SpO2, peripheral capillary oxygen saturation.
Patients with respiratory failure had significantly lower median lymphocyte count (0.7*109/l vs. 1.3*109/l, P < 0.0001), median platelet count (155.5*109/l vs. 186.5*109/l, P = 0.0008) and median albumin levels (40.0 vs. 41.9 g/l, P = 0.0005). In addition, patients with respiratory failure had significantly higher median levels of C-reactive protein (40.6 vs. 10.0 mg/l, P < 0.0001), median fibrinogen levels (4.3 vs. 3.4 g/l, P < 0.0001) and median D-dimer levels (0.3 vs. 0.2 mg/l, P = 0.0004) (Table 2).
IQR, inter-quartile range; WBC, white blood cell; APTT, activated partial thromboplastin time.
For visually evaluated CT features at hospital admission, patients with respiratory failure had significantly greater median of CT quadrant score (4.00 vs. 2.00, P < 0.0001) and median of pulmonary opacity score (52.50% vs. 20.00%, P < 0.0001) (Table 3).
CT, computer tomography; IQR, inter-quartile range.
a The total number of cases is 484 (43 respiratory failures and 441 no respiratory failures).
b The total number of cases is 496 (47 respiratory failures and 449 no respiratory failures).
Oxygen was delivered to patients with respiratory failure via nasal cannula (80.4%), simple face masks (23.2%), HFNC (39.3%), non-invasive mechanical ventilation (NIV) (53.6%), IMV (8.9%) and in prone position (30.4%) (Table 4). Patients with respiratory failure were more likely to receive supportive treatments and medications (all P < 0.05), except for the interferon.
Twenty variables at admission were found to be related to the occurrence of respiratory failure in univariate logistic regression analysis (Table 5). When they were included in the multivariate logistic regression model simultaneously, four variables were independently related to the occurrence of respiratory failure: age (in years) (OR, 1.07; 95% CI: 1.03–1.10; P = 0.0002), respiratory rate (breaths/min) (OR, 1.23; 95% CI: 1.08–1.40; P = 0.0020), lymphocyte count (109/l) (OR, 0.18; 95% CI: 0.05–0.69; P = 0.0157) and pulmonary opacity score (per 5%) (OR, 1.38; 95% CI: 1.19–1.61; P < 0.0001).
OR, odds ratio; CI, confidence interval; HR, heart rate.
a Univariate analysis is based on the complete cases without missing value.
b Multivariate analysis is based on imputed values for missing data in lymphocyte, platelet, albumin, C-reactive protein, fibrinogen, D-dimer, quadrant score and pulmonary opacity using multiple imputation method.
The sensitivity logistic regression model with only above four significant variables was estimated on complete cases (without missing data), and these variables remained statistically significant.
Discussion
This study assessed the level and risk factors of respiratory failure among patients with COVID-19. During the 14-day follow-up, 9% (56 out of 625) of patients with COVID-19 suffered from respiratory failure (mainly type I, hypoxaemic) in Jiangsu province, China. At the end of the study (15th March 2020), no patients died and all patients were discharged from hospitals. Among many factors explored in this study, four of them (older age, increased respiratory rate, decreased lymphocyte count and greater pulmonary opacity score) were identified as independent risk factors of respiratory failure after controlling for other confounders. Therefore, patients with such factors need to be carefully and thoroughly managed by physicians.
Jiangsu province is non-neighbouring with and geographically distant around 600 km from Hubei province, where Wuhan city (the epicentre of COVID-19 pandemic in China) is located in. The rate of respiratory failure in Jiangsu province is similar to the figure in the national report from China [Reference Guan20] but this rate is much lower than reported earlier in Wuhan city (26%–32%) [Reference Huang5, Reference Liu21]. This may be due to the early responses and measures adopted by the Jiangsu provincial health authorities to deal with the disease during the pandemic, including more adequate medical resources, deeper understanding and better management of respiratory failure for patients with COVID-19 [Reference Sun11]. The rate of respiratory failure in this study is also lower than reported in other countries, including Italy (29%–40%) [Reference Bartoletti7] and the United States (~14%) [Reference Richardson17], which may result from variable population demographics [Reference Richardson17], and the early identification and early treatment of patients at high risk of respiratory failure in Jiangsu province [Reference Sun11].
Our study found that age was associated with respiratory failure. This is consistent with previous studies reporting that middle-aged and elderly patients with COVID-19 were susceptible to respiratory failure [Reference Liu4, Reference Bartoletti7], and is also similar to the results for patients with severe acute respiratory syndrome (SARS) [Reference Lee22–Reference Tsui24]. A previous study showed that the older population have a higher incidence of comorbidities and hence possible poorer immune response to COVID-19 and poorer clinical outcomes [Reference Luo25]. This study also found that patients with respiratory failure tended to have more comorbidities including hypertension, coronary heart disease and diabetes.
Increased respiratory rate as a risk factor of respiratory failure identified in this study is simple and convenient to apply for clinicians in their practice. This result is consistent with previous studies [Reference Bartoletti7, Reference Murphy26].
In our study, reduced lymphocyte count at admission was strongly associated with respiratory failure. Previous studies also show that the lymphocyte count of COVID-19 patients admitted to the ICU continued to decline [Reference Wang27, Reference Yang28], and the reduced lymphocyte count was a risk factor for respiratory failure [Reference Bartoletti7]. Reduced lymphocytes may be part of the pathogenesis of COVID-19 [Reference Rothan and Byrareddy29, Reference Chan30]. Therefore, more attention needs to be paid to patients whose lymphocyte count decline more severely. Laboratory abnormalities including lymphopenia have been previously reported in severe cases of other respiratory viral diseases, including SARS, Middle East respiratory syndrome (MERS) and influenza [Reference Lee22, Reference Ko31–Reference Chu34].
All radiologic images collected in our study were CT which have the advantage of high-resolution transversal imaging and accurate display of the extent and range of lung lesions. The potential measurement bias and misclassification bias resulted from visual assessment of CT images were controlled by double assessments of CT images by two independent radiologists with more than 5 years’ experience in pulmonary imaging, and by discussion and reaching agreements on different assessment results. This study showed that the pulmonary opacity score, one of the radiologic features, was strongly and independently associated with respiratory failure in patients with COVID-19, indicating that the more severe abnormality of lung function is an important factor to identify patients at high risks of respiratory failure. Former studies demonstrate that CT lung lesions can predict death and ICU admission in patients with COVID-19 [Reference Ruch13, Reference Davarpanah35]. On the other hand, this study shows that between patients with and without respiratory failure, there were no significant differences in the distribution of lung lesions (both were more likely to involve the outer third of lungs), lesion density (both were more likely to have less than 20% of consolidation of lungs) and lesion boundary definition (both were more likely to have the ill-defined border of lungs).
In Jiangsu province with adequate medical material and human resources, all COVID-19 patients with respiratory failure received ICU monitoring. In comparison, due to limited resources at disease outbreak in Spain, some comorbid patients with respiratory failure requiring mechanical ventilation and ICU treatment did not receive these treatments which were reserved for non-comorbid patients, hence a large number of patients with respiratory failure finally died [Reference Yarza36]. For similar settings with limited resources, the experience in Jiangsu province may be beneficial to COVID-19 patients with respiratory failure: in Jiangsu, a large proportion of patients with respiratory failure conducted prone position (30%) and received HFNC (40%) or NIV (~50%), which help reduce the further use of IMV (~10%) and mortality (no death occurred and all patients were discharged at the end of the study). This confirms the findings from several previous small studies [Reference Patel37, Reference Sivaloganathan38].
This study addressed several limitations of previous studies by (1) adding radiologic features and several epidemiologic, clinical and laboratory features into analysis to reduce residual confounding; (2) using multiple imputation method to provide unbiased estimates of levels and risk factors and (3) conducting a sensitivity analysis to confirm that results from multivariate logistic regression analysis, from which our main conclusion was drawn, were insensitivity to missing data and hence were robust and credible. This study included 625 patients from 24 hospitals, i.e. nearly all patients in Jiangsu province, China (a province with a population of 80 million; only six cases were excluded due to missing medical records), to make the findings from this study be subject to less selection bias and be more generalisable to populations in similar settings.
Our study has several limitations. Firstly, this is a retrospective observational study and its results may be subject to measurement bias and information bias, and some unobserved confounders (e.g. obesity, gene cluster) [Reference Rottoli16, Reference Ellinghaus39]. Secondly, some data on laboratory test were missing and hence fewer laboratory parameters were included in the analysis. Thirdly, we were unable to analyse the impact of medical management on respiratory failure including supportive treatments and medication, because of the chronological order and treatment information collected before, when, or after respiratory failure occurred.
Conclusions
In conclusion, this large cohort study based on a representative sample of 625 patients with COVID-19 shows that the rate of respiratory failure in Jiangsu province, China (9%), was similar to the national level in China, but much lower than in Wuhan city (the epicentre of COVID-19 pandemic in China) and some other countries. The study has also identified four independent risk factors of respiratory failure in patients with COVID-19 including older age, increased respiratory rate, decreased lymphocyte count and greater pulmonary opacity score at admission. For successful control of mortality related to COVID-19, patients with COVID-19 having these risk factors need to be intensively managed during hospitalisation.
Author contributions
SJ, YY, DW, YW and HL conceived and designed the study. YW, HL, SJ and SL contributed to the literature search and had full access to all of the data and took responsibility for the integrity of the data. SL, SJ, YY, DW, YW and HL contributed to data collection, quality checks and data management. DW, YW, HL, YY, SJ and SL contributed to data analysis and results presentation. YW, HL, SL, KM, DW, YY and SJ were responsible for results interpretation. YW, HL, SL, KM, DW, YY and SJ contributed in drafting the paper. All authors provided critical review of the paper and approved the final draft for publication.
Financial support
This work was supported, in part, by the Ministry of Science and Technology of the People's Republic of China [grant number 2020YFC0843700 67].
Conflict of interest
None.
Ethical standards
The Ethics Committee of Zhongda Hospital Affiliated to Southeast University approved this study (2020ZDSYLL013–P01 and 2020ZDSYLL019–P01) and waived patient informed consent due to the retrospective study design.
Data
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.