Volume 20, Number 6—June 2014
Short-Term Malaria Reduction by Single-Dose Azithromycin during Mass Drug Administration for Trachoma, Tanzania
Single-dose mass drug administration of azithromycin (AZT) is underway to eliminate trachoma worldwide. Studies in Ethiopia showed a reduction in all-cause childhood deaths after administration. To examine the effect of single-dose AZ MDA on prevalent malaria infections in a large prospective cohort of children and parents in Dodoma Province, Tanzania, we quantified the temporal prevalence of malaria parasitemia by real-time PCR for 6 months after single-dose AZT. In the first month after treatment but not in subsequent months, Plasmodium falciparum infections were reduced by 73% (95% CI 43%–89%) in treatment versus control villages and differences remained significant (p = 0.00497) in multivariate models with village-level random effects. Genetic sequencing of P. falciparum ribosomal L4 protein showed no mutations associated with AZT resistance. AZT mass drug administration caused a transient, 1-month antimalarial effect without selecting for P. falciparum ribosomal L4 resistance mutations in a region with a 10-year history of treating trachoma with this drug.
Malaria can be treated or prevented with the broad-spectrum antimicrobial drugs tetracycline or azithromycin (AZT) (1). In vitro, AZT interferes with malarial parasite replication by targeting the unique apicoplast organelle of the parasite (1). AZT inhibits malarial parasite growth 10–fold every 48 hours, and the pharmacokinetics of AZT predict that it remains at concentrations high enough to limit parasite growth for >1 week (2). AZT might interfere with transmission by exoerythrocytic inhibition of parasite liver stages in humans and mice (3) and by interference with ookinete and sporozoite production in mosquitoes (4). Monotherapy with AZT is not typically used to treat malaria. However AZT is highly effective against Chlamydia trachomatis, the causative agent of trachoma, which causes blindness (5,6). Persons with malaria who live in trachoma-endemic regions may undergo repeated AZT therapy as part of the World Health Organization–sponsored global trachoma eradication program (7,8).
Data regarding the effects of AZ mass drug administration (MDA) on malaria are limited. In a cluster randomized trial in The Gambia, AZ MDA given in 3 doses (20 mg/kg) 7 days apart reduced malaria rates by half when measured at 1 time point in children 5–14 years of age (9). More recently, a 49% reduction in the odds of death (95% CI 29%–90%) was reported for children 1–9 years of age in an AZ MDA treatment group compared with controls in Ethiopia (10,11). Porco et al. suggested that reductions in malaria prevalence associated with AZT MDA might have contributed to observed decreases in overall deaths. (10).
For malaria treatment, a randomized clinical trial that compared AZT/artesunate with artemether/lumefantrine in Muheza, Tanzania, reported that the odds of treatment failure were 5 times greater (95% CI 3.3–11.4) in the group that received AZT/artesunate (12). The authors postulated that the AZT/artesunate showed treatment failure because MDA for trachoma in Tanzania could have led to localized Plasmodium spp. resistance to AZT (12). Plasmodium spp. drug resistance to AZT has not been documented in the field. However, in vitro selection for AZT resistance identified a G76V mutation among conserved active site amino acids at position 71–79 in P. falciparum apicoplast-encoded ribosomal protein L4 (PfRpL4) (PFC10_API0043) (1.) A Cochrane meta-analysis report stated that “azithromycin’s future for the treatment of malaria does not look promising” (13) and cited studies in which AZT, although well tolerated, was inferior to tetracycline for malaria prophylaxis (14–17).
In the current study, we evaluated malaria prevalence in a cohort of 2,053 children and adults in central Tanzania to examine the effect of single dose AZ MDA on prevalent malaria infections. We also searched for PfRpL4 mutations that might confer P. falciparum AZT resistance.
Study Site and Sampling
Study participants were from 8 rural agricultural villages in Dodoma Province, central Tanzania (18–20). The study period was January 12, 2009–July 21, 2009 and was coincidental with a period in which rainfall was <60% of the average amount (21,22). AZT was offered to all the residents of the treatment villages, including residents not part of the follow-up investigation. Four treatment villages were selected on the basis of trachoma prevalence >10% in children 1–9 years of age (trachoma prevalence 12%, 18%, 18%, and 14%), rather than by randomization (18,19), which was consistent with World Health Organization AZT MDA guidelines for trachoma control (4). Four control villages with lower trachoma rates (8%–10%) were chosen on the basis of geographic proximity.
A complete census was conducted in the villages> In each village, 130 families with children <5 years of age were randomly selected for follow-up prevalence blood sampling; 1 child and 1 adult were randomly selected from each family. Sample size estimations were based on prestudy estimates of a malaria prevalence of 10% in children <5 years of age in the Dodoma region (23). Follow-up fingerprick blood sampling was performed at baseline and at 1, 3, 4, and 6 months later. Children were divided into 3 subgroups that were sampled once at weeks 2 (group 1), 6 (group 2), and 8 (group 3). Contemporaneous to the monthly blood sampling, staff members from villages conducted an active surveillance program in which children were visited weekly and screened for axillary temperature and questionnaires were answered regarding fever, diarrhea, and respiratory disease (18,19).
Blood samples were collected on ProteinSaver903 (GE Healthcare, Pittsburgh, PA, USA) for analysis by real-time PCR for all participants at each scheduled monthly visit (20). When a fever (temperature ≥37.5°C) in adults and children was observed or when a caregiver reported a history of fevers in children during the weekly surveillance, a rapid diagnostic test (RDT) (Paracheck-Pf; Orchid Biomedical Systems, Goa, India) and thick and thin blood films were prepared. RDTs were performed by trained study staff and interpreted in the field at the time of sampling. Two experienced microscopists at the Amani Laboratory in Muheza, Tanzania, read the slides blinded to the PCR or RDT results. Discordant results were read by a third microscopist. Persons with positive results by RDT and all febrile children <5 years of age were treated with artemether/lumefantrine according to national guidelines and removed from later analysis.
Information on bednet ownership was obtained through standardized participant interviews. Latitude, longitude, and altitude of home locations were measured by using the GPSMAP 76 unit (Garmin, Olathe, KS, USA).
DNA extraction and real-time PCR are described elsewhere (20). Five 3-mm diameter punches, equivalent to 25 μL of whole blood, were removed from the filter papers. DNA was extracted by using a commercial 96-well kit (Promega, Fitchburg, WI, USA). The DNA was concentrated by glycogen acetate and acetate/ethanol precipitation and low-speed (3,000 × g) centrifugation for 30 min.
Multiplex real-time PCR was used to amplify the 18S P. falciparum ribosomal gene with a Cy5-labeled probe (20). Samples were processed in duplicate on manually loaded 384-well plates. A 40-cycle standard PCR protocol was used in a CFX 384 real-time PCR Detection System Thermocycler (BioRad, Hercules, CA, USA). Baseline relative fluorescence units (RFUs) were readjusted by using Bio-Rad CFX manager software. The real-time PCR system also detected Borrelia spp. reported by Reller et al. (24). The sensitivity and specificity of the real-time PCR versus that of RDT or microscopy were reported by Schachterele et al. (20). The real-time PCR detected 1–100 parasites/μL and showed the highest sensitivity in latent class analysis on febrile study patients. For the monthly surveys, microscopic analysis was not performed for the 8,711 samples for which real-time PCR was conducted. The RFU cutoff used was greater than that for >50 negative blood samples from control patients at Johns Hopkins University. The cutoff of 650 RFUs for the last cycle of the highest real-time PCR replicate had a specificity of 100% for samples from the control group at Johns Hopkins University and a specificity of 94% for samples from febrile patients in Tanzania; microscopy was used as a reference method (20).
Single-Nucleotide Polymorphism Analysis of PfRpL4 Gene for AZT Resistance
Blood samples from treatment and control villages that had higher parasite densities by real-time PCR and microscopy were subjected to PCR amplification and whole gene fragment cloning into the pCR2.1 plasmid. Amplification and cloning were followed by full PfRpL4 gene DNA sequencing of multiple bacterial clones for each patient and examination of mutations anywhere in this gene (1).
Data were analyzed from the perspective of an intention to treat in which participants from AZ MDA–treated villages who refused treatment were classified with those who accepted AZT as members of treatment villages. Univariate ratios were used to compare prevalence proportions in treated versus untreated villages. CIs for ratios that compared proportions of prevalent infection in treated versus untreated villages were estimated by using 2-sided exact binomial tests.
Maps were used to display geospatial patterns in infection prevalence and malaria clustering over time (Google Maps, Mountain View, CA, USA, and Quantum Geographic Information System Open Source Geospatial Foundation Project (www.qgis.org/en/site/). Geographic coordinates of P. falciparum–positive and –negative blood samples were projected into kilometers, and locations were smoothed by using quadratic kernel intensity estimation. Kernel intensity estimators predict spatial intensity at unsampled points from nearby sampled points by using a quadratic function to heavily weight the nearest measured points and de-emphasize points sampled at greater distances. Weighted integrals were used to estimate intensity smoothing parameters for quadratic functions, assuming a 2-dimensional (i.e., spatial), stationary, isotropic point pattern process. Spatial odds ratios (ratio of kernel intensities from positive blood samples to kernel intensities from negative blood samples) (25), were determined for P. falciparum.
To examine the effect of AZT MDA while adjusting for malaria clusters observed on maps, we used multivariate logistic models with random effects to compare odds of prevalent malaria infection for AZT MDA treatment villages with odds for control villages. Random intercepts were fit to account for the multilevel or hierarchical design that nested persons within villages. This model enabled accurate statistical inference by adjusting CIs and p values for P. falciparum clusters in villages, and this inference was apparent in prevalence maps (25). Residual spatial autocorrelation was assessed by using variograms, and residual temporal autocorrelation was eliminated by restricting models to 1 sampling interval (26). Multiple models that examined the effect of AZT MDA on prevalent malaria were fit to control for confounding by other drivers of malaria prevalence and to address the robustness of standard errors to potential spatial residual autocorrelation. Altitude and self-reported bednet ownership were held constant to control for confounding in the multivariate model. A priori, we believed any AZT MDA effect would be most perceptible between month 1 and baseline.
To distinguish between the effect of AZT on prevalence and incident infections, we performed a sensitivity analysis for only P. falciparum infections that had been preceded by a negative P. falciparum test result in the previous sampling interval. These data also underwent a second sensitivity analysis that included participants treated with artemether/lumefantrine who had been removed from the previous analyses.
Data management and analysis were conducted by using R statistical software (27). R with the SPLANCS package was used for kernel intensity estimation (28), and the LME4 package was used for the random effects model (29). All R code is available upon request.
Characteristics of Study Population
The study site in the Kongwa District of Tanzania included 8 villages and 12,898 persons. Four villages were chosen on the basis of increased prevalence of trachoma (>10%) in children 1–9 years of age. Excluded were 34 of 66 villages that had received AZT MDA in the previous year. In the 4 villages that received AZT MDA, 6,252 persons received AZT and 642 did not receive AZT. A census in untreated control villages identified 5,991 persons. These villages were chosen because of geographic proximity to treatment villages. In each village, pairs of a parent and a child <5 years of age were randomly chosen for inclusion (1,045 persons in the AZT MDA treatment group and 1,008 persons in the control group for follow-up) (Figure 1). Blood samples for real-time PCR testing were collected from participants at baseline and at months 1, 3, 4, and 6 (4,437 patient-time samples from AZT MDA treatment villages and 4,274 patient-time samples from control villages). At any follow-up period after baseline, 1,010 (97%) study participants in AZT MDA villages who received AZT and 976 (97%) study controls were sampled for real-time PCR at least once. At baseline, treatment and control villages reported antimalarial drug use, latrine access, education (based on the highest level of education attained by the father), and home elevations. However, self-reported bednet ownership and fever histories were higher in treatment villages than in control villages (p<0.05) (Table 1).
Univariate Analysis by Time
Overall, the proportion of P. falciparum–infected participants was highest (6%) at baseline. Infections decreased sharply in treated villages and gradually in control villages throughout follow-up period (Figure 2). At the baseline evaluation, 6% (53/854) of participants from AZT MDA treatment villages were positive for P. falciparum compared with 6% (54/894) of participants from AZT MDA control villages. No differences in odds of infection were observed between control and treatment villages (odds ratio 1.03, 95% CI 0.68–1.55) (Figure 3, Table 2).
By month 1, P. falciparum prevalence in AZT MDA treatment villages decreased to 2% (14/851) and P. falciparum prevalence in AZT MDA control villages decreased to 5% (37/779). The odds ratio for AZ MDA treated villages compared with control villages was 0.34 (95% CI 0.17–0.64), which is consistent with a 66% reduction in the odds of P. falciparum infection in AZT MDA treatment participants compared with AZT MDA control participants. Beyond month 1, the association between of AZT MDA and reduced prevalence of malaria infection decreased. However, overall rates of P. falciparum infection decreased below levels for which reliable inferences could be made.
A village level comparison of changes in malaria prevalence for AZT MDA treatment villages versus prevalence in control villages indicated that 1 village in each group had no change in malaria prevalence, and 3 AZ MDA treatment villages had a significant (p<0.05, by Fisher exact test) decrease in malaria prevalence (p = 0.489, p = 0.0002, and p = 0.0005) and 2 in the control group had significant (p = 0.0008 and p = 0.001) increases. Village 5 in the control group had a significant (p = 0.0001) decrease (Figure 4; Technical Appendix Table 1).
In a subgroup of 200 children examined for malaria infection at weeks 2, 6, and 8, the malaria prevalence of 2%–5% was lower than the 10% expected (23) Thus, reliable prevalence comparisons were not relevant. During the study, febrile patients were tested by RDT and microscopy. During the critical period of 5–42 days after AZT MDA, 46 and 48 children in the treatment and control villages, respectively, were evaluated for fever. Among febrile children, only 5 children from treatment villages and 3 children from control villages were positive by real-time PCR.
At month 1, reductions in P. falciparum malaria prevalence were observed by logistic regression models that controlled for differences in self-reported bednet ownership and home altitude (Table 3) and included village level random effects. The odds of P. falciparum infection were 63% (95% CI 28%–81%) but less in treated villages than in control villages after adjusting for bednet ownership and home altitude in a model with village-level random intercepts. Because self-reported history of fever, malaria medication, latrine access, and education of the head of the household were not associated with malaria prevalence and did not appreciably alter the interpretations of the multivariate models, they were not included in the reported estimates. As was evident in univariate analysis, no differences in malaria rates existed at baseline, and the association between AZT MDA and reduced malaria rates decreased after the first month of follow up.
Spatial information on home locations enabled a thorough examination of residual spatial autocorrelation. A total of 51% (19/37) of P. falciparum infections at month 1 occurred in control village 8 (Figure 4), which raised concerns that residual spatial autocorrelation artificially decreased CIs from the reported models despite use of the random effects model. However, a variogram of standardized residuals from the model suggested that residual spatial autocorrelation was appropriately controlled by the village level random effects. The data were reanalyzed with incident infections that had been preceded by a negative real-time PCR result for P. falciparum; results did not change appreciably. We also observed an age-independent decrease in malaria prevalence at month 1 for persons 1–10 years of age (univariate analysis) and for persons >10 years of age (univariate analysis and multivariate analysis) (Technical Appendix Table 2).
Analysis of PfRpL4 Mutations Associated with Azithromycin Resistance
Sequencing of full-length P. falciparum ribosomal L4 protein was performed for samples from 12 patients. We did not find evidence of single-nucleotide polymorphisms (SNPs) conferring AZT resistance in amino acid region 71–79 (1). A synonymous SNP at position K36 was found in all bacterial plasmid clones from samples of 1 patient and in a P. falciparum–positive control patient who was not from Tanzania. A single plasmid clone from a sample of participant from an AZT MDA control village contained a nonsynonymous SNP in the active site of the PfRpL4 gene at A78S. However, 17 other plasmid clones from the same blood sample did not contain the mutation, which suggested the aberrant sequence resulted from a Taq polymerase error during plasmid cloning of the PfRpL4 gene. Moreover, restriction enzyme HindIII, which was specific for the nonsynonymous SNP at A78S, did not digest multiple PCR amplicons from the original patient blood sample, which further suggested that the SNP on the PfRpL4 gene was present at a low concentration or may have been an artifact of the Taq polymerase amplification. A78S has not been identified in other mutant bacteria associated with AZT resistance (1). Eighteen other nonsynonymous mutations were found in samples from participants from AZT MDA treatment villages, but these mutations were not in conserved amino acid regions for PfRpL4 associated with AZT resistance in bacteria or the selected P. falciparum (except for an I39S were single occurrences among multiple clones from a single patient). Four synonymous SNPs were also detected (Table 4).
AZT MDA was associated with a minimum reduction of 66% in odds of P. falciparum infection compared with odds of reduction for controls in the first month after drug administration, as shown by univariate and multivariate models on an individual level. Beyond 1 month, prevalence in treatment and control groups decreased, and there was no difference in prevalence between treatment and control groups. In addition, no mutations were detected in samples from treated persons who had parasitemias after drug treatment in the only gene associated with in vitro selection of AZT resistance. However, these results show more modest effects than those of other studies of AZT MDA, which reported large reductions in malaria (9) and that AZT MDA protection against malaria contributed to broader decreases in illness and death (10).
Because our data were analyzed from the perspective of an intention to treat, participants who did not receive AZT but lived in villages that received AZT MDA were classified with those who received AZT. The intention-to-treat analysis could have biased reported results toward the null because persons in AZT MDA treatment villages who benefited from AZT might have been more likely to become infected (Technical Appendix Tables 3, 4). Because intention-to-treat analysis would probably bias results toward the null, we believe that it was a conservative analytic assumption.
An additional potential confounder is that all febrile children <5 years of age and adults who were positive for P. falciparum by RDT were treated according to national and Integrated Management of Childhood Illness guidelines with artemether/lumefantrine, which has a lingering 4-week prophylactic effect on parasitemia. However, inclusion of artemether/lumefantrine–treated participants did not affect the relevant proportions of malaria infection at month 1 (Technical Appendix Table 5).
AZT treatment occurred at the village level, which might have exacerbated similarities in malaria risk among participants in the same village and created village-level clustering (25). Clustering within villages can bias results if it is not addressed in the statistical models (30). Village-level random effects were used to adjust variance estimates for within-village clustering and provide valid multilevel statistical inferences, given the clustering inherent in the study design (20,31). We argue that a relevant comparison is between the proportion of AZT MDA–treated and control villages, rather than a grouping of 4 villages in each AZT MDA intervention and control villages. Grouping findings into only 8 units on a village level does not produce relevant results and was not the design of the study.
AZT MDA was conducted occurred according to the WHO SAFE (surgery, antibiotics, facial cleanliness, and environmental changes) trachoma control strategy, and assigning treatment to villages by randomization would have been unethical because of the established effectiveness of this strategy in the study region. AZT MDA was examined from an observational perspective with the longitudinal cohort. We sought to minimize lack of comparability by selecting villages from a similar geographic region and by controlling for bednet ownership and altitude in multivariate analysis. We believe that residual confounding and bias are minimal because similar P. falciparum infection prevalence rates were observed at baseline and at months 3–6 after the effect of the AZT MDA had decreased.
In Ethiopia, Porco et al. reported a 49% (95% CI 10%–71%) decrease in odds of death in a randomized trial of AZ MDA villages compared with control villages that received a single dose of AZ (20 mg/kg to ≤1 g) (10). AZT protection from P. falciparum infection might have contributed to the overall decrease in deaths. Although the effect of AZT MDA on prevalent infections was strong, most infections observed were asymptomatic. The study in Ethiopia also had 3 dosing strategies (annual, biannual, and quarterly), which resulted in two thirds of study participants receiving more annual doses than used in the study in Tanzania. Our data provide only weak evidence that a reduction in malaria infection contributed to the reduction in deaths reported in Ethiopia.
At the study site, we did not find evidence of P. falciparum AZT resistance markers on a gene previously implicated in AZT resistance in vitro. We also observed a short-lived but major reduction in malaria. AZT MDA does not produce Plasmodium spp. drug resistance probably because other drugs with different mechanisms of action are used to treat malaria. This study might provide an optimistic note for AZ MDA planners and those living in areas in developing nations to which malaria and trachoma are endemic.
Dr Schachterle is an associate director at Pfizer, Inc., New York, NY. His primary research interest is spatial analysis of infectious diseases.
We thank Harran Mkocha, Diane Stare, and the PretPlus study team for providing tireless efforts during this study; Edward Sambu and Fikirini Msuya for assisting with training of local field staff and reading peripheral blood smears; Stephen Magesa for providing technical support; and Natalie Sanders for assisting with DNA sequencing.
This study was supported by the Bill and Melinda Gates Foundation (grant 48027); Research to Prevent Blindness (unrestricted grant); and the National Institute of Allergy and Infectious Diseases (grant U01 AI068613: (HIV Prevention Trials Network–Laboratory Network). Pfizer, Inc. donated the azithromycin for the study. S.W. was supported by a Senior Scientific Award from Research to Prevent Blindness.
- Sidhu AB, Sun Q, Nkrumah LJ, Dunne MW, Sacchettini JC, Fidock DA. In vitro efficacy, resistance selection, and structural modeling studies implicate the malarial parasite apicoplast as the target of azithromycin. J Biol Chem. 2007;282:2494–504.
- Duran JM, Amsden GW. Azithromycin: indications for the future? Expert Opin Pharmacother. 2000;1:489–505.
- Friesen J, Silvie O, Putrianti ED, Hafalla JC, Matuschewski K, Borrmann S. Natural immunization against malaria: causal prophylaxis with antibiotics. Sci Transl Med. 2010;2:40ra49 .
- Shimizu S, Osada Y, Kanazawa T, Tanaka Y, Arai M. Suppressive effect of azithromycin on Plasmodium berghei mosquito stage development and apicoplast replication. Malar J. 2010;9:73 .
- Burton MJ, Mabey DC. The global burden of trachoma: a review. PLoS Negl Trop Dis. 2009;3:e460.
- Mathew AA, Turner A, Taylor HR. Strategies to control trachoma. Drugs. 2009;69:953–70.
- West SK, Munoz B, Mkocha H, Gaydos CA, Quinn TC. Number of years of annual mass treatment with azithromycin needed to control trachoma in hyper-endemic communities in Tanzania. J Infect Dis. 2011;204:268–73.
- Polack S, Brooker S, Kuper H, Mariotti S, Mabey D, Foster A. Mapping the global distribution of trachoma. Bull World Health Organ. 2005;83:913–9 .
- Sadiq ST, Glasgow KW, Drakeley CJ, Muller O, Greenwood BM, Mabey DC, Effects of azithromycin on malariometric indices in The Gambia. Lancet. 1995;346:881–2.
- Porco TC, Gebre T, Ayele B, House J, Keenan J, Zhou Z, Effect of mass distribution of azithromycin for trachoma control on overall mortality in Ethiopian children: a randomized trial. JAMA. 2009;302:962–8.
- Keenan JD, Ayele B, Gebre T, Zerihun M, Zhou Z, House JI, Childhood mortality in a cohort treated with mass azithromycin for trachoma. Clin Infect Dis. 2011;52:883–8.
- Sykes A, Hendriksen I, Mtove G, Mandea V, Mrema H, Rutta B, Azithromycin plus artesunate versus artemether-lumefantrine for treatment of uncomplicated malaria in Tanzanian children: a randomized, controlled trial. Clin Infect Dis. 2009;49:1195–201.
- van Eijk AM, Terlouw DJ. Azithromycin for treating uncomplicated malaria. Cochrane Database Syst Rev. 2011;2:CD006688 .
- Andersen SL, Oloo AJ, Gordon DM, Ragama OB, Aleman GM, Berman JD, Successful double-blinded, randomized, placebo-controlled field trial of azithromycin and doxycycline as prophylaxis for malaria in western Kenya. Clin Infect Dis. 1998;26:146–50.
- Heppner DG Jr, Walsh DS, Uthaimongkol N, Tang DB, Tulyayon S, Permpanich B, Randomized, controlled, double-blind trial of daily oral azithromycin in adults for the prophylaxis of Plasmodium vivax malaria in western Thailand. Am J Trop Med Hyg. 2005;73:842–9 .
- Taylor WR, Richie TL, Fryauff DJ, Picarima H, Ohrt C, Tang D, Malaria prophylaxis using azithromycin: a double-blind, placebo-controlled trial in Irian Jaya, Indonesia. Clin Infect Dis. 1999;28:74–81.
- Taylor WR, Richie TL, Fryauff DJ, Ohrt C, Picarima H, Tang D, Tolerability of azithromycin as malaria prophylaxis in adults in northeast Papua, Indonesia. Antimicrob Agents Chemother. 2003;47:2199–203.
- Coles CL, Levens J, Seidman JC, Mkocha H, Munoz B, West S. Mass distribution of azithromycin for trachoma control is associated with short-term reduction in risk of acute lower respiratory infection in young children. Pediatr Infect Dis J. 2012;31:341–6.
- Coles CL, Seidman JC, Levens J, Mkocha H, Munoz B, West S. Association of mass treatment with azithromycin in trachoma-endemic communities with short-term reduced risk of diarrhea in young children. Am J Trop Med Hyg. 2011;85:691–6.
- Schachterle SE, Mtove G, Levens JP, Clemens EG, Shi L, Raj A, Prevalence and density-related concordance of three diagnostic tests for malaria in a region of Tanzania with hypoendemic malaria. J Clin Microbiol. 2011;49:3885–91.
- Taylor R, Todd M, Kongloa L, Maurice L, Nahozya E, Sanga H, Evidence of the dependence of groundwater resources on extreme rainfall in east Africa. Nature Climate Change. 2013;3:374–8.
- World AgroMeteorological Information Service, Southern African Development Community (SADC). SADC Agromet Bulletin Jan 21–31, 2009 [cited 2014 Feb 24]. http://wamis.gmu.edu/countries/sadc_rev.php
- National Bureau of Statistics. Tanzania HIV/AIDS and malaria indicator survey: key findings. Calverton (MD): The Bureau; 2009.
- Reller ME, Clemens EG, Schachterle SE, Mtove GA, Sullivan DJ, Dumler JS. Multiplex 5′ nuclease-quantitative PCR for diagnosis of relapsing fever in a large Tanzanian cohort. J Clin Microbiol. 2011;49:3245–9.
- Waller L, Gotway C. Applied spatial statistics for public health data. Hoboken (NJ): Wiley-IEEE; 2004.
- Diggle P, Ribero P Jr, Christensen O. An introduction to model-based geostatistics. Spatial Statistics and Computational Methods. 2003;173:43–86.
- R Development Core Team. R: a language and environment for statistical computing. Vienna (Austria): R Foundation for Statistical Computing; 2008.
- Rowlingson B, Diggle P. SPLANCS: spatial and space-time point pattern analysis. R package version 2.01–27; 2010 [cited 2014 Feb 24]. https://r-forge.r-project.org/scm/viewvc.php/pkg/splancs/DESCRIPTION?pathrev=1194&sortby=date&root=rspatial&sortdir=down&view=diff&r1=1111&r2=1112
- Bates D, Maechler M, Dai D. lme4:Linear mixed-effects models using S4 classes. R package version 0.999375–42; 2008 [cited 2014 Feb 24]. http://lme4.r-forge.r-project.org/
- von Seidlein L, Greenwood BM. Mass administrations of antimalarial drugs. Trends Parasitol. 2003;19:452–60.
- Diggle P, Heagerty P, Liang K, Zeger S. Analysis of longitudinal data. 2nd ed. New York: Oxford University Press; 2002.
TablesCite This Article
1Current affiliation: Pfizer, Inc., New York, New York, USA.