Europe PMC

This website requires cookies, and the limited processing of your personal data in order to function. By using the site you are agreeing to this as outlined in our privacy notice and cookie policy.

A longitudinal characterization of the Non-Cystic Fibrosis Bronchiectasis airway microbiome.

Author information

Affiliations

  • 1. Department of Biological Sciences, University of Calgary, Calgary, AB, Canada.
      Authors
    • Woo TE 1, 2
    • Storey DG 1, 3
    (2 authors)
  • 2. Department of Medicine, University of Calgary, Calgary, AB, Canada.
      Authors
    • Woo TE 1, 2
    • Lim R 2
    • Rabin HR 2, 3
    • Mody CH 2, 3
    • Somayaji R 2, 3
    • Parkins MD 2, 3
    (6 authors)
  • 3. Department of Microbiology, Immunology and Infectious Diseases, University of Calgary, Calgary, AB, Canada.
      Authors
    • Heirali AA 3
    • Acosta N 3
    • Rabin HR 2, 3
    • Mody CH 2, 3
    • Somayaji R 2, 3
    • Surette MG 3
    • Sibley CD 3
    • Storey DG 1, 3
    • Parkins MD 2, 3
    (9 authors)

ORCIDs linked to this article

Scientific Reports, 03 May 2019, 9(1):6871
https://doi.org/10.1038/s41598-019-42862-y PMID: 31053725 PMCID: PMC6499777

This article is in the Europe PMC Open access subset. Refer to the copyright information in the article for licensing details.
Free full text in Europe PMC

Abstract 

A diverse microbiota exists within the airways of individuals with non-cystic fibrosis bronchiectasis (nCFB). How the lung microbiome evolves over time, and whether changes within the microbiome correlate with future disease progression is not yet known. We assessed the microbial community structure of 133 serial sputa and subsequent disease course of 29 nCFB patients collected over a span of 4-16 years using 16S rRNA paired-end sequencing. Interestingly, no significant shifts in the microbial community of individuals were observed during extended follow-up suggesting the microbiome remains relatively stable over prolonged periods. Samples that were Pseudomonas aeruginosa culture positive displayed markedly different microbial community structures compared to those that were positive for Haemophilus influenzae. Importantly, patients with sputum of lower microbial community diversity were more likely to experience subsequent lung function decline as defined by annual change in ≥-1 FEV1% predicted. Shannon diversity values <1 were more prevalent in patients with FEV1 decline (P = 0.002). However, the relative abundance of particular core microbiota constituents did not associate with risk of decline. Here we present data confirming that the microbiome of nCFB individuals is generally stable, and that microbiome-based measurements may have a prognostic role as biomarkers for nCFB.

Free full text 

Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
Sci Rep. 2019; 9: 6871.
Published online 2019 May 3. https://doi.org/10.1038/s41598-019-42862-y
PMCID: PMC6499777
PMID: 31053725

A longitudinal characterization of the Non-Cystic Fibrosis Bronchiectasis airway microbiome

T. E. Woo,1,2 R. Lim,2 A. A. Heirali,3 N. Acosta,3 H. R. Rabin,2,3 C. H. Mody,2,3 R. Somayaji,2,3 M. G. Surette,3,4 C. D. Sibley,3 D. G. Storey,1,3 and M. D. Parkinscorresponding author2,3
Author information Article notes Copyright and License information Disclaimer

Associated Data

Supplementary Materials
Data Availability Statement

Abstract

A diverse microbiota exists within the airways of individuals with non-cystic fibrosis bronchiectasis (nCFB). How the lung microbiome evolves over time, and whether changes within the microbiome correlate with future disease progression is not yet known. We assessed the microbial community structure of 133 serial sputa and subsequent disease course of 29 nCFB patients collected over a span of 4–16 years using 16S rRNA paired-end sequencing. Interestingly, no significant shifts in the microbial community of individuals were observed during extended follow-up suggesting the microbiome remains relatively stable over prolonged periods. Samples that were Pseudomonas aeruginosa culture positive displayed markedly different microbial community structures compared to those that were positive for Haemophilus influenzae. Importantly, patients with sputum of lower microbial community diversity were more likely to experience subsequent lung function decline as defined by annual change in ≥−1 FEV1% predicted. Shannon diversity values <1 were more prevalent in patients with FEV1 decline (P = 0.002). However, the relative abundance of particular core microbiota constituents did not associate with risk of decline. Here we present data confirming that the microbiome of nCFB individuals is generally stable, and that microbiome-based measurements may have a prognostic role as biomarkers for nCFB.

Subject terms: Infectious-disease diagnostics, Microbiome

Introduction

Bronchiectasis manifests as the abnormal dilation and thickening of the airways secondary to impaired mucocilliary clearance1. Studies of bronchiectasis have focused on infections caused by Pseudomonas aeruginosa and Haemophilus influenzae due to the prevalence of these pathogens in cultured respiratory samples, relationship with accelerated lung function decline, increased exacerbation frequency, and associated morbidity and mortality2,3. More recently, however, our understanding of chronic respiratory infections has evolved from a one-host-one-pathogen model to a polymicrobial model of infection4. A community of microorganisms, including bacteria, fungi, and viruses, collectively refered to as the lung microbiome, have been demonstrated in airways5–7. However, the role these organisms play in health and disease is not well understood8,9.

