Skip to main content

New approaches to measuring anthelminthic drug efficacy: parasitological responses of childhood schistosome infections to treatment with praziquantel



By 2020, the global health community aims to control and eliminate human helminthiases, including schistosomiasis in selected African countries, principally by preventive chemotherapy (PCT) through mass drug administration (MDA) of anthelminthics. Quantitative monitoring of anthelminthic responses is crucial for promptly detecting changes in efficacy, potentially indicative of emerging drug resistance. Statistical models offer a powerful means to delineate and compare efficacy among individuals, among groups of individuals and among populations.


We illustrate a variety of statistical frameworks that offer different levels of inference by analysing data from nine previous studies on egg counts collected from African children before and after administration of praziquantel.


We quantify responses to praziquantel as egg reduction rates (ERRs), using different frameworks to estimate ERRs among population strata, as average responses, and within strata, as individual responses. We compare our model-based average ERRs to corresponding model-free estimates, using as reference the World Health Organization (WHO) 90 % threshold of optimal efficacy. We estimate distributions of individual responses and summarize the variation among these responses as the fraction of ERRs falling below the WHO threshold.


Generic models for evaluating responses to anthelminthics deepen our understanding of variation among populations, sub-populations and individuals. We discuss the future application of statistical modelling approaches for monitoring and evaluation of PCT programmes targeting human helminthiases in the context of the WHO 2020 control and elimination goals.


Human helminthiases comprise six of the seven most prevalent neglected tropical diseases (NTDs) [1] that are targeted by the World Health Organization (WHO) for elimination, where feasible, by 2020 [2]. Schistosomiasis infects over 250 million people with an estimated global burden of 3.31 million disability-adjusted life years [3, 4]. The main strategy to control and eliminate human helminthiases is preventative chemotherapy (PCT) by mass drug administration (MDA) using a handful of safe and efficacious anthelminthic drugs [57]. The cornerstone of schistosomiasis control and elimination efforts is praziquantel which is effective in killing adult Schistosoma mansoni, S. haematobium and S. japonicum, the most globally important causes of intestinal, urogenital and Asian intestinal schistosomiasis, respectively.

The reliance on single drugs with essentially no available alternatives makes the long-term effectiveness of the PCT strategy susceptible to the potentially devastating consequences of emerging anthelminthic resistance. Although examples of praziquantel-resistant isolates of S. mansoni are currently scant and virtually non-existent for S. haematobium or S. japonicum [8], there is broad consensus that the efficacy of praziquantel—and other anthelminthics used for MDA—should be monitored to detect atypical responses that may indicate dwindling efficacy, possibly caused by emerging drug resistance [6, 912].

The efficacy of anthelminthics is typically expressed as either a cure rate (CR), or an intensity reduction rate (IRR; Table 1), calculated using data on parasite transmission stages collected before and after treatment. Intensity reduction rates are recommended by the WHO for monitoring the efficacy of both praziquantel and the benzimidazoles, which are widely used in PCT targeting soil-transmitted helminthiasis [13]. In this context, IRRs are calculated using data on helminth egg counts and so are referred to as egg reduction rates (ERRs). The methods of estimating ERRs can be divided into model-free and model-based approaches.

Table 1 Definitions

Model-free approaches calculate ERRs directly from data using simple arithmetic operations, without invocation of distributional (modelling) assumptions. These so-called sample estimates are easy to calculate and straightforward to interpret as population averages [14] and are the most commonly reported estimates of efficacy [15, 16]. Although egg count values are not normally distributed even after log transformation, the WHO recommends arithmetic means over geometric means [13] as the former are more sensitive to outliers and thus more apt to identify suboptimal group responses [14]. However, they are not readily compatible with exploring associations between efficacy and covariates, nor do they permit inference on the underlying distribution of drug responses among individuals.

Model-based approaches are seldom used to estimate efficacy, despite offering a powerful means to conduct multivariate analyses of longitudinal data [17] on egg counts to delineate and compare efficacy among individuals, among groups of individuals and among populations. Furthermore, modelling approaches, particularly in conjunction with Bayesian techniques, can be used to define distributions of responses to anthelminthics among individuals within demographic strata and among populations [16].

Here, we illustrate two distinct modelling approaches by analysing data on schistosome egg counts collected from children infected with S. mansoni or S. haematobium before and up to six weeks after administration of praziquantel. The data are from several past studies on the efficacy of praziquantel from communities in Côte d’Ivoire, Kenya and Uganda, predominantly naïve to MDA, or having received only a few rounds of MDA. We illustrate how two classes of statistical model can be used to (i) identify geographical, demographic and drug regimen covariates associated with ERRs and (ii) explore distributions of individual responses to praziquantel in key population demographics and evaluate the frequency of nominally optimal and sub-optimal responders. We discuss the context in which each modelling approach is most appropriate, depending on the goal of the analysis and the level of inference sought [18]. We also describe how distributions of drug responses among individuals infected with predominantly drug-naïve and maximally susceptible parasites could facilitate identification of sub-optimally or atypically responding individuals, ultimately providing a practical tool for the monitoring and evaluation (M&E) of anthelmintic efficacy during PCT programmes.

Models and methods

Ethics, consent and permissions

This paper reports a secondary analysis of data collected in past studies, all of which obtained the necessary ethical approvals from relevant institutional review boards and local and national ethics committees. All data were collected in accordance with international ethical standards. Data received were completely and irreversibly anonymised. Additional information can be found in the original publications and in Additional file 1: Supplementary Tables, Table S1.

Data selection criteria

We obtained the datasets presented in Olliaro et al. [14] that comprise individual-level data on schistosome egg counts measured before and after administration of praziquantel, collected from 13 studies (Additional file 1: Supplementary Tables, Table S1). The data also comprise individual covariates including: age; sex; dose of praziquantel, and days of follow-up after treatment, and the population covariate, country. Studies were largely non-controlled, non-blinded public health interventions, with one exception [19], which was placebo-controlled and double blinded. Placebo-controlled randomised controlled trials have been rare for helminthiases as it is now generally considered unethical to withhold effective treatment from infected individuals. Following the criteria outlined in Fig. 1, we selected three studies with data on S. haematobium [1921], and six studies with data on S. mansoni [19, 2226] infections in children. Key features of these datasets are summarised in Table 2 with further details given in Additional file 1: Supplementary Tables, Table S1. Raw egg counts and child-specific mean egg counts before and after treatment with praziquantel are depicted in Fig. 2.

