Skip to main content

Characteristics of the sputum microbiome in COPD exacerbations and correlations between clinical indices



Chronic obstructive pulmonary disease (COPD) is a prevalent, progressive respiratory disease, and acute exacerbations of COPD (AECOPD) can accelerate the deterioration of the disease. Increasing evidence suggests that airway bacterial dysbiosis is associated with AECOPD. However, the exact relationship between changes in the sputum microbiome during AECOPD and clinical indices remains unclear.


In this study, a total of 76 sputum samples were collected from patients with AECOPD (n = 28), stable COPD (n = 23), recovery (n = 15) and healthy controls (HCs; n = 10). The sputum microbiome profile was analysed by sequencing the V3‑V4 amplicon of the 16S rRNA (ribosomal RNA) gene.


The bacterial diversity (Shannon and Simpson’s index) was found to be significantly decreased in the AECOPD and recovery groups when compared to that in the stable COPD and HC groups. The most dominant phylum identified in the sputum samples of AECOPD patients was Proteobacteria, accounting for 30% of the microbiome. Compared to the stable COPD groups, the relative abundances of Firmicutes and Bacteroidetes were decreased, whereas those of Proteobacteria and Actinobacteria were increased in AECOPD patients. Furthermore, discriminative bacteria, such as Haemophilus, were identified as being specific taxa in AECOPD patients. Functional analysis showed that genes involved in membrane transport and signal transduction metabolism were enriched in the AECOPD group. Importantly, the proportions of Veillonella were positively correlated with lung function, and Staphylococcus was positively correlated with inflammatory indices.


Our study revealed variations in the sputum microbiome of AECOPD (based on composition and function) in a Chinese cohort and highlighted its correlation to clinical indices. These results indicated that microbial dysbiosis may contribute to disease progression and provide microbial biomarkers for the diagnosis of AECOPD.


Chronic obstructive pulmonary disease (COPD), one of the leading causes of morbidity and mortality worldwide, is characterized by irreversible airflow limitation due to airway and/or alveolar abnormalities [1]. Acute exacerbation of COPD (AECOPD) is defined as an acute worsening of respiratory symptoms, such as increased cough, purulent sputum production, and dyspnoea, which requires a change in treatment [2]. Frequent acute exacerbations lead to an increase in mortality as well as socioeconomic loss and health expenditures [3].

Bacteria-induced respiratory infection is the major cause of acute exacerbations of COPD [4]. Alterations in the airway microbiome are associated with decreased lung function and enhanced airway inflammation [5, 6]. Clinic indices such as spirometry indices are the indicators of the severity of airflow limitation and key to the diagnosis of COPD [7]. Inflammatory biomarkers can indicate bacterial exacerbation in COPD patients [8]. These clinical indices can be utilized to accurately profile disease severity and prognosis in COPD patients and may be influenced by microorganisms [9,10,11]. Recent evidence shows that different bacterial taxa identified in sputum samples from AECOPD patients are associated with different clinical outcomes [12]. However, the manner in which these detected bacteria are involved in pathogenesis and how they ultimately affect the clinical indices of the host remain unclear. Thus, it is necessary to investigate the correlation between the relating microbiota composition and the pathogenesis of acute exacerbation and clinical indices [13].

Traditional culture techniques have significant limitations regarding unculturable bacteria, as well as a poor sensitivity to detecting low-abundance bacteria [14]. In some situations, which mainly occurs for patients chronically colonized by Pseudomonas aeruginosa, chronic colonizing bacteria (rather than pathogenic bacteria) grow easily from the analysed sample by culturing, which can interfere with the potential cause of the exacerbations [15]. In recent years, advances in diagnostic technology, such as the 16S ribosomal RNA (rRNA) gene sequencing, have now allowed for the assessment of microbiome features with unprecedented personalization and precision, thus demonstrating alterations in the composition and relative abundance of bacteria related to COPD patients [16,17,18]. However, the exact relationship between changes in the sputum microbiome during AECOPD and clinical indices has not yet been elucidated by sequencing the 16S rRNA gene.

In this study, we investigated the correlation between key bacteria and certain clinical indices in COPD exacerbations by conducting 16S rRNA gene sequencing. We revealed that Streptococcus is the most dominant genus of AECOPD and found that Veillonella and Staphylococcus correlated with FEV1/FVC and C-reactive protein, respectively. These bacteria may serve as potential biomarkers and therapeutic targets for AECOPD.


Subject enrollment

From March 2019 to March 2021, 76 AECOPD patients who were admitted to the Respiratory Department of Peking Union Medical College Hospital and its affiliated hospitals were recruited for this study. Stable COPD patients and healthy controls were enrolled from outpatient clinics. Inclusion criteria were as follows: (a) diagnosed with COPD or AECOPD according to the GOLD criteria [1]; (b) age ≥ 40 years; (c) stable COPD patients were defined as clinically stable for at least 3 months free from acute exacerbation. Study exclusion criteria were: (a) a history of any other respiratory illnesses, such as asthma, pulmonary tuberculosis, concomitant bronchiectasis, sleep apnea syndrome, lung cancer, or interstitial lung disease; (b) underwent chronic treatment with oral corticosteroids or immunosuppressive drugs; (c) history of any antibiotics usage within the previous 4 weeks before enrollment. All AECOPD patients received antibiotics and inhaled corticosteroids (ICS) according to guidelines [19]. Patients in the recovery group were defined as 2 weeks post-therapy. They should meet discharge criteria and have discontinued antibiotics at the same time. All participants signed an informed written consent agreement.

