Skip to main content

Trends of insecticide resistance monitoring in mainland Tanzania, 2004–2020

Abstract

Background

Insecticide resistance is a serious threat to the continued effectiveness of insecticide-based malaria vector control measures, such as long-lasting insecticidal nets (LLINs) and indoor residual spraying (IRS). This paper describes trends and dynamics of insecticide resistance and its underlying mechanisms from annual resistance monitoring surveys on Anopheles gambiae sensu lato (s.l.) populations conducted across mainland Tanzania from 2004 to 2020.

Methods

The World Health Organization (WHO) standard protocols were used to assess susceptibility of the wild female An. gambiae s.l. mosquitoes to insecticides, with mosquitoes exposed to diagnostic concentrations of permethrin, deltamethrin, lambdacyhalothrin, bendiocarb, and pirimiphos-methyl. WHO test papers at 5× and 10× the diagnostic concentrations were used to assess the intensity of resistance to pyrethroids; synergist tests using piperonyl butoxide (PBO) were carried out in sites where mosquitoes were found to be resistant to pyrethroids. To estimate insecticide resistance trends from 2004 to 2020, percentage mortalities from each site and time point were aggregated and regression analysis of mortality versus the Julian dates of bioassays was performed.

Results

Percentage of sites with pyrethroid resistance increased from 0% in 2004 to more than 80% in the 2020, suggesting resistance has been spreading geographically. Results indicate a strong negative association (p = 0.0001) between pyrethroids susceptibility status and survey year. The regression model shows that by 2020 over 40% of An. gambiae mosquitoes survived exposure to pyrethroids at their respective diagnostic doses. A decreasing trend of An. gambiae susceptibility to bendiocarb was observed over time, but this was not statistically significant (p = 0.8413). Anopheles gambiae exhibited high level of susceptibility to the pirimiphos-methyl in sampled sites.

Conclusions

Anopheles gambiae Tanzania’s major malaria vector, is now resistant to pyrethroids across the country with resistance increasing in prevalence and intensity and has been spreading geographically. This calls for urgent action for efficient malaria vector control tools to sustain the gains obtained in malaria control. Strengthening insecticide resistance monitoring is important for its management through evidence generation for effective malaria vector control decision.

Background

Long-lasting insecticidal nets (LLINs) and indoor residual spraying (IRS) have significantly contributed to the reduction of malaria globally, and remain important interventions for countries to progress towards malaria elimination [1, 2]. In order to achieve 80% nationwide LLINs coverage of at-risk populations in Tanzania, LLINs have been provided since 2009 through health facilities to pregnant women attending antenatal care and children receiving their first measles vaccine [3], to school-age children under the School Net Programme (SNP), under-five catch-up campaigns (U5CC) and universal coverage campaigns (UCC) [4]; the U5CC and UCC alone cumulatively distributed over 60 million LLINs between 2009 and 2020 [5]. As a result, the coverage of LLINs increased nationally from 38.0% in 2008 to 91.5% in 2012 [6]; in 2017, coverage was 77.9% [7].

IRS was introduced in 2007 in Muleba and Karagwe districts in northwest Tanzania [8], with subsequent gradual scale-up across the Lake Victoria regions, covering approximately 6.5 million people in 18 districts by 2012. In 2020 targeted IRS was conducted in six districts with the highest malaria burden in the Lake Victoria regions. These districts were also covered by LLINs.

The scaling-up of LLIN and IRS together with the provision of prompt malaria diagnosis and effective treatment has resulted in the reduction of malaria prevalence in Tanzania, from 14% in 2016 to 7.5% in 2017 for children aged 6–59 months [7]. The effectiveness of these malaria vector control interventions relies on a limited number of insecticides available for public health use and the continued susceptibility of Anopheles vectors to these insecticides. Unfortunately, across sub-Saharan Africa and indeed Tanzania, Anopheles mosquitoes have been developing resistance to all the classes of available insecticides used for their control [9,10,11]. The presence of insecticide resistance among malaria vectors presents a challenge to the effectiveness and long-term implementation of insecticide-based vector control interventions in malaria-endemic countries.

Metabolic and target site resistance are the two broad mechanisms by which mosquitoes develop resistance to insecticides. With metabolic resistance, mosquitoes produce increased quantities of enzymes, which either metabolize or sequester the insecticide molecules before reaching the target sites hence rendering them ineffective. A common mosquito insecticide resistance, particularly against pyrethroids and DDT, is caused by non-synonymous mutations to the knockdown resistance (kdr) gene, which is conferred by the alteration in the target site of action, i.e., the voltage-gated sodium channel. This genetic resistance mechanism has been widely studied in various insect species and linked to phenotypic resistance [12,13,14].

To safeguard the efficacy of insecticide-based vector control interventions, evidence-based strategies for preventing and managing insecticide resistance are needed. As per the Global Plan for Insecticide Resistance Management (GPIRM), the WHO urges all malaria endemic countries to develop strategies for monitoring and managing insecticide resistance for the purpose of ensuring that the limited number of insecticides available for vector control remain effective [15].

The National Institute for Medical Research (NIMR) in collaboration with the National Malaria Control Programme (NMCP) and other stakeholder began monitoring insecticide resistance in 1999 with annual monitoring surveys conducted at sentinel sites in more than 15 regions of mainland Tanzania. Formal insecticide resistance monitoring and management strategy has been implemented in the country in 2016 [16]. The strategy provides guidelines on monitoring the susceptibility of malaria vectors to insecticides and how to manage resistance.

While longitudinal Anopheles insecticide resistance data is available for Tanzania, these data have never been comprehensively analysed in full, a missed opportunity to offer important information to guide decisions on future malaria vector control efforts in the country.

