Persistently high proportions of plasmodium-infected Anopheles funestus mosquitoes in two villages in the Kilombero valley, South-Eastern Tanzania

Distribution and insecticide resistance profile of the major malaria vector Anopheles funestus group across the African continent

There has been significant progress in malaria control in the last 2 decades, with a decline in mortality and morbidity. However, these gains are jeopardised by insecticide resistance, which negatively impacts the core interventions, such as insecticide-treated nets (ITN) and indoor residual spraying (IRS). While most malaria control and research efforts are still focused on Anopheles gambiae complex mosquitoes, Anopheles funestus remains an important vector in many countries and, in some cases, contributes to most of the local transmission. As countries move towards malaria elimination, it is important to ensure that all dominant vector species, including An. funestus, an important vector in some countries, are targeted. The objective of this review is to compile and discuss information related to A. funestus populations' resistance to insecticides and the mechanisms involved across Africa, emphasising the sibling species and their resistance profiles in relation to malaria elimination goals. Data on insecticide resistance in An. funestus malaria vectors in Africa were extracted from published studies. Online bibliographic databases, including Google Scholar and PubMed, were used to search for relevant studies. Articles published between 2000 and May 2023 reporting resistance of An. funestus to insecticides and associated mechanisms were included. Those reporting only bionomics were excluded. Spatial variation in species distribution and resistance to insecticides was recorded from 174 articles that met the selection criteria. It was found that An. funestus was increasingly resistant to the four classes of insecticides recommended by the World Health Organisation for malaria vector control; however, this varied by country. Insecticide resistance appears to reduce the effectiveness of vector control methods, particularly IRS and ITN. Biochemical resistance due to detoxification enzymes (P450s and glutathione-S-transferases [GSTs]) in An. funestus was widely recorded. However, An. funestus in Africa remains susceptible to other insecticide classes, such as organophosphates and neonicotinoids. This review highlights the increasing insecticide resistance of An. funestus mosquitoes, which are important malaria vectors in Africa, posing a significant challenge to malaria control efforts. While An. funestus has shown resistance to the recommended insecticide classes, notably pyrethroids and, in some cases, organochlorides and carbamates, it remains susceptible to other classes of insecticides such as organophosphates and neonicotinoids, providing potential alternative options for vector control strategies. The study underscores the need for targeted interventions that consider the population structure and geographical distribution of An. funestus, including its sibling species and their insecticide resistance profiles, to effectively achieve malaria elimination goals.

Reagent-free detection of Plasmodium falciparum malaria infections in field-collected mosquitoes using mid-infrared spectroscopy and machine learning

Field-derived metrics are critical for effective control of malaria, particularly in sub-Saharan Africa where the disease kills over half a million people yearly. One key metric is entomological inoculation rate, a direct measure of transmission intensities, computed as a product of human biting rates and prevalence of Plasmodium sporozoites in mosquitoes. Unfortunately, current methods for identifying infectious mosquitoes are laborious, time-consuming, and may require expensive reagents that are not always readily available. Here, we demonstrate the first field-application of mid-infrared spectroscopy and machine learning (MIRS-ML) to swiftly and accurately detect Plasmodium falciparum sporozoites in wild-caught Anopheles funestus, a major Afro-tropical malaria vector, without requiring any laboratory reagents. We collected 7178 female An. funestus from rural Tanzanian households using CDC-light traps, then desiccated and scanned their heads and thoraces using an FT-IR spectrometer. The sporozoite infections were confirmed using enzyme-linked immunosorbent assay (ELISA) and polymerase chain reaction (PCR), to establish references for training supervised algorithms. The XGBoost model was used to detect sporozoite-infectious specimen, accurately predicting ELISA and PCR outcomes with 92% and 93% accuracies respectively. These findings suggest that MIRS-ML can rapidly detect P. falciparum in field-collected mosquitoes, with potential for enhancing surveillance in malaria-endemic regions. The technique is both fast, scanning 60-100 mosquitoes per hour, and cost-efficient, requiring no biochemical reactions and therefore no reagents. Given its previously proven capability in monitoring key entomological indicators like mosquito age, human blood index, and identities of vector species, we conclude that MIRS-ML could constitute a low-cost multi-functional toolkit for monitoring malaria risk and evaluating interventions.

Rapid assessment of the blood-feeding histories of wild-caught malaria mosquitoes using mid-infrared spectroscopy and machine learning

BACKGROUND

The degree to which Anopheles mosquitoes prefer biting humans over other vertebrate hosts, i.e. the human blood index (HBI), is a crucial parameter for assessing malaria transmission risk. However, existing techniques for identifying mosquito blood meals are demanding in terms of time and effort, involve costly reagents, and are prone to inaccuracies due to factors such as cross-reactivity with other antigens or partially digested blood meals in the mosquito gut. This study demonstrates the first field application of mid-infrared spectroscopy and machine learning (MIRS-ML), to rapidly assess the blood-feeding histories of malaria vectors, with direct comparison to PCR assays.

METHODS AND RESULTS

Female Anopheles funestus mosquitoes (N = 1854) were collected from rural Tanzania and desiccated then scanned with an attenuated total reflectance Fourier-transform Infrared (ATR-FTIR) spectrometer. Blood meals were confirmed by PCR, establishing the 'ground truth' for machine learning algorithms. Logistic regression and multi-layer perceptron classifiers were employed to identify blood meal sources, achieving accuracies of 88%-90%, respectively, as well as HBI estimates aligning well with the PCR-based standard HBI.

CONCLUSIONS

This research provides evidence of MIRS-ML effectiveness in classifying blood meals in wild Anopheles funestus, as a potential complementary surveillance tool in settings where conventional molecular techniques are impractical. The cost-effectiveness, simplicity, and scalability of MIRS-ML, along with its generalizability, outweigh minor gaps in HBI estimation. Since this approach has already been demonstrated for measuring other entomological and parasitological indicators of malaria, the validation in this study broadens its range of use cases, positioning it as an integrated system for estimating pathogen transmission risk and evaluating the impact of interventions.

Rapid classification of epidemiologically relevant age categories of the malaria vector, Anopheles funestus

BACKGROUND

Accurately determining the age and survival probabilities of adult mosquitoes is crucial for understanding parasite transmission, evaluating the effectiveness of control interventions and assessing disease risk in communities. This study was aimed at demonstrating the rapid identification of epidemiologically relevant age categories of Anopheles funestus, a major Afro-tropical malaria vector, through the innovative combination of infrared spectroscopy and machine learning, instead of the cumbersome practice of dissecting mosquito ovaries to estimate age based on parity status.

METHODS

Anopheles funestus larvae were collected in rural south-eastern Tanzania and reared in an insectary. Emerging adult females were sorted by age (1-16 days old) and preserved using silica gel. Polymerase chain reaction (PCR) confirmation was conducted using DNA extracted from mosquito legs to verify the presence of An. funestus and to eliminate undesired mosquitoes. Mid-infrared spectra were obtained by scanning the heads and thoraces of the mosquitoes using an attenuated total reflection-Fourier transform infrared (ATR-FT-IR) spectrometer. The spectra (N = 2084) were divided into two epidemiologically relevant age groups: 1-9 days (young, non-infectious) and 10-16 days (old, potentially infectious). The dimensionality of the spectra was reduced using principal component analysis, and then a set of machine learning and multi-layer perceptron (MLP) models were trained using the spectra to predict the mosquito age categories.

RESULTS

The best-performing model, XGBoost, achieved overall accuracy of 87%, with classification accuracy of 89% for young and 84% for old An. funestus. When the most important spectral features influencing the model performance were selected to train a new model, the overall accuracy increased slightly to 89%. The MLP model, utilizing the significant spectral features, achieved higher classification accuracy of 95% and 94% for the young and old An. funestus, respectively. After dimensionality reduction, the MLP achieved 93% accuracy for both age categories.

CONCLUSIONS

This study shows how machine learning can quickly classify epidemiologically relevant age groups of An. funestus based on their mid-infrared spectra. Having been previously applied to An. gambiae, An. arabiensis and An. coluzzii, this demonstration on An. funestus underscores the potential of this low-cost, reagent-free technique for widespread use on all the major Afro-tropical malaria vectors. Future research should demonstrate how such machine-derived age classifications in field-collected mosquitoes correlate with malaria in human populations.

Non-Coding RNAs Potentially Involved in Pyrethroid Resistance of Anopheles funestus Population in Western Kenya

