Impact of prophylactic and ‘rescue pack’ antibiotics on the airway microbiome in chronic lung disease
•,,,,,,,.
...
Abstract
The management of many chronic lung diseases involves multiple antibiotic prescriptions either to treat acute exacerbations or as prophylactic therapy to reduce the frequency of exacerbations and improve patients’ quality of life.
Aim To investigate the effects of antibiotics on the homeostasis of bacterial communities in the airways, and how this may contribute to antimicrobial resistance (AMR) among respiratory pathogens and microbiota.
Methods Within an observational cohort study, sputum was collected from 84 patients with chronic obstructive pulmonary disease and/or bronchiectasis at stable state: 47 were receiving antibiotic prophylaxis therapy. V3-V4 16S-rRNA sequencing on Illumina MiSeq, quantitative PCR for typical respiratory pathogens, bacteriology cultures and antimicrobial susceptibility testing of sputum isolates, resistome analysis on a subset of 17 sputum samples using MinION metagenomics sequencing were performed.
Finding The phylogenetic α-diversity and the total bacterial density in sputum were significantly lower in patients receiving prophylactic antibiotics (p=0.014 and 0.029, respectively). Antibiotic prophylaxis was associated with significantly lower relative abundance of respiratory pathogens such as Pseudomonas aeruginosa, Moraxella catarrhalis and members of family Enterobacteriaceae in the airway microbiome, but not Haemophilus influenzae and Streptococcus pneumoniae. No major definite directional shifts in the microbiota composition were identified with prophylactic antibiotic use at the cohort level. Surveillance of AMR and resistome analysis revealed a high frequency of resistance to macrolide and tetracycline in the cohort. AMR expressed by pathogenic bacterial isolates was associated with antibiotics prescribed as ‘rescue packs’ for prompt initiation of self-treatment of exacerbations (Spearman’s rho=0.408, p=0.02).
Conclusions Antibiotic prophylactic therapy suppresses recognised pathogenic bacteria in the sputum of patients with chronic lung disease. The use of antibiotic rescue packs may be driving AMR in this cohort rather than prophylactic antibiotics.
What is already known on this topic
In the management of chronic lung diseases, substantial amounts of antibiotics are prescribed either as prophylaxis or treatment of acute exacerbations.
Substantial scientific evidence supports the clinical usefulness of prophylactic antibiotic use, especially macrolides which can reduce the rate of exacerbations.
A key question is which of these strategies leads to a greater risk of altering the microbiome and/or the development of antimicrobial resistance.
What this study adds
We found that antibiotic prophylaxis suppressed common pathogenic bacteria such as Pseudomonas aeruginosa in sputum from patients with chronic lung diseases; however, use of antibiotic rescue packs to self-treat exacerbations may be associated with antimicrobial resistance in our cohort.
How this study might affect research, practice or policy
The potential impact of prescribed antibiotic rescue packs is concerning and requires review, especially in patients with frequent exacerbations who might benefit more from antibiotic prophylaxis.
Introduction
Chronic lung diseases such as bronchiectasis and chronic obstructive pulmonary disease (COPD) are associated with both a considerable socioeconomic burden and impacts on patients’ lives.1 Many people with chronic lung conditions also periodically experience intermittent acute deteriorations in respiratory health (exacerbations) which cause significant morbidity, impact significantly on quality of life and necessitate a change in regular medication.2 3 Frequent exacerbations have been associated with progressive lung damage, faster decline in lung function and worse quality of life.4 5
Guidelines for bronchiectasis and COPD recommend that patients have a self-management plan (SMP) in place for acute exacerbations. Within these SMPs, antibiotics are often prescribed and kept at home by patients as a ‘rescue pack’ to allow for prompt start of antibiotic treatment.2 3
Reducing the rate of exacerbations has been a major aim of patient management.6 Evidence from systematic reviews and meta-analyses of randomised placebo-controlled clinical trials has shown the clinical usefulness of prophylactic antibiotics in reducing the frequency of exacerbations in COPD7 8 and bronchiectasis.9 10 Azithromycin specifically has very good penetration into respiratory tissues, and a long serum half-life (up to 60 hours), permitting once-daily and even thrice-weekly dosing, with a characteristic postantibiotic effect11—making it the macrolide of choice.
Nevertheless, the adverse effects of long-term use of antimicrobials in chronic respiratory disease has not been fully investigated. Concern exists that the selection pressure on microbiota from a prolonged antibiotic regimen may alter host microbial homeostasis. In addition, there is the risk of AMR emergence in human microbiota, which may act as a reservoir for antimicrobial resistance dissemination into the wider population.12 Given that one aim of prophylactic antibiotics is to reduce the need and use of acute ‘rescue’ antibiotics in people with chronic lung diseases, a key question is which of these strategies leads to a greater risk of altering the microbiome and/or the development of antimicrobial resistance?
Here, we investigated the impact of prolonged antibiotic prophylactic therapy on the airway microbiome in patients with chronic lung disease, and the influence of prophylactic antibiotics and acute ‘rescue packs’ on the resistome and AMR in sputum bacterial isolates.
Methods
Study design and population
Participants with diagnosis of bronchiectasis (on CT) or COPD (on spirometry) were recruited between February 2017 and February 2018, participants were followed-up for 12 months. This was a convenience sample. The inclusion criteria were: age ≥20 years, confirmed diagnosis of COPD and/or bronchiectasis, history of frequent exacerbations (≥2/year), ability to spontaneously expectorate sputum and consent to participate. The exclusion criteria were diagnosis of cystic fibrosis or lung cancer, history of lung transplantation, known tuberculosis or HIV infections at the time of recruitment and clinical instability (sputum collected within 1 week before and 2 weeks after antibiotic treatment of an exacerbation or any other acute infection was excluded in the presented work). Clinical data collection included interviewing participants, and reviewing medical records and daily diary cards filled by participants reporting changes in the symptoms and treatment.
Patient and public involvement
Our research proposal was discussed with academics, physicians and students at University College London and with a sample of patients and their families/carers in pulmonary rehabilitation programme held in Peckwater Centre, London prior to submitting the ethics application. We adjusted the design in the light of their comments. Feedback on the study’s procedure was obtained from participants during clinic visits. Any concerning findings in sputum microbiology results were communicated with the treating physician.
Sputum processing and bacteriology cultures
Details on the methods are provided in online supplemental data. Sputum samples were processed with Sputasol (Oxoid, UK), comprehensive bacteriology cultures were obtained by plating sputum on Columbia agar with chocolate horse blood, Columbia colistin-nalidixic acid agar and cystine-lactose-electrolyte-deficient agar. All morphologically distinct colonies on the three cultures were isolated, purified and identified using matrix-assisted laser desorption/ionisation-time-of-flight (MALDI-TOF) mass spectrometry. Antimicrobial susceptibility of bacterial isolates was determined using the Kirby-Bauer agar diffusion technique according to the standards of the European Committee on Antimicrobial Susceptibility Testing (EUCAST).13
DNA extraction
DNA was then extracted from 1 mL sputum samples that had been centrifuged, pellet heated at 95°C for 30 min, then mechanically disrupted by bead-beating step, on the automated LIAISON Ixt extraction platform using DiaSorin Arrow-DNA extraction kit.
The total bacterial density in sputum indicated by the number of copies of the 16S rRNA gene per µL; in addition, the bacterial loads of Pseudomonas aeruginosa,14 and Haemophilus influenzae, Moraxella catarrhalis and Streptococcus pneumoniae15 were quantified using two multiplex quantitative PCR (qPCR) Taq-Man assays on a Qiagen Rotor-gene 6000 machine (online supplemental table S1).16 A mean of technical triplicates was taken for each sample. An internal amplification control, SPUD A, was used to test for PCR inhibition.
16S rRNA profiling
A sequence library was created by amplification of V3-V4 regions of the bacterial 16S rRNA gene using 341 forward-primer and 805 reverse-primer. The PCR products (584 bp amplicons) were cleaned using Agencourt AMPure XP beads (Beckman Coulter, UK). The samples were pooled in an equimolar ratio at 8 nM into one library which was checked by TapeStation (Agilent, USA). Sequencing was performed using Illumina MiSeq Platform using costume sequencing primers, and MiSeq Reagent Kit v2 (500 cycles).17
An extraction negative control and a no-template PCR control (water) were run alongside each batch of samples throughout the process as negative controls. A laboratory-prepared mock community was run as a positive control.
We adopted the workflow established by Microbiome helper18 using QIIME pipeline V.1.9.1.19 The appropriate statistical significance tests were calculated using QIIME wrapper scripts. STAMP (V.2.1.3)20 was used to visualise the results and testing the differential abundances of taxa between the two study groups using White’s non-parametric t-test.21 All p values were corrected using Benjamini-Hochberg false discovery rate method for multiple comparisons. The statistical analysis was performed using IBM SPSS V.25.0.22
Resistome analysis
Metagenomic sequencing was carried out on a subset of 17 sputum samples: 8 were from patients with bronchiectasis on prophylactic antibiotic therapy and 9 from the comparator group, using Oxford Nanopore MinION system. Human DNA was depleted as per the published method by Charalampous et al.23 Metagenomics libraries of six multiplexed samples were prepared using the Rapid Barcoding Kit (Ref: SQK-RLB004, ONT). Sequencing was performed on MinION flow cell (R.9.4.1) for 48 hours. The genomes were assembled using the miniasm/minimap pipeline.24 A Basic Local Alignment Search Tool (BLAST) search against the ResFinder database performed25 to detect AMR genes. The alignments with accuracy <90% have been excluded. The prevalence of AMR genes within the samples was measured in parts per million reads (ppm), that is, the number of sequence reads identified as AMR genes relative to the total number of reads representing the sample.
Results
The clinical and demographic characteristics of patients who were and were not on prescribed antibiotic prophylaxis therapy are shown in table 1. Forty-seven of 84 participants were prescribed antibiotic prophylaxis therapy for over a year, in the context of their routine medical care; in 66% these were macrolides, specifically, 29 patients were on azithromycin 250 mg thrice weekly and 8 on 500 mg thrice weekly; 2 were on clarithromycin (250 mg twice daily). Other prescribed antibiotics were co-trimoxazole (960 mg once daily or 1920 mg, divided into two doses, thrice weekly), ciprofloxacin (250 mg twice daily), doxycycline (100 mg once daily), amoxicillin (250 mg twice daily), cephalexin (500 mg twice daily) and phenoxy-methyl penicillin (250 mg twice daily).
Table 1
|
Demographic and clinical characteristics of participants
There were no significant differences in the proportions of patients with COPD and bronchiectasis in the two groups (p=0.157). All participants were regularly receiving the influenza vaccine and had pneumococcal vaccine. Lung function as indicated by the spirometry results was comparable between the two study groups (table 1).
A significantly higher proportion of patients with common variable immunodeficiency (CVID) were on antibiotic prophylaxis therapy (66%), compared with seven in the comparator group (19%) (p<0.0001); 90% of the patients with CVID were on immunoglobulin replacement therapy.
Seventy-nine per cent of the participants were prescribed a rescue pack of antibiotics to keep at home as an SMP; 47% of the rescue packs contained amoxicillin/clavulanic acid, 20% doxycycline, 14% quinolones (either ciprofloxacin, levofloxacin or moxifloxacin), 15% amoxicillin and 7.5% had macrolides either azithromycin or clarithromycin.
For comparing the microbiome of patients who were and were not receiving antibiotic prophylaxis therapy, one sputum sample per participant was selected to represent the microbiome profile at stable state. To avoid biases due to exacerbations and antibiotic treatment, sputum samples collected within 1 week before and 2 weeks after exacerbation and/or breakthrough antibiotic treatment were excluded.
The total bacteria density and the phylogenetic α-diversity were significantly lower in the sputum collected from patients on prophylactic antibiotics in comparison with those who were not (figure 1). The median 16S rRNA gene copies (IQR) was 6.20 (5.54–6.73) log10copies/µL in the antibiotic prophylaxis group vs 6.60 (6.17–6.92) log10copies/µL in the comparator group (p=0.029) (figure 1A). The difference of means of the phylogenetic α-diversity index PD whole tree in the two populations (SD) was 1.4 (0.56), p=0.014. This was despite the greater proportion of patients with bronchiectasis in the prophylactic antibiotic group (table 1) and the fact that the total bacterial density was significantly higher in the patients with bronchiectasis compared with those with COPD, where the median 16S rRNA gene copies (IQR) were 6.63 (6.18–6.96) vs 6.18 (5.41–6.88) log10copies/µL (p=0.042), respectively (figure 1B).
The total bacterial density, α-diversity and richness of airway bacterial community in patients with prophylactic antibiotics (AB) (n=47, navy blue) and those who were not (n=37, yellow). (A) Total number of 16S rRNA copies (a proxy to total bacterial density in sputum) was significantly lower in the prophylactic antibiotic group (p=0.029); however, (B) it was significantly higher in the patients with bronchiectasis (sky blue, n=50) compared with patients with chronic obstructive pulmonary disease (COPD) (red, n=23) (p=0.042). (C) α-Diversity expressed by phylogenetic diversity index: PD whole tree (p=0.014) and (D) richness expressed by total number of observed Operational Taxonomic Units (OTUs) (p=0.078).
Antibiotic prophylaxis was the only significant covariate in weighted-UniFrac β-diversity index (p=0.038 by permutational multivariate analysis of variance (PERMANOVA)) and unweighted-UniFrac β-diversity index (p=0.012 by PERMANOVA) (online supplemental data,online supplemental figure S1). Chronic lung disease (whether bronchiectasis or COPD) was a significant covariate in unweighted-UniFrac β-diversity index (p=0.01 by PERMANOVA) but not in weighted-UniFrac β-diversity index (p=0.152 by PERMANOVA). Primary immunodeficiency status (CVID) was not a significant covariate in both weighted and unweighted-UniFrac β-diversity indices (p>0.05 by PERMANOVA) (online supplemental data, online supplemental figure S2).
Comparison between the average airway microbiome profile of patients not receiving antibiotic (AB) prophylaxis (comparator group) (A) and those on prophylactic AB therapy (B) at the phylum level (p>0.05); (C) the relative abundance of phylum Synergistetes was significantly lower in patients receiving antibiotic prophylaxis therapy (p=0.00247); yellow: comparator group, navy blue: AB prophylactic group.
The bacterial community structure was shifted towards more Proteobacteria (p=0.041) in patients with bronchiectasis while in patients with COPD it was shifted towards more Firmicutes (p=0.048). Bacteroidetes was slightly lower in patients with COPD and was significantly associated with airflow obstruction (p=0.049) (online supplemental figures S3 and S4).
The microbiome profiles demonstrated that the microbiota composition in sputum was similar in patients who did or did not use prophylactic antibiotics at the phylum level. Only the phylum Synergistetes, which is represented by a minor taxon present at relative abundance (RA) of <0.5% in all samples, was significantly less abundant in the patients on prophylactic antibiotic therapy (p=0.002) (figure 2). At the genus level, the major taxa constituting the microbiome profiles and representing the core respiratory microbiota were similar in the patients who did or did not use prophylactic antibiotics (figure 3). Nevertheless, potentially pathogenic taxa were significantly less abundant in patients on prophylactic antibiotic therapy including: Pseudomonas (p=0.027), Enterobacteriaceae (p=0.021), Klebsiella (p=0.046), Pasteurella (p=0.012) and Morganella (p<0.0001). The genus Moraxella tended to be less abundant in patients on prophylactic antibiotic therapy; however, the observed difference did not reach statistical significance (p=0.087) (figure 4). Some minor taxa with RA <1% such as Bacteroidetes (p=0.02), Schwartzia (p=0.001) and Sphaerochaeta (p=0.01) were also significantly lower in the antibiotic prophylaxis group.
Airway microbiome profiles of patients with chronic lung disease who were receiving antibiotic (AB) prophylaxis therapy in comparison with patients in the comparator group; not receiving antibiotic prophylaxis; profiles of the Mock community resulted from two sequencing runs R3 and R4. f, family; g, genus; NC, no template negative control; NCext, extraction negative control.
Suppression of potential pathogenic taxa in the sputum of patients on prophylactic antibiotic therapy (navy blue) compared with others in the comparator group (yellow): (A) genus Pseudomonas (p=0.027), (B) family Enterobacteriaceae (p=0.021), (C) genus Moraxella (p=0.087).
qPCR confirmed that P. aeruginosa load was significantly supressed in the sputum of patients on antibiotic prophylaxis; the median load (IQR) was 2.89 (2.62–3.64) log10CFU/mL vs 7.23 (3.44–7.95) log10CFU/mL in those who were not (p=0.001), even though P. aeruginosa was more frequently detected in the antibiotic prophylaxis group; 61% vs 44% in the comparator group (figure 5D). M. catarrhalis was significantly less prevalent in the sputum of patients on antibiotic prophylaxis (4%) compared with 20% in those who were not (p=0.039). However, when detected, the load of M. catarrhalis was no different between the two groups (figure 5A,C). The prevalence and loads of both H. influenzae and S. pneumoniae were similar in both groups.
Prevalence, co-existence and bacterial loads of pathogenic airway bacteria in patients with chronic lung disease who were receiving antibiotic (AB) prophylaxis therapy (navy blue) and those not (comparator group, yellow) as determined by the multiplex quantitative PCRs. (A) Percentage of samples tested positive for the tested pathogens, (B) number of organisms tested positive within the same sputum samples (p=0.203), (C) mean bacterial loads of Haemophilus influenzae (Hi), Moraxella catarrhalis (Mc), Streptococcus pneumoniae (Spn) were similar in both groups (p=0.189, multivariate analysis of variance). (D) Pseudomonas aeruginosa (Pa) load was significantly lower in the prophylactic AB group (p=0.001). Bacterial load was measured in CFU/mL of original sputum sample for each bacterium. Samples which gave negative results for a given bacteria were excluded from the analysis. Error bars show ±1 SD. *P<0.05, **p≤0.01.
Co-existence of two respiratory pathogens was detected in 39% of the examined sputum samples (figure 5B). The most common combinations were P. aeruginosa+H. influenzae (12%) and H. influenzae+S. pneumoniae (9%). The co-existence of these three pathogens was detected in 7% of samples. Four pathogens were detected in one patient with bronchiectasis who was not receiving antibiotic prophylaxis. However, there were no significant differences in the co-existence of these respiratory pathogens in the sputum of patients who were or were not on antibiotic prophylaxis (p=0.203).
In bacterial cultures, a total of 187 distinct bacterial isolates were isolated from sputum cultures and identified by MALDI-TOF. The prevalence of the 68 sputum isolates; 36 from patients on prophylactic antibiotics and 32 from the comparator group are shown in figure 6. The variety and frequency of bacterial pathogens isolated from the sputum of patients using prophylactic antibiotics were significantly less than those not receiving prophylactic antibiotics in the comparator group. Twenty-six per cent of the Gram-negative isolates were isolated from sputum cultures of patients on prophylactic antibiotics compared with 48% from the rest of the cohort (p=0.011).
The taxa isolated from the sputum of patients with chronic lung disease and their prevalence among the patients who were receiving antibiotic (AB) prophylaxis therapy (n=36, navy blue) and those not (comparator group, yellow, n=32).
Antibiotic susceptibility testing revealed high proportions of resistance to macrolide (64% and 29%), tetracycline (45% and 34%) and ampicillin (10% and 67%) of the tested Gram-positive and Gram-negative isolates, respectively in the whole cohort. Ciprofloxacin and cefotaxime resistance was detected in 21% and 17% of the Gram-negative isolates, respectively (online supplemental figures S5 and S6). The observed AMR in the bacterial isolates was associated with the corresponding prescribed antibiotics in the rescue packs for self-treatment of acute exacerbations (Spearman’s correlation ρ=0.408, p=0.02) but not with the prophylactic antibiotics used (Spearman’s correlation ρ=−0.419, p=0.136) (figure 7B).
(A) Comparison between the antimicrobial resistance (AMR) in the commensal viridans streptococci group isolated from the sputum of patients receiving antibiotic (AB) prophylaxis therapy (group 1) and those not (comparator (group 0) (p<0.0001), (B) association between the AMR in pathogenic sputum isolates and the prescribed ABs in the rescue packs for the self-treatment of acute exacerbations (p=0.04) prevalence: percentage of the resistant isolates. AZM, azithromycin; AMP, ampicillin; CTX, cefotaxime; DO, doxycycline, FQ, fluoroquinolone; LEV, levofloxacin; R, resistant (red); S, sensitive (green); TE, tetracycline; VA, vancomycin.
Among 119 tested viridans streptococci isolates, resistance to multiple antibiotics was common in both groups such as: azithromycin (95% in the antibiotic prophylaxis group vs 76% in the comparator group), tetracycline (74% vs 62%, respectively), ampicillin (65% vs 76%, respectively), cefotaxime (61% vs 54%, respectively), levofloxacin (48% vs 32%, respectively) and vancomycin (29% vs 6%, respectively). Nevertheless, resistance was significantly more frequently detected in the antibiotic prophylaxis group compared with the comparator group (p<0.0001) (figure 7A).
In the resistome analysis, no significant differences in the prevalence of AMR genes were found in the sputum samples of patients with bronchiectasis who were on prophylactic antibiotics (seven were on 250 mg azithromycin thrice weekly and one on ciprofloxacin 250 mg twice daily) compared with those who were not (n=9) (p=0.48). Although not statistically significant due to the small sample sizes, a trend was observed in which the patients who were not receiving antibiotic prophylaxis therapy but had frequent exacerbations, median frequency of exacerbations/year 4 (IQR 2–5), exhibited a broader range of AMR gene prevalence: the median (IQR) was 21 ppm (6–121 ppm, n=9) compared with those on prophylactic antibiotics which was 28 ppm (IQR 14–43 ppm, n=8) (figure 8A). Also, a greater diversity of AMR genes conferring resistance to multiple antibiotic classes was observed in the frequent exacerbator patients using rescue packs and not on prophylactic antibiotics (figure 8B).
(A) Comparison between the relative abundance of the total antimicrobial resistance (AMR) genes (measured in ppm sequence reads) in sputum samples from patients with bronchiectasis (BTX) on antibiotic (AB) prophylaxis (blue, n=8) and other patients with BTX having frequent exacerbations treated with AB rescue packs but not receiving AB prophylaxis therapy (yellow, n=9) (p=0.48 by Mann–Whitney). (B) The prevalence of AMR genes in the patients on antibiotic prophylaxis and those not (*p<0.05 by MW). AG, aminoglycosides resistance; FQ, fluoroquinolones resistance; Rf, rifamycin resistance; Rp, ribosomal protection proteins; Sulph/trimeth, sulfonamides/trimethoprim.
Discussion
Antibiotic prophylaxis is a common approach in the management of chronic lung diseases including COPD and bronchiectasis, especially in advanced cases with frequent and/or severe exacerbations. Here, we report the impact of prolonged antibiotic prophylaxis therapy on the airway microbiome within a cohort of 84 patients with chronic respiratory conditions: COPD and/or bronchiectasis, 47 received prophylactic antibiotics prescribed for at least 1 year prior to joining the study as part of their routine clinical care. We also relate AMR in sputum-resistome and sputum-derived bacterial isolates to antibiotics: prophylactic and acute ‘rescue’ prescription.
Antibiotic prophylaxis was associated with lower α-diversity of the airway bacterial community, as indicated by the PD whole tree (a quantitative measure of the phylogenetic diversity within an ecosystem)26 and decline in the total bacterial burden in sputum. Previous studies have also reported a decline in richness and α-diversity in bronchoalveolar lavage (BAL) samples from patients with COPD and moderate-to-severe asthma who were receiving azithromycin prophylaxis therapy.27 28
The microbiome profiles at the phylum level were similar between those prescribed and not prescribed prophylactic antibiotics. Synergistetes (a minor phylum) was the only phylum significantly less abundant in the antibiotic prophylaxis group. Synergistetes is a recently recognised bacterial phylum which has been detected in various body microbiomes such as the oral cavity, gut, umbilicus and vagina. Since its presence has been associated with disease in sites such as periodontitis, abscesses and cysts, Synergistetes are best considered opportunistic pathogens.29
There were no significant differences in the relative abundances of any of the genera that would normally be defined as components of the healthy microbiota, between the patients receiving antibiotic prophylaxis and those who did not. Similarly, Rogers et al reported that there were no significant differences in microbiota composition between erythromycin treatment and placebo arms in the bronchiectasis and low-dose erythroymycin study (BLESS).30
Some pathogenic genera such as Pseudomonas, Enterobacteriaceae, Klebsiella, Pasteurella and Morganella were significantly less abundant in the patients receiving antibiotic prophylaxis compared with those who did not. Sequencing data suggested that Moraxella was also less abundant within the antibiotic prophylaxis group, although statistical significance was not achieved. Nevertheless, this observation was confirmed by specific qPCR (p=0.039).
qPCR is more sensitive and specific compared with 16S rRNA sequencing, data on the sensitivity and specificity of the methods used is presented in the online supplemental data, online supplemental figure S7, online supplemental tables S2 and S3. Apart from S. pneumoniae which could not be distinguished by V3-V4 16S rRNA sequencing from the viridans streptococci that are abundant in sputum, a strong significant correlation between the qPCR results of H. influenzae (Spearman’s ρ=0.798 p=8.4E-39) and M. catarrhalis (ρ=0.621, p=3.2E-15) and the corresponding taxa in the sequencing results of V3-V4 variable regions of 16S rRNA gene was found (online supplemental data, online supplemental figure S8). This correlation was less for P. aeruginosa (ρ=0.238 p=2.9E-4). Although P. aeruginosa was more frequently detected by qPCR in the antibiotic prophylaxis group, as a significantly higher proportion of CVID patients were in this study group and P. aeruginosa was significantly more prevalent in the patients with CVID (68%, p=0.01) compared with the rest of the cohort, the load of P. aeruginosa was in fact significantly lower in patients receiving prophylactic antibiotics (p=0.001). The same trend was observed after excluding the patients with CVID, although the sample size in this case was too small to assess true statistical significance (online supplemental figure S9).
The data from molecular methods were reflected in the bacterial culture results where pathogenic bacteria were isolated significantly less frequently from patients’ sputum using prophylactic antibiotic therapy compared with the other participants. This is in line with previous clinical trials which evaluated the efficacy of antibiotic prophylaxis therapy in various chronic lung diseases, and reported that respiratory colonisation with typical respiratory bacterial pathogens was either eliminated or inhibited following introduction of the prophylactic antibiotic.31–33 It is noteworthy that microbiology was a secondary outcome in these trials and the results in most cases were based on culture-dependent techniques only, which are less holistic compared with the microbiome approach and can be insensitive when detecting bacteria present at low loads.
This finding may explain how prophylactic antibiotics reduce the risk of exacerbations in chronic lung disease, where the mechanism would be long-term suppression of pathogenic bacteria within the airway microbiome rather than ‘snapshot’ elimination of pathogens at acute exacerbation events. This would support the use of antibiotic prophylaxis in the management of chronic lung disease.
Most participants (62%) in the antibiotic prophylaxis group in our cohort were on an intermittent azithromycin regimen (thrice weekly). Macrolides are reported to have additional anti-inflammatory and immunomodulatory properties34; therefore, the mechanism of action may not be mediated solely through antimicrobial activity. For example, macrolides have poor antibacterial activity against P. aeruginosa, but azithromycin was found to inhibit its biofilm formation by interfering with essential quorum sensing pathways.35 36 Segal et al have demonstrated that azithromycin could modulate the metabolome of the microbiome which in turn mediated the anti-inflammatory and immunomodulatory activity in the BAL samples of patients with COPD receiving azithromycin prophylaxis therapy.27
Surveillance of AMR revealed a high frequency of resistance across the whole cohort; however, it was significantly greater in the antibiotic prophylaxis group. For example, azithromycin resistance was detected in 95% of the viridans streptococci isolates in the antibiotic prophylaxis group vs 75% in the comparator group. Previous studies confirmed that the oropharyngeal carriage of macrolide resistant viridans streptococci is common in general populations; it was estimated as 71% in a Belgian cohort of healthy adults and 94% in another Spanish cohort. In the Belgian study, co-resistance to tetracycline was identified in 73% of the isolates.37 38
The resistome analysis revealed similar trends in which AMR genes were slightly richer (in terms of frequency of detection) in the samples from patients with bronchiectasis on prophylactic antibiotics, however, the broadest diversity of AMR genes was observed in those patients who had frequent exacerbations events treated with antibiotic rescue packs but were not receiving antibiotic prophylaxis. Nevertheless, due to the small sample sizes, the study lacked the power to demonstrate that any of the observed trends were statistically significant. The acquired and/or inherent resistance of the microbiota might explain the resilience of the airway microbial community in response to antimicrobial treatment.
The relationship between the development of AMR and use of antibiotic prophylaxis is unclear within the current literature. Many studies report a significant rise in the acquisition of AMR, or a significant elevation in the minimum inhibitory concentrations of clinical isolates, in response to antibiotic prophylaxis.31 39 40 The viridans streptococci develop resistance rapidly, especially to macrolides and tetracyclines.41 42 Whole genome sequencing of commensal streptococci in one study revealed that macrolide transmissible resistance genes were carried with tetracycline-resistance determinants on transposable elements.42 Commensal streptococci were also reported to carry the same macrolide resistance genes as pathogenic streptococci.43 Therefore, microbiota may be regarded as a reservoir for AMR in the microbial ecosystem.
On the other hand, some studies reported no significant differences in the detection rates of the antibiotic resistance between patients receiving antibiotic prophylaxis therapy and those who did not.32 44 Other studies have found that the resistance acquired during the antibiotic prophylaxis therapy was temporary, this suggested that acquired resistance is at the cost of fitness and does not persist.45 46
In current clinical practice, courses of antibiotics are often prescribed in ‘rescue packs’ kept by patients to initiate treatment of acute exacerbations promptly. The decision on antibiotic choice is based on clinician preference and may consider the individual’s clinical condition, drug tolerance and allergies, previous sputum bacteriology culture results and the antimicrobial susceptibility profile (antibiogram) of previous sputum pathogenic bacterial isolates. Nevertheless, the benefit of antibiotics in the treatment of acute exacerbation, especially in COPD, remains controversial, and the evidence supporting the universal treatment of exacerbation with antibiotics is limited.2 In our results, AMR detected in sputum isolates was more closely related to antibiotic rescue packs rather than to prophylactic antibiotics. Therefore, the practice of prescribing antibiotic rescue packs might be more concerning in this sector rather than antibiotic prophylaxis therapy, especially in patients who suffer frequent exacerbations. This may justify the benefit of antibiotic prophylaxis in these patients. Patients’ education on the rational use of antibiotics during exacerbations is crucial (such as having clear criteria when to start and stop them). Also, rotating prescriptions between different classes of antibiotics may theoretically help to mitigate AMR in this sector, although the available options may be limited and our data cannot inform on this.
There are several limitations to our results. As this was an observational study, participants were not randomly assigned to receive prophylactic antibiotics, therefore, our results demonstrate association but cannot prove causation. The decision to place a patient on antibiotic prophylaxis therapy was a clinical decision taken by the treating clinicians prior to joining the study. Antimicrobial prescription behaviour can be subjective and clinicians have different thresholds for prescribing antibiotic prophylaxis, although in most cases this decision is reserved for severe cases. This may explain why the frequency of exacerbations was comparable in both study groups despite the evidence that antibiotic prophylaxis therapy reduces the rate of exacerbations. Our cohort included patients with COPD, bronchiectasis and 45% of participants had an underlying diagnosis of CVID (the majority of whom were receiving immunoglobulin replacement). A relatively higher proportion of participants had bronchiectasis compared with COPD in the whole cohort; nevertheless, similar proportions of patients with bronchiectasis and COPD were present in both study groups. Apart from azithromycin, the small number of participants in the other six antibiotic prophylaxis regimes did not allow for comparisons between different antibiotics. The residual antimicrobial activity in the samples could have biased the bacteriology cultures in the antibiotic prophylaxis group but the molecular methods confirmed the findings of the culture-based approach. All these sources of variability might have masked some significant trends in the data.
In conclusion, antibiotic prophylaxis therapy was associated with reduced phylogenetic α-diversity of the bacterial communities and lower bacterial density in sputum. It selectively suppressed specific taxa that represent bacterial respiratory pathogens without disrupting the homeostasis of the respiratory microbiota. It did not induce a definitive compositional shift in the airway microbiota composition at the cohort level. In general, macrolide resistance was high in the whole cohort, nevertheless, it was significantly greater in patients receiving antibiotic prophylaxis. The practice of antibiotic rescue packs may be driving AMR in this cohort since the detected AMR expressed by sputum bacterial isolates was associated with prescribed antibiotics in the rescue packs kept by patients for self-treatment of exacerbations. Therefore, the clinical decision of antibiotic prophylaxis in the management of chronic lung disease should be carefully considered on an individual basis by weighing the benefits of suppressing pathogenic bacteria and reducing the rate of exacerbations, and hence need for acute antibiotic treatment courses, against the risk of enriching resistance to the prescribed antibiotic among bacterial populations and the considerable adverse effects that can result from the long-term use of these antimicrobial chemotherapeutic agents.
Contributors: SR, JRH, TMcH, ML, DS and JB designed the study. SR, SQ, ML,DML and JB recruited participants, conducted the study and collected clinical data and samples. SR processed the samples, did the laboratory work, performed the bioinformatics and data analyses. SR, JRH TMcH, ML, DML and DS interpreted the data. TMcH is the principal investigator and acts a guarantor of the study. SR wrote the first draft of the manuscript with contributions from all other authors.
Funding: The authors have not declared a specific grant for this research from any funding agency in the public, commercial or not-for-profit sectors.
Disclaimer: The funder had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Competing interests: Newton Mosharafa funds by British Council Egypt and Egyptian Ministry of Higher Education sponsored the PhD of SR.
Patient and public involvement: Patients and/or the public were involved in the design, or conduct, or reporting, or dissemination plans of this research. Refer to the 'Methods' section for further details.
Provenance and peer review: Not commissioned; externally peer reviewed.
Supplemental material: This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.
Data availability statement
Data are available on reasonable request.
Ethics statements
Patient consent for publication:
Not applicable.
Ethics approval:
The prospective observational cohort study conducted in the Royal Free London NHS Foundation Trust was approved by Harrow Research Ethics Committee (16/LO/1490). Participants gave informed consent to participate in the study before taking part.
Acknowledgements
We would like to thank Sary Workman and Fernando Moreira in the Department of Clinical Immunology Royal Free London NHS Foundation Trust for their help and support in the recruitment and samples' collection. The first author would like to thank the Missions sector in the Egyptian Ministry of Higher Education and the British Council Egypt for funding her PhD as part of the Newton-Mosharafa scheme. The funder had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
WHO. Forum of international respiratory societies. The global impact of respiratory disease. Sheffield, European Respiratory Society2017; Google Scholar
GOLD. Global strategy for prevention, diagnosis and management of COPD: 2022 report. 2022; Google Scholar
Donaldson GC, Seemungal TAR, Bhowmik A, et al. Relationship between exacerbation frequency and lung function decline in chronic obstructive pulmonary disease. Thorax2002; 57:847–52. doi:10.1136/thorax.57.10.847•Google Scholar
Martínez-García MA, Soler-Cataluña J-J, Perpiñá-Tordera M, et al. Factors associated with lung function decline in adult patients with stable non-cystic fibrosis bronchiectasis. Chest2007; 132:1565–72. doi:10.1378/chest.07-0490•Google Scholar
Ni W, Shao X, Cai X, et al. Prophylactic use of macrolide antibiotics for the prevention of chronic obstructive pulmonary disease exacerbation: a meta-analysis. PLoS ONE2015; 10. doi:10.1371/journal.pone.0121257•Google Scholar
Hnin K, Nguyen C, Carson KV, et al. Prolonged antibiotics for non-cystic fibrosis bronchiectasis in children and adults. Cochrane Database Syst Rev2015; 2015. doi:10.1002/14651858.CD001392.pub3•Google Scholar
Girard D, Finegan SM, Dunne MW, et al. Enhanced efficacy of single-dose versus multi-dose azithromycin regimens in preclinical infection models. J Antimicrob Chemother2005; 56:365–71. doi:10.1093/jac/dki241•Google Scholar
Wenzel RP, Fowler AA, Edmond MB, et al. Antibiotic prevention of acute exacerbations of COPD. N Engl J Med2012; 367:340–7. doi:10.1056/NEJMct1115170•Google Scholar
EUCAST. The European committee on antimicrobial susceptibility testing. Breakpoint tables for interpretation of mics and zone diameter. Version 7.0. 2017;
Stressmann FA, Rogers GB, Marsh P, et al. Does bacterial density in cystic fibrosis sputum increase prior to pulmonary exacerbation? J Cyst Fibros2011; 10:357–65. doi:10.1016/j.jcf.2011.05.002•Google Scholar
Garcha DS, Thurston SJ, Patel ARC, et al. Changes in prevalence and load of airway bacteria using quantitative PCR in stable and exacerbated COPD. Thorax2012; 67:1075–80. doi:10.1136/thoraxjnl-2012-201924•Google Scholar
Rofael SAD. Airway microbiome in chronic lung disease. UCL (University College London)2021; Google Scholar
Rofael SAD, Brown J, Pickett E, et al. Enrichment of the airway microbiome in people living with HIV with potential pathogenic bacteria despite antiretroviral therapy. EClinicalMedicine2020; 24:100427. doi:10.1016/j.eclinm.2020.100427•Google Scholar
Comeau AM, Douglas GM, Langille MGI, et al. Microbiome helper: a custom and streamlined workflow for microbiome research. MSystems2017; 2. doi:10.1128/mSystems.00127-16•Google Scholar
Caporaso JG, Kuczynski J, Stombaugh J, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods2010; 7:335–6. doi:10.1038/nmeth.f.303•Google Scholar
White JR, Nagarajan N, Pop M, et al. Statistical methods for detecting differentially abundant features in clinical metagenomic samples. PLoS Comput Biol2009; 5. doi:10.1371/journal.pcbi.1000352•Google Scholar
Petrie A. Medical statistics at a glance. Newark, John Wiley & Sons, Incorporated2019; Google Scholar
Charalampous T, Kay GL, Richardson H, et al. Nanopore metagenomics enables rapid clinical diagnosis of bacterial lower respiratory infection. Nat Biotechnol2019; 37:783–92. doi:10.1038/s41587-019-0156-5•Google Scholar
Bhandari V, Gupta RS. Molecular signatures for the phylum Synergistetes and some of its subclades. Antonie Van Leeuwenhoek2012; 102:517–40. doi:10.1007/s10482-012-9759-2•Google Scholar
Rogers GB, Bruce KD, Martin ML, et al. 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 Med2014; 2:988–96. doi:10.1016/S2213-2600(14)70213-9•Google Scholar
Albert RK, Connett J, Bailey WC, et al. Azithromycin for prevention of exacerbations of COPD. N Engl J Med2011; 365:689–98. doi:10.1056/NEJMoa1104623•Google Scholar
Uzun S, Djamin RS, Kluytmans J, et al. Azithromycin maintenance treatment in patients with frequent exacerbations of chronic obstructive pulmonary disease (COLUMBUS): a randomised, double-blind, placebo-controlled trial. Lancet Respir Med2014; 2:361–8. doi:10.1016/S2213-2600(14)70019-0•Google Scholar
Berkhof FF, Doornewaard-ten Hertog NE, Uil SM, et al. Azithromycin and cough-specific health status in patients with chronic obstructive pulmonary disease and chronic cough: a randomised controlled trial. Respir Res2013; 14. doi:10.1186/1465-9921-14-125•Google Scholar
Venditto VJ, Haydar D, Abdel-Latif A, et al. Immunomodulatory effects of azithromycin revisited: potential applications to COVID-19. Front Immunol2021; 12. doi:10.3389/fimmu.2021.574425•Google Scholar
Köhler T, Perron GG, Buckling A, et al. Quorum sensing inhibition selects for virulence and cooperation in pseudomonas aeruginosa. PLoS Pathog2010; 6. doi:10.1371/journal.ppat.1000883•Google Scholar
Burr LD, Rogers GB, Chen A-H, et al. Macrolide treatment inhibits pseudomonas aeruginosa quorum sensing in non-cystic fibrosis bronchiectasis. An analysis from the bronchiectasis and low-dose erythromycin study trial. Ann Am Thorac Soc2016; 13:1697–703. doi:10.1513/AnnalsATS.201601-044OC•Google Scholar
Malhotra-Kumar S, Lammens C, Martel A, et al. Oropharyngeal carriage of macrolide-resistant viridans group streptococci: a prevalence study among healthy adults in Belgium. J Antimicrob Chemother2004; 53:271–6. doi:10.1093/jac/dkh026•Google Scholar
Perez-Trallero E. High proportion of pharyngeal carriers of commensal streptococci resistant to erythromycin in Spanish adults. J Antimicrob Chemother2001; 48:225–9. doi:10.1093/jac/48.2.225•Google Scholar
Altenburg J, de Graaff CS, Stienstra Y, et al. Effect of azithromycin maintenance treatment on infectious exacerbations among patients with non-cystic fibrosis bronchiectasis: the bat randomized controlled trial. JAMA2013; 309:1251–9. doi:10.1001/jama.2013.1937•Google Scholar
Brill SE, Law M, El-Emir E, et al. Effects of different antibiotic classes on airway bacteria in stable COPD using culture and molecular techniques: a randomised controlled trial. Thorax2015; 70:930–8. doi:10.1136/thoraxjnl-2015-207194•Google Scholar
Serisier DJ, Martin ML, McGuckin MA, et al. Effect of long-term, low-dose erythromycin on pulmonary exacerbations among patients with non-cystic fibrosis bronchiectasis: the BLESS randomized controlled trial. JAMA2013; 309:1260–7. doi:10.1001/jama.2013.2290•Google Scholar
Taylor S, Leong L, Mobegi F, et al. The impact of azithromycin therapy on the lung microbiota and resistome in the AMAZES trial cohort. Respirology2019; 24. Google Scholar
Stadler C, Teuber M. The macrolide efflux genetic assembly of Streptococcus pneumoniae is present in erythromycin-resistant streptococcus salivarius. Antimicrob Agents Chemother2002; 46:3690–1. doi:10.1128/AAC.46.11.3690-3691.2002•Google Scholar
Seemungal TAR, Wilkinson TMA, Hurst JR, et al. Long-term erythromycin therapy is associated with decreased chronic obstructive pulmonary disease exacerbations. Am J Respir Crit Care Med2008; 178:1139–47. doi:10.1164/rccm.200801-145OC•Google Scholar
Sethi S, Jones PW, Theron MS, et al. Pulsed moxifloxacin for the prevention of exacerbations of chronic obstructive pulmonary disease: a randomized controlled trial. Respir Res2010; 11. doi:10.1186/1465-9921-11-10•Google Scholar
Hare KM, Grimwood K, Chang AB, et al. Nasopharyngeal carriage and macrolide resistance in Indigenous children with bronchiectasis randomized to long-term azithromycin or placebo. Eur J Clin Microbiol Infect Dis2015; 34:2275–85. doi:10.1007/s10096-015-2480-0•Google Scholar