The aim of this analysis was to describe the evolution and dynamics of insecticide resistance in Tanzania using available results from surveys conducted in surveillance sentinel sites across mainland Tanzania on Anopheles mosquito population susceptibility to insecticides and resistance mechanisms to guide future vector control interventions.

Methods

Study sites and design

This was an analysis of data generated from cross-sectional surveys conducted in selected sentinel surveillance sites to detect and monitor insecticide resistance in Anopheles mosquitoes (Fig. 1a and b). Data sources from published and un-published reports are summarized in Table 1. Selection of these sentinel sites was based on WHO recommended criteria, which included: a history of insecticide use in the sampled areas; malaria endemicity; high coverage of vector control interventions, i.e. LLINs and IRS; demographic setting (urban/rural); and site accessibility. The characteristics of the sentinel sites have been previously described by Kabula et al. [9, 17, 18] and in unpublished survey reports (Table 1). Nine sentinel sites were established in 1999 with funding from the Ministry of Health Gender and Elderly and Children (MoHGEC). As LLIN distribution was scaled up, and in response to the concerns around insecticide resistance within the global malaria community, the number of sentinel sites was increased to 11 in 2004, to 14 in 2012, and then to 22 in 2015. Due to changes in malaria epidemiology and geographic targeting of IRS operations, different sites have been sampled over the years, with 71 total sentinel districts having been surveyed across all of the country’s six different agro-ecological zones at one point between 2004 and 2020 (Additional file 1: Table S1). Data prior to 2004 were limited and, therefore, excluded from the analysis.

Fig. 1
figure 1

Maps of Tanzania showing the distribution of sentinel sites used for insecticide resistance monitoring by agro-ecological zone (a) and annual rainfall (b) between 2004 and 2020

Table 1 Data sources from published and unpublished reports by survey year

Mosquito collections

Between 2004 and 2008, indoor-resting adult Anopheles mosquitoes were collected from human habitations during the rainy season using mechanical aspirators between 0600–0900 h according to a WHO standard protocol for insecticide resistance monitoring [19]. However, because using adult mosquitoes collected through indoor-resting yielded mosquitoes of unknown ages and physiological status which could confound test results, from 2009 onwards Anopheles mosquitoes were reared from larvae collected from sentinel sites as described elsewhere [9]. The geographical coordinates of each larvae sampling site were recorded. Larvae collected from several different breeding sites and/or villages in the same district were pooled together for rearing and testing.

WHO susceptibility bioassays

Two to five-day old female (F1 generation) mosquitoes were tested using standard WHO insecticide susceptibility procedures [19, 20], with 20–25 mosquitoes being placed inside holding tubes before being exposed to WHO test papers treated with either 0.05% deltamethrin, 0.75% permethrin, 0.05% lambdacyhalothrin, 4% DDT, 0.1% bendiocarb or 0.25% pirimiphos-methyl; test papers were prepared at the WHO Collaborating Centre at Universiti Sains, Malaysia [20]. Prior to 2009, insecticide susceptibility tests were conducted to adult mosquitoes collected using indoor resting catches as previously described by Kabula et al. [18].

Intensity of resistance

Monitoring of insecticide resistance intensity was introduced in 2018. The WHO permethrin and deltamethrin treated papers at 5× and 10× of the discriminating concentrations (DC) were tested to assess the intensity of resistance in 17 out of 22 sentinel sites that recorded resistance to these insecticides. Mosquitoes were first tested to determine their insecticide resistance status in each district (sentinel site), and only those with confirmed resistance were assessed for intensity. Testing procedures and the determination of intensity resistance categories (i.e., low, moderate and high intensity) across locations followed the standard WHO protocol [20]. Briefly, low-intensity resistance means mosquito mortality is < 90% after exposure to 1× DC and ≥ 98% after exposure to 5× DC. Moderate-intensity resistance means mosquito mortality is < 90% after exposure to 1× DC, < 98% after exposure to 5× DC, and is ≥ 98% at 10× DC. High-intensity resistance means mosquito mortality is < 90% after exposure to the 1× DC and < 98% after exposure both to 5× and to 10× DC [9].

Synergy bioassays

To ascertain involvement of mixed function oxidases in the observed phenotypic resistance, synergy tests using piperonyl butoxide (PBO) were conducted when mosquitoes were found to be resistant to permethrin and/or deltamethrin. Additional bioassays using An. gambiae s.l. were done in surveys between 2017 and 2020 using pre-exposure to PBO, a synergist that interferes with the activity of oxidase enzymes that are involved in pyrethroid resistance. Tests were conducted in accordance with standard WHO procedures [20]. Briefly, 2–5 days old F1 adult An. gambiae s.l. mosquitoes were pre-exposed to 4% PBO paper for 1 h and immediately exposed to 0.75% permethrin, or 0.05% deltamethrin for another 1 h. Two controls were used during this experiment; the first control constituted mosquitoes exposed to clean papers neither with insecticides nor with PBO, while the second control constituted mosquitoes exposed to papers treated with PBO only [20]. The number of mosquitoes tested for each insecticide varied between 40 and 80. Mortalities recorded 24 h post-exposure were compared between the PBO synergized and the un-synergized (insecticide only) groups. The comparison was used to establish the potential role of cytochrome P450 genes in the observed resistance.

Molecular identification of members of the Anopheles gambiae complex

For each sentinel site, DNA samples from 50 randomly selected adult mosquitoes per insecticide per survey were extracted individually. Genomic DNA was extracted according to the method described by Collins et al. [21]. Molecular identification of sibling species of respective An. gambiae s.l. were conducted based on PCR methods described by Scott and colleagues [22].

Detection of knockdown resistance (kdr) alleles

Mutations associated with knockdown resistance (kdr) L1014S (east mutation) and L1014F (west mutation) in An. gambiae s.l. to pyrethroids were assayed using standard TaqMan® quantitative real-time polymerase chain reaction technique (qPCR) [13, 22], using aliquots from the DNA samples extracted for species identification.

Data analysis

Data from the same sentinel site location district and village and date of sample collection were combined for each species across all bioassays. Percent mortality and 95% confidence intervals for the WHO susceptibility tests were calculated by the binomial exact method using STATA software 13 (Stata Corp LP, College Station, TX, USA). Trend analyses were performed for all tested insecticides. To estimate insecticide resistance trends from 2004 to 2020, aggregated mosquito mortalities from each sentinel site and time were performed using regression analysis of mortality versus the Julian dates of bioassays using STATA software. Only insecticides that were tested across multiple survey years have been included in the analysis. A linear regression model was fitted for species composition as the dependent variable against time (years).

The association between resistance markers and phenotypic resistance was analysed using Chi-square tests on a subsample of An. gambiae and An. arabiensis. The subsample was weighted by the inverse of the sampling fraction i.e., subsample/total collected to represent the relative proportion in the total population. The kdr genotype frequencies among dead and live An. gambiae s.l. in susceptibility test across different survey years and clusters were compared using the Genepop software, version 4.0.

Ethical clearance

Approval was obtained from the ethics committee of the Tanzanian National Institute of Medical Research (Ref: NIMR/HQ/R.8a/VOL.IX/3604). In addition, verbal consent was sought from community leaders of each sentinel site. For mosquito larvae collection, oral consent was obtained from field owners in each location. These locations were not protected land, and the field studies did not involve endangered or protected species.

Results

Phenotypic resistance

In 2004, vectors from all surveyed sentinel districts were fully susceptible (100% mortality) to permethrin and deltamethrin (pyrethroids). Thereafter, resistance was reported over time and space. Considering each sentinel site and survey year (time) as a single data point. there has been a statistically significant negative correlation between mortality and time (p = 0.0033) (Fig. 2). Regression analysis showed that by the year 2020 over 45% and 60% of An. gambiae s.l. mosquitoes survived exposure to permethrin and deltamethrin. In 2017 An. gambiae s.l. from 10 of the 22 sentinel sites were resistant to deltamethrin i.e., Mortality < 90% (Fig. 3a and b). Anopheles gambiae s.l. showed reduced susceptibility to lambda-cyhalothrin (another pyrethroid insecticide) over time with a significant negative correlation between mortality and time (p = 0.033). The proportion of sentinel sites with evidence of resistance increased from zero in 2004 to more than 80% in 2020 (Fig. 3a and b). Generally, results indicate a strong negative association (p < 0.0001) between pyrethroids susceptibility and survey year. The estimated odds of pyrethroid susceptibility were significantly less from one survey year to the next. An. gambiae s.l. showed decreasing susceptibility levels to DDT over time and there was a negative correlation between DDT mortality and time; however, this was not significant (p = 0.4341) (Fig. 4). Although there was considerable variation among sites and years in the mortality of An. gambiae s.l. exposed to bendiocarb (Fig. 5), the decreasing trend of susceptibility over time was not statistically significant (p = 0.8413). However, between 2012 and 2014, when bendiocarb was replaced with an alternative insecticide i.e., pirimiphos-methyl in areas where IRS was being implemented susceptibility to bendiocarb increased significantly but never achieved full susceptibility (F = 23.4, p = 0.04) (Fig. 6). Anopheles gambiae s.l. exhibited high levels of susceptibility to pirimiphos-methyl at all sites tested over time (Fig. 5).

Fig. 2
figure 2

Mortality of An. gambiae s.l. exposed to permethrin and deltamethrin over time. Each point denotes the average mortality per site for that year at 95% confidence interval. The line graph shows median value for that year

Fig. 3
figure 3

Spatiotemporal trend of insecticide resistance. a Spatiotemporal trend of permethrin phenotypic resistance. b Spatiotemporal trend of deltamethrin phenotypic resistance. The solid bars indicate number of sites with resistance and the line graph indicates the percentage of sites with resistance

Fig. 4
figure 4

Mortality of An. gambiae s.l. exposed to DDT over time. Each point denotes the mortality for a single population with 95% confidence interval. The line is the linear regression best fitting line

Fig. 5
figure 5

Mortality of An. gambiae s.l. exposed to bendiocarb and pirimiphos-methyl over time. Each point denotes the mean mortality for a population in particular a sentinel site with 95% confidence interval. The line is the linear regression best fitting line

Fig. 6
figure 6

Mortality of An. gambiae s.l. from Muleba district exposed to bendiocarb 3 years after it was rotated to pirimiphos-methyl

Intensity of resistance

Overall, An. gambiae s.l. reared from eight, six, and three of the 17 sentinel sites had high, moderate, and low intensity resistance to permethrin, respectively (Fig. 7a). Similarly, An. gambiae s.l. reared from seven, five, and five of the 17 sentinel sites had high, moderate and low intensity resistance to deltamethrin (Fig. 7b).

Fig. 7
figure 7

Percentage sites per insecticides resistance intensity after exposure to different concentrations {1×, 5× and 10×} of insecticide. a Concentrations of permethrin (0.75%). b Concentrations of deltamethrin (0.05%)

Species identification

Out of An. gambiae s.l. samples analysed annually between 2010 and 2020, An. arabiensis represented 71.7% (95% CI 64.6–77.4) and An. gambiae sensu stricto (s.s.) represented 28.3% (95% CI 22–34.6%). There was a significant positive association over time with An. arabiensis (F = 6.7, p < 0.001) and negative non-significant association with An. gambiae s.s. (Fig. 8).

Fig. 8
figure 8

Species composition of An. gambiae s.s. and An. arabiensis over time. Each point denotes a percentage composition for that survey year. The dashed line is the linear regression best fitting line

Knockdown resistance (kdr) genotyping

In 2012, the kdr east-mutation (L1014S) was detected in An. gambiae s.s. in two sentinel sites at allelic frequencies of 20.8% (Geita District) and 88.9% (Ngara District). In 2013, the kdr east-mutation was detected in An. arabiensis at one sentinel site (Kinondoni District) at an allelic frequency of 17.0%. In 2014, the kdr east-mutation was detected in An. arabiensis from sentinel sites in Muheza, Kinondoni, Musoma, Ngara and Kyela Districts at allelic frequencies ranging from 6.0 to 25.0%. In 2015, kdr east-mutation was detected in An. gambiae s.s. in Muleba District at a frequency of 0.6%. In 2016, the kdr east-mutation was detected in An. gambiae s.s. in two sentinel sites with allelic frequencies of 83.0% (Muleba District) and 53.0% (Ruangwa District). In 2012, kdr west-mutation (L1014F) was detected in both An. gambiae s.s. and An. arabiensis in Muleba District at allelic frequencies of 100% and in An. arabiensis in Arumeru District at an allelic frequency of 0.6%. In 2019, kdr east-mutation was detected in An. gambiae s.s. mosquitoes in 9 (41%) of the 22 sampled sentinel sites, with allelic frequencies ranging from 3.0% in Kilwa District to 93.0% in Misenyi District (Table 2). In 2019, kdr east-mutation was detected in An. arabiensis mosquitoes in 11 (50%) of the 22 sampled sentinel sites, with allelic frequencies ranging from 2.0% in Kilwa and Geita Districts to 63.0% in Tandahimba District (Table 2). In 2020, kdr east-mutation was also detected in An. gambiae s.s. mosquitoes with allelic frequencies ranging from 15.0% in Bukoba District to 50.0% in Kigoma District Council (Table 3). The kdr east-mutation was detected in An. arabiensis with allelic frequencies ranging from 1.0% in Ngara and Ruangwa Districts to 8.0% in Misenyi District (Table 4).

Table 2 Distribution and allelic frequencies of L1014S (kdr east-mutation genotypes) in 2019
Table 3 Distribution and allelic frequencies of kdr mutations in Anopheles gambiae s.s. among the six sentinel districts registered with kdr resistance in 2020
Table 4 Distribution and allelic frequencies of kdr mutations in Anopheles arabiensis among the six sentinel districts registered with kdr resistance in 2020

Synergy bioassays

The effect of the PBO synergist was evident in sentinel sites where synergy tests were carried out with mortality rate range 85–100% for the permethrin and deltamethrin, confirming involvement of monooxygenase enzyme activity in the pyrethroid detoxification. For permethrin alone, mortality ranged from 12 to 85% (Fig. 9a). At the same sites, mortality to permethrin after exposure to PBO ranged from 67.5 to 100% (Fig. 9a). For deltamethrin alone, mortality in mosquitoes ranged from 4 to 90%. After exposure to PBO, mortality among An. gambiae s.l. exposed to deltamethrin ranged from 92 to 100% (Fig. 9b).

Fig. 9
figure 9

Mortality of An. gambiae s.l. exposed to insecticides. a Mortality of An. gambiae s.l. exposed to permethrin and permethrin-PBO. b Mortality of An. gambiae s.l. exposed to deltamethrin. The blue bars indicate mortality of An. gambiae s.l. without PBO pre-exposure, while the orange bars indicate mortality after pre-exposure to PBO

Discussion

Between 2009 and 2020, insecticide resistance in An. gambiae s.l. in Tanzania was not only increasing in frequency, it was as well spreading geographically. However, results from entomological surveys conducted between 2004 and 2009 demonstrate that An. gambiae s.l. populations were susceptible to all tested insecticides, including pyrethroids. Nevertheless, the West African leucine-phenylalanine (L1014F kdr) mutations wad detected in this time period in one sentinel site [23]. In 2010 and 2011, pyrethroid resistance was detected in An. gambiae s.l. in 34% of all the sites surveyed. By 2014, resistance to permethrin and deltamethrin was detected in An. gambiae s.l. in 64% and 70% of sampled sentinel sites, respectively. By 2020, resistance to permethrin, deltamethrin, and lambda-cyhalothrin was detected in over 90% of sampled sentinel sites. In nearly one quarter of the sentinel sites in 2020, the observed pyrethroid resistance was of high intensity, meaning that mosquitoes were able to survive the exposure to insecticides at both 5× and 10× DC. Results from these studies also show that both target site voltage gated sodium channel (vgsc) and metabolic mechanisms are associated with the reported resistance. The target sites kdr east-mutation-(L1014S) and kdr west-mutations (L1014F) associated with DDT and pyrethroids were detected within the An. gambiae s.l. populations. Furthermore, synergist test results indicate that metabolic resistance is partially or fully involved in the observed phenotypic resistance. This indicates that controlling mosquitoes using pyrethroid-based vector control interventions without factoring in a resistance management plan will be difficult.

Results of the analyses presented here confirm what was earlier reported in a series of previous studies in Tanzania. The first published report of insecticide resistance in wild An. gambiae s.l. was recorded against permethrin in Moshi in 2006, with West African leucine-phenylalanine (L1014F kdr) mutations being identified [23], followed by a report of permethrin resistance in the same area due to elevated levels of both mixed function oxidases and β-esterases [24]. Furthermore, between 2009 and 2010, resistance to DDT was observed in Ilala District in Dar es Salaam, and resistance to lambdacyhalothrin and permethrin were observed in Moshi District [9]. Resistance to other pyrethroids (permethrin and deltamethrin) was also observed in Muleba District, with bioassay mortality of anopheline mosquitoes exposed to these pyrethroids consistently less than 35% [25]. There was a sharp decrease in mortality of An.gambiae s.l. exposed to bendiocarb after two rounds of bendiocarb‐based IRS in Muleba district [26]. An. gambiae s.l. from other Districts was reported to have developed insecticide resistance to pyrethroids, DDT, and bendiocarb [10].

The rapid increase in pyrethroid resistance across Tanzania is most likely attributed to the cumulative effects of increased use of insecticide-treated nets, that began in the 1990s, followed by national scaling up of mass distribution of LLINs to achieve universal coverage in 2011 and consistently maintained to date. The contribution of LLINs in the development of pyrethroid resistance in malaria vectors has been reported throughout Africa [27,28,29]. In Tanzania, coverage of LLINs has been scaled-up in the country since 2011 when the initial 8.7 million LLINs were distributed [30]. Most of the LLINs that have been distributed since that time and are in use are treated with permethrin [7] or deltamethrin. The presence of millions of pyrethroid-treated LLINs in the field presents selection pressure for pyrethroids resistance.

IRS operations, which previously have been implicated in the selection of pyrethroid resistance [31], also may have compounded selection pressure in locations where IRS with lambdacyhalothrin was applied, such as Muleba District. Since being launched in 2007 in Muleba and Karagwe Districts, IRS was expanded in 2009 to cover up to 18 districts of Kagera, Mwanza, and Mara Region located around Lake Victoria. From 2012 onwards, IRS discontinued use of pyrethroids and only has been applied in a few targeted districts in the same Lake Victoria zone. Though not part of the NMCP’s strategy for malaria control, the frequent use of aerosol spray and coils, some of which include pyrethroids, in urban areas could have contributed to selective pyrethroid pressure in mosquito populations. Lastly, extensive use of pyrethroids in agriculture might also have contributed to the observed insecticide resistance pattern.

Given that carbamates are used for veterinary and public health purposes, the observed resistance to this class of insecticides in An. gambiae s.l. is not surprising. It is possible that the observed carbamate resistance in the selected areas is due to the use of carbamates for IRS and other insecticides as herbicides in agriculture, as reported in previous studies showing that the use of this product contributes to resistance [32,33,34]. The rise of carbamate resistance has precluded the continued use of bendiocarb for IRS in Tanzania, although use of this insecticide in neighboring Mozambique seemed to be effective despite high levels of bendiocarb resistance observed there [35].

Despite high levels of resistance to pyrethroids and carbamates, An. gambiae s.l. populations have remained susceptible to pirimiphos-methyl (organophosphate) in all sentinel sites across the country between 2015 and 2020. This may be due to the fact that organophosphates are less widely used compared to pyrethroids. Currently, pirimiphos-methyl is used for IRS in six districts only. Similar results have been reported in Malawi [36]. Earlier results obtained from Dar es Salaam showed unexpected high resistance of An. gambiae to DDT [9], and it was speculated that there could be selection due to cross-resistance associated with pyrethroids since both share the same mode of action [37].

The development of pyrethroid resistance has important implications for malaria vector control [38,39,40]. Pyrethroids have been used for both LLINs and IRS and loss of effectiveness against mosquitoes has serious consequences. The development of pyrethroid resistance in An. funestus in South Africa resulted in a sharp increase in malaria cases and consequently necessitated a change to DDT in combination with pyrethroids for IRS [41, 42]. Although IRS may be implemented with several non-pyrethroid insecticides, most of these are significantly more expensive than pyrethroids. As control programmes have been switching to non-pyrethroids as pyrethroid resistance expanded, the programmatic coverage of IRS decreased and, in some cases, NMCPs abandoned IRS altogether. Unfortunately, most available LLINs are treated with pyrethroids and resistance is a serious concern given that it has been estimated that 68% of the malaria cases averted since 2000 are attributable to LLINs [1]. There are several previous studies that shows pyrethroid resistance is undermining the effectiveness of LLINs [38,39,40, 43]. Laboratory experiments in which wild, resistant Anopheles strains directly exposed to net samples have indicated reduced efficacy of LLINs against resistant strains, while both experimental hut studies and field studies have demonstrated that resistant mosquitoes are more likely to enter LLINs, feed upon the occupants, and survive [44, 45]. A case–control study carried out in Malawi showed that among sick children aged < 5 years old reporting to the main referral hospital in the area, there was no protective effect of LLIN use [31]. The observed high pyrethroid resistance intensity might have epidemiological significance by compromising the effectiveness of LLINs in the community, and resulting in LLIN programmatic failure [46].

The spread of insecticide resistance in Tanzania described in this analysis suggests a need for insecticide resistance management for malaria vector control. WHO GPIRM provides recommendations for managing resistance in malaria endemic countries that includes the use of alternative insecticides in rotation or combination delivered through LLINs or IRS. Modelling has shown that the most promising way to delay the selection of resistance is to apply mixtures of unrelated insecticides [47]. Rather than use a non-pyrethroid insecticide to overcome resistance, an equally valid approach is to deploy a chemical synergist on the fibers. PBO is one of the synergists that overcomes resistance by inhibiting the enzymes responsible for certain types of resistance, essentially reversing resistance to pyrethroids. In both the laboratory [44, 48] and field [49] An. gambiae s.l. mortality was shown to be higher when PBO LLINs are used compared to pyrethroid-only LLINs; this study’s data showed synergy with PBO managed to restore the efficacy of permethrin and deltamethrin in killing mosquitoes in all sentinel sites. Furthermore, cluster randomized trials have shown a 44% reduction in morbidity in the area where PBO LLINs were distributed and used as compared to areas with standard pyrethroid-only LLINs [46]. These results suggest that PBO LLINs [50] may be a potential resistance management option in Tanzania.

Similarly, insecticide rotation for IRS is also a GPIRM recommended option for managing insecticide resistance. Results from our analysis indicate a significant increase in the susceptibility of bendiocarb in Muleba district after 3 years of pirimiphos-methyl use for IRS. The availability of the newly WHO prequalified insecticide products, Sumishield™ (clothianidin, a neonicotinoid) and Fludora Fusion, currently in use for IRS in Tanzania, provides more options that allow national malaria control programmes to implement suitable resistance management strategies. In addition, new insecticides and vector control tools are expected to enter the market in the near future, like new generation nets, such as Interceptor G2 [45], that are expected to add additional effective insecticide resistance management options.

In all surveys and sites, An. arabiensis was the predominant species, followed by An. gambiae s.s. The role of An. arabiensis in outdoor biting, particularly early in the evening hours (before bed time), may result in outdoor residual malaria transmission as described in parts of Tanzania [51]. Ecologically, An. arabiensis is known to have a wide range of behaviour patterns, making it more difficult to control compared to An. gambiae s.s.

One potential limitation of this analysis is that the mosquitoes used in the susceptibility testing were collected using different methods. Initially from 2004 to 2008, indoor-resting Anopheles mosquitoes of unknown ages and physiological status were collected from human habitations and tested for their susceptibility status. From 2009 onwards, Anopheles mosquitoes were reared from larvae collected from sentinel sites. In some sentinel sites, larvae were collected from only a few habitats; therefore, resistance results for the larger district area may be misrepresented.

Conclusion

The major malaria vector in Tanzania, Anopheles gambiae s.l. is currently highly resistant to pyrethroids across the country, with increasing prevalence and intensity spreading geographically since 2009. However, the efficacy of permethrin and deltamethrin was partially of fully restored with pre-exposure to PBO in all sentinel districts in which synergy study was carried out, suggesting a significant role of cytochrome P450 detoxification enzymes in the observed pyrethroid resistance. Thus, the use of PBO LLINs and other next-generation LLINs, as well as rotation of insecticides with different modes of action for IRS guided by the ongoing insecticide resistance monitoring, might be effective strategies to mitigate resistance and ensure effectiveness of these vital interventions.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

References

  1. WHO. World malaria report 2020. Geneva: World Health Organization; 2020. http://www.who.int/malaria/publications/world-malaria-report-2020/wmr2020-profiles.pdf. Accessed Nov 2021.

  2. Bhatt S, Weiss DJ, Cameron E, Bisanzio D, Mappin B, Dalrymple U, et al. The effect of malaria control on Plasmodium falciparum in Africa between 2000 and 2015. Nature. 2015;526:207–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Mushi AK, Schellenberg J, Mponda H, Lengeler C. Targeted subsidy for malaria control with treated nets using a discount voucher system in Tanzania. Health Policy Plan. 2003;18:163–71.

    Article  PubMed  Google Scholar 

  4. Universal Coverage Campaign. Final impact. http://medatanzania.org/?p=ucc_final.

  5. Yukich J, Stuck L, Scates S, Wisniewski J, Chacky F, Festo C, et al. Sustaining LLIN coverage with continuous distribution: the school net programme in Tanzania. Malar J. 2020;19:158.

    Article  PubMed  PubMed Central  Google Scholar 

  6. National Bureau of Statistics. Malaria atlas: the 2011–12 Tanzania HIV/AIDS and malaria indicator survey. Dar es Salaam, Tanzania. 2012. pp. 12.

  7. National Bureau of Statistics. Tanzania Malaria Indicator Survey. Ministry of Health, Community Development, Gender, Elderly and Children (Tanzania, Mainland); MoH Zanzibar, NBS, Office of the Chief Government Statistician, ICF. Dar es Salaam, Tanzania, and Rockville, USA; 2017.

  8. WHO. RBM-progress and impact series: country reports: focus on mainland Tanzania. Geneva: World Health Organization; 2012. p. 72.

    Google Scholar 

  9. Kabula B, Tungu P, Malima R, Rowland M, Minja J, Wililo R, et al. Distribution and spread of pyrethroid and DDT resistance among the Anopheles gambiae complex in Tanzania. Med Vet Entomol. 2014;28:244–52.

    Article  CAS  PubMed  Google Scholar 

  10. Kisinza WN, Nkya TE, Kabula B, Overgaard HJ, Massue DJ, Mageni Z, et al. Multiple insecticide resistance in Anopheles gambiae from Tanzania: a major concern for malaria vector control. Malar J. 2017;16:439.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Nkya TE, Akhouayri I, Poupardin R, Batengana B, Mosha F, Magesa S, et al. Insecticide resistance mechanisms associated with different environments in the malaria vector Anopheles gambiae: a case study in Tanzania. Malar J. 2014;13:28.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Hemingway J, Hawkes NJ, McCarroll L, Ranson H. The molecular basis of insecticide resistance in mosquitoes. Insect Biochem Mol Biol. 2004;34:653–65.

    Article  CAS  PubMed  Google Scholar 

  13. Martinez-Torres D, Chandre F, Williamson MS, Darriet F, Berge JB, Devonshire AL, et al. Molecular characterization of pyrethroid knockdown resistance (kdr) in the major malaria vector Anopheles gambiae s.s. Insect Mol Biol. 1998;7:179–84.

    Article  CAS  PubMed  Google Scholar 

  14. Ranson H, Jensen B, Vulule JM, Wang X, Hemingway J, Collins FH. Identification of a point mutation in the voltage-gated sodium channel gene of Kenyan Anopheles gambiae associated with resistance to DDT and pyrethroids. Insect Mol Biol. 2000;9:491–7.

    Article  CAS  PubMed  Google Scholar 

  15. WHO. Global Programmeplan for insecticide resistance management in malaria vectors. Geneva: World Health Organization; 2012. https://apps.who.int/iris/handle/10665/44846. Accessed Nov 2021.

  16. NMCP. Tanzania National Malaria Strategic Plan 2014–2020. 2014. https://www.out.ac.tz/wp-content/uploads/2019/2010/Malaria-Strategic-Plan-2015-2020-2011.pdf. Accessed Nov 2021.

  17. Abbott WS. A method of computing the effectiveness of an insecticide. J Econ Entomol. 1925;18:265–7.

    Article  CAS  Google Scholar 

  18. Kabula B, Tungu P, Matowo J, Kitau J, Mweya C, Emidi B, et al. Susceptibility status of malaria vectors to insecticides commonly used for malaria control in Tanzania. Trop Med Int Health. 2012;17:742–50.

    Article  PubMed  Google Scholar 

  19. WHO. Test procedures for insecticide resistance monitoring in malaria vectors, Bio-efficancy and persistance of insecticides on treated surfaces. Geneva: World Health Organization; 1998.

    Google Scholar 

  20. WHO. Test procedures for insecticide resistance monitoring in malaria vector mosquitoes. Geneva: World Health Organization; 2016.

    Google Scholar 

  21. Collins FH, Mendez MA, Rasmussen MO, Mehaffey PC, Besansky NJ, Finnerty V. A ribosomal RNA gene probe differentiates member species of the Anopheles gambiae complex. Am J Trop Med Hyg. 1987;37:37–41.

    Article  CAS  PubMed  Google Scholar 

  22. 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.

    Article  CAS  PubMed  Google Scholar 

  23. Kulkarni MA, Rowland M, Alifrangis M, Mosha FW, Matowo J, Malima R, et al. Occurrence of the leucine-to-phenylalanine knockdown resistance (kdr) mutation in Anopheles arabiensis populations in Tanzania, detected by a simplified high-throughput SSOP-ELISA method. Malar J. 2006;5:56.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Matowo J, Kulkarni MA, Mosha FW, et al. (2010) Biochemical basis of permethrin resistance in Anopheles arabiensis from Lower Moshi, north-eastern Tanzania. Malar J. 2010;9:193.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Protopopoff N, Matowo J, Malima R, Kavishe R, Kaaya R, Wright A, et al. High level of resistance in the mosquito Anopheles gambiae to pyrethroid insecticides and reduced susceptibility to bendiocarb in north-western Tanzania. Malar J. 2013;12:149.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Matowo J, Kitau J, Kaaya R, Kavishe R, Wright A, Kisinza W, et al. Trends in the selection of insecticide resistance in Anopheles gambiae s.l. mosquitoes in northwest Tanzania during a community randomized trial of long lasting insecticidal nets and indoor residual spraying. Med Vet Entomol. 2015;29:51–9.

    Article  CAS  PubMed  Google Scholar 

  27. Czeher C, Labbo R, Arzika I, Duchemin JB. Evidence of increasing Leu-Phe knockdown resistance mutation in Anopheles gambiae from Niger following a nationwide long-lasting insecticide-treated nets implementation. Malar J. 2008;7:189.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Lindblade KA, Mwandama D, Mzilahowa T, Steinhardt L, Gimnig J, Shah M, et al. A cohort study of the effectiveness of insecticide-treated bed nets to prevent malaria in an area of moderate pyrethroid resistance, Malawi. Malar J. 2015;4:31.

    Article  Google Scholar 

  29. Ndiath MO, Sougoufara S, Gaye A, Mazenot C, Konate L, Faye O, et al. Resistance to DDT and pyrethroids and increased kdr mutation frequency in An. gambiae after the implementation of permethrin-treated nets in Senegal. PLoS ONE. 2012;7:e31943.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Bonner K, Mwita A, McElroy PD, Omari S, Mzava A, Lengeler C, et al. Design, implementation and evaluation of a national campaign to distribute nine million free LLINs to children under five years of age in Tanzania. Malar J. 2011;10:73.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Sharp BL, Ridl FC, Govender D, Kuklinski J, Kleinschmidt I. Malaria vector control by indoor residual insecticide spraying on the tropical island of Bioko, Equatorial Guinea. Malar J. 2007;6:52.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Diabate A, Baldet T, Chandre F, Akogbeto M, Guiguemde TR, Darriet F, et al. The role of agricultural use of insecticides in resistance to pyrethroids in Anopheles gambiae s.l. in Burkina Faso. Am J Trop Med Hyg. 2002;67:617–22.

    Article  CAS  PubMed  Google Scholar 

  33. Nkya TE, Poupardin R, Laporte F, Akhouayri I, Mosha F, Magesa S, et al. Impact of agriculture on the selection of insecticide resistance in the malaria vector Anopheles gambiae: a multigenerational study in controlled conditions. Parasit Vectors. 2014;7:480.

    PubMed  PubMed Central  Google Scholar 

  34. Reid MC, McKenzie FE. The contribution of agricultural insecticide use to increasing insecticide resistance in African malaria vectors. Malar J. 2016;15:107.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Kloke RG, Nhamahanga E, Hunt RH, Coetzee M. Vectorial status and insecticide resistance of Anopheles funestus from a sugar estate in southern Mozambique. Parasit Vectors. 2011;4:16.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Mzilahowa T, Chiumia M, Mbewe RB, Uzalili VT, Luka-Banda M, Kutengule A, et al. Increasing insecticide resistance in Anopheles funestus and Anopheles arabiensis in Malawi, 2011–2015. Malar J. 2016;15:563.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Brogdon WG, McAllister JC. Insecticide resistance and vector control. Emerg Infect Dis. 1998;4:605–13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Churcher TS, Lissenden N, Griffin JT, Worrall E, Ranson H. The impact of pyrethroid resistance on the efficacy and effectiveness of bednets for malaria control in Africa. Elife. 2016;5:e16090.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Ranson H, N’Guessan R, Lines J, Moiroux N, Nkuni Z, Corbel V. Pyrethroid resistance in African anopheline mosquitoes: what are the implications for malaria control? Trends Parasitol. 2011;27:91–8.

    Article  CAS  PubMed  Google Scholar 

  40. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Hargreaves K, Hunt RH, Brooke BD, Mthembu J, Weeto MM, Awolola TS, et al. Anopheles arabiensis and An. quadriannulatus resistance to DDT in South Africa. Med Vet Entomol. 2003;17:417–22.

    Article  CAS  PubMed  Google Scholar 

  42. Hargreaves K, Koekemoer LL, Brooke BD, Hunt RH, Mthembu J, Coetzee M. Anopheles funestus resistant to pyrethroid insecticides in South Africa. Med Vet Entomol. 2000;14:181–9.

    Article  CAS  PubMed  Google Scholar 

  43. Kleinschmidt I, Bradley J, Knox TB, Mnzava AP, Kafy HT, Mbogo C, et al. Implications of insecticide resistance for malaria vector control with long-lasting insecticidal nets: a WHO-coordinated, prospective, international, observational cohort study. Lancet Infect Dis. 2018;18:640–9.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Tungu P, Magesa S, Maxwell C, Malima R, Masue D, Sudi W, et al. Evaluation of PermaNet 3.0 a deltamethrin-PBO combination net against Anopheles gambiae and pyrethroid resistant Culex quinquefasciatus mosquitoes: an experimental hut trial in Tanzania. Malar J. 2010;9:21.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Tungu PK, Michael E, Sudi W, Kisinza WW, Rowland M. Efficacy of interceptor(R) G2, a long-lasting insecticide mixture net treated with chlorfenapyr and alpha-cypermethrin against Anopheles funestus: experimental hut trials in north-eastern Tanzania. Malar J. 2021;20:180.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. WHO. World malaria report 2016. Geneva: World Health Organization; 2016. https://www.who.int/publications/i/item/9789241511711. Accessed Nov 2021.

  47. Tabashnik BE. Managing resistance with multiple pesticide tactics: theory, evidence, and recommendations. J Econ Entomol. 1989;82:1263–9.

    Article  CAS  PubMed  Google Scholar 

  48. Kweka EJ, Lyaruu LJ, Mahande AM. Efficacy of PermaNet(R) 3.0 and PermaNet(R) 2.0 nets against laboratory-reared and wild Anopheles gambiae sensu lato populations in northern Tanzania. Infect Dis Poverty. 2017;6:11.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Protopopoff N, Mosha JF, Lukole E, Charlwood JD, Wright A, Mwalimu CD, et al. Effectiveness of a long-lasting piperonyl butoxide-treated insecticidal net and indoor residual spray interventions, separately and together, against malaria transmitted by pyrethroid-resistant mosquitoes: a cluster, randomised controlled, two-by-two factorial design trial. Lancet. 2018;391:1577–88.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. WHO. Report of the fifteenth WHOPES Working Group Meeting 18–22 June 2012 Review Of: Olyset Plus, Interceptor LN, Malathion 440 Ew and Vectobac Gr. Geneva: World Health Organization; 2012.

    Google Scholar 

  51. 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.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We are grateful to the invaluable contributions and support from various people and institutions. We thank the Tanzanian National Malaria Control Programme and Tropical Pesticides Research Institute for their fruitful collaboration in conducting the insecticide resistance monitoring. We also wish to express our sincerely gratitude to the District Executive Directors and District Medical Officers from the respective study districts for their invaluable co-operation in allowing their District Vector Control Officers/Health Officers to participate and support the resistance monitoring activities.

Disclaimer

The opinions expressed herein are those of the authors and do not necessarily reflect the views of the President’s Malaria Initiative, the U.S. Agency for International Development, the U.S. Centers for Disease Control and Prevention, or other employing organizations or sources of funding.

Funding

The resistance surveys were initially supported by the Government of the United Republic of Tanzania followed by Bill and Melinda Gates Foundation through WHO, and thereafter by the generous support of the U.S. President’s Malaria Initiative (PMI) through the U.S. Agency for International Development’s Tanzania Vector Control Scale-Up Project, VectorLink, and RTI-Okoa Maisha Dhibiti Malaria. This analysis was supported by RTI International through Okoa Maisha Dhibiti Malaria (OMDM) activity (Cooperative agreement number: 72062118CA-00002) in collaboration with the President's Malaria Initiative, the U.S. Agency for International Development and the U.S. Centers for Disease Control and Prevention.

Author information

Authors and Affiliations

Authors

Contributions

WNK, SM, and BK conceived and designed the study. PT, BK and TN managed the study. RM, VM, DJM, PM, and PT supervised the field collection and laboratory analyses and ensured data quality. DM, SN, HM, JJ, AC, EE, ER, RR, LN, NA and NS advised on the setting and design of the study. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Patrick Tungu.

Ethics declarations

Ethics approval and consent to participate

Ethical clearance was obtained from the ethics committees of NIMR Tanzania (Ref: NIMR/HQ/R.8a/VOL.IX/3604). In addition, verbal consent was sought from community leaders of each sentinel site. For mosquito larvae collection, oral consent was obtained from field owners in each location. These locations were not protected land, and the field studies did not involve endangered or protected species. This paper is published with permission from the Director General, National Institute for Medical Research Tanzania.

Consent for publication

Not applicable.

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.

Supplementary Information

Additional file 1: Table S1.

List of sentinel sites and their respective agro-ecological zones, 2004–2020.

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 http://creativecommons.org/licenses/by/4.0/. 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 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

Tungu, P., Kabula, B., Nkya, T. et al. Trends of insecticide resistance monitoring in mainland Tanzania, 2004–2020. Malar J 22, 100 (2023). https://doi.org/10.1186/s12936-023-04508-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12936-023-04508-3

Keywords