Backgrounds The resurgence of Anopheles funestus , a dominant vector of human malaria in western Kenya was partly attributed to insecticide resistance. However, evidence on the molecular basis of pyrethroid resistance in western Kenya is limited. Noncoding RNAs (ncRNAs) form a vast class of RNAs that do not code for proteins and are ubiquitous in the insect genome. Here, we demonstrated that multiple ncRNAs could play a potential role in An. funestus resistance to pyrethroid in western Kenya. Materials and Methods Anopheles funestus mosquitoes were sampled by aspiration methods in Bungoma, Teso, Siaya, Port Victoria and Kombewa in western Kenya. The F1 progenies were exposed to deltamethrin (0.05%), permethrin (0.75%), DDT (4%) and pirimiphos-methyl (0.25%) following WHO test guidelines. A synergist assay using piperonyl butoxide (PBO) (4%) was conducted to determine cytochrome P450s' role in pyrethroid resistance. RNA-seq was conducted on a combined pool of specimens that were resistant and unexposed, and the results were compared with those of the FANG susceptible strain. This approach aimed to uncover the molecular mechanisms underlying pyrethroid resistance. Results Pyrethroid resistance was observed in all the sites with an average mortality rate of 57.6%. Port Victoria had the highest level of resistance to permethrin (MR=53%) and deltamethrin (MR=11%) pyrethroids. Teso had the lowest level of resistance to permethrin (MR=70%) and deltamethrin (MR=87%). Resistance to DDT was observed only in Kombewa (MR=89%) and Port Victoria (MR=85%). A full susceptibility to P-methyl (0.25%) was observed in all the sites. PBO synergist assay revealed high susceptibility (>98%) to the pyrethroids in all the sites except for Port Victoria (MR=96%, n=100). Whole transcriptomic analysis showed that most of the gene families associated with pyrethroid resistance comprised non-coding RNAs (67%), followed by imipenemase (10%),cytochrome P450s (6%), cuticular proteins (5%), olfactory proteins (4%), glutathione S-transferases (3%), UDP-glycosyltransferases (2%), ATP-binding cassettes (2%) and carboxylesterases(1%). Conclusions This study unveils the molecular basis of insecticide resistance in An. funestus in western Kenya, highlighting for the first time the potential role of non-coding RNAs in pyrethroid resistance. Targeting non-coding RNAs for intervention development could help in insecticide resistance management.

Evaluation of the DN-Mini (miniaturized double net) trap for sampling host-seeking Anopheles mosquitoes in malaria-endemic villages of southern Tanzania

BACKGROUND

Surveillance of malaria vectors is crucial for assessing the transmission risk and impact of control measures. Human landing catches (HLC) directly estimate the biting rates but raise ethical concerns due to the exposure of volunteers to mosquito-borne pathogens. A common alternative is the CDC-light trap, which is effective for catching host-seeking mosquitoes indoors but not outdoors. New, exposure-free methods are needed for sampling mosquitoes indoors and outdoors in ways that reflect their natural risk profiles. The aim of this study was therefore to evaluate the efficacy of the miniaturized double net trap (DN-Mini) for sampling host-seeking mosquitoes in south-eastern Tanzania, where malaria transmission is dominated by Anopheles funestus.

METHODS

Adult mosquitoes were collected from 222 randomly selected houses across three villages (74 per village) in Ulanga district, south-eastern Tanzania, using the DN-Mini traps, CDC-Light traps, and Prokopack aspirators. First, we compared CDC-light and DN-Mini traps for collecting indoor host-seeking mosquitoes, while Prokopack aspirators were used for indoor-resting mosquitoes. Second, we deployed the DN-Mini and Prokopack aspirators to collect host-seeking and resting mosquitoes indoors and outdoors. Generalized linear mixed models (GLMM) with a negative binomial distribution were used to compare the effectiveness of the traps for catching different mosquito species.

RESULTS

The DN-Mini was 1.53 times more efficient in collecting An. funestus indoors (RR = 1.53, 95% CI: 1.190-1.98) compared to the CDC-Light trap. However, for Anopheles arabiensis, the DN-Mini caught only 0.32 times as many mosquitoes indoors as the CDC-Light traps (RR = 0.32, 95% CI: 0.183-0.567). Both An. funestus and An. arabiensis were found to be more abundant indoors than outdoors when collected using the DN-Mini trap. Similarly, the Prokopack aspirator was greater indoors than outdoors for both An. funestus and An. arabiensis.

CONCLUSION