Fig. 1
figure 1

Data selection criteria

Table 2 Summary of data included in the analysis
Fig. 2
figure 2

Schistosome egg counts by number of days after treatment with praziquantel. Panel a depicts Schistosoma haematobium egg counts measured by urine filtration. Panel b depicts S. mansoni egg counts measured by Kato-Katz technique. Each data point represents a single count (i.e. not an average of multiple counts). Panels c and d depict the arithmetic mean egg counts per person connected by a line. Treatment with praziquantel occurred following the counts made at day zero

Model-free approach

We calculated model-free sample ERRs [13] using

$$ \mathrm{Sample}\ \mathrm{ERR} = 1\ \hbox{-} \frac{\mathrm{mean}\ \mathrm{egg}\ \mathrm{count}\ \mathrm{after}\ \mathrm{treatment}}{\mathrm{mean}\ \mathrm{egg}\ \mathrm{count}\ \mathrm{before}\ \mathrm{treatment}}, $$

which we compared with model-based estimates. We used a non-parametric percentile block bootstrap method (Table 1) to calculate associated 95 % confidence intervals (CIs). Block bootstrap methods [27] account for correlation among observations (egg counts) from the same individual by randomly sampling (with replacement) blocks of data; in this case, all of an individual’s egg counts before and after treatment. Details are given in Additional file 1: Supplementary Methods S1 Percentile block bootstrap.

Modelling approach

We employed marginal models and conditional mixed models [17] to estimate ERRs [16]. These are two distinct classes of statistical regression model suitable for analysing dependent (clustered/correlated) data, which here arises because egg counts are measured repeatedly from the same individual and individuals are sampled repeatedly within the same study. We defined two variants of the latter, one in a classical (frequentist) manner, and the other in a Bayesian framework. Full mathematical details of the models are given in the Additional file 1: Supplementary Methods sections S2 Marginal models, S3 Conditional mixed models, and S4 Bayesian conditional mixed models. Here, we give a brief synopsis of the model classes and a description of the salient features for estimating ERRs. Key distinctions are summarised in Table 3.

Table 3 Summary of approaches used to estimate egg reduction rates among children infected with schistosomes following treatment with praziquatel

Marginal models offer population average (marginal) inference, empirically accounting for the dependence of the data using a postulated correlation matrix. In combination with sandwich estimators (Table 1) of coefficient standard errors, marginal models yield robust estimates of uncertainty. Conditional mixed models offer inference at the level of the individual by modelling explicitly the conditional dependence of the data using fixed and random effects. This permits estimation of individual ERRs and the degree of variation among them. By casting conditional mixed models in a Bayesian architecture— defining the necessary parameter prior distributions—one can fully integrate uncertainty into the estimated posteriors and hence derive robust indices of uncertainty, including those associated with the estimated distribution of ERRs among individuals.

The essential ingredient of both model classes (marginal and conditional mixed models) is a log-linear regression structure that describes the change in egg counts after treatment, x = 1, compared to before treatment, x = 0, in a multiplicative fashion. Hence, the accompanying regression coefficient β quantifies the risk ratio (RR) of egg counts after treatment compared to before treatment, and the ERR is given (generically) by 1 – exp(βx). Covariates enter the regression structure as interacting with x. In marginal models this permits ERRs to vary among strata. In conditional mixed models, this permits ERRs to vary also among strata, via fixed effects, and additionally among individuals, via random effects.


We defined marginal and conditional mixed models separately for the S. haematobium and S. mansoni datasets, including the covariates of ERRs listed in Table 4. In the conditional mixed models, these covariates were treated as exerting both fixed and random effects permitting variation among population strata and among individuals within strata. We fitted the models in R [28] using: (a) generalized estimating equation techniques, implemented with the geepack package (marginal models) [29]; (b) restricted maximum likelihood estimation by Laplace approximation, implemented with lme4 (conditional mixed models) [30]; and (c) Markov chain Monte Carlo (MCMC) methods, implemented with MCMCglmm (Bayesian conditional mixed models) [31] (see Table 1 for descriptions of these statistical techniques). We ran three MCMC chains for the Bayesian models, monitoring for convergence and checking that our final conclusions were not dependent on the choice of initial values [32]. In general, 5,000 iterations were discarded as burn-in and an additional 20,000 were sufficient to estimate parameter posterior distributions.

Table 4 Covariates included in the regression models used to estimate egg reduction rates among children infected with schistosomes following treatment with praziquatel


Average egg reduction rates

The model-free sample estimates of S. haematobium and S. mansoni average ERRs following treatment with praziquantel, aggregated across studies, are—with 95 % CIs given in parentheses—99.3 % (98.7 %, 99.7 %) and 83.8 % (77.7 %, 88.9 %), respectively. The corresponding marginal-model estimates (excluding covariates) are 99.6 % (98.1 %, 99.9 %) and 77.9 % (72.7 %, 82.0 %), respectively. The notable difference between the model-free and model-based estimates for S. mansoni is because the marginal model accounts for the correlation among the repeated measures. This is compounded by the high variation among the number of observations per individual, which ranged from 4 to 37 [22].

Stratum average egg reduction rates

The average ERRs and their accompanying CIs estimated for each stratum (defined by the covariates listed in Table 4) using the model-free and marginal model approaches are depicted in Fig. 3. The two approaches yield similar estimates, albeit the assumptions of the modelling approach stabilise estimates in some poorly populated strata, and generally reduce variation. Moreover, model-free bootstrap CIs cannot be constructed in some strata because no eggs were counted in any samples after treatment. Therefore, in these strata, the estimated ERRs are 100 % with no associated uncertainty (grey circles, Fig. 3b). Some of the model-based average ERRs among children infected with S. mansoni fall below the WHO’s 90 % threshold of ‘optimal’ praziquantel efficacy albeit less so at 21 days, which is the WHO-recommended maximum follow-up time (Fig. 3b) [13].

Fig. 3
figure 3

