Introduction

Patients with coronavirus disease 2019 (COVID-19) have occasionally developed severe pneumonia, and some of these patients progress to life-threatening respiratory failure, acute respiratory distress syndrome (ARDS) and multiple organ failure1,2. Although the mechanisms of COVID-19–induced lung injury and multiple organ failure are still being elucidated, patients with severe COVID-19 are reported to have higher serum cytokine levels than those with mild to moderate COVID-19, suggesting that a “cytokine storm” may be one of the etiological factors2,3,4,5. This condition is associated with rapid deterioration in the severe acute respiratory syndrome coronavirus (SARS-CoV)-1 and Middle East respiratory syndrome coronavirus (MERS-CoV)6,7,8. If the COVID-19–induced lung injury worsens to the degree that invasive mechanical ventilation or extracorporeal membrane oxygen therapy (ECMO) is required, the mortality is very high9,10. Therefore, appropriate anti-inflammatory therapy to suppress the cytokine storm is considered crucial to prevent progression to irreversible ARDS and multiple organ failure11,12.

Corticosteroid therapy is expected not only to suppress the cytokine storm but also to prevent the progression to pulmonary fibrosis associated with COVID-19 pneumonia, and has been widely used to treat previously prevalent SARS-CoV-1 and MERS-CoV13,14. Recent results of several randomized trials of corticosteroids against COVID-19 have been reported15,16,17,18, in which the therapy reduced the 28-day mortality and increased the number of ventilator-free days in critically ill patients with COVID-1915,16. In light of these results, the latest World Health Organization (WHO) guidance recommends corticosteroids for severe and critical patients19. On the other hands, most of the randomized trials reported so far did not include non-severe patients who did not require invasive or noninvasive ventilation. Only in the Randomized Evaluation of COVid-19 thERapY (RECOVERY) trial has the efficacy of corticosteroids for non-severe patients been validated, and corticosteroids failed to show a survival benefit for patients not receiving respiratory support, and might even be harmful15. Based on this result alone, WHO guidance suggested not to use corticosteroids for the treatment of non-severe patients.

Therefore, the usefulness and necessity of corticosteroids for COVID-19 remains controversial, especially for the patients who do not require invasive or noninvasive ventilation. Moreover, many factors remain unclear regarding the appropriate use of corticosteroids for COVID-19, such as initial dose, administration period, and timing of initiation. In this context, this multicenter, retrospective, propensity score–matched study was launched to evaluate the efficacy of systemic corticosteroid administration for hospitalized patients with COVID-19 ranging in the degree of severity from mild to critically-ill disease. In addition, various subgroup analyses were performed to examine in detail the appropriate use of corticosteroids for COVID-19.

Methods

Study design and participants

This multicenter, retrospective study was conducted at 30 institutions in Japan. The study enrolled all consecutive patients who met the following inclusion criteria: (1) SARS-CoV-2 infection confirmed by polymerase chain reaction (PCR) test; (2) diagnosed between January 23–April 30, 2020; (3) required hospitalization for COVID-19; and (4) did not require home oxygen therapy before infection with COVID-19. Clinical and laboratory data were retrieved from patient medical records. Clinical outcomes for COVID-19 patients who received systemic corticosteroids (corticosteroid group) were compared with those who did not receive this therapy (non-corticosteroid group), after adjusting for propensity scores. The case registration period was from May 1–June 30, 2020.

Propensity score matching

The method of propensity score matching was used to minimize the bias due to confounding factors, assuming that an imbalance in patient background between the corticosteroid and non-corticosteroid groups may exist. The propensity score for each patient was calculated as a probability from a logistic regression model, including all covariates that were considered clinically important and had an impact on the patient's prognosis: (1) gender; (2) age; (3) body mass index; (4) smoking history; (5) comorbid hypertension; (6) comorbid diabetes mellitus; (7) time from symptom onset to admission; (8) score of 7-point ordinal scale on Day 1; (9) oxygen saturation (SpO2)/fraction of inspired oxygen (FiO2) on Day 1; (10) dyspnea; (11) pneumonia on initial chest X-ray or computed tomography (CT); (12) C-reactive protein (CRP); (13) concomitant use of favipiravir, and (14) concomitant use of any non-steroidal treatment for COVID-19. In mild to severe patients who do not require invasive/noninvasive ventilation, FiO2 was estimated from the delivery system and flow rate using a commonly used conversion table.

Endpoints

The primary endpoint was the odds ratio (OR) for improvement of the score on a 7-point ordinal scale on Day 15, with the first day of hospitalization as Day 1. The ordinal scale is an assessment of the clinical status on a given day. The 7-point scale is as follows: (1) death; (2) hospitalized, on invasive mechanical ventilation or ECMO; (3) hospitalized, on noninvasive positive pressure ventilation (NIPPV) or high-flow nasal cannula (HFNC); (4) hospitalized, requiring low flow supplemental oxygen; (5) hospitalized, not requiring supplemental oxygen, requiring ongoing medical care; (6) hospitalized, not requiring supplemental oxygen, no longer required ongoing medical care; and (7) discharged/not hospitalized.

The key secondary endpoints were as follows: (1) time to PCR negativity of the swab solution; (2) duration of fever; (3) percentage of improvement in radiological findings; (4) time to improvement in radiological findings; (5) proportion of patients requiring invasive mechanical ventilation with tracheal intubation/ECMO; (6) time to requiring invasive mechanical ventilation with tracheal intubation; (7) duration of invasive mechanical ventilation with tracheal intubation; (8) hospitalization period, and (9) survival period.

Statistical analysis

In the primary analysis, ordinal variables were compared between groups using a proportional odds model. In the secondary and exploratory analysis, time to event was estimated using the Kaplan–Meier method. The Cox proportional hazards model was used to calculate the hazard ratio (HR) and its 95% confidence interval (CI) for the treatment effect between groups. Categorical variables were presented as numbers (percentages), and compared using chi square test or Fisher exact test. Normally distributed continuous variables were presented as mean and standard deviation (SD), and compared using t test. Continuous variables related to time were presented as median (interquartile ranges) and compared using t test. A p value < 0.05 was considered statistically significant. All statistical analyses were performed using statistical software package SAS (version 9.4, SAS Institute).

Ethics approval and participant consent

This study was performed in accordance with the Declaration of Helsinki. The study protocol was approved by the Ethics Committee of the Kanagawa Cardiovascular and Respiratory Center (approval date: April 21, 2020, approved number: KCRC-20–0004), and the Institutional Review Board or Ethics Committee of other participating facilities. According to the Ethical Guidelines for Medical Research on Human Subjects in Japan, this research falls under the category of research, which does not involve intervention and does not use samples obtained from the human body. The need for patient consent was waived because this was a retrospective study and anonymity was secured. For this reason, the Institutional Review Boards or Ethics Committees of all participating facilities approved that we applied opt-out method by publishing this study on either the participating facility's website or on a bulletin board. This study was registered in the University hospital Medical Information Network Clinical Trials Registry on April 21, 2020 (ID: UMIN000040211).

Results

Patient disposition

Of 1141 consecutive hospitalized patients with COVID-19 registered by June 30, 2020, 49 patients were excluded from this study based on these criteria: (1) diagnosis after May 2020 (25 patients); (2) duplicate registrations from 2 facilities (11 patients); (3) no clinical status information for Day 1 due to transfer from another hospital (9 patients); (4) no need to be hospitalized (2 patients), (5) PCR testing only performed on spinal fluid (1 patient), and (6) negative results on PCR, but clinical diagnosis (1 patient). Thus, 1092 patients were included in the final analysis (Fig. 1).

Figure 1
figure 1

Patient Disposition. COVID-19, coronavirus disease 2019; PCR, polymerase chain reaction.

Clinical characteristics and prognosis of the 1092 patients analyzed are shown in Supplemental Table 1 and 2. The mortality was 2.1% on Day 14, and 3.8% on Day 28. Of the 235 patients who received corticosteroids for COVID-19, 163 (69.4%) received early corticosteroids within 3 days after admission. The remaining 72 patients (30.6%) who started corticosteroid > 4 days after admission had a greater decline in SpO2/FiO2 from admission to just before corticosteroid initiation and a significantly worse score on the 7-point ordinal scale on Day 15 compared with the 163 early-treatment patients (Supplemental Table 3 and 4). We considered that the more patients in the corticosteroid group who started delayed corticosteroids after their respiratory status had already deteriorated, the more difficult it would be to assess the primary endpoint (improvement in clinical status on Day 15) in comparison with the non-corticosteroid group, even using propensity score matching; moreover, early corticosteroids have been reported to be effective against COVID-19.20 Therefore, the 163 patients who received early corticosteroids within 3 days after admission were designated as the unmatched corticosteroid group, and the 857 patients who did not receive corticosteroids for COVID-19 were categorized as the unmatched non-corticosteroid group. After propensity score matching, 118 patients were assigned to either the corticosteroid and non-corticosteroid groups.

Baseline characteristics before/after propensity score matching

The distribution of the patients’ baseline characteristics according to corticosteroid exposure is shown in Table 1, both in the unmatched and matched samples. The unmatched samples included a significantly higher number of male patients and those who were older age, had a higher weight and body mass index, and had more comorbidities (hypertension and diabetes) in the corticosteroid versus non-corticosteroid group. In addition, clinical and laboratory data for the corticosteroid versus non-corticosteroid group showed significantly poorer clinical status in a 7-point ordinal scale on Day 1, lower SpO2/FiO2, higher rates of fever and dyspnea, higher CRP concentrations, and lower lymphocyte counts.

Table 1 Baseline characteristics before and after propensity score matching.

Standardized mean differences for each covariate before and after propensity score matching are shown in Fig. 2. The differences between corticosteroid and pretreatment variables were attenuated in the matched versus unmatched samples for propensity score. In fact, baseline characteristics were well balanced between the corticosteroid and non-corticosteroid groups after propensity score matching (Table 1). Regarding the baseline score on the 7-point ordinal scale in the matched samples, 4 was the most common score for both the corticosteroid and non-corticosteroid groups (44.9% vs. 50.8%), followed by a score of 5 (40.7% vs. 39.0%).

Figure 2
figure 2

Standardized Mean Difference Before and After Propensity Score Matching. * Covariates used for propensity score matching. COVID-19, coronavirus disease 2019; CT, computed tomography.

Regarding the specific COVID-19 treatment administered both in the propensity score–unmatched and matched corticosteroid groups, nearly 90% of the corticosteroids administered for COVID-19 were methylprednisolone, with a median starting dose of 80 mg/day and a mean administration period of 11.0 days (Table2).

Table 2 Treatment for coronavirus disease 2019.

Primary outcome

The odds of improvement in a 7-point ordinal scale on Day 15 were significantly lower in the corticosteroid versus non-corticosteroid group (OR, 0.611; 95% CI, 0.388–0.962; p = 0.034) (Table 3). In critically ill patients with a baseline 7-point ordinal score of 2 or 3, the clinical status on Day 15 was similar in both groups (OR, 0.953; 95% CI, 0.215–4.224; p = 0.950). In contrast, for patients with mild to severe disease with a baseline score of 4 or 5, the odds of improvement were lower in the corticosteroid group than in the non-corticosteroid group.

Table 3 Primary outcome.

Key secondary outcomes

The key secondary outcomes are shown in Table 4. No significant differences were observed between the two groups with respect to time to PCR negativity or duration of hospitalization. The duration of fever was significantly longer in the corticosteroid group (HR, 0.746; 95% CI, 0.560–0.994; p = 0.045). The time to improvement in radiological findings was significantly shorter in the corticosteroid versus non-corticosteroid group (HR, 1.758; 95% CI, 1.323–2.337; p < 0.001), regardless of baseline score of 7-point ordinal scale (Fig. 3). The number of patients requiring invasive mechanical ventilation was higher in the corticosteroid versus non-corticosteroid group (33.9% vs. 17.8%; p = 0.0072), with median time from admission to tracheal intubation of 2 days for both groups (Supplemental Fig. 1). The duration of invasive mechanical ventilation was shorter in the corticosteroid versus non-corticosteroid group (HR, 1.466; 95% CI, 0.841–2.554; p = 0.177) (Fig. 4A). Mortality on Day 28 tended to be higher in the corticosteroid versus non-corticosteroid group (10.2% vs. 4.2%; p = 0.1289), and the HR was 2.417 (95% CI, 0.868–6.733; p = 0.091) (Supplemental Fig. 2A).

Table 4 Secondary outcomes.
Figure 3
figure 3

Time to Improvement in Radiological Findings. Kaplan–Meier curves for the time to improvement in radiological findings. Cox proportional hazards model was used to calculate the hazard ratio and its 95% confidence interval for the treatment effect between groups. IQR, interquartile ranges; CI, confidence interval.

Figure 4
figure 4

Duration of Invasive Mechanical Ventilation. (A) Kaplan–Meier curves for the duration of invasive mechanical ventilation comparing corticosteroid group and non-corticosteroid group. (B) Kaplan–Meier curves for the duration of invasive mechanical ventilation comparing the pulse/semi-pulse group (initial dose ≥ 250 mg/day) and the standard dose group (initial dose < 250 mg/day) among patients receiving methylprednisolone. Cox proportional hazards model was used to calculate the hazard ratio and its 95% confidence interval for the treatment effect between groups. IQR, interquartile ranges; CI, confidence interval.

Subgroup analysis based on initial dose, administration period and timing of corticosteroids

Subgroup analysis was performed based on initial dose, administration period, and timing of corticosteroids (Table 5). Of the 106 patients who received methylprednisolone, the duration of invasive mechanical ventilation was significantly shorter in the pulse/semi-pulse group (initial dose ≥ 250 mg/day) than in the standard dose group (initial dose < 250 mg/day) (median, 8 days vs. 15 days; HR, 2.831; 95% CI, 1.347–5.950; p = 0.006) (Fig. 4B). In the patients receiving corticosteroids for ≤ 10 days, the time to PCR negativity of the swab solution tended to be shorter (HR, 1.437; 95% CI, 0.968–2.132; p = 0.072) compared with the patients receiving corticosteroids for > 11 days.

Table 5 Subgroup analysis in the propensity-score matched corticosteroid group.

Safety outcome

Safety outcomes for both the corticosteroid and non-corticosteroid groups were also analyzed. Results showed no significant difference in the frequency of thromboembolism between the corticosteroid and non-corticosteroid groups (2.5% vs. 3.4%).

Discussion

This study demonstrated the following 3 important clinical observations. First, corticosteroids did not lead to avoidance of tracheal intubation or lower mortality in patients with mild to severe COVID-19. Second, for the critically ill patients, corticosteroid therapy reduced not only the time to improvement in radiological findings, but also the duration of invasive mechanical ventilation. Third, methylprednisolone pulse/semi-pulse therapy significantly shortened the duration of invasive mechanical ventilation compared with the standard dose.

In both the corticosteroid and non-corticosteroid groups, > 85% of patients did not require invasive or noninvasive ventilation at baseline. Because data are lacking on the benefit of corticosteroids for patients with mild to severe COVD-19 at baseline, the results of this study deserve to be noteworthy. In mild to severe patients with a baseline 7-point ordinal score of 4 or 5, the clinical status on Day 15 tended to be worse in the corticosteroid group than in the non-corticosteroid group. Administration of corticosteroids to patients with mild to severe COVD-19 was expected to prevent the progression to critical conditions that would require ECMO or invasive mechanical ventilation. However, a high proportion of patients in the corticosteroid group required invasive mechanical ventilation, despite the fact that corticosteroids reduced the time to improvement in radiological findings and suggested some benefit. This finding may have been due to the fact that the time from admission to tracheal intubation in patients who required invasive mechanical ventilation (median 2 days) was shorter than the time to improvement in imaging findings (median 8 days). Because corticosteroid treatment takes a certain amount of time to show benefit, it may not be expected to improve short-term outcomes, such as the avoidance of tracheal intubation.

This study also suggested that corticosteroids may have a negative impact on survival as assessed by Day 14, Day 28 mortality, and the HR in mild to severe patients with a baseline 7-point ordinal score of 4 or 5 (Supplemental Fig. 2C and 2D). Similarly, the RECOVERY study suggested that dexamethasone may rather worsen the prognosis among patients who were not receiving any respiratory support at randomization15. Although there is concern that corticosteroid administration within 7 days of onset may inhibit antibody production, the median time from symptom onset to admission in the corticosteroid group in this study was 8.4 days. One possible explanation is that adverse effects of corticosteroids may have affected the prognosis. In the present study, no increase in thromboembolism was observed with corticosteroids. Although it cannot be ruled out that impaired hyperglycemic control and secondary infections may have an impact on prognosis, these data were not collected in this study and are a limitation for the application of study findings. However, it is questionable whether such events really have a significant impact on prognosis. Among the studies of COVID-19 patients requiring hospitalization and treatment, the patient population included in this study had a clearly better prognosis with a lower mortality than that described in many previous reports21,22. Although the cause is unclear, the mortality is low not only in Japan, but also in most of the countries in East and Southeast Asia23. In light of these findings, it may be difficult to draw definitive conclusions about the survival endpoint based on the results of this study, and caution should be used when comparing and interpreting data from this study with data from previous studies of populations in Europe and in North and South America.

Meanwhile, for the critically ill patients with a baseline ordinal score of 2 or 3, clinical status on Day 15 assessed by the 7-point ordinal scale was similar between the corticosteroid and non-corticosteroid groups. The present study showed 2 positive effects of corticosteroids as expected for critically ill COVID-19 patients. First, the reduction in the time to improvement in radiological findings may suggest the effectiveness of corticosteroids. Second, corticosteroids tended to reduce the duration of invasive mechanical ventilation, which is consistent with previous reports from randomized clinical trials. In the CoDEX trial, dexamethasone significantly increased the number of ventilator-free days in patients with moderate to severe ARDS who required intubation and ventilation17. Although this study had a relatively small proportion of critical ill cases at baseline, early administration of corticosteroids to patients who develop severe respiratory failure requiring invasive or noninvasive ventilation may be beneficial.

It is also noteworthy that a subgroup analysis in this study showed that methylprednisolone pulse/semi-pulse therapy shortened the duration of mechanical intubation compared with the usual dose regimen. A small, single-blind, randomized, controlled, clinical trial in Iran reported that methylprednisolone pulse therapy (intravenous injection, 250 mg/day for 3 days) reduced the time of clinical improvement and discharge from the hospital or death in severe hospitalized patients compared to the standard of care24. However, to date, no previous reports have examined the differences in efficacy of different starting dose of corticosteroids. In contrast, subgroup analysis by administration period in this study suggested that prolonged corticosteroid administration over 11 days may prolong the time to PCR negativity. This result suggests that prolonged corticosteroid administration may delay the elimination of the virus from the body. Therefore, high-dose, short-term corticosteroid therapy should be considered in critically ill patients with COVID-19 pneumonia.

As a limitation of this study, even among the groups matched for propensity score, the corticosteroid group may still have included more rapidly deteriorating patients than the non-corticosteroid group. In fact, a higher rate of invasive mechanical ventilation (many of which cases are introduced within a few days), longer-lasting fevers despite corticosteroid therapy, and worsening of SpO2/FiO2 even in a short period of time from baseline to just before corticosteroid initiation were observed in the matched corticosteroid group. As a possible reason, although the impact on severity and prognosis of COVID-19 patients differed greatly among the covariates used in the propensity score matching method in this study, all of these covariates were treated as having equal weight. In addition, the presence or absence of each variable implies the equivalence between multiple variables in terms of severity, but this is not appropriate because being dyspneic is clearly not the same as being diabetic, at least for COVID-19. Therefore, the methodology used for matching was one of the limitations of this study. The biomarker to identify a rapidly deteriorating population among COVID-19 is not well established and may have been difficult, at least with the factors used for propensity score matching in this study. Although serum ferritin concentration could not be used for matching in this study because of the large number of deficiencies, this variable may be useful as a biomarker (serum ferritin concentration was measured for only 26 patients in the corticosteroid group [mean 1534.9 μg/dL] and 51 in the non-corticosteroid group [mean 774.1 μg/dL]). As for radiological findings, it was difficult to make a uniform and detailed evaluation because some patients did not have CT scans, so we only collected information on the presence of pneumonia as determined by the investigator. Detailed analysis of baseline CT images may also be a useful biomarker based on reports that the degree of extension of lung opacities and lung volume loss on CT had an impact on prognosis25,26,27,28. Because the clinical data in this study were collected retrospectively from the electronic medical records of each 30 participating institutions, and because we were afraid to overburden the investigators in the current situation where COVID-19 is still raging, subjective symptoms such as dyspnea and comorbidities such as diabetes were assessed only for presence or absence and not for severity of illness. As an additional limitation, it is necessary to discuss whether the endpoint using an ordinal scale was valid as a primary endpoint. None of the previously reported studies evaluating the efficacy of corticosteroids for COVID-19 have met the endpoint on an ordinal scale17. Many infectious disease studies have relatively short-term measures set as primary endpoints, but the primary endpoint for the study of COVID-19 may need to be established with a longer-term perspective.

Conclusions

Corticosteroids for hospitalized patients with COVID-19 did not improve clinical status on Day 15. However, corticosteroids reduced not only the time to improvement in radiological findings in all patients regardless of disease severity, but also the duration of invasive mechanical ventilation in the critically ill patients. Methylprednisolone pulse/semi-pulse therapy significantly shortened the duration of invasive mechanical ventilation compared with the standard dose.