The DN-Mini outperformed the CDC-light trap in sampling the dominant malaria vector, An. funestus species, but was less effective in capturing An. arabiensis, and for both vector species, the biting risk was greater indoors than outdoors when measured using the DN-Mini trap. These findings highlight the importance of selecting appropriate trapping methods based on mosquito species and behaviors.

Geospatial modelling of dry season habitats of the malaria vector, Anopheles funestus, in south-eastern Tanzania

BACKGROUND

Anopheles funestus is a major malaria vector in Eastern and Southern Africa and is currently the dominant malaria-transmitting vector in many parts of Tanzania. Previous research has identified its preference for specific aquatic habitats, especially those that persist in dry months. This observation suggests the potential for targeted control through precise habitat mapping and characterization. In this study, we investigated the influence of habitat characteristics, land cover and human population densities on An. funestus distribution during dry seasons. Based on the results, we developed a habitat suitability model for this vector species in south-eastern Tanzania.

METHODS

Eighteen villages in south-eastern Tanzania were surveyed during the dry season from September-December 2021. Water bodies were systematically inspected for mosquito larvae and characterized by their physico-chemical characteristics and surrounding environmental features. A generalized linear model was used to assess the presence of An. funestus larvae as a function of the physico-chemical characteristics, land use and human population densities. The results obtained from this model were used to generate spatially explicit predictions of habitat suitability in the study districts.

RESULTS

Of the 1466 aquatic habitats surveyed, 440 were positive for An. funestus, with river streams having the highest positivity (74%; n = 322) followed by ground pools (15%; n = 67). The final model had an 83% accuracy in predicting positive An. funestus habitats, with the most important characteristics being permanent waters, clear waters with or without vegetation or movement and shading over the habitats. There was also a positive association of An. funestus presence with forested areas and a negative association with built-up areas. Human population densities had no influence on An. funestus distribution.

CONCLUSIONS

The results of this study underscore the crucial role of both the specific habitat characteristics and key environmental factors, notably land cover, in the distribution of An. funestus. In this study area, An. funestus predominantly inhabits river streams and ground pools, with a preference for clear, perennial waters with shading. The strong positive association with more pristine environments with tree covers and the negative association with built-up areas underscore the importance of ecological transitions in vector distribution and malaria transmission risk. Such spatially explicit predictions could enable more precise interventions, particularly larval source management, to accelerate malaria control.

Influence of larval growth and habitat shading on retreatment frequencies of biolarvicides against malaria vectors

Effective larviciding for malaria control requires detailed studies of larvicide efficacies, aquatic habitat characteristics, and life history traits of target vectors. Mosquitoes with brief larval phases present narrower timeframes for biolarvicidal effects than mosquitoes with extended periods. We evaluated two biolarvicides, VectoBac (Bacillus thuringiensis israelensis (Bti)) and VectoMax (Bti and Bacillus sphaericus) against Anopheles funestus and Anopheles arabiensis in shaded and unshaded habitats; and explored how larval development might influence retreatment intervals. These tests were done in semi-natural habitats using field-collected larvae, with untreated habitats as controls. Additionally, larval development was assessed in semi-natural and natural habitats in rural Tanzania, by sampling daily and recording larval developmental stages. Both biolarvicides reduced larval densities of both species by >98% within 72 h. Efficacy lasted one week in sun-exposed habitats but remained >50% for two weeks in shaded habitats. An. funestus spent up to two weeks before pupating (13.2(10.4-16.0) days in semi-natural; 10.0(6.6-13.5) in natural habitats), while An. arabiensis required slightly over one week (8.2 (5.8-10.6) days in semi-natural; 8.3 (5.0-11.6) in natural habitats). The findings suggest that weekly larviciding, which is essential for An. arabiensis might be more effective for An. funestus whose prolonged aquatic growth allows for repeated exposures. Additionally, the longer residual effect of biolarvicides in shaded habitats indicates they may require less frequent treatments compared to sun-exposed areas.

Entomological longitudinal surveys in two contrasted eco-climatic settings in Cameroon reveal a high malaria transmission from Anopheles funestus associated with GSTe2 metabolic resistance

BACKGROUND

The impact of metabolic resistance to insecticides on malaria transmission remains poorly characterised notably through application of entomological parameters. The lack of resistance markers has been one of the limiting factors preventing a robust assessment of such impact. To this end, the present study sought to investigate how the L119F-Gste2 metabolic gene influences entomological parameters underpinning mosquitos' propensity to transmit Plasmodium spp.

