key: cord-0776070-65w2tgcq
authors: Li Bassi, Gianluigi; Gibbons, Kristen; Suen, Jacky Y.; Dalton, Heidi J.; White, Nicole; Corley, Amanda; Shrapnel, Sally; Hinton, Samuel; Forsyth, Simon; Laffey, John G.; Fan, Eddy; Fanning, Jonathon P.; Panigada, Mauro; Bartlett, Robert; Brodie, Daniel; Burrell, Aidan; Chiumello, Davide; Elhazmi, Alyaa; Esperatti, Mariano; Grasselli, Giacomo; Hodgson, Carol; Ichiba, Shingo; Luna, Carlos; Marwali, Eva; Merson, Laura; Murthy, Srinivas; Nichol, Alistair; Ogino, Mark; Pelosi, Paolo; Torres, Antoni; Ng, Pauline Yeung; Fraser, John F.
title: Early short course of neuromuscular blocking agents in patients with COVID-19 ARDS: a propensity score analysis
date: 2022-05-17
journal: Crit Care
DOI: 10.1186/s13054-022-03983-5
sha: cda71971c5681c25e742e1e055d4d2821a3a6ff6
doc_id: 776070
cord_uid: 65w2tgcq

BACKGROUND: The role of neuromuscular blocking agents (NMBAs) in coronavirus disease 2019 (COVID-19) acute respiratory distress syndrome (ARDS) is not fully elucidated. Therefore, we aimed to investigate in COVID-19 patients with moderate-to-severe ARDS the impact of early use of NMBAs on 90-day mortality, through propensity score (PS) matching analysis. METHODS: We analyzed a convenience sample of patients with COVID-19 and moderate-to-severe ARDS, admitted to 244 intensive care units within the COVID-19 Critical Care Consortium, from February 1, 2020, through October 31, 2021. Patients undergoing at least 2 days and up to 3 consecutive days of NMBAs (NMBA treatment), within 48 h from commencement of IMV were compared with subjects who did not receive NMBAs or only upon commencement of IMV (control). The primary objective in the PS-matched cohort was comparison between groups in 90-day in-hospital mortality, assessed through Cox proportional hazard modeling. Secondary objectives were comparisons in the numbers of ventilator-free days (VFD) between day 1 and day 28 and between day 1 and 90 through competing risk regression. RESULTS: Data from 1953 patients were included. After propensity score matching, 210 cases from each group were well matched. In the PS-matched cohort, mean (± SD) age was 60.3 ± 13.2 years and 296 (70.5%) were male and the most common comorbidities were hypertension (56.9%), obesity (41.1%), and diabetes (30.0%). The unadjusted hazard ratio (HR) for death at 90 days in the NMBA treatment vs control group was 1.12 (95% CI 0.79, 1.59, p = 0.534). After adjustment for smoking habit and critical therapeutic covariates, the HR was 1.07 (95% CI 0.72, 1.61, p = 0.729). At 28 days, VFD were 16 (IQR 0–25) and 25 (IQR 7–26) in the NMBA treatment and control groups, respectively (sub-hazard ratio 0.82, 95% CI 0.67, 1.00, p = 0.055). At 90 days, VFD were 77 (IQR 0–87) and 87 (IQR 0–88) (sub-hazard ratio 0.86 (95% CI 0.69, 1.07; p = 0.177). CONCLUSIONS: In patients with COVID-19 and moderate-to-severe ARDS, short course of NMBA treatment, applied early, did not significantly improve 90-day mortality and VFD. In the absence of definitive data from clinical trials, NMBAs should be indicated cautiously in this setting. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1186/s13054-022-03983-5.

Since early 2020, SARS-CoV-2 infections have placed tremendous burden on patients and international healthcare services [1] . A high proportion of diseased patients require hospitalization, and a small subset with severe coronavirus disease 2019 (COVID-19) become critically ill and require invasive mechanical ventilation (IMV) for life-threatening respiratory failure [2] [3] [4] [5] . High mortality has been reported in this subpopulation [4, [6] [7] [8] [9] , irrespective of survival benefits from established treatments, such as corticosteroids and IL-6 receptor antagonists [10, 11] .

Neuromuscular blocking agents (NMBAs) have been commonly used for acute respiratory distress syndrome (ARDS) [12, 13] to reduce patient-ventilator asynchrony and-in the context of severely damaged and poorly compliant lungs-to improve oxygenation, while minimizing the work of breathing and risks of barotrauma. International studies have confirmed the use of NMBAs in up to 26% of ARDS patients [14] . Yet, conflicting results were provided by large randomized clinical trials, such as the ACURASYS [15] and ROSE trials [16] , and indication, efficacy, and safety on the use of NMBAs in such patients remain uncertain.

The severe respiratory failure associated with COVID-19 has been described as a form of ARDS, and the potential benefits of supportive treatments have largely been extrapolated from evidence in non-COVID-19-related ARDS. Thus, many COVID-19 patients have been treated with NMBAs [17, 18] , and in some areas even shortage of these medications has been reported [19] , especially during the early phase of the pandemic. To the best of our knowledge, no international studies have clearly elucidated the effects of NMBAs on mortality in COVID-19 patients. Interestingly, in a multicenter observational study, Courcelle et al. [20] reported that the duration of NMBA treatment in this population was often longer than 48 h and not associated with shorter duration of IMV.

To further delineate the role of NMBA in ventilated COVID-19 patients, data extracted from the multicenter registry of the international COVID-19 Critical Care Consortium incorporating the ExtraCorporeal Membrane Oxygenation for 2019 novel Coronavirus Acute Respiratory Disease (COVID-19-CCC/ECMOCARD) were examined. Our hypothesis was that in patients with microbiologically confirmed COVID-19 and moderateto-severe ARDS, no significant difference in 90-day hospital mortality was to be found between populations who received or not early NMBAs. Thus, the primary goal of this study was to identify difference in 90-day mortality, through propensity score (PS)-adjusted analysis, between patients who received or not a short course of NMBAs, within 48 h from commencement of IMV.

This was a comparative study in which mechanically ventilated COVID-19 patients with moderate-to-severe ARDS, who received an early short course of NMBAs, were compared with patients who did not receive NMBA or underwent NMBA treatment only on the day of commencement of IMV to explore the impact on in-hospital mortality during a follow-up period up to 90 days. Of note, the STROBE (Strengthening the Reporting of Observational Studies in Epidemiology) guidelines were used to ensure the reporting of this observational study [21] .

We studied a convenience sampling of patients (≥ 16 years old) with microbiologically confirmed SARS-CoV-2 infection, via rapid nucleic acid amplification test (NAAT) or antigen tests, who were admitted to an intensive care unit (ICU) from February 1, 2020, through October, 31, 2021, at any of the COVID-19 Critical Care Consortiums participating hospitals. In addition, patients were selected if they presented with moderate-to-severe ARDS, defined by a ratio of the partial pressure of arterial oxygen to the fraction of inspired oxygen of < 150 mm Hg [16] , within 48 h from commencement of IMV.

We excluded patients with unrecorded date of commencement of IMV or who were transferred from other institutions, after commencement of IMV. In addition, we excluded from the primary analysis all patients with continuous use of NMBA for periods longer than 3 days, and intermittent use of NMBA within a 3-day period of treatment.

Data on demographics, comorbidities, clinical symptoms, and laboratory results were collected by clinical and research staff of the participating ICUs in an electronic case report form [22] . Details of respiratory and hemodynamic support, physiological variables, and laboratory results were collected daily, with worst daily values recorded preferentially. Duration of IMV, length of ICU and hospital stays, and hospital mortality were also recorded.

In line with previous therapeutic protocols applied in patients with moderate-to-severe ARDS [15, 16] , we aimed to appraise a short course of NMBAs, applied early during the course of IMV. Thus, patients were assigned to the following study groups.

Use of NMBA treatment for at least 48 h and up to 3 consecutive days, initiated early, within 48 h from commencement of IMV.

No use of NMBAs or administration only on the day IMV was commenced.

Of note, the time of first NMBAs administration was not recorded on the case report form; thus, the treatment group was designed to include up to 3 consecutive days of NMBAs exposure to ensure that at least 2 days of full treatment were achieved in most of the patients, even if treatment commenced during the night or a single NMBA dose was given upon endotracheal intubation. In addition, patients who were tracheally intubated and transferred from institutions not collaborating with the COVID-19-CCC/ECMOCARD were excluded from the analysis to avoid enrolment of subjects who might have received NMBA prior to study monitoring. Recording of NMBA use was censored at 28 days from ICU admission, or upon discharge from the ICU or death, whichever occurred first.

The primary objective was comparison between groups in 90-day in-hospital mortality, through PS matching analysis. Secondary objectives were the numbers of ventilator-free days (days since successful weaning from mechanical ventilation) between day 1 and day 28 and between day 1 and day 90 in the PS-matched population.

VFD was calculated as follows: VFDs = 0 if subject died in hospital within 28 or 90 days of mechanical ventilation; VFDs = 28 or 90 − x if successfully liberated from ventilation x days after initiation; VFDs = 0 if the subject was mechanically ventilated for more than 28 or 90 days.

We analyzed the COVID-19-CCC/ECMOCARD study dataset. COVID-19-CCC/ECMOCARD is an international, multicenter, observational cohort study ongoing in 354 hospitals across 54 countries (Additional file 1: Appendix). The full study protocol has been published elsewhere [22] . The COVID-19-CCC/ECMOCARD observational study was reviewed under the National Mutual Acceptance scheme by the Alfred Health Human Research Ethics Committee on February 27, 2020. The ethics committee certified that the study protocol met the requirements of the National Statement on Ethical Conduct in Human Research (2007), granting approval on March 2, 2020. In addition, in accordance with the Office of the Health Services Commissioners Statutory Guidelines on Research issued for the purposes of Health Privacy Principles, the Alfred Health Human Research Ethics Committee granted a waiver of consent for the collection, use, and disclosure of participants health and personal information. Subsequently, the protocol was approved by all international participating hospitals prior to data collection and waiver of consent was granted in all centers.

The COVID-19 Consortium collaborates with the International Severe Acute Respiratory and Emerging Infection Consortium (ISARIC) group [23] and their Short PeRiod IncideNce sTudy of Severe Acute Respiratory Infection (SPRINT-SARI) [24] project. De-identified patient data are recorded by data collectors at each collaborating site via the REDCap (Vanderbilt/NIH/NCATS UL1 TR000445 v.10.0.23) electronic data capture tool, using instances hosted at the University of Oxford (UK), University College Dublin (Ireland), Monash University (Australia), and The University of Queensland (Australia). A detailed data dictionary was provided to all sites to assist in data collection. In addition, biweekly drop-in data sessions were scheduled to assist with all queries that data collectors might have. Importantly, the database quality audits of the COVID-19 Consortium dataset are a continuing and intensive process encompassing (1) data cleaning rules; (2) checks for outliers; (3) filtering rules setup during the initial development of the case report form, which was periodically monitored/adjusted; and (4) data completeness checks. Finally, in the case any issue was detected during monitoring of data quality, or statistical analysis, these matters were pursued to address any data collection/process limitation in a timely manner, often including follow-up with the site that entered the data for value verification and correction where appropriate.

Categorical data are presented as frequency and percentage, while continuous data as mean and standard deviation (SD; normally distributed) and median and interquartile range (IQR; non-normally distributed). Normality of continuous data was reviewed via visual inspection of histograms. Bivariate comparisons between patients receiving NMBA treatment and controls were compared using Fishers exact test, Students unpaired t test, or the Mann-Whitney U test for categorical, normally, and non-normally distributed continuous variables, respectively.

Analyses were performed to examine the relationship between the primary outcome-mortality censored at 90 days post-commencement of IMV-and NMBA treatment. Patients entered the study analysis once they commenced IMV. As an imbalance in several characteristics at ICU admission between the study groups was noted, propensity score (PS) matching analysis was undertaken to balance these patient characteristics. Matching was undertaken, in preference to other techniques involving propensity scores, to best mimic the results of a randomized controlled trial. Baseline features considered potential variables for inclusion in the PS calculation were demographic and clinical characteristics with either a documented relationship with mortality or a baseline imbalance. Variables with more than 5% missing data were not considered. Based on the available literature on risk factors associated with mortality in severe COVID-19, the following covariates were included during PS modeling: age, sex, region (reference group: North America), time from symptom onset to hospital admission and baseline comorbidities (PaO 2 /FiO 2 within 48 h from commencement of IMV, hypertension, chronic cardiac disease, chronic kidney disease, obesity). Matched cohorts were constructed based on the logit of the PS using nearest-neighbor matching with a caliper width of 0.2 and replacement to promote better balance, and accounting for the estimation of the standard errors. Rubins B [25] (a summation measure of bias with value < 25% indicating adequate balance) and Rubins R (0.5-2.0 indicating appropriate balancing) statistics were initially calculated to assess the success of the matching (i.e., if the resulting matched cohort had balanced characteristics). Standardized mean differences for all covariates, before and after matching, were then estimated, with an absolute difference of 10% or greater considered indicative of imbalance [26, 27] .

Following PS matching, survival curves were generated for the matched sample to compare patients undergoing NMBA treatment against control. Finally, multivariable Cox regression models were constructed to assess the effect of NMBA treatment on outcome-incorporating robust estimators of variance to account for matched observations-as well as investigating effect modification on the association between mortality and NMBA by concomitant use of corticosteroids and respiratory failure severity. Covariates considered for inclusion in the multivariable Cox regression were laboratory results upon admission, treatments received and ventilator settings; covariates with known relationship with mortality in COVID-19 patients available in this dataset with less than 5% missing data were included in the final model. Where potential covariates were correlated, the most relevant clinical variable was chosen for inclusion in the model to avoid multicollinearity and the violation of model assumptions. Continuous covariates were standardized by subtracting the mean and dividing by the standard deviation to facilitate meaningful interpretation of the hazard ratios. Multiple imputation was not undertaken, and complete case analysis was performed [28] . Assumptions of the Cox regression model, particularly the proportional hazards assumption, were evaluated using test of the Schoenfeld residuals, as well as graphical exploration of the log-log plot; the proportional hazards assumption was met. Competing risk regression was performed to assess the secondary outcomes of VFDs at 28 and 90 days post-starting of IMV [29] .

Considering the observational nature of our report and potential limitations related to PS matching, we also balanced baseline patient characteristics by weighting each individual in the analysis, by the inverse probability of receiving exposure to NMBAs (inverse probability of treatment weighting [IPTW]) [30] . In addition, as the dataset did not report exact NMBA administration practice, the following sensitivity analyses were conducted to confirm the results of the primary analysis: (1) in patients with clinically suspected and microbiologically confirmed COVID-19, with the treatment group defined as per the primary analysis; and in the treatment group only including patients who received NMBAs (2) continuously for 2 days; (3) continuously for 3 days; and (4) continuously for more than 3 days.

Of note, a post hoc power calculation was undertaken, as a convenience sample was being used. In the unmatched cohort, assuming type I error of 0.05, as well as the rate of primary outcome and study group size as reported in the results, the analysis had 98% power, adequate to undertake and report on the analyses of interest. Data were analyzed using StataSE version 17.0 (StataCorp Pty Ltd., College Station, Texas). Any two-tailed p value less than 0.05 was considered significant. No correction was made for multiple comparisons.

There were 11,873 patients with microbiologically confirmed, or suspected, COVID-19 admitted to 244 hospitals between February 1, 2020, and October 31, 2021. After excluding 3244 patients because of lack of microbiological confirmation of SARS-CoV-2 infection and 1844 patients who were transferred from non-collaborating centers, or transferred out while still on mechanical ventilation, 4616 patients who were on IMV were included in the primary analysis ( Fig. 1 ). Figure 2 shows enrolment rate of those mechanically ventilated patients by date of ICU admission. A matched cohort based on propensity scores was generated, with 210 patients who underwent NMBA treatment and 210 patients as control group. Overall, baseline characteristics were balanced in the propensity score-matched cohort. The mean age was 60.3 years (standard deviation [SD] 13.2 years) ( Table 1) and 296 (70.5%) were male. The mean acute physiology and chronic health evaluation (APACHE) II score was 19.5 (SD 11.6; N = 117). In the NMBA treatment and control group, respiratory system compliance was 32.5 mL/ cmH 2 O (SD 12.6) and 33.9 (SD 12.2), respectively, and driving pressure was 22.7 cmH 2 O (SD 7.5) and 22.4 (SD 8.6). The most common comorbidities among patients in the PS-matched cohort were hypertension (239, 56.9%), obesity (172, 41.1%), and diabetes (125, 30.0%).

In 3 days. In the matched control group, patients never received NMBAs, not even upon the day IMV started. NMBA was provided on average 1 day after ICU admission (N = 207; IQR 0-2 days). The median length from commencement of IMV to initiation of NMBA therapy was 0 days (IQR 0-0 days).

In the PS-matched cohort, as summarized in Table 2 , PaO 2 /FiO 2 was 88.6 (SD 29.7) vs 86.0 (SD 30.7), in NMBA treatment and control group, respectively. Upon IMV commencement, patients who received NMBA treatment presented PaCO 2 of 51.0 mmHg (SD 13.7) vs 48.0 mmHg (SD 15.5) in those who did not. Positive endexpiratory pressure in those undergoing NMBA treatment or not was 12.8 cmH 2 O (SD 3.3) vs 11.9 cmH 2 O (SD 3.1), respectively. As for adjunctive therapies, extracorporeal membrane oxygenation, renal replacement therapy, and recruitment maneuvers were provided more frequently to patients who received NMBA treatment (Table 3) . Pneumothorax occurred in 21 (10.4%) and 19 (9.6%) of the patients undergoing NMBA treatment or not, respectively.

The Kaplan-Meier (Fig. 3A ) (log rank test p < 0.001) and the Cox regression model (unadjusted HR 1.78; 95% CI 1.38, 2.30; p < 0.001) confirmed that NMBA treatment was associated with increased mortality risk. The PS-matched cohort analysis comprised 210 patients in the NMBA treatment group vs 210 controls (using replacements). The summaries of balance for unmatched and matched critical parameters are depicted in Fig. 4 . Kaplan-Meier curves (Fig. 3B ) (log rank p = 0.537) showed no effect of NMBA treatment on mortality, which was also corroborated by the Cox regression model (unadjusted HR = 1.12, 95% CI 0.79, 1.59, p = 0.534). The lack of association with 90-day mortality was consistent after adjusting for smoking habit and critical therapeutic covariates (adjusted HR 1.07; 95% CI 0.72, 1.61, p = 0.729) ( Table 4 ). In subgroup analyses, PaO 2 /FiO 2 and the concomitant use of corticosteroids did not alter negative NMBA association with mortality.

Sensitivity analyses confirmed the lack of association between NMBA treatment and 90-day mortality, when inverse probability weighting analysis was applied (adjusted HR 1.28; 95% CI 0.89, 1.84; p = 0.187) and in Tidal volume (ml/PBW) mean (SD) 7. of 90-day mortality (adjusted HR 1.73, 95% CI 1.22, 2.37, p = 0.001). In this context, the reported median duration of NMBA treatment was 6 days (IQR 5-10).

In the PS-matched cohort at 28 days after 

To the best of our knowledge, this is the largest, international observational study of patients with COVID-19 requiring mechanical ventilation to assess the impact of short-course NMBA treatment, commenced early during IMV, on 90-day in-hospital mortality. We found that NMBA use was common, specifically at European sites, and frequently applied in patients who presented hypercapnic and required higher levels of PEEP upon commencement of IMV. PS-matching analysis confirmed that early NMBA treatment did not result in lower mortality, while post hoc sensitivity analysis found increased mortality risk when NMBAs use was extended beyond 3 days. In addition, at 28 and 90 days, there were no between-group differences in days free of mechanical ventilation. NMBAs are commonly used in critically ill patients who require IMV, but this practice has considerably changed throughout the years. In the 1980s, a survey from Great Britain reported NMBA administered in 90% of the patients on IMV [13] , while in 2005 data from an international large cohort of mechanically [16] . Meta-analyses emphasized potential limitations of those previous trials, i.e., heterogeneity in the use of sedation and prone positioning, selection bias, and crossover between treatment groups. Thus, to date, available evidence supports the use of NMBAs to reduce risks of barotrauma and to improve oxygenation [34] [35] [36] , but without clear advantage in survival rate. ARDS caused by COVID-19 has broad similarities with historic ARDS caused by other etiologies [18, [37] [38] [39] , although pulmonary blood flow derangement and resulting pulmonary shunt seem to play a primary role in COVID-19 ARDS [40] . As a result, various interventions previously employed for ARDS, such as prone position [41] , ECMO [42] , and NMBAs [20] , were extensively applied in the COVID-19 population.

In the current study, we describe COVID-19 patients who required mechanical ventilation and who presented with moderate-to-severe hypoxemia. Various small case reports have demonstrated substantial ventilatory asynchrony in these patients [43] [44] [45] and adverse sequelae, such as pneumothorax or pneumomediastinum [46] [47] [48] . Patients receiving NMBA also frequently required adjunctive therapies for ARDS, suggesting that the severity of disease was a main driver in the use of these drugs. The COVID-19 population subset investigated in our analyses provides further value to our study; indeed, we identified patients with moderate-to-severe ARDS who received, early, a short course of NMBAs. This is in line with protocols applied in previous large randomized trials [15, 16] . Furthermore, in comparison with the ACU-RASYS and ROSE trials, we found that the COVID-19 population presented similar baseline impairment in respiratory function, and as in the ROSE trial, low tidal volume and high PEEP were applied and resulted in comparable airway plateau pressure. However, prone position was used much more frequently in the COVID-19 population and the improvements associated with its use [49, 50] might have offset any additional benefit related to NMBAs.

Courcelle and collaborators specifically investigated the effects of NMBAs in COVID-19 patients enrolled in French/Belgian ICUs. After PS matching, they did not find a significant difference in 28-day mortality [20] . Similar to these preliminary findings, our international, larger multicenter investigation also found no improvement in 90-day mortality in patients receiving a short course of NMBAs. Moreover, Courcelle found a median duration of NMBA use of 5 days; conversely, our study appraised only patients who received continuous NMBAs up to 3 days, and interestingly, we found higher hazard ratio of 90-day mortality when the use of NMBAs was extended more than 3 days. This finding provides further evidence on the risks associated with NMBAs in COVID-19 patients. Relative to Courcelles PS-matched cohort, our matched population was similar in age, gender proportion, and BMI and initial ventilatory settings. Yet, in the French/Belgian study, prone positioning was applied in over 90% and ECMO in 15% of the patients, in line with French expert consensus guidelines recommendations [51] and evidence in favor of the aforementioned interventions in COVID-19 patients [42, 50, 52] . Conversely, our dataset encompassed a global population of critically ill COVID-19 patients, including ICU centers both with and without access to modalities such as prone position and ECMO, which were applied in up to 69.6% and 16.7%, respectively. In the early waves of the pandemic, centers with the highest patient numbers acknowledged that treatments, such as prone positioning, could sometimes not be accomplished, due to staffing limitations and patient safety concerns [11] . Similarly, NMBA use at centers with overwhelming ICU surge may also be greater than in less resource-constrained environments.

Non-depolarizing NMBAs inhibit the acetylcholine (Ach) receptor on the motor endplate and are available as aminosteroid or benzylisoquinolinium compound, the latter being the first choice for ARDS patients [53] . Renal and hepatic disease can drastically prolong the clearance of aminosteroid NMBAs, and in these conditions benzylisoquinolinium agents are preferred, since they undergo spontaneous degradation via Hofmann elimination. The NMBA class used was not recorded in our dataset, but severe liver and chronic kidney diseases were present in only 2.9 and 5.7% of the patients treated with NMBAs, respectively. Of note, COVID-19 can cause a broad variety of neurological symptoms and sequelae [54] and has been associated with the development of anti-Ach receptors antibodies [55] and myasthenia gravis [56] . Hence, comprehensive research is needed in this field to elucidate whether underlying neurological mechanisms associated with COVID-19 lead to an increased risk of death when NMBA are administered, and caution should be applied.

To date, this is the first international report of a large group of mechanically ventilated COVID-19 patients with severe disease to investigate the effects of NMBA. Employing comprehensive PS matching analysis, we demonstrated no improvement in mortality with NMBA use. These findings were consistently corroborated in sensitivity analysis applying inverse probability of treatment weighting. In addition, we provided inferences not limited by clinical practice specific to single-country studies. Nonetheless, certain limitations must be highlighted. First, although we conducted a comprehensive comparative PS-matched analysis that resulted in significant homogeneity in the evaluated sub-cohorts, the observational nature of our study and risks of drawing causal inferences should be highlighted. While immortal time bias could be considered as a potential limitation of this study design, given the short duration between cohort entry and exposure to NMBA therapy, the magnitude of the bias is expected to be limited [57] . In addition, similar times from ICU admission to death were found in both groups. Yet, the risk of immortal time bias should also be considered for several interventions potentially associated with survival-such as prone positioning, ECMO, and corticosteroids-since only patients who survived long enough to receive these interventions were analyzed. Second, the analyzed dataset did not provide any information on the type of NMBA, doses, or adequacy of the block achieved; hence, we were unable to extrapolate the appropriateness of their clinical use. In addition, the time of commencement of NMBA treatment was not recorded; thus, although we included patients who received NMBAs up to 3 days, and we explored different treatment regimens in sensitivity analyses, potential discrepancy in the clinically indicated 2-day full course of NMBA treatment [15, 16] should be considered when interpreting the results. Third, a proportion of the studied population was admitted to the ICU early during the pandemic. Consequently, potential logistical limitations related to managing critically ill patients treated in newly developed and understaffed ICUs, due to shortages in ICU beds or healthcare providers, must be considered. Fourth, as shown in Fig. 4 , residual discrepancy for Asia likely persisted post-matching. Further assessment on the potential reasons for this divergence suggested possible increased mortality in control patients from Asia. This could be related to the specific regional differences in ventilatory practice or application of adjunctive therapies, which should be further explored in future studies. Irrespective, potential random effect related to heterogenous practice among collaborating hospitals should be considered while interpreting our results. Fifth, the voluntary nature of site participation in the study must be emphasized, especially during a pandemic, as data may be skewed toward centers with enough resources to enter data. Lastly, as our observational study acquired data from routine clinical records, missing data could have biased our estimates.

Our results derived by PS-matching analysis suggest that among patients with COVID-19 and moderate-to-severe ARDS, early administration of a short course of NMBAs did not result in improved 90-day in-hospital mortality. In addition, NMBA treatment did not impact ventilatorfree days. Thus, the use of NMBAs should be cautiously assessed in this population, pending further studies that could elucidate specific indications for NMBAs in COVID-19. 

The online version contains supplementary material available at https:// doi. org/ 10. 1186/ s13054-022-03983-5.

Additional file 1. Appendix.

14 Department of Medicine, Columbia College of Physicians and Surgeons, and Center for Acute Respiratory Failure

27 Department of Intensive Care, The Alfred Hospital

31 Hospital Clinic of Barcelona

35 Avera McKennan Hospital & University Health Centre

44 Carilion Clinic, Roanoke, USA. 45 Chiba University Graduate School of Medicine

54 Medical ICU, Columbia College of Physicians and Surgeons

59 Fondazione Policlinico Universitario Agostino Gemelli IRCCS

68 Hamad General Hospital -Weill Cornell Medical College in Qatar

79 Hospital Emergencia Ate Vitarte, Ate, Peru. 80 Hospital for Tropical Diseases

São Paulo, Brazil. 90 INOVA Fairfax Hospital

Lankenau Institute of Medical Research (Main Line Health), Wynnewood, USA. 104 Launceston General Hospital

110 Manipal Hospital Whitefield, Bengaluru, India. 111 Mar del Plata Medical Foundation Private Community Hospital

121 Mount Sinai Medical Centre, Paranthal, India. 122 Nagoya University Hospital

Padova, Italy. 137 Pamela Youde Nethersole Eastern Hospital, Chai Wan, Hong Kong. 138 Papworth Hospitals NHS Foundation Trust

141 Perth Childrens Hospital

148 Dr Sulianti Saroso Hospital, Jakarta, Indonesia. 149 Providence Saint Johns Health Centre

167 Sao Camilo Cura Dars, Fortaleza, Brazil. 168 São João Hospital Centre

172 Severance Hospital, Seongnam-si, South Korea. 173 Shizuoka Childrens Hospital

178 St Christophers Hospital for Children

Tufts Medical Centre (and Floating Hospital for Children

218 COVID-19 Critical Care Consortium

227 Azienda Ospedaliero Universitaria Parma

232 Brigham and Womens Hospital

234 Brighton and Sussex Medical School

237 Cairo University Hospital, Cairo, Egypt. 238 Cardio VID, Medellín, Colombia. 239 Careggi Hospital

242 Cengkareng Hospital

246 Childrens Health Ireland (CHI) at Crumlin, Dublin, Ireland. 247 Childrens Hospital Bambino Gesù

259 Dignity Health Medical Group-Dominican

264 Suat Seren Chest Diseases and Surgery Practice and Training Centre

272 Giesinger Medical Centre, Philadelphia, USA. 273 Gunma University Graduate School of Medicine

277 Hamilton General Hospital

280 Hervey Bay Hospital, Pialba, Australia. 281 Hospital Clinico de La Pontificia Universidad Catolica

Porto Alegre, Brazil. 284 Hospital Felicio Rocho, Porto Alegre, Brazil. 285 Hospital Israelita Albert Einstein

287 Hospital Pro Cardíaco

290 Hospital Samaritano Paulista

296 Hospital Universitario Lucus Augusti, Lugo, Spain. 297 Hyogo Prefectural Kobe Childrens Hospital, Minatojima, Japan

309 King Abdullah Medical City, Manama, Bahrain. 310 King Abdullah Medical City Specialist Hospital, Mecca, Saudi Arabia. 311 King Abdullah Medical Complex

Saudi Arabia. 314 KK Womens and Childrens Hospital

324 Medanta Hospital, Gurgaon, India. 325 Medica Superspecialty Hospital, Kolkata, India. 326 Medical University of South Carolina

344 Novant Health (NH) Presbyterian Medical Centre

351 Parirenyatwa General Hospital, Harare, Zimbabwe. 352 Pediatric and Neonatal Cardiac Intensive Care at the American University

367 Ramón y Cajal University Hospital

372 Royal Brompton & Harefield NHS Foundation Trust, London, UK. 373 Royal Columbian Hospital

376 Saiseikai Yokohamashi Tobu Hospital

380 Sardjito Hospital, Yogyakarta, Indonesia. 381 Scripps Memorial Hospital La Jolla

392 The Catholic University of Seoul St Mary Hospital

397 Tuen Mun Hospital

409 University of Calgary and Alberta Health Services, Calgary, Canada. 410 University of California at San Diego

420 University of Wisconsin & American Family Childrens Hospital

COVID-19 Map-Johns Hopkins Coronavirus Resource Center

Baseline characteristics and outcomes of 1591 patients infected with SARS-CoV-2 admitted to ICUs of the Lombardy Region, Italy

Clinical characteristics of coronavirus disease 2019 in China

Presenting characteristics, comorbidities, and outcomes among 5700 patients hospitalized with COVID-19 in the New York City Area

Case characteristics, resource use, and outcomes of 10 021 patients with COVID-19 admitted to 920 German hospitals: an observational study

Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study

Quality of life and lung function in survivors of extracorporeal membrane oxygenation for acute respiratory distress syndrome

Clinical characteristics of 138 hospitalized patients with 2019 novel coronavirus-infected pneumonia in Wuhan

Case fatality rates for patients with COVID-19 requiring invasive mechanical ventilation

Effect of hydrocortisone on mortality and organ support in patients with severe COVID-19 the REMAP-CAP COVID-19 corticosteroid domain randomized clinical trial the members of the REMAP-CAP investigators appear in supplement 2

Interleukin-6 receptor antagonists in critically ill patients with Covid-19

Sedating drugs and neuromuscular blockade during mechanical ventilation

The techniques used to sedate ventilated patients. Intensive Care Med

Epidemiology, patterns of care, and mortality for patients with acute respiratory distress syndrome in intensive care units in 50 countries

Neuromuscular blockers in early acute respiratory distress syndrome

Early Neuromuscular Blockade in the Acute Respiratory Distress Syndrome

Epidemiology, clinical course, and outcomes of critically ill adults with COVID-19 in New York City: a prospective cohort study

Clinical features, ventilatory management, and outcome of ARDS caused by COVID-19 are similar to other causes of ARDS

Neuromuscular monitoring and neuromuscular blocking agent shortages when treating critically ill COVID-19 patients: a multicentre retrospective analysis

Neuromuscular blocking agents (NMBA) for COVID-19 acute respiratory distress syndrome: a multicenter observational study

The strengthening the reporting of observational studies in epidemiology (STROBE) statement: guidelines for reporting observational studies

Open source clinical science for emerging infections

Using research to prepare for outbreaks of severe acute respiratory infection

Using propensity scores to help design observational studies: application to the tobacco litigation

Constructing a control group using multivariate matched sampling methods that incorporate the propensity score

Balance diagnostics for comparing the distribution of baseline covariates between treatment groups in propensity-score matched samples

Accounting for missing data in statistical analyses: multiple imputation is not always the answer

Reappraisal of ventilator-free days in critical care research

Variance reduction in randomised trials by inverse probability weighting using the propensity score

Use of sedatives and neuromuscular blockers in a cohort of patients receiving mechanical ventilation

convenient online submission • thorough peer review by experienced researchers in your field • rapid publication on acceptance • support for research data, including large and complex data types • gold Open Access which fosters wider collaboration and increased citations maximum visibility for your research: over 100M website views per year your research ?

Clinical practice guidelines for sustained neuromuscular blockade in the adult critically ill patient

Effect of neuromuscular blocking agents on gas exchange in patients presenting with acute respiratory distress syndrome

Neuromuscular blockers in the acute respiratory distress syndrome: a meta-analysis

Validation of neuromuscular blocking agent use in acute respiratory distress syndrome: a meta-analysis of randomized trials

Neuromuscular blockade in acute respiratory distress syndrome: a systematic review and meta-analysis of randomized controlled trials

COVID-19-associated acute respiratory distress syndrome: is a different approach to management warranted?

COVID-19 acute respiratory distress syndrome (ARDS): clinical features and differences from typical pre-COVID-19 ARDS

Respiratory physiology of COVID-19-induced respiratory failure compared to ARDS of other etiologies

Covid-19 does not lead to a typical" acute respiratory distress syndrome

Prone position in intubated, mechanically ventilated patients with COVID-19: a multi-centric study of more than 1000 patients

Extracorporeal membrane oxygenation support in COVID-19: an international cohort study of the Extracorporeal Life Support Organization registry

Lung mechanics of mechanically ventilated patients with covid-19: analytics with highgranularity ventilator waveform data

Ventilation management and clinical outcomes in invasively ventilated patients with COVID-19 (PRoVENT-COVID): a national, multicentre, observational cohort study

Impact of asynchronies in acute respiratory distress syndrome due to coronavirus disease 2019

Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: a descriptive study

COVID-19 with spontaneous pneumomediastinum

Pneumothorax in COVID-19 disease-incidence and clinical characteristics

Prone positioning in severe acute respiratory distress syndrome

Feasibility and physiological effects of prone positioning in non-intubated patients with acute respiratory failure due to COVID-19 (PRON-COVID): a prospective cohort study

Formal guidelines: management of acute respiratory distress syndrome

Extracorporeal membrane oxygenation for severe acute respiratory distress syndrome associated with COVID-19: a retrospective cohort study

An observational study of the efficacy of cisatracurium compared with vecuronium in patients with or at risk for acute respiratory distress syndrome

Neurologic manifestations of hospitalized patients with coronavirus disease 2019 in Wuhan, China

SARS-CoV-2 infection associated with diplopia and anti-acetylcholine receptor antibodies

Myasthenia gravis associated with SARS-CoV-2 infection

Problem of immortal time bias in cohort studies: example using statins for preventing progression of diabetes

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations

We recognize the crucial importance of the ISARIC and SPRINT-SARI networks in developing and expanding the global COVID-19 Critical Care Consortium. We thank the generous support we received from the Extracorporeal Life Support Organization (ELSO) and the International ECMO Network (ECMONet). We greatly acknowledge Adrian Barnett for his invaluable contribution to the analysis plan and statistical insights. We greatly acknowledge Gabriele Fior (University of Milan) for his input and contribution to the critical revision of the manuscript. We owe Li Wenliang, MD, from the Wuhan Central Hospital, an eternal debt of gratitude for reminding the world that doctors should never be censored during a pandemic. Finally, we acknowledge all members of the

COVID-19 Critical Care Consortium and various collaborators as reported in the Additional file 1: Appendix.Author contributions GLB, JF, and JYS conceived of the study, participated in its design and coordination, and helped to draft the manuscript; KG and HD participated in the design of the study and coordination, conducted the statistical analysis, and helped to draft the manuscript; NW performed the statistical analysis and helped to draft the manuscript; AC, SS, SH, SF, JGL, EF, JoF, RB, DB, AB, DC, AA, ME, GG, CH, SH, CL, EM, LM, SM, AN, MO, PP, AT, and PY participated in collection of data and reviewed the initial draft of the manuscript; MP participated in collection of data, helped in the statistical analysis and reviewed the initial draft of the manuscript, Kristen Gibbons contributed to statistical analysis. All authors read and approved the final manuscript. 

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.