Corticosteroid therapy for coronavirus disease 2019-related acute respiratory distress syndrome: a cohort study with propensity score analysis

Background The impact of corticosteroid therapy on outcomes of patients with coronavirus disease 2019 (COVID-19) is highly controversial. We aimed to compare the risk of death between COVID-19-related ARDS patients with corticosteroid treatment and those without. Methods In this single-center retrospective observational study, patients with ARDS caused by COVID-19 between January 20, 2020, and February 24, 2020, were enrolled. The primary outcome was 60-day in-hospital death. The exposure was prescribed systemic corticosteroids or not. Time-dependent Cox regression models were used to calculate hazard ratios (HRs) and 95% confidence intervals (CIs) for 60-day in-hospital mortality. Results A total of 382 patients [60.7 ± 14.1 years old (mean ± SD), 61.3% males] were analyzed. The median of sequential organ failure assessment (SOFA) score was 2.0 (IQR 2.0–3.0). Of these cases, 94 (24.6%) patients had invasive mechanical ventilation. The number of patients received systemic corticosteroids was 226 (59.2%), and 156 (40.8%) received standard treatment. The maximum dose of corticosteroids was 80.0 (IQR 40.0–80.0) mg equivalent methylprednisolone per day, and duration of corticosteroid treatment was 7.0 (4.0–12.0) days in total. In Cox regression analysis using corticosteroid treatment as a time-varying variable, corticosteroid treatment was associated with a significant reduction in risk of in-hospital death within 60 days after adjusting for age, sex, SOFA score at hospital admission, propensity score of corticosteroid treatment, comorbidities, antiviral treatment, and respiratory supports (HR 0.42; 95% CI 0.21, 0.85; p = 0.0160). Corticosteroids were not associated with delayed viral RNA clearance in our cohort. Conclusion In this clinical practice setting, low-dose corticosteroid treatment was associated with reduced risk of in-hospital death within 60 days in COVID-19 patients who developed ARDS.


Background
The World Health Organization (WHO) declared that the outbreak of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) constitutes a pandemic [1]. Since the first confirmed case on January 22, 2020, the virus has been emerged in 216 countries. More than 14,765,000 laboratory-confirmed cases were reported, with an average mortality approaching 4.1% as of July 22, 2020 [2]. The spread of SARS-CoV-2 has led to serious socioeconomic consequences worldwide.
Currently, there is no specific treatment or vaccine for coronavirus disease 2019 . Up to 29% of the COVID-19 patients developed acute respiratory distress syndrome (ARDS) [3][4][5] as a consequence of cytokine storm. ARDS was the major cause of morbidity [6]. Adjunctive corticosteroids may be theoretically beneficial [7], and has been widely used by clinicians to suppression of hyperinflammation in COVID-19 patients, especially those with critical illness [3,5,8]. However, there was comprehensive controversy on its efficacy [9,10], due to the results of observational studies that showed corticosteroid treatment was associated with increased mortality and nosocomial infections for influenza and delayed virus clearance for severe acute respiratory syndrome (SARS), Middle East respiratory syndrome (MERS) [11].
In the early pandemic of COVID-19, lower mortality was reported among the critically ill subgroup of SARS patients treated with corticosteroids in a retrospective study [12]. Since then, evidence is growing that corticosteroid treatment is beneficial for some COVID-19 patients. The Randomized Evaluation of COVID-19 Therapy (RECOVERY) trial from the UK reported reduced mortality in patients treated with oral or intravenous dexamethasone 6 mg/d for up to 10 days [13]. The efficacy was significant only in severe patients who receiving oxygen or invasive mechanical ventilation. A prospective meta-analysis of 7 RCTs of critically ill patients with COVID-19 also showed the association between systemic corticosteroids and lower all-cause mortality [14]. Other two randomized controlled trials (RCTs) did not show benefit on mortality from intravenous dexamethasone or methylprednisolone treatment, but corticosteroid treatment increased ventilator-free days and number of days alive in patients with moderate to severe ARDS [15,16]. A meta-analysis of corticosteroids on the mortality of ARDS showed beneficial effects of corticosteroids on short-term mortality may be counteracted by the delayed onset of adverse effects [17], such as secondary infection due to immunosuppression and altered tissue repair. A follow-up period of 28 days in these randomized trials may underestimate the late adverse effects of corticosteroids on all-cause mortality. In addition, most of the aforementioned trials studied dexamethasone, and evidence for methylprednisolone was limited.
In this observational study, we thoroughly examined the associations of corticosteroid treatment with 60-day in-hospital mortality among a population of COVID-19 patients who have developed ARDS.

