key: cord-0701771-nlhpou4k
authors: Parish, Austin J.; West, Jason R.; Caputo, Nicholas D.; Janus, Trevor M.; Yuan, Denley; Zhang, John; Singer, Daniel J.
title: Early Intubation and Increased Coronavirus Disease 2019 Mortality: A Propensity Score–Matched Retrospective Cohort Study
date: 2021-06-15
journal: Crit Care Explor
DOI: 10.1097/cce.0000000000000452
sha: 53200f5c57182b2b6e59fa2165ce9a26fc52381b
doc_id: 701771
cord_uid: nlhpou4k

OBJECTIVES: There has been controversy about the timing and indications for intubation and mechanical ventilation in novel coronavirus disease 2019. This study assessed the effect of early intubation and mechanical ventilation on all-cause, inhospital mortality for coronavirus disease 2019 patients. DESIGN: Multicenter retrospective cohort study. SETTING: Eleven municipal hospitals in New York City from March 1, 2020, to December 1, 2020. PATIENTS: Adult patients who tested positive for coronavirus disease 2019 in the emergency department were subsequently admitted. Patients with do-not-intubate orders at admission were excluded. INTERVENTIONS: Intubation within 48 hours of triage and intubation at any point during hospital stay. MEASUREMENTS AND MAIN RESULTS: Data from 7,597 coronavirus disease 2019 patients were included; of these, 1,628 (21%) were intubated overall and 807 (11%) were intubated within 48 hours of triage. After controlling for available confounders, intubation rates for coronavirus disease 2019 patients varied significantly across hospitals and decreased steadily as the pandemic progressed. After nearest neighbor propensity score matching, intubation within 48 hours of triage was associated with higher all-cause mortality (hazard ratio, 1.30 [1.15–1.48]; p < 0.0001), as was intubation at any time point (hazard ratio, 1.62 [1.45–1.80]; p < 0.0001). Among intubated patients, intubation within 48 hours of triage was not significantly associated with differences in mortality (hazard ratio, 1.09 [0.94–1.26]; p = 0.26). These results remained robust to multiple sensitivity analyses. CONCLUSIONS: Intubation within 48 hours of triage, as well as at any time point in the hospital course, was associated with increased mortality in coronavirus disease 2019 patients in this observational study.

hours (9) . The optimal treatment of respiratory failure from COVID-19 has yet to be established, and it has been suggested that COVID-19 patients should have an earlier threshold for receiving invasive mechanical ventilation due to high inspiratory efforts and transpulmonary pressures causing self-inflicted lung injury (8, 10) . However, there are no robust data on the timing of tracheal intubation and its relationship to mortality for COVID-19 patients, and the existing studies have inconsistent findings (11) (12) (13) , suggesting either no relationship between time to invasive mechanical ventilation and mortality (11) or that early intubation may affect mortality (12, 13) .

In this observational study using data from 11 public hospitals in New York City gathered from March to December 2020, we examined the relationship between invasive mechanical ventilation and inhospital mortality, using propensity score matching and adjustment to reduce confounding. We looked at the effect on mortality of intubation within 48 hours of hospital triage, as well as the effect of intubation at any time point.

We conducted a retrospective chart review of 11 New York City Health and Hospitals (H + H), public hospitals for all adult patients, seen in emergency departments (EDs) between March 1, 2020, and December 1, 2020, who were tested with a polymerase chain reaction test for SARS-CoV-2 (COVID-19) during their time in the ED and subsequently admitted. Patients with negative, discontinued, or indeterminate tests were excluded, as were patients that were transferred to hospitals outside of the New York City H + H system. As the same patient may have presented to the ED multiple times, we used only the earliest visit resulting in admission so that each patient contributed unique, noncorrelated data. We obtained institutional review board (IRB) approval for this study both from the Lincoln Medical Center IRB and from the New York City H + H IRB (IRB approval number: 20-013).

We extracted a range of demographic and clinical data for each patient, including initial labs obtained within 24 hours of triage. For each patient, we extracted their inhospital exposure to invasive ventilation (endotracheal intubation and mechanical ventilation).

Data were extracted automatically from the Epic electronic medical record system (14) . For a subset of 300 patients, all variables were extracted manually by three authors (A.J.P., D.Y., T.M.J.), and Cohen kappa was calculated to the assess accuracy of the automatic extraction. For continuous variables where 30% or fewer patients had missing values, we used multiple imputation with predictive mean matching to estimate the missing numerical data (15) ; variables with more than 30% missing values were not used in any models.

We excluded patients that died within 48 hours of triage. We also excluded all patients that had donot-intubate (DNI) orders placed within the first 48 hours of triage. To reduce the effect of immortal time bias, in the primary analysis, we compared patients that were intubated within 48 hours of ED triage with all other patients, including those intubated after this time point (16, 17) .

The primary outcome was all-cause, inhospital mortality; patients that were discharged, left against medical advice, or were still in the hospital as of December 1, 2020, were censored. We calculated the hazard ratio (HR) for mortality using the Cox proportional hazards model (18) and estimated median survival times using Kaplan-Meier analysis (19) . We calculated propensity scores for intubation with multivariate logistic regression and used 1:1 nearest neighbor propensity score matching without replacement (20) to balance covariates between those that were intubated and those that were not (21) . Additionally, we used propensity score adjustment (22) and overlap propensity score weighting (23) as alternative methods to estimate the primary outcome.

We included in the propensity score model variables considered likely to be related to the intervention (intubation) and primary outcome (mortality) (24) . The variables that were included were as follows: age, sex, body mass index, ever-smoker, ethnicity, race, hospital that the patient was admitted to, admission date, total hospital census (number of inpatients) at admission date, whether the patient was a Medicare or Medicaid recipient, history of diabetes mellitus, hypertension, asthma, chronic obstructive pulmonary disease, chronic kidney disease, end-stage renal disease, congestive heart failure, coronary artery disease or cancer, triage vitals (including initial Spo 2 , heart rate, respiratory rate, mean arterial pressure, and temperature), emergency severity index, altered mental status at triage, initial labs (ferritin, lactate dehydrogenase, d-dimer, procalcitonin, anion gap, absolute lymphocyte count, absolute neutrophil count, creatinine, blood urea nitrogen, hemoglobin, platelet count, and Pco 2 in venous or arterial blood gas), presence of bilateral opacities or infiltrate on initial chest radiograph, and whether the patient received each of the following supplemental oxygen therapies: nasal cannula, nonrebreather mask, high-flow nasal cannula (HFNC), or noninvasive positive-pressure ventilation (NIPPV).

We carried out two main analyses: 1) the effect of intubation within 48 hours of triage on inhospital mortality, with the control group consisting of both patients that were never intubated and those that were intubated after 48 hours and 2) the effect of intubation at any point in time on mortality, with the control group consisting only of patients that were never intubated. We also carried out a sensitivity analysis comparing those that were intubated within 48 hours with those that were intubated at a later point in their hospital course, excluding patients that were never intubated.

We analyzed the rate of intubation and mortality over time by plotting the data and by calculating the per-week odds ratio (OR) for mortality and for intubation. We attempted to isolate the effect of time on intubation and mortality by adjusting this per-week OR for age, sex, race, ethnicity, hospital location, smoking status, and comorbidities. We analyzed the effect of admission hospital on rate of intubation by using mixed-effects logistic regression and treating admission hospital as a random effect (25) .

