# Macrofilaricidal Efficacy of Repeated Doses of Ivermectin for the Treatment of River Blindness

Martin Walker, Sébastien Pion, Hanwei Fang, J. Gardon, Joseph Kamgno,

Maria-Gloria Basáñez, Michel Boussinesq

► To cite this version:

Martin Walker, Sébastien Pion, Hanwei Fang, J. Gardon, Joseph Kamgno, et al. Macrofilaricidal Efficacy of Repeated Doses of Ivermectin for the Treatment of River Blindness. Clinical Infectious Diseases, 2017, 65 (12), pp.2026-2034. 10.1093/cid/cix616. hal-02429016

HAL Id: hal-02429016

https://hal.science/hal-02429016

Submitted on 16 Feb 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Macrofilaricidal Efficacy of Repeated Doses of Ivermectin for the Treatment of River Blindness

Martin Walker,¹² Sébastien D. S. Pion,³ Hanwei Fang,⁴ Jacques Gardon,⁵ Joseph Kamgno,⁶,⁷ Maria-Gloria Basáñez,² and Michel Boussinesq³

**Background.** Mass drug administration (MDA) with ivermectin is the cornerstone of efforts to eliminate human onchocerciasis by 2020 or 2025. The feasibility of elimination crucially depends on the effects of multiple ivermectin doses on *Onchocerca volvulus*. A single ivermectin (standard) dose clears the skin-dwelling microfilarial progeny of adult worms (macrofilariae) and temporarily impedes the release of such progeny by female macrofilariae, but a macrofilaricidal effect has been deemed minimal. Multiple doses of ivermectin may cumulatively and permanently reduce the fertility and shorten the lifespan of adult females. However, rigorous quantification of these effects necessitates interrogating longitudinal data on macrofilariae with suitably powerful analytical techniques.

**Methods.** Using a novel mathematical modeling approach, we analyzed, at an individual participant level, longitudinal data on viability and fertility of female worms from the single most comprehensive multiple-dose clinical trial of ivermectin, comparing 3-monthly with annual treatments administered for 3 years in Cameroon.

**Results.** Multiple doses of ivermectin have a partial macrofilaricidal and a modest permanent sterilizing effect after 4 or more consecutive treatments, even at routine MDA doses (150 µg/kg) and frequencies (annual). The life expectancy of adult *O. volvulus* is reduced by approximately 50% and 70% after 3 years of annual or 3-monthly (quarterly) exposures to ivermectin.

**Conclusions.** Our quantification of macrofilaricidal and sterilizing effects of ivermectin should be incorporated into transmission models to inform onchocerciasis elimination efforts in Africa and residual foci in Latin America. It also provides a framework to assess macrofilaricidal candidate drugs currently under development.

**Keywords.** onchocerciasis; river blindness; ivermectin; macrofilaricide; multiple dose.

In 1987 Merck & Co. Inc. committed to donate ivermectin for as long as needed to control human onchocerciasis. Subsequently, more than 2 billion treatments have been distributed in 29 countries throughout sub-Saharan Africa and in 6 countries in Latin America (http://www.merck.com/about/featured-stories/mectizan.html) by mass drug administration (MDA) programs. The Onchocerciasis Elimination Program for the Americas (OEPA) was launched in 1992 and the African Programme for Onchocerciasis Control (APOC) was launched in 1995 to coordinate these efforts, maintain high MDA coverage, and control onchocerciasis as a public health problem. Since release of the World Health Organization's Roadmap on Neglected Tropical Diseases

(NTDs) [1], the objective (in Africa) has been refocused toward elimination of the infection (OEPA aimed at elimination from the outset).

Barring the Amazonian onchocerciasis focus [2], elimination in Latin America is broadly on course [3]. In sub-Saharan Africa, where 99% of the infection cases occur, and notwithstanding some notable elimination successes [4–6], transmission is ongoing in many communities despite more than 15–20 rounds of annual (or 6-monthly) ivermectin treatments [7–10]. Addressing the question of whether elimination in Africa is feasible using current intervention strategies is a key issue. Before its closure in December 2015, APOC released guidance on alternative treatment strategies to accelerate progress toward elimination [11].

Mathematical modeling has demonstrated that the effectiveness and cost effectiveness of ivermectin MDA and the timeframes for elimination crucially depend on the long-term antifilarial action of ivermectin [12]. Ivermectin rapidly clears the skin-dwelling microfilarial progeny of *Onchocerca volvulus*, which is the filarial nematode that causes onchocerciasis, and temporarily prevents the release of new microfilariae by female worms. The dynamics of these so-called microfilaricidal and temporary sterilizing effects are well understood following a

<footnote n="1">¹Department of Pathobiology and Population Sciences and London Centre for Neglected Tropical Disease Research, Royal Veterinary College, Hatfield, and ²Department of Infectious Disease Epidemiology and London Centre for Neglected Tropical Disease Research, Imperial College London, United Kingdom; ³UMI233-TransVIHMI, Institut de Recherche pour le Développement, INSERM U1175, University of Montpellier, France; ⁴Department of Microbiology and Molecular Medicine, University of Geneva Medical School, Switzerland; ⁵Hydrosciences Montpellier, Institut de Recherche pour le Développement, France; and ⁶Centre for Research on Filariasis and Other Tropical Diseases, and ⁷Faculty of Medicine and Biomedical Sciences University of Yaoundé I, Yaoundé, Cameroon.</footnote>

<footnote>Received 16 March 2017; editorial decision 23 June 2017; accepted 14 July 2017; published online July 19, 2017.</footnote>

<footnote>M. W. and S. D. S. P. contributed equally to this work. M.-G. B. and M. B. contributed equally to this work.</footnote>

<footnote>Correspondence: M. Walker, Department of Pathobiology and Population Health Sciences & London Centre for Neglected Tropical Disease Research, Royal Veterinary College, Hawkshead Lane, Hatfield, Hertfordshire, AL9 7TA, UK (mwalker@rvc.ac.uk).</footnote>

<footnote>**Clinical Infectious Diseases®** 2017;65(12):2026–34</footnote>

<footnote>© The Author 2017. Published by Oxford University Press for the Infectious Diseases Society of America. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited. DOI: 10.1093/cid/cix616</footnote>
single standard ivermectin dose (150 µg/kg) [13]. However, because adult worms are long lived, with a life span of approximately 10 years [14], they are exposed to multiple ivermectin doses during an MDA program. Therefore, quantification of cumulative antifilarial effects on repeatedly exposed worms is essential for understanding the long-term epidemiological impact of MDA.

To date, cumulative antifilarial effects have been documented in an indirect manner (estimating their magnitude by fitting models to microfilarial trends following several treatment rounds), yielding conflicting results [15, 16]. However, direct observations on adult worm viability and reproductive status have not been linked to models for the underlying population biology of *O. volvulus*, a powerful approach that can translate clinical trial data into measures of drug efficacy [17]. Macrofilaricidal activity has been inferred from the smaller proportion of live worms extracted from participants given high-intensity regimens (frequency of more than twice per year; dose >150 µg/kg) compared to those given standard regimens [18] or not treated [19]. Similarly, fertility indicators from female *O. volvulus* exposed to high-intensity regimens are lower than those from females exposed to standard regimens, yet such indicators have not been translated into cumulative and potentially permanent effects of incremental exposures to ivermectin.

Here, we use our recently developed mathematical framework [17] to analyze longitudinal data on the viability and fertility of female adult *O. volvulus* collected from participants of the single largest and most recent randomized control trial, which was conducted in an endemic area of Cameroon [18] (a summary of other smaller and older trials is provided in the Supplementary Materials, Supplementary Figures 1–4 and Supplementary Table 1). We estimate the macrofilaricidal and antifertility activity of multiple ivermectin doses given at standard (annual) and high (3-monthly) frequencies and relate our findings to observations from previous studies on the antifilarial activity of multiple ivermectin doses and other macrofilaricidal therapies. We consider how these effects should be incorporated into transmission models to support onchocerciasis control and elimination and discuss how our findings can be used to assess the performance of candidate macrofilaricidal drugs.

# METHODS

## Trial Design and Data Collection

The data were collected between 1994 and 1998 during a clinical trial undertaken in a hyperendemic onchocerciasis focus with no history of ivermectin distribution or vector control and located in the Mbam River Valley, Central Region, Cameroon. The trial design is represented schematically in Figure 1; details are given in Gardon et al [18]. Briefly, 657 consenting males (to avoid contraindications associated with ivermectin in pregnancy), aged 18 to 60 years and with at least 2 palpable onchocercomas (subcutaneous nodules), were randomly assigned to 1 of 4 ivermectin treatment groups (Table 1). A single, randomly selected nodule was surgically excised from each participant on up to 3 occasions: prior to ivermectin treatment in May 1994; 3 years after the first trial dose in August 1997; and in either June 1998 or November 1998, 6 or 12 months after the last ivermectin treatment (Figure 1, Table 1). Following guidelines provided by Gardon et al [18], female worms within excised nodules were histologically classified (see Supplementary Materials) as nonfertile, *N*: live, nonfertile worms including potentially fertile females not currently inseminated, producing oocytes which transform into unfertilized ova and then degenerate as well as empty, senescent worms no longer producing oocytes; fertile, *F*: inseminated or re-inseminated and producing embryos of any stage up to live microfilariae; and moribund or dead, *D*: with few or no internal organs, collapsed body wall, and partial or total calcification. Data were analyzed per protocol rather than intent-to-treat; individuals who missed a treatment were systematically excluded from the cohort. We did not analyze data on male worms because their smaller size hinders their detection in histological sections and because very few were identified as dead or moribund [18]. This is common [20], probably because males tend to leave nodules to seek females [21] and may die outside of nodules.

## Modeling Approach and Statistical Inference

We used the treatment submodel embedded within our transmission model EPIONCHO [22, 23], originally developed to describe the population dynamics of adult female *O. volvulus* exposed to a single standard dose (150 µg/kg) of ivermectin [13]. We extended this model to track mean numbers of live and dead worms per host, in addition to the fertility (fertile or nonfertile) status of live worms, through different multiple-dose ivermectin treatment regimens. It was assumed that

<table><caption>**Table 1. Summary of Trialed Ivermectin Treatment Regimens**</caption><thead><tr><td>Group</td><td>Dose</td><td>Frequency</td><td>Description<sup>a</sup></td></tr></thead><tbody><tr><td>1, Standard</td><td>150 µg/kg</td><td>Annual</td><td>3 × 150 µg/kg annually for 3 years</td></tr><tr><td>2, High dose</td><td>400–800 µg/kg</td><td>Annual</td><td>1 × 400 µg/kg and then 2 × 800 µg/kg for 2 years</td></tr><tr><td>3, High frequency</td><td>150 µg/kg</td><td>3-monthly</td><td>12 × 150 µg/kg 3-monthly for 3 years</td></tr><tr><td>4, High dose, high frequency</td><td>400–800 µg/kg</td><td>3-monthly</td><td>2 × 400 µg/kg and then 10 × 800 µg/kg for 2.5 years</td></tr></tbody></table>

<footnote>ªAll individuals were additionally treated with a clearing dose of 150 µg/kg 3 months prior to the start of the study (immediately after the baseline nodulectomy) and with a final dose of 150 µg/kg 6 months after the final treatment.</footnote>
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Schematic representation of the study design.</figcaption></figure>

the (prior) life expectancy of (untreated) adult female *O. volvulus* was 10 years, with 95% of worms dead by age 14 years [14]. Mathematical details are given in the Supplementary Materials, a schematic in Supplementary Figure 5 and a list of parameter definitions in Supplementary Table 2.

The framework captures macrofilaricidal and sterilizing (temporary and/or permanent) effects, referred to collectively as antifilarial effects, on worms differentially exposed to ivermectin via parameters that increase the relevant population dynamics rates above that of unexposed worms (eg, an excess mortality rate; Supplementary Figure 6). These effects can vary with dose and frequency and can change dynamically

such that the cumulative effects on viability or fertility can accrue in a constant decreasing or increasing manner with incremental exposures. Because of ongoing transmission, previously unexposed worms could have been acquired at any time during the trial. Therefore, we modeled explicitly the antifilarial effects on differentially exposed worms (Supplementary Figure 6) [13].

We modeled the temporally dynamic proportion of live female worms and the proportion of fertile (live) worms within hosts given each treatment regimen. By introducing host-specific random effects, we reconciled the modeled group-based dynamics with participant-specific longitudinal data, producing a
<table><caption>**Table 2. Models and Variants Incorporating Different Antifilarial Effects and Different Assumptions on the Transmission Dynamics During the Trial**</caption><thead><tr><th>Model</th><th>Permanent Antifilarial Activity</th><th>Variant</th><th>Assumptions on Short-term and Long-term Changes in Transmission</th></tr></thead><tbody><tr><td>0</td><td><ul><li>No macrofilaricidal activity</li><li>No sterilizing activity</li></ul></td><td>A</td><td><ul><li>Constant short-term transmission</li><li>Constant long-term transmission</li></ul></td></tr><tr><td>1</td><td><ul><li>No macrofilaricidal activity</li><li>Sterilizing activity</li></ul></td><td>B</td><td><ul><li>Seasonal short-term transmission</li><li>Constant long-term transmission</li></ul></td></tr><tr><td>2</td><td><ul><li>Macrofilaricidal activity</li><li>No sterilizing activity</li></ul></td><td>C</td><td><ul><li>Constant short-term transmission</li><li>Declining long-term transmission</li></ul></td></tr><tr><td>3</td><td><ul><li>Macrofilaricidal activity</li><li>Sterilizing activity</li></ul></td><td>D</td><td><ul><li>Seasonal short-term transmission</li><li>Declining long-term transmission</li></ul></td></tr></tbody></table>

Nota bene: All variants A to D of models 0 to 3 were fitted to the data separately (ie, 16 separate fits in total).

statistical model to account for correlation among worm counts made repeatedly on the same individual, permitting (random) nonspecific variation among participants. We modeled state probabilities (proportions of live and fertile worms) to nullify the effects of variation in the incidence of new infections among participants (see Supplementary Materials for details). We did not use data on skin microfilarial densities [18] to avoid invoking additional assumptions on the production of microfilariae by female worms and on the competing dynamics of ivermectin's microfilaricidal and sterilizing effects [13]. Hence, because the focus here is on antimacrofilarial activity, we used data on macrofilariae only.

We considered 4 models, model 0 to model 3, covering all combinations of antifilarial effects, with model 0 being the null model of no macrofilaricidal or permanent sterilizing activity of ivermectin but including a temporary sterilizing effect [13] (Table 2). We assumed that the first exposure to ivermectin does not exert any macrofilaricidal or permanent sterilizing effect [13], and we did not consider postulated prophylactic effects (Supplementary, Prophylactic effects of ivermectin). We fitted these models to the data in a Bayesian framework using Markov chain Monte Carlo techniques (see Supplementary Materials for details). We reflected a priori uncertainty in parameter values using prior distributions (priors). Parameters with information available in the literature were assigned informative priors, reflecting published uncertainties. These include the parameters that define the temporary sterilizing effect of a single standard dose of ivermectin [13], permitting permanent sterilizing effects to be resolved. Parameters without published information were assigned uninformative priors and include all parameters that govern the macrofilaricidal and permanent sterilizing activity (Supplementary Table 2).

### Transmission Intensity

Without data on transmission intensity during the trial (eg, numbers of infective larvae per blackfly) [18], we tested the robustness of the models to different assumptions on short-term (within year) and long-term (between year) changes in transmission. We fitted 4 variants (A to D) of models 0 to 3 (Table 2) covering combinations of short-term seasonal

variation in transmission (mediated by fluctuations in the abundance of blackflies) and a long-term decline in transmission, perhaps induced by the trial (although this was unlikely; see Supplementary Materials for discussion). To capture seasonal variation, we fitted a sinusoidal function to data on blackfly numbers collected in the vicinity just before the trial (see Supplementary Materials for details). To capture long-term reductions in transmission, we modeled an exponential decline such that during the last year of the trial, the intensity of transmission was approximately half of its initial value (see Supplementary Materials and Supplementary Figure 7).

## RESULTS

### Antifilarial Indications and Population Dynamics

The best fitting models under all assumptions of seasonal variation and long-term declines in transmission include both macrofilaricidal activity and permanent sterilization (Supplementary Table 3), indicating that both processes operate. Model 3A had the lowest deviance information criterion (DIC) [24], indicating adequacy of a constant transmission assumption (Table 2). Parameter estimates are only modestly affected by different assumptions on seasonal variation and long-term declines in transmission (Figure 2; Supplementary Table 4).

Figure 3 shows the dynamics, estimated from model 3A, in the proportions of live and fertile female *O. volvulus* through the 4 multiple-dose treatment regimens (Table 1). Model 3A includes macrofilaricidal activity, permanent sterilization, and constant (seasonal and long-term) transmission (parameter estimates are depicted in Supplementary Figure 8). Depicted are both the mean dynamics among all trial participants given a particular regimen and the individual (host) dynamics. The observed (grouped) data are shown for visual appraisal of model fit.

### Macrofilaricidal Activity

The decline in the proportion of live worms (Figure 3A and 3B) is protracted and driven by the permanent excess (additional) mortality of live worms exposed to ivermectin compared to unexposed worms (Figure 2A). Macrofilaricidal activity is not
<figure><img src="image_3.png" /><figcaption>**Figure 2.** Summary of estimated parameters that define the antifilarial activity of multiple doses of ivermectin under different assumptions on the transmission dynamics during the trial. In each panel, the data points and vertical lines indicate the means and 95% Bayesian credible intervals (BCIs) of the posterior distributions estimated from variants A (constant short-term and long-term transmission), B (seasonal short-term transmission; constant long-term transmission), C (constant short-term transmission; declining long-term transmission), and D (seasonal short-term transmission; declining long-term transmission) of the best fitting model 3 (see Supplementary Table 3 for deviance information criteria). The horizontal dotted lines in each panel indicate the null (zero) effect of each parameter. Hence, parameters with BCIs that cross the dotted line can be interpreted as not statistically significantly different from zero. For example, in panels (A) and (D) none of the BCIs include 0, indicating statistically significant (permanent) macrofilaricidal and sterilizing activity of multiple-dose ivermectin regimens. By contrast, in panels (B) and (E) all BCIs cross 0, indicating no statistically significant effect of dose on either macrofilaricidal or permanent sterilizing activity. The parameter posteriors estimated from the model variants are generally similar and therefore robust to different assumptions on the transmission dynamics during the trial.</figcaption></figure>

significantly associated with dose (150 µg/kg vs 800 µg/kg; **Figure 2B**). The pace at which the macrofilaricidal effect accumulates declines with incremental exposures (**Figure 2C** negative incremental exposure adjustment for all model variants).

Figure 4 shows how the estimated life expectancy of female *O. volvulus* is reduced with increasing exposure to ivermectin (see Supplementary Materials for details of this calculation). Worms that establish in participants before the first treatment with ivermectin (and therefore are maximally exposed) have a life expectancy of approximately 4.8 years for the standard-dose annual regimen and 3.3 years for the 3-monthly regimen. This corresponds to 50% and 67% reductions compared to hypothetically unexposed worms, which have an estimated life expectancy of 9.6 years based on the life expectancy of *O. volvulus* in West African savannah [14].

## Sterilizing Activity

The modeled transient effects of the temporary sterilizing effect (Supplementary Figure 6) drives the oscillating proportion of fertile (live) worms (**Figure 3C** and **3D**). The more protracted systematic decline is caused predominantly by the permanent sterilization of worms that are exposed to ivermectin on multiple occasions. Although this effect is estimated with substantial uncertainty (**Figure 2D**, models that included a permanent sterilizing effect (**Table 2**) were superior (according to the DIC)

to those that omitted it (Supplementary Table 3). There is no statistically significant effect of dose on permanent sterilization (**Figure 2E**), but there is a nonsignificant tendency for the higher dose regimens to elicit a stronger effect. The pace at which the permanent sterilizing effect accumulates remains constant with incremental exposures (**Figure 2F**).

## DISCUSSION

Multiple doses of ivermectin have a partial macrofilaricidal effect and a modest permanent and irreversible sterilizing effect on female *O. volvulus* after 4 or more consecutive treatments, even at the annual frequency and 150 µg/kg dose used for MDA. A high 800-µg/kg dose does not significantly modify either effect. We used the scaffold of our newly developed modeling framework [17], linking longitudinal data collected at different follow-up times to the effects of ivermectin on fundamental parasite population processes. This was done in order to estimate the pronounced macrofilaricidal effect of 3 years of annual or 3-monthly multiple-dose regimens as, respectively, a 50% and 67% reduction in the life expectancy of female *O. volvulus*.

In their original analysis of these data, Gardon et al [18] identified a statistically significant higher proportion of dead female worms in patients who received the 3-monthly treatment compared to those treated annually. Raised proportions
<figure><img src="image_4.png" /><figcaption>**Figure 3.** Dynamics of *Onchocerca volvulus* viability and fertility through multiple rounds of ivermectin treatment. Each panel depicts either the percentage of live female *O. volvulus* (panels A and B) or the percentage of fertile (live) female *O. volvulus* (panels C and D) exposed to ivermectin treatments given annually (panels A and C) or 3-monthly (panels B and D) at standard (150 µg/kg) and high doses (800 µg/kg). The thick gray and black solid lines indicate, respectively, the marginal mean (averaged across all hosts) dynamics for standard or high-dose regimens. The thin gray lines indicate the individual host dynamics mediated by the estimated random effects terms within the statistical model (see Supplementary Materials for statistical details). For clarity, no graphical distinction is given among individual hosts given either standard or high-dose regimens. In each panel, the triangular and circular markers denote, respectively, the observed percentages of live or fertile (live) females extracted from participants given either standard or high-dose regimens. Vertical error bars denote 95% confidence intervals that were calculated using a nonparametric weighted bootstrapping procedure to account for the variable number of female *O. volvulus* extracted from each host. Horizontal error bars (often narrower than the data point) indicate the range of nodulectomy dates at each sampling round.</figcaption></figure>

of dead worms from treated compared to untreated controls were also reported from trials in Liberia, Guatemala, and Sierra Leone (see Supplementary Materials and Supplementary Table 1 for a summary of these trials). The data in these studies were analyzed using traditional statistical techniques, making it impossible to translate proportions of dead worms (observed at variable follow-up times) into estimates of worm life expectancy.

Moreover, with the exception of 9 treated participants from a trial in Guatemala [25], (Supplementary Table 1), none of the previous trials followed participants longitudinally.

Similarly, the original conclusions on the sterilizing activity of the multiple-dose regimens were based on comparisons of proportions of fertile worms among differently treated participants [18]. A lower proportion among 3-monthly compared
<figure><img src="image_5.png" /><figcaption>**Figure 4.** The life expectancy of adult *Onchocerca volvulus* multiply exposed to ivermectin. The dashed and dotted lines in panel A show the estimated life expectancy (in years) of adult female *O. volvulus* exposed to annual and 3-monthly treatments with ivermectin at a standard dose (150 µg/kg) (there is no substantive difference in the corresponding estimates for a high dose of 800 µg/kg dose, Figures 2 and 3). The solid line in this panel (A) depicts the estimated life expectancy of hypothetically unexposed worms. In panel (B), the dashed and dotted lines express the estimated life expectancies as percentage reductions compared to hypothetically unexposed worms. Estimates are plotted against the date of infection, that is, the time of establishment of female worms. Treatment dates with ivermectin are represented by the triangles: 3-monthly treatments on the top row and annual treatments below. Hence, worms that establish at an earlier date are exposed to more treatments and incur greater reductions in life expectancy. In each panel, gray shaded areas represent 95% Bayesian credible intervals. Life expectancies were calculated from the estimated parameter posteriors of model 3A (the best fitting model, see Supplementary Table 3 for deviance information criteria). Details of these calculations are given in the Supplementary Materials.</figcaption></figure>

to annually treated participants was reported; however, it was not possible to determine whether this was consistent with the (since well-described [13]) temporary sterilizing effect of ivermectin or indicative of irreversible impairment (a conundrum discussed elsewhere [26]). Our modeling approach resolves this and suggests a modest, permanent sterilization effect, albeit bounded by high uncertainty because both the macrofilaricidal and sterilizing effects influence observed proportions of fertile worms. The macrofilaricidal activity of ivermectin kills an increasing proportion of nonfertile worms before they become fertile for the first time, altering the worm population structure to comprise relatively more juvenile worms.

We recently used our modeling framework to quantify the macrofilaricidal activity of doxycycline, a tetracycline antibiotic that depletes the endosymbiotic *Wolbachia* bacteria essential to *O. volvulus* viability and fertility [17]. We reported an 80% reduction in the life expectancy of *O. volvulus* treated with 4–6 weeks of 100 mg or 200 mg daily doxycycline, which is higher than the macrofilaricidal efficacy of multiple doses of ivermectin estimated here. Other macrofilaricidal therapies include intravenous suramin, which is poorly tolerated with significant side effects and is therefore seldom used, and flubendazole, which was originally associated with serious intramuscular injection-site reactions [27] and subsequent development of an oral reformulation has been recently discontinued over toxicology concerns [28]. Hence, doxycycline is the de facto gold standard macrofilaricidal treatment for onchocerciasis, despite possible challenges associated with adherence to a relatively long treatment course.

This paucity of macrofilaricides emphasizes the need for new drugs to improve treatment of onchocerciasis in clinical settings and to bolster the armory of treatment options to assist

with global control and elimination efforts. There is substantial and ongoing research into new drugs and drug combinations [29]. The analysis presented here, like that used to evaluate the efficacy of doxycycline, shows that the interpretation of clinical trial data is facilitated by population dynamics modeling techniques. These approaches can be used in a retrospective, analytical manner (as presented here) and prospectively to inform clinical trial design and the formulation of so-called target product profiles by simulating (parasitological or other) outcomes using postulated antifilarial activities.

The specific structure of the model used here is based on the treatment submodel embedded within EPIONCHO [22, 23], a transmission model that is used to explore the effectiveness of current and alternative intervention strategies that target onchocerciasis elimination in Africa [23, 30, 31]. In the future, to give our efficacy estimates more explicit programmatic interpretation, we will explore how partial macrofilaricidal and sterilizing effects of repeated rounds of ivermectin given during MDA alter current projections on the feasibility of achieving elimination targets [1, 32]. This will include revising projections on annual and 6-monthly treatment [30] and modeling the potential benefit of 3-monthly MDA regimens; the latter showing a significantly higher macrofilaricidal efficacy than annual MDA. Irrespective of our new understanding of the partial macrofilaricidal activity of repeated doses of ivermectin, maintaining high coverage and adherence to treatment, particularly in highly endemic areas, will remain extremely important to the effectiveness of MDA [22, 30].

Some epidemiological projections that incorporate an instantaneous macrofilaricidal effect of ivermectin have been made
using the ONCHOSIM transmission model [22, 33]. These analyses attribute long-term reductions in microfilarial densities to either macrofilaricidal or permanent sterilizing activity, but not to both mechanisms acting simultaneously. Furthermore, their estimates of macrofilaricidal efficacy, based on clinical trial data from Guatemala [34] of 3-monthly treatments (Supplementary Table 1), rested on an assumption of no ongoing transmission during the trial and no resorption of dead worms [17, 35], risking underestimation of the true macrofilaricidal effect. Our approach untangles the permanent sterilizing and macrofilaricidal effects for both annual and 3-monthly ivermectin regimens at different doses and adjusts for the diluting effects of ongoing transmission.

Notwithstanding the partial macrofilaricidal activity of ivermectin, the search for new drugs to assist with the control and elimination of onchocerciasis remains extremely important [29]. Ivermectin cannot be used safely in loiasis-onchocerciasis co-endemic areas of West and Central Africa because of the risk of severe adverse events associated with high-grade *Loa loa* microfilaremia [36, 37]. Many communities in these areas are currently not served with effective onchocerciasis control. In such areas, macrofilaricidal (and permanently sterilizing) anti-*Wolbachia* therapies that do not affect *L. loa* (devoid of *Wolbachia*) are particularly promising. Anti-*Wolbachia* therapies may also impede parasite development in the blackfly vector [38], exerting a potential transmission-blocking effect. Elsewhere in Africa, reports of suboptimal responses to ivermectin continue [39] and should be considered an impetus to seeking alternative drugs to mitigate potentially emerging ivermectin resistance. Moxidectin is a promising candidate for MDA in these circumstances [40].

## CONCLUSIONS

We modeled longitudinal data from the most comprehensive clinical trial of multiple doses of ivermectin treatment using a powerful and contemporary population dynamics and statistical framework. We show that 4 or more consecutive treatments with ivermectin are partially macrofilaricidal and exert a permanent sterilizing effect on female *O. volvulus*, even at the standard dose and frequency used in routine MDA in Africa. It is essential to incorporate these effects into mathematical transmission models to reevaluate existing projections on the feasibility of eliminating onchocerciasis.

## Supplementary Data

Supplementary materials are available at *The Journal of Infectious Diseases* online. Consisting of data provided by the authors to benefit the reader, the posted materials are not copyedited and are the sole responsibility of the authors, so questions or comments should be addressed to the corresponding author.

## Notes

**Acknowledgments.** We thank J.A.N. Filipe for a preliminary modeling analysis of the data.

**Disclaimer.** The funders had no influence on study design, collection, analysis, and interpretation of the data or the decision to submit the paper for publication.

**Financial support.** This work was supported by the Wellcome Trust (grant 092677/Z/10/Z), the Neglected Tropical Disease Modelling Consortium (supported by the Bill & Melinda Gates Foundation in partnership with the Task Force for Global Health), and the United Nations Children's Fund/United Nations Development Programme/World Bank World Health Organization Special Programme for Research and Training in Tropical Diseases. The data were originally collected using funds from the River Blindness Foundation.

**Potential conflicts of interest.** All authors: No reported conflicts of interest. All authors have submitted the ICMJE Form for Disclosure of Potential Conflicts of Interest. Conflicts that the editors consider relevant to the content of the manuscript have been disclosed.

## References

1. World Health Organization. Accelerating work to overcome the global impact of neglected tropical diseases: a roadmap for implementation. Geneva: World Health Organization, 2012. Available at: http://www.who.int/neglected_diseases/NTD_RoadMap_2012_Fullversion.pdf. Accessed 7 March 2017.

2. Botto C, Basañez MG, Escalona M, et al. Evidence of suppression of onchocerciasis transmission in the Venezuelan Amazonian focus. Parasit Vectors **2016**, 9:40.

3. Progress towards eliminating onchocerciasis in the WHO Region of the Americas: verification of elimination of transmission in Mexico. Wkly Epidemiol Rec, **2015**, 90: 577–81.

4. Traore MO, Sarr MD, Badji A, et al. Proof-of-principle of onchocerciasis elimination with ivermectin treatment in endemic foci in Africa: final results of a study in Mali and Senegal. PLoS Negl Trop Dis **2012**, 6:e1825.

5. Tekle AH, Elhassan E, Isiyaku S, et al. Impact of long-term treatment of onchocerciasis with ivermectin in Kaduna State, Nigeria: first evidence of the potential for elimination in the operational area of the African Programme for Onchocerciasis Control. Parasit Vectors **2012**, 5:28.

6. Tekle A, Zouré H, Noma M, et al. Progress towards elimination in the participating countries of African Programme for Onchocerciasis Control: epidemiological evaluation results. Infect Dis Pov, **2016**, 5:66.

7. Katabarwa MN, Eyamba A, Nwane P, et al. Seventeen years of annual distribution of ivermectin has not interrupted onchocerciasis transmission in North Region, Cameroon. Am J Trop Med Hyg **2011**, 85:1041–9.

8. Katabarwa MN, Lakwo T, Habomugisha P, et al. Transmission of *Onchocerca volvulus* continues in Nyagak-Bondo focus of northwestern Uganda after 18 years of a single dose of annual treatment with ivermectin. Am J Trop Med Hyg **2013**, 89:293–300.

9. Wanji S, Kengne-Ouafo JA, Esum ME, et al. Situation analysis of parasitological and entomological indices of onchocerciasis transmission in three drainage basins of the rain forest of South West Cameroon after a decade of ivermectin treatment. Parasit Vectors **2015**, 8:202.

10. Eisenbarth A, Achukwi MD, Renz A. Ongoing transmission of *Onchocerca volvulus* after 25 years of annual ivermectin mass treatments in the Vina du Nord River Valley, in North Cameroon. PLoS Negl Trop Dis **2016**, 10:e0004392.

11. African Programme for Onchocerciasis Control. Report of the consultative meetings on strategic options and alternative treatment strategies for accelerating onchocerciasis elimination in Africa. Geneva: World Health Organization, 2015. Available at: http://www.who.int/apoc/ATS_Report_2015.12.pdf. Accessed 7 March 2017.

12. Turner HC, Churcher TS, Walker M, Osei-Atweneboana MY, Prichard RK, Basáñez MG. Uncertainty surrounding projections of the long-term impact of ivermectin treatment on human onchocerciasis. PLoS Negl Trop Dis **2013**, 7:e2169.

13. Basáñez MG, Pion SD, Boakes E, Filipe JA, Churcher TS, Boussinesq M. Effect of single-dose ivermectin on *Onchocerca volvulus*: a systematic review and meta-analysis. Lancet Infect Dis **2008**, 8:310–22.

14. Plaisier AP, van Oortmarssen GJ, Remme J, Habbema JD. The reproductive lifespan of *Onchocerca volvulus* in West African savanna. Acta Trop **1991**, 48:271–84.

15. Plaisier AP, Alley ES, Boatin BA, 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.
16. Bottomley C, Isham V, Collins RC, Basáñez MG. Rates of microfilarial produc-
tion by *Onchocerca volvulus* are not cumulatively reduced by multiple ivermectin
treatments. *Parasitology* **2008**, 135:1571–81.

17. Walker M, Specht S, Churcher TS, Hoerauf A, Taylor MJ, Basáñez MG. Therapeutic efficacy and macrofilaricidal activity of doxycycline for the treatment of river blindness. *Clin Infect Dis* **2015**, 60:1199–207.

18. Gardon J, Boussinesq M, Kamgno J, Gardon-Wendel N, Demanga-Ngangue, Duke BO. Effects of standard and high doses of ivermectin on adult worms of *Onchocerca volvulus*: a randomised controlled trial. *Lancet* **2002**, 360:203–10.

19. Cupp EW, Cupp MS. Short report: impact of ivermectin community-level treat-
ments on elimination of adult *Onchocerca volvulus* when individuals receive mul-
tiple treatments per year. Am J Trop Med Hyg **2005**, 73:1159–61.

20. Karam M, Schulz-Key H, Remme J. Population dynamics of *Onchocerca volvulus* after 7 to 8 years of vector control in West Africa. *Acta Trop* **1987**, 44:445–57.

21. Schulz-Key H, Jean B, Albiez EJ. Investigations on female *Onchocerca volvulus* for the evaluation of drug trials. *Tropenmed Parasitol* **1980**, 31:34–40.

22. Basáñez MG, Walker M, Turner HC, Coffeng LE, de Vlas SJ, Stolk WA. River blindness: mathematical models for control and elimination. *Adv Parasitol* **2016**, 94:247–341.

23. Walker M, Stolk WA, Dixon MA, et al. Modelling the elimination of river blindness using long-term epidemiological and programmatic data from Mali and Senegal. *Epidemics* **2017**, 18:4–15.

24. Spiegelhalter DJ, Best NG, Carlin BP, van der Linde A. Bayesian measures of com-
plexity and fit. *J R Statist Soc B* **2002**, 64:583–639.

25. Duke BOL, Zea-Flores G, Castro J, Cupp EW, Muñoz B. Comparison of the effects of a single dose and of 4 6-monthly doses of ivermectin on adult *Onchocerca volvulus*. Am J Trop Med Hyg **1991**, 45:132–7.

26. Kläger SL, Whitworth JA, Downham MD. Viability and fertility of adult *Onchocerca volvulus* after 6 years of treatment with ivermectin. *Trop Med Int Health* **1996**, 1:581–9.

27. Mackenzie CD, Geary TG. Flubendazole: a candidate macrofilaricide for lymph-
atic filariasis and onchocerciasis field programs. Expert Rev Anti Infect Ther
**2011**; 9:497–501.

28. Johnson & Johnson. Janssen discontinues development of flubendazole formu-
lation to treat onchocerciasis. Available at: https://www.jnj.com/media-center/
press-releases/janssen-discontinues-development-of-flubendazole-formula-
tion-to-treat-onchocerciasis. Accessed 19 June 2017.

29. Kuesel AC. Research for new drugs for elimination of onchocerciasis in Africa. *Int J Parasitol* Drugs Drug Resist **2016**, 6:272–86.

30. Stolk WA, Walker M, Coffeng LE, Basáñez MG, de Vlas SJ. Required duration of mass ivermectin treatment for onchocerciasis elimination in Africa: a comparative modelling analysis. *Parasit Vectors* **2015**, 8:552.

31. Turner HC, Walker M, Churcher TS, et al. Reaching the London Declaration on Neglected Tropical Diseases goals for onchocerciasis: an economic evaluation of increasing the frequency of ivermectin treatment in Africa. *Clin Infect Dis* **2014**, 59:923–32.

32. African Programme for Onchocerciasis Control. 2012 Eighteenth Session of the Joint Action Forum. Available at: http://www.who.int/apoc/about/structure/jaf/FinalCommunique_JAF_18_English_final_with_annexes.pdf. Accessed 7 March 2017.

33. Coffeng LE, Stolk WA, Hoerauf A, et al. Elimination of African onchocerciasis: modeling the impact of increasing the frequency of ivermectin mass treatment. *PLoS One* **2014**, 9:e115886

34. Duke BOL, Zea-Flores G, Castro J, Cupp EW, Muñoz B. Effects of 3-month doses of ivermectin on adult *Onchocerca volvulus*. Am J Trop Med Hyg **1992**, 46:189–94.

35. Albeiz EJ. Calcification in *Onchocerca volvulus*. Trop Med Parasitol **1985**; 36:180-1.

36. Gardon J, Gardon-Wendel N, Demanga-Ngangue, Kamgno J, Chippaux JP, Boussinesq M. Serious reactions after mass treatment of onchocerciasis with ivermectin in an area endemic for *Loa loa* infection. *Lancet* **1997**, 350: 18–22.

37. Zouré HG, Wanji S, Noma M, et al. The geographic distribution of *Loa loa* in Africa: results of large-scale implementation of the rapid assessment procedure for loiasis (RAPLOA). *PLoS Negl Trop Dis* **2011**, 5:e1210.

38. Albers A, Esum ME, Tendongfor N, et al. Retarded *Onchocerca volvulus* L1 to L3 larval development in *Simulium damnosum* vector after anti-wolbachial treatment in human host. *Parasit Vectors* **2012**, 5:12.

39. Frempong KK, Walker M, Cheke RA, et al. Does increasing treatment frequency address suboptimal responses to ivermectin for the control and elimination of river blindness? *Clin Infect Dis* **2016**, 62:1338–47.

40. Turner HC, Walker M, Attah SK, et al. The potential impact of moxidectin on onchocerciasis elimination in Africa: an economic evaluation based on the phase II clinical trial data. *Parasit Vectors* **2015**; 8:167.