Evaluating the impact of larviciding with Bti and community education and mobilization as supplementary integrated vector management interventions for malaria control in Kenya and Ethiopia
Malaria Journal volume 19, Article number: 390 (2020)
Malaria prevention in Africa is mainly through the use of long-lasting insecticide treated nets (LLINs). The objective of the study was to assess the effect of supplementing LLINs with either larviciding with Bacillus thuringiensis israelensis (Bti) or community education and mobilization (CEM), or with both interventions in the context of integrated vector management (IVM).
The study involved a factorial, cluster-randomized, controlled trial conducted in Malindi and Nyabondo sites in Kenya and Tolay site in Ethiopia, to assess the impact of the following four intervention options on mosquitoes and malaria prevalence: LLINs only (arm 1); LLINs and Bti (arm 2); LLINs and CEM (arm 3); and, LLINs combined with Bti and CEM (arm 4). Between January 2013 and December 2015, CDC light traps were used to sample adult mosquitoes during the second, third and fourth quarter of each year in 10 houses in each of 16 villages at each of the three study sites. Larvae were sampled once a fortnight from potential mosquito-breeding habitats using standard plastic dippers. Cross-sectional malaria parasite prevalence surveys were conducted involving a total of 11,846 primary school children during the 3-year period, including 4800 children in Tolay, 3000 in Malindi and 4046 in Nyabondo study sites.
Baseline relative indoor anopheline density was 0.11, 0.05 and 0.02 mosquitoes per house per night in Malindi, Tolay and Nyabondo sites, respectively. Nyabondo had the highest recorded overall average malaria prevalence among school children at 32.4%, followed by Malindi with 5.7% and Tolay 1.7%. There was no significant reduction in adult anopheline density at each of the three sites, which could be attributed to adding of the supplementary interventions to the usage of LLINs. Malaria prevalence was significantly reduced by 50% in Tolay when using LLINs coupled with application of Bti, community education and mobilization. The two other sites did not reveal significant reduction of prevalence as a result of combining LLINs with any of the other supplementary interventions.
Combining LLINs with larviciding with Bti and CEM further reduced malaria infection in a low prevalence setting in Ethiopia, but not at sites with relatively higher prevalence in Kenya. More research is necessary at the selected sites in Kenya to periodically determine the suite of vector control interventions and broader disease management strategies, which when integrated would further reduce adult anopheline populations and malaria prevalence beyond what is achieved with LLINs.
After more than two decades of malaria being associated with an estimated one million global deaths annually prior to 2000, the disease cases declined by 40–60% in many countries between 2000 and 2015 . Malaria deaths stood at 435,000 globally in 2017, far below the 2000 levels despite a slowdown in the rate of decline reported from 2015 onwards . The decrease in malaria burden in nearly all the endemic countries since 2000 is commonly attributed to a multi-pronged strategy involving the scaling up of the use of long-lasting insecticide-treated nets (LLINs), selective indoor residual spraying (IRS) and better access to prompt diagnosis and effective treatment of the disease [3,4,5], among other programmatic interventions. Other factors that have been speculated as generally having contributed to a decline of malaria vectors and malaria itself in endemic countries include climatic and environmental changes and general development leading to better socio-economic conditions including improved housing [6,7,8].
Unfortunately, as different countries pursue the goal of malaria elimination, encouraged by the gains experienced since 2000, serious concerns have been expressed at national and international level regarding sustainability of the primary interventions, particularly those targeting vector control. In the first instance, there has been a looming threat posed to malaria control programmes by widespread vector resistance to the insecticides commonly used in IRS and LLINs [9, 10]. Secondly, the ability of governments in many malaria-affected countries to finance prevention and control activities outside the ambit of external donor-funded programmes is often in doubt . In recognition of these, among other technical, community and health system challenges, the World Health Organization (WHO) has repeatedly recommended the adoption of integrated vector management (IVM) among the strategies that could lead to sustainable malaria control and, ultimately, its elimination [12, 13]. Key elements of IVM include evidence-based decision making, the integration of non-chemical and chemical vector control methods, advocacy and social mobilization, and inter-sectoral collaboration and action. The key elements have recently been re-emphasised and endorsed by the World Health Assembly (WHA) in the form of pillars of action in the Global Vector Control Response framework for 2017–2030 .
However, research evidence is largely unavailable to support the perception of greater impact of IVM on malaria when compared to relatively non-integrated vector control in different eco-epidemiological settings. Towards this end, several studies have been conducted to assess whether adding a secondary anti-vector intervention in a setting where the primary vector control method is the usage of LLINs can lead to further site-specific decline of malaria prevalence towards elimination. The secondary methods most generally evaluated include IRS, screening of houses with fine wire mesh to reduce mosquito entry and consequently reducing vector-host contact, and larval source management, particularly involving biolarviciding with Bacillus thuringiensis israelensis (Bti) [15,16,17,18]. The relevance and urgency of such research has received new impetus following a stalling in the general decline of malaria noted from around 2015 .
IVM for malaria control has previously been promoted in areas where the present studies were conducted in Kenya and Ethiopia through collaborative research and non-research developmental activities . The IVM strategy adopted in those previous studies involved community education and mobilization geared towards improved environmental management by the local community and other stakeholders, routine application of the common biopesticide Bti in potential mosquito breeding sites, and the promotion of proper use of LLINs [20,21,22,23,24]. An external review conducted in 2012 to evaluate this earlier work noted that while participatory IVM may have led to significant reductions in malaria in the study sites, it was nevertheless not possible to tease out and quantify the incremental effects of the IVM effort, over and above those ordinarily due to the conventional singular usage of LLINs [19, 25].
The objective of this study was therefore to determine the impact of adding larviciding using Bti and community education and mobilization (CEM) onto the frontline intervention of LLINs. It specifically aimed at verifying under field conditions if adding the former two interventions, herein, also referred to as supplementary interventions, either individually or combined can further reduce indoor malaria vector density and malaria prevalence in the community compared to when LLINs alone are used.
Study areas and context
The studies were conducted from January 2013 to December 2015 simultaneously in three different geographic locations: the Nyabondo plateau in the western part of Kenya, Malindi Sub-County at the Kenyan coast, and the Tolay region in southwestern Ethiopia (Fig. 1). The three study sites were selected to enable comparisons of the likely effectiveness and impacts of IVM to be made between situations of low and moderate to high malaria prevalence, represented, respectively, by Tolay on the one hand, and Malindi and Nyabondo on the other [3, 4].
Nyabondo is a rural setting in a plateau located in Upper Nyakach Division of Kisumu County (0° 23′ S; 34° 58′ E), about 30 km northeast of Lake Victoria, at an altitude between 1520 and 1670 m above sea level. An estimated 34,000 people live in homesteads spread out in the area but administratively grouped into villages, each of up to about 200 houses. The main livelihood activities include subsistence farming, mainly of maize, small-scale livestock rearing and brick-making. Water accumulation in shallow ground pits from the brick-making activity contributes to mosquito breeding . Other breeding sites commonly recorded in Nyabondo include abandoned and poorly managed fish ponds . Malaria is endemic in the Lake Victoria region, with a reported average prevalence of 27% in 2015 . This is the highest average parasite rate among the different malaria-endemic regions in Kenya. Malaria peak seasons follow the long and short rains but this may vary from year to year. Vectors of malaria in Nyabondo mainly belong to the species Anopheles arabiensis . Local houses are constructed of mud walls, iron sheet roofs and have narrow open eaves .
Malindi Sub-County is located in Kenya’s coastal area about 108 km north of Mombasa (3° 13′ S; 40° 7′ E) at an altitude of between 40 and 400 m above sea level. This study was conducted in a rural area of Malindi with approximately 50,000 people. Subsistence farming, fishing and trading are the major economic activities in this area. Malaria is one of the key causes of morbidity and mortality within the Sub-County with an average prevalence of 4–8% . The main vectors of malaria in Malindi are An. arabiensis, Anopheles merus and Anopheles funestus . There are two main rainy seasons in the Kenyan study areas, i.e., the long (March to June) and short (September to November) rains. Common mosquito breeding sites include shallow ponds, drainage channels and wells for domestic water use. Houses in the study area are constructed of mud walls, palm thatched roofs and have wide open eaves.
Tolay is a rural semi-arid, agro-climatic zone located in the southwestern part of Ethiopia (8° 14′ N; 37° 35′ E) at an altitude of 1100–1600 m above sea level. The community is mainly engaged in small-scale mixed agriculture and trade at village level. The major food crops in the region are maize, sorghum and teff. Malaria is unstable and seasonal, with a prevalence of about 2% or less . Peak malaria transmission and epidemic season includes the period from October to November, coinciding with the peak rain season. The main vector of malaria in Tolay is An. arabiensis . Common mosquito-breeding sites include semi-permanent and permanent ponds, pools at the edge of streams and drainage ditches. Most houses in Tolay are constructed of wooden walls plastered with mud and cow dung, and grass thatch or corrugated iron sheet roofs. The houses do not have open eaves.
LLINs (experimental control)
All the three project sites were in areas with universal coverage with LLINs, distributed by the countries’ national malaria control programmes (NMCPs) via multiple channels, including free mass and continuous distributions mainly through antenatal and immunization services [3, 4, 28]. Mass distribution of LLINs had been conducted in the respective study sites in Kenya and Ethiopia in 2012 [3, 4]. Universal coverage implies that people at risk of malaria own LLINs at a ratio of one net for every two persons sleeping in a household. Thus, the existing situation of LLIN coverage constituted the default experimental control treatment to which one other mosquito control (larviciding with Bti) and one social intervention (community education and mobilization) were added either individually or in combination. Survey information available retrospectively regarding household LLIN ownership at the beginning of the baseline data collection in 2013 indicates household net ownership of 83% in Nyabondo area, 84% in Malindi and 64% in Tolay [4, 23, 29].
Larviciding with Bti
Application of Bti granules (VectoBac G; Valent Biosciences Corp, Libertyville, IL, USA) was done once every 2 weeks to all potential Anopheles breeding habitats, identified in advance at a particular study site where this intervention was to be implemented. All habitats found to have anopheline mosquito larvae within a study village were treated. Also included were those habitats within a 200-m buffer zone from the last house in the village in order to mitigate potential spillover effects in entomological assessments due to mosquito flight range . The Bti was applied by hand, by a team of 6–8 people, half of whom were project field staff and the other half members of the community, commonly referred to as mosquito scouts [19, 25], who had been trained locally by the project staff on how to apply the product. The Bti granules were broadcast by hand within the manufacturer’s recommended dose of about 3.0 kg/ha .
Community education and mobilization
The community education and mobilization intervention entailed conducting intensive community educational activities deliberately planned to surpass the awareness creation that normally accompanies LLIN distribution by the Ministry of Health (MOH) of either country. In all three project sites, resident project field staff of about 2–3 persons conducted the CEM activities on regular monthly basis assisted by about 6–10 mosquito scouts. The scouts had previously been trained by the research team in basic aspects of IVM including: visually identifying presence of mosquito larvae in aquatic habitats, larval source management like draining of stagnant water and filling up of temporary pools of water with soil or levelling them wherever appropriate, proper care and usage of mosquito nets, adult mosquito identification, and community participation in malaria prevention and control. The following methods of community education and mobilization were applied: door to door campaigns, engagement of primary school pupils through school health and environment clubs, distribution of information, education and communication (IEC) materials, neighbourhood clean-up campaigns with emphasis on elimination of mosquito breeding and resting habitats, and community meetings on mosquito and malaria control activities.
A cluster-randomized, controlled trial with a factorial design [32, 33] was implemented between January 2013 and December 2015 to study the effects of four intervention treatments (arms) at each of the three study sites (Fig. 2). In each site there were 16 villages (clusters), with four villages randomly assigned to each of the four study arms: Arm 1 comprised LLIN use only; Arm 2 comprised LLIN use and application of Bti; Arm 3 comprised LLIN use and community education and mobilization; and, Arm 4 comprised LLIN use, application of Bti and community education and mobilization. In each village at each site, 10 houses were randomly included in entomological assessment using indoor adult anopheline mosquito relative density as the proxy indicator of malaria transmission risk. A total of 38 schools located within the study sites (10 schools in Malindi, 12 in Nyabondo, 16 in Tolay) were randomly sampled for cross-sectional malaria parasitological surveys among school children, using parasite prevalence as the measure of epidemiological impact. The overriding hypothesis was that integrating usage of LLINs with CEM and larviciding with Bti (arm 4) has greater impact in reducing populations of anopheline mosquitoes and malaria prevalence than when LLINs are used alone or integrated with only either one of the other two supplementary interventions.
The study villages at each study site were randomly selected from a larger pool of villages, which were purposively chosen from the entire list of villages in each study site, and fulfilled most or all of the following criteria: ecologically similar including, among other features, the presence of natural or man-made potential mosquito breeding habitats; location in an area with reported universal MOH-LLIN coverage; accessibility throughout the year; availability of a local primary school; existence of a local health facility such as dispensary or hospital. During selection of the 16 villages, a minimum distance of approximately 0.5 km was maintained between study households in one village and those in the next in order to mitigate potential spillover of intervention effects [30, 34].
Larval mosquito sampling
Mosquito larvae were sampled once every 2 weeks from all the study villages. All stagnant water bodies within and about 200 m from the last house in a village were inspected for mosquito larvae using a standard dipper (11.5 cm diameter and 350 ml capacity) . Larval sampling was carried out between 09:00 and 12:00 by mosquito scouts from the study villages and trained by the research team. The number of scoops varied according to the size of micro-habitat and ranged from one to 10 scoops per habitat for small and large habitats, respectively. The collected larvae were identified as either anopheline or culicine genera using morphological identification keys  and their number recorded. Information on mosquito breeding habitat types was recorded.
Adult mosquito sampling
Adult mosquitoes were sampled indoors once every 3 months from a total of 160 houses at each of the three project sites. The total number of houses comprised of 10 houses from each of the 16 villages, initially selected using simple random sampling. These same houses continued to be used for the rest of the study for ease of field logistics. Consent of the house owners was sought in order to make the houses accessible for the duration of the study. Mosquitoes were collected indoors using standard Centers for Disease Control (CDC) light traps (Model 512; John W. Hock Co, Gainesville, FL, USA). In each selected house, the trap was hung about 1.5 m from the floor, at the rear end of a regularly used bed, with the occupant protected with a bed net . Mosquito scouts deployed the traps at 18:00 and collected them the following morning at 06:00. Mosquitoes collected were sorted in the field laboratory and morphologically identified as belonging to either anopheline or culicine genera. Anophelines were further identified to distinguish between the main malaria vector, Anopheles gambiae sensu lato (s.l.) species complex and other Anopheles species . Finally, sub-samples of An. gambiae s.l. were analysed for sibling species composition using recombinant deoxyribonucleic acid (rDNA)-polymerase chain reaction (PCR) technique [38,39,40]. The relative density of mosquitoes in a village was calculated by dividing the total number collected by the number of houses sampled in order to obtain the average number per house per night.
Cross-sectional malaria surveys
Cross-sectional malaria parasite surveys in children aged between 6 and 12 years were conducted at each project site to assess the local prevalence of malaria. A total of 38 public primary schools were used in the study, including 12 schools in Nyabondo, 10 in Malindi and 16 in Tolay. The schools were intended to be matched with the study villages at a ratio of one school per village. However, this was possible in Tolay but not in Malindi and Nyabondo. For this reason, some of the study arms at the latter two sites had fewer than four schools included in the prevalence assessments in spite of them having four study villages. In each school, 20 children of equal number of males and females were randomly selected from each of classes 2–6 based on the children present that day, and using computer-generated random number tables. This sample size of 100 children per school sought to estimate a 5% change in prevalence of infection across the years, assuming a power of 80% and test size of 5%. Oversampling of 20% was considered to take care of any contingencies such as drop-outs, data entry error, non-response, and damaged samples.
Plasmodium falciparum prevalence (PfPR)  was established through rapid diagnostic tests (RDT, Paracheck Ver.3, Orchid Biomedical Systems, Mumbai, India). The selected children were asked to provide a finger-prick blood sample, which was used to assess Plasmodium infection in the peripheral blood. The sample collection was done by trained MOH technicians. All children who were found to be positive were treated with artemisinin-based combination therapy (ACT), artemether-lumefantrine, by the MOH staff according to the national guidelines for the diagnosis, treatment and prevention of malaria in Kenya and Ethiopia [3, 4].
The primary objective of the statistical analysis was to: (1) determine the malaria transmission risk as measured by the relative vector density each year of the survey; (2) determine the prevalence of P. falciparum infection in school children as assessed by the school-level cross-sectional surveys; and, (3) estimate the magnitude of the intervention effects over time for the four intervention groups on malaria transmission risk and P. falciparum infection prevalence.
Entomological analysis for vector characteristics and treatment effects
The following key entomological indicators for measuring vector characteristics were calculated: (1) adult mosquito relative density defined as the number of adult mosquitoes per house per night; and (2) mosquito larval density defined as the number of mosquito larvae per dip. Vector density was estimated using negative binomial regression and adjusted for the number of habitats (for larval mosquitoes) and household clusters (for adult mosquitoes) with the number of trap-nights included as an offset to account for differences in collection effort in the mosquito collection method. The effect of the treatment interventions was estimated using generalized estimating equations (GEE) on vector density, allowing for within-subject correlation using robust variance estimator to calculate standard errors (SEs). From the GEE model, the incidence rate ratios (IRRs) were reported and herein described as density ratios (DR), the control intervention (LLINs only) was used as the reference against the other interventions.
Plasmodium falciparum prevalence and treatment effects
Plasmodium falciparum infection was defined as a positive RDT result [41,42,43]. Proportions were calculated for variables of interest at school level with a 95% confidence interval (CI) using generalized linear models that accounted for school clusters. Overall, analysis of the treatment effect of the various interventions on infection prevalence was initially analysed by fitting GEE model assuming within-subject correlation and binomial family function to test whether the proportion positive for P. falciparum infection varied significantly among the treatment arms. The covariates included in the model were treatment interventions, age, gender, and study site. Additionally, multivariable mixed effects logistic regression model at three levels was separately fitted, individuals nested within schools selected within treatment groups and finally within study site with age and gender retained as fixed terms in the model and reporting the adjusted odds ratio (aOR). The results from the first and second model were compared and found no significant difference between the two models. Hence, the results of the latter model only were reported here.
All statistical analyses were performed using STATA version 14.1 (STATA Corporation, College Station, TX, USA). Graphs were developed using the ggplot package implemented in R .
Entomological findings from adult and larval mosquito collection
Table 1 provides the baseline household and entomological characteristics. The overall baseline mean number of adult anopheline (χ2 = 101.85, p = 0.090) and culicine (χ2 = 146.84, p = 0.178) mosquito densities did not vary significantly in the treatment intervention groups. Similarly, the baseline larval anopheline and culicine mosquito densities did not vary significantly in the treatment intervention groups.
Overall, 24,763 adult anopheline and culicine mosquitoes were collected over the three years in all the three study sites with majority being culicines (83.4%, n = 20,640). The Anopheles captured included An. gambiae (14.1%, n = 3501), An. funestus (1.9%, n = 471), and other Anopheles species (0.6%, n = 151). Over half (69.8%, n = 14,398) of all the culicine mosquitoes collected were from Nyabondo site, followed by Malindi (28.4%, n = 5859) and Tolay (1.9%, n = 383). On the other hand, the majority of all the anopheline mosquitoes collected were from Malindi site (56.6%, n = 2334), followed by Nyabondo (37.5%, n = 1547) and Tolay (5.9%, n = 242).
Anopheles arabiensis was the predominant An. gambiae sibling species in all the sites being more than 99% in Tolay and Nyabondo where the number of specimens identified by PCR were 94 and 100, respectively. The composition of An. gambiae s.l. in Malindi as determined from 207 specimens identified by PCR was An. arabiensis (64%) and Anopheles merus (36%). The sibling species An. gambiae sensu stricto of the An. gambiae complex was virtually absent at the three project sites.
The overall adult anopheline mosquito density was 0.07 with site-specific anopheline density of 0.20, 0.03 and 0.02 for Malindi, Tolay and Nyabondo sites, respectively. The least adult anopheline density was observed in arm 4 (0.04) and arm 3 (0.06) of the study. Notably, arm 4 of the study showed least adult anopheline density in all the three sites (Fig. 3). Similarly, the overall adult culicine density was 0.14 with highest density observed in Malindi study site. Least adult culicines density was observed in arm 4 (0.08) and arm 3 (0.11) of the study. In two study sites, Malindi and Tolay, least adult culicines density was seen in arm 4 of the study while in Nyabondo site it was seen in arm 3.
Analysis of the effect of the treatment interventions on adult vector density showed that overall arm 4 of the study significantly reduced the adult anophelines density by nearly half (DR = 0.55, p = 0.012), and also significantly reduced adult culicines density by nearly two-thirds (DR = 0.38, p < 0.001). In all the three study sites, arm 4 non-significantly reduced both the adult anophelines and culicines density (Table 2).
Mosquito larvae were collected in different habitats around the sampled houses in all the 16 villages in the three study sites. Characteristics of potential aquatic breeding habitats generally fell in the following three categories based on their estimated lifespan: temporary habitats lasting from 2 weeks up to about 3 months including rain puddles, shallow drainage pools especially by the sides of dirt roads, shallow soil excavation sites left behind after brick making activities, and car tracks; semi-permanent habitats lasting up to 6 months including ground pools, seepage pools next to streams, brick-making pits, shallow ponds by the edge of streams; and, permanent habitats lasting more than 6 months including ponds used for domestic water and watering of livestock, fish ponds, wells. The combined semi-permanent and permanent habitats in a village in either Nyabondo or Malindi were about twice as many as those in a village in Tolay.
The overall larval anopheline density was 0.12 with highest anopheline larval density observed in Tolay site (0.66) followed by Malindi (0.55) and Nyabondo (0.08) (Fig. 4). Similarly, the overall larval culicines density was 0.18 with the highest larval culicines density observed in Tolay (1.20), followed by Malindi (0.42) and Nyabondo (0.13) sites. In overall, arm 2 had the least larval vector density for both anophelines and culicines compared to other study arms. At individual site level, arm 2 had the least larval density for anophelines and culicines in Nyabondo and Tolay study sites.
Analysis of the effect of the treatment interventions on larval vector density showed that in overall arm 2 of the study significantly reduced the larval anophelines density by nearly a third (DR = 0.65, p < 0.001), and also significantly reduced the larval culicines density by nearly 12% (DR = 0.88, p = 0.0016). In two study sites, Malindi and Nyabondo, arm 2 significantly reduced the larval anophelines density by 85 and 20% in the two sites, respectively. However, the larval culicines density was significantly reduced in arm 2 and both arms 2 and 3 only in Nyabondo and Malindi study sites, respectively (Table 3).
Epidemiological findings from the school-level surveys
Baseline demographic and epidemiological characteristics for the school children is shown in Table 1. A total of 3599 children (38 schools) with median age of 9 years (range: 3–19 years) were included in the baseline cross-sectional survey in all the three sites. The reported household membership was five people per household with high household membership observed in Malindi study site. Baseline overall malaria prevalence was low 8.8% (95% CI 5.5–14.0%) with the baseline prevalence being significantly different in each of the four study arms, χ2 = 216.06, p < 0.001 and ranged from 5.1 to 12.6%.
Table 4 outlines the number of schools and children surveyed, malaria prevalence and the overall treatment interventions effect on malaria prevalence. The overall malaria prevalence was 13.3% (95% CI 8.8–20.0%) with the highest malaria prevalence observed in Nyabondo site 32.4% followed by Malindi 5.7% and Tolay 1.7%. In overall, highest malaria prevalence was observed in the third year of the survey (the year 2015) compared to the other 2 years.
Overall, the least malaria prevalence was observed in the control arm (7.7%) and highest prevalence in arm 4 of the study (19.3%). However, arm 3 of the study showed least prevalence in two study sites, Malindi and Nyabondo, with arm 4 showing the least prevalence in Tolay site only.
Figure 5 provides the site-specific malaria prevalence by the different treatment interventions. At global level, analysis of the effect of the treatment interventions on malaria prevalence showed that only arm 3 reduced malaria prevalence by 10% (aOR = 0.9, p = 0.356), although this was not significant. The other intervention combinations did not show any prevalence reduction when compared to the control arm of the study. However, at individual site level Tolay showed significant reduction in malaria prevalence in arm 4 by 50% (aOR = 0.5, p = 0.032), while Malindi and Nyabondo showed no prevalence reduction in any of the study arms.
The results demonstrated that children under 5 years of age (aOR = 2.3, p < 0.001) as well as those aged 6 to 10 years (aOR = 1.6, p < 0.001) were at significant risk of malaria infection when compared to older children (over 10 years). However, male children were at non-significant risk of malaria infection compared to females, aOR = 1.1, p = 0.635.
The study assessed the impact of integrating the usage of LLINs with larviciding with Bti, and intensive CEM, on malaria vector populations and malaria prevalence at three rural sites, namely Nyabondo and Malindi in Kenya, and Tolay in Ethiopia. While LLINs and Bti are the actual mosquito control methods, CEM is a social intervention which is necessary for their implementation and sustainability. By analysing the individual and combined effects of two technical and one social interventions, the study aimed at a better understanding of their respective contributions in the set-up broadly defined as IVM strategy [12, 13, 45]. This is the first time in the study areas when a direct comparison of such functionally different technical and social interventions has been attempted, to evaluate their individual and synergistic effects as key elements of IVM.
The relative density of anopheline mosquitoes indoors was found to be generally very low at all the three sites, being less than one mosquito per house per night. The overall average malaria prevalence varied considerably across the three sites with the order of magnitude being 19:3:1 for Nyabondo:Malindi:Tolay, respectively. The assessment of the adult anopheline density at the combined three-site (global level) showed that compared to the control arm (use of LLINs only), LLINs coupled with Bti application and CEM (arm 4) significantly reduced by half the adult anopheline density. Nonetheless, at individual site level, no significant reduction in adult anopheline density was observed in any of the three sites.
For the adult culicines density, the global analysis revealed that all three study arms (i.e., arms 2, 3 and 4) were significantly effective in reducing adult culicines density, with arm 4 showing up to two-thirds reduction. At site-specific level, two sites (i.e., Nyabondo and Tolay) did not show significant reduction of adult culicines density in any of the interventions, however, Malindi site showed that arms 3 and 4 were significantly effective in reducing adult culicines density.
Similarly, assessment of larval anopheline density at global level revealed that arm 2 was significantly effective in reducing larval anopheline density by up to a third. At site-specific level, all the study arms showed significant reduction of larval anopheline density in Malindi by at least 70%, while in Nyabondo only arm 2 reduced density by 20%, and no study arm showed significant reduction of larval anopheline density in Tolay.
Study arm 2 was similarly effective in reducing the overall larval culicines density by about 12%, with the same study arm reducing larval culicines density by nearly a third in Nyabondo site. In Tolay site, both arms 2 and 3 were significantly effective in reducing larval culicines density by at least 30%. However, no significant reductions were observed in any of the interventions in Malindi site.
At global level, arm 3 reduced malaria prevalence by circa 10% though this was not statistically significant. However, it is important to point out that in Tolay site, malaria prevalence was significantly reduced by 50% in arm 4. The two other sites did not reveal significant reduction of prevalence as a result of the use of any of the interventions.
The results of parasite prevalence obtained during this study were similar to what has previously been reported for the three respective project sites [3, 4]. In this regard, the results confirmed that Tolay and Malindi had significantly lower malaria prevalence than Nyabondo as presumed during the selection of the study sites. The prevalence results at the three project sites importantly favoured the original intention of the study, i.e., to have the IVM interventions assessed among situations with low and moderate to high malaria prevalence. In this regard, it was interesting to note that compared to the control (use of LLINs only), the adopted relatively comprehensive strategy in arm 4 significantly reduced malaria prevalence in the setting with the lowest malaria prevalence, i.e., Tolay. The same strategy reduced the relative densities of anopheline and culicine mosquitoes in the settings with higher malaria prevalence. However, in these latter settings neither the relatively comprehensive intervention strategy, nor the strategies separately having either CEM or Bti, had a significant reduction of malaria prevalence when compared to the control arm.
A logical inference based on the above evidence is that bringing on board larviciding with Bti can help to reduce malaria further in settings where prevalence is already low, either naturally or due to conventional implementation of a primary vector control method such as LLINs. As regards the study sites with moderate to high malaria prevalence and more extensive potential anopheline breeding habitats, the results of Bti application corroborated previous assertions of low or no additional impact of larviciding in malaria control, as expressed in a recent systematic review of past studies on the intervention .
According to the results, low adult anopheline densities recorded at all three sites during the 2013–2015 period of the study did not necessarily translate to equally low malaria prevalence across the sites. Thus, malaria prevalence in Nyabondo and Malindi sites in Kenya remained higher than in Tolay, in spite of adult anopheline densities being similarly low at less than one mosquito collected per trap per night indoors. Interestingly, this situation contrasts with that in a follow-up study from 2017 to 2018 in Nyabondo alone, where indoor anopheline densities were higher, and screening of house eaves was found to lower indoor anopheline density as well as malaria parasite prevalence . Further research may therefore be needed to understand how mosquito relative density and other entomological variables of the vectorial capacity  of anopheline populations in an area might vary over time, and the implications of any such variation for malaria prevalence in the human population. Other factors which may require further investigation as part of understanding malaria dynamics in an area include certain human behaviour, such as sleeping, working or socializing outdoors at night, which are associated with outdoor malaria transmission .
In Malindi, the study demonstrated significantly lower anopheline and culicine larval populations in the experimental arms that had Bti, i.e., arm 2 and arm 4, compared to those that did not have Bti. The results support previous findings of Bti’s known high efficacy against mosquito larvae under controlled conditions , However, the failure to also significantly reduce adult anopheline mosquitoes at the same site in the same study arms involving the application of Bti likely point to operational gaps in the way larviciding was conducted. Thus, some of the potential anopheline breeding sites might have gone untreated due to failure to detect them, especially during the wet season when they are many and spread over a wide area. WHO policy recommendation in response to this particular challenge as it relates to malaria vector control is that larviciding should only be considered, with or without other interventions, when the breeding sites in the targeted area are few, fixed and findable . Consequently, WHO recommends that prior to adding more interventions in an actual vector control situation, an assessment should be made to determine if the ineffectiveness of what is currently being applied is due to operational shortcomings that may require resolving through improved management rather than adding of new vector control methods. Operational problems may arise due to inadequate human resources, poor supervision, or lack of tools for mapping and tracking mosquito habitats.
In the case of Tolay study site in Ethiopia, it should be noted that vector control also often involves the use of IRS in addition to LLINs, and this could generally have also contributed to low indoor anopheline vector densities observed in the study site [4, 51]. However, combining LLINs with IRS has been shown not to further reduce malaria incidence beyond what is achieved with LLINs in Ethiopia in spite of its having an impact on indoor vector density . Finally, low indoor anopheline densities in Tolay could also have been due to local houses not having open eaves. On the other hand, low densities in Nyabondo and Malindi could have been due to other environmental, socio-economic or past vector control factors since there was no IRS at the two sites during the study period, and house eaves were generally open.
CEM in the study was expected to influence communities’ behaviour leading to a range of anti-vector activities, including improved usage of LLINs and general environmental management. The lack of significant additional impact on adult anopheline density and malaria prevalence in the experimental arm where CEM was added to LLINs could suggest that, in spite of household net ownership being relatively high and meeting universal coverage, net usage, although not monitored during the study, might have been below the desired optimal levels anticipated in universal coverage . Results on the assessment of how CEM changes people’s knowledge and practice towards vector control will be published separately. However, such an assessment conducted in other research in Tolay area of Ethiopia has recently confirmed that malaria education conveyed through primary school pupils significantly improves community knowledge and leads to behaviour change towards malaria and malaria vector control .
Certain challenges were encountered in the study methodology which require that the results are interpreted with caution, but which may also provide a useful basis for refinement of future similar studies. For instance, larviciding with Bti and CEM were interventions that essentially operated at community level, while the usage of LLINs is a self-protection measure at personal and household level, in spite of having a communal dimension when scaled up in a village setting. Thus, implementing the former two interventions could potentially have led to unintended spillover effects among the study clusters in spite of randomization efforts. For instance, a school enrolled in the cross-sectional malaria surveys and associated with a particular IVM intervention might have had some of its pupils who were tested for malaria resident in clusters with different interventions, but not recorded as such. Furthermore, the CEM intervention may have been confounded by intermingling of residents of different clusters. The problem of spillover effects is well known and has been a subject of discussion in other studies on educational interventions . Moreover, it is also likely that certain information about mosquito and malaria control was generally available at all the study sites due to their being located in areas where malaria research and other anti-malaria activities have been going on for many years, especially in Kenya [19, 55,56,57].
Furthermore, information on universal net coverage was not verified through actual household surveys. Nevertheless, congruence of data from other detailed studies and MIS reports lends credibility to net coverage as cited for the different study sites in “Methods” section of this paper. Lack of more detailed information about larval breeding sites and unavailability of current insecticide resistance test data for the respective study sites were also important limitations of the study.
The study demonstrated that integration of conventional usage of LLINs with larviciding with Bti and CEM led to further reduction in malaria in a setting of low disease prevalence in Tolay in Ethiopia, but not in settings with relatively higher prevalence in Kenya. The results suggest that reducing malaria beyond what has been achieved with LLINs in the latter settings would require a different IVM configuration, possibly including other anti-vector measures as well as integration with prompt diagnosis and effective treatment of malaria cases. The observations underscore the need to view IVM as an adaptive approach which may need to be calibrated differently for different settings based on prior knowledge of the socio-ecological situation, disease epidemiology and effectiveness and practicality of interventions as is often reiterated by WHO .
Availability of data and materials
The data supporting the conclusions of this article are provided within the article. The datasets analyzed are available upon request.
Artemisinin-based combination therapy
Adjusted odds ratio
- Bti :
Bacillus thuringiensis israelensis
- CEM :
Community education and mobilization
Centers for Disease Control
Generalized estimating equations
Information, education and communication
Incidence rate ratios
Indoor residual spraying
Integrated vector management
Long-lasting insecticide treated nets
Ministry of Health
National Malaria Control Programmes
Polymerase chain reaction
Plasmodium falciparum prevalence
Rapid diagnostic test
World Health Assembly
World Health Organization
WHO/UNICEF. Achieving the Malaria MDG Target: reversing the incidence of malaria 2000–2015. Geneva: World Health Organization; 2015.
WHO. World Malaria Report 2018. Geneva: World Health Organization; 2018.
Kenya MOH. Kenya Malaria Indicator Survey 2015. Nairobi: National Malaria Control Programme, Ministry of Health; 2016.
Ethiopia MOH. Ethiopia National Malaria Indicator Survey 2015. Addis Ababa: Ethiopian Public Health Institute, Ministry of Health; 2016.
Snow RW, Kibuchi E, Karuri SW, Sang G, Gitonga CW, Mwandawiro C, et al. Changing malaria prevalence on the Kenyan Coast since 1974: climate, drugs and vector control. PLoS ONE. 2015;10:e0128792.
Tusting L, Willey B, Lucas H, Thompson J, Kafy HT, Smith R, et al. Socioeconomic development as an intervention against malaria: a systematic review and meta-analysis. Lancet. 2013;382:963–72.
Meyrowitsch DW, Pedersen EM, Alifrangis M, Scheike TH, Malecela MN, Magesa SM, et al. Is the current decline in malaria burden in sub-Saharan Africa due to a decrease in vector population? Malar J. 2011;10:188.
Mwangangi JM, Mbogo CM, Orindi BO, Muturi EJ, Midega JT, Nzovu J, et al. Shifts in malaria vector species composition and transmission dynamics along the Kenyan coast over the past 20 years. Malar J. 2013;12:13.
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.
Hemingway J. The role of vector control in stopping the transmission of malaria: threats and opportunities. Phil Trans R Soc B. 2014;369:20130431.
Leach-Kemon K, Chou DP, Schneider MT, Tardif A, Dielman JL, Brooks BPC, et al. The global financial crisis has led to a slowdown in growth of funding to improve health in many developing countries. Health Aff. 2012;31:228–35.
WHO. Global Strategic Framework for Integrated Vector Management. Geneva: World Health Organization; 2004.
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.
WHO. Global Vector Control Response 2017–2030. Geneva: World Health Organization; 2017.
Kleinschmidt I, Schwabe C, Shiva M, Segura JL, Sima V, Mabunda SJ, et al. Combining indoor residual spraying and insecticide-treated net interventions. Am J Trop Med Hyg. 2009;81:519–24.
Okumu FO, Moore SJ. Combining indoor residual spraying and insecticide-treated nets for malaria control in Africa: a review of possible outcomes and an outline of suggestions for the future. Malar J. 2011;10:208.
Kirby MJ, Ameh D, Green C, Jawara M, Milligan PJ, Bottomley C, Snell PC, Conway DJ, Lindsay SW. Efficacy of two different house screening interventions against exposure to malaria and anaemia in children in The Gambia: a randomized controlled trial. Lancet. 2009;374:998–1009.
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.
Mutero C, Mbogo C, Mwangangi J, Imbahale S, Kibe L, Orindi B, et al. An assessment of participatory integrated vector management for malaria control in Kenya. Environ Health Perspect. 2015;123:1145–51.
Imbahale S, Mukabana WR. Efficacy of neem chippings for mosquito larval control under field conditions. BMC Ecol. 2015;25:8.
Demissew A, Balkew M, Girma M. Larvicidal activity of chinaberry, neem and Bti to an insecticide resistant population of Anopheles arabiensis from Tolay, South West Ethiopia. Asian Pac J Trop Biomed. 2016;6:554–61.
Macintyre K, Keating J, Sosler S, Kibe L, Mbogo CM, Githeko AK, et al. Examining the determinants of mosquito-avoidance practices in two Kenyan cities. Malar J. 2002;1:14.
Kibe LW, Mbogo CM, Keating J, Molyneux S, Githure JI, Beier JC. Community based vector control in Malindi, Kenya. Afr Health Sci. 2006;6:240–6.
Mwangangi JM, Kahindi SC, Kibe LW, Nzovu JG, Luethy P, Githure JI, Mbogo CM. 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.
ICIPE 2012. Report of the comprehensive evaluation of the Biovision-ICIPE IVM projects in Kenya and Ethiopia. Nairobi. International Centre of Insect Physiology and Ecology; 2012.
Imbahale SS, Abonyo OK, Aduogo OP, Githure IJ, Mukabana WR. Conflict between the need for income and the necessity of controlling endemic malaria. J Ecosys Ecograph. 2013;3:129.
Ng’ang’a PN, Okoyo C, Mbogo C, Mutero CM. Evaluating effectiveness of screening house eaves as an intervention for integrated vector management for malaria control in a rural area of Western Kenya. Malar J. 2020;19:341.
WHO. Methods for achieving universal coverage with long-lasting insecticidal nets in malaria control. Report to MPAC. Geneva. Vector Control Technical Expert Group. World Health Organization; 2013.
Zhou G, Li JS, Ototo EN, Atieli HE, Githeko AK, Yan G. Evaluation of universal coverage of insecticide-treated nets in western Kenya: field surveys. Malar J. 2014;13:351.
Thomas CJ, Cross DE, Bøgh C. Landscape movements of Anopheles gambiae malaria vector mosquitoes in rural Gambia. PLoS ONE. 2013;8:e68679.
Valent Biosciences. Technical Use Sheet. VectorBac Granule. Biological larvicide. Accessed online June 2019 at: https://www.valentbiosciences.com/publichealth/wp-content/uploads/sites/4/2019/04/VectoBac-Granule-Technical-Use-Sheet-sp.pdf.
WHO. How to design vector control efficacy trials, guidance on phase III vector control field trial design. Geneva: World Health Organization; 2017.
Kramer RA, Mboera LEG, Senkoro K, Lesser A, Shayo EH, Paul CJ, et al. A randomized longitudinal factorial design to assess malaria vector control and disease management interventions in rural Tanzania. Int J Environ Res Public Health. 2014;11:5317–32.
Benjamin-Chung J, Arnold BF, Berger D, Luby SP, Miguel E, Colfor JM Jr, et al. Spillover effects in epidemiology: parameters, study designs and methodological considerations. Int J Epidemiol. 2018;47:332–47.
Service MW. Mosquito ecology. Field sampling methods. 2nd ed. London: Elsevier Chapman and Hall; 1993.
Gillies MT, De Meillon B. The Anopheline of Africa South of the Sahara (Ethiopian Zoogeographical region). Publications of the South African Institute for Medical Research; 1968.
Mboera LEG, Kihonda J, Braks MA, Knols BGJ. Influence of Centers for Disease Control light trap position, relative to a human-baited bednet, on catches of Anopheles gambiae and Culex quinquefasciatus in Tanzania. Am J Trop Med Hyg. 1998;59:595–6.
Paskewitz SM, Collins FH. Use of the polymerase chain reaction to identify mosquito species of the Anopheles gambiae complex. Med Vet Entomol. 1990;4:367–73.
Scott JA, Brodgon WG, Collins FH. Identification of single specimens of Anopheles gambiae complex by 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 (PCR) assay to identify the Anopheles funestus (Diptera: Culicidae) group. Am J Trop Med Hyg. 2002;6:78–83.
Drakeley C, Abdulla S, Agnandji ST, Fernandes JF, Kremsner P, et al. Longitudinal estimation of Plasmodium falciparum prevalence in relation to malaria prevention measures in six sub-Saharan African countries. Malar J. 2017;16:433.
Gitonga CW, Karanja PN, Kihara J, Mwanje M, Juma E, Snow RW, et al. Implementing school malaria surveys in Kenya: towards a national surveillance system. Malar J. 2010;9:306.
Gitonga CW, Kihara JH, Njenga SM, Awuondo K, Noor AM, Snow RW, et al. Use of rapid diagnostic tests in malaria school surveys in Kenya: does their under-performance matter for planning malaria control? Am J Trop Med Hyg. 2012;87:1004–11.
Wickham H. Ggplot2: elegant graphics for data analysis. Berlin: Springer; 2009.
Mutero CM, Schlodder D, Kabatereine N, Kramer R. Integrated vector management for malaria control in Uganda: knowledge, perceptions and policy development. Malar J. 2012;11:21.
Choi L, Majambere S, Wilson AL. Larviciding to prevent malaria transmission. Cochrane Database Syst Rev. 2019;8:CD012736.
Brady OJ, Godfray HCJ, Tatem AJ, Gething PW. Vectorial capacity and vector control: reconsidering sensitivity to parameters for malaria elimination. Trans R Soc Trop Med Hyg. 2016;110:107–17.
Monroe A, Moore S, Koenker H, Lynch M, Ricotta E. Measuring and characterizing night time human behaviour as it relates to residual malaria transmission in sub-Saharan Africa: a review of the published literature. Malar J. 2019;18:6.
Dambachet P, Louis VR, Kaiser A, Quedraogo S, Sié A, Sauerbron R, et al. Efficacy of Bacillus thuringiensis var. israelensis against malaria mosquitoes in north-western Burkina Faso. Parasit Vectors. 2014;7:371.
WHO. Larval source management: a supplementary measure for malaria vector control. An operational manual. Geneva: World Health Organization; 2013.
Kenea O, Balkew M, Tekie H, Deressa W, Loha E, Lindtjørn B, et al. Impact of combining indoor residual spraying and long-lasting insecticidal nets on Anopheles arabiensis in Ethiopia: results from a cluster randomized controlled trial. Malar J. 2019;18:182.
Loha E, Deressa W, Gari T, Balkew M, Kenea O, Solomon T, et al. Long-lasting insecticidal nets and indoor residual spraying may not be sufficient to eliminate malaria in a low malaria incidence area: results from a cluster randomized controlled trial in Ethiopia. Malar J. 2019;18:141.
Kebede Y, Abebe L, Alemayehu G, Sudhakar M, Birhanu Z. Effectiveness of peer-learning assisted primary school students educating the rural community on insecticide-treated nets utilization in Jimma-zone Ethiopia. Malar J. 2020;19:331.
VanderWeele TJ, Hong G, Jones SM, Brown JL. Mediation and spillover effects in group-randomized trials: a case study of the 4Rs educational intervention. J Am Stat Assoc. 2013;108:469–82.
USAID. Malaria in Kenya. Washington: USAID Knowledge Services Centre; 2008.
Jenkins R, Omollo R, Ongecha M, Sifuna P, Othieno C, Ongeri L, et al. Prevalence of malaria parasites in adults and its determinants in malaria endemic area of Kisumu County, Kenya. Malar J. 2015;4:263.
Idris ZM, Chan CW, Kongere J, Gitaka J, Logedi J, Omar A, et al. High and heterogeneous prevalence of asymptomatic and sub-microscopic malaria infections on islands in Lake Victoria, Kenya. Sci Rep. 2016;6:36958.
WHO. A tool kit for integrated vector management in sub-Saharan Africa. Geneva: World Health Organization; 2016.
The authors wish to thank icipe and KEMRI field staff, members of the community and schools in Malindi and Nyabondo in Kenya, and Tolay in Ethiopia, for their general support and participation in the study.
The work was conducted with funding from Biovision Foundation (Grant No. BV HH-07/2013-2015), UK’s Foreign, Commonwealth & Development Office (FCDO); the Swedish International Development Cooperation Agency (Sida); the Swiss Agency for Development and Cooperation (SDC); the Federal Democratic Republic of Ethiopia; and the Government of the Republic of Kenya. The views expressed herein do not necessarily reflect the official opinion of the donors.
Ethics approval and consent to participate
The study protocol received ethical approval from the Kenya Medical Research Institute (KEMRI) Scientific and Ethics Review Unit (SERU) (SSC Protocol No. 2675), as part of Integrated Vector Management (IVM) for Sustainable Malaria Control in Eastern Africa (Ref: KEMRI/RES/7/3/1). Additional approval was provided by the appropriate ministries of health in Kenya and Ethiopia. At household level, approval was provided by the house owner. At school, parental consent was based on passive opt-out consent rather than written opt-in consent owing to the low risk and routine nature of the survey procedure. Child individual assent was obtained before participation in the surveys.
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.
About this article
Cite this article
Mutero, C.M., Okoyo, C., Girma, M. et al. Evaluating the impact of larviciding with Bti and community education and mobilization as supplementary integrated vector management interventions for malaria control in Kenya and Ethiopia. Malar J 19, 390 (2020). https://doi.org/10.1186/s12936-020-03464-6