For all statistical tests, nominal statistical significance was set at α = 0.05. A post hoc power analysis was carried out, and power was estimated using the Schoenfeld model (26) . All statistical analyses were carried out using Version 4.0.0 of the R programming language (R Core Team, 2020, https://www.R-project. org/). Propensity score matching was carried out using the MatchIt package in R (Vienna, Austria) (27) .

Between March 1, 2020, and December 1, 2020, a total of 82,578 adult patients were tested for COVID-19 in the ED; of these tests, 12,902 (15.6%) were positive. Of these patients, 3,377 (26.2%) were discharged home, 487 (3.8%) were transferred to another facility outside the hospital system, 129 (1.0%) left against medical advice, 272 (2.1%) died before admission, and 8,637 (66.9%) were admitted. Of these, 8,510 (98.5%) represented a unique patient admission. Of these unique patients, 263 (3.1%) died within 48 hours of triage and were excluded. Of the remaining 8,247 patients, 650 (7.9%) had DNI orders placed within 48 hours of triage and were excluded, leaving 7,597 patients in the primary data set for this study. The agreement between the manually and automatically extracted data was substantial (see Supplementary Table 1 , http://links.lww.com/CCX/ A673). The overall flowchart for data acquisition and inclusion is summarized in Supplementary Figure 1 (http://links.lww.com/CCX/A674).

As of December 1, 2020, 1,628 of these patients (21.4%) had been intubated and 1,898 had died (25.0%). Of the 1,628 intubations, 807 (49.6%) occurred within 48 hours of triage. The distribution of days until death or discharge and days between triage and intubation are summarized in Figure 1 .

Intubation rates varied significantly from hospital to hospital, with as few as 11.4% of admitted patients being intubated to as many as 35.3%. After treating location as a random effect and treating the variables summarized above for the propensity score model as fixed effects, location was significantly associated with intubation in multilevel logistic regression (median OR, 1.60; p = 0.0013).

The dynamics of mortality and intubation over the course of the pandemic from March 1, 2020, to December 1, 2020, are summarized in Figure 2 , and the associated data are included in Supplementary 

The association between intubation and mortality was robust to several sensitivity analyses (summarized in Supplementary Table 4 , http://links.lww.com/CCX/ A673). Notably, not excluding patients with DNI orders and restricting only to patients with significant oxygen requirements (those receiving HFNC, NIPPV, or mechanical ventilation) and/or those admitted to the ICU did not change the direction of the effect. In addition to the effect of intubation within 48 hours, we found intubation at any time point to still be associated with increased mortality after matching (HR, 1.62 [1.45-1.80]; p < 0.0001).

Notably, we also found that when restricting to intubated patients (n = 1628), intubation within 48 hours of triage (n = 807) was not associated with significantly increased mortality compared with those intubated later in their hospital course (HR, 1.09 [0.94-1.26]; p = 0.26). 

In this work, we examined patterns of invasive ventilation in COVID-19 patients admitted during the New York City coronavirus outbreak. Previous studies have focused on recommendations for safe intubations (28) , reported data regarding intubation success and safety (29, 30) , or analyzed the effect of timing of intubation on mortality in COVID-19 patients (13). To our knowledge, this is the first study that has directly compared outcomes for intubated and nonintubated COVID-19 patients using propensity score matching.

We found rates of intubation to vary significantly from hospital to hospital, with providers at certain hospitals having intubated significantly more patients than others, even after adjusting for a variety of available confounders such as age or markers of disease severity. This difference may be due to local practice standards (such as lower clinical thresholds to intubate) that may have been amplified by the uncertainty of a novel viral pandemic. We found intubation to have decreased significantly over time, alongside mortality, and relationships that again remained stable after controlling for potential confounding variables. A portion of this change may be due to less severe cases presenting later in the pandemic; some amount, however, may be due to changes in clinical practice as understanding of the disease and treatment options grew. In the early stages of the pandemic, it was uncertain which patients would benefit from early intubation and mechanical ventilation. Some researchers and groups argued for early intubation (3, 13, 31) , whereas others cautioned against it (2, 7, 32) . The potential benefit of adjunctive airway treatments such as HFNC or NIPPV (33, 34) or awake proning (35, 36) was uncertain. Hence, even within one hospital system and after controlling for confounders, there was significant variation in rates of intubation.

After propensity score matching, we found intubation to be associated with increased mortality for COVID-19 patients, and this relationship remained robust to several sensitivity analyses. Intubation and mechanical ventilation carry a variety of risks. Early complications of tracheal intubation include cardiovascular collapse (37) , hypoxemia (38) , and aspiration (39) . Longer term complications include ventilatorinduced lung injury (40) and ventilator-associated pneumonia (41) . Mortality in intubated COVID-19 patients is remarkably high (42) . There have been reports of increased rates of barotrauma (43) , high rates of ventilator-associated pneumonia (44) , as well as severe acute kidney injury resulting from intubation in these patients (45) . Additionally, mechanically ventilated patients require substantially increased nursing and provider time, resources that may become increasingly limited during a pandemic surge (46) . A recent study found that greater ICU patient load was associated with increased mortality for COVID-19 patients (47) .

Notably, we did not find a significant association between mortality and timing of intubation (within 48 hr of triage vs later in hospital course), when restricting only to intubated patients. This differs from the results of Hyman et al (13) , who found that earlier intubation was associated with significantly reduced mortality, with each additional day increasing the hazard of death by 3%. The differences between our results and theirs may be due to baseline differences in populations, differing clinical decision thresholds for intubation, the inclusion of early laboratory markers of disease severity in our propensity score model (such as ferritin and d-dimer), or different study timeframes (from March 1, 2020, to December 1, 2020, in our study vs from January 30, 2020, to April 30, 2020, in theirs). Their study also focused on timing of intubation and did not compare intubated patients with nonintubated patients. Further studies from other populations and locations investigating this issue are needed.

Our study has several advantages. We drew data from a large, heterogeneous population across multiple centers in a city with significant racial, ethnic, and socioeconomic diversity. We used nearest neighbor propensity score matching and adjustment, which has theoretical advantages over multivariate regression (48) . We carefully selected the covariates included in the propensity score model and achieved good covariate balance. We attempted to reduce the effect of immortal time bias (16, 17) by treating as exposure patients that were intubated within 48 hours of triage. We performed several sensitivity analyses and found the primary effect to be robust to them.

This work also has multiple limitations. As a retrospective observational study, unmeasured confounding cannot be eliminated, and the design precludes us from drawing causal conclusions about intubation and mortality. The association between intubation and mortality after propensity matching/adjustment may still be due to more severely ill patients being intubated. Although there was substantial agreement between automatically and manually extracted chart data, there was little available information on whether patients were receiving oxygen therapy at the exact time the initial vitals were recorded. Similarly, we did not have access to exact information about these therapies, such as the effective Fio 2 being delivered. By restricting to patients that were intubated within 48 hours of triage, the power of our study was reduced. However, one of the sensitivity analyses we performed did not restrict to intubation within 48 hours, and in post hoc power analysis, we still had sufficient power to detect the effects identified in this study. Our analysis did not use timevarying covariates, as we only used data available at the time of triage, or the exposure (intubation) and outcome (mortality). 

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

Ventilation techniques and risk for transmission of coronavirus disease, including COVID-19: A living systematic review of multiple streams of evidence

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

Clinical course and outcomes of critically ill patients with SARS-CoV-2 pneumonia in Wuhan, China: A single-centered, retrospective, observational study

Rethinking the early intubation paradigm of COVID-19: Time to change gears?

Basing respiratory management of COVID-19 on physiological principles

In defense of evidence-based medicine for the treatment of COVID-19 acute respiratory distress syndrome

COVID-19 pneumonia: Different respiratory treatments for different phenotypes?

Timing of intubation and clinical outcomes in adults with acute respiratory distress syndrome

Management of COVID-19 respiratory distress

et al; and the Emory COVID-19 Quality and Clinical Research Collaborative: Timing of intubation and mortality among critically ill coronavirus disease 2019 patients: A single-center cohort study

Timing of invasive mechanic ventilation in critically ill patients with coronavirus disease 2019

Timing of intubation and in-hospital mortality in patients with coronavirus disease 2019

Epic Systems Corporation: Epic Electronic Medical Records

Multiple imputation for missing data in epidemiological and clinical research: Potential and pitfalls

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

Immortal time bias in pharmaco-epidemiology

Regression models and life-tables

Logistic regression, survival analysis, and the Kaplan-Meier curve

A comparison of 12 algorithms for matching on the propensity score

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

On regression adjustment for the propensity score

Overlap weighting: A propensity score method that mimics attributes of a randomized clinical trial

Variable selection for propensity score models

A brief conceptual tutorial of multilevel analysis in social epidemiology: Using measures of clustering in multilevel logistic regression to investigate contextual phenomena

Sample-size formula for the proportional-hazards regression model

MatchIt: Nonparametric preprocessing for parametric causal inference

Asthma prolongs intubation in COVID-19

A prospective, observational, cohort study of airway management of patients with Covid-19 by specialist tracheal intubation teams

collaborators: Emergency tracheal intubation in 202 patients with COVID-19 in Wuhan, China: Lessons learnt and international expert recommendations

Consensus statement: Safe Airway Society principles of airway management and tracheal intubation specific to the COVID-19 adult patient group

Clinical significance of timing of intubation in critically ill patients with Covid-19: A multi-center retrospective study

Surviving sepsis campaign: Guidelines on the management of critically ill adults with coronavirus disease 2019 (Covid-19)

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

Early, awake proning in emergency department patients with Covid-19

Proseva Investigators: Prone position-induced improvement in gas exchange does not predict improved survival in the acute respiratory distress syndrome

Cardiac arrest and mortality related to intubation procedure in critically ill adult patients: A multicenter cohort study

Incidence and duration of continuously measured oxygen desaturation during emergency department intubation

The occurrence of aspiration pneumonia after emergency endotracheal intubation

The Acute Respiratory Distress Syndrome Network: Ventilation with lower tidal volumes as compared with traditional tidal volumes for acute lung injury and the acute respiratory distress syndrome

Ventilator-associated pneumonia

Patient characteristics and outcomes of 11 721 patients with coronavirus disease 2019 (COVID-19) hospitalized across the United States

Increased incidence of barotrauma in patients with COVID-19 on invasive mechanical ventilation

Factors influencing liberation from mechanical ventilation in coronavirus disease 2019: Multicenter observational study in fifteen Italian ICUs

Severe acute kidney injury in patients with COVID-19 and acute respiratory distress syndrome

Fair allocation of scarce medical resources in the time of Covid-19

Association of intensive care unit patient load and demand with mortality rates in US Department of Veterans Affairs Hospitals during the COVID-19 pandemic

Statistical primer: Propensity score matching and its alternatives

The authors have disclosed that they do not have any potential conflicts of interest.For information regarding this article, E-mail: parisha@nychhc. org; auparish@gmail.com