METHODS

Longitudinal studies were carried out in Mibellon and Elende, two different eco-climatic settings in Cameroon and mosquitoes were collected using Human Landing Catch (HLC), Centre for Disease Control Light Trap (CDC-LT) and Pyrethrum Spray Catch (PSC) technics. Plasmodium sporozoite parasites were detected by TaqMan and Nested PCR, and blood meal origin by ELISA. The allele-specific PCR (AS-PCR) method was used to genotype the L119F-GSTe2 marker and association with malaria transmission was established by comparing key transmission parameters such as the Entomological Inoculation Rate (EIR) between individuals with different L119F-GSTe2 genotypes.

RESULTS

An. funestus s.l was the predominant malaria vector collected during the entomological survey in both sites (86.6% and 96.4% in Elende and Mibellon, respectively) followed by An. gambiae s.l (7.5% and 2.4%, respectively). Sporozoite infection rates were very high in both collection sites (8.7% and 11% in Elende and Mibellon, respectively). An. funestus s.s exhibited a very high entomological inoculation rate (EIR) (66 ib/h/month and 792 ib/h/year) and was responsible for 98.6% of all malaria transmission events occurring in both sites. The Human Blood Index was also high in both locations (HBI = 94%). An. funestus s.s. mosquitoes with both 119 F/F (RR) and L119F (RS) genotypes had a significantly higher transmission intensity than their susceptible L/L119 (SS) counterparts (IRR = 2.2, 95%CI (1.1-5.2), p = 0.03; IRR = 2.5, 95% CI (1.2-5.8), p = 0.01 respectively).

CONCLUSION

This study highlights the major role that An. funestus s.s plays in malaria transmission in Cameroon with an aggravation from GSTe2-based metabolic resistance.

Common predators and factors influencing their abundance in Anopheles funestus aquatic habitats in rural south-eastern Tanzania

BACKGROUND

The role of larval predators in regulating the Anopheles funestus population in various malaria-endemic countries remains relatively unknown. This study aimed to investigate the common predators that co-exist with Anopheles funestus group larvae and evaluate factors that influence their abundance in rural south-eastern Tanzania.

METHODS

Mosquito larvae and predators were sampled concurrently using standard dipper (350 ml) or 10 L bucket in previously identified aquatic habitats in selected villages in southern Tanzania. Predators and mosquito larvae were identified using standard identification keys. All positive habitats were geo-located and their physical features characterized. Water physicochemical parameters such as dissolved oxygen (DO), pH, electrical conductivity (EC), total dissolved solids (TDS) and temperature were also recorded.

RESULTS

A total of 85 previously identified An. funestus aquatic habitats in nine villages were sampled for larvae and potential predators. A total of 8,295 predators were sampled. Of these Coenagrionidae 57.7% (n = 4785), Corixidae 12.8% (n = 1,060), Notonectidae 9.9% (n = 822), Aeshnidae 4.9% (n = 405), Amphibian 4.5% (n = 370), Dytiscidae 3.8% (n = 313) were common. A total of 5,260 mosquito larvae were sampled, whereby Anopheles funestus group were 60.3% (n = 3,170), Culex spp. 24.3% (n = 1,279), An. gambie s.l. 8.3% (n = 438) and other anophelines 7.1% (n = 373). Permanent and aquatic habitats larger than 100m2 were positively associated with An. funestus group larvae (P<0.05) and predator abundance (P<0.05). Habitats with submerged vegetation were negatively associated with An. funestus group larvae (P<0.05). Only dissolved oxygen (DO) was positively and significantly affect the abundance of An. funestus group larvae (P<0.05). While predators' abundance was not impacted by all physicochemical parameters.

CONCLUSION

Six potential predator families were common in aquatic habitats of An. funestus group larvae. Additional studies are needed to demonstrate the efficacy of different predators on larval density and adult fitness traits. Interventions leveraging the interaction between mosquitoes and predators can be established to disrupt the transmission potential and survival of the An. funestus mosquitoes.

CDC light traps underestimate the protective efficacy of an indoor spatial repellent against bites from wild Anopheles arabiensis mosquitoes in Tanzania

BACKGROUND

Methods for evaluating efficacy of core malaria interventions in experimental and operational settings are well established but gaps exist for spatial repellents (SR). The objective of this study was to compare three different techniques: (1) collection of blood-fed mosquitoes (feeding), (2) human landing catch (HLC), and (3) CDC light trap (CDC-LT) collections for measuring the indoor protective efficacy (PE) of the volatile pyrethroid SR product Mosquito Shield^™ METHODS: The PE of Mosquito Shield^™ against a wild population of pyrethroid-resistant Anopheles arabiensis mosquitoes was determined via feeding, HLC, or CDC-LT using four simultaneous 3 by 3 Latin squares (LS) run using 12 experimental huts in Tanzania. On any given night each technique was assigned to two huts with control and two huts with treatment. The LS were run twice over 18 nights to give a sample size of 72 replicates for each technique. Data were analysed by negative binomial regression.

RESULTS

The PE of Mosquito Shield^™ measured as feeding inhibition was 84% (95% confidence interval (CI) 58-94% [Incidence Rate Ratio (IRR) 0.16 (0.06-0.42), p < 0.001]; landing inhibition 77% [64-86%, (IRR 0.23 (0.14-0.36) p < 0.001]; and reduction in numbers collected by CDC-LT 30% (0-56%) [IRR 0.70 (0.44-1.0) p = 0.160]. Analysis of the agreement of the PE measured by each technique relative to HLC indicated no statistical difference in PE measured by feeding inhibition and landing inhibition [IRR 0.73 (0.25-2.12) p = 0.568], but a significant difference in PE measured by CDC-LT and landing inhibition [IRR 3.13 (1.57-6.26) p = 0.001].

CONCLUSION

HLC gave a similar estimate of PE of Mosquito Shield^™ against An. arabiensis mosquitoes when compared to measuring blood-feeding directly, while CDC-LT underestimated PE relative to the other techniques. The results of this study indicate that CDC-LT could not effectively estimate PE of the indoor spatial repellent in this setting. It is critical to first evaluate the use of CDC-LT (and other tools) in local settings prior to their use in entomological studies when evaluating the impact of indoor SR to ensure that they reflect the true PE of the intervention.

Relationships between biological age, distance from aquatic habitats and pyrethroid resistance status of Anopheles funestus mosquitoes in south-eastern Tanzania

BACKGROUND

Malaria transmission can be highly heterogeneous between and within localities, and is influenced by factors such as survival and biting frequencies of Anopheles mosquitoes. This study investigated the relationships between the biological age, distance from aquatic habitats and pyrethroid resistance status of Anopheles funestus mosquitoes, which currently dominate malaria transmission in south-east Tanzania. The study also examined how such relationships may influence malaria transmission and control.

METHODS

Female An. funestus were collected in houses located 50-100 m, 150-200 m or over 200 m from the nearest known aquatic habitats. The mosquitoes were exposed to 1×, 5× and 10× the diagnostic doses of deltamethrin or permethrin, or to the synergist, piperonyl butoxide (PBO) followed by the pyrethroids, then monitored for 24 h-mortality. Ovaries of exposed and non-exposed mosquitoes were dissected to assess parity as a proxy for biological age. Adults emerging from larval collections in the same villages were tested against the same insecticides at 3-5, 8-11 or 17-20 days old.

FINDINGS

Mosquitoes collected nearest to the aquatic habitats (50-100 m) had the lowest mortalities compared to other distances, with a maximum of 51% mortality at 10× permethrin. For the age-synchronized mosquitoes collected as larvae, the insecticide-induced mortality assessed at both the diagnostic and multiplicative doses (1×, 5× and 10×) increased with mosquito age. The highest mortalities at 1× doses were observed among the oldest mosquitoes (17-20 days). At 10× doses, mortalities were 99% (permethrin) and 76% (deltamethrin) among 8-11 day-olds compared to 80% (permethrin) and 58% (deltamethrin) among 3-5 day-olds. Pre-exposure to PBO increased the potency of both pyrethroids. The proportion of parous females was highest among mosquitoes collected farthest from the habitats.

CONCLUSION

In this specific setting, older An. funestus and those collected farthest from the aquatic habitats (near the centre of the village) were more susceptible to pyrethroids than the younger ones and those caught nearest to the habitats. These findings suggest that pyrethroid-based interventions may remain at least moderately effective despite widespread pyrethroid-resistance, by killing the older, less-resistant and potentially-infective mosquitoes. Further studies should investigate how and whether these observations could be exploited to optimize malaria control in different settings.