Skip to main content

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



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.


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.


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.


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.


The risk of malaria transmission in endemic communities is influenced by multiple environmental and biological factors, often leading to a heterogeneous pattern, which can be further influenced by intervention [1,2,3,4]. These factors include the distribution of mosquito breeding habitats and human dwellings [5,6,7], as well as mosquito-related factors such as blood-feeding behaviours, dispersal range and non-random host selection [8,9,10,11]. When human settlements are close to the breeding sites, the range at which mosquitoes can disperse is limited [12], and vector densities tend to be highest where the settlements are most concentrated [11]. These variations can occur at all levels of malaria endemicity, and at different geographical scales. However, they tend to be greater in low-transmission compared to high-transmission areas and can be significant over fine geographic scales such as between or within villages [1].

The landscape-level variations, especially as influenced by non-random mosquito-biting frequencies, have been elucidated broadly both in situ and in silico, and are known to have significant epidemiological implications [1, 9, 13]. Understanding these variations of malaria transmission risk can enable more strategic approaches and better resource allocation for malaria control, even though the actual impact of such targeting remains unclear and poorly understood [14]. It is well known that female Anopheles mosquitoes can travel for varied distances [12, 15], but final aggregations of these vectors tend to be highest where human densities and household occupancy are greatest [11]. Also, while mosquito densities tend to be highest near their aquatic habitats, it is the females found far from their larval habitations that tend to be older and more likely to be infected with malaria parasites [7]. Furthermore, the infectivity of wild-caught mosquitoes depends on their biological age and the density of the parasite at ingestion [16]. This implies differential infectivity since malaria vectors generally require more than ten days to ensure complete parasite development and maturation [17].

Beyond the environmental factors and their influence on mosquito dispersal, experimental observations have also shown that the susceptibility of mosquitoes to insecticides can increase with age; and that older mosquitoes are generally more likely to be killed than younger ones [18,19,20]. Moreover, resistance in major malaria vectors can vary seasonally and spatially, in some cases at fine geographical scales [21, 22]. The observations may also have direct implications on the effectiveness of the insecticidal interventions commonly used against malaria vectors, namely insecticide-treated nets (ITNs) and indoor residual spraying (IRS). Unfortunately, there is currently limited research on how these interactions play out in communities such as rural south-eastern Tanzania, where ITNs have been widely used for more than two decades; and malaria transmission risk has become highly fragmented on a spatial scale, with greater than 10 fold differences in transmission intensity over short distances less than 50 km [23,24,25,26,27]. In these areas, the dominant vector species is An. funestus, which mediates over 80% of all malaria transmission events [25, 26]. In particular, there has not been any detailed analysis that integrates observations of the environmental factors (e.g. distance from habitats) and biological factors (e.g. mosquito infectivity and dispersal rates) with intervention-related factors (e.g. insecticide resistance levels in the dominant malaria vectors); and how such correlations may influence disease transmission and control.

This study, therefore, assessed, at fine-scale, the interplay between key biological and environmental factors that influence malaria transmission by An. funestus in rural south-eastern Tanzania. The key factors investigated were: (a) age of An. funestus mosquitoes, (b) the susceptibility of these mosquitoes to pyrethroid insecticides, and (c) the spatial distances from nearest aquatic habitats.


Study site

Adult female An. funestus mosquitoes were collected from the spatially isolated hamlet in Ikwambi village in Kilombero district, south-eastern Tanzania (location: 7.98033°S, 36.81701°E; altitude: 400 m above the sea level; annual rainfall: 1200–1800 mm; temperatures: 20–33 °C) [28] (Fig. 1). Anopheles funestus dominates this village in both densities and malaria transmission activity. The collections were done in the dry months from August to early December 2020.

Fig. 1
figure 1

Map of the study area as a representative for adult mosquito collections at different distances from aquatic habitats

Mosquito collections

Initial surveys were undertaken to identify and mark all possible aquatic habitats of An. funestus, based on attributes and methods previously described by Nambunga et al. [29]. The habitats generally had clear waters, were formed by the river stream, had emergent vegetation and remained with water even during the dry season. The distance from confirmed habitats of An. funestus to human habitations were calculated based on geo-locations data obtained using a handled GPS receiver (eTrex, Vista, Garmin, USA). Three distance ranges were defined from the aquatic habitat as follows; (a) 50–100 m (edge of the village), (b) 150–200 m (intermediate distance), and (c) over 200 m (at the centre of the village) (Fig. 1).

Centre for Disease Control and Prevention (CDC) light traps placed beside volunteer-occupied bed nets were used to collect host-seeking female mosquitoes indoors from 18:00 to 07:00 h each night. The collections were grouped by distance and acclimatized for 24 h in an insectary where they were supplied with 10% glucose, and only the survivors were used for subsequent tests. Additional collections were done using volunteer-occupied miniaturized double net (DN-Mini) traps placed near human dwellings [30]. The identification of the collected mosquitoes was done using dichotomous keys of African Anopheles [31].

Besides the adult collections, An. funestus larvae were also collected using procedures described by Nambunga et al., [29] and transported in their natural water for onward rearing in the insectary (at 80 ± 10% relative humidity and 27 ± 2 °C) at the Ifakara Health Institute vector biology laboratory, the VectorSphere. Tetramin® fish food (Tetra GmbH, Melle, Germany) was provided to supplement the aquatic diet, and the pupae collected were placed in separate cages supplied with 10% glucose. The emergent, mosquitoes were used to assess variations of susceptibility in three age groups, as follows; (i) 3–5 days, (ii) 8–11 days and (iii) 17–20 days.

Insecticide resistance bioassays

Two pyrethroids (deltamethrin and permethrin) were evaluated using WHO guidelines [32], starting with standard diagnostic insecticide doses (1×) for phenotypic resistance evaluation, followed by 5× and 10× the diagnostic doses to assess intensities of resistance. Since malaria vectors in the area are resistant to pyrethroids, additional tests were done using the synergist, piperonyl butoxide (PBO), to investigate the role of cytochrome P450 enzymes, and the potency of PBO in reversing resistance [21, 22].

Each assay consisted of six replicates with 20 individual mosquitoes per replicate, totalling 120 mosquitoes per bioassay. In the first round of bioassays (using adults collected from houses), the tests were done separately for each distance range from aquatic habitats (50–100 m, 150–200 m and over 200 m); and in the second round (using age-synchronized adults collected as larvae) the tests were done separately for the different age classes (3–5, 8–11 and 17–20 days). Mosquitoes were exposed to the insecticides for 60 min and moved to non-insecticidal holding tubes with 10% glucose, then monitored for 24 h post-exposure mortalities. Since the mosquitoes were resistant to the diagnostic concentrations of both pyrethroids, the additional tests for resistance intensity were completed as prescribed by WHO [32].

For synergy tests with PBO four cohorts of mosquitoes were exposed to: (i) deltamethrin or permethrin only at diagnostic concentration, (ii) 4% PBO followed by diagnostic dose of deltamethrin or permethrin, (iii) 4% PBO only, or (iv) silicone oil coated papers (control group). Each test was repeated three times, and the mosquitoes were provided 10% glucose during the holding period. Mortality was recorded 24 h post-exposure.

Dissections to assess parity status and estimate the number of gonotrophic cycles

Following insecticide bioassays, live and dead mosquitoes 24 h post-exposure, were immediately dissected (n = 1577). The dissections were done under the stereo-microscope and the dissected ovaries were observed under a compound microscope for the presence or absence of the coiled tracheolar skeins indicating nulliparous or parous status respectively [33, 34]. Since studies on landscape distribution of the biological age of mosquitoes can be affected by the methods used for mosquito collection, an additional cohort of mosquitoes collected using the CDC light traps (n = 560) and the DN-Mini trap (n = 78), and not tested for resistance were dissected for more detailed identification of the ovariole dilations to determine how many times the mosquito had laid eggs (i.e. the number of gonotrophic cycles) [34, 35]. The results of the dissections were recorded by distance from the aquatic habitat(s) and method of collection.

Molecular identification of the mosquitoes

A sub-sample of the tested mosquitoes (at least 10% from each replicate, n = 409) was packed in micro-centrifuge tubes containing silica desiccant. Sibling species of An. funestus were identified by PCR using nucleic acid material extracted from the legs, to screen for species-specific nucleotide sequences (internal transcribed spacer 2) in the ribosomal DNA (rDNA) [36].

Data analysis

The percentage mortality of mosquitoes was calculated as a fraction of the total number exposed and interpreted according to the WHO guidelines [32]. Since no control mortalities exceeded 5%, observed but not corrected mortalities are reported. The mosquitoes were considered susceptible if mortality was \(\ge\)98%, resistant, if mortality was < 90% or possibly resistant and requiring additional tests if mortality ranged between 90% and 97% [32]. The mean proportions of parous mosquitoes at different distances were compared using the analysis of variance (ANOVA). To determine which specific distance pairs were significantly different, Tukey’s post hoc test was applied.

