Impact Factor 4.019

The world's most-cited Microbiology journal

Original Research ARTICLE

Front. Microbiol., 10 July 2018 | https://doi.org/10.3389/fmicb.2018.01413

Respiratory Microbiome Profiling for Etiologic Diagnosis of Pneumonia in Mechanically Ventilated Patients

Georgios D. Kitsios1,2*, Adam Fitch2, Dimitris V. Manatakis3, Sarah F. Rapport1, Kelvin Li2, Shulin Qin1,2, Joseph Huwe1, Yingze Zhang1, Yohei Doi4, John Evankovich1, William Bain1, Janet S. Lee1, Barbara Methé1,2, Panayiotis V. Benos3, Alison Morris1,2,5† and Bryan J. McVerry1,2†
  • 1Division of Pulmonary, Allergy and Critical Care Medicine, Department of Medicine, University of Pittsburgh School of Medicine and University of Pittsburgh Medical Center, Pittsburgh, PA, United States
  • 2Center for Medicine and the Microbiome, University of Pittsburgh, Pittsburgh, PA, United States
  • 3Department of Computational and Systems Biology, University of Pittsburgh, Pittsburgh, PA, United States
  • 4Division of Infectious Diseases, University of Pittsburgh Medical Center, Pittsburgh, PA, United States
  • 5Department of Immunology, University of Pittsburgh School of Medicine, Pittsburgh, PA, United States

Etiologic diagnosis of bacterial pneumonia relies on identification of causative pathogens by cultures, which require extended incubation periods and have limited sensitivity. Next-generation sequencing of microbial DNA directly from patient samples may improve diagnostic accuracy for guiding antibiotic prescriptions. In this study, we hypothesized that enhanced pathogen detection using sequencing can improve upon culture-based diagnosis and that certain sequencing profiles correlate with host response. We prospectively collected endotracheal aspirates and plasma within 72 h of intubation from patients with acute respiratory failure. We performed 16S rRNA gene sequencing to determine pathogen abundance in lung samples and measured plasma biomarkers to assess host responses to detected pathogens. Among 56 patients, 12 patients (21%) had positive respiratory cultures. Sequencing revealed lung communities with low diversity (p < 0.02) dominated by taxa (>50% relative abundance) corresponding to clinically isolated pathogens (concordance p = 0.009). Importantly, sequencing detected dominant pathogens in 20% of the culture-negative patients exposed to broad-spectrum empiric antibiotics. Regardless of culture results, pathogen dominance correlated with increased plasma markers of host injury (receptor of advanced glycation end-products-RAGE) and inflammation (interleukin-6, tumor necrosis factor receptor 1-TNFR1) (p < 0.05), compared to subjects without dominant pathogens in their lung communities. Machine-learning algorithms identified pathogen abundance by sequencing as the most informative predictor of culture positivity. Thus, enhanced detection of pathogenic bacteria by sequencing improves etiologic diagnosis of pneumonia, correlates with host responses, and offers substantial opportunity for individualized therapeutic targeting and antimicrobial stewardship. Clinical translation will require validation with rapid whole meta-genome sequencing approaches to guide real-time antibiotic prescriptions.

Introduction

Severe pneumonia is a leading cause of hospitalization and death among adults in the US, often requiring admission to an intensive care unit (ICU) (Chalmers et al., 2011; Barrett et al., 2014; Jain et al., 2015; Valley et al., 2015). While appropriate antibiotic therapy is the cornerstone of pneumonia management, etiologic pathogen diagnosis with current culture-based microbiologic tests is often negative in patients with a clinical picture of pneumonia (Jain et al., 2015) or requires long incubation periods (∼3 days) to provide actionable results (Zumla et al., 2014; Jain et al., 2015; Vincent et al., 2015). Consequently, antibiotic prescriptions for severe pneumonia are empiric and typically include two or three broad-spectrum agents prescribed for seven or more days (Mandell et al., 2007; Kalil et al., 2016). This “one-size-fits-all” practice is hazardous for individual patients, who may receive insufficient or disproportionately intense antibiotics, and further contributes to antibiotic resistance, a global health threat (Laxminarayan et al., 2013; Modi et al., 2014; Kitsios et al., 2017).

Next-generation sequencing (NGS) of microbial DNA extracted directly from patient samples without the need for ex vivo organismal growth may overcome shortcomings of culture-based diagnosis. By sequencing either amplified bacterial marker genes (typically the 16S rRNA gene [16S sequencing]) or whole metagenomes, NGS provides comprehensive profiling of resident microbial communities with relative abundance information of the constituent bacteria, regardless of whether they are alive, dead or fastidious. With improvements in fidelity and accessibility of whole metagenome sequencing, identification of non-bacterial microbes (viruses, fungi and parasites) (Kuczynski et al., 2011) is becoming increasingly feasible and holds promise for clinical utility in the near future (Judge et al., 2015). Despite the theoretical advantages of NGS, the technology has not yet been validated as a diagnostic tool to guide antimicrobial prescriptions in the ICU.

To examine the clinical validity (Khoury et al., 2003) of 16S sequencing for etiologic diagnosis of bacterial pneumonia in patients requiring mechanical ventilation, we conducted this proof-of-concept study with the Microbiome Cohort in Acute Lung Injury Registry (MICALIR) at the University of Pittsburgh Medical Center (UPMC), assessing the upper and lower respiratory tract microbiome composition and its association with clinical diagnoses, host responses, and clinical outcomes.

Materials and Methods

This study is reported in compliance with the Strengthening the Reporting of Observational studies in Epidemiology (STROBE) statement (details presented in Supplementary Table 1) (Elm et al., 2007). Extensive methods regarding clinical data recording, sample collection, experimental protocols and statistical analyses are presented in the Supplement. Primary data and statistical code for replication of our findings are also provided in the online Supplement and are also available for download at https://github.com/MicrobiomeALIR/Resp_Microbiome_Profiles_Pneumonia.

Study Design and Participants

We conducted a prospective cohort study from June 2015 – March 2017 enrolling consecutive adult patients in the Medical ICU with acute respiratory failure within 72 h of intubation. Eligible patients were 18 years or older with acute respiratory failure requiring mechanical ventilation via endotracheal intubation. Exclusion criteria included inability to obtain informed consent, presence of tracheostomy, or mechanical ventilation for more than 72 h. Given our focus on etiologic diagnosis of bacterial pneumonia, we divided patients into culture-positive and culture-negative cases, based on clinical microbiologic results of respiratory specimens. We considered specimens obtained within 48 h of research sample timing acquisition so that such specimens would be reflective of the same infectious process being studied by NGS. Clinical cultures were obtained at the discretion of the treating physicians who were not involved in the MICALIR study. We considered microbiologic cultures of respiratory specimens [sputum, endo-tracheal aspirates (ETAs), or bronchoalveolar lavage – (BAL)] as positive when pathogenic bacterial species had been isolated by the clinical laboratory and treating physicians covered these bacteria with antibiotics (i.e., clinically these bacteria were not considered to be airway colonizers). Culture-negative cases were defined as those in which no organismal growth was observed or presence of “normal respiratory flora” was reported, as per standard clinical practices. We also recorded results of respiratory viral panel testing performed in nasopharyngeal swabs or respiratory specimens. The University of Pittsburgh Institutional Review Board approved the MICALIR study and written informed consent was provided by all participants or their surrogates in accordance with the Declaration of Helsinki.

