Skip to main content

Respiratory microbiome profiles differ by recent hospitalization and nursing home residence in patients on mechanical ventilation



Healthcare-associated pneumonia (HCAP) is a heterogeneous disease. We redefined nursing-home- and hospital-associated infections (NHAI) group by revising existing HCAP risk factors. The NHAI group comprised nursing home residents with a poor functional status, or recent (past 90 days) hospitalization or recent (past 180 days) antibiotic therapy. Our aim was to determine whether respiratory microbiota profiles are related to newly defined NHAI group in critically ill patients on mechanical ventilation.


The 180 endotracheal aspirates (ETAs) from 60 mechanically ventilated ICU patients (NHAI group, n = 24; non-NHAI group, n = 36) were prospectively collected on days 1, 3 and 7 in a university hospital. The bacterial community profiles of the ETAs were explored by 16S rRNA gene sequencing. A phylogenetic-tree-based microbiome association test (TMAT), generalized linear mixed models (GLMMs), the Wilcoxon test and the reference frame method were used to analyze the association between microbiome abundance and disease phenotype.


The relative abundance of the genus Corynebacterium was significantly higher in the pneumonia than in the non-pneumonia group. The microbiome analysis revealed significantly lower α-diversity in the NHAI group than in the non-NHAI group. In the analysis of β-diversity, the structure of the microbiome also differed significantly between the two groups (weighted UniFrac distance, Adonis, p < 0.001). The abundance of Corynebacterium was significantly higher, and the relative abundances of Granulicatella, Staphylococcus, Streptococcus and Veillonella were significantly lower, in the NHAI group than in the non-NHAI group.


The microbiota signature of the ETAs distinguished between patients with and without risk factors for NHAI. The lung microbiome may serve as a therapeutic target for NHAI group.


Among the recommendations of the 2016 and 2019 American Thoracic Society/Infectious Diseases Society of America (ATS/IDSA) guidelines was the removal of the category of healthcare-associated pneumonia (HCAP) [1, 2]. This is in contrast to the 2005 guidelines, which advocated the use of broad-spectrum antibiotic therapy in HCAP and nosocomial pneumonia [3]. The change was based on several studies showing that, while patient characteristics such as comorbidities and age are determinants of the risk of multidrug-resistant (MDR) infection, this is not the case for all of the HCAP criteria [4,5,6]. Jones et al. showed that acceptance of the HCAP criteria lowered the threshold for broad-spectrum antibiotics for infection with nosocomial pathogens [7]. However, some HCAP criteria, such as admission from a nursing home or recent (preceding 90 days) hospitalization, are independent predictors of MDR infection [5, 8, 9]. For example, the risk of MDR infections in patients with HCAP requiring mechanical ventilation is similarly high to that of nosocomial infections [10, 11].

Based on previous studies, three risk factors may be better than the previous HCAP criteria for predicting the risk of infection with nosocomial pathogens: recent hospitalization, nursing home residence and antibiotic use in the past 90 days [5, 8, 12]. The application of new, culture-independent methods of characterizing microbial communities, such as 16S rRNA, may lead to better characterization of patient risk groups and more appropriate therapeutic strategies [13]. We hypothesized that the composition of the bacterial microbiome in the respiratory tract during mechanical ventilation would be related to the newly defined risk factors for nursing-home and hospital-associated infections (NHAI), and to the 28-day all-cause and final hospital mortality rates. Therefore, we analyzed the microbial profiles of the endotracheal tube aspirates (ETAs) of intubated patients and determined the relationship between the ETA microbiome and risk factors for NCAI.


Patients and study design

Patients admitted to the intensive care unit (ICU) of Chuncheon Sacred Hospital (South Korea) between July 2017 and October 2018 were prospectively recruited. All patients placed on mechanical ventilation at the time of ICU admission were included in the study. The exclusion criteria were age < 18 years, initiation of mechanical ventilation > 48 h after ICU admission and a duration of mechanical ventilation < 7 days.

In 2005, the ATS/IDSA guidelines introduced HCAP as a new category of pneumonia requiring therapy similar to that prescribed for nosocomial pneumonia, namely broad-spectrum antibiotics [3]. However, several studies showed that HCAP is a heterogeneous disease, and that broad-spectrum antibiotic therapy may not be necessary for all patients with HCAP [6, 8, 14]. Therefore, attempts were made to refine the definition of HCAP based on risk factors for MDR infections. On the basis of the available data [5, 8, 10, 12, 14, 15], we redefined NHAI patients as those meeting at least one of the following criteria: (1) nursing home resident with poor functional status; (2) recent (past 90 days) hospitalization; or (3) recent (past 180 days) antibiotic therapy.

The 60 patients (41 with and 19 without pneumonia) enrolled in this study were divided into NHAI and non-NHAI groups and prospectively followed. This study was approved by the Institutional Review Board of Chuncheon Sacred Heart Hospital (IRB approval number: 2017‐47).

Data collection and clinical outcome measures

Data on demographic characteristics (age, sex, preexisting comorbidities), the indications for intubation, the presence of acute respiratory distress syndrome (ARDS), the arterial oxygen tension/fraction of inspired oxygen (SaO2/FiO2) ratio, and the Glasgow Coma Scale (GCS), Acute Physiology and Chronic Health Evaluation (APACHE) II) and Sequential Organ Failure Assessment (SOFA) scores were recorded. The Charlson Comorbidity Index (CCI) score was calculated as described previously [16]. The clinical outcomes of interest in this study were the 28-day all-cause and final hospital mortality rates. Respiratory samples were acquired by endotracheal aspiration, with serial samples collected 1, 3 and 7 days after the initiation of mechanical ventilation.

DNA extraction, PCR and sequencing

Genomic DNA was extracted using a commercial microbial DNA isolation kit (MP Biomedicals, USA) and then amplified using primers targeting the V3–V4 region of the prokaryotic 16S rRNA gene. The amplicons were purified using XSEP MagBead (CELEMICS), and then subjected to PCR using the Nextera Index Kit (Illumina, USA) following the manufacturer’s instructions. The resultant product was further purified using XSEP MagBead (CELEMICS). The prepared bacterial amplicon library was then quantified, mixed with multiple libraries and sequenced using the MiSeq v3 platform (Illumina).

The raw sequencing data files were preprocessed before downstream data analysis. To remove sequences with low-quality scores, the raw sequence reads were pre-filtered using PRINSEQ [17]. Adapter sequences were removed using CUTADAPT [18]. Paired-end reads were merged using PEAR [19] and then filtered with PRINSEQ. Chimeric sequences and singletons were screened and reduced using USEARCH [20]. Finally, 180 samples were analyzed and an average of 133,219 reads per sample were obtained (minimum, 29,247; maximum, 587,337). The downstream data analysis was carried out using QIIME [21] with the EzBioCloud 16S rRNA gene sequence database [22]. Operational taxonomic units (OTUs) were defined as clusters of sequences with ≥ 97% identity.

Data processing and statistical analysis

Various Microbial community comparisons were performed using Quantitative Insights in Microbial Ecology software (QIIME, v. 1.8) and the R package Vegan. In the analysis, two period samples (day 1 and day 7) were selected to clarify the differences between groups and to confirm the changes over time. The microbial sequences from the ETAs represented 554 genera and 2009 species. Absolute OTU of each subject for both time points was combined, and species were filtered by two criteria; (1) eliminating strains that appear only in a few patients, (2) removing strains that account for less than 1000 OTU. Finally, 11 predominant genera (Acinetobacter, Streptococcus, Corynebacterium, Staphylococcus, Prevotella, Neisseria, Veillonella, Mycoplasma, Granulicatella, Actinomyces, Campylobacter) and 8 predominant species (Acinetobacter baumannii, Streptococcus mitis, Corynebacterium ulcerans, Staphylococcus caprae, Veillonell adispar, Granulicatella adiacens, Streptococcus parasanguinis, Streptococcus lactarius) were included in the statistical analysis.

Ace, Chao1, the Shannon index and the Simpson index were used to express α-diversity; β‐diversity (inter-sample diversity) was defined as the extent of the similarity between microbial communities based on the degree of structural overlap. Nonmetric multidimensional scaling plots based on weighted UniFrac distances were used to visualize the differences between groups in microbial community structure. Principle coordinate analysis was conducted using a permutational multivariate ANOVA (PERMANOVA), performed via the Adonis function of the R package vegan (1000 permutations).

We conducted statistical analysis to identify the association of pneumonia and NHAI risk factor on the abundance with each OTU. Each subject was measured at two different time points (day 1 and 7) and generalized linear mixed effects models (GLMM) was utilized to handle the repeatedly observed measurement. Pneumonia and NHAI risk factor were considered as response variables. Relative abundances of OTUs at each time point were Log2 transformed, and its effect was evaluated by including it as a covariate. We also considered the paired and longitudinal distance-based approach (Pldist) to detect the longitudinal changes in the microbiome over time with various outcome types [23]. It can also calculate differences in non-phylogenetic dissimilarities such as Gower’s distance, Bray–Curtis dissimilarities, Jaccard distance and Kulczynski distance between time points. Last we found that pneumonia and NHAI risk factor are time-invariant and thus considered four methods for cross-sectional data; phylogenetic tree-based microbiome association test (TMAT), OMiAT (version 5.1), Wilcoxon test and reference frame method [24, 25]. For the cross-sectional methods, means of log-transformed read count per million (CPM) at two different point (day 1 and 7) were considered as the response variables. Among all approaches, other tests except for reference frame method can adjust the effect of the other covariates and age, sex, APACHE score and CCI were included as covariates for both. GLMM was conducted with SAS 9.4 (Cary, NC), and the other analyses were done with R software. All test results were considered significant using a p value of ≤ 0.05.


Characteristics of the study participants

The study population consisted of 60 patients who had been placed on mechanical ventilation, including 41 patients diagnosed with pneumonia and 24 who had risk factors for NHAI. Table 1 shows the patient characteristics. The proportions of ARDS, a high C-reactive protein (CRP) level and an abnormal paO2/FiO2 ratio were higher in patients with than without pneumonia, as expected. Patients in the NHAI and non-NHAI groups did not differ in terms of sex, severity index (APACHE and SOFA) scores, or the GCS score. However, patients in the NHAI group were older and had more comorbidities than the non-NHAI group. There were no significant differences between the two groups in either the 28-day all-cause or final hospital mortality rate.

Table 1 Baseline characteristics of the study participants

Differences in the microbial composition of the ETAs between patient groups (pneumonia vs non-pneumonia, NHAI vs. non-NHAI) are shown in the taxonomic summary (Fig. 1). The average abundance of the genera Corynebacterium, Staphylococcus and Pseudomonas was higher in the pneumonia group than in the non-pneumonia group. The relative abundances of Streptococcus and Prevotella were higher in the non-NHAI group, whereas the relative abundance of Corynebacterium was higher in the NHAI group.

Fig. 1

Relative abundance of bacterial communities in the endotracheal aspirates of the study participants over time. Taxonomic summaries of the 180 samples are shown for pneumonia and non-pneumonia patients (a) and for NHAI and non-NHAI patients (b)

Respiratory microbiota during mechanical ventilation in patients with and without pneumonia

Neither α-diversity nor β-diversity differed significantly between patients with and without pneumonia (not shown). The association between pneumonia and microbial taxa was analyzed, adjusting for the effects of age, sex, and APACHE and CCI scores in the TMAT, optimal microbiome-based association test (OMiAT), and GLMM, but not in the Wilcoxon test. A correlation between Corynebacterium and pneumonia was demonstrated by the TMAT and Wilcoxon test (Table 2). Especially, the GLMMs identified a significant difference in the abundance of C. ulcerans between patients with and without pneumonia (FDR corrected p = 0.035), although it did not find the correlation between Corynebacterium and pneumonia at the genus level.

Table 2 Bacterial genera associated with pneumonia

Figure 2a shows the significantly increased log-transformed CPM of the genus Corynebacterium in the pneumonia than in the non-pneumonia group by TMAT with false discovery rate using Benjamini–Hochberg correction. The difference of this genus abundance was mainly caused by a species named C. ulcerans. No other Corynebacterium species detected in the samples showed meaningful difference between the pneumonia and non-pneumonia groups. Figure 2b shows the larger increase in C. ulcerans in the pneumonia than in the non-pneumonia group, both on day 1 and on day 7 (Day 1 p-value = 0.003, Day 7 p-value = 0.007). A trend towards a larger difference in the OTU abundance of Corynebacterium in the pneumonia group compared to the non-pneumonia group on day 7 was also identified.

Fig. 2

Increased abundance of Corynebacterium ulcerans in patients with pneumonia. a Box and whisker plot showing the relationship between Corynebacterium OTUs and pneumonia group. For the lines in a box and whisker plot: the extreme bars are the overall range, the bottom and top of the box are the 25th and 75th percentiles and the line inside the box is the 50th percentile (median). The significantly increased log-transformed CPM of the genus Corynebacterium in the pneumonia than in the non-pneumonia group was shown by TMAT with false discovery rate using Benjamini–Hochberg correction (FDR corrected p-value = 0.0226). b Longitudinal patterns of relative abundance of C. ulcerans OTUs in the pneumonia and non-pneumonia groups. Wilcoxon test was used (Day 1 p-value = 0.003, Day 7 p-value = 0.007). **p < 0.01

Respiratory microbiota during mechanical ventilation in patients with and without risk factors for NHAI

An analysis of α-diversity revealed significant differences between the NHAI and non-NHAI groups (expressed as median [interquartile range]; Fig. 3a). Ace (67.5 [48.2, 108.9] vs. 139.9 [106.2, 215.7]; p < 0.001), Chao (64.4 [45.3, 112.0] vs. 137.5 [106.2, 230.6]; p < 0.001), Shannon index (2.3 [0.7, 2.9] vs. 3.8 [2.3, 4.7]; p < 0.001) and Simpson index (0.7 [0.2, 0.8] vs. 0.9 [0.6, 0.9]; p < 0.001) values were lower in the NHAI group than in the non-NHAI group. The microbiome structure of the NHAI group could also be clearly distinguished from that of the non-NHAI group based on the β-diversity (PERMANOVA, adjusted p < 0.001, Fig. 3b).

Fig. 3