Additional analysis was done using statistical software, R-Software version 3.6.0 [37]. A generalized linear mixed effect model (GLMM) with binomial distribution and logit link function was used to assess the proportion of mosquitoes that died 24 h post-insecticide exposure at different age groups and distances. Age and distance were added as fixed factors, whereas replicate was added as a random factor. The nearest distance to the habitats (50–100 m) and the youngest age class of mosquitoes (3–5 days) were considered as references for these analyses. Odd ratios with their corresponding 95% confidence intervals were reported. Statistical significance was considered when the p-value is less than 0.05.


Relationship between insecticide susceptibility and distance from aquatic habitats

Mosquitoes collected in different houses were resistant to both deltamethrin and permethrin at the diagnostic doses (1×) but the percentage of mortality varied by distance from the aquatic habitats (Fig. 2). The mosquitoes collected farthest from the habitats showed significantly higher mortality than those collected nearest to the habitats (Table 1). A similar trend was observed when the mosquitoes were exposed to 5× and 10× the diagnostic doses of permethrin but not deltamethrin, with mosquitoes from the nearest distances showing the lowest mortality and those at the farthest distance showing the highest mortality. Overall, there was lower mortality in tests against deltamethrin compared to permethrin (Fig. 2). Further analysis showed that the association of 24-h mortality with distance was statistically significant for deltamethrin only at the diagnostic dose, and for permethrin at all doses (Table 1).

Fig. 2
figure 2

Percentage mortality in An. funestus mosquitoes, collected at different distances from the aquatic habitats and exposed to deltamethrin or permethrin. The red-dotted lines represent 90% mortality and the blue-dotted lines represent 98% mortality

Table 1 Summary of 24-hour mortality of Anopheles funestus mosquitoes in tests of adults collected at different distances from aquatic habitats and exposed to deltamethrin or permethrin with or without the synergist, piperonyl butoxide (PBO)

Pre-exposure to the synergist, PBO, significantly increased the potency of the candidate pyrethroids against the resistant mosquitoes, achieving more than 80% mortality at all distances when tested using the diagnostic doses (Fig. 3). However, the association of 24-h mortality with distance in these PBO tests was statistically significant only for deltamethrin (Table 1).

Fig. 3
figure 3

Percentage mortality in An. funestus mosquitoes, collected at different distances from the aquatic habitats and exposed to either the candidate pyrethroids alone or PBO followed by pyrethroids. The red-dotted lines represent 90% mortality and the blue-dotted lines represent 98% mortality

Relationship between biological age of mosquitoes and distance from aquatic habitats

A total of 1397 female An. funestus were dissected from the first cohort. Although the overall difference was statistically significant (F = 4.38, p = 0.015), only the mosquitoes collected by CDC light traps farthest from aquatic habitats (over 200 m) had significantly greater proportions of parous mosquitoes compared to those collected nearest to the habitats (50–100 m) (p = 0.017, Fig. 4). Some cohorts of mosquitoes collected beyond 200 m were 100% parous. When the additional cohort of mosquitoes collected by the light traps was assessed for the number of gonotrophic cycles (n = 560), the majority were found to be nulliparous; and those with multiple gonotrophic cycles were equally distributed over the distances (50–100 m, 150–200 m, over 200 m) (Fig. 5). On the other hand, among the mosquitoes collected by DN-Mini (n = 78), the proportion that was nulliparous dropped while the proportion with multiple gonotrophic cycles was slightly higher at 150–200 m and \(>\)200 m distances compared to 50–100 m (Fig. 5).

Fig. 4
figure 4

The parity percentages of An. funestus mosquitoes, collected at different distances and pre-exposed to insecticide. Red dots indicate the means

Fig. 5
figure 5

Proportions of mosquitoes with different numbers of gonotrophic cycles, among those collected by either CDC light traps (top panel) or the miniaturized double net (DN-Mini) trap (bottom panel) at different distances

Relationship between insecticide susceptibility and chronological age of mosquitoes

For both pyrethroids, at the diagnostic dose (1×), the percentage of mortality increased with age, suggesting that the older An. funestus were more susceptible than younger ones from the same locations (Fig. 6, Table 2). Similarly, in the intensity bioassays conducted with 5× or 10× the diagnostic doses, the levels of resistance varied by age, but mortality significantly increased, reaching 99% for permethrin and 81% for deltamethrin (Fig. 6). Further analysis showed that the association of 24-hour mortality with mosquito age was statistically significant for both deltamethrin and permethrin, and at all doses (Table 2).