Our knowledge of infections in non-cystic fibrosis bronchiectasis (nCFB) lags behind that in cystic fibrosis (CF). Whereas the lung microbiome of paediatric CF patients is dynamic and changes throughout their adolescence10,11, by adulthood a core climax community evolves which remains relatively stable thereafter11–13. The observation that microbial diversity within the airways of CF patients inversely correlates with lung function in cross-sectional studies suggests that these communities may influence disease progression6,9,13. However, how the microbiota changes over time and associates with disease progression remains to be elucidated. With this study, we sought to characterize the natural history of airways communities within individuals with nCFB using a longitudinal collection of samples. Furthermore, we sought to determine if the nCFB microbiota correlated with subsequent long-term pulmonary outcomes and might thus serve as a novel biomarker to predict future clinical course.

Results

Patient Demographics

One hundred and thirty-three sputum samples from 29 patients were included. Demographics and baseline patient characteristics are listed in Table 1. The median duration of follow-up was 5.8 years (IQR 5.0–7.9, range 4.3–16). The majority of patients were never-smokers (93%; 27/29). Etiologies of bronchiectasis included post-infectious (62%), idiopathic (34%) and ciliary dyskinesia (4%). Lung function was generally stable during follow up with a median change in FEV1% percent predicted of −0.29%/year (IQR: −1.3%–+0.65%). Twelve exacerbation events were captured in the selected sputa, from 10 of the 29 (34%) patients. Two patients ultimately died and no patients received lung transplantation during the study period.

Table 1

Demographics and treatments used in nCFB patients at enrolment.

Patient Demographics (n = 29)Decliners (n = 11)Stables (n = 18)Total (n = 29)
Median Age at Study Enrolment (IQR)63 (39–70)57 (50–68)60 (48–69)
Median FEV1 at enrolment (% predicted) (IQR)71 (45–93)53 (39–67)55% (40.5–74)
Respiratory Comorbidities
Sinusitis5 (45%)10 (56%)15 (52%)
Asthma1 (9%)3 (17%)4 (14%)
Chronic Obstructive Pulmonary Disease (COPD)02 (11%)2 (7%)
Treatments Recorded During Sample Collection
Inhaled corticosteroids7 (64%)14 (78%)21 (72%)
Long-acting beta-agonists (LABA)6 (55%)14 (78%)20 (69%)
Long-acting muscarinic antagonist (LAMA)1 (9%)8 (44%)9 (31%)
Short-acting beta-agonists (SABA)6 (55%)12 (67%)18 (62%)
Ipratropium bromide3 (27%)7 (39%)10 (34%)
Saline2 (18%)4 (22%)6 (21%)
Proton-pump inhibitors1 (9%)2 (11%)3 (10%)
H2-receptor antagonists1 (9%)1 (6%)2 (7%)
Cultured Pathogens
Pseudomonas aeruginosa 6 (55%)10 (56%)16 (55%)
Haemophilus influenzae 4 (36%)4 (22%)8 (28%)
Staphylococcus aureus 3 (27%)5 (28%)8 (28%)
Moraxella catarrhalis 2 (18%)1 (6%)3 (10%)
Streptococcus pneumoniae 1 (9%)3 (17%)4 (14%)

Inter-quartile Range (IQR) and percent are indicated by the parenthesis.

The nCFB Microbial Community Displays Inter- and Intra-patient Variability

A total of 7,698,972 sequences were generated (median of 49,437 sequences/sample, IQR: 33,655–73,464, range: 5136–170,812) with 3202 total observed OTUs (median of 138 OTU per sample, IQR: 90–215, range: 29–817). Across the sample population, the majority of OTUs were classified as Streptococcus (30%), Pseudomonas (28%), and Haemophilus (17%). The longitudinal lung microbiota of taxa found in >1% relative abundance in the dataset was visualized using taxonomic summaries (Fig. 1). Notably, we observed varying degrees of inter- and intra-patient heterogeneity. Whereas patients 22, 23 and 27 had remarkably stable lung microbiomes, patients 4, 11, 15 and 26 demonstrated considerable temporal variation. When samples collected during exacerbations were excluded, significant clustering of samples by patient was observed (PERMANOVA, p = 0.03) (Fig. 2). In order to understand the natural history of the nCFB airways microbiota, we compared the composition of samples by first binning them into distinct time intervals (Fig. 3A–C). No significant changes in alpha-diversity was observed in any of these periods relative to the initial samples and PERMANOVA analysis did not reveal any community-wide differences (p = 0.54). In particular, Bray-Curtis analysis did not reveal any significant inter- and intra-patient changes between the initial, intermediate or final sample intervals - although certain individuals had dynamic ranges (Fig. 3D).

An external file that holds a picture, illustration, etc. Object name is 41598_2019_42862_Fig1_HTML.jpg

Taxonomic summaries of OTU present in >1% relative abundance of all samples (n = 133) obtained from a cohort of 29 nCFB patients as a function of the years from first sample. Samples for each patient are grouped together by a unique patient identifier – as indicated by the grey heading. Patient’s numbers from the Decliner group are indicated in magenta. Analysis of the taxonomic summary reveals a unique microbiome associated with each patient and varying degrees of inter and intra-patient diversity. Samples collected under antibiotic pressure are indicated by red text indidcating year of collection.

An external file that holds a picture, illustration, etc. Object name is 41598_2019_42862_Fig2_HTML.jpg

Clustering of samples by patient during stable periods (excluding exacerbation events) was observed (n = 121). Samples obtained from each patient are designated by the patient number adjacent to each point. Samples were categorized by color based on the time since baseline sample. Patients classified as Decliners are indicated by a circle, whereas patients with a Stable course are represented by triangles. Significant clustering was observed using PERMANOVA analysis (p = 0.03).

