Relationship between the respiratory microbiome and the severity of airflow limitation, history of exacerbations and circulating eosinophils in COPD patients
BMC Pulmonary Medicine volume 19, Article number: 112 (2019)
The respiratory microbiome is altered in COPD patients but its relationship with core components of the disease, such as the severity of airflow limitation, the frequency of exacerbations or the circulating levels of eosinophils, is unclear.
Cross-sectional study comprising 72 clinically stable COPD patients (mean age 68 [SD 7.9] years; FEV1 48.7 [SD 20.1]% of reference) who provided spontaneous sputum samples for 16S rRNA gene amplification and sequencing. The microbiome composition was analysed with QIIME.
We observed that: (1) more severe airflow limitation was associated with reduced relative abundance (RA) of Treponema and an increase in Pseudomonas; (2) patients with ≥2 exacerbations the previous year showed a significantly different bacterial community with respect to non-exacerbators (p = 0.014), with changes in 13 genera, including an increase of Pseudomonas, and finally, (3) peripheral eosinophils levels ≥2% were associated with more diverse microbiome [Chao1 224.51 (74.88) vs 277.39 (78.92) p = 0.006; Shannon 3.94 (1.05) vs 4.54 (1.06) p = 0.020], and a significant increase in the RAs of 20 genera.
The respiratory microbiome in clinically stable COPD patients varies significantly according to the severity of airflow limitation, previous history of exacerbations and circulating eosinophils levels.
Summary at a glance
Core components of COPD such as airflow limitation, history of previous exacerbations and level of circulating eosinophils have an impact in the bronchial respiratory microbiome of clinically stable COPD patients.
Chronic obstructive pulmonary disease (COPD) is a heterogeneous disease [1,2,3]. The study of the respiratory microbiome in COPD has revealed a specific bacterial community composition in these patients [4, 5]. However, the relationship between this microbiome and core components of the disease, such as, the severity of the airflow limitation and the type of treatment received remains unclear. In addition, changes in the microbiome have been described in COPD exacerbations [6, 7], but it is not known if differences between patients who suffer two or more exacerbations per year, who are considered frequent exacerbators [8, 9], and non-exacerbators can be detected during clinical stability. Likewise, levels of circulating eosinophils ≥2% in clinically stable patients identifies a subgroup of COPD patients who are prone to recurrent exacerbations and are more responsive to treatment [10,11,12], but it is not clear if this is associated with a different respiratory microbiome. This work sought to investigate these questions.
Methods are detailed in the Additional file 1 and summarized below.
Study design and ethics
This is a cross-sectional, prospective, uncontrolled, multicentre, observational study. The study protocol was approved by the Ethics Committees of the participating hospitals (IMIM-Hospital del Mar, Hospital Universitari Parc Taulí, Hospital Clinic, Hospital 12 Octubre, Fundación Jimenez Díaz and Hospital Son Espases), and all patients included signed their informed consent.
Current or former smokers (≥ 10 pack-year) with stable COPD, attending the outpatients’ clinics of five Spanish hospitals between 2014 and 2016 were included in this study. The diagnosis and severity staging of COPD was established in accordance with GOLD criteria . Exclusion criteria were: age less than 40 years; a lifetime diagnosis of asthma, cystic fibrosis, bronchiectasis or cancer; patients receiving long-term treatment with oral corticosteroids or immunosuppressants; any comorbidity limiting cognitive capabilities and ≥ 3 admissions or 1 episode severe enough to require more than 30 days in hospital the previous year. Patients who had been treated with short-term antibiotics and/or corticosteroids at any time during the previous three months were considered unstable and not considered for the study.
Variables and measurements
Sociodemographic data were recorded by specific questionnaires. Lung function values during stability were obtained from the most recently available forced spirometry with reversibility testing performed according to standard techniques the previous year . Peripheral blood cell counts were obtained at enrolment and used to identify patients with ≥2% circulating eosinophils . Episodes of increased dyspnoea, sputum production and/or purulence during the previous year were identified and considered as exacerbations when treated with antibiotics and/or corticosteroids [15, 16]. Participants were considered as frequent exacerbators (FE) when they reported ≥2 exacerbations the previous year.
Sample collection, DNA extraction, PCR amplification and 16S sequencing
Spontaneous sputum samples were collected and processed within 60 min on the day of the visit. Sputum quality was assessed according to Murray-Washington criteria  and only samples with > 25 leucocytes per field (M-W ≥ 3) were considered for the study. Sputum samples were frozen until processing, which was carried out in a certified BSL2 hood with appropriate laminar flow.16S rRNA gene was amplified following the 16S Metagenomic Sequencing Library Preparation Illumina protocol (Part # 15044223 Rev. A, Illumina, CA, USA). Details are provided in the Additional file 1.
The Quantitative Insights Into Microbial Ecology (QIIME) pipeline 1.9.0  was used for sequence processing to obtain taxonomic information. Further technical details are provided in the Additional file 1.
Details are provided in the Additional file 1. In brief, categorical variables are expressed as absolute and relative frequencies, and continuous variables as means and standard deviations (SD) when the distribution was normal, or as medians and interquartile range (IQR) otherwise. Linear discriminant analysis Effect Size (LEfSe) was used to identify the differentially abundant taxa that explained the differences between the groups of participants. The threshold value of the logarithmic LDA score for discriminative features was 2.0. Bacterial α-diversity was assessed through the Chao1 estimator  and the Shannon index , calculating both indexes after subsampling with QIIME so as to avoid sequencing effort bias. Principal Coordinates Analysis (PCoA) with Bray-Curtis dissimilarity index  was used to study community composition, assessing the statistical significance of the differences in sample groupings through Adonis testing. Interaction between independent variables was assessed through stratification and multivariate analyses with α-diversity as dependent variable. Statistical tests used in the study were two-sided, and a p value of 0.05 or less was reported as statistically significant. Statistical analyses were performed using the SPSS statistical software package version 18 (SPSS Inc., Chicago, IL, USA).
Table 1 summarizes the main demographic and clinical characteristics of the 72 patients included. They were mostly men (88.9%), with a mean age of 68 (SD 7.9) years and FEV1 of 48.7 (SD 20.1)% of reference.
16S rRNA analysis
At phylum level, 13 different phyla were identified, six of them with median relative abundance (RA) above 0.1% (Additional file 1: Table S1). At genus level, 190 different genera were identified and, after removing the genera present in only one sample, 171 remained for subsequent analyses, 26 of them with RA above 0.1% (Table 2).
Alpha-diversity parameters showed a negative relationship with age (R2 = 0.075 p = 0.020 and R2 = 0.074 p = 0.020 respectively), but β-diversity analysis did not show significant differences in relation with this variable (p = 0.389).
We found a significant progressive increase in the RA of Pseudomonas genus and a decrease in the RA of Treponema in patients with more severe airflow limitation (Fig. 1). Regarding bacterial diversity, neither α-diversity parameters nor β-diversity analysis showed significant differences between GOLD grades of airflow limitation. Of note, airflow limitation severity was not related to age (p = 0.245).
Forty-nine COPD patients had not modified their inhaled maintenance treatment during the previous year; thirty-six of them (73.5%) used a combination of LAB/ICS, 9 (18.4%) were treated with LAB as monotherapy and 4 (8.2%) were not receiving COPD treatment. LAB/ICS treatment did not have any effect on either α-diversity (p = 0.365) or bacterial community composition in the patients studied (p = 0.963), when compared with patients not receiving this treatment. Similarly, the continuous use of LAB as monotherapy was not associated with significant changes in the respiratory microbiome (p = 0.854).
In the previous year, 31 patients (43.1%) did not report any acute episodes, 18 (25%) referred only one and 23 suffered two or more (31.9%), and were considered FE. Demographic and clinical characteristics of these three groups only showed statistically significant differences in lung function, with lower values in COPD patients reporting one or more exacerbations the previous year (Table 1) Comparisons between their respiratory microbiomes were made in pairs using patients without exacerbations as the reference. Patients with one exacerbation had significantly lower RA of the phylum TM7 (Additional file 1: Figure S1) and lower RAs of 13 different genera (Additional file 1: Figure S2). However, α-diversity parameters did not show significant differences between the groups, and β-diversity analysis did not demonstrate bacterial communities with a different composition (p = 0.081). FE showed a significant decrease in the RA of TM7 and Spirochaetes at phylum level (Additional file 1: Figure S3). At genus level, the RAs of Pseudomonas, Selenomonas and Anaerococcus increased, while 10 different genera decreased (Fig. 2). Alpha-diversity analysis did not show significant differences between groups, but β-diversity analysis demonstrated that the bacterial communities of COPD patients with frequent exacerbations differed significantly (p = 0.014).
Forty-two of the participants (58.3%) had ≥2% blood eosinophils. There were no significant differences in age (p = 0.368), sex (p = 1.00) and number of exacerbations the previous year (p = 0.080) between patients with ≥2% circulating eosinophils or less. The bacterial community in the former had significantly higher RAs of the phyla Bacteroidetes and Spirochaetes (Additional file 1: Figure S4). At genus level, 20 genera showed significantly higher RA and one genus, Peptostreptococcus, had lower RA in these patients (Fig. 3). Alpha-diversity was significantly higher in patients with ≥2% circulating eosinophils [Chao1 index: 224.51 (74.88) vs 277.39 (78.92), p = 0.006; and Shannon index: 3.94 (1.05) vs 4.54 (1.06), p = 0.020] (Fig. 4). Pearson's correlation coefficients were r = 0.282 (p = 0.016) for Chao1 and r = 0.231 (p = 0.051) for Shannon. β-diversity analysis showed a trend towards different bacterial communities (p = 0.072).
Multivariate analyses were performed with α-diversity as dependent variable and eosinophils levels as predictive factor, including age and lung function as covariates. Eosinophils in blood, expressed as percentage, kept a statistically significant relationship with Chao1 in this analysis (p = 0.026) and a borderline significance for Shannon (p = 0.051), a finding confirming that the bronchial microbiome was related to blood eosinophils independently of the functional limitations suffered by the patient.
To explore potential interactions between the previous history of exacerbations and eosinophils levels, we compared the microbiome in COPD patients with and without circulating eosinophils ≥2% stratified by the frequency of exacerbations. We found that the significant differences in the microbial composition related to patients with eosinophils ≥2% were maintained in the subsample of patients with no exacerbations or only one episode (p = 0.033), but this effect disappeared in FE (p = 0.995).
The main findings of this study were that the diversity and composition of the respiratory microbiome in clinically stable COPD patients change in relation to age, the severity of airflow limitation, exacerbation frequency and eosinophils in peripheral blood.
In our study, older age was significantly associated with a loss of diversity, which has been also found in the gut microbiome . Besides, in patients with severe asthma, an inverse correlation between α-diversity and age has been also reported . A previous work has shown less microbial diversity of the respiratory microbiome in younger COPD patients using bronchoalveolar lavage , but this sample targets the peripheral airway of the lung and it is not representative of the bronchial tree mainly sampled by sputum .
We found that patients with more severe airflow limitation had a significant decrease in the RA of Treponema and a progressive increase in the RA of Pseudomonas. These results suggest that severity-related changes in the respiratory microbiome are based on a decrease in specific genera, which are partially substituted by Pseudomonas. This change may be partly related to recurrent antibiotic exposure in previous years, considering the antibiotic sensitivity of the microorganisms part of Treponema genus. Previous cross-sectional studies evaluating the relation between bacterial diversity and more severe airway limitation have mostly showed a decline in advanced stages [26,27,28], associated with changes in the RAs of specific genera such as Haemophilus [28, 29]. These partly discordant results may be due to patient selection, considering that most of the previous studies have focused on a restricted number of patients with moderate or severe disease [26, 27] or an overrepresentation of patients with moderate disease  whereas we studied a wider range of disease severity (GOLD 1–4). Our results, therefore, support a significant role for Pseudomonas as the severity of the disease increases to higher lung function impairment.
We also found that the respiratory microbiome was significantly different in FE. Previous studies have investigated the characteristics of the respiratory microbiome during exacerbations [7, 30], and recently, like we do in this study, Mayhew and cols.  reported specific characteristics in the bronchial microbiome recovered from FE patients during clinical stability. Both studies show that FE have a different respiratory microbiome during clinical stability, suggesting that the microbial changes during exacerbations in FE may be a mixture of the dysbiosis found in stability and specific exacerbation-related perturbations of the lung bacteria community composition [28, 31].
Circulating eosinophils ≥2% were associated with higher microbial diversity in the population studied. Patients with ≥2% blood eosinophils have been reported to have more frequent exacerbations and a better response to ICS preventive therapy  Previous studies have demonstrated a different bronchial microbiome in eosinophilic COPD exacerbations [7, 31], which seems related to Th2 inflammation in both COPD and asthma . In our study we observed that in patients with ≥2% blood eosinophils higher microbial diversity is already present in stability, with an increase in the RA of 20 genera. Similar results have been reported in stable asthmatic patients, who showed a good correlation between the percentage of eosinophils in bronchoalveolar lavage and bacterial diversity . Higher bacterial diversity may have a protective role in patients with ≥2% blood eosinophils avoiding the presence of pathogenic bacteria such as Haemophilus influenzae and Streptococcus pneumonia which has been reported to be overrepresented in patients with eosinophils counts below 2% treated with ICS . Yet, when we stratified the COPD patients included according to both the level of circulating eosinophils and the frequency of exacerbations, we observed that the differences related to blood eosinophils disappeared in FE, likely highlighting a higher impact of frequent exacerbations on the respiratory microbiome in these patients.
This study has some potential limitations. First, we do not have a wide representation of the respiratory microbiome, which has been shown to be heterogeneous throughout the airway, because only sputum samples were analysed. Second, although the patients included had not taken antibiotics three months before their inclusion, we lack information on previous antibiotic treatments, which may have had an effect on their microbial communities. Finally, we analysed only bacterial communities, fungi and virus may also have an effect on these patients, either directly or through interactions with other microorganisms and the host.
This study shows that the respiratory microbiome in clinically stable COPD patients changes in relation to age, severity of airflow limitation, history of previous exacerbations and level of circulating eosinophils. These factors need to be considered when interpreting respiratory microbiome changes in patients with COPD.
Availability of data and materials
Bacterial 16SrRNA datasets from this study are accessible in the European Nucleotide Archive under the study PRJEB26773 with the sample numbers ERS2486515–609.
Chronic obstructive pulmonary disease
Forced expiratory volume the first second
Global Initiative for Chronic Obstructive Lung Disease
Linear discriminant analysis
Linear discriminant analysis Effect Size
Principal Coordinates Analysis
Polymerase chain reaction
Quantitative Insights Into Microbial Ecology
Ribosomal ribonucleic acid
Han MK, Agusti A, Calverley PM, Celli BR, Criner G, Curtis JL, et al. Chronic obstructive pulmonary disease phenotypes: the future of COPD. Am J Respir Crit Care Med. 2010;182:598–604.
Vestbo J. COPD: definition and phenotypes. Clin Chest Med. 2014;35:1–6.
Agustí A. Phenotypes and disease characterization in chronic obstructive pulmonary disease. Toward the extinction of phenotypes? Ann Am Thorac Soc. 2013;10(Suppl):S125–30.
Hilty M, Burke C, Pedro H, Cardenas P, Bush A, Bossley C, et al. Disordered microbial communities in asthmatic airways. PLoS One. 2010;5:e8578.
Faner R, Sibila O, Agustí A, Bernasconi E, Chalmers JD, Huffnagle GB, et al. The microbiome in respiratory medicine: current challenges and future perspectives. Eur Respir J. 2017;49(4):1602086.
Millares L, Ferrari R, Gallego M, Garcia-Nuñez M, Pérez-Brocal V, Espasa M, et al. Bronchial microbiome of severe COPD patients colonised by Pseudomonas aeruginosa. Eur J Clin Microbiol Infect Dis Off Publ Eur Soc Clin Microbiol. 2014;33:1101–11.
Wang Z, Bafadhel M, Haldar K, Spivak A, Mayhew D, Miller BE, et al. Lung microbiome dynamics in COPD exacerbations. Eur Respir J. 2016;47:1082–92.
From the Global Strategy for the Diagnosis, Management and Prevention of COPD, Global Initiative for Chronic Obstructive Lung Disease (GOLD) 2017. Available from: https://goldcopd.org.
Miravitlles M, Soler-Cataluña JJ, Calle M, Molina J, Almagro P, Quintano JA, et al. Spanish guidelines for Management of Chronic Obstructive Pulmonary Disease (GesEPOC) 2017. Pharmacological treatment of stable phase. Arch Bronconeumol. 2017;53:324–35.
Ho J, He W, Chan MTV, Tse G, Liu T, Wong SH, et al. Eosinophilia and clinical outcome of chronic obstructive pulmonary disease: a meta-analysis. Sci Rep. 2017;7:13451.
Vedel-Krogh S, Nielsen SF, Lange P, Vestbo J, Nordestgaard BG. Blood eosinophils and exacerbations in chronic obstructive pulmonary disease. The Copenhagen general population study. Am J Respir Crit Care Med. 2016;193:965–74.
Contoli M, Pauletti A, Rossi MR, Spanevello A, Casolari P, Marcellini A, et al. Long-term effects of inhaled corticosteroids on sputum bacterial and viral loads in COPD. Eur Respir J. 2017;50:1700451.
Standardization of Spirometry, 1994 Update. American Thoracic Society. Am J Respir Crit Care Med. 1995;152:1107–36.
Moermans C, Heinen V, Nguyen M, Henket M, Sele J, Manise M, et al. Local and systemic cellular inflammation and cytokine release in chronic obstructive pulmonary disease. Cytokine. 2011;56:298–304.
Rodriguez-Roisin R. Toward a consensus definition for COPD exacerbations. Chest. 2000;117(5 Suppl 2):398S–401S.
Anthonisen NR, Manfreda J, Warren CP, Hershfield ES, Harding GK, Nelson NA. Antibiotic therapy in exacerbations of chronic obstructive pulmonary disease. Ann Intern Med. 1987;106:196–204.
Murray PR, Washington JA. Microscopic and baceriologic analysis of expectorated sputum. Mayo Clin Proc. 1975;50:339–44.
Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.
Chao A. Nonparametric estimation of the number of classes in a population. Scand J Stat. 1984;11:265–70.
Shannon CE. The mathematical theory of communication. 1963. MD Comput Comput Med Pract. 1997;14:306–17.
Bray JR, Curtis JT. An ordination of the upland forest communities of southern Wisconsin. Ecol Monogr. 1957;27:325–49.
Maffei VJ, Kim S, Blanchard E, Luo M, Jazwinski SM, Taylor CM, et al. Biological aging and the human gut microbiota. J Gerontol A Biol Sci Med Sci. 2017;72:1474–82.
Taylor SL, Leong LEX, Choo JM, Wesselingh S, Yang IA, Upham JW, et al. Inflammatory phenotypes in patients with severe asthma are associated with distinct airway microbiology. J Allergy Clin Immunol. 2018;141:94–103.e15.
Pragman AA, Kim HB, Reilly CS, Wendt C, Isaacson RE. The lung microbiome in moderate and severe chronic obstructive pulmonary disease. PLoS One. 2012;7:e47305.
Cabrera-Rubio R, Garcia-Núñez M, Setó L, Antó JM, Moya A, Monsó E, et al. Microbiome diversity in the bronchial tracts of patients with chronic obstructive pulmonary disease. J Clin Microbiol. 2012;50:3562–8.
Garcia-Nuñez M, Millares L, Pomares X, Ferrari R, Pérez-Brocal V, Gallego M, et al. Severity-related changes of bronchial microbiome in chronic obstructive pulmonary disease. J Clin Microbiol. 2014;52:4217–23.
Galiana A, Aguirre E, Rodriguez JC, Mira A, Santibañez M, Candela I, et al. Sputum microbiota in moderate versus severe patients with COPD. Eur Respir J. 2014;43:1787–90.
Mayhew D, Devos N, Lambert C, Brown JR, Clarke SC, Kim VL, 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.
Sze MA, Dimitriu PA, Suzuki M, McDonough JE, Campbell JD, Brothers JF, et al. Host response to the lung microbiome in chronic obstructive pulmonary disease. Am J Respir Crit Care Med. 2015;192:438–45.
Bafadhel M, McKenna S, Terry S, Mistry V, Reid C, Haldar P, et al. Acute exacerbations of chronic obstructive pulmonary disease: identification of biologic clusters and their biomarkers. Am J Respir Crit Care Med. 2011;184:662–71.
Wang Z, Singh R, Miller BE, Tal-Singer R, Van Horn S, Tomsho L, et al. Sputum microbiome temporal variability and dysbiosis in chronic obstructive pulmonary disease exacerbations: an analysis of the COPDMAP study. Thorax. 2018;73:331–8.
Pascoe S, Locantore N, Dransfield MT, Barnes NC, Pavord ID. Blood eosinophil counts, exacerbations, and response to the addition of inhaled fluticasone furoate to vilanterol in patients with chronic obstructive pulmonary disease: a secondary analysis of data from two parallel randomised controlled trials. Lancet Respir Med. 2015;3:435–42.
Ghebre MA, Pang PH, Diver S, Desai D, Bafadhel M, Haldar K, et al. Biological exacerbation clusters demonstrate asthma and chronic obstructive pulmonary disease overlap with distinct mediator and microbiome profiles. J Allergy Clin Immunol. 2018;141:2027–2036.e12.
Denner DR, Sangwan N, Becker JB, Hogarth DK, Oldham J, Castillo J, et al. Corticosteroid therapy and airflow obstruction influence the bronchial microbiome, which is distinct from that of bronchoalveolar lavage in asthmatic airways. J Allergy Clin Immunol. 2016;137:1398–1405.e3.
Authors thank participants for their willingness to contribute to medical research and Michael Maudsley for providing an outline for this manuscript.
This study was funded by FIS PI15/00167 and PI15/02042 from Instituto de Salud Carlos III (ISCIII), SAF 2014–54371 (FEDER) from Ministerio de Economía y Competitividad, Sociedad Española de Neumología y Cirugía Torácica (SEPAR), Laboratorios Menarini S. A, Fundació Pla Armengol and Centro de Investigación Biomédica en Red de Enfermedades Respiratorias (CIBERES). CIBERES is an initiative of the ISCIII.
This work was also supported by grants from Agència de Gestió d’Ajuts Universitaris i de Recerca—Generalitat de Catalunya (2017 GRC-716).
Ethics approval and consent to participate
This is a cross-sectional, prospective, uncontrolled, multicentre, observational study. The study protocol was approved by the Ethics committees of the participating hospitals (IMIM-Hospital del Mar, Hospital Universitari Parc Taulí, Hospital Clinic, Hospital 12 Octubre, Fundación Jimenez Díaz and Hospital Son Espases), and all patients included signed their informed consent.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Relative abundances of the phyla detected. Figure S1. The TM7 phylum had significantly lower relative abundance in patients with one exacerbation than patients without exacerbations the previous year (dotted line = median). Figure S2. Thirteen genera with significantly lower relative abundances in COPD patients with one exacerbation the previous year compared to non-exacerbators. Figure S3. A significant reduction in the RA of phyla TM7 and Spirochaetes in patients with ≥2 exacerbations the previous year, using patients without exacerbations as the reference (dotted line = median). Figure S4. Phyla with significantly higher relative abundances in COPD patients showing circulating eosinophils ≥2%. (DOCX 725 kb)
About this article
Cite this article
Millares, L., Pascual, S., Montón, C. et al. Relationship between the respiratory microbiome and the severity of airflow limitation, history of exacerbations and circulating eosinophils in COPD patients. BMC Pulm Med 19, 112 (2019). https://doi.org/10.1186/s12890-019-0867-x