- Open Access
Host-seeking activity of a Tanzanian population of Anopheles arabiensis at an insecticide treated bed net
© The Author(s) 2017
Received: 23 February 2017
Accepted: 26 June 2017
Published: 4 July 2017
Understanding how mosquitoes respond to long lasting insecticide treated nets (LLINs) is fundamental to sustaining the effectiveness of this essential control tool. We report on studies with a tracking system to investigate behaviour of wild anophelines at an LLIN, in an experimental hut at a rural site in Mwanza, Tanzania.
Groups of adult female mosquitoes (n = 10 per replicate) reared from larvae of a local population, identified as predominantly (95%) Anopheles arabiensis, were released in the hut. An infrared video tracking system recorded flight and net contact activity over 1 h as the mosquitoes attempted to reach a supine human volunteer within a bed net (either a deltamethrin-treated LLIN or an untreated control net). A range of activities, including flight path, position in relation to the bed net and duration of net contact, were quantified and compared between treatments.
The total time that female An. arabiensis spent in flight around LLINs was significantly lower than at untreated nets [F(1,10) = 9.26, p = 0.012], primarily due to a substantial reduction in the time mosquitoes spent in persistent ‘bouncing’ flight [F(1,10) = 18.48, p = 0.002]. Most activity occurred at the net roof but significantly less so with LLINs (56.8% of total) than untreated nets [85.0%; Χ2 (15) = 234.69, p < 0.001]. Activity levels at the bed net directly above the host torso were significantly higher with untreated nets (74.2%) than LLINs [38.4%; Χ2 (15) = 33.54, p = 0.004]. ‘Visiting’ and ‘bouncing’ rates were highest above the volunteer’s chest in untreated nets (39.9 and 50.4%, respectively) and LLINs [29.9 and 42.4%; Χ2 (13) = 89.91, p < 0.001; Χ2 (9) = 45.73, p < 0.001]. Highest resting rates were above the torso in untreated nets [77%; Χ2 (9) = 63.12, p < 0.001], but in LLINs only 33.2% of resting occurred here [Χ2 (9) = 27.59, p = 0.001], with resting times spread between the short vertical side of the net adjacent to the volunteer’s head (21.8%) and feet (16.2%). Duration of net contact by a single mosquito was estimated at 204–290 s on untreated nets and 46–82 s on LLINs. While latency to net contact was similar in both treatments, the reduction in activity over 60 min was significantly more rapid for LLINs [F(1,10) = 6.81, p = 0.026], reiterating an ‘attract and kill’ rather than a repellent mode of action.
The study has demonstrated the potential for detailed investigations of behaviour of wild mosquito populations under field conditions. The results validate the findings of earlier laboratory studies on mosquito activity at LLINs, and reinforce the key role of multiple brief contacts at the net roof as the critical LLIN mode of action.
The effectiveness of long lasting insecticide-treated nets (LLINs) in malaria prevention , and their contribution to reductions in malaria incidence and morbidity in Africa , underscore their central role in current and planned malaria prevention and elimination programmes [3–5]. Yet, even as such achievements are being reported, the future of LLINs is under threat from behavioural and physiological changes in mosquito populations; behavioural resistance, such as shifts to feeding outdoors or on animals, could reduce mosquito exposure to LLINs, and the emergence and rapid spread of resistance to the pyrethroid insecticides used on LLINs could decrease the impact of nets [6–8]. Possible solutions to this challenge include using novel insecticides or other treatments delivered via nets, including two-in-one or combination nets [9–12] and novel fabric technologies . Whatever solutions are found to safeguard LLIN effectiveness against pyrethroid-resistant populations, the importance of identifying their mode of action and effect on mosquito behaviour should not be underestimated.
A number of test methods have been used to investigate the behaviour of Anopheles sp. when host seeking at LLINs. Baited WHO tunnel tests have been in use for some years to test effects of insecticide on blood-feeding and mortality of approaching mosquitoes . Modified tunnel tests have been adapted for filming, to allow precise quantification of contact with insecticide, and investigation of the host seeking flight of approaching mosquitoes [15–18]. Other approaches have used techniques of tracking and sticky net traps to investigate mosquito activity around human baited LLINs [19–22]. The optical imaging and flight-tracking system used by Parker et al.  permitted remote tracking, recording and quantitative analysis of the space around an entire bed net, with multiple mosquitoes simultaneously flying without restriction, over long periods, and responding to the human host ‘bait’ within the bed net. In that report, Anopheles gambiae s.s. responses at a human-occupied net were classified into four distinct behavioural modes, termed swooping, visiting, bouncing and resting. Net contact largely comprised multiple brief ‘visits’ or ‘bounces’ on the bed net roof, and the majority of flight paths were above the roof. Behaviour at LLINs was similar to untreated nets inasmuch as insecticide treatment did not repel mosquitoes prior to net contact, although the duration of contact was significantly lower at the LLIN (maximum contact duration in 60 min tests = 96 and 334 s at LLIN and untreated nets, respectively ).
The tracking system of Parker et al.  was installed in an experimental hut in rural Tanzania and used to describe the flight behaviour and nature of the contact at LLINs for Culex quinquefasciatus . Though field tests cannot feasibly control temperature and humidity as is done in laboratory work, a tracking system that produces such high quality information in a natural field setting has many advantages; studies can be carried out on wild rather than colonised populations, minimising the influence of behavioural changes resulting from genetic bottlenecks or unknown selection pressures that might occur during colonisation and the problems associated with laboratory studies of behavioural responses to insecticides can be overcome, e.g. laboratory data are not always consistent with field data [24, 25] and colonisation can result in altered host responses [26, 27] or changes in mosquito responses to insecticides [28–32].
Additional explanations for discrepancies between laboratory and field studies are numerous, but the smaller scale of laboratory test arenas is likely to be a major constraint [33–35], as are the sealed draught-proof laboratory environments that are unlikely to mimic the complex air movements, odour plumes and microclimatic changes that occur in field or more realistic settings [36–39].
The establishment of the field-based tracking system in Tanzania  allowed wild mosquito populations to be studied and compared with laboratory studies of colonised An. gambiae s.s. For this study, we tracked behaviour of Anopheles arabiensis, the dominant species in the area. Although they are closely related members within the same species complex and both are primary malaria vectors in Africa, An. gambiae s.s. and An. arabiensis differ in a number of key bloodfeeding behaviours. Typically, An. gambiae s.s. populations are highly anthropophagic, with strong endophagic and endophilic tendencies, while An. arabiensis populations feed outdoors as well as indoors and will feed on cattle and other hosts in addition to humans [40–42]. Various studies report species-specific preferences for biting certain areas of the human body, and these preferences may vary according to whether the host is supine, seated or standing [43–45]. This raises the possibility that the arrival site and contact rates of An. gambiae s.l. species at an LLIN could be influenced by these preferences.
Finally, responses to insecticides may differ between the two species. Anopheles arabiensis were reported to be more inhibited by various insecticides during blood-feeding than An. gambiae s.s. , while another study controlling for insecticide resistance, showed that both species exhibited similar behavioural responses to permethrin in contact irritancy assays .
We report here on a study using a large scale tracking system  at a rural field site in Mwanza, northern Tanzania to describe and quantify the effects of an LLIN on host-seeking flight of adult female An. arabiensis. Results were used to validate the results of previous laboratory tests, through qualitative comparison with existing data on behaviour of colonised An. gambiae s.s. mosquitoes .
During July and August 2014, anopheline larvae were collected from natural ground water pools in Magu, Tanzania (2°33′41″S 33°18′11″E) and reared to adults in the insectaries at the National Institute of Medical Research in Mwanza. Members of the An. gambiae s.l. species complex were identified at adult emergence using morphological keys [40, 48] and maintained under a 12:12 light: dark cycle (approximating the natural cycle) and provided 10% sugar solution ad libitum. Only the first generation of emerged mosquitoes (i.e. F0) were used in the tests.
Initial trials attempted to track wild mosquitoes entering the hut through the eaves or open door after they had emerged naturally from the surrounding area. However, a mixed population of numerous mosquito species within the Anopheles, Culex and Mansonia genera entered the hut, along with numerous other insects from taxa including Diptera, particularly Chironomidae and Chaoboridae, and Lepidoptera. It has not been possible to reliably classify an individual mosquito track as belonging to a specific mosquito genus based on quantified mosquito flight behaviours . Therefore, the approach to use a mixed population of wild mosquitoes was abandoned.
Species identification and insecticide susceptibility status
The identity of species within the An. gambiae complex was confirmed periodically throughout the testing period using the standard protocol .
In late July (mid-way through the testing period) WHO insecticide susceptibility bioassays were performed on 3–5 day old mosquitoes (n = 22) using 0.05% deltamethrin treated papers, with an exposure time of 1 h following standard methodology . Mortality outcomes at 24 h were corrected according to Abbott’s formula.
Experimental hut and tracking system
The experimental hut was erected in a rice growing area near the village of Kayenze (2°23′43″S, 33°0′5″E), approximately 13 km north-east of Mwanza city, Tanzania. The hut was situated 70 m from the nearest house, and 240 m from a cattle enclosure.
Ten female An. gambiae s.l. mosquitoes, aged 3–5 days post eclosion, were released in each test. Mosquitoes were selected for experiments 2 h prior to testing by placing an arm against a cage and aspirating those that attempted to feed into a paper cup. Mosquitoes were sugar-starved during transport from the insectary to the field site (typically 1–1.5 h). Tests were conducted between the hours of 21:00 and 00:30, to coincide with peak biting periods of An. gambiae complex mosquitoes in this region [51–55].
The volunteer entered the bed net 1 h prior to recording, and the paper cup of mosquitoes was hung at eave height on the wall of the experimental hut. At the start of tests, a cord was pulled from outside the experimental hut to remove the net cover from the paper cup, inverting it and releasing mosquitoes into the room. Activity was recorded for 60 min after which the mosquitoes were collected using a prokopack aspirator .
Mosquito responses were tracked at Permanet 2.0® LLINs (75 denier polyester net; 55 mg/m2 deltamethrin; Vestergaard, Lausanne, Switzerland) and untreated nets (assembled from untreated polyester net of similar mesh), both of which had been tailored slightly to fit the field of view. Ten tests were conducted with each net (i.e. 20 tests in total, using a total of 200 mosquitoes); order of recording mosquitoes at treated and untreated nets was randomised. Up to two tests could be conducted per night, although on days when two tests were conducted the same net was used for both tests.
Members of the local community volunteered to act as human baits; four female and seven male volunteers participated in the ten tests. Half of the volunteers lay with their heads on the left hand side of the field of view, and half with their heads on the right, to control for possible effects of the position of the mosquito release point on net approach.
Tracking mosquito flights
Mosquito activity was recorded using the tracking system described in detail previously [21, 23]. Two cameras filming at 50 fps recorded a total 1.2 × 2.4 m field of view (1.2 × 1.2 m per camera), which was illuminated by two infrared LEDs, with four (two pairs) of large Fresnel lenses (1 × 1.4 m) to collimate the light through the region around the bed net and focus that light into the cameras (Fig. 1b). Each Fresnel lens pair was positioned 1.96 m apart, on either side of the bed. Flight tracking was restricted to the total volume between the two lens pairs combined (Fig. 1b), resulting in a total tracking volume of 1.2 × 2.4 × 1.96 m.
Videos were recorded to sequence (.seq) files using StreamPix software. Mosquito tracks were identified and analysed using custom written Matlab applications .
Mosquito activity and behavioural modes
Mean total activity time (minutes) of female mosquitoes spent in different behavioural modes over 60 min tests in the field hut
The effects of net treatment were investigated statistically only for swooping and bouncing behavioural modes, to avoid analysing multiple non-independent variables. As these behavioural categories are not independent, increases in activity in one mode necessarily reduces activity in the others; hence only two categories were selected for analysis: swooping (non-contact flight) and bouncing (flight involving high levels of net contact). Effects of net treatment on time spent in either behavioural mode were assessed using cluster adjusted regression analysis (StataCorp, 2013), as described above. Activity times spent in all behavioural modes are presented in qualitative summaries in graphs, tables and text.
Flight speed and tortuosity
Flight speed and tortuosity were calculated as described previously  with speed values described only for swooping tracks, because tracks that included net contact were assumed to be slower and would not be equally represented in different net treatments. Tortuosity values were drawn from swooping tracks and the initial sections of net-contacting tracks, up to their point of first contact. Note that since recordings were made in 2D, speed and tortuosity values could potentially underestimate the true values, as track information did not include movement in the z axis. The effect of the explanatory variable of net treatment on the dependent variables of speed and tortuosity was analysed using cluster adjusted regression analysis (StataCorp, 2013), adjusting for repeated tests by the same volunteers.
Distribution of activity on the bed net
The field of view was sub-divided into 16 regions (Fig. 4), and activity allocated to different regions for spatial analysis of movement. As recordings were made in 2D, there was a degree of uncertainty over the precise location of mosquitoes within a region. This ambiguity affected whether tracks were assigned to regions 1–10 (the net surface), or regions 15 and 16 (the space surrounding the net). Contact with the net was identified using physical contact definitions detailed below, and used to allocate tracks to the appropriate region.
In spatial analysis of total activity time, and time spent in the four different behavioural modes (Fig. 4b, d–g), time spent in each region was scaled by region size (s/m2) to account for the different sizes of the various regions. Physical contact times were analysed by unscaled time (seconds).
Spatial preferences were assessed using a generalised linear model that included terms for net treatment and region, and an interaction term between net treatment and region to assess whether activity was evenly distributed across the field of view (SPSS Statistics, version 21, IBM).
Physical contact with the net surface
In this 2D system, mosquito position relative to the net was unknown, and therefore was inferred by examination of track characteristics. Tracks were analysed to find points where mosquitoes made physical contact with the net. Three different track elements were classed as indicating contact with the net; sharp angle turns in tracks (angles of 80° or more [15, 16]); the troughs of bouncing tracks, in which a mosquito changed flight direction at intervals of 0.4 s or less; resting periods, in which a mosquito remained static, or moved at a speed of less than 1.33 mm/s for more than 0.75 s. Examples of all three types of contact can be seen in the Additional file 1: video S1. The total duration of physical contact a mosquito made with the net surface incorporated every contact accrued through all three of these track types.
Effects of net treatment on time spent in physical contact with the net were assessed by cluster adjusted regression analysis (StataCorp, 2013), adjusting for any potential variations resulting from differences in the attractiveness of different volunteers (STATA). As the limits of the camera field of view meant that it was not possible to track individual mosquitoes for the full hour’s test, the range of times that an individual mosquito may have spent in contact with the LLIN was calculated as described in Parker et al. . The maximum of the range makes the assumption that not all mosquitoes released into the room responded to the human bait. The minimum of the range was based on a 100% response rate where all mosquitoes released contacted the net. The maximum of this range was based on a low response rate, calculated by dividing the total contact time by the maximum number of mosquitoes observed responding to the net at one time (i.e. the number of insects we know responded to, and were active around the net). The minimum contact time per individual was calculated by dividing the total duration of physical contact amassed by all mosquitoes by 10 (the total number of mosquitoes released).
Activity rates over time
Repellency was quantified using measures of time lag from mosquito release to first appearance in the field of view, and time elapsed until the first mosquito contacted the net. This value was scored for the first mosquito entry only (i.e. one mosquito per test). These values were evaluated using a Log Rank Mantel-Cox survival analysis in SPSS version 21 (IBM). Two tests were conducted to assess the effect of the explanatory variable of net treatment against the outcome variables of time lag to first appearance, and the time between mosquito release and net contact. If nets were repellent, mosquitoes in treated tests would be expected to take longer to appear in the camera field of view, and take longer to make first contact with the net.
Mosquito activity over the entire recording period (60 min) was assessed for suitability for exponential decay modelling, but many of the tests violated that equation’s constraints. Instead, tests were used to assess the difference between activity in the first and the final 5-min intervals (0–5 and 55–60 min). Total activity recorded in the first interval was subtracted from total activity recorded in the final interval: if negative, the value indicated activity decay, while a positive value indicated that activity had increased between mosquito release and termination of the test. These values were compared using cluster adjusted regression analysis (StataCorp, 2013), adjusting for volunteer effects (STATA) to investigate effect of net treatment on attack persistence.
A total of 142 individuals, morphologically identified as Anopheles gambiae s.l., were processed by PCR, and 129 produced a PCR band. Of these, 95.3% (n = 123) were identified as An. arabiensis, 4.7% (n = 6) as An. gambiae s.s. Hence, during the study period, adult mosquitoes reared from immature stages from this site were referred to as An. arabiensis.
Adult female An. gambiae s.l. from the same population were tested in WHO insecticide bioassays with 0.05% deltamethrin treated papers. Abbott’s correction was applied to 24 h mortality data, as control mortality at this time point was 9.5%. The knockdown rate after 1 h was 95%, the mortality rate after 24 h was 100%, and the population was classified as susceptible to pyrethroids.
Responses of mosquitoes to LLINs and untreated nets
A segment of a video recording from the study, showing a sequence of An. arabiensis flight around a human-occupied LLIN, and including flights tracks on the roof of the net that exhibit characteristic bouncing and visiting patterns, is provided in the Additional file 1: Video S1 (available online). This clip shows mosquito tracks superimposed over a still reference image of the volunteer baited bednet.
Flight speed and tortuosity
Swooping mosquitoes in untreated tests flew at a mean instantaneous speed of 327 mm/s (95% CI 306–348). In LLIN tests, the mean swooping speed was 353 mm/s (95% CI 318–388) and was not significantly different from untreated nets [F(1,10) = 3.09, p = 0.109].
Track tortuosity was not significantly affected by net treatment [F(1,10) = 0.22, p = 0.650]: the mean tortuosity of tracks at an untreated net was 1.35 (95% CI 1.20–1.50), and 1.40 (95% CI 1.27–1.53) at LLINs.
Location of activity at the bed net interface
Spatial analysis was conducted to investigate effect of region, insecticide treatment, and potential interaction between region and insecticide on the outcome activity density (i.e. activity scaled by region size). In these models, insecticide treatment significantly reduced total activity levels [Χ2 (1) = 17.81, p < 0.001]. This decrease was seen in bouncing [Χ2 (1) = 16.01, p < 0.001] and resting [Χ2 (1) = 21.96, p < 0.001] modes, but not swooping [Χ2 (1) = 3.77, p = 0.052] or visiting [Χ2 (1) = 0.92, p = 0.337].
Only 10.4 and 10.9% of all swooping at untreated nets and LLINs respectively occurred in regions 15 and 16 (the spaces in front of the bed net, see Fig. 4a), the lowest levels of this behaviour mode in any region [Χ2 (5) = 66.77, p < 0.001]. Net treatment did not affect distribution of swooping flights [Χ2 (5) = 3.71, p = 0.592; Fig. 4d].
Most visiting activity occurred on the roof of the net above the volunteer’s torso regions 1–3 [39.9% on untreated nets; 29.9% on LLINs; Χ2 (13) = 89.91, p < 0.001]. Regions 7 and 10 accounted for 16.1 and 10.2% respectively of visiting activity on untreated nets, and 13.7 and 20.8% at LLINs. Net treatment did not significantly affect the distribution of visiting behaviour [Χ2 (13) = 10.42, p = 0.659; Fig. 4e].
The majority of bouncing activity occurred in region 2 on both untreated nets (50.4%) and LLINs [42.4%; Χ2 (9) = 45.73, p < 0.001; Fig. 4f]. In contrast, very low levels of bouncing occurred at the lower body portion of the net regions 4–10 (16.1% of untreated bouncing, 31.5% of LLIN bouncing occurred in these seven regions). Net treatment affected distribution of bouncing [Χ2 (9) = 28.14, p = 0.001], with significantly higher bouncing rates at region 2 of the roof of untreated nets (Fig. 4f).
Activity in resting mode was unevenly distributed between net regions [Χ2 (9) = 63.12, p < 0.001]. High levels of resting were observed on region 2 above the volunteer’s chest in both untreated nets (38.6%) and LLINs (21.0%; Fig. 4g). However, there were significant differences [Χ2 (9) = 27.59, p = 0.001] in distribution of resting events according to net treatment: in LLINs, the highest density of resting (21.8%) was recorded on the vertical surface of the net adjacent to the head (region 7), where only 5.6% of resting occurred on untreated nets. Thus, the majority of resting events occurred on regions 1–3 (77% of total) at untreated nets, but at the LLIN only 33.2% of resting occurred here (Fig. 4g).
Quantifying duration of net contact
Since the levels of physical contact with the net are calculated from the combination of visiting, bouncing and resting activity, the distribution and duration of net contact mirrors the relative preferences for those behavioural modes (Fig. 4c). Hence, the highest level of physical contact occurred in region 2 on both net types, where the mean total duration (by ten mosquitoes) of net contact per test was 774 s in untreated nets and 126 s in LLINs [equivalent to 37.7 and 26.6% of total contact time; Χ2 (9) = 30.09, p < 0.001]. The frequencies of net contact at different regions were influenced by net treatment [Χ2 (1) = 20.00, p = 0.011]: in untreated nets, the majority of net contact (76.7%) occurred at roof regions 1–3, whereas in LLINs, net contact occurred at 1–3 (47.0%) and at roof regions 5 and 6, above the volunteer’s feet (21.2%).
Duration of Anopheles arabiensis physical contact with a bed net, during a 60-min test
Duration of physical contact with the bed net surface (60 min test)
Mean total time (all contacts)a (min)
Mean time/mosquito (10 mosquitoes)b (s)
Mean time/mosquito (observed max number)c (s)
Responses to the bed net over time
The results did not provide any evidence for any repellent effect of the LLIN prior to net contact. The delay prior to the first mosquito’s appearance in the field of view was not significantly affected by net treatment [Χ2 (1) = 0.60, p = 0.438], with a geometric mean delay from release to first appearance of 8 s (95% CI 4–14) in untreated nets, and 16 s (95% CI 1–39) in LLINs.
In untreated nets, mosquitoes first contacted the net at a geometric mean of 36 s (95% CI 7–89) after release, compared to 46 s (95% CI 9–119) in LLINs, times that were not significantly different [Χ2 (1) = 0.89, p = 0.766].
The study investigated the behaviour of a wild Anopheles females responding to a human bait within an insecticide-treated bed net, in the field in Africa. The mosquitoes were a pyrethroid susceptible population comprising predominantly An. arabiensis and the results show detailed flight data on host seeking by wild mosquitoes under natural conditions.
Following widespread use of IRS and LLINs, which target endophily and endophagy respectively, An. arabiensis is becoming increasingly important as a vector of residual malaria in Africa [53, 57], a consequence of its high levels of zoophilic, exophagic and exophilic behaviours. Nonetheless, some levels of endophilic activity still occur as evidence indicates that LLINs are still effective against An. arabiensis [58–60].
Comparing responses at LLINs with untreated control nets, insecticide treatment reduced net attack (as measured by the number or frequency of flights), chiefly affecting bouncing and resting, the behavioural modes involving highest levels of net contact. The majority of flight activity occurred on the bed net roof in the area above the volunteer’s torso (regions 1–3; Fig. 4b), though this preference was more pronounced in untreated nets than LLINs. Net contact estimates indicated that an individual mosquito accumulated at least 46 s of physical contact with the LLIN during a 60 min test. Since no evidence was found for insecticide repellency, the LLIN appeared to exert its effect only after direct contact with the insecticide, resulting in reduced activity.
These findings are broadly similar to those recently reported using the same tracking system and a long-established colony of An. gambiae s.s. under laboratory conditions . Considering that fewer An. arabiensis were released in this field study (10 mosquitoes per test) than in the laboratory (25 An. gambiae s.s. per test), activity per mosquito was higher in the An. arabiensis field population: total activity time at untreated tests was 91.9 min for An. arabiensis, compared to 124.6 min for An. gambiae s.s. Insecticide treatment reduced activity to 32% of untreated net values with An. arabiensis in the field, compared to a reduction to 17% with An. gambiae s.s. Whether these differences resulted from colonisation, innate differences between these species, or differences in test conditions and setting, is impossible to determine. However, earlier work has suggested that some mosquito species are more active than others during host-seeking and in behavioural bioassays [20, 23, 61, 62].
As seen in previous laboratory reports with other Anopheles spp., total activity, bouncing, and resting flight densities were highest at the net roof above the volunteer’s torso in all treatments, although insecticide treatment significantly lowered this preference relative to other areas of the net (see track examples in Additional file 1: video S1) [20–22]. Notably, in both the present study and the earlier An. gambiae s.s. study , the proportions of activity on the roof (Fig. 4; regions 1–6) and foot end of untreated nets (region 10) were comparable: An. arabiensis, 85.0 and 4.6%, An. gambiae s.s., 74.7 and 10.9% at roof and foot regions, respectively. At LLINs however, the equivalent values were markedly different between the two studies with roof and foot end activity at 56.8 and 2.0% in the An. arabiensis field study, and 78.3 and 8.8% in the An. gambiae s.s. laboratory study, respectively. Nevertheless, the majority of mosquito activity was focused on the roof of a human-occupied bed net, both here and in the earlier study .
Little activity was observed at the sides of the net, an important observation as this is the region where nets are most likely to be damaged, particularly so at the bottom of the net close to the mattress, where holes are most commonly found . We acknowledge that the 2D nature of our tracking system could have underestimated net contacts in regions 8 and 9, as the sharp angled movements towards and away from the net that comprise ‘visiting’ might not have been visible in every case and some events potentially could have been misclassified as ‘swooping’. However, since total activity (Fig. 4b), which includes swooping, upholds the mosquito preference for the net roof observed when measuring total physical contact times (Fig. 4c), misclassifications would not have been significant. Moreover, sticky net studies, which trap mosquitoes that touch the net, also showed similar attacking patterns [20, 22].
The use of 10 mosquitoes per test, compared with 25 in previous laboratory studies, provided a more precise estimate of net contact times per mosquito. In laboratory tests, a single An. gambiae s.s. accumulated 18–96 s of net contact in a 60 min test , whereas in the field tests with An. arabiensis, the range was 46–82 s. However, in the present study, mortality was not recorded as it was not possible to recapture all mosquitoes in the experimental hut following tests, and those that were collected were often physically damaged by the prokopak aspirator, resulting in elevated mortality rates. Elsewhere, in flight tracking experiments that precisely quantified LLIN contact, An. gambiae that contacted deltamethrin for 40 s or more were knocked down at 1 h, and dead after 24 h, although the number of mosquitoes tested was low (6/35 tested were knocked down) . The insecticide concentration used was comparable with that in the present study and as susceptible An. arabiensis and An. gambiae s.s. display similar knock-down times , it is reasonable to assume that mosquitoes in our study would have been affected 1–24 h after net attack.
Here, and in a previous study , the nature and duration of insecticide exposure at an LLIN have been revealed as being markedly lower than that used in the standard bioassays for assessing KD50 , where mosquitoes are forced into almost continuous contact with insecticides for a constant 3 min. The validity of such tests should be re-evaluated, while investigating the consequence for mosquitoes of more realistic or natural types of exposure should be prioritised, aiming for more accurate measurement of net effectiveness and mosquito susceptibility.
Decay in mosquito activity over the hour’s test was studied to assess whether mosquitoes desisted attacking LLINs following initial contact with insecticide. There was little evidence for such striking temporal decay in activity by An. arabiensis at LLINs as had been reported in laboratory tests with An. gambiae s.s., where high initial activity fell rapidly to negligible levels within 30 min of release . Though less obvious here, where activity of An. arabiensis at LLINs remained at relatively low levels throughout, the reduction in activity over the 60 min of testing was greater at LLINs than with untreated nets (Fig. 5). As previously stated, whether this is the result of differences between a relatively homogeneous colonised strain and a heterogeneous wild population, differences in experimental conditions or between the two species, remains to be determined. No volatile or pre-contact repellent effects of the LLINs had been observed in either study, suggesting that observed behaviours were responses by mosquitoes after LLIN contact. This is supported by a recent report that LLINs did not affect house entry rates . A limitation of our study was that the closed test room design prevented the detection of exiting behaviour following net contact. Clearly, further work on the insecticides used on bed nets is needed to characterise the range of possible effects (e.g. toxic, irritant, sensory impairment) and their impact on mosquitoes (e.g. immediate/delayed, reversible, lethal, consequences post-blood meal, etc.), many of which might be subtler than previously recognised.
Neither this study nor the earlier laboratory study  detected any effects of net treatment on track speed and tortuosity. Speeds of An. arabiensis in all baited tests were between 322 and 355 mm/s, and tortuosity ranged from 1.36 to 1.66. In the laboratory, An. gambiae s.s. speeds were comparable : mean speeds ranged from 321 to 327 mm/s, and mean tortuosity values were between 1.63 and 1.66. Tracking in 2D will lead to underestimation of true speeds as the system does not record movement in the z axis. Nonetheless, flight speeds recorded in the present study with An. arabiensis were comparable to those observed in 3D tracked wind-tunnel host seeking experiments, where upwind velocities ranged from 50 to 444 mm/s, and the highest average downwind flight velocity recorded was 272 mm/s [67, 68].
Conditions in the field hut used in the present study differed to those of the previous laboratory study  in a number of ways. Field tests were conducted in a slightly larger room with a higher ceiling and some air movement, with wild mosquitoes, at night time and with naturally fluctuating temperature, humidity and lunar illumination. However the similarity between these two studies, and indeed a third report using the tracking system , indicate that data obtained in the laboratory can be considered a reliable representation of LLIN performance in a natural context. On this basis, the key role of the bed net roof for targeting Anopheles sp. and Culex sp. and delivering insecticide via multiple brief contacts, is reinforced by the results reported here.
The study has demonstrated the potential for detailed investigations of mosquito behaviour under semi-field conditions. The results validate the findings of earlier laboratory studies on mosquito activity at LLINs, and reinforce the evidence for ‘attract and kill’ via multiple brief contacts at the net roof as the principal mode of action of LLINs.
PJM conceived the study. DPT, NAJ and CET designed and built the tracking system; NAJ and DPT designed segmentation and tracking algorithms; NAJ implemented and tested algorithms; PJM, JEAP and MA designed the experiments; JEAP and JM carried out the experiments; FM supervised field experiments; KG assisted with data extraction and analyses; JEAP analysed the data and drafted the manuscript. All authors read and approved the final manuscript.
We thank Prof Brian Faragher for statistical advice, and John Changalucha for institutional support. We are grateful to David Giesbrecht for assistance in construction of the experimental hut and to Joseph & Sara Mberere for permission to carry out the field studies on their farm. We acknowledge with sincere thanks, the contribution of the late Mr. Benedict Wareba as the project’s driver and translator in Tanzania throughout the study.
The authors declare that they have no competing interests.
This research was supported by the European Union Seventh Framework Programme FP7 (2007–2013) under Grant Agreement No. 265660 ‘AvecNET’ and co-funded by MRC Grant No. MR/M011941/1.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Lengeler C. Insecticide-treated bed nets and curtains for preventing malaria. Cochrane Database Syst Rev. 2004;2:CD000363.Google Scholar
- Bhatt S, Weiss DJ, Cameron E, Bisanzio D, Mappin B, Dalrymple U, Battle KE, Moyes CL, Henry A, Eckhoff PA, et al. The effect of malaria control on Plasmodium falciparum in Africa between 2000 and 2015. Nature. 2015;526:207–11.View ArticlePubMedPubMed CentralGoogle Scholar
- WHO. World Malaria Report 2014. Geneva: World Health Organisation; 2014. p. 2014.Google Scholar
- WHO. World malaria report 2015. Geneva: World Health Organization; 2015.Google Scholar
- Roll Back Malaria. The global malaria action plan. Geneva: Roll Back Malaria Partnership; 2008.Google Scholar
- Strode C, Donegan S, Garner P, Enayati AA, Hemingway J. The impact of pyrethroid resistance on the efficacy of insecticide-treated bed nets against African anopheline mosquitoes: systematic review and meta-analysis. Plos Med. 2014;11:e1001619.View ArticlePubMedPubMed CentralGoogle Scholar
- Ranson H, Lissenden N. Insecticide resistance in African Anopheles mosquitos: a worsening situation that needs urgent action to maintain malaria control. Trends Parasitol. 2016;32:187–96.View ArticlePubMedGoogle Scholar
- Gatton ML, Chitnis N, Churcher T, Donnelly MJ, Ghani AC, Godfray HCJ, Gould F, Hastings I, Marshall J, Ranson H, et al. The importance of mosquito behavioural adaptations to malaria control in Africa. Evol Int J Org Evol. 2013;67:1218–30.View ArticleGoogle Scholar
- Abílio AP, Marrune P, de Deus N, Mbofana F, Muianga P, Kampango A. Bio-efficacy of new long-lasting insecticide-treated bed nets against Anopheles funestus and Anopheles gambiae from central and northern Mozambique. Malar J. 2015;14:352.View ArticlePubMedPubMed CentralGoogle Scholar
- Ngufor C, N’Guessan R, Fagbohoun J, Odjo A, Malone D, Akogbeto M, Rowland MW. Olyset Duo® (a pyriproxyfen and permethrin mixture net): an experimental hut trial against pyrethroid resistant Anopheles gambiae and Culex quinquefasciatus in Southern Benin. PLoS ONE. 2014;9:e93603.View ArticlePubMedPubMed CentralGoogle Scholar
- Pennetier C, Bouraima A, Chandre F, Piameu M, Etang J, Rossignol M, Sidick I, Zogo B, Lacroix MN, Yadav R, et al. Efficacy of Olyset Plus(R), a new long-lasting insecticidal net incorporating permethrin and piperonyl-butoxide against multi-resistant malaria vectors [corrected]. PLoS ONE. 2013;8:e75134.View ArticlePubMedPubMed CentralGoogle Scholar
- Koffi AA, Ahoua Alou LP, Djenontin A, Kabran JP, Dosso Y, Kone A, Moiroux N, Pennetier C. Efficacy of Olyset(R) Duo, a permethrin and pyriproxyfen mixture net against wild pyrethroid-resistant Anopheles gambiae s.s. from Cote d’Ivoire: an experimental hut trial. Parasite. 2015;22:28.View ArticlePubMedPubMed CentralGoogle Scholar
- Andriessen R, Snetselaar J, Suer RA, Osinga AJ, Deschietere J, Lyimo IN, Mnyone LL, Brooke BD, Ranson H, Knols BGJ, Farenhorst M. Electrostatic coating enhances bioavailability of insecticides and breaks pyrethroid resistance in mosquitoes. Proc Natl Acad Sci. 2015;112:12081–6.View ArticlePubMedPubMed CentralGoogle Scholar
- World Health Organisation. Guidelines for laboratory and field testing of long-lasting insecticidal mosquito nets. Geneva: World Health Organisation; 2005.Google Scholar
- Hawkes F, Gibson G. Seeing is believing: the nocturnal malarial mosquito Anopheles coluzzii responds to visual host-cues when odour indicates a host is nearby. Parasit Vectors. 2016;9:320.View ArticlePubMedPubMed CentralGoogle Scholar
- Spitzen J, Ponzio C, Koenraadt CJM, Pates Jamet HV, Takken W. Absence of close-range excitorepellent effects in malaria mosquitoes exposed to deltamethrin-treated bed nets. Am J Trop Med Hyg. 2014;90:1124–32.View ArticlePubMedPubMed CentralGoogle Scholar
- Webster B, Lacey ES, Cardé RT. Waiting with bated breath: opportunistic orientation to human odor in the malaria mosquito, Anopheles gambiae, is modulated by minute changes in carbon dioxide concentration. J Chem Ecol. 2015;41:59–66.View ArticlePubMedGoogle Scholar
- Miller JE, Gibson G. Behavioral response of host-seeking mosquitoes (Diptera: Culicidae) to insecticide-impregnated bed netting: A new approach to insecticide bioassays. J Med Entomol. 1994;31:114–22.View ArticlePubMedGoogle Scholar
- Sutcliffe J, Colborn KL. Video studies of passage by Anopheles gambiae mosquitoes through holes in a simulated bed net: effects of hole size, hole orientation and net environment. Malar J. 2015;14:199.View ArticlePubMedPubMed CentralGoogle Scholar
- Sutcliffe JF, Yin S. Behavioural responses of females of two anopheline mosquito species to human-occupied, insecticide-treated and untreated bed nets. Malar J. 2014;13:294.View ArticlePubMedPubMed CentralGoogle Scholar
- Parker JEA, Angarita-Jaimes N, Abe M, Towers CE, Towers D, McCall PJ. Infrared video tracking of Anopheles gambiae at insecticide-treated bed nets reveals rapid decisive impact after brief localised net contact. Sci Rep. 2015;5:1–14.View ArticleGoogle Scholar
- Lynd A, McCall PJ. Clustering of host-seeking activity of Anopheles gambiae mosquitoes at the top surface of a human-baited bed net. Malar J. 2013;12:267.View ArticlePubMedPubMed CentralGoogle Scholar
- Angarita-Jaimes NC, Parker JEA, Abe M, Mashauri F, Martine J, Towers CE, McCall PJ, Towers DP. A novel video-tracking system to quantify the behaviour of nocturnal mosquitoes attacking human hosts in the field. J R Soc Interface. 2016;13:20150974.View ArticlePubMedPubMed CentralGoogle Scholar
- Grieco JP, Achee NL, Chareonviriyaphap T, Suwonkerd W, Chauhan K, Sardelis MR, Roberts DR. A new classification system for the actions of IRS chemicals traditionally used for malaria control. PLoS ONE. 2007;2:e716.View ArticlePubMedPubMed CentralGoogle Scholar
- N’Guessan R, Asidi A, Boko P, Odjo A, Akogbeto M, Pigeon O, Rowland MW. An experimental hut evaluation of PermaNet® 3.0, a deltamethrin-piperonyl butoxide combination net, against pyrethroid-resistant Anopheles gambiae and Culex quinquefasciatus mosquitoes in southern Benin. Trans R Soc Trop Med Hyg. 2010;104:758–65.View ArticlePubMedGoogle Scholar
- Laarman JJ. The host-seeking behaviour of anopheline mosquitoes. Trop Geogr Med. 1958;10:293–305.PubMedGoogle Scholar
- Clark GG, Bernier UR, Allan SA, Kline DL, Golden FV. Changes in host-seeking behavior of Puerto Rican Aedes aegypti after colonization. J Med Entomol. 2011;48:533–7.View ArticlePubMedGoogle Scholar
- Chareonviriyaphap T, Roberts DR, Andre RG, Harlan HJ, Manguin S, Bangs MJ. Pesticide avoidance behavior in Anopheles albimanus, a malaria vector in the Americas. J Am Mosq Control Assoc. 1997;13:171–83.PubMedGoogle Scholar
- Sathantriphop S, Ketavan C, Prabaripai A, Visetson S, Bangs MJ, Akratanakul P, Chareonviriyaphap T. Susceptibility and avoidance behavior by Culex quinquefasciatus say to three classes of residual insecticides. J Vector Ecol. 2006;31:266–74.View ArticlePubMedGoogle Scholar
- Sathantriphop S, Thanispong K, Sanguanpong U, Achee NL, Bangs MJ, Chareonviriyaphap T. Comparative behavioral responses of pyrethroid—susceptible and—resistant Aedes aegypti (Diptera: Culicidae) populations to citronella and eucalyptus oils. J Med Entomol. 2014;51:1182–91.View ArticlePubMedGoogle Scholar
- Thanispong K, Achee NL, Bangs M, Grieco JP. Irritancy and repellency behavioral responses of three strains of Aedes aegypti exposed to DDT and α-cypermethrin. J Med Entomol. 2009;46:1407–14.View ArticlePubMedGoogle Scholar
- Malaithong N, Tisgratog R, Tainchum K, Prabaripai A, Juntarajumnong W, Bangs MJ, Chareonviriyaphap T. Locomotor behavioral responses of Anopheles minimus and Anopheles harrisoni to alpha-cypermethrin in Thailand. J Am Mosq Control Assoc. 2011;27:217–26.View ArticlePubMedGoogle Scholar
- Hossain MI, Curtis CF. Permethrin-impregnated bed nets: behavioural and killing effects on mosquitoes. Med Vet Entomol. 1989;3:367–76.View ArticlePubMedGoogle Scholar
- Barnard DR, Posey KH, Smith D, Schreck CE. Mosquito density, biting rate and cage size effects on repellent tests. Med Vet Entomol. 1998;12:39–45.View ArticlePubMedGoogle Scholar
- McMeniman Conor J, Corfas Román A, Matthews Benjamin J, Ritchie Scott A, Vosshall Leslie B. Multimodal integration of carbon dioxide and other sensory cues drives mosquito attraction to humans. Cell. 2014;156:1060–71.View ArticlePubMedPubMed CentralGoogle Scholar
- Murlis J, Elkinton JS, Cardé RT. Odor plumes and how insects use them. Annu Rev Entomol. 1992;37:505–32.View ArticleGoogle Scholar
- Geier M, Bosch O, Boeckh J. Influence of odour plume structure on upwind flight of mosquitoes towards hosts. J Exp Biol. 1999;202:1639–48.PubMedGoogle Scholar
- Ferguson HM, Ng’habi K, Walder T, Kadungula D, Moore SJ, Lyimo IN, Russell TL, Urassa H, Mshinda H, Killeen GF, Knols BGJ. Establishment of a large semi-field system for experimental study of African malaria vector ecology and control in Tanzania. Malar J. 2008;7:158.View ArticlePubMedPubMed CentralGoogle Scholar
- Lacey ES, Cardé RT. Location of and landing on a source of human body odour by female Culex quinquefasciatus in still and moving air. Physiol Entomol. 2012;37:153–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Gillies MT, De Meillon B. The Anophelinae of Africa South of the Sahara (Ethiopian zoogeographical region). 2nd ed. Johannesburg: South African Institute for Medical Research; 1968.Google Scholar
- Pates H, Curtis CF. Mosquito behavior and vector control. Annu Rev Entomol. 2005;50:53–70.View ArticlePubMedGoogle Scholar
- Takken W, Verhulst NO. Host preferences of blood-feeding mosquitoes. Annu Rev Entomol. 2013;58:433–53.View ArticlePubMedGoogle Scholar
- De Jong R, Knols BGJ. Selection of biting sites on man by two malaria mosquito species. Experientia. 1995;51:80–4.View ArticlePubMedGoogle Scholar
- Dekker T, Takken W, Knols BGJ, Bouman E, Laak S, Bever A, Huisman PW. Selection of biting sites on a human host by Anopheles gambiae s.s., An. arabiensis and A. quadriannulatus. Entomol Exp Appl. 1998;87:295–300.View ArticleGoogle Scholar
- Braack L, Hunt R, Koekemoer LL, Gericke A, Munhenga G, Haddow AD, Becker P, Okia M, Kimera I, Coetzee M. Biting behaviour of African malaria vectors: 1. Where do the main vector species bite on the human body? Parasit Vectors. 2015;8:76.View ArticlePubMedPubMed CentralGoogle Scholar
- Malima RC, Oxborough RM, Tungu PK, Maxwell C, Lyimo IN, Mwingira V, Mosha FW, Matowo J, Magesa SM, Rowland MW. Behavioural and insecticidal effects of organophosphate-, carbamate- and pyrethroid-treated mosquito nets against African malaria vectors. Med Vet Entomol. 2009;23:317–25.View ArticlePubMedGoogle Scholar
- Kawada H, Ohashi K, Dida GO, Sonye G, Njenga SM, Mwandawiro C, Minakawa N. Insecticidal and repellent activities of pyrethroids to the three major pyrethroid-resistant malaria vectors in western Kenya. Parasit vectors. 2014;7:208.View ArticlePubMedPubMed CentralGoogle Scholar
- Gillies MT, Coetzee M. A supplement to the Anophelinae of Africa South of the Sahara (Afrotropical Region). Johannesburg: South African Institute for Medical Research; 1987.Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- World Health Organisation. Test procedures for insecticide resistance monitoring in malaria vectors, bio-efficacy and persistence of insecticides on treated surfaces. Geneva: World Health Organisation; 1998.Google Scholar
- Mathenge EM, Gimnig JE, Kolczak M, Ombok M, Irungu LW, Hawley WA. Effect of permethrin-impregnated nets on exiting behavior, blood feeding success, and time of feeding of malaria mosquitoes (Diptera: Culicidae) in western Kenya. J Med Entomol. 2001;38:531–6.View ArticlePubMedGoogle Scholar
- Fornadel CM, Norris LC, Glass GE, Norris DE. Analysis of Anopheles arabiensis blood feeding behavior in southern Zambia during the 2 years after introduction of insecticide-treated bed nets. Am J Trop Med Hyg. 2010;83:848–53.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Yohannes M, Boelee E. Early biting rhythm in the afro-tropical vector of malaria, Anopheles arabiensis, and challenges for its control in Ethiopia. Med Vet Entomol. 2012;26:103–5.View ArticlePubMedGoogle Scholar
- Bayoh MN, Walker ED, Kosgei J, Ombok M, Olang GB, Githeko AK, Killeen GF, Otieno P, Desai M, Lobo NF, et al. Persistently high estimates of late night, indoor exposure to malaria vectors despite high coverage of insecticide treated nets. Parasit Vectors. 2014;7:380.View ArticlePubMedGoogle Scholar
- Vazquez-Prokopec GM, Galvin WA, Kelly R, Kitron U. A new, cost-effective, battery-powered aspirator for adult mosquito collections. J Med Entomol. 2009;46:1256–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Derua YA, Alifrangis M, Hosea KM, Meyrowitsch DW, Magesa SM, Pedersen EM, Simonsen PE. Change in composition of the Anopheles gambiae complex and its possible implications for the transmission of malaria and lymphatic filariasis in north-eastern Tanzania. Malar J. 2012;11:188.View ArticlePubMedPubMed CentralGoogle Scholar
- Okumu F, Mbeyela E, Lingamba G, Moore J, Ntamatungiro AJ, Kavishe DR, Kenward MG, Turner EL, Lorenz L, Moore SJ. Comparative field evaluation of combinations of long-lasting insecticide treated nets and indoor residual spraying, relative to either method alone, for malaria prevention in an area where the main vector is Anopheles arabiensis. Parasit Vectors. 2013;6:46.View ArticlePubMedPubMed CentralGoogle Scholar
- Kawada H, Ohashi K, Dida GO, Sonye G, Njenga SM, Mwandawiro C, Minakawa N. Preventive effect of permethrin-impregnated long-lasting insecticidal nets on the blood feeding of three major pyrethroid-resistant malaria vectors in western Kenya. Parasit Vectors. 2014;7:383.View ArticlePubMedPubMed CentralGoogle Scholar
- Killeen GF, Govella NJ, Lwetoijera DW, Okumu FO. Most outdoor malaria transmission by behaviourally-resistant Anopheles arabiensis is mediated by mosquitoes that have previously been inside houses. Malar J. 2016;15:225.View ArticlePubMedPubMed CentralGoogle Scholar
- Chareonviriyaphap T, Prabaripai A, Bangs MJ. Excito-repellency of deltamethrin on the malaria vectors, Anopheles minimus, Anopheles dirus, Anopheles swadiwongporni, and Anopheles maculatus, in Thailand. J Am Mosq Control Assoc. 2004;20:45–54.PubMedGoogle Scholar
- Cooperband MF, Allan SA. Effects of different pyrethroids on landing behavior of female Aedes aegypti, Anopheles quadrimaculatus, and Culex quinquefasciatus mosquitoes (Diptera: Culicidae). J Med Entomol. 2009;46:292–306.View ArticlePubMedGoogle Scholar
- World Health Organisation. Report of the Twelfth WHOPES working group meeting, vol. WHO/HTM/NTD/WHOPES/2009.1. Geneva: WHO; 2008.Google Scholar
- Vanden Eng JL, Chan A, Abilio AP, Wolkon A, Ponce de Leon G, Gimnig J, Morgan J. Bed net durability assessments: exploring a composite measure of net damage. PLoS ONE. 2015;10:e0128499.View ArticlePubMedPubMed CentralGoogle Scholar
- World Health Organization. Test procedures for insecticide resistance monitoring in malaria vector mosquitoes. Geneva: World Health Organization; 2013.Google Scholar
- Spitzen J, Koelewijn T, Mukabana WR, Takken W. Effect of insecticide-treated bed nets on house-entry by malaria mosquitoes: the flight response recorded in a semi-field study in Kenya. Acta Trop. 2017;172:180–5.View ArticlePubMedGoogle Scholar
- Takken W, Dekker T, Wijnholds YG. Odor-mediated flight behavior of Anopheles gambiae giles Sensu Stricto and An. stephensi liston in response to CO2, acetone, and 1-octen-3-ol (Diptera: Culicidae). J Insect Behavior. 1997;10:395–407.View ArticleGoogle Scholar
- Spitzen J, Spoor CW, Grieco F, Ter Braak C, Beeuwkes J, van Brugge SP, Kranenbarg S, Noldus LPJJ, van Leeuwen JL, Takken W. A 3D analysis of flight behavior of Anopheles gambiae sensu stricto malaria mosquitoes in response to human odor and heat. PLoS ONE. 2013;8:e62995.View ArticlePubMedPubMed CentralGoogle Scholar