Skip to main content

Is outdoor vector control needed for malaria elimination? An individual-based modelling study

Abstract

Background

Residual malaria transmission has been reported in many areas even with adequate indoor vector control coverage, such as long-lasting insecticidal nets (LLINs). The increased insecticide resistance in Anopheles mosquitoes has resulted in reduced efficacy of the widely used indoor tools and has been linked with an increase in outdoor malaria transmission. There are considerations of incorporating outdoor interventions into integrated vector management (IVM) to achieve malaria elimination; however, more information on the combination of tools for effective control is needed to determine their utilization.

Methods

A spatial individual-based model was modified to simulate the environment and malaria transmission activities in a hypothetical, isolated African village setting. LLINs and outdoor attractive toxic sugar bait (ATSB) stations were used as examples of indoor and outdoor interventions, respectively. Different interventions and lengths of efficacy periods were tested. Simulations continued for 420 days, and each simulation scenario was repeated 50 times. Mosquito populations, entomologic inoculation rates (EIRs), probabilities of local mosquito extinction, and proportion of time when the annual EIR was reduced below one were compared between different intervention types and efficacy periods.

Results

In the village setting with clustered houses, the combinational intervention of 50% LLINs plus outdoor ATSBs significantly reduced mosquito population and EIR in short term, increased the probability of local mosquito extinction, and increased the time when annual EIR is less than one per person compared to 50% LLINs alone; outdoor ATSBs alone significantly reduced mosquito population in short term, increased the probability of mosquito extinction, and increased the time when annual EIR is less than one compared to 50% LLINs alone, but there was no significant difference in EIR in short term between 50% LLINs and outdoor ATSBs. In the village setting with dispersed houses, the combinational intervention of 50% LLINs plus outdoor ATSBs significantly reduced mosquito population in short term, increased the probability of mosquito extinction, and increased the time when annual EIR is less than one per person compared to 50% LLINs alone; outdoor ATSBs alone significantly reduced mosquito population in short term, but there were no significant difference in the probability of mosquito extinction and the time when annual EIR is less than one between 50% LLIN and outdoor ATSBs; and there was no significant difference in EIR between all three interventions. A minimum of 2 months of efficacy period is needed to bring out the best possible effect of the vector control tools, and to achieve long-term mosquito reduction, a minimum of 3 months of efficacy period is needed.

Conclusions

The results highlight the value of incorporating outdoor vector control into IVM as a supplement to traditional indoor practices for malaria elimination in Africa, especially in village settings of clustered houses where LLINs alone is far from sufficient.

Background

Huge progress has been made toward malaria elimination following the scale up of long-lasting insecticidal nets (LLINs) and indoor residual spraying (IRS) [1,2,3,4,5], the two major vector control tools recommended by the World Health Organization (WHO). However, outdoor residual malaria transmission has been consistently reported in many areas where these interventions are in place [6, 7]. The wide use of LLINs and IRS has resulted in increased insecticide resistance in the malaria vector mosquitoes both chemically [8,9,10] and behaviourally [7, 11,12,13,14,15,16,17], reducing the efficacy of these tools and increasing outdoor malaria transmission [18]. Studies have suggested that indoor interventions alone are not sufficient to achieve malaria elimination, especially in places with high malaria transmission [7, 11, 19, 20]. Hence, it is of interest to evaluate the necessity of using outdoor vector control tools to complement the indoor control by LLINs and IRS to achieve malaria elimination.

A major benefit of outdoor vector control is that it can target and kill mosquitoes in their natural outdoor habitats and is not limited only to the vicinity of residential houses. In addition, it can target the vectors displaying more outdoor human biting activity in response to the indoor interventions [7, 11, 12, 18, 21]. There are several outdoor adult mosquito control methods, such as thermal fogging [22, 23], insecticide-treated cattle [24,25,26], male swarm spraying [27], “push–pull system” (indoor spatial repellent and outdoor traps) [28, 29], and mosquito landing box [30, 31]. There are also outdoor larval source management tools which are recommended by the WHO as supplementary malaria vector control methods [32, 33], but the feasibility is limited because of the large number of small and temporary larval habitats of malaria vectors [34]. One of the most promising new tools for outdoor vector control is attractive toxic sugar baits (ATSBs) [35, 36], which target the sugar-feeding behaviour of mosquitoes [37] and can be used outdoors either sprayed on vegetation or as bait stations. Several field trials have demonstrated the efficacy of both outdoor and indoor ATSBs against malaria vector mosquitoes [35, 36, 38,39,40]. An individual-based modelling (IBM) study has estimated and compared the effectiveness of different spatial configurations of ATSB treatment on malaria control [41]. In addition, insecticide resistance will not affect the efficacy of this tool because there are already plenty of oral toxins available from different chemical groups and having different modes of actions, which can be used/exchanged in ATSBs [38].

Both LLINs and ATSBs have limited efficacy periods. Studies in multiple locations have reported that LLINs were not fully used because of discarding [42, 43], giving away [42], or using for purposes other than indoor vector control [42]. In addition, the efficacy period of LLIN products varies, and they gradually lose efficacy due to reduced insecticidal activity, holes, damage, and improper use [42,43,44,45,46]. Similar for ATSBs, baits sprayed on vegetation may be washed out by rain, and although bait stations can resist rains, they may lose efficacy due to damage or expiration. Hence, predicting the recovery pattern of vector populations and malaria transmission after the efficacy period of these interventions can provide information on the durability of the interventions to guide the tool replacement schedule and further cost-effectiveness evaluations.

Integrated vector management (IVM) is important to extend the malaria elimination process. IVM suggests evidence-based decision-making and integrated approaches [47]. Considering the option of outdoor vector control, a basic question pertinent to developing IVM programmes is the optimal combination of indoor and outdoor interventions that are required to achieve desired reductions in malaria transmission given financial, technical, and logistical constraints. An IBM study has compared different indoor intervention strategies and suggested outdoor intervention should be added in settings with exophilic vectors; however, outdoor intervention was not evaluated directly in this study [48]. A mathematical modelling study has found that the use of LLINs in combination with ATSBs can achieve significantly better control of anopheline mosquitoes than LLINs combined with IRS, and the ATSB–LLIN combination was relatively effective against outdoor-biting Anopheles arabiensis mosquitoes [49]. However, the ATSB intervention was not specified as indoor or outdoor due to lack of a spatial component in the model.

Community trials are appropriate for these evaluations. However, it can be difficult to identify villages as comparable sites in both demographics and environmental factors. In addition, a randomized controlled trial (RCT) is very time and resource-consuming. To avoid ethical consequences, more complicated designs such as stepwise design may be needed. Therefore, a spatial IBM that simulates the life cycle of Anopheles gambiae mosquitoes and their interactions with the environment, humans, and the vector control interventions may be the most pragmatic way to initially predict the consequences of combining outdoor and indoor interventions to guide further empirical studies. The objectives of this study are to: (1) compare the immediate (during efficacy period) and long-term (after losing efficacy) reductions of malaria vectors and entomological inoculation rate (EIR) between indoor and outdoor interventions to evaluate the necessity of incorporating outdoor vector control into IVM programme; and (2) estimate the effect of efficacy periods to guide schedule of maintenance and replacement of the intervention tools.

Methods

Model design

For this study, a previous IBM that simulated the interactions between An. gambiae mosquitoes and humans, environment, and ATSBs was updated with the features of LLINs. The details of the basic model design were introduced in two previous studies, one evaluating the impact of environmental resources on the survival and biting behaviour of An. gambiae [50] and the other comparing the effectiveness of different spatial configurations of ATSB treatment [41]. A summary of the model used in this study is as follows:

Environment

A hypothetical isolated African village of 600 × 600 m with 25 houses was simulated. To simulate an environment representing the village configuration of clustered houses like those is Mali, the 25 houses were randomly located in the 100 × 100 m center of the simulated environment; to simulate the other type of environment which has dispersed settlement of houses, the 25 houses were randomly located in the whole simulated environment. Fifty natural sugar sources, 50 outdoor resting sites, and 50 larval habitats were randomly placed in the whole village to simulate a “resource rich” environment.

Human

A hundred humans were simulated, and they walked randomly outdoors from 7:00 to 19:59, went back home and walked inside their houses from 20:00 to 22:59, after 23:00 they went to bed and were protected by LLINs if covered.

Female Anopheles gambiae

Female An. gambiae were active during night (19:00 to 05:00). Each mosquito had its physiological status recorded in the model: age, energy level, need for blood, and gravid status. In each step, each mosquito performed different behaviours: fly, sugar feed, blood feed, oviposit, or rest, depending on the environmental factors and their own physiological status. The maximum lifespan was 30 days [51,52,53]. The extrinsic incubation period was 10 days [54, 55]. The mosquitoes were assumed to have a 20% chance of being infected by biting a human once [56,57,58,59,60,61,62,63], natural history of malaria in human was not considered here because of the complexity and uncertainty of malaria immunity [64]. More details of the density-dependent development of mosquito aquatic stages were described in a previous study [41].

Intervention

Long-lasting insecticidal nets were used as the example of the indoor intervention since LLINs are the WHO-recommended and commonly used indoor tool. When a LLIN was efficacious, it protected a human sleeping inside it from mosquito bites and killed mosquitoes that contacted with it; no insecticide resistance was simulated in the model because limited data is available on it, additional assumptions could lead to greater uncertainty. Beyond the efficacy period, both physical and chemical protection was removed to mimic holes on the nets and the loss of insecticide efficacy. ATSB is considered one of the best outdoor tools available currently, and so bait stations were placed in 7 × 7 grid configuration [41] over the whole site as the example of the outdoor intervention. The bait stations attracted and killed sugar-feeding mosquitoes. The efficacy for both LLINs and ATSBs were simulated to be binomial (on/off) instead of gradual degrading because data is limited on the degrading pattern of the two tools, and physical damage could be common in both tools, which is consistent with the binomial pattern. In addition, the efficacy periods were increased by month, increasing the resolution to day or second may not have significant impact on the comparisons between different interventions. Three interventions were simulated: 50% LLIN coverage, ATSBs alone, and 50% LLIN coverage plus ATSBs. The three interventions were introduced in different simulations at the beginning of day 1 (the initial equilibration of mosquito population), and efficacy periods (time when interventions remain efficacious) of 30, 60, 90, 120, 150, and 180 days were simulated. Two controls for each of the three interventions were used: one was the negative control, in which no intervention was applied; the other was the positive control, in which continuous efficacy period until the end of the simulation (day 420) was simulated. Each simulation scenario was repeated 50 times.

Data analysis

Four measures were calculated and analysed to examine the effects of different interventions and efficacy periods. To describe the best possible effects while the interventions were efficacious, means of mosquito population sizes and annual EIRs at the last day of the efficacy periods (e.g. day 60 for ATSB intervention with an efficacy period of 60 days, day 90 for 50% LLIN coverage with an efficacy period of 90 days) were compared between different interventions and efficacy periods. To control for the over dispersion and fit the data distribution, negative binomial regression models with a square root transformation of the outcomes were selected. Intervention type, efficacy period, and their interaction were included in the model. The annual EIRs were derived as daily EIR times 365 days, daily EIRs were calculated as the total number of infectious bites per human per day. A value of 0.0001 was added to EIRs to eliminate cells with zero variance and enable model fitting. To describe the long-term effects after the interventions lost efficacy, the probabilities of mosquito population extinction at the end of simulation (day 420) were compared between different interventions and efficacy periods using a logistic regression model. Intervention type, efficacy period, and their interaction were included in the model. Interventions that had none of the 50 trials ending up with extinction were excluded from the analysis. To describe the overall effects, the total number of days while EIR was reduced to the level below one per person per year for each intervention and efficacy period were compared using a negative binomial regression model.

Results

The environment simulated in this model represented a malaria endemic village: in the 600 × 600 m area, female An. gambiae mosquito populations equilibrated at around 400 without any intervention. Annual EIR equilibrated at around 200 infectious bites per person per year.

Figure 1 shows the dynamics of female mosquito population with different intervention types and efficacy periods in the two types of village configurations. Mosquito populations dropped sharply in the first 30 days following the treatments and then decreased slowly toward equilibration in the next 30 days. Mosquito populations began to reestablish immediately after the efficacy period of the interventions. However, with different types of interventions and efficacy periods, the re-equilibrated mosquito populations averaged at different levels. The mosquito populations were the means (averages) of 50 repeated simulations. In a proportion of the 50 repetitions, mosquito populations were reduced to very low levels and reestablished to the negative control level. In the other repetitions, mosquito populations were annihilated, so smaller means of re-equilibrated population sizes represent lower chances of recovery to a negative control level. Interventions of or with outdoor ATSB treatment led to better mosquito control results, especially in village setting of clustered houses. The detailed comparison of immediate and long-term effects of different interventions and efficacy periods in the two village configurations are presented in Tables 1, 2 and 3.

Fig. 1
figure 1

Dynamics of female mosquito population with different intervention types and efficacy periods in clustered and dispersed village settings. The two columns of figures show the dynamics of female mosquito populations in village setting of clustered houses (on the left) and dispersed houses (on the right). The first row of figures shows results of 50% LLIN intervention, the second row shows results of outdoor ATSB intervention, and the third row shows results of 50% LLIN plus ATSB intervention. In all six figures, the y-axis represents the mean number of female An. gambiae mosquitoes at each day, and the x-axis shows the time in days. Day 0 is the day that the intervention was applied. The lines in different colors represent different simulated efficacy periods for each intervention

Table 1 Average population sizes and EIRs at the last day of efficacy periods for different intervention types and efficacy periods
Table 2 Proportions (%) of trials ending up with local mosquito extinction for different intervention types and efficacy periods
Table 3 Proportions (%) of time when EIR was below 1 per person per year for different intervention types and efficacy periods

Table 1 shows the means of mosquito population sizes and annual EIRs at the last day of the efficacy periods for different intervention types and efficacy periods. In the village setting with clustered houses, female population size was reduced from over 400 in negative control to a range of 2–193, and annual EIR was reduced from over 200 to a range of 0–4 with different intervention types and efficacy periods. The results of the negative binomial regression model show that all three interventions reduced mosquito population and EIR to significantly lower levels than negative control (Ps < 0.0001), outdoor ATSB intervention reduced mosquito population and EIR to significantly lower levels than 50% LLIN intervention (P < 0.0001), the intervention of 50% LLINs plus outdoor ATSBs reduced mosquito population to significantly lower level than the intervention of ATSBs alone (P < 0.0001), however, EIR was not significantly different between the combinational intervention and outdoor ATSB alone (P = 0.3178). Interventions with efficacy period of 2 months or longer reduced mosquito population to significantly lower levels than efficacy period of 1 month (Ps < 0.0001), and there was no significant difference in EIR between different efficacy periods (Ps > 0.05).

In the village setting with dispersed houses, female population size was reduced from over 400 in negative control to a range of 1–36, and annual EIR was reduced from over 200 to a range of 0–0.4 with different intervention types and efficacy periods. All the three interventions reduced population size and EIR to significantly lower levels than negative control (P < 0.0001). Although the differences of population size between each intervention type were still statistically significant (Ps < 0.0001), the absolute differences between 50% LLIN and outdoor ATSB were much smaller than those in the village setting with clustered houses. There was no significant difference in EIR between the three interventions (P > 0.05). In addition, there was no significant difference in EIR between different efficacy periods (P > 0.05).

Table 2 shows the probabilities (proportions of trials) of mosquito extinction for each intervention type and efficacy period after the interventions lost efficacy. In the village setting of clustered houses, the probability of mosquito extinction was increased from 0 to a range of 0–90% by the three interventions. Outdoor ATSB intervention increased the probability to significantly higher levels than 50% LLIN intervention (P < 0.0001), and the intervention of 50% LLINs plus outdoor ATSBs increased the probability to significantly higher level than the intervention of outdoor ATSBs alone (P < 0.0001). Longer efficacy periods significantly increased the probability of mosquito extinction (Ps < 0.0001), except that there was no significant difference between efficacy periods of 3–6 months (Ps > 0.05).

In the village setting of dispersed houses, the probability of mosquito extinction was increased from 0 to a range of 4–90% by the three interventions. The intervention of 50% LLINs plus outdoor ATSBs increased the probability to significantly higher level than the other two interventions (Ps < 0.0001), but there was no significant difference between the intervention of 50% LLINs and outdoor ATSB alone (P = 0.2286), except that with efficacy period of 1 month, 50% LLINs intervention led to a significant higher probability than outdoor ATSB intervention (P = 0.0108). The comparison between different efficacy periods was similar with that in the village setting of clustered houses: longer efficacy periods significantly increased the probability of mosquito extinction (Ps < 0.0001), except that there was no significant difference between efficacy periods of 3–6 months (Ps > 0.05).

Table 3 shows the proportion of time that EIR was reduced to below one per person per year for different intervention types and efficacy periods. In the village setting of clustered houses, the proportion of time that EIR was less than one was increased from 13.68% to a range of 20.31–98.77% by the three interventions. Outdoor ATSB intervention increased the proportion of time to significantly higher levels than 50% LLIN intervention (P < 0.0001), and the intervention of 50% LLINs plus outdoor ATSBs increased the probability to significantly higher level than the intervention of outdoor ATSBs alone (P = 0.0002). Longer efficacy periods significantly increased the proportion of time most of times (Ps < 0.05), except that there was no significant difference between efficacy periods of 3 and 4 months (P = 0.2666), 4 and 5 months (P = 0.3015), 4 and 6 months (P = 0.0940), and 5 and 6 months (P = 0.5207).