Fig. 6
figure 6

Percentage mortality in different age cohorts of emergent An. funestus mosquitoes, collected as larvae and exposed to deltamethrin or permethrin. The red-dotted line represents 90% mortality, the blue-dotted line represents 98% mortality

Table 2 Summary of 24-hour mortality of Anopheles funestus mosquitoes collected as larvae and exposed to pyrethroids (deltamethrin or permethrin) at different age classes (3–5, 8–11 or 17–20 days old)

Molecular identification of the mosquitoes

A total of 409 An. funestus mosquitoes were tested by PCR to identify sibling species in the An. funestus group, including 60% (n = 245) from mosquitoes collected as adults; and 40% (n = 164) of the adults emerging from larval collections. The amplification rate for the adult collected mosquitoes was 95% (n = 123), and all of them (100%) were identified as An. funestus sensu stricto (s.s.). On the other hand, the amplification rate for those collected as larvae was 90% (n = 148), of which 99.3% (n = 147) were An. funestus s.s, the remaining one specimen being Anopheles rivulorum.


To improve surveillance and control of malaria in endemic Africa, there is need for greater understanding of both landscape and local variations in the risk of transmission [7, 38]. Data on population distribution, age-structure and susceptibility of malaria vectors are particularly crucial for analysing the local risk patterns, designing fine-scale targets of interventions and selection of appropriate interventions. This study demonstrated a positive correlation between the mortality of An. funestus females exposed to the pyrethroid insecticides and the distance between the sites where these mosquitoes were collected and the nearest aquatic habitats. A positive correlation was observed in tests with the diagnostic dose of insecticides as well as in tests with the five and ten-fold increases in the diagnostic dose. The highest mortalities were obtained among mosquitoes farthest from the aquatic habitats, nearest to the centre of the village, and among the oldest mosquito cohorts. This suggests that insecticidal interventions may have differential impact dependent on distance from the main mosquito habitats.

Mathematical predictions by Smith et al. initially showed that human-biting densities would be greatest near breeding sites where adult mosquitoes emerge, even as the proportions old enough to transmit malaria would be found far from the habitats [7]. An important question, therefore, was how the insecticide resistance profiles influence this system, and whether these potentially infectious adults would also be the most susceptible as suggested by the first experiments in this study. This study has provided a unique opportunity for a concurrent analysis of the age distribution and insecticide susceptibility of malaria vectors over the same distances, thus enabling improved understanding of these variations and their potential epidemiological implications. The proportion of parous females was greatest among the mosquitoes collected farthest from the aquatic habitats, compared to those collected nearest to the aquatic habitats. These findings validate previous studies using the mark-release and recapture techniques, where older female Anopheles mosquitoes were collected at the farthest distance, over 600 m from a point of release and survived over 17 days in the wild [15]. Such non-random age distribution of mosquito populations could have a direct impact on the disease transmission risk [6, 7, 9].

Less than 30% of the mosquitoes caught by the CDC had two or more gonotrophic cycles, and the maximum number of cycles was three. The low proportions of old mosquitoes in light traps may reflect the tendency of CDC light traps to catch mostly nulliparous mosquitos in such settings [39, 40]. Secondly it may be due to the fact that only live or freshly dead mosquitoes were dissected, yet parous mosquitoes can be more likely to die during collection [39]. Interestingly, when the DN-Mini trap was used, the proportion of nulliparous mosquitoes decreased (relative to the CDC light trap collections), and there were significant proportions with at least one or two gonotrophic cycles; and a slight increase in proportions with multiple gonotrophic cycles at distances far from the habitats. A related lesson here is that trapping methods may yield biased catches requiring careful interpretation in a broader context [30, 41].

In the tests using mosquitoes collected as larvae from the field sites, the age-synchronized females also showed positive correlations between the chronological age classes and the insecticide potency. The youngest An. funestus (3–5 days old) were the most strongly resistant to the candidate insecticides, while the oldest mosquitoes (17–20 days old) were the least resistant. These findings suggest that despite widespread pyrethroid resistance, older mosquitoes may remain substantially susceptible and that insecticidal interventions may retain some degree of efficacy in the communities. The findings by Collins et al., using the wild population, the resistant sample were more likely to be nulliparous compared to the susceptible samples [42]. This is important since it is the older mosquitoes that are responsible for parasite transmission as mosquitoes need approximately 2 weeks for parasite maturity in their bodies [7, 17]. Future research should, therefore, investigate whether this demographic of mosquitoes can or should be selectively targeted and how such a strategy could impact malaria control. Several studies have indeed demonstrated age-related variations in insecticide resistance, in mosquito species other than An. funestus [18,19,20, 43,44,45], but practical applications of this knowledge remain limited. This might be explained by the fact that insecticide metabolism by mosquitoes decreases with age due to reduced active enzymes [18, 46,47,48]. Other factors like blood-feeding and larval nutrition may also play a role in insecticide-induced mortality [45, 49].

The tests where mosquitoes were pre-exposed to PBO showed potential significant involvement of mixed-function oxidases in the resistance. Overall, the potency of pyrethroids and mortalities significantly increased after the PBO exposures, reaching between 82% and 98%. The lack of full recovery in some mosquitoes may suggest that additional mechanisms may be partially involved in a resistance profile of the local An. funestus mosquitoes. It is reasonable to assume that since the pyrethroid resistance in these An. funestus populations were broadly metabolic, the associated detoxifying enzymes may also be differentially expressed by age. However, when all the PBO test data were analysed together, there was no significant correlation between the distances (a proxy for population age structure) and the performance of PBO. These differences have been observed in previous studies and the involvement of the enzymes in these age-related differences remain poorly understood. For example, Christian et al., who analysed mosquitoes aged 3 to 30 days, concluded that two resistance-related genes, CYP6P9 and CYP6P13 were most highly expressed at 10 days of age, but that there was generally no correlation between the expression of these genes and mosquito susceptibility to permethrin [43]. On the other hand, Hazelton et al., who measured the activity of glutathione-s-transferases in mosquitoes during growth, maturity and senescence phases concluded that the enzyme activities generally declined with age [48]. Taken together, the results of this study and the previous studies indicate the need to further investigate the differential expression of resistance genes in different mosquito species and in different settings.

Despite the broad success of this study in analysing the correlations between biological ages, distance from aquatic habitats and pyrethroid resistance status of An. funestus mosquitoes, there were a few limitations. Firstly, as previously described by Charlwood et al., the older mosquitoes are more likely to die during collection when the CDC-LT is used during specimen collection [39]. This might have marginally confounded the age evaluation of both insecticide pre-exposed and non-exposed cohorts. However, age synchronised, larval collected mosquitoes showed the same trend as the mortality of mosquitoes increased with the increased age and concentration. Thus, these results could still be considered informative. Secondly, fewer mosquitoes from the farthest distance (over 200 m) were dissected compared to the other distances as most of them were susceptible to the insecticide, and were therefore dead and too dry to be dissected. The effect of insecticide pre-exposure on dissection was rectified by the additional collections done using both CDC light traps and DN-Mini. Lastly, due to the difficulties of collecting An. funestus at the larval stage, the PBO pre-exposure tests were done only on the adult collections but not on the age-synchronized mosquitoes collected as larvae.


In south-eastern Tanzania, where An. funestus predominates malaria transmission, the older mosquitoes are more likely to be found farthest from the aquatic habitats and near the centre of the villages, and are also less resistant to pyrethroid insecticides. Similarly, in the age synchronized mosquitoes, the oldest cohorts were the least resistant to the candidate pyrethroids. These findings provide new insights into how the insecticide-based interventions could still provide modest levels of protection against malaria in the communities irrespective of the strong insecticide resistance currently reported. Additionally, the findings also imply that greater impact may be achieved if vector control methods were spatially targeted at a fine-scale, in ways that exploit the observed variations in mosquito age and susceptibility to insecticides. While the actual configuration of such a strategy remains to be determined, it could enable a more effective onslaught against the older populations of malaria vectors, which drive pathogen transmission and are more readily controlled by the insecticides.

Availability of data and materials

Data are available under the Ifakara Health Institute data sharing policy upon request.



internal transcribed spacer 2 region


Ifakara Health Institute


Institutional Review Board


Long-lasting insecticidal nets


National Institute for Medical Research


Centre for disease control and prevention


Piperonyl butoxide


Indoor residual spraying


Insecticide-treated bed nets


World Health Organizations


Polymerase chain reaction


Deoxyribonucleic acid


