Azithromycin alters spatial and temporal dynamics of airway microbiota in idiopathic pulmonary fibrosis ======================================================================================================= * Pieter-Jan Gijs * Cécile Daccord * Eric Bernasconi * Martin Brutsche * Christian Clarenbach * Katrin Hostettler * Sabina A. Guler * Louis Mercier * Niki Ubags * Manuela Funke-Chambour * Christophe von Garnier ## Abstract **Background** High bacterial burden in lung microbiota predicts progression of idiopathic pulmonary fibrosis (IPF). Azithromycin is a macrolide antibiotic known to alter the lung microbiota in several chronic pulmonary diseases and observational studies have shown a positive effect of azithromycin on mortality and hospitalization rate in IPF. However, the effect of AZT on lung microbiota in IPF remain unknown. **Methods** We sought to determine the impact of a three-month course of azithromycin on lung microbiota in IPF. We assessed sputum and oropharyngeal swab specimens from 24 adults with IPF included in a randomized controlled cross-over trial of a thrice-weekly 500 mg oral azithromycin. 16S rRNA sequencing and quantitative polymerase chain reaction (qPCR) were performed to assess bacterial communities. Antibiotic resistance genes (ARG) were assessed using real-time qPCR. **Results** Azithromycin significantly decreased community diversity with a stronger and more persistent effect in lower airways. During treatment, turnover of airway microbiota decreased in upper and lower airways, resulting in greater similarity between microbiota of the two sites persisting one month after macrolide cessation. Patients with increased expression of ARG had a lower bacterial load and an enrichment of the genus *Streptococcus*. In contrast, patients without increased in ARG expression had a higher bacterial load and an enrichment in *Prevotella*. **Conclusions** We observed that AZT caused sustained changes in the diversity and composition of the upper and lower airway microbiota in IPF, with effects on the temporal and spatial dynamics between the two sites. Keywords * Microbiota * idiopathic pulmonary fibrosis * azithromycin * antibiotic resistance gene ## Introduction Idiopathic pulmonary fibrosis (IPF) is a progressive and fatal interstitial lung disease of unknown origin and the most common and severe form of the idiopathic interstitial pneumonias1,2. Active infection plays a role in IPF progression3. In healthy subjects, the composition of the microbiota in the upper respiratory tract (URT) and lower respiratory tract (LRT) has considerable similarities4. In the LRT, the influx of bacteria by microaspiration is counterbalanced by mucociliary clearance, resulting in a physiological turnover. Respiratory health is therefore associated with a dynamic turnover of the microbiota between the URT and the LRT. Converging evidence suggests that in IPF the spatial and temporal dynamics of the airway microbiota are disturbed. LRT microbiota is more abundant5,6 and less diverse7, suggesting either greater local growth, accumulation due to impaired clearance, or both (**Figure 1**). One consequence is greater dissimilarity between URT and LRT microbiota, with LRT microbiota carrying a majority of genes absent in the URT8. This, alongside with evidence from longitudinal data on persistence of LRT microbiota disturbances9, provides additional evidence for a lower microbiota turnover in IPF. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/09/15/2022.09.14.22279909/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2022/09/15/2022.09.14.22279909/F1) Figure 1. In health, the LRT microbiota is a community of transiently present microorganisms originating from the URT, whereas in IPF the LRT microbiota is a thriving community with a higher bacterial load and a decreased community diversity. The spatial dynamics between the URT and LRT microbiota in IPF, their turnover rates, and how these are affected by AZT remains unknown. *Abbreviations:* LRT = lower respiratory tract; URT = upper respiratory tract; IPF = idiopathic pulmonary fibrosis; AZT = azithromycin Preclinical data suggest that disruption of LRT microbiota precedes chronic lung epithelial injury and repair7,10, potentially perpetuating inflammation9. Bacterial load is currently the LRT microbiota feature most consistently associated with disease, in terms of progression, exacerbations and mortality. Bacterial load is higher in IPF patients than in healthy subjects, making it a potential therapeutic target11. A single-centre retrospective study showed good tolerance to prophylactic AZT in IPF patients, with fewer non-elective hospitalizations12. Also, treatment with AZT during acute IPF exacerbations improved survival rates compared to fluoroquinolones in a retrospective single-center study13. The impact of AZT on the respiratory microbiota in IPF is unknown. In patients with severe asthma and emphysema, multiple effects on the composition and structure of the LRT microbiota, as well as increased anti-inflammatory bacterial metabolites have been reported14,15. This effect of AZT on LRT microbiota was independent of a decrease in bacterial load. Acquired macrolide resistance is a global health concern and macrolide administration increases the carriage of macrolide-resistant bacteria in URT16,17 and LRT18. Higher carriage was observed after AZT for seven genes involved in antibiotic resistance, five of which associated with macrolide resistance, and two with tetracycline resistance14. Additionally, the airway resistome positively correlates with bacterial load in chronic obstructive pulmonary disease (COPD)19. This study examines the impact of a three-month treatment of AZT on the airway microbiota of IPF patients. ## Materials and Methods ### Study population and sample collection This study is a post hoc analysis of samples collected in the “*Azithromycin for the Treatment of Chronic Cough in Idiopathic Pulmonary Fibrosis”* study (Clinical trial identifier: [NCT02173145](http://medrxiv.org/lookup/external-ref?link_type=CLINTRIALGOV&access_num=NCT02173145&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F15%2F2022.09.14.22279909.atom))20. This previous study was a multi-centre, double-blind, randomised, placebo-controlled cross-over trial to determine the effect of AZT on chronic cough in IPF patients (see Supplementary materials for key inclusion and exclusion criteria). Patients underwent a 12-week treatment period with oral AZT 500 mg 3 times per week, and a 12-week treatment period with placebo 3 times per week in random order. The two periods were separated by a 4-week washout period. Expectorated sputum and oropharyngeal swab (OPS) samples were collected both before and after AZT and placebo periods, and an additional specimen was collected after the 4-week washout period following the second treatment period. Samples were stored at -80°C until analysis. The study design and detailed sample collection are presented in Figure E1 in the Supplementary data. ### DNA extraction and 16S rRNA amplicon quantification Sputum and OPS samples were treated with dithiothreitol to homogenise the mucosal phase and DNA was extracted using the DNeasy UltraClean microbial kit (Qiagen, Hilden, Germany), with inclusion of a lysozyme digestion step (see details in Supplementary data). Negative controls (N = 11) underwent the same procedure and included blank swabs, blank sputum collection tubes, as well as blank DNA extractions (reagent control). To obtain a proxy of bacterial load, we determined the copy numbers of the 16S rRNA gene by qPCR using previously reported primers specific to panbacteria21 (see **Supplementary Table E1** for full sequences). Standard curves were obtained using purified amplicon products. ### Bacterial 16S rRNA amplicon sequencing Amplicon sequencing targeted the V1-V2 region of the 16S rRNA gene with primers F-27 and R-338 (see Supplementary Table E1 for full sequences and Supplementary Data for details). Amplification was performed using the Accuprime Taq DNA Polymerase High Fidelity kit (Invitrogen, Waltham, MA). No-template PCR reaction controls (N = 2) were included. Libraries were loaded onto an Illumina MiSeq using pairwise chemistry, generating 250 × 2 read lengths (Lausanne Genomic technologies facility, University of Lausanne, Switzerland). ### Analysis of antibiotic resistance gene carriage Quantification of antibiotic resistance genes (ARG) carriage targeting 23S ribosomal RNA methyltransferases (*erm*(B) and *erm*(F)), ATP-binding cassette ribosomal protection protein (*mel* and *msr*[E]), major facilitator superfamily antibiotic efflux pump (*mef*), and tetracycline-resistant ribosomal protection proteins (*tet*[M] and *tet*[W]) was performed on sputum specimens using dye-based (SsoAdvanced Universal SYBR Green, Bio-Rad) or probe-based real-time PCR assays, using primer pairs, probes and conditions previously described14. The copy numbers of resistance genes per sample were normalised relative to the copy numbers of the 16S rRNA gene. To obtain a synthetic picture of ARG carriage per sample, the counts obtained for each individual gene were scaled from 0 to 1 to provide equal importance to each gene, and the cumulative counts were reported. ### Bioinformatics and statistical analysis All analyses were performed in R version 4.1.022. Demultiplexing, removal of chimeric and short reads, single-base resolution of reads into amplicon sequence variants (ASVs) using the Divisive Amplicon Denoising Algorithm 2 (DADA2) algorithm23, and taxonomic annotation using the SILVA database24 were performed using a dedicated pipeline available at [https://github.com/chuvpne/dada2-pipeline](https://github.com/chuvpne/dada2-pipeline). Analyses were performed on a rarefied dataset at a sequencing depth of 11,741 (**Figure E2;** see Supplementary materials for details on processing and quality control). Multiple group comparisons were made using Kruskal Wallis test with Dunn’s post hoc test and Holm’s adjustment. The Wilcoxon signed-rank test was used to compare paired data. In all tests, we considered an alpha significance level of 0.05. To assess the difference in bacterial community composition between OPS and sputum samples during each phase of AZT treatment, we computed the Aitchison distance, corresponding to the Euclidean distance between the samples after centered log-ratio transformation25. Comparisons of proportions of ASVs represented only in one condition or shared between two conditions were assessed using the chi-square goodness-of-fit test. To test for differences in microbiota composition between the different phases of AZT treatment, we used a matrix built on unweighted UniFrac distance and conducted a Permutational Multivariate Analysis of Variance (PERMANOVA) with 999 permutations. ## Results ### Study cohort Of the 25 patients who completed randomization, 12 were initially randomized to intervention and 13 to placebo. 20 patients completed the study. The Consolidated Standards of Reporting Trials diagram is shown in **Figure E3** in Supplementary materials. 24 patients had at least one specimen, with a total of 67 sputum and 90 OPS specimens. **Table 1** shows patient demographics and clinical characteristics. View this table: [Table 1.](http://medrxiv.org/content/early/2022/09/15/2022.09.14.22279909/T1) Table 1. Demographics and clinical characteristics of study cohort20 *Abbreviations*: UIP = usual interstitial pneumonia; HRCT = high-resolution computed tomography; TLC = total lung capacity; FVC = forced vital capacity; FEV1 = forced vital capacity in 1 second; DLCO = diffusing capacity of the lung for carbon monoxide. Values are presented as frequency and percentage or mean and standard deviation (SD) Bacterial quantification, 16S rRNA gene sequencing and ARG analysis were performed on all specimens, except on 1 OPS specimen with insufficient DNA. 10 OPS specimens not meeting the 11’741 read threshold were excluded from further sequencing analysis (**Figure E4** in Supplementary materials). Finally, 67 sputum and 79 OPS specimens were analyzed. To investigate temporal changes in respiratory microbiota linked to AZT, we distinguished five treatment phases: *PreAZT* for specimens collected before the start of AZT treatment, available only in patients who started with placebo (N=35; sputum=14; OPS=21), *StartAZT* for specimens collected at the beginning of AZT (N=31; sputum=15; OPS=16), *EndAZT* for specimens collected at the end of AZT (N=23; sputum=12; OPS=11), *PostAZT_1 month* for specimens collected 1 month after the end of treatment (N=32; sputum=14; OPS=18) and *PostAZT>=3 months* for specimens collected 3 months or more after the end of AZT (N=25; sputum=12; OPS=13). ### The respiratory microbiota of IPF patients is distinct from environmental noise As our study used low volume specimens from body sites with low microbial biomass, it was vulnerable to environmental noise26. Nevertheless, we detected a markedly higher load in both LRT (*p* < 0.001; **Figure E5a**) and URT (*p* < 0.001; **Figure E5b**) specimens compared to procedural controls. Rank abundance analysis further showed that the dominant bacteria were mostly different in patient samples versus controls (**Figure E5c**) and allowed us to exclude from further analysis two ASVs dominant in controls (ASV6_*Cutibacterium*; ASV26_*Pseudomonas*). ### The upper and lower airway microbiota of IPF patients is altered after AZT treatment There was no difference in bacterial load between the different treatment phases in either LRT (*p* = 0.95) or URT (*p* = 0.22; **Figure 2**), but we observed a decrease in community richness after AZT, with a stronger and more persistent effect in LRT (*p* = <0.001) compared to URT (*p* = 0.008; **Figure 3a**). A decrease in bacterial phylogenetic diversity was also observed in both LRT (*p* = <0.001) and URT (*p* = 0.0093), without reverting to the pre-treatment level until the end of the observation period for LRT and more transiently for URT (**Figure 3b**). In contrast, evenness increased transiently after treatment in LRT only (*p* = 0.0021; **Figure 3c**). The Shannon diversity index reflects the effects of the AZT on richness (Figure 3d). In addition, we observed a transient increase in the dominance of core community ASVs in LRT (*p* = 0.014), with a non-significant trend in URT (*p* = 0.22; **Figure E6** Principal coordinate analysis (PCoA) based on unweighted UniFrac distance, which accounts for phylogeny showed a significant shift between pre- and post-treatment communities, both in LRT (*p* = 0.008; **Figure E7a**) and URT (*p* = 0.009; **Figure E7b**). ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/09/15/2022.09.14.22279909/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2022/09/15/2022.09.14.22279909/F2) Figure 2. qPCR determination of 16S rRNA gene copy number per microlitre of sputum or total OPS samples showing no variation in bacterial load during the different phases of the study. *Abbreviations:* qPCR = quantitative polymerase chain reaction; OPS = oropharyngeal swab ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/09/15/2022.09.14.22279909/F3/graphic-4.medium.gif) [](http://medrxiv.org/content/early/2022/09/15/2022.09.14.22279909/F3/graphic-4) ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/09/15/2022.09.14.22279909/F3/graphic-5.medium.gif) [](http://medrxiv.org/content/early/2022/09/15/2022.09.14.22279909/F3/graphic-5) Figure 3. Decrease in community diversity after AZT treatment. (*a* to *d*) Alpha diversity metrics showing a decrease in Chao1 richness (*a*) and Faith’s diversity (*b*) during AZT treatment in LRT and URT, an increase in Camargo evenness in LRT (*c*), and a decrease in Shannon diversity in LRT (*d*). *Abbreviations*: AZT = azithromycin; LRT = lower respiratory tract; URT = upper respiratory tract Together, these results show a broad impact of AZT on the respiratory microbiota of IPF patients. ### AZT treatment reduces spatial dissimilarity between URT and LRT microbiota and affects their temporal dynamics Based on unweighted UniFrac distance, bacterial communities in URT and LRT were distinct (**Figure 4a**). Per-patient comparisons showed greater similarity between URT and LRT samples at the end of AZT treatment, with the effect persisting one month after treatment with subsequent attenuation thereafter (**Figure 4b**). Divergence of the microbiota between the two sites was driven by ASVs present only in LRT, more numerous (*p* < 0.001) than those present only in URT or shared between the two sites, in each treatment phase. The transient decrease in the divergence between URT and LRT microbiota after AZT was associated with a non-significant decrease in the proportion of ASVs represented in LRT only, and a marginal increase in the number of ASVs shared by both sites (**Figure 4c**). ![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/09/15/2022.09.14.22279909/F4.medium.gif) [Figure 4.](http://medrxiv.org/content/early/2022/09/15/2022.09.14.22279909/F4) Figure 4. AZT treatment reduces dissimilarly between LRT and URT microbiota. (*a*) PCoA based on unweighted UniFrac distance showing an overall difference in bacterial community composition between URT and LRT. (*b*) Aitchison distance showing a decrease in dissimilarity between URT and LRT microbiota between the start and end of treatment. For each treatment phase, Aitchison distance was calculated between each OPS sample and the centroid obtained from the corresponding set of sputum samples. (*c*) Venn diagrams of the numbers and proportions of ASVs represented in either the LRT or URT alone, or in both sites, showing that the decrease in dissimilarity between the LRT and URT microbiota during treatment was in part driven by a decrease in the proportion of ASVs present in the LRT alone.(*d*) Venn diagrams of the numbers and proportions of ASVs represented either at the start, at the end or retained over a 3-month time window, showing a decrease in ASV turnover after AZT treatment in both LRT and URT. *Abbreviations*: AZT = azithromycin; LRT = lower respiratory tract; URT = upper respiratory tract; PCoA = principal coordinate analysis; OPS = oropharyngeal swab; ASVs = amplicon sequence variants Analysis of samples from IPF patients first receiving placebo allowed quantification of bacterial turnover over a three-month period in the absence of AZT. We found a mean of 40% ASVs retained in URT, and 45% in LRT (**Figure 4d, left diagrams**). During AZT treatment, patients had a slower turnover in airway microbiota, resulting in fewer ASVs present at the end compared to the start of AZT treatment in LRT (*p* < 0.001) and URT (*p* = 0.008). In LRT, the number of newly acquired ASVs was also lower than the number of retained ASVs (*p* = 0.004), in contrast to the observations in patients receiving placebo first (**Figure 4d, middle diagrams**). However, the bacterial turnover after treatment was increased compared to the placebo first group, with numbers of newly acquired ASVs exceeding those of cleared ASVs in both the LRT (*p* = 0.045) and URT (*p* = 0.026; **Figure 4d, right diagrams**). We observed that bacteria cleared from LRT during treatment, which accounted for 13 to 47% of the local community at the start of treatment, were already present in this site (between 7 and 24% relative abundance) at least four months before treatment started. The effect of treatment on these bacteria was long-lasting, with only partial resilience five months after the end of treatment. These same bacteria followed similar kinetics in URT, with the relative abundance ranging from 2-42% in the 4 months prior to the start of treatment and then decreasing to less than 0.05% (3 patients) or below detection (3 patients) (**Figure 5a**). ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/09/15/2022.09.14.22279909/F5/graphic-7.medium.gif) [](http://medrxiv.org/content/early/2022/09/15/2022.09.14.22279909/F5/graphic-7) ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/09/15/2022.09.14.22279909/F5/graphic-8.medium.gif) [](http://medrxiv.org/content/early/2022/09/15/2022.09.14.22279909/F5/graphic-8) Figure 5. Longitudinal changes in the relative abundance of bacteria cleared (*a*) or acquired (*b*) in the LRT during AZT treatment, showing different magnitudes of change in relative abundance between cleared and acquired bacteria, and similar kinetics between LRT and URT. *Abbreviations:* LRT = lower respiratory tract; AZT = azithromycin; URT = upper respiratory tract Bacteria acquired in LRT during treatment were represented in some cases at low levels during the four months prior to treatment. Their local relative abundance was highest at the end of treatment, without reaching the levels of bacteria cleared by the treatment, before decreasing during the 5 months following the end of treatment. These same bacteria were also poorly represented in URT prior to treatment (less than 1% in 7 of the 8 patients with data available). Their local relative abundance also increased during treatment, and, for some patients, these bacteria remained detectable in URT for up to 5 months after the end of treatment (**Figure 5b**). Taken together, these observations indicate marked changes in the spatial distribution and temporal dynamics of the airway microbiota of IPF patients secondary to AZT treatment. ### Changes in the composition of LRT microbiota correlate with antibiotic resistance gene carriage We next investigated whether alterations in respiratory microbiota during AZT treatment were related to ARG acquisition, focusing on LRT and previously described target genes14. We observed that ARG carriage was limited to a minority of LRT samples and ASVs (**Figure 6a)** and that almost all samples with ARG were collected at the end of AZT treatment or later (**Figure 6b**). This was confirmed by longitudinal within-patient analysis, which showed a peak in total ARG carriage at the end of AZT treatment, with substantial inter-individual differences in post-treatment kinetics (**Figure 6c**). ![Figure 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/09/15/2022.09.14.22279909/F6.medium.gif) [Figure 6.](http://medrxiv.org/content/early/2022/09/15/2022.09.14.22279909/F6) Figure 6. Acquisition of ARG during AZT treatment in LRT. (*a* and *b*) Canonical correlation analysis showing that ARG carriage was limited to a minority of sputum samples (circles) and ASVs (crosses) (*a*), and stratification by treatment phase showing that almost all samples with ARG were taken at the end of AZT treatment or later (*b*). (*c*) Within-patient monitoring of the cumulative sum of normalised copy numbers of the full set of ARG genes, focusing on patients with more than one sputum sample available, showing the increase in ARG carriage after the start of AZT treatment. *Abbreviations:* ARG = antibiotic resistance genes; AZT = azithromycin; LRT = lower respiratory tract; ASVs = amplicon sequence variants To investigate the implications for LRT microbiota of an increase in ARG carriage, we separated patients with LRT samples available at the start and end of AZT treatment (N = 10) into two groups of five patients according to the median change in the cumulative sum of seven pooled ARG (median fold-change = 3.8) (**Figure E8**). This revealed an association between ARG carriage and bacterial load, with a higher load in patients with stable resistance compared with those with increased resistance during AZT treatment (*p* < 0.001) (**Figure 7a**). However, there was no change in bacterial load between the different treatment phases, whether ARG carriage status was stable or increased (**Figure E9**). ![Figure 7.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/09/15/2022.09.14.22279909/F7.medium.gif) [Figure 7.](http://medrxiv.org/content/early/2022/09/15/2022.09.14.22279909/F7) Figure 7. Relationship between the difference in ARG carriage during AZT treatment, bacterial load and the 15 dominant genera of the LRT microbiota. (a) qPCR determination of bacterial load showing lower levels in patients with increased resistance. (b and c) Double PCoA showing the influence of *Prevotella_7* abundance in the compositional changes occurring during AZT treatment in patients with stable resistance, respectively *Streptococcus* abundance in patients with increased resistance. In panel c, the arrows link LRT samples collected from a single patient at the start and end of AZT treatment. *Abbreviations:* ARG = antibiotic resistance genes; AZT = azithromycin; LRT = lower respiratory tract; qPCR = quantitative polymerase chain reaction; PCoA = principal coordinate analysis We observed changes in the composition of LRT microbiota between the start and end of treatment, depending on resistance status. Specifically, double PCoA focused on the 15 most abundant genera showed that AZT treatment was associated with a shift in composition influenced by *Prevotella_7* abundance, in all patients with stable resistance and only in one patient with increased resistance. In contrast, the other four patients with increased ARG carriage showed compositional changes mainly driven by *Streptococcus* (**Figure 7b and c**). Genus abundance analysis during treatment confirmed an increase in *Prevotella*_7 and a decrease in *Streptococcus* in patients with stable resistance, and an increase in *Streptococcus* in 4 out of the 5 patients with increased resistance (**Figure E10a and b**). Finally, enrichment of *Prevotella*_7 in patients with stable resistance during treatment (*p* = 0.006; **Figure E10c**), and of *Streptococcus* in those with increased resistance (*p* = 0.002; **Figure E10d**), was also observed when all samples of these patients were considered, including those taken before the start or after the end of treatment. Together, these observations link the resistome to the composition of LRT microbiota, with a direct correlation with local bacterial load, but without change in lung function during the study period (**Figure E11**)28 ## Discussion This study provides valuable insight on the spatial and temporal distribution of the respiratory microbiota in IPF by comparing URT and LRT, microbial turnover, and the impact of AZT on these dynamics. Our findings suggest an attenuating effect of AZT on airway ecological disruption but also disturbance in LRT through a decrease in community diversity and an increase of potential pathogenic community members. These effects persisted up to five months after the end of treatment, possibly due to AZT retention in lung macrophages27. No changes in clinical status were observed during the study period20. URT and LRT communities in our cohort of IPF patients differed significantly, as in other chronic respiratory diseases28,29, with microbial richness in LRT exceeding that in URT, which contrasts with the situation in healthy subjects. Greater decrease in microbial richness in LRT versus URT during AZT treatment, paralleled by a diminishing phylogenetic dissimilarity between microbiota of the two sites, suggests that treatment decreases airway ecological disturbance. The absence of a concomitant decrease in bacterial load implies that the numerous taxa cleared were replaced by a smaller number of taxa already present and becoming more dominant. As increased dominance may alter the bacterial impact on host, it is often interpreted as being harmful, but little is known about the threshold for detrimental dominance, which depends on bacteria and the clinical context. Further evidence of AZT treatment effects was investigated by assessing bacterial temporal dynamics. In the absence of AZT, bacterial turnover was lower in the LRT compared to the URT. Combined with excessive microbial richness in LRT, this suggests a local persistence of bacteria potentially harmful to the host. AZT treatment had two effects on this turnover, namely by eliminating taxa (decrease in richness) and by preventing the acquisition of new taxa. However, these effects vanished progressively during the months following treatment. The shift in LRT microbiota composition associated with AZT, and the decrease in richness and phylogenetic diversity, mirrors the findings reported in asthma and emphysema14,15. These reported alterations were mainly driven by a mutually exclusive increase in relative abundance of the two predominant genera of LRT microbiota in IPF, *Streptococcus* (Gram-positive *Firmicutes*) and *Prevotella* (Gram-negative *Bacteroidetes*). Enrichment of these genera was consistent with the development of macrolide resistance previously described in cystic fibrosis30,31. The difference we observed with the increase in ARG carriage during AZT treatment when either of these genera predominated suggests that different mechanisms may be present. The increase in ARG carriage (between 0.5- and 11-fold for 12 weeks of AZT exposure) across the cohort was consistent with findings reported in asthma14. While *Streptococcus* enrichment was associated with a marked increase in ARG carriage (in 4 out of 4 patients), *Prevotella* enrichment was linked to more stable ARG carriage (in 5 out of 6 patients). *Prevotella* acquisition of AZT resistance could be due to carriage of genes not listed in our targets, although we included the most likely candidates based on a previous metagenomic study of AZT treatment in severe asthma14. A larger sample size will be necessary to strengthen conclusions and to perform in-depth analysis to determine whether carriage of genes involved in tetracycline versus macrolide resistance is associated with different bacteria. Although patients remained clinically stable during the study period, it is difficult to predict the impact of even a relatively transient predominance of *Streptococcus* or *Prevotella* on lung health. Both genera include commensal species that can act as opportunistic pathogens under permissive conditions. Some members of *Streptococcus* have been associated with IPF progression6. *Prevotella* are maintained at low levels in the healthy lung, where they establish a subclinical level of inflammation favoring immune surveillance4. However, their abundance is increased in IPF where they can become predominant32. Furthermore, *Prevotella* are involved in several inflammatory diseases33, including periodontitis, with a possible impact on the frequency of exacerbations in COPD34. Accordingly, preclinical evidence suggests that *Prevotella* predominance in the airways promotes pulmonary fibrosis through a mechanism involving IL-17B35. Decreased microbial diversity in IPF, such as observed with the predominance of *Streptococcus* or *Prevotella*, has previously been associated with increased alveolar concentrations of pro-inflammatory cytokines and profibrotic growth factors10. Further studies are required to determine the relationship between alterations of microbiota and gene expression in the airways during AZT treatment. Our study has limitations. First, we report an association between AZT and ARG carriage, but separate analysis of microbiota by amplicon sequencing and ARG carriage by qPCR does not allow us to identify bacteria that acquire resistance during treatment. Genome-wide metagenomic analyses would answer this question. Second, the absence of a healthy control group prevented us from assessing whether disruption of the airway microbiota by AZT is specific for IPF. Third, although active infection or antibiotic therapy 4 weeks prior to enrolment was an exclusion criterion for the study, a preconditioning of the respiratory microbiota that might have impacted the subsequent AZT effect in some patients cannot be excluded. In conclusion, we found that AZT alters the spatial and temporal dynamics of airway microbiota in IPF patients with a decline in richness and phylogenetic diversity in URT and LRT, without impact on bacterial load. The AZT impact was characterized by increased *Prevotella* or *Streptococcus* abundance and persisted with partial resilience five months after treatment. Also, increased ARG carriage was present in half of patients in the presence of *Streptococcus* predominance. This study provides novel insights into airway ecology disturbances in IPF by longitudinal sampling of URT and LRT and expands on previous knowledge from studies in distal LRT. From a clinical care perspective this approach also encourages non-invasive sampling of more than one respiratory tract site, particularly in patients with IPF for whom bronchoscopic sampling is frequently associated with an unfavorable risk-benefit ratio. ## Supporting information Supplemental Material [[supplements/279909_file03.docx]](pending:yes) ## Data Availability All data produced in the present study are available upon reasonable request to the authors [https://github.com/chuvpne/dada2-pipeline](https://github.com/chuvpne/dada2-pipeline). [https://github.com/CHUVpulmonology/Airway\_microbiota-Lung\_fibrosis-Azithromycin](https://github.com/CHUVpulmonology/Airway_microbiota-Lung_fibrosis-Azithromycin) ## Acknowledgements We thank the patients who participated in the study and the study nurses from Bern, Basel, St. Gallen and Zurich for their work and data acquisition. We also thank the Lung Association of the Canton of Vaud (Ligue Pulmonaire Vaudoise, LPV) and the Research Fund of the Swiss Lung Association, Bern for financial support. 16S rRNA sequencing was performed at the Lausanne Genomic Technologies Facility, University of Lausanne, Switzerland ([https://www.unil.ch/gtf/en/home.html](https://www.unil.ch/gtf/en/home.html)). ## Footnotes * * Shared first authorship * + Shared last authorship * Received September 14, 2022. * Revision received September 14, 2022. * Accepted September 15, 2022. * © 2022, Posted by Cold Spring Harbor Laboratory The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission. ## References 1. 1.Maher TM, Wells AU, Laurent GJ. Idiopathic pulmonary fibrosis: multiple causes and multiple mechanisms? Eur Respir J 30, 835–839 (2007). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiZXJqIjtzOjU6InJlc2lkIjtzOjg6IjMwLzUvODM1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDkvMTUvMjAyMi4wOS4xNC4yMjI3OTkwOS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 2. 2.Raghu, G, Remy-Jardin M, Myers JL et al. Diagnosis of Idiopathic Pulmonary Fibrosis An Official ATS/ERS/JRS/ALAT Clinical Practice Guideline. Am J Respir Crit Care Med 198, 44–68 (2018). 3. 3.Invernizzi R, Molyneaux PL. The contribution of infection and the respiratory microbiome in acute exacerbations of idiopathic pulmonary fibrosis. Eur Respir Rev 28, (2019). 4. 4.Bassis CM, Erb-Downward JR, Dickson RP et al. Analysis of the upper respiratory tract microbiotas as the source of the lung and gastric microbiotas in healthy individuals. mBio 6, (2015). 5. 5.Molyneaux PL, Cox MJ, Willis-Owen SA et al. The role of bacteria in the pathogenesis and progression of idiopathic pulmonary fibrosis. Am J Respir Crit Care Med 190, 906–913 (2014). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1164/rccm.201403-0541OC&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25184687&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F15%2F2022.09.14.22279909.atom) 6. 6.Han MK, Zhou Y, Murray S et al. Lung microbiome and disease progression in idiopathic pulmonary fibrosis: An analysis of the COMET study. Lancet Respir Med 2, 548–556 (2014). 7. 7.Takahashi Y, Saito A, Chiba H et al. Impaired diversity of the lung microbiome predicts progression of idiopathic pulmonary fibrosis. Respir Res 19, 1–10 (2018). 8. 8.Tong X, Su F, Xu X et al. Alterations to the lung microbiome in idiopathic pulmonary fibrosis patients. Front Cell Infect Microbiol 9, 149 (2019). 9. 9.Molyneaux PL, Willis-Owen SAG, Cox MJ et al. Host-Microbial Interactions in Idiopathic Pulmonary Fibrosis. Am J Respir Crit Care Med 195, 1640–1650 (2017). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1164/rccm.201607-1408OC&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F15%2F2022.09.14.22279909.atom) 10. 10.O’Dwyer DN, Ashley SL, Gurczynski SJ et al. Lung Microbiota Contribute to Pulmonary Inflammation and Disease Progression in Pulmonary Fibrosis. Am J Respir Crit Care Med 199, 1127–1138 (2019). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1164/rccm.201809-1650OC&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F15%2F2022.09.14.22279909.atom) 11. 11.Invernizzi R, Barnett J, Rawal B et al. Bacterial burden in the lower airways predicts disease progression in idiopathic pulmonary fibrosis and is independent of radiological disease extent. Eur Respir J 55, (2020). 12. 12.Macaluso C, Maritano Furcada J, Alzaher O et al. The potential impact of azithromycin in idiopathic pulmonary fibrosis. Eur Respir J 53, (2019). 13. 13.Kawamura K, Ichikado K, Yasuda Y et al. Azithromycin for idiopathic acute exacerbation of idiopathic pulmonary fibrosis: A retrospective single-center study. BMC Pulm Med 17, (2017). 14. 14.Taylor SL, Leong LEX, Mobegi FM et al. Long-Term Azithromycin Reduces Haemophilus influenzae and Increases Antibiotic Resistance in Severe Asthma. Am J Respir Crit Care Med. 200(3):309–317 (2019). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1164/rccm.201809-1739OC&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F15%2F2022.09.14.22279909.atom) 15. 15.Segal LN, Clemente JC, Wu BG et al. Randomised, double-blind, placebo-controlled trial with azithromycin selects for anti-inflammatory microbial metabolites in the emphysematous lung. Thorax 72, 13–22 (2017). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToidGhvcmF4am5sIjtzOjU6InJlc2lkIjtzOjc6IjcyLzEvMTMiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wOS8xNS8yMDIyLjA5LjE0LjIyMjc5OTA5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 16. 16.Kastner U, Guggenbichler JP. Influence of macrolide antibiotics on promotion of resistance in the oral flora of children. Infection 29, 251–256 (2001). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s15010-001-1072-3&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11688901&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F15%2F2022.09.14.22279909.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000171439500002&link_type=ISI) 17. 17.Choo JM, Abell GCJ, Thomson R et al. Impact of Long-Term Erythromycin Therapy on the Oropharyngeal Microbiome and Resistance Gene Reservoir in Non-Cystic Fibrosis Bronchiectasis. mSphere 3, (2018). 18. 18.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. JAMA 309, 1251–1259 (2013). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2013.1937&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23532241&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F15%2F2022.09.14.22279909.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000316625200025&link_type=ISI) 19. 19.Ramsheh MY, Haldar K, Bafadhel M et al. Resistome analyses of sputum from COPD and healthy subjects reveals bacterial load-related prevalence of target genes. Thorax 75, 8–16 (2020). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToidGhvcmF4am5sIjtzOjU6InJlc2lkIjtzOjY6Ijc1LzEvOCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzA5LzE1LzIwMjIuMDkuMTQuMjIyNzk5MDkuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 20. 20.Guler SA, Clarenbach C, Brutsche M et al. Azithromycin for the Treatment of Chronic Cough in Idiopathic Pulmonary Fibrosis: A Randomized Controlled Crossover Trial. Ann Am Thorac Soc 18, 2018–2026 (2021). 21. 21.Bacchetti De Gregoris T, Aldred N, Clare AS et al. Improvement of phylum-and class-specific primers for real-time PCR quantification of bacterial taxa. J Microbiol Methods 86, 351–356 (2011). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.mimet.2011.06.010&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21704084&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F15%2F2022.09.14.22279909.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000294507400013&link_type=ISI) 22. 22.R Core Team (2021). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 23. 23.Callahan BJ, McMurdie PJ, Rosen MJ et al. DADA2: High-resolution sample inference from Illumina amplicon data. Nature Methods 2016 13:7 13, 581–583 (2016). 24. 24.Quast C, Pruesse E, Yilmaz P et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Research 41, D590 (2013). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gks1219&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23193283&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F15%2F2022.09.14.22279909.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000312893300084&link_type=ISI) 25. 25.Aitchison, J, Barceló-Vidal C, Martín-Fernández JA et al. Log ratio Analysis and Compositional Distance. Mathematical Geology 2000 32:3 32, 271–275 (2000). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1023/A:1007529726302&link_type=DOI) [GeoRef](http://medrxiv.org/lookup/external-ref?access_num=2009075610&link_type=GEOREF) 26. 26.Salter SJ, Cox MJ, Turek EM et al. Reagent and laboratory contamination can critically impact sequence-based microbiome analyses. BMC Biol 12, 1–12 (2014). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1741-7007-12-1&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24417977&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F15%2F2022.09.14.22279909.atom) 27. 27.Aubert JD, Juillerat-Jeanneret L, Fioroni P et al. Function of human alveolar macrophages after a 3-day course of azithromycin in healthy volunteers. Pulm Pharmacol Ther 11, 263–269 (1998). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1006/pupt.1998.0123&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10101742&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F15%2F2022.09.14.22279909.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000078462600003&link_type=ISI) 28. 28.Willner D, Haynes MR, Furlan M et al. Spatial distribution of microbial communities in the cystic fibrosis lung. ISME J 6, 471–474 (2012). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ismej.2011.104&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21796216&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F15%2F2022.09.14.22279909.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000300984200023&link_type=ISI) 29. 29.Venkataraman A, Bassis CM, Beck JM et al. Application of a neutral community model to assess structuring of the human lung microbiome. mBio 6, (2015). 30. 30.Webb KA, Olagoke O, Baird T et al. Genomic diversity and antimicrobial resistance of Prevotella species isolated from chronic lung disease airways. Microb Genom 8, 000754 (2022). 31. 31.Tazumi A, Maeda Y, Goldsmith CE et al. Molecular characterization of macrolide resistance determinants [erm(B) and mef(A)] in Streptococcus pneumoniae and viridans group streptococci (VGS) isolated from adult patients with cystic fibrosis (CF). J Antimicrob Chemother 64, 501–506 (2009). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/jac/dkp213&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19584106&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F15%2F2022.09.14.22279909.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000273181400006&link_type=ISI) 32. 32.Dickson RP, Huffnagle GB, Flaherty KR et al. Radiographic honeycombing and altered lung microbiota in patients with idiopathic pulmonary fibrosis. Am J Respir Crit Care Med 200, 1544–1547 (2019). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F15%2F2022.09.14.22279909.atom) 33. 33.Larsen JM. The immune response to Prevotella bacteria in chronic inflammatory disease. Immunology 151, 363–374 (2017). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/imm.12760&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28542929&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F15%2F2022.09.14.22279909.atom) 34. 34.Kelly N, Winning L, Irwin C et al. Periodontal status and chronic obstructive pulmonary disease (COPD) exacerbations: a systematic review. BMC Oral Health 21, 1–11 (2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12903-020-01373-3&link_type=DOI) 35. 35.Yang D, Chen X, Wang J et al. Dysregulated Lung Commensal Bacteria Drive Interleukin-17B Production to Promote Pulmonary Fibrosis through Their Outer Membrane Vesicles. Immunity 50, 692-706.e7 (2019).