1 Introduction
Community-acquired pneumonia (CAP) is the leading reason of hospitalization worldwide, which claims thousands of lives each year [
1]. Severe community-acquired pneumonia (SCAP) refers to patients complicated with severe hypoxemia, sepsis/septic shock, and acute respiratory distress syndrome (ARDS), who often require intensive care and mechanical ventilation [
2,
3]. Although broad-spectrum antibiotics and advanced oxygen support are used commonly as the therapeutic regimen, the mortality among SCAP patients remains as high as 25%–50% [
4,
5].
Next-generation sequencing has revealed that the human lung contains complex and dynamic microbiota community [
6]. The composition of lung microbiota in SCAP patients is complex and susceptible to changes during disease progression. The overgrowth of the invasive pneumonia-associated pathogens in patients’ lower respiratory tract may cause a decline in lung microbiota diversity, and the microbiota might even be dominated by single species [
7]. Based on the new conceptual models of respiratory microbiology proposed by Robert P Dickson, the lung microbiota is determined by three factors, namely, microbial immigration, elimination, and the relative reproduction rates of the members, which are affected by the regional growth conditions such as oxygen tension, pH, alveolar ventilation, and temperature in the lung [
8,
9]. Critical illness alters the internal environment and pathophysiology of the respiratory tract in patients with SCAP. Laryngeal dysfunction, supine positioning, aspiration or confusion, and intensive care unit (ICU) management, including antibiotic use, mechanical ventilation, and vasopressors, affect the balance of the immigration and elimination of airway microorganisms [
10]. Collectively, these clinical parameters have the potential to change the lung microbial composition of SCAP patients.
The relationship between respiratory microbial composition and respiratory infections has attracted increasing research and clinical interests. Recently, several studies have suggested that the lung microbiota is associated with respiratory infection and disease outcomes. The murine model suggests that the intact microbiota contributes to the protection against bacterial pathogens with GM-CSF signaling [
11]. The sputum microbial composition is associated with the length of stay and ICU admission in children hospitalized for CAP [
12]. Moreover, the enrichment of specific taxon in the airway microbiota, such as Enterobacteriaceae, might be related to high plasma inflammatory cytokine level and the development of ARDS, which indicate fewer ventilator-free days in patients with critical illness [
13,
14]. Although previous studies have identified that factors such as elder age, underlying diseases, and immunocompromised status predicted poor prognosis of SCAP patients [
15], the prognostic role of lung microbiota in SCAP patients remains unclear.
As few studies have described the key features and prognostic roles of the lung microbiota in patients with SCAP, we conducted a study to explore the clinical factors that may be associated with the lung microbiota in patients with SCAP and identify the important taxa that may predict clinical prognosis.
2 Materials and methods
2.1 Study populations
All patients with SCAP admitted in ICU between March 2018 and March 2019 were prospectively enrolled in this observational study. The inclusion criteria were as follows: (1) age≥18 years; (2) diagnosed with SCAP according to the 2007 Infectious Disease Society of America/American Thoracic Society guidelines (Table S1) [
16]; (3) time from illness onset to ICU admission≤7 days; for patients transferred to ICU from normal wards in our hospital because of disease exacerbation within 7 days of hospital admission, the time from illness onset to ICU admission can be extended to 14 days (Fig. 1). Patients were excluded if they met one of the following criteria: (1) having a history of hospitalization within 14 days before illness onset; (2) bronchoscopy couldn’t be performed within 48 h after admission to ICU; (3) peripheral blood specimens were not available within 24 h; (4) pregnant or breastfeeding women; (5) had an alternative diagnosis at the end of the study, including lung cancer, pulmonary tuberculosis, or pulmonary embolism.
2.2 Sample collection
Bronchoscopy was performed at bedside within 48 h of ICU admission using the standard clinical protocol [
17]. The bronchoscope was inserted through the nose or orotracheal tube, and the bronchoalveolar lavage fluid (BALF) samples were collected. Approximately 15 mL of BALF specimen was immediately sent to the microbiology laboratory for routine bacterial, fungal, and viral examinations, and the remaining BALF was stored at −80 °C until further processing. Blood samples were obtained within 24 h of ICU admission. The blood samples were centrifuged, and the remaining plasma was stored at −80 °C.
2.3 16S rRNA gene sequencing
DNA was extracted from BALF and negative control using the Maxwell® RSC Whole Blood DNA Kit (Promega, USA). The V3–V4 hypervariable region of the 16S rRNA gene of all samples was amplified by polymerase chain reaction (PCR). The pooled library was sequenced on the Illumina Miseq platform (Illumina, San Diego, CA, USA) using pair-end sequencing (2 × 300 bp). Details are shown in the supplementary material.
2.4 Data collection
The following data of included patients were collected using a standard case report form: demographic data, underlying diseases, the time of illness onset, clinical symptoms, laboratory findings, microbiology results, radiographic data, antimicrobial use, glucocorticoid use, and mechanical ventilation use. The included patients were followed up until discharged or until the patient died, whichever happened first.
2.5 Cytokine measurement
The plasma interleukin-4 (IL-4), IL-6, IL-8, macrophage inflammatory protein-1β (MIP-1β), vascular endothelial growth factor-A (VEGF-A), and matrix metalloproteinase-9 (MMP-9) of the included SCAP patients were detected using a magnetic bead-based multiplex immunoassay and read on a Bio-Plex 200 suspension array system (Bio-Rad, Hercules, CA, USA) according to the manufacturer’s instructions.
2.6 Definitions and outcomes
The primary outcome was the rate of clinical improvements, which was defined as a decrease of 2 categories and above on a 7-category ordinal scale within 14 days following bronchoscopy. The ordinal scales, which were also used in our influenza and COVID-19 studies, consisted of the following categories: (1) not hospitalized with resumption of normal activities; (2) not hospitalized, but unable to resume normal activities; (3) hospitalization, not requiring supplemental oxygen; (4) hospitalization, requiring supplemental oxygen; (5) hospitalization, requiring nasal high-flow oxygen therapy and/or noninvasive mechanical ventilation; (6) hospitalization, requiring extracorporeal membrane oxygenation and/or invasive mechanical ventilation; (7) death [
18–
20]. In our study, immunocompromised status was defined as meeting any of the following conditions: having a history of cancer with neutropenia (neutrophil count<0.5 × 10
9/L), hematological malignancies, solid malignancies receiving chemotherapy during the previous 3 months, solid organ or bone-marrow transplant, active graft-versus-host disease, bronchiolitis obliterans, human immunodeficiency virus infection, immunoglobulin deficiency, using immunosuppressive agents, or current treatment with systemic corticosteroids (≥20 mg of prednisone per day or its equivalent) for>30 continuous days before illness onset [
21]. The presence of ARDS on admission was diagnosed according to the Berlin definition [
22]. For nonventilated subjects who had a history of acute respiratory failure within 7 days because of a known respiratory event and bilateral pulmonary infiltration on chest X-ray with PaO
2/FiO
2 below 300 mmHg, ARDS was also diagnosed. Sepsis and septic shock were diagnosed on the basis of the third international consensus [
23]. Acute kidney injury (AKI) was diagnosed on the basis of the Kidney Disease Improving Global Guidelines (KDIGO) [
24]. Acute cardiac failure was diagnosed by cardiologists on the basis of clinical vitals, laboratory findings, and echocardiography. Pneumonia severity was assessed by CURB-65, APACHE-II, and PSI risk class. For patients admitted to ICU on the same day, the order of performing bronchoscopy was arranged according to the time of signing the bronchoscopy informed consent.
2.7 Microbial laboratory tests
All the specimens for diagnostic microbiology were collected within 48 h after ICU admission. Virus laboratory tests were considered positive if the respiratory virus was detected in sputum, endotracheal aspirates (ETA), BALF, or nasopharyngeal (NP) swabs by real-time PCR (Zhijiang, Shanghai, China), including respiratory syncytial virus, influenza virus types A and B, parainfluenza virus, rhinovirus, coronavirus, human metapneumovirus, and adenovirus. Positive bacteria or atypical pathogen laboratory test results were defined as meeting one of the following criteria: (1) positive bacterial culture in the blood or pleural fluid sample; (2) positive antigen test for
Legionella pneumophila (Binax Now; Trinity Biotech, Bray, Ireland) or
Streptococcus pneumoniae (Binax Now; Emergo Europe, Amsterdam, The Netherlands) in urine sample; (3) identification of
Mycoplasma pneumoniae or
Chlamydia pneumoniae in sputum, BALF, ETA, or NP swabs by real-time PCR (Zhijiang, Shanghai, China); (4) identification of
L. pneumophila in sputum, BALF, or ETA by real-time PCR (Zhijiang, Shanghai, China); (5) bacteria with moderate-to-heavy growth (>3+ growth) in qualified sputum or ETA or quantified culture in BALF of≥10
4 CFU/mL [
21]. The diagnosis of invasive fungal disease (IFD) was based on the revision and update of the consensus definitions of IFD [
25].
2.8 Statistical analysis
The software VSEARCH (version 2.7.1) and USEARCH (version 10.0) were used to process the sequencing data. Reads were denoised into Zero-radius Operational Taxonomic Units (ZOTUs) with UNOISE3. Contaminants and ZOTUs were identified by decontam package or comparison with the control samples (Table S2). A total of 490 ZOTUs were analyzed after removal of contaminants and ZOTUs whose relative abundance was less than 0.01%. Statistical analysis was performed in R version 3.6.2 via the Rstudio interface. PERMANOVA (vegan R-package) based on the Bray–Curtis distance was performed to assess the association between the clinical factors and the lung microbiota. A random forest learning approach (randomForest R-package) was used to identify the clinical factors associated with taxon. The Wilcoxon rank sum test and generalized linear models were performed to compare the relative abundance of the species. Multivariable Cox regression was performed to assess the association between lung microbiota and clinical improvements. The Cox regression model was adjusted for sampling season, plasma IL-8 level, CURB-65, APACHEII scores, presence of shock at sampling, oxygen index on admission, creatinine level, and microbiology results. Detailed statistical analysis was described in the supplementary material.
2.9 Data access
The sequence data have been deposited in the NCBI Sequence Read Archive under accession PRJNA616057. Feature, taxonomy, metadata tables, and a reproducible workflow of the analysis are available for download at related website.
3 Results
3.1 Patients’ characteristics and bacterial profiles
A total of 67 patients were included (Table 1), and the study flowchart is depicted in Fig. 1. The median age of study patients was 64 years. Fifty (74.63%) patients in the study population were male. Sixteen (23.88%) patients were current smokers, and 10 patients (14.93%) were in immunocompromised status. The sampling season was winter (27/67, 40.3%). Twenty-four (35.82%) patients had a decrease of 2 categories and above on a 7-category ordinal scale within 14 days, and the average number of days was 9.54 days. Forty (59.7%) patients had a decrease of 1 category and above on a 7-category ordinal scale within 14 days, and the average number of days was 6.8 days. The distribution of patients in each changed category ordinal from day 1 to 14 is shown in Fig. S1. A total of 67 BALF samples from the patients were analyzed. Proteobacteria was the most abundant phylum in the lung microbiota, accounting for 43.7% of the total filtered reads, followed by Firmicutes (30.94%), Actinobacteria (7.93%), Bacteroidetes (4.6%), and Fusobacteria (1.17%). At the family level, the most abundant taxa were Pseudomonadaceae (23.22%), Streptococcaceae (13.71%), Pasteurellaceae (11.33%), Corynebacteriaceae (5.72%), Moraxellaceae (5.15%), Staphylococcaceae (4.94%), Bacillaceae_1 (4.5%), Veillonellaceae (2.98%), Prevotellaceae (2.11%), and Enterobacteriaceae (1.74%).
3.2 Clinical factors associated with lung microbiota on ICU admission
We first performed a univariate PERMANOVA on the association between clinical factors and lung microbiota. The plasma MIP-1β level (Bray–Curtis distance; R2 = 0.029; P = 0.03), procalcitonin (PCT) level (R2 = 0.029; P = 0.026), complicated with AKI (R2 = 0.033; P = 0.012), positive bacteria laboratory test results (R2 = 0.033; P = 0.013), and the time from admission to sampling (R2 = 0.049; P = 0.018; Fig. 2A–2E) were significantly associated with the lung microbiota based on the univariate PERMANOVA. By contrast, age, baseline chronic respiratory diseases, antibiotic use, complicated with ARDS, and septic shock (Fig. S2, Table S3) were not significantly associated with the lung microbiota.
We performed a multivariable analysis to assess the independent association of the clinical factors. Smoking status (R2 = 0.036; P = 0.185), time on mechanical ventilator before sampling (R2 = 0.04; P = 0.097), and sputum production (R2 = 0.021; P = 0.118; Fig. 2F) whose P value was less than 0.2 were also included. κ test did not show a severe multicollinearity problem among the eight factors (K= 4.8; Fig. 2F). The positive bacteria laboratory test results showed the strongest independent association with the lung microbiota composition (multivariable PERMANOVA, R2 = 0.033; P = 0.018), followed by AKI (R2 = 0.032; P = 0.011) and plasma MIP-1β level (R2 = 0.027; P = 0.044).
3.3 Bacteria laboratory test results and lung microbial taxa
The Prevotellaceae family played the most important role in distinguishing bacterial samples with positive bacteria laboratory test results from those with negative results, followed by the pathogenic taxa, such as, the families Moraxellaceae, Staphylococcaceae, and Streptococcaceae (Fig. 3A). In patients with positive bacteria laboratory test results, the mean relative abundance of the families Prevotellaceae, Microbacteriaceae, and Actinomycetaceae decreased, whereas the families Moraxellaceae, Staphylococcaceae, and Streptococcaceae increased (Fig. 3A and 3B), which were in line with their high isolation rate of typical bacterial culture (Table 2). 16S rRNA gene sequencing could produce the same results as culture methods in 61% of patients (Table 2). However, no significant difference in the relative abundance of the family Pseudomonadaceae was observed between the two groups of patients (Fig. 3A). The samples with bacteria detected had lower α diversity, and the difference was not significant (Fig. 3C).
3.4 Lung microbiota and clinical improvements
The lung microbial community composition of the patients on admission was not associated with clinical improvements (PERMANOVA; Bray–Curtis distance; R2 = 0.011; P = 0.76). The increased Shannon index predicted a faster decrease of 2 categories, which was analyzed continuously (adjusted HR 1.43, 95% CI 1.04–1.98, P = 0.03) or by organizing their value into ranges (adjusted HR 5.73, 95% CI 1.41–23.26, P = 0.01; Fig. 4B, Table 3). Although positive bacteria laboratory test results were not associated with clinical improvement, we found that every 1% increase in the relative abundance of the families Prevotellaceae and Actinomycetaceae in the lung microbiota of the patients increased the probability for clinical improvements by 14% and 10%, respectively (95% CI 1.04–1.25, P = 0.006; 95% CI 1.02–1.18, P = 0.01; Table 3, Fig. 4C and 4D). However, the families Moraxellaceae, Staphylococcaceae, and Streptococcaceae were not associated with clinical improvement (Table 3, Fig. 4C and 4D). The proportions of patients with ARDS, septic shock, CURB-65≥3, and PSI≥IV were similar irrespective of the different levels of α diversity or relative abundance of taxa (Table S4).
4 Discussion
This prospective cohort study described the composition of lung microbiota and the important taxa predicting clinical prognosis in patients with SCAP. The core findings of our study indicated that positive bacteria laboratory test results had the strongest independent association with the lung microbiota composition on patients’ ICU admission, followed by AKI and plasma MIP-1β level. The increased α diversity and the enrichment of the families Prevotellaceae and Actinomycetaceae in the lung microbiota were associated with clinical improvements.
In this study, the positive bacteria laboratory test results on admission were the most important factor associated with patients’ lung microbiota. The bacteria that could be identified by the laboratory tests used in this study (primarily culture-based methods) might be more abundant in the lower respiratory tract [
26]. This finding might explain why these taxa (e.g.,
Acinetobacter baumannii and
Staphylococcus aureus) dominated the lung microbiota in patients with positive bacteria laboratory results as revealed by 16S rRNA gene sequencing. Although the results of 16S rRNA sequencing and culture-based methods were consistent with regard to dominant bacteria, the results related to
Pseudomonas aeruginosa differed in these two technologies. Although culture-based methods had high isolation rate (33.33%) of
Pseudomonas aeruginosa, 16S rRNA sequencing did not show a significant difference in the relative abundance of the family Pseudomonadaceae between the positive bacteria laboratory test and negative bacteria laboratory test groups. The ICU environment might be the reason behind the phenomenon.
P. aeruginosa is a primary nosocomial pathogen in ICU. The risk factors for
P. aeruginosa acquisition consisted of antibiotic pressure, mechanical ventilation use, length of hospitalization, cross-transmission among patients and medical staffs, and environmental factors, such as contamination of tap water [
27–
29]. All of our patients had received antimicrobial treatments, and more than half were exposed to mechanical ventilation before bronchoscopy. Thus, these patients might have the same risk of obtaining
P. aeruginosa colonization under the ICU status. Although 16S rRNA gene sequencing had a low accuracy for pathogen identification, it could capture all the information related to bacterial species, including fastidious bacteria in the respiratory tract, which might explain the different results between culture-dependent and culture-independent analysis in our study [
30,
31].
AKI was common in SCAP patients [
32]. Many studies had suggested a complex cross-talk between the lung and the kidney [
33]. Cigarette smoking could induce endothelial cell injury with increased AGEs-RAGE levels in the lung and kidney [
34]. In critically ill patients, mechanical ventilation was proven to be a risk factor for AKI by a meta-analysis [
34]. A study in mice with sepsis established by cecal ligation and puncture showed that ventilator-induced lung injury increased the renal expression of VEGF, VCAM-1, and angiopoietin-2 [
35]. In our study, we observed that the presence of AKI on admission was associated with lung microbiota of the patients. However, few studies described the effect of AKI on lung microbiota. The plasma MIP-1β level was significantly associated with the abundance of the Streptococcaceae family.
In vitro experiments confirmed that
S. pneumoniae D39 could induce the upregulation of MIP-1β gene in human monocytes [
36]. Similar conclusions, that is, the increased level of MIP-1β in BALF was related to the enrichment of
S.
pneumoniae in the lung microbiota, were drawn from patients with asthma [
37]. The enrichment of
Streptococcus spp. in the lung microbiota was associated with high plasma MIP-1β level in our patients. The ZOTU4, belonging to the family Streptococcaceae, was significantly higher in the high MIP-1β level group (Fig. S3).
In our study, after adjustment for microbiology results and severity of illness, multivariate analysis showed that clinical improvements were associated with increased relative abundance of the family Prevotellaceae and Actinomycetaceae. The results of our study were supported by several previous studies. Based on previous studies, the relative abundance of
Prevotella was inversely associated with the extent of airway inflammation in patients with cystic fibrosis [
38]. An increase in the abundance of
Prevotella in nose/throat microbiota was related to lower susceptibility to influenza A infection [
39]. As the families Prevotellaceae and Actinomycetaceae usually resided in the healthy respiratory tract or oral cavity [
40–
42], the protective role of these families might stem from its ability to inhibit the growth of pathogenic bacteria. For example, coculture experiments observed that
Prevotella spp. could inhibit
Hemophilus influenzae-induced IL-12p70 in dendritic cells [
43]. Nevertheless, another set of research showed that the enrichment of lung microbiota with
Prevotella enhanced the level of BALF inflammatory cytokines, which was associated with the development of asthma in children and ARDS in severe patients [
13,
40,
44]. The findings in current studies appeared to be inconsistent probably because of the observational nature and potential confounders in these studies. In addition, the disease type and population in these studies were different, and the result of current studies had limited generalizability. Therefore, the actual role of these taxa in disease severity and prognosis of patients with respiratory diseases should be further explored.
Different from some respiratory microbiome studies in critically ill patients, our study found that the composition of lung microbiota at sampling was not associated with the presence of ARDS or septic shock. For example, unlike the report in two previous published studies [
13,
14], the relative abundance of the family Enterobacteriaceae was not associated with disease severity and clinical outcomes in our study probably because these studies primarily focused on patients with ARDS caused by severe trauma, cardiac arrest, or cerebral vascular events but not severe infection. In Robert P. Dickson
et al.’s study, patients with pneumonia or sepsis merely accounted for 13% (12/91), and in Ariane R. Panzer’s study, the ARDS in patients was due to trauma [
13,
14]. Thus, the conclusions of the studies might not be applied to SCAP patients.
Our study has several limitations. First, the lung microbiota analysis was only performed on ICU admission. The variation of lung microbiota during patients’ hospitalization and the association of lung microbiota with clinical improvements could not be analyzed. Second, the sample size of the included patients was small compared with studies for other diseases. Further studies are needed to validate the conclusions.
5 Conclusions
In this study on respiratory microbiome of SCAP patients, we found that the positive bacteria laboratory test results, AKI, and plasma MIP-1β level were associated with the composition of lung microbiota on patients’ ICU admission. The increased α diversity and the enrichment of the families Prevotellaceae and Actinomycetaceae in the lung microbiota were associated with clinical improvements.