In the village setting of dispersed houses, the proportion of time that EIR was less than one was increased from 4.35% to a range of 20.68–99.00% by the three interventions. The intervention of 50% LLINs plus outdoor ATSBs increased the probability to significantly higher level than the other two interventions (Ps < 0.0001), but there was no significant difference between the intervention of 50% LLINs and outdoor ATSB alone (P = 0.6493), except that with efficacy period of 1 month, 50% LLINs intervention led to a significant higher proportion than outdoor ATSB intervention (P < 0.0001). Longer efficacy periods significantly increased the proportion of time most of times (Ps < 0.05), except that there was no significant difference between efficacy periods of 3–6 months (Ps > 0.05).

Discussion

In this study, the potential benefits of incorporating outdoor vector control into the widely used indoor vector control strategy (LLINs) in a malaria endemic village were evaluated by simulating three types of indoor and outdoor control interventions with different efficacy periods. According to WHO recommendations and current research findings, LLINs with 50% coverage were used as the example of indoor mosquito control; a pattern of ATSB stations placed in a 7 × 7 grid configuration over the whole area was used as the example of outdoor mosquito control. Incorporating outdoor ATSB intervention significantly improved both immediate and long-term effects of mosquito and malaria control, especially in the village setting with clustered houses.

The mosquito population size and EIR at the last day of efficacy periods is the best possible effect of an intervention with a certain efficacy period, which represents the immediate effect. The results show that all three types of interventions are effective at reducing mosquito population and EIR, however, to achieve a better mosquito control result, outdoor ATSBs alone or the combination of 50% LLINs and outdoor ATSBs is suggested, especially in village setting of clustered houses where even outdoor ATSBs alone achieved substantially lower population size than 50% LLINs. The difference in EIR between different interventions was smaller, thus to only achieve EIR reduction, if resources are limited, outdoor ATSBs alone in village setting of clustered houses is acceptable, and either 50% LLINs or outdoor ATSBs alone in village setting of dispersed houses is acceptable. In addition, for any interventions, a minimum of 2 months of efficacy period is needed to bring out the best possible effectiveness. This provides guidance for schedule of tool maintenance and education for the human compliance.

The probability of local mosquito extinction after the interventions lose efficacy represents the long-term effect. The explanation for the long-term mosquito reduction may be that a good combination of interventions can rapidly knock down the mosquito populations to very low levels; and after a certain time, as the existing eggs develop to adults and are killed and as overall reproduction decreases, the population can be locally eliminated. In both village settings, the combinational intervention of 50% LLIN and outdoor ATSBs was suggested to achieve the best long-term mosquito control results. Even with very limited resources, intervention of 50% LLINs alone is not suggested in village setting of clustered houses. In addition, a minimum of 3 months of efficacy period is needed to achieve the best long-term mosquito reduction.

There is a study showing that malaria elimination can be achieved when annual EIR is below one [65]. The results of the proportion of time when annual EIR is less than one represents the overall protection. In both village settings, incorporating outdoor ATSBs into the LLIN intervention with an efficacy period of 3 months or longer achieved the same level of protection as the other two interventions with continuous efficacy until the end of simulation. This suggests that an investment of the combinational intervention short-term may save the work of long-time maintenance of the interventions.

The underlying mechanism for the increase in effectiveness by the incorporation of the outdoor ATSB intervention is that the outdoor ATSBs not only killed mosquitoes that entered houses, but also killed those searching for sugar meals outdoors. LLINs only target the mosquitoes searching for a blood meal indoors, and these mosquitoes constitute a proportion of the mosquito population since females would blood feed mainly once per gonotrophic cycle (3 days in average) [66,67,68]. Whereas outdoor ATSBs target the mosquitoes searching for sugar meals, which constitute a larger proportion of the mosquito population since they require at least one sugar meal per night for energy and more activity would increase the needs for sugar meals [37, 69]. In addition, ATSBs can well target male An. gambiae mosquitoes because males solely feed on sugar for energy and they require at least two sugar meals per night [69]. Male mosquitoes were simulated in ten repetitions and followed similar population dynamic patterns as females, so they were not included in the results. However, as outdoor ATSBs killed more male mosquitoes, it would reduce the mating opportunities for females, which would either lead to lower mosquito reproduction or additional flights for females and hence increased mortality. Successful mating was assumed in the model, if mating opportunity was taken into account, outdoor ATSB intervention may have even better vector control results. In the current model, a sugar-rich village was simulated. Whereas in villages that lack natural sugar resources, or during dry seasons, the outdoor ATSB intervention could achieve even better results [41].

A LLIN coverage of 50% was used in this study, this is due to the fact that according to the estimation by WHO, the average bed net coverage increased from less than 2% in 2000 to an estimated 55% in 2015 [70]. Except for the considerations of incorporating outdoor vector control, there are also considerations of scaling up LLIN coverage for malaria elimination. However, most studies evaluating the effectiveness of scaling up LLIN coverage investigated the effect where the coverage rate increased from very low to intermediate/high level and showed significant reduction of malaria transmission [71,72,73,74], but there is no empirical evidence showing that further increasing the coverage can reduce malaria transmission: only one study evaluated locations with this magnitude of increase, and the results showed that asexual parasite prevalence decreased significantly, but gametocyte prevalence did not [75]. In addition, results from another RCT showed that universal LLIN coverage did not reduce malaria transmission than targeted LLIN coverage for pregnant women and children under six [76]. These are consistent with the results of 100% LLIN intervention in this study, which are not reported here. The results of the simulations show that there was no significant difference between 100% coverage and 50% coverage. Another explanation for this phenomenon may be that no insecticide resistance was simulated for LLIN. However, with 50% LLIN coverage, if behavioural resistance in mosquitoes was simulated, the mosquitoes could avoid LLINs and select those unprotected human for blood-feeding. Since those unprotected humans could then compensate the reduced bites for those protected under LLINs, the EIR could be higher in 50% LLINs in reality, and thus the difference between 50 and 100% coverages could be greater. Because of the simplification of not including insecticide resistance in this model, the relevant results can be inaccurate and are not included. Further modelling studies with these details and empirical studies are needed for the confirmation of the effectiveness of scaling up LLIN coverage.

