Abstract
Background The fungal cell-wall constituent (1,3)-β-d-glucan (BDG) is a pathogen-associated molecular pattern (PAMP) that can stimulate innate immunity. We hypothesized that BDG from colonizing fungi in critically-ill patients may translocate into the systemic circulation and thus be associated with host inflammatory responses and outcomes.
Methods We enrolled 453 mechanically-ventilated patients with acute respiratory failure with no evidence of invasive fungal infection (IFI). From serial plasma samples, we measured BDG, innate immunity and epithelial permeability biomarkers. From lower respiratory tract and stool samples we quantified bacterial and fungal DNA load using culture-independent techniques.
Results A positive BDG test (>60pg/ml) at baseline was detected in 19% of patients. BDG levels were significantly associated with markers of innate immunity (interleukin-6, tumor necrosis factor receptor-1 and procalcitonin), epithelial barrier disruption (receptor for advanced glycation end-products and fatty-acid binding protein-2, for lung and gut respectively) and with higher probability of classification in an adverse prognosis hyperinflammatory subphenotype (all p<0.05). No differences in fungal or bacterial DNA load were found by BDG test positivity. Positive BDG testing was associated with higher incidence of acute kidney injury, fewer ventilator free days and worse 30-day survival (adjusted p<0.05). Patients with positive BDG test on follow-up sampling (>3 days from intubation) had higher mortality than patients with persistently negative test on follow-up (p<0.05).
Conclusions This is the first study to demonstrate the prognostic role of BDG in critically ill patients with no evidence of IFI. Translocation of BDG into systemic circulation may contribute to inflammation and clinical outcomes.
Funding support National Institutes of Health [K23 HL139987 (GDK); U01 HL098962 (AM); P01 HL114453 (BJM); R01 HL097376 (BJM); K24 HL123342 (AM); U01 HL137159 (DVM, PVB); R01 LM012087 (DVM, PVB); R01 HL142084 (JSL); R01 HL136143 (JSL); F32 HL137258 (JWE); F32 HL142172 (WB); K08 HS025455 (IJB); K23 GM122069 (FS)].
Introduction
Dysregulation of innate immunity is a central pathogenetic feature of the heterogeneous syndromes of sepsis and the acute respiratory distress syndrome (ARDS) (1, 2). In several independent patient populations, patient subgroups with elevated plasma levels of inflammatory cytokines and biomarkers of tissue injury (classified as a hyperinflammatory subphenotype) exhibit worse organ dysfunction and clinical outcomes compared to their hypoinflammatory counterparts (3–5). Such biologic subphenotyping of critically-ill patients offers improved risk prediction over stratification with clinical variables alone, and presents new opportunities for targeted pharmacologic interventions (6). However, the molecular mechanisms of aberrant innate immunity stimulation and resultant tissue injury remain largely undefined.
Inter-patient variability of endogenous microbiota may explain the observed differences in host-responses. Many of the predictive biomarkers for the detrimental hyperinflammatory subphenotype (e.g. interleukins [IL]-6 and −8, soluble tumor necrosis factor receptor-1 [TNFR1], procalcitonin and receptor for advanced glycation end-products [RAGE]) (3, 4) reflect canonical pathways of immune cell stimulation from pathogen and damage associated molecular patterns (PAMPs and DAMPs, respectively) (7, 8). PAMPs include microbial nucleic acids as well as cell wall components, such as lipopolysaccharide of gram-negative bacteria and 1-3-beta-D-glucan (BDG) from fungi (8, 9). Recent culture-independent studies with bacterial DNA sequencing have uncovered important roles of intestinal and respiratory bacterial communities in the evolution and outcome of critical illness (10–13). However, very few studies have performed in depth examinations of the fungal communities (mycobiome) in critical illness (14) and little is known about how fungal PAMPs may interact with the innate immunity system in the critically-ill host.
Epidemiologic evidence suggests that host-fungal interactions at different mucosal surfaces are clinically important (15). In mechanically-ventilated immunocompromised patients, airway colonization by Candida species portends worse outcome (16), and in patients with suspected ventilator-associated pneumonia, Candida species in respiratory tract secretions have been associated with increased mortality and levels of IL-6 and procalcitonin (17, 18). C.albicans overgrowth in the human gut occurs in the setting of antibiotic-treated sepsis, raising the risk for secondary candidemia (19, 20). The respiratory and intestinal tracts in critical illness contain high levels of the fungal PAMP BDG. BDG is recognized by mucosal innate immune cells through the dectin-1 receptor and stimulates local inflammatory responses (21). It can also exert systemic effects on host inflammation via translocation into the bloodstream, in the setting of disruption of mucosal epithelial integrity from various insults, such as hypoxemia, ischemia or direct pathogen invasion (22). Preliminary observations have suggested that higher levels of circulating BDG are associated with worse clinical outcomes in patients in the intensive care unit (ICU) (23), but evidence has been limited. BDG levels are routinely checked in the diagnostic work-up of critically-ill patients with clinical suspicion of invasive fungal infection (IFI) (24), but whether BDG could be used as a clinical biomarker of fungal PAMP translocation and its contribution to the development and perpetuation of hyperinflammatory host-responses in critical illness is unknown.
In this study, we demonstrate that approximately 20% of critically-ill patients without IFI have elevated plasma BDG levels (as defined by a positive BDG test). Patients with a positive plasma BDG test had higher levels of circulating biomarkers of innate immunity and tissue injury, were more likely to be classified as the adversely prognostic hyperinflammatory subphenotype, and had worse clinical outcomes independent of other prognostic factors. Exploratory analyses for the potential source of BDG e.g. the intestinal or respiratory tract, suggested the potential presence of heterogeneous mechanisms underlying fungal PAMP translocation in critical illness.
Results
Study population
We analyzed data from 453 mechanically-ventilated patients with acute respiratory failure, who were consecutively enrolled from ICUs at the University of Pittsburgh Medical Center (UPMC) (3, 25). We excluded 11 patients with a clinical diagnosis of IFI (24) because we wanted to examine BDG as a fungal PAMP in a broad ICU population and not in the context of IFI as a diagnostic test. However, inclusion of patients with IFI in sensitivity analyses did not impact the observed associations and outcomes. Characteristics of included patients are shown in Table 1.
Baseline characteristics of enrolled subjects, stratified by positive and negative BDG testing. P-values from nonparametric tests comparing patients with a negative (≤60 pg/ml) and positive (>60 pg/ml) BDG result are shown in bold when significant (p<0.05).
Distribution of plasma BDG levels in critically-ill patients
At baseline, 85 patients (19%) had positive BDG defined as greater or equal to the conventional threshold of 60 pg/ml (Figure 1) (24). Patients with positive BDG levels (median 95, interquartile range [IQR] 76-227 pg/ml) had similar distribution of baseline clinical characteristics as patients with negative BDG (median 26, IQR 15-49 pg/ml), with the exception of higher sequential organ failure assessment (SOFA) scores in the positive group (p=0.01, Table 1). Similar results were obtained when we used a more stringent threshold of BDG positivity (80 pg/ml, Table S1). Although beta-lactam antibiotics have been previously implicated in false-positive BDG detection (26, 27), we did not find any differences in BDG levels between patients receiving beta-lactam antibiotics, or not, at the time of sampling. A small proportion of patients (2.9%) receiving topical antifungals for either oral thrush or skin candidiasis (clotrimazole troche or nystatin) had significantly higher BDG levels compared to those not receiving these medications (n = 13 [2.9%], p=0.03).
19% of patients had BDG levels >60pg/ml which is the conventional threshold for a positive BDG test in the diagnostic workup of invasive fungal infection.
Plasma BDG correlates with host-response biomarkers and the hyperinflammatory subphenotype
We systematically profiled host-responses by measuring plasma biomarkers in pathways of innate immunity (IL-6, IL-8, IL-10, TNFR1, fractalkine, suppression of tumorigenicity-2 [ST-2]), alveolar epithelial injury (RAGE), endothelial injury (angiopoeitin-2) and response to bacterial infection (procalcitonin and pentraxin-3) (3). BDG levels at baseline positively correlated with levels of all measured biomarkers (p<0.05), with seven of these correlations remaining significant after adjustments for multiple testing (Figure 2A). We synthesized these host-response profiles into a parsimonious regression model for classification into a hyperinflammatory (28%) vs. a hypoinflammatory subphenotype (72%) based on latent class analyses as previously described (3). Hyperinflammatory patients had significantly higher mean BDG levels than hypoinflammatory ones (39 [23-68] vs. 23 [14-39] pg/ml, p<10−6, Figure 2B).
A. Correlogram of 10 host-response biomarkers and BDG levels. Effect sizes of positive correlations are shown with variations of red, and statistical significance by size of circles. Comparisons were adjusted for multiple testing with the Benjamini-Hochberg method. Non-significant correlations are marked with X. B. Patients assigned to the hyperinflammatory subphenotype had significantly higher levels of BDG (p<0.0001). C. Probabilistic Graphical Model analysis demonstrating first and second neighbors of the BDG variable, when considered in conjunction with 32 clinical and biomarker variables. Two first neighbors were identified: TNFR1 and ST-2. Abbreviations: IL-6: Interleukin-6; IL-10: Interleukin-10; TNFR-1: Tumor necrosis factor receptor-1; Ang-2: Angiopoietin-2; ST-2: Suppression of tumorigenicity-2; RAGE: Receptor for advanced glycation end products.
Network analyses reveal association of BDG with key elements of the hyperinflammatory response
To comprehensively examine for causal associations between baseline clinical variables, biomarkers and BDG, we utilized the Probabilistic Graphical Models (PGM) framework (3, 28, 29). From a total of 32 variables, only two were directly linked with BDG (i.e. first neighbors in PGM): TNFR1, a predictive biomarker for the hyperinflammatory subphenotype, and ST-2, the soluble receptor of IL-33 that functions as an alarmin released by respiratory epithelia during acute lung injury (30). Both TNFR1 and ST-2 were further linked to additional components of the hyperinflammatory response (Figure 2C).
Patients with higher BDG levels had worse clinical outcomes
BDG levels analyzed as a continuous variable were associated in regression models with higher incidence of acute kidney injury (AKI), fewer ventilator-free days (VFDs) and higher 30-day mortality (all p<0.05, Table 2). These associations with clinical endpoints remained significant after adjustments for other important predictors of outcome in critical illness (age, SOFA scores and hyperinflammatory subphenotype) as well as potential confounders of BDG measurement (beta-lactam antibiotics that have been associated with false positive BDG and batch of BDG measurement). Patients with positive BDG testing had worse 30-day survival (adjusted hazard ratio: 1.37, 95% confidence interval 1.13-1.67, p<0.001) compared to patients with negative BDG (Figure 3).
Associations between BDG levels and clinical outcomes. Outcomes of BDG-positive and negative patients are shown in the first two columns. P-values for the statistical significance of the association between BDG levels and each outcome (Logistic regression models for acute kidney injury and 30-day mortality; zero-inflated negative binomial regression for ventilator-free days) are shown for unadjusted and adjusted analyses (adjustments for age, sequential organ failure assessment score, hyperinflammatory subphenotype assignment, beta-lactam antibiotics and batch of BDG measurement). Statistically-significant associations are highlighted in bold.
Kaplan-Meier curve for 30-day survival. Hazard ratio adjusted for age, sequential organ failure assessment score, hyperinflammatory subphenotype assignment, beta-lactam antibiotics and batch of BDG measurement.
BDG is associated with markers of epithelial barrier disruption in the respiratory and intestinal tract
Given the clinically significant prognostic associations of BDG, we sought to understand potential sources of plasma BDG and whether it is derived from intestinal or respiratory tract fungi. To that end, we examined associations between plasma BDG with biomarkers of epithelial barrier disruption and differences in microbial load (bacterial and fungal) in each body compartment studied by culture-independent techniques.
For the respiratory tract, plasma BDG levels were significantly associated with plasma levels of RAGE (p<0.005), a biomarker of alveolar epithelial injury. This association remained significant after adjustment for clinical variables potentially related to fungal colonization in the airways (systemic antibiotics and steroid administration on day of sampling). We directly quantified bacterial and fungal load in lower respiratory tract communities (endotracheal aspirate samples) by quantitative polymerase chain reaction (qPCR) of the bacterial 16S rRNA gene and by absolute number of reads (sequences) of the Internal Transcribed Spacer (ITS) region sequencing of the fungal genome, respectively. BDG was not significantly associated with bacterial (16S qPCR) or fungal load (number of ITS reads) in lower respiratory tract samples. However, BDG remained significantly associated with RAGE (p<0.005) after adjustment for bacterial and fungal load, indicating that epithelial permeability may be a more important contributor to BDG translocation rather than the absolute amount of fungi in the airways.
To examine the potential origin of BDG from the intestinal tract, we measured plasma levels of fatty acid binding protein-2 (FABP-2), a validated biomarker of intestinal barrier integrity (31), in a random subset of 220 patients. BDG levels were significantly associated with FABP-2 (p<0.05, Figure 3), an association that remained significant after adjustment for variables potentially related to intestinal tissue ischemia and fungal overgrowth (administration of vasopressors and systemic antibiotics at time of sampling, respectively). We did not find any significant association between plasma BDG and measurements of fungal and bacterial load in stool samples in a small subset with available data (n=16). Thus, integrative analyses of biomarkers of epithelial disruption and microbial load suggested that the possibility that BDG may translocate to the systemic circulation from both the respiratory and gastrointestinal tract.
Longitudinal evolution of BDG levels and association with outcomes
For 156 patients (34% of the cohort), follow-up plasma samples were obtained in the time window from day 3 to 5 from baseline while the patients remained in the ICU. Patients with a negative baseline test who developed a positive test on follow-up (N = 10, 6%) and those with an initial positive test who had a subsequent negative test (N = 17, 11%) had worse 30-day mortality compared to patients who remained negative on follow-up (p<0.05, Figure 5).
P-values of linear regression models are shown (A: effect of RAGE adjusted for systemic antibiotics and systemic steroids expressed as prednisone equivalent dose at the time of sampling; B: effect of FABP2 adjusted for receipt of vasopressors and systemic antibiotics).
Patients with a negative baseline test who developed a positive test on follow-up (N = 10, 6%) and those with an initial positive test who had a subsequent negative test (N = 17, 11%) had worse 30-day mortality compared to patients who remained negative on follow-up (60% vs 23% and 47% vs 23% respectively, p<0.05).
Discussion
We demonstrated that BDG is detectable in the plasma in 19% of mechanically ventilated patients from a general ICU population. BDG levels correlated with host-response biomarkers of injury and inflammation and independently predicted adverse clinical outcomes. Over time, patients with incident BDG positivity had higher mortality than patients with persistently negative BDG testing. Significant associations with markers of intestinal and respiratory tract epithelial permeability suggested the possibility of fungal PAMP translocation from either vascular site, independent of fungal or bacterial loads.
BDG represents a highly immunostimulatory molecule that is likely to mediate the effects of colonizing fungi in critically ill patients. As a constituent of the cell-wall of many potentially pathogenic fungi, BDG is recognized by pattern recognition receptors such as complement receptor 3 and dectin-1 present on the surface of macrophages, dendritic cells, neutrophils and pulmonary epithelial cells (21, 32–34). Co-stimulation of dectin-1 with Toll-like receptors augments cytokine production via NF-κB-dependent (nuclear factor kappa-light-chain-enhancer of activated B cells) signaling (32, 34). Such innate immune responses, albeit necessary for the recognition and clearance of invading fungi, may also lead to exacerbated mucosal inflammation and injury upon perpetual and aberrant stimulation.
The role of BDG in the pathogenesis of mucosal injury and local inflammation, however, may not be limited necessarily to that of a direct stimulant of the innate immune system. Fungal PAMPs in the form of heat inactivated Candida were shown to reduce the killing activity of bone-marrow derived macrophages while enhancing cytokine production (35), while intact Candida has been shown to reduce the phagocytic abilities of macrophages (36), potentially by a soluble molecule (37) among other poorly understood mechanisms. Such mechanisms can lead to impaired first-line defenses of mucosal immunity with concurrent precipitation of host injury and inflammation either by fungi or other pathogenic microbes at mucosal sites. Therefore, complex interactions between colonizing fungi, non-fungal members of the microbial community and the critically-ill host need to be considered.
Dysregulation of host-responses with high levels of circulating inflammatory cytokines and immune cell activation is a key pathogenetic feature in the evolution of critical illness (38). Unsupervised classification analyses in patients with sepsis, severe pneumonia or ARDS have discovered the presence of distinct subphenotypes, characterized by a hyperinflammatory vs. a hypoinflammatory cluster of host-responses(3–5). The markedly worse outcomes and differential treatment responses associated with the hyperinflammatory subphenotype offer a new framework for predictive enrichment and the delivery of targeted efficacious therapies (e.g. anti-inflammatory drugs to be given only to patients with the hyperinflammatory subphenotype who are likely to benefit) (6). However, the proximal, molecular determinants of ARDS and sepsis subphenotypes have not been defined, and, similarly, clinical stratifications by types of insult (e.g. pneumonia or extrapulmonary infection etc.) have not been informative for distinguishing host-responses (3, 39). Thus, we posit that at least some of the observed host-response heterogeneity may be accounted by host-microbial interactions. Beyond the well-recognized importance of bacterial pathogens, translocating fungal PAMPS and their interactions with innate immune cells may also significantly contribute to the evolution and outcome of critical illness.
Epidemiologic observations support the role of fungi in critical illness, even outside the context of conventionally defined IFI. Evidence of Candida colonization in the airways or fungal overgrowth in the gastrointestinal tract has been associated with worse clinical outcomes, including but not limited to the development of serious bloodstream infections due to systemic invasion of fungi from different mucosal surfaces (16, 17, 19, 20). Factors that have been associated with fungal colonization in critically-ill patients include widespread antibiotic use, which allow fungi to grow and overtake environmental niches (20), as well as primary or iatrogenic immunosuppression, such as with inhaled corticosteroid therapies leading to upper respiratory tract fungal colonization (40). Our study utilized culture-independent methods (BDG measurement and next-generation sequencing approaches) to examine for translocation of fungal elements into the blood stream, their interactions with innate immunity and prognostic importance in a broad ICU population with no clinical evidence of IFI.
We observed associations of circulating BDG with innate immune responses as well as biomarkers of epithelial permeability in the respiratory and intestinal tract. BDG may be stimulating innate immunity not only through mucosal interactions, but also in the bloodstream and in tissue locations within the reticuloendothelial system (35, 41). Therefore, circulating BDG molecules may not simply represent a marker of mucosal barrier injury that has found its way into the bloodstream, but also an active mediator of inflammatory organ dysfunction and illness severity in sepsis (9, 35). Our integrative analyses of permeability biomarkers and culture-independent assays of bacterial/fungal load did not reveal a specific primary source of circulating BDG from either the lungs or the gut. However, we identified significant independent associations of BDG with corresponding markers of epithelial permeability in both compartments. With regards to the lungs, soluble RAGE, an extensively validated biomarker of ARDS radiographic and clinical severity (25, 42, 43), was strongly associated with plasma BDG levels, potentially reflecting the degree of pulmonary epithelial/endothelial disruption and ability of BDG from airway colonizing fungi to enter the systemic circulation. Similarly for the gut, BDG levels were also associated with FABP2, a biomarker of intestinal barrier disruption (31). In animal models of intestinal leakage via direct physical or chemical injury, serum BDG testing showed low sensitivity but high specificity for a permeable gut (44). Direct measurement and study of the phenomenon of microbial translocation in human translational studies is challenging, and no gold-standard method exists for in vivo study (45). Thus, blood levels of BDG may offer a non-invasive and clinically available surrogate measure of gut/lung permeability in critically ill patients, with important prognostic implications.
A positive BDG test (defined as >60pg/ml) conferred an adverse prognosis in our inclusive cohort of subjects with acute respiratory failure. BDG-positive patients demonstrated higher severity of illness by SOFA score, accumulated more severe organ dysfunction (incident AKI and fewer VFDs) and had significantly worse 30-day survival. These associations with adverse outcomes remained significant after multivariate model adjustments for age, severity of illness and host-response subphenotypes, suggesting that BDG testing may capture prognostic information not offered by clinically available and emerging stratification tools based on biomarkers. Similar observations of poor prognosis associated with BDG had been reported in smaller cohorts of patients with suspected VAP as well as in patients with sepsis (23, 44). Our findings in a much larger and inclusive cohort suggest that BDG may offer additional prognostic information in the course of critical illness. Further validation of the prognostic utility of BDG in independent cohorts is needed.
Our study has several limitations. It is a single-center study, and the generalizability of our findings in critically ill populations beyond our tertiary care ICU requires external validation. Our study is also limited by the available sample size. Despite being the largest study of BDG measurement in acute respiratory failure, results from analyses for subgroups and specific biomarkers (e.g. FABP2) require cautious interpretation, as the effective sample size for certain analyses was smaller. Although we aimed to control for potential confounders of BDG measurement from available variables (batch of measurement, beta-lactam antibiotics or bacteremia), there may be other sources of BDG false positivity, such as hemodialysis or blood product preparation filter processing that were not possible to account for in our analyses (46, 47). Finally, the extent of microbiologic workup for diagnosis of bacterial infection or IFI was directed by the treating physicians and, as such, microbiologic testing was not uniform. Thus, it is possible that some of the BDG positive cases may represent clinically undiagnosed or preclinical IFI. Nonetheless, inclusion of patients with IFI in sensitivity analyses did not impact the observed associations.
In summary, we demonstrate that circulating BDG is an independent predictor of a hyperinflammatory host-response profile and adverse clinical outcomes in mechanically ventilated patients with acute respiratory failure. These associations highlight a novel, underappreciated mechanism of biological heterogeneity in critical illness that involves interactions between translocating fungal PAMPs and innate immunity, which is independent of fungal load at mucosal sites and points towards complex interactions between the host, fungi and other microbes. Our findings call for further mechanistic interrogation in animal studies and independent prospective clinical study to validate the biologic relevance and prognostic information of circulating BDG.
Methods
Clinical cohort and sample collection
From October 2011 – June 2019, we prospectively enrolled a convenience sample of consecutive, adult patients with acute respiratory failure, who were intubated and mechanically-ventilated in the Medical or Cardiac Intensive Care Units (ICU) at the University of Pittsburgh Medical Center. Exclusion criteria included inability to obtain informed consent, presence of tracheostomy, or mechanical ventilation for more than 72 hours prior to enrollment. The study was approved by the University of Pittsburgh Institutional Review Board (protocol PRO10110387), and written informed consent was provided by all participants or their surrogates.
From enrolled subjects, we collected blood samples for centrifugation and separation of plasma and serum. We also collected noninvasive biospecimens for study of lower respiratory tract microbiota with endotracheal aspirates (ETA) and we analyzed stool samples for study of the intestinal microbiota, as previously described (48, 49).
Laboratory analyses
We measured plasma BDG levels using the commercially available Fungitell® Limulus Amebocyte Lysate (LAL) assay (Associates of Cape Cod, Inc, East Falmouth, MA, USA) at the manufacturer’s facility (50). As a biomarker of intestinal permeability, we measured levels FABP2 with the Quantikine® Human FABP2/I-FABP Immunoassay (R&D Systems, Minneapolis, USA). For profiling host-responses, we constructed a custom Luminex multi-analyte panel (R&D Systems, Minneapolis) targeting biomarkers associated with ARDS outcomes (RAGE, IL-6, IL-8, IL-10, TNFR1, ST-2, fractalkine, angiopoietin-2, procalcitonin and pentraxin-3), as previously described (3). We extracted microbial DNA directly from samples and quantified bacterial load profiling by amplifying via quantitative polymerase chain reaction the V3-V4 hypervariable region of the bacterial 16S rRNA gene as previously described (48). For fungal community profiling, we sequenced amplified regions 1 and 2 of the ITS rRNA gene on the Illumina MiSeq platform (Supplement).
Clinical classifications
A consensus committee reviewed clinical and radiographic data and performed retrospective classifications of the etiology and severity of acute respiratory failure. Classifications were performed without knowledge of biomarker data. We retrospectively classified subjects as having ARDS per established criteria (49), at risk for ARDS due to the presence of direct (pneumonia or aspiration) or indirect (e.g. extra-pulmonary sepsis or acute pancreatitis) lung injury risk factors (51) but lacking ARDS diagnostic criteria, acute respiratory failure without risk factors for ARDS, or acute on chronic respiratory failure. We recorded clinical microbiologic results of respiratory and blood culture specimens obtained within 48hrs of research sample acquisition. We excluded patients with diagnosis of IFI based on microbiologic evidence and consistent clinical picture (24). We recorded antibiotic prescriptions, including beta-lactams associated with false-positive BDG tests, as well as topical antifungals for skin candidiasis or oral thrush. We followed patients prospectively for incidence of AKI, VFD at 30 days, and 30-day mortality.
Statistics
From baseline plasma BDG levels, we classified patients into those with a negative (≤60 pg/ml) or positive (>60 pg/ml) result. In sensitivity analyses, we also classified BDG levels using an 80pg/ml threshold(52). Biomarker values were log-transformed for use in regression models. From available biomarker and clinical variables, we classified patients into hyper-inflammatory vs. hypo-inflammatory subphenotype based on an internally derived and validated parsimonious logistic regression model, following application of latent class analyses (53). Baseline characteristics and outcomes between groups were compared using Wilcoxon rank-sum for non-normally distributed continuous variables and Student’s t-test for normally distributed continuous variables. Fisher’s exact was used for categorical variables. We used linear regression models to examine for associations between BDG levels and biomarkers of host responses. Multiple testing adjustments were implemented with the Benjamini-Hochberg method. We examined for differences in BDG levels between hyper-inflammatory and hypo-inflammatory subphenotypes using Wilcoxon rank sum tests. We examined for associations between BDG levels and clinical outcomes with logistic regression models for AKI and 30-day mortality, zero-inflated negative binomial models for VFD, and Cox-proportional hazard models for 30-day survival, adjusted for age, SOFA score, subphenotype classification, beta-lactam antibiotics and batch of BDG measurement. We examined for association between BDG levels and RAGE (marker of alveolar epithelial damage) using a linear regression adjusted for variables potentially associated with fungal overgrowth in the airways (systemic antibiotics and systemic steroids expressed as prednisone equivalent dose at the time of sampling). For the association between BDG and FABP2 (marker of intestinal permeability), we built a linear regression model adjusted for variables that may impact the integrity of the intestinal epithelium (receipt of vasopressors) and fungal overgrowth in the gut (systemic antibiotics). For bacterial and fungal microbiota analyses, we performed comparisons of microbial load (16S qPCR and N of fungal ITS reads) by BDG test results. Statistical significance was defined as p<0.05. All statistical analyses were performed with the R statistical software, version 3.6.0 (54).
Study approval
The University of Pittsburgh Institutional Review Board approved the study (STUDY19050099) and written informed consent was provided by all participants or their surrogates in accordance with the Declaration of Helsinki.
Data Availability
The data that support the findings will be available on Github at https://github.com/MicrobiomeALIR following acceptance of this manuscript for publication.
Author contributions
Georgios D. Kitsios: Conceptualization, Methodology, Investigation, Software, Formal analysis, Writing - Original Draft, Writing - Review & Editing, Visualization, Supervision
Daniel Kotok: Conceptualization, Methodology, Investigation, Software, Formal analysis, Writing - Original Draft, Writing - Review & Editing, Visualization, Supervision
Haopu Yang: Methodology, Investigation, Writing - Review & Editing, Visualization
Malcolm Finkelman: Methodology, Investigation, Writing - Review & Editing
Yonglong Zhang: Investigation, Writing - Review & Editing
Noel Britton: Investigation, Writing - Review & Editing
Rui Guo: Investigation, Writing - Review & Editing
John W. Evankovich: Investigation, Writing - Review & Editing
William Bain: Investigation, Writing - Review & Editing
Faraaz Shah: Investigation, Writing - Review & Editing
Yingze Zhang: Investigation, Writing - Review & Editing
Panayiotis V. Benos: Methodology, Investigation, Software, Writing - Review & Editing, Visualization
Bryan J. McVerry: Conceptualization, Methodology, Investigation, Writing - Original Draft, Writing - Review & Editing, Supervision
Alison Morris: Conceptualization, Methodology, Investigation, Writing - Original Draft, Writing - Review & Editing, Supervision
Conflict of interest
GDK receives research funding from Karius, Inc. The other authors have no competing interests to declare. MF and YZ are employed by the manufacturer of the BDG assay (Associates of Cape Cod, Inc, East Falmouth, MA, USA) used in this study (“Fungitell®”).
Acknowledgements
We would like to thank all members of the research team of the Acute Lung Injury Registry (ALIR) and Biospecimen Repository at the University of Pittsburgh, the medical and nursing staff in the Medical Intensive Care Unit at the University of Pittsburgh Medical Center, and all patients and their families for participating in this research project. Off-label use disclaimer: “Fungitell® is cleared for use as an adjunct test in the diagnosis of invasive fungal disease. The data presented herein are for research purposes only.”
Footnotes
↵* co-first authors
List of abbreviations
- ARDS
- Acute Respiratory Distress Syndrome
- ALIR
- Acute Lung Injury and Biospecimen Repository
- AKI
- Acute kidney injury
- CHF
- Congestive heart failure
- ICU
- Intensive care unit
- IL
- Interleukin
- ITS
- Internal transcribed spacer
- IFI
- Invasive fungal infection
- P:F ratio
- Ratio of partial pressure of arterial oxygen [P] and fraction of inspired oxygen [F]
- RAGE
- Receptor of advanced glycation end-products
- SOFA
- Sequential Organ Failure Assessment
- ST-2
- Suppression of tumorigenicity-2
- TNFR1
- Tumor necrosis factor receptor 1
- UPMC
- University of Pittsburgh Medical Center
- VFD
- Ventilator-free days