Comparison of egg reduction rates among children infected with schistosomes following treatment with praziquantel estimated by model-free and marginal-model methods. Panels a and b depict, respectively, estimates from individuals infected with Schistosoma haematobium and S. mansoni. Subplots within each panel are stratified according to the different covariate combinations defined by the marginal model; some strata are unpopulated and therefore have no data points. Marginal model and model-free estimates are plotted at each follow-up time for ease of visual comparison. Error bars represent 95 % confidence intervals, calculated using bootstrap methods for model-free sample estimates and using robust sandwich estimators of the standard error for marginal-model estimates. Circular data points (depicting model-free estimates) that are coloured grey do not have an associated uncertainty interval since, in the corresponding strata, all egg counts after treatment were zero, and hence are incompatible with the bootstrap approach. The dashed lines in panel b highlight the decreasing trend in efficacy for increasing follow-up times as estimated by the marginal model fitted to the S. mansoni data (see Fig. 5 for coefficient estimates)

Underlying the marginal-model estimates shown in Fig. 3 are the estimated covariate coefficients. These are presented in the form of RRs in Fig. 4, alongside accompanying 95 % CIs. We also present the ERRs corresponding to these RRs in Table 5. The estimates in Fig. 4a indicate that average S. haematobium egg counts from older SAC (Fig. 5a) after treatment relative to before treatment are approximately 12 times greater than those from younger SAC (P-value = 0.016). However, this seemingly pronounced difference corresponds to an average ERR of 95.9 % (85.8 %, 98.8 %) compared to 99.7 % (99.4 %, 100 %), respectively (with other covariates set to their baseline values, i.e. males from Côte d’Ivoire, Table 5).

Fig. 4
figure 4

Coefficient estimates of covariates associated with average egg reduction rates among children infected with schistosomes following treatment with praziquantel. Panels a and b depict coefficients estimated from the marginal models fitted to the data on, respectively, Schistosoma haematobium and S. mansoni egg counts measured from children before and after treatment with praziquantel. The coefficient point estimates (black circles) indicate the multiplicative change (risk ratio, RR) in egg counts after treatment in a particular covariate group compared to the change after treatment in the reference group. Hence, a RR <1 is associated with an increased efficacy and a RR >1 is associated with a decreased efficacy (compared with the reference group). Error bars depict 95 % confidence intervals (CIs). A covariate is deemed to exert a statistically significant effect only when its CI does not cross the vertical grey line at RR = 1. For example, older school-aged children (SAC) infected with S. haematobium are associated with a statistically significant decrease in efficacy (RR >1) compared to younger SAC

Table 5 The effect of covariates on average egg reduction rates among children infected with schistosomes following treatment with praziquantel
Fig. 5
figure 5

Egg reduction rates among children infected with schistosomes following treatment with praziquantel. Panels a and b depict, respectively, estimates from children infected with Schistosoma haematobium and S. mansoni. Egg reduction rates are calculated from the empirical best linear unbiased predictors (see Table 1 for definition) estimated from the classical (frequentist) conditional mixed models. Negative estimates of ERRs (a: n = 1, 0.59 %; b: n = 24, 2.4 %), which correspond to an increase in egg counts after treatment compared to before treatment, are not shown

The estimates in Fig. 4b highlight the increasing trend in the RRs from 21 days to 28 days to 42 days, corresponding to a decreasing trend in the average ERR (as also evident by the trend lines in Fig. 3b). Egg counts made at 42 days after treatment relative to before treatment were 5.77 times greater than those made at 21 days (P-value = 0.012), corresponding to ERRs of 73.6 and 95.4 % respectively (Table 5). Also noteworthy is that average egg counts following a 60 mg/kg oral dose of praziquantel are 66.8 % lower compared with a 40 mg/kg dose (P-value = 0.064), corresponding to a (not statistically significant) increase in the ERRs from 95.4 to 98.5 %, (with other covariates held at their baseline values, i.e. male younger SAC in Côte d’Ivoire, Table 5).

Individual egg reduction rates

In Fig. 5 we show individual ERRs across all studies, adjusted for covariate fixed effects and estimated by the classical (frequentist) and Bayesian conditional mixed models for S. haematobium (Fig. 5a) and S. mansoni (Fig. 5b). The point estimated ERRs (so-called empirical best linear unbiased predictors, EBLUPs, Table 1) from the classical conditional mixed models (denoted by stars in Fig. 5) indicate that the percentage of individuals with an ERR of greater than 90 % is 97.7 and 80.7 % for S. haematobium and S. mansoni respectively. The corresponding percentages calculated using the Bayesian posterior medians (denoted by dots in Fig. 5) are 96.5 and 80.8 %. However, when the uncertainty in the estimated ERRs is taken into account—by calculating the fraction of individuals with an ERR > 90 % for each draw from the estimated parameter posterior—the median percentages and associated Bayesian credible intervals (BCIs, given in parentheses) are 97.1 % (94.2 %, 98.8 %) and 75.9 % (67.0 %, 81.0 %) for S. haematobium and S. mansoni respectively.

The cumulative distributions (percentiles) of individual ERRs within strata, estimated from the Bayesian conditional mixed models, are depicted in Figs. 6 and 7 for S. haematobium and S. mansoni respectively. The corresponding fractions of individual responses greater than 90 % are given in Table 6. The distributions of ERRs, like the point estimates shown in Fig. 5, show that praziquantel is highly efficacious in the majority of children (ERRs > 90 %) but that a substantial minority have ERRs below the 90 % threshold. In particular, the distributions estimated from children infected with S. mansoni in Uganda, or from those followed-up after 42 days, have longer left tails and greater uncertainty than those estimated from children in Côte d’Ivoire or followed-up after a shorter duration. In Uganda, the median and 95 % BCI associated with the proportion of children with an ERR greater than 90 % is 75.9 % (59.7 %, 86.4 %) compared to 94.4 % (85.3 %, 98.3 %) in Côte d’Ivoire. The corresponding median and 95 % BCI associated with a 42-day follow up is 59.1 % (39.7 %, 76.4 %) compared to 94.4 % (85.3 %, 98.3 %) with a 21-day follow-up.

Fig. 6
figure 6

Cumulative distributions of egg reduction rates among children infected with Schistosoma haematobium following treatment with praziquantel. Cumulative distributions (black lines) are constructed from the posterior distributions of the fixed and random effects components of egg reduction rates estimated from the Bayesian conditional mixed models. Distributions are depicted by country, age group and sex in panels a, b and c respectively. In all panels, covariates not indicated in the legend are set to their baseline levels, i.e. male younger school-aged children from Côte d’Ivoire followed up after 21 days, see Table 4. Grey shaded areas depict 95 % Bayesian credible intervals

Fig. 7
figure 7

Cumulative distributions of egg reduction rates among children infected with Schistosoma mansoni following treatment with praziquantel. Cumulative distributions (black lines) are constructed from the posterior distributions of the fixed and random effects components of egg reduction rates estimated from the Bayesian conditional mixed models. Distributions are depicted by country, age group, sex, dose and follow-up days in panels a, b, c, d and e respectively. In all panels, covariates not indicated in the legend are set to their baseline levels, i.e. male younger school-aged children given 40 mg/kg praziquantel from Côte d’Ivoire followed up after 28 days, see Table 4. Grey shaded areas depict 95 % Bayesian credible intervals

Table 6 The effect of covariates on the percentage of egg reduction rates greater than 90 % among children infected with schistosomes following treatment with praziquantel


Modelling methods are currently underused for assessing anthelminthic drug efficacy. By analysing longitudinal data on the intensity of schistosome infections before and after treatment with praziquantel, we show that marginal and conditional mixed models can be used to obtain robust estimates of both population- and individual-level efficacies, while concurrently evaluating the effects of covariates. While a small number of anthelminthic drug studies have employed various somewhat unconventional approaches [3335], statistical modelling techniques have not translated into more general use in either the veterinary or human fields (but see [36]). Beyond the realm of estimating anthelminthic drug efficacy, longitudinal and hierarchical modelling techniques are often used incorrectly across a variety of disciplines in ecology and evolutionary biology [37], indicating a general lack of applied understanding about these powerful analytical tools.

The efficacy of praziquantel within populations and among individuals

Model-free sample-based methods for estimating anthelminthic drug efficacy are hindered by their inherent inefficiency in handling covariates (although we note that the WHO protocol for measuring ERRs recommends evaluating ERRs at a standardized 21-day follow-up in SAC using a single test on a single sample before and after treatment) [16]. Moreover, they can produce biased estimates when individuals are assessed with different numbers of repeated measures (e.g. different numbers of Kato-Katz counts per stool sample, or different numbers of stool samples, either before or after treatment). This is because all observations are weighted equally and, therefore, individuals contributing more observations disproportionately influence the estimated statistic [38], i.e. the sample ERR. This probably explains some of the differences between the model-free and marginal-model estimates of ERR presented here. Confidence intervals associated with average ERRs estimated using the marginal models capture the effect of correlated repeated measures using robust sandwich estimators of the standard errors (Table 1). Although the block bootstrap method used to construct CIs for the model-free approach accounts for correlation in a more rudimentary fashion, it is inappropriate in strata when no egg counts are observed after treatment, generating a statistically invalid CI with a width of zero and a 100 % point-estimated ERR (grey circles in Fig. 3b).

The estimated ERRs of praziquantel against S. haematobium from both the model-free and marginal modelling approaches were noticeably higher than those for S. mansoni (Fig. 4). It is possible that this difference is driven by a mechanistic superiority of praziquantel against S. haematobium, for which there is some supportive in vitro evidence [39]. It is also possible that the difference is driven by a differential distribution of praziquantel to the parts of the perivesical venous plexus and the mesenteric/rectal veins where adult S. haematobium and S. mansoni flukes reside, respectively [40, 41]. Another possibility is that the discrepancy in ERRs reflects the presence of praziquantel-tolerant S. mansoni, which albeit rare, has been previously observed in the field [12]. By contrast, tolerant S. haematobium specimens have not been [8], notwithstanding the reports of individuals requiring multiple doses of praziquantel to clear infection [42, 43].

Perhaps most pertinent to the inter-species difference is that the modelling approaches illustrated here do not account for the sensitivity of the diagnostic method (Kato-Katz thick smear for S. mansoni and urine filtration for S. haematobium diagnosis) and, in particular, the manner in which sensitivity declines with decreasing infection intensity, before and after praziquantel administration [44, 45]. This phenomenon can cause overestimation of drug efficacy [23, 46] with differential bias between species if the severity of this effect is different between the urine filtration technique and the Kato-Katz method. Moreover, it is assumed that the specificity of the diagnostic methods are 100 %, with no incorrect identification of schistosome eggs in truly negative samples. In particular, as the intensity of infection declines after treatment, the number of true negative samples increases, increasing the probability of false positives, even with very good specificity. The next generation of modelling frameworks must account for the sensitivity and specificity of the diagnostic method and capture how it varies with the underlying intensity of infection.

While marginal models are powerful tools for assessing stratum-level average efficacy (average ERRs), conditional mixed models offer additional insight at the level of the individual. The individual estimates of ERRs reveal numerous so-called sub-optimally responding individuals, with an efficacy lower than the WHO’s empirical 90 % threshold for group average. A similar result is presented in the recent sample-based analysis of the full 13-study database (see Additional file 1: Supplementary Tables, Table S1) [14]. A substantive 3.7 % of children infected with S. mansoni had ERRs of less than 50 % (Fig. 6) and more than twice as many children infected with S. mansoni had an ERR below 90 % compared with those infected with S. haematobium (10.4 % versus 4.0 %). Suboptimal or atypical responses to praziquantel have been previously described in terms of the number of doses necessary to clear the parasite, with any result above one dose considered suboptimal [42, 43, 4749]. Compared to this measure of sub-optimality, conditional mixed models offer a more nuanced insight into the distribution of drug responses among individuals.

When does the efficacy of praziquantel appear greatest?

The observed efficacy of an anthelminthic crucially depends on when after treatment measurements of infection intensity (or presence/absence of transmission stages) are made. This is because efficacy is a snapshot of the competing dynamics of parasite clearance, followed by reinfection [50] or repopulation initiated by surviving parasites, perhaps subjected to temporarily reduced or inhibited fertility or, in the case of praziquantel, maturation of juvenile parasites [45].

Our results indicate that, among children infected with S. mansoni, average ERRs decreased from follow-up times of 21 days to 28 days, and from 28 days to 42 days. This mirrors the findings of a meta-analysis of 11 studies of the efficacy of praziquantel against S. mansoni in Africa where a similar relationship between longer follow-ups and lower CRs was interpreted as being probably due to reinfection [48]. Repopulation by surviving juvenile parasites is another likely cause. Schistosomes exhibit a biphasic susceptibility to praziquantel throughout their lifecycle; early-stage migrating larvae are susceptible to praziquantel, but after about one month of growth, susceptibility drops precipitously, and is only regained after another two months [51]. Therefore, by 42 days after treatment, many of the juvenile schistosomes that survived treatment may have matured or repopulated the organs of preferred location as fertile egg-producing adult schistosomes. Moreover, recent studies using circulating cathodic antigen (CCA) tests indicate that 50–100 % of praziquantel-treated children retain their antigenemia/antigenuria when retested 1–7 weeks after treatment [52, 53]. Hence, the effective praziquantel efficacy is probably much lower than previously thought and surviving parasites are also likely contributors to the recovery in egg counts after treatment. We did not have access to data collected at shorter follow-up times, to perhaps model when ERRs are at a maximum, although analyses presented elsewhere [23] have suggested that this occurs 2–3 weeks after treatment, which is the timeframe currently recommended by WHO [13].

Why does age affect the efficacy of praziquantel?

The decreasing (weak and not statistically significant) trend in estimated ERRs with increasing age (pre-SAC to younger SAC to older SAC, see Figs. 5b and 7c) of children infected with S. mansoni is in accordance with several existing hypotheses. Older children are likely to have stronger acquired immunity to schistosomes than younger children [54], as flukes killed or damaged from exposure to praziquantel release previously ‘unseen’ antigens [5558] eliciting protective immune responses thought to enhance (but not always [59]) the efficacy of subsequent treatments [60, 61]. Consequently, one might expect older children to be more amenable to treatment and exhibit higher ERRs than their younger counterparts. However, and crucially, the data analysed here were collected from communities unexposed or minimally exposed to praziquantel MDA (see selection criteria in Fig. 1). Hence, older children in these communities might be more difficult to treat and respond less well to praziquantel having been left to grow older with untreated schistosome infections. This might explain the estimates from the marginal model indicating that the average ERR from older SAC infected with S. haematobium is markedly lower than that from younger SAC. Alternatively, this result may reflect a general limitation of the analysis. For example, most of the studies under consideration encompassed multiple villages or study areas within a single country, but village-level identifiers were not available, and small-scale variation in praziquantel responses among villages can be quite large [62].

How does infection intensity affect the efficacy of praziquantel?

Numerous studies have demonstrated a negative association between the infection intensity before treatment and estimated CRs following treatment with praziquantel; the higher the intensity, the lower the CR [20, 48, 63]. This is because, if adult schistosomes die with a fixed probability when exposed to praziquantel, cure will be less likely in heavily infected individuals than in lightly infected individuals. Indeed, the WHO no longer recommends using CRs for the monitoring and evaluation (M&E) of anthelmintic efficacy [13] because it is impossible to observe incremental reductions in parasite burden using a binary measure of cure and, therefore, CRs fail to capture the impact of multiple doses of anthelminthic drugs over the course of PCT programmes [64]. In this work, we constructed the conditional mixed models to estimate the association between an individual’s egg count before treatment (the random intercept term) and their ERR (governed by the random ‘gradient’ term, see Additional file 1Supplementary Methods, S3 Conditional mixed models), but we found no statistically significant relationship. Intuitively, again invoking the assumption of a constant probability of death by praziquantel, the percentage reduction in intensity will be constant, explaining why there is no association between egg counts before treatment and the estimated ERR.

Modelling for M&E of MDA interventions

Modelling has an important role in the M&E of anthelmintic drug efficacy in the context of human helminth PC programmes. Marginal models offer a robust method of estimating (sub-) population-level drug efficacy, which would allow disease control managers to identify whether or not target ERRs are being met. Such targets could be defined using data from populations predominantly naïve to MDA, in different demographic groups, which may respond in a systematically differential manner to drug treatment. Individual-level ERRs estimated using conditional mixed models offer an additional depth of insight, permitting characterisation of the distribution of drug responses among individuals. This is important for the rapid identification of changing responses to anthelminthic drugs that may be indicative of declining drug efficacy, potentially caused by emerging drug-resistant parasites [6, 1012, 49]. By comparing the observed distribution of drug responses to a reference distribution of expected responses estimated before MDA—ideally from the same community, but otherwise using data from demographically and geographically matched or partially matched populations—it would be possible to identify whether or not individuals are responding aberrantly to the drug. That is, one could quantify how atypical an observed response is compared to the usual or expected distribution of responses. Individuals responding suspiciously could be investigated, perhaps using in vitro drug sensitivity tests on the infecting parasites. Moreover, over multiple rounds of MDA, one could identify shifts in the distribution of responses from the original reference distribution; shifts towards decreased efficacy triggering further programmatic and parasitological investigation.


Marginal and conditional mixed models are robust approaches for calculating population- and individual-level estimates of anthelminthic drug efficacy. We illustrate these techniques by analysing data collated from nine previous studies on schistome egg counts from children before and after administration of praziquantel. We show that model-based analyses: (a) offer more stable and robust estimates of average ERRs compared to traditional sample-based methods, especially when sample sizes are small; (b) can be used to evaluate how and to what degree drug responses vary among population strata, in terms of an average response, and among individuals within a stratum, in terms of the distribution of individual responses. We show that Bayesian methods are particularly useful in quantifying uncertainties, and permit creation of prototype ‘reference’ distributions describing the range of drug responses expected in communities predominantly naïve to MDA. These distributions have potentially important applications to the M&E of anthelmintic efficacy in helminthiasis PCT programmes, particularly for identifying individual atypical responses and distributional shifts, potentially indicative of emerging drug resistance. Therfore, the approaches illustrated in this paper have an important role in supporting the control and elimination of human helminthiases.



Bayesian credible interval


circulating cathodic antigen


confidence interval


cure rate


egg reduction rate


generalized estimating equation


generalized linear model


generalized linear mixed model


intensity reduction rate


monitoring & evaluation


Markov chain Monte Carlo


mass drug administration


neglected tropical disease


preventative chemotherapy


school-aged children


World Health Organization


  1. Reddy M, Gill SS, Kalkar SR, Wu W, Anderson PJ, Rochon PA. Oral drug therapy for multiple neglected tropical diseases: a systematic review. JAMA. 2007;298:1911–24.

    Article  CAS  PubMed  Google Scholar 

  2. World Health Organization. Accelerating Work to Overcome the Global Impact of Neglected Tropical Diseases: A Roadmap for Implementation. Geneva: World Health Organization; 2012.

    Google Scholar 

  3. Hotez PJ, Alvarado M, Basáñez MG, Bolliger I, Bourne R, Boussinesq M, et al. The global burden of disease study 2010: interpretation and implication for the neglected tropical diseases. PLoS Negl Trop Dis. 2014;8, e2865.

    Article  PubMed Central  PubMed  Google Scholar 

  4. Murray CJL, Vos T, Lozano R, Naghavi M, Flaxman AD, Michaud C, et al. Disability-adjusted life years (DALYs) for 291 diseases and injuries in 21 regions, 1990–2010: a systematic analysis for the Global burden of disease study 2010. Lancet. 2012;380:2197–223.

    Article  PubMed  Google Scholar 

  5. World Health Organization. Preventative Chemotherapy in Human Helminthiasis: Coordinated use of Anthelmintic Drugs in Control Interventions: A Manual for Health Professionals and Programme Managers. Geneva: World Health Organization; 2006. p. 1–62.

    Google Scholar 

  6. Bockarie MJ, Kelly-Hope LA, Rebollo M, Molyneux DH. Preventive chemotherapy as a strategy for elimination of neglected tropical parasitic diseases: endgame challenges. Phil Trans R Soc B. 2013;368:20120144.

    Article  PubMed Central  PubMed  Google Scholar 

  7. Utzinger J, Becker SL, Knopp S, Blum J, Neumayr AL, Keiser J, et al. Neglected tropical diseases: diagnosis, clinical management, treatment and control. Swiss Med Wkly. 2012;142:w13727.

    PubMed  Google Scholar 

  8. Wang W, Wang L, Liang Y. Susceptibility or resistance of praziquantel in human schistosomiasis: a review. Parasitol Res. 2012;111:1871–7.

    Article  PubMed  Google Scholar 

  9. Albonico M, Engels D, Savioli L. Monitoring drug efficacy and early detection of drug resistance in human soil-transmitted nematodes: a pressing public health agenda. Int J Parasitol. 2004;34:1205–10.

    Article  CAS  PubMed  Google Scholar 

  10. Fenwick A, Webster JP, Bosque-Oliva E, Blair L, Fleming F, Zhang Y, et al. The schistosomiasis control initiative (SCI): rationale, development and implementation from 2002–2008. Parasitology. 2009;136:1719–30.

    Article  CAS  PubMed  Google Scholar 

  11. Basáñez MG, McCarthy JS, French MD, Yang GJ, Walker M, Gambhir M, et al. A research agenda for helminth diseases of humans: modelling for control and elimination. PLoS Negl Trop Dis. 2012;6:e1548.

    Article  PubMed Central  PubMed  Google Scholar 

  12. Webster JP, Molyneux DH, Hotez PJ, Fenwick A. The contribution of mass drug administration to global health: past, present and future. Phil Trans R Soc B. 2014;369:2013034.

    Article  Google Scholar 

  13. World Health Organization. Assessing the Efficacy of Anthelminthic Drugs Against Schistosomiasis and Soil-Transmitted Helminthiases. Geneva: World Health Organization; 2013.

    Google Scholar 

  14. Olliaro PL, Vaillant M, Diawara A, Coulibaly JT, Garba A, Keiser J, et al. Toward measuring Schistosoma response to praziquantel treatment with appropriate descriptors of egg excretion. PLoS Negl Trop Dis. 2015;9:e0003821.

    Article  PubMed Central  PubMed  Google Scholar 

  15. Cabaret J, Berrag B. Faecal egg count reduction test for assessing anthelmintic efficacy: average versus individually based estimations. Vet Parasitol. 2004;121:105–13.

    Article  CAS  PubMed  Google Scholar 

  16. Walker M, Churcher TS, Basáñez MG. Models for measuring anthelmintic drug efficacy for parasitologists. Trends Parasitol. 2014;30:528–37.

    Article  CAS  PubMed  Google Scholar 

  17. Diggle PJ, Heagerty PJ, Liang KY, Zeger SL. Analysis of Longitudinal Data. Oxford: Oxford University Press; 2002.

    Google Scholar 

  18. Subramanian SV, O’Malley AJ. Modeling neighborhood effects: the futility of comparing mixed and marginal approaches. Epidemiol. 2010;21:475–8.

    Article  CAS  Google Scholar 

  19. Olds GR, King C, Hewlett J, Olveda R, Wu G, Ouma J, et al. Double-blind placebo-controlled study of concurrent administration of albendazole and praziquantel in schoolchildren with schistosomiasis and geohelminths. J Infect Dis. 1999;179:996–1003.

    Article  CAS  PubMed  Google Scholar 

  20. Coulibaly JT, Fürst T, Silué KD, Knopp S, Hauri D, Ouattara M, et al. Intestinal parasitic infections in schoolchildren in different settings of Côte d’Ivoire: effect of diagnostic approach and implications for control. Parasit Vectors. 2012;5:135.

    Article  PubMed Central  PubMed  Google Scholar 

  21. Stete K, Krauth SJ, Coulibaly JT, Knopp S, Hattendorf J, Müller I, et al. Dynamics of Schistosoma haematobium egg output and associated infection parameters following treatment with praziquantel in school-aged children. Parasit Vectors. 2012;5:298.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. Utzinger J, N’Goran EK, N’dri A, Lengeler C, Tanner M. Efficacy of praziquantel against Schistosoma mansoni with particular consideration for intensity of infection. Trop Med Int Health. 2000;5:771–8.

    Article  CAS  PubMed  Google Scholar 

  23. Utzinger J, Booth M, N’Goran EK, Müller I, Tanner M, Lengeler C. Relative contribution of day-to-day and intra-specimen variation in faecal egg counts of Schistosoma mansoni before and after treatment with praziquantel. Parasitology. 2001;122:537–44.

    Article  CAS  PubMed  Google Scholar 

  24. Raso G, N’Goran EK, Toty A, Luginbühl A, Adjoua CA, Tian-Bi NT, et al. Efficacy and side effects of praziquantel against Schistosoma mansoni in a community of western Côte d’Ivoire. Trans R Soc Trop Med Hyg. 2004;98:18–27.

    Article  CAS  PubMed  Google Scholar 

  25. Scherrer AU, Sjöberg MK, Allangba A, Traoré M, Lohourignon LK, Tschannen AB, et al. Sequential analysis of helminth egg output in human stool samples following albendazole and praziquantel administration. Acta Trop. 2009;109:226–31.

    Article  CAS  PubMed  Google Scholar 

  26. Sousa-Figueiredo JC, Betson M, Atuhaire A, Arinaitwe M, Navaratnam AM, Kabatereine NB, et al. Performance and safety of praziquantel for treatment of intestinal schistosomiasis in infants and preschool children. PLoS Negl Trop Dis. 2012;6, e1864.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  27. Lahiri SN. Resampling Methods for Dependent Data. New York: Springer; 2003.

    Book  Google Scholar 

  28. R Core Team. R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing; 2015.

    Google Scholar 

  29. Halekoh U, Højsgaard S, Yan J. The R package geepack for generalized estimating equations. J Stat Softw. 2006;15:1–11.

    Article  Google Scholar 

  30. Bates D, Maechler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67:1–48.

    Article  Google Scholar 

  31. Hadfield JD. MCMC methods for multi-response generalized linear mixed models: the MCMCglmm R package. J Stat Softw. 2010;33:1–22.

    Article  Google Scholar 

  32. Gelman AC, Stern J, Rubin DB. Bayesian Data Analysis. London: Chapman & Hall; 2004.

    Google Scholar 

  33. Mejía ME, Igartúa BMF, Schmidt EE, Cabaret J. Multispecies and multiple anthelmintic resistance on cattle nematodes in a farm in Argentina: the beginning of high resistance? Vet Res. 2003;34:461–7.

    Article  PubMed  Google Scholar 

  34. Nielsen M, Vidyashankar A, Hanlon B, Diao G, Petersen S, Kaplan R. Hierarchical model for evaluating pyrantel efficacy against strongyle parasites in horses. Vet Parasitol. 2013;197:614–22.

    Article  CAS  PubMed  Google Scholar 

  35. Vidyashankar A, Kaplan R, Chan S. Statistical approach to measure the efficacy of anthelmintic treatment on horse farms. Parasitology. 2007;134:2027–39.

    Article  CAS  PubMed  Google Scholar 

  36. Levecke B, Anderson RM, Berkvens D, Charlier J, Devleesschauwer B, Speybroeck N. Mathematical inference in helminth egg counts in stool and its applications in mass drug administration programmes to control soil-transmitted helminthiases in public health. Adv Parasitol. 2015;87:193–247.

    Article  PubMed  Google Scholar 

  37. Bolker BM, Brooks ME, Clark CJ, Geange SW, Poulsen JR, Stevens MHH, et al. Generalized linear mixed models: a practical guide for ecology and evolution. Trends Ecol Evol. 2009;24:127–35.

    Article  PubMed  Google Scholar 

  38. Hanley JA, Negassa A, Edwardes MD, Forrester JE. Statistical analysis of correlated data using generalized estimating equations: an orientation. Am J Epidemiol. 2003;157:364–75.

    Article  PubMed  Google Scholar 

  39. Andrews P. A summary of the efficacy of praziquantel against schistosomes in animal experiments and notes on its mode of action. Arzneimittelforschung. 1981;31:538–41.

    CAS  PubMed  Google Scholar 

  40. Anayi S, Al-Nasiri N. Acute mesenteric ischaemia caused by Schistosoma mansoni infection. Br Med J (Clin Res Ed). 1987;294:1197.

    Article  CAS  Google Scholar 

  41. Smith J, Christie J. The pathology of human Schistosoma haematobium infections. In: Miller MJ, Love EJ, editors. Parasitic Diseases: Treatment and Control. Boca-Raton: CRC Press; 1989.

    Google Scholar 

  42. Silva IM, Thiengo R, Conceição MJ, Rey L, Lenzi HL, Pereira Filho E, et al. Therapeutic failure of praziquantel in the treatment of Schistosoma haematobium infection in Brazilians returning from Africa. Mem Inst Oswaldo Cruz. 2005;100:445–9.

    PubMed  Google Scholar 

  43. Alonso D, Munoz J, Gascon J, Valls ME, Corachan M. Failure of standard treatment with praziquantel in two returned travelers with Schistosoma haematobium infection. Am J Trop Med Hyg. 2006;74:342–4.

    PubMed  Google Scholar 

  44. Booth M, Vounatsou P, N’Goran EK, Tanner M, Utzinger J. The influence if sampling effort and the performance of the Kato-Katz technique in diagnosing Schistosoma mansoni and hookworm co-infections in rural Côte d’Ivoire. Parasitology. 2003;127:525–31.

    Article  CAS  PubMed  Google Scholar 

  45. Lamberton PH, Kabatereine NB, Oguttu DW, Fenwick A, Webster JP. Sensitivity and specificity of multiple Kato-Katz thick smears and a circulating cathodic antigen test for Schistosoma mansoni diagnosis pre- and post-repeated-praziquantel treatment. PLoS Negl Trop Dis. 2014;8:e3139.

    Article  PubMed Central  PubMed  Google Scholar 

  46. de Vlas SJ, Gryseels B. Underestimation of Schistosoma mansoni prevalences. Parasitol Today. 1992;8:274–7.

    Article  PubMed  Google Scholar 

  47. Ismail M, Metwally A, Farghaly A, Bruce J, Tao LF, Bennett JL. Characterization of isolates of Schistosoma mansoni from Egyptian villagers that tolerate high doses of praziquantel. Am J Trop Med Hyg. 1996;55:214–8.

    CAS  PubMed  Google Scholar 

  48. Danso-Appiah A, de Vlas SJ. Interpreting low praziquantel cure rates of Schistosoma mansoni infections in Senegal. Trends Parasitol. 2002;18:125–9.

    Article  PubMed  Google Scholar 

  49. Melman SD, Steinauer ML, Cunningham C, Kubatko LS, Mwangi IN, Wynn NB, et al. Reduced susceptibility to praziquantel among naturally occurring Kenyan isolates of Schistosoma mansoni. PLoS Negl Trop Dis. 2009;3, e504.

    Article  PubMed Central  PubMed  Google Scholar 

  50. Jia TW, Melville S, Utzinger J, King CH, Zhou XN. Soil-transmitted helminth reinfection after drug treatment: a systematic review and meta-analysis. PLoS Negl Trop Dis. 2012;6, e1621.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  51. Doenhoff MJ, Cioli D, Utzinger J. Praziquantel: mechanisms of action, resistance and new derivatives for schistosomiasis. Curr Opin Infect Dis. 2008;21:659–67.

    Article  CAS  PubMed  Google Scholar 

  52. Mwinzi PNM, Kittur N, Ochola E, Copper PJ, Campbell Jr CH, King CH, et al. Additional evaluation of the point-of-contact circulating cathodic antigen assay for Schistosoma mansoni infection. Front Public Health. 2015;3:1–8.

    Article  Google Scholar 

  53. Bustinduy AL. Pharmacokinetics and pharmacodynamics study of praziquantel in at 40 mg/kg or 60 mg/kg dosing in young Ugandan children infected with Schistosoma mansoni. In: The Burden of Schistosomiasis Morbidity in African Children: Validating Novel, low Cost and Scalable Detection Tools and Optimizing Praziquantel Treatment, PhD thesis. Amsterdam: University of Amsterdam; 2015.

    Google Scholar 

  54. Butterworth AE. Immunological aspects of human schistosomiasis. Br Med Bull. 1998;54:357–68.

    Article  CAS  PubMed  Google Scholar 

  55. Harnett W, Kusel J. Increased exposure of parasite antigens at the surface of adult male Schistosoma mansoni exposed to praziquantel in vitro. Parasitology. 1986;93:401–5.

    Article  CAS  PubMed  Google Scholar 

  56. Doenhoff MJ, Sabah AA, Fletcher C, Webbe G, Bain J. Evidence for an immune-dependent action of praziquantel on Schistosoma mansoni in mice. Trans R Soc Trop Med Hyg. 1987;81:947–51.

    Article  CAS  PubMed  Google Scholar 

  57. Brindley PJ, Sher A. Immunological involvement in the efficacy of praziquantel. Exp Parasitol. 1990;71:245–8.

    Article  CAS  PubMed  Google Scholar 

  58. Modha J, Lambertucci J, Doenhoff M, McLaren DJ. Immune dependence of schistosomicidal chemotherapy: an ultrastructural study of Schistosoma mansoni adult worms exposed to praziquantel and immune serum in vivo. Parasite Immunol. 1990;12:321–34.

    Article  CAS  PubMed  Google Scholar 

  59. King CH, Muchiri EM, Ouma JH. Evidence against rapid emergence of praziquantel resistance in Schistosoma haematobium, Kenya. Emerg Infect Dis. 2000;6:585–94.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  60. Woolhouse MEJ, Hagan P. Seeking the ghost of worms past. Nat Med. 1999;5:1225–7.

    Article  CAS  PubMed  Google Scholar 

  61. Black CL, Steinauer ML, Mwinzi PN, Evan Secor W, Karanja D, Colley DG. Impact of intense, longitudinal retreatment with praziquantel on cure rates of schistosomiasis mansoni in a cohort of occupationally exposed adults in western Kenya. Trop Med Int Health. 2009;14:450–7.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  62. Thiongo F, Mbugua G, Ouma J, Sturrock R. Efficacy of oxamniquine and praziquantel in school children from two Schistosoma mansoni endemic areas. East Afr Med J. 2002;79:29–33.

    Article  CAS  Google Scholar 

  63. Fallon P, Hamilton J, Doenhoff M. Efficacy of treatment of murine Schistosoma mansoni infections with praziquantel and oxamniquine correlates with infection intensity: role of host antibody. Parasitology. 1995;111:59–66.

    Article  CAS  PubMed  Google Scholar 

  64. Montresor A. Cure rate is not a valid indicator for assessing drug efficacy and impact of preventive chemotherapy interventions against schistosomiasis and soil-transmitted helminthiasis. Trans R Soc Trop Med Hyg. 2011;105:361–3.

    Article  PubMed  Google Scholar 

Download references


MW and MGB acknowledge funding from the Wellcome Trust (, grant 092677/Z/10/Z) and also acknowledge that this investigation received financial support from the UNICEF/UNDP/World Bank WHO Special Programme for Research and Training in Tropical Disease (TDR). TSM acknowledges funding from the Thouron Award (, a graduate exchange programme between British universities and the University of Pennsylvania, to pursue an MSc in epidemiology at Imperial College London. CHK acknowledges funding from a WHO Project Development Unit award (ID-910357) and JRS acknowledges funding from the Wellcome Trust (WT085440MA). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.


The opinions expressed in this paper are those of the authors and may not reflect those of their employing organizations. PLO is a staff member of the World Health Organization (WHO); the authors alone are responsible for the views expressed in this publication and it does not necessarily represent the decisions, policy, or views of their organizations.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Martin Walker.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

JTC, CHK, GR, AUS, JRS, JCSF, KS and JU coordinated fieldwork and data collection for the original studies. MW, PLO and MGB provided overall project supervision for the analyses. MW, TSM and MGB developed the statistical approach. MW and TSM conducted the statistical analyses. MW and TSM wrote the first draft of the manuscript. All authors read and approved the final version of the manuscript.

Martin Walker and Tarub S. Mabud contributed equally to this work.

Additional file

Additional file 1:

Supplementary Methods, Figures and Tables. (DOCX 108 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, 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 ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Walker, M., Mabud, T.S., Olliaro, P.L. et al. New approaches to measuring anthelminthic drug efficacy: parasitological responses of childhood schistosome infections to treatment with praziquantel. Parasites Vectors 9, 41 (2016).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: