A crossover study to evaluate the diversion of malaria vectors in a community with incomplete coverage of spatial repellents in the Kilombero Valley, Tanzania
Parasites & Vectors volume 9, Article number: 451 (2016)
Malaria elimination is unlikely to occur if vector control efforts focus entirely on transmission occurring indoors without addressing vectors that bite outdoors and outside sleeping hours. Additional control tools such as spatial repellents may provide the personal protection required to fill this gap. However, since repellents do not kill mosquitoes it is unclear if vectors will be diverted from households that use spatial repellents to those that do not.
A crossover study was performed over 24 weeks in Kilombero, Tanzania. The density of resting and blood-engorged mosquitoes and human blood index (HBI) of malaria vector species per household was measured among 90 households using or not using 0.03 % transfluthrin coils burned outdoors under three coverage scenarios: (i) no coverage (blank coils); (ii) complete coverage of repellent coils; and (iii) incomplete coverage of repellent and blank coils. Mosquitoes were collected three days a week for 24 weeks from the inside and outside of all participating households using mosquito aspirators. Paired indoor and outdoor human landing collections were performed in three random households for six consecutive nights to confirm repellent efficacy of the coils and local vector biting times.
The main vectors were Anopheles arabiensis and Anopheles funestus (sensu stricto), which fed outdoors, outside sleeping hours, on humans as well as animals. Anopheles arabiensis landings were reduced by 80 % by the spatial repellent although household densities were not reduced. The HBI for An. arabiensis was significantly higher among households without repellents in the incomplete coverage scenario compared to houses in the no coverage scenario (Odds ratio 1.71; 95 % CI: 1.04–2.83; P = 0.03). This indicated that An. arabiensis mosquitoes seeking a human blood meal were diverted from repellent users to non-users. The repellent coils did not affect An. funestus densities or HBI.
Substantial malaria vector activity is occurring outside sleeping hours in the Kilombero valley. Repellent coils provided some protection against local An. arabiensis but did not protect against local (and potentially pyrethroid-resistant) An. funestus. Pyrethroid-based spatial repellents may offer a degree of personal protection, however the overall public health benefit is doubtful and potentially iniquitous as their use may divert malaria vectors to those who do not use them.
Insecticide-treated nets (ITNs) protect people from mosquito bites and kill host seeking mosquitoes during sleeping hours and consequently reduce malaria transmission . However, the protective efficacy of ITNs is dependent on a series of factors: the time and place where malaria vectors bite , the users’ sleeping hours, correct maintenance of the net, and consistent use of the bed-net . Because the nets prevent mosquitoes from obtaining a blood meal from their human hosts during the late hours of the night, these mosquitoes have been shown to adapt their behavior to host-seeking periods earlier in the evening when people are still active [4–7]. In this case, malaria transmission occurs when people are not yet sleeping under a protective bed net. Rural communities in Tanzania commonly sit outside their homes during and after dusk. In the Kilombero Valley, the majority of the population retires to bed between 9 and 10 pm, until which an estimated 20 % of the total malaria transmission occurs due to early evening biting . Repellents may be a method suitable to provide protection when people are sitting outdoors in the evening.
Repellents are compounds that interfere with the mosquito’s olfactory system preventing them from identifying their hosts and succeeding in taking a blood-meal . These compounds may be applied directly to the skin as a topical formulation; or dispersed into the air creating a mosquito bite free “space” that protects all people within that space. In the latter case they are called spatial repellents. Examples of spatial repellents include mosquito coils, emanators, as well as traditional practices such as burning or smoldering repellent plants . The most common and best-characterized spatial repellent is the mosquito coil. These have been used for centuries in Asia and are traditionally made from a dried paste of pyrethrum powder shaped into a spiral. Currently mosquito coils can be found throughout the world containing volatile pyrethroids such as D-alethrin and transfluthrin. These compounds act as insect neurotoxins and are dispersed into the surrounding air by the smoke of the burning coil. When the mosquito is exposed to the volatilized active ingredient it becomes unable to find a host and feed .
Repellents reduce human vector contact but do not kill mosquitoes. Therefore there is a possibility of mosquitoes being redirected (diverted) from repellent users to those individuals who do not use repellents. A previous study conducted in Tanzania investigated the degree of mosquito diversion occurring in villages with incomplete coverage of topical repellents between households that used the topical repellent 15 %-DEET (N,N-Diethyl-meta-toluamide) to households that did not use a repellent . Households not assigned to the repellent group and surrounded by households that did use repellents had three-fold more mosquitoes in and around their dwellings opposed to households in an area where repellents were not used at all. Thus, if complete coverage is not reached topical repellents will divert mosquitoes to households of individuals who do not use them and could potentially increase their exposure to mosquito-borne diseases such as malaria and filariasis. This would lead to health inequity between community members but also would allow pathogens to persist within the community through the vulnerable pockets of population that do not use topical repellents. On the other hand, spatial repellent pyrethroids such as transfluthrin and metofluthrin do not necessarily force the mosquito to move away from the repellent source. Unlike DEET, their main mode of action is neural excitation that results inhibition of host seeking and feeding rather than to repel . However, the host-seeking mosquitoes remain in the environment and may recover from their exposure to the volatile pyrethroid and become ready to host-seek once again. Also, exposure to volatile chemicals outdoors may be insufficient to disarm host-seeking mosquitoes if they simply move away from the repellent source before picking up a sufficient dose. The current study was designed to investigate if diversion of malaria vectors occurs between households that use mosquito coils to those that do not. Health equity implications and the role of spatial repellents as a tool to reduce outdoor residual malaria transmission in the Kilombero Valley are discussed.
The project was conducted for 24 weeks from December 2012 to June 2013 in Mbingu, which is approximately 40 km west of Ifakara at 8.21°S and 36.24°E, in Kilombero Valley, in south-eastern Tanzania. The site is characterized by typical rural houses surrounded by rice, maize and banana fields close to the Londo River and the slopes of the Udzungwa Mountains. The main mosquito disease vectors are Anopheles arabiensis, An. funestus (sensu stricto) (s.s.), Mansonia africanus, Mansonia uniformis, Coquillettidia aureus and Culex univittatus .
Three villages within the area of Mbingu were selected as study site. A total of 90 households, 30 from each of three villages, Matete, Uwata and Igima, were enrolled in the project. These villages were selected because they had 30 households enabling the crossover design. Geographically, the villages were nearly aligned in a straight line with Igima in the south, Uwata in the north and Matete in between. The distance between Matete and Igima was around 2.5 km and between Matete and Uwata 1 km. Households in the villages were no further than 100 m apart from each other whilst never being directly next to each other (<10 m). All households from these villages agreed to participate in the study; a total 90 households and 381 individuals were enrolled. Each recruited household was given one Olyset® long-lasting insecticide-treated net per sleeping space before the start of the study. Households were numbered and their structural characteristics recorded: wall structure, presence or absence of eaves, roof type, screened or unscreened windows, number of rooms as well as number of occupants and ownership of domestic animals.
The study was designed to test if mosquitoes are diverted from households that use treated mosquito coils to those that do not. A crossover design was used because of the heterogeneous nature of the data generated from resting mosquito collections from households. This design accounted for differences between households by allowing each household to act as its own control. Three coverage scenarios were rotated between the three villages on a two weekly basis for a period of 24 weeks from December 2012 to May 2013 (Fig. 1). The objective of the study was to measure mosquito densities inside and outside households under three coverage scenarios: (i) Complete coverage; (ii) No coverage; and (iii) Incomplete coverage with 80 % coverage of repellents where 24 households were given treated coils and 6 households were randomly selected using a lottery system and assigned to the repellent non-user group. These households were then removed from the lottery in the following weeks and new houses were randomly selected in order to allow all households to be assigned to the repellent non-user group in the incomplete coverage scenario. The coverage level of 80 % was selected based on a previous study conducted in the same area investigating mosquito diversion caused by topical repellents at 80 % coverage level .
Mosquito coils were chosen as the spatial repellent intervention because they are simple to use and well characterized . During the first six weeks households assigned to the no coverage scenario as well as repellent non-users of the incomplete coverage scenario were not given any type of treatment. After week 7 and until the end of the 24 weeks completely identical blank coils were introduced as a placebo. This was done to assess the effect of smoke produced by the blank coils on mosquito densities. Both blank coils and treated coils were produced for this study by SC Johnson and had identical packages, shape, color and lavender smell. Participants were not aware if they were assigned a treated or an untreated coil. Given the nature of a crossover design researchers and the head of the field technicians were not blinded to the intervention. In order to obtain a protective “bubble” effect, household leaders were asked to burn three coils each evening around the space where they sat outdoors in the evenings . Coils were used every day of the week, they were lit at sunset and extinguished once householders retired to bed. Compliance was assessed each morning by inspecting the ashes produced by the mosquito coils during the previous night. Mosquito coils were distributed on a daily basis by the field team and were only sufficient for one night; this was done in order to avoid sharing between households. In addition, the householders were asked what time they had gone to bed the previous evening.
Resting mosquitoes were collected every week on three consecutive weekdays for 24 weeks. Collections were done in all 90 households from all three villages on each experimental day. The team collecting the mosquitoes rotated each between villages to reduce collection bias. Coverage scenarios were assigned to each village on a 2-weekly basis. A total of 72 village/week blocks (3 villages over 24 weeks) were divided between coverage scenarios: (i) Complete coverage: 18 weeks (54 days of resting collections); (ii) No coverage scenario: 18 weeks (54 days of resting collections); and (iii) Incomplete coverage scenario: 36 weeks (108 days of resting collections).
Resting mosquitoes were collected early morning using Prokopack aspirating devices from under outdoor kitchen thatch roofs as well as indoors . All households were surveyed within a 30-minute window period to avoid differences in resting mosquito densities caused by mosquitoes naturally leaving the households in the mornings. The outdoor kitchens, also known as “kibanda” were selected as previous evaluations had shown them to be a common outdoor resting site for mosquitoes . These structures are a basic thatch roof supported by wooden poles.
Mosquitoes were morphologically identified to genus level and the anophelines were identified to species complex level using morphological keys . A sample of Anopheles gambiae (sensu lato) (s.l.) and Anopheles funestus (s.l.) was pooled and sent to the IHI laboratory for PCR species identification [17, 18]. A subsample of the PCR confirmed specimens was sent to a laboratory at the University of Notre Dame for quality control. The blood meals of engorged mosquitoes were analyzed using ELISA (enzyme linked immunosorbent assay) to detect the source of blood in engorged mosquitoes: human or animal and allow the calculation of human blood index (HBI).
Human landing catches
All night indoor-outdoor human landing catches were performed on 6 consecutive nights from dusk till dawn (18:30 to 06:30) during the month of May 2013. The methodology used aimed at reducing all possible risks to collectors and supervisors as described by Gimnig et al. . A total of 6 volunteers were grouped in pairs and each pair was assigned to a randomly selected village household. Each day 3 new households from the study group were randomly selected using a lottery system. A total of 18-paired human landing collections were performed of which 14 were done in households using blank coils and 4 were done in houses using treated coils. This was designed to evaluate the effect of the treatments as well as the biting patterns of the vectors without the spatial repellent. One individual was asked to collect mosquitoes indoors and the other outdoors. After each hour the two individuals swapped positions. All volunteers were males between 18 and 45 years old and recruited under written informed consent. Before starting the mosquito collection each volunteer was seen by a clinical officer who assessed the individual’s general fitness and their malaria status by malaria rapid diagnostic test (mRDT) SD Bioline Pf Pan. In order to protect the volunteers from malaria exposure, mefloquine prophylaxis was provided to each individual. In addition, volunteers were asked to wear closed shoes and short trousers and were provided with “bug-jackets” made of untreated netting that completely covered their upper body. Mosquitoes landing on the volunteers bare lower legs were siphoned into a paper cup labeled with time and location. The volunteers were given a hot meal before beginning collections and were allowed 15 minutes break each hour to drink water, tea or use the restroom. The chief field technician supervised the collections and volunteer’s alertness was checked on an hourly basis. The following day a trained technician morphologically identified the mosquitos’ genus or species complex. A subsample of An. gambiae (s.l.) and An. funestus (s.l.) was pooled and sent to the IHI laboratory for PCR species identification.
The analyses were carried out in R 2.10.0  using the package lmer and glmm ADMB  and in STATA (version 11), (StataCorp) . The densities of resting and fed mosquitoes were compared between the three coverage scenarios using Generalized Linear Mixed Models , in which differences between coverage scenarios were assessed whilst adjusting for village, household and coverage scenario. Mosquito abundance data are typically over-dispersed [24, 25], consequently a negative binomial distribution was used to model the data after assessing their fit to this distribution. Model fit was assessed using a likelihood ratio test and visually by examining the residuals versus fitted plot for deviation from the assumptions of linearity and homoscedasticity. Mosquito densities were compared for the following: (i) household in a no coverage scenario (reference); (ii) household in a complete coverage scenario; (iii) household using repellent coils in an incomplete coverage scenario; and (iv) households not using repellent coils in an incomplete coverage scenario. Coverage scenario was included as a fixed effect together with week to adjust for overall conditions at the time of collection, while village and household were included as random effects. Households assigned to treated coils but that did not comply were re-assigned to “incomplete coverage scenario - repellent non-users” as they did not use a repellent intervention. First, the impact of coverage scenario on resting mosquito and fed mosquito abundance was estimated from the models. We estimated the incidence rate ratios (IRR) and 95 % confidence intervals (95 % CI) of mosquito density per household per night by comparing each coverage scenarios to the coverage scenario where everyone was given a blank coil. Human blood index (HBI) was calculated for each coverage scenario by dividing the total number of mosquitoes fed on humans by the total number of blood-fed mosquitoes. The data for human blood fed as a proportion of blood-fed mosquitoes was analysed with a generalised mixed effects model using a binary distribution with coverage scenario and week included as a fixed effects and village and household as random effects and the odds ratios (OR) and 95 % confidence intervals (95 % CI) were estimated.
Data from the human landing collections was analysed to determine differences between landing rates of mosquitoes when exposed to treated coils opposed to blank coils. A mixed effects Poisson regression with a random intercept to account for overdispersion was carried out in STATA version 11 (Stata Corp) . The presence or absence of the coil was a fixed effect and household and day were random effects. We estimated the IRR of mosquitoes landing on a human per hour when sitting surrounded by three blank coils relative to three transfluthrin 0.03 % coils.
Study villages and household participants
Most households of Matete, Uwata and Igima had brick walls with roof made of thatch or metal sheets; eaves were present in virtually all households. The median number of household members was four individuals (Inter-quartile range IQR = 2–6). The median number of rooms was three (IQR = 2–4). Almost every household owned chickens or ducks as domestic animals; goats were common but cattle were found only in one household. Villages were statistically different in terms of roof, wall type and number of rooms per household (Table 1) underlining the importance of the crossover design used. All the villagers cooked outside and ate their meals outdoors before going to bed; the average bedtime was around 21:00.
Malaria vector species composition and biting patterns
Biting patterns were plotted using human landing catch data collected from households using blank coils only. A median number of 19 mosquitoes per household were collected indoors (IQR = 6–40) and a median number of 16 mosquitoes per household were collected outdoors (IQR = 9–27). Sixty-two percent of the collected mosquitoes were Culex spp., 24 % An. arabiensis, 8 % An. funestus and 5 % Mansonia spp.
Results show that An. arabiensis were mostly host-seeking during the early evening hours with peak biting time between 20.00 and 22.00 h as well as activity in the early morning (Fig. 2). Anopheles funestus preferred biting indoors. Early morning biting was observed for both species. Approximately 23 % of exposure to An. arabiensis and 11 % of exposure to An. funestus bites occurred outdoors outside sleeping hours.
A total of 1,625 (7.5 % of the total catch) Anopheles gambiae (s.l.); 1,757 (8.1 % of the total catch) Anopheles funestus (s.l.); 17,793 (83 % of the total catch) Culex spp. (mainly Cx. quinquefasciatus and Cx. univittatus); 323 (1.4 % of the total catch) Mansonia africana/ M. uniformis; 13 (0.006 % of the total catch) Aedes spp.; and 104 (0.04 % of the total catch) Coquillettidia aureus/ C. versicolor were collected from resting catches (Table 2). Human landing catches collected 1,113 (24 %) An. gambiae (s.l.); 369 (8 %) An. funestus (s.l.); 2,860 (62.5 %) Culex spp.; 211 (4.6 %) Mansonia spp.; and 17 (0.4 %) Coquilettidia spp. PCR was performed on 1,412 An. gambiae (s.l.) of which 724 (100 % successful amplifications) amplified as An. An arabiensis; 1,393 Anopheles funestus (s.l.) of which 345 (91 %) amplified as An. funestus (s.s.); 21 (6 %) as An. leesoni; and 11 (3 %) as An. rivolorum. Because the amplification rate was low, further PCR amplifications were done for quality assurance at the University of Notre Dame: 163 An. gambiae (s.l.) were analysed, all amplified, of which 160 (98 %) were An arabiensis, 3 (1.5 %) An. quadriannulatus and one (0.5 %) An. gambiae (s.s.); 58 An. funestus (s.l.) were analysed, 57 (98 %) amplified as An. funestus (s.s.) and one (2 %) as An. maculipennis.
Assessment of treatment efficacy
The 0.03 % transfluthrin-treated coils repelled approximately 80 % of the An. arabiensis mosquitoes attempting to bite a human outdoors (IRR = 0.22, 95 % CI: 0.052–0.91, P = 0.037) (Table 2). In contrast, both An. funestus and Culex spp. were not significantly repelled by the treated coils (An. funestus: IRR = 0.69, 95 % CI: 0.21–2.27, P = 0.54; Culex spp.: IRR = 1.27, 95 % CI: 0.54–2.97, P = 0.59).
The effect of smoke from blank coils did not significantly reduce the number of mosquitoes resting in the households (An. arabiensis: IRR = 0.84, 95 % CI: 0.58–1.20, P = 0.33; An. funestus: IRR = 1.10, 95 % CI: 0.71–1.73, P = 0.66; Culex spp.: IRR = 0.90, 95 % CI: 0.68– 1.19, P = 0.45), but it did elicit a significant reduction in fed mosquitoes of all types (An. arabiensis, IRR = 0.53, 95 % CI: 0.33–0.84, P = 0.007; An. funestus, IRR = 0.56, 95 % CI: 0.34–0.95, P = 0.031; Culex spp., IRR = 0.16, 95 % CI: 0.10–0.24, P < 0.0001) (Table 3).
Resting collections and diversion
Diversion of An. arabiensis from households using repellent coils to non-users was seen in the incomplete coverage scenario. Human blood index (HBI), for An. arabiensis was significantly higher among those households without repellents in the incomplete coverage scenario with HBI increasing from 0.20 to 0.30 (Odds ratio = 1.71, 95 % CI: 1.04–2.83, P = 0.03) (Table 4). However, the use of coils had no impact on the probability of finding a resting or a blood-fed An. arabiensis in the household (resting: IRR = 0.87, 95 % CI: 0.70–1.09, P = 0.21; blood-fed: IRR = 0.91, 95 % CI: 0.68–1.21, P = 0.59) (Tables 5 and 6). Nightly An. arabiensis densities were highest in households that did not use repellent coils regardless of coverage scenario (households in no coverage scenario IRR = 1; repellent coil non-users in the incomplete coverage scenario IRR = 0.91, P = 0.46). The same was observed for blood-engorged females, with highest densities found in households allocated to the no coverage scenario (IRR = 0.86, 95 % CI: 0.62– 1.19, P = 0.32) (Table 6).
No diversion was seen between coil users and non-users with An. funestus in the incomplete coverage scenario. The HBI was similar in all coverage scenarios ranging between 0.66 and 0.71 and the use of mosquito coils had no effect on the proportion of An. funestus fed on humans (Table 4). Unexpectedly, the use of mosquito coils increased the probability of finding An. funestus in the household (repellent users in complete coverage scenario: IRR = 1.44, 95 % CI: 1.16–1.79, P < 0.001; repellent users in the incomplete coverage scenario: IRR = 1.63, 95 % CI: 1.31–2.02, P < 0.001; both compared to no coverage.). Those houses that did not use spatial repellents in the incomplete coverage scenario (where one would expect diversion) also had more An. funestus in and around their dwellings compared to when nobody was given a repellent coil (IRR = 1.56, 95 % CI: 1.22–1.98, P = 0.003) but this was not different from the densities measured among repellent users and so not indicative of diversion. The number of blood-fed mosquitoes was also 40 % higher among repellent users in both the complete coverage (IRR = 1.35, 95 % CI: 1.01–1.80, P = 0.004) and incomplete coverage scenario (IRR = 1.39, 95 % CI: 1.04–1.86, P = 0.02). Surprisingly, 30 % of the An. funestus (s.s.) population had fed on animals (Table 4) indicating that high LLIN use is forcing An. funestus to feed on alternate hosts to humans. In addition, very little An. funestus activity occurred in the early evening (Fig. 2) so the use of spatial repellents would have a limited impact in preventing man-vector contact with this vector.
No diversion was measured with Culex species through change in the HBI, with HBI ranging between 0.30 for no intervention and 0.50 for houses with no repellent in the incomplete coverage scenario, and this was not statistically significant (OR = 1.58, 95 % CI: 0.96–2.62, P = 0.07) (Table 4). Households not using repellents in the incomplete coverage scenario had 19 % more Culex spp. resting in their dwellings (IRR = 1.19, 95 % CI: 1.06–1.34, P = 0.003) but the number of blood-fed Culex spp. was not significantly different between coverage scenarios. Coils reduced household densities of Culex spp. (repellent users in the complete coverage scenario: IRR = 0.74, 95 % CI: 0.66–1.77, P < 0.001; repellent users in the incomplete coverage scenario: IRR = 0.70, 95 % CI: 0.63–0.77, P < 0.001).
There was no withdrawal of study participants throughout the project. Total compliance was measured at a 92 % level and was similar in all coverage scenario groups (no coverage with blank coils: 93 % compliance; complete coverage with treated coils: 91 % compliance; incomplete coverage blank coil users: 90 % compliance; and incomplete coverage treated coil users: 90 % compliance). Reasons given for not complying with lighting the coils included: (i) not being at home that evening; (ii) family had gone to bed earlier than sunset; (iii) very windy weather; or (iv) because some individuals felt unwell for reasons unrelated to the project and the smoke disturbed them, such as “mafua” (Kiswahili word to describe congestion of the nose). Generally the community was pleased to use the mosquito coils and no adverse effects were reported.
Transfluthin mosquito coils reduced mosquito bites from An. arabiensis by 80 % among people outdoors although they did not affect densities of An. arabiensis. However, households that did not use spatial repellents but were surrounded by households using treated coils had a greater probability of mosquitoes containing a blood meal from a human (Odds ratio 1.71), indicating a slight diversion of bites of Anopheles arabiensis to repellent non-users. These findings agree with those of Maia et al.  conducted in the same field site where topical repellents diverted mosquitoes to repellent non-using households and increased mosquito densities by over four-fold, although the HBI was not measured in that study and so the results are not directly comparable. The data may also explain to some extent the findings of Wilson et al.  demonstrating that mosquito repellents do not reduce malaria at a community scale and that this effect is particularly apparent among communities where compliance with personal protection is lower. Complete coverage of repellent usage at the community level can never realistically be reached and this may lead to diversion of mosquitoes to those who are not protected. In terms of public health benefit it is possible that repellents in general could lead to inequity by increasing the exposure of less informed or poorer individuals to bites of malaria vectors.
The transfluthrin coils did not repel An. funestus. This was consistently observed in the data from the human landing collections as well as from the resting collections. Diversion cannot occur if the treatment does not produce a repellent effect that would cause mosquitoes to preferentially feed on repellent non-users. Although the number of An. funestus found in households that did not use repellents but were surrounded by repellent users was high, this should not be recognised as diversion. These results indicate that the volatile pyrethroids did not affect the host-seeking behaviour of wild An. funestus and are not providing protection against this malaria vector. It is a shortfall of the study that pyrethroid resistance, including susceptibility to transfluthrin was not established. For future studies it is recommended to measure mosquito susceptibility to the repellents used when investigating the effect of volatile pyrethroids on mosquito behaviours, as it has been observed that pyrethroid resistance is a crucial factor that influences the behavioural susceptibility of Aedes aegypti mosquitoes to transfluthrin .
Households using transfluthrin-treated mosquito coils had fewer indoor and outdoor resting Culex spp. Households of repellent non-users in the incomplete coverage scenario had 19 % more Culex spp. mosquitoes in their dwellings although the HBI was not significantly different. The variety of Culex spp. mosquitoes that occur in this area  with potentially different behavioural responses to transfluthrin spatial repellents make it unwise to draw conclusions on diversion for this particular genus as a whole. Unfortunately the study did not identify the culicines to the species level. It would have been interesting to quantify diversion for Culex quinquefasciatus in particular, as it is an important disease vector and nuisance agent. It is useful for laboratory and semi-field studies to characterise the behavioural responses of Cx. quinqefasciatus to spatial repellents because the nuisance caused by this mosquito may compromise the acceptability of the spatial repellent in the community if mosquito bites are still being perceived.
Results from the human landing collections show that a considerable proportion of exposure to malaria vectors is occurring during the early evening hours. Anopheles arabiensis biting densities were highest at 21:00 h, which coincides with the mean bedtime in the area. During recent years there have been consistent reports of vector biting activity shifting towards early evening hours in Kenya [7, 28], Mozambique , Zambia  as well as Tanzania [31, 32]. Classically malaria vectors in sub-Saharan Africa were mostly host-seeking in the middle of the night . The shift to early evening biting is thought to occur due to high levels of LLIN coverage which has pressured the vector population to modify biting patterns in order to feed on a human host when he or she is available . Data from previous studies in the nearby village of Namwawala, show that since the 1990s there is an increasing trend in outdoor biting and early evening activity of malaria vectors . The village of Mbingu lies approximately 20 km southwest of Namwawala and although retrospective data are not specifically available for Mbingu it is likely that a similar phenomenon has occurred. Also, changes in species composition of the main malaria vectors in Mbingu have taken place. A study conducted in the same villages in 2010 reported collecting 83 % An. gambiae (s.s.), 17 % An. arabiensis and absence of An. funestus . Since then An. gambiae (s.s.) has nearly disappeared and been replaced by An. arabiensis. It is also remarkable how An. funestus progressed from nearly absent to one of the main vector species of the Kilombero valley within two years despite universal bed-net coverage being in place [11, 35, 36]. Given the classical anthropophilic and endophilic character of An. funestus it would have been expected that LLINs reduce the vector densities as happened with An. gambiae (s.s.). It is unlikely that LLINs alone are the reason why An. gambiae (s.s.) has been eliminated from the Kilombero Valley, environmental factors such as competition for breeding sites and weather may have also played a role. The present study observed a large proportion of An. funestus feeding on animals with HBI ranging between 0.66 and 0.71. Classically An. funestus is an almost strictly anthropophilic vector with HBI around 0.95 [37, 38]. Between 1990 and 1992 resting catches were performed just a few kilometres north of Mbingu, hundreds of blood engorged An. funestus were collected and virtually all had fed on humans . The present study reports a variation in the host-feeding patterns of An. funestus in the Kilombero Valley. Given the selective pressure exerted by LLINs An. funestus mosquitoes may have been forced to adapt their host choice to the more available animals, and therewith avoided elimination unlike the highly anthropophilic An. gambiae (s.s.) . This adaptation to non-human hosts has also been observed in anthropophilic An. gambiae (s.s.) in Kenya . Human biting rates for An. funestus were relatively similar indoors and outdoors. Comparable biting patterns have been reported in Uganda . Changes in biting time preferences as well as increases in zoophily show how the host-seeking behaviour of An. funestus is dynamic and may change if pressured. Host-seeking activity was highest in the late hours of the night onto early morning hours. It has been reported that most rural families in the Kilombero Valley start their day around 5–6 am [8, 36] and so may become exposed to early morning biting vectors. Early morning feeding by An. funestus has also been reported in Benin , where up to 70 % of exposure to An. funestus bites occur after 6 am during daylight. In addition, in nearby villages insecticide susceptibility tests indicated high levels of resistance in An. funestus against deltamethrin (87 %), permethrin (65 %), lambda cyhalothrin (74 %), bendiocarb (65 %) and DDT (66 %). Similarly, An. arabiensis showed insecticide resistance to deltamethrin (64 %), permethrin (77 %) and lambda cyhalothrin (42 %) in 2014 , which may be associated with the reduced sensitivity of An. funestus to the pyrethroid transfluthrin. It is certainly necessary to understand this relationship before the consideration of a new pyrethroid-based spatial repellent into vector control.
The epidemiological significance of insecticide resistance on malaria transmission is unclear . Ideally vector control should target the vector in multiple locations and physiological stages. Current strategies only control vectors indoors. Vectors are given sufficient room to maneuver their behavior outside the control zone. If elimination is to be reached vector control has to be played like a game of chess, where the opponent must be fought strategically with all available means targeting the vector at different stages of its life. Successes achieved through LLIN campaigns may be followed by a sharp increase in malaria transmission if strategies do not anticipate changes in the vectors behavior. It is likely that malaria will significantly increase in the Kilombero valley if the current situation is not re-evaluated and new tools implemented within the near future.
A considerable degree of malaria transmission in the Kilombero Valley is occurring outside sleeping hours when people are unprotected. This gap in personal protection from disease vectors needs to be addressed. In areas where the main vector is An. arabiensis, transfluthrin spatial repellents may offer some complimentary control aimed at early evening outdoor biting vectors (23 % of exposure as measured in this study) that are precluded by LLINs since HLC demonstrated that the mosquito coils showed good efficacy in preventing mosquito landings. However, the data did not demonstrate a large decline in either mosquito density or HBI in the clusters of houses where the transfluthrin repellents were used despite showing a good efficacy in preventing mosquito landings as measured by human landing catches. Pyrethroid-based spatial repellents may provide a degree of personal protection; however the overall public health benefit is doubtful and potentially iniquitous as their use may divert malaria vectors to those who are not wealthy or wise enough to use them.
An. funestus proved to be insensitive to volatile pyrethroids, this is quite worrisome given its increasing presence in the region. The behavioural responses of susceptible and resistant An. funestus mosquitoes to pyrethroid volatiles needs to be rigorously characterised in order to understand if pyrethroid-based spatial repellents can play a role in malaria endemic areas where the main vectors are pyrethroid-resistant, or if there is an intrinsic species specific variability of response to transfluthrin. Successful implementation of integrated vector control strategies using volatile pyrethroid spatial repellents will be highly dependent on local vector composition, resistance to pyrethroids, feeding times and host preference. Entomological surveys must be made to areas where the implementation of spatial repellents is being considered in order to assess their suitability before implementation.
CI, Confidence interval; DEET, N,N-Diethyl-meta-toluamide; ELISA, enzyme linked immunosorbent assay; HBI, human blood index; HLC, human landing catch; IHI, Ifakara Health Institute; IHI-IRB, Ifakara Health Institute Institutional Review Board; IQR, inter-quartile range; IRR, incidence rate ratio; ITNs, insecticide-treated nets; LLINs, long-lasting insecticide-treated nets; mRDT, malaria rapid diagnostic test; NIMR, National Institute of Medical Research; OR, Odds ratio; PCR, Polymerase chain reaction; RR, Rate ratio
Lengeler C. Insecticide-treated bednets and curtains for preventing malaria. Cochrane Database Syst Rev. 2004;2:CD000363.
Pates H, Curtis CF. Mosquito behaviour and vector control. Ann Rev Entomol. 2005;50:53–70.
Alaii JA, van den Borne HW, Kachur SP, Shelley K, Mwenesi H, Vulule JM, et al. Community reactions to the introduction of permethrin-treated bed nets for malaria control during a randomized controlled trial in western Kenya. Am J Trop Med Hyg. 2003;68:128–36.
Charlwood JD. Biological variation in Anopheles darlingi (Root). Mem Inst Oswaldo Cruz. 1996;91:391–8.
Takken W. Do insecticide-treated bednets have an effect on malaria vectors? Trop Med Int Health. 2002;7:1022–30.
Magesa SM, Wilkes TJ, Mnzava AE, Njunwa KJ, Myamba J, Kivuyo MD, et al. Trial of pyrethroid impregnated bednets in an area of Tanzania holoendemic for malaria. Part 2. Effects on the malaria vector population. Acta Tropica. 1991;49:97–108.
Mbogo CN, Baya NM, Ofulla AV, Githure JI, Snow RW. The impact of permethrin-impregnated bednets on malaria vectors of the Kenyan coast. Med Vet Entomol. 1996;10:251–9.
Killeen GF, Kihonda J, Lyimo E, Oketch FR, Kotas ME, Mathenge E, et al. Quantifying behavioural interactions between humans and mosquitoes: evaluating the protective efficacy of insecticidal nets against malaria transmission in rural Tanzania. BMC Infect Dis. 2006;6:161.
White GB, Moore SJ. Terminology of insect repellents. In: Debboun M, Frances S, Strickman D, editors. Insect Repellent: Principles, Methods, and Uses. Second edition. Boca Raton: CRC Press; 2015. p. 3.
Ogoma SB, Moore SJ, Maia MF. A systematic review of mosquito coils and passive emanators: defining recommendations for spatial repellency testing methodologies. Parasit Vectors. 2012;5:287.
Ogoma S, Lorenz LM, Ngonyani H, Sangusangu R, Kitumbukile M, Kilalangongono M, et al. An experimental hut study to quantify the effect of DDT and airborne pyrethroids on entomological parameters of malaria transmission. Malar J. 2014;13:131.
Maia MF, Onyango SP, Thele M, Simfukwe ET, Turner EL, Moore SJ. Do topical repellents divert mosquitoes within a community? Health equity implications of topical repellents as a mosquito bite prevention tool. PLoS One. 2013;8, e84875.
Onyango SP, Lwetiojera DW, Simfukwe E, Ngonyani H, Mbeleyela E, Lugiko D, et al. Use of a semi-field system to evaluate the efficacy of topical repellents under user conditions provides a disease exposure free technique comparable with field data. Malar J. 2014;13:159.
Ogoma SB, Ngonyani H, Simfukwe ET, Mseka A, Moore J, Maia MF, et al. The mode of action of spatial repellents and their impact on vectorial capacity of Anopheles gambiae sensu stricto. PLoS One. 2014;9, e110433.
Maia MF, Robinson A, John A, Mgando J, Simfukwe E, Moore SJ. Comparison of the CDC Backpack aspirator and the Prokopack aspirator for sampling indoor- and outdoor-resting mosquitoes in southern Tanzania. Parasit Vectors. 2011;4:124.
Gillies MT, Coetzee M. A Supplement to the Anophelinae of Africa South of the Sahara (Afrotropical region). Publ S Afr Inst Med Res. 1987;55:143.
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.
Cohuet A, Simard F, Toto JC, Kengne P, Coetzee M, Fontenille D. Species identification within the Anopheles funestus group of malaria vectors in Cameroon and evidence for a new species. Am J Trop Med Hyg. 2003;69:200–5.
Gimnig JE, Walker ED, Otieno P, Kosgei J, Olang G, Ombok M, et al. Incidence of malaria among mosquito collectors conducting human landing catches in western Kenya. Am J Trop Med Hyg. 2013;88:301–8.
R_Core_Team. R: A language and environment for statistical computing. In: Book R: A language and environment for statistical computing; 2016. https://cran.r-project.org/doc/manuals/r-release/fullrefman.pdf
lme4. Linear mixed-effects models using S4 classes [http://CRAN.R-project.org/package=lme4]
StataCorp. Stata Statistical Software: Release 12. College Station, TX: StataCorp LP. 2011.
Bolker BM, Brooks ME, Clark CJ, Geange SW, Poulsen JR, Stevens MH, et al. Generalized linear mixed models: a practical guide for ecology and evolution. Trends Ecol Evol. 2009;24:127–35.
O'Hara RB, Kotze DJ. Do not log-transform count data. Methods Ecol Evol. 2010;1:118–22.
Eberhardt LL. Appraising variability in population studies. J Wildl Manage. 1978;42:207–38.
Wilson AL, Chen-Hussey V, Logan JG, Lindsay SW. Are topical insect repellents effective against malaria in endemic populations? A systematic review and meta-analysis. Malar J. 2014;13:446.
Wagman JM, Achee NL, Grieco JP. Insensitivity to the spatial repellent action of transfluthrin in Aedes aegypti: a heritable trait associated with decreased insecticide susceptibility. PLoS Negl Trop Dis. 2015;9:e0003726.
Githeko AK, Adungo NI, Karanja DM, Hawley WA, Vulule JM, Seroney IK, et al. Some observations on the biting behavior of Anopheles gambiae s.s. Anopheles arabiensis, and Anopheles funestus and their implications for malaria control. Exp Parasitol. 1996;82:306–15.
Mendis C, Jacobsen JL, Gamage-Mendis A, Bule E, Dgedge M, Thompson R, et al. Anopheles arabiensis and An. funestus are equally important vectors of malaria in Matola coastal suburb of Maputo, southern Mozambique. Med Vet Entomol. 2000;14:171–80.
Fornadel CM, Norris LC, Glass GE, Norris DE. Analysis of Anopheles arabiensis blood feeding behavior in southern Zambia during the two years after introduction of insecticide-treated bed nets. Am J Trop Med Hyg. 2010;83:848–53.
Geissbuhler Y, Chaki P, Emidi B, Govella NJ, Shirima R, Mayagaya V, et al. Interdependence of domestic malaria prevention measures and mosquito-human interactions in urban Dar es Salaam, Tanzania. Malar J. 2007;6:126.
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.
Gillies MT, De Meillon B. Anophelinae of Africa South of Sahara (Ethiopian Zoogeographical region). 2nd ed. Johannesburg: The South African Institute for Medical Research; 1968.
Durnez L, Coosemans M. Residual transmission of malaria: an old issue for new approaches. In: Manguin S, editor. Anopheles mosquitoes - new insights into malaria vectors. Rijeka: Intech; 2013.
Lwetoijera DW, Harris C, Kiware SS, Dongus S, Devine GJ, McCall PJ, et al. Increasing role of Anopheles funestus and Anopheles arabiensis in malaria transmission in the Kilombero Valley, Tanzania. Malar J. 2014;13:331.
Matowo NS, Moore J, Mapua S, Madumla EP, Moshi IR, Kaindoa EW, et al. Using a new odour-baited device to explore options for luring and killing outdoor-biting malaria vectors: a report on design and field evaluation of the Mosquito Landing Box. Parasit Vectors. 2013;6:137.
De Meillon B. The Anophelini of the Ethiopian geographical region. Johannesburg: The South African Institute for Medical Research; 1947. p. 270.
Chandler JA, Boreham PF, Highton RB, Hill MN. A study of the host selection patterns of the mosquitoes of the Kisumu area of Kenya. Trans R Soc Trop Med Hyg. 1975;69:415–25.
Charlwood JD, Smith T, Kihonda J, Heiz B, Billingly PF, Takken W. Density independent feeding success of malaria vectors (Diptera: Culicidae) in Tanzania. Bull Ent Res. 1995;85:29–35.
Mwangangi JM, Mbogo CM, Nzovu JG, Githure JI, Yan G, Beier JC. Blood-meal analysis for anopheline mosquitoes sampled along the Kenyan coast. J Am Mosq Control Assoc. 2003;19:371–5.
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–10.
Kabbale FG, Akol AM, Kaddu JB, Onapa AW. Biting patterns and seasonality of Anopheles gambiae sensu lato and Anopheles funestus mosquitoes in Kamuli District, Uganda. Parasit Vectors. 2013;6:340.
Moiroux N, Damien GB, Egrot M, Djenontin A, Chandre F, Corbel V, et al. Human exposure to early morning Anopheles funestus biting behavior and personal protection provided by long-lasting insecticidal nets. PLoS One. 2014;9, e104967.
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.
The authors thank SC Johnson for donating the coils used in this study. In addition, the authors would like to thank the people of Mbingu village, particularly from Matete, Uwata and Igima for their support and understanding. We also kindly thank Fulbright volunteer Michael Henry for his contribution to linking the field data and the lab results. We specially thank the Mbingu Sisters Convent and Catholic Health Centre for admitting and treating two children of Igima who suffered severe malaria during the course of this study and survived.
We thank the Bill and Melinda Gates Foundation for funding this project (award number 51431).
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article. Raw data are available from the authors upon request.
MFM and SJM designed the experiment. MFM coordinated logistics, staff and managed all the data collected in the field and in the lab. EM supervised the field technicians and mosquito collections. DR managed all the lab samples and carried out the PCRs. VM carried out the ELISAs and handling of blood fed mosquitoes. NFL carried out the PCR for the mosquitoes pooled from the HLC. The data was analyzed by MFM, KK, SJM and AR. The manuscript was prepared by MFM and reviewed by all other authors who read and approved the final version of the manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Ethical approval was obtained from the Ifakara Institutional Review Board (IHI-IRB) Certificate number: IHI/IRB/No:01-2012, from the Tanzanian National Institute for Medical Research (NIMR) Certificate number: NIMR/HQ/R.8a/Vol.IX/1428 and from the Interventions Research Ethics Committee of the London School of Hygiene and Tropical Medicine, ethics reference 6159. All villagers were invited to several community meetings held by the village heads and the research team where the outline of the study and their participating role was explained in the local language. All residents of the participating households were enrolled upon written informed consent; parental consent was obtained for children under the age of 18. Withdrawal of participation without justification was permitted at any time. Each household was given LLINs for all sleeping spaces and a medical officer was employed by the project to offer free health assistance to all participants during the study period. In the case of severe illness (unrelated to the study) participants were taken to a local Catholic Mission where appropriate medical care was provided.
About this article
Cite this article
Maia, M.F., Kreppel, K., Mbeyela, E. et al. A crossover study to evaluate the diversion of malaria vectors in a community with incomplete coverage of spatial repellents in the Kilombero Valley, Tanzania. Parasites Vectors 9, 451 (2016). https://doi.org/10.1186/s13071-016-1738-4