An external file that holds a picture, illustration, etc. Object name is 41598_2019_42862_Fig3_HTML.jpg

The microbial community was observed to remain relatively constant over the study period as measured by Observed OTUs (A), Shannon-Wiener (B) and Simpson diversity (C), and Bray-Curtis dissimilarity (D). Samples were categorized based on the time interval relative to first sample as: 0 (n = 29), 0–2 years (n = 36), 2–4 (n = 22), and 4 + (n = 46) years since initial sample collection. The median and interquartile ranges (IQR) are represented by the middle, top, and bottom lines of each boxplot. No statistically significant comparisons were observed.

The nCFB microbiome associates with disease trajectory

We investigated whether the baseline microbiota - as established from the first sputum sample on record - correlated with either baseline disease severity or subsequent lung function decline using our prospective collection. Notably, baseline demographics were not significantly associated with changes in microbial diversity (data not shown). In our small sample population, we did not confirm that disease severity at baseline associated with either alpha or beta-diversity (Fig. S1). We did, however, observe that baseline microbial diversity correlated with subsequent rate of lung function decline (Fig. 4). Eleven (38%) patients met our a priori definition for Decliner (>−1%/year FEV1). Median rates of annual FEV1 decline between Decliners were −1.8% (IQR: −1.3–−1) vs. +0.4% (IQR: −0.07–1.25) in Stable patients (p = 0.004), respectively. Baseline measurements of Observed OTU (p = 0.90), Shannon diversity (p = 0.51), or Simpson diversity (p = 0.64) were not significantly different between both groups. When comparing these two groups, we identified that Decliners generally had lower alpha-diversity using the Observed OTU (p = 0.05), but not Shannon (p = 0.097) and Simpson (p = 0.06) indices although similar trends were evident. When rarefied, significant differences in Shannon (p = 0.035) and Simpson diversity (p = 0.033), but not Observed OTUs (p = 0.07) were observed. However, no community-wide differences were observed between the two patient groups (p = 1.0). When using alternate definitions for Decliners of −0.5%/year and −1.5%/year the same trends were evident, but did not reach significance (data not shown).

An external file that holds a picture, illustration, etc. Object name is 41598_2019_42862_Fig4_HTML.jpg

Sputum samples collected from patients (n = 44 samples, 11 patients) who were categorized as Decliners (≥1 FEV1% per year) generally displayed lower levels of alpha diversity as compared to Stables (n = 89, 18 patients). No significant differences were detected using the Observed OTU (A), Shannon (B), and Simpson (C) diversity indices. The median and interquartile ranges of each group are represented by the middle, top, and bottom lines of each boxplot. Statistical analysis was performed using the Wilcoxon-signed rank test.

As lower microbial diversity was observed to associate with Decliner status, we investigated whether the initial Observed OTU, Simpson, and Shannon diversity measurement was correlated with change in lung function. No significant correlation between the initial diversity measures for each patient and annual rate of lung function decline was found using a generalized linear model using the Shannon (p = 0.61), Observed (p = 0.23) and Simpson diversity (p = 0.61) indices (Fig. S2). Likewise, we found that no association between the abundance of specific genera (found within the overall community at >1% of reads) and subsequent lung function decline. The initial abundance (t = 0) of Pseudomonas (rs = −0.21), Haemophilus (rs = 0.34), Streptococcus (rs = −0.18), Prevotella (rs = −0.29), Fusobacterium (rs = −0.18), Gemella (rs = −0.10), Granicutella (rs = −0.06), Neisseria (rs = 0.08), Rothia (rs = −0.08), and Staphyloccocus (rs = 0.17) within each sample did not correlate strongly with lung function for each patient. Lastly, we examined whether Shannon diversity measurements associated with the clinical prognosis of patients within our study. Using a Shannon diversity index cut-off value of 1, used for biomarker studies in CF14, we observed that the majority of samples from Stable patients had a Shannon diversity value of >1 (73/93 samples) as compared to Decliners (20/40 samples) (p = 0.002). When the analysis was constrained by patient, similar trends were observed but did not reach significance (p = 0.39). Samples with a Shannon diversity index <1 were associated with changes in the microbial community (p = 0.04). In particular the abundance of seven genera had relative lower abundance; Gemella, Lactobacillus, Nicoletella, Haemophilus, Bordetella, Staphylococcus, and Achromobacter (Fig. S3) whereas Streptococcus and Moraxella were increased.

Pseudomonas and Haemophilus display an antagonistic relationship in nCFB airways

We observed a strong correlation between airway pathogens obtained during routine real-time clinical laboratory culture and microbial community diversity determined retrospectively from these same samples. P. aeruginosa was isolated in 60 (45%) of the total samples (16/29 patients), whereas H. influenzae was only found in 17 samples (13%) (8/29 patients). Samples which were culture positive for P. aeruginosa trended towards coming from individuals with worse baseline lung disease relative to H. influenzae positive individuals (median FEV1 of 53% (IQR: 36–62) versus 65% (IQR: 39–75) (p = 0.10). We observed that the real-time isolation of P. aeruginosa was exclusive to H. influenzae (and vice-versa) such that these two organisms were never recovered together in culture from the same sample. Further differences arose when assessing the community microbiota in samples obtained in association with P. aeruginosa, H. influenzae, and in samples with neither (Fig. 5). When comparing these three groups, we found that P. aeruginosa culture positive samples had relatively higher Simpson (p = 0.01) and Shannon (p = 0.05) diversity, but not Observed OTU diversity (p = 0.85) as compared to H. influenzae culture positive samples. No significant differences in the alpha diversity were observed between either P. aeruginosa or H. influenzae culture positive samples and samples where neither of these organisms were identified. Community-wide differences were observed using PERMANOVA analysis (p = 0.001). In particular, the relative abundance of 10 genera was found to be significantly different between samples (Fig. S4). Interestingly, the increased log-fold relative abundance of Pseudomonas was observed in tandem with a >5 fold decrease change in relative abundance of Haemophilus (Fig. S5).

An external file that holds a picture, illustration, etc. Object name is 41598_2019_42862_Fig5_HTML.jpg

Observed OTU (A), Shannon (B), and Simpson (C) diversity, and Bray-Curtis dissimilarity (D) illustrate clustering of P. aeruginosa and H. influenzae positive samples. Differences in the sputum microbial diversity of samples collected in relation to the growth of Pseudomonas aeruginosa (PA, n = 60) and Haemophilus influenzae (HI, n = 17) and in samples with neither (Neither, n = 56). The ellipses for each group represent the 90% confidence interval. The median and interquartile ranges (IQR) are represented by the middle, top, and bottom lines of each boxplot. Statistically significant comparisons were shown using the p-value.

Acute antibacterial therapies introduce transient shifts within the microbiome

When we examined the influence of acute antimicrobial therapy with microbial changes in community structure, only fluoroquinolone use was associated with changes in diversity. Observed OTU richness (p = 0.004) was reduced under acute antibacterial pressures, but not Shannon (p = 0.19) or Simpson (p = 0.55) diversity indices (Fig. 6). Interestingly, community-wide differences were observed using PERMANOVA analysis (p = 0.01). In particular, samples obtained during fluoroquinolone use were found to have significant changes in the abundance of 12 genera. Two of the twelve genera, Escherichia and Moraxella, had increased relative abundance and the remaining had decreased relative abundance including: Veillonella, Dialister, Azorhizophilus, Bulleidia, Capnocytophaga, Lactobacillus, Prevotella, Treponema, Haemophilus, and Achromobacter (Fig. S6).

An external file that holds a picture, illustration, etc. Object name is 41598_2019_42862_Fig6_HTML.jpg

Acute fluoroquinolones use was associated with reduced community richness as represented by the Observed OTU. Samples collected during acute fluoroquinolone use (FQ, n = 55) were compared to samples collected in their absence (No-FQ, n = 78). Significant differences were seen in Observed OTUs (A), but not for Shannon (B) or Simpson diversity (C). Significant community-wide differences were found using PERMANOVA when permutations were constrained by patient and visualized using PCoA (p = 0.01).

Discussion

The lung microbiome has been suggested to play an important role in modifying disease progression for individuals with CF and nCFB6,11,12. Whereas recent studies have alluded to the relative stability of the nCFB microbiome over short periods15, we believe our ability to assess serial samples collected up to 16 years apart, provides unique insight into the longer-term adaptations of the lung microbiome. Interestingly, we observed no significant changes in the alpha and beta-diversity of the lung microbiome. This is all the more surprising given the extended duration of observation (for some patients up to 16 years). Stability was not influenced by baseline clinical characteristics or disease severity. Indeed, even over protracted periods of time, we observed, as others have over short time periods, microbiota tended to cluster by patient, alluding to unique patterns of colonization15. In comparison, studies in CF suggest that the childhood microbiota fluctuates and changes over time, while the adult lung microbiome remains relatively stable10,12,16.

Given the limited capacity by which cultured pathogens within nCFB airways predict clinical outcomes, the potential of the microbiome to serve as a biomarker for clinical management is of great interest. For instance, early studies have alluded to the predominance of P. aeruginosa and Veillona species as predictors for future exacerbations17. In contrast to early point prevalence studies, we observed no significant relationship between alpha diversity and static measures of lung function6,18. However, these trends are not ubiquitous amongst all studies, with other studies showing no association between diversity and lung function19. Rather, we observed a general relationship between alpha diversity and the rate at which lung function declined in these patients. It is likely that the slow rate of progressive lung disease was prohibitive in establishing community-wide differences as no significant decline was observed between baseline and end-of-study lung function measures within our cohort. Furthermore, we observed that samples with Shannon diversity values <1 tended to be associated with Decliner status and may warrant further interest in its use as a biomarker.

Amongst the pathogens isolated from the airways, P. aeruginosa and H. influenzae play an important role. The presence of these organisms within the airways has been associated with decreased microbial diversity, worsened clinical outcomes, and increased morbidity and mortality4,11,20,21. Importantly, the interactions between these airway pathogens is of great interest, with an antagonistic relationship observed in both cultured and culture-independent microbiome4,8. Similarly, we believe that these trends may reflect both direct and indirect competition between these two species. It is possible that it is the accessory microbiome associated with Pseudomonas-dominant patients that may exclude the colonization of H. influenzae and vice-versa8. Further, this competition may extend to other organisms and may explain the change in genera observed within our study. Though our V3/V4 sequencing of the 16S rRNA allowed only for genus level classification, our complimentary analysis of real-time culture data revealed that the only Pseudomonas and Haemophilus species identified from these samples were P. aeruginosa and H. influenzae.

The impact of acute exacerbations on microbial communities is not well understood – with no significant changes in the microbiome associated with the onset of these events19. Antibiotic use is the mainstay of management of exacerbations in patients with CF and nCFB. Indeed, antibiotic therapy may influence the respiratory microbiota and play an important role in clinical prognosis4,11,20,21. Within our study and in congruence with prior findings15, we suggest that the broad-spectrum effects of fluoroquinolone decrease the overall richness without affecting the evenness of the microbial community.

Importantly, this retrospective study of prospectively collected samples has a number of limitations. While our sample size of 29 patients is limited, it remains comparable in size to similar microbiome studies in nCFB6,19, and even longitudinal studies in CF16,21. Furthermore, we believe that the repeated longitudinal measurements analysis – assessing a time period vastly greater than prior works - provides unique insight into the natural history of the disease well beyond cross-sectional studies. Notably, as the primary focus of our study was to assess longitudinal changes in the microbial community, our analysis may have been insensitive to the influence of antimicrobial therapies or the incidence of pulmonary exacerbations on microbial diversity – as seen in previous studies15,16,22,23. Given our use of “opportunistically collected samples” during clinical encounters, it may be that our analysis is insensitive to detect changes in the microbial diversity prior to the onset of an exacerbation. As such, studies characterizing samples in their longitudinal relation to such acute events are important in highlighting other factors which may lead to onset of exacerbation21. Furthermore, as the vast majority of samples were collected during stable periods it is possible that little antibacterial selective pressures persisted to drive changes in the microbial community. Previous studies in CF have observed that microbial diversity did not change significantly during periods of stability16, which we suspect also extends to nCFB. Lastly, our collection of samples is derived from a regional referral clinic and patients are not exclusively seen by nCFB clinic providers and therefore may not be all encompassing.

The Calgary nCFB Biobank represents a unique opportunity in which to characterize the longitudinal nCFB lung microbiome. To our knowledge, this study represents the longest examination of the nCFB microbiome to-date suggesting that the lung microbiome is relatively stable over time, and is highly individualized. We suggest that the lung microbial diversity may be an important contributor to clinical course – although this is likely but one of many host, organism/community and environmental factors. Furthermore, we confirmed the reciprocal relationship between P. aeruginosa and H. influenzae in nCFB airways4,8. This study has provided further insight into the longitudinal microbiome of individuals with nCFB, and suggests that further studies exploring the microbiome’s association with subsequent clinical outcomes are warranted in our attempt to develop biomarkers for patient prognostication.

Methods

Patients and samples

The Calgary Bronchiectasis Biobank (1998–current) is a prospectively inventoried collection of sputum obtained from nCFB adults (>18 years) as part as their routine clinical care and stored at −80 °C. Patients followed at the Calgary Bronchiectasis Clinic have symptoms consistent with and radiographic evidence of bronchiectasis24. For inclusion in this study patients had to have ≥4 sputum samples spanning ≥4 years and a diagnosis of nCFB. Individuals with bronchiectasis due to CF were excluded.

Clinical Information and Definitions

Patient demographic and clinical information was collected through detailed chart review. Pathogens detected through traditional culture-based approaches from these samples were reported in real-time from the clinical microbiology laboratory and these data were extracted during chart review. Baseline patient characteristics including bronchiectasis aetiology, gender, respiratory therapies, antibacterial therapies and co-morbidities were analyzed for their relationship to microbial diversity. Bronchiectasis aetiology was classified as post-infectious, idiopathic, or other. Lung disease stages were classified by spirometry as determined using the Knudson calculation as advanced; Forced expiratory volume in one second (FEV1) < 40% predicted, moderate; 40–70%, mild; ≥70%18. Rates of lung function decline were calculated by dividing the FEV1% by the total number of years followed for each patient. Patients with an annual decline in FEV1% predicted of ≥1% were a priori considered Decliners (D)25,26 whereas those whose decline was less than this were considered Stables (S). Samples were separately coded and considered to be collected during acute antimicrobial therapy if obtained within two weeks of a new antibiotic administration (as acute parenteral antibiotic therapies have been associated with transient changes in the microbiota of suppurative lung disease)27–29. In order to assess differences in community composition over time, samples were categorized based on years from collection relative to the initial sample (0), >0.25–2 years, >2–4 years, and >4 years for each patient over the study period. The study was were performed in accordance with relevant guidelines and regulations of the University of Calgary and consistent with those required by Scientific Reports. The Conjoint Health Research Board of the University of Calgary has approved the ongoing collection and maintenance of the biobank and the samples maintained in it (REB16-0035). Furthermore, the CHREB has granted the investigators permission to evaluate patient outcomes associated with biobank derived samples REB16-0854.

DNA Isolation, 16s rRNA Variable 3–4 (V3/V4) Amplification, and Sequence Processing

Genomic DNA was isolated from sputum samples as previously described30,31. Barcoded universal primers adapted from Bartam et al. were used to amplify the V3 and V4 region of the 16s rRNA prior to MiSeq Illumina sequencing32. Reagent blanks were run for each set of DNA extractions and DNA amplification. If contamination was seen in the technical controls, samples were excluded and the DNA extraction and amplification was subsequently repeated. Technical replicates were not sequenced. Sequence analysis was carried out using previously published Perl scripts33. Low quality reads and primer sequences were removed using Cutadapt34 and paired-end reads were merged using PANDAseq35. Operational taxonomic units (OTUs) were then generated using AbundantOTU+ based on ≥97% similarity of sequences and classified taxonomically using the Ribosomal Database Project Classifier and the Greengenes database36–38. No samples were excluded based on our a priori requirement cut-off of >3500 reads. Bacterial OTU tables were generated by the removal of singletons (sequences with single reads) and non-bacterial OTU’s (Supplementary Table 1). Singletons were removed in an effort to more accurately evaluate community diversity by reducing the impact of base substitutions, low-quality reads, variable read lengths, non-target amplification, and undetected chimeric sequences39,40. The removal of singletons, while debatable, has been shown to reduce alpha-diversity but does significantly alter beta-diversity and the removal of singletons likely reflects a more accurate community39.

Microbial Community and Statistical Analysis

Analysis of the resulting OTU tables was performed in R (V 3.99) using the phyloseq package41. For alpha-diversity and beta-diversity calculations, we included OTUs that were present in ≥1% of total relative abundance. Alpha-diversity was calculated using Observed OTUs as measure of richness, and Shannon and Simpson’s diversity indices which include richness and evenness were all calcuated6,42,43. Unless otherwise stated, all analysis was performed without the use of rarefication as it is suggested to avoid introducing false positives44. However, a secondary analysis of alpha-diversity was performed after rarefying the samples to a read depth of 15,000. This analysis was not included within the manuscript unless a discrepancy with the non-rarfied data was observed. If discrepancies arose, both the non-rarefied and rarefied data are presented. Alpha diversity indices were compared using Wilcoxon Signed Rank for paired non-parametric factors, and Mann-Whitney U test for unpaired non-normally distributed variables. Beta-diversity differences in the microbiome were analyzed using Bray-Curtis distances and visualized using principal coordinate analysis. Community-wide OTU level differences were assessed using the permutational multivariate analysis of variance (PERMANOVA) using distance matrices using 999 permutations following proportional normalization of the input data13,45. If community-wide differences were identified, OTUs belonging to the same genus were analysed using the DeSeq2 package using the test = “wald” and fitType = “parametric” settings46. For multiple testing, we utilized the Benjamin Hochberg multiple test correction. Log abundance plots were visualized using ggplot247. A generalized linear model was constructed in R to examine the relationship between microbial diversity and repeated lung function measurements. Each analyzed variable was tested independently. The correlation between the abundance of individual genera and FEV1 was examined using the Spearman rank correlation coefficient (rs). Analysis of clinical characteristic and dichotomous variables was done using the Chi-Squared test and Fischer’s exact test using Prism 5.0 (GraphPad). Asymmetrically distributed variables were represented using median and interquartile ranges (IQR).

Supplementary information

Figs S1-S6(2.6M, docx)

Acknowledgements

We would like to thank our clinical partners at the Calgary Bronchiectasis Clinic and Calgary Laboratory Services for their assistance in procuring and maintaining samples in the biobank. This study was funded by a grant from the Lung Association of Alberta and Northwest Territories and a grant from the Natural Sciences and Engineering Research Council of Canada (RGPIN 203595-2013). T.E.W. was supported by a scholarship from Natural Sciences and Engineering Research Council of Canada and Alberta Innovates Technology Futures.

Author Contributions

T.E.W. was responsible for the collection strains. T.E.W. and A.H. performed the genomic extractions and DNA amplification. T.E.W., R.L., H.R.R. and C.H.M. were responsible for patient data collection. T.E.W., N.A., M.G.S., R.S., M.D.P. and C.S. contributed to the statistical analysis. All authors contributed to the preparation of the manuscript. The Calgary Non-CF Bronchiectasis Biobank is maintained by H.R.R. and M.D.P. D.G.S. and M.D.P. envisioned the project together and serve as guarantors of the work.

Data Availability

The data from which this study derives is available for analysis as per Scientific Reports requirements (PRJNA514329).

Notes

Competing Interests

T.E.W., R.L., A.A.H., N.A., C.H.M., R.S. and C.D. declare no competing interests. M.D.P., H.R.R. and M.G.S. have received research support from Gilead Sciences.

Footnotes

Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Supplementary information accompanies this paper at 10.1038/s41598-019-42862-y.

References

1. Weycker D, Edelsberg J, Oster G, Tino G. Prevalence and economic Burden of Bronchiectasis. Clin. Pulm. Med. 2005;12:205–209. 10.1097/01.cpm.0000171422.98696.ed. [CrossRef] [Google Scholar]
2. Davies G, Wells AU, Doffman S, Watanabe S, Wilson R. The effect of Pseudomonas aeruginosa on pulmonary function in patients with bronchiectasis. Eur. Respir. J. 2006;28:974–9. 10.1183/09031936.06.00074605. [Abstract] [CrossRef] [Google Scholar]
3. Ho PL, et al. The effect of Pseudomonas aeruginosa infection on clinical parameters in steady-state bronchiectasis. Chest. 1998;114:1594–1598. 10.1378/chest.114.6.1594. [Abstract] [CrossRef] [Google Scholar]
4. Purcell P, et al. Polymicrobial airway bacterial communities in adult bronchiectasis patients. BMC Microbiol. 2014;14:130. 10.1186/1471-2180-14-130. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
5. Sibley CD, et al. Culture enriched molecular profiling of the cystic fibrosis airway microbiome. PLoS One. 2011;6:25–27. [Europe PMC free article] [Abstract] [Google Scholar]
6. Rogers GB, et al. Clinical measures of disease in adult non-CF bronchiectasis correlate with airway microbiota composition. Thorax. 2013;68:731–7. 10.1136/thoraxjnl-2012-203105. [Abstract] [CrossRef] [Google Scholar]
7. Green H, Jones AM. The Microbiome and Emerging Pathogens in Cystic Fibrosis and Non–Cystic Fibrosis Bronchiectasis. Semin. Respir. Crit. Care Med. 2015;36:225–235. 10.1055/s-0035-1546752. [Abstract] [CrossRef] [Google Scholar]
8. Rogers GB, van der Gast CJ, Serisier DJ. Predominant pathogen competition and core microbiota divergence in chronic airway infection. ISME J. 2015;9:217–225. 10.1038/ismej.2014.124. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
9. van der Gast CJ, et al. Partitioning core and satellite taxa from within cystic fibrosis lung bacterial communities. ISME J. 2011;5:780–791. 10.1038/ismej.2010.175. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
10. Cox MJ, et al. Airway Microbiota and Pathogen Abundance in Age-Stratified Cystic Fibrosis Patients. PLoS One. 2010;5:e11044. 10.1371/journal.pone.0011044. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
11. Coburn B, et al. Lung microbiota across age and disease stage in cystic fibrosis. Sci. Rep. 2015;5:10241. 10.1038/srep10241. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
12. Carmody LA, et al. The daily dynamics of cystic fibrosis airway microbiota during clinical stability and at exacerbation. Microbiome. 2015;3:12. 10.1186/s40168-015-0074-9. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
13. Acosta, N. et al. The Evolving Cystic Fibrosis Microbiome: A Comparitive Cohort Study Spanning Sixteen Years. Ann. Am. Thorac. Soc. 14  1288–1297 (2017). [Abstract]
14. Acosta N, et al. Sputum microbiota is predictive of long-term clinical outcomes in young adults with cystic fibrosis. Thorax. 2018;73:1016–1025. 10.1136/thoraxjnl-2018-211510. [Abstract] [CrossRef] [Google Scholar]
15. Cox Michael J., Turek Elena M., Hennessy Catherine, Mirza Ghazala K., James Phillip L., Coleman Meg, Jones Andrew, Wilson Robert, Bilton Diana, Cookson William O. C., Moffatt Miriam F., Loebinger Michael R. Longitudinal assessment of sputum microbiome by sequencing of the 16S rRNA gene in non-cystic fibrosis bronchiectasis patients. PLOS ONE. 2017;12(2):e0170622. 10.1371/journal.pone.0170622. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
16. Zhao J, et al. Decade-long bacterial community dynamics in cystic fibrosis airways. Proc. Natl. Acad. Sci. USA. 2012;109:5809–14. 10.1073/pnas.1120577109. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
17. Rogers GB, et al. A Novel Microbiota Stratification System Predicts Future Exacerbations in Bronchiectasis. Ann. Am. Thorac. Soc. 2014;11:496–503. 10.1513/AnnalsATS.201310-335OC. [Abstract] [CrossRef] [Google Scholar]
18. Paganin P, et al. Changes in cystic fibrosis airway microbial community associated with a severe decline in lung function. PLoS One. 2015;10:1–19. 10.1371/journal.pone.0124348. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
19. Tunney MM, et al. Lung microbiota and bacterial abundance in patients with bronchiectasis when clinically stable and during exacerbation. Am. J. Respir. Crit. Care Med. 2013;187:1118–1126. 10.1164/rccm.201210-1937OC. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
20. Rogers GB, Bruce KD, Martin ML, Burr LD, Serisier DJ. The effect of long-term macrolide treatment on respiratory microbiota composition in non-cystic fibrosis bronchiectasis: an analysis from the randomised, double-blind, placebo-controlled BLESS trial. Lancet. Respir. Med. 2014;2:988–996. 10.1016/S2213-2600(14)70213-9. [Abstract] [CrossRef] [Google Scholar]
21. Whelan FJ, et al. Longitudinal sampling of the lung microbiota in individuals with cystic fibrosis. PLoS One. 2017;12:e0172811. 10.1371/journal.pone.0172811. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
22. Carmody LA, et al. Changes in cystic fibrosis airway microbiota at pulmonary exacerbation. Ann. Am. Thorac. Soc. 2013;10:179–187. 10.1513/AnnalsATS.201211-107OC. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
23. Collie D, et al. Lung Microbiota Changes Associated with Chronic Pseudomonas aeruginosa Lung Infection and the Impact of Intravenous Colistimethate Sodium. PLoS One. 2015;10:e0142097. 10.1371/journal.pone.0142097. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
24. Chalmers JD, et al. The EMBARC European Bronchiectasis Registry: protocol for an international observational study. ERJ Open Res. 2016;2:00081–2015. 10.1183/23120541.00081-2015. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
25. Fardon, T. et al. Identifying Modifiable Risk Factors for Rapid Lung Function Decline in Bronchiectasis. In American journal of respiratory and critical care medicine A2878 (2016).
26. Martínez-García MA. Factors Associated With Lung Function Decline in Adult Patients With Stable Non-Cystic Fibrosis Bronchiectasis. Chest. 2007;132:1565. 10.1378/chest.07-0490. [Abstract] [CrossRef] [Google Scholar]
27. Smith DJ, et al. Pyrosequencing reveals transient cystic fibrosis lung microbiome changes with intravenous antibiotics. Eur. Respir. J. 2014;44:922–930. 10.1183/09031936.00203013. [Abstract] [CrossRef] [Google Scholar]
28. Stressmann FA, et al. Long-term cultivation-independent microbial diversity analysis demonstrates that bacterial communities infecting the adult cystic fibrosis lung show stability and resilience. Thorax. 2012;67:867–873. 10.1136/thoraxjnl-2011-200932. [Abstract] [CrossRef] [Google Scholar]
29. Nielsen S, et al. Disrupted progression of the intestinal microbiota with age in children with cystic fibrosis. Sci. Rep. 2016;6:1–11. 10.1038/s41598-016-0001-8. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
30. Thevaranjan N, et al. Streptococcus pneumoniae colonization disrupts the microbial community within the upper respiratory tract of aging mice. Infect. Immun. 2016;84:IAI.01275-15. 10.1128/IAI.01275-15. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
31. Whelan FJ, et al. The loss of topography in the microbial communities of the upper respiratory tract in the elderly. Ann. Am. Thorac. Soc. 2014;11:513–521. 10.1513/AnnalsATS.201310-351OC. [Abstract] [CrossRef] [Google Scholar]
32. Bartram AK, Lynch MDJ, Stearns JC, Moreno-Hagelsieb G, Neufeld JD. Generation of Multimillion-Sequence 16S rRNA Gene Libraries from Complex Microbial Communities by Assembling Paired-End Illumina Reads. Appl. Environ. Microbiol. 2011;77:3846–3852. 10.1128/AEM.02772-10. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
33. Whelan FJ, Surette MG. A comprehensive evaluation of the sl1p pipeline for 16S rRNA gene sequencing analysis. Microbiome. 2017;5:100. 10.1186/s40168-017-0314-2. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
34. Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet.journal. 2011;17:10–12. 10.14806/ej.17.1.200. [CrossRef] [Google Scholar]
35. Masella AP, Bartram AK, Truszkowski JM, Brown DG, Neufeld JD. PANDAseq: paired-end assembler for illumina sequences. BMC Bioinformatics. 2012;13:31. 10.1186/1471-2105-13-31. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
36. DeSantis TZ, et al. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl. Environ. Microbiol. 2006;72:5069–5072. 10.1128/AEM.03006-05. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
37. Ye Y. Identification and Quantification of Abundant Species from Pyrosequences of 16S rRNA by Consensus Alignment. Proceedings. (IEEE. Int. Conf. Bioinformatics Biomed). 2011;2010:153–157. [Europe PMC free article] [Abstract] [Google Scholar]
38. Wang Q, Garrity GM, Tiedje JM, Cole JR. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 2007;73:5261–7. 10.1128/AEM.00062-07. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
39. Auer L, Mariadassou M, O’Donohue M, Klopp C, Hernandez-Raquet G. Analysis of large 16S rRNA Illumina data sets: Impact of singleton read filtering on microbial community description. Mol. Ecol. Resour. 2017;17:e122–e132. 10.1111/1755-0998.12700. [Abstract] [CrossRef] [Google Scholar]
40. Schirmer Melanie, Ijaz Umer Z., D'Amore Rosalinda, Hall Neil, Sloan William T., Quince Christopher. Insight into biases and sequencing errors for amplicon sequencing with the Illumina MiSeq platform. Nucleic Acids Research. 2015;43(6):e37–e37. 10.1093/nar/gku1341. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
41. McMurdie, P. J. & Holmes, S. Phyloseq: An R Package for Reproducible Interactive Analysis and Graphics of Microbiome Census Data. PLoS One8, e61217 (2013). [Europe PMC free article] [Abstract]
42. Rogers GB, et al. Reducing bias in bacterial community analysis of lower respiratory infections. ISME J. 2013;7:697–706. 10.1038/ismej.2012.145. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
43. Gihring TM, Green SJ, Schadt CW. Massively parallel rRNA gene sequencing exacerbates the potential for biased community diversity comparisons due to variable library sizes. Environ. Microbiol. 2012;14:285–290. 10.1111/j.1462-2920.2011.02550.x. [Abstract] [CrossRef] [Google Scholar]
44. McMurdie PJ, Holmes S. Waste Not, Want Not: Why Rarefying Microbiome Data Is Inadmissible. PLoS Comput. Biol. 2014;10:e1003531. 10.1371/journal.pcbi.1003531. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
45. Anderson M. A new method for non parametric multivariate analysis of variance. Austral Ecol. 2001;26:32–46. [Google Scholar]
46. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq. 2. Genome Biol. 2014;15:550. 10.1186/s13059-014-0550-8. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
47. Wickham H. ggplot2: Elegant Graphics for Data Analysis. New York: Springer-Verlag; 2009. [Google Scholar]
Articles from Scientific Reports are provided here courtesy of Nature Publishing Group

Full text links 

Read article at publisher's site: https://doi.org/10.1038/s41598-019-42862-y

Read article for free, from open access legal sources, via Unpaywall: https://www.nature.com/articles/s41598-019-42862-y.pdfUnpaywall

Citations & impact 

Impact metrics

28
Citations
Jump to Citations

Citations of article over time

Alternative metrics

Altmetric item for https://www.altmetric.com/details/59986204
Altmetric
Discover the attention surrounding your research
https://www.altmetric.com/details/59986204

Smart citations by scite.ai
Smart citations by scite.ai include citation statements extracted from the full text of the citing article. The number of the statements may be higher than the number of citations provided by EuropePMC if one paper cites another multiple times or lower if scite has not yet processed some of the citing articles.
Explore citation contexts and check if this article has been supported or disputed.
https://scite.ai/reports/10.1038/s41598-019-42862-y

Supporting
Mentioning
Contrasting
4
32
1

Article citations

Go to all (28) article citations

Data 

Data behind the article

This data has been text mined from the article, or deposited into data resources.

BioStudies: supplemental material and supporting data

Similar Articles 

To arrive at the top five similar articles we use a word-weighted algorithm to compare words from the Title and Abstract of each citation.

Funding 

Funders who supported this work.

CIHR, Cystic Fibrosis Canada, Lung Association of Alberta and NWT

    CIHR, Cystic Fibrosis Canada, Lung Association of Alberta and NWT, Crohn&apos;s and Colitis Foundation of Canada

      NA

        NSERC, Lung Association of Alberta and NWT