key: cord-0010973-v6tj958h authors: Peñuelas, Oscar; Muriel, Alfonso; Abraira, Victor; Frutos-Vivar, Fernando; Mancebo, Jordi; Raymondos, Konstantinos; Du, Bin; Thille, Arnaud W.; Ríos, Fernando; González, Marco; del-Sorbo, Lorenzo; Ferguson, Niall D.; del Carmen Marín, Maria; Pinheiro, Bruno Valle; Soares, Marco Antonio; Nin, Nicolas; Maggiore, Salvatore M.; Bersten, Andrew; Amin, Pravin; Cakar, Nahit; Suh, Gee Young; Abroug, Fekri; Jibaja, Manuel; Matamis, Dimitros; Zeggwagh, Amine Ali; Sutherasan, Yuda; Anzueto, Antonio; Esteban, Andrés title: Inter-country variability over time in the mortality of mechanically ventilated patients date: 2020-01-07 journal: Intensive Care Med DOI: 10.1007/s00134-019-05867-9 sha: fb8effe32fb79abcbb332960089ef6d3c5d3eead doc_id: 10973 cord_uid: v6tj958h PURPOSE: Variations in clinical characteristics and management and in the mortality of mechanically ventilated patients have not been sufficiently evaluated. We hypothesized that mortality shows a variability associated with country after adjustment for clinical characteristics and management. METHODS: Analysis of four studies carried out at 6-year intervals over an 18-year period. The studies included 26,024 patients (5183 in 1998, 4968 in 2004, 8108 in 2010, and 7765 in 2016) admitted to 1253 units from 38 countries. The primary outcome was 28-day mortality. We performed analyses using multilevel logistic modeling with mixed-random effects, including country as a random variable. To evaluate the effect of management strategies on mortality, a mediation analysis was performed. RESULTS: Adjusted 28-day mortality decreased significantly over time (first study as reference): 2004: odds ratio 0.82 (95% confidence interval [CI] 0.72–0.93); 2010: 0.63 (95% CI 0.53–0.75); 2016: 0.49 (95% CI 0.39–0.61). A protective ventilatory strategy and the use of continuous sedation mediated a moderate fraction of the effect of time on mortality in patients with moderate hypoxemia and without hypoxemia, respectively. Logistic multilevel modeling showed a significant effect of country on mortality: median odds ratio (MOR) in 1998: 2.02 (95% CI 1.57–2.48); in 2004: 1.76 (95% CI 1.47–2.06); in 2010: 1.55 (95% CI 1.37–1.74), and in 2016: 1.39 (95% CI 1.25–1.54). CONCLUSIONS: These findings suggest that country could contribute, independently of confounder variables, to outcome. The magnitude of the effect of country decreased over time. Clinical trials registered with http://www.clinicaltrials.gov (NCT02731898). ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (10.1007/s00134-019-05867-9) contains supplementary material, which is available to authorized users. The management of mechanically ventilated patients has improved in recent years with the incorporation of strategies aimed at preventing or reducing ventilatorinduced lung injury: use of non-invasive positive-pressure ventilation (NPPV) [1] [2] [3] , a lung-protective strategy [4] [5] [6] , titration of positive end-expiratory pressure [7] , prone positioning [8] , and neuromuscular blocking agents [9] , and at decreasing the duration of mechanical ventilation through the reduced use of sedatives [10, 11] and the early identification of patient readiness to breathe unassisted [12] . Some of these interventions which initially focused on the management of patients with acute respiratory distress syndrome now appear to be substantially applicable to all mechanically ventilated patients [6, 13] . Previous observational data [14] [15] [16] have shown that short-term mortality of mechanically ventilated patients has decreased over time. Nevertheless, several aspects in the management of critically ill patients have been shown to differ across intensive care units (ICUs) and countries, leading to considerable clinical variability [17] . Variability in clinical practice is a current and controversial topic, which is of interest not only to clinicians, but also to health administrators and patient advocates. Some variations may be explained by the characteristics of individual patients or by differences in the ability of doctors who rely too much on subjective judgement. However, unwarranted clinical variability may exist, defined as variations in patient outcomes that cannot be explained by the patient's underlying illness or patient care based on evidence-based medicine [18] . The existence of unwarranted variability may have a negative impact on quality of care and clinical outcomes [19] . Geographical variability in outcome has been reported in several conditions such as ischemic heart disease [20, 21] , cancer [22] [23] [24] , and cardiac arrest [25] . However, the geographical or inter-country variability in mortality in critically ill patients has been little evaluated. In this study, we estimate the inter-country variability in the outcomes of mechanically ventilated patients through the analysis of four international observational studies carried out at 6-year intervals from 1998 to 2016. We hypothesized that, after adjustment for variables related to patient management and clinical characteristics, mortality would be found to have decreased over time and to show a variability related to country. We carried out four prospective, observational, multicenter international studies, in 1998 [14] , 2004 [15] , 2010 [16] , and 2016 (first report), in adults who had need of invasive mechanical ventilation for more than 12 h, or received NPPV for more than 1 h from ICU admission. The decision to use NPPV was based on the local practices of the units participating in the study. We followed the same methodology for all studies (descriptions of the variables collected in each study are provided in eTable 1). Briefly, in all studies, we collected data on baseline characteristics [age, sex and severity at admission estimated by Simplified Acute Physiology Score (SAPS II) which ranges from 0 (lower severity) to 163 (higher severity)], daily gas exchange, variables related to management ventilator settings, sedation, neuromuscular blockers, and complications [acute respiratory distress syndrome (ARDS), sepsis, ventilator-associated pneumonia]. Evaluation of organ function (cardiovascular, renal, hepatic, hematologic) was performed according to the sequential organ failure assessment (SOFA) score and organ failure was defined as a SOFA sub-score higher than two points for the organ in question. These data were collected while patients were ventilated or until day 28. Ventilator-free days were calculated according to a previously published formula [26] . Patients were followed until discharge from hospital. In the event of discharge from hospital before day 28 after starting mechanical ventilation, the status on this day was registered. National coordinators recruited local investigators from eligible ICUs (see full list of Investigators in the Electronic Supplementary Material). Only the research team members at each site were aware of the purpose and the precise timing of the study. Ethics Committees of each participating institution approved the protocol; a waiver of informed consent was obtained according to local regulations. For the 1998 and 2004 studies, all investigators and research coordinators were provided with a manual describing data collection and definitions, and case report forms on paper. Investigators filled out a case report form for each included patient and sent it to the Coordinating Center to be included in an electronic database. Three study coordinators checked each case report form to identify errors, and inconsistent data were corrected. In the 2010 and 2016 studies, detailed instructions were available through a secure website. Validity checks were made at the time of electronic data entry, including plausibility checks within each variable and between variables. At the time of data entry, the site investigators were required to answer all queries raised by the case report form before they could electronically complete each individual patient data set. Patient data sets that were not complete were not included in the analysis. In addition, prior to analysis, all data were screened for potentially erroneous data and outliers. These data were verified and corrected by site investigators. There was no on-site data monitoring. We followed the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement guidelines for observational cohort studies [27] . The main objective of the analysis was to evaluate, over time, the inter-country variability in 28-day mortality The mortality of mechanically ventilated patients could show country-related variability, albeit decreasing over time. after the start of mechanical ventilation. The first step was to estimate the crude temporal trends of mortality over the four periods by creating a logistic regression model with study year, coded as a dummy variable, as the only independent variable, and 28-day mortality as the dependent variable. To adjust this estimate, we then fitted a multilevel (patients and countries) logistic regression model accounting for baseline demographic variables (age, sex, SAPS II, and reason for initiation of mechanical ventilation); variables related to the evolution of patients over the course of mechanical ventilation [development of organ failure (cardiovascular, renal, hepatic, and hematologic) defined as a SOFA score higher than two points]; complications during the course of mechanical ventilation, such as acute respiratory distress syndrome, sepsis, and ventilator-associated pneumonia; variables related to the ventilatory setting [use of NPPV, a lung-protective strategy (arbitrarily defined as tidal volume below 6 ml/kg actual body weight or tidal volume below 8 ml/ kg actual body weight plus plateau or peak inspiratory pressure less than 30 cmH 2 O)]; and variables related to ventilator management (use of sedation, neuromuscular blocking, and performance of tracheotomy). The random part of the model included two effects and their possible correlation: a random effect in the intercept modeling heterogeneity of mortality among countries in the first study, and a random effect in the slope modeling heterogeneity of changes in mortality over time among countries. The model was estimated by maximum likelihood using the adaptive Gaussian quadrature approximation (with seven quadrature points) [28] . The likelihood ratio test was used to assess the goodness of fit and to compare alternative models. The measure of fixed effects was the odds ratio with 95% confidence interval. Second, to estimate the random variability of outcomes between countries, we calculated the median odds ratio (MOR) defined as the median value of the odds ratio between the country at highest risk of mortality and the area at lowest risk when randomly picking out two countries. The MOR can be conceptualized as the median increased risk of mortality that a patient would have if he or she moved to a country with a higher risk [29] . For the purpose of this analysis, countries that included fewer than 30 patients were excluded to avoid spurious findings in reliabilityadjusted models. We performed two sensitivity analyses: first, using data from the countries that participated in all four studies, and second, including only patients with invasive ventilation from the beginning of ventilatory support. A secondary analysis was performed to estimate changes, over time, in the practice of three management strategies: use of NPPV as the first mode of ventilatory support, use of a lung-protective strategy on the first day of mechanical ventilation, and use of sedation (defined as a daily continuous infusion for more than three consecutive hours of any sedative drug) over the course of mechanical ventilation. For this purpose, a multivariate logistic regression model was performed for each management strategy. Variables entered in these models were: age, sex, SAPS II, and reason for starting mechanical ventilation. We calculated and plotted the margin probability for each strategy in each study year with the estimators obtained from the model. To estimate how these changes influence mortality, a mediation analysis was performed. The mediation analysis assumes a causal relationship between the exposure (period of time) and mediator variables (NPPV, sedation, lung-protective strategy) which in turn may affect the outcome (28-day mortality), and it clarifies the nature of the relationship between the exposure and mediator to better understand the pathways by which the exposure potentially affects the outcome [30] . Therefore, we performed the mediation analysis using the inverse odds ratio weighting method to estimate the direct effect of time on 28-day mortality and the effect, on 28-day mortality, that is mediated by the use of NPPV, sedation and the lung-protective strategy (indirect effect). We calculated direct and indirect effects of each period using inverse odds weighting, which is a semi-parametric weight-based approach to mediation analysis, implemented with standard software. This method can accommodate multiple mediators [31, 32] . Standard errors for the indirect effect were estimated using bootstrapping. For this mediation analysis, we considered two clinical approaches: first, in the overall cohort of mechanically ventilated patients; second, taking into account that the application of the above strategies cannot be uniform, since they depend basically on the severity of respiratory failure at initiation of mechanical ventilation. For this reason, we classified the overall cohort into four sub-groups: (a) patients with a PaO 2 /FiO 2 ratio higher than 300 at the beginning of mechanical ventilation; (b) patients with a PaO 2 /FiO 2 ratio of between 200 and 300; (c) patients with a PaO 2 /FiO 2 ratio of between 100 and 200, and (d) patients with a PaO 2 /FiO 2 ratio lower than 100. Statistical analyses were conducted using Stata Software 14.0 (StataCorp LP, College Station, Texas). Overall 27,467 of 63,313 patients (43%) admitted to participating ICUs met the inclusion criteria in the four periods of study. After excluding patients from countries with fewer than 30 patients, full data for analysis were available for 26,024 patients: 5183 patients in 1998, 4968 patients in 2004, 8108 patients in 2010, and 7765 patients in 2016 (eFigure 1). We included patients from 38 countries (eTable 2). The sample size of enrolled patients for each participating country in each study is shown in eTable 3. Fifteen countries (Argentina, Canada, Chile, Colombia, Ecuador, France, Greece, Italy, Mexico, Peru, Portugal, Spain, Tunisia, USA, and Uruguay) had patients in all four studies and accounted for 63.5% of the overall cohort (16,537 of 26,024 patients). Comparisons of baseline patient characteristics across studies are shown in Table 1 . Over time, the most common patient type was consistently found to be male, in the sixth decade of life, with a severity at admission that predicted a hospital mortality around 35%. Changes in the reason for mechanical ventilation were observed, with a decrease in patients with chronic obstructive pulmonary disease (COPD) (from 10% in 1998 to 7% in 2016; p < 0.001) and cardiac failure (from 10% in 1998 to 7% in 2016; p < 0.001) and an increase in patients with neurological disease (from 17% in 1998 to 21% in 2016; p < 0.001). Data on evolution and outcomes over the course of mechanical ventilation are shown in Table 2 . It emerged that there was, over time, a significant increase in the use of the three management strategies, which were found to influence the outcome of ventilated patients: NPPV, lung-protective strategy, and continuous sedation (Table 3 and eFigure 2). The implementation of NPPV as first ventilatory support increased from 5% in 1998 to 16% in 2016 (p < 0.001) although without significant changes in the failure rate (around 30% of patients required intubation in all studies). The most important change in the management of mechanically ventilated patients was a decrease in the tidal volume (from a mean of 8.8 ml per kg of actual body weight in 1998 to 6.7 ml per kg of actual body weight in 2016; p < 0001) and an increment in the applied positive end-expiratory pressure (PEEP) (from a mean of 5 cm of water in 1998 to 7 cm of water in 2016; p < 0.001). These changes were observed in both patients with ARDS and patients without ARDS (Fig. 1) . Also, this evolution implied changes in other relevant variables such as plateau pressure and driving pressure (eFigure 3). A mixed-effects logistic regression showed the magnitude of these changes, taking as reference the first study, and the inter-country variability (Table 4 ). Raw 28-day mortality by country and study year is shown in eFigure 4. On average, inter-country variability decreased over time, with countries more distant from the mean changing more in the opposite direction, as shown in Fig. 2 . The bivariate analysis and fitted multilevel logistic regression are shown in eTable 4. The results of the mixed-effect analysis (eTable 5) showed an independent effect associated with country as estimated by the MOR: 2.02 in 1998 (95% CI 1.57-2.47); 1.76 (95% CI 1.47-2.05) in 2004; 1.55 in 2010 (95% CI 1.36-1.73); and 1.39 (95% CI 1.24-1.53) in 2016 (Fig. 3) . Also, the sensitivity analysis, including the 15 countries that participated in all four studies, showed similar results in the adjusted model: MOR (Fig. 3) . The sensitivity analysis including only patients with invasive ventilation from beginning of ventilatory support showed following results: MOR (Fig. 3) . No effect of mediation variables was found in the overall cohort of patients (eTable 6). However, NPPV, a lung-protective strategy, and the use of sedation showed different effects on mortality according to severity of hypoxemia (eTable 7). In patients with hypoxemia, implementation of a lung-protective strategy had an indirect effect over Table 2 We showed that 28-day mortality in an overall cohort of mechanically ventilated patients decreased significantly over time. Our analysis suggests that mortality could show an effect of country, independently of patient and management factors. However, some management strategies (lungprotective strategy and sedation) mediated an important fraction of the effect of temporal evolution on mortality. Over the past 20 years, an extensive body of research in different areas has led to advances in mechanical ventilation to improve oxygenation and prevent ventilatorinduced lung injury [33] . The question is whether or Tidal Volume (ml/kg actual body weight) Applied PEEP (cm of water) Fig. 1 Changes over time in tidal volume (blue box plot) and positive end-expiratory pressure (red box plot) set at initiation of mechanical ventilation in patients without criteria of acute respiratory distress syndrome (ARDS) and patients with ARDS. For each box plot, the middle line represents the median, the lower hinge represents the first quartile, the upper hinge represents the third quartile, the whiskers extend to 1.5 times interquartile range, and the outliers are values outside the whiskers' range not the implementation of better practices in the management of mechanically ventilated patients [1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] has spread homogeneously across the world and whether these strategies themselves have an impact on mortality. We think that a reasonable way to address these issues could be through a proper estimate of the influence of country on mortality rates in mechanically ventilated patients. Our objective was to evaluate whether the observed variations in mortality of mechanically ventilated patients have an inter-country effect after proper adjustment for variables associated with clinical management. However, our aim was not to make an explicit comparison between individual countries or geographical areas, as other authors have done [17] . Interestingly, we found that the variability in 28-day mortality may be partially explained by a random effect of country and that this effect showed a significant reduction over time. Our results also showed a significant geographical variability in 28-day mortality among countries, with the risk of death found to differ up to twofold between any two randomly selected countries. The observed random effect of the international variability in mortality could, itself, be influenced by possible unmeasured variables such as ICU organization and utilization [34, 35] , intensive-care resources [36] , intensivecare admissions policy, intensivist staffing schedules [37] , Table 4 Odds ratio for the implementation of management strategies (the 1998 study was taken as reference) and the inter-country variability a Defined as a tidal volume below 6 ml/kg actual body weight or a tidal volume below 8 ml/kg actual body weight plus plateau or peak inspiratory pressure less than 30 2 Regression lines of the evolution of the 28-day mortality by countries (thin lines) and globally (thick line). For simplicity of the graph, only the fifteen countries that participated in the four studies are shown. Overall 28-day mortality decreased over time (thick line). The variability between countries declined (the dispersion is lower in the last year than in the first year). The rate of change in mortality in each country is inversely related to the initial dispersion: countries farthest from the average changed faster in the opposite direction. Countries with initial higher mortality tended to see this decline more rapidly (blue lines), while those with initial lower mortality tended to raise it slowly (red lines) or variability in the withdrawal of life-sustaining treatment across countries [38, 39] . Additional studies are needed to evaluate whether all these variables, together with the variables assessed in our study, could contribute to that effect. We speculate that progressive and overall implementation of the best available clinical practice may have led to this variability in mortality [40] . We could argue that the countries with worse results in the past have improved their clinical practice through higher rates of adherence to best practices more experience [41] , or an economic improvement in the case of the low-or middle-income countries if we consider the geo-economy a variable associated with the outcome of the overall cohort of mechanically ventilated patients as described in patients with ARDS [17] . However, it is possible to argue that the implementation of certain management strategies, such as NPPV, protective ventilation and sedation, could, per se, be the reason for the observed changes in mortality. A metaanalysis that included individual patient data from nine randomized clinical trials showed that reductions in the daily fluid balance, plateau pressure and tidal volume, and the increase in PEEP decreased the mortality rate in ARDS patients over a period of 17 years [42] . We performed a mediation analysis to evaluate the effect of the aforementioned management strategies. Despite the methodological complexity, we have found a moderate (lower than 25%) mediated effect of lung-protective strategy and sedation in a subgroup of patients with moderate hypoxemia. Our study has limitations. Despite enrolling a large number of ICUs from around the world, our sample may not have avoided a selection bias liable to limit the generalizability of its findings. However, given that 60% of our sample was collected from countries that participated in all four studies, as well as the fact that it covered a time span of almost 20 years, its results can be considered very consistent and reliable. It may also be questioned whether the participating ICUs are sufficiently representative of the countries. We are confident that the availability and use of health resources within each country, including the socioeconomic and demographic aspects of patients, showed similar random distribution. Another limitation is that we did not have all the same ICUs participating in all four studies. To some extent, this was unavoidable due to closure of some previously participating hospitals, changes in resources available for data collection at participating centers, and our goal to expand the representability of the results by including ICUs from more countries. Finally, we had missing data for some relevant variables; however, we believe that the size of the database should have minimized the influence of the missing information in the results. In this sense, we analyzed the evolution over time of tidal volume expressed as the tidal volume in milliliters per kilogram The median odds ratio can be conceptualized as the increased risk of mortality that a patient would have if he/she moved to a country with different mortality rate. Adjusted inter-country variability estimated by mean odds ratio (95% confidence interval) at each study year in all the participating countries (red) and in the 15 countries that took part in all four studies (blue). The Median Odds Ratio (MOR) is defined as the median value of the odds ratio between the area at highest risk of mortality and the area at lowest risk when randomly picking out two areas of actual body weight, because height was not registered in the first study. Some may argue that the data on tidal volumes are difficult to interpret given the use of actual body weight instead ideal body weight. However, the differences between actual body weight and ideal body weight were found to decrease slightly over time (13. 69 kg in 2004, 12.91 kg in 2010, and 11.97 kg in 2016), which resulted in only minor differences between tidal volume estimated by actual or predicted body weight (1.46 ml in 2004, 1.14 in 2010, and 1.13 in 2016). Therefore, we consider that this possible limitation does not invalidate our observation about the evolution of tidal volume. In conclusion, we have provided clinical evidence on the observed decrease in mortality of mechanically ventilated patients over a period of 18 years. We observed a variability associated with country, which also decreased over time, but could be partially mediated by the implementation of some management strategies in patients with moderate hypoxemia. Effect of noninvasive ventilation on tracheal reintubation among patients with hypoxemic respiratory failure following abdominal surgery: a randomized clinical trial Effect of noninvasive ventilation vs oxygen therapy on mortality among immunocompromised patients with acute respiratory failure: a randomized clinical Trial Changing use of noninvasive ventilation in critically ill patients: trends over 15 years in francophone countries Association between use of lung-protective ventilation with lower tidal volumes and clinical outcomes among patients without acute respiratory distress syndrome: a meta-analysis Effect of a low vs intermediate tidal volume strategy on ventilator-free days in intensive care unit patients without ARDS: a randomized clinical trial Lung protective ventilation strategy for the acute respiratory distress syndrome Higher vs lower positive end-expiratory pressure in patients with acute lung injury and acute respiratory distress syndrome: systematic review and meta-analysis Prone positioning in severe acute respiratory distress syndrome Neuromuscular blockers in early acute respiratory distress syndrome Early intensive care sedation predicts long-term mortality in ventilated critically ill patients Daily sedation interruption in mechanically ventilated critically ill patients cared for with a sedation protocol: a randomized controlled trial Efficacy and safety of a paired sedation and ventilator weaning protocol for mechanically ventilated patients in intensive care (Awakening and Breathing Controlled trial): a randomised controlled trial The future of mechanical ventilation: lessons from the present and the past Characteristics and outcomes in adult patients receiving mechanical ventilation: a 28-day international study Evolution of mechanical ventilation in response to clinical research Evolution of mortality over time in patients receiving mechanical ventilation Geo-economic variations in epidemiology, patterns of care, and outcomes in patients with acute respiratory distress syndrome: insights from the LUNG SAFE prospective cohort study Disseminating innovations in health care Unwarranted variations in healthcare delivery: implications for academic medical centres Mechanical ventilation and clinical practice heterogeneity in intensive care units: a multicenter case-vignette study Geographic variability in patient characteristics, treatment and outcome in an International Trial of Magnesium in acute myocardial infarction Assessment of the magnitude of geographical variations and socioeconomic contextual effects on ischaemic heart disease mortality: a multilevel survival analysis of a large Swedish cohort Hospital and geographic variability in two colorectal cancer surgery outcomes: complications and mortality after complications Geographic variations in lung cancer lobectomy outcomes: the general thoracic surgery database Hospital variation in survival after in-hospital cardiac arrest Statistical evaluation of ventilator-free days as an efficacy measure in clinical trials of treatments for acute respiratory distress syndrome Strengthening the reporting of observational studies in epidemiology (STROBE) statement: guidelines for reporting observational studies Multilevel and longitudinal modelling using stata A brief conceptual tutorial of multilevel analysis in social epidemiology: using measures of clustering in multilevel logistic regression to investigate contextual phenomena Mediation analysis allowing for exposuremediator interactions and causal interpretation: theoretical assumptions and implementation with SAS and SPSS macros Inverse odds ratio-weighted estimation for causal mediation analysis Practical guidance for conducting mediation analysis with multiple mediators using inverse odds Ratio weighting Mechanical ventilation: state of the art The impact of hospital and ICU organizational factors on outcome in critically ill patients: results from the extended prevalence of infection in intensive care study Hospital volume and the outcomes of mechanical ventilation Hospital mortality in relation to staff workload: a 4-year study in an adult intensive-care unit The effect of nighttime intensivist staffing on mortality and length of stay among intensive care unit patients Withholding and withdrawal of life-sustaining treatments in intensive care units in Asia Global variability in withholding and withdrawal of life-sustaining treatment in the intensive care unit: a systematic review Five-year trends of critical care practice and outcomes Acute respiratory distress syndrome cases volume and ICU mortality in medical patients Declining mortality in patients with acute respiratory distress syndrome: an analysis of the acute respiratory distress syndrome network trials To all national coordinators and participants in the four International Studies on Mechanical Ventilation listed in Electronic Supplementary Material. The authors declare there is no conflict of interest. Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.