Abstract
Background
The World Health Organization recommends periodic mass antibiotic distributions to reduce the ocular strains of chlamydia that cause trachoma, the world's leading cause of infectious blindness. Their stated goal is to control infection, not to completely eliminate it. A single mass distribution can dramatically reduce the prevalence of infection. However, if infection is not eliminated in every individual in the community, it may gradually return back into the community, so often repeated treatments are necessary. Since public health groups are reluctant to distribute antibiotics indefinitely, we are still in need of a proven longterm rationale. Here we use mathematical models to demonstrate that repeated antibiotic distributions can eliminate infection in a reasonable time period.
Methods
We fit parameters of a stochastic epidemiological transmission model to data collected before and 6 months after a mass antibiotic distribution in a region of Ethiopia that is one of the most severely affected areas in the world. We validate the model by comparing our predicted results to Ethiopian data which was collected biannually for two years past the initial mass antibiotic distribution. We use the model to simulate the effect of different treatment programs in terms of local elimination of infection.
Results
Simulations show that the average prevalence of infection across all villages progressively decreases after each treatment, as long as the frequency and coverage of antibiotics are high enough. Infection can be eliminated in more villages with each round of treatment. However, in the communities where infection is not eliminated, it returns to the same average level, forming the same stationary distribution. This phenomenon is also seen in subsequent epidemiological data from Ethiopia. Simulations suggest that a biannual treatment plan implemented for 5 years will lead to elimination in 95% of all villages.
Conclusion
Local elimination from a community is theoretically possible, even in the most severely infected communities. However, elimination from larger areas may require repeated biannual treatments and prevention of reintroduction from outside to treated areas.
Background
Trachoma remains the major cause of infectious blindness in the world [1,2]. Repeated infection with the ocular strains of Chlamydia trachomatis can lead to a cascade of conjunctival scarring, inturned eyelids and eyelashes, and eventually blindness due to corneal opacity. To reduce the prevalence of infection, the World Health Organization (WHO) has recommended communitywide distributions of oral azithromycin as part of its strategy to control blinding trachoma by the year 2020[3]. They believe that infection can be reduced to a level low enough that resulting blindness will no longer be a major public health concern, but they do not believe that infection can be completely eliminated from an area.
A single dose of oral azithromycin is clearly effective in eliminating ocular chlamydial infection in an individual [4,5]. Antibiotic distributed simultaneously to an entire community has been shown to reduce the prevalence of infection to low levels [611]. The WHO recommends mass treatments be administered with an observed clinical prevalence greater than 10 percent[3]. Unfortunately, infection returns back into communities over time [7,1215]. Repeated treatments have progressively reduced infection, at least in two areas with moderate disease. After three annual treatments, only a single infection could be identified in a village in western Nepal[8]. Similarly, after a single mass oral azithromycin distribution and three biannual topical tetracycline distributions, only a single infection was found in a Tanzanian village [9]. However, infection has never been eliminated from a hyperendemic area, and if antibiotics are discontinued infection will presumably return. We have no proven longterm rationale for repeated mass antibiotic distributions in severely affected areas. Randomised controlled trials are the gold standard for assessing the effect of different treatment strategies, and several are currently underway. However, it is difficult with trials to explore the long term effects several years in the future for a range of strategies. Here, here we use short term data from Ethiopia to mathematically model longer term outcomes.
In previous studies, we have used deterministic models of difference [16] and differential equations [7] to model infection in communities. Mass antibiotic treatment is incorporated into the model by lowering infection by an amount proportional to the coverage of the community and the efficacy of antibiotic in an individual [7,16]. In the model, infection returns into a community between treatments according to logistic growth at a rate estimated from baseline data[16] or determined from posttreatment Ethiopian data [7]. The models suggest that periodic mass distributions can progressively reduce chlamydial infection in a community, as long as infection is reduced more by each treatment than it returns between treatments (Figure 1a) [7,16]. In this report, we use data collected from 16 villages in Ethiopia to estimate parameters for a mathematical model that incorporates the effects of chance. We then determine whether local elimination is possible, and if so, in what time span.
Figure 1. a. Deterministic Model of Time vs. Prevalence with biannual treatments: Results from a differential equation based model demonstrating that biannual coverage of 80% of the population should progressively reduce ocular chlamydial infection (blue curve). The deterministic model is an excellent approximation for the expectation of the stochastic model (mean of 1000 simulations, red curve). b. Stochastic Model of Time vs. Prevalence with biannual treatments: The mean of 1000 simulations of a stochastic model, assuming biannual treatments with 80% coverage (again, red curve) vs. the average prevalence of only those villages which still harbor infection (green curve). After the third treatment, the average prevalence of infection in these villages returns to approximately the same level with each subsequent treatment (green curve).
Methods
Clinical data
Sixteen villages in the Gurage region of southern Ethiopia were enrolled in a mass antibiotic treatment program for trachoma as previously described [7]. These villages are very remote with approximately 2 kilometres or more between them, thus we consider them to be closed communities. Briefly, 1–5 year old children, those most likely to harbor infection, were monitored at baseline and 2, 6, 12, 18, and 24 months after baseline. Every member of all sixteen villages>= 1 year old in the study were offered antibiotic treatments 1 week after baseline, 6,12, and 18 month visits. Pregnant women and those allergic to macrolides were offered topical tetracycline. The right upper conjunctiva of each child was swabbed and tested for the presence of chlamydial DNA using Amplicor PCR (Roche Diagnostics, Branchburg, NJ). Since posttreatment prevalence was expected to be relatively low, swabs from the same village were randomly pooled into groups of 5 with any remainder pooled into a final tube, as previously described [7,17]. The number of positive individual samples in a village most likely to have resulted in the observed pooled results was determined by maximum likelihood estimation.
Models
Previously, we constructed a simple SIS (susceptible, infected, susceptible) model of ocular chlamydial infection in a core group of children using differential equations [18,19].
where I is number of infectious cases, t is time, γ is the rate of recovery from infection, β is a transmission parameter, and N is the effective population size. Treatment was simulated by reducing the number of individuals infected by the product of the coverage, c, and the efficacy, f, of antibiotic in an individual (Figure 1a). For this study, a stochastic model (continuous time, discrete individuals) analogous to equation 1 having the transitions shown in Table 1 and using the parameters shown in Table 2 and was used.
Note that treatment is not given at a constant rate, but is given periodically to a certain proportion of the population. Each member of the community has a chance of being treated based on the coverage level and the efficacy of antibiotic in an individual.
We constructed the analogous Markov model by letting p_{i}(t) denote the probability that there are i infected individuals in the population at time t (where i varies from 0 to N). Assume that the mass treatments happen at time 0,T,2T,...,kT,...; between periodic mass treatments, the model is a standard continuous time Markov process with constant coefficients. Assuming that infected individuals recover naturally at rate γ, and uninfected individuals become infected at rate βI/N leads us to the Kolmogorov forward equations, which are
between the kth mass treatment and the next (δ_{ij }is the Kronecker delta, which equals 1 if i = j and zero otherwise). At the mass treatment times, however, we assume that each infected individual has a probability c of being treated (the coverage), times a probability f of being cured if treated (the efficacy).
Parameters for this stochastic model were fit to the 6 month Ethiopian data using maximum likelihood estimation. We started out simulations at the average prevalence, and ran for 24 months to allow the distribution of prevalence to approximate the pretreatment distribution at the time point zero We approximated the joint probability distribution function of the 3 points (baseline, 2, and 6 month posttreatment prevalence) by fitting a standard kernel density estimator (a 3dimensional productEpanechnikov kernel), to the results of 1000 simulations. The bandwidth of the kernel in each dimension was set using standard techniques [20]. We used the estimated likelihood to determine the log likelihood of the observed data. The values of the parameters β and γ which maximized the estimated likelihood were determined iteratively and uncertainty was estimated by inverting the Hessian of the loglikelihood function. Sensitivity analyses were performed by varying a single parameter at a time, keeping other parameters set at β = 0.044 and γ = 0.017, N = 100, effective coverage of 90%, and treatment frequency biannual.
For the sensitivity analysis, we varied γ from 0.005 to 0.022, showing (Figure 7). We varied β from 0.01 to 0.1, illustrated in Figure 6, with the corresponding baseline prevalence estimated from the deterministic model
(
Figure 6. Baseline prevalence vs. Years until elimination: Here we vary β in the stochastic model, while keeping other parameters the same. Other parameters are biannual treatment with 90% effective coverage, an effective population size of 100 children, and γ = 0.017.
Figure 7. Recovery rate vs. Years until elimination: Here we vary the recovery rate in the stochastic model, while keeping other parameters the same. Other parameters are biannual treatment with 90% effective coverage, an effective population size of 100 children, and β = 0.044.
To incorporate the effects of treatment, we used the analogous discrete time, Markov process where time is incremented every period of mass antibiotic distribution (e.g. 26 weeks). The quasistationary distribution (Figure 3a, violet curve) was determined by constructing a vector q, with q_{i }being the probability of having i infectious cases in the community (from i = 1 to the number of individuals in the community, n) and an n × n matrix M (determined from simulations), such that,
Figure 3. a. Probability Density of infection prevalence found in biannually treated villages where "village level" elimination has not yet occurred: Probability distribution of 1000 simulations at baseline, 6, 12, 18, and 24 months. The prevalence of infection in a simulated community pretreatment varies in a normal distribution [21]. Each mass treatment eliminates infection in some villages, but in those that it does not, the distribution is shifted to the left, rapidly approaching a quasistationary distribution 3b. Probability density graphs using Ethiopian data: Distribution of the prevalence of infection in preschool children in 16 Ethiopian villages. Baseline, 2, and 6 month data were used to fit the parameters of the stochastic model. Subsequent data from 12, 18, and 24 months confirm that the distribution of infection in the villages also approaches the quasistationary distribution.
M·q^{t }= q^{t+treatment period}
Note that q is similar to p above, except that q_{0, }the state that would represent elimination, is specifically not included in this formulation, so only communities where infection has not been eliminated are followed. The largest eigenvalue of M is the reciprocal of the expected time to extinction, and the corresponding eigenvector provides the quasistationary distribution (violet curve, Figure 3a) [21]. All models were built using Mathematica Version 5.1 using the statistical packages Discrete Distributions and Continuous Distributions.
Results and Discussion
Parameters for a stochastic model were estimated from longitudinal data from approximately 5000 people within 16 separate Ethiopian villages during the first 6 months after a mass distribution. The MLE estimates for our model gave β = 0.044 and γ = 0.017. Gamma can be interpreted as the reciprocal of the duration of infection, and here would be approximately 1 year.
With estimated parameters, our mathematical models did indeed show that elimination is possible with repeated treatments. The stochastic element of the model recreates what we see happening in real life. In some villages we see infection returning, while in others it disappears. As more treatments are administered, we see more and more villages achieving elimination. Observing results from a 1000 simulations using the stochastic model, the average prevalence of infection is progressively reduced with each periodic treatment as in the deterministic model (Figure 1a and 1b, red curve). However, this reduction does not necessarily occur in all of the villages. In some simulated villages infection is eliminated with a single mass treatment. In others, it is reduced to a low level and fades out by chance (Figure 2a, red curve). In still others, infection may return quite rapidly (Figure 2b, blue curve). Simulations of 1000 communities demonstrate that infection is eliminated in more villages with each subsequent treatment. If not eliminated, it returns on average to the same level before the next scheduled treatment (green curve, Figure 1b). In fact, if we set aside the villages that have achieved elimination, the distribution of the prevalence of infection after the 3^{rd }treatment is similar to that after the 4^{th }treatment, and prevalence after the 4^{th }even more similar to that after the 5^{th }(Figure 3a). After two treatments, simulated communities quickly approach a distribution which, if conditioned on nonelimination, is stationary (termed a quasistationary distribution) [2123]. This quasistationary distribution can be determined analytically (violet curve in Figure 3a, see methods).
Figure 2. a. Simulation data after one treatment: In identical simulated communities, infection responds to identical treatment in different ways. It may return after a single mass antibiotic treatment relatively rapidly (blue curve) or fade out (red curve) due to the effects of chance. b. Real Ethiopian villages after one treatment: In Ethiopian communities with similar pretreatment prevalence of infection and similar antibiotic coverage levels, infection may return relatively rapidly (blue curve), or fade out (red curve).
Subsequent data from the same sixteen Ethiopian villages at 12, 18, and 24 months appear to confirm that elimination is possible. In one village, 50% of the children were infected at baseline, but no infection was identified in any of the 5 visits after the first treatment. In another, 45% of children were infected at baseline, but no infection could be found after the second treatment. If infection was not eliminated in a village then the distribution to which it returned 6 months after the 2^{nd }treatment approaches that to which it returns after the 3rd treatment (Figure 3b). There is no statistical difference between distributions 6 months after the 3^{rd }and 4^{th }treatments (KolmogorovSmirnov test, P = 0.95 comparing 18 and 24 month data).
This stochastic model does not incorporate reintroduction of infection from neighbouring villages. The rate that infection is reintroduced is not known and is difficult to estimate. Reintroduction appears not to be frequent, as communities do well in the short term, whether or not neighbouring villages were treated (unpublished data). On the other hand, it may well be occurring [13]. Some villages in which no infection could be found in children for several visits had infection identified at a subsequent visit. Even occasional transmission from other communities would prevent sustained elimination. Interestingly, the quasistationary distribution can be estimated from models by only considering villages where infection has occurred (the limiting case of the curves in Figure 3), by preventing the last infectious case from disappearing (violet curve in Figure 3a), or by artificially introducing a single infection into a community whenever infection has been eliminated, much as might be expected with reintroduction[21]. Since each of these methods mimics reintroduction, then the quasistationary distribution may be a reasonable estimate of the distribution of prevalence of infection.
Investigators are struggling to determine why infection returns in some villages and not others. This search is important, but may not be fruitful, since a great deal of variation is expected even in otherwise similar villages. Simulations allow us to construct communities that are absolutely identical, and variation is still observed due to the vagaries of who infects whom. Infection will be eliminated in some fortunate villages. In other, essentially identical villages, it will continue to return to the same average level even after multiple treatments (Figure 2a). In practice, we have also observed a varied response to treatment in similar villages (Figure 2b) [12].
Often overall prevalence is reported, but we believe there needs to be a shift in the way we evaluate progress. The average prevalence of infection across several villages can be misleading. A regional prevalence of 5% does not mean that 1 in 20 children in each village are infected. More likely, infection has been eliminated in most villages but has returned in a few villages to a level that may far exceed 5%. A better measure of success of a trachoma program may be the proportion of villages in which infection has been eliminated. While the power of sampling allows the mean prevalence to be obtained relatively easily, it should be interpreted in terms of what it reveals about the underlying process; after several rounds of treatment, the mean represents some villages where infection has been completely eliminated and others where infection has returned to a level chosen from the quasistationary distribution. Thus the decision for stopping distributions may need to depend on whether infection has been eliminated locally.
Periodic mass antibiotic treatments can locally eliminate the ocular chlamydia that cause blinding trachoma, if given frequently and to a large portion of the population [7,16]. Models suggest it will return to the same distribution before the next scheduled treatment. If infection is not eliminated from a community after several mass treatments, then stochastic Ethiopian data demonstrates that this quasistationary distribution is closely approached after only two or three treatments. The WHO's strategy requires only for antibiotics to reduce the level of infection and expects that other measures will be necessary for a permanent solution. Their current recommendation of three treatments before reevaluation is reasonable. If constant reinfection of communities prevents elimination, then subsequent treatments may only serve to maintain the prevalence in a stationary distribution, not lower it. However, if infection between neighbouring communities is rare or can be reduced by treating large areas within a narrow time frame, then elimination over larger areas than a village may be achieved. There is some concern that communities from which trachoma has been partially eliminated will lose much of their immunity to chlamydia, only to have it return with a vengeance after treatments have been discontinued [2426]. This has yet to be demonstrated convincingly in practice [26], but potential loss of immunity could easily be included in future models.
Partners in the WHO's trachoma program have distributed over 30 million doses of oral azithromycin so far, and over 1 million in Ethiopia. These programs are expanding rapidly. There has been evidence of subsequent reemergence of infection in some recent trials of azithromycin, suggesting coverage and dosing intervals in mass therapy need careful consideration [7,14,24,27]. Biannual treatments are costly, require more resources, and may contribute to antibiotic resistance in chlamydia and other pathogens such as Streptococcus pneumoniae [26,2830]. The benefits of mass therapy thus need to be weighed against the potential emergence of antibiotic resistance which will require long term surveillance. But if necessary for local elimination biannual distributions may be more cost effective in the long run in severely affected areas. The elimination of a bacterial disease from a large area with mass antibiotics would be major medical breakthrough.
Simulations were run to estimate the number of rounds of treatment necessary to achieve elimination in 50%, 75%, and 95% of the villages. Elimination depends on the efficacy and coverage of treatment, as well as the effective population size in a community. It should be noted that these models were developed using data from a hyperendemic region, and they are not generalizable to all areas.
The stochastic model suggests that treatment given to 90% of the population biannually would eliminate infection in 95% of the villages after 5 years. Lower coverage or less frequent treatments would not be as successful (Figures 4 and 5). For example, the WHO's recommended strategy of annual treatment of 80% of the population would eliminate infection in an estimated 95% of hyperendemic villages in 12 years. Not surprisingly, areas with a lower baseline prevalence of infection should achieve elimination more rapidly (Figure 6). If the rate that an individual recovers from an infection is more rapid than estimated here, then it will be more difficult to eliminate infection with repeated treatments (Figure 7). Infection can be eliminated from smaller communities more rapidly than larger ones, in part because stochastic fade out is more likely with fewer cases (Figure 8). Other factors, such as the cost of medication and of distributions, were not analyzed in this report, but clearly play a role in determining the optimal strategy for a given program.
Figure 4. Coverage vs. Years until elimination: Here we vary coverage in the stochastic model, while keeping other parameters the same. Other parameters are biannual treatment an effective population size of 100 children, γ = 0.017, and β = 0.044.
Figure 5. Number of months between treatments vs. Years until elimination: Here we vary treatment frequency in the stochastic model, while keeping other parameters the same. Other parameters are 90% effective coverage, an effective population size of 100 children, γ = 0.017, and β = 0.044.
Figure 8. Population vs. Years until elimination: Here we vary the effective population size of children in the stochastic model, while keeping other parameters the same. Other parameters are biannual treatment with 90% effective coverage, γ = 0.017, and β = 0.044.
Conclusion
The local elimination of infectious trachoma is possible with repeated mass antibiotic distributions. A stochastic mathematical model suggests that elimination (reduction of the prevalence of infection to zero in a community) can occur within a reasonable time period even in the most severely affected areas, as long as the coverage and frequency of distributions are high enough. For example, in hyperendemic areas of Ethiopia, we estimate that biannual treatment with 90% effective coverage would result in elimination in more than half of communities in 3 years, and in 95% within 5 years. Most areas have far less trachoma than the 16 villages in Ethiopia presented here, so elimination should prove far easier. The WHO's current recommendation of three annual treatments with at least 80% coverage will likely eliminate infection in a large number of communities worldwide. The time to elimination is also dependent on whether transmission of infection can be effectively reduced by any other means, such as hygiene (clean faces are associated with less clinically active trachoma) and fly control (Musca sorbens is thought to play a role in the transmission of ocular chlamydia), both of which are major parts of the WHO's overall trachoma strategy [3134]. Researchers have noted that trachoma is disappearing in the absence of active trachoma programs [3538]; such a secular trend would also make infection easier to eliminate.
Competing interests
The author(s) declare that they have no competing interests.
Authors' contributions
KR acquisition of data, mathematical modelling, analysis and interpretation of data, drafting of the manuscript, critical revision of the manuscript for important intellectual content, statistical analysis.
TP study concept and design, analysis and interpretation of data, drafting of the manuscript, critical revision of the manuscript for important intellectual content, statistical analysis.
KH acquisition of data, analysis and interpretation of data, critical revision of the manuscript for important intellectual content, administrative, technical, or material support, study supervision.
DL mathematical modeling, analysis and interpretation of data, critical revision of the manuscript for important intellectual content, statistical analysis.
WA study concept and design, acquisition of data, critical revision of the manuscript for important intellectual content, obtained funding.
MM study concept and design, acquisition of data, critical revision of the manuscript for important intellectual content, obtained funding.
TL study concept and design, acquisition of data, critical revision of the manuscript for important intellectual content, obtained funding.
EY acquisition of data, critical revision of the manuscript for important intellectual content, administrative, technical, or material support, study supervision.
JH acquisition of data, critical revision of the manuscript for important intellectual content, administrative, technical, or material support, study supervision.
JC acquisition of data, analysis and interpretation of data, critical revision of the manuscript for important intellectual content.
JW study concept and design, acquisition of data, drafting of the manuscript, critical revision of the manuscript for important intellectual content, obtained funding.
BG study concept and design, acquisition of data, critical revision of the manuscript for important intellectual content.
TL conceived of the study concept and design, acquisition of data, mathematical modelling, drafting of the manuscript, critical revision of the manuscript for important intellectual content, statistical analysis, obtained funding, study supervision.
Acknowledgements
We would like to thank the Bernard Osher Foundation, the International Trachoma Initiative, the Bodri Foundation, the South Asia Research Fund, the Harper Inglis Trust, Research to Prevent Blindness, That Man May See, and the National Institutes of Health (Grant U10EY016214) for financial support, and Dr. Arek Goetz and Dr. Mohammad Kafai for guidance.
References

Thylefors B, Négrel AD, Pararajasegaram R, Dadzie KY: Global data on blindness.
Bull World Health Organ 1995, 73(1):115121. PubMed Abstract

Resnikoff S, Pascolini D, Etya'ale D, Kocur I, Pararajasegaram R, Pokharel GP, Mariotti SP: Global data on visual impairment in the year 2002.
Bull World Health Organ 2004, 82(11):844851. PubMed Abstract  Publisher Full Text

Report of the Eighth meeting of the W.H.O. Alliance for the Global Elimination of Blinding Trachoma. Geneva
2004.

Bailey RL, Arullendran P, Whittle HC, Mabey DC: Randomised controlled trial of singledose azithromycin in treatment of trachoma.
Lancet 1993, 342(8869):453456. PubMed Abstract  Publisher Full Text

Dawson CR, Schachter J, Sallam S, Sheta A, Rubinstein RA, Washton H: A comparison of oral azithromycin with topical oxytetracycline/polymyxin for the treatment of trachoma in children.
Clin Infect Dis 1997, 24(3):363368. PubMed Abstract

Schachter J, West SK, Mabey D, Dawson CR, Bobo L, Bailey R, Vitale S, Quinn TC, Sheta A, Sallam S, Mkocha H, Faal H: Azithromycin in control of trachoma.
Lancet 1999, 354(9179):630635. PubMed Abstract  Publisher Full Text

Melese M, Chidambaram JD, Alemayehu W, Lee DC, Yi EH, Cevallos V, Zhou Z, Donnellan C, Saidel M, Whitcher JP, Gaynor BD, Lietman TM: Feasibility of eliminating ocular Chlamydia trachomatis with repeat mass antibiotic treatments.
JAMA 2004, 292(6):721725. PubMed Abstract  Publisher Full Text

Gaynor BD, Miao Y, Cevallos V, Jha H, Chaudary JS, Bhatta R, OsakiHolm S, Yi E, Schachter J, Whitcher JP, Lietman T: Eliminating trachoma in areas with limited disease.
Emerg Infect Dis 2003, 9(5):596598. PubMed Abstract  Publisher Full Text

Solomon AW, Holland MJ, Alexander ND, Massae PA, Aguirre A, NatividadSancho A, Molina S, Safari S, Shao JF, Courtright P, Peeling RW, West SK, Bailey RL, Foster A, Mabey DC: Mass treatment with singledose azithromycin for trachoma.
N Engl J Med 2004, 351(19):19621971. PubMed Abstract  Publisher Full Text

Chidambaram JD, Melese M, Alemayehu W, Yi E, Prabriputaloong T, Lee DC, Cevallos V, Zhou Z, Whitcher JP, Gaynor BD, Lietman TM: Mass antibiotic treatment and community protection in trachoma control programs.
Clin Infect Dis 2004, 39(9):e957. PubMed Abstract  Publisher Full Text

West ES, Munoz B, Mkocha H, Holland MJ, Aguirre A, Solomon AW, Bailey R, Foster A, Mabey D, West SK: Mass treatment and the effect on the load of Chlamydia trachomatis infection in a trachomahyperendemic community.
Invest Ophthalmol Vis Sci 2005, 46(1):8387. PubMed Abstract  Publisher Full Text

Chidambaram JD, Lee DC, Porco TC, Lietman TM: Mass antibiotics for trachoma and the Allee effect.
Lancet Infect Dis 2005, 5(4):194196. PubMed Abstract  Publisher Full Text

Burton MJ, Holland MJ, Makalo P, Aryee EA, Alexander ND, Sillah A, Faal H, West SK, Foster A, Johnson GJ, Mabey DC, Bailey RL: Reemergence of Chlamydia trachomatis infection after mass antibiotic treatment of a trachomaendemic Gambian community: a longitudinal study.
Lancet 2005, 365(9467):13211328. PubMed Abstract  Publisher Full Text

West SK, Munoz B, Mkocha H, Holland MJ, Aguirre A, Solomon AW, Foster A, Bailey RL, Mabey DC: Infection with Chlamydia trachomatis after mass treatment of a trachoma hyperendemic community in Tanzania: a longitudinal study.
Lancet 2005, 366(9493):12961300. PubMed Abstract  Publisher Full Text

Chidambaram JD, Alemayehu W, Melese M, Lakew T, Yi E, Cevallos V, Zhou Z, Maxey K, Shapiro B, Srinivasan M, Whitcher JP, Gaynor BD, Lietman TM: Return of Infectious Trachoma after a Single Mass Antibiotic Distribution.
JAMA 2006., 295(10) PubMed Abstract  Publisher Full Text

Lietman T, Porco T, Dawson C, Blower S: Global elimination of trachoma: how frequently should we administer mass chemotherapy?
Nat Med 1999, 5(5):572576. PubMed Abstract  Publisher Full Text

Diamont J, Moncada J, Pang F, Jha H, Benes R, Schachter J, Lietman. T: Pooling of Chlamydial Laboratory Tests to Determine Prevalence of Trachoma.
Ophthalmic Epidemiology 2001, 8(23):109117. PubMed Abstract  Publisher Full Text

Hethcote HW, Yorke JA: Gonorrhea transmission dynamics and control. In Series Lecture notes in biomathematics. Volume 56. Berlin , SpringerVerlag; 1984.

Lee DC, Chidambaram JD, Porco TC, Lietman TM: Seasonal effects in the elimination of trachoma.
Am J Trop Med Hyg 2005, 72(4):468470. PubMed Abstract  Publisher Full Text

Scott DW, Wand MP: Feasibility of Multivariate Density Estimates.
Biometrika 1991, 78(1):197205. Publisher Full Text

Nasell I: On the quasistationary distribution of the stochastic logistic epidemic.
Mathematical Biosciences 1999, 156(12):2140. PubMed Abstract  Publisher Full Text

Jacquez JA, Simon CP: The stochastic SI model with recruitment and deaths. I. Comparison with the closed SIS model.
Math Biosci 1993, 117(12):77125. PubMed Abstract  Publisher Full Text

Nasell I: On the quasistationary distribution of the Ross malaria model.
Mathematical Biosciences 1991, 107(2):187207. PubMed Abstract  Publisher Full Text

Atik B, Thanh TT, Luong VQ, Lagree S, Dean D: Impact of annual targeted treatment on infectious trachoma and susceptibility to reinfection.
Jama 2006, 296(12):14881497. PubMed Abstract  Publisher Full Text

Brunham RC, Pourbohloul B, Mak S, White R, Rekart ML: The unexpected impact of a Chlamydia trachomatis infection control program on susceptibility to reinfection.
J Infect Dis 2005, 192(10):18361844. PubMed Abstract  Publisher Full Text

Yang JL Lietman, T.L.: The return of trachoma in partially treated communities.
Archives of Ophthalmology 2007, 125(7):989991. PubMed Abstract  Publisher Full Text

Chidambaram JD, Alemayehu W, Melese M, Lakew T, Yi E, House J, Cevallos V, Zhou Z, Maxey K, Lee DC, Shapiro BL, Srinivasan M, Porco T, Whitcher JP, Gaynor BD, Lietman TM: Effect of a single mass antibiotic distribution on the prevalence of infectious trachoma.
Jama 2006, 295(10):11421146. PubMed Abstract  Publisher Full Text

Leach AJ, ShelbyJames TM, Mayo M, Gratten M, Laming AC, Currie BJ, Mathews JD: A prospective study of the impact of communitybased azithromycin treatment of trachoma on carriage and resistance of Streptococcus pneumoniae.
Clin Infect Dis 1997, 24(3):356362. PubMed Abstract

Solomon AW, Mohammed Z, Massae PA, Shao JF, Foster A, Mabey DC, Peeling RW: Impact of mass distribution of azithromycin on the antibiotic susceptibilities of ocular Chlamydia trachomatis.
Antimicrob Agents Chemother 2005, 49(11):48044806. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Fry AM, Jha HC, Lietman TM, Chaudhary JS, Bhatta RC, Elliott J, Hyde T, Schuchat A, Gaynor B, Dowell SF: Adverse and beneficial secondary effects of mass treatment with azithromycin to eliminate blindness due to trachoma in Nepal.
Clin Infect Dis 2002, 35(4):395402. PubMed Abstract  Publisher Full Text

West S, Munoz B, Lynch M, Kayongoya A, Chilangwa Z, Mmbaga BB, Taylor HR: Impact of facewashing on trachoma in Kongwa, Tanzania.
Lancet 1995, 345(8943):155158. PubMed Abstract  Publisher Full Text

Emerson PM, Cairncross S, Bailey RL, Mabey DC: Review of the evidence base for the 'F' and 'E' components of the SAFE strategy for trachoma control.
Trop Med Int Health 2000, 5(8):515527. PubMed Abstract  Publisher Full Text

Emerson PM, Lindsay SW, Alexander N, Bah M, Dibba SM, Faal HB, Lowe KO, McAdam KP, Ratcliffe AA, Walraven GE, Bailey RL: Role of flies and provision of latrines in trachoma control: clusterrandomised controlled trial.
Lancet 2004, 363(9415):10931098. PubMed Abstract  Publisher Full Text

Miller K, Pakpour N, Yi E, Melese M, Alemayehu W, Bird M, Schmidt G, Cevallos V, Olinger L, Chidambaram J, Gaynor B, Whitcher J, Lietman T: Pesky trachoma suspect finally caught.
The British journal of ophthalmology 2004, 88(6):750751. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Chidambaram JD, Bird M, Schiedler V, Fry AM, Porco T, Bhatta RC, Jha H, Chaudary JSP, Gaynor B, Yi E, Whitcher JP, OsakiHolm S, Lietman TM: Trachoma Decline and Widespread Use of Antimicrobial Drugs.
Emerg Infect Dis 2004, 10(11):18951899. PubMed Abstract  Publisher Full Text

Dolin PJ, Faal H, Johnson GJ, Minassian D, Sowa S, Day S, Ajewole J, Mohamed AA, Foster A: Reduction of trachoma in a subSaharan village in absence of a disease control programme.
Lancet 1997, 349(9064):15111512. PubMed Abstract  Publisher Full Text

Hoechsmann A, Metcalfe N, Kanjaloti S, Godia H, Mtambo O, Chiopeta T, Barrows J, Witte C, Courtright P: Reduction of trachoma in the absence of antibiotic treatment: Evidence from a populationbased survey in Malawi.
Ophthalmic Epidemiology 2001, 8(23):145153. PubMed Abstract  Publisher Full Text

Jha H, Chaudary JSP, Bhatta R, Miao Y, OsakiHolm S, Gaynor B, Zegans M, Bird M, Yi E, Holbrook K, Whitcher JP, Lietman. T: Disappearance of trachoma in western Nepal.
Clinical Infectious Diseases 2002, 35(6):765768. PubMed Abstract  Publisher Full Text
Prepublication history
The prepublication history for this paper can be accessed here: