The Fitness Cost of Antibiotic Resistance in Streptococcus pneumoniae: Insight from the Field

doi:10.1371/journal.pone.0029407

PLoS One. 2012; 7(1): e29407.

Published online 2012 January 17. doi: 10.1371/journal.pone.0029407

PMCID: PMC3260144

Copyright Maher et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

The Fitness Cost of Antibiotic Resistance in Streptococcus pneumoniae: Insight from the Field

M. Cyrus Maher,¹ Wondu Alemayehu,² Takele Lakew,² Bruce D. Gaynor,^3,⁴ Sara Haug,³ Vicky Cevallos,³ Jeremy D. Keenan,^3,⁴ Thomas M. Lietman,^1,^3,^4,⁵ and Travis C. Porco^1,^3,^4,⁵^*

¹Department of Epidemiology & Biostatistics, University of California San Francisco, San Francisco, California, United States of America

²ORBIS International, New York, New York, United States of America

³F.I. Proctor Foundation, University of California San Francisco, San Francisco, California, United States of America

⁴Department of Ophthalmology, University of California San Francisco, San Francisco, California, United States of America

⁵Institute for Global Health, University of California San Francisco, San Francisco, California, United States of America

Bernard Beall, Editor

Centers for Disease Control & Prevention, United States of America

* E-mail: travis.porco/at/ucsf.edu

Conceived and designed the experiments: TCP TML. Performed the experiments: BDG JDK WA TL SH VC. Analyzed the data: MCM TCP. Wrote the paper: MCM TCP TML BDG JDK.

Received August 13, 2011; Accepted November 28, 2011.

This article has been cited by other articles in PMC.

Abstract

Background

Laboratory studies have suggested that antibiotic resistance may result in decreased fitness in the bacteria that harbor it. Observational studies have supported this, but due to ethical and practical considerations, it is rare to have experimental control over antibiotic prescription rates.

Methods and Findings

We analyze data from a 54-month longitudinal trial that monitored pneumococcal drug resistance during and after biannual mass distribution of azithromycin for the elimination of the blinding eye disease, trachoma. Prescription of azithromycin and antibiotics that can create cross-resistance to it is rare in this part of the world. As a result, we were able to follow trends in resistance with minimal influence from unmeasured antibiotic use. Using these data, we fit a probabilistic disease transmission model that included two resistant strains, corresponding to the two dominant modes of resistance to macrolide antibiotics. We estimated the relative fitness of these two strains to be 0.86 (95% CI 0.80 to 0.90), and 0.88 (95% CI 0.82 to 0.93), relative to antibiotic-sensitive strains. We then used these estimates to predict that, within 5 years of the last antibiotic treatment, there would be a 95% chance of elimination of macrolide resistance by intra-species competition alone.

Conclusions

Although it is quite possible that the fitness cost of macrolide resistance is sufficient to ensure its eventual elimination in the absence of antibiotic selection, this process takes time, and prevention is likely the best policy in the fight against resistance.

Introduction

Streptococcus pneumoniae is the leading cause of serious illness in children and adults worldwide [1]. Despite the availability of a vaccine to combat the disease, rates of colonization remain high, as vaccine-induced immunity often results in replacement of targeted strains as other pneumococcal serotypes fill the newly opened ecological niche [2]. This antigenic diversity renders eradication nearly impossible, emphasizing the importance of antibiotic treatment for control of invasive pneumococcal disease, and illuminating the necessity for understanding and predicting long-term trends in resistance.

The diversity and adaptability of S. pneumoniae is facilitated in large part by active DNA import and extensive genomic repeats that greatly increase the likelihood of intra- and interspecific homologous recombination [3]. Considering its high genomic plasticity, it is not surprising that substantial levels of drug resistance are observed in this organism in response to antibiotic treatment [4]. Yet, recent research has suggested that antibiotic resistance may come at a cost to bacteria harboring the trait, limiting fitness in the absence of antimicrobial drug selection [5].

Competition studies both in vitro [6], [7], [8] and in animal models [9], [10], [11] have found evidence for a competitive disadvantage from drug resistance in a range of bacterial species. At the epidemiologic level, several observational studies have demonstrated that lower antibiotic prescription rates are associated with reduced antibiotic resistance [12], [13]. Others have used correlations between resistance and prevalence to estimate fitness costs, either directly [14], [15] or by fitting a mathematical transmission model [16]. To our knowledge, however, no prior study has observed the development of antibiotic resistance upon introduction of a previously unused antibiotic class, and followed the decay of that resistance after antibiotic pressure is removed.