Differences in respiratory microbiome in NHAI group and non-NHAI group. a Comparison of the different metrics of α-diversity between the NHAI group (group 1) and non-NHAI group (group 0). For the lines in a box and whisker plot: the extreme bars are the overall range, the bottom and top of the box are the 25th and 75th percentiles and the line inside the box is the 50th percentile (median). b Microbial community structure in the endotracheal aspirates according to risk factors for NHAI based on weighted UniFrac distances. A PERMANOVA, performed using the Adonis function in the R package vegan) was conducted to compare the microbial community structure between the NHAI and non-NHAI groups. The x- and y-axis represent the first and second principal coordinates with the proportion of variance. The ellipses represent the 95% confidence interval for each group assuming a multivariate normal distribution. The analysis was adjusted for age, sex, APACHE score and Charlson Comorbidity Index score. ***p < 0.001 ****p < 0.0001

The microbial taxa associated with NHAI were identified using the methods shown in Table 3. Adjustments were made for age, sex, and the APACHE and CCI scores in the TMAT, OMiAT, and GLMMs, but not in the Wilcoxon test. The TMAT and Wilcoxon analyses revealed that Granulicatella, Streptococcus and Veillonella were negatively correlated with NHAI, while Corynebacterium was positively correlated (Table 3).

Table 3 Bacterial genera associated with NHAI

Similarly, the GLMMs revealed that the genera Granulicatella and Prevotella were more abundant in the non-NHAI group, while at the species level C. ulcerans was more abundant in the NHAI group and Granulicatella adiacens was more abundant in the non-NHAI group (Additional file 1: Table S1).

Additional file 2: Figure S1 shows the longitudinal trends in microbial taxa according to the risk factors for NHAI. In ETAs from days 1 and 7, the OTU abundances of Veillonella, Granulicatella and Streptococcus were significantly higher in the non-NHAI than in the NHAI group, whereas Corynebacterium was more abundant in the

NHAI group. However, none of the longitudinal changes in microbial taxa differed significantly by disease phenotype.

Longitudinal change in the respiratory microbiome according to clinical outcome

The associations between non-phylogenetic dissimilarities (based on pldist) and the clinical outcomes (the 28-day all-cause and final hospital mortality rates) were analyzed using the microbiome regression-based kernel association test (MiRKAT). The results showed that the variability in the respiratory microbiome was not associated with the clinical outcome (Table 4). The GLMM analysis yielded similar results. There was no association of microbial OTU abundance with the 28-day all-cause or final hospital mortality rate, in either the pneumonia or non-pneumonia group (Additional file 1: Table S2).

Table 4 Association between microbiome variability and clinical outcome


Patients with HCAP are a heterogeneous group [14]. Because the existing HCAP criteria did not reflect the risk of nosocomial pathogen infection and led to overuse of broad-spectrum antibiotics, they were omitted from the 2016 hospital-acquired pneumonia/ventilator-associated pneumonia (HAP/VAP) guidelines [1, 7]. Attempts to develop alternative criteria for HCAP, and a specific therapy [8, 11, 14], included our concept of NHAI, which may better predict the occurrence of MDR infections and consequences of continued exposure to hospital and other healthcare environments.

Recent studies using conventional microbiologic culture tests have shown different patterns of infection among countries and geographic regions [26,27,28,29]. Studies from the USA and the Asia/Pacific region reported that HCAP is characterized by high frequencies of infections with drug-resistant pathogens [26,27,28, 30], whereas other studies showed similarities in the involved pathogens between HCAP and community-acquired pneumonia (CAP) [29].

Culture-independent 16S rRNA gene sequencing is useful for detecting fastidious bacteria, the abundance of which tends to be underestimated; this method is more sensitive than conventional culture methods [31]. To our knowledge, this is the first study to use microbiome sequencing to determine the associations of healthcare/hospital-related factors with the composition of the respiratory microbiome of intubated patients.

The most important finding of this study was that the ETA microbiome profile of intubated patients was significantly associated with newly defined risk factors for NHAI. Generally, microbial diversity was higher in patients without than with risk factors for NHAI. The analysis of β-diversity, based on weighted UniFrac distances, showed distinct clusters of samples, identified as NHAI and non-NHAI groups. Beyond microbial diversity, our study showed the association of specific bacterial genera with continued exposure to healthcare and hospital environments.

