- Open Access
Efficacy and persistence of long-lasting microbial larvicides against malaria vectors in western Kenya highlands
Parasites & Vectorsvolume 11, Article number: 438 (2018)
Chemical-based malaria vector control interventions are threatened by the development of insecticide resistance and changes in the behavior of the vectors, and thus require the development of alternative control methods. Bacterial-based larvicides have the potential to target both insecticide resistant and outdoor-biting mosquitoes and are safe to use in the environment. However, the currently available microbial larvicide formulations have a short duration of activity requiring frequent re-applications which increase the cost of control interventions. This study was designed to evaluate the efficacy and duration of activity of two long-lasting formulations of Bacillus thuringiensis var. israelensis (Bti) and Bacillus sphaericus (Bs) (LL3 and FourStar®) under field conditions in western Kenya highlands.
Three sites were selected for this study in the highlands of western Kenya. In each site, one hundred anopheline larval habitats were selected and assigned to one of three arms: (i) LL3; (ii) FourStar®; and (iii) untreated control larval habitats. Four types of larval habitats were surveyed: abandoned gold mines, drainage canals, fish ponds and non-fish ponds. The habitats were sampled for mosquito larvae by using a standard dipping technique and collected larvae were recorded according to the larval stages of the different Anopheles species. The larvicides were applied at manufacturers’ recommended dosage of 1 briquette per 100 square feet. Both treatment and control habitats were sampled for mosquito larvae immediately before treatment (day 0), and then at 24 hours, 3 days and weekly post-treatment for 5 months.
Overall larval density in treatment habitats was significantly reduced after application of the two microbial larvicides as compared to the control habitats. Post-intervention reduction in anopheline larval density by LL3 was 65, 71 and 84% for 1 day, 2 weeks and 4 weeks, respectively. FourStar® reduced anopheline larval density by 60, 66 and 80% for 1 day, 2 weeks and 4 weeks, respectively. Comparisons between the treatments reveal that LL3 and FourStar® were similar in efficacy. A higher reduction in Anopheles larval density was observed in the abandoned goldmines, while drainage canals had the lowest reduction.
Both LL3 and FourStar® long-lasting microbial larvicides were effective in reducing immature stages of An. gambiae complex and An. funestus group species, with significant reductions lasting for three months post-application.
In recent years, the western Kenya highlands have experienced increased frequencies of malaria epidemics and high transmission rates despite prolonged use of long-lasting insecticidal nets (LLINs) and indoor residual spraying (IRS) for malaria vector control [1,2,3]. Human population increase and changes in land use such as deforestation and swamp cultivation have been linked with this increase in malaria transmission and local epidemics [2, 4]. Moreover, Anopheles arabiensis is increasingly becoming an important malaria vector in the highlands of western Kenya in addition to the well-known An. gambiae (s.s.) and An. funestus. Anopheles arabiensis is difficult to control due to its exophilic and zoophilic tendencies allowing it to survive chemical control better than the highly anthropophagic and endophilic An. gambiae (s.s.) and An. funestus [5, 6]. Moreover, the current chemical-based malaria vector control interventions (LLINs and IRS) are threatened by the development of insecticide resistance and changes in vector biting and resting behavior that require the development of alternative control interventions [7,8,9,10].
Mosquito larval control is once again gaining importance owing to the need for integrated vector management approaches to complement measures against adult mosquitoes [11,12,13,14]. Biolarvicides based on the Gram-positive spore-forming bacteria, Bacillus thuringiensis var. israelensis (Bti) and Bacillus sphaericus (Bs) can target both insecticide resistant and outdoor biting malaria vectors [15,16,17]. Several studies in Africa have demonstrated the effectiveness of these bacterial larvicides in granular, powder and tablet formulations in reducing the density of malaria vectors [15, 18,19,20,21,22,23,24,25,26,27] and malaria transmission [23, 28]. However, these formulations are of short duration of activity and hence require regular applications which are associated with high operational costs as compared to the longer-lasting chemical-based interventions; hence, their use has not been widely embraced by the national malaria control programs in sub-Saharan Africa [29, 30].
To overcome the problem of low residual activity of the microbial larvicides, new long-lasting formulations that release effective levels of Bti and Bs at the water surface over prolonged periods of time have been developed. However, before deployment of these products for malaria vector control, their effectiveness and persistence have to be verified in natural anopheline breeding habitats in different ecological settings. This study was designed to evaluate two new long-lasting microbial larvicide formulations based on Bti and Bs against different stages of malaria vectors in different aquatic larval habitats in western Kenya highlands.
The study was conducted in three villages in the western Kenya highlands, Iguhu (00°09'42"N, 34°45'42"E) in Kakamega County, Emutete (00°01'35"N, 34°37'00"E) and Emakakha (00°06'32"N, 34°39'12"E) in Vihiga County (Fig. 1). All study sites are located in-between Kisumu town and Kakamega town, the two largest cities in western Kenya. Iguhu is located about 15 km south of Kakamega town along the Kisumu-Kakamega road, Emakakha is about 15 km southwest of Iguhu, and Emutete is about 10 km south of Emakakha (Fig. 1). Briefly, the study villages were located in an area of hills and valleys. The hillside mostly comprised comprises maize plantations and a few patches of indigenous forests located along the valleys. The principal occupation of inhabitants of the study villages includes subsistence farmers of maize and some vegetables, and livestock rearing such as cattle, goats, sheep and chickens. The climate consists of a bimodal pattern of rainfall, with a long rainy season from April to June and a short rainy season in October and December. There is no clear dry season in the villages, but usually there is less rainfall from July to September. January and February are the hottest and driest months.
Larval habitat characterization, selection and randomization
The study villages were surveyed for the presence of anopheline larval habitats and habitats were characterized by type, size, permanence, vegetation coverage and land use types as previously classified . All potential larval habitats identified were enumerated and mapped using a hand held global positioning system (GPS) device before being sampled for mosquito larvae and pupae (Fig. 1). Four types of larval habitats were surveyed, namely abandoned gold mines, drainage canals, non-fish ponds and fish ponds. Abandoned gold mines were pits left after the cessation of gold mining activities while drainage canals were canals used to drain water from farms. Non-fish ponds were natural or man-made relative large water bodies with or without emergent aquatic vegetation. Fish ponds were man-made ponds used for fish farming, and some had fish present during the period of survey.
In each of the three villages, 100 anopheline larval breeding habitats were selected for the study and randomized equally into three arms: (i) FourStar®; (ii) LL3; and (iii) untreated control larval habitats. FourStar® is a slow-release briquet formulation with combination of Bti and Bs that the manufacturer claims is effective from 90 to 180 days (Adapco Inc., Sanford, FL, USA). The 30 g formulation consists of 6% by weight of Bacillus sphaericus serotype H5a5b strain 2362, 1% by weight of Bacillus thuringiensis subspecies israelensis, strain BMP 144 and other ingredients that slowly release the bacterial toxins. LL3 has the same Bti/Bs contents as FourStar®, the only difference being that LL3 floats in the water (density approximately 0.99 g/cm3) and FourStar® sinks. The briquets were applied at the manufacturers’ recommended dosage of 1 briquet per 100 square feet of larval habitat though hand broadcasting. For habitats with a likelihood of overflow or flow of water during heavy rains, the briquets were fastened to a long loose thin string attached to a pole at the shallow water margins.
Pre- and post-intervention larval sampling
Selected larval habitats were sampled for mosquito immatures using World Health Organization (WHO) 350 ml standard mosquito dippers. Depending on habitat size, 5 to 20 dips were taken from each larval habitat: 5 dips were undertaken in small larval habitats of ≤ 1 m2, 10 dips for medium-sized habitats (2–15 m2) and 20 dips for relatively large habitats (> 15 m2). The immatures were classified into early instars (1st- and 2nd-stage larvae), late instars (3rd- and 4th-stage larvae) and pupae of different species of Anopheles. Both treatment and control larval habitats were sampled on a weekly basis for five weeks before the intervention to obtain baseline information on abundance of mosquito larvae and pupae. The larval habitats were monitored immediately prior to application (day 0) and then on days 1, 3 and weekly after 7 days post-treatment for a period of 5 months.
Identification of An. gambiae complex and An. funestus group
In the laboratory, mosquito specimens from each habitat were morphologically identified into An. gambiae complex and An. funestus group and stored in 80% ethanol. These were later identified as respective sibling species by polymerase chain reaction (PCR) as described previously [32, 33]. In brief, for members of An. gambiae complex, PCR reactions were conducted in a final volume of 20 μl consisting of 0.25 μM of each of the five primers, 1:1 DreamTaq Green PCR Master Mix (ThermoScientific, Waltham, USA) and 2 μl of DNA extract. The samples were amplified in Bio-Rad T100 thermo cycler (Bio-Rad, Hercules, CA, USA) and cycling conditions were 95 °C for 5 min followed by 30 cycles of denaturation at 94 °C for 30 s, annealing at 50 °C for 30 s, extension at 72 °C for 30 s and final extension at 72 °C for 10 min.
For An. funestus sibling species identification, each PCR run was conducted in a final volume of 25 μl consisting of 0.4 μM of each of the six primers, 1:1 DreamTaq Green PCR Master Mix and 3 μl of extracted DNA. The samples were amplified in a Bio-Rad T100 thermocycler and cycling conditions were 95 °C for 5 min followed by 40 cycles of denaturation at 94 °C for 30 s, annealing at 50 °C for 30 s, extension at 72 °C for 40 s and final extension at 72 °C for 10 min. The amplified DNA for both sibling species of An. gambiae and An. funestus species complexes were separated based on their fragment size by gel electrophoresis.
The percentage reduction in the larval and pupal densities was calculated using Mulla’s formula  as follows: % reduction = 100 – (C1/T1 × T2/C2) × 100, where C1 and C2 are the counts in control habitats before and after treatment, respectively, and T1 and T2 are the counts in treated habitats before and after treatment, respectively. The differences in abundance of various stages of immature Anopheles observed were compared using generalized estimating equations (GEE) based on the Poisson distribution assumption in which baseline and observation time (in weeks) were treated as covariates [34, 35]. The values of the covariates were constant for the repeated elementary observations at each habitat. GEE is considered as a suitable model that handles longitudinal data which is potentially correlated between subsequent observations and thus distorts the independence assumption of the ordinary linear regression model. The correlation between longitudinal observations was tested against four assumptions, i.e. independent (not correlated), exchangeable (fixed correlation), lag 1 autoregression and unstructured (correlations are all different). The Akaike information criterion was used for model parameter estimation and selection. The significant level was P = 0.05. If 0.05 < P < 0.10, then we defined it as marginally significant, and insignificant otherwise. The models were first run using the interventions against control to evaluate the impact of interventions on the abundance of different species of Anopheles immatures, then the two interventions against each other (LL3 versus FourStar®) to determine the difference between the two larvicides. Relative reduction in abundance of immature anophelines was calculated as reduction in intervention habitats against control habitats based on pre-intervention observations. Data analysis was conducted by the use of open source language R 3.3.1. For the GEE analysis, the geeglm function in the geepack package was used.
During the 25-week study period, 7,896 samples (1,486 pre-intervention and 6,410 post-intervention) were collected. A total of 43,147 immature anophelines were recorded. Among these, 29,221 were An. gambiae (s.l.) and 13,926 were An. funestus (s.l.). During the 5-week pre-intervention period, 49% of habitats had immature An. gambiae with an average density of 5.1 individuals per habitat per sampling occasion, and 37% of habitats had immature An. funestus with an average density of 2.0 individuals per habitat per sampling occasion. During the 20-week post-intervention period, 27% of habitats had immature An. gambiae with an average density of 3.4 individuals per habitat per sampling occasion, 32% of habitats had immature An. funestus with an average density of 1.7 individuals per habitat per sampling occasion.
Impact of FourStar® and LL3 on general Anopheles larval density
During the pre-intervention period, the mean density of Anopheles larvae was 8.05, 7.56 and 7.44 larvae per dip for LL3, FourStar® and control arms, respectively (Fig. 2). The mean density of Anopheles larvae per dip in treatment habitats (LL3 and FourStar®) and control larval habitats was not significantly different during the pre-intervention period (GEE, P > 0.1, Table 1). However, during 20 weeks of post-intervention observation, the mean density of Anopheles larvae was significantly reduced to 2.88 and 3.09 for LL3- and FourStar®-treated habitats (Fig. 2), respectively, when compared to the baseline pre-treatment values (GEE, P < 0.05, Table 1). On the other hand, Anopheles larvae density in control habitats was slightly increased to 7.77 larvae per dip. Post-intervention reduction in anopheline larval density by LL3 was 65, 71 and 84% for 24 hours, 1 week and 4 weeks, respectively. At the same time, FourStar® reduced anopheline larval density by 60, 66 and 80% (Fig. 3, Table 2). By the 5–7 weeks post-application, reduction in larval density averaged 60 and 44% for LL3 and FourStar®, respectively. The reduction in larval density by LL3 remained above 50% for up to 12 weeks (Fig. 3). Overall, the two interventions, LL3 and FourStar® did not show significant difference in reducing immature mosquito densities (Table 2).
Impact of FourStar® and LL3 on density of different larval stages of An. gambiae complex and An. funestus group
During the baseline, the mean density of An. gambiae complex 1st- and 2nd-instar larvae was 3.17, 3.07 and 3.18 larvae per dip for LL3-treated, FourStar®-treated and control habitats, respectively (Fig. 4). The mean larval density for early and late instars of An. gambiae complex in the three experimental arms did not differ significantly during this period (GEE, P > 0.1, Table 1). However, within the 20-week post-intervention period, the mean density of An. gambiae complex early instars was significantly reduced to 1.35 and 1.47 for LL3- and FourStar®-treated habitats, respectively, whereas in control habitats the mean density did not change significantly (GEE, P > 0.1). On the other hand, the mean density of late instars of An. gambiae complex in the baseline survey was 1.95, 2.12 and 1.95 larvae per dip, for LL3-treated, FourStar®-treated and control habitats, respectively. However, during the 20-week post-intervention period this mean density was significantly reduced to 0.33 and 0.51 larvae per dip (GEE, P < 0.05), whereas in control habitats, the mean density did not change significantly (GEE, P > 0.1). Likewise, the mean pupal density for An. gambiae complex in the pre-intervention period was 0.28, 0.26 and 0.26 pupae per dip in LL3-treated, FourStar®-treated and control larval habitats, respectively. During this period, the mean pupal density between the control and treatment habitats was not significantly different (GEE, P > 0.1). However, during the post-treatment period, the mean pupal density for An. gambiae complex was significantly reduced to 0.03 and 0.03 pupae per dip in LL3- and FourStar®-treated habitats, respectively (GEE, P < 0.05). In all cases, the two interventions, LL3 and FourStar® did not show significant differences in immature mosquito densities (Table 2).
For the An. funestus group, the mean larval density for early instars during the baseline period was 1.40, 1.00 and 1.11 larvae per dip in LL3-treated, FourStar®-treated and control habitats, respectively, while for late instars it was 1.04, 0.88 and 0.78 larvae per dip, respectively (Fig. 4). During the 20-week post-intervention period, the mean larval density of An. funestus group was significantly reduced in treatment habitats to 0.84 and 0.77 larva per dip for early instars and 0.31 and 0.29 larvae per dip for late instars in LL3- and FourStar®-treated larval habitats, respectively (GEE, P < 0.05). The mean larval density for early and late instars of An. funestus group in control habitats did not change significantly from the baseline values, and was maintained at 1.38 and 0.82 larvae per dip, respectively (GEE, P > 0.1). The mean pupal density for An. funestus group species in the pre-intervention period was 0.21, 0.25 and 0.17 pupae per dip in LL3-treated, FourStar®-treated and control habitats, respectively. The mean pupal density in treated and control habitats pre-treatment did not differ significantly (GEE, P > 0.1). However, during the post treatment period, there was only a marginal reduction in An. funestus group mean pupal density in both the LL3 and FourStar® treatments (Fig. 4). In all cases, the two interventions, LL3 and FourStar® did not show significant difference in reducing immature mosquito densities (Table 2).
Impact of FourStar® and LL3 in different larval habitat types
During the study, four types of aquatic habitats were found to be common breeding sites of Anopheles larvae: abandoned goldmines, drainage canals, non-fish ponds and fish ponds (Fig. 5). In the abandoned goldmines, mean Anopheles larval density in the control habitats during the pre-intervention period was 12.25 and increased to 15.17 larvae per dip during the intervention period. In the treated habitats, the baseline mean Anopheles larvae density was 12.69 and 12.76 for LL3- and FourStar®-treated habitats, respectively (Fig. 5). However, during the post-intervention period, this decreased significantly to an average of 5.61 and 5.73 larvae per dip for LL3- and FourStar®-treated habitats, respectively (GEE, P < 0.05). This reduction in the mean Anopheles larval density remained significantly low compared to the baseline pre-intervention values for up to 12 weeks.
For drainage canals, the Anopheles larval productivity was reduced significantly after larviciding from the baseline pre-intervention values of 4.87 to 1.57 larvae per dip after treatment in LL3 and 5.46 to 1.55 larvae per dip after treatment in FourStar®-treated habitats (GEE, P < 0.05). The trend was similar in the fish ponds, where the pre-intervention period Anopheles larvae density was 6.79 and 6.05 larvae per dip for LL3- and FourStar®-treated habitats, respectively, but decreased significantly to 1.75 and 2.34 larvae per dip, respectively, during the post-intervention period (GEE, P < 0.05). In the non-fish ponds, mean Anopheles larval density during the pre-intervention period was 5.56, 6.84 and 6.77 in control, LL3- and FourStar®-treated habitats, respectively. Twenty-four hours post-intervention period mean Anopheles larval density decreased significantly to 1.84 and 2.31 larvae per dip in LL3- and FourStar®-treated habitats, respectively (GEE, P < 0.05), while in the control the density was not significantly different from the baseline value (Fig. 5). Overall, the highest reduction in Anopheles larval density was observed in the abandoned goldmines, while drainage canals had the lowest reduction.
Species identity of larval specimens of An. gambiae complex and An. funestus group
A total of 366 An. gambiae complex and 237 An. funestus group larval specimens were morphologically identified and then processed for sibling species identity using PCR technique. In the An. gambiae complex, An. gambiae (s.s.) were the predominant sibling species, accounting for 70.8% of the identified specimens. On the other hand, An. funestus (s.s.) accounted for 73.6% of the identified sibling species of the An. funestus group (Table 3).
Mosquito larval control has been found to be effective in lowering malaria transmission as it targets both indoor and outdoor biting mosquito species [23, 36]. Larviciding with microbial larvicides has proved to be effective in mosquito control but with low residual activity requiring repeated applications with high cost implications [29, 30]. Recently, long-lasting microbial larvicide formulations with potential for sustained larval control have been developed to overcome the problem of low persistence, thereby improving cost-effectiveness. The current study was designed to test two long-lasting formulations of microbial larvicides (FourStar® and LL3) against different stages of malaria vectors in different aquatic larval habitats in the highlands of western Kenya. In the study area, the malaria vectors An. gambiae (s.s.) and An. funestus (s.s.) were the predominant sibling species of the An. gambiae complex and An. funestus group, respectively.
Results showed that both FourStar® and LL3 long-lasting microbial larvicides were effective in reducing the density of Anopheles larval populations, as shown by the significant differences observed between the treated and control larval habitats during the post-intervention period. This was not the case during the baseline survey period where the larval density between the treatment and control habitats was broadly similar, but this trend changed markedly after application of FourStar® and LL3. The observed fluctuations in Anopheles larval densities after application of microbial larvicides could be due to the effect of heavy rains experienced during the study period. The heavy rains results in flushing of the floating toxins and thereby reducing the impact of the larvicides. The impact of both larvicides in reducing Anopheles larval densities was clearly observed during the first twelve weeks post-application and this is consistent with another study performed previously in the same area .
The findings have shown that the two long-lasting microbial larvicides tested were equally effective in reducing the density of both the early and late larval stages of An. gambiae complex and An. funestus group. This reduction was, however, higher in species of An. gambiae complex as compared to An. funestus group. This could be due to the fact that An. gambiae complex were found in open sunlit temporary to semi-permanent habitats with less vegetation cover where the larvicidal toxins can disperse better with high chances of being consumed by the target larvae. For the case of An. funestus group, larval habitats were more permanent with high levels of vegetation which may interfere with the spread of the toxins and hence reduce their availability to the target larvae.
In the current study, reductions in pupal densities were significant in An. gambiae complex while in An. funestus group they were only marginal. Reduced activity of long-lasting microbial larvicides against pupae of An. funestus group observed in the current study agreed with findings of a recent previous study performed in the same area . Since a reduction in pupal densities is a good indicator of reduced adult emergence, low reductions in pupal of An. funestus group may indicate low activity of the long-lasting microbial larvicides against this important malaria vector. Analysis of larval and pupal reductions in different larval habitat types studied indicated a reduced persistence of long-lasting microbial larvicides in some larval habitats. The larvicides were more effective in reducing larval densities in the abandoned goldmines, non-fish ponds and fish ponds, which were mainly stagnant water breeding habitats. Efficacy was, however, relatively low in drainage canals and this may be due to the fact that drainage may allow water to flow during the rains and wash away floating larvicide toxins and hence reduce the effectiveness of the intervention. This finding is in line with a study at the Kenyan coast  which showed reduced persistence of Bti/Bs in drainage canals. Despite the observed effectiveness of LL3 and FourStar® larvicides in reducing larvae of malaria vectors, their activity varied with vector species and habitat types. Thus, based on our findings, habitat type is an important criterion when considering LL3 and FourStar® for integrated vector management operations. The study demonstrates that the two larvicides are very target specific and environmentally benign and can be applied in various water bodies without fear of detrimental effects on non-target organisms.
This study showed that both LL3 and FourStar® long-lasting microbial larvicides were effective in reducing larvae and pupae of An. gambiae complex and An. funestus group mosquitoes for more than three months. Based on marginal levels of immature reduction observed, these long-lasting microbial larvicides are not suitable for use as stand-alone interventions, and should be integrated with the existing tools for malaria vector control.
Bacillus thuringiensis var. israelensis
Generalized estimating equations
Global positioning system
Indoor residual spraying
Long-lasting insecticidal nets
Polymerase chain reaction
World Health Organization
Githeko AK, Ndegwa W. Predicting malaria epidemics in the Kenyan highlands using climate data: a tool for decision makers. Glob Chang Hum Health. 2001;2:54–63.
Himeidan YE, Kweka EJ. Malaria in East African highlands during the past 30 years: impact of environmental changes. Front Physiol. 2012;3:315.
Zhou G, Afrane YA, Vardo-Zalik AM, Atieli H, Zhong D, Wamae P, et al. Changing patterns of malaria epidemiology between 2002 and 2010 in western Kenya: the fall and rise of malaria. PLoS One. 2011;6:e20318.
Munga S, Minakawa N, Zhou G, Githeko AK, Yan AG. Survivorship of immature stages of Anopheles gambiae s.l. (Diptera: Culicidae) in natural habitats in western Kenya highlands. J Med Entomol. 2007;44:758–64.
Govella NJ, Ferguson H. Why use of interventions targeting outdoor biting mosquitoes will be necessary to achieve malaria elimination. Front Physiol. 2012;3:199.
Githeko AK, Ototo EN, Guiyun Y. Progress towards understanding the ecology and epidemiology of malaria in the western Kenya highlands: opportunities and challenges for control under climate change risk. Acta Trop. 2012;121:19–25.
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.
Cooke MK, Kahindi SC, Oriango RM, Owaga C, Ayoma E, Mabuka D, et al. “A bite before bed”: exposure to malaria vectors outside the times of net use in the highlands of western Kenya. Malar J. 2015;14:259.
Mathias DK, Ochomo E, Atieli F, Ombok M, Bayoh NM, Olang G, et al. Spatial and temporal variation in the kdr allele L1014S in Anopheles gambiae s.s. and phenotypic variability in susceptibility to insecticides in western Kenya. Malar J. 2011;10:10.
Ochomo E, Bayoh NM, Kamau L, Atieli F, Vulule J, Ouma C, et al. Pyrethroid susceptibility of malaria vectors in four districts of western Kenya. Parasit Vectors. 2014;7:310.
Bukhari T, Takken W, Koenraadt CJM. Biological tools for control of larval stages of malaria vectors - a review. Biocontrol Sci Technol. 2013;23:987–1023.
Fillinger U, Lindsay SW. Larval source management for malaria control in Africa: myths and reality. Malar J. 2011;10:353.
Gu W, Novak RJ. Habitat-based modeling of impacts of mosquito larval interventions on entomological inoculation rates, incidence, and prevalence of malaria. Am J Trop Med Hyg. 2005;73:546–52.
Gu W, Utzinger J, Novak RJ. Habitat-based larval interventions: a new perspective for malaria control. Am J Trop Med Hyg. 2008;78:2–6.
Fillinger U, Knols BGJ, Becker N. Efficacy and efficiency of new Bacillus thuringiensis var. israelensis and Bacillus sphaericus formulations against Afrotropical anophelines in Western Kenya. Trop Med Int Health. 2003;8:37–47.
Lacey LA, Merritt RW. The safety of bacterial microbial agents used for black fly and mosquito control in aquatic environments. In: Hokkanen HMT, Hajek AE, editors. Environmental impact of microbial insecticides: Need and methods for risk assessment. Dordrecht, Netherlands: Kluwer Academic Publishers; 2003. p. 151–68.
Mittal PK. Biolarvicides in vector control: challenges and prospects. J Vector Borne Dis. 2003;40:20–32.
Djènontin A, Pennetier C, Zogo B, Soukou KB, Ole-Sangba M, Akogbéto M, et al. Field efficacy of vectobac GR as a mosquito larvicide for the control of anopheline and culicine mosquitoes in natural habitats in Benin, West Africa. PLoS One. 2014;9:e87934.
Dambach P, Louis V, Kaiser A, Ouedraogo S, Sie A, Sauerborn R, et al. Efficacy of Bacillus thuringiensis var. israelensis against malaria mosquitoes in northwestern Burkina Faso. Parasit Vectors. 2014;7:371.
Majambere S, Lindsay SW, Green C, Kandeh B, Fillinger U. Microbial larvicides for malaria control in the Gambia. Malar J. 2007;6:76.
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.
Kandyata A, Mbata KJ, Shinondo CJ, Katongo C, Kamuliwo RM, Nyirenda F, et al. Impacts of Bacillus thuringiensis var. israelensis and Bacillus sphaericus insect larvicides on mosquito larval densities in Lusaka, Zambia. Med J Zambia. 2012;39:33–8.
Geissbühler Y, Kannady K, Chaki PP, Emidi B, Govella NJ, Mayagaya V, et al. Microbial larvicide application by a large-scale, community-based program reduces malaria infection prevalence in urban Dar Es Salaam, Tanzania. PLoS One. 2009;4:e5107.
Kahindi SC, Midega JT, Mwangangi JM, Kibe LW, Nzovu J, Luethy P, et al. Efficacy of vectobac DT and culinexcombi against mosquito larvae in unused swimming pools in Malindi, Kenya. J Am Mosq Control Assoc. 2008;24:538–42.
Fillinger U, Lindsay SW. Suppression of exposure to malaria vectors by an order of magnitude using microbial larvicides in rural Kenya. Trop Med Int Health. 2006;11:1629–42.
Mwangangi JM, Kahindi SC, Kibe LW, Nzovu JG, Luethy P, Githure JI, et al. Wide-scale application of Bti/Bs biolarvicide in different aquatic habitat types in urban and peri-urban Malindi, Kenya. Parasitol Res. 2011;108:1355–63.
Mpofu M, Becker P, Mudambo K, De Jager C. Field effectiveness of microbial larvicides on mosquito larvae in malaria areas of Botswana and Zimbabwe. Malar J. 2016;15:586.
Zhou G, Afrane YA, Dixit A, Atieli HE, Lee MC, Wanjala CL, et al. Modest additive effects of integrated vector control measures on malaria prevalence and transmission in western Kenya. Malar J. 2013;12:256.
Fillinger U, Ndenga B, Githeko A, Lindsay SW. Integrated malaria vector control with microbial larvicides and insecticide-treated nets in western Kenya: a controlled trial. Bull World Health Organ. 2009;87:655–65.
Worrall E, Fillinger U. Large-scale use of mosquito larval source management for malaria control in Africa: a cost analysis. Malar J. 2011;10:338.
Kweka EJ, Zhou G, Lee M, Gilbreath TM, Mosha F, Munga S, et al. Evaluation of two methods of estimating larval habitat productivity in western Kenya highlands. Parasit Vectors. 2011;4:110.
Scott JA, Brogdon WG, Collins FH. Identification of single specimens of the Anopheles gambiae complex by the polymerase chain reaction. Am J Trop Med Hyg. 1993;49:520–9.
Koekemoer LL, Kamau L, Hunt RH, Coetzee M. A cocktail polymerase chain reaction assay to identify members of the Anopheles funestus (Diptera: Culicidae) group. Am J Trop Med Hyg. 2002;66:804–11.
Gosho M. Criteria to select a working correlation structure for the generalized estimating equations method in SAS. J Stat Software. 2014;57:1–10.
Li C-S. Generalized estimating equations in longitudinal data analysis: a review and recent developments. Adv Stat. 2014;2014:303728.
Afrane YA, Mweresa NG, Wanjala CL, Gilbreath TM, Zhou G, Lee MC, et al. Evaluation of long-lasting microbial larvicide for malaria vector control in Kenya. Malar J. 2016;15:577.
The help of the field assistants from Iguhu, Emakakha and Emutete villages and the highland malaria program field team is highly appreciated for data collection throughout the study. This paper is published with the permission of the Director General, Kenya Medical Research Institute.
This work was supported by grants from the National Institutes of Health (R01 A1050243, U19 AI129326 and D43 TW001505). The funding body had no role in the design, data collection, analysis and interpretation of data and in writing the manuscript.
Availability of data and materials
All the relevant data for this article are included in the text, tables and figures.
Ethics approval and consent to participate
The study was approved by the Scientific and Ethical Unit of the Kenya Medical Research Institute (Ref: KEMRI/RES/7/3/1). Before field work, meetings were held with County leaders of the study sites to inform them about the study and to seek their cooperation. Oral informed consent was sought and obtained from farm owners before beginning larval habitat surveys and application of the microbial larvicides.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.