Short-term increase in prevalence of nasopharyngeal carriage of macrolide-resistant Staphylococcus aureus following mass drug administration with azithromycin for trachoma control

Background Mass drug administration (MDA) with azithromycin is a corner-stone of trachoma control however it may drive the emergence of antimicrobial resistance. In a cluster-randomized trial (Clinical trial gov NCT00792922), we compared the reduction in the prevalence of active trachoma in communities that received three annual rounds of MDA to that in communities that received a single treatment round. We used the framework of this trial to carry out an opportunistic study to investigate if the increased rounds of treatment resulted in increased prevalence of nasopharyngeal carriage of macrolide-resistant Staphylococcus aureus. Three cross-sectional surveys were conducted in two villages receiving three annual rounds of MDA (3 × treatment arm). Surveys were conducted immediately before the third round of MDA (CSS-1) and at one (CSS-2) and six (CSS-3) months after MDA. The final survey also included six villages that had received only one round of MDA 30 months previously (1 × treatment arm). Results In the 3 × treatment arm, a short-term increase in prevalence of S. aureus carriage was seen following MDA from 24.6% at CSS-1 to 38.6% at CSS-2 (p < 0.001). Prevalence fell to 8.8% at CSS-3 (p < 0.001). A transient increase was also seen in prevalence of carriage of azithromycin resistant (AzmR) strains from 8.9% at CSS-1 to 34.1% (p < 0.001) in CSS-2 and down to 7.3% (p = 0.417) in CSS-3. A similar trend was observed for prevalence of carriage of macrolide-inducible-clindamycin resistant (iMLSB) strains. In CSS-3, prevalence of carriage of resistant strains was higher in the 3 × treatment arm than in the 1 × treatment (AzmR 7.3% vs. 1.6%, p = 0.010; iMLSB 5.8% vs. 0.8%, p < 0.001). Macrolide resistance was attributed to the presence of msr and erm genes. Conclusions Three annual rounds of MDA with azithromycin were associated with a short-term increase in both the prevalence of nasopharyngeal carriage of S. aureus and prevalence of carriage of AzmR and iMLSB S. aureus. Trial registration This study was ancillary to the Partnership for the Rapid Elimination of Trachoma, ClinicalTrials.gov NCT00792922, registration date November 17, 2008.


Background
Trachoma, caused by ocular infection with the intracellular bacterium Chlamydia trachomatis, is the leading infectious cause of blindness worldwide. Mass drug administration (MDA) with the broad-spectrum antibiotic azithromycin is an important part of the World Health Organization's integrated strategy for trachoma control [1,2]. This treatment serves to decrease the reservoir of infection, thereby reducing transmission.
There has been increased interest in MDA with azithromycin following the publication of studies conducted in Ethiopia suggesting treatment is associated with a significant reduction in childhood mortality [3,4]. Calls to expand azithromycin distribution beyond trachoma-endemic countries [5] and a large-scale clinical trial to evaluate the effect of treatment on mortality that is underway in three African countries [6] bring a greater need to document unintended effects of treatment, including the emergence of antimicrobial resistance, which is a global public health concern.
There is, as yet, no evidence to suggest MDA of azithromycin at the community-level leads to increased azithromycin resistance in ocular Chlamydia trachomatis infection [7][8][9]. However, there are data supporting an association of MDA with the emergence of macrolide-resistant Streptococcus pneumoniae isolated from the nasopahrnyx, at least in some settings. While studies carried out in Tanzania, Nepal and The Gambia have shown no evidence of such resistance following a single treatment round [10][11][12], other studies in Tanzania, Nepal and Australia suggest resistance does emerge after just one or two annual rounds of mass treatment [13][14][15]. Further studies in Ethiopia have documented increased macrolide resistant pneumococci isolated following four rounds of MDA given at 3 month intervals [16] and following six biannual rounds over a period of 3 years [17]. To date, little work has been carried out to assess the effect of MDA with azithromycin on other bacterial pathogens.
Staphylococcus aureus colonization is a risk factor for many conditions ranging from skin and soft tissue infections in children to invasive disease such as neonatal sepsis, bacteraemia and endocarditis [18][19][20][21][22]. In West Africa, it has been shown to be a significant cause of invasive disease in young children [23,24]. However, the potential effect of azithromycin MDA on prevalence of carriage of S. aureus including macrolide resistant strains has not yet been explored. In the present study, we used the framework of a clinical trial, which compared the prevalence of active trachoma in Gambian communities that received three annual rounds of azithromycin MDA to that of communities that received a single MDA round [25,26] to explore whether three rounds of MDA drove increased nasopharyngeal carriage of azithromycin resistant S. aureus.

Study design
The Partnership for the Rapid Elimination of Trachoma (PRET) study (ClinicalTrials.gov NCT00792922) was a cluster randomized controlled trial, the design of which has been described elsewhere [25,26]. Briefly, the study compared the effectiveness of three versus one round of azithromycin MDA in reducing the prevalence of active trachoma and ocular C. trachomatis infection. Treatment was a single oral dose of 20 mg azithromycin per kg to a maximum of 1 g and height was used as a proxy for weight. A pneumococcal carriage study was nested within PRET [12] and was carried out in eight villages that were a part of the larger trial. This included two villages that had been randomized, by the underlying PRET trial, to three annual rounds of MDA (3 × treatment arm) and six villages that received a single treatment round (1 × treatment arm) (Figs. 1 and 2). All villages had also participated in a trial of pneumococcal conjugate vaccine and were part of that study's control arm, where children under 5 years of age received PCV-7 [27].
Three cross-sectional surveys (CSS) were conducted in the 3 × treatment arm: CSS-1, 1 month prior to the third round of MDA; CSS-2, 1 month following the third round of MDA; and CSS-3, 6 months following the third round (Figs. 1 and 2). CSS-3 also included villages in the 1 × treatment arm that had received one round of MDA 30 months previously (Figs. 1 and 2).
Census data were gathered in the week before CSS-1. All censused children under the age of 15 years and present at the time of sampling were invited to participate. For individuals 15 years of age or older, 150 individuals were randomly selected at each of CSS-1 and CSS-2 and in both treatment arms in CSS-3. Random selection was independent at each CSS. In the 3 × treatment arm, participation in CSS-2 and 3 was restricted to those who were documented to have received azithromycin during the third round of MDA in July 2010. At each CSS, a nasopharyngeal swab (NPS) was collected from each participant and a questionnaire, pertaining to demographic characteristics, socioeconomic status and recent antibiotic use (within the last 30 days), was conducted. In the present study, we analysed all NPS samples collected as a part of the pneumococcal carriage study [12] to determine the prevalence of S. aureus carriage, including carriage of macrolide resistant strains.

Specimen collection
NPS samples were collected using calcium alginate swabs and inoculated into sterile vials containing skimmilk-tryptone-glucose glycerol transport medium, which has been shown to preserve the viability of respiratory pathogens for up to 12 years when stored at −70°C [28]. Samples were kept on wet-ice in the field, transferred to a 4°C refrigerator within 8 h of collection and moved to long-term storage at −80°C within 24 h of collection. Samples remained in long-term storage at −80°C with back-up power and twice daily temperature monitoring for 18 months before being cultured for S. aureus.

Laboratory methods
NPS samples were thawed at room temperature and 10 μl transport medium was inoculated onto Mannitol Salt Agar (Oxoid Ltd, Basingstoke, UK) and incubated for 18-24 h at 35°C. Presumptive S. aureus colonies were subcultured for purity on a Blood Agar Base No. 2 (Oxoid Ltd, Basingstoke, UK) supplemented with 5% sheep blood (TCS Biosciences Ltd., Botolph Clayton, UK) and incubated under the same conditions. A coagulase agglutination test was performed on well-isolated colonies to confirm the identification of S. aureus using the SLIDEX Staph-kit (Biomerieux UK Limited, Basingstoke, UK).
Sensitivity results were interpreted according to the Clinical and Laboratory Standards Institute guidelines [29]. Azithromycin resistance (Azm R ) was defined by a zone size ≤13 mm. Macrolide-inducible-clindamycin resistance, also referred to as the inducible Macrolide-Lincosamide-Streptogramin B (iMLS B ) phenotype, was defined as flattening of the circular zone of inhibition around the clindamycin disc on the side facing the erythromycin disc (D-test) [30]. Constitutive clindamycin resistance (Cli R ) was defined as a circular zone of inhibition ≤ 14 mm.
Twenty-three S. aureus isolates obtained from the 3 × treatment arm at CSS-1 and CSS-2 were randomly selected for analysis of the presence of 117 of the most prevalent and transferable antibiotic resistance determinants found in Gram-positive bacteria using micro-array, as previously described [31].
Laboratory technicians were blind to the treatment arm and the timing of sample collection.

Data management and statistical analysis
Data were double-entered into an OpenClinica database and the validated data imported into Stata version 12 (StataCorp LP, College Station, Texas, USA) for statistical analyses. Proportions of overall carriage and carriage of antimicrobial resistant S. aureus were compared using Chi-Square or Fisher's exact tests. Logistic regression analyses were further used to identify risk factors for carriage, control for confounders and test for interactions, reporting odds ratios with 95% confidence intervals (CI). CIs and pvalues were estimated using clustered robust standard errors to account for within village correlation of participants. Analysis included overall prevalence of nasopharyngeal carriage, prevalence of carriage of Azm R and iMLS B strains and proportions of Azm R and iMLS B strains isolated.

Ethical review
This study was ancillary to the PRET trial and was approved by The Gambia Government/Medical Research Council Unit, The Gambia Joint Ethics Committee. Written, informed consent was obtained from all participants of the PRET trial and additional written consent was obtained for this ancillary study. In the case of minors, informed consent was obtained from the parent or guardian.

Study participants
In the 3 × treatment arm, NPS were collected from 415, 417 and 343 individuals at CSS-1, CSS-2 and CSS-3, respectively. Four hundred participants in the 1 × treatment arm were also sampled at CSS-3. Overall, 25 NPS from the 1 × treatment arm (6.2% of 400 samples) had insufficient volume remaining to conduct the laboratory assays and were excluded from the analysis. Demographic and epidemiological characteristics of the participants, including S. pneumoniae carriage [12], are given in Table 1.

Antibiotic resistance determinants
In a subset of S. aureus isolates (N = 23), the presence of antibiotic resistance determinants was assayed using a DNA microarray [31]. Results are shown in Table 5. No macrolide resistant determinants were found in six isolates that were sensitive to both azithromycin and clindamycin. Of five isolates displaying resistance to azithromycin but sensitivity to clindamycin, all were positive for the msr gene, which conveys resistance to macrolides and streptogramin B. Twelve isolates were azithromycin resistant and had the iMLS B phenotype and all of these carried erm genes [11 erm(C) and 1 erm(T)], confirming their resistance to the MLS B antibiotics.

Discussion
In order to explore the effect of repeated MDA with azithromycin on the prevalence of carriage of macrolideresistant S. aureus and the proportion of resistant strains, we compared communities receiving one or three annual treatment rounds. Our results indicate that MDA was associated with a significant increase in the prevalence of carriage of Azm R and iMLS B S. aureus strains immediately following treatment, which returned to lower levels 6 months later. When comparing treatment arms at CSS-3, the prevalence of carriage of resistant S. aureus and proportion of resistant strains were higher in those communities that received three rounds of MDA as compared to a single treatment round.
The observation that the prevalence of resistant strains, in the 3 × treatment arm, decreased between CSS-2 and CSS-3 is consistent with research of macrolide resistant Streptococcus pneumoniae following azithromycin MDA that has demonstrated resistant phenotypes have decreased fitness when antibiotic pressure is relieved [32,33]. It may be that, had we collected additional samples at longer time points following treatment, the prevalence in the 3 × treatment arm may have eventually reached that seen in villages that received only a single round of treatment.
When we examine our results as the proportion of S. aureus isolates that display macrolide resistance, rather than the prevalence of carriage of resistant strains,  [34,35] and this is consistent with the timing of our surveys (CSS-1 and CSS2 were conducted in the wet season while CSS-3 was conducted in the dry season). However, other external factors may also have played a role. For example, conducting surveys in the wet season while crops are being planted may result in under representation of able-bodied adults amongst those surveyed. As the proportion of isolates displaying antimicrobial resistance was high, we chose to validate our findings on a subset of isolates using a microarray designed to detect antimicrobial resistance determinants common to Grampositive bacteria [31]. The results confirmed the presence of determinants encoding resistance to macrolide, lincosamide and streptogramin B antibiotics. While macrolides are not first-line therapy for S. aureus infection in Gambia (treatment would usually be cloxacillin or chloramphenicol), their use is indicated in respiratory disease in the case of penicillin allergy and recurrent skin infection, also in the case of penicillin allergy [36]. They would also be considered in the case of suspected atypical pneumonia. The presence of macrolide resistance therefore, while not a cause for immediate concern, is worth monitoring, especially as 91% of the isolates examined by microarray also carried the blaZ gene, suggesting concurrent resistance to penicillins in the population.
The majority of the resistance to macrolides was attributed to the presence of either msr or erm(C) genes however, one strain contained an erm(T) gene. To date, erm(T) has been primarily reported in Streptococcus species [37][38][39] and has rarely been identified in S. aureus isolates [40,41] suggesting it may have been acquired under selective pressure. The trimethoprim resistance gene dfr(G), which was detected in one third of our isolates is reported to be widespread among S. aureus isolates in Africa [42]. Almost half of the isolates that were tested by microarray carried the gene encoding Panton-Valentine leukocidin, a pore-forming cytotoxin that has been associated with skin and soft tissue infections and with community-acquired, necrotising haemorrhagic pneumonia [43,44].
S. pneumoniae colonization in the nasopharynx is thought, by many, to be negatively associated with S. aureus colonization and interventions to reduce pneumococcal carriage have been associated with an increase in S. aureus carriage and disease in some populations [45,46]. Plausible molecular mechanisms driving competition between the two bacteria are the pneumococcal pilus, which may allow better adherence of S. pneumoniae [47] and hydrogen peroxidase production by S. pneumoniae, which inhibits S. aureus growth [48]. In our study, S. aureus carriage was inversely associated with pneumococcal carriage at CSS-2, 1 month following MDA. One possible explanation for this is that the decrease in S. pneumoniae carriage immediately following MDA (Table 1) provided S. aureus the opportunity to expand its presence in the nasopharyngeal niche.
This was an opportunistic study that utilised the framework of the PRET trial [25,26] and the pneumococcal carriage study that was nested within PRET [12] to explore associations between carriage of macrolideresistant S. aureus and azithromycin MDA. As such, it has a number of limitations that could have been avoided had this been a prospective study of S. aureus carriage. One of these is the lack of baseline data, collected before any MDA. In communities that received just a single round of MDA (1 × treatment arm), 24% of strains isolated were resistant to azithromycin 30 months following treatment (at CSS-3). Data on carriage of macrolide resistant S. aureus in The Gambia are scarce, however, samples collected in 2003-2004 from the same region of country as our study (Western Division, now named Brikama Local Government Area) showed 64% erythromycin susceptibility [49,50]. This suggests macrolide resistance was no higher in our 1 × treatment arm at CSS-3 that it was in the region prior to azithromycin MDA. A survey of school-going children in Cameroon has also reported 75% susceptibility to erythromycin [51] suggesting our results are similar to levels of resistance in other West African countries. Another limitation was the use of nasopharyngeal swabs, rather than oropharyngeal or nasal swabs, which may have yielded higher numbers of S. aureus. However, there is no evidence to suggest the dynamics of macrolide resistance would differ between these respiratory sites. Sampling at additional time points post-MDA would also

Conclusions
Three rounds of MDA for trachoma control were associated an increase in the prevalence of carriage of Azm R and iMLS B S. aureus and in the proportion of isolated strains that were resistant to these antibiotics. While the increase in prevalence of carriage of resistant strains was transient, the increase in proportion of resistant strains was maintained for at least 6 months following the final round of MDA. Both the prevalence and the proportion of resistant strains was higher in the 3 × treatment arm than in the 1 × treatment arm. These findings highlight the need for continued antimicrobial resistance monitoring in communities receiving azithromycin treatment at the community-level.

Abbreviations
Azi R : Azithromycin resistance; CI: Confidence interval; Cli R : Clindamycin resistance; cMLS B : Constitutive Macrolide-Lincosamide-Streptogramin B phenotype; CSS: Cross sectional survey; iMLS B : Inducible Macrolide-Lincosamide-Streptogramin B phenotype; MDA: Mass drug administration; NPS: Nasopharyngeal swab; PRET: Partnership for the Rapid Elimination of Trachoma Bacteriology, University of Bern) for assistance with the microarray work and the community leaders and villagers for their participation in the study.

Funding
This work was supported by the Bill and Melinda Gates Foundation (grant number 48027) and by the Wellcome Trust (grant number WT093368MA). The funders had no role in the study design, data collection, analysis and interpretation of the data, decision to publish, or preparation of the manuscript.

Availability of data and materials
Data supporting the findings is contained within the manuscript.
Authors' contributions EB carried out the bacterial culture and antibiotic susceptibility testing and drafted the manuscript, JJ performed the statisitical analysis, VP carried out the microarray work and assisted in writing the manuscript, JH provided supervision for sample collection, EMH-E provided supervision in the field and assisted in writing the manuscript, AS provided supervision in the field, DCWM obtained funding and assisted in writing the manuscript, MJH obtained funding and assisted in writing the manuscript, RB obtained funding and assisted in writing the manuscript, AR was involved in study design, supervision in the field, interpretation of the data and writing of the manuscript, SEB was involved in study design, supervision in the field and laboratory, interpretation of the data and drafting of the manuscript. All authors read and approved the final manuscript.

Competing interests
The authors declare that they have no competing interests.

Consent to publish
Not applicable.

Ethics approval and consent to participate
This study was approved by The Gambia Government/Medical Research Council Unit, The Gambia Joint Ethics Committee under study numbers L2010.34, L2011.65 and L2013.32. Written, informed consent was obtained from all adult participants. In the case of minors, informed consent was obtained from the parent or guardian.

Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Author details