- Open Access
Feasibility of controlling hookworm infection through preventive chemotherapy: a simulation study using the individual-based WORMSIM modelling framework
Parasites & Vectorsvolume 8, Article number: 541 (2015)
Globally, hookworms infect 440 million people in developing countries. Especially children and women of childbearing age are at risk of developing anaemia as a result of infection. To control hookworm infection and disease (i.e. reduce the prevalence of medium and heavy infection to <1 %), the World Health Organization has set the target to provide annual or semi-annual preventive chemotherapy (PC) with albendazole (ALB) or mebendazole (MEB) to at least 75 % of all children and women of childbearing age in endemic areas by 2020. Here, we predict the feasibility of achieving <1 % prevalence of medium and heavy infection, based on simulations with an individual-based model.
We developed WORMSIM, a new generalized individual-based modelling framework for transmission and control of helminths, and quantified it for hookworm transmission based on published data. We simulated the impact of standard and more intense PC strategies on trends in hookworm infection, and explored the potential additional impact of interventions that improve access to water, sanitation, and hygiene (WASH). The individual-based framework allowed us to take account of inter-individual heterogeneities in exposure and contribution to transmission of infection, as well as in participation in successive PC rounds.
We predict that in low and medium endemic areas, current PC strategies (including targeting of WCBA) will achieve control of hookworm infection (i.e. the parasitological target) within 2 years. In highly endemic areas, control can be achieved with semi-annual PC with ALB at 90 % coverage, combined with interventions that reduce host contributions to the environmental reservoir of infection by 50 %. More intense PC strategies (high frequency and coverage) can help speed up control of hookworm infection, and may be necessary in some extremely highly endemic settings, but are not a panacea against systematic non-participation to PC.
Control of hookworm infection by 2020 is feasible with current PC strategies (including targeting of WCBA). In highly endemic areas, PC should be combined with health education and/or WASH interventions.
Globally, over 1 billion people in developing countries are infected with soil-transmitted helminths (STH), of which about 440 million people are infected with at least one type of hookworm (Necator americanus or Ancylostoma duodenale) . Adult hookworms attach themselves to the intestinal mucosae to feed on host blood, causing leakage of intestinal blood and thus contributing to the development of iron-deficiency anaemia. The risk of anaemia is highest in heavily infected individuals, as well as children and women of childbearing age (WCBA), given their naturally low iron reserves [2, 3]. To control the global disease burden of hookworm and other STH infections, the World Health Organization (WHO) has set the operational target to provide regular preventive chemotherapy (PC) to at least 75 % of the population at highest risk for hookworm morbidity by 2020, i.e. pre-school (preSAC) and school-age children (SAC), and WCBA. The associated parasitological goal is to reduce the prevalence of medium and heavy infection (≥2000 eggs per gram (epg) faeces) to levels under 1 % among preSAC, SAC, and WCBA by 2020 . PC targeting preSAC and SAC is typically implemented at the level of schools, its frequency (annual or semi-annual) depending on pre-control STH infection levels . PC targeting WCBA is not yet widely implemented, but planned to be rolled out over the coming years. Mass drug administration (MDA; i.e. PC targeting the whole population) against STH is not officially implemented, but practically taking place in areas where The Global Program for Elimination of Lymphatic Filariasis  is fighting lymphatic filariasis by means of MDA using a combination of albendazole and diethylcarbamazine or ivermectin. Similarly, mass drug administration with ivermectin by the African Programme for Onchocerciasis Control will most likely have already had a significant impact on the STH burden .
Hookworm epidemiology differs from that of other STH, bringing with it a particular challenge for control. The intensity and prevalence of hookworm infection typically rise during childhood and reach a plateau in adult persons, whereas in ascariasis and trichuriasis, infection levels typically peak in children and then decline with age [2, 3]. An explanation for this difference is that ascariasis and trichuriasis are transmitted through ingestion of worm eggs, a mechanism that is strongly related to hygienic practices, which are typically poorer in children than in adults. In contrast, hookworm infection is acquired through larval skin penetration, a mechanism that is related to footwear practices and movement patterns, which are relatively stable over different ages. Further, although not practiced in all endemic regions, the use of human excrements as fertiliser (night soil) provides an additional mechanism by which adults are exposed to hookworm infection. Given that current PC programs mostly target preSAC and SAC, while adult hosts harbour most hookworms and therefore probably contribute most to transmission, the feasibility of controlling hookworm infection by 2020 with current PC strategies can be questioned .
Here, we predict the feasibility of achieving <1 % prevalence of medium and heavy hookworm infection by 2020 with the currently recommended PC strategies (annual or semi-annual PC at 75 % coverage, targeting preSAC, SAC, and WCBA). We further predict the impact of more intense PC strategies (higher frequency and coverage) and targeting of the entire population of age two and above. We performed simulations in WORMSIM, a newly developed generalised modelling framework for transmission and control of helminth infections. The individual-based nature of WORMSIM allows us to take account of important sources of heterogeneity at the individual level, notably exposure and contribution to transmission of infection, as well as participation in successive PC rounds. We evaluated the feasibility of control for various scenarios pertaining to pre-control infection levels and patterns in individual participation in PC, including varying levels of systematic (non-)participation of a subgroup of individuals.
General outline of the WORMSIM modelling framework
WORMSIM is a generalised individual-based modelling framework for transmission and control of helminth infections in humans (Fig. 1), and is based on earlier individual-based models for onchocerciasis, schistosomiasis, and lymphatic filariasis [9–11]. Here we describe the general outline of the framework, partly based on an earlier description of ONCHOSIM . Additional file 1 provides more technical details and the mathematical formulae underlying the model. Additional file 2 contains a zip-compressed version of the WORMSIM modelling framework, including example input files.
WORMSIM simulates the life histories of individual humans and individual worms within humans. Simulated humans are born and die in a stochastic fashion, based on user-specified fertility and life tables. The size of the simulated human population depends on the specified fertility and survival rates, as well as an optionally defined maximum population size. When the maximum population size is exceeded, a random fraction of the population is removed from the simulation. Other than that, WORMSIM simulates a closed population, meaning that there is no migration of humans into or out of the population.
Acquisition of infection and worm reproduction
Simulated humans are exposed and contribute to a central reservoir of infection, which can be taken to represent a cloud of vectors that transmit infection (e.g. as for onchocerciasis), or an environmental reservoir of infective material (more details below). The rate at which each simulated individual is exposed to the reservoir of infection may depend on the season of the year, age, sex, and random personal factors defined in terms of relative individual exposure (following some distribution with mean one, for which we here use a gamma distribution with equal shape and rate α Exi ). The overall exposure rate in the simulation (i.e. for a person with relative exposure 1.0) is defined in terms of parameter ζ. When infective material from the reservoir is successfully transmitted to an individual, it may develop into a male or female adult worm. Once a female worm lives beyond a user-specified prepatent period, she starts producing eggs or larvae, as long as she is regularly inseminated by a patent male worm present in the same host. The egg or larval production can be specified to depend on the age of the female worm. The lifespan of adult worms is simulated in a stochastic fashion, given some user-defined, positively bounded continuous distribution.
Transmission of infection to environmental reservoir or vector
Humans containing reproductive adult worms contribute infective material (larvae or eggs) to the central reservoir of infection. An individual’s contribution rate may depend on the season of the year, age, sex, and random personal factors. When the reservoir is taken to represent a cloud of vectors that transmit larvae (as for onchocerciasis), infective material passes through the reservoir instantaneously (a reasonable assumption when passage through the vector is short relative to the discrete simulation time steps of 1 month), such that the force of infection acting on the human population is always proportional to the amount of infective material contributed to the central reservoir by the human population. When the reservoir is taken to represent an environmental reservoir of infection (as for soil-transmitted helminths), infective material is considered to accumulate and decay within the reservoir, given the total contribution of the human population and some exponential decay rate for infective material in the reservoir. In this case, the force of infection acting on the human population is proportional to the amount of infective material currently present in the environmental reservoir.
Density dependence in transmission
Density dependence in transmission can be specified at several points in the transmission cycle: uptake of infectious material by the central reservoir of infection (e.g. due to limited vectorial capacity), worm fecundity (e.g. due to host immune response and/or competition for nutrients), and worm establishment (e.g. worms already present in the host may trigger partial immunity).
Drug treatment can be specified to temporarily and/or permanently reduce the reproductive capacity of female worms, and to kill adult worms and/or infective material (larvae or eggs) present in the host. The probability that an individual participates in a PC programme is determined by the overall PC coverage level, and the relative probability of participating, given an individual’s age, sex (taking into account treatment eligibility), and/or a lifelong compliance factor. Individual participation to PC is assumed to be either random (given age and sex), fully systematic (given the lifelong compliance factor), or a mix of random and systematic participation (see Additional file 1 for details). In addition, the user can specify that treatment fails in a random fraction of people (e.g. due to malabsorption).
WORMSIM quantification for hookworm transmission
We set the average lifespan of hookworms in the human host to 3 years [3, 13–15], and assumed that the worm mortality rate increases linearly with worm age (i.e. worm lifespan follows a Weibull distribution with mean 3.0 and shape 2.0). Based on literature, we assumed that hookworms can reproduce after a fixed prepatent period of 7 weeks [2, 3, 13, 16]. Fecundity-related parameters were set such that egg production was independent of female worm age, and female worms could produce eggs as long as at least one male worm was also present in the host. The total egg output of all female worms in a host was assumed to be negatively density-dependent on the total number of female worms x. This density dependence was defined in terms of the hyperbolic saturating function αx/(1 + αx/β), where α is the average egg production per female worm in absence of density dependence (α = 200 epg ), and β is the average maximum egg output for a host (see section Endemicity scenarios below for details about quantification of β).
Host suitability for infection
To capture inter-individual variation in host suitability for worm infection (e.g. due to genetic factors [18–20], nutrition status, and/or immunocompetence), each individual’s saturation level for total egg output β i was assumed to be a random lifelong trait with inter-individual variation proportional to a gamma distribution with mean 1.0. As there simply is no field data on this particular host characteristic, we first assumed that variation in relative host susceptibility was low with 95 %-CI: 0.74–1.30 (gamma distribution with shape and rate equal to 50), such that predictions were very similar to the assumption of “no variation in susceptibility” (i.e. relative susceptibility = 1.0, which has been typically employed so far). Next, we repeated the whole analysis with an arbitrary, much more extreme assumption about high variation in host susceptibility (95 %-CI: 0.12–2.29; shape and rate equal to two, such that the value of zero relative susceptibility still has zero density).
Larvae in the environment
The average lifespan of larvae in the environmental reservoir was set to 2 weeks, assuming an exponential distribution [2, 16, 21]. Because the survival of larvae may vary geographically with environmental conditions (average humidity and temperature), we alternatively assumed that average larval lifespan is 4 weeks (implying higher reinfection rates between PC rounds). Host contribution and exposure to the environmental reservoir were assumed to be perennial, with age-patterns increasing linearly from zero to one between ages zero and ten, such that we achieve the typical age pattern of infection levels increasing with age until they reach a plateau at age 20 (Fig. 2) . In addition, we assumed that individual exposure and contribution rates vary randomly due to personal factors (for quantification, see section Endemicity scenarios). Contribution and exposure were assumed to be perfectly correlated for each individual.
Effect of preventive chemotherapy
We assumed that treatment with albendazole (ALB) or mebendazole (MEB) kills a given fraction of prepatent and adult worms, and has no effect on worm fecundity. We further assumed that the proportion of worms killed by ALB or MEB is equal to observed reductions in mean egg counts (95 and 80 % respectively) . This is a reasonable assumption, as in the cited study reductions in hookworm egg counts were not correlated with pre-treatment egg counts. Drug efficacy was assumed to be equal for all host and worm ages.
Simulation output on infection levels was defined in terms of prevalence of no, light, medium, and heavy infection (cut-offs: 1, 2000, 4000 epg), based on single Kato-Katz slides of 41.7 mg. Kato-Katz slides are wet mounts of faecal samples, which are systematically examined under a microscope to count worm eggs . The sensitivity of Kato-Katz slides to detect infection increases with the number of eggs in the sample. To simulate this, we assumed that sampling error in Kato-katz slides follows a negative binomial distribution with mean egg count as predicted by WORMSIM for a given person, and aggregation parameter k Kato-Katz = 0.40, based on an analysis of repeated slides from 2083 Ugandan individuals (see Additional file 3). We assumed perfect specificity of Kato-Katz testing (i.e. as if performed by a trained and experienced laboratory technician). Simulation output on prevalence of infection was stratified for the following sub-populations: infants (age <2), pre-SAC (age 2–5), SAC (age 5–15), WCBA (females of age 15–45), and other (females of age 45 and above, and males of age 15 and above).
Comparison to field data
First, we compared WORMSIM predictions for the impact of PC to field trial data on trends in hookworm burden in a population of Vietnamese WCBA who were offered 4-monthly treatment with albendazole for 1 year, and 6-monthly from then onwards [24, 25]. For this setting, we assumed that the average saturation level for egg output was 1500 epg (between 1113 and 1943 for 95 % of individuals) or 2000 epg (1484–2591), and that the coverage of mass treatment targeting WCBA was as reported during a single cross-sectional survey (85.8 %) and remained stable for the entire 54-month study period. Because we had no exact information on the timing and coverage of PC targeting preSAC and SAC in the study area, we assumed that these were treated at the same time as the WCBA, and at equal coverage.
For the purpose of predicting the impact of PC, we defined three endemicity levels (high, medium, low) representative of field conditions (Fig. 3), based on an analysis of literature data [22, 26–31] (previously collated and described elsewhere ). Each endemicity level was defined in terms of mean egg counts and a negative binomial aggregation k, from which we derived target values for the distribution of intensity of infection (none, light, medium, heavy) to reproduce in WORMSIM (Table 1). The endemicity scenarios were arbitrarily chosen such that each scenario had a pre-control prevalence of medium and heavy infection (≥2000 epg) in 10–14 year old children of at least 1 % (such that preventive chemotherapy is still indicated), and such that together, the endemicity scenarios spanned the range of infection levels observed in literature (Fig. 3). For technical details on how these endemicity scenarios were arrived at, see Additional file 3.
The average saturation level for host egg output was set to β = 1500 epg such that WORMSIM could reproduce the distribution of intensity of infection for the highly endemic scenario, and such that it could also still simulate sustained transmission for the low endemicity scenario. We also investigated the alternative assumption that average saturation level in egg output is higher (β = 2000 epg, i.e. weaker density dependence in transmission). Similarly, for low and medium endemic areas, we investigated the alternative assumption that average saturation level is lower (β = 1000 epg, i.e. stronger density dependence). Next, we performed a grid search to quantify the remaining free transmission parameters ζ and α Exi in WORMSIM, so as to reproduce the distribution of light, medium, and heavy infection as expected for each endemicity level (see in Additional file 3: Table A3–2).
Preventive chemotherapy programmes
Next, for each endemicity scenario (high, medium, low), we performed 250 repeated simulations to predict the average impact of PC on hookworm transmission in a closed population of about 400 individuals, given the estimated efficacy of a single drug treatment. Simulations that resulted in interruption of transmission before start of control were dropped from the analyses (this occurred in about 40 % of the simulations for the low endemic scenario with β = 2000). We varied assumptions about PC frequency (annual, semi-annual, 4-monthly, or quarterly), target population (preSAC and SAC; preSAC, SAC and WCBA; or total population of age two and higher), coverage of target population (75 %, the WHO operational nation-level target, or 90 %, the level of coverage typically achieved in individual schools), and patterns in participation (random, mixed, or systematic, see in Additional file 1: Figure A1-2 for illustration) as also previously used in mathematical modelling of onchocerciasis . We assumed that PC coverage is stable over time, once PC is implemented, even though national coverage rates for preSAC and SAC published by WHO increase over time . However, these national figures are based on a mix of local coverage rates of mostly zeroes (no PC) and high rates (75–95 % wherever PC is successfully implemented), and therefore mainly reflect geographical scaling up of PC. Because transmission takes place locally, we chose to simulate stable coverage only in our main analysis. Only to compare our model predictions with those by Truscott et al. , who assume that local PC coverage levels do follow the national average trend, we also provide predictions for this scenario (i.e. assuming a linear increase in coverage of preSAC and SAC, starting at 0 % in 2002, and increasing up to 75 % in 2020, and stable from then onwards).
Health education and WASH
To explore the potential effects of health education and improved access to water, sanitation, and hygiene (WASH), we simulated the impact of a 50 % reduction in contribution of all individuals to the environmental reservoir (WORMSIM does not yet support an effect of WASH on exposure, e.g. as a result of improved flooring in dwellings).
Comparison of prediction to field data
WORMSIM-predicted trends for light, medium, and heavy infection were in good agreement with field data on WCBA from Vietnam at 3 and 12 months after the first treatment rounds (Fig. 4). However, at 30 months after start of the PC programme, the reported prevalence of heavy infection was higher than predicted by WORMSIM. At 54 months, the reported prevalence of light and medium infection was lower than predicted by WORMSIM. These differences may be explained by sampling issues at 30 months and changes in behaviour over time (see Discussion).
Predicted impact of preventive chemotherapy in low and medium endemic areas
For the low and medium endemic scenarios, current PC strategies targeting only children (annual or semi-annual PC at 75 % coverage) are predicted to achieve control of hookworm infection (prevalence of medium and heavy infection <1 %) in preSAC and SAC within a few years (Fig. 5, first two columns). However, PC specifically also targeting WCBA (third and fourth column) is required to also achieve control in WCBA (blue line) within the same time span (low endemic scenario), or at all (medium endemic scenario). For women over 45 years of age and men over 15 years of age (pink line; not targeted by PC) PC targeting preSAC, SAC, and WCBA provides some benefit through indirect transmission effects, but this decreases with higher pre-control infection levels. These findings are independent of the choice of drug: ALB (Fig. 5) or MEB (Additional file 4).
Predicted impact of preventive chemotherapy in highly endemic areas
For the highly endemic scenario, semi-annual PC at 75 % coverage is not sufficient to achieve control (Fig. 5). However, we predict that control in preSAC, SAC, and WCBA can be achieved in highly endemic areas by implementing PC with ALB either quarterly at 75 % coverage, 4-monthly at 90 % coverage, or semi-annually targeting the whole population (under our standard assumption of density-dependent fecundity with β = 1500; Fig. 6). In case PC with MEB is implemented, more intensive strategies are required (90 % coverage and/or higher frequency; Additional file 1). Alternatively, control of hookworm infection in highly endemic settings may also be achieved by a combination of the standard strategy of semi-annual PC at 90 % coverage and interventions that reduce the contribution of infectious material by each host by 50 % (Fig. 7). However, under the alternative assumption of lower density dependence in worm fecundity (β = 2000), we predict that control can even be achieved with semi-annual PC with ALB alone, when implemented at 90 % coverage and targeting preSAC, SAC, and WCBA.
Impact of systematic (non-)participation
Figure 8 illustrates how systematic (non-)participation of individuals to PC with albendazole can dramatically decrease the impact of PC, regardless of treatment frequency. However, as long as all individuals have a chance to be treated at some point (mixed participation), the impact of PC was predicted to be very similar to that in the scenario of completely random participation. This finding was the same for PC with MEB (Additional file 4).
All aforementioned findings with regard to achieving control were robust to alternative assumptions about the level of density dependence in transmission, inter-individual variation in host suitability for infection, and the lifespan of larvae in the environmental reservoir (Additional file 4).
Impact of PC during scaling up
Figure 9 illustrates the impact of annual PC targeting children, with coverage scaling up from 0 to 75 % between 2002 and 2020, based on trends in national coverage as reported by WHO.
In low and medium endemic settings, achieving control of hookworm infection in preSAC and SAC (i.e. prevalence of medium and heavy infection <1 %) can most likely be achieved with current PC strategies (annual or semi-annual PC at 75 % or 90 % coverage) within as little as 1 to 2 years. Although adults benefit indirectly from PC programs targeting children, this effect is not enough to control of hookworm infection in WCBA in a timely fashion; the latter is best achieved by also actively targeting WCBA. This will be especially important in Asia, where in some localities infection levels keep on rising with age into late adulthood. For low and medium endemic settings, control can be achieved equally well with MEB and ALB.
For highly endemic areas, we predict that control of hookworm infection is feasible by means of semi-annual PC with ALB at 90 % coverage, targeting preSAC, SAC, and WCBA, combined with concomitant interventions that reduce host contributions to the environmental reservoir of infection by, say, about half (e.g. health education or WASH), if implemented successfully. Furthermore, if infection levels in highly endemic areas are (partly) driven by lower density-dependent worm fecundity (e.g. resulting from lower host immunocompetence, such that people actually do not have excessively many worms, but just very productive worms), control may be even achieved with semi-annual PC alone, implemented at 90 % coverage. More intense PC strategies (higher frequency and coverage) can help speed up control of hookworm infection, and may be necessary in some extremely highly endemic settings. Control programs should however keep in mind that more intense PC is not a panacea against systematic non-compliance to PC, which may occur among adults in particular (who are not necessarily treated centrally at school, like children), as observed in mass drug administration programs against lymphatic filariasis . Further, in settings where transmission is seasonal (rather than perennial, as assumed for our predictions), timing of PC right before the transmission season (warm and humid weather) may provide an additional means of maximizing the impact of control efforts.
WORMSIM could reasonably well reproduce trends in light, medium, and heavy infection in Vietnamese WCBA during 1 year of 4-monthly treatment with albendazole. The high prevalence of heavy infection reported at 30 months after the start of the PC programme is most likely caused by issues related to sampling (only 187 WCBA were tested) and suboptimal PC coverage preceding the 30-month survey. Further, our long-term predictions for prevalence of light infection (at 30 and 54 months) were relatively pessimistic compared to reported prevalence figures. The most likely explanation is that the health education component of the field trial helped improve human behaviour related to contribution and/or exposure to the environmental reservoir of infection (nurses and health workers were actively trained and provided with educational material to use during their repeated contacts with the women). Another explanation may be that albendazole not only kills worms, but also somehow cumulatively impairs reproduction of surviving worms, as has been proposed to be an effect of ivermectin on Onchocerca volvulus worms .
In our exploratory simulations for the potential effects of health education and WASH, we assumed that a hypothetical intervention reduces the host contribution of infective material to the environmental reservoir by 50 %, and equally so for all hosts (WORMSIM does not yet support interventions that impact host exposure). In reality, the effect of health education and WASH interventions most likely varies between individuals due to behavioural factors, and may therefore be relatively lower at the population level than predicted here. However, while we only considered an impact on host contribution to the environmental reservoir of infection, the impact of health education and WASH on exposure to hookworm infection is probably important as well. For instance, natural dwelling floor types (e.g., earth, sand, dung, or mud) greatly contribute to host exposure to infection, and have been estimated to be responsible for as much as 86 % of all hookworm infections . Therefore, we expect that the overall impact of successfully implemented health education and WASH interventions may be even larger than predicted here. Still, some may argue that a 50% reduction is currently too optimistic, given that so far, large-scale community-based WASH trials in India have shown little impact on latrine use and child health [38, 39]. However, a recent meta-analysis did show strong associations between availability of sanitation and STH infections (odds ratio around 0.5), and use of sanitation and hookworm infection (odds ratio 0.6) . Further, a school-based WASH RCT in Kenya did show a reduction in reinfection rate of about 50 % . Also, the Magic Glasses project, a video-based hygiene education program implemented in Chinese schools, was able to reduce incidence of STH infection in children (after albendazole treatment) by 50 % . It may be feasible to achieve similar impacts with WASH in community settings, if implemented successfully, which will most likely depend more on behavioural factors than bricks and mortar .
The literature data underlying our predictions are probably most representative of Necator americanus, the most prevalent human hookworm species globally, but whose eggs are morphologically indistinguishable from Ancylostoma duodenale eggs. Because the egg production rate of A. duodenale is believed to be two to three times higher than that of N. americanus , our model predictions in terms of absolute egg counts may not apply directly to localities where A. duodenale is the prominent hookworm species. Still, assuming that density dependence in transmission of A. duodenale and N. americanus is comparable (in terms of the ratio of eggs per worms and maximum host egg output), our model predictions also apply to A. duodenale in a qualitative sense.
Defining the endemicity scenarios for this study, we assumed that in highly endemic areas, heterogeneity in individual exposure and contribution to the environmental reservoir is lower. This clearly resulted from our analysis of literature data on distribution of intensity of infection (Fig. 3), and seems plausible, as in such areas everybody can be thought to walk barefoot and/or defecate in the same area. It may also be that in low endemic areas, density dependence in worm fecundity is stronger because people are less susceptible to heavy infections, e.g. due to higher immunocompetence. However, the finding that at low and medium intensity levels of hookworm infection, faecal egg reduction rates of albendazole and mebendazole are independent of pre-treatment egg counts (in contrast to ascariasis and trichuriasis, for which negative correlations were observed), suggests that at these intensities of infection, density dependence only plays a marginal role . Furthermore, our sensitivity analyses show that our predictions for achieving control of hookworm infection in low and medium endemic areas are robust to alternative assumptions about the mean level of and variation in host suitability for infection.
As previously suggested , our findings confirm that school-based deworming may not always be enough to control hookworm infection in population at high risk for morbidity, and that women of childbearing age should also be targeted with PC. Still, the relative importance of human subpopulations of different ages in hookworm transmission is not well known. Anderson et al. explored different assumptions about the relative contribution and exposure of children and adults to the environmental reservoir . Their findings suggest that if children contaminate the environment relatively more often than adults, the impact of PC targeting only children would be relatively larger and also significantly impact infection levels in adults. However, we argue that this is an unlikely scenario, as children’s contributions would have to be higher than adults’ because of their open defecation practices. However, these practices are most likely also correlated with exposure to infection, while typically, observed infection levels are lower in children. In the current study, we therefore assumed that host exposure and contribution to the environmental reservoir of infection are proportional to each other, and that both increase with age (up to the age of 10 years) as a result of open defecation practices. If after all, children do indeed contaminate the environment relatively more often (but are not exposed relatively more often), PC targeting children will have a larger impact on transmission of hookworm infection than predicted here.
A set of previous STH modelling studies focussing on interruption of STH transmission highlight several issues that may also be important for STH control [13–15]. First of all, to reduce hookworm transmission, high PC coverage of adults is more important than for ascariasis transmission due to different age-patterns in infection levels [13, 14]. This is supported by our conclusion that PC should also target women of childbearing age to control hookworm. Second, helminth mating processes and the dynamics of the environmental reservoir play an important role in elimination of STH . In the current study we account for both by explicitly simulating mating events between male and female worms and decay of larvae in the environment. Mating processes are probably less important for control than elimination of infection as in a control situation, relatively many infected individuals will still harbour multiple worms. However, the lifespan of infective material in the environment may weigh into the speed at which control or elimination is achieved. In our simulations, an average larval lifespan of 2 or 4 weeks did not matter much for the speed at which hookworm control is achieved. However, for ascariasis and trichuriasis, the average lifespan of eggs in the environment (which is in the order of months) will probably play a more important role in control and elimination of infection. We will further examine both aspects (PC coverage of age-groups and egg lifespan) in future modelling studies with WORMSIM.
In a similar modelling study, Truscott et al. take a more global perspective on the impact of current PC strategies on ascariasis, trichuriasis, and hookworm, using national trends on PC coverage . Although Truscott et al. predict trends in average worm burdens, their predictions are qualitatively very similar to ours, in that for hookworm, PC targeting children will not benefit adults very much. An important difference in study design however is that Truscott et al. used data national trends in PC coverage, while in our main analysis, we assumed PC coverage is stable over time. We believe our approach is more representative of local situations where STH control is actually taking place and where PC coverage rates are relatively high compared to the national average. When we employed axxproimately the same assumptions about scaling up of PC as Truscott et al. (PC coverage of preSAC and SAC increases linearly from 0 to 75 % between 2002 and 2020), we get very similar results: the impact of PC only becomes really noticeable when coverage levels approach 50 % (between 2010 and 2015). Still, we believe that these predictions do not realistically represent national or global trends in infection levels, but only a local setting (i.e. community-level) in which PC coverage is scaled up.
We predict that control of hookworm infection in low and medium endemic areas by 2020 is feasible with current PC strategies, which we consider to include targeting of WCBA. Control of hookworm infection in highly endemic areas may be achieved with a minimum of semi-annual PC with ALB at 90 % coverage, combined with health education and/or WASH interventions. More intense PC strategies (high frequency and coverage) may help speed up control of hookworm infection, and may be necessary in some extremely highly endemic settings.
Eggs per gram
Mass drug administration
Negative binomial (distribution)
Women of childbearing age
World Health Organization
Pullan RL, Smith JL, Jasrasaria R, Brooker SJ. Global numbers of infection and disease burden of soil transmitted helminth infections in 2010. Parasit Vectors. 2014;7:37.
Hotez PJ, Brooker S, Bethony JM, Bottazzi ME, Loukas A, Xiao S. Hookworm infection. N Engl J Med. 2004;351:799–807.
Bethony J, Brooker S, Albonico M, Geiger SM, Loukas A, Diemert D, et al. Soil-transmitted helminth infections: ascariasis, trichuriasis, and hookworm. Lancet. 2006;367:1521–32.
World Health Organization. Sustaining the drive to overcome the global impact of neglected tropical diseases. Geneva: World Health Organization. 2013.
Gabrielli A-F, Montresor A, Chitsulo L, Engels D, Savioli L. Preventive chemotherapy in human helminthiasis: theoretical and operational aspects. Trans R Soc Trop Med Hyg. 2011;105:683–93.
Ottesen EA, Hooper PJ, Bradley M, Biswas G. The global programme to eliminate lymphatic filariasis: health impact after 8 years. PLoS Negl Trop Dis. 2008;2:e317.
Krotneva SP, Coffeng LE, Noma M, Zouré HGM, Bakoné L, Amazigo UV, et al. African Program for Onchocerciasis Control 1995-2010: impact of annual ivermectin mass treatment on off-target infectious diseases. PLoS Negl Trop Dis. 2015;9:e0004051.
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:e2027.
Plaisier AP, van Oortmarssen GJ, Habbema JD, Remme J, Alley ES. ONCHOSIM: a model and computer simulation program for the transmission and control of onchocerciasis. Comput Methods Programs Biomed. 1990;31:43–56.
De Vlas SJ, Van Oortmarssen GJ, Gryseels B, Polderman AM, Plaisier AP, Habbema JD. SCHISTOSIM: a microsimulation model for the epidemiology and control of schistosomiasis. Am J Trop Med Hyg. 1996;55:170–5.
Plaisier AP, Subramanian S, Das PK, Souza W, Lapa T, Furtado AF, et al. The LYMFASIM simulation program for modeling lymphatic filariasis and its control. Method Inf Med. 1998;37:97–108.
Coffeng LE, Stolk WA, Hoerauf A, Habbema D, Bakker R, Hopkins AD, et al. Elimination of African onchocerciasis: modeling the impact of increasing the frequency of ivermectin mass treatment. PLoS One. 2014;9:e115886.
Anderson RM, Truscott J, Hollingsworth TD. The coverage and frequency of mass drug administration required to eliminate persistent transmission of soil-transmitted helminths. Philos Trans R Soc L B Biol Sci. 2014;369:20130435.
Truscott JE, Hollingsworth TD, Brooker SJ, Anderson RM. Can chemotherapy alone eliminate the transmission of soil transmitted helminths? Parasit Vectors. 2014;7:266.
Truscott J, Hollingsworth TD, Anderson RM. Modeling the interruption of the transmission of soil-transmitted helminths by repeated mass chemotherapy of school-age children. PLoS Negl Trop Dis. 2014;8:e3323.
Brooker S, Bethony J, Hotez PJ. Human Hookworm Infection in the 21st Century. Adv Parasitol. 2004;58:197–288.
Anderson RM, Schad GA. Hookworm burdens and faecal egg counts: an analysis of the biological basis of variation. Trans R Soc Trop Med Hyg. 1985;79:812–25.
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:e713.
Williams-Blangero S, Blangero J, Bradley M. Quantitative genetic analysis of susceptibility to hookworm infection in a population from rural Zimbabwe. Hum Biol. 1997;69:201–8.
Quinnell RJ, Pullan RL, Breitling LP, Geiger SM, Cundill B, Correa‐Oliveira R, et al. Genetic and household determinants of predisposition to human hookworm infection in a Brazilian Community. J Infect Dis. 2010;202:954–61.
Augustine DL. Investigations on the control of hookworm disease. XVI. Length of life of hookworm larvae from the stools of different individuals. Am J Epidemiol. 1923;3:127–36.
Levecke B, Montresor A, Albonico M, Ame SM, Behnke JM, Bethony JM, et al. Assessment of anthelmintic efficacy of mebendazole in school children in six countries where soil-transmitted helminths are endemic. PLoS Negl Trop Dis. 2014;8:e3204.
World Health Organization. Basic laboratory methods in medical parasitology. Geneva; 1991.
Casey GJ, Phuc TQ, Macgregor L, Montresor A, Mihrshahi S, Thach TD, et al. A free weekly iron-folic acid supplementation and regular deworming program is associated with improved hemoglobin and iron status indicators in Vietnamese women. BMC Public Health. 2009;9:261.
Casey GJ, Montresor A, Cavalli-Sforza LT, Thu H, Phu LB, Tinh TT, et al. Elimination of iron deficiency anemia and soil transmitted helminth infection: evidence from a fifty-four month iron-folic acid and de-worming program. PLoS Negl Trop Dis. 2013;7:e2146.
Montresor A, Urbani C, Camara B, Bha AB, Albonico M, Savioli L. Preliminary survey of a school health program implementation in Guinea. Med Trop. 1997;57:294–8.
Albonico M, Bickle Q, Ramsan M, Montresor A, Savioli L, Taylor M. Efficacy of mebendazole and levamisole alone or in combination against intestinal nematode infections after repeated targeted mebendazole treatment in Zanzibar. Bull World Heal Organ. 2003;81:343–52.
Gabrielli AF, Ramsan M, Naumann C, Tsogzolmaa D, Bojang B, Khoshal MH, et al. Soil-transmitted helminths and haemoglobin status among Afghan children in World Food Programme assisted schools. J Helminthol. 2005;79:381–4.
Phommasack B, Saklokham K, Chanthavisouk C, Nakhonesid-Fish V, Strandgaard H, Montresor A, et al. Coverage and costs of a school deworming programme in 2007 targeting all primary schools in Lao PDR. Trans R Soc Trop Med Hyg. 2008;102:1201–6.
Pasricha S-R, Caruana SR, Phuc TQ, Casey GJ, Jolley D, Kingsland S, et al. Anemia, iron deficiency, meat consumption, and hookworm infection in women of reproductive age in northwest Vietnam. Am J Trop Med Hyg. 2008;78:375–81.
Vercruysse J, Behnke JM, Albonico M, Ame SM, Angebault C, Bethony JM, et al. Assessment of the anthelmintic efficacy of albendazole in school children in seven countries where soil-transmitted helminths are endemic. PLoS Negl Trop Dis. 2011;5:e948.
Montresor A, À Porta N, Albonico M, Gabrielli AF, Jankovic D. Soil-transmitted helminthiasis: the relationship between prevalence and classes of intensity of infection. Trans R Soc Trop Med Hyg. 2015;109:262–7.
World Health Organization (n.d.) PCT Databank. www.who.int/neglected_diseases/preventive_chemotherapy. Accessed 6 Aug 2015.
Truscott JE, Turner HC, Anderson RM. What impact will the achievement of the current World Health Organisation targets for anthelmintic treatment coverage in children have on the intensity of soil-transmitted helminth infections? Parasit Vectors 2015. In press.
Krentel A, Fischer PU, Weil GJ. A review of factors that influence individual compliance with mass drug administration for elimination of lymphatic filariasis. PLoS Negl Trop Dis. 2013;7:e2447.
Plaisier AP, Alley ES, Boatin BA, Van Oortmarssen GJ, Remme H, De Vlas SJ, et al. Irreversible effects of ivermectin on adult parasites in onchocerciasis patients in the Onchocerciasis Control Programme in West Africa. J Infect Dis. 1995;172:204–10.
Soares Magalhaes RJ, Barnett AG, Clements ACA. Geographical analysis of the role of water supply and sanitation in the risk of helminth infections of children in West Africa. Proc Natl Acad Sci U S A. 2011;108:20084–9.
Clasen T, Boisson S, Routray P, Torondel B, Bell M, Cumming O, et al. Effectiveness of a rural sanitation programme on diarrhoea, soil-transmitted helminth infection, and child malnutrition in Odisha, India: a cluster-randomised trial. Lancet Glob Heal. 2014;2:e645–53.
Patil SR, Arnold BF, Salvatore AL, Briceno B, Ganguly S, Colford JM, et al. The effect of India’s total sanitation campaign on defecation behaviors and child health in rural Madhya Pradesh: a cluster randomized controlled trial. PLoS Med. 2014;11:e1001709.
Ziegelbauer K, Speich B, Mäusezahl D, Bos R, Keiser J, Utzinger J. Effect of sanitation on soil-transmitted helminth infection: systematic review and meta-analysis. PLoS Med. 2012;9:e1001162.
Freeman MC, Clasen T, Brooker SJ, Akoko DO, Rheingans R. The impact of a school-based hygiene, water quality and intervention on soil-transmitted helminth reinfection: a cluster-randomized trial. Am J Trop Med Hyg. 2013;89:875–83.
Bieri FA, Gray DJ, Williams GM, Raso G, Li Y-S, Yuan L, et al. Health-education package to prevent worm infections in Chinese Schoolchildren. N Engl J Med. 2013;368:1603–12.
Brocklehurst C. Scaling up rural sanitation in India. PLoS Med. 2014;11:e1001710.
LEC, RB, and SJdV gratefully acknowledge funding of the NTD Modelling Consortium by the Bill and Melinda Gates Foundation in partnership with the Task Force for Global Health. The views, opinions, assumptions or any other information set out in this article are solely those of the authors. We would like to thank Simon Brooker (Londen School of Hygiene and Tropical Medicine) for sharing Ugandan data on repeated egg counts.
The authors declare that they have no competing interests.
LEC designed the extension of WORMSIM for soil-transmitted helminths, analysed and interpreted the data, performed model simulations, and drafted the manuscript. RB programmed the WORMSIM extension for soil-transmitted helminths, interpreted the results of the analysis, and critically revised the manuscript. AM collated literature data on the association between mean level and distribution of intensity of infection, interpreted the results of the analysis, and critically revised the manuscript. SJdV conceived of the study, interpreted the results of the analysis, and critically revised the manuscript. All authors read and approved the final manuscript.
LEC is a postdoctoral researcher in infectious disease control working on mathematical modelling of transmission and control of soil-transmitted helminths, onchocerciasis, lymphatic filariasis, and visceral leishmaniasis. RB is a biologist and senior scientific programmer of individual-based models for infectious disease transmission. AM is a medical officer and focal point for STH control at World Health Organization. SJdV is an Associate Professor in infectious disease control working on mathematical modelling of transmission and control of both microparasitic and macroparasitic infectious diseases.
This document provides a description of the WORMSIM model structure, the default parameter quantification as used in the simulations, an annotated input and output files, and instructions on how to run WORMSIM. (PDF 11191 kb)
Zip archive containing the JAVA-based WORMSIM modelling framework (version 2.58Ap9), including example input files (XML) for hookworm transmission and detailed documentation of the XML schema for input files. The contents of this zip archive are licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License. To view a copy of this license, visit http://creativecommons.org/licenses/by-nc-nd/4.0/ or send a letter to Creative Commons, PO Box 1866, Mountain View, CA 94042, USA. By opening Additional File X, you agree to the aforementioned license. You are free to use and share (copy and redistribute the material in any medium or format) the material contained within Additional File 4 under the following terms: Attribution — You must give appropriate credit, provide a link to the license, and indicate if changes were made. You may do so in any reasonable manner, but not in any way that suggests the licensor endorses you or your use. NonCommercial — You may not use the material for commercial purposes. NoDerivatives — If you remix, transform, or build upon the material, you may not distribute the modified material. (ZIP 5.69 MB)
Methods and results of analyses for sampling error in Kato-Katz slides and the association between mean and aggregation of egg counts in populations. (PDF 1238 kb)