What is required in terms of mass drug administration to interrupt the transmission of schistosome parasites in regions of endemic infection?
© Anderson et al. 2015
Received: 15 September 2015
Accepted: 7 October 2015
Published: 22 October 2015
Schistosomiasis is endemic in 54 countries, but has one of the lowest coverages by mass drug administration of all helminth diseases. However, with increasing drug availability through donation, the World Health Organisation has set a goal of increasing coverage to 75 % of at-risk children in endemic countries and elimination in some regions. In this paper, we assess the impact on schistosomiasis of the WHO goals in terms of control and elimination.
We use an age-structured deterministic model of schistosome transmission in a human community and the effect of mass drug administration. The model is fitted to baseline data from a longitudinal re-infection study in Kenya and validated against the subsequent re-infection data. We examine the impact on host worm burden of the current treatment trend, extrapolated to meet the WHO goals, and its sensitivity to uncertainty in important parameters. We assess the feasibility of achieving elimination.
Model results show that the current treatment trend, extrapolated to the WHO goals, is able to greatly reduce host worm burdens. If coverage is continued at the same level beyond 2020, elimination is possible for low to moderate transmission settings, where transmission intensity is defined by the basic reproduction number, R0. Low levels of adult coverage have a significant impact on worm burden in all settings. Model validation against the re-infection survey demonstrates that the age-structured model is able to match post-treatment data well in terms of egg output, but that some details of re-infection among school children and young adults are not currently well represented.
Our work suggests that the current WHO treatment goals should be successful in bringing about a major reduction in schistosome infection in treated communities. If continued over a 15 year period, they are likely to result in elimination, at least in areas with lower transmission.
The third progress report of the London Declaration records that schistosomiasis remains red in the progress score card chart (recently developed by Uniting to Combat NTDs) as it has the lowest coverage of all helminth diseases treatable by mass drug administration (MDA) at 14.4 % in 2012 and 15.6 % in 2013 . In addition, new mapping of schistosomiasis in affected countries is increasing the number of identified districts with endemic infection. Only twenty-six countries (50 %) of 52 endemic countries reported MDA in 2013 . The report argues that significant improvements could be made in the coming years as drug supply of praziquantel is expected to increase and the launching of the new Global Schistosomiasis Alliance will increase collaboration within this disease community to help countries scale up MDA of children and adults .
Global treatment numbers for schistosomiasis have slightly improved in recent years. In 2013, 47.3 million people were treated with praziquantel, while the comparative figure was 42 million in 2012. However, donations of praziquantel are increasing significantly. In 2014, pharmaceutical donor Merck increased its contribution to nearly 75 million tablets, and of 41 African countries requiring treatment, 36 (88 %) were reached. Merck has increased its donation in 2015 to 100 million tablets. For 2016, Merck has committed to increasing its donation of praziquantel up to 250 million tablets, equivalent to 100 million treatments. As reported in the third progress report, DFID, USAID, the World Bank, and World Vision also purchased over 100 million of praziquantel tablets in 2015 .
In this paper, we address what might be required to eliminate transmission for Schistosoma mansoni using annual MDA coverage in pre-SAC, SAC and adults. We also assess the impact on host worm burden of historic coverage levels, projected forward to 75 % in SAC by 2020. Only in six of 31 reporting countries was the target of reaching at least 75 % of school children achieved .
To provide quantitative predictions of MDA impact we employ deterministic age structured mathematical models that incorporating sexual reproduction, aggregated distributions of worm numbers per host (described by the negative binomial distribution with aggregation parameter k), density dependent egg production and variable contact with, and production of, infective stages to continue the life cycle, by different age groupings of the host population, defined as pre-school aged children (pre-SAC, 0–4 years of age), school aged children (SAC, 5 to 14 years of age), and adults (> 15 years of age). The model template is described in Anderson & May  and modifications to this structure to include variable exposure to infection and production of eggs by host age groupings are detailed in Truscott et al. .
Details of the mathematical model employed are described in previous publications [6–8] and also in Additional file 1. Most research on model development for schistosome transmission is based on the pioneering studies of Macdonald in the 1960s . He was the first to point out that sexual reproduction in dioecious species created a breakpoint in transmission (an unstable equilibrium separating the two stable points of endemic infection and parasite extinction) below which mating frequency is insufficient to sustain transmission. More recent modifications include the inclusion of negative binomial distributions of parasite number per host [10, 11]. Worm counts based on autopsy studies show clearly that parasites are highly aggregated within the host population where most hosts harbour few parasites and a few harbour many . These autopsy results are also confirmed by plots of the prevalence of infection versus mean epg counts which also show high degrees of aggregation [13, 14].
Other important modifications in past work include the demography of the human host population , the impact of MDA [16–18], acquired immunity (dependent on past exposure to infection) , spatial structure and host gender , density dependence in egg production , predisposition to heavy infection , stochastic individual based models  and the generation of parasite aggregation through heterogeneity in exposure amongst the host population [20, 21].
In this paper we employ a deterministic hybrid model which assumes the parasite is dioecious and monogamous, has density dependent egg production  and a degree of parasite aggregation defined by the negative binomial probability distribution with a fixed value of the aggregation parameter k. The details of the model can be found in the supplementary information. The within-host section of the model describes the evolution of the worm burden in individuals as a function of age. It has full age structure, but outputs are grouped into three age groupings for pre-SAC (0–4 years of age), SAC (5–14 years of age) and adults (>15 years of age), reflecting the key age distinctions made within MDA programs. We include variability in host exposure to infection and contributions to egg production dependent on host age group (described by a parameter β i which describes the intensity of infectious contact of age group i with the infectious material in the environment. Because of the relative fast time scales of turnover in the larval stages and snail intermediate host of schistosomes (see review by Anderson & May ), relative to the long life expectancy of the adult worm in the human host, the details of the dynamics of the free living and snail based stages of the parasite are collapsed into an equation for infectious material in the environment.
Parameter values used in the numerical evaluations of model predictions for S mansoni. The age intensity profile from the cited sources based on faecal egg counts were used to derive estimate of R0 and βi using maximum likelihood methods. The age groupings for the contact parameter βi are defined as the intervals between 0, 5, 10, 16 and 80 years
Basic Reproductive number, R0
Adult worm life expectancy
Negative binomial clumping parameter, k
Density dependence fecundityparameter, ɣ
Beta, βi for intensity of age group i’s contact with infectious material
Drug efficacy as a proportion of worms killed by praziquantel
Egg output per female worm
Source of age - mean intensity data
Ietune village, Machakos District, Kenya 
The issue of whether or not the convex age intensity curves recorded for the schistosomes are due to acquired immunity or age related exposure remains uncertain. The problem of observed epidemiological pattern interpretation was first enunciated by Warren in 1973 who described the problem succinctly as one of differentiating between ecology and immunology . The conundrum in interpretation remains today, although there is abundant clinical evidence of the mounting of strong and specific immunological responses to schistosome infections by the human host . In our analyses of the impact of MDA, we assume that the convex age intensity profiles are generated by age dependent exposure. Future work will assess the impact of weakening this assumption to include the more realistic premise that both ecology and immunology act as determinants.
Treatment coverage is assumed to be at random at each round of MDA in each age group. This is unlikely to be true in reality and a future publication will examine adherence issues within a stochastic framework. If those not treated are predisposed to this state for social or behavioural factors, the impact of MDA at 75 % coverage will be less than predicted in the results discussed in this paper.
Host demographic details (age dependent death rates and population pyramids) were taken from published data on sub-Saharan countries with the profile in Uganda chosen to represent a typical profile for sub-Saharan Africa where S. mansoni is endemic .
Validation of model predictions is essential and can be performed in variety of ways. Qualitative comparisons of model prediction with observed pattern of reinfection post mass treatment is one approach. Ideally, however, quantitative comparisons of predicted outcome of the impact of MDA require parameter estimation using data from a defined population with little or no recent past treatment and then the comparison of expected and observed patterns of reinfection at subsequent rounds. Given the age structured nature of the current model, a detailed validation requires precise information on who has been treated, at what age, how often and their representation within the final reinfection survey. For schistosomes , worm loads are usually measured indirectly by eggs per gram of faeces (epg) or per unit volume of urine in the case of S haematobium. This measure is known to be unreliable, typically underestimating the true prevalence. Furthermore, in estimating intensity, it has a high variance as recorded in repeated samples from the same stool and in stools taken on consecutive days . New tools, such as quantitative PCR, offer some hope for the future in providing precise scores for prevalence, and perhaps even egg burdens in stools. A further complication arises from density dependence in egg production by female worms which can result in higher per capita egg output post the initial round of treatment during the early phases of reinfection.
Model validation has rarely been attempted for helminth infections, but it is a central need for ongoing and future epidemiological studies of MDA impact, if predictions on coverage requirements to interrupt transmission are to carry weight amongst public health policy makers. A previous attempt by French et al. (2010) used pre and post treatment epg counts from a study of S mansoni in Uganda recording reinfection over 3 years of follow up, but validation was based on estimating parameter values from the trends in reinfection . Ideally, it should be done by estimating these parameters first, and then predicting reinfection pattern to compare with observed outcomes given data on treatment coverage.
We report a validation exercise using a reinfection study of S mansoni recording prevalence and mean epgs in a set of villages in Kenya that had not previously been treated prior to the initiation of mass drug administration . We also examine the sensitivity of model predictions of the impact of achieving 75 % coverage in SAC to changes of R0, to different values of the contact parameter βi and to the level of drug coverage especially in adults.
Sensitivity analyses of model projections to changes in the values of key parameters. The numerical values record the percentage reduction in the predicted mean intensity of infection in each age group from a rising treatment coverage to 75 %, and to 27.4 % in adults, in 2020 as compared to baseline in 2003
Pre-SAC mean burden
SAC mean burden
Adults mean burden
Overall population mean burden
Fitted parameter values (Table 1)
Higher and lower transmission setting (R0 + 1.0, −0.25)
49 %, 81 %
78 %, 94 %
76 %, 86 %
76 %, 90 %
Coverage level achieved in 2020 (60 %, 90 %)
77 %, 67 %
94 %, 85 %
87 %, 79 %
90 %, 81 %
Degree of exposure and contribution of children to the infectious reservoir (5–16 years) (β2 and β2 + 50 %, −50 %)
73 %, 71 %
91 %, 90 %
84 %, 82 %
87 %, 85 %
We performed a basic sensitivity analysis on the impact of the treatment coverage trend with respect to some of the more important or uncertain parameters, the results of which are summarised in Table 2. Transmission intensity, as captured by R0 in our model, is encapsulates a large number of features of the transmission cycle and is likely to vary from community to community. The value obtained from the current fit (R0 = 1.7) is relatively low with respect to previous estimates [5, 31]. An increase of R0 by 1 decreases the impact of treatment by 2020 considerably. However, examination of the time series suggests that elimination is still possible, but delayed beyond 2030.
The coverage trend shown in Fig. 2 is an average across many countries in which schistosomiasis is endemic. Across individual communities, however, levels of coverage are likely to vary considerably. Table 2 shows that the impact of the 2020 level of coverage of children is relatively small and also linear across all age groups. This indicates that, even though coverage may vary across treated communities in a country, average coverage levels will give a good estimate of the average effect on worm burden in the population as a whole.
Sensitivity of model projections (percentage reduction in the predicted mean intensity of infection) to changes in the level of coverage of adults achieved in 2020 (linear trend from last data point in 2013)
Adult coverage by 2020
Pre-SAC mean burden
SAC mean burden
Adults mean burden
Overall population mean burden
Much of the shape of the elimination surface can be understood with reference to the baseline age profile shown in Fig. 4. For example, egg output from the pre-SAC age group is -quite low and this corresponds with the low sensitivity of elimination to coverage in pre-SAC (the surface is very flat in the direction of changing pre-SAC coverage). Treatment of pre-SAC is not currently part of the WHO goals for schistosomiasis. The drug is not authorised for use in this age group (see Methods section), although efforts are being made to develop a suitable formulation . However, our results suggest that they would be of limited use in bringing about elimination, although they could be important in reducing morbidity. Coverage of adults is important at low levels, enabling elimination to be achieved within the WHO goals for SAC. Further increase in adult coverage has little additional effect due to the relatively low worm burdens among adults. The impact of treating adults is reduced if a longer period of treatment is recovered, showing that the treating adults reduces the time to elimination rather than the long-term possibility of elimination.
The close association between the baseline age-profile and the elimination surface suggests that if the baseline varies between populations, the conditions under which elimination can be achieved will vary too. For example, other age profiles have a stronger presence in worms within the adult population and this will translate into a greater sensitivity to adult coverage.
The predictions reported in this paper are for S mansoni, but given that the life expectancy of S heamatobium and its age-profile within the host population are similar, the predictions should apply for that species too. In the case of S. japonicum, the added complexity of multiple zoonotic reservoirs will likely mean that conclusions drawn from S. mansoni cannot be directly transferred to the other species that infest humans.
As outlined in the Methods section, validation of model predictions requires comparison of observed outcomes with corresponding post-treatment predictions, given baseline initial conditions and parameter estimates and data on who is treated, how frequently and who is surveyed. Given the nature of our model, ideal data would be individual longitudinal egg count data with matching data on age and treatment history.
Few published epidemiological studies on schistosomes have this level of detail, at least in their published form. The study by Butterworth et al. of S mansoni transmission and treatment in Machakos district, Kenya offers a good compromise in that it reports age group-stratified baseline and re-infection data across several study arms with data on MDA coverage and timing over a number of years [24, 26, 34]. We focus on the study arm in Iietune village. In this population, an initial survey in 1983 was followed immediately by a round of MDA. A second round occurred in 1985 and a re-infection survey was completed in 1987. A coverage of 93 % was achieved in both rounds, assumed for the purpose of validation to be coverage of those eligible for treatment (5 years or over). We have extracted epg data from the baseline and the re-infection survey, with the baseline data being used to parameterize the model using maximum likelihood techniques.
The general conclusion generated by the model predictions is that an annual 75 % coverage in SAC and a modest level of coverage in adults should reduce mean intensities of infection to very low levels by 2020. In low to moderate transmission settings continuation of this coverage for a decade or so should result in transmission elimination. For high transmission settings, however, somewhat higher levels of coverage will be required. These predictions seem to be relatively insensitive to changes in the parameter estimates explored in Table 2.
The key question, however, is whether coverage can be raised to 75 % in SAC and modest levels in adults? Concomitantly, can drug production levels also be raised to meet such demands? These questions are beyond the scope of this paper and the resolution of them lies with WHO, policy makers, implementers within countries and the pharmaceutical companies who makes the major donations of praziquantel.
Model validation against reinfection data indicate that, while broadly in agreement with the data, there are clear discrepancies. These arise from two equally important sources. On the one hand, the failure to capture the behaviour of the population in the 12–25 year age range suggests that a more sophisticated model of infectious contact among children is required. On the other, there is a lot of missing and uncertain material in the data being used in this validation, such as past history of treatment (which we have assumed no past history). The availability of high quality data is a major issue for model testing and validation. Far too few published studies record data on intensity and prevalence by age group (ideally individual person epg counts) electronically and make such data freely available to other researchers. This is an important need in the coming years as MDA coverage is hopefully increased.
As it currently stands, our model does not include the effects of immunity. Immune responses are not well understood, but possible mechanisms include acquired immunity from past exposure to infection and a long-lived immune response generated by death of adult worms [35, 36]. Without a clear understanding of the mechanisms of acquired immunity, it is very hard to disentangle the effects of immunity from age-dependent exposure to infection (β). However, re-infection data, stratified by age, may offer a chance to distinguish the two processes. Our current model overestimates reinfection in an age-group that was heavily infected at baseline, perhaps suggesting that an immune response triggered by treatment may be responsible. Future work will examine this possibility.
Another potential source of error in our model is the treatment of host demography. The data used for both fitting and validation is from studies conducted in the 1980s, but our modelling concerns the past decade and near future. Within this interval, significant changes in demography have certainly taken place. We use a frequency-dependent model of infection whose behaviour will not be effected by the absolute host population (whether absolute population size effects disease dynamics in reality will require detailed analysis of datasets from different time periods). More significant are changes in the age distribution of the host population which will affect the mean force of infection experienced by individuals and their contributions to the infectious reservoir. Such effects can be investigated with our model, but will require high quality infection data from a range of different demographies.
Broader scientific issues relate to the uncertainties on biological questions, the quality of data available for parameter estimation, data on drug coverage levels and the impact of heterogeneities in transmission. Biological uncertainties are key in two areas. The first is the question of monogamy for adult parasites. The calculations recorded in this paper will be extended in future work to look how the assumption of polygamy influences predictions. In rough terms it will make elimination of transmission a little more difficult (raising the surface in Fig. 7).
The second biological question relates to the importance of acquired immunity. Again, this requires further exploration. Previous studies have outline various model frameworks for helminths to address this issue using infection functions that depend on past exposure to infection and these will be employed in future analyses. However, it is possible to speculate on what influence acquired immunity may have . Broadly speaking, if immunity is important and builds up through early life exposure, then repeated MDA will act to reduce herd immunity in a manner that will increase reinfection rates after many rounds of treatment as children lack past experience of infection. This would act to lessen the long term impact of MDA and require the elimination of transmission to make sure that, if treatment stopped (or decreased), the population was not more susceptible to infection and associated morbidity.
Data availability to test predictions against observations is good in one context and poor in others. Epidemiological data, well stratified by age and sex, on reinfection post rounds of MDA is good, as well illustrated in Fig. 3 by the studies of Butterworth and colleagues [17, 26]. Such data facilitates model validation.
Data is poor for the estimation of key population dynamic parameters such as density dependent fecundity and the relationship between epg measures and worm counts because direct measurements of host worm populations in vivo is impossible. Estimates can be obtained from experimental animal models (mice and some primate species [37, 38]) as well as from autopsy studies, although these are very rare and probably constitute a biased sample . The alternative source for estimates for these parameters is through fitting complex models to extensive datasets . In this sense, model development and parameterisation for schistosome species is much more difficult than for soil transmitted helminths, where many worm expulsion studies have been published [40–42].
Drug coverage data for praziquantel treatment stratified by the age group treated is poor at present. This could be easily remedied by requesting countries to record this information as treatment coverage expands. This data is essential for making accurate predictions, given that for all schistosome species, a significant fraction of the total adult parasite population is harboured by those over 15 years of age .
The transmission of schistosome species relies on human contact with water, and as such much heterogeneity exists between endemic settings in the exposure of individuals by age and gender as a consequence of the prevailing cultural, social and behavioural conditions. Our predictions were based on parameter values derived from one medium intensity transmission setting, but in the villages studied by Fulford et al. in Kenya much variation in exposure to infection was apparent village by village, relating to the type of water contact and what sort of aquatic habitat was present near each village . Addressing such heterogeneity requires the use of individual based stochastic models incorporating spatial transmission, and these will be a priority in future research.
A stochastic framework will also facilitate addressing other important sources of heterogeneity, such as predisposition to heavy infection , and variability in, or predisposition to, treatment at each round of MDA. Data requirements in this context include the keeping of accurate records of precisely who receives treatment at each round of MDA.
In conclusion, we return to the main prediction of the model calculations. S. mansoni control, and even the interruption of transmission, is possible if high levels of annual treatment coverage in SAC are reached and maintained, plus some modest coverage in adults of around 30 %. Reaching the 75 % coverage in SAC in all endemic regions is the challenge for the coming years.
Research for this article was funded by The Children’s Investment Fund Foundation (UK) (“CIFF”). The views, opinions, assumptions or any other information set out in this paper are solely those of the authors and should not be attributed to CIFF or any person connected with CIFF.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Country leadership and collaboration on neglected tropical diseases: third progress report of the London Declaration [http://unitingtocombatntds.org/sites/default/files/document/UTCNTD%20FULL%20REPORT.pdf]
- Global Schistosomiasis Alliance [http://www.eliminateschisto.org/]
- WHA resolution 54.19 Schistosomiasis and soil-transmitted helminth infections [http://www.who.int/neglected_diseases/mediacentre/WHA_54.19_Eng.pdf]
- WHO: Investing to overcome the global impact of neglected tropical diseases: third WHO report on neglected tropical diseases. In. Edited by Holmes P; 2015: 191.Google Scholar
- Anderson RM, May RM. Infectious diseases of humans : dynamics and control. Oxford. New York: Oxford University Press; 1991.Google Scholar
- Truscott JE, Hollingsworth TD, Brooker SJ, Anderson RM. Can chemotherapy alone eliminate the transmission of soil transmitted helminths? Parasites & Vectors. 2014;7(1):266.View ArticleGoogle Scholar
- Anderson RM, Truscott JE, Hollingsworth TD. The coverage and frequency of mass drug administration required to eliminate persistent transmission of soil-transmitted helminths. Philos Trans R Soc Lond B Biol Sci. 2014;369(1645):20130435.PubMed CentralView ArticlePubMedGoogle Scholar
- Truscott J, Hollingsworth TD, Anderson R. Modeling the interruption of the transmission of soil-transmitted helminths by repeated mass chemotherapy of school-age children. PLoS Negl Trop Dis. 2014;8(12), e3323.PubMed CentralView ArticlePubMedGoogle Scholar
- Macdonald G. The dynamics of helminth infections, with special reference to schistosomes. Trans R Soc Trop Med Hyg. 1965;59(5):489–506.View ArticlePubMedGoogle Scholar
- Anderson RM, May RM. Herd immunity to helminth infection and implications for parasite control. Nature. 1985;315(6019):493–6.View ArticlePubMedGoogle Scholar
- May RM. Togetherness among Schistosomes: its effects on the dynamics of the infection. Math Biosci. 1977;35:301–43.View ArticleGoogle Scholar
- Cheever AW, Kamel IA, Elwi AM, Mosimann JE, Danner R. Schistosoma mansoni and S. haematobium infections in Egypt. II. Quantitative parasitological findings at necropsy. Am J Trop Med Hyg. 1977;26(4):702–16.PubMedGoogle Scholar
- Bensted-Smith R, Anderson RM, Butterworth AE, Dalton PR, Kariuki HC, Koech D, et al. Evidence for predisposition of individual patients to reinfection with Schistosoma mansoni after treatment. Trans R Soc Trop Med Hyg. 1987;81(4):651–4.View ArticlePubMedGoogle Scholar
- Guyatt HL, Smith T, Gryseels B, Lengeler C, Mshinda H, Siziya S, et al. Aggregation in schistosomiasis: comparison of the relationships between prevalence and intensity in different endemic areas. Parasitology. 1994;109(Pt 1):45–55.View ArticlePubMedGoogle Scholar
- Anderson RM, May RM. Population dynamics of human helminth infections: control by chemotherapy. Nature. 1982;297(5867):557–63.View ArticlePubMedGoogle Scholar
- Anderson RM, May RM. Helminth infections of humans: mathematical models, population dynamics, and control. Adv Parasitol. 1985;24:1–101.View ArticlePubMedGoogle Scholar
- French MD, Churcher TS, Gambhir M, Fenwick A, Webster JP, Kabatereine NB, et al. Observed reductions in Schistosoma mansoni transmission from large-scale administration of praziquantel in Uganda: a mathematical modelling study. PLoS Negl Trop Dis. 2010;4(11), e897.PubMed CentralView ArticlePubMedGoogle Scholar
- Tingley GA, Butterworth AE, Anderson RM, Kariuki HC, Koech D, Mugambi M, et al. Predisposition of humans to infection with Schistosoma mansoni: evidence from the reinfection of individuals following chemotherapy. Trans R Soc Trop Med Hyg. 1988;82(3):448–52.View ArticlePubMedGoogle Scholar
- Gurarie D, King CH. Heterogeneous model of schistosomiasis transmission and long-term control: the combined influence of spatial variation and age-dependent factors on optimal allocation of drug therapy. Parasitology. 2005;130(Pt 1):49–65.View ArticlePubMedGoogle Scholar
- Anderson RM, Medley GF. Community control of helminth infections of man by mass and selective chemotherapy. Parasitology. 1985;90:629–60.View ArticlePubMedGoogle Scholar
- Gurarie D, King CH. Population biology of Schistosoma mating, aggregation, and transmission breakpoints: more reliable model analysis for the end-game in communities at risk. PLoS One. 2014;9(12), e115875.PubMed CentralView ArticlePubMedGoogle Scholar
- Number of children requiring PCT for STH [http://apps.who.int/neglected_diseases/ntddata/sth/sth.html]
- WHO: Preventive Chemotheraphy in Human Helminthiasis. In. Geneva, Switzerland: World Health Organization; 2006.Google Scholar
- Butterworth AE, Sturrock RF, Ouma JH, Mbugua GG, Fulford AJ, Kariuki HC, et al. Comparison of different chemotherapy strategies against Schistosoma mansoni in Machakos District, Kenya: effects on human infection and morbidity. Parasitology. 1991;103(Pt 3):339–55.View ArticlePubMedGoogle Scholar
- Fulford AJ, Butterworth AE, Sturrock RF, Ouma JH. On the use of age-intensity data to detect immunity to parasitic infections, with special reference to Schistosoma mansoni in Kenya. Parasitology. 1992;105(Pt 2):219–27.View ArticlePubMedGoogle Scholar
- Fulford AJ, Butterworth AE, Ouma JH, Sturrock RF. A statistical approach to schistosome population dynamics and estimation of the life-span of Schistosoma mansoni in man. Parasitology. 1995;110(Pt 3):307–16.View ArticlePubMedGoogle Scholar
- Warren KS. Regulation of the prevalence and intensity of schistosomiasis in man: immunology or ecology? J Infect Dis. 1973;127(5):595–609.View ArticlePubMedGoogle Scholar
- Siddiqui AA, Siddiqui BA, Ganley-Leal L. Schistosomiasis vaccines. Hum Vaccin. 2011;7(11):1192–7.PubMed CentralView ArticlePubMedGoogle Scholar
- Pullan RL, Kabatereine NB, Quinnell RJ, Brooker S. Spatial and genetic epidemiology of hookworm in a rural community in Uganda. PLoS Negl Trop Dis. 2010;4(6), e713.PubMed CentralView ArticlePubMedGoogle Scholar
- McCullough FS, Bradley DJ. Egg output stability and the epidemiology of Schistosoma haematobium. I. Variation and stability in Schistosoma haematobium egg counts. Trans R Soc Trop Med Hyg. 1973;67(4):475–90.View ArticlePubMedGoogle Scholar
- Woolhouse ME, Hasibeder G, Chandiwana SK. On estimating the basic reproduction number for Schistosoma haematobium. Trop Med Int Health. 1996;1(4):456–63.View ArticlePubMedGoogle Scholar
- Truscott JE, Turner HC, Farrell SH, Anderson RM. Soil Transmitted Helminths: mathematical models of transmission, the impact of mass drug administration and transmission elimination criteria. In: Anderson RM, Basanez MG, editors. Mathematical Models for Neglected Tropical Diseases: Essential Tools for Control and Elimination, vol. Volume B. 2015.Google Scholar
- Pediatric praziquantel consortium [http://www.pediatricpraziquantelconsortium.org/]
- Chan MS, Guyatt HL, Bundy DA, Booth M, Fulford AJ, Medley GF. The development of an age structured model for schistosomiasis transmission dynamics and control and its validation for Schistosoma mansoni. Epidemiol Infect. 1995;115(2):325–44.PubMed CentralView ArticlePubMedGoogle Scholar
- Chan MS, Anderson RM, Medley GF, Bundy DA. Dynamic aspects of morbidity and acquired immunity in schistosomiasis control. Acta Trop. 1996;62(2):105–17.View ArticlePubMedGoogle Scholar
- Mitchell KM, Mutapi F, Savill NJ, Woolhouse ME. Protective immunity to Schistosoma haematobium infection is primarily an anti-fecundity response stimulated by the death of adult worms. Proc Natl Acad Sci U S A. 2012;109(33):13347–52.PubMed CentralView ArticlePubMedGoogle Scholar
- Crombie JA, Anderson RM. Population dynamics of Schistosoma mansoni in mice repeatedly exposed to infection. Nature. 1985;315(6019):491–3.View ArticlePubMedGoogle Scholar
- Zhang W, Ahmad G, Torben W, Noor Z, Le L, Damian RT, et al. Sm-p80-based DNA vaccine provides baboons with levels of protection against Schistosoma mansoni infection comparable to those achieved by the irradiated cercarial vaccine. J Infect Dis. 2010;201(7):1105–12.PubMed CentralView ArticlePubMedGoogle Scholar
- Cheever AW. A quantitative post-mortem study of Schistosomiasis mansoni in man. Am J Trop Med Hyg. 1968;17(1):38–64.PubMedGoogle Scholar
- Bradley M, Chandiwana SK, Bundy DA, Medley GF. The epidemiology and population biology of Necator americanus infection in a rural community in Zimbabwe. Trans R Soc Trop Med Hyg. 1992;86(1):73–6.View ArticlePubMedGoogle Scholar
- Bundy DA, Thompson DE, Cooper ES, Golden MH, Anderson RM. Population dynamics and chemotherapeutic control of Trichuris trichiura infection of children in Jamaica and St. Lucia. Trans R Soc Trop Med Hyg. 1985;79(6):759–64.View ArticlePubMedGoogle Scholar
- Elkins DB, Haswell-Elkins M, Anderson RM. The epidemiology and control of intestinal helminths in the Pulicat Lake region of Southern India. I. Study design and pre- and post-treatment observations on Ascaris lumbricoides infection. Trans R Soc Trop Med Hyg. 1986;80(5):774–92.View ArticlePubMedGoogle Scholar
- Anderson RM, Truscott JE, Pullan RL, Brooker SJ, Hollingsworth TD. How effective is school-based deworming for the community-wide control of soil-transmitted helminths? PLoS Negl Trop Dis. 2013;7(2), e2027.PubMed CentralView ArticlePubMedGoogle Scholar
- Zwang J, Olliaro PL. Clinical efficacy and tolerability of praziquantel for intestinal and urinary schistosomiasis-a meta-analysis of comparative and non-comparative clinical trials. PLoS Negl Trop Dis. 2014;8(11), e3286.PubMed CentralView ArticlePubMedGoogle Scholar
- de Vlas SJ, Nagelkerke NJ, Habbema JD, van Oortmarssen GJ. Statistical models for estimating prevalence and incidence of parasitic diseases. Stat Methods Med Res. 1993;2(1):3–21.View ArticlePubMedGoogle Scholar