Analysis of variance


  1. Greenwood BM. The microepidemiology of malaria and its importance to malaria control. Trans R Soc Trop Med Hyg. 1989;83(Suppl 1):25–9.

    Article  PubMed  Google Scholar 

  2. Hagenlocher M, Castro MC. Mapping malaria risk and vulnerability in the United Republic of Tanzania: a spatial explicit model. Popul Health Metr. 2015;13:2.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Homan T, Maire N, Hiscox A, Di Pasquale A, Kiche I, Onoka K, et al. Spatially variable risk factors for malaria in a geographically heterogeneous landscape, western Kenya: an explorative study. Malar J. 2016;15:1.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Sadoine ML, Smargiassi A, Ridde V, Tusting LS, Zinszer K. The associations between malaria, interventions, and the environment: a systematic review and meta-analysis. Malar J. 2018;17:73.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Carter R, Mendis KN, Roberts D. Spatial targeting of interventions against malaria. Bull World Health Organ. 2000;78:1401–11.

    CAS  PubMed  PubMed Central  Google Scholar 

  6. Thomas CJ, Lindsay SW. Local-scale variation in malaria infection amongst rural gambian children estimated by satellite remote sensing. Trans R Soc Trop Med Hyg. 2000;94:159–63.

    Article  CAS  PubMed  Google Scholar 

  7. Smith DL, Dushoff J, McKenzie FE. The risk of a mosquito-borne infection in a heterogeneous environment. PLoS Biol. 2004;2:e368.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Takken W, Verhulst NO. Host preferences of blood-feeding Mosquitoes. Annu Rev Entomol. 2013;58:433–53.

    Article  CAS  PubMed  Google Scholar 

  9. Burkot TR. Non-random host selection by anopheline mosquitoes. Parasitol Today. 1988;4:156–62.

    Article  CAS  PubMed  Google Scholar 

  10. Guelbéogo WM, Gonçalves BP, Grignard L, Bradley J, Serme SS, Hellewell J, et al. Variation in natural exposure to Anopheles mosquitoes and its effects on malaria transmission. Elife. 2018;7:e32625.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Kaindoa EW, Mkandawile G, Ligamba G, Kelly-Hope LA, Okumu FO. Correlations between household occupancy and malaria vector biting risk in rural Tanzanian villages: implications for high-resolution spatial targeting of control interventions. Malar J. 2016;15:199.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Takken W, Charlwood JD, Billingsley PF, Gort G. Dispersal and survival of Anopheles funestus and A. gambiae s.l (Diptera: Culicidae) during the rainy season in southeast Tanzania. Bull Entomol Res. 1998;88:561–6.

    Article  Google Scholar 

  13. Woolhouse MEJ, Dye C, Etard JF, Smith T, Charlwood JD, Garnett GP, et al. Heterogeneities in the transmission of infectious agents: implications for the design of control programs. Proc Natl Acad Sci USA. 1997;94:338–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Stresman G, Bousema T, Cook J. Malaria hotspots: is there epidemiological evidence for fine-scale spatial targeting of interventions? Trends Parasitol. 2019;35:822–34.

    Article  PubMed  Google Scholar 

  15. Midega JT, Mbogo CM, Mwambi H, Wilson MD, Ojwang G, Mwangangi JM, et al. Estimating dispersal and survival of Anopheles gambiae and Anopheles funestus along the kenyan coast by using Mark–Release–Recapture methods. J Med Entomol. 2007;44:923–9.

    Article  PubMed  Google Scholar 

  16. Dawes EJ, Churcher TS, Zhuang S, Sinden RE, Basáñez M-G. Anopheles mortality is both age- and plasmodium-density dependent: implications for malaria transmission. Malar J. 2009;8:228.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Ohm JR, Baldini F, Barreaux P, Lefevre T, Lynch PA, Suh E, et al. Rethinking the extrinsic incubation period of malaria parasites. Parasit Vectors. 2018;11:178.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Rajatileka S, Burhani J, Ranson H. Mosquito age and susceptibility to insecticides. Trans R Soc Trop Med Hyg. 2011;105:247–53.

    Article  CAS  PubMed  Google Scholar 

  19. Lines JD, Nassor NS. DDT resistance in Anopheles gambiae declines with mosquito age. Med Vet Entomol. 1991;5:261–5.

    Article  CAS  PubMed  Google Scholar 

  20. Chouaibou MS, Chabi J, Bingham GV, Knox TB, N’Dri L, Kesse NB, et al. Increase in susceptibility to insecticides with aging of wild Anopheles gambiae mosquitoes from Côte d’Ivoire. BMC Infect Dis. 2012;12:214.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Matowo NS, Munhenga G, Tanner M, Coetzee M, Feringa WF, Ngowo HS, et al. Fine-scale spatial and temporal heterogeneities in insecticide resistance profiles of the malaria vector, Anopheles arabiensis in rural south-eastern Tanzania. Wellcome Open Res. 2017;2:96.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Pinda PG, Eichenberger C, Ngowo HS, Msaky DS, Abbasi S, Kihonda J, et al. Comparative assessment of insecticide resistance phenotypes in two major malaria vectors, Anopheles funestus and Anopheles arabiensis in south-eastern Tanzania. Malar J. 2020;19:408.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Finda MF, Limwagu AJ, Ngowo HS, Matowo NS, Swai JK, Kaindoa E, et al. Dramatic decreases of malaria transmission intensities in Ifakara, south-eastern Tanzania since early 2000s. Malar J. 2018;17:362.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Swai JK, Mmbando AS, Ngowo HS, Odufuwa OG, Finda MF, Mponzi W, et al. Protecting migratory farmers in rural Tanzania using eave ribbons treated with the spatial mosquito repellent, transfluthrin. Malar J. 2019;18:414.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Kaindoa EW, Matowo NS, Ngowo HS, Mkandawile G, Mmbando A, Finda M, et al. Interventions that effectively target Anopheles funestus mosquitoes could significantly improve control of persistent malaria transmission in south – eastern Tanzania. PLoS ONE. 2017;12:e0177807.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Mapua SA, Hape EE, Kihonda J, Bwanary H, Kifungo K, Kilalangongono M, et al. Persistently high proportions of Plasmodium-infected Anopheles funestus mosquitoes in two villages in the Kilombero valley, South-Eastern Tanzania. Parasite Epidemiol Control. 2022;18:e00264.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Okumu F, Finda M. Key characteristics of residual malaria transmission in two districts in South-Eastern Tanzania—Implications for Improved Control. J Infect Dis. 2021;223:143–54.

    Article  Google Scholar 

  28. Spark W. Average Weather in Ifakara, Tanzania, Year Round - Weather Spark.

  29. Nambunga IH, Ngowo HS, Mapua SA, Hape EE, Msugupakulya BJ, Msaky DS, et al. Aquatic habitats of the malaria vector Anopheles funestus in rural south-eastern Tanzania. Malar J. 2020;19:219.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Limwagu AJ, Kaindoa EW, Ngowo HS, Hape E, Finda M, Mkandawile G, et al. Using a miniaturized double-net trap (DN-Mini) to assess relationships between indoor–outdoor biting preferences and physiological ages of two malaria vectors, Anopheles arabiensis and Anopheles funestus. Malar J. 2019;18:282.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Coetzee M. Key to the females of Afrotropical Anopheles mosquitoes (Diptera: Culicidae). Malar J. 2020;19:70.

    Article  PubMed  PubMed Central  Google Scholar 

  32. WHO. Test procedures for insecticide resistance monitoring in malaria vector mosquitoes 2nd Edn. Geneva, World Health Organization; 2018.

  33. Detinova T, Sergeevna DS, Bertram, WHO. Age-grouping methods in Diptera of medical importance, with special reference to some vectors of malaria. Geneva: World Health Organization; 1962.

    Book  Google Scholar 

  34. Silver JB. Methods of age-grading adults and estimation of adult Survival Rates. In: Mosquito ecology field sampling methods. New York, Springer; 2008. pp. 1161–271.

  35. Polovodova V. [The determination of the physiological age of female Anopheles by the number of gonotrophic cycles completed] (in rusian). Med Parasitol. 1949;28:352–5.

    Google Scholar 

  36. Koekemoer LL, Kamau L, Hunt RH, Coetzee M. A cocktail polymerase chain reaction assay to identify members of the Anopheles funestus (Diptera: Culicidae) group. Am J Trop Med Hyg. 2002;66:804–11.

    Article  CAS  PubMed  Google Scholar 

  37. R Core Team. A language and environment for statistical computing. R Found. Stat. Comput. Vienna, Austria. http// 2019.

  38. Cotter C, Sturrock HJW, Hsiang MS, Liu J, Phillips AA, Hwang J, et al. The changing epidemiology of malaria elimination: new strategies for new challenges. Lancet. 2013;382:900–11.

    Article  PubMed  Google Scholar 

  39. Charlwood JD, Tomás EVE, Andegiorgish AK, Mihreteab S, LeClair C. “We like it wet”: a comparison between dissection techniques for the assessment of parity in Anopheles arabiensis and determination of sac stage in mosquitoes alive or dead on collection. PeerJ. 2018;2018:e5155.

    Article  Google Scholar 

  40. Drakeley C, Schellenberg D, Kihonda J, Sousa CA, Arez AP, Lopes D, et al. An estimation of the entomological inoculation rate for Ifakara: a semi-urban area in a region of intense malaria transmission in Tanzania. Trop Med Int Heal. 2003;8:767–74.

    Article  CAS  Google Scholar 

  41. Mboera LE. Sampling techniques for adult afrotropical malaria vectors and their reliability in the estimation of entomological inoculation rate. Tanzan Health Res Bull. 2005;7:117–24.

    CAS  PubMed  Google Scholar 

  42. Collins E, Vaselli NM, Sylla M, Beavogui AH, Orsborne J, Lawrence G, et al. The relationship between insecticide resistance, mosquito age and malaria prevalence in Anopheles gambiae s.l from Guinea. Sci Rep. 2019;9:8846.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Christian R, Matambo T, Spillings B, Brooke B, Coetzee M, Koekemoer L, et al. Age-related pyrethroid resistance is not a function of P450 gene expression in the major african malaria vector, Anopheles funestus (Diptera: Culicidae). Genet Mol Res. 2011;10:3220–9.

    Article  CAS  PubMed  Google Scholar 

  44. Mbepera S, Nkwengulila G, Peter R, Mausa EA, Mahande AM, Coetzee M, et al. The influence of age on insecticide susceptibility of Anopheles arabiensis during dry and rainy seasons in rice irrigation schemes of Northern Tanzania. Malar J. 2017;16:364.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Kulma K, Saddler A, Koella JC. Effects of Age and Larval Nutrition on phenotypic expression of insecticide-resistance in Anopheles mosquitoes. PLoS ONE. 2013;8:e58322.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Rowland M, Hemingway J. Changes in malathion resistance with age in Anopheles stephensi from Pakistan. Pestic Biochem Physiol. 1987;28:239–47.

    Article  CAS  Google Scholar 

  47. Hunt RH, Brooke BD, Pillay C, Koekemoer LL, Coetzee M. Laboratory selection for and characteristics of pyrethroid resistance in the malaria vector Anopheles funestus. Med Vet Entomol. 2005;19:271–5.

    Article  CAS  PubMed  Google Scholar 

  48. Hazelton GA, Lang CA. Glutathione S-transferase activities in the yellow-fever mosquito aedes aegypti (Louisville) during growth and aging. Biochem J. 1983;210:281–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Oliver SV, Brooke BD. The effect of multiple blood-feeding on the longevity and insecticide resistant phenotype in the major malaria vector Anopheles arabiensis (Diptera: Culicidae). Parasit Vectors. 2014;7:390.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


We thank the Ifakara Health Institute research team for their great support and for Community members who volunteered their houses and environment for mosquito collections.


The activities in this work were supported by Anopheles funestus 100, a Grant from the Bill & Melinda Gates Foundation (INV-002138). PP and FOO were supported by Howard Hughes International Scholarship (OPP1099295). LLK was supported by a DST/NRF South African Research Chairs Initiative Grant (UID 64763).

Author information

Authors and Affiliations



PGP, EWK, FOO and LLK conceptualized the idea, designed the study and supervised the experiments; PGP analysed the data and drafted the manuscript; DSM, LLM, IHM, RMN and HSN helped in the data analysis and drafting of the manuscript; JK and HB helped with the field mosquito collection and morphological identification. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Polius G. Pinda or Fredros O. Okumu.

Ethics declarations

Ethics approval and consent to participate

This study was permitted by the Institute Review Board of Ifakara Health Institute; IHI/IRB/No: 26-2020 and the Medical Research Coordinated Committee of the National Institute for Medical Research of the United Republic of Tanzania; NIMR/HQ/R.8a/Vol. IX/3495. With these approvals, an ethical approval waiver was obtained from the Animal Research Ethics Committee of the University of the Witwatersrand, where the lead author is a post-graduate student (Waiver 11-10-2021-O).

Consent for publication

This manuscript has been approved for publication by the National Institute for Medical Research of the United Republic of Tanzania (NIMR/HQ/P.12 VOL XXXIV/78).

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Pinda, P.G., Msaky, D.S., Muyaga, L.L. et al. Relationships between biological age, distance from aquatic habitats and pyrethroid resistance status of Anopheles funestus mosquitoes in south-eastern Tanzania. Malar J 21, 365 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: