- Open Access
Microbial burden and viral exacerbations in a longitudinal multicenter COPD cohort
Respiratory Research volume 21, Article number: 77 (2020)
Chronic obstructive pulmonary disease (COPD) is a heterogeneous disease characterized by frequent exacerbation phenotypes independent of disease stage. Increasing evidence shows that the microbiota plays a role in disease progression and severity, but long-term and international multicenter assessment of the variations in viral and bacterial communities as drivers of exacerbations are lacking.
Two-hundred severe COPD patients from Europe and North America were followed longitudinally for 3 years. We performed nucleic acid detection for 20 respiratory viruses and 16S ribosomal RNA gene sequencing to evaluate the bacterial microbiota in 1179 sputum samples collected at stable, acute exacerbation and follow-up visits.
Similar viral and bacterial taxa were found in patients from the USA compared to Bulgaria and Czech Republic but their microbiome diversity was significantly different (P < 0.001) and did not impact exacerbation rates. Virus infection was strongly associated with exacerbation events (P < 5E-20). Human rhinovirus (13.1%), coronavirus (5.1%) and influenza virus (3.6%) constitute the top viral pathogens in triggering exacerbation. Moraxella and Haemophilus were 5-fold and 1.6-fold more likely to be the dominating microbiota during an exacerbation event. Presence of Proteobacteria such as Pseudomonas or Staphylococcus amongst others, were associated with exacerbation events (OR > 0.17; P < 0.02) but more strongly associated with exacerbation frequency (OR > 0.39; P < 4E-10), as confirmed by longitudinal variations and biotyping of the bacterial microbiota, and suggesting a role of the microbiota in sensitizing the lung.
This study highlights bacterial taxa in lung sensitization and viral triggers in COPD exacerbations. It provides a global overview of the diverse targets for drug development and explores new microbiome analysis methods to guide future patient management applications.
Chronic obstructive pulmonary disease (COPD) is defined by airflow limitation but encompasses several lung diseases. This heterogeneity includes differences in clinical characteristics, source of inflammation, response to therapies and causes of exacerbation . As COPD progresses, exacerbations become more frequent and more severe. Exacerbation rates reflect an independent susceptibility phenotype , which could be mediated by host factors , environmental factors , viral infections and/or the bacterial microbiome [5, 6].
Infections are predominant causes of COPD exacerbations, with approximately half reported to be caused by bacterial infections including non typeable Haemophilus influenzae (NTHi), Moraxella catarrhalis, Streptococcus pneumoniae, or Pseudomonas aeruginosa, and the other half by viral infections, primarily human Rhinovirus (HRV), but also Influenza virus, Coronavirus and Respiratory syncytial virus (RSV) to name a few [5, 6]. Bacteria and viruses are also frequently isolated in the airways of stable COPD patients [6,7,8]. The advent of culture-independent testing has suggested viral persistence  and colonization of the lower airways with a resident bacterial microbiota , implicating a role for the microbiota in disease pathogenesis, progression and treatment outcome of lung diseases .
Studies of the microbiome provide a new framework to understand host-pathogen interactions, which can also yield new markers for patient diagnosis and management. Microbiota diversity is seen as a potential biomarker in cases where a single pathogenic organism reduces community complexity such as in bacterial vaginosis , or Crohn’s disease . Lung microbial dysbiosis in COPD is characterized by decreased diversity [10,11,12], which may contribute to altered immune response to environmental insults . Dysbiosis at the time of COPD exacerbation contributes to increased disease severity  and higher 1-year mortality rates .
Geography could also be a potential covariate in COPD patients microbiota. The gut microbiota has been shown to be geographically variable . Previous studies in different conntries have evaluated the COPD microbiota [16, 17], but the effect of geographical variations has not yet been evaluated in a single study in relation to disease severity. The present cohort stems from a study on the incidence of viral infections in COPD . Patients were enrolled in Europe and North America and followed up for up to 3 RSV seasons, with scheduled wellness visits and unscheduled illness visits. To further our understanding of COPD exacerbation dynamics, we retrospectively evaluated the sputum bacterial microbiota from a subset of this study. The goals were to identify differences in patients with higher rates of exacerbations, to assess geographical differences in the microbiota between Europe and the USA, and to determine the influence of viral infections on microbiota diversity and the frequency of exacerbations.
Materials and methods
The patient cohort is part of an observational study on the incidence of acute respiratory illness (ARI) or events leading to the worsening of cardiorespiratory status in COPD (ClinicalTrials.gov, NCT01455402) . The protocol was approved by independent institutional review boards, and all subjects signed written informed consent at enrollment. The study population included 200 adults ≥50 years of age with COPD, recruited at 6 sites in Bulgaria, 5 sites in Czech Republic, and 14 sites across the USA from fall 2011 to spring 2014 (Fig. 1a). Cohort from Bulgaria and Czech Republic showed similar characteristics, and are analyzed jointly as Europe for ease of representation. Subjects had scheduled wellness visits in May and October each year to obtain sputum and clinical data (Fig. 1a). Unscheduled illness visits to collect sputum and clinical data were performed when a subject experienced an ARI or acute exacerbation of COPD symptoms, and during follow-up illness visits. Samples were considered stable if collected at least 30 days from the last day of hospitalization or from the last ARI or acute exacerbation event if it did not require hospitalization. Acute exacerbation samples were collected within 72 h of an event, additional samples were collected 4–30 days after an acute exacerbation event during unscheduled visits (Fig. 1a, c). Patients were characterized with a frequent exacerbator phenotype if they experienced 2 or more exacerbation events per year. No investigational drug was administered in the study. The subject’s physician prescribed and recorded all treatment deemed necessary to provide adequate supportive care. Samples collected 1 to 7 days following treatment were considered treatment associated (Fig. 1d).
Sputum sample collection
Spontaneous rather than induced sputum collection was possible in this severe COPD study subset following a standardized collection visit . Subjects were asked to gargle with water immediately prior to sputum collection to reduce the number of oral bacteria . Subjects were asked to cough deeply and expectorate into a cup, mixed 1:1 with cold Transport Media and kept at -60 °C or below .
The GenMark respiratory virus panel (GenMark Diagnostics, Inc. Carlsbad, CA) was used to detect common respiratory viruses from all sputum samples. RSV detection was confirmed by RT-PCR as previously described . Additionally, RT-PCR using primers against HRV VP1-VP4  was used to detect human rhinovirus subtypes.
16S rRNA gene sequencing
Bacterial genomic DNA was extracted at a single central lab using the Zymobiomics 96 DNA kit (Zymo, California, USA) following manufacturer’s instructions. The V4 hypervariable region of the 16S rRNA gene was PCR amplified  and sequenced using Illumina Miseq platform along with negative controls, Zymobiomics Microbial community and DNA standards, and a PhiX library for quality control.
16S rRNA gene sequences were analysed using QIIME2 version 2018.11 . DADA2 software package , wrapped in QIIME2, was used for correcting sequences and obtained 9853 annotated sequence variants (ASVs). Taxonomy was assigned using a Naïve Bayes classifier [25, 26] that was trained on the Greengenes database version 13.8 clustered at 99% identity  V4 sequences. Alignment was performed with MAFFT , masked and used in FastTree  to build the phylogenetic tree. Alpha-diversity metrics , beta diversity metrics , and Principal coordinate analysis (PCoA) were estimated after samples were rarefied to 1000 sequences per sample. Significant features of interest were re-tested using non-rarefying alpha  and beta diversity estimates . Analysis of composition (ANCOM) was used for differential taxa abundance calculations  on non-rarefied data. ANCOM account for the structure of the data and controls for the false discovery rate. Microbiota profiles were clustered into biotypes using the biotypeR package in R . Longitudinal microbiota variations at stable time points were assessed by calculating the median weighted UniFrac distance  for patients with at least 3 stable samples. Patients were categorized into the bottom and top quartile of within-patient stable samples median UniFrac distances, respectively.
Associations of the microbiota and viral components with demographics and clinical data were calculated in R using Fisher’s test for categorical variables, Welsh’s T-test for comparison of means in continuous variables, and chi-squared test to compare expected frequencies. Odds ratio were calculated using questionr package in R . Associations of microbial diversity with demographics and clinical data were calculated with Qiime2’s diversity plugin using grouped and pairwise Kruskal-Wallis test corrected for false discovery rate for analysis of alpha diversity, and PERMANOVA following 999 permutations for analysis of beta diversity distances.
The cohort was composed of 200 patients, 101 from Europe (Bulgaria and Czech Republic) and 99 from the USA (Fig. 1 and Table 1). All but 5 patients presented with severe or very severe COPD. Patients from Europe and the USA were matched by age, sex, and GOLD stage. At enrollment, their forced expiratory volume as a percent of predicted (FEV1%) and their comorbidities (congestive heart failure, diabetes, hypertension and malignancy) were similar between Europe and the USA. The USA patients had significantly higher smoking history in pack-years and longer COPD duration. The number of exacerbations per year and the number of frequent exacerbator phenotype, defined as 2 or more exacerbations/year, was higher in the USA compared to Europe, but not statistically significant (Table 1). A total of 1179 sputum samples were collected from these patients over 3 years, with similar proportions of samples collected at acute exacerbation in Europe (33.3%) and in the USA (33.9%) (Supplementary Table 1). Approximately 11 and 16% of samples were associated with antibiotic and corticosteroid treatment, respectively. Significantly more samples collected in the USA were associated with antibiotic treatment (Supplementary Table 1).
Overview of 16S microbiota and viral detection in COPD sputum
Analysis of bacterial taxa in 1179 sputum samples shows phyla commonly observed in the lung microbiota, with Firmicutes, Proteobacteria and Bacteroidetes representing a majority (> 80%) of the phyla identified (Fig. 2a). Prevotella, Veillonella, Streptococcus and Haemophilus represented the most prevalent (> 65%) bacterial genera (Fig. 2a). Samples from European patients had more Bacteroidetes and less Proteobacteria overall than samples from the USA patients. Within the phylum Firmicutes, the USA patients had more Streptococcus than Veillonella compared to European patients (Fig. 2a). Microbiota predominant with Prevotella, Streptococcus and Veillonella were found in a majority of samples (1012/1179; 85.8%) (Fig. 2b). Prevotella, Streptococcus and Veillonella represented the majority of the microbiota in samples collected at stable states, Haemophilus and Moraxella were predominant in acute exacerbation samples, and Pseudomonas was prevalent in exacerbation follow-up samples (Fig. 2b).
Viral testing on 1179 sputum samples showed that 25.6% of sputum samples were positive for at least one virus, with HRV found in the largest proportion of samples (13.6%), followed by Coronavirus (5.1%) and Influenza virus (3.6%) (Fig. 2c). All viruses, except for adenoviruses, were detected more frequently at acute exacerbation (14.4%), than at follow-up visits (7.2%) or stable (5.6%) (Fig. 2c). There were no significant differences in viral incidence between Europe and the USA, with the exception of Coronavirus HKU1, Influenza B virus and RSV A (Supplementary Table 1).
Alpha diversity metrics such as the number of observed ASV, Shannon evenness, or Faith phylogenetic diversity (PD) indexes, represent the mean number of taxons in a sample, sometimes weighted by the phylogenetic relatedness. Overall, samples from the USA patients had significantly lower shannon diversity index than those from Europe (Fig. 2a), even when corrected for antibiotic use or clinical sites (Supplementary Fig. 1 and 2). Years of severe COPD and smoking pack-years were associated with differences in microbiota diversity but not as significantly as geography (Supplementary Table 2). No significant differences in the number of observed ASVs and in Faith PD index were observed at exacerbation compared to stable samples (Supplementary Table 2). Prevotella-predominant microbiota showed the highest alpha diversity, while samples predominant with Escherichia, Haemophilus, and Pseudomonas had the lowest alpha diversity (Fig. 2b). There were no significant differences in microbial diversity between samples infected or not with a virus, or between samples infected with different viruses (Fig. 2c). Beta diversity are metrics such as weighted UniFrac and robust Aitchison PCA used for comparing microbiota communities resulting in distance matrices. Here, principal coordinate visualization of weighted UniFrac distance supports the above conclusions, showing some differences in geography and predominant bacteria but not affected by the type of viral infection (Supplementary Fig. 3). These results were also supported by testing non-rarefied alpha and beta-diversity estimates, with geography (P < 0.044) and dominant bacterial genus (P < 0.001), but not viruses (P > 0.094) being associated with significant changes in diversity (data not shown).
Odds ratio of exacerbation
The odds ratio of an acute exacerbation event and frequent exacerbations (≥2 events/ year) was calculated for demographic and clinical data, viral infections and abundance of certain bacterial taxa in the lung microbiota (Fig. 3). Viral infections were more strongly associated with an exacerbation event than with frequent exacerbations. Parainfluenzaviruses (PIV), Influenza B virus and RSV B had the highest odds ratio of an exacerbation event (Fig. 3a). Interestingly, Influenza B virus was negatively correlated with frequent exacerbations, as it was only detected in patients that exacerbated infrequently.
Bacteria were more strongly associated with exacerbation frequency than with exacerbation events. Presence or higher abundance of Enterococcus, Lactobacillus, Moraxella, Pseudomonas, Staphylococcus and Streptococcus was correlated with frequent exacerbations. Neisseria, Prevotella and Veillonella were significantly associated with a lower exacerbation frequency (Fig. 3b). Interestingly, top quartile abundance odds ratio of the genus Haemophilus was not associated with higher exacerbation rate. This effect, however, seemed mediated by H.parainfluenzae. We noted that differential abundance of only 2 bacterial taxa were significantly associated with exacerbation events, while 19 taxa were associated with exacerbation frequency (Supplementary Fig. 4), with high E.coli, Lactobacillus, and Staphylococcus in stable samples as potential predictors of frequent exacerbation (Supplementary Fig. 5).
Patients with hypertension had a significantly higher odds ratio of being frequent COPD exacerbators (Fig. 3c). Other comorbidities did not significantly influence exacerbation frequency.
Longitudinal variations of stable state sputum microbiota and viruses
We quantified temporal variability of the sputum microbiota at stable state within individual subjects. For 84 patients with more than 3 longitudinal stable samples, we calculated their median weighted UniFrac distances and categorize the top and bottom quartile patients into consistent and variable microbiota over time (Fig. 4a). Patients with a more variable sputum microbiota (median weighted UniFrac > 0.22) were more likely to have a higher relative abundance of Bacillus, Escherichia, Lactobacillus, Moraxella, and Staphylococcus (Fig. 4 and Supplementary Fig. 6). Microbiota variability in seemingly stable disease state were associated with higher exacerbation frequency and frequent viral infections (Fig. 4b).
Longitudinal sampling also enabled the assessment of recurrent viral infections by the same viral species, strain or subtype in consecutive samples. Viruses were detected in 2 to 3 consecutive samples from 14 patients (Fig. 5). The time beween 2 consecutive samples ranged from 3 days to 1 year. Viruses that were detected within 2–3 weeks of sampling corresponded most probably to a typical single infection, and these included the detection of Coronavirus OC43, HRV A, PIV4, RSV A and B. Viruses that were detected in the same patient over 60 days apart may correspond to chronic or recurrent infections, with Adenovirus C, Coronavirus HKU1, HRV A, HRV B and HRV C detected. Eighteen out of 30 (60%) samples from these 14 patients corresponded to acute exacerbation events, and 6/14 (42%) patients were frequent exacerbators (Fig. 5).
Biotyping of COPD sputum microbiota
We sought to define common microbiota clusters and their association with clinical characteristics. Microbiota in stable samples could be separated into 2 biotypes, as indicated by the highest Calinski-Harabasz (CH) index following iterative partitioning-around-medoids clustering analysis over the Jensen-Shannon distance. Between-class analysis showed 2 major clusters, biotype 1 represented by Neisseria and Veillonella, and biotype 2 represented by Streptococcus and Rothia (Fig. 6a). Samples with high relative abundance of Streptococcus/Rothia (biotype 2) were found in greater proportion in USA patients, which was associated with longer history of COPD and less frequent detection of viruses at stable visits (Fig. 6a). Samples at acute exacerbation visits could be separated into 3 biotypes, characterized by a high relative abundance of either Prevotella (biotype 1), Streptococcus (biotype 2), or Haemophilus (biotype 3) (Fig. 6b). Streptococcus and Haemophilus were found in a majority of USA samples, and associated with longer COPD duration, higher exacerbation frequency, antibiotics and corticosteroid use, but did not significantly correlate with higher viral infections (Fig. 6b). Exacerbation follow-up samples were more diverse and could be optimally clustered into 6 biotypes. Biotype 6 was characterized by a high relative abundance of Pseudomonas and significantly associated with longer COPD duration and antibiotic use. Biotype 4 was characterized by a high relative abundance of Streptococcus or Rothia and was significantly associated with higher exacerbation frequency. Other biotypes were not significantly associated with clinical characteristics (Fig. 6c). Biotyping shows that microbiota profile diversity is dynamically dependent on disease state.
Understanding of the presence and role of both bacterial and viral pathogens over time in the heterogeneous and dynamic COPD disease  is needed for patient treatment and management. The characterization of the 16S rRNA gene microbiota and respiratory viruses from a longitudinal and international cohort of severe COPD patients described in this study provides the largest survey to date on their complex associations with geography, exacerbation frequency and other demographic and clinical history.
COPD patients are particularly susceptible to respiratory infections . HRV was identified as the most prominent agent in respiratory tract infections in this cohort of COPD patients. Viral characterization is of particular importance as few reports exist on the diversity of respiratory viral agents in COPD  and viral diversity should be accounted for when designing new treatments. In particular, the newly described HRV C was detected in 23 of 11,179 samples (21 of 200 patients) in the present cohort, which was reported only once previously .
We detected viruses in 15% of stable samples. Asymptomatic viral infections in COPD patients are common . However, the role of these asymptomatic infections in disease progression is unclear. Most respiratory viruses tested here showed highly significant positive odds ratio with exacerbation events, but lower significance in regards to exacerbation frequency. This may indicate that viral infections alone do not sensitize the lung to repeat exacerbations as much as the bacterial microbiota. Interestingly, similar trends have been shown in asthma where respiratory viral infections in early life resulted in microbiome changes and hypersensitivity predisposition that can lead to asthma [41, 42].
Repeat viral detection were more frequent in patients with frequent exacerbator phenotype, but the number was small, and contradicting reports exist on the repeated detection of a single virus in COPD . Further complete viral genomic characterization will be needed to understand the nature of viral infection.
With the advent of culture-independent techniques, it appears that all microbiomes harbor potential bacterial pathogens as characterized here and elsewhere , but that only a portion of them will develop exacerbation-prone phenotypes, a phenotype that appears independent of disease GOLD stage yet linked to microbiota diversity . It was previously observed that microbiota diversity in the COPD lung correlated with disease severity but not disease state . Here, microbiota diversity alone was not correlated with frequent exacerbations, but was highly correlated with certain bacterial taxa dominating the microbiota. Microbiota predominant with Escherichia, Pseudomonas or Streptococcus, showed significantly lower alpha diversity and significant positive odds ratio with the frequent exacerbation phenotype, suggesting a role of the microbiota in sensitizing the COPD lung to acute exacerbations. This study was limited in disease severity metrics with the exception of baseline evaluation and sampling of events. Longitudinal monitoring of symptoms scales would help to better understand the relation of certain bacteria to not only exacerbation frequency, but also the symptom severity and COPD progression.
The use of biotyping has been seldom used in respiratory microbiota research  and not yet explored in COPD. A complex resident bacterial community could be identified in all COPD sputum samples and categorized into 2, 3 and 6 biotypes at stable, acute exacerbation and exacerbation follow-up visits, respectively. Biotype 1 at stable state was associated with higher viral infections, while biotypes 2 and 3 at acute exacerbation were associated with high exacerbation frequency. These findings are interesting because they mirror another study showing the partitioning of COPD exacerbation samples into 3 cytokine profile clusters  with associations to specific ratios of Proteobacteria, Firmicutes and Bacteroidetes that highlighted the heterogeneity of exacerbation profiles in COPD patients. During exacerbation follow-up visits, biotype 6 with a high relative abundance of Pseudomonas was found over-represented in samples associated with antibiotics use. Antibiotic treatment inadequacy is the cause for secondary infection or the emergence of multi-drug resistant P.aeruginosa . New targeted treatments, such as monoclonal antibodies, could be useful in such settings .
The principal novelty of this study cohort was the long term patient follow-up. We were able to collect several sputum samples per patient at stable state over the course of 3 years and studied the COPD microbiota longitudinally. The lung microbiome is inherently variable, shaped by a process of inhalation and elimination . The lung microbiome is also personal, with large inter-patient variability . Previously, it was shown that microbial dysbiosis from stable to exacerbated state correlated with greater exacerbation severity . Here, patients with greater microbiota variability at stable state correlated with higher exacerbation frequency. Proteobacteria such as Pseudomonas and Moraxella were more abundant in patients with more variable microbiota at stable state. Interestingly, P.aeruginosa and M.cattharalis are prominent causes of exacerbations , but their role in stable disease is less clear . It was previously shown that chronic colonization with P. aeruginosa occurs more frequently in more severe COPD patients  and that M.catarrhalis asymptomatic colonization was associated with a greater frequency of a sputum IgA response than exacerbation . Our results suggest that dysbiotic burden at stable state by Pseudomonas, Moraxella and others might sensitize the lung to further exacerbations and viral infections. Pseudomonas and Moraxella, like many opportunistic Proteobacteria, are pro-inflammatory [52, 53]. Imbalanced inflammation can improve P.aeruginosa’s fitness , allow the acquisition of new M.cattharalis strains , leading to exacerbation and possibly infections from other pathogens in a coupled cycle of inflammation and dysbiosis . Microbiology clinical testing in COPD patients is most often performed at exacerbation or follow-up visits. Patients might benefit from clinical monitoring of these bacteria at stable state to assess their presence and/or growth which could lead to potential future exacerbations.
We also noted geographical differences in COPD lung microbiota. Geographical differences in gut microbiota have previously been noted , but not yet in the lung. There were significant differences in alpha and beta diversity between the USA and Europe, but not within countries or sites. Microbiota diversity in the USA was lower and although frequent exacerbator phenotypes were more common than in Europe, the difference was not significant. USA patients tended to have samples with high relative abundance of Streptococcus (biotype 2) and Haemophilus (biotype 3) associated with the frequent exacerbator phenotype. S.pneumoniae and H.influenzae are commonly associated with exacerbations , and should also be considered as potential risk factors in the frequent exacerbator phenotype. Significantly more samples collected in the USA were associated with antibiotics use, but this alone did not explain differences in diversity. This observational study included a variety of standard-of-care medications, doses and timings precluding precise treatment effect modeling on the microbiota. Clinical trials exploring current and novel treatment modalities will lead to better patient management and antibiotics stewardship as reviewed elsewhere .
Predominance of Haemophilus was over-represented in acute exacerbation samples (Figs. 2 and 6), as noted in previous studies [6, 14, 17, 43, 56]. However, interestingly, using odds ratio (Fig. 3) or ANCOM (Suppplementary Fig. 4a) over Haemophilus abundance, H.influenzae was high but not significantly associated with acute exacerbation event, whereas H.parainfluenzae was significantly higher at stable state. Differences in Haemophilus abundance compared to other studies might be due to the type of sputum collected, transportation media, extraction protocol or an effect of the sample size. Although H.parainfluenzae can be the cause of respiratory infection in healthy subjects, it has not been associated with exacerbation in COPD . H.parainfluenzae could compete for niche resources leading to overgrowth of the more pathogenic H.influenzae. Previous work has shown that H.influenzae competes with S.pneumoniae  and that patients colonized by NTHi and acquiring HRV have more frequent and severe exacerbations . Here, Streptococcus species could not be resolved using 16S rRNA V4 region, and although speciation of Haemophilus was attempted, further validation using targeted PCR or whole genome sequencing will be necessary to ensure correct discrimination. To note, constant improvements in 16S databases can also affects taxonomic resolution. SILVA database  version 132 updated in 2017 classified reads into more genera (n = 562) compared to Greengenes version 13.8 updated in 2013 (n = 395). However, bacterial taxa discussed in this study showed less than 1% variations in read classification between the 2 database classifiers (data not shown), and conclusions were unchanged. Speciations and typing of bacteria and viruses are critical to understanding their pathogenicity and complex relationships. Greater taxonomic resolution will be achieved using updated databases and more comprehensive techniques like shotgun metagenomics.
In summary, our study provides a broad survey of viruses and bacteria colonizing severe COPD patients, providing clinicians with potential targets for clinical testing and patient treatment. It demonstrates that viral infections are strongly associated with acute exacerbation events, and that particular components of the microbiota are associated with higher exacerbation frequency. Geographic and longitudinal differences in the lung COPD microbiota exist and were correlated with exacerbation outcomes. Stable state longitudinal microbiota monitoring and biotyping could lead to the identification of potential biomarkers indicative of future exacerbations from bacterial sensitization. Comprehensive microbiota profiling and respiratory viral detection will be useful in the development of anti-microbial agents for therapeutic intervention or for better patient management.
Availability of data and materials
Data underlying the findings described in this manuscript may be obtained in accordance with AstraZeneca’s data sharing policy described at https://astrazenecagrouptrials.pharmacm.com/ST/Submission/Disclosure
Analysis of composition
Acute respiratory illness
Annotated sequence variants
Chronic obstructive pulmonary disease
Forced expiratory volume in 1 s as a percent of predicted
Global initiative for chronic obstructive lung disease
Non typeable Haemophilus influenzae
Respiratory syncytial virus
Agusti A, Calverley PMA, Celli B, Coxson HO, Edwards LD, Lomas DA, MacNee W, Miller BE, Rennard S, Silverman EK, Tal-Singer R, Wouters E, Yates JC, Vestbo J. Evaluation of COPD longitudinally to identify predictive surrogate endpoints (ECLIPSE) investigators. Characterisation of COPD heterogeneity in the ECLIPSE cohort. Respir Res. 2010;11:122.
Hurst JR, Vestbo J, Anzueto A, Locantore NW, Müllerova H, Tal-Singer R, Miller BE, Lomas DA, Agusti A, MacNee W, Calverley PMA, Rennard SI, Wouters EFM, Wedzicha JA. Susceptibility to exacerbation in chronic obstructive pulmonary disease. N Engl J Med. 2010;363:1128–38.
Ghebre MA, Pang PH, Diver S, Desai D, Bafadhel M, Haldar K, Kebadze T, Cohen S, Newbold P, Rapley L, Woods J, Rugman P, Pavord ID, Johnston SL, Barer M, May RD, Brightling CE. 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.
Lee J, Jung HM, Kim SK, Yoo KH, Jung K-S, Lee SH, Rhee CK. Factors associated with chronic obstructive pulmonary disease exacerbation, based on big data analysis. Sci Rep. 2019;9(1):6679.
Falsey AR, Walsh EE, Esser MT, Shoemaker K, Yu L, Griffin MP. Respiratory syncytial virus-associated illness in adults with advanced chronic obstructive pulmonary disease and/or congestive heart failure. J Med Virol. 2019;91:65–71.
Sethi S, Murphy TF. Infection in the pathogenesis and course of chronic obstructive pulmonary disease. N Engl J Med. 2008;359:2355–65.
Leung JM, Tiew PY, Mac Aogáin M, Budden KF, Yong VFL, Thomas SS, Pethe K, Hansbro PM, Chotirmall SH. The role of acute and chronic respiratory colonization and infections in the pathogenesis of COPD. Respirol Carlton Vic. 2017;22:634–50.
Wilkinson TMA, Donaldson GC, Johnston SL, Openshaw PJM, Wedzicha JA. Respiratory syncytial virus, airway inflammation, and FEV1 decline in patients with chronic obstructive pulmonary disease. Am J Respir Crit Care Med. 2006;173:871–6.
Sikkel MB, Quint JK, Mallia P, Wedzicha JA, Johnston SL. Respiratory syncytial virus persistence in chronic obstructive pulmonary disease. Pediatr Infect Dis J. 2008;27:S63–70.
Mammen MJ, Sethi S. COPD and the microbiome. Respirol Carlton Vic. 2016;21:590–9.
Oakley BB, Fiedler TL, Marrazzo JM, Fredricks DN. Diversity of human vaginal bacterial communities and associations with clinically defined bacterial vaginosis. Appl Environ Microbiol. 2008;74:4898–909.
Manichanh C, Rigottier-Gois L, Bonnaud E, Gloux K, Pelletier E, Frangeul L, Nalin R, Jarrin C, Chardon P, Marteau P, Roca J, Dore J. Reduced diversity of faecal microbiota in Crohn’s disease revealed by a metagenomic approach. Gut. 2006;55:205–11.
Yatsunenko T, Rey FE, Manary MJ, Trehan I, Dominguez-Bello MG, Contreras M, Magris M, Hidalgo G, Baldassano RN, Anokhin AP, Heath AC, Warner B, Reeder J, Kuczynski J, Caporaso JG, Lozupone CA, Lauber C, Clemente JC, Knights D, Knight R, Gordon JI. Human gut microbiome viewed across age and geography. Nature. 2012;486:222–7.
Wang Z, Singh R, Miller BE, Tal-Singer R, Van Horn S, Tomsho L, Mackay A, Allinson JP, Webb AJ, Brookes AJ, George LM, Barker B, Kolsum U, Donnelly LE, Belchamber K, Barnes PJ, Singh D, Brightling CE, Donaldson GC, Wedzicha JA, Brown JR, COPDMAP. Sputum microbiome temporal variability and dysbiosis in chronic obstructive pulmonary disease exacerbations: an analysis of the COPDMAP study. Thorax. 2018;73:331–8.
Leitao Filho FS, Alotaibi NM, Ngan D, Tam S, Yang J, Hollander Z, Chen V, FitzGerald JM, Nislow C, Leung JM, Man SFP, Sin DD. Sputum microbiome is associated with 1-year mortality after chronic obstructive pulmonary disease hospitalizations. Am J Respir Crit Care Med. 2019;199:1205–13.
Pragman AA, Gould TJ, Knutson K, Hodgson S, Isaacson RE, Reilly CS, Wendt CH. The chronic obstructive pulmonary disease frequent Exacerbator phenotype is associated with decreased upper airway microbiota α diversity. Ann Am Thorac Soc. 2018;15:S286.
Wang Z, Bafadhel M, Haldar K, Spivak A, Mayhew D, Miller BE, Tal-Singer R, Johnston SL, Ramsheh MY, Barer MR, Brightling CE, Brown JR. Lung microbiome dynamics in COPD exacerbations. Eur Respir J. 2016;47:1082–92.
Tangedal S, Aanerud M, Persson LJP, Brokstad KA, Bakke PS, Eagan TM. Comparison of inflammatory markers in induced and spontaneous sputum in a cohort of COPD patients. Respir Res. 2014;15:138.
Muzanye G, Morgan K, Johnson J, Mayanja-Kizza H. Impact of mouth rinsing before sputum collection on culture contamination. Afr Health Sci. 2009;9:200.
Kimball AM, Foy HM, Cooney MK, Allan ID, Matlock M, Plorde JJ. Isolation of respiratory syncytial and influenza viruses from the sputum of patients hospitalized with pneumonia. J Infect Dis. 1983;147:181–4.
McIntyre CL, Knowles NJ, Simmonds P. Proposals for the classification of human rhinovirus species a, B and C into genotypically assigned types. J Gen Virol. 2013;94:1791–806.
Kozich JJ, Westcott SL, Baxter NT, Highlander SK, Schloss PD. Development of a dual-index sequencing strategy and Curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform. Appl Environ Microbiol. 2013;79:5112–20.
Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet C, Al-Ghalith GA, Alexander H, Alm EJ, Arumugam M, Asnicar F, Bai Y, Bisanz JE, Bittinger K, Brejnrod A, Brislawn CJ, Brown CT, Callahan BJ, Caraballo-Rodríguez AM, Chase J, Cope E, Silva RD, Dorrestein PC, Douglas GM, Durall DM, Duvallet C, Edwardson CF, Ernst M, Estaki M, Fouquier J, et al. QIIME 2: reproducible, interactive, scalable, and extensible microbiome data science: PeerJ Inc.; 2018. https://doi.org/10.7287/peerj.preprints.27295v2.
Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2: high-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13:581–3.
Werner JJ, Koren O, Hugenholtz P, DeSantis TZ, Walters WA, Caporaso JG, Angenent LT, Knight R, Ley RE. Impact of training sets on classification of high-throughput bacterial 16s rRNA gene surveys. ISME J. 2012;6:94–103.
Bokulich NA, Kaehler BD, Rideout JR, Dillon M, Bolyen E, Knight R, Huttley GA, Gregory CJ. Optimizing taxonomic classification of marker-gene amplicon sequences with QIIME 2’s q2-feature-classifier plugin. Microbiome. 2018;6:90.
McDonald D, Price MN, Goodrich J, Nawrocki EP, DeSantis TZ, Probst A, Andersen GL, Knight R, Hugenholtz P. An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J. 2012;6:610–8.
Katoh K, Misawa K, Kuma K, Miyata T. MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002;30:3059–66.
Price MN, Dehal PS, Arkin AP. FastTree: computing large minimum evolution trees with profiles instead of a distance matrix. Mol Biol Evol. 2009;26:1641–50.
Faith DP. Conservation evaluation and phylogenetic diversity. Biol Conserv. 1992;61:1–10.
Lozupone CA, Hamady M, Kelley ST, Knight R. Quantitative and qualitative beta diversity measures lead to different insights into factors that structure microbial communities. Appl Environ Microbiol. 2007;73:1576–85.
Willis A, Bunge J. Estimating diversity via frequency ratios. Biometrics. 2015;71:1042–9.
Martino C, Morton JT, Marotz CA, Thompson LR, Tripathi A, Knight R, Zengler K. A novel sparse compositional technique reveals microbial perturbations. mSystems. 2019;4(1).
Mandal S, Van Treuren W, White RA, Eggesbø M, Knight R, Peddada SD. Analysis of composition of microbiomes: a novel method for studying microbial composition. Microb Ecol Health Dis. 2015;26:27663.
Arumugam M, Raes J, Pelletier E, Le Paslier D, Yamada T, Mende DR, Fernandes GR, Tap J, Bruls T, Batto J-M, Bertalan M, Borruel N, Casellas F, Fernandez L, Gautier L, Hansen T, Hattori M, Hayashi T, Kleerebezem M, Kurokawa K, Leclerc M, Levenez F, Manichanh C, Nielsen HB, Nielsen T, Pons N, Poulain J, Qin J, Sicheritz-Ponten T, et al. Enterotypes of the human gut microbiome. Nature. 2011;473:174–80.
Flores GE, Caporaso JG, Henley JB, Rideout JR, Domogala D, Chase J, Leff JW, Vázquez-Baeza Y, Gonzalez A, Knight R, Dunn RR, Fierer N. Temporal variability is a personalized feature of the human microbiome. Genome Biol. 2014;15(12):531.
Barnier J, Briatte F, Larmarange J. questionr: Functions to Make Surveys Processing Easier. 2018. at <https://CRAN.R-project.org/package=questionr>.
Faner R, Agustí Á. Multilevel, dynamic chronic obstructive pulmonary disease heterogeneity. A challenge for personalized medicine. Ann Am Thorac Soc. 2016;13:S466–70.
Lau SKP, Yip CCY, Lin AWC, Lee RA, So L-Y, Lau Y-L, Chan K-H, Woo PCY, Yuen K-Y. Clinical and molecular epidemiology of human rhinovirus C in children and adults in Hong Kong reveals a possible distinct human rhinovirus C subgroup. J Infect Dis. 2009;200:1096–103.
Gandhi A, Walsh EE, Formica MA, Hennessey PA, Criddle MM, Peterson DR, Baran A, Falsey AR. Factors associated with symptomatic rhinovirus infection in patients with COPD. J Clin Virol Off Publ Pan Am Soc Clin Virol. 2012;55:343–7.
Kusel MMH, de Klerk NH, Kebadze T, Vohma V, Holt PG, Johnston SL, Sly PD. Early-life respiratory viral infections, atopic sensitization, and risk of subsequent development of persistent asthma. J Allergy Clin Immunol. 2007;119:1105–10.
Dinwiddie DL, Denson JL, Kennedy JL. Role of the airway microbiome in respiratory infections and asthma in children. Pediatr Allergy Immunol Pulmonol. 2018;31:236–40.
Mayhew D, Devos N, Lambert C, Brown JR, Clarke SC, Kim VL, Magid-Slav M, Miller BE, Ostridge KK, Patel R, Sathe G, Simola DF, Staples KJ, Sung R, Tal-Singer R, Tuck AC, Horn SV, Weynants V, Williams NP, Devaster J-M, Wilkinson TMA. Longitudinal profiling of the lung microbiome in the AERIS study demonstrates repeatability of bacterial and eosinophilic COPD exacerbations. Thorax. 2018;73:422–30.
Kelly MS, Surette MG, Smieja M, Pernica JM, Rossi L, Luinstra K, Steenhoff AP, Feemster KA, Goldfarb DM, Arscott-Mills T, Boiditswe S, Rulaganyang I, Muthoga C, Gaofiwe L, Mazhani T, Rawls JF, Cunningham CK, Shah SS, Seed PC. The nasopharyngeal microbiota of children with respiratory infections in Botswana. Pediatr Infect Dis J. 2017;36:e211–8.
Eklöf J, Gliese KM, Ingebrigtsen TS, Bodtger U, Jensen J-US. Antibiotic treatment adequacy and death among patients with Pseudomonas aeruginosa airway infection. PLoS One. 2019;14:e0226935.
Le HN, Tran VG, Vu TTT, Gras E, Le VTM, Pinheiro MG, Aguiar-Alves F, Schneider-Smith E, Carter HC, Sellman BR, Stover CK, DiGiandomenico A, Diep BA. Treatment Efficacy of MEDI3902 in Pseudomonas aeruginosa blood stream infection and acute pneumonia rabbit models. Antimicrob Agents Chemother. 2019. https://doi.org/10.1128/AAC.00710-19.
Dickson RP, Martinez FJ, Huffnagle GB. The role of the microbiome in exacerbations of chronic lung diseases. Lancet Lond Engl. 2014;384:691–702.
Erb-Downward JR, Thompson DL, Han MK, Freeman CM, McCloskey L, Schmidt LA, Young VB, Toews GB, Curtis JL, Sundaram B, Martinez FJ, Huffnagle GB. Analysis of the lung microbiome in the “healthy” smoker and in COPD. PLoS One. 2011. https://doi.org/10.1371/journal.pone.0016384.
Beasley V, Joshi PV, Singanayagam A, Molyneaux PL, Johnston SL, Mallia P. Lung microbiology and exacerbations in COPD. Int J Chron Obstruct Pulmon Dis. 2012;7:555–69.
Millares L, Ferrari R, Gallego M, Garcia-Nuñez M, Pérez-Brocal V, Espasa M, Pomares X, Monton C, Moya A, Monsó E. 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.
Murphy TF, Brauer AL, Grant BJB, Sethi S. Moraxella catarrhalis in chronic obstructive pulmonary disease: burden of disease and immune response. Am J Respir Crit Care Med. 2005;172:195–9.
Lin CK, Kazmierczak BI. Inflammation: a double-edged sword in the response to Pseudomonas aeruginosa infection. J Innate Immun. 2017;9:250–61.
Parameswaran GI, Wrona CT, Murphy TF, Sethi S. Moraxella catarrhalisacquisition, airway inflammation and protease-antiprotease balance in chronic obstructive pulmonary disease. BMC Infect Dis. 2009;9:178.
Jacobs DM, Ochs-Balcom HM, Zhao J, Murphy TF, Sethi S. Lower airway bacterial colonization patterns and species-specific interactions in chronic obstructive pulmonary disease. J Clin Microbiol. 2018;56(10).
Toraldo DM, Conte L. Influence of the lung microbiota Dysbiosis in chronic obstructive pulmonary disease exacerbations: the controversial use of corticosteroid and antibiotic treatments and the role of Eosinophils as a disease marker. J Clin Med Res. 2019;11:667–75.
Smith CB, Kanner RE, Golden CA, Renzetti AD. Hæmophilus influenzæ and hæmophilus parainfluenzæ in chronic obstructive pulmonary disease. Lancet. 1976;307:1253–5.
Pericone CD, Overweg K, Hermans PWM, Weiser JN. Inhibitory and bactericidal effects of hydrogen peroxide production by Streptococcus pneumoniae on other inhabitants of the upper respiratory tract. Infect Immun. 2000;68:3990–7.
Wilkinson TMA, Aris E, Bourne S, Clarke SC, Peeters M, Pascal TG, Schoonbroodt S, Tuck AC, Kim V, Ostridge K, Staples KJ, Williams N, Williams A, Wootton S, Devaster J-M, AERIS Study Group. A prospective, observational cohort study of the seasonal dynamics of airway pathogens in the aetiology of exacerbations in COPD. Thorax. 2017;72:919–27.
Yilmaz P, Parfrey LW, Yarza P, Gerken J, Pruesse E, Quast C, Schweer T, Peplies J, Ludwig W, Glöckner FO. The SILVA and “all-species living tree project (LTP)” taxonomic frameworks. Nucleic Acids Res. 2014;42:D643–8.
We would like to thank the study participants and the site investigators listed as follows: Andreea Antonescu-Turcu/USA, Wesley Bray/USA, Ana Dancheva/Bulgaria, Stefan Denchev/Bulgaria, Anthony Dimarco/USA, Mark Dransfield/USA, Ali El Solh/USA, Ann Falsey/USA, Jonathan Ilowite/USA, Ivan Kiselov/Bulgaria, Petr Kolman/Czech Republic, Daniela Kopecka/Czech Republic, Radovan Kozel/Czech Republic, Peter Krumpe/USA, Camil Kreit/USA, DavidLaman/USA, Arturo Meade/USA, Frank Sciurba/USA, Milkana Simeonova/Bulgaria, Amir Sharafkhaneh/USA, Sotir Sotirov/Bulgaria, John Southard/USA, Anelia Stoyanova/Bulgaria, Patrick Sturm/USA, Keipp Talbot/USA, Josef Veverka/Czech Republic, and Jiri Vytiska/Czech Republic. We would also like to thank Alison Humbles, Paul Warrener, Taylor Cohen, Omari Jones, Fiona Fernandes and Outi Vaarala for helpful discussions.
This study was fully funded by AstraZeneca.
Ethics approval and consent to participate
This observational study is registered under NCT01455402 at ClinicalTrials.gov. The protocol was approved by independent institutional review boards, and all subjects signed written informed consent at enrollment.
All authors are employees and share holders of AstraZeneca.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. Sample characteristics. Table S2. Statistical significance of Alpha and Beta microbiota diversity metrics against demographic and clinical variables using Qiime2’s diversity plugin. Figure S1. Shannon diversity following antibiotic treatment in Europe and USA patients. Figure S2. Microbiota composition and diversity across study sites. Figure S3. Principal coordinate analysis of weighted UniFrac distances of the sputum microbiota colored by (A) geography and sample type, (B) most predominant bacterial taxa, and (C) viral infection. Figure S4. Bacteria associated with COPD exacerbation or frequency of exacerbation. (A) Abundance of Moraxella and H.parainfluenzae identified by analysis of composition (ANCOM) of microbiota between stable and exacerbated samples. (B) Cladogram of bacterial taxa identified by ANCOM comparing Frequent (> 2 exacerbation event/ year) to Infrequent exacerbator. Figure S5. Adjusted odds ratio of bacterial abundance (top/bottom quartile) in stable samples only to be associated with frequent exacerbations. Figure S6. Cladogram of 55 Differentially abundant bacterial taxa identified by ANCOM comparing consistent to variable longitudinal microbiota patient profiles at stable state.
About this article
Cite this article
Bouquet, J., Tabor, D.E., Silver, J.S. et al. Microbial burden and viral exacerbations in a longitudinal multicenter COPD cohort. Respir Res 21, 77 (2020). https://doi.org/10.1186/s12931-020-01340-0