Study design and patients
This single-center, retrospective, cohort study was conducted at the Jin Yin-tan Hospital, Wuhan, China. We identified all adult patients with confirmed COVID-19 according to WHO interim guidance [18] and patients were admitted between January 20, 2020, and February 24, 2020. Then, we identified those who developed ARDS according to the WHO definition for analysis [18]. Among these patients, 84 have been described previously by Wu et al. [8], and 91 participated in the open-label trial of Lopinavir-Ritonavir [19]. To avoid the influence of early mortality before cortical steroid presenting treatment efficacy, patients who died or discharged within 2 days on hospital admission were excluded. Other exclusion criteria were: (1) participating in any double-blind clinical trial, (2) under long-term corticosteroid therapy for at least 1 month as part of treatment for chronic underline diseases, or (3) could not provide valid medical history because of mental disease. The Jin Yin-tan Hospital Ethics Committee approved the study (No. KY-2020-44.01) and granted a waiver of informed consent from study participants.
In Jin Yin-tan hospital, systemic corticosteroids were considered if patients had progressive respiratory failure or laboratory findings indicated the presence of hyperinflammatory response. In patients receiving mechanical ventilation, preventive ventilation strategy of tidal volume 4-8 mL/kg of predicted body weight, inhale positive airway pressure < 30 cmH 2 O, and PEEP > 5 cmH 2 O was followed.

Data collection
Administration of corticosteroids was defined as systemic use (oral or intravenous) of corticosteroids, including methylprednisolone, dexamethasone, hydrocortisone, and prednisone. The primary outcome was Keywords: Corticosteroids, Coronavirus disease 2019, Severe acute respiratory syndrome coronavirus 2, Mortality, Propensity score, Methylprednisolone 60-day in-hospital mortality. Patients were followed to death or discharge from hospital up to 60 days since hospital admission (last clinical outcome was observed on March 21, 2020). The secondary outcomes were time to SARS-CoV-2 viral clearance since symptom onset. Data on demographics, medical history, laboratory findings, chest radiology, medication use, and clinical outcomes were extracted retrospectively from electronic medical records using a standardized data collection form. All data were checked independently by two physicians (DH and XC). From January 11, 2020, to monitor the clearance of viral RNA, SARS-CoV-2 RNA was tested using polymerase chain reaction from throat-swab specimens for every other day after clinical remission of symptoms, including fever, cough, and dyspnea [8]. Viral clearance was defined as two consecutive negative results. The definitions of ARDS and other diseases were described in eMethods (Additional file 1).

Statistical analysis
Baseline characteristics were compared between patients with and without corticosteroid treatment. Data were reported as percentage for categorical variables and as mean ± standard deviation (SD) or median with interquartile range (IQR, 25-75%) for continuous variables. Categorical variables were compared by Fisher's exact test or Pearson chi-square test, as appropriate, and continuous variables were compared by Mann-Whitney U test or Student's t test.

Propensity score adjustment
To reduce the effect of steroids treatment bias and potential confounding factors, we performed propensity score analysis [20] to adjust the differences in baseline characteristics. For each patient, a propensity score indicating the likelihood of receiving systemic corticosteroid treatment was calculated by a logistic regression model. The model included 10 pre-selected baseline variables based on the clinical guidelines from National Health Commission of China, which recommended corticosteroids for patients with progressive respiratory failure and hyperinflammation response. Specifically, SpO 2 /FiO 2 ratio and respiratory rate were included for indicating severe respiratory failure; temperature, heart rate, SOFA score, blood lymphocyte count, blood neutrophil count, and level of CRP at hospital admission were included for indicating systemic inflammatory response syndrome; age and sex were included as basic characteristics of each patient. The outcome variable was whether or not the patient received corticosteroid therapy in the current hospital stay. Goodness of fit was evaluated by the c-statistic and the Hosmer-Lemeshow test.

Cox proportional-hazard regression model
The effect of corticosteroid treatment on risk of 60-day in-hospital all-cause death was analyzed using a series of Cox proportional-hazard regression models. First, we constructed a univariable Cox regression model on hospital death by 60 days since hospital admission with corticosteroid treatment treated as a time-varying covariate. Then, we constructed a multivariable Cox model of 60-day hospital death with corticosteroid treatment as time-varying covariate and incorporated the individual propensity score into the model as a covariable to calculate the propensity adjusted hazard ratio (HR). In the final model, the effects of corticosteroids on 60-day in-hospital death were adjusted for propensity score of corticosteroid treatment, as well as the following pre-selected covariates: age, sex, sequential organ failure assessment (SOFA) score at hospital admission, comorbidities (diabetes, hypertension, coronary artery disease, chronic pulmonary disease, chronic renal or liver disease, solid malignant tumor, hematologic malignancy, and immunosuppressive status), antiviral treatment (Lopinavir-Ritonavir, oseltamivir, and ganciclovir), and respiratory supports (high-flow oscillation oxygen, noninvasive mechanical ventilation, and invasive mechanical ventilation) in hospital [9].

Sensitivity analysis
Several sensitivity analyses were performed to assess the robustness of our findings. To test whether the findings were influenced by the time point of baseline variables used in propensity score analysis, Cox proportional-hazard regression models were repeated in sensitivity analyses: (1) comparing patients receiving corticosteroids within 2 days after hospital admission versus no corticosteroids; (2) comparing patients receiving corticosteroids versus no corticosteroids using ARDS onset date as baseline, where values of SpO 2 /FiO 2 ratio, respiratory rate, temperature, heart rate, respiratory rate, SOFA score, blood lymphocyte count, blood neutrophil count, and level of CRP at ARDS onset were used; (3) comparing patients receiving corticosteroids within 2 days after ARDS diagnosis versus no corticosteroids with ARDS onset as baseline. To test whether the findings might be influenced by ARDS definition, we conducted survival analysis using the same model among patients diagnosed with ARDS by Berlin definition [21]. The difference between Berlin Definition and WHO definition was that the latter included patients with SpO 2 /FiO 2 ≤ 315 when PaO 2 is not available, and did not restrict to ventilated patients.

Viral shedding and inflammation
Differences in the time to SARS-CoV-2 RNA clearance were analyzed using Cox proportional-hazard regression adjusted for the same covariables with corticosteroid therapy as a time-varying covariate. Patients died without viral shedding or discharged alive before they had two consecutive negative SARS-CoV-2 RNA tests were censored.
The associations between c-reactive protein with corticosteroids treatment after ARDS onset were analyzed using the interaction between corticosteroids and days after ARDS onset based on a linear regression model.
Results were analyzed with SAS (version 9.4, SAS Institute, Cary, NC). Unadjusted and adjusted hazard ratios and their 95% confidence intervals (Cis) were reported. Two-sided P values less than 0.05 were considered statistically significant.

Results
A total of 1147 patients with COVID-19 were screened for the study. Forty patients were excluded for participating in any double-blind clinical trial (n = 15), death or discharge from the hospital within 2 days after hospital admission (n = 21), underwent long-term corticosteroid therapy for chronic kidney disease or rheumatic disease (n = 3), or no valid medical history provided because of mental disease (n = 1). From 1107 patients remained, 382 patients were identified as ARDS (Additional file 1: Fig.  S1).
Baseline characteristics of the ARDS patients at hospital admission by receiving systemic corticosteroid treatment are shown in Table 1. In the entire cohort, the mean age was 60.7 ± 14.1 years, and 234 (61.3%) patients were male. 147 (38.5%) were treated with NIMV, 94 (24.6%) with IMV, and 11 (2.9%) with ECMO. All but one of the patients reached end point of decease (53.1%) or discharge (46.6%) during the follow-up period of 60 days. The median duration of follow-up was 12.0 (IQR 7.0-18.0) days. A total of 226 (59.2%) ARDS patients had a prescription of systemic corticosteroids. Corticosteroids were more likely prescribed to the younger (p = 0.0077) and males (p = 0.0135). Corticosteroids group had lower lymphocyte count and higher levels of CRP and lactate dehydrogenase at hospital admission than non-corticosteroids group, indicating a propensity in prescribing corticosteroids to patients with more severe immune dysfunction and inflammatory response ( Table 1). The 60-day hospital death in patients who ever used corticosteroids was higher than the patients who did not use corticosteroids Among patients prescribed corticosteroids, methylprednisolone was the most frequently administered corticosteroids (213/226, 94.2%) ( Table 2) In the logistic regression model generating propensity score, the pre-selected variables most closely correlated with prescription of systemic corticosteroids included age, blood lymphocyte count, heart rate and CRP (eTable 2, Additional file 1). The multivariable regression model of propensity for corticosteroid treatment had area under the receiver operating characteristic curve (ROC) of 0.71.
In sensitivity analysis, narrowing to patients meet the Berlin definition of ARDS did not alter the association between corticosteroids and lower risk of death (HR 0.43; 95% CI 0.21, 0.88; p = 0.0208). The HR of corticosteroids on risk of death was constant when using hospital admission as baseline and excluding patients received corticosteroids 2 days after hospital admission (HR 0.37; 95% CI 0.18, 0.76; p = 0.0072). When using ARDS onset date as baseline, the associations between corticosteroids and the risk of death were also significant ( Blood CRP level decreased among corticosteroids group on the first 4 days after ARDS onset (Fig. 2),

Table 3 Estimated effects of corticosteroid treatment on 60-day mortality in patients with ARDS associated with COVID-19
All of the models assessed the effects of corticosteroids as a time-varying covariate ARDS acute respiratory distress syndrome, FIO 2 fraction of inspired oxygen, SOFA sequential organ failure assessment, SpO 2 pulse oxygen saturation a Adjusted for age, sex, SOFA score at hospital admission, propensity score of corticosteroid treatment, comorbidities (diabetes, hypertension, chronic pulmonary disease, chronic renal or liver disease, solid malignant tumor, hematologic malignancy, and immunosuppressive status), antiviral treatment (Lopinavir-Ritonavir, oseltamivir, and ganciclovir), and respiratory supports (high-flow oscillation oxygen, noninvasive mechanical ventilation, and invasive mechanical ventilation). Propensity score was calculated by a non-parsimonious logistic regression model that included: age, sex, SOFA score, temperature, respiratory rate, SpO 2 /FiO 2 ratio, blood lymphocyte count, blood neutrophil count, and level of c-reactive protein at hospital admission b Using values of SpO 2 /FiO 2 ratio, respiratory rate, temperature, heart rate, respiratory rate, SOFA score, blood lymphocyte count, blood neutrophil count, and level of CRP at ARDS onset while an increase was found in non-corticosteroids group. CRP levels were significantly lower in corticosteroids-treated group after 2 days of ARDS onset (p for interaction = 0.0434).

Discussion
In this observational study, prescription of low-to-moderate dose systemic corticosteroids was associated with lower risk of 60-day in-hospital death among COVID-19 patients who developed ARDS. The efficiency of corticosteroids was further supported by the reduction of CRP, as the marker for suppressed systemic inflammation responses. No associations between corticosteroid treatment with viral shedding were found in our study.
Our study demonstrates the association between corticosteroid treatment with long-term (60 days) risk of death in severe COVID-19 patients. It is biologically plausible that suppression of inflammatory response by corticosteroids may be beneficial for patients with ARDS, which was caused by dysregulated systemic inflammation [3][4][5] and proved the main cause of death. The RECOV-ERY trial and a metanalysis of ongoing RCTs showed reduced 28-day mortality and longer ventilation-free days in patients with corticosteroid treatment [13]. We used 60-day in-hospital death as primary endpoint. To our knowledge, it was longer than previously reported RCTs. The results showed significant association between corticosteroids and risk of death, which further provide evidence on the long-term benefit of corticosteroids. It was compatible with studies that indicated corticosteroid treatment was not associated with increases in secondary infections in COVID-19-related ARDS patients [15,16].
Type, dosage, and duration of corticosteroids therapy were fundamental variables of corticosteroid treatment regimens. Different from the most of the published RCTs using dexamethasone, most of the patients in our study were treated with methylprednisolone. Methylprednisolone is a rapid onset glucocorticoid with shorter half-life and less mineralocorticoid effects than dexamethasone, which indicate shorter effects on systemic immunity and preventing corticosteroids-related fluid retention. Meduri et al. [22] firstly promoted the early use of corticosteroids in ARDS. They found that methylprednisolone reduced the duration of mechanical ventilation, ICU stay, and ICU mortality in early severe ARDS patients. In our study, a similar maximum dose (equivalent to methylprednisolone of 1-2 mg/kg) of corticosteroids was used, which was also close to the dose in RECOVERY trial and recommended by Society of Critical Care Medicine (SCCM) and European Society of Intensive Care Medicine (ESICM) for ARDS patients [23,24]. Our results were in line with the RCTs of COVID-19 [13,14] and previous studies of other ARDS patients, which showed low-dose corticosteroid treatment (equivalent to methylprednisolone of 1-2 mg/kg) accelerates the resolution of ARDS [23,25,26], indicating low-dose methylprednisolone as an alternative to dexamethasone in COVID-19-related ARDS. Higher dose may increase risks of Cox regression model with corticosteroid treatment was time-varying variable, adjusting for age, sex, SOFA score at hospital admission, propensity score of corticosteroid treatment, comorbidities (diabetes, hypertension, chronic pulmonary disease, chronic renal or liver disease, solid malignant tumor, hematologic malignancy, and immunosuppressive status), antiviral treatment (Lopinavir-Ritonavir, oseltamivir, and ganciclovir), and respiratory supports (high-flow oscillation oxygen, noninvasive mechanical ventilation, and invasive mechanical ventilation). Propensity score was calculated by a non-parsimonious logistic regression model that included: age, sex, SOFA score, temperature, respiratory rate, SpO 2 /FiO 2 ratio, blood lymphocyte count, blood neutrophil count, and level of c-reactive protein at hospital admission. ARDS, acute respiratory distress syndrome; SOFA, sequential organ failure assessment immune-suppression and corticosteroid-induced complications [27][28][29]. In this cohort, tapering strategy was performed as has been suggested by the guidelines for the ARDS-related corticosteroid insufficiency (CIRCI) to reduce deterioration from the development of a reconstituted inflammatory response and febrile response. Of note, a randomized trial included mild COVID-19 patients using a similar dose and shorter course of methylprednisolone (0.5 mg/kg twice daily for 5 days) than in the regimens of our study found no benefit on mortality. More research is needed to determine the best duration of corticosteroid therapy. A recent analysis of four trials showed prolonged corticosteroids therapy reduced mortality [25]. However, chronic side effects of corticosteroids including secondary infection and osteoporosis may occur in prolonged course of treatment.
Delayed virus clearance was reported in corticosteroidtreated patients with both SARS, MERS, and influenza [30][31][32], which was a major concern for the immune suppressive effects of corticosteroids, albeit its uncertain clinical relevance. We found no difference in viral shedding duration from symptom onset between corticosteroid and non-corticosteroid groups, which may explain the heterogeneity in efficacy of corticosteroids between COVID-19 and other virus infections. Of note, positive SARS-CoV-2 test results have been reported after two consecutive negative results [33]. Viral tests of throat swabs were not monitored after two consecutive negative tests in our cohort. More evidence is needed for assessing the effects of corticosteroids on clearance of SARS-CoV-2 RNA.
We used rigorous statistical method to control for survival and indication bias. Survivors-treated bias exists in observational studies that assess exposure after the start of follow-up, where only patients survived long enough had an opportunity to receive the intervention. Therefore, the patients died early are more likely to be misclassified to the no-treatment group, leading to overestimation of the effectiveness of medicine [34]. This study was specifically designed to address survivors-treated bias of corticosteroid treatment, by using a time-dependent variable for corticosteroids initiation to define corticosteroids group and non-corticosteroids group [35]. In addition, there was a propensity of clinicians to give corticosteroids to patients who were critically ill in non-randomized clinical condition. The imbalance in baseline characteristics may introduce confounders in comparison of mortality between corticosteroids and non-corticosteroids group. In this cohort, lymphopenia and elevation of CRP and lactate dehydrogenase levels were more severe in patients who received corticosteroids therapy. Propensity score is a validated method to account for baseline confounding and control selection bias in this case [36]. In this study, we performed a rigorous propensity adjustment analysis accounting for the baseline variables related to propensity of corticosteroid treatment. These added to the strength of our results that found associations between corticosteroid treatment and risk of death.
Our study had some limitations. First, unlike randomized controlled trials, the selection bias and potential confounding effects might exist. We used propensity analysis rather than standard multivariable analysis to rigorously adjust for selection bias, and time-dependent model to avoid survivors-treated bias. Nonetheless, only measured factors were controlled for due to the nature of observational study design. Second, secondary infections were not monitored in this study, because microbiological culture results needed for definite diagnosis of secondary infection were possibly affected by antibiotic treatment the patients received simultaneously. To include the delayed effects of secondary infections on mortality, a longer follow-up period of 60 days was used. Third, this study was single center and patients were sicker and transferred from other hospital, so might lacking of generality. Forth, some of the patients were treated with Lopinavir-Ritonavir; it might be a confounder because the efficacy of Lopinavir-Ritonavir in COVID-19 was unclear. Fifth, our cohort was collected in the early outbreak of COVID-19; thus, the mortality was relatively higher than other studies, which limited generalization of our results.

Conclusion
Our findings suggest administration of low dose of corticosteroids might reduce the risk of death in COVID-19 patients who developed ARDS.