One limitation of the study was that the simulated village was isolated/closed to influx of mosquitoes from the outside. Thus, after the mosquito population was knocked down, the effect could be sustained even after the loss of interventions efficacy. However, previous studies have shown that malaria can be locally eliminated when vector density is below a certain threshold [77, 78]. If the influx of mosquitoes is below a certain level and not able to reestablish the population, the effect of control may still be sustained. In addition, anopheline mosquitoes have a maximum dispersal range from several 100 m to almost 3000 m [79]. With limited human migration between villages, if the villages are separated by distances exceeding the maximum dispersal range, especially in arid environments, or if several villages are treated at the same time in a big area exceeding the maximum dispersal range, then “isolation” can be achieved and local extinction may be sustained.

In this study, ATSB stations placed outdoor in 7 × 7 grid configuration over the whole area was used to represent outdoor adult mosquito control in general. ATSB attracts sugar-seeking mosquitoes from a distance and are not limited by constraints of insecticide resistance; in fact, ATSB is being considered as a new tool for breaking insecticide resistance. Generalization from the study results from this study to other available outdoor vector control tools should be careful because of these advantages of ATSB tool. No insecticide resistance was simulated for the indoor LLIN treatments, and no damage was simulated for the outdoor ATSB treatments during its efficacy period. However, in reality, these perfect conditions are never achieved. Therefore, this study may have over-estimated the effectiveness of both tool. Nevertheless, even with the simplifications in this modelling study, it provides insights of how expanding mosquito control from inside houses to a larger scale outdoors can benefit overall vector and malaria control.

Conclusions

In summary, incorporating outdoor vector control interventions such as ATSBs into IVM strategies is suggested for malaria elimination, especially in village setting of clustered houses where indoor LLINs alone are far from sufficient. With the combination of indoor and outdoor vector control interventions, it is likely that mosquito population will become locally extinct even after the efficacy period of the interventions in certain conditions, leading to long-term mosquito control and EIR reduction. Further confirmation from empirical studies are recommended to provide additional evidence for updating the malaria elimination policy.

Abbreviations

LLIN:

long-lasting insecticidal net

ATSB:

attractive toxic sugar bait

IVM:

integrated vector management

EIR:

entomological inoculation rate

WHO:

World Health Organization

RCT:

randomized controlled trial

References

  1. Mabaso ML, Sharp B, Lengeler C. Historical review of malarial control in southern African with emphasis on the use of indoor residual house-spraying. Trop Med Int Health. 2004;9:846–56.

    Article  PubMed  Google Scholar 

  2. Sharp BL, Kleinschmidt I, Streat E, Maharaj R, Barnes KI, Durrheim DN, et al. Seven years of regional malaria control collaboration—Mozambique, South Africa, and Swaziland. Am J Trop Med Hyg. 2007;76:42–7.

    PubMed  PubMed Central  Google Scholar 

  3. Kleinschmidt I, Schwabe C, Benavente L, Torrez M, Ridl FC, Segura JL, et al. Marked increase in child survival after four years of intensive malaria control. Am J Trop Med Hyg. 2009;80:882–8.

    PubMed  PubMed Central  Google Scholar 

  4. Chizema-Kawesha E, Miller JM, Steketee RW, Mukonka VM, Mukuka C, Mohamed AD, et al. Scaling up malaria control in Zambia: progress and impact 2005–2008. Am J Trop Med Hyg. 2010;83:480–8.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Steketee RW, Campbell CC. Impact of national malaria control scale-up programmes in Africa: magnitude and attribution of effects. Malar J. 2010;9:299.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Killeen GF. Characterizing, controlling and eliminating residual malaria transmission. Malar J. 2014;13:330.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Russell TL, Govella NJ, Azizi S, Drakeley CJ, Kachur SP, Killeen GF. Increased proportions of outdoor feeding among residual malaria vector populations following increased use of insecticide-treated nets in rural Tanzania. Malar J. 2011;10:80.

    Article  PubMed  PubMed Central  Google Scholar 

  8. N’Guessan R, Corbel V, Akogbéto M, Rowland M. Reduced efficacy of insecticide-treated nets and indoor residual spraying for malaria control in pyrethroid resistance area, Benin. Emerg Infect Dis. 2007;13:199.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Ranson H, N’Guessan R, Lines J, Moiroux N, Nkuni Z, Corbel V. Pyrethroid resistance in African anopheline mosquitoes: what are the implications for malaria control? Trends Parasitol. 2011;27:91–8.

    Article  CAS  PubMed  Google Scholar 

  10. Padonou GG, Sezonlin M, Ossé R, Aizoun N, Oké-Agbo F, Oussou O, et al. Impact of three years of large scale Indoor Residual Spraying (IRS) and Insecticide Treated Nets (ITNs) interventions on insecticide resistance in Anopheles gambiae s.l. in Benin. Parasit Vectors. 2012;5:72.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Reddy MR, Overgaard HJ, Abaga S, Reddy VP, Caccone A, Kiszewski AE, et al. Outdoor host seeking behaviour of Anopheles gambiae mosquitoes following initiation of malaria vector control on Bioko Island, Equatorial Guinea. Malar J. 2011;10:184.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Moiroux N, Gomez MB, Pennetier C, Elanga E, Djènontin A, Chandre F, et al. Changes in Anopheles funestus biting behavior following universal coverage of long-lasting insecticidal nets in Benin. J Infect Dis. 2012;206:1622–9.

    Article  CAS  PubMed  Google Scholar 

  13. Durnez L, Mao S, Denis L, Roelants P, Sochantha T, Coosemans M. Outdoor malaria transmission in forested villages of Cambodia. Malar J. 2013;12:329.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Sougoufara S, Diédhiou SM, Doucouré S, Diagne N, Sembène PM, Harry M, et al. Biting by Anopheles funestus in broad daylight after use of long-lasting insecticidal nets: a new challenge to malaria elimination. Malar J. 2014;13:125.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Russell TL, Beebe NW, Bugoro H, Apairamo A, Chow WK, Cooper RD, et al. Frequent blood feeding enables insecticide-treated nets to reduce transmission by mosquitoes that bite predominately outdoors. Malar J. 2016;15:156.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Ndenga BA, Mulaya NL, Musaki SK, Shiroko JN, Dongus S, Fillinger U. Malaria vectors and their blood-meal sources in an area of high bed net ownership in the western Kenya highlands. Malar J. 2016;15:1.

    Article  Google Scholar 

  17. Killeen GF, Govella NJ, Lwetoijera DW, Okumu FO. Most outdoor malaria transmission by behaviourally-resistant Anopheles arabiensis is mediated by mosquitoes that have previously been inside houses. Malar J. 2016;15:76.

    Article  Google Scholar 

  18. Govella NJ, Ferguson H. Why use of interventions targeting outdoor biting mosquitoes will be necessary to achieve malaria elimination. Front Physiol. 2012;3:199.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Bugoro H, Cooper RD, Butafa C, Iro’ofa C, Mackenzie DO, Chen C-C, Russell TL. Bionomics of the malaria vector Anopheles farauti in Temotu Province, Solomon Islands: issues for malaria elimination. Malar J. 2011;10:133.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Monroe A, Asamoah O, Lam Y, Koenker H, Psychas P, Lynch M, et al. Outdoor-sleeping and other night-time activities in northern Ghana: implications for residual transmission and malaria prevention. Malar J. 2015;14:35.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Meyers JI, Pathikonda S, Popkin-Hall ZR, Medeiros MC, Fuseini G, Matias A, et al. Increasing outdoor host-seeking in Anopheles gambiae over 6 years of vector control on Bioko Island. Malar J. 2016;15:239.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Xue R-D, Smith ML, Qualls WA. Field evaluation of truck-mounted thermal fog of DUET™ and Aqualuer® using different solvents against Aedes albopictus. J Am Mosq Control Assoc. 2013;29:301–3.

    Article  PubMed  Google Scholar 

  23. Fulcher A, Farooq M, Richardson AG, Smith ML, Scott JM, Gaines MK, Xue R-D. Characteristics and efficacy of three commercial handheld thermal foggers with pyrethroid insecticides against three species of mosquitoes. J Am Mosq Control Assoc. 2016;32:44–50.

    Article  PubMed  Google Scholar 

  24. Rowland M, Durrani N, Kenward M, Mohammed N, Urahman H, Hewitt S. Control of malaria in Pakistan by applying deltamethrin insecticide to cattle: a community-randomised trial. Lancet. 2001;357:1837–41.

    Article  CAS  PubMed  Google Scholar 

  25. Chaccour CJ, Kobylinski KC, Bassat Q, Bousema T, Drakeley C, Alonso P, et al. Ivermectin to reduce malaria transmission: a research agenda for a promising new tool for elimination. Malar J. 2013;12:153.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Poché RM, Burruss D, Polyakova L, Poché DM, Garlapati RB. Treatment of livestock with systemic insecticides for control of Anopheles arabiensis in western Kenya. Malar J. 2015;14:351.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Diabate A, Tripet F. Targeting male mosquito mating behaviour for malaria control. Parasit Vectors. 2015;8:347.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Wagman JM, Grieco JP, Bautista K, Polanco J, Briceño I, King R, Achee NL. The field evaluation of a push–pull system to control malaria vectors in Northern Belize, Central America. Malar J. 2015;14:184.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Menger DJ, Omusula P, Wouters K, Oketch C, Carreira AS, Durka M, et al. Eave screening and push–pull tactics to reduce house entry by vectors of malaria. Am J Trop Med Hyg. 2016;94:868–78.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Matowo NS, Koekemoer LL, Moore SJ, Mmbando AS, Mapua SA, Coetzee M, et al. Combining synthetic human odours and low-cost electrocuting grids to attract and kill outdoor-biting mosquitoes: field and semi-field evaluation of an improved mosquito landing box. PLoS ONE. 2016;11:e0145653.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Mmbando AS, Okumu FO, Mgando JP, Sumaye RD, Matowo NS, Madumla E, et al. Effects of a new outdoor mosquito control device, the mosquito landing box, on densities and survival of the malaria vector, Anopheles arabiensis, inside controlled semi-field settings. Malar J. 2015;14:494.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Smith DL, Perkins TA, Tusting LS, Scott TW, Lindsay SW. Mosquito population regulation and larval source management in heterogeneous environments. PLoS ONE. 2013;8:e71247.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. WHO. Larval source management: a supplementary malaria vector control measure: an operational manual. Geneva: World Health Organization; 2013.

    Google Scholar 

  34. Majambere S, Pinder M, Fillinger U, Ameh D, Conway DJ, Green C, et al. Is mosquito larval source management appropriate for reducing malaria in areas of extensive flooding in The Gambia? A cross-over intervention trial. Am J Trop Med Hyg. 2010;82:176–84.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Müller GC, Kravchenko VD, Schlein Y. Decline of Anopheles sergentii and Aedes caspius populations following presentation of attractive toxic (spinosad) sugar bait stations in an oasis. J Am Mosq Control Assoc. 2008;24:147–9.

    Article  PubMed  Google Scholar 

  36. Müller GC, Schlein Y. Efficacy of toxic sugar baits against adult cistern-dwelling Anopheles claviger. Trans R Soc Trop Med Hyg. 2008;102:480–4.

    Article  PubMed  Google Scholar 

  37. Foster WA. Mosquito sugar feeding and reproductive energetics. Annu Rev Entomol. 1995;40:443–74.

    Article  CAS  PubMed  Google Scholar 

  38. Müller GC, Beier JC, Traore SF, Toure MB, Traore MM, Bah S, et al. Successful field trial of attractive toxic sugar bait (ATSB) plant-spraying methods against malaria vectors in the Anopheles gambiae complex in Mali, West Africa. Malar J. 2010;9:210.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Beier JC, Müller GC, Gu W, Arheart KL, Schlein Y. Attractive toxic sugar bait (ATSB) methods decimate populations of Anopheles malaria vectors in arid environments regardless of the local availability of favoured sugar-source blossoms. Malar J. 2012;11:31.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Qualls WA, Müller GC, Traore SF, Traore MM, Arheart KL, Doumbia S, et al. Indoor use of attractive toxic sugar bait (ATSB) to effectively control malaria vectors in Mali, West Africa. Malar J. 2015;14:301.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Zhu L, Marshall JM, Qualls WA, Schlein Y, McManus JW, Arheart KL, et al. Modelling optimum use of attractive toxic sugar bait stations for effective malaria vector control in Africa. Malar J. 2015;14:492.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Koenker H, Kilian A, de Beyl CZ, Onyefunafoa EO, Selby RA, Abeku T, et al. What happens to lost nets: a multi-country analysis of reasons for LLIN attrition using 14 household surveys in four countries. Malar J. 2014;13:464.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Batisso E, Habte T, Tesfaye G, Getachew D, Tekalegne A, Kilian A, et al. A stitch in time: a cross-sectional survey looking at long lasting insecticide-treated bed net ownership, utilization and attrition in SNNPR, Ethiopia. Malar J. 2012;11:183.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Hakizimana E, Cyubahiro B, Rukundo A, Kabayiza A, Mutabazi A, Beach R, et al. Monitoring long-lasting insecticidal net (LLIN) durability to validate net serviceable life assumptions, in Rwanda. Malar J. 2014;13:344.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Morgan J, Abílio AP, do Rosario Pondja M, Marrenjo D, Luciano J, Fernandes G, et al. Physical durability of two types of long-lasting insecticidal nets (LLINs) three years after a mass LLIN distribution campaign in Mozambique, 2008–2011. Am J Trop Med Hyg. 2015;92:286–93.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Tan KR, Coleman J, Smith B, Hamainza B, Katebe-Sakala C, Kean C, et al. A longitudinal study of the durability of long-lasting insecticidal nets in Zambia. Malar J. 2016;15:106.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Beier JC, Keating J, Githure JI, Macdonald MB, Impoinvil DE, Novak RJ. Integrated vector management for malaria control. Malar J. 2008;7(Suppl 1):S4.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Griffin JT, Hollingsworth TD, Okell LC, Churcher TS, White M, Hinsley W, et al. Reducing Plasmodium falciparum malaria transmission in Africa: a model-based evaluation of intervention strategies. PLoS Med. 2010;7:1000324.

    Article  Google Scholar 

  49. Marshall JM, White MT, Ghani AC, Schlein Y, Müller GC, Beier JC. Quantifying the mosquito’s sweet tooth: modelling the effectiveness of attractive toxic sugar baits (ATSB) for malaria vector control. Malar J. 2013;12:291.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Zhu L, Qualls WA, Marshall JM, Arheart KL, DeAngelis DL, McManus JW, et al. A spatial individual-based model predicting a great impact of copious sugar sources and resting sites on survival of Anopheles gambiae and malaria parasite transmission. Malar J. 2015;14:59.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Okech BA, Gouagna LC, Killeen GF, Knols BG, Kabiru EW, Beier JC, et al. Influence of sugar availability and indoor microclimate on survival of Anopheles gambiae (Diptera: Culicidae) under semifield conditions in western Kenya. J Med Entomol. 2003;40:657–63.

    Article  PubMed  Google Scholar 

  52. Afrane YA, Zhou G, Lawson BW, Githeko AK, Yan G. Effects of microclimatic changes caused by deforestation on the survivorship and reproductive fitness of Anopheles gambiae in western Kenya highlands. Am J Trop Med Hyg. 2006;74:772–8.

    PubMed  Google Scholar 

  53. Olayemi IKA, Ande AT. Life table analysis of Anopheles gambiae (Diptera: Culicidae) in relation to malaria transmission. J Vector Borne Dis. 2009;46:295–8.

    CAS  PubMed  Google Scholar 

  54. Lines J, Wilkes T, Lyimo E. Human malaria infectiousness measured by age-specific sporozoite rates in Anopheles gambiae in Tanzania. Parasitology. 1991;102:167–77.

    Article  PubMed  Google Scholar 

  55. Gary R, Foster W. Effects of available sugar on the reproductive fitness and vectorial capacity of the malaria vector. J Med Entomol. 2001;38:22–8.

    Article  PubMed  Google Scholar 

  56. Githeko A, Brandling-Bennett A, Beier M, Atieli F, Owaga M, Collins F. The reservoir of Plasmodium falciparum malaria in a holoendemic area of western Kenya. Trans R Soc Trop Med Hyg. 1992;86:355–8.

    Article  CAS  PubMed  Google Scholar 

  57. Shililu J, Maier W, Seitz H, Orago A. Seasonal density, sporozoite rates and entomological inoculation rates of Anopheles gambiae and Anopheles funestus in a high altitude sugarcane growing zone in western Kenya. Trop Med Int Health. 1998;3:706–10.

    Article  CAS  PubMed  Google Scholar 

  58. Sattabongkot J, Maneechai N, Phunkitchar V, Eikarat N, Khuntirat B, Sirichaisinthop J, et al. Comparison of artificial membrane feeding with direct skin feeding to estimate the infectiousness of Plasmodium vivax gametocyte carriers to mosquitoes. Am J Trop Med Hyg. 2003;69:529–35.

    PubMed  Google Scholar 

  59. Drakeley C, Sutherland C, Bousema JT, Sauerwein RW, Targett GA. The epidemiology of Plasmodium falciparum gametocytes: weapons of mass dispersion. Trends Parasitol. 2006;22:424–30.

    Article  PubMed  Google Scholar 

  60. Schneider P, Bousema JT, Gouagna LC, Otieno S, Van de Vegte-Bolmer M, Omar SA, et al. Submicroscopic Plasmodium falciparum gametocyte densities frequently result in mosquito infection. Am J Trop Med Hyg. 2007;76:470–4.

    PubMed  Google Scholar 

  61. Ouédraogo AL, Bousema T, Schneider P, De Vlas SJ, Ilboudo-Sanogo E, Cuzin-Ouattara N, et al. Substantial contribution of submicroscopical Plasmodium falciparum gametocyte carriage to the infectious reservoir in an area of seasonal transmission. PLoS ONE. 2009;4:e8410.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Bousema T, Drakeley C. Epidemiology and infectivity of Plasmodium falciparum and Plasmodium vivax gametocytes in relation to malaria control and elimination. Clin Microbiol Rev. 2011;24:377–410.

    Article  PubMed  PubMed Central  Google Scholar 

  63. Da DF, Churcher TS, Yerbanga RS, Yaméogo B, Sangaré I, Ouedraogo JB, et al. Experimental study of the relationship between Plasmodium gametocyte density and infection success in mosquitoes; implications for the evaluation of malaria transmission-reducing interventions. Exp Parasitol. 2015;149:74–83.

    Article  PubMed  Google Scholar 

  64. Langhorne J, Ndungu FM, Sponaas A-M, Marsh K. Immunity to malaria: more questions than answers. Nat Immunol. 2008;9:725–32.

    Article  CAS  PubMed  Google Scholar 

  65. Ulrich JN, Naranjo DP, Alimi TO, Müller GC, Beier JC. How much vector control is needed to achieve malaria elimination? Trends Parasitol. 2013;29:104–9.

    Article  PubMed  PubMed Central  Google Scholar 

  66. Gillies MT, Wilkes TJ. A study of the age-composition of populations of Anopheles gambiae Giles and A. funestus Giles in North-Eastern Tanzania. Bull Entomol Res. 1965;56:237–62.

    Article  CAS  PubMed  Google Scholar 

  67. Straif SC, Beier JC. Effects of sugar availability on the blood-feeding behavior of Anopheles gambiae (Diptera: Culicidae). J Med Entomol. 1996;33:608–12.

    Article  CAS  PubMed  Google Scholar 

  68. Takken W, van Loon JJ, Adam W. Inhibition of host-seeking response and olfactory responsiveness in Anopheles gambiae following blood feeding. J Insect Physiol. 2001;47:303–10.

    Article  CAS  PubMed  Google Scholar 

  69. Gary R, Foster W. Diel timing and frequency of sugar feeding in the mosquito Anopheles gambiae, depending on sex, gonotrophic state and resource availability. Med Vet Entomol. 2006;20:308–16.

    Article  PubMed  Google Scholar 

  70. WHO. World Malaria report 2015. Geneva: World Health Organization; 2015.

    Google Scholar 

  71. Otten M, Aregawi M, Were W, Karema C, Medin A, Bekele W, et al. Initial evidence of reduction of malaria cases and deaths in Rwanda and Ethiopia due to rapid scale-up of malaria prevention and treatment. Malar J. 2009;8:14.

    Article  PubMed  PubMed Central  Google Scholar 

  72. Okiro EA, Alegana VA, Noor AM, Snow RW. Changing malaria intervention coverage, transmission and hospitalization in Kenya. Malar J. 2010;9:285.

    Article  PubMed  PubMed Central  Google Scholar 

  73. van Eijk AM, Hill J, Larsen DA, Webster J, Steketee RW, Eisele TP, et al. Coverage of intermittent preventive treatment and insecticide-treated nets for the control of malaria during pregnancy in sub-Saharan Africa: a synthesis and meta-analysis of national survey data, 2009–2011. Lancet Infect Dis. 2013;13:1029–42.

    Article  PubMed  Google Scholar 

  74. Ouattara AF, Dagnogo M, Constant EA, Koné M, Raso G, Tanner M, et al. Transmission of malaria in relation to distribution and coverage of long-lasting insecticidal nets in central Côte d’Ivoire. Malar J. 2014;13:109.

    Article  PubMed  PubMed Central  Google Scholar 

  75. Bhattarai A, Ali AS, Kachur SP, Mårtensson A, Abbas AK, Khatib R, et al. Impact of artemisinin-based combination therapy and insecticide-treated nets on malaria burden in Zanzibar. PLoS Med. 2007;4:e309.

    Article  PubMed  PubMed Central  Google Scholar 

  76. Corbel V, Akogbeto M, Damien GB, Djenontin A, Chandre F, Rogier C, et al. Combination of malaria vector control interventions in pyrethroid resistance area in Benin: a cluster randomised controlled trial. Lancet Infect Dis. 2012;12:617–26.

    Article  PubMed  Google Scholar 

  77. Ross R. Studies on malaria. London: John Murray; 1928.

    Google Scholar 

  78. McKenzie FE, Killeen GF, Beier JC, Bossert WH. Seasonality, parasite diversity, and local extinctions in Plasmodium falciparum malaria. Ecology. 2001;82:2673–81.

    Article  PubMed  PubMed Central  Google Scholar 

  79. Kligler I. Flight of anopheles mosquitoes. Trans R Soc Trop Med Hyg. 1924;18:199–202.

    Article  Google Scholar 

Download references

Authors’ contributions

LZ developed the individual-based model, designed the study, and analysed the results. GCM, WAQ, YS, SFT, and SD contributed to the literature review. GCM, JMM, WMH, and JCB contributed to the study design. KLA contributed to data analysis. All authors contributed to the interpretation. LZ wrote the first draft, and all authors contributed to manuscript writing. All authors read and approved the final manuscript.

Acknowledgements

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

The datasets generated and/or analysed during the current study available from the corresponding author on reasonable request.

Funding

Research reported in this publication was supported by the National Institute of Allergy And Infectious Diseases of the National Institutes of Health under Award Number R01AI100968. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Lin Zhu.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhu, L., Müller, G.C., Marshall, J.M. et al. Is outdoor vector control needed for malaria elimination? An individual-based modelling study. Malar J 16, 266 (2017). https://doi.org/10.1186/s12936-017-1920-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12936-017-1920-y

Keywords