Ivermectin susceptibility and sporontocidal effect in Greater Mekong Subregion Anopheles

Novel vector control methods that can directly target outdoor malaria transmission are urgently needed in the Greater Mekong Subregion (GMS) to accelerate malaria elimination and artemisinin resistance containment efforts. Ivermectin mass drug administration (MDA) to humans has been shown to effectively kill wild Anopheles and suppress malaria transmission in West Africa. Preliminary laboratory investigations were performed to determine ivermectin susceptibility and sporontocidal effect in GMS Anopheles malaria vectors coupled with pharmacokinetic models of ivermectin at escalating doses. A population-based pharmacokinetic model of ivermectin was developed using pre-existing data from a clinical trial conducted in Thai volunteers at the 200 µg/kg dose. To assess ivermectin susceptibility, various concentrations of ivermectin compound were mixed in human blood meals and blood-fed to Anopheles dirus, Anopheles minimus, Anopheles sawadwongporni, and Anopheles campestris. Mosquito survival was monitored daily for 7 days and a non-linear mixed effects model with probit analyses was used to calculate concentrations of ivermectin that killed 50% (LC50) of mosquitoes for each species. Blood samples were collected from Plasmodium vivax positive patients and offered to mosquitoes with or without ivermectin at the ivermectin LC25 or LC5 for An. dirus and An. minimus. The GMS Anopheles displayed a range of susceptibility to ivermectin with species listed from most to least susceptible being An. minimus (LC50 = 16.3 ng/ml) > An. campestris (LC50 = 26.4 ng/ml) = An. sawadwongporni (LC50 = 26.9 ng/ml) > An. dirus (LC50 = 55.6 ng/ml). Mosquito survivorship results, the pharmacokinetic model, and extensive safety data indicated that ivermectin 400 µg/kg is the ideal minimal dose for MDA in the GMS for malaria parasite transmission control. Ivermectin compound was sporontocidal to P. vivax in both An. dirus and An. minimus at the LC25 and LC5 concentrations. Ivermectin is lethal to dominant GMS Anopheles malaria vectors and inhibits sporogony of P. vivax at safe human relevant concentrations. The data suggest that ivermectin MDA has potential in the GMS as a vector and transmission blocking control tool to aid malaria elimination efforts.


Background
Artemisinin-resistant Plasmodium falciparum isolates have been identified in the Greater Mekong Subregion (GMS) which threatens to undermine malaria elimination goals [1]. The primary GMS Anopheles vectors of Plasmodium frequently feed outdoors and before people go to sleep [2,3], rendering classical vector control measures such as insecticide-treated bed nets and indoor residual spraying with insecticides less effective. There are currently no vector control tools in the GMS that specifically target outdoor feeding Anopheles which is a major impediment for malaria vector control in the region. Novel vector control tools that target outdoor malaria transmission are urgently needed to support artemisinin resistance containment and malaria elimination efforts in the GMS.
Ivermectin has been shown to be lethal to over a dozen Anopheles species worldwide [4]. Ivermectin mass drug administration (MDA) has been suggested as a possible malaria parasite transmission control tool as it directly targets the vector at the point of human blood feeding. It is one of the few vector control measures that targets outdoor malaria transmission. Laboratory studies [5,6], clinical trials [7][8][9], and field studies [10,11] demonstrate that ivermectin is lethal to Anopheles gambiae at human relevant concentrations. Ivermectin MDA campaigns in Senegal, Liberia, and Burkina Faso have demonstrated that ivermectin can suppress P. falciparum transmission by wild An. gambiae s.l. [11,12]. In addition to direct mosquito-lethal effects, ivermectin suppresses development of P. falciparum in An. gambiae [13,14].
Plasmodium transmission in the GMS is complex with numerous Anopheles species serving as primary and secondary vectors throughout the region. Anopheles dirus s.s., found in the GMS east of Myanmar, is a primary malaria vector in forested areas and feeds predominantly outdoors on humans [2,15]. Anopheles minimus s.s., found throughout the GMS and parts of mainland Asia, is a primary malaria vector outside of forested areas, feeds both indoors and outdoors, and displays a variable feeding preference on humans and cattle across its range [3,16]. There are several secondary malaria vectors that may facilitate malaria transmission in the GMS. Anopheles sawadwongporni, a member of the Anopheles maculatus group, has been incriminated as a malaria vector in Thailand [17] and Anopheles campestris, predominantly found in rice paddies throughout the GMS is suspected to be a secondary malaria vector in Thailand [18].
Ivermectin MDA campaigns in West Africa using a 150-200 µg/kg dose were shown to be effective at reducing P. falciparum transmission by An. gambiae [11,12]. However, not all Anopheles may be equally susceptible to ivermectin [4]. This implies that higher doses of ivermectin may be required during MDA to effectively target all Anopheles in a given region. Doses of ivermectin up to 2000 µg/kg were safe and well tolerated in healthy volunteers [19], which is ten times the amount approved for strongyloidiasis treatment in Thailand [20]. The ivermectin dose of 800 µg/kg has been assessed in onchocerciasis-infected patients in an extended trial in Ghana [21,22] and repeatedly every 3 months in a trial in Cameroon [23]. Adverse events in the ivermectin 800 µg/kg trials may be correlated with the immune response to dead microfilariae and not necessarily linked directly to ivermectin treatment. Trials in healthy and malaria infected patients without concomitant onchocerciasis infection at the ivermectin 800 μg/kg dose are warranted. The lack of adverse events in healthy volunteers at up to ivermectin doses of 2000 μg/kg [19] would support the notion that the adverse events observed in Ghana [21,22] and Cameroon [23] trials may be linked to the death and clearance of Onchocerca volvulus parasites and not ivermectin toxicity. Several clinical trials have investigated the safety and tolerability of ivermectin at 400 and 800 µg/ kg, however, to the best of our knowledge, none have assessed the pharmacokinetic properties of ivermectin at these doses, particularly in an Asian population, which is the population of interest in this study.
The mosquito-lethal effects of ivermectin on the GMS malaria vectors An. dirus s.s., An. minimus s.s., An. sawadwongporni, and An. campestris, and the sporontocidal effects of ivermectin on Plasmodium vivax in An. dirus s.s., and An. minimus s.s. were investigated. A population-based pharmacokinetic model of ivermectin and simulated concentration-time profiles after 200, 400 and 800 µg/kg doses was developed, and correlated with mosquito survivorship results to rationally select ivermectin doses for MDA use in the GMS.

Pharmacokinetic modelling
Frequent ivermectin plasma concentrations were collected in 23 healthy Thai volunteers after a standard 200 µg/kg dose [24]. Ivermectin concentrations were transformed into their natural logarithms and characterized using nonlinear mixed-effects modelling in NON-MEM version.7.3 (Icon Development Solution, Ellicott City, MD, USA). The first-order conditional estimation method with interaction was used throughout model development. Model diagnostics and automation were performed using Xpose version 4.0 [25], Pirana [26] and Pearl-speaks-NONMEM (PsN; version 3.6.0) [27]. The objective function value (OFV; calculated by NONMEM as proportional to −2 × loglikelihood of data) was used to discriminate between hierarchical models. A difference between two models (ΔOFV) of >3.84 and >10.83 were considered statistically significant at a P value of <0.05 and <0.001, respectively. Models were assessed using standard goodness-of-fit and simulation-based diagnostics (i.e. numerical and visual predictive checks using 2000 simulations). The robustness of the final model was evaluated by bootstrapping (n = 1000).
Different disposition models (one-, two-and three-compartment models) and absorption models (first-order and transit models) were evaluated to describe the pharmacokinetic structural model of ivermectin. Pharmacokinetic parameters were assumed to be log-normally distributed and inter-individual variability was implemented as exponential functions. The effect of body size was evaluated using an allometric function (centred on the medium trial weight of 55 kg), scaling clearances with an exponent of ¾ and volumes with an exponent of 1. The influence of other individual covariates including gender, age, haemoglobin, haematocrit, total bilirubin, direct bilirubin, alkaline phosphatase, aspartate aminotransferase, and alanine aminotransferase, serum creatinine, blood urea nitrogen, and albumin level were investigated using a stepwise covariate approach (P value <0.05 and <0.001 for forward and backward addition and deletion, respectively). The developed final population pharmacokinetic model of ivermectin was used to simulate population mean ivermectin pharmacokinetic concentration-time profiles after a single oral dose of 200, 400 and 800 µg/kg using Berkeley Madonna [28].

Mosquitoes
All mosquitoes were reared at the Armed Forces Research Institute of Medical Sciences Department of Entomology in Bangkok, Thailand. Anopheles dirus s.s., An. minimus s.s., An. sawadwongporni, and An. campestris were produced as described previously [29]. Adult mosquitoes used for experiments were provided 10% sucrose solution ad libitum. Mosquitoes were reared at 25 ± 2 °C and 80 ± 10% relative humidity, and 12 h light:12 h dark photoperiod. Mosquitoes were between 5 and 8 days post emergence at time of first blood feed, and mosquitoes were sugar starved with access to water from 12 to 18 h prior to their first blood meal.

Drug
Powdered formulation of ivermectin compound was obtained from Sigma-Aldrich (St. Louis, USA). Ivermectin was dissolved in dimethylsulfoxide (DMSO) to concentrations of 10 mg/ml and 20 µl aliquots were frozen at −20 °C. Ivermectin was thawed and serial dilutions were made in phosphate buffered saline (PBS) and 10 μl was added to 990 μl of blood to reach final concentration desired for mosquito membrane feeding assays. Control blood meals consisted of previously frozen DMSO diluted in PBS to match the ratio of DMSO and PBS fed to mosquitoes in the ivermectin-containing blood meals.

Blood and serum
Pooled whole blood collected in 350-ml citrate phosphate dextrose adenine (CPDA-1) anticoagulant bags (National Blood Centre, Thai Red Cross Society, Bangkok, Thailand) was used to maintain colonies and perform ivermectin lethal concentration experiments. Blood was never more than 2 weeks post collection at the time of mosquito blood feeds. Plasmodium vivax-infected blood was drawn from individual volunteers into 10-ml sodium heparin tubes (NH) (158 USP units, BD Vacutainer, Franklin Lakes, NJ, USA) for the ivermectin sporogony experiments. Serum (type AB+) from malaria-naïve persons (Interstate Blood Bank Inc, Memphis, TN, USA) was processed and stored frozen at −20 °C until ready for use in ivermectin sporogony experiments.

Ivermectin lethal concentration calculations
Anopheles dirus, An. minimus, An. sawadwongporni, and An. campestris were blood fed multiple concentrations of ivermectin to determine the lethal concentration that killed 50% (LC 50 ), 25% (LC 25 ) and 5% (LC 5 ) of the mosquitoes at day 7 post blood meal, following previous methods [5,13]. Control blood meals consisted of DMSO diluted in PBS to match the concentration found in the highest ivermectin treatment group in each replicate. After blood feeding via a membrane feeder, unfed mosquitoes were removed and discarded. Blood-fed mosquitoes were held in 3.5-l plastic containers with access to 10% sucrose and kept in an incubator at 25 ± 1 °C and 70 ± 10% relative humidity, on a 12 h light:12 h dark photoperiod. Mosquito survivorship was monitored for 7 days, every 24 h dead mosquitoes were removed and recorded, and on day 7 all remaining mosquitoes were frozen and counted as alive.
The time above the LC 50 -values for different species of mosquitoes was determined from the simulated pharmacokinetic concentration-time profiles after a single oral dose of 200, 400 and 800 µg/kg, as described above.