The abundance of Corynebacterium increased, whereas the abundances of Granulicatella, Streptococcus, Staphylococcus and Veillonella decreased, in intubated patients with risk factors for NHAI. In particular, members of the genus Corynebacterium were significantly more abundant in intubated patients with than without pneumonia, and the abundance of C. ulcerans in the pneumonia group tended to increase over time. Although Corynebacterium spp. are among the commensal flora of the nasopharynx [32], a causative role in pneumonia was recently described [33,34,35,36]. The identification of Corynebacteria subspecies is now possible using new diagnostic technologies, such as matrix-assisted laser desorption ionization time of flight mass spectrometry and 16S rRNA gene sequencing [13, 33,34,35, 37]. Based on such studies, Clariot et al. suggested that Corynebacterium spp. may be responsible for pneumonia in mechanically ventilated patients, similar to our results [33], while Yasuda et al. reported that C. ulcerans produces diphtheria toxin and causes severe pneumonia complicated by diffuse pseudomembrane formation in the central airways [36]. Corynebacterium spp. may also be responsible for HAP. An analysis of bronchoalveolar lavage fluid (BALF) obtained from Japanese patients revealed Corynebacterium spp. in the samples of 11.8% of the patients with HAP, and in 4.9% of those from patients with HCAP; however, Corynebacterium spp. was present in only 1.6% of the samples from patients with CAP [37]. Our results also demonstrate the need to consider Corynebacterium spp. as a possible healthcare/hospital-associated pathogen.

Several studies using molecular techniques have revealed the presence of unexpected microorganisms in specific phenotypes of pneumonia. For example, Steenhuisen et al. found an increase in Lactobacillus and Rothia in CAP, while Zarkharkina et al. reported a decrease in these microorganisms, in patients with VAP [38, 39]. Similar to our own findings, a study of adult ICU patients identified Prevotella as being dominant in controls, whereas Corynebacteria species were uniquely present in VAP cohorts [40].

Segal et al. reported that a lung microbiome containing taxa from the oral cavity (i.e., Streptococcus, Veillonella, Granulicatella, and Prevotella) induced host cellular mucosal immunity of the Th17/neutrophilic phenotype and blunted TLR4 responses [41]. In the study of Kitsios et al., the dominance of oral taxa within lung communities was strongly associated with culture negativity and weaker host immune responses [42]. From these data, it can be inferred that non-NHAI patients in whom oral taxa dominate the lung microbiome will have a weaker immune response than NHAI patients, which would also explain why patients with HCAP present with more severe disease and have clinical outcomes resembling those of nosocomial pneumonia [26, 28, 43].

Our data showed no difference in microbial community according to clinical outcome; this is in contrast to Dickson et al., who reported that the lung microbiome predicts outcomes in critically ill patients [44]. This discordance may be due to differences in study population and sampling method (ETA vs. BALF).

Our study had several limitations. First, it used a single-center design and the sample size was small. Multicenter studies with larger populations are needed to validate our results. Second, we did not obtain BALF samples, although this was in accordance with ATS guidelines recommending non-invasive testing [1]. Nevertheless, analyzing BAL and ETA samples may enable comparison with the results of other studies, such that more robust conclusions could be drawn. Third, although our findings support the concept of an NHAI risk factor-associated microbiome, its clinical significance should be investigated via functional analyses, for example. Such studies would allow conclusions to be drawn regarding the association between an NHAI-specific microbiome and outcome/immune metabolic responses.

The strengths of our study included the homogeneity of the study population in terms of age (all elderly patients), and the multiple serial samples acquired from each patient. The analysis showed that the composition of the NHAI-related microbiome was similar regardless of antibiotic use or time of day.


In this prospective observational cohort study of mechanically ventilated patients, the loss of diversity and dysbiosis of the respiratory microbiome were more profound in patients with than without risk factors for NHAI, which were in turn positively associated with the presence of Corynebacterium, and negatively associated with that of Granulicatella, Streptococcus, Staphylococcus and Veillonella. Further studies are needed to determine the clinical relevance of microbial community profiling, and to identify novel therapeutic targets in critically ill elderly patients with risk factors for NHAI.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request. We submitted our raw sequence reads to NCBI. The SRA submission is in process and the BioProject accession number is PRJNA678854.



Healthcare-associated pneumonia


Nursing home and hospital associated infection


Intensive care unit


Tree-based microbiome association test


Generalized linear mixed effects models


Endotracheal tube aspirates


Operational taxonomic units


Community acquired pneumonia


Ventilator associated pneumonia


Multidrug resistant


Acute Physiology, Age, Chronic Health Evaluation II


Sequential Organ Failure Assessment


Charlson comorbidity index


  1. 1.

    Kalil AC, Metersky ML, Klompas M, Muscedere J, Sweeney DA, Palmer LB, Napolitano LM, O’Grady NP, Bartlett JG, Carratala J, et al. 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. 2016;63:e61–111.

    PubMed  PubMed Central  Article  Google Scholar 

  2. 2.

    Metlay JP, Waterer GW, Long AC, Anzueto A, Brozek J, Crothers K, Cooley LA, Dean NC, Fine MJ, Flanders SA, et al. Diagnosis and treatment of adults with community-acquired pneumonia. An official clinical practice guideline of the American Thoracic Society and Infectious Diseases Society of America. Am J Respir Crit Care Med. 2019;200:e45–67.

    PubMed  PubMed Central  Article  Google Scholar 

  3. 3.

    American Thoracic S, Infectious Diseases Society of A. Guidelines for the management of adults with hospital-acquired, ventilator-associated, and healthcare-associated pneumonia. Am J Respir Crit Care Med. 2005;171:388–416.

    Article  Google Scholar 

  4. 4.

    Chalmers JD, Rother C, Salih W, Ewig S. Healthcare-associated pneumonia does not accurately identify potentially resistant pathogens: a systematic review and meta-analysis. Clin Infect Dis. 2014;58:330–9.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  5. 5.

    Gross AE, Van Schooneveld TC, Olsen KM, Rupp ME, Bui TH, Forsung E, Kalil AC. Epidemiology and predictors of multidrug-resistant community-acquired and health care-associated pneumonia. Antimicrob Agents Chemother. 2014;58:5262–8.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  6. 6.

    Yap V, Datta D, Metersky ML. Is the present definition of health care-associated pneumonia the best way to define risk of infection with antibiotic-resistant pathogens? Infect Dis Clin N Am. 2013;27:1.

    Article  Google Scholar 

  7. 7.

    Jones BE, Jones MM, Huttner B, Stoddard G, Brown KA, Stevens VW, Greene T, Sauer B, Madaras-Kelly K, Rubin M, et al. Trends in antibiotic use and nosocomial pathogens in hospitalized veterans with pneumonia at 128 medical centers, 2006–2010. Clin Infect Dis. 2015;61:1403–10.

    PubMed  PubMed Central  Article  Google Scholar 

  8. 8.

    Park SC, Kang YA, Park BH, Kim EY, Park MS, Kim YS, Kim SK, Chang J, Jung JY. Poor prediction of potentially drug-resistant pathogens using current criteria of health care-associated pneumonia. Respir Med. 2012;106:1311–9.

    PubMed  Article  PubMed Central  Google Scholar 

  9. 9.

    Schreiber MP, Chan CM, Shorr AF. Resistant pathogens in nonnosocomial pneumonia and respiratory failure: is it time to refine the definition of health-care-associated pneumonia? Chest. 2010;137:1283–8.

    PubMed  Article  PubMed Central  Google Scholar 

  10. 10.

    El-Solh AA, Sikka P, Ramadan F, Davies J. Etiology of severe pneumonia in the very elderly. Am J Respir Crit Care Med. 2001;163:645–51.

    CAS  PubMed  Article  Google Scholar 

  11. 11.

    El Solh AA, Pietrantoni C, Bhat A, Bhora M, Berbary E. Indicators of potentially drug-resistant bacteria in severe nursing home-acquired pneumonia. Clin Infect Dis. 2004;39:474–80.

    PubMed  Article  Google Scholar 

  12. 12.

    Shorr AF, Zilberberg MD, Micek ST, Kollef MH. Prediction of infection due to antibiotic-resistant bacteria by select risk factors for health care-associated pneumonia. Arch Intern Med. 2008;168:2205–10.

    PubMed  Article  Google Scholar 

  13. 13.

    Xia LP, Bian LY, Xu M, Liu Y, Tang AL, Ye WQ. 16S rRNA gene sequencing is a non-culture method of defining the specific bacterial etiology of ventilator-associated pneumonia. Int J Clin Exp Med. 2015;8:18560–70.

    CAS  PubMed  PubMed Central  Google Scholar 

  14. 14.

    Brito V, Niederman MS. Healthcare-associated pneumonia is a heterogeneous disease, and all patients do not need the same broad-spectrum antibiotic therapy as complex nosocomial pneumonia. Curr Opin Infect Dis. 2009;22:316–25.

    PubMed  Article  Google Scholar 

  15. 15.

    Umeki K, Tokimatsu I, Yasuda C, Iwata A, Yoshioka D, Ishii H, Shirai R, Kishi K, Hiramatsu K, Matsumoto B, Kadota J. Clinical features of healthcare-associated pneumonia (HCAP) in a Japanese community hospital: comparisons among nursing home-acquired pneumonia (NHAP), HCAP other than NHAP, and community-acquired pneumonia. Respirology. 2011;16:856–61.

    PubMed  Article  Google Scholar 

  16. 16.

    Charlson ME, Pompei P, Ales KL, MacKenzie CR. A new method of classifying prognostic comorbidity in longitudinal studies: development and validation. J Chronic Dis. 1987;40:373–83.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  17. 17.

    Schmieder R, Edwards R. Quality control and preprocessing of metagenomic datasets. Bioinformatics. 2011;27:863–4.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  18. 18.

    Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnry J. 2011;17:10–2.

    Article  Google Scholar 

  19. 19.

    Zhang J, Kobert K, Flouri T, Stamatakis A. PEAR: a fast and accurate Illumina Paired-End reAd mergeR. Bioinformatics. 2014;30:614–20.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  20. 20.

    Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26:2460–1.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  21. 21.

    Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Pena AG, Goodrich JK, Gordon JI, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  22. 22.

    Yoon SH, Ha SM, Kwon S, Lim J, Kim Y, Seo H, Chun J. Introducing EzBioCloud: a taxonomically united database of 16S rRNA gene sequences and whole-genome assemblies. Int J Syst Evol Microbiol. 2017;67:1613–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  23. 23.

    Plantinga AM, Chen J, Jenq RR, Wu MC. pldist: ecological dissimilarities for paired and longitudinal microbiome association analysis. Bioinformatics. 2019;35:3567–75.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  24. 24.

    Kim KJ, Park J, Park SC, Won S. Phylogenetic tree-based microbiome association test. Bioinformatics. 2020;36:1000–6.

    PubMed  Google Scholar 

  25. 25.

    Morton JT, Marotz C, Washburne A, Silverman J, Zaramela LS, Edlund A, Zengler K, Knight R. Establishing microbial composition measurement standards with reference frames. Nat Commun. 2019;10:2719.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  26. 26.

    Kollef MH, Shorr A, Tabak YP, Gupta V, Liu LZ, Johannes RS. Epidemiology and outcomes of health-care-associated pneumonia: results from a large US database of culture-positive pneumonia. Chest. 2005;128:3854–62.

    PubMed  Article  Google Scholar 

  27. 27.

    Micek ST, Kollef KE, Reichley RM, Roubinian N, Kollef MH. Health care-associated pneumonia and community-acquired pneumonia: a single-center experience. Antimicrob Agents Chemother. 2007;51:3568–73.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  28. 28.

    Shindo Y, Sato S, Maruyama E, Ohashi T, Ogawa M, Hashimoto N, Imaizumi K, Sato T, Hasegawa Y. Health-care-associated pneumonia among hospitalized patients in a Japanese community hospital. Chest. 2009;135:633–40.

    PubMed  Article  PubMed Central  Google Scholar 

  29. 29.

    Chalmers JD, Taylor JK, Singanayagam A, Fleming GB, Akram AR, Mandal P, Choudhury G, Hill AT. Epidemiology, antibiotic therapy, and clinical outcomes in health care-associated pneumonia: a UK cohort study. Clin Infect Dis. 2011;53:107–13.

    PubMed  Article  PubMed Central  Google Scholar 

  30. 30.

    Park HK, Song JU, Um SW, Koh WJ, Suh GY, Chung MP, Kim H, Kwon OJ, Jeon K. Clinical characteristics of health care-associated pneumonia in a Korean teaching hospital. Respir Med. 2010;104:1729–35.

    PubMed  Article  PubMed Central  Google Scholar 

  31. 31.

    Kelly BJ, Imai I, Bittinger K, Laughlin A, Fuchs BD, Bushman FD, Collman RG. Composition and dynamics of the respiratory tract microbiome in intubated patients. Microbiome. 2016;4:7.

    PubMed  PubMed Central  Article  Google Scholar 

  32. 32.

    Lipsky BA, Goldberger AC, Tompkins LS, Plorde JJ. Infections caused by nondiphtheria corynebacteria. Rev Infect Dis. 1982;4:1220–35.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  33. 33.

    Clariot S, Constant O, Lepeule R, Fihman V, Razazi K, Cook F, Attias A, Merle JC, Hemery F, Levesque E, et al. Clinical relevance and impact of Corynebacterium isolation in lower respiratory tract of critically ill patients requiring mechanical ventilation. Infection. 2020;48:413–20.

    CAS  PubMed  Article  Google Scholar 

  34. 34.

    Yatera K, Mukae H. Corynebacterium species as one of the major causative pathogens of bacterial pneumonia. Respir Investig. 2020;58:131–3.

    PubMed  Article  Google Scholar 

  35. 35.

    Yang K, Kruse RL, Lin WV, Musher DM. Corynebacteria as a cause of pulmonary infection: a case series and literature review. Pneumonia (Nathan). 2018;10:10.

    Article  Google Scholar 

  36. 36.

    Yasuda I, Matsuyama H, Ishifuji T, Yamashita Y, Takaki M, Morimoto K, Sekino M, Yanagihara K, Fujii T, Iwaki M, et al. Severe pneumonia caused by toxigenic Corynebacterium ulcerans infection, Japan. Emerg Infect Dis. 2018;24:588–91.

    PubMed  PubMed Central  Article  Google Scholar 

  37. 37.

    Yatera K, Noguchi S, Yamasaki K, Kawanami T, Fukuda K, Naito K, Akata K, Kido T, Ishimoto H, Sakamoto N, et al. Determining the possible etiology of hospital-acquired pneumonia using a clone library analysis in Japan. Tohoku J Exp Med. 2017;242:9–17.

    PubMed  Article  PubMed Central  Google Scholar 

  38. 38.

    de Steenhuijsen Piters WA, Huijskens EG, Wyllie AL, Biesbroek G, van den Bergh MR, Veenhoven RH, Wang X, Trzcinski K, Bonten MJ, Rossen JW, et al. Dysbiosis of upper respiratory tract microbiota in elderly pneumonia patients. ISME J. 2016;10:97–108.

    PubMed  Article  CAS  Google Scholar 

  39. 39.

    Zakharkina T, Martin-Loeches I, Matamoros S, Povoa P, Torres A, Kastelijn JB, Hofstra JJ, de Wever B, de Jong M, Schultz MJ, et al. The dynamics of the pulmonary microbiome during mechanical ventilation in the intensive care unit and the association with occurrence of pneumonia. Thorax. 2017;72:803–10.

    PubMed  Article  PubMed Central  Google Scholar 

  40. 40.

    Bousbia S, Papazian L, Saux P, Forel JM, Auffray JP, Martin C, Raoult D, La Scola B. Repertoire of intensive care unit pneumonia microbiota. PLoS ONE. 2012;7:e32486.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  41. 41.

    Segal LN, Clemente JC, Tsay JC, Koralov SB, Keller BC, Wu BG, Li Y, Shen N, Ghedin E, Morris A, et al. Enrichment of the lung microbiome with oral taxa is associated with lung inflammation of a Th17 phenotype. Nat Microbiol. 2016;1:16031.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  42. 42.

    Kitsios GD, Fitch A, Manatakis DV, Rapport SF, Li K, Qin S, Huwe J, Zhang Y, Doi Y, Evankovich J, et al. Respiratory microbiome profiling for etiologic diagnosis of pneumonia in mechanically ventilated patients. Front Microbiol. 2018;9:1413.

    PubMed  PubMed Central  Article  Google Scholar 

  43. 43.

    Carratala J, Mykietiuk A, Fernandez-Sabe N, Suarez C, Dorca J, Verdaguer R, Manresa F, Gudiol F. Health care-associated pneumonia requiring hospital admission: epidemiology, antibiotic therapy, and clinical outcomes. Arch Intern Med. 2007;167:1393–9.

    PubMed  Article  Google Scholar 

  44. 44.

    Dickson RP, Schultz MJ, van der Poll T, Schouten LR, Falkowski NR, Luth JE, Sjoding MW, Brown CA, Chanderraj R, Huffnagle GB, et al. Lung microbiota predict clinical outcomes in critically ill patients. Am J Respir Crit Care Med. 2020;201:555–63.

    CAS  PubMed  Article  Google Scholar 

Download references


The authors thank all of study participants and the medical staff from Chuncheon Sacred Heart Hospital for their contributions to this study.


This study was supported by the Bio & Medical Technology Development Program of the National Research Foundation (NRF) funded by the Korean government (MSIT) (NRF-2017M3A9E8033225) and by the National Research Foundation of Korea Grant funded by Korean Government (NRF 2020R1A2C1011455). This research was supported by Hallym University Research Fund. Funding agencies had no role in the study design, data collection and analysis, decision to publish, manuscript preparation, or decision to submit the manuscript for publication. The remaining authors have no conflicts of interest to report.

Author information




MB, CB, SW, YK and HY had the idea for and designed the study and had full access to all data in the study and take responsibility for the integrity of the data analysis. JYH and JJL wrote the first full draft of the report. MB, NEK, CB, SW, YK, JJL, HY and JYH contributed to critical revision of the report. NEK, SW and HY contributed to the statistical analysis. All authors contributed to data acquisition, data analysis, or data interpretation, and reviewed the final version. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Hana Yi or Ji Young Hong.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Institutional Review Board of Chuncheon Sacred Heart Hospital (IRB approval number: 2017‐47). Informed consent was obtained from all patients or their legal representatives. This research was conducted in accordance with the Declaration of Helsinki.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1.

Statistics for the relative abundances of bacterial taxa in ETAs of NHAI group based on GLMM analysis. Table S2. Association between the relative abundances of the bacterial taxa in the ETAs and clinical outcomes: results from GLMM analyses.

Additional file 2: Figure S1.

Longitudinal patterns of the relative OTU abundances of the dominant genera in the NHAI and non-NHAI groups. For the lines in a box and whisker plot: the extreme bars are the overall range, the bottom and top of the box are the 25th and 75th percentiles and the line inside the box is the 50th percentile (median). *P < 0.05; **P < 0.01; ***P < 0.001; **** P < 0.0001.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Baek, Mg., Woo, S.J., Kim, N.E. et al. Respiratory microbiome profiles differ by recent hospitalization and nursing home residence in patients on mechanical ventilation. J Transl Med 18, 464 (2020).

Download citation


  • HCAP
  • Microbiome
  • Pneumonia
  • Mechanical ventilation