This study was approved by the Ethical Committee of West China Hospital (WCH2018-409). It was reported according to the REporting of studies Conducted using Observational Routinely-collected Data (RECORD) Checklist18.
Data source
We conducted a retrospective cohort study based on two large-scale validated ICU databases: the Registry of Healthcare-associated Infections in ICUs in China (ICU-HAI) and the Medical Information Mart for Intensive Care (MIMIC) -IV database.
The ICU-HAI registry database included all patients admitted to six ICUs (general, surgical, respiratory, neurological, thoracic, and pediatric) in the West China Hospital (WCH) system. The WCH system is a national critical care centre. There were over 8300 ICU admissions per year, and approximately 93% of them had undergone MV. The ICU-HAI registry was established by integrating data from three systems (electronic medical record (EMR) system, ICU-HAI system, and ICU system) using the patient’s unique ID. This registry contained over 1000 variables with detailed information on demographics, laboratory and microbiology results, diagnosis, treatment, outcome, hospital costs, etc. Among these, many variables were recorded hourly by well-trained specialist nurses, such as the Richmond Agitation-Sedation Scale (RASS), vital signs, fluid intake and output volume, and life support (oxygen, mechanical ventilation, positive end-expiratory pressure, and fraction of inspired oxygen). The ICU-HAI registry has been validated and proven to be of high quality. It has supported various ICU-related studies 19–22. A detailed database profile of the ICU-HAI registry was described in our previously published study 23. In this study, we used data from April 2015 to December 2018, which contained information on 28,848 ICU patients.
MIMIC-IV is sourced from two data systems (an EHR system and an ICU system) of the Beth Israel Deaconess Medical Center in Boston, Massachusetts. The data profile of MIMIC-IV is available elsewhere24. In this study, MIMIC-IV v2.0 was used, containing 76,943 ICU admissions of 69,639 patients from 2008 to 2019, with abundant clinical and intensive care information.
Study population and sedation strategies
We identified ICU patients who consecutively received MV for at least three days from the ICU-HAI registry and MIMIC-IV. We excluded patients who were missing basic information such as age, sex and diagnosis at discharge, or who were under 18 years of age. To reduce indication bias, we also excluded patients who underwent MV for less than three days, stayed in the ICU for more than 30 days, or were in a coma (RASS = -5) before MV treatment.
We defined three sedation strategies according to the daily RASS within 7 days after the initiation of MV, and the strategies included continuous light sedation (CLS), light to deep sedation (LTDS) and continuous deep sedation (CDS). CLS referred to light sedation within 2 days after initiation of MV followed by nondeep sedation; LTDS referred to deep sedation within 2 days after initiation of MV followed by nondeep sedation; and CDS was defined as continuous deep sedation after MV initiation. We categorised the depth of sedation according to the mean daily value of the RASS: deep sedation (RASS≦-3), light sedation (-3 < RASS≦--1), calm (-1 < RASS < 1) and restless/agitated (RASS ≧ 1) 1,25.
Study outcomes
Our outcomes included ventilator mortality, ICU mortality and hospital mortality. Patients who died one day after extubation were defined as ventilator mortality, indicating extubation failure. Similarly, patients who died one day after ICU discharge were defined as ICU mortality.
Covariates
Potential confounding covariates included demographic information (age and sex), ICU type, Acute Physiology and Chronic Health Evaluation (APACHE) II score at ICU admission, comorbidities (diabetes, ischaemic heart disease (IHD), heart failure, kidney failure, renal failure, pulmonary vascular disease, hypertension), acute conditions (acute respiratory distress syndrome (ARDS), gastrointestinal bleeding, shock, pneumonia, sepsis, trauma), and treatments. Treatments included prescriptions (sedatives, neuroleptics, opioids, antithrombotics, neuromuscular blockers, acid inhibitors, antibiotics, intestinal probiotics, immunosuppressants, vasopressors), surgery (cranial or cardiac surgery) and other interventions (fibreoptic bronchoscopy, tracheotomy, enteral nutrition, gastrointestinal decompression and head-of-bed elevation).
For the MIMIC dataset, ethnicity (white, Asian, black, Latin American and other) was also taken into account. The MIMIC dataset did not collect the APACHE II score, instead, we used two other clinical scores, the Oxford Acute Severity of Illness Score (OASIS) and Acute Physiology Score III (APSIII), to assess the health status of patients on ICU admission.
Demographic information, type of ICU, clinical scores on ICU admission and comorbidity/acute conditions were measured as time-fixed variables. Treatments were measured as time-dependent variables, which were recorded daily from admission to discharge.
Statistical analyses
All analyses were conducted separately in the ICU-HAI and MIMIC cohorts, and the adjusted hazard ratios (HRs) of three sedation strategies were pooled through meta-analysis.
We described the sedation pattern of patients who received consecutive MV for at least 3 days by plotting the individual trajectory of the sedation depth within 30 days of MV initiation. The patients were classified into 3 categories according to the clinical outcome within 30 days: dead, successfully extubated and still requiring ventilatory support.
We compared the characteristics and clinical outcomes among patients who underwent three sedation strategies of interest using descriptive analyses and univariate tests. We depicted numeric variables as the means (standard deviation) and categorical variables as frequencies (percentages). The Student’s t test was applied to assess differences in the means. The \({\chi }^{2}\) test or Fisher’s exact test was used to assess differences in percentages.
We assessed the effects of the three sedation strategies on ventilator mortality, ICU mortality and hospital mortality. To control the time-dependent confounding affected by time-dependent sedation status, we conducted a weighted survival analysis 26. This analysis involved three steps. First, we calculated the daily probability of patients being exposed to different sedation depths by estimating the propensity score (PS) for multivalued treatments. As the ICU data were highly complex with nonlinearity, interactions, etc., the XGBoost algorithm was used to generate the PS. Second, we constructed the inverse probability of treatment weights based on the cumulative probability of the sedation trajectory during MV 27. Finally, a Fine-Gray subdistribution hazard competing risk model combined with inverse probability weighting (IPTW) was adopted to analyse the competing risk of extubation alive versus ventilator mortality, ICU discharge alive versus ICU mortality, and hospital discharge alive versus hospital mortality28. Based on the Fine-Gray model weighted by IPTW, marginal structural models (MSMs) were used to assess the risk of sedation strategies on ventilator mortality, and doubly robust models were applied to analyse ICU mortality and hospital mortality with adjustment for confounders both during MV and after extubation 29. The sedation effects for each mortality outcome estimated in the ICU-HAI and MIMIC cohorts were pooled using fixed-effect meta-analysis. Missing data were imputed by the last observation carried forwards.
Sensitivity analyses
We conducted six sensitivity analyses to test the robustness of our results: 1) we combined the estimated effects for each outcomes through random effect meta analysis; 2)we adopted multiple imputation to tackle missing data; 3) we assessed the association using complete cases without imputation; 4) we used an alternative popular boosting algorithm, namely the gradient boosting machine (gbm), for weighting; 5) we employed a time-dependent Cox model regardless of the competing risk; and 6) we used an alternative definition of deep to light sedation strategy, that is deep sedation within 1 day of MV initiation and then switching to nondeep sedation. In the last five sensitivity analyses, the estimated effects of ICU-HAI cohort and MIMIC cohort were combined by fix-effect meta analysis same as the main analysis.
We used R packages twang, RISCA and survival, version 4.2.2. A two-sided P < 0.05 was considered to be statistically significant. Data analysis was conducted from August 2022 to April 2023.