Effect of ivermectin on Plasmodium vivax sporogony
Previously it was demonstrated that the ivermectin LC 25 inhibited the development of P. falciparum in An. gambiae [13], and a similar experimental design was repeated here. Mosquitoes reared in Bangkok were transported to a field laboratory at the Mae Sot Malaria Medical Clinic in Mae Sot, Tak Province. Plasmodium vivaxinfected blood was collected from malaria-infected patients reporting to malaria clinics operated by the Thailand Ministry of Public Health in Tak Province (protocol WRAIR#1949A). Patient plasma was removed and replaced by centrifuging blood samples at 3000 rpm for 5 min, plasma removed and blood washed with RPMI twice, and malaria-naïve AB+ serum was added at a 50:50 ratio [30]. One millilitre of blood was fed to each group of 100 mosquitoes in 0.5-l cardboard containers and unfed mosquitoes were removed. Mosquitoes were maintained on 10% sucrose ad libitum at the Mae Sot field insectary (25 ± 1 °C and 70 ± 10% relative humidity, light 12 h:12 h dark photoperiod). Once infected with P. vivax, blood-fed mosquitoes were securely transported back to Bangkok and maintained in a separate insectary (25 ± 1 °C and 70 ± 10% relative humidity, light 12 h:12 h dark photoperiod).
For An. dirus, four experimental feeding schedules performed where ivermectin LC 25 and LC 5 was blood fed concomitantly with P. vivax at day post infection zero (DPI 0), ivermectin LC 5 three days before (DPI −3) P. vivax, and ivermectin LC 25 six (DPI 6) or nine (DPI 9) after P. vivax (Fig. 1). For An. minimus only one experimental feeding schedule was performed when ivermectin LC 25 and LC 5 was blood fed concomitantly with P. vivax (DPI 0). For each mosquito species and feeding schedule there was an ivermectin (LC 25 or LC 5 ) and control group, which had matching concentrations of DMSO and PBS. When mosquitoes were fed ivermectin or control meals at DPI −3, 6 or 9, uninfected donor blood at a 50% haematocrit was used for the blood feed. Unfed mosquitoes were removed from the feeding containers after each blood meal. Mosquitoes from the DPI 0 and DPI −3 feeding schedules were dissected 7 days post parasite ingestion to enumerate oocysts. Midguts were dissected with minuten pins into saline on a microscope slide and stained with 0.4% mercurochrome and viewed at 40× magnification with a compound microscope to determine oocyst prevalence and intensity. Approximately 25 mosquitoes were dissected from each control and treatment group at each dissection time point. Dissectors were blinded to treatment groups. Due to lack of initial mosquito blood feeding and mosquitoes lost to experimental mortality and oocyst dissections at day 7 post parasite feeding, not enough mosquitoes remained at day 14 post parasite feeding to determine sporozoite prevalence for the DPI −3 and 0 experiments. Mosquitoes from the DPI 6 and 9 experiments were dissected 14 days post parasite ingestion to determine sporozoite prevalence. Salivary glands were dissected with minuten pins into saline on a microscope slide and viewed at 40× magnification to determine sporozoite prevalence.

Effect of Plasmodium vivax infection on Anopheles dirus survivorship response to ivermectin
Previously it was shown that An. gambiae infected with P. falciparum were slightly but significantly more susceptible to an ivermectin blood meal at day 14 post parasite ingestion [13]. Anopheles dirus were fed a P. vivax blood meal or a non-infectious blood meal prepared from uninfected donor blood and serum was replaced as described above. Both blood meals were spiked with ivermectin at the LC 25 . After blood feeding, 50 blood fed females from each treatment group were transferred to a new 0.5-l cardboard carton. Mosquito survivorship was monitored daily and dead mosquitoes were removed from the containers until day 14 post blood meal. Any live mosquitoes at day 14 were frozen and counted as alive. A separate group of mosquitoes were fed a P. vivax control blood meal to ensure that the P. vivax sample was infectious to mosquitoes and only isolates that successfully infected at least 80% of An. dirus were included in the results.

Statistical analysis
A non-linear mixed effects model with probit analysis was used to calculate LC 50 , LC 25 and LC 5 values with Statistical Analysis Software (SAS Institute Inc, Cary, NC, USA) as described previously [5]. Hazard ratios for mosquito mortality at day 7 post blood meal were calculated using Poisson regression analysis with STATA version 12.1 (StataCorp LLC, College Station, TX, USA). Oocyst and sporozoite prevalence (i.e. proportion of infected mosquitoes) was compared by Fishers Exact test. Oocyst intensity (i.e. number of oocysts per infected mosquito) was compared by the Mann-Whitney U test. To determine if P. vivax infection altered mosquito survivorship post-ivermectin ingestion survival data were analysed by the Mantel-Cox method. The Fishers Exact, Mann-Whitney U, and survival analysis data were analysed with Prism 7 (GraphPad Software Inc, San Diego, CA, USA).

Pharmacokinetic modelling
A total of 534 venous blood samples were collected from 23 healthy volunteers (12 males and 11 females) for up to 7 days after a single oral dose of 200 µg/kg of ivermectin. There were no observations reported to be below the lower limit of quantification. Observed ivermectin concentrations were best described by a two-compartment disposition model, with no further improvement of an additional disposition compartment (∆OFV = 0). A flexible transit absorption model (n = 2) improved the model fit substantially compared to the first-order absorption (∆OFV = −378). The implementation of body weight as an allometric function had a negligible impact on model fit (∆OFV = 0.516), but it was retained in the final model based on the strong biological prior. Implementing body weight as a covariate in the model also allowed for translation and interpretation to other study populations.
No additional covariates had a significant impact on the model fit. The final model showed satisfactory goodness of fits (Additional file 1: Figure S1) and simulation-based predictive performance (Fig. 2). A numerical predictive check resulted in 2.62% (95% CI 0.37-6.18%) and 0.56% (95% CI 0.56-5.43%) of ivermectin observations below and above the simulated 95% prediction interval, respectively. Ivermectin population pharmacokinetic parameter estimates and secondary parameters from final model are summarized in Table 1. The simulated concentration-time profiles after a single oral dose of 200, 400 and 800 µg/kg of ivermectin are illustrated in Fig. 3.
The GMS Anopheles displayed a range of susceptibility to ivermectin with species listed from most to least susceptible being An. minimus > An. campestris = An. sawadwongporni > An. dirus ( Table 2). The amount of time that ivermectin compound was predicted to be above each Anopheles species 7-day-LC 50 concentration in human plasma for each 200, 400 and 800 µg/kg dose is reported in Table 3. For An. dirus all ivermectin concentrations had a significantly increased hazard of mortality compared to the control group; for An. minimus all ivermectin concentrations had significantly increased hazard of mortality compared to the control group except for 6, 3 and 2 ng/ml; for An. sawadwongporni all ivermectin concentrations had significantly increased hazard of mortality compared to the control group except for 12, 10, 8, 6, 5 and 4 ng/ml; and for An. campestris all ivermectin concentrations had significantly increased hazard of mortality compared to the control group except for 10, 6, 5 and 3 ng/ml (Additional file 1: Table S1).

Discussion
All four Anopheles species investigated were susceptible to ivermectin at concentrations predicted to be present in humans following oral administration. This suggests that ivermectin MDA has a potential role in malaria elimination in the GMS as it is a novel vector control tool that could directly combat outdoor malaria transmission. Anopheles dirus and An. minimus are the two most important primary vectors of malaria in the GMS while An. campestris and An. sawadwongporni are potential secondary vectors [31]. Since An. dirus is arguably the most    (Table 3). Based on the An. dirus in vitro survivorship results, pharmacokinetic model output and extensive safety data discussed below, the ivermectin 400 µg/kg dose appears to be the ideal minimal dose to administer during MDA for malaria control in the GMS. The simulated time above the An. dirus ivermectin 7-day-LC 50 at ivermectin doses of 400 µg/kg (0.4 days) or 800 µg/kg (1.1 days) would not seem to deliver a substantial mosquito-killing window following ivermectin   MDA. However, it should be noted that the 200 µg/ kg dose would only provide a simulated 0.6 days above the An. gambiae 5-day-LC 50 value of 22.4 ng/ml [5], and yet ivermectin MDAs in Senegal at doses of 150 µg/ kg substantially reduced An. gambiae 5-day survivorship by 43.6% for up to 6 days post MDA [10]. Furthermore, the pharmacokinetic model predicts that only 0.432 ng/ml of ivermectin parent compound would be present in a typical patient of 55 kg body weight 6 days post MDA at the 150 µg/kg dose, which is well below the concentration capable of killing An. gambiae. This clearly illustrates that in vitro mosquito membrane feeds and pharmacokinetic predictions of parent compound likely underestimate the full mosquito-lethal potential of ivermectin-treated humans. One possible explanation may be that ivermectin produces active mosquito-lethal in vivo metabolites with different pharmacokinetic properties that extend the duration of mosquito-lethal effect beyond that of the parent compound. Human liver microsomes have been used to characterize some of the ivermectin metabolites [32]. A small (n = 4) mass balance study in humans determined that mean peak plasma concentration of metabolites was 2.5-fold greater than that of the parent compound and the effective half-lives of the metabolites was approximately 2.9 days while the parent compound half-life was 11.8 h [33]. Further attention to the characterization of ivermectin mosquitolethal metabolites is warranted, especially in light of the novel long-lasting ivermectin formulations in development for human use [34][35][36]. Plasmodium vivax oocyst prevalence (a) and intensity (b) in Anopheles minimus when ivermectin co-ingested with parasites. Plasmodium vivax oocyst prevalence (a) and intensity (b) in An. minimus when ivermectin LC 25 (11.3 ng/ml) and LC 5 (6.7 ng/ml) co-ingested with parasites at DPI 0. Oocyst prevalence was significantly reduced at the LC 25 and LC 5 concentrations as determined by the Fishers Exact test. Oocyst intensity was significantly reduced at the LC 5 but not the LC 25 concentration as determined by the by the Mann-Whitney U test. Prevalence error bars represent standard error

Table 1 Parameter estimates from final population pharmacokinetic model of ivermectin in healthy volunteers
The final population pharmacokinetic model developed in this study described ivermectin concentrations in healthy volunteers satisfactorily. Overall parameter estimates were similar to that reported previously in healthy volunteers and patients with onchocerciasis [37,38]. Simulation-based diagnostics demonstrated a good predictive performance, which suggests that the developed model is suitable to use for simulations. However, the model was developed on data from a single oral dose of 200 µg/kg and extrapolations beyond that (i.e. simulations of 400 and 800 µg/kg doses) assumes that ivermectin shows dose linearity at this dosing range [19]. Body weight, implemented as allometric function, produces a biologically plausible covariate relationship [39,40] between ivermectin exposure and body weight so that the developed pharmacokinetic model can be used to simulate other populations at risk, such as children. However, the pharmacokinetic properties of a drug can be very different in children and adults due to the rapid change in body size, organ function, body composition, and enzyme maturation, which occur during the early years of life. Prospective clinical trials are urgently needed in children since there are currently no pharmacokinetic assessments of ivermectin in children and adolescents below 17 years of age. A previous study [8] reported that female participants and participants with higher body mass index (BMI) had higher day-7 ivermectin concentrations. However, BMI was not found to be a significant covariate on any pharmacokinetic parameters in the current analysis, perhaps due to the relatively narrow range of BMI (17.8-22.8 kg/m 2 ) studied here. Simulation of ivermectin concentration-time profiles at different dosing regimens were also performed in a previous study [41]. The authors obtained pharmacokinetic parameters from the literature, based on American healthy volunteers and assuming a 30% inter-individual variability in each parameter. Their simulated concentration-time profile after a single dose of 800 µg/kg of ivermectin resulted in an average peak concentration of 108.1 ng/ml [41], which is somewhat lower than the value simulated using the model developed here (174 ng/ml).
The sporogony experiments conducted here clearly indicate a sporontocidal effect of ivermectin against P. vivax in both An. dirus and An. minimus (Figs. 5, 6). Both P. vivax oocyst prevalence and intensity were significantly reduced when ivermectin compound at the LC 25 and LC 5 when ingested concomitantly (DPI 0) with parasites in An. dirus and An. minimus. A previous study found that ivermectin compound significantly reduced P. falciparum prevalence in An. gambiae when ingested concomitantly (DPI 0) at the ivermectin LC 25 but not the LC 5 , and had no effect on oocyst intensity at either concentration [13]. A sporontocidal effect would reduce onward transmission from infected persons that received ivermectin MDA while gametocytaemic. This may be more relevant for transmission suppression for P. vivax than P. falciparum as P. vivax gametocytes mature much more rapidly and are therefore present before people become ill enough to seek treatment [67]. These studies only assess the sporontocidal effect of ivermectin in the mosquito, future studies should investigate the potential gametocytocidal action of ivermectin.
There was no effect of ivermectin LC 5 compound on P. vivax oocyst prevalence (Fig. 7) or intensity (Additional file 1: Figure S2) when ingested by An. dirus at DPI −3. This is in direct contrast to previous findings which showed ivermectin compound LC 5 at DPI −3 reduced oocyst prevalence of P. falciparum in An. gambiae [13]. Ingestion of blood from ivermectin-treated cattle 4 days before ingestion of field isolates of P. falciparum from infected patients did not have a sporontocidal effect [68]. Ivermectin LC 25 compound failed to reduce P. vivax sporozoite prevalence in An. dirus when ingested at DPI 6 or 9 (Fig. 7), while there was a significant sporontocidal effect at the ivermectin LC 25 at DPI 6 and 9 for P. falciparum in An. gambiae [13]. This indicates that ivermectin can have differential sporontocidal impact with different Anopheles and Plasmodium species combinations.
Previously it was suggested that ivermectin sporontocidal effect may be due to direct effects on the mosquito, specifically ivermectin delaying or altering formation of the peritrophic matrix [13]. Indeed, ivermectin has been shown to delay and alter peritrophic matrix formation in Aedes aegypti [69], delay blood meal digestion in An. gambiae [5], and upregulates peritrophic matrix gene expression following ivermectin ingestion in An. gambiae [70]. However, limited investigations by two other laboratories were unable to identify a sporontocidal effect of ivermectin against P. falciparum NF54 strain in An. gambiae or Anopheles stephensi [8] or with field isolates of P. falciparum in An. gambiae [71]. This suggests more complex interaction of ivermectin on mosquito and parasite interaction, possibly by influencing mosquito midgut microbiota. The mosquito midgut microbiota composition has recently been shown to dramatically alter Anopheles immune response and thus Plasmodium infection [72]. Original investigations of avermectin, the biological precursor of ivermectin, suggested no direct effect on a range of bacteria species [73]. However, more recent evidence indicates that avermectin can inhibit growth of Staphylococcus aureus [74], and ivermectin can inhibit growth of Chlamydia trachomatis [75] and Mycobacterium tuberculosis [76]. This recent antibacterial evidence suggests that ivermectin may influence mosquito midgut microbiota composition, which could in turn alter Plasmodium infection. Recently, the midgut microbiota composition was shown to significantly alter formation of the peritrophic matrix [77], thus possible alteration of the midgut microbiota by ivermectin may influence formation of the peritrophic matrix and in turn alter Plasmodium infection. Since different insectaries around the world likely have different midgut bacteria microbiota compositions in colonized mosquitoes, this may explain the differences observed by different laboratories when investigating ivermectin sporontocidal effects. Much remains to be explored to determine the sporontocidal mode of action of ivermectin.
Concomitant ingestion of P. vivax and ivermectin LC 25 failed to alter An. dirus survivorship compared to ivermectin LC 25 alone (Fig. 8). This in contrast to the finding that P. falciparum-infectious An. gambiae that ingested ivermectin LC 25 at DPI 14 were significantly, albeit modestly, more susceptible to ivermectin compared to uninfected An. gambiae [13]. There were some caveats to the P. vivax and An. dirus survival study in that mosquitoes had to be transported between and housed in two different insectaries and P. vivax blood was collected freshly from a single donor in sodium heparin tubes and uninfected blood from multiple donors was stored briefly in CPDA-1 bags. There was no significant difference in An. dirus survivorship when fed a control or ivermectin LC 25 meal mixed with blood collected in CPDA-1 or sodium heparin tubes (Additional file 1: Figure S3, additional information text). Future investigations are warranted to determine whether Plasmodium infection in Anopheles alters susceptibility to ivermectin.
During ivermectin MDAs for malaria, numerous neglected tropical diseases (NTDs) in the GMS would be affected, including lymphatic filariasis, scabies, lice, gnathostomiasis, and soil-transmitted helminths (STHs) such as strongyloidiasis, ascariasis and trichuriasis. It has been estimated that 50% of persons in resource-poor communities in the GMS have one or more STH [78], and strongyloidiasis afflicts between 40 to 60% of persons in rural Cambodia [79] and Laos [80]. Ivermectin was found to be very effective at treating strongyloidiasis and repeated MDAs would further benefit afflicted communities as re-infection rates in Cambodia can be quite high [81]. Years of experience with ivermectin MDA for onchocerciasis in Africa have demonstrated that treated persons clearly recognize and appreciate the secondary benefits that ivermectin treatment has on NTDs [82].
MDA with artemisinin combination therapy (ACT) is being piloted in the GMS [1] and Africa [83]. One problem facing ACT MDA is that asymptomatic persons may not perceive a direct personal benefit of clearing malaria parasites. However, if ivermectin were incorporated into ACT MDA then this may improve compliance as persons could observe a direct personal benefit to MDA participation by reducing NTD burdens. Ivermectin and ACT MDA for malaria control, when rationally deployed in targeted hotspot areas of active Plasmodium transmission, would act in concert by clearing infected persons of their malaria parasites while reducing the transmission potential of the extant mosquito population. This could essentially reset malaria transmission when the next wave of naïve Anopheles emerge from the larval habitat and feed on a population cleared of their malaria parasites. Modelling suggests that if ivermectin is added to anti-malarial drug MDA this will reduce the number of MDAs and time required to achieve elimination [84]. The combination of ivermectin and artemether-lumefantrine was shown to be very safe [8]. Two clinical trials are currently being conducted to assess the safety, tolerability, pharmacokinetics and mosquito-lethal efficacy of ivermectin and dihydroartemisinin-piperaquine in Thailand (NCT02568098) and Kenya (NCT02511353). If ivermectin can integrate with ACT MDA, then this could become a powerful new tool for malaria elimination.

Conclusions
Key malaria vectors in the GMS are susceptible to ivermectin at human relevant concentrations. Ivermectin is sporontocidal to P. vivax in An. dirus and An. minimus. A population pharmacokinetic model was developed for ivermectin after a single oral dose administration of 200 µg/ kg in healthy volunteers. The developed model described observed data adequately and was used to simulate population mean concentration-time profiles after 200, 400 and 800 µg/kg doses of ivermectin. Mosquito in vitro survivorship results and pharmacokinetic model results indicate that the ivermectin dose of 400 µg/kg is the ideal minimal dose to administer during MDA in the GMS. There may be active ivermectin metabolites that extend the duration of mosquito-lethal effect beyond what is predicted by parent compound concentrations. Ivermectin MDA could be a powerful new tool to combat outdoor malaria transmission and assist malaria elimination efforts in the GMS.