Lung fluid biomarkers for acute respiratory distress syndrome: a systematic review and meta-analysis

Background With the development of new techniques to easily obtain lower respiratory tract specimens, bronchoalveolar lavage fluid and other lung fluids are gaining importance in pulmonary disease diagnosis. We aimed to review and summarize lung fluid biomarkers associated with acute respiratory distress syndrome diagnosis and mortality. Methods After searching PubMed, Embase, Web of Science, and the Cochrane Library for articles published prior to January 11, 2018, we performed a meta-analysis on biomarkers for acute respiratory distress syndrome diagnosis in at-risk patients and those related to disease mortality. From the included studies, we then extracted the mean and standard deviation of the biomarker concentrations measured in the lung fluid, acute respiratory distress syndrome etiologies, sample size, demographic variables, diagnostic criteria, mortality, and protocol for obtaining the lung fluid. The effect size was measured by the ratio of means, which was then synthesized by the inverse-variance method using its natural logarithm form and transformed to obtain a pooled ratio and 95% confidence interval. Results In total, 1156 articles were identified, and 49 studies were included. Increases in total phospholipases A2 activity, total protein, albumin, plasminogen activator inhibitor-1, soluble receptor for advanced glycation end products, and platelet activating factor-acetyl choline were most strongly associated with acute respiratory distress syndrome diagnosis. As for biomarkers associated with acute respiratory distress syndrome mortality, interleukin-1β, interleukin-6, interleukin-8, Kerbs von Lungren-6, and plasminogen activator inhibitor-1 were significantly increased in the lung fluid of patients who died. Decreased levels of Club cell protein and matrix metalloproteinases-9 were associated with increased odds for acute respiratory distress syndrome diagnosis, whereas decreased levels of Club cell protein and interleukin-2 were associated with increased odds for acute respiratory distress syndrome mortality. Conclusions This meta-analysis provides a ranking system for lung fluid biomarkers, according to their association with diagnosis or mortality of acute respiratory distress syndrome. The performance of biomarkers among studies shown in this article may help to improve acute respiratory distress syndrome diagnosis and outcome prediction. Electronic supplementary material The online version of this article (10.1186/s13054-019-2336-6) contains supplementary material, which is available to authorized users.


Background
Acute respiratory distress syndrome (ARDS) is a clinical syndrome comprising a rapid onset of respiratory failure in patients with risk factors, such as refractory arterial hypoxemia with low reaction to supplemental oxygen and the presence of bilateral infiltrates on radiographic imaging [1]. To date, the diagnosis of ARDS and acute lung injury (ALI) is mostly based on clinical characterization. Frequently-used criteria are the American European Consensus Conference (AECC) criteria [2] and the Berlin definition [3].
As the accuracy of a diagnosis of ARDS based only on the clinical syndrome has been questioned, countless studies have focused on the identification of biomarkers for ARDS. Terpstra et al. conducted a meta-analysis in 2014 focused on plasma biomarkers for ARDS in humans and provided a ranking system for distinguishing the disease from at-risk patients and determining the prognosis [4]. They reviewed multiple plasma biomarkers for ARDS, ranked by pooled odds ratio (OR). However, they only summarized biomarkers for ARDS in plasma; biomarkers in other fluids, such as bronchoalveolar lavage fluid (BALF), were not evaluated.
BALF and other lung fluids, such as pulmonary edema fluid (PEF), epithelial lining fluid (ELF), and lung aspirational fluid (LAF), are definitive in respiratory disease diagnosis. Since BALF provides a sample closest to the site of the disease process, it reflects the local lung environment directly. In 2017, García-Laorden et al. reported that biomarkers representing epithelial apoptosis, such as Fas and FasL, as well as biomarkers reflecting extracellular matrix injury, such as procollagen peptide III (PCP III) and procollagen peptide I (PCP I), were elevated in ARDS BALF samples [5]. The aim of the present study was to compare biomarker levels in lung fluid samples among patients with ARDS and the at-risk controls, as well as those of non-survivors versus survivors of ARDS.

Data source and study selection
We manually searched PubMed, Embase, Wed of Science, and the Cochrane Library for studies on biomarkers for ARDS in lung fluid samples published prior to January 11, 2018. Details of the search strategy are listed in Additional file 1. We also searched the references of included studies. Two researchers screened and evaluated the eligibility of all studies independently, and a third reviewer intervened whenever there was a disagreement. The inclusion criteria were (1) original research report of adult with or at-risk of ARDS, (2) report of exact values of biomarker concentration in lung fluid related to a clinical outcome (diagnosis of ARDS in at-risk patients and/or mortality of ARDS), (3) description of demographic variables, and (4) written in English. The exclusion criteria were (1) written in languages other than English, (2) not related to ARDS/ALI, (3) not an original research, (4) in vivo/in vitro studies, (5) pediatric studies, (6) biomarker not measured in lung fluid, (7) biomarker used for treatment monitoring, and (8) only one article available for a specific biomarker for no mergeable effect size and low reliability.

Data extraction and quality assessment
We built Excel spreadsheets (Microsoft Corp., Redmond, WA) to extract data from the included studies, and the two researchers finished data extraction independently. The ARDS etiology and the mean or median level and standard deviation (SD) of the biomarker in the lung fluid were obtained. When a biomarker was measured sequentially, only the day 1 measurement was extracted. We extracted lung fluid biomarker levels from different subgroups as follows: patients with ARDS versus critically ill non-ARDS controls and survivors versus non-survivors in patients with ARDS. The mean value of a biomarker's concentration was equal to the median level in this study, and standard error (SE) was converted to SD using an Excel formula. In addition, demographic variables (age, sex, and number of participants for each subgroup), diagnostic criteria for ARDS, ARDS mortality, the moment the lung fluid sample was retrieved, the sample type (BALF/other than BALF), sample retrieval location, and volume of BALF irrigation solution used were recorded. The recovery rate of BALF was also recorded, if provided.
All studies were assessed for quality according to the Quality Assessment of Diagnostic Accuracy Studies Score-2 (QUADAS-2), and the content was tailored according to the guideline of QUADAS-2 [6]. Details of the tailored QUADAS-2 are listed in Additional file 2. Risk of bias and an applicability concerns graph/summary was conducted using Review Manager version 5.3 (Cochrane Collaboration, Oxford, UK).

Data synthesis and data analysis
Meta-analysis was performed with Stata 13.1 (StataCorp LLC, College Station, TX). The ratio of means (RoM) was employed to assess the effect size [7][8][9]. RoM is the mean value of a biomarker in the ARDS group divided by the mean value in the at-risk group (mean ARDS /mean at-risk ) or the mean value of a biomarker in the non-survivors group divided by the mean value in the survivors group (mean non-survivor /mean survivor ). RoM of each study was log transformed and pooled using the inverse-variance method to gain a pooled, transformed RoM, which was then back-transformed to determine the pooled RoM and 95% confidence interval, using the fixed effect model of the Stata software. The significance level for this meta-analysis model was set at p < .05. Forest plots were provided for biomarkers of which four or more studies were included in this meta-analysis. Biomarkers were ranked according to pooled RoM and statistical significance. We used the Q statistic to test the existence of heterogeneity; a p value of less than 0.10 was considered significant for heterogeneity. I 2 was employed to assess the proportion of total variability due to heterogeneity. An I 2 value of approximately 25% was regarded as low heterogeneity, 50% as medium, and 75% as high heterogeneity. Publication bias was assessed with Egger's regression test [10], where a p value of less than 0.10 was considered significant for publication bias. Duval and Tweedie's trim and fill was then conducted [11].
For the biomarkers with a significant RoM and existence of heterogeneity, we performed a subgroup meta-analysis on study type (case-control study versus another study type) or sample type (BALF versus other lung fluid), when three or more studies were included.

Literature search
The total literature search yielded 1156 articles from the databases as follows: PubMed, 340 articles; Web of Science, 522 articles; Embase, 279 articles; Cochrane Library, 12 articles; and 3 articles from the reference lists of included studies. By reviewing the titles and abstracts, studies were mainly excluded due to the following: in vitro/animal studies (n = 434), duplication (n = 423), not an original research (reviews, editorials, or case reports, n = 90), and biomarkers not related with occurrence or mortality of ARDS (n = 30). After the initial screening, 95 articles remained for full-text review. Of these, 25 articles only reported on a specific biomarker, 16 articles contained insufficient data, and 4 articles had no full-text copy available, despite attempts to contact the authors. The remaining 49 articles were used for the meta-analysis   (Fig. 1).

Study characteristics and quality assessment
Demographic variables of the included studies are summarized in (Table 1). A total of 49 articles involving 2189 patients were included in this meta-analysis. ARDS/ALI was diagnosed according to the AECC criteria in 71% of the studies. Other criteria, such as edema fluid/plasma protein ratio [61], lung injury score [62], Fowler criteria [63], and clinical criteria, were used along with the AECC criteria. The mean age ranged from 37 to 70 years, mortality rate ranged from 15 to 77%, and lung fluid was collected between 30 min of intubation and 72 h of ARDS diagnosis. As for sample retrieval location, the right middle lobe and lingular lobe were the most common. Other locations were based on abnormal areas identified on chest radiographs, and a blind sampling of BALF was performed in two studies. In regard to sample type, 69% of the studies measured biomarkers in BALF with a certain volume of irrigation solution.  Ten articles used pulmonary edema fluid, and four articles measured a biomarker in epithelial lining fluid. Only three articles provided the recovery rate of the irrigation solution; therefore, we could only assume a stable recovery rate between subgroups for this study.
The quality assessment is displayed in Additional file 4, including the risk of bias and applicability of studies to the review question.

Biomarkers associated with ARDS diagnosis
We performed a meta-analysis on 22 biomarkers in lung fluid associated with the diagnosis of ARDS in the at-risk population ( We performed an influence analysis to examine the sensitivity of the results. Influence analysis showed that the heterogeneity was possibly caused by the limited Numbers within the square brackets were reference numbers RoM ratio of means, CI confident interval number of studies. By removing the studies with extreme RoM, we observed a robust effect on the biomarkers. The outcome of the influence analysis is displayed in Additional file 5. Subgroup analysis was performed for biomarkers IL-8, total protein, albumin, sRAGE, PAF-AcH, and IL-6. Since most of the included studies were case-control studies, we excluded studies with other design types. Heterogeneity for total protein, albumin, and IL-6 was partly explained. However, heterogeneity for IL-8 was not clarified. Only one article was case-control study; therefore, the source of heterogeneity was not determined for PAF-AcH because of a limited number of study. We also excluded studies with ARDS that was not diagnosed using AECC criteria, which significantly reduced the heterogeneity for IL-6, but not for albumin. In addition, we assumed sample type may be a variable between studies because biomarker measurement in BALF was influenced by recovery rate and dilution. IL-8 remained significantly increased when BALF studies were excluded. As only one article measured biomarkers in lung fluid other than BALF, it was not evaluated. Results of the subgroup analysis are presented in Additional file 6.

Biomarkers associated with ARDS mortality
We performed a meta-analysis on 11 biomarkers in lung fluid associated with ARDS mortality (Table 3) were associated with a high mortality rate. Heterogeneity was displayed for many of the biomarkers, and influence analysis indicated that the heterogeneities were not likely caused by extreme RoM values. Due to the small number of studies, subgroup analysis based on design type or sample type could not be performed. Subgroup analysis for tumor necrosis factor-α (TNF-α) was performed when excluding patients with ARDS not diagnosed using AECC criteria. Results of the subgroup analysis are presented in Additional file 7.

Publication bias
Among the biomarkers associated with ARDS, Egger's regression test demonstrated a p value of less than 0.10 for IL-6; furthermore, when we adjusted for possible publication bias by Duval and Tweedie's trim and fill, the RoM remained significant for IL-6. Among the biomarkers associated with mortality, no publication bias was noted. The results of the publication bias analysis are presented in Additional file 8.

Discussion
In this systematic review and meta-analysis, we summarized the biomarkers related to ARDS diagnosis in the at-risk population and those related to ARDS mortality.
By searching several databases and screening for related articles, 49 studies involving 2189 patients were identified.
We discovered that total protein, albumin, PAI-1, PAF-ACH, sTNFα-RII, HGF, IL-8, PCP I, PCP III, soluble receptor for advanced glycation end products (SRAGE), and IL-6 were significantly increased in the lung fluid of patients with ARDS. Although total PLA2 activity, soluble intercellular adhesion molecule-1 (SICAM-1), IL-2, CC16, and MMP-9 were also significantly different between patients with ARDS and at-risk patients, few studies were included on each of these biomarkers, so the results are unreliable.
Pioneering work by Terpstra et al. reported on plasma biomarkers for ARDS diagnosis and prognosis. They reported that KL-6, lactate dehydrogenase (LDH), SRAGE, von Willebrand factor (vWF), and IL-8 displayed the highest effect size for ARDS diagnosis, and interleukin-4 (IL-4), IL-2, Angiopoietin-2 (Ang-2), and KL-6 had the highest effect size for ARDS prognosis (assessed by pooled odds ratio). These biomarkers represent pathophysiological processes, which led to the hypothesis that ARDS diagnosis is correlated with tissue damage, whereas ARDS mortality is correlated with systemic inflammation (4). In our meta-analysis, biomarkers for ARDS diagnosis were related to inflammation (IL-8 and IL-6), endothelial injury (SICAM), epithelial injury (SRAGE and HGF), lung fibroproliferation (PCPI and PCPIII), and coagulopathy (PAF-ACH). With regard to ARDS mortality, biomarkers related to inflammation (IL-8, IL-6, and IL-1β), epithelial injury (KL-6), and lung fibroproliferation (PCPIII) were elevated in the lung fluid of patients with ARDS presenting the worst outcomes. Therefore, we assume that, aside from tissue damage and systemic inflammation, lung fibroproliferation is vital in both ARDS diagnosis and prognosis. Several studies on lung biopsy of patients with ARDS showed a strong relationship between fibrosis activity and ARDS mortality [64,65]. To our knowledge, this is the first meta-analysis of biomarkers in lung fluid for ARDS. Since lower respiratory tract specimens are now easily obtained, BALF and other lung fluids have become frequently used clinical samples for pulmonary disease diagnosis, secondary to plasma/serum. Since the most intensive physiological processes in ARDS occur in the lung, theoretically, lung fluid can reflect the pathophysiological process differently from other body fluids.
In this analysis, we applied RoM to assess the effect size and to attempt to eliminate the bias caused by dilution of different kinds of lung fluid. For example, edema fluid was completely undiluted, whereas BALF might be quite diluted. This methodology has been proved to be robust and is widely used [7-9, 66, 67].
This systematic review and meta-analysis may prompt further research on ARDS diagnosis and prognosis in many different ways. First, it demonstrates the research priorities and indicates that research on ARDS biomarkers in lung fluid and other compartments is needed. Second, it establishes ARDS biomarkers and their performance in lung fluid as an innovative research field. Third, it identifies numerous novel translational approaches for biomarker measurement in different compartments, such as chromatography for metabolomics separation and mass spectrometry or nuclear magnetic resonance spectroscopy for biomarker detection [68]. Although some non-quantitative methods were not included in this meta-analysis, they are still worth exploring.
There were limitations in this meta-analysis as well. First, although we performed a subgroup analysis of the biomarkers related to ARDS diagnosis and mortality, heterogeneity was not explainable for every biomarker. We assume that this could be related to the different etiologies of ARDS, variation in BALF procedures between studies, multiple control types used in the studies, wide range of intervals between study inclusion and biomarker measurement, and different treatments for ARDS. We were not able to conduct further analysis due to the limited information.
Second, a limited number of studies were included for each biomarker, which impedes the reproducibility of the results. The number of studies for each biomarker should be taken into consideration while assessing the performance in the ranking system. Third, only the biomarkers addressed by two or more studies were included. Because of this, promising biomarkers evaluated in a single study were not considered, which may limit the view on lung fluid biomarker research as a whole.
Finally, the use of lung fluid as a study object itself had some limitations. Due to the lack of information on specific BALF procedures, we could only assume the recovery rates between subgroups in one study were within an acceptable range, which may have caused some of the heterogeneity between studies.