key: cord-0865610-nz3rwwbb authors: Cusacovich, Ivan; Aparisi, Álvaro; Marcos, Miguel; Ybarra-Falcón, Cristina; Iglesias-Echevarria, Carolina; Lopez-Veloso, Maria; Barraza-Vengoechea, Julio; Dueñas, Carlos; Juarros Martínez, Santiago Antonio; Rodríguez-Alonso, Beatriz; Martín-Oterino, José-Ángel; Montero-Baladia, Miguel; Moralejo, Leticia; Andaluz-Ojeda, David; Gonzalez-Fuentes, Roberto title: Corticosteroid Pulses for Hospitalized Patients with COVID-19: Effects on Mortality date: 2021-03-12 journal: Mediators Inflamm DOI: 10.1155/2021/6637227 sha: 00dc9435e19f3c8bb02d896af6aa51ab904b46c1 doc_id: 865610 cord_uid: nz3rwwbb OBJECTIVES: To assess the influence of corticosteroid pulses on 60-day mortality in hospitalized patients with severe COVID-19. METHODS: We designed a multicenter retrospective cohort study in three teaching hospitals of Castilla y León, Spain (865,096 people). We selected patients with confirmed COVID-19 and lung involvement with a pO2/FiO2<300, excluding those exposed to immunosuppressors before or during hospitalization, patients terminally ill at admission, or those who died in the first 24 hours. We performed a propensity score matching (PSM) adjusting covariates that modify the probability of being treated. Then, we used a Cox regression model in the PSM group to consider factors affecting mortality. RESULTS: From 2933 patients, 257 fulfilled the inclusion and exclusion criteria. 124 patients were on corticosteroid pulses (250 mg of methylprednisolone for three days), and 133 were not. 30.3% (37/122) of patients died in the corticosteroid pulse group and 42.9% (57/133) in the nonexposed cohort. These differences (12.6%, 95% CI [8·54-16.65]) were statically significant (log-rank 4.72, p = 0, 03). We performed PSM using the exact method. Mortality differences remained in the PSM group (log-rank 5.31, p = 0.021) and were still significant after a Cox regression model (HR for corticosteroid pulses 0.561; p = 0.039). CONCLUSIONS: This study provides evidence about treatment with corticosteroid pulses in severe COVID-19 that might significantly reduce mortality. Strict inclusion and exclusion criteria with that selection process set a reliable frame to compare mortality in both the exposed and nonexposed groups. In December 2019, a new betacoronavirus called SARS-CoV-2 induced severe bilateral pneumonia similar to severe acute respiratory syndrome (SARS), described in 2003. This coronavirus disease (COVID-19) had lower mortality than SARS-CoV-1 infection but higher infective capacity. The epidemic began in Wuhan, mainland China, but in a few months became pandemic. Spain was one of the world's most affected countries, especially in Madrid, Catalonia, and Castilla y León regions [1] . We analyzed patients with COVID-19 admitted between March 12th and May 20th to three tertiary teaching hospitals in Castilla y León, Spain: Hospital Clínico Universitario de Valladolid (HCUV), Hospital Universitario de Salamanca (HUSA), and Hospital Universitario de Burgos (HUBU). The three hospitals cover all hospital admissions in a geographical area corresponding to 865 096 people. The treating team decided on the prescription of all drugs without any intervention from investigators. We obtained the local ethics committee (CEIC) permission to perform the study. Informed consent was obtained. We designed a retrospective cohort study and compared a cohort of patients exposed to corticosteroid pulses and an unexposed one. 2.1. Data Source. We analyzed paper and electronic records in all hospitals. We recorded variables related to clinical outcomes and corticosteroids exposure (supplementary material section A (available here)). We included patients older than 18 years, testing positive on SARS-CoV-2 PCR (nasopharyngeal or oropharyngeal swab specimens). Patients with positive ELISA serology and consistent clinical symptoms were also considered confirmed cases. All included patients had a significant lung involvement, defined as a pO2/FiO2<300, maintained for 24 hours or repeated for three days. We measured pO2/FiO2 in arterial gasometry or estimated it from pulse oximetry data (nonlinear estimate model) [41, 42] . We excluded patients receiving classic immunosuppressors or cytokine blockers (as cyclosporine, tocilizumab, or anakinra). Concomitant drugs allowed were hydroxychloroquine, azithromycin, remdesivir, lopinavir/ritonavir, and colchicine. We excluded patients who died in the first 24 hours of admission. Patients on corticosteroid treatment in a different regimen than the one described in this study were excluded. We also excluded pregnant women, terminally ill patients, and patients under a limitation of therapeutic efforts during the first 24 hours of admission. We considered exposure to corticosteroid pulses if administered at a daily dose of 125 to 500 mg of intravenous methylprednisolone for two to five days. We did not include patients with repeated corticosteroid pulses nor treatments longer than five days. About timing, we considered corticosteroid pulses in the ±3 days, respecting the inclusion criterion date. The primary endpoint was 60-day mortality in exposed versus nonexposed patients. Secondary endpoints were 30-day mortality, intensive care unit (ICU) admission, in-hospital stay, viral shedding until negative PCR, and serious adverse events, including infections. 2.5. Statistical Analysis. We expressed continuous variables with the median and interquartile range (mean and standard deviation if they had normal distribution). We used chisquare to compare qualitative variables and the t-test (if normal distribution) or the Mann Withney test to compare two quantitative variables. We performed the Kolmogorov-Smirnov test to examine normal distribution. We performed a propensity score matching to balance the difference of covariates related to exposure to corticosteroid pulses. To select suitable covariates to control, we set biologically plausible variables related to the probability of being treated with corticosteroids pulses (propensity score). First, we analyzed variables associated with the propensity score in univariate analysis. Variables found significant were dichotomized, and then, we performed a binary logistic regression to evaluate independent variables associated with the propensity score. We chose three matching methods (propensity score matching) to preprocess the sample: the nearest neighbor, the nearest neighbor with a caliper (at a distance of 0.05, 0.1, 0.2, and 0.3), and exact matching. We performed the propensity score matching using the R software, with the MatchIt and Cobalt libraries. We checked the balance of the propensity score matching through the "difference of means" to ensure that the distribution of covariates was similar in the treated and control group, and we picked the best-matched model. Once we completed the propensity score matching, we performed a Cox proportional hazard regression analysis to evaluate mortality and intensive care admission. From 2933 patients in our cohort, 257 fulfilled the inclusion and exclusion criteria. 767 fulfilled the inclusion criterion, and 546 had any exclusion criteria (see Figure 1 ). We diagnosed with COVID-19 in 243 patients based on SARS-CoV-2 PCR in the nasopharyngeal or oropharyngeal swabbing and 14 patients based on positive serology with compatible symptoms. 124 patients were on corticosteroid pulses, and 133 were not. The most used corticosteroid pulses dose was 250 mg daily for three days (92%,114/124). 3.2. Propensity Score Matching. We calculated the propensity score (probability of being treated with corticosteroid pulses) in each participant from a binary logistic regression. Variables statically significant in the binary logistic regression were as follows: epidemiological week, presence of bilateral infiltrates or not, Center in Castilla y León, Ferritin, and COVID-gram score (see supplementary material section B (available here)). After performing the propensity score, we tried several preprocessing methods for matching (see Methods), and we selected the exact matching method as it got the minimum difference between groups with the minimum sample loss (28 controls and 2 treated patients). The difference of means was zero in the treated versus the controls because we used the exact matching method (Figure 2 and supplementary material section B (available here)). After the propensity score matching, the sample consisted of 207 patients (119 treated and 88 controls). Comparing patients exposed to corticosteroid pulses and the nonexposed ones, we found that age and comorbidities were similar in both groups, without significant differences. The COVID-gram score [43] was 155.8 in the bolus group and 152.3 in the control group, even decreasing these differences after matching. The classic Charlson score was significantly higher (0.7 points) in the control group (p = 0:012). These differences disappeared after dichotomizing the variable in the matched group (≤2 or >2 more than 2 points, p = 0:171) (see Table 1 ). LDH and ferritin at admission were higher in the pulse group, but these differences disappeared after matching. Peak ferritin and peak LDH during hospitalization were significantly higher in the pulse group, even after matching (see Table 1 ). Concomitant treatments with colchicine, interferon beta-1b, lopinavir/ritonavir, and azithromycin were more common in the corticosteroid pulse group than in the controls, both before and after matching (see Table 1 ). We did not find differences in pO2/FiO2 between the pulses and control group (p = 0:183 in all participants and p = 0:69 in the matched group), but bilateral lung infiltrates were more frequent in the corticosteroid pulse group 3.4.1. Primary Endpoint. Ninety-four patients died during the 60 days after admission, representing 36.9% (94/255) of the sample. 30.3% (37/122) of patients died in the corticosteroid pulse group and 42.9% (57/133) in the nonexposed cohort. These differences (12.6%) were statically significant in the Kaplan Meier curve (log-rank 4.72, p = 0:03) (see Figure 3 ). We carried out a propensity score matching and calculated the mortality in the matched group. The corticosteroid pulse group had a 60 days mortality of 29.6% (34/115), while mortality in the control group was 44.3% (39/88). Again, these differences were statically significant (log-rank 5.31, p = 0:021) (see Figure 3 and Table 2 ). We performed a multivariate analysis using a Cox regression model in the propensity score matching group, after dichotomizing variables, and those independently related to 60-day mortality were as follows: corticosteroid pulses (HR 0.561, p = 0:039), age older than 80 years (HR 7.3, p < 0:001 ), CPR at admission higher than 200 mg/dL (HR 3.35, p < 0:001), neutrophil/lymphocyte index > 7:4 (HR 2.15, p = 0:010), Charlson index higher than 2 points (HR 2.15, p = 0:018), and LDH at admission > 372 UI/mL (HR 2.29, p = 0:008) (see Figure 4 ). In the equation, we did not include other variables (like pO2/FiO2, lung infiltrates, or D-dimer) that were not significant in the multivariate model, although important in the univariate analysis. We studied the possible association of colchicine, azithromycin, and lopinavir/ritonavir on mortality, as those treatments were more frequently used in the corticosteroid pulse group. The three treatments were associated with lower mortality in the Kaplan-Meier curves (p = 0:008 for colchicine, p = 0:032 for lopinavir/ritonavir, and p < 0:001 for azithromycin). These differences disappear after adjusting for age higher than 80 years in both three drugs (p = 0:1 for colchicine, p = 0:794 for lopinavir/ritonavir, and p = 0:378 for azithromycin). Patients on corticosteroid pulses had a higher, but not significant, anticoagulation rate (86.9% vs. 83.9%, p = 0:49). Patients on pulses also used significantly lower prophylactic dose heparin and higher intermediate and anticoagulant dose. Still, we did not find in our study an association between the absence of anticoagulation and mortality (logrank 0.9, p = 0:338) (see supplementary material Section D). ) ). After using the same Cox regression model used in the 60-days mortality analysis, corticosteroid pulses remain a protective factor for 30-days mortality (p = 0:049) (see supplementary material section C1 (available here) and Figure 3 ). Data are shown as the median (IQR) or n (%). * Significant (p < 0, 05) difference between methylprednisolone and usual care population in both groups, the overall population, and propensity score match. It was calculated using the χ 2 test or Kruskal-Wallis test as appropriate. PaO2/FiO2 = ratio of arterial oxygen partial pressure to fractional inspired oxygen. Forty-nine patients were admitted to ICU during this period. ICU admission was 18.5% (23/124) in the corticosteroid pulses cohort and 19.5% (26/133) in the nonexposed group (p = 0:838). We found similar results in the propensity score matching group (ICU admission 16.2% in the corticosteroid group and 23.9% in the control group, p = 0:173). The time from hospital admission to ICU admission was zero to four days, and 83.6% of patients were admitted in the first 24 hours. Mortality among the ICU admitted patients was 17.4% (4/23) in the exposed group and 61.5% (16/26) in the nonexposed group. This difference was statistically significant for 60-day mortality (log-rank = 9:7, p = 0:002). This difference remains in the propensity score-matched group (15.8% vs. 61.9%, p = 0:003). After adjusting for several variables (peak LDH, peak CRP, number of comorbidities, D-dimer at admission, SaO2/FiO2, and age) in a Cox regression, those differences in mortality were not significant (see supplementary material section C2 (available here)). The ICU average stay was 19.2 and 21.67 days for both the exposed and nonexposed cohorts, respectively (p = 0:750). The in-hospital median stay was 12 [7.25-19 .75] days in the treated cohort and 8 [5] [6] [7] [8] [9] [10] [11] [12] [13] [14] [15] days in the nonexposed group (p < 0:001). These differences remain in the matched group (p = 0:001). However, differences were not statistically significant if we analyzed them in the survivor's group (difference of means 1.2 days, p = 0:619). Viral shedding until negative PCR was shorter (but not significant, p = 0:279) in the corticosteroid cohort (25.02 days in the corticosteroid group and 30.65 days in the nonexposed). We reported serious adverse events in 17 patients in the exposed cohort and 15 patients in the control group (p = 0:133). Hemorrhage happened in 9 patients, without difference between groups (p = 0:053). In-hospital infections were not higher in the corticosteroid bolus group (29/124) than in the nonexposed cohort (32/129) (p = 0:792). We found 27 patients with pulmonary embolism by CT scan (8 in the exposed group and 19 in nonexposed, p = 0:066). We also reported two ischemic strokes and one acute myocardial infarction. Corticosteroid pulses have been widely used in Spain, especially in the Castilla y León region, but not in other countries to treat COVID-19. The rationale of its use is to stop the systemic inflammation process [3] that develop in some patients with severe COVID-19. Some studies [44] [45] [46] have described the positive effects of corticosteroid pulses on mortality in patients with severe COVID-19. We found a significant improvement in survival in patients treated with corticosteroid pulses. We designed a retrospective cohort study to confirm this statement, which is the main limitation of our work. As there is no randomization, unknown confounders might be unattended. We carried out a multicenter study with three teaching hospitals in the Castilla y León region in Spain. This fact is one of the main strengths of the study. On the one hand, we summarized various treatment protocols in each center, showing a wider specter of treatment options for severe COVID-19. This protocol variety determines different probabilities of being treated with corticosteroids in each center and enables us to adjust them in the latter propensity score matching. It also considers different hospital admission criteria and different extrahospital resources that may change the baseline features of hospitalized patients. On the other hand, it represents all hospital admissions in a geographic area in Castilla y León with more than 865 000 people that have similar epidemiological features. They also had the same timing of lockdown and the same mobility restrictions over time. We used strict inclusion and exclusion criteria to avoid mixing the effects of other immunosuppressive treatments in mortality. They were also useful to find an adequate patient profile who, a priori, should benefit on an antiinflammatory treatment as corticosteroid pulses. We selected patients, at the inclusion time, with the onset of an acute respiratory distress syndrome. To limit potential biases, we performed a propensity score matching (PSM) using the exact method. Thus, we obtained a more homogeneous sample with baseline features that were similar in both groups. Both the exposed and nonexposed cohorts in the propensity score matching group had comparable age, a similar pretest probability of dying (through the COVID-gram score), probability of being treated with corticosteroids, and uniform comorbidities. We again found, in this PSM group, the same mortality decrease in the corticosteroid pulse arm. Then, to adjust other possible mortality causes, we performed a Cox regression multivariate model on the PSM group, once more finding a significant protective role of corticosteroid pulses in mortality. Altogether, joining these strict inclusion and exclusion criteria with all this selection process sets a reliable frame to compare mortality in both the exposed and nonexposed groups. Thus, corticosteroid pulses might be a good option for the treatment of severe COVID-19, as they have been shown to be effective in reducing mortality in our cohort and they are inexpensive and highly available worldwide. Our results can only extrapolate to patients with severe COVID-19, with a pO2/FiO2 lower than 300, not exposed to any kind of immunosuppression, and in the absence of a terminally ill situation at admission. Some recent studies have confirmed that oral or intravenous low-dose corticosteroids positively affect mortality [29] . There was a decrease in mortality between 3.1 and 12.1% (in the ICU admitted group), lower than the 12.6% of global mortality reduction (even higher in the ICU subset) that we found. We hypothesized that the effect of cor-ticosteroid pulses might be higher than the low-dose corticosteroids because they act in different pathways (genomic vs. nongenomic) and behave, in fact, as different drugs [37] . Our study is not powered to compare both low-dose and pulse corticosteroid treatments, so we cannot assure this statement. Future studies must investigate this topic. Pulse corticosteroids did not reduce the ICU admission rate in our study. Most patients were moved to the ICU in the first 24 hours of hospitalization (83.6%), so we understand that those patients were critically ill at admission time, requiring ICU in any case. The in-hospital stay was significantly longer in the pulses corticosteroid arm, but those differences disappear in the survivor's group. Thus, this difference in the hospital stay was due to higher survival in the corticosteroid pulse group. The rate of adverse events and serious adverse events declared was similar in both groups. In the same way, inhospital infection and viral shedding time were similar in both groups, but some of these adverse events and infections might be underreported. The study was not powered to detect these adverse events because they were not always reported in the medical record in all patients. In conclusion, this study provides evidence about the treatment with corticosteroid pulses in severe COVID-19 that might significantly reduce mortality. This data must be confirmed in prospective randomized studies. Acute respiratory distress syndrome COVID-19: Coronavirus disease 2019 ICU: Intensive care unit PSM: Propensity score matching SARS-CoV-2: Severe acute respiratory syndrome coronavirus 2 PCR: Polymerase chain reaction IFN-γ: Interferon-γ TNF-α: Tumor necrosis factor α LDH: Lactate dehydrogenase CRP: C-reactive protein. We recorded all data in a database. These data could be requested if there is an important reason. Bullet Points. (i) Corticosteroid pulses can improve mortality in severe COVID-19 patients. (ii) The Kaplan Meier curve with the propensity score matching and multivariate analysis in this group is reliable. (iii) There are no substantial changes in ICU admission and hospital stay. We sent a preprint to MedRxiv, with the following doi:10 .1101/2020.09.30.20204719. Centro de coordinación y alertas sanitarias Predictors of mortality for patients with COVID-19 pneumonia caused by SARS-CoV-2: a prospective cohort study COVID-19: consider cytokine storm syndromes and immunosuppression Interleukin-1 blockade with high-dose anakinra in patients with COVID-19, acute respiratory distress syndrome, and hyperinflammation: a retrospective cohort study Silencing the cytokine storm: the use of intravenous anakinra in haemophagocytic lymphohistiocytosis or macrophage activation syndrome COVID-19: an update on the epidemiological, clinical, preventive and therapeutic evidence and guidelines of integrative Chinese-Western medicine for the management of 2019 novel coronavirus disease The Italian coronavirus disease 2019 outbreak: recommendations from clinical practice Ineffectiveness of high-dose methylprednisolone in preventing parenchymal lung injury and improving mortality in patients with septic shock Clinical evidence does not support corticosteroid treatment for 2019-nCoV lung injury Efficacy and safety of corticosteroids for persistent acute respiratory distress syndrome Corticosteroid treatment in critically ill patients with severe influenza pneumonia: a propensity score matching study Adjuvant corticosteroid treatment in adults with influenza a (H7N9) viral pneumonia The effect of corticosteroids on mortality of patients with influenza pneumonia: a systematic review and meta-analysis Corticosteroids in the prevention and treatment of acute respiratory distress syndrome (ARDS) in adults: meta-analysis Effect of prolonged methylprednisolone therapy in unresolving acute respiratory distress syndrome: a randomized controlled trial Established acute respiratory distress syndrome: benefit of corticosteroid rescue therapy Early methylprednisolone treatment for septic syndrome and the adult respiratory distress syndrome High-Dose corticosteroids in patients with the adult respiratory distress syndrome Lung function and organ dysfunctions in 178 patients requiring mechanical ventilation during the 2009 influenza a (H1N1) pandemic Effect of low-to-moderate-dose corticosteroids on mortality of hospitalized adolescents and adults with influenza A(H1N1)pdm09 viral pneumonia Early corticosteroids in severe influenza A/H1N1 pneumonia and acute respiratory distress syndrome Corticosteroid therapy in patients with primary viral pneumonia due to pandemic (H1N1) 2009 influenza Effect of immunomodulatory therapies in patients with pandemic influenza A (H1N1) 2009 complicated by pneumonia Methylprednisolone infusion in early severe ards: results of a randomized controlled trial Early steroid therapy for respiratory failure Dexamethasone treatment for the acute respiratory distress syndrome: a multicentre, randomised controlled trial Risk factors associated with acute respiratory distress syndrome and death in patients with coronavirus disease 2019 pneumonia in Wuhan, China Association between administration of systemic corticosteroids and mortality among critically ill patients with COVID-19 Dexamethasone in hospitalized patients with Covid-19 -preliminary report Endogenous glucocorticoids control neutrophil mobilization from bone marrow to blood and tissues in non-inflammatory conditions Glucocorticoids inhibit apoptosis of human neutrophils Effect of dexamethasone on reactive oxygen species generation by leukocytes and plasma interleukin-10 concentrations: a pharmacodynamic study Effect of fluticasone propionate on neutrophil chemotaxis, superoxide generation, and extracellular proteolytic activity in vitro Glucocorticoids severely impair differentiation and antigen presenting function of dendritic cells despite upregulation of toll-like receptors Glucocorticoids affect human dendritic cell differentiation and maturation Effects of antirejection drugs on innate immune cells after kidney transplantation Genomic and nongenomic effects of glucocorticoids Effectiveness of prolonged glucocorticoid treatment in acute respiratory distress syndrome: the right drug, the right way? Khamashta M. Glucocorticoid use and abuse in SLE Pulse glucocorticoid therapy Nonlinear imputation of PaO2/FIO2 from SpO2/FIO2 among patients with acute respiratory distress syndrome Nonlinear imputation of PaO2/FIO2 from SpO2/FIO2 among mechanically ventilated patients in the ICU Development and validation of a clinical risk score to predict the occurrence of critical illness in hospitalized patients with COVID-19 Second week methyl-prednisolone pulses improve prognosis in patients with severe coronavirus disease 2019 pneumonia: an observational comparative study using routine care data intravenous methylprednisolone pulse as a treatment for hospitalised severe COVID-19 patients: results from a randomised controlled clinical trial Effectiveness of corticoid pulses in patients with cytokine storm syndrome induced by SARS-CoV-2 infection Roche granted the present study through the IECSCYL foundation. Ainhoa Cusacovich contributed to the paper translation. Section A: Variable Description and Coding in the Study. Section B: Construction of the Propensity Score and Matching Process. Section C: Secondary Endpoints. Section D: Anticoagulation and Survival Analysis. (Supplementary Materials) The authors declare that there are no conflicts of interest. All authors contributed to data collection. Ivan Cusacovich, Álvaro Aparisi, Miguel Marcos, David Andaluz-Ojeda, Roberto Gonzalez-Fuentes, and Carlos Dueñas contributed to the study design. Ivan Cusacovich and David Andaluz-Ojeda performed the statistical analysis. Ivan Cusacovich, Álvaro Aparisi, David Andaluz-Ojeda, Roberto Gonzalez-Fuentes, Miguel Marcos, and Carlos Dueñas contributed to writing the draft. David Andaluz-Ojeda and Roberto Gonzalez-Fuentes contributed equally to this work.