Around 5% of coronavirus disease-2019 (COVID-19) patients develop a critical form of the disease leading to ICU admission (1 , 2 ). The main manifestation of COVID-19 is acute hypoxemic respiratory failure (AHRF) requiring respiratory support. The optimal supportive treatment for AHRF related to COVID-19 (COVID-AHRF) is not yet established.
Initial guidelines advised against the use of noninvasive positive pressure ventilation (NIPPV) or high-flow nasal cannula (HFNC) as they are aerosol-generating procedures and, thus, a risk of infection among healthcare workers (3 ). However, some concerns promptly emerged regarding poor prognosis of patients with COVID-19 pneumonia related to invasive mechanical ventilation (IMV) use (4 ). In late March 2020, guidelines of the surviving sepsis campaign on the management of critically ill adults with COVID-AHRF advocated using HFNC first over other noninvasive techniques such as NIPPV and over IMV (5 ).
There is no clear-cut evidence of the greater effect of noninvasive oxygen support (NI-OS), including NIPPV, HFNC, continuous positive airway pressure (CPAP), and oxygen delivered through a nonrebreather face mask, than that of early IMV for COVID-AHRF management. NI-OS can improve patient outcome by avoiding intubation (6–8 ) and ventilator-associated complications (9 ) but can also worsen it by delaying intubation and increasing the risk of cardiac arrest before intubation (10 ). In addition, debate continues about the role of early intubation in preventing the risk of patient self-inflicted lung injury (P-SILI) (11–13 ). COVID-AHRF encompasses different patient phenotypes and thus, as reported by Gattinoni et al (12 ) (L = low elastance and H = high elastance), requires different ventilator supports. However, assessing patients to determine their phenotype at an early stage of ICU management is difficult, especially in the event of a massive influx of patients (13 ). Finally, the COVID-19 pandemic has underlined the need for a comprehensive national approach to ventilatory support management, given that the unprecedented influx of patients leads to a shortage of ICU capacity and a lack of resources such as ventilators (14 ). In such a context, HFNC could be administered in intermediate care units, thereby reserving ventilators for most needy patients. These complex issues have fueled incentives to conduct studies to provide evidence for improving decision-making processes.
In the absence of randomized clinical trials to assess the causal effect of IMV on mortality during the COVID pandemic, analysis of observational longitudinal studies is a suitable alternative (15 , 16 ).
The aim of this study was to assess the effect of early IMV in COVID-AHRF on 60-day mortality, using a large high-quality database and applying an inverse probability of treatment weight (IPTW) weighted Cox-survival model.
MATERIALS AND METHODS
Data Source
The study used data from the French prospective multicenter OutcomeRea database (n = 11 ICUs). The methods of data collection and quality of the database have been described in detail elsewhere (17 ). In accordance with French law, the OutcomeRea database was approved by the French Advisory Committee for Data Processing in Health Research and the French National Commission for Data Protection and Liberties (registration number 8999262). The database protocol was submitted to the Institutional Review Board of the Clermont-Ferrand University Hospital, France, which waived the need for informed consent (Institutional Review Board number 5891).
Study Population
Patients over 18 years were eligible for inclusion in the analysis if they were admitted to one of the participant ICUs belonging to the OutcomeRea network and had at admission an AHRF related to severe COVID-19 pneumonia defined as the combination of: 1) radiological features compatible with this diagnosis, 2) Pao 2 /Fio 2 ratio ≤ 300 mm Hg, and 3) a positive severe acute respiratory syndrome coronavirus-2 test using reverse-transcriptase polymerase chain reaction.
Patients were excluded if they were referred from another ICU or intermediate are unit, when a decision was made to discontinue life-sustaining treatments during the first 2 calendar days after ICU admission, if ICU length of stay was less than or equal to 2 days and if they had a shock or a Glasgow Coma Scale (GCS) less than or equal to 12 on ICU admission.
Data Collection
All data were prospectively collected and comprised those recorded on ICU admission, several variables collected throughout the ICU stay and outcomes (demographics, chronic disease/comorbidities as assessed by the Knaus Scale (18 ), baseline severity indexes, Simplified Acute Physiology Score [SAPS] II (19 ), and Sequential Organ Failure Assessment [SOFA] (20 ) scores; daily throughout ICU stay: clinical and biological parameters, requirement for NI-OS and IMV, other organ support [vasopressors and renal replacement therapy], and occurrence of ICU-acquired pneumonia and bacteremia; ICU and hospital length of stay; and vital status at ICU and hospital discharge, and at day 60 after ICU admission).
Definitions, Group Assignment, and Respiratory Support Modalities
COVID-AHRF was classified into three categories based on the severity of hypoxemia at admission: mild (200 < Pao 2 /Fio 2 ≤ 300 mm Hg), moderate (100 < Pao 2 /Fio 2 < 200 mm Hg), and severe (Pao 2 /Fio 2 ≤ 100 mm Hg) (21 ).
The early IMV group comprised all patients who received IMV within the first 2 calendar days in after ICU admission and the nonearly IMV group comprised all other patients who received at least one of these NI-OSs: NIPPV, HFNC, CPAP, and oxygen delivered through a nonrebreather face mask and not earlier than the third day after ICU admission.
Strategies for IMV and descriptions of the NI-OS techniques are given in the Online Data Supplement (https://links.lww.com/CCX/A495 ).
Statistical Analysis
Patient characteristics were expressed as n (%) for categorical variables and median (interquartile range [IQR]) for continuous variables. Comparisons were made with exact Fisher tests for categorical variables and Wilcoxon tests for continuous variables.
To estimate the average causal effect of early IMV on day-60 mortality, on longitudinal observational data, we used an IPTW estimator, which is the inverse of the patients’ predicted probability of being in the early IMV group on the basis of their baseline covariates.
IPTW estimator is a statistical tool that allows causal inference on observational data. It was developed nearly 20 years ago (22–26 ) and has been widely used in medical fields. The methodology of IPTW has been fully described elsewhere. It is detailed in the Online Data Supplement (https://links.lww.com/CCX/A495 ) and briefly presented here. When randomized clinical trials are impossible or difficult to be implemented in the presence of baseline confounding factors, IPTW estimators are a suitable alternative for observational data to estimate the real causal effect of a treatment on outcome. By weighting all patients with their IPTW, two pseudopopulations are created, similar with regard to baseline confounding factors, and their outcomes are compared with survival models.
The impact of early IMV on day-60 mortality was estimated by a two-step process: 1) weight estimation by the IPTW estimator and 2) estimation of the impact of early IMV on day-60 mortality by a weighted Cox model. In the first step, the weight model, a nonparsimonious multivariate logistic regression model, was constructed to estimate each patient’s predicted probability of being in the early IMV group. We included in the weight model the following baseline covariates, recorded at admission and not affected by study groups: period of admission, time between symptom onset and ICU admission, time between hospital and ICU admission, age, gender, body mass index, comorbidities including presence of chronic liver failure, cardiovascular, respiratory and kidney chronic diseases, immunosuppression, clinical and laboratory features at admission, T greater than 39°C, renal SOFA item (> 2), GCS < 15, Pao 2 /Fio 2 ratio, respiratory rate, lactatemia, lymphocyte count, ferritin, d -dimer plasma level, C-reactive protein serum level, and treatments received at admission including lopinavir ritonavir, hydroxychloroquine, tocilizumab, anakinra, and corticosteroids received at admission. In the model, continuous variables were kept linearly unless in the absence of log linearity. All variables included in the weight model reflected knowledge available at baseline. To avoid extreme weights, we used stabilized weights, and to ensure the positivity, assumption was respected, weights were truncated at the 1–99th percentile (27 ) (online data supplement, https://links.lww.com/CCX/A495 ). In the second step, we used a weighted Cox proportional-hazard model to estimate the risk of death within the first 60 days of the ICU stay of early IMV. A hazard ratio greater than 1 indicated an increased risk of death. The proportionality of hazard risk for IMV was tested with martingale residuals. Further analyses were performed to confirm the results obtained with the IPTW: 1) a raw (nonweighted) multivariable Cox model with adjustment on severity of the patients at admission and 2) a case-control analysis for which cases (patients who had died at day 60) were matched with controls (survivors) based on age, SOFA without respiratory item, and Pao 2 /Fio 2 ratio at admission. Characteristics of the population by centers are reported in Table E1 (https://links.lww.com/CCX/A495 ). All models were stratified by centers. To assess the potential effect of early IMV on subgroups of patients, sensitivity analyses were performed in the subgroup of the patients: 1) admitted to hospital within less than 5 days, 2) with a Pao 2 /Fio 2 less than 150 mm Hg at admission, and 3) with a Pao 2 /Fio 2 greater than 150 mm Hg at admission. These groups were defined prior to the statistical analyses. For all tests, a two-sided α of 0.05 was considered to be significant. Missing baseline variables were handled by median. All statistical analyses were performed with the SAS software, Version 9.4 (SAS Institute, Cary, NC).
RESULTS
Patients and Baseline Characteristics
From February 15, 2020, to May 1, 2020, 400 patients with laboratory confirmed COVID-19 were admitted to ICUs of the OutcomeRea network. Of the 245 included in the study (Fig. 1 ) 187 (76.4%) were male. The overall median (IQR) was 61 years (52–69 yr). The sex distribution and median age of included and excluded patients were similar (Table E2 , https://links.lww.com/CCX/A495 ). One or more comorbidities were present in 157 patients (64%), with obesity and cardiovascular disease being the most frequently coexisting medical conditions, confirmed in 89 patients (36.4%) and 61 patients (24.8%), respectively. The median duration from the onset of symptoms to ICU admission was 10 days (7–12 d), and between hospital and ICU admission, 2 days (1–4 d). At admission, the median (IQR) SAPS II score was 34 (26–44). Median (IQR) Pao 2 /Fio 2 ratio was 121 mm Hg (90–174) and 206 patients (84%) had moderate-to-severe AHRF. Overall, 117 patients (47.8%) received IMV (early IMV group), oxygen by mask or nasal prongs 16 (6.6%), HFNC 85 (34.6%), CPAP 18 (7.4%), and NIPPV 9 (3.6%). After intubation, their median (IQR) tidal volume was 6 mL/kg (5.9–6.4 mL/kg) with a median (IQR) compliance of 36.7 mL/mm Hg (27.6–53.3 mL/mm Hg). Ninety-five patients (38.8%) received lopinavir-ritonavir, 21 patients (8.6%) hydroxychloroquine, 22 patients (9%) tocilizumab, 22 patients (9%) anakinra, and 68 patients (28%) corticosteroids. The comparison of baseline characteristics between the early and the nonearly IMV groups is shown in Table E3 (https://links.lww.com/CCX/A495 ). The rates of ICU-acquired pneumonia, bacteremia, and the ICU length of stay were higher in the early than in the nonearly IMV group (p < 0.01, p < 0.01, and p < 0.01, respectively).
Figure 1.: Flowchart. AHRF = acute hypoxemic respiratory failure, COVID-19 = coronavirus disease 2019, GCS = Glasgow Coma Scale, IMV = invasive mechanical ventilation, LOS = length of stay.
Propensity Score Model Development
Propensity scores ranged from 0.01 to 0.92 and from 0.02 to 0.97 in the no-early IMV and in the early IMV groups, respectively, with 93.8% in the region of common support (0.02–0.92) (Fig. E1 , https://links.lww.com/CCX/A495 ). All the covariates in the planned propensity score were kept in the final model. After applying IPTW, all covariates in the planned propensity score had weighted standardized differences below 10%, which is in favor of an equilibration of the covariates between the subgroups and ensure the exchangeability at baseline for these confounders (Table E4 (https://links.lww.com/CCX/A495 ); and Fig. 2 ).
Figure 2.: Standardized differences of variables used to generate the propensity score before and after proportional-hazard model weighted by inverse probability (IPTW). Propensity scores ranged from 0.01 to 0.92 and from 0.02 to 0.97 in the no early invasive mechanical ventilation (IMV) and in the early IMV groups, respectively, with 93.8% in the region of common support (0.02–0.92). All the covariates in the planned propensity score were kept in the final model. After applying IPTW, all covariates in the planned propensity score had weighted standardized differences below 10%, which is in favor of an equilibration of the covariates between the subgroups and ensures the exchangeability at baseline for these confounders. ACE = angiotensin-converting-enzyme inhibitors, SOFA = Sequential Organ Failure Assessment.
Follow-Up and Outcomes
The day-60-mortality in the whole study population was 31.8%. It was higher in the early than in the nonearly IMV group: 42.7% versus 21.9% (pval < 0.01), respectively.
After weighted Cox model analysis, the risk of death at day 60 was higher in patients receiving early IMV (HRw = 1.74, 95% CI, 1.07–2.83; p = 0.03). Similar results were observed in a sensitivity analysis using truncated HRw (Table E5 , https://links.lww.com/CCX/A495 ). In all subgroup analyses, early IMV was or tended to be associated with an increased risk of day-60 mortality (Fig. 3 ).
Figure 3.: Effect of early invasive mechanical ventilation on ICU death of patients in the main cohort and in different subgroups (sensitivity analyses).
In other sensitivity analyses, early IMV was also associated with day-60-mortality (Tables E6 and E7 , and Fig. E2 , https://links.lww.com/CCX/A495 ).
In addition, we found that in the case of delayed intubation, that is, intubation after the first 2 calendar days after ICU admission, the patients in the late-IMV group (n = 45, 18.4%) had a similar outcome to those in the early IMV group, with a day-60-mortality of 42.2% and 42.7%, respectively. Patients without any IMV (n = 83, 33.9%) had the best survival rate during the ICU stay, with a day-60 mortality of 10.8% (Table E8 , https://links.lww.com/CCX/A495 ).
DISCUSSION
Our study adds to current knowledge on the management of COVID-19 patients. It shows that using IMV during the first 2 calendar days after ICU admission in critically ill COVID-19 patients with AHRF was associated with increased day-60 mortality as compared to initial use of NI-OS. It also shows that patients intubated later because of failure of the noninvasive oxygenation strategy had a similar outcome to those ventilated early. The day-60 survival of patients with this successful conservative noninvasive strategy was better.
Deciding when to initiate IMV in critically ill patients with AHRF is challenging. The benefits of IMV must be weighed against its inherent risks (9 ). Since IMV is associated with adverse events entailing substantial morbidity and mortality, physicians have developed multiple means of NI-OS that avoid endotracheal intubation, such as NIPPV, HFNC, and CPAP. In de novo AHRF, the safety of NI-OS is debated. European Guidelines do not recommend CPAP and NIPPV, because they can postpone endotracheal intubation and increase the risk of hypoxic cardiac arrest (28 ). However, HFNC has recently shown clinical benefits in de novo AHRF (29 ). In COVID-19 patients with AHRF or ARDS, the question is even more complex, especially in the early phase of the illness, because these patients may have normal lung mechanics and elastance (11 , 12 , 30 ). Thus, several experts have argued that protective IMV with effective sedation and paralysis should be implemented early to prevent delays in endotracheal intubation with the subsequent risk of increased mortality, barotrauma, volotrauma, and P-SILI due to large respiratory effort during noninvasive ventilatory support (12 ). The risk of aerosolization of viral particles and contamination of healthcare workers has also been cited to warrant early intubation instead of using NI-OS (31 , 32 ). To date, no studies related to COVID-19 have directly evidenced this risk. Furthermore, findings with regard to droplet dispersion and aerosol generation with HFNC are uncertain and the increased risks of aerosolization with NIPPV compared with HFNC are largely unknown (33–35 ).
The need for IMV varied widely, from 29.1% to 89.9%, among patients admitted to the ICU with COVID-AHRF. IMV is invariably associated with high mortality ranging from 16% to 78% (2 , 4 , 31 , 36 , 37 ). Very few studies have assessed the impact on the mortality of the mode of respiratory support in severe forms of COVID-19 pneumonia. In a cohort from Lombardy (36 ) of 3,940 ICU COVID-19 patients, 2,929 (87.3%) underwent intubation. The patients initially treated with NI-OS had a lower risk for mortality (heart rate [HR], 0.62; 95% CI, 0.52–0.75; p < 0.01) than those treated with IMV. Two recent retrospective cohort studies were not able to show a deleterious effect of intubation and subsequent IMV. Demoule et al (38 ) showed that HFNC significantly reduced intubation and subsequent IMV—55% (95% CI, 46–63) in the HFNC patient group (n = 146) versus 72% (95% CI, 64–79%) in patients without HFNC (n = 233)—but did not affect case fatality. Hernandez-Romieu et al (39 ) showed that mortality among 231 patients was neither associated with time from ICU admission to intubation nor with HFNC use.
Hospital capacity in France was exceeded during the pandemic. Some patients were intubated early to be transferred safely to the ICU and consequently did not undergo HFNC or other techniques to delay intubation. Furthermore, the recommendations, published at the beginning of the pandemic (3 ), advised against NIPPV and HFNC to avoid virus transmission. Some of our patients, therefore, were intubated before or at admission without having the option of noninvasive procedures.
In our study, higher mortality was observed in the early IMV group. Several factors could have contributed to this result. First, most patients in the nonearly IMV group received HFNC for oxygenation and respiratory support (62.4%) and very few received NIPPV. In de novo AHRF, the impact of HFNC on survival is uncertain (32 , 40 ), but several studies have reported a lower risk of mortality associated with its use (41 , 42 ) and recent meta-analyses have concluded that HFNC reduces the risk for endotracheal intubation (29 , 40 ). The benefit of HFNC over standard oxygen delivered through a facemask or by NIPPV can be explained by its physiologic effects. HFNC provides a high and stable Fio 2 , delivers a positive end-expiratory pressure value maintenance around 4 mm Hg, and can adjust independently the amount of Fio 2 and the flow rate in an attempt to match the patient’s oxygenation requirements and inspiratory flow demands. Thus, HFNC seems to improve oxygenation while reducing excessive respiratory efforts, pleural pressure swings, and exacerbation of lung injury (P-SILI) (41 ).
Second, NI-OS significantly reduced nosocomial infections. Several studies have already shown that NIPPV in critically ill patients with acute exacerbation of chronic obstructive pulmonary disease or severe cardiogenic pulmonary edema was associated with improved survival and reduction of nosocomial infections (43 , 44 ).
One important finding of our study was that patients who ultimately did need IMV according to physicians had a high mortality rate. This underlines again the complexity of the timing of intubation in patients suffering from COVID-AHRF. This finding is consistent with the results from the Lombardy cohort (36 ), in which the 151 patients who were initially treated noninvasively and subsequently underwent intubation had a lower chance of survival than the 199 patients who underwent NI-OS throughout the ICU stay (HR, 1.69; 95% CI, 1.43–1.98; p < 0.01). The mortality of the patients undergoing subsequent intubation was similar to that of patients who were treated with mechanical ventilation at ICU admission: HR for IMV versus NIV failure, 1.20; 95% CI, 0.95–1.53; p = 0.12). Tools have been developed to avoid a deleterious delayed intubation in patients suffering from COVID-AHRF (45 , 46 ). In pneumonia patients with AHRF under HFNC, the ratio of oxygen saturation as measured by pulse oximetry/Fio 2 to respiratory rate is an index that can help to identify those with low and those with high risk for intubation (47 ).
Our study has certain limitations. First, although the question of whether noninvasive oxygenation techniques can be administered outside ICUs is crucial in the case of limited ICU capacity and lack of ventilators, we focused on the impact of early IMV in ICU patients. Consequently, the study does not address the potential role of NI-OS in intermediate care units or wards. However, if hospitals and physicians cannot do otherwise, NI-OS could be an alternative in such a pandemic context, provided that strict protocols are discussed beforehand among the various hospital stakeholders involved. Second, most of our patients were managed with HFNC, and none were under Helmet, which limits the external validity of our results. Third, we pooled all the noninvasive ventilation procedures and could not compare HNFC with either NIPPV or CPAP. Fourth, several data could not be collected or determined retrospectively, including the phenotype (L or H) of the COVID-AHRF, the tidal volume in patients under noninvasive oxygenation techniques, and the medical reasons leading physicians to intubate or not the patients and the timing of the procedure. Finally, our study was only observational. As a result, although we tried to consider most confounding factors, our model could still be biased mostly because of unmeasured confounding. Then, to avoid selection bias and bias due to immortal time, we excluded all the patients transferred from another ICU or intermediate care unit who might have benefitted from NI-OS before ICU admission and also all patients with shock or neurologic failure at admission, for whom intubation is recommended. In addition, all our sensitivity analyses yielded the same results.
CONCLUSIONS
In ICU patients admitted with COVID-19-induced acute hypoxemic respiratory failure, early IMV was associated with an increased day-60 mortality. Those results deserved to be confirmed by randomized controlled trials.
ACKNOWLEDGMENT
The OutcomeRea study group members thank Celine Feger, MD, from Emibiotech for assistance in medical writing.
REFERENCES
1. Wu Z, McGoogan JM. Characteristics of and important lessons from the coronavirus disease 2019 (COVID-19) outbreak in China: Summary of a report of 72 314 cases from the Chinese Center for Disease Control and Prevention. JAMA. 2020; 323:1239–1242
2. Xie J, Wu W, Li S, et al. Clinical characteristics and outcomes of critically ill patients with novel coronavirus infectious disease (COVID-19) in China: A retrospective multicenter study. Intensive Care Med. 2020; 46:1863–1872
3. Bouadma L, Lescure FX, Lucet JC, et al. Severe SARS-CoV-2 infections: Practical considerations and management strategy for intensivists. Intensive Care Med. 2020; 46:579–582
4. Wunsch H. Mechanical ventilation in COVID-19: Interpreting the current epidemiology. Am J Respir Crit Care Med. 2020; 202:1–4
5. Alhazzani W, Møller MH, Arabi YM, et al. Surviving sepsis campaign: Guidelines on the management of critically ill adults with coronavirus disease 2019 (COVID-19). Crit Care Med. 2020; 48:e440–e469
6. Lemiale V, Mokart D, Resche-Rigon M, et al.; Groupe de Recherche en Réanimation Respiratoire du patient d’Onco-Hématologie (GRRR-OH). Effect of noninvasive ventilation vs oxygen therapy on mortality among immunocompromised patients with acute respiratory failure: A randomized clinical trial. JAMA. 2015; 314:1711–1719
7. Lemiale V, Mokart D, Azoulay E; Groupe de Recherche en Réanimation Respiratoire Onco-Hématologique (GRRR-OH). Noninvasive ventilation and outcomes among immunocompromised patients–reply. JAMA. 2016; 315:1902–1903
8. Demoule A, Chevret S, Carlucci A, et al.; oVNI Study Group; REVA Network (Research Network in Mechanical Ventilation). Changing use of noninvasive ventilation in critically ill patients: Trends over 15 years in francophone countries. Intensive Care Med. 2016; 42:82–92
9. Schnell D, Timsit JF, Darmon M, et al. Noninvasive mechanical ventilation in acute respiratory failure: Trends in use and outcomes. Intensive Care Med. 2014; 40:582–591
10. Thille AW, Contou D, Fragnoli C, et al. Non-invasive ventilation for acute hypoxemic respiratory failure: Intubation rate and risk factors. Crit Care. 2013; 17:R269
11. Marini JJ, Gattinoni L. Management of COVID-19 respiratory distress. JAMA. 2020; 323:2329–2330
12. Gattinoni L, Chiumello D, Caironi P, et al. COVID-19 pneumonia: Different respiratory treatments for different phenotypes? Intensive Care Med. 2020; 46:1099–1102
13. Tobin MJ, Laghi F, Jubran A. Caution about early intubation and mechanical ventilation in COVID-19. Ann Intensive Care. 2020; 10:78
14. Tobin MJ. Basing respiratory management of COVID-19 on physiological principles. Am J Respir Crit Care Med. 2020; 201:1319–1320
15. Ross JS. Randomized clinical trials and observational studies are more often alike than unlike. JAMA Intern Med. 2014; 174:1557
16. Anglemyer A, Horvath HT, Bero L; The Cochrane Collaboration. Healthcare outcomes assessed with observational study designs compared with those assessed in randomized trials. In: Cochrane Database of Systematic Reviews. 2014, Chichester, UK: John Wiley & Sons.
17. Zahar JR, Timsit JF, Garrouste-Orgeas M, et al. Outcomes in severe sepsis and patients with septic shock: Pathogen species and infection sites are not associated with mortality. Crit Care Med. 2011; 39:1886–1895
18. Knaus WA, Wagner DP, Draper EA, et al. The APACHE III prognostic system. Risk prediction of hospital mortality for critically ill hospitalized adults. Chest. 1991; 100:1619–1636
19. Le Gall JR, Lemeshow S, Saulnier F. A new simplified acute physiology score (SAPS II) based on a European/North American multicenter study. JAMA. 1993; 270:2957–2963
20. Vincent JL, de Mendonça A, Cantraine F, et al. Use of the SOFA score to assess the incidence of organ dysfunction/failure in intensive care units: Results of a multicenter, prospective study. Working group on “sepsis-related problems” of the European Society of Intensive Care Medicine. Crit Care Med. 1998; 26:1793–1800
21. Acute Respiratory Distress Syndrome. The Berlin definition. JAMA. 2012; 307:2526–2533
22. Hernán MA, Robins JM. Causal Inference: What If. Boca Raton, Chapman & Hall/CRC, 2020
23. Hernán MA, Robins JM. Estimating causal effects from epidemiological data. J Epidemiol Community Health. 2006; 60:578–586
24. Austin PC, Mamdani MM. A comparison of propensity score methods: A case-study estimating the effectiveness of post-AMI statin use. Stat Med. 2006; 25:2084–2106
25. Austin PC, Schuster T. The performance of different propensity score methods for estimating absolute effects of treatments on survival outcomes: A simulation study. Stat Methods Med Res. 2014; 25:2214–2237
26. Heinze G, Jüni P. An overview of the objectives of and the approaches to propensity score analyses. Eur Heart J. 2011; 32:1704–1708
27. Cole SR, Hernán MA. Constructing inverse probability weights for marginal structural models. Am J Epidemiol. 2008; 168:656–664
28. Rochwerg B, Brochard L, Elliott MW, et al. Official ERS/ATS clinical practice guidelines: Noninvasive ventilation for acute respiratory failure. Eur Respir J. 2017; 50:1602426
29. Rochwerg B, Granton D, Wang DX, et al. High flow nasal cannula compared with conventional oxygen therapy for acute hypoxemic respiratory failure: A systematic review and meta-analysis. Intensive Care Med. 2019; 45:563–572
30. Pfeifer M, Ewig S, Voshaar T, et al. Position paper for the state-of-the-art application of respiratory support in patients with COVID-19. Respiration. 2020; 99:521–542
31. Schünemann HJ, Khabsa J, Solo K, et al. Ventilation techniques and risk for transmission of coronavirus disease, including COVID-19: A living systematic review of multiple streams of evidence. Ann Intern Med. 2020; 173:204–216
32. Agarwal A, Basmaji J, Muttalib F, et al. High-flow nasal cannula for acute hypoxemic respiratory failure in patients with COVID-19: Systematic reviews of effectiveness and its risks of aerosolization, dispersion, and infection transmission. Can J Anesth Can Anesth. 2020; 67:1217–1248
33. Remy KE, Lin JC, Verhoef PA. High-flow nasal cannula may be no safer than non-invasive positive pressure ventilation for COVID-19 patients. Crit Care. 2020; 24:169
34. Li J, Fink JB, Ehrmann S. High-flow nasal cannula for COVID-19 patients: Low risk of bio-aerosol dispersion. Eur Respir J. 2020; 55:2000892
35. Honore PM, Barreto Gutierrez L, Kugener L, et al. Compared to NIPPV, HFNC is more dangerous regarding aerosol dispersion and contamination of healthcare personnel: We are not sure. Crit Care. 2020; 24:482
36. Grasselli G, Greco M, Zanella A, et al. Risk factors associated with mortality among patients with COVID-19 in intensive care units in Lombardy, Italy. JAMA Intern Med. 2020; 180:1345–1355
37. Chang R, Elhusseiny KM, Yeh YC, et al. COVID-19 ICU and mechanical ventilation patient characteristics and outcomes - a systematic review and meta-analysis. medRxiv. 202008.16.20035691
38. Demoule A, Vieillard Baron A, Darmon M, et al. High flow nasal canula in critically ill severe COVID-19 patients. Am J Respir Crit Care Med. 2020; 202:1039–1042
39. Hernandez-Romieu AC, Adelman MW, Hockstein MA, et al. Timing of intubation and mortality among critically ill coronavirus disease 2019 patients: A single-center cohort study. Crit Care Med. 2020; 48:e1045–e1053
40. Ferreyro BL, Angriman F, Munshi L, et al. Association of noninvasive oxygenation strategies with all-cause mortality in adults with acute hypoxemic respiratory failure: A systematic review and meta-analysis. JAMA. 2020; 324:57–67
41. Frat JP, Thille AW, Mercat A, et al.; FLORALI Study Group; REVA Network. High-flow oxygen through nasal cannula in acute hypoxemic respiratory failure. N Engl J Med. 2015; 372:2185–2196
42. Macé J, Marjanovic N, Faranpour F, et al. Early high-flow nasal cannula oxygen therapy in adults with acute hypoxemic respiratory failure in the ED: A before-after study. Am J Emerg Med. 2019; 37:2091–2096
43. Girou E, Brun-Buisson C, Taillé S, et al. Secular trends in nosocomial infections and mortality associated with noninvasive ventilation in patients with exacerbation of COPD and pulmonary edema. JAMA. 2003; 290:2985–2991
44. Girou E, Schortgen F, Delclaux C, et al. Association of noninvasive ventilation with nosocomial infections and survival in critically ill patients. JAMA. 2000; 284:2361–2367
45. Apigo M, Schechtman J, Dhliwayo N, et al. Development of a work of breathing scale and monitoring need of intubation in COVID-19 pneumonia. Crit Care. 2020; 24:477
46. Zucman N, Mullaert J, Roux D, et al.; Contributors. Prediction of outcome of nasal high flow use during COVID-19-related acute hypoxemic respiratory failure. Intensive Care Med. 2020; 46:1924–1926
47. Roca O, Caralt B, Messika J, et al. An index combining respiratory rate and oxygenation to predict outcome of nasal high-flow therapy. Am J Respir Crit Care Med. 2019; 199:1368–1376
APPENDIX. OUTCOMEREA STUDY GROUP
Jean-François Timsit, Elie Azoulay, Maïté Garrouste-Orgeas, Jean-Ralph Zahar, Bruno Mourvillier, Christophe Clec’h, Corinne Alberti, Stephane Ruckly, Sébastien Bailly, Aurélien Vannieuwenhuyze, Christophe Adrie, Carole Agasse, Bernard Allaouchiche, Olivier Andremont, Pascal Andreu, Laurent Argaud, Claire Ara-Somohano, Elie Azoulay, Francois Barbier, Déborah Boyer, Jean-Pierre Bedos, Thomas Baudry, Jérome Bedel, Julien Bohé, Lila Bouadma, Jeremy Bourenne, Noel Brule, Cédric Brétonnière, Frank Chemouni, Christine Cheval, Julien Carvelli, Elisabeth Coupez, Martin Cour, Claire Dupuis, Etienne de Montmollin, Loa Dopeux, Anne-Sylvie Dumenil, Jean-Marc Forel, Marc Gainnier, Charlotte Garret, Dany Goldgran-Tonedano, Steven Grangé, Antoine Gros, Hédia Hammed, Akim Haouache, Romain Hernu, Tarik Hissem, Vivien Hong Tuan Ha, Sébastien Jochmans, Jean-Baptiste Joffredo, Hatem Kallel, Guillaume Lacave, Virgine Laurent, Alexandre Lautrette, Clément Le bihan Eric Magalhaes, Virgine Lemiale, Guillaume Marcotte, Jordane Lebut, Maxime Lugosi, Sibylle Merceron, Benoît Misset, Mathild Neuville, Laurent Nicolet, Johanna Oziel, Laurent Papazian, Juliette Patrier, Benjamin Planquette, Aguila Radjou, Marie Simon, Romain Sonneville, Jean Reignier, Bertrand Souweine, Carole Schwebel, Shidasp Siami, Romain Sonneville, Nicolas Terzi, Gilles Troché, Marie Thuong, Guillaume Thierry, Marion Venot, Sondes Yaacoubi, Olivier Zambon, Julien Fournier, Stéphanie Bagur, Mireille Adda, Vanessa Vindrieux, Sylvie de la Salle, Pauline Enguerrand, Vincent Gobert , Stéphane Guessens, Helene Merle, Nadira Kaddour, Boris Berthe, Samir Bekkhouche, Kaouttar Mellouk, Mélaine Lebrazic, Carole Ouisse, Diane Maugars, Christelle Aparicio, Igor Theodose, Manal Nouacer, Veronique Deiler, Fariza Lamara, Myriam Moussa, Atika Mouaci, and Nassima Viguier.