Results

Longitudinal trial data

Presented in Figure 1A is the village-level prevalence of azithromycin-resistant strains during the course of the study, with average prevalence superimposed. These average values have been presented previously by our group [17]. Briefly, average macrolide resistance rose from 28.2% after 4 treatments to 76.8% after the 6th treatment. Resistance then declined sharply to 30.6% at 12 months after the last antibiotic treatment, and finally to 20.8% one year after that.

Figure 1

The prevalence of pneumococcal azithromycin resistance by village in the Ethiopian trachoma elimination study.

To estimate baseline antibiotic resistance, we analyzed baseline data from eight of the control villages that were within a 5 mile radius of the intervention sites. Laboratory testing found only one resistant case among 112 pneumococcal isolates, a rate of approximately 0.9%. Combined with survey data showing that macrolides constituted approximately 0.5% of the antibiotic doses dispensed locally (B. Ayele, unpublished data), it is likely that baseline rates of azithromycin resistance were comparably low in intervention villages. For illustrative purposes, this assumption is included in Figure 1A with an error bar to represent its uncertainty, and an asterisk to distinguish it from directly observed data.

Figures 1B and 1C, respectively, show the proportion of azithromycin resistance due to the mefA/E and ermB phenotypes. We observed that pneumococcal strains harboring the mefA/E gene constituted a larger share of azithromycin-resistant samples than strains with the ermB phenotype. In addition, we observed that from month 42 to month 54, the overall decrease in azithromycin resistance was mostly due to the decrease in prevalence of mefA/E strains. In fact, in one village, we saw the measured prevalence of ermB strains increase from 0% to over 20% in one year. This suggests that reintroduction events played a role in long-term transmission dynamics, and in future studies, should be considered when feasible.

Fitness cost estimation

Our model calculated the fitness of the mefA/E phenotype to be 0.86 (95% CI 0.80 to 0.90) relative to the azithromycin-sensitive strain, and similarly to be 0.88 (95% CI 0.82 to 0.93) for the ermB phenotype (log likelihood An external file that holds a picture, illustration, etc.
Object name is pone.0029407.e001.jpg

An external file that holds a picture, illustration, etc.
Object name is pone.0029407.e001.jpg

). Table 1 demonstrates the sensitivity of these estimates to model inputs. Relative fitness estimates were insensitive to nearly all fixed parameters, with the exception of the duration of infection (1/ An external file that holds a picture, illustration, etc.
Object name is pone.0029407.e002.jpg

An external file that holds a picture, illustration, etc.
Object name is pone.0029407.e002.jpg

) and the antibiotic efficacy against drug sensitive strains. In both cases, however, relative fitness estimates remained significantly below one across the range of input values tested.

Table 1

Parameter ranges tested during univariate sensitivity analysis, and their influence on calculated relative fitness.

Probability of Strain Survival

To examine the effect of fitness cost on predicted long-term trends in resistance, we performed simulations extrapolating the probability of survival for each azithromycin-resistant strain, out to 8 years after the last treatment. Figure 2 demonstrates that, under the basecase point estimate for the relative fitness, we would expect less than a 5% chance of survival for both drug resistant strains 5 years after the last antibiotic administration. Using the less optimistic upper bound for this estimate, we find the probability of survival at this time point to be only slightly higher: 6% for the mefA/E strain, and 5.4% for the ermB strain.

Figure 2

The probability of survival versus time under the basecase model for A.) the mefA/E strain and B.) the ermB strain.

Discussion

Previous research has demonstrated that antibiotic resistance often comes with a fitness cost, the magnitude of which may depend on growth conditions, mechanism of resistance, genetic background, and the presence of compensatory mutations [5]. In vitro and animal model experiments are capable of analyzing the influence of these variables, but in humans, even quantifying an average fitness cost of drug resistance is made difficult by the increased complexity and decreased control introduced by working with field data. As a result, most epidemiological studies interested in the fitness cost of drug resistance in S. pneumoniae have looked at correlations between naturally occurring changes in prescription rates and the prevalence of corresponding antibiotic resistances [14], [18]. Although this has demonstrated a significant relationship between the two factors, lack of an underlying transmission model and relatively small changes in antibiotic consumption have hindered the estimation of relative fitness values.

In this study, we have used data gathered as part of a longitudinal trachoma elimination study to examine the fitness cost of macrolide resistance in S. pneumoniae. Observed macrolide resistance rates rose from 28.2% following four rounds of biannual mass antibiotic distribution, to 76.8% after two subsequent rounds, followed by a drop to 20.6% twenty-four months after the conclusion of the program. Fitting a mathematical model to these data, we estimated the relative fitness for strains harboring two different mechanisms of macrolide resistance, and showed how these estimates may depend on assumed transmission parameters and antibiotic effectiveness. Calculated fitness costs increased significantly with higher antibiotic efficacy against drug-sensitive strains, and longer durations of infection (1/ An external file that holds a picture, illustration, etc.
Object name is pone.0029407.e006.jpg

An external file that holds a picture, illustration, etc.
Object name is pone.0029407.e006.jpg

). These results agree with theory. Respectively, a larger available ecological niche, or fewer infection cycles, would both be expected to magnify the competitive disadvantage of drug resistance given identical data.

The clinical significance of these differing fitness costs can be measured, in part, by antibiotic-resistant strains' probability of survival versus time. Our model predicts that, in the base scenario, it may take up to 5 years for macrolide resistant strains of S. pneumoniae to reach high probabilities of elimination by competition alone. In fact, this number may be an underestimate of the true value, as our model does not allow for the stable coexistence of pneumococcal strains. Previous modeling studies have shown that coexistence can be explained by simultaneous carriage of multiple strains [19], a phenomenon we were unable to model due to a lack of data on superinfection. However, the decreased complexity of the model allowed us to follow the probability of all community-level infection states over time, better capturing the stochastic effects that are important in small populations.

It should also be noted that due to bias inherent in our phenotyping method, any individual strains with both mechanisms of resistance would be classified as ermB by our protocol. A study from a different area of Ethiopia observed these strains to constitute 6.6% of resistant samples after 4 rounds of azithromycin treatment [20]. This misclassification would most likely increase the estimated fitness cost for the ermB strain, although low prevalence of mefA/E+, ermB+ strains would attenuate this effect.

In addition, our model does not allow for strain reintroduction, or for the amelioration of fitness costs by compensatory mutation, an adaptation that has been observed to develop quickly in other organisms [5]. The rate of compensatory mutation is only one part of the picture, however. In pneumococcal, both mechanisms of macrolide resistance are mediated by entire genes rather than single nucleotide polymorphisms (SNPs). As a result, the genetic target for mutations that attenuate or even eliminate resistance is comparatively much larger, consequently making these events far more likely. In other words, in the absence of antibiotic selection against pneumococcus, we would expect to see two opposing effects not captured in our model. Mutation could reduce antibiotic resistance, as well as the fitness cost associated with it.

Notwithstanding uncertainty in the long-term prediction of resistance rates, the existence of an epidemiologically verified fitness cost of macrolide resistance in S. pneumoniae suggests that intraspecific competition may be exploited to at least partially reclaim antibiotic efficacy after the development of high levels of resistance. In practice, however, this process may be highly variable and take years even under the most optimistic scenarios. Although it is quite possible that the fitness cost of macrolide resistance is sufficient to ensure its eventual elimination in the absence of antibiotic selection, this process takes time, and prevention is likely the best policy in the fight against resistance.

Methods

The Trachoma Elimination Follow-up Study

The Trachoma Elimination Follow-up (TEF) study [21] is a longitudinal cohort study conducted in the Gurage Zone in Ethiopia, in which 40 villages were chosen from the Gurage Zone in Ethiopia, assessed for ocular chlamydial infection, randomized to one of 5 treatment protocols, and monitored over a 54 month period. Individuals aged 1 year or older were offered directly observed treatment with single-dose oral azithromycin (1 g in adults; 20 mg/kg in children). Pregnant women or those allergic to macrolides were ordered a six-week course of topical 1% tetracycline ointment to be applied twice daily to both eyes (not directly observed).

A subset of 8 villages was selected randomly selected from the biannual treatment arm for a more detailed pneumococcal antibiotic resistance investigation. During the follow-up, we collected pneumococcal samples from fifteen randomly chosen one-to-five year olds in each village (a total of 120 children were studied at each time point). Nasopharyngeal swabs were taken at 24, 36, 42, and 54 months. To estimate pneumococcal prevalence prior to treatment, 240 children (also aged 1–5 years old) were randomly selected and sampled from sixteen nearby villages prior to the enrollment of the villages in the treatment program.

Verbal consent was obtained from the guardians of all children participating in the study. Written consent was not sought due to the low literacy rate in the area. Verbal consent was approved by both the University of California, San Francisco Institutional Review Board (IRB), and the Ethiopian IRB (The Ethiopian Science and Technology Commission). Consent was documented on the encounter form by a check mark. Approval for the study as a whole was also obtained from the UCSF IRB and the Ethiopian IRB. The study was carried out in accordance with the Declaration of Helsinki.

Determining Strain Assignments

Fifteen samples were collected from each village at each time point, tested for resistance to a variety of antibiotics, and then classified as (1) azithromycin-sensitive (2) azithromycin-resistant, with mefA/E phenotype or (3) azithromycin-resistant, with ermB phenotype. All samples with mean inhibitory concentrations (MIC) of azithromycin less than 1 mg/L where categorized as drug-sensitive. For samples with MICs over this threshold, individual-level antibiotic resistance profiles were used to determine the mechanism of macrolide resistance.

This inference was possible because more than 97% of pneumococcal resistance to macrolides is mediated by one of two genotypes [22]. Azithromycin resistance arising from the ermB genotype also gives rise to resistance to lincosamides, such as clindamycin [23]. In contrast, mefA/E offers no such cross-protection, allowing for discernment between the two mechanisms by examining whether we also observe clindamycin resistance [24], [25].

It is important to note here that our phenotyping method would classify mefA/E+, ermB+ strains as being ermB+ only. Data from a different area of Ethiopia showed that mefA/E+, ermB+ strains constituted 6.6% of resistant strains after four rounds of azithromycin treatment [20]. It is likely that strains carrying both resistance alleles would have increased antibiotic resistance, but also lowered fitness. Although our model is insensitive to the assumption of antibiotic resistance rates for drug-resistant strains (see Table 1), it is possible that this misclassification could amplify the estimated fitness cost of ermB strains.

Constructing the model

We employed a dynamic SIS (susceptible-infective-susceptible) model including mefA/E, ermB, and antibiotic-sensitive pneumococcal strains, so that within each village, there was a given number of individuals who were uninfected, infected with drug susceptible strains, infected by mefA/E strains, and infected by ermB strains. Between mass treatments, these numbers changed as individuals recovered or became newly infected. During a mass treatment, each individual's probability of cure was the product of the chance of treatment and the chance of cure given the individual's infection status. All equations relevant to this model are included in Supporting Information S1.

The model focused on children between the ages of one and five, as this group has the highest rates of acquisition, and new cases within the age group appear to occur primarily as the result of interaction with fellow members [26], [27]. Based on the distances between villages and the lack of motorized transportation in the area, transmission between villages was assumed to be negligible.

A mixed discrete-continuous time system was used to model this process, with ordinary differential equations (Equation A2) governing the probability that a village has given numbers of individuals of each status between mass treatments. The probability of a given village state following a mass treatment was the sum of the products of two probabilities: the probability of being in a given state, and the probability of transitioning from that state to the state of interest (Equation A2).

We initialized the system of equations using the first set of available village-specific pneumococcal sample results (taken at month 24 after the beginning of the study). Assuming a priori no preference for any given number of individuals infected with each of the three possible strains, we used the observed sampling results to compute a posterior probability distribution for the number of infected individuals in each village (see Equation A1). This posterior distribution was used as the initial condition for data-fitting simulations. Similarly, extrapolation of the probability of extinction versus time was achieved by initializing the maximum likelihood model with our observations at month 54.

In both cases, all strains were assumed to have the same rate of recovery ( An external file that holds a picture, illustration, etc.
Object name is pone.0029407.e007.jpg

An external file that holds a picture, illustration, etc.
Object name is pone.0029407.e007.jpg

), which was extracted from the literature [28]. It should be noted that younger children may exhibit a slower average rate of recovery than our baseline value. As we see in the sensitivity analysis (Table 1), the resulting longer durations of infection imply larger proportional fitness costs, since the reduction in observed prevalence must take place over fewer generations of transmission. The infection rate constant for the antibiotic-sensitive strain, An external file that holds a picture, illustration, etc.
Object name is pone.0029407.e008.jpg

An external file that holds a picture, illustration, etc.
Object name is pone.0029407.e008.jpg

, was calculated from the equation: equilibrium prevalence An external file that holds a picture, illustration, etc.
Object name is pone.0029407.e009.jpg

An external file that holds a picture, illustration, etc.
Object name is pone.0029407.e009.jpg

, where

An external file that holds a picture, illustration, etc.
Object name is pone.0029407.e010.jpg

is the basic reproduction number. Infection rate constants for the antibiotic-resistant strains were allowed to vary during the optimization. In several cases, values examined in the sensitivity analysis improved the likelihood of the model. For the basecase, we chose to keep the parameters that we judged to be the closest match to values from the literature.

For fitness cost estimation, cure rate parameters were also selected for each strain to probabilistically model treatment at months 24 and 30 (Equation A2). These values were chosen based on survey data (B. Ayele. unpublished data), and pneumococcus-specific pharmacodynamics simulations [29]. Comparing these data with observed MIC values allowed us to calculate an expected cure rate for mefA/E of 5%. The cure rate for ermB strains was estimated to be near-zero by this method, but was also set to 5% to avoid introducing an ab initio difference between the two strains. The cure rate was estimated to be 88% for the antibiotic-sensitive strain, by extrapolating estimated treatment efficacies of around 93% from studies that used higher doses of azithromycin, but in populations that presumably harbored resistance [30], [31].

Between treatments, we integrated the ordinary differential equations for the probability of each village state (equation A2) using the GNU Scientific Library [32]. From these probabilities, we used Equation A4 to calculate the likelihood function for observing the village-level strain type data over time for each village (taking the number of children per village to be 50, consistent with survey results). In short, for each village we summed over the probability of each An external file that holds a picture, illustration, etc.
Object name is pone.0029407.e011.jpg

An external file that holds a picture, illustration, etc.
Object name is pone.0029407.e011.jpg

infection state multiplied by the probability of our observations at that time point, given that state. Assuming the villages to be independent, we then multiplied these village-level probabilities together to get an overall probability of the data at that time point given the parameters. These time point probabilities were then multiplied together to get the overall likelihood of the data given the parameters.

We optimized this likelihood with respect to the mefA/E and ermB transmission rate constants, An external file that holds a picture, illustration, etc.
Object name is pone.0029407.e012.jpg

An external file that holds a picture, illustration, etc.
Object name is pone.0029407.e012.jpg

and

An external file that holds a picture, illustration, etc.
Object name is pone.0029407.e013.jpg

, using the optim function in R (http://www.r-project.org). All other parameters were held fixed, and subsequently examined by univariate sensitivity analysis. Asymptotic standard errors for these estimates were found using the observed Fisher information. Finally, the upper bound, lower bound, and point estimates for the relative fitness values of the mefA/E and ermB strains were used to extrapolate each strain's probability of survival versus time.

Supporting Information

Supporting Information S1

This file describes the equations relevant to the construction of our model.

(PDF)

Click here for additional data file.^{(116K, pdf)}

Acknowledgments

We would like to thank Dr. Phillip Spector and the Statistical Computing Facility at the University of California, Berkeley for the use of the center's computing resources.

Footnotes

Competing Interests: The authors have declared that no competing interests exist.

Funding: This work was supported by the International Trachoma Initiative, the Bernard Osher Foundation, That Man May See, the Peierls Foundation, the Bodri Foundation, the Harper Inglis Trust, the South Asia Research Fund, Research to Prevent Blindness, grants U10 EY016214 and R21 AI 55752 from the National Institutes of Health, the Gates Exploration Program (OPP1015714) and the Modeling Infectious Disease Agents Study (NIH GM 1-R01-GM08778). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

References

1. Advisory Committee On Immunization Practices. Preventing pneumococcal disease among infants and young children. Recommendations of the Advisory Committee on Immunization Practices (ACIP). MMWR. Recommendations and Reports. 2000;49(RR-9):1–35. [PubMed]

2. Mera R, Miller LA, Fritsche TR, Jones RN. Serotype replacement and multiple resistance in Streptococcus pneumoniae after the introduction of the conjugate pneumococcal vaccine. Microb Drug Resist. 2008;14(2):101–107. [PubMed]

3. Hanage WP, Fraser C, Tang J, Connor TR, Corander J. Hyper-recombination, diversity, and antibiotic resistance in pneumococcus. Science. 2009;324(5933):1454–1457. [PubMed]

4. Appelbaum PC. Resistance among Streptococcus pneumoniae: Implications for drug selection. Clin Infect Dis. 2002;34(12):1613–1620. [PubMed]

5. Andersson DI, Hughes D. Antibiotic resistance and its cost: is it possible to reverse resistance? Nat Rev Microbiol. 2010;8(4):260–271. [PubMed]

6. Rozen DE, McGee L, Levin BR, Klugman KP. Fitness costs of uoroquinolone resistance in Streptococcus pneumoniae. Antimicrob Agents Chemother. 2007;51(2):412–416. [PMC free article] [PubMed]

7. Marcusson LL, Frimodt-Moller N, Hughes D. Interplay in the selection of uoroquinolone resistance and bacterial fitness. PLoS Pathog. 2009;5(8):e1000541. [PMC free article] [PubMed]

8. Levin BR, Perrot V, Walker N. Compensatory mutations, antibiotic resistance and the popu- lation genetics of adaptive evolution in bacteria. Genetics. 2000;154(3):985–997. [PubMed]

9. Trzcinski K, Thompson CM, Gilbey AM, Dowson CG, Lipsitch M. Incremental increase in fitness cost with increased beta-lactam resistance in pneumococci evaluated by competition in an infant rat nasal colonization model. J Infect Dis. 2006;193:1296–1303. [PubMed]

10. Zhang Q, Sahin O, McDermott PF, Payot S. Fitness of antimicrobial-resistant Campylobacter and Salmonella. Microbes Infect. 2006;8(7):1972–1978. [PubMed]

11. McCallum N, Karauzum H, Getzmann R, Bischoff M, Majcherczyk P, et al. In vivo survival of teicoplanin-resistant Staphylococcus aureus and fitness cost of teicoplanin resistance. Antimicrob Agents Chemother. 2006;50(7):2352–2360. [PMC free article] [PubMed]

12. Dagan R, Barkai G, Givon-Lavi N, Sharf AZ, Vardy D, et al. Seasonality of antibiotic-resistant streptococcus pneumoniae that causes acute otitis media: a clue for an antibiotic-restriction policy? J Infect Dis. 2008;197(8):1094–1102. [PMC free article] [PubMed]

13. Kristinsson KG. Effect of antimicrobial use and other risk factors on antimicrobial resistance in pneumococci. Microb Drug Resist. 1997;3(2):117–123. [PubMed]

14. Riedel S, Beekmann SE, Heilmann KP, Richter SS, Garcia-de Lomas J, et al. Antimicrobial use in Europe and antimicrobial resistance in Streptococcus pneumoniae. Euro J Clin Microb Infect Dis. 2007;26(7):485–490. [PubMed]

15. Dagan R, Klugman K. Impact of conjugate pneumococcal vaccines on antibiotic resistance. Lancet Infect Dis. 2008;8(12):785–795. [PubMed]

16. Luciani F, Sisson SA, Jiang H, Francis AR, Tanaka MM. The epidemiological fitness cost of drug resistance in Mycobacterium tuberculosis. Proc Natl Acad Sci U S A. 2009;106(34):14711–14715. [PubMed]

17. Haug S, Lakew T, Habtemariam G, Alemayehu W, Cevallos V, et al. The Decline of Pneu- mococcal Resistance after Cessation of Mass Antibiotic Distributions for Trachoma. Clin Infect Dis. 2010;51:5–8. [PubMed]

18. Garcia-Rey C, Aguilar L, Baquero F, Casal J, Dal-Re R. Importance of local variations in antibiotic consumption and geographical differences of erythromycin and penicillin resistance in Strep- tococcus pneumoniae. Journal of clinical microbiology. 2002;40(1):159. [PMC free article] [PubMed]

19. Colijn C, Cohen T, Fraser C, Hanage W, Goldstein E, et al. What is the mechanism for persistent coexistence of drug-susceptible and drug-resistant strains of Streptococcus pneumoniae? J R Soc Interface. 2010;7(47):905–919. [PMC free article] [PubMed]

20. Skalet AH, Cevallos V, Ayele B, Gebre T, Zhou Z, et al. Antibiotic selection pressure and macrolide resistance in nasopharyngeal Streptococcus pneumoniae: a cluster-randomized clinical trial. PLoS Med. 2010;7(12):e1000377. [PMC free article] [PubMed]

21. Muluken Melese M, Chidambaram JD, Alemayehu W, Lee DC, Yi EH, et al. Feasibility of eliminating ocular Chlamydia trachomatis with repeat mass antibiotic treatments. JAMA. 2004;292(6):721–725. [PubMed]

22. Isozumi R, Ito Y, Ishida T, Osawa M, Hirai T, et al. Genotypes and related factors reecting macrolide resistance in pneumococcal pneumonia infections in Japan. J Clin Microbiol. 2007;45(5):1440–1446. [PMC free article] [PubMed]

23. Spízek J, Rezanka T. Lincomycin, clindamycin and their applications. Appl Microbiol Biotech- nol. 2004;64(4):455–464. [PubMed]

24. Johnston NJ, De Azavedo JC, Kellner JD, Low DE. Prevalence and characterization of the mechanisms of macrolide, lincosamide, and streptogramin resistance in isolates of Streptococcus pneu- moniae. Antimicrob Agents Chemother. 1998;42(9):2425–2426. [PMC free article] [PubMed]

25. K Waites K, C Johnson C, B Gray B, K Edwards K, M Crain M, et al. Use of clindamycin disks To detect macrolide resistance mediated by ermB and mefE in Streptococcus pneumoniae isolates from adults and children. J Clin Microbiol. 2000;38(5):1731–1734. [PMC free article] [PubMed]

26. Melegaro A, Choi Y, Pebody R, Gay N. Pneumococcal carriage in United Kingdom fami- lies: estimating serotype-specific transmission parameters from longitudinal data. Am J Epidemiol. 2007;166(2):228–235. [PubMed]

27. Hill PC, Townend J, Antonio M, Akisanya B, Ebruke C, et al. Transmission of Streptococcus pneumoniae in rural Gambian villages: a longitudinal study. Clin Infect Dis. 2010;50(11):1468–1476. [PubMed]

28. Melegaro A, Choi Y, George R, Edmunds WJ, Miller E, Gay N. Dynamic models of pneu- mococcal carriage and the impact of the Heptavalent Pneumococcal Conjugate Vaccine on invasive pneumococcal disease. BMC Infect Dis. 2010;10:90. [PMC free article] [PubMed]

29. Sevillano D, Alou L, Aguilar L, Echevarría O, Giménez MJ, et al. Azithromycin iv pharma- codynamic parameters predicting Streptococcus pneumoniae killing in epithelial lining uid versus serum: an in vitro pharmacodynamic simulation. J Antimicrob Chemother. 2006;57(6):1128–1133. [PubMed]

30. Paris R, Confalonieri M, Dal Negro R, Ligia GP, Mos L, et al. Efficacy and safety of azithromycin 1 g once daily for 3 days in the treatment of community-acquired pneumonia: an open-label randomised comparison with amoxicillin-clavulanate 875/125 mg twice daily for 7 days. J Chemother. 2008;20(1):77–86. [PubMed]

31. D'Ignazio J, Camere ME, Lewis DE, Jorgensen D, Breen JD. Novel, single-dose microsphere formulation of azithromycin versus 7-day levooxacin therapy for treatment of mild to moderate community-acquired Pneumonia in adults. Antimicrob Agents Chemother. 2005;49(10):4035–4041. [PMC free article] [PubMed]

32. Rossi F, Galassi M, Davies J, Theiler J, Gough B, et al. GNU Scientific Library Reference Manual. Network Theory Ltd., v1.12 edition; 2009.

Articles from PLoS ONE are provided here courtesy of
Public Library of Science