Sample collection and DNA extraction

Fresh sputum samples were collected from AECOPD patients at admission before any therapy and recovery patients discontinued antibiotics. We collected either spontaneous or induced sputa from each patient. Participants were asked to rinse their mouths with a 3% hypertonic saline solution before sampling. Induced sputum samples were then collected via inhalation of 3% hypertonic saline and expectoration of sputum into a sputum cup within 30 min [20]. A total of 76 specimens were assigned to four clinical states: 1. AECOPD (n = 28); 2. Stable (n = 23); 3. Recovery (n = 15); 4. Healthy controls (n = 10). Then the samples were divided into 0.5 ml aliquots and stored at − 80 ℃ for DNA extraction and the remaining part of the samples were submitted for routine bacterial culture.

Total genome DNA from all sputum samples was extracted using the CTAB method following the manufacturer’s instructions. DNA concentration and purity were monitored on 1% agarose gels. According to the concentration, DNA was diluted to 1 ng/µL using sterile water.

Library preparation and sequencing

16S rRNA genes of distinct regions were amplified by a specific primer with the barcode. PCR reactions were performed using Phusion® High-Fidelity PCR Master Mix (NEB) as manufacturer’s instructions. Thermal cycling consisted of an initial denaturation (1 min at 98 ℃), followed by 30 cycles of denaturation (10 s at 98 ℃), annealing (30 s at 50 ℃), elongation (30 s at 72 ℃) and finally 72 ℃ for 5 min. Mix PCR products with the same volume of 1X loading buffer (SYB green) and operate electrophoresis on 2% agarose gel for detection. PCR products were then mixed in equidensity ratios and purified using Qiagen Gel Extraction Kit (Qiagen, Germany). Followed by the manufacturer's recommendations, sequencing libraries were generated by TruSeq® DNA PCR-Free Sample Preparation Kit (Illumina, USA), and then index codes were added. Library quality was assessed using the Qubit@ 2.0 Fluorometer (Thermo Scientific) and Agilent Bioanalyzer 2100 system. Finally, the library was then sequenced on the Illumina NovaSeq platform and subsequently generate 250 bp paired-end reads.

Bioinformatics analysis

The 16S rRNA sequencing data of 76 samples were examined for paired-end reads assembly and quality control. Sequences analysis was performed by Uparse software (Uparse v7.0.1001) [21]. Sequences were assigned to the same OTUs with similarity ≥ 97%. A representative sequence from each OTU was selected for further annotation. To annotate taxonomic information, each representative sequence uses the Silva Database based on the Mother algorithm [22]. To study the difference of the dominant species in different samples and the phylogenetic relationship of different OTUs, multiple sequence alignment analysis was carried out by the MUSCLE software (Version 3.8.31) [23]. OTUs abundance data were normalized by a standard of sequence number consistent with the sample with the least sequences. Further analysis of alpha and beta diversity analyses were all calculated based on the output normalized data using QIIME software (Version 1.9.1). Alpha diversity analysis is applied in assessing the complexity of species diversity within a sample and included observed species, Chao1, Shannon, and Simpson indices. Beta diversity analysis was applied to explore the differences of bacterial communities between COPD patients at different stages and healthy controls in species complexity. PCoA analysis was carried out by stat packages, WGCNA package, and ggplot2 packages. Microbial metabolic function profiles were predicted using Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt) software to generate the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways based on 16S rRNA gene sequence [24]. Finally, we calculated the Spearman correlation coefficient and corresponding p-value between each module and clinical traits.

Statistical analyses

SPSS 21.0 software was used for basic statistical analysis of data, GraphPad Prism9 and R software (Version 2.15.3) were used for chart production (p < 0.05 was considered statistically significant). Normally distributed variables were presented as mean ± standard deviation and compared through one-way analysis of variance (ANOVA) followed by Tukey’s multiple comparison test. The nonparametric data were presented as median (interquartile range) and were analyzed by Kruskal–Wallis tests. Categorical variables were expressed as count and percentages (N, %) and compared using the χ2 test. Specifically, alpha diversity was calculated by the Wilcox rank-sum tests, and beta diversity was acquired by Adonis analysis. Use linear discriminate analysis effect size (LEfSe) to obtain different species between groups. In this study, the linear discriminant analysis (LDA) score threshold used was > 4.0. Spearman correlation analysis was performed to evaluate the relationship between the sputum microbiome and clinical indices and then use the pheatmap function in the pheatmap package for visualization.


Characteristics of the study participants

A total of 76 sputum samples were collected from COPD patients at different stages of disease and healthy controls. The demographics and clinical features of the patients are summarized in Table 1. No statistical differences in age, gender, BMI, current smoker, comorbid conditions, medications, or exacerbation frequency were observed among the groups. The results of sputum culture showed no statistical difference may be due to the low positive rate. The smoking index was low in healthy people, but there was no significant difference among the COPD groups. As expected, lung function (FEV1%pred and FEV1/FVC values) was significantly decreased in all COPD patients when compared to that of healthy controls (p < 0.001). In addition, patients in the AECOPD group had a higher modified Medical Research Council (mMRC) dyspnoea score than stable COPD patients (p = 0.023). Blood routine parameters, glycemia, arterial blood gas analysis, inflammation index values (PCT and CRP), and GOLD classification were not found to be significantly different among the four groups.

Table 1 Major clinical characteristics of participants in this study

16S gene sequence‑based characterization

The samples yielded a total of 4053 operational taxonomic units (OTUs) at 97% sequence identity (2677, 1894, 2127, and 988 in the AECOPD, recovery, stable COPD, and healthy control groups, respectively). We measured alpha diversity to evaluate the diversity and richness of the bacterial community (Fig. 1a). The Shannon index and Simpson’s index represent microbial diversity, while the Chao1 index and observed species index reflect microbial richness. AECOPD samples had significantly decreased indices of diversity compared to those in the stable COPD and healthy control groups (p < 0.05). In addition, the AECOPD samples tended to have lower bacterial richness than the stable COPD and healthy control samples, but the differences were not statistically significant (p > 0.05). Samples collected from the recovery group had the lowest alpha diversity measures (p < 0.05).

Fig. 1
figure 1

Bacterial community structure in healthy controls and COPD patients at different stages (AECOPD, recovery, and stable COPD, respectively). a Comparison of sputum microbiome alpha richness and diversity (Shannon index, Simpson’s index, Chao1, and Observed species) by Wilcoxon rank-sum test. *p < 0.05, **p < 0.01. b Principal coordinates analysis (PCoA) based on Bray–Curtis’s distance is represented. Each point represents one subject and different colors represent different groups. ADONIS analysis showed that the separation of bacterial communities was significant (p < 0.05) as described in the text, except for stable vs controls (p = 0.066), AECOPD vs recovery (p = 0.057). c Venn diagram showing different OTUs among groups

In addition, we calculated the Bray–Curtis distance to investigate the clustering of microbial taxa in the four groups. Principal coordinates analysis (PCoA) showed that the separation of microbial communities was significant in the four groups (p < 0.05). Significant differences in microbial community structure were observed between the AECOPD and stable groups (p = 0.001), AECOPD and healthy controls (p = 0.035), the healthy control and recovery groups (p = 0.004), the stable COPD and recovery groups (p = 0.001) (Fig. 1b).

Venn diagrams distinguished between the sputum microbiomes of the four groups. As shown in Fig. 1c, there were 3467 OTUs in the AECOPD and stable COPD samples, of which 1337 (38.56%) were shared by the two groups; there were 2909 OTUs in the AECOPD and healthy controls samples, of which 756 (25.99%) were shared by the two groups; there were 3256 OTUs in the AECOPD and recovery samples, of which 1315 (40.39%) were shared by the two groups; there were 2830 OTUs in the recovery and stable COPD samples, of which 1191 (42.08%) were shared by the two groups; there were 2159 OTUs in the recovery and healthy controls samples, of which 723 (33.49%) were shared by the two groups; and there were 2371 OTUs in the stable COPD and healthy controls samples, of which 744 (31.38%) were shared by the two groups.

Abundance of bacterial taxa changes among groups

At the phylum level, the numbers of phyla detected in the AECOPD, recovery, stable COPD, and healthy control groups were 35, 29, 24, and 14, respectively. The dominant bacterial phylum in the AECOPD group was Proteobacteria (30.29%), followed by Firmicutes (29.85%) and Bacteroidetes (14.02%). In the stable COPD group, the major phylum was Firmicutes (31.63%), followed by Bacteroidetes (28.94%), and Proteobacteria (19.68%). In the recovery group, the major phyla were Firmicutes (44.04%), Proteobacteria (21.94%), and Bacteroidetes (13.35%). Lastly, in the healthy controls group, Firmicutes (34.01%), Bacteroidetes (26.01%), and Proteobacteria (23.09%) were predominant (Fig. 2a).

Fig. 2
figure 2

Relative abundance of the most prevalent bacterial phyla (a) and genera (b) in healthy controls and COPD patients at different stages (AECOPD, recovery, and stable COPD, respectively)

At the genus level, we detected 506, 399, 397, and 270 genera in the AECOPD, recovery, stable COPD, and healthy control groups, respectively (Fig. 2b). The dominant genera in the AECOPD group were Streptococcus (14.31%), Neisseria (11.60%), unidentified_Prevotellaceae (8.90%), Haemophilus (7.49%), and Veillonella (6.37%). In the stable COPD group, the dominant genera were unidentified_Prevotellaceae (15.69%), Streptococcus (14.31%), Neisseria (12.13%), Veillonella (7.38%), and Haemophilus (3.62%). In the recovery group, Streptococcus (27.32%) was the most highly represented genus, followed by Neisseria (8.69%), unidentified_Prevotellaceae (6.68%), Veillonella (4.35%), unidentified_Corynebacteriaceae (3.98%), and Staphylococcus (3.79%). Streptococcus was the most common genus in the healthy controls (18.41%), followed by unidentified_Prevotellaceae (14.63%), Neisseria (13.21%), Veillonella (9.00%), and Haemophilus (5.93%).

Sputum microbiota alterations in AECOPD patients at the level of phylum and genus

At the phylum level, the relative abundance of Bacteroidetes was significantly decreased while those of Actinobacteria were increased in AECOPD patients compared to stable COPD and healthy controls (Fig. 3a, b). Although there were no statistical differences presented, we observed a decrease in Firmicutes and an increase in Proteobacteria in the AECOPD group compared to the stable COPD group. At the genus level, we observed significantly increased proportions of Rothia, unidentified_Corynebacteriaceae, Stenotrophomonas, and significantly decreased proportions of Prevotella, Alloprevotella, Porphyromonas, unidentified_Prevotellaceae in AECOPD patients compared to stable COPD(Fig. 3c, d).

Fig. 3
figure 3

Heatmap and barplot showing the statistic differences of relative abundance of taxa among healthy controls, and COPD patients at different stages (AECOPD, recovery, and stable COPD, respectively) at phylum (a, b) and genus (c, d) levels. *, ** means p < 0.05, p < 0.01, respectively

Bacterial taxonomic differences among the four groups

To identify biomarkers from the disease groups, we performed a LEfSe analysis at the genera level. As expected, discriminative bacteria were identified in the four groups. We selected a linear discriminant analysis (LDA) to score higher than 4 to represent a significantly enriched genus in each group (Fig. 4a, b). The sputum microbiome of the AECOPD group was characterized by a dominance of Haemophilus, Pasteurellaceae, etc., whereas the microbiome in the stable COPD group was dominated by the genus Prevotella, Clostridiales, etc. (p < 0.05).

Fig. 4
figure 4

LEfSe analysis revealed the discriminative microbes in AECOPD, stable COPD, control, and recovery patient groups. a The histogram indicates the LDA score. These taxa showed statistically significant differences between the four groups (p < 0.05 by the Wilcoxon test). We selected the LDA threshold value of > 4 as significant in each group. b The cladogram of major differentially abundant taxa is based on LEfSe analysis. The circles radiating from inside to outside represent the taxonomic level from phylum to genus. Each small circle at a different classification level represents a classification at that level, and the diameter of the small circle is proportional to relative abundance. Different species biomarker follows the group for coloring. For example, the red node indicates the group of microorganisms that play an important role in the AECOPD group. The diameter of each circle represents the relative abundance of each taxon

Functional analysis of the microbiome in AECOPD patients by PICRUSt analysis

We performed a PICRUSt analysis based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) database to evaluate sputum microbiome functions across groups in our study cohort [24]. Specifically, we found that the levels of metabolism, such as replication and repair, translation, nucleotide metabolism, glycan biosynthesis and metabolism, cell growth and death, and biosynthesis of other secondary metabolites, were decreased in the AECOPD group than in the stable COPD group. In contrast, metabolism related to membrane transport and signal transduction, including ABC transporters and secretion system were enriched in the AECOPD group (p < 0.05; Fig. 5a, b). The above results indicate that the sputum microbiome affects genetic information processing, nucleotide, and sugar metabolism.

Fig. 5
figure 5

The differences in metabolic pathways by PICRUSt functional analysis. a Related KEGG pathways are plotted as a bar graph at KEGG level 2 between two groups. b Comparison of functional analysis in healthy controls and COPD patients at different stages at KEGG level 3

The sputum microbiome of AECOPD patients was associated with clinical indices

We used Spearman’s correlation analysis to evaluate the correlation between each group’s sputum microbiome and clinical index parameters, including the demographics and clinical features of the AECOPD patients, as shown in Fig. 6. Strong correlations (correction |r|> 0.4, p < 0.05) were found among 5 taxa and the 3 clinical indices in COPD subjects (Fig. 6a). Among the most abundant genera, the relative abundance of Veillonella exhibited a significant positive correlation with FEV1/FVC (Fig. 6b). Staphylococcus was positively correlated with inflammation index value CRP, whereas Alloprevotella was negatively correlated with CRP (Fig. 6c and Additional file 1: Fig. S1). Another finding was that Haemophilus and Prevotella were negatively correlated with the mMRC Dyspnoea scale (Additional file 1: Fig. S1).

Fig. 6
figure 6

