Abstract
ARDS is severe form of respiratory failure with significant impact on the morbidity and mortality of critical care patients. Epidemiological data are crucial for evaluating the efficacy of therapeutic interventions, designing studies, and optimizing resource distribution. The goal of this review is to present general aspects of mortality data published over the past decades. A systematic search of the MEDLINE/PubMed was performed. The articles were divided according to their methodology, type of reported mortality, and time. The main outcome was mortality. Extracted data included study duration, number of patients, and number of centers. The mortality trends and current mortality were calculated for subgroups consisting of in-hospital, ICU, 28/30-d, and 60-d mortality over 3 time periods (A, before 1995; B, 1995–2000; C, after 2000). The retrospectivity and prospectivity were also taken into account. Moreover, we present the most recent mortality rates since 2010. One hundred seventy-seven articles were included in the final analysis. General mortality rates ranged from 11 to 87% in studies including subjects with ARDS of all etiologies (mixed group). Linear regression revealed that the study design (28/30-d or 60-d) significantly influenced the mortality rate. Reported mortality rates were higher in prospective studies, such as randomized controlled trials and prospective observational studies compared with retrospective observational studies. Mortality rates exhibited a linear decrease in relation to time period (P < .001). The number of centers showed a significant negative correlation with mortality rates. The prospective observational studies did not have consistently higher mortality rates compared with randomized controlled trials. The mortality trends over 3 time periods (before 1995, 1995–2000, and after 2000) yielded variable results in general ARDS populations. However, a mortality decrease was present mostly in prospective studies. Since 2010, the overall rates of in-hospital, ICU, and 28/30-d and 60-d mortality were 45, 38, 30, and 32%, respectively.
Introduction
ARDS is a highly heterogeneous life-threatening severe form of acute respiratory failure. It is characterized by an acute onset of hypoxemia refractory to oxygen treatment and the presence of bilateral pulmonary infiltrates on chest radiograph.1–4 The syndrome is based on alveolar-capillary membrane injury followed by the development of pulmonary edema, which cannot be fully explained by heart failure.1–3 ARDS is usually accompanied by multiple- organ system failure and therefore substantially contributes to high morbidity and mortality in critical care.5,6 It increases the cost of treatment of critical care patients by prolonging ICU and in-hospital stay. Reliable epidemiologic data concerning ARDS are important for all critical care physicians. These data facilitate clinical decision making, with implications for the efficacy of current therapeutic interventions, and suggest ways for their future development.7 Moreover, the data help us to better assess the present clinical and economic impact of the condition and may lead to more appropriate resource distribution.
The history of ARDS should be divided into 3 main periods: (1) before 1995 (before publication of the American-European Consensus Conference [AECC] criteria), (2) from 1995 to 2000 (implementation of the AECC criteria and development of lung-protective mechanical ventilation), and (3) from 2001 to the present (the era of lung-protective mechanical ventilation).
The main goals of this systematic review of the literature are to highlight the complexity of ARDS mortality reporting, to outline the trends over the last 2 decades, and to estimate current mortality rates of ARDS, suggesting benchmarks in various clinical trial designs. Moreover, uncertainties and controversies regarding the topic are assessed, and possible explanations are presented.
Methods
We used the PRISMA (Preferred Reporting Items for Systematic Reviews and Meta-Analyses) statement as a guide.8 A computerized search through the online web database MEDLINE/PubMed was performed. The search was performed by the MeSH (Medical Subject Headings) system using the key words “acute lung injury” AND “acute respiratory distress syndrome” AND “mortality” AND “outcome” AND ”incidence” AND “clinical trial.” The bibliographies of the retrieved articles were also reviewed. Human studies with >20 subjects published from January 1994 to January 2015 were taken into consideration. Non-English articles were included, provided that they contained all required data. Two independent reviewers (MJ and JO) screened the titles, abstracts, and methodological quality of the articles with a high degree of agreement (93.2%). All eligible articles were retrieved in full text (see online supplement I at http://www.rcjournal.com).
The data collected in a data sheet were as follows: study group characteristics (specific vs mixed ARDS population). We defined a mixed population as a group of subjects with ARDS irrespective of etiology. A group of subjects enrolled into the study primarily due to the presence of a unique diagnosis (ie, multiple trauma, sepsis, burns, military injury) or intervention (ie, extracorporeal oxygenation, high-frequency ventilation) was defined as the specific ARDS population. The specific population was omitted from final analysis to improve generalization of the results. Additional extracted data included number of centers, number of enrolled patients with ARDS, time of enrollment (in months), time period when the study was started and performed, type of mortality, and absolute value of the reported mortality rate. The characteristics of the specific ARDS populations are described in Table 1.
Studies with mixed ARDS populations were divided into 3 groups according to the study's methodology: (1) prospective observational, (2) retrospective observational, and (3) randomized controlled trials (RCTs). In prospective and retrospective observational studies, we considered all enrolled subjects with ARDS; in RCTs, we screened only the mortality rates in the control (non-interventional) group to obtain a set of subjects comparable with the observational ones and to avoid possible influence of the intervention.
The reported mortality rates were also categorized by type: in-hospital, ICU, 28-d, 30-d, 60-d, 90-d, 365-d, etc. Some papers reported more than one type of mortality rate. For further analysis, only the most frequent mortality types were selected, and they were analyzed separately: in-hospital, ICU, 28/30-d mortality (note: 28-d mortality and 30-d mortality were considered together), and 60-d mortality. The rest of the mortality types that were reported rarely in the papers were omitted from further analysis.
To assess the aspect of time, we subdivided the main groups (prospective observational studies, retrospective observational studies, and RCTs) into 3 subgroups in relation to the time of the initiation and performance of the study: (1) before 1995, before the publication of the AECC criteria (subgroup A), (2) from 1995 to 2000, when the highlighted ARDS network study9 was published (subgroup B), and (3) from 2001 to the present (subgroup C). Overall mortality rates (weighted by number of subjects) and their trends over time were calculated for all subgroups. Moreover, we separately calculated the most recent cumulative mortality rates since 2010.
To resolve the bias stemming from different sample sizes among the studies, aggregated mortality rates were computed as a weighted average of the mortality rate reported in each study, with the weight being the number of subjects involved. To assess simultaneous influences of the independent variables on weighted mortality, a binomial family generalized linear model was fitted to the data. Ordered categories A, B, and C, which represent the time periods, were encoded using orthogonal polynomial coding, which allows trends to be captured in a linear or quadratic manner.
The terminology and classification of acute respiratory failure used in this review and during the search of databases respected the usual nomenclature in the respective period of time (eg, the AECC criteria).1,2 The authors are aware of the current recommendations to use only the term ARDS according to the recently published Berlin definition.3
The certified statistician (MB) used R 3.2.0 (R Foundation for Statistical Computing, Vienna, Austria) for statistical analysis and for creating the figures and tables presented in this study.10 The data set obtained underwent statistical assessment for quality. Missing values, outliers, and other suspicious values were noted and double-checked. The analysis of the data set dealt mainly with mortality rates (ie percentage rates of occurrence of some event, which is a random variable with binomial distribution). Therefore, to resolve the bias stemming from different sample sizes among the studies, aggregated mortality rates were obtained by computing arithmetic means of reported mortality rates weighted by numbers of subjects included in the studies. To assess simultaneous influences of multiple variables on weighted mortality, a binomial family generalized linear model was fitted to the data.
Results
Study Selection and Characteristics
The initial search identified 1,265 potentially eligible articles. A total of 1,096 papers were excluded for the following reasons: 843 were considered to be not relevant mainly because of their topic, focus, and methodological quality; 85 involved only pediatric subjects; 20 were preventive; 25 used an inappropriate definition; 6 were missing required data; 94 were reviews/meta-analyses/post hoc analyses; 5 were congress proceedings/reports; 16 were case reports or included fewer than the required number of subjects; and 2 were currently recruiting for trials. The bibliographies of the retrieved articles contained 8 articles that met the inclusion criteria and were additionally considered for this systematic review. In total, the resulting 177 papers reported 189 enrollment intervals (historical cohorts), which we treated as separate studies (see Fig. 1).
Mortality Rates in Relation to Study Methodology and Time Period
In general, mortality rates ranged from 13 to 80% in the specific group and from 11 to 87% in the mixed group. More studies screened mixed ARDS populations than specific ARDS populations (n = 139 vs n = 50, respectively). Many articles report multiple mortality rates. The list of specific ARDS populations is displayed in Table 1. The total number of subjects in all selected studies was 38,351.
Among the mixed ARDS population, the most frequently reported type of mortality rate was in-hospital (no. = 50), followed by ICU mortality (no. = 34), 28/30-d mortality (no. = 31), and 60-d mortality (no. = 16). Only studies with mixed population were selected for further analysis (no. = 106). The total number of subjects in the mixed group was 25,966. Figure 2 depicts the number of studies and the total number of subjects involved in studies categorized by mortality type and study methodology. Table 2 shows in detail the number of studies, number of subjects, and average mortality rates related to study methodology, mortality type, and time period.
The generalized linear model with the binomial family was used to assess simultaneous influences of various variables (number of centers, number of months, and study design) on mortality rates. A logarithm of odds ratios (logits) was predicted using the linear regression (see online supplement II). The mortality is significantly influenced by the study design. If the type of mortality is 28/30-d or 60-d, the logits (as surrogate for the mortality rate) decrease significantly by 0.44 or 0.45 (both P < .001, respectively). If the study is prospective (ie, prospective observational studies or RCTs), then the reported mortality increases by 0.2 logits (P < .001).
The orthogonal polynomial contrasts allowed us to evaluate the time period linearity and quadraticity, assessing the influence of time period on the reported mortality rates. The mortality rate has a linear decreasing tendency by logit 0.15 in relation to time period (P < .001). Quadratic influences of time period were not identified as statistically significant.
The mortality rates showed significant dependence on the number of centers. With an increasing number of centers included in the studies, the reported mortality rate decreases by 0.003 logits (P < .001). A dependence of mortality rates on the length of the studies (months) was not found (P = .57).
Table 3 shows the detailed mortality rates in defined time periods (A <1995; B 1995–2000; C >2000) grouped by prospectivity or retrospectivity. Lower mortality was reported in retrospective mode in almost all calculable subgroups (9 of 10). A significant difference was found for 28/30-d mortality rates in period A (P < .001), in-hospital mortality in periods B and C (P < .001, and P < .001, respectively), and ICU mortality rates in period B (P = .018) (4 of 10).
The differences between prospective modes of study subgroups (prospective observational studies and RCTs) are described in Table 4. We found higher mortality in prospective observational studies compared with RCTs in 7 of 9 calculable subgroups with significance for in-hospital morality in periods B and C (both P < .001, respectively), in 28/30-d mortality in the same periods (P < .001 and P = .046, respectively), and in 60-d mortality in period C (P < .001) (5 of 9).
Trends and Recent Mortality Rates
A detailed analysis of trends (see Fig. 3) shows that the mortality rate decline is statistically significant for in-hospital mortality in retrospective observational studies, ICU mortality in RCTs, 28/30-d mortality in prospective observational studies and RCTs, and 60-d mortality in RCTs. There was a mild to moderate decrease tendency for ICU mortality in prospective observational studies and for 60-d mortality in RCTs. There was no significant change in mortality for in-hospital RCTs, ICU mortality in retrospective observational studies, and 28/30-d mortality in retrospective observational studies. The significant increase in mortality was observed for in-hospital mortality in prospective observational studies. The data for calculation of the 60-d mortality trend in retrospective observational studies was missing due to a lack of studies reporting the particular type of mortality. Overall, 7 of the 12 subgroups showed a decrease (58%) and 3 showed stagnation (25%) of the mortality rate trend. Moreover, the mortality rate decrease was present in 6 prospective and one retrospective subgroup.
Table 5 summarizes the absolute values of the most recent mortality rates since 2010. Currently, the overall rates of reported in-hospital, ICU, 28/30-d, and 60-d mortality are 45, 38, 30, and 32%, respectively.
Discussion
Reporting of ARDS epidemiology is extremely heterogeneous. Its variability is based not only on differences between causal factors (type and intensity) and route of lung injury (pulmonary vs extrapulmonary ARDS) but also on genetic background, the presence of various complications (ie, secondary infections), and either acute or chronic comorbidities. The cause of death in patients with ARDS is usually not lung pathology per se but mostly concurrent extrapulmonary organ dysfunction. Early mortality is often related to the causal factor of the lung injury, and late mortality is more associated with complications (eg, sepsis and multiple-organ system failure).11,12
From the first clinically relevant description of acute lung injury/ARDS in 196713 to the early 1990s, the epidemiological data often yielded ambiguous and biased results due to the lack of generally accepted diagnostic criteria. The studies also included subjects with infantile respiratory distress syndrome, acute kidney injury, fluid overload, drug intoxication, disseminated intravascular coagulation, and burns.14 AECC criteria established in 1994 enhanced the diagnostic accuracy of the syndrome and facilitated research in the field. However, many articles published in this period reported mortality rates in a surprisingly wide range from 15 to 72%.6,15–25 Our findings confirmed the wide span of mortality rates during the time period (ie, from 13 to 80% in the specific group and from 11 to 87% in the mixed group). The third period (since 2001) is characterized by the implementation of lung-protective mechanical ventilation positively affecting the mortality of ARDS.9
The mentioned mortality differences could be explained by a type of the causal pathologic condition (eg, sepsis and hospital-acquired/ventilator-associated pneumonia are associated with higher mortality, whereas trauma-related injury or aspiration is associated with a lower mortality).26,27 The incidence of nosocomial infections varies among regions,28 with the type of pathogen (viruses vs bacteria, antibiotic-sensitive vs multidrug-resistant bacteria) also playing a significant role.29,30 The other explanations are as follows: population-based features (eg, composition of genetic factors, age, and sex); environmental features (eg, level of economic development, prevailing type of industry, local environment, cultural habits); health-care system differences (eg, health-care policy and area coverage, type of health insurance, and financial resources); research-related factors (eg, variability in the process of diagnosis; time factors [length of the enrollment, years vs months]; number of centers and ICUs involved; number of ICU beds; and type of data collection, processing, and storage). More reasons for the mortality-related differences could be (1) the data acquisition in various regions during specific time periods (ie, the changes in etiological patterns during year [incidence of viral infections in late winter, incidence of multiple traumas during warm seasons]; (2) the incomplete application of favorable interventions for patients with ARDS in routine practice; (3) the absence of discrimination between the terms acute lung injury and ARDS in diagnostic process or using the inappropriate term “acute respiratory failure”; (4) changes in the frequency of specific forms of ARDS (eg, due to outbreaks of various infectious diseases [severe acute respiratory syndrome or influenza A H1N1 2009]); (5) gradual increases in mean age and the number of serious comorbidities of subjects; and (6) differences in the general hospital management of subjects after discharge from ICU, which might contribute up to 3–15% of the overall mortality.31,32 These mentioned factors might also limit our analysis and findings. Moreover, we observed that the types of mortality rates reported by studies often vary (eg, overall, 28-d, 30-d, 60-d, 90-d, in-hospital, ICU, 8-week), which significantly contributed to the inconsistency among clinical studies. The 28/30-d mortality and 60-d mortality are qualitatively different from ICU and in-hospital mortality. We suggest that due to the extended hospital stay of patients with ARDS, the ICU and especially in-hospital mortality are close to or may even be longer than 28/30 d. The 60-d mortality is described in a limited number of studies mostly reporting just one type of mortality. One can see in Figure 3 that in the corresponding time periods, in almost all cases, overall 28/30-d mortality is lower than 60-d mortality, as could be expected. The same applies also for ICU and in-hospital mortality. The particular study methodology could have a significant influence on reported mortality values and might be a source of bias.
Given the known problems of static mortality data, an analysis of the kinetics (trends) of mortality rates over time33–38 seems more valuable. In the period before 1994, several studies reported steady mortality rates (51 ± 19%) from 1967 to 1994.39 However, Milberg et al40 reported different results, describing a steady state of mortality from 1983 to 1987 and a subsequent decline until 1993. In another single-center study, the cohorts of subjects in specific years (1982, 1990, 1994, and 1998) presented a decrease in overall mortality, from 68% in 1982 to 29% in 1996, with a plateau in the mid-1990s.41 In 2008, Zambon and Vincent15 analyzed 72 studies, including a total number of 11,426 subjects, and described a tendency for mortality to decrease from 1994 to 2006. They reported a wide range of mortality rates, from 15 to 72%, with an overall pooled mortality of 43%. They observed reductions in mortality in both interventional and non-interventional studies (trials with and without inclusion/exclusion criteria, respectively). The authors suggested that the improvement in survival may not have only been due to better respiratory management, but several other non-pulmonary factors could also have played a role (eg, better timing and rationalization of therapeutic interventions [antibiotic administration, fluid therapy, blood transfusions], glucose control, hygienic measures, improved infection source control and overall sepsis management).15 The second study challenged the decrease in mortality trends. It analyzed data from 89 studies conducted between 1984 and 2006, including 53 observational and 36 randomized trials with a total of 18,900 subjects. The reported overall pooled mortality rate was 44%. The observational studies showed a decrease in mortality until 1993, but the overall mortality did not decrease from 1994 to 2006. The mortality value was lower in RCTs (36%) than it was in observational studies (44%).16 We subdivided the studies in greater detail according to time, study methodology, and reported type of mortality rate. In contrast to Phua et al,16 we found a significant decreasing tendency for in-hospital mortality in retrospective observational studies, ICU mortality in RCTs, and 28/30-d mortality in prospective observational studies and RCTs over all 3 followed time periods. Furthermore, Phua et al16 found a significant difference between observational and randomized trials (including prospective studies and both interventional and non-interventional arms of RCTs). We did not conclusively confirm their results, because only 5 of 9 of calculable subgroups (56%) were significantly different. The reason for the difference might be that we analyzed a different number of studies (139 vs 89, respectively) including more subjects (25,996 vs 18,900). Moreover, we subdivided the studies according to mortality type and deliberately did not include interventional arms of RCTs due to possible influence of particular intervention on mortality results.
Moreover, we performed a comparison based on prospectivity versus retrospectivity. We found that prospective studies (prospective observational studies and non-interventional arms of RCTs considered together) report significantly higher mortality rates only in 4 of 10 calculable subgroups (40%) compared with retrospective observational studies. Thus, we did not confirm the presumed disadvantage of retrospective data collection: the inability to enroll all eligible patients and missing opportunity to collect all required data.
In the regression model, we found that the number of centers negatively correlated with reported mortality. The reason might be that the studies including more centers enrolled fewer patients per center with more focused management considered as a modified Hawthorne effect within medical staff. The length of the time interval of subjects' enrollment has no significant influence on reported mortality. We suggest that the high variability of the length of studies prevents the demonstration of significance. The linear tendency of mortality to decrease over time might be considered as the most important finding.
Finally, to present the most recent mortality data, we additionally calculated the overall rates for in-hospital, ICU and 28/30-d and 60-d mortality (45, 38, 30, and 32%, respectively) since 2010. We suggest these values as possible benchmarks for further clinical trials reporting the particular types of mortality.
Due to our experience with adult patients only, we omitted results from pediatric studies, and age <18 y of the subjects was one of the exclusion criteria. However, we are aware of the fact that children of adolescent age are somatically almost identical to adults and could have been included. It might be considered as a limitation of this review. However, we were afraid of problems that might have been caused by including the data of school-age and younger children with different physiology and ARDS risk factors compared with adults.
Conclusions
The reported mortality rates are highly heterogeneous and often biased, and their interpretation has many limitations. We advise against the unqualified adoption of reported mortality data from different studies performed in different time periods. Comparing these data with routine practice could be misleading, and using them in designing new clinical trials could diminish the value of the research results. We suggest carefully taking into account the study methodology (observational vs randomized and retrospective vs prospective, specific vs mixed ARDS population), with preference given to prospective modes of data collection and critical consideration of the age and length of the study. The mortality trends suggest a decreasing tendency in prospectively designed studies over time; however, mortality remains high.
Acknowledgments
We thank Dr E David G McIntosh (Imperial College, London, United Kingdom) for the final correction of the manuscript.
Footnotes
- Correspondence: Jan Máca MD, Department of Anesthesiology and Intensive Care Medicine, University Hospital, Ostrava, 17. listopadu 1790, 708 52 Ostrava-Poruba, Czech Republic. E-mail: jan.maca{at}fno.cz.
This study was supported by Grants SVV 260 026 and MO 1012, Charles University Prague. The authors have disclosed no conflicts of interest.
Supplementary material related to this paper is available at http://www.rcjournal.com.
- Copyright © 2017 by Daedalus Enterprises