Clinical Data Collection

We collected prospective clinical data on participants including age, gender, body mass index and history of smoking, comorbid conditions, such as diabetes and chronic obstructive pulmonary disease (COPD), on the day of enrollment. Physiological and laboratory variables such as PaO2 to FiO2 ratio, levels of positive end-expiratory pressure, plateau pressure, systolic blood pressure, white blood cell count etc. were obtained from the medical record by recording the physiologically worse value within a 24 h period on the day of enrollment (e.g., lowest blood pressure or highest creatinine value). We measured modified sequential organ failure assessment (SOFA) scores (we did not include the neurologic component of SOFA score because all patients were intubated and receiving sedative medications) by using the physiologically worse values within 24 h of enrollment. We recorded antibiotics and vasopressors administered during the first week of ICU course from intubation. A consensus committee of clinical experts (GDK, JE, WB, JSL, AM, BJM) reviewed all available data for the enrolled patients and performed retrospective classifications of the etiology and severity of their acute respiratory failure. Classifications were performed without knowledge of microbiome sequencing or host biomarker data and included sepsis, acute respiratory distress syndrome (ARDS), pneumonia or aspiration, and intubation for airway protection without risk factors for ARDS per established criteria (Gong et al., 2005; ARDS Definition Task Force et al., 2012; Singer et al., 2016; Neto et al., 2016). Prospective clinical outcomes included 30-day mortality, duration of ICU stay, acute kidney injury (Mehta et al., 2007), incident shock (defined as need for vasopressors) and ventilator-free days (VFD) (Kalanuria et al., 2014).

Sample Collection

Immediately after enrollment, we obtained baseline samples for microbiome analyses of the oral and lung communities, with swabs of the base of the tongue and ETA, respectively. Simultaneous blood samples (10cc) were collected for centrifugation and separation of plasma. All samples were frozen and stored at -80°C.

Laboratory Analyses

We extracted bacterial DNA directly from oral swabs and ETAs and amplified the V4 hypervariable region of the 16S rRNA gene for sequencing on the Illumina MiSeq platform (Morris et al., 2013). We also performed qPCR of the V3-V4 region of the 16S rRNA gene to obtain absolute bacterial loads in each sample (Liu et al., 2012). For plasma biomarker analyses, we constructed a custom Luminex multi-analyte panel (R&D Systems, Minneapolis, MI, United States) (McKay et al., 2017) targeting biomarkers associated with pneumonia diagnosis (procalcitonin) and acute lung injury outcomes (RAGE: receptor of advanced glycation end-products, IL-6: interleukin-6, IL-8: interleukin-8, sTNFR1: soluble tumor necrosis factor receptor 1) (Narvaez-Rivera et al., 2012; Calfee et al., 2014; Famous et al., 2017). Clinical microbiologic cultures were processed per standard procedures as described in the Supplement.

Data Processing and Statistical Analysis Methods

From derived 16S sequences, we applied a custom pipeline for Operational Taxonomic Units (OTUs-taxa) classification (Supplement). We calculated descriptive statistics of baseline characteristics and performed non-parametric comparisons using the R software (R Foundation for Statistical Computing, 2016). We performed ecological analyses of alpha diversity (richness-Shannon and evenness-Dominance), beta diversity (Bray–Curtis dissimilarity index), and taxonomic descriptions between culture-positive and negative cases with the Quantitative Insights in Microbial Ecology software (QIIME) and the R vegan package (Dixon, 2003; Caporaso et al., 2010). Beta-diversity comparisons with permutation analysis of variance (Permanova at 1000 permutations) were visualized with non-metric multidimensional scaling (NMDS). From the taxonomic composition of the reference-standard culture-positive cases and the available literature on the composition of the healthy lung microbiome (Morris et al., 2013; Segal et al., 2016), we operationally defined taxa as “pathogenic” (when corresponding to clinically relevant bacterial species isolated in cultures) vs. “oral-origin” taxa, for those taxa that have been included in the supraglottic pneumotype of the lung microbiome created by microaspiration of oral bacteria (Segal et al., 2016). From observed taxonomic profiles and distribution of pathogen and oral taxa abundances, we defined pathogen dominance or oral taxa dominance as relative abundance of >50%, respectively. We compared log-transformed concentrations of host-response biomarkers and pathogen or oral taxa dominance with linear regression models, adjusted for culture results and history of COPD.

Integrative Analysis and Modeling of Microbiome and Clinical Data

To move beyond simple correlations and comprehensively examine which microbiome variables are directly linked to clinical variables, we used Probabilistic Graphical Models (PGMs). PGMs can estimate and graphically represent the complex relationships of large numbers of variables that interact with each other, allowing for the discovery of direct links between variables based on their conditional dependencies. We used the CausalMGM (Causal Mixed Graphical Model) R package1, a novel algorithm that can accurately identify the underlying graphical model structure over mixed data types (continuous and discrete) (Sedgewick et al., 2016; Raghu et al., 2017).

Results

Cohort Descriptive Data

Fifty six patients were enrolled (mean age 56 years, 61% men), 12 (21%) with positive cultures for common respiratory pathogens (Table 1). Culture-positive and negative patients had similar distribution of comorbid conditions, severity of illness scores, mechanical ventilation parameters, laboratory values and clinical outcomes (Table 1). Empiric antibiotics were prescribed for 54/56 (96%) of patients at the time of enrollment, with a median exposure of two different classes of antibiotics (gram-positive, gram-negative, and atypical coverage).

TABLE 1
www.frontiersin.org

TABLE 1. Baseline characteristics and clinical outcomes of patients enrolled in the Microbiome Cohort of the Acute Lung Injury Registry (MICALIR) study, categorized as patients with positive or negative respiratory cultures.

16S Sequencing Results

We analyzed a total of 110 clinical microbiome samples (56 ETAs and 54 oral swabs) and 71 experimental control samples. Clinical samples produced 3835 (standard deviation [SD] = 1329) reads (high-quality 16S sequences) on average, whereas negative experimental controls generated relatively few reads (mean = 136 [SD = 146]) (p < 10-16) (Figure 1A) and were compositionally dissimilar to clinical samples by Bray–Curtis indices (Figure 1B).

FIGURE 1
www.frontiersin.org

FIGURE 1. Comparisons of sequencing results between clinical samples and experimental controls. (A) N of reads by sample type in the MICALIR study. Clinical samples (lung and oral) produced about 30 times more reads than negative control samples (p < 10-16) [ExtractionNeg: negative control for DNA extraction experiments; TrachAspControls: sterile left-over saline from the one used to instill into the endotracheal tube for suctioning endotracheal aspirates; PCRNeg: PCR negative controls]. A rarefaction level of 850 reads was selected for alpha diversity analyses, which excluded small numbers of clinical samples for these analyses. (B) Non-metric multidimensional scaling (NMDS) plot of Bray–Curtis dissimilarity indices between clinical samples (lung and oral) versus negative controls (extraction, PCR and tracheal sampling procedure controls) and positive controls. Experimental control samples were compositionally markedly dissimilar to clinical samples by Bray–Curtis indices (Permanova p-value = 0.001). No taxa detected in negative controls were filtered from downstream analyses.

Lung Microbial Community Profiles

Lung communities in culture-positive subjects had significantly lower alpha diversity (richness [p = 0.02] and evenness [p = 0.04]) compared to culture-negative subjects (Figure 2A). Culture-negative communities had a wide distribution of alpha diversity, ranging from very low richness as in culture-positive cases (Shannon = 0-1) to high alpha diversity (Shannon > 2.9) in the range of the healthy lung microbiome, despite the fact that these patients were exposed to broad-spectrum antibiotics (Morris et al., 2013). Culture-positive communities had overall significantly different taxonomic composition by Bray–Curtis indices compared to culture-negative communities (Permanova p = 0.003) (Figure 2B). On the other hand, lung communities from patients with a history of COPD (a disease process that is known to affect lung microbiome composition) (Sze et al., 2014) had modest differences in alpha and beta diversity compared to patients without COPD (Supplementary Figure 3).

FIGURE 2
www.frontiersin.org

FIGURE 2. Culture-positive lung communities have significantly lower alpha diversity and are compositionally distinct from culture-negative communities. (A) Alpha diversity comparisons between culture-positive and culture-negative cases showed statistically significant differences in Shannon and Dominance indices, representing richness and evenness, respectively, indicating that culture-positive communities are being dominated by fewer bacterial species. Culture-negative communities demonstrated a wide range of alpha diversity. (B) Beta diversity (Bray–Curtis dissimilarity indices) comparisons between respiratory culture-positive (Resp Cx Positive) and negative (Resp Cx negative) cases showed significant differences (Permanova p-value = 0.003), but certain culture-negative samples overlapped with the culture-positive cluster indicating underlying compositional similarity.

By examining the taxonomic composition of the 12 culture-positive samples (Figure 3A), we found that 16S sequencing detected taxa concordant to the clinically isolated bacterial pathogens in 11/12 samples (e.g., Staphylococcus genera for S. aureus and Enterobacteriaceae for Klebsiella pneumoniae). In 9 (75%) samples, the concordant taxa were the most abundant organisms in their respective communities. In two other cases (cases 10 and 11, Figure 3A), sequencing revealed that taxa corresponding to clinical isolates (Klebsiella and Staphylococcus, respectively) had low abundance in their communities, which were dominated by potentially pathogenic taxa undetected by clinical cultures (Enterococcus and Fusobacterium, respectively) (Grupper et al., 2009; Johannesen et al., 2016; Kelly et al., 2016). Although definitive causal inferences about the clinical impact of these bacteria are not possible, in the case of Fusobacterium dominance (case 11), vancomycin monotherapy targeted against the cultured methicillin-resistant S. aureus (MRSA) failed to produce a clinical response, whereas improvement ensued after empiric addition of piperacillin/tazobactam which in fact is effective against the fastidious anaerobe Fusobacterium (Johannesen et al., 2016).

FIGURE 3
www.frontiersin.org

FIGURE 3. Pathogen dominance detection in culture-positive (A) and negative lung samples (B). Taxonomic composition is shown as stacked bar-graphs, with each bar representing a patient’s community, with taxa colored individually and heights of component bars corresponding to relative abundance of each taxon. In culture-positive samples (A), the clinically isolated organisms by routine microbiologic cultures are spelled out vertically in each bar (Methicillin-resistant S. aureus in cases 1, 2, 11; Methicillin-sensitive S. aureus in cases 3, 4, 5; Haemophilus influenza in case 5; Pseudomonas aeruginosa in case 6, Klebsiella pneumoniae in cases 7, 8, and 10, Escherichia coli in case 9; Serratia marcescens in case 12. In cases 1–9, the most abundant taxon corresponded to the clinically isolated pathogen (culture-concordance). In three cases (10–12), there was discordance between cultures and sequencing (i.e., the most abundant organism was not the one isolated by cultures. In cases 10 and 11, the clinically isolated Klebsiella pneumoniae and S. aureus corresponded to a minority of concordant reads in these communities that were dominated by Enterococcus and Fusobacterium taxa, respectively. In case 12, sequencing showed dominance by Haemophilus taxa whereas cultures isolated Serratia marcescens. Among culture-negative samples (B), 20% were dominated by pathogenic taxa similar to the ones detected in culture-positive cases, and the remaining samples showed high abundance of oral bacteria. In six cases highlighted with the “#” symbol, respiratory viral panels of the nasopharynx or respiratory specimens were positive (for influenza, respiratory syncytial virus, metapneumovirus or parainfluenza virus). The “other” taxonomic assignment corresponds to multiple genera not corresponding to “pathogens” or “oral taxa” lumped together for display purposes. H.Flu, Haemophilus Influenza.

Based on the overall concordance between culture-isolated bacteria and taxonomic abundance by sequencing, we operationally defined the most abundant taxa in culture-positive communities as “pathogens” (Supplementary Table 2 and Supplementary Figure 4). Looking then into the culture-negative cases that represent a diagnostic “black-box,” we found that 9/44 (20%) communities were dominated by pathogenic taxa (e.g., Staphylococcus or Pseudomonas genera) suggesting a specific etiology that was undefined by standard culture-based methods (Figure 3B). The remaining culture-negative samples were populated by the most common members of the supraglottic pneumotype of the lung microbiome, such as Prevotella, Veillonella, and Streptococcus taxa (“oral-origin” taxa) (Segal et al., 2016; Dickson et al., 2017). In four culture-negative cases, clinical viral panel testing was positive for respiratory viruses (influenza, parainfluenza, metapneumovirus, and respiratory syncytial virus) without identification of common bacterial pathogens (Figure 2B), and thus no indication of the commonly suspected bacterial super-infection that leads to empiric antibiotic courses in patients with viral pneumonia (van der Sluijs et al., 2010).

After examining different thresholds of pathogen abundance (Supplementary Figure 5), we found that a relative abundance of ∼50% (community dominance) was strongly associated with bacterial culture results with an odds-ratio (OR) of 48.2 (95% CI: 2.6-898.9, p = 0.009) (Figure 4). In contrast, oral taxa dominance was strongly “protective” against respiratory culture positivity (OR = 0.01 [95% CI: 0.0008–0.26, p = 0.004]).

FIGURE 4
www.frontiersin.org

FIGURE 4. Dominance of lung communities by pathogens or oral taxa was strongly associated with respiratory culture results. (A) Pathogen dominance (>50% abundance) was strongly associated with concordant pathogen culture-positivity [Fisher’s odds ratio (OR) with continuity correction and associated 95% confidence interval shown]. The reference standard here was chosen to be concordant pathogen positivity, as the sequencing diagnostic test would be clinically valid if able to detect the same organism as our current reference standard of cultures. (B) Oral taxa dominance (>50% abundance) practically eliminated the odds of culture-positivity by any pathogen (OR = 0.01). The reference standard here was defined as “any pathogen,” given that oral taxa are not generally considered as pathogens or speciated by routine microbiologic cultures, and this comparison aims to assess the negative predictive value of high oral abundance in a lung community for ruling out culture positivity by any pathogenic bacteria.

We also quantified the bacterial load in lung communities by 16S qPCR and found no significant differences in absolute number of 16S rRNA gene copies between culture-positive and negative samples, underscoring the fact that several culture-negative samples were carrying high bacterial loads that were clinically undetected. There was a wide range of 16S rRNA copies across samples (range from 34 to 50,390,355 copies) (Figure 5), and we did not identify a diagnostic threshold of absolute pathogen abundance for association with culture positivity. Thus, in our dataset relative abundance of pathogens, which is reflective of within community microbial dynamics, appeared as a stronger predictor of culture positivity compared to absolute pathogen abundance as quantified by 16S qPCR copies.

FIGURE 5
www.frontiersin.org

FIGURE 5. Bacterial load in culture-positive and negative lung communities as quantified by 16S rRNA gene qPCR. No significant differences between culture-positive and negative samples were found. The Y-axis showing number of 16S rRNA gene copies per sample is square-root transformed for visualization purposes.

Oral Microbial Community Profiles

Oral microbiome profiles closely reflected the patterns observed with the lung microbiome, with lower alpha diversity in subjects with positive respiratory cultures compared to culture-negative subjects (p = 0.002) (Figure 6). We further examined the taxonomic composition of the oral communities in subjects with pathogen dominance in their lung communities (n = 19) and found that the oral communities were dominated by the same pathogen (Supplementary Figure 6) in seven subjects (36%), implicating colonization of the oral cavity as a potential source of pneumonia pathogens and suggesting the potential utility of oral sample sequencing for pneumonia diagnosis.

FIGURE 6
www.frontiersin.org

FIGURE 6. Alpha and beta diversity comparisons in oral and lung communities, stratified by respiratory specimen culture positivity. (A) Alpha diversity comparisons in oral communities, showing statistically significantly lower richness and evenness in culture-positive samples compared to culture-negative ones (p = 0.002). (B) Bray–Curtis dissimilarity indices comparison in 4 groups: red circles for oral communities of culture-positive samples, blue circles for oral communities of culture-negative samples, red triangles for lung communities of culture-positive samples, blue triangles for lung communities of culture-negative samples. Permanova indicates significant differences overall, but oral and lung communities are overlapping when stratified by respiratory sample culture positivity, indicating that oral communities in culture-positive cases were taxonomically more similar to their corresponding culture-positive lung communities, rather than the oral communities of culture-negative cases.

Host Responses to Lung Microbiota and Clinical Outcomes

Lung community pathogen dominance was significantly associated with higher levels of circulating inflammatory cytokines (IL-6: p = 0.007; sTNFR1: p = 0.03) and epithelial injury biomarkers (RAGE: p = 0.02) (Figure 7). For IL-6, the association remained significant after adjusting for culture results (p = 0.03), suggesting that such pathogenic bacteria induce host inflammation regardless of their viability or ability to grow in cultures at the time of sample acquisition (Table 2). These associations remained statistically significant and with larger effect sizes when adjusted for history of COPD, a disease process that could also confound host-microbiome associations (Table 2). We did not find any significant associations between important patient-centered outcomes (mortality, shock, VFD and length of ICU stay) and pathogen dominance in the lung communities.

FIGURE 7
www.frontiersin.org

FIGURE 7. Pathogen abundance in lung communities was associated with higher levels of lung epithelial injury and inflammation. Associations between plasma biomarkers of injury (RAGE) and inflammation (IL-6, IL-8, and sTNFR1) with pathogen dominance in lung communities (relative abundance > 50%) are shown for the entire cohort. Biomarker concentrations (pg/ml) are shown in logarithmic scale. Culture-positive and negative samples are shown with filled and open circles, respectively. Statistically significant p-values are shown in boxes. No significant association with procalcitonin levels was found.

TABLE 2
www.frontiersin.org

TABLE 2. Associations between pathogen abundance and host-response biomarkers in unadjusted linear regression models, and in adjusted models for culture positivity or history of COPD.

Network Analyses

Probabilistic graph interrogation of our dataset using CausalMGM identified five microbial taxa and the levels of hemoglobin to be the most informative values for respiratory culture positivity, as well as a composite variable measuring pathogen taxa abundance (Figure 8). These findings extended our univariate taxonomic analyses in that they identified the taxa that are directly linked to positive cultures (not simply correlated) and highlighted sequenced bacterial taxa as the strongest explanatory variables of culture positivity. To examine the ability of 16S taxonomic data alone to predict culture positivity, we used the Markov blanket around the culture-positivity variable as a feature selection method (Huang et al., 2015). The taxonomy-based classifier yielded mean accuracy of 82.3% (SD = 7%) (Table 3), indicating proof-of-concept utility for use of sequencing in clinical practice for predicting culture results, if sequencing results were available real-time.

FIGURE 8
www.frontiersin.org

FIGURE 8. Pathogen abundance by sequencing was the strongest correlate of culture-positivity by probabilistic graphical modeling. Network analyses included clinical (yellow), individual 16S taxa (blue), composite pathogen taxa abundance (orange) and biomarker (purple spheres) variables. First and second neighbors around the clinically important variable of respiratory culture positivity (highlighted by a dashed square) are shown. Edges (links) between variables are shown in red for positive and in green for negative correlations. The thickness of the edges is proportional to the stability metric of the detected associations. Respiratory culture positivity was positively associated with the composite variable of pathogen abundance, Enterobacteriaceae, Haemophilus, Escherichia, and Enterococcus taxa, and negatively associated with Prevotella abundance and hemoglobin levels. COPD, chronic obstructive pulmonary disease; Hgb, hemoglobin; P/F Ratio, Pa02/Fi02 ratio; BMI, body mass index.

TABLE 3
www.frontiersin.org

TABLE 3. Predictive taxa of culture positivity along with their numbers of appearance in the 10 Markov blankets.

Discussion

Our prospective study in mechanically ventilated patients provides the largest examination to date of the clinical validity of NGS for etiologic diagnosis of index severe pneumonia. By utilizing minimally invasive samples, we detected high abundance of pathogenic bacteria in both culture-positive and negative cases of pneumonia, which could allow for timely etiologic pathogen identification and antibiotic adjustments once technology evolves to allow rapid NGS. Our analyses highlighted overall community structure differences associated with culture-positivity and uncovered distinct differences in the host immune responses to dominant bacterial taxa in the lungs. Probabilistic graphical models that inclusively consider the wide range of variables in our dataset provided insights on the taxonomic and clinical variables directly linked to respiratory culture positivity. Our results indicate that respiratory microbiome profiles may provide clinically meaningful and actionable information above that currently afforded by standard microbiologic cultures.

16S sequencing provided important insights into the etiology of both culture-positive and negative bacterial pneumonias. In culture-positive cases, taxa concordant to clinically isolated pathogens dominated the respective communities in the vast majority of cases, consistent with the concept of community collapse during infection (Dickson et al., 2014). In the few exceptions of discordance between dominant pathogen and culture results (cases 10–12, Figure 3A), sequencing did find the cultured pathogen although it was not dominant and may have uncovered co-primary or alternative pathogens not clinically considered. Thus, even in cases of supposed diagnostic certainty, cultures may only capture a small portion of the underlying bacterial community and miss dominant bacteria that may have important therapeutic implications.

Culture-negative samples represent a diagnostic “black-box” in clinical practice, and in up to 75% of pneumonia cases no bacterial pathogen is ever isolated (Jain et al., 2015), similar to the culture-negative rate seen in our cohort. We found that approximately 20% of our culture-negative samples actually contained a predominance of pathogenic taxa similar to those observed in culture-positive specimens (e.g., Staphylococcus or Pseudomonas). In contrast, eighty percent of culture-negative samples without a dominant pathogen had high oral bacteria abundance, similar to what is often found in the lung microbiome and suggesting ongoing micro-aspiration of oral secretions around endotracheal tube cuffs in these patients (Kitsios and McVerry, 2018). In general, oral-origin taxa were associated with lower levels of host immune responses compared to pathogenic taxa. Nevertheless, oral taxa dominance cannot be uniformly considered as clinically innocuous, given significant interindividual variability in host responses and absolute bacterial load, making context-specific integration of clinical and microbiome data necessary for clarifying the clinical importance of individual oral bacteria.

Improvements in the etiologic diagnosis of pneumonia offered by NGS can translate into measurable benefits as defined by appropriate antibiotic targeting. For example, when a culprit pathogen is identified to dominate a community, early antibiotic tailoring against this pathogen becomes feasible and avoids the hazards of empiric broad-spectrum antibiotic courses. Furthermore, in up to 40% of culture-negative communities, no pathogenic sequence abundance was found and communities had high alpha diversity metrics in the range of the healthy human lung microbiome. With no indication of bacterial pneumonia in these patients by 16S sequencing, early discontinuation of antibiotics could substantially decrease their cumulative antibiotic exposure. Nonetheless, the 16S sequencing technique we utilized cannot yet be applied clinically and does not provide species-level or antibiotic resistance information. The advent of ultra-rapid sequencing devices and bioinformatics pipelines offer the capacity for whole metagenome sequencing in a matter of hours, theoretically allowing for bedside pathogen identification including viruses and fungi and antibiotic resistance prediction (Naccache et al., 2014; Judge et al., 2015; Schmidt et al., 2017). In a recent proof-of-concept study, sequencing with a portable, point-of-care device (Oxford Nanopore Technologies) was able to identify the culprit pathogens before a clinical microbiology lab, underlining that NGS may alleviate the time-consumption problem with traditional culture techniques (Pendleton et al., 2017). Whole metagenome sequencing techniques remain to be optimized for pathogen detection from biological samples containing substantial amounts of human DNA.

Apart from the technological challenges outlined above, clinical research in the field also has to overcome the formidable challenge of comparing a new sensitive test (NGS) against a standard-of-care (microbiologic cultures) that is not a gold-standard (Rutjes et al., 2007). Given that clinical cultures can effectively grow only a subset of cultivable bacteria (Venkataraman et al., 2015), conventional sensitivity/specificity analyses of NGS become obsolete. To overcome such diagnostic challenges, the clinical reference standard has to be refined, with synthesis of multi-level data (clinical, radiographic, biomarkers etc.) to be combined in a “construct gold-standard” for pneumonia diagnosis with supervised (involving expert input) or unsupervised classification methods (Rutjes et al., 2007). Our PGM analyses highlighted sequencing variables as the strongest predictors of culture-positivity. Through iterative training of machine-learning algorithms comparing sequencing profiles to “construct gold-standard” cases of pneumonia, a diagnostic, sequencing-based classifier can be developed for clinical use (Kitsios, 2018). Ultimately, demonstration of clinical efficacy of NGS-based testing for improving patient outcomes in randomized clinical trials will provide the evidentiary support necessary for clinical adoption.

The striking similarities between the oral and lung communities provided us with further understanding of the bacterial topography of the intubated respiratory tract (Kelly et al., 2016; Dickson et al., 2017; Kitsios et al., 2017). Although not part of routine clinical screening, oral cavity colonization by pathogenic gram-negative bacteria is a well-known harbinger of pneumonia in hospitalized patients (Dickson, 2016). Our analyses uncovered oral communities with dysbiosis (very low alpha diversity) in patients with positive sputum or BAL cultures. With cross-sectional examination of the respiratory microbiome, we could not assess for temporal bacterial immigration along the oro-tracheal tract. However, the observed patterns of oral and lung co-dominance by the same pathogens strongly indict the mouth as the originating pool of pathogens, which notably occurred despite routine oropharyngeal decontamination with chlorhexidine in our ICU (Price et al., 2014). Thus, sequencing of non-invasive oral swabs offers an attractive option for screening for pathogen colonization and for plausible pathogen detection in patients with pneumonia when lower respiratory specimens are unavailable.

The significant associations between sequencing-detected abundance of pathogens and host biomarkers of injury and inflammation (RAGE, IL-6, TNFR1) provided further validation of the biological and clinical relevance of sequencing profiles for pneumonia diagnosis. Consistent with other recent reports correlating lung microbiota composition with concurrently measured plasma biomarkers in critically ill patients (Panzer et al., 2017), our findings indicate that inter-individual heterogeneity in patients with severe pneumonia and ARDS may be explained on the basis of host-microbiome interactions, offering a new avenue for identifying clinical subphenotypes for personalized therapies (Calfee et al., 2014). Our linear regression models identified history of COPD as a significant confounder of host–microbiome associations, and thus COPD should be taken into account in future statistical modeling of these associations.

Our study has several limitations. It is a single-center design and is limited by available sample size, despite being the largest study of its kind. Consequently, the panel of pathogenic bacteria detected by sequencing is limited (Srinivasan et al., 2015), and generalizing to other critically ill populations and bacteria should be cautious (Calfee et al., 2014). Streptococcus taxa were commonly found in our patients with oral taxa-predominant bacterial communities, and without species-level information, we cannot exclude that certain 16S sequences belonged to pathogenic Streptococcus pneumoniae (albeit not recovered in any culture). However, in sensitivity analyses in which we coded Streptococcus sequences as “pathogens,” the associations between pathogen abundance and host responses were attenuated and no longer statistically significant (data not shown), suggesting that Streptococcus taxa detected in our samples are less likely to induce host inflammation and injury.

Since we performed 16S sequencing in extracted DNA, we did not consider viruses (either DNA or RNA) or fungi in our analyses, which can be important pathogens in subgroups of hospitalized patients (Jain et al., 2015). Sequencing of the bacteriome is nonetheless of paramount importance, even when a viral pathogen is identified, as the concern of bacterial super- or co-infection is pervasive among clinicians. For practical and ethical purposes, we did not perform BAL for microbiome analyses, thus we could not assess for regional variability of communities. Our approach to rely on ETAs was informed by the recent clinical practice guidelines by the American Thoracic Society and the Infectious Disease Society of America, recommending non-invasive testing via ETAs over invasive testing (with bronchoscopy or blind bronchial sampling) for hospital-acquired or ventilator-associated pneumonia diagnosis (Kalil et al., 2016). These recommendations are based on evidence from five randomized clinical trials (Berton et al., 2014), including a multicenter clinical trial (Canadian Critical Care Trials Group, 2006) showing that there is no clinical advantage (in terms of mortality, length of ICU stay, duration of mechanical ventilation or antibiotic management changes) between invasive or non-invasive sampling practices. For research purposes, ETAs are considered minimal-risk for enrolled patients, and represent a generally accepted “summary statement” of the pulmonary microbiome in intubated patients (Kelly et al., 2016; Kitsios et al., 2017; Panzer et al., 2017). Thus, we considered a priori that our comparisons between clinical BAL samples and research ETAs are internally valid for assessing concordance between cultures and sequencing. Our results showing striking concordance between dominance pathogen by sequencing with culture positivity (Figure 4) further validated the comparability of ETAs vs. BAL for microbiologic diagnosis. Timing differences of sample acquisition between clinical cultures and research ETAs may account for some of the observed sequencing-culture discordances, although research samples were obtained within 24hrs of clinical samples in >75% cases. Finally, we did not compare NGS against evolving techniques for rapid pathogen identification (Zumla et al., 2014), because cultures remain the current standard-of-care and allowed us to perform pragmatic comparisons in clinically challenging cases.

In summary, our study provides proof-of-concept evidence that as NGS technologies develop further, they will become useful as pneumonia diagnostic tools in the ICU to allow for fast and reliable etiologic pathogen identification. Our results demonstrate the clinical relevance of comprehensive microbial community profiling to provide information beyond what is currently available in clinical practice by microbiologic cultures and to directly impact clinical decision-making. Further prospective study in larger patient cohorts will allow for meaningful integration of sequencing output in culture-independent definitions of pneumonia.

Data Availability Statement

All datasets generated and analyzed for this study are included in the Supplementary Files and are also publicly available for download at https://github.com/MicrobiomeALIR/Resp_Microbiome_Profiles_Pneumonia.

Author Contributions

GK, AM, and BJM: conception and design. GK, AF, SR, DM, KL, SQ, JH, YZ, YD, JE, WB, JL, BM, PB, AM, and BJM: acquisition, analysis or interpretation of data. GK, JE, WB, JL, AM, and BJM: clinical cohort phenotyping. GK, AF, SR, DM, KL, SQ, JH, YZ, YD, JE, WB, JL, BM, PB, AM, and BJM: drafting of work and/or revising for important intellectual content. GK, AF, SR, DM, KL, SQ, JH, YZ, YD, JE, WB, JL, BM, PB, AM, and BJM: final approval of version to be published; agreement to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.

Funding

This work was supported by National Institutes of Health [T32 HL007562 (GK); K23 HL139987 (GK); U01 HL098962 (AM); P01 HL114453 (BJM); R01 HL097376 (BJM); K24 HL123342 (AM); U01 HL137159 (DM and PB); R01 LM012087 (DM and PB); R01 HL086884 (JL and WB); HL136143 (JL and WB); F32 HL137258-01 (JE)].

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.01413/full#supplementary-material

DATA SHEET S1 | Metadata file.

DATA SHEET S2 | Combined metadata and taxonomic table.

DATA SHEET S3 | Taxa table at genus level.

DATA SHEET S4 | Extended Methods and Results.

Footnotes

  1. ^ https://github.com/benoslab/causalMGM

References

ARDS Definition Task Force, Ranieri, V. M., Rubenfeld, G. D., Thompson, B. T., Ferguson, N. D., Caldwell, E., et al. (2012). Acute respiratory distress syndrome: the Berlin definition. JAMA 307, 2526–2533. doi: 10.1001/jama.2012.5669

PubMed Abstract | CrossRef Full Text | Google Scholar

Barrett, M. L., Smith, M. W., Elixhauser, A., Honigman, L. S., and Pines, J. M. (2014). Utilization of Intensive Care Services, 2011. HCUP Statistical Brief # 185. Rockville, MD: Agency for Healthcare Research and Quality.

Google Scholar

Berton, D. C., Kalil, A. C., and Teixeira, P. J. Z. (2014). Quantitative versus qualitative cultures of respiratory secretions for clinical outcomes in patients with ventilator-associated pneumonia. Cochrane Database Syst. Rev. 10:CD006482. doi: 10.1002/14651858.CD006482.pub4

PubMed Abstract | CrossRef Full Text | Google Scholar

Calfee, C. S., Delucchi, K., Parsons, P. E., Thompson, B. T., Ware, L. B., Matthay, M. A., et al. (2014). Subphenotypes in acute respiratory distress syndrome: latent class analysis of data from two randomised controlled trials. Lancet Respir. Med. 2, 611–620. doi: 10.1016/S2213-2600(14)70097-9

CrossRef Full Text | Google Scholar

Canadian Critical Care Trials Group (2006). A randomized trial of diagnostic techniques for ventilator-associated pneumonia. N. Engl. J. Med. 355, 2619–2630. doi: 10.1056/NEJMoa052904

PubMed Abstract | CrossRef Full Text | Google Scholar

Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303

PubMed Abstract | CrossRef Full Text | Google Scholar

Chalmers, J. D., Taylor, J. K., Mandal, P., Choudhury, G., Singanayagam, A., Akram, A. R., et al. (2011). Validation of the infectious diseases society of America/American thoratic society minor criteria for intensive care unit admission in community-acquired pneumonia patients without major criteria or contraindications to intensive care unit care. Clin. Infect. Dis. 53, 503–511. doi: 10.1093/cid/cir463

PubMed Abstract | CrossRef Full Text | Google Scholar

Dickson, R. P. (2016). The microbiome and critical illness. Lancet Respir. Med. 4, 59–72. doi: 10.1016/S2213-2600(15)00427-0

CrossRef Full Text | Google Scholar

Dickson, R. P., Erb-Downward, J. R., Freeman, C. M., McCloskey, L., Falkowski, N. R., Huffnagle, G. B., et al. (2017). Bacterial topography of the healthy human lower respiratory tract. MBio 8:e02287-16. doi: 10.1128/mBio.02287-16

PubMed Abstract | CrossRef Full Text | Google Scholar

Dickson, R. P., Erb-Downward, J. R., and Huffnagle, G. B. (2014). Towards an ecology of the lung: new conceptual models of pulmonary microbiology and pneumonia pathogenesis. Lancet Respir. Med. 2, 238–246. doi: 10.1016/S2213-2600(14)70028-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Dixon, P. (2003). VEGAN, a package of R functions for community ecology. J. Veg. Sci. 14, 927–930. doi: 10.1111/j.1654-1103.2003.tb02228.x

CrossRef Full Text | Google Scholar

Elm, E., von Altman, D. G., Egger, M., Pocock, S. J., Gøtzsche, P. C., Vandenbroucke, J. P., et al. (2007). The strengthening the reporting of observational studies in epidemiology (STROBE) statement: guidelines for reporting observational studies. Lancet 370, 1453–1457. doi: 10.1016/S0140-6736(07)61602-X

CrossRef Full Text | Google Scholar

Famous, K. R., Delucchi, K., Ware, L. B., Kangelaris, K. N., Liu, K. D., Thompson, B. T., et al. (2017). Acute respiratory distress syndrome subphenotypes respond differently to randomized fluid management strategy. Am. J. Respir. Crit. Care Med. 195, 331–338. doi: 10.1164/rccm.201603-0645OC

PubMed Abstract | CrossRef Full Text | Google Scholar

Gong, M. N., Thompson, B. T., Williams, P., Pothier, L., Boyce, P. D., and Christiani, D. C. (2005). Clinical predictors of and mortality in acute respiratory distress syndrome: potential role of red cell transfusion. Crit. Care Med. 33, 1191–1198. doi: 10.1097/01.CCM.0000165566.82925.14

PubMed Abstract | CrossRef Full Text | Google Scholar

Grupper, M., Kravtsov, A., and Potasman, I. (2009). Enterococcal-associated lower respiratory tract infections: a case report and literature review. Infection 37, 60–64. doi: 10.1007/s15010-007-7123-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, G. T., Tsamardinos, I., Raghu, V., Kaminski, N., and Benos, P. V. (2015). T-ReCS: stable selection of dynamically formed groups of features with application to prediction of clinical outcomes. Pac. Symp. Biocomput. 20, 431–442.

PubMed Abstract | Google Scholar

Jain, S., Self, W. H., Wunderink, R. G., Fakhran, S., Balk, R., Bramley, A. M., et al. (2015). Community-acquired pneumonia requiring hospitalization among U.S. Adults. N. Engl. J. Med. 373, 415–427. doi: 10.1056/NEJMoa1500245

PubMed Abstract | CrossRef Full Text | Google Scholar

Johannesen, K., Dessau, R., Heltberg, O., and Bodtger, U. (2016). Bad news itself or just the messenger? The high mortality of Fusobacterium spp. infections is related to disseminated malignancy and other comorbidities. Eur. Clin. Respir. J. 3, 30287. doi: 10.3402/ecrj.v3.30287

PubMed Abstract | CrossRef Full Text | Google Scholar

Judge, K., Harris, S. R., Reuter, S., Parkhill, J., and Peacock, S. J. (2015). Early insights into the potential of the Oxford Nanopore MinION for the detection of antimicrobial resistance genes. J. Antimicrob. Chemother. 70, 2775–2778. doi: 10.1093/jac/dkv206

PubMed Abstract | CrossRef Full Text | Google Scholar

Kalanuria, A. A., Ziai, W., and Mirski, M. (2014). Ventilator-associated pneumonia in the ICU. Crit. Care 18:208. doi: 10.1186/cc13775

PubMed Abstract | CrossRef Full Text | Google Scholar

Kalil, A. C., Metersky, M. L., Klompas, M., Muscedere, J., Sweeney, D. A., Palmer, L. B., et al. (2016). Executive summary: management of adults with hospital-acquired and ventilator-associated pneumonia: 2016 clinical practice guidelines by the infectious diseases society of America and the American thoracic society. Clin. Infect. Dis. 63, 575–582. doi: 10.1093/cid/ciw504

PubMed Abstract | CrossRef Full Text | Google Scholar

Kelly, B. J., Imai, I., Bittinger, K., Laughlin, A., Fuchs, B. D., Bushman, F. D., et al. (2016). Composition and dynamics of the respiratory tract microbiome in intubated patients. Microbiome 4:7. doi: 10.1186/s40168-016-0151-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Khoury, M., Little, J., and Burke, W. (2003). “ACCE: a model process for evaluating data on emerging genetic tests,” in Human Genome Epidemiology: A Scientific Foundation for Using Genetic Information to Improve Health and Prevent Disease, eds M. Khoury, J. Little, and W. Burke (Oxford: Oxford University Press), 217–233.

PubMed Abstract | Google Scholar

Kitsios, G. D. (2018). Translating advancements in lung microbiome analytics into the next-generation diagnosticgold-standard for pneumonia: a physician-investigator’s perspective. mSystems 13:e00153-17.

Kitsios, G. D., and McVerry, B. J. (2018). Host-microbiome interactions in the subglottic space: bacteria ante Portas! Am. J. Respir. Crit. Care Med. doi: 10.1164/rccm.201802-0276ED [Epub ahead of print].

PubMed Abstract | CrossRef Full Text

Kitsios, G. D., Morowitz, M. J., Dickson, R. P., Huffnagle, G. B., McVerry, B. J., and Morris, A. (2017). Dysbiosis in the intensive care unit: microbiome science coming to the bedside. J. Crit. Care 38, 84–91. doi: 10.1016/j.jcrc.2016.09.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuczynski, J., Lauber, C. L., Walters, W. A., Parfrey, L. W., Clemente, J. C., Gevers, D., et al. (2011). Experimental and analytical tools for studying the human microbiome. Nat. Rev. Genet. 13, 47–58. doi: 10.1038/nrg3129

PubMed Abstract | CrossRef Full Text | Google Scholar

Laxminarayan, R., Duse, A., Wattal, C., Zaidi, A. K. M., Wertheim, H. F. L., Sumpradit, N., et al. (2013). Antibiotic resistance-the need for global solutions. Lancet Infect. Dis. 13, 1057–1098. doi: 10.1016/S1473-3099(13)70318-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, C. M., Aziz, M., Kachur, S., Hsueh, P.-R., Huang, Y.-T., Keim, P., et al. (2012). BactQuant: an enhanced broad-coverage bacterial quantitative real-time PCR assay. BMC Microbiol. 12:56. doi: 10.1186/1471-2180-12-56

PubMed Abstract | CrossRef Full Text | Google Scholar

Mandell, L. A., Wunderink, R. G., Anzueto, A., Bartlett, J. G., Campbell, G. D., Dean, N. C., et al. (2007). Infectious diseases society of America/American thoracic society consensus guidelines on the management of community-acquired pneumonia in adults. Clin. Infect. Dis. 44(Suppl. 2), S27–S72. doi: 10.1086/511159

PubMed Abstract | CrossRef Full Text | Google Scholar

McKay, H. S., Margolick, J. B., Martínez-Maza, O., Lopez, J., Phair, J., Rappocciolo, G., et al. (2017). Multiplex assay reliability and long-term intra-individual variation of serologic inflammatory biomarkers. Cytokine 90, 185–192. doi: 10.1016/j.cyto.2016.09.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Mehta, R. L., Kellum, J. A., Shah, S. V., Molitoris, B. A., Ronco, C., Warnock, D. G., et al. (2007). Acute kidney injury network: report of an initiative to improve outcomes in acute kidney injury. Crit. Care 11:R31. doi: 10.1186/cc5713

PubMed Abstract | CrossRef Full Text | Google Scholar

Modi, S. R., Collins, J. J., and Relman, D. A. (2014). Antibiotics and the gut microbiota. J. Clin. Invest. 124, 4212–4218. doi: 10.1172/JCI72333

PubMed Abstract | CrossRef Full Text | Google Scholar

Morris, A., Beck, J. M., Schloss, P. D., Campbell, T. B., Crothers, K., Curtis, J. L., et al. (2013). Comparison of the respiratory microbiome in healthy nonsmokers and smokers. Am. J. Respir. Crit. Care Med. 187, 1067–1075. doi: 10.1164/rccm.201210-1913OC

PubMed Abstract | CrossRef Full Text | Google Scholar

Naccache, S. N., Federman, S., Veeraraghavan, N., Zaharia, M., Lee, D., Samayoa, E., et al. (2014). A cloud-compatible bioinformatics pipeline for ultrarapid pathogen identification from next-generation sequencing of clinical samples. Genome Res. 24, 1180–1192. doi: 10.1101/gr.171934.113

PubMed Abstract | CrossRef Full Text | Google Scholar

Narvaez-Rivera, R. M., Rendon, A., Salinas-Carmona, M. C., and Rosas-Taraco, A. G. (2012). Soluble RAGE as a severity marker in community acquired pneumonia associated sepsis. BMC Infect. Dis. 12:15. doi: 10.1186/1471-2334-12-15

PubMed Abstract | CrossRef Full Text | Google Scholar

Neto, A. S., Barbas, C. S. V., Simonis, F. D., Artigas-Raventós, A., Canet, J., Determann, R. M., et al. (2016). Epidemiological characteristics, practice of ventilation, and clinical outcome in patients at risk of acute respiratory distress syndrome in intensive care units from 16 countries (PRoVENT): an international, multicentre, prospective study. Lancet Respir. Med. 4, 882–893. doi: 10.1016/S2213-2600(16)30305-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Panzer, A. R., Lynch, S. V., Langelier, C., Christie, J. D., McCauley, K., Nelson, M., et al. (2017). Lung microbiota is related to smoking status and to development of ARDS in critically ill trauma patients. Am. J. Respir. Crit. Care Med. 197, 621–631. doi: 10.1164/rccm.201702-0441OC

PubMed Abstract | CrossRef Full Text

Pendleton, K. M., Erb-Downward, J. R., Bao, Y., Branton, W. R., Falkowski, N. R., Newton, D. W., et al. (2017). Rapid pathogen identification in bacterial pneumonia using real-time metagenomics. Am. J. Respir. Crit. Care Med. 196, 1610–1612. doi: 10.1164/rccm.201703-0537LE

PubMed Abstract | CrossRef Full Text | Google Scholar

Price, R., MacLennan, G., Glen, J., and SuDDICU Collaboration. (2014). Selective digestive or oropharyngeal decontamination and topical oropharyngeal chlorhexidine for prevention of death in general intensive care: systematic review and network meta-analysis. BMJ 348:g2197. doi: 10.1136/bmj.g2197

PubMed Abstract | CrossRef Full Text | Google Scholar

R Foundation for Statistical Computing (2016). R: A Language and Environment for Statistical Computing. Vienna: CRAN.

Google Scholar

Raghu, V. K., Ramsey, J. D., Morris, A., Manatakis, D. V., Sprites, P., Chrysanthis, P. K., et al. (2017). Comparison of Strategies for Scalable Causal Discovery of Latent Variable Models from Mixed Data. Halifax, NS: Narcotics Anonymous.

Google Scholar

Rutjes, A. W. S., Reitsma, J. B., Coomarasamy, A., Khan, K. S., and Bossuyt, P. M. M. (2007). Evaluation of diagnostic tests when there is no gold standard. A review of methods. Health Technol. Assess. 11, iii, ix–51. doi: 10.3310/hta11500

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmidt, K., Mwaigwisya, S., Crossman, L. C., Doumith, M., Munroe, D., Pires, C., et al. (2017). Identification of bacterial pathogens and antimicrobial resistance directly from clinical urines by nanopore-based metagenomic sequencing. J. Antimicrob. Chemother. 72, 104–114. doi: 10.1093/jac/dkw397

PubMed Abstract | CrossRef Full Text | Google Scholar

Sedgewick, A. J., Shi, I., Donovan, R. M., and Benos, P. V. (2016). Learning mixed graphical models with separate sparsity parameters and stability-based model selection. BMC Bioinformatics 17:175. doi: 10.1186/s12859-016-1039-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Segal, L. N., Clemente, J. C., Tsay, J.-C. J., Koralov, S. B., Keller, B. C., Wu, B. G., et al. (2016). Enrichment of the lung microbiome with oral taxa is associated with lung inflammation of a Th17 phenotype. Nat. Microbiol. 1:16031. doi: 10.1038/nmicrobiol.2016.31

PubMed Abstract | CrossRef Full Text | Google Scholar

Singer, M., Deutschman, C. S., Seymour, C. W., Shankar-Hari, M., Annane, D., Bauer, M., et al. (2016). The third international consensus definitions for sepsis and septic shock (Sepsis-3). JAMA 315, 801–810. doi: 10.1001/jama.2016.0287

PubMed Abstract | CrossRef Full Text | Google Scholar

Srinivasan, R., Karaoz, U., Volegova, M., MacKichan, J., Kato-Maeda, M., Miller, S., et al. (2015). Use of 16S rRNA gene for identification of a broad range of clinically relevant bacterial pathogens. PLoS One 10:e0117617. doi: 10.1371/journal.pone.0117617

PubMed Abstract | CrossRef Full Text | Google Scholar

Sze, M. A., Hogg, J. C., and Sin, D. D. (2014). Bacterial microbiome of lungs in COPD. Int. J. Chron. Obstruct. Pulmon. Dis. 9, 229–238. doi: 10.2147/COPD.S38932

PubMed Abstract | CrossRef Full Text | Google Scholar

Valley, T. S., Sjoding, M. W., Ryan, A. M., Iwashyna, T. J., and Cooke, C. R. (2015). Association of intensive care unit admission with mortality among older patients with pneumonia. JAMA 314, 1272–1279. doi: 10.1001/jama.2015.11068

PubMed Abstract | CrossRef Full Text | Google Scholar

van der Sluijs, K. F., van der Poll, T., Lutter, R., Juffermans, N. P., and Schultz, M. J. (2010). Bench-to-bedside review: bacterial pneumonia with influenza - pathogenesis and clinical implications. Crit. Care 14:219. doi: 10.1186/cc8893

PubMed Abstract | CrossRef Full Text | Google Scholar

Venkataraman, A., Bassis, C. M., Beck, J. M., Young, V. B., Curtis, J. L., Huffnagle, G. B., et al. (2015). Application of a neutral community model to assess structuring of the human lung microbiome. MBio 6, e2284–14. doi: 10.1128/mBio.02284-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Vincent, J.-L., Brealey, D., Libert, N., Abidi, N. E., O’Dwyer, M., Zacharowski, K., et al. (2015). Rapid diagnosis of infection in the critically Ill, a multicenter study of molecular detection in bloodstream infections, pneumonia, and sterile site infections. Crit. Care Med. 43, 2283–2291. doi: 10.1097/CCM.0000000000001249

PubMed Abstract | CrossRef Full Text | Google Scholar

Zumla, A., Tawfiq, J. A., Al-Enne, V. I., Kidd, M., Drosten, C., Breuer, J., et al. (2014). Rapid point of care diagnostic tests for viral and bacterial respiratory tract infections–needs, advances, and future prospects. Lancet Infect. Dis. 14, 1123–1135. doi: 10.1016/S1473-3099(14)70827-8

CrossRef Full Text | Google Scholar

Keywords: microbiome, pneumonia, 16S rRNA gene sequencing, respiratory failure, antibiotic stewardship

Citation: Kitsios GD, Fitch A, Manatakis DV, Rapport SF, Li K, Qin S, Huwe J, Zhang Y, Doi Y, Evankovich J, Bain W, Lee JS, Methé B, Benos PV, Morris A and McVerry BJ (2018) Respiratory Microbiome Profiling for Etiologic Diagnosis of Pneumonia in Mechanically Ventilated Patients. Front. Microbiol. 9:1413. doi: 10.3389/fmicb.2018.01413

Received: 20 February 2018; Accepted: 08 June 2018;
Published: 10 July 2018.

Edited by:

Leonard Peruski, Centers for Disease Control and Prevention (CDC), United States

Reviewed by:

Mirjam Kooistra-Smid, Certe and University of Groningen, Netherlands
Brian Gregory George Oliver, University of Technology Sydney, Australia

Copyright © 2018 Kitsios, Fitch, Manatakis, Rapport, Li, Qin, Huwe, Zhang, Doi, Evankovich, Bain, Lee, Methé, Benos, Morris and McVerry. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Georgios D. Kitsios, kitsiosg@upmc.edu

Co-senior authors