Spearman correlation between major sputum microbiome and clinical indices in COPD patients. a Heatmap of Spearman correlation analysis between the relative abundances of sputum microbiome and the clinical indices. *p < 0.05, **p < 0.01. b Veillonella shows a positive correlation with FEV1/FVC. c Staphylococcus shows a positive correlation with CRP


In this study, we investigated the characteristics of sputum microbiota in patients with acute exacerbation and other disease states of COPD patients via cross-sectional observation by using high-throughput 16S rRNA sequencing. It was found that the lung microbiome of AECOPD and recovery patients demonstrated lower bacterial richness, as well as significantly different compositions, compared to those of the stable COPD and healthy control group subjects. We also found alterations in metabolic pathways, and common clinical indices were associated with AECOPD microbial communities. Collectively, these differences in our findings indicate that microbial dysbiosis may contribute to AECOPD and highlight its correlation to clinical indices, which could be considered microbial biomarkers used in potential targets for therapeutic interventions.

Our results demonstrated that microbiome diversity was reduced in AECOPD patients, compared with stable COPD patients and healthy controls. Samples collected from the recovery group had the lowest alpha diversity measures, which may likely be due to antibiotic exposure during clinical treatments. A large longitudinal study showed an overall reduction in microbial diversity during COPD exacerbations compared to samples from stable patients [25]. These results were consistent with the present study, thus suggesting that microbial diversity may be a biological indicator of AECOPD.

In our study, the most dominant phylum in sputum samples of AECOPD patients was Proteobacteria, which was consistent with the previous reports that the occurrence of AECOPD was related to the reduction in microbiome diversity and increased proportion of Proteobacteria [25, 26]. Proteobacteria and its members such as Haemophilus were associated with the onset of AECOPD [27]. In the present study, the most dominant phylum in sputum samples from the stable COPD group was Firmicutes, consistent with previous findings [28, 29].

At the genus level, the top 5 dominant genera present in all sputum samples were Streptococcus, Neisseria, Prevotellaceae, Haemophilus, and Veillonella, which has been reported in a previous study [28] except Prevotellaceae. In the present study, Streptococcus and Prevotellaceae were the most dominant genus in AECOPD and stable COPD groups respectively. Streptococcus has been reported as the most common genus among all samples in COPD patients [29, 30]. Currently, Prevotellaceae and Streptococcus, which both belong to the Firmicutes, are considered members of the core pulmonary microbiome [31].

Interestingly, compared to the stable COPD group, the AECOPD group had a higher level of Actinobacteria and the associated genera Corynebacteriaceae and Rothia. The AECOPD group also had a higher level of Proteobacteria and the associated genera Stenotrophomonas and Haemophilus, which are common COPD-related pathogens reported in the previous study [32, 33]. The level of Bacteroidetes and the associated genus Prevotella and Porphyromonas were markedly lower in the AECOPD group than in the stable COPD and healthy controls groups. It has been reported that the proportion of Prevotella of Bacteroidetes in the airways of asthma patients was lower than that of the normal population [34]. These results may indicate that alteration of normal flora distribution may lead to enhanced inflammation and increased exacerbation risk.

Functional prediction showed that the transporters, peptidases, purine/pyrimidine, and nucleotide metabolic pathways altered in the AECOPD group. For example, ABC transporters and the secretion system were enriched in AECOPD groups compared to the controls. ABC transporters and their role in substrate transport across the bacterial membrane have been associated with antibiotic resistance [35]. The bacterial secretion system contributes to secrete virulence factors for host invasion [36]. These pathways are important for the survival of pathogenic bacteria. Furthermore, AECOPD groups exhibited decreased glycan biosynthesis and metabolism relative to the stable COPD group and controls. Similar to previous findings, impaired glucose metabolism was also found in COPD patients [37]. The above results further indicate that the change in microbiome function may influence the occurrence or development of AECOPD.

In this study, our results indicated that the sputum microbiome was associated with disease indices. Veillonella was found to have a significant positive association with FEV1% (p < 0.05) and FEV1/FVC (p < 0.001). Additionally, Staphylococcus had a highly significant correlation with the inflammatory index CRP (p < 0.01). The proportion of Veillonella was slightly decreased in the AECOPD group (6.37%) when compared to that of the stable COPD group (7.38%) and healthy controls (9.00%) in our cohort. Simultaneously, the proportion of Staphylococcus was slightly increased in the AECOPD group (0.79%) compared to that of the stable COPD group (0.03%) and healthy controls (0.04%) in our cohort. Interestingly, Filho et al. found that a higher abundance of Staphylococcus was associated with decreased 1-year mortality while the survivors were abundant in Veillonella [12]. The results of our study may explain the reason for this result. Veillonella is a gram-negative anaerobic bacterium that belongs to normal oral commensals [38]. Veillonella species were observed to be greater in the healthy populations than in the COPD groups, thus suggesting a beneficial role for symbiosis in health and/or disease states [39]. Staphylococcus, a group of gram-positive pathogens that cause infections in the lungs [40], has been shown to result in the need for more antibiotics and longer hospitalizations [41]. Staphylococcus aureus in particular has been reported to directly trigger the formation of neutrophil extracellular traps (NETs), which may consequently influence the cycle of inflammation in COPD [42]. Taken together, these results indicate that Veillonella and Staphylococcus may be involved in COPD pathogenesis through their effects on lung function and local inflammatory responses, respectively. AECOPD patients may suffer a lung microecological imbalance, thereby aggravating the inflammatory responses and airflow limitations, which eventually increases the risk of readmission and mortality. Therefore, the sputum microbiome may be used to identify the clinical outcome and prognosis of COPD patients.

There were several limitations to our study. First, samples were cross-sectional and collected from different patients with different disease states, and the interference of confounding factors was unavoidable. Second, in addition to bacteria, viruses and fungi are also involved in the pathogenesis of acute exacerbations of COPD [33]. However, due to the lack of these data, we did not take them into account in the present study. Finally, our results were mainly restricted to the urban area of Beijing. Future large-scale multicentre studies with different biogeographical backgrounds are required to validate our findings.


In conclusion, this study demonstrates that alterations in the abundance and composition of the sputum microbiome in COPD patients progress from stable to exacerbated states. Additionally, there are metabolic pathway alterations in COPD patients at different stages and the sputum microbiome is associated with common clinical indicators. These findings may provide new insight into the potential use of the sputum microbiome in AECOPD patients.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.



Chronic obstructive pulmonary disease


Acute exacerbations of chronic obstructive pulmonary disease


Healthy controls


Body mass index


Forced expiratory volume in 1 s


Forced vital capacity




C-reactive protein


Global initiative for obstructive lung disease


Modified Medical Research Council Dyspnoea scale


White blood cell count


Inhaled corticosteroid


Serum creatinine


Principal coordinates analysis


Linear discriminant analysis effect size


Linear discriminant analysis


  1. Vestbo J, Hurd SS, Agusti AG, Jones PW, Vogelmeier C, Anzueto A, Barnes PJ, Fabbri LM, Martinez FJ, Nishimura M, et al. Global strategy for the diagnosis, management, and prevention of chronic obstructive pulmonary disease: GOLD executive summary. Am J Respir Crit Care Med. 2013;187:347–65.

    CAS  PubMed  Google Scholar 

  2. Viniol C, Vogelmeier CF. Exacerbations of COPD. Eur Respir Rev. 2018;27:170103.

    PubMed  Google Scholar 

  3. McGhan R, Radcliff T, Fish R, Sutherland E, Welsh C, Make B. Predictors of rehospitalization and death after a severe exacerbation of COPD. Chest. 2007;132:1748–55.

    PubMed  Google Scholar 

  4. Ko FW, Chan KP, Hui DS, Goddard JR, Shaw JG, Reid DW, Yang IA. Acute exacerbation of COPD. Respirology. 2016;21:1152–65.

    PubMed  PubMed Central  Google Scholar 

  5. Wilkinson T, Patel I, Wilks M, Donaldson G, Wedzicha J. Airway bacterial load and FEV1 decline in patients with chronic obstructive pulmonary disease. Am J Respir Crit Care Med. 2003;167:1090–5.

    PubMed  Google Scholar 

  6. Dicker AJ, Huang JTJ, Lonergan M, Keir HR, Fong CJ, Tan B, Cassidy AJ, Finch S, Mullerova H, Miller BE, et al. The sputum microbiome, airway inflammation, and mortality in chronic obstructive pulmonary disease. J Allergy Clin Immunol. 2021;147:158–67.

    CAS  PubMed  Google Scholar 

  7. Pellegrino R, Viegi G, Brusasco V, Crapo R, Burgos F, Casaburi R, Coates A, van der Grinten C, Gustafsson P, Hankinson J, et al. Interpretative strategies for lung function tests. Eur Respir J. 2005;26:948–68.

    CAS  PubMed  Google Scholar 

  8. Peng C, Tian C, Zhang Y, Yang X, Feng Y, Fan H. C-reactive protein levels predict bacterial exacerbation in patients with chronic obstructive pulmonary disease. Am J Med Sci. 2013;345:190–4.

    PubMed  Google Scholar 

  9. Shaw JG, Vaughan A, Dent AG, O’Hare PE, Goh F, Bowman RV, Fong KM, Yang IA. Biomarkers of progression of chronic obstructive pulmonary disease (COPD). J Thorac Dis. 2014;6:1532–47.

    PubMed  PubMed Central  Google Scholar 

  10. Soumagne T, Levesque F, Milot J, Godbout K, Lacasse Y, Maltais F. Significance of Stenotrophomonas maltophilia when detected in sputum of ambulatory patients with COPD. Int J Chron Obstruct Pulmon Dis. 2021;16:2895–900.

    PubMed  PubMed Central  Google Scholar 

  11. Jones TPW, Brown J, Hurst JR, Vancheeswaran R, Brill S. COPD exacerbation phenotypes in a real-world five year hospitalisation cohort. Respir Med. 2020;167:105979.

    PubMed  Google Scholar 

  12. Leitao Filho FS, Alotaibi NM, Ngan D, Tam S, Yang J, Hollander Z, Chen V, FitzGerald JM, Nislow C, Leung JM, et al. Sputum microbiome is associated with 1-year mortality after chronic obstructive pulmonary disease hospitalizations. Am J Respir Crit Care Med. 2019;199:1205–13.

    CAS  PubMed  Google Scholar 

  13. Sethi S. Chronic obstructive pulmonary disease and infection. Disruption of the microbiome? Ann Am Thorac Soc. 2014;11:S43-47.

    PubMed  Google Scholar 

  14. Han MK, Huang YJ, Lipuma JJ, Boushey HA, Boucher RC, Cookson WO, Curtis JL, Erb-Downward J, Lynch SV, Sethi S, et al. Significance of the microbiome in obstructive lung disease. Thorax. 2012;67:456–63.

    PubMed  Google Scholar 

  15. Monso E. Microbiome in chronic obstructive pulmonary disease. Ann Transl Med. 2017;5:251.

    PubMed  PubMed Central  Google Scholar 

  16. Noecker C, McNally CP, Eng A, Borenstein E. High-resolution characterization of the human microbiome. Transl Res. 2017;179:7–23.

    CAS  PubMed  Google Scholar 

  17. Galiana A, Aguirre E, Rodriguez JC, Mira A, Santibanez M, Candela I, Llavero J, Garcinuno P, Lopez F, Ruiz M, et al. Sputum microbiota in moderate versus severe patients with COPD. Eur Respir J. 2014;43:1787–90.

    PubMed  Google Scholar 

  18. Sze MA, Dimitriu PA, Suzuki M, McDonough JE, Campbell JD, Brothers JF, Erb-Downward JR, Huffnagle GB, Hayashi S, Elliott WM, et al. Host response to the lung microbiome in chronic obstructive pulmonary disease. Am J Respir Crit Care Med. 2015;192:438–45.

    CAS  PubMed  PubMed Central  Google Scholar 

  19. Jacobs D, Pandit U, Sethi S. Acute exacerbations in chronic obstructive pulmonary disease: should we use antibiotics and if so, which ones? Curr Opin Infect Dis. 2019;32:143–51.

    PubMed  Google Scholar 

  20. Pin I, Gibson PG, Kolendowicz R, Girgis-Gabardo A, Denburg JA, Hargreave FE, Dolovich J. Use of induced sputum cell counts to investigate airway inflammation in asthma. Thorax. 1992;47:25–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  21. Edgar R. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods. 2013;10:996–8.

    CAS  PubMed  Google Scholar 

  22. Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, Peplies J, Glöckner F. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41:D590-596.

    CAS  PubMed  Google Scholar 

  23. Edgar R. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792–7.

    CAS  PubMed  PubMed Central  Google Scholar 

  24. Langille M, Zaneveld J, Caporaso J, McDonald D, Knights D, Reyes J, Clemente J, Burkepile D, Vega Thurber R, Knight R, et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat Biotechnol. 2013;31:814–21.

    CAS  PubMed  PubMed Central  Google Scholar 

  25. Wang Z, Bafadhel M, Haldar K, Spivak A, Mayhew D, Miller BE, Tal-Singer R, Johnston SL, Ramsheh MY, Barer MR, et al. Lung microbiome dynamics in COPD exacerbations. Eur Respir J. 2016;47:1082–92.

    PubMed  Google Scholar 

  26. Ma N, Qi Y, Liang X, Bai J, Deng J, Li M, He Z. Compare the effect of inhaled corticosteroids and systemic corticosteroids on sputum microbiome of AECOPD. Front Med. 2021;8:637246.

    Google Scholar 

  27. Mayhew D, Devos N, Lambert C, Brown JR, Clarke SC, Kim VL, Magid-Slav M, Miller BE, Ostridge KK, Patel R, et al. Longitudinal profiling of the lung microbiome in the AERIS study demonstrates repeatability of bacterial and eosinophilic COPD exacerbations. Thorax. 2018;73:422–30.

    PubMed  Google Scholar 

  28. Yang CY, Li SW, Chin CY, Hsu CW, Lee CC, Yeh YM, Wu KA. Association of exacerbation phenotype with the sputum microbiome in chronic obstructive pulmonary disease patients during the clinically stable state. J Transl Med. 2021;19:121.

    CAS  PubMed  PubMed Central  Google Scholar 

  29. Wang J, Chai J, Sun L, Zhao J, Chang C. The sputum microbiome associated with different sub-types of AECOPD in a Chinese cohort. BMC Infect Dis. 2020;20:610.

    CAS  PubMed  PubMed Central  Google Scholar 

  30. Pragman AA, Lyu T, Baller JA, Gould TJ, Kelly RF, Reilly CS, Isaacson RE, Wendt CH. The lung tissue microbiota of mild and moderate chronic obstructive pulmonary disease. Microbiome. 2018;6:7.

    PubMed  PubMed Central  Google Scholar 

  31. Diao W, Shen N, Du Y, Qian K, He B. Characterization of throat microbial flora in smokers with or without COPD. Int J Chron Obstruct Pulmon Dis. 2017;12:1933–46.

    CAS  PubMed  PubMed Central  Google Scholar 

  32. Oladunjoye O, Oladunjoye A, Oladiran O, Donato A. Stenotrophomonas maltophilia infection in a patient with acute exacerbation of chronic obstructive pulmonary disease (COPD): a colonizer or true infection? Am J Case Rep. 2020;21:e924577.

    PubMed  PubMed Central  Google Scholar 

  33. Leung JM, Tiew PY, Mac Aogain M, Budden KF, Yong VF, Thomas SS, Pethe K, Hansbro PM, Chotirmall SH. The role of acute and chronic respiratory colonization and infections in the pathogenesis of COPD. Respirology. 2017;22:634–50.

    PubMed  PubMed Central  Google Scholar 

  34. Guo MY, Chen HK, Ying HZ, Qiu FS, Wu JQ. The role of respiratory flora in the pathogenesis of chronic respiratory diseases. Biomed Res Int. 2021;2021:6431862.

    PubMed  PubMed Central  Google Scholar 

  35. Wilson D. The ABC of ribosome-related antibiotic resistance. MBio. 2016;7:e00598.

    CAS  PubMed  PubMed Central  Google Scholar 

  36. Depluverez S, Devos S, Devreese B. The role of bacterial secretion systems in the virulence of gram-negative airway pathogens associated with cystic fibrosis. Front Microbiol. 2016;7:1336.

    PubMed  PubMed Central  Google Scholar 

  37. Agarwal AR, Kadam S, Brahme A, Agrawal M, Apte K, Narke G, Kekan K, Madas S, Salvi S. Systemic immuno-metabolic alterations in chronic obstructive pulmonary disease (COPD). Respir Res. 2019;20:171.

    PubMed  PubMed Central  Google Scholar 

  38. Man W, de Steenhuijsen PW, Bogaert D. The microbiota of the respiratory tract: gatekeeper to respiratory health. Nat Rev Microbiol. 2017;15:259–70.

    CAS  PubMed  PubMed Central  Google Scholar 

  39. Einarsson GG, Comer DM, McIlreavey L, Parkhill J, Ennis M, Tunney MM, Elborn JS. Community dynamics and the lower airway microbiota in stable chronic obstructive pulmonary disease, smokers and healthy non-smokers. Thorax. 2016;71:795–803.

    CAS  PubMed  Google Scholar 

  40. Tong SYC, Davis JS, Eichenberger E, Holland TL, Fowler VG. Staphylococcus aureus infections: epidemiology, pathophysiology, clinical manifestations, and management. Clin Microbiol Rev. 2015;28:603–61.

    CAS  PubMed  PubMed Central  Google Scholar 

  41. Narewski ER, Kim V, Marchetti N, Jacobs MR, Criner GJ. Is methicillin-resistant Staphylococcus aureus colonization associated with worse outcomes in COPD hospitalizations? Chronic Obstr Pulm Dis. 2015;2:252–8.

    PubMed  PubMed Central  Google Scholar 

  42. Uddin M, Watz H, Malmgren A, Pedersen F. NETopathic inflammation in chronic obstructive pulmonary disease and severe asthma. Front Immunol. 2019;10:47.

    CAS  PubMed  PubMed Central  Google Scholar 

Download references


We gratefully acknowledge all of the study participants and the medical staff who helped collect the samples in this study.


This work was supported by the National Key Research and Development Program (No. 2018YFC1315103) from the Ministry of Science and Technology of China and the National Key Technology R&D Program of China (No. 2013BAI09B10). The funding body had no role in the design of the study, the collection, analysis, or interpretation of data, or in writing the manuscript.

Author information

Authors and Affiliations



YX, HBZ, HMZ, JW, and YXQ conceived of and designed the experimental protocol; LFS, YXQ, JML, and RH contributed to patients’ recruitment and clinical management; LFS and HMZ were involved in data processing and analysis; LFS drafted the manuscript; ZWL assisted the statistical analysis; LFS, HMZ, ZWL and YX revised and edited the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Hongmei Zhao, Jing Wang or Yi Xiao.

Ethics declarations

Ethics approval and consent to participate

All experimental protocols in this study were approved by the ethics committee of Peking Union Medical College Hospital (no.B306).

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: Fig. S1.

Spearman correlations between major sputum microbiome with clinical indices in COPD patients. a. Alloprevotella shows negative correlation with CRP. Prevotella (b) and Haemophilus (c) show a negative correlation with the MRC Dyspnoea scale.

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

Su, L., Qiao, Y., Luo, J. et al. Characteristics of the sputum microbiome in COPD exacerbations and correlations between clinical indices. J Transl Med 20, 76 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Chronic obstructive pulmonary disease
  • 16S ribosomal RNA gene sequencing
  • Sputum microbiome