High level of resistance in Anopheles arabiensis mosquito to pyrethroid insecticides from low malaria transmission zone of Moroto district, Karamoja region, Uganda: Implication for malaria vector control

Abstract

Background: Karamoja region of Uganda previously classified as low malaria transmission zone is currently experiencing significant upsurge of malaria incidences. Long lasting insecticidal nets (LLINs) impregnated with pyrethroids constitute a major tool for malaria control in this region. Efficacy of this tool can be hampered by resistance to the pyrethroids in the Anopheles mosquito vectors. Resistance status of these mosquitoes in this region is poorly understood, effectively hampering better understanding of the impact of LLINs in the malaria control initiative. Here, we assessed susceptibility of the Anopheles arabiensis from the region to deltamethrin, permethrin (pyrethroids) and pirirmiphos-methyl (organophosphate) insecticides.

Method: We collected anopheline mosquito larvae from their natural habitats and reared them to adult emergence in situ field insectary in Karamoja region. We then identified them morphological to species level and exposed 513 emerge adult female An gambiae s.l., mosquitoes to diagnostic dosages of deltamethrin (0.05%), permethrin (0.75%) and pirimiphos-methyl (0.25%) pyrethroids exposure using the standard WHO insecticide susceptibility test assay. Synergic assays using piperonyl butoxide (PBO) were done to check for the involvement of detoxification enzymes in pyrethroid resistant populations. We then screened for knockdown resistance (KDR) and mosquito species diversity using Polymerase Chain Reaction (PCR).

Results: Majority (96%) of the mosquitoes we sampled were identified as An. arabiensis and 4% as An. gambiae sensu stricto. We observed cross-resistance to both deltamethrin (11.9%) and permethrin (47%) but susceptibility (100% mortality) to pirimiphos-methyl in An. arabiensis. The pre-exposure to PBO ameliorated the resistance to both pyrethroids. We detected homozygous KDR -eastern variant in 1.8 and 50% of the An. arabiensis and An. gambiae s.s. respectively.

Conclusion: Anopheles arabiensis and An. gambiae s.s. are the malaria vector in Karamoja region with An. arabiensis predominating. Both species are susceptible to pirimiphos-methyl but resistant to both deltamethrin and permethrin, through a metabolic process (phenotype). Mosquotoes with genetic (kdr) mutations for resistance were minimal and hence have minimal contribution to the pyrethroid resistance profile. An. arabiensis can thus be controled in Karamoja region using deltamethrin and/or permethrin impregnated mosquito nets integrated with PBO and/or through indoor residual spraying of sprayable human dwellings with pirimiphos-methyl.

Background

Malaria is a leading cause of morbidity and mortality in the tropics despite reduction in global malaria burden [1,2] The decline has been attributed to use of long-lasting insecticidal nets (LLINs), indoor residual spraying (IRS), introduction of artemisin-based combination therapy (ACT) for malaria treatment and intermittent preventive treatment (IPTp) during pregnancy [2]. Despite these significant gains, long-term effectiveness of malaria vector controls using LLINs and IRS is currently being challenged by widespread insecticide resistance in mosquito populations [3].

World Health Organization (WHO) recommends expansion of insecticide resistance monitoring and surveillance within national and regional control programs [4]. This initiative can 1) provide insights into mechanisms that can drive resistance, 2) characterize effectiveness of control efforts and 3) identify causes of any failures within the controlled areas for effective implementation and sustainable vector control [5,6]. WHO recommend organochlorines, organophosphates, pyrethroids and carbamates for control of malaria vectors, among which pyrethroids are the only WHO prequalified insecticides for LLINs, and are also extensively used for IRS due to their low cost, human safety and long duration of residual action [4]. However, there are other novel pyrolle and neonicotinoid insecticide compounds like pyrolle chlorfenapyr and neonicotinoid clothianidin that are used in LLINS and IRS in combination with pyrethroids and all have different mode of actions [7–9] The LLINs and IRS are widely promoted and implemented as public health intervention tool for malaria control in most malaria-endemic countries including Uganda [10]. In most of Uganda, LLINs is selectively implemented in the different regions, with an average operational coverage of over 95% [11]. The control programmes are concentrated in high/medium transmission zones while low malaria transmission zones, including Karamoja region have largely been neglected [12,13] Despite extensive LLINs implementation, over the recent years Karamoja experienced significant (> 60%) malaria incidences between 2015 to 2017 [11,14] with Moroto district registering 334.5 cases per 1000 children under 5 years [15]. This translates into about 33 % disease prevalence and may be linked to factors such as pyrethroid resistance that might be impeding efficacy of the LLINs [16–18]

Pyrethroids exert their insecticidal effect on the voltage-gated sodium channel (VGSC) located on the membrane of neurons [19]. When pyrethroids bind an open channel, they block its closure, thus extending the action potential and resulting in the insect’s rapid paralysis, a phenomenon known as “knockdown resistance” (kdr) [20]. However, mosquito resistance to pyrethroids is an emerging challenge to efficacy of otherwise successful insecticide-treated nets (ITN) based malaria control intervention against adult vector [21]. This is particularly of great concern to sub-Saharan Africa, with high malaria transmission levels and major vectors (An. funestus, An. gambiae s.s. and An. arabiensis) resistant to pyrethroids [22]. The resistance is mediated through knockdown resistance (kdr) mutations, enhanced detoxification of pyrethroids by the mosquito enzymes (esterases, monooxygenases and glutathione S-transferases [23,24] and probably other unknown mechanisms such as behavioral and penetration[25,26] known to occur in other vectors [27,28]

This present study was initiated to establish distribution of major malaria vectors and their resistance status to deltamethrin, permethrin (pyrethroids) and pirirmiphos-methyl (organophosphate) insecticides in Karamoja region of Uganda. Our findings, reported herein, will contribute to designing suitable control interventions and improve implementation of resistance management strategies across Karamoja region of Uganda.

Methods

Study area

Karamoja lies between 1°30′ to 4° 06′ N and 33° 30′ to 35° covering an area of 29,430 km² forming the north-eastern part of Uganda with a human population of about 1.1million [29]. This region encompasses Abim, Amudat, Kaabong, Karenga, Kotido, Moroto, Nabilatuk, Napak and Nakapiripit districts. The region has 53.8 to 63.5% Human Poverty Indices (HPI) compared to 37.5 % Uganda national average [30]. Most of the populations are nomadic pastoralist. The region is characterized by savannah woodlands and semi-desert vegetation with 1200 and 1500m plateaus and Kadam, Akisim and Napak Mountains and Moroto mountain (about 3083m above sea level) [31]. The region is drained by Turkwel River, part of internal drainage basin of Lake Turkana in Kenya. Karamoja experiences only one rain season which occurs from May to July with malaria transmission peak picking up immediately after the rains. The mean annual rainfall level is between 350 and 750 mm [32]. Ambient temperatures range from 16ºC to 30ºC. The species of Anopheles mosquitoes in the and Anopheles gambiae sensu stricto Giles 1902 [33]. Malaria transmission in the region is absolutely (100%) attributed to Plasmodium falciparum parasite [34]. Malaria prevalence in the region has increased by 30% with average incidence ranging from 166 in 2015 to 295 in 2018 per 1000 people [11]. The LLINs which is the only malaria vector control strategy employed in the region with 57% households owning more than one ITN with ratio of use and access of 0.95 [13,35] This indicate low ratio of people sleeping inside nets.

Mosquito collections

We collected mosquito larvae stages 3, 4 and pupae from four aquatic breeding habitats in Moroto (N02°32.0′ E34°40.0′) in the rainy month of July 2019 using dipping technique. The sites included brick pits and puddles (N2°35'4" E34°41'50"), road side drains (N2°31'52" E34°39'22") and excavations (N2°31'13" E34°31'21"), sand pits and open gardens (N2°28'57" E34°37'7"). Larvae were collected once within the three days in the month of July 2019. Immediately after sampling, the larvae were transferred to a field insectary where they were maintained at a density of 500 larvae per three liters in their natural water that contained dissolved natural foods. We kept the emergent adult mosquitoes in standard 30×30×30 cm cages in under the ambient insectary environmental conditions. We identified An. gambiae s. l among emerged mosquito adults using morphological keys [36] offered them 10% glucose solution ad libitum and tested the females for insecticide resistance 2-5 days post emergence. We obtained ethical clearance for mosquito collection from Uganda National Council for Science and Technology (authorization No SS 4610).

Insecticide susceptibility bioassay

We conducted insecticide susceptibility tests on 2-5 days old non-fed adult female An. gambiae s.l. mosquitoes using the standard WHO tube bioassay protocol for assessing potential insecticide resistance [37] at Gulu University mosquitoes insectary. To simulate natural conditions, we maintained an average temperature in the insectary at 29 ± 2°C (day) and 24 ± 2°C (night), with relative humidity (RH) ranging from 57 to 70 and a photoperiod of 12:12 h (L–D). Our overall maintenance of the colony followed standard operating procedure for rearing Anopheles mosquitoes [38,39] We exposed 492 female An. arabiensis to 0.05% deltamethrin, 0.75% permethrin or 0.25% pirimiphos-methyl concentrations on insecticide-impregnated papers or control (oil-treated) (Vector Control Research Unit, Universiti Sains Malaysia) for 60 minutes. During this exposure duration, we recorded the number of mosquitoes knocked-down at 10, 15, 20, 30, 40, 50 and 60 minutes intervals post exposure. After the exposure period, mosquitoes were transferred to holding tubes and maintained on 10% glucose solution. The final mortalities were determined after 24 h post exposure. After the bioassays, we collected dead (susceptible) and alive (resistant) mosquitoes, and stored them individually in separate Eppendorf tubes with silica gel for subsequent molecular laboratory analysis.

Evaluation of potential impact metabolism on resistance to the pyrethroids in the mosquitoes

To evaluate if metabolism of the pyrethroids by the mosquitoes was responsible for the insecticide resistance we observed, we assessed relative mortality of the mosquitoes exposed to the pyrethroids with or without piperonyl butoxide (PBO) that prevents pyrethroid catabolism by Cytochrome P450 oxidase that can detoxify the active ingredient before an insecticidal effect can occur [40]. This inhibition produces a synergistic effect [41]. Briefly, we assessed this effect using WHO protocol [37] where we pre-exposed adult female An. gambiae s.l to WHO papers impregnated with 4% PBO an oxidase inhibitor for 1 hour. This was followed immediately with exposure to discriminating doses of 0.05% deltamethrin and 0.75% permethrin on impregnated papers for 60 minutes. Five replicated were performed for each exposure set. We recorded the knock-down rates at 10, 15, 20, 30, 40, 50- and 60-minutes intervals during the one-hour exposure to synergist and to insecticides (deltamethrin or permethrin). We immediately assessed the resultant mortality of the flies from the three treatments (PBO + deltamethrin or permethrin exposed or unexposed mosquitoes). The live mosquitoes were fed on 10% glucose solution. Mortality rates from assays conducted with and without exposure to PBO (synergist) were scored after 24 hours of exposure. Mortality rates were compared between PBO deltamethrin or permethrin exposed and unexposed mosquitoes. Resistance was also compared with pirimiphos-methyl to ascertain which test mosquitoes were fully susceptible. All mosquitoes tested were identified to species level by PCR as described in mosquitoes species identification subsection below.

PCR identification of memebers of An. gambiae s.l. and An. funestus s.l.

We identified our mosquitoes (post exposure) to their An. gambiae s.l sibling species status using PCR methods of Scott et al. [42], and for members of the An. funestus s.l. group using a method of Koekemoer et al. [43]. We extracted total genomic DNA (gDNA) from whole flies using the DNeasy blood and tissue kits (Qiagen, Valencia, CA) following the manufacturer’s protocols and assessed the quality of the DNA on 1.5% agarose gel as visualized on Gel Doc Imaging System (UVITEC, Cambridge). Our primers included those specific to An. gambiaes.s, An. arabiensis as well as the universal An. gambiae s.l. complex primer [42]. Similarly, specific primers for An. funestus s.s. and universal primer for the An. funestus group were also used for the An. funestus s.l. [43]. These two Ribosomal DNA polymerase chain reaction methods for An. gambiae s.l. and An. funestus complex are based on species-specific nucleotide sequences in the ribosomal DNA internal transcribed spacer region 2 (ITS2) and are used to identify both species and interspecies hybrids in mosquitoes regardless of the life stage using extracted DNA. The primer sequences are derived from the 28S coding region and the 5’ end of the internecine spacers (IGS) of the members of the An. gambiae and An. funestus complexes. In the PCR reaction, we amplified 2.5µl gDNA with 1 unit of GoTaq Green Master Mix (Promega, Madison, MO) in the buffer in a total volume 12.5 µl. For An. funestus complex, conditions remained the same except that we added 0.5 mM of MgCl₂. We run the reactions in touch screen thermal cycler (SimpliAmp, Applied Biosystems, Life Technologies, Singapore). The first cycle included five minutes at 95°C, 30 seconds at 50°C, and 30 seconds at 72°C. Subsequent cycles involved 1 minute at 94°C, 30 seconds at 50°C, and 1 minute at 72 °C for 30 cycles. We also run positive and no-sample negative controls. We loaded The PCR products onto 1X SYBR safe (Invitrogen, 5791 Allen Way Carlsbad CA 92008, USA) 2 % agarose gels (AppliChem GmbH Ottoweg Damstadt Germany) in a TBE (40mM Tris-HCl pH 8.3, 45 mM boric acid, and 1mM EDTA) buffer and run a 100 bp DNA ladder molecular weight marker (Life Technologies, Rockville, MD) to confirm expected molecular weights of the amplification products. We documented our PCR products using GelDoc Imaging System (UVITEC, Cambridge) and scored expected PCR product sizes for An. gambiae s.s at 390bp, An. arabiensis 315 bp and 505bp An. funestus [42,43].

Detection of East African kdr resistance mutations in An. gambiae s.l.

We also used the methods described by Ranson et al. [24] to detect single base pair specific single nucleotide polymorphisms (SNPs) leucine to serine substitution TTA/TCA mutation in the voltage-gated sodium channel known as knockdown resistance to DDT and pyrethroids in East Africa. Briefly, we amplified 5 µl gDNA, 0.2μM of the specific primers (AgD1, AgD2, AgD4 and AgD5) with 1 unit of GoTaq Green Master Mix (Promega, Madison, MO) in the buffer in a total volume 25µl. We run the reactions in touch screen thermal cycler (SimpliAmp, Applied Biosystems, Life Technologies, Singapore). We used touch down PCR conditions with the initial denaturation at 95ºC for 5 minutes followed by 10 cycles of denaturation at 94ºC for 1 minute, annealing at 54ºC for 30 seconds and extension of 72ºC for 30 seconds. This was followed by 30 cycles of denaturation at 94ºC for 1 minute, annealing at 47ºC for 30 seconds and extension at 72ºC for 30 seconds with the final extension of 72ºC for 10 minutes and holding at 4ºC until when collected. We also run positive and no-sample negative controls. PCR products analysis was carried out as described in mosquito species identification subsection above.

Data analysis

We corrected the knockdown rates for testing the toxicity of each insecticide, where the knockdown in our control treatment was greater than 5% but less than 20% using Abbott’s formula [44]and then transformed them to Probits [44] for linear regression analysis and the determination of 50 % knockdown (KD_t50). For Probit analysis, we used GraphPad Prism version 7.00 for Mac (GraphPad Software, La Jolla California USA). We used weighted mean summarize knockdown due to different insecticides and adopted WHO criteria of interpretation of our results. Consequently, we considered mosquito population to be 1) resistant (confirmed) to a particular insecticide if mortality rate was 0-79%, 2) resistant suspected if mortality was 80-97% and more investigations are required, and 3) susceptible when the mortality was 98-100%. We summarized data for the metabolic tests as mean percent mortality of the five replicates, and calculated 95% confidence intervals around those means.

Results

Susceptibility of An. arabiensis to pyrethroids

Results of bioassay of adult female An. gambiae s.l. (N= 513) for susceptibility to the three insecticides are summarized in Figure 1, 2 and Table 1. Based on WHO criteria, [37], resistance (< 95% mortality) of An. arabiensis to deltamethrin and permethrin was observed in all the samples, with mortality varying from 12% to 91% in Moroto (Table 2). Susceptibility of An. arabiensis to pirimiphos-methyl was observed in Moroto (Table 2).

Impact of PBO on toxicity permethrin and deltamethrin to An. arabiensis   

An hour pre-exposure of An. arabiensis to PBO enhanced toxicity of deltamethrin (from 11.8 to 90%) or permethrin (from 47 to 89.2%) in Moroto, compared to cohorts directly exposed to each of the two candidate insecticides without PBO pre-exposure  (Figure 1 and Table 2), suggesting potential cytochrome P450 oxidase mediated resistance [45] to these insecticides in the mosquito. However, our observations that the pre-exposure to PBO partially abolished resistance to deltamethrin and permethrin in Moroto still suggests partial role of the cytochrome P450 oxidase in the resistance phenotype while other mechanisms might also play a role (Figure 1 and Table 2).

In the absence of PBO, the LT₅₀ for deltamethrin and permethrin on An. arabiensis were similar (Table 3 and Figure 3). However, the PBO significantly reduced the LT₅₀ for both deltamethrin and permethrin on An. arabiensis relative to their respective native formulations (without PBO). Additionally, incorporation of PBO reduced the LT₅₀ in permethrin more than it did with deltamethrin (Table 3).

An. gambiae s.l. species identification

Our subsequent PCR of sibling species on most of the mosquitoes (N = 342) post exposure revealed that most of the mosquitoes were An. arabiensis (96%) with the rest being An. gambiae s.s. (4%). No An. funestus species was detected.

Table 1: Knockdown summary during 60 minutes exposure

*KT50: Time after which 50% of the Anopheles tested are knocked down.

*KD60: Proportion of Anopheles knocked down after 60 minutes.

Table 2: Mean percentage mortalities following exposure of 2-5 day old An. arabiensis to pirirmiphos-methyl, deltamethrin and permethrin. Mortalities were recorded 24 hours post exposure (for 60-minutes)

Table 3: Median Lethal Time for knockdown in Adult female An. arabiensis mosquitoes (from Moroto district, Uganda) by various formulations of pyrethroids.

PBO – Piperonyl Butoxide, (+) - with PBO, (-) - Without PBO, CI- Confidence Interval, LT - Median lethal Time.

Prevalence of East African  (L1014S) knockdown resistance (kdr) point mutations in An. arabiensis and An, gambiae s.s. in Moroto district

We have summarized our distribution of L1014S mutations from genotyping of An. arabiensis (n=328) and An. gambiae s.s. (15) in Table 4. Of the 328 An. arabiensis tested for kdr east allele, 66.2% (n=217) were homozygous for wild type alleles (SS) and 1.2% (n=4) were homozygous for mutation alleles (RR). Of the 15 An. gambiae s.s, 33.3% (n=5) were homozygous for the susceptible wild type (SS), 13.3% (n=2) were homozygote mutation alleles (L1014S) and only 6.7% (n=1) was heterozygous (HR).

Table 4: KDR genotypes

KDR - knockdown resistance, PBO – Piperonyl Butoxide, # - number. SS represents homozygous for wild type alleles, RR homozygous mutation alleles and HR heterozygous mutation alleles.

Discussions

We report the first results of resistance in An. arabiensis against pyrethroid in Karamoja region. We established An. arabiensis as the major malaria vector in Karamoja areas followed by An. gambiae s.s. As in other studies, the major malaria vectors in Uganda include An. gambiae s.s. An. arabiensis, An. Bwambae, and An. funestus, with 90% of malaria transmission in the country being attributed to An. gambiae s.l [46–48]. An. arabiensis and An. gambiae s.s. reported in our study have different genetics and behavior that can have important implications for the epidemiology of malaria and their control. The current use of LLIN in Karamoja region for malaria vector control might be less effective against An. arabiensis. This is mainly due to resistance of the mosquito to the insecticides (deltamethrin and permethrin), more exophagic, exophilic and zoophilic nature of An. arabiensis.

An. arabiensis being less killed by they two insecticides deltamethrin and permethrin commonly used for the treatment of LLIN. Similar studies have been reported in Kenya and Tanzania [49,50] This happens as a result of the more exophagic, exophilic and zoophilic tendencies of An. arabiensis compared to An. gambiae s.s. behaviors. As shown in other studies, An. arabiensis have peak biting activities happening during the early parts of the night when most people have not gone to bed [51,52] . What this could means is that the early evening and outdoor biting of An. arabiensis can compromise on the efficacy of the LLINS in Karamoja region as seen in Ethiopia [52] and Tanzania [49]. The exophagic and exophilic nature of An. arabiensis similarly can affect the performance of IRS since these mosquitoes tend to rest and feed outdoors [53].

The high prevalence of malaria seen in Karamoja region as describe in other studies [11,14] may be linked to population dynamics of An. arabiensis in the region. Exophagic and exophilic nature of An. arabiensis, increases the outdoor biting and resting of these mosquitoes leading to high malaria transmission in the area. Besides, An. arabiensis, bite early evening or morning when the population is not protected by LLIN potentially sustaining residual malaria in area with high coverage of net [54,55]. Additionally communities in Karamoja region thrive on pastoralism looking for pasture and water for their livestock which pre-expose them to outdoors activities when keeping the livestock and increasing the chances of more outdoor mosquitoes’ bites and thus more malaria in the area.

In this study, we confirm existence of both mechanisms of insecticide resistances: knock down resistance (kdr) and metabolic resistance in the An. gambiaes.s. and An. arabiensis in Karamoja region. We confirm the restoration of An. arabiensis susceptibility following the pre-exposure to PBO indicating increased activity of detoxification of enyzymes which are contributing to pyrethroid resistance in Karamoja region. However, future further studies are needed in characterizing these mechanisms of insecticide resistances. Besides, future planned roll out of LLINs enhanced with PBO might not be a magic bullet since more than one resistance mechanisms were observed in An. gambiae s.l. in Karamonja region.

While the two pyrethroids seems to be having similar performance, PBO significantly enhances the performance of both deltamethrin and permethrin and further enhances that of permethrin over that of deltamethrin (Table 3). This suggests that permethrin incorporated with PBO would have a much better performance in addressing the issues of pyrethroids resistance in Moroto than deltamethrin. These insecticide resistances could as well have important implications for the epidemiology of malaria as well as malaria vector control in the low malaria transmission zones of Karamoja areas. Although pyrethoids are currently used in all LLINS as recommended by the World Health Organization, there are also used in combination with either synergists (PBO) [55] or other active ingredient like clothianidin, chlorfenapyr in Interceptor G2 [7,8]. However, pyrethroid resistance of malaria vectors is still widespread in Africa as well as other classes of insecticides [23,37,56]. Increased resistance has been attributed to selection pressure from the scale-up of LLINs and IRS [23,57] and use of similar classes of insecticides in agriculture [58], although the relative contribution of these mechanisms varies by area [59].

Given that resistance has been reported in An. gambiae s.s and An. funestus [57,60–66] that are the major malaria vector in Uganda and now An. arabiensis in this current study in Karamoja region, this is a major threat to LLINs and IRS use in malaria vector control programmes in the country and the whole of the neighbouring countries. This appears to reflect the need to monitor malaria vector resistance. As demonstrated in other studies on the host feeding and insecticide resistance, the fitness costs associated with insecticide resistance can have influence on malaria transmission directly by altering host seeking feeding and mating behaviours [67–69] or fecundity [70] and reducing mosquito life span [71,72] It can also indirectly influence malaria transmission by impairing parasite development inside the mosquitoes [70,73,74].

Current WHO guidelines recommend combining ITNs and IRS in various transmission settings, especially in areas with holoendemic and epidemic malaria [5] . LLINs and IRS could be employed together in the same households in Karamoja region. However, due to unique bionomics of An. arabiensis, local anthropological and variable environmental factors in Karamoja region, the combined entomological outcomes of IRS and LLIN interventional trails can be carried out. Besides, insecticide resistance management should rely on tactical deployment of the active ingredients used for IRS and on LLINs in rotation, combinations (particularly LLINs), mosaics and mixtures as recommended in other studies [5]. What we observed would be important is, if LLINs are to be combined with IRS for malaria prevention and control, the selection of appropriate LLIN types and design should put in consideration the housing structures and cultural issues concerning in-house sleeping patterns and IRS chemicals should be done with caution to avoid further exacerbating existing resistance.

Implication for future malaria vector control

Nevertheless, the increasing intensity of resistance currently seen in An. arabiensis mosquito populations in Karamoja region will reduce the efficacy of pyrethroid-based interventions. The right course of action for malaria vector control program is therefore to adopt a proactive approach and modify current practices so as to delay the spread of resistance and preserve the effectiveness of deltamethrin and permethrin insecticides.

The high frequencies or intensities of resistance can lead to failure of IRS and can thereby have an epidemiologically significant effect on malaria incidence [75]. There is also cross-resistance between different classes of insecticides that share the same mode of action. The existence of cross- resistance and multiple resistance restricts the choice of alternative insecticides in situations where resistance has been detected. Besides, the impact of the observed spread of resistance will have effectiveness of current vector-control programmes. Therefore, a key element of effective resistance management is the use of alternations, rotations, or sequences of different insecticide mode of action classes. Other studies have shown that even in the presence of pyrethroid resistance, insecticide- treated mosquito nets (ITNs) perform better than untreated nets in terms of protection against blood-feeding, and ITNs can induce significant mosquito mortality [76].

We also recommend the use of Larval Source Management (LSM) in order to prevent the completion of development of the immature stages. Mosquitoes breeding sites in Karamoja region can easily be mapped out during the dry season. LSM can be done through 1) Larviciding; the regular application of biological or chemical insecticides to water bodies and 2) biological control: the introduction of natural predators into water bodies.

Limitations

We acknowledge the limitations of the current study including:

• The time constraints of conducting this research during.
• The sites sampled were predominantly suited for gambiae s.l breeding, thus findings don’t rule out possibility existence of An. funestus s.l in the area.

Conclusions

Our findings show that Anopheles arabiensis and An. gambiae s.s. are the malaria vectors in Karamoja region with An. arabiensis predominating. Both species are still susceptible to pirimiphos-methyl but resistant to both deltamethrin and permethrin, through a metabolic process (phenotype). An. arabiensis with genetic (kdr) mutations for resistance were minimal and we think this is having minimal contribution to the pyrethroid resistance profile seen in the region. However, we confirm restoration of An. arabiensis susceptibility following the pre-exposure to PBO indicating elevated levels of detoxifying enzymes in these vectors in Karamoja region. Therefore gudielines for management and controlling An. arabiensis insecticide resistance in Karamoja region are needed. An. arabiensis in Karamoja region can be controled using deltamethrin and/or permethrin impregnated mosquito nets integrated with PBO and/or through indoor residual spraying of sprayable human dwellings with pirimiphos-methyl.

List Of Abbreviations

IRS: indoor residual spraying

LLINs: long lasting insecticide-treated nets

WHO: World Health Organization

KDR: knockdown resistance

PCR: Polymerase chain reaction assay

PBO: Piperonyl butoxide

MFOs: mixed-function oxidases

Declarations

Authors’ contributions

RE, JI and DO conceived, contributed design of the study, field collections, performed laboratory work, analyzed the data, and drafted an initial version of the manuscript. JA, PM, GMM and TI performed laboratory work and analyzed the data. EAO conceived, designed the study, coordinated fieldwork and provided guidance. All authors read and approved the final manuscript.

Author details

¹Faculty of Science, Gulu University, P. O. Box 166,Gulu, Uganda,

²Gulu University Biosciences Research Laboratories, P. O. Box 166,Gulu, Uganda.

³Department of Biochemistry, Biotechnology Research Institute - Kenya Agricultural and Livestock Research Organization, Kikuyu, Kenya

Acknowledgements

We are grateful to the communities in in Moroto district in Karamoja region for participating in this research.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

The authors declare that all the main data supporting the findings of this study are available within the article. Any additional data sets are available from the corresponding author upon reasonable request.

Consent for publication

Not applicable.

Ethics approval and consent to participate

This study was approved by Gulu University Ethical Review Committee. Formal approval to conduct the study was granted by the Uganda National Council for Science and Technology and the Office of the Ugandan president (SS4610). Community leaders in Moroto district provided written informed consent.

Funding

The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported through the DELTAS Africa Initiative grant # DEL-15-011 to THRiVE-2 awarded through Career Development Award to Dr. Richard Echodu. The DELTAS Africa Initiative is an independent funding scheme of the African Academy of Sciences (AAS)’s Alliance for Accelerating Excellence in Science in Africa (AESA) and supported by the New Partnership for African’s Development Planning and Coordinating Agency (NEPAD Agency) with funding from the Wellcome Trust grant # 107742/Z/15/Z and the UK government. The views expressed in this publication are those of the author(s) and not necessarily those of AAS, NEPAD Agency, Wellcome Trust or the UK government.

References

[1]      World Health Organization, World Health Statistics 2015, 2015. doi:10.1016/j.jbankfin.2017.09.006.

[2]      World Health Organization, World malaria report 2018, 2018. www.who.int/malaria%0Ahttps://apps.who.int/iris/bitstream/handle/10665/275867/9789241565653-eng.pdf?ua=1%0Ahttps://www.who.int/malaria/publications/world-malaria-report-2018/en/; consulté le 22/03/2019%0Ahttps://www.who.int/malaria/media/world-malaria-rep.

[3]      World Health Organization., Global report on insecticide resistnace in malaria vectors: 2010-2016, World Heal. Organ. (2018) 2010–2016.

[4]      WHO, Conditions for deployment of mosquito nets treated with a pyrethroid and piperonyl butoxide, Glob. Malar. Program. World Heal. Organ. 2017 (2017) 1–5. doi:10.1007/s40122-017-0070-9.

[5]      WHO, The technical basis for coordinated action against insecticide resistance:, World Heal. Organ. Tech. Rep. Ser. (2010) 1–42.

[6]      I. Dusfour, J. Vontas, J.P. David, D. Weetman, D.M. Fonseca, V. Corbel, K. Raghavendra, M.B. Coulibaly, A.J. Martins, S. Kasai, F. Chandre, Management of insecticide resistance in the major Aedes vectors of arboviruses: Advances and challenges, PLoS Negl. Trop. Dis. 13 (2019) 1–22. doi:10.1371/journal.pntd.0007615.

[7]      S. Camara, L.P. Ahoua Alou, A.A. Koffi, Y.C.M. Clegban, J.P. Kabran, F.M. Koffi, K. Koffi, C. Pennetier, Efficacy of Interceptor® G2, a new long-lasting insecticidal net against wild pyrethroid-resistant Anopheles gambiae s.s. from Côte d’Ivoire: A semi-field trial, Parasite. 25 (2018) 1–11. doi:10.1051/parasite/2018042.

[8]      K. Dagg, S. Irish, R.E. Wiegand, J. Shililu, D. Yewhalaw, L.A. Messenger, Evaluation of toxicity of clothianidin (neonicotinoid) and chlorfenapyr (pyrrole) insecticides and cross-resistance to other public health insecticides in Anopheles arabiensis from Ethiopia, Malar. J. 18 (2019) 1–11. doi:10.1186/s12936-019-2685-2.

[9]      A. Sovi, C. Keita, Y. Sinaba, A. Dicko, I. Traore, M.B.M. Cisse, O. Koita, D. Dengela, C. Flatley, E. Bankineza, J. Mihigo, A. Belemvire, J. Carlson, C. Fornadel, R.M. Oxborough, Anopheles gambiae (s.l.) exhibit high intensity pyrethroid resistance throughout Southern and Central Mali (2016-2018): PBO or next generation LLINs may provide greater control, Parasites and Vectors. 13 (2020) 1–16. doi:10.1186/s13071-020-04100-7.

[10]    S.G. Staedke, M.R. Kamya, G. Dorsey, C. Maiteki-Sebuguzi, S. Gonahasa, A. Yeka, A. Lynd, J. Opigo, J. Hemingway, M.J. Donnelly, LLIN Evaluation in Uganda Project (LLINEUP) - Impact of long-lasting insecticidal nets with, and without, piperonyl butoxide on malaria indicators in Uganda: Study protocol for a cluster-randomised trial, BMC Trials. 20 (2019) 1–13. doi:10.1186/s13063-019-3382-8.

[11]    N.M.C.D. Uganda Ministry of Health, National Malaria Annual Report 2017-2018, Uganda Minist. Heal. (2018) 2–72. file:///Users/walusimbidennis/Downloads/Malaria Annual Report July 2017 web (2) (3).pdf.

[12]    Uganda Bureau of Statistics, Environmental, Demographic, Socio-economic, Production and Macroeconomic sectors, Uganda Bur. Stat. (2015) 1–353. doi:10.1017/CBO9781107415324.004.

[13]    U.P.M.I. USAID, President ’ S Malaria Initiative Angola, USAID. (2018) 1–111.

[14]    Uganda Ministry of Health, Uganda Malaria Bulletin Issue 9, January - March 2015, (2015) 1–10. https://www.pmi.gov/docs/default-source/default-document-library/implementing-partner-reports/uganda-malaria-bulletin-issue-9.pdf.

[15]    J. Ssempiira, J. Kissa, B. Nambuusi, C. Kyozira, D. Rutazaana, E. Mukooyo, J. Opigo, F. Makumbi, S. Kasasa, P. Vounatsou, The effect of case management and vector-control interventions on space-time patterns of malaria incidence in Uganda, Malar. J. 17 (2018) 1–11. doi:10.1186/s12936-018-2312-7.

[16]    K.H. Toé, C.M. Jones, S. N’fale, H.M. Ismai, R.K. Dabiré, H. Ranson, Increased pyrethroid resistance in malaria vectors and decreased bed net effectiveness Burkina Faso, Emerg. Infect. Dis. 20 (2014) 1691–1696. doi:10.3201/eid2010.140619.

[17]    R. Pwalia, J. Joannides, A. Iddrisu, C. Addae, D. Acquah-Baidoo, D. Obuobi, G. Amlalo, S. Akporh, S. Gbagba, S.K. Dadzie, D.K. Athinya, M.P. Hadi, H.P. Jamet, J. Chabi, High insecticide resistance intensity of Anopheles gambiae (s.l.) and low efficacy of pyrethroid LLINs in Accra, Ghana, Parasites and Vectors. 12 (2019) 1–9. doi:10.1186/s13071-019-3556-y.

[18]    J.M. Riveron, F. Watsenga, H. Irving, S.R. Irish, C.S. Wondji, High Plasmodium Infection Rate and Reduced Bed Net Efficacy in Multiple Insecticide-Resistant Malaria Vectors in Kinshasa, Democratic Republic of Congo, J. Infect. Dis. 217 (2018) 320–328. doi:10.1093/infdis/jix570.

[19]    L.M. Field, T.G. Emyr Davies, A.O. O’Reilly, M.S. Williamson, B.A. Wallace, Voltage-gated sodium channels as targets for pyrethroid insecticides, Eur. Biophys. J. 46 (2017) 675–679. doi:10.1007/s00249-016-1195-1.

[20]    N.M. Bowman, K. Akialis, G. Cave, R. Barrera, C.S. Apperson, S.R. Meshnick, Pyrethroid insecticides maintain repellent effect on knock-down resistant populations of Aedes aegypti mosquitoes, PLoS One. 13 (2018) 1–14. doi:10.1371/journal.pone.0196410.

[21]    J. Etang, F. Chandre, P. Guillet, L. Manga, Reduced bio-efficacy of permethrin EC impregnated bednets against an Anopheles gambiae strain with oxidase-based pyrethroid tolerance, BMC Malar. J. 3 (2004) 1–7. doi:10.1186/1475-2875-3-46.

[22]    A.P. Mnzava, T.B. Knox, E.A. Temu, A. Trett, C. Fornadel, J. Hemingway, M. Renshaw, Implementation of the global plan for insecticide resistance management in malaria vectors: Progress, challenges and the way forward, Malar. J. 14 (2015) 1–9. doi:10.1186/s12936-015-0693-4.

[23]    H. Ranson, R. N’Guessan, J. Lines, N. Moiroux, Z. Nkuni, V. Corbel, Pyrethroid resistance in African anopheline mosquitoes: What are the implications for malaria control?, Trends Parasitol. 27 (2011) 91–98. doi:10.1016/j.pt.2010.08.004.

[24]    J. Hemingway, H. Ranson, Insecticide Resistance in Insect Vectors of Human Disease, Annu. Rev. Entomol. 45 (2000) 371–391. doi:10.1146/annurev.ento.45.1.371.

[25]    A.F. Spíndola, C.S.A. Silva-Torres, A.R.S. Rodrigues, J.B. Torres, Survival and behavioural responses of the predatory ladybird beetle, Eriopis connexa populations susceptible and resistant to a pyrethroid insecticide, Bull. Entomol. Res. 103 (2013) 485–494. doi:10.1017/S0007485313000072.

[26]    G. Rt, G. Lm, S. Jm, C. Robin, Behavioural response to combined insecticide and temperature stress in natural populations of Drosophila melanogaster ., J. Evol. Biol. 29 (2016) 1030–1044. doi:10.1111/jeb.12844.

[27]    C. Nansen, O. Baissac, M. Nansen, K. Powis, G. Baker, Behavioral avoidance - will physiological insecticide resistance level of insect strains affect their oviposition and movement responses?, PLoS One. 11 (2016) 1–12. doi:10.1371/journal.pone.0149994.

[28]    D.G. Lilly, S.L. Latham, C.E. Webb, S.L. Doggett, Cuticle thickening in a pyrethroid-resistant strain of the common bed bug, cimex lectularius l. (hemiptera: Cimicidae), PLoS One. 11 (2016) 6–16. doi:10.1371/journal.pone.0153302.

[29]    Uganda Bureau of Statistics, National Population and Housing Census 2014, Uganda Bur. Stat. (2016) 1–105. doi:10.1017/CBO9781107415324.004.

[30]    World Vision, Addressing food security and household income challenges, World Vis. (2018) 2018.

[31]    FAO Forestry Division, Forests , livelihoods and poverty alleviation : the case of Uganda, FAO. (2013) 1–87.

[32]    P.O. Otim, Baseline Study for the Ugandan side of the Karamoja Cluster, CEWARN/ICAD, Addis Abeba Ethiop. (2004) 1–29.

[33]    Onor E, Malaria in Karamoja District, Uganda, World Heal. Organ. (1969) 1–20.

[34]    National Malaria Control Programme, Abt Associates, the INFORM Project, An epidemiological profile of malaria and its control in Uganda, Natl. Malar. Control Program. Program. (2013).

[35]    UBOS, Uganda Demograpic and Health Survey 2011, Uganda Bur. Stat. (2012) 1–461.

[36]    C.M.A. Gillies MT, A Supplement to the Anophelinae of Africa South of the Sahara (Ethiopian zoogeographical region), South African Inst. Med. Res. (1987) 1–146.

[37]    WHO/ GMP, Test procedures for insecticide resistance monitoring in malaria vector mosquitoes: Second edition, World Heal. Organ. Tech. Rep. Ser. (2016) 1–54. https://www.who.int/malaria/publications/atoz/9789241511575/en/.

[38]    G.E. Ford H.R, Laboratory rearing of Anopheles albimanus Wiedemann, Mosq. News. 32 (1972) 509–513. doi:10.1017/CBO9781107415324.004.

[39]    N.I. of H. ATCC, Bei Resources, Methods in Anopheles Research, MR4. (2015).

[40]    J. Hemingway, N.J. Hawkes, L. McCarroll, H. Ranson, The molecular basis of insecticide resistance in mosquitoes, Insect Biochem. Mol. Biol. 34 (2004) 653–665. doi:10.1016/j.ibmb.2004.03.018.

[41]    H.A.A. Khan, W. Akram, S.A. Shad, J.J. Lee, Insecticide Mixtures Could Enhance the Toxicity of Insecticides in a Resistant Dairy Population of Musca domestica L, PLoS One. 8 (2013) 1–8. doi:10.1371/journal.pone.0060929.

[42]    J.A. Scott, W.G. Brogdon, F.H. Collins, Identification of single specimens of the Anopheles gambiae complex by the polymerase chain reaction, Am. J. Trop. Med. Hyg. 49 (1993) 520–529. doi:10.4269/ajtmh.1993.49.520.

[43]    L.L. Koekemoer, L. Kamau, R.H. Hunt, M. Coetzee, A cocktail polymerase chain reaction assay to identify members of the Anopheles funestus (Diptera: Culicidae) group, Am. J. Trop. Med. Hyg. 66 (2002) 804–811. doi:10.4269/ajtmh.2002.66.804.

[44]    J.R. Busvine, Critical review of the techniques for testing insecticides, Commonw. Inst. Entomol. London. (1958).

[45]    S. Horstmann, R. Sonneck, Contact bioassays with phenoxybenzyl and tetrafluorobenzyl pyrethroids against target- site and metabolic resistant mosquitoes, PLoS One. 11 (2016) 1–16. doi:10.1371/journal.pone.0149738.

[46]    P.E. Okello, W. Van Bortel, A.M. Byaruhanga, A. Correwyn, P. Roelants, A. Talisuna, U. D’Alessandro, M. Coosemans, Variation in malaria transmission intensity in seven sites throughout Uganda, Am. J. Trop. Med. Hyg. 75 (2006) 219–225. doi:10.4269/ajtmh.2006.75.219.

[47]    R. Echodu, J.J. Lutwama, J. Enyaru, R. Ocan, R.B. Asaba, F. Ajuga, D. Bradley, C. Mutero, C. Kabonesa, J. Olobo, Heterogeneity of Anopheles Mosquitoes in Nyabushozi County , Kiruhura district , Uganda, J. Parasitol. Vector Biol. 2 (2010) 28–34.

[48]    A. Lynd, S. Gonahasa, S.G. Staedke, A. Oruni, C. Maiteki-Sebuguzi, G. Dorsey, J. Opigo, A. Yeka, A. Katureebe, M. Kyohere, J. Hemingway, M.R. Kamya, M.J. Donnelly, LLIN Evaluation in Uganda Project (LLINEUP): A cross-sectional survey of species diversity and insecticide resistance in 48 districts of Uganda, Parasites and Vectors. 12 (2019) 1–10. doi:10.1186/s13071-019-3353-7.

[49]    J. Kitau, R.M. Oxborough, P.K. Tungu, J. Matowo, R.C. Malima, S.M. Magesa, J. Bruce, F.W. Mosha, M.W. Rowland, Species shifts in the anopheles gambiae complex: Do LLINs successfully control anopheles arabiensis?, PLoS One. 7 (2012) 1–7. doi:10.1371/journal.pone.0031481.

[50]    R.S. McCann, E. Ochomo, M.N. Bayoh, J.M. Vulule, M.J. Hamel, J.E. Gimnig, W.A. Hawley, E.D. Walker, Reemergence of Anopheles funestus as a vector of Plasmodium falciparum in Western Kenya after long-term implementation of insecticide-treated bed nets, Am. J. Trop. Med. Hyg. 90 (2014) 597–604. doi:10.4269/ajtmh.13-0614.

[51]    O. Kenea, M. Balkew, H. Tekie, W. Deressa, E. Loha, B. Lindtjørn, H.J. Overgaard, Impact of combining indoor residual spraying and long-lasting insecticidal nets on Anopheles arabiensis in Ethiopia: Results from a cluster randomized controlled trial, Malar. J. 18 (2019) 1–11. doi:10.1186/s12936-019-2811-1.

[52]    B.E. Yohannes M, Early biting rhythm in the afro-tropical vector of malaria, Anopheles arabiensis, and challenges for its control in Ethiopia, Med. Vet. Entomol. 26 (2011) 103–105.

[53]    B. Abong’o, J.E. Gimnig, S.J. Torr, B. Longman, D. Omoke, M. Muchoki, F. ter Kuile, E. Ochomo, S. Munga, A.M. Samuels, K. Njagi, J. Maas, R.T. Perry, C. Fornadel, M.J. Donnelly, R.M. Oxborough, Impact of indoor residual spraying with pirimiphos-methyl (Actellic 300CS) on entomological indicators of transmission and malaria case burden in Migori County, western Kenya, Sci. Rep. 10 (2020) 1–14. doi:10.1038/s41598-020-61350-2.

[54]    F.M. Nafisah Braimah, Chris Drakeley, Eliningaya Kweka, Tests of bednet traps (Mbita traps) for monitoring mosquito populations and time of biting in Tanzania and possible impact of prolonged insecticide treated net use, Int. Jounral Trop. Insect Sci. 25 (2007) 208–213. doi:https://doi.org/10.1079/IJT200576.

[55]    T.L. Russell, D.W. Lwetoijera, D. Maliti, B. Chipwaza, J. Kihonda, J.D. Charlwood, T.A. Smith, C. Lengeler, M.A. Mwanyangala, R. Nathan, B.G. Knols, W. Takken, G.F. Killeen, Impact of promoting longer-lasting insecticide treatment of bed nets upon malaria transmission in a rural Tanzanian setting with pre-existing high coverage of untreated nets, BMC Malar. J. 9 (2010) 1–14. doi:10.1186/1475-2875-9-187.

[56]    R. N’Guessan, V. Corbel, M. Akogbéto, M. Rowland, Reduced efficacy of insecticide-treated nets and indoor residual spraying for malaria control in pyrethroid resistance area, Benin, Emerg. Infect. Dis. 13 (2007) 199–206. doi:10.3201/eid1302.060631.

[57]    T.A. Abeku, M.E.H. Helinski, M.J. Kirby, J. Ssekitooleko, C. Bass, I. Kyomuhangi, M. Okia, G. Magumba, S.R. Meek, Insecticide resistance patterns in Uganda and the effect of indoor residual spraying with bendiocarb on kdr L1014S frequencies in Anopheles gambiae s.s., Malar. J. 16 (2017) 1–11. doi:10.1186/s12936-017-1799-7.

[58]    R. Maharaj, Global trends in insecticide resistance and impact on disease vector control measures, Open Access Insect Physiol. 3 (2011) 27–33.

[59]    W.G. Brogdon, J.C. McAllister, Insecticide resistance and vector control, Emerg. Infect. Dis. 4 (1998) 605–613. doi:10.1300/J096v09n02_21.

[60]    R. Echodu, J. Iga, W.S. Oyet, P. Mireji, J. Anena, D. Onanyang, T. Iwiru, J.J. Lutwama, E.A. Opiyo, High insecticide resistances levels in Anopheles gambiaes s.l. In northern Uganda and its relevance for future malaria control, BMC Res. Notes. 13 (2020) 1–6. doi:10.1186/s13104-020-05193-0.

[61]    M. Okia, R. Ndyomugyenyi, J. Kirunda, A. Byaruhanga, S. Adibaku, D.K. Lwamafa, F. Kironde, Bioefficacy of long-lasting insecticidal nets against pyrethroid-resistant populations of Anopheles gambiae s.s. from different malaria transmission zones in Uganda, Parasites and Vectors. 6 (2013) 1–10. doi:10.1186/1756-3305-6-130.

[62]    M. Okia, D.F. Hoel, J. Kirunda, J.B. Rwakimari, B. Mpeka, D. Ambayo, A. Price, D.W. Oguttu, A.P. Okui, J. Govere, Insecticide resistance status of the malaria mosquitoes: Anopheles gambiae and Anopheles funestus in eastern and northern Uganda, BMC Malar. J. 17 (2018) 1–12. doi:10.1186/s12936-018-2293-6.

[63]    R. John, T. Ephraim, A. Andrew, Reduced susceptibility to pyrethroid insecticide treated nets by the malaria vector Anopheles gambiae s.l. in western Uganda, BMC Malar. JournalMalaria J. 7 (2008) 1–5. doi:10.1186/1475-2875-7-92.

[64]    U. Ramphul, T. Boase, C. Bass, L.M. Okedi, M.J. Donnelly, P. Müller, Insecticide resistance and its association with target-site mutations in natural populations of Anopheles gambiae from eastern Uganda, Trans. R. Soc. Trop. Med. Hyg. 103 (2009) 1121–1126. doi:10.1016/j.trstmh.2009.02.014.

[65]    and M.J.D. Henry D. Mawejje, Craig S. Wilding, Emily J. Rippon, Angela Hughes, David Weetman, Insecticide resistance monitoring of field-collected Anopheles gambiae s.l. populations from Jinja, eastern Uganda, identifies high levels of pyrethroid resistance, Med. Vet. Entomol. 23 (2013) 276–283. doi:10.1038/jid.2014.371.

[66]    K. Verhaeghen, W. Van Bortel, P. Roelants, P.E. Okello, A. Talisuna, M. Coosemans, Spatio-temporal patterns in kdr frequency in permethrin and DDT resistant Anopheles gambiae s.s. from Uganda, Am. J. Trop. Med. Hyg. 82 (2010) 566–573. doi:10.4269/ajtmh.2010.08-0668.

[67]    P.Y. Siegert, E. Walker, J.R. Miller, D. Behavioral, Differential Behavioral Responses of Anopheles gambiae (Diptera: Culicidae) Modulate Mortality Caused by Pyrethroid-Treated Bednets, J. Econ. Entomol. 102 (2009) 2061–2071.

[68]    C.P. Malal M. Diop, Nicolas Moiroux, Fabrice Chandre, Hadrien Martin-Herrou, Pascal Milesi, Olayidé Boussari, Angélique Porciani, Stéphane Duchon, Pierrick Labbé, Behavioral cost & overdominance in anopheles gambiae, PLoS One. 10 (2015) 1–12. doi:10.1371/journal.pone.0121755.

[69]    N. Platt, R.M. Kwiatkowska, H. Irving, A. Diabaté, R. Dabire, C.S. Wondji, Target-site resistance mutations (kdr and RDL), but not metabolic resistance, negatively impact male mating competiveness in the malaria vector Anopheles gambiae, Heredity (Edinb). 115 (2015) 243–252. doi:10.1038/hdy.2015.33.

[70]    H. Alout, R.K. Dabiré, L.S. Djogbénou, L. Abate, V. Corbel, F. Chandre, A. Cohuet, Interactive cost of Plasmodium infection and insecticide resistance in the malaria vector Anopheles gambiae, Sci. Rep. 6 (2016) 1–11. doi:10.1038/srep29755.

[71]    E. Collins, N.M. Vaselli, M. Sylla, A.H. Beavogui, J. Orsborne, G. Lawrence, R.E. Wiegand, S.R. Irish, T. Walker, L.A. Messenger, The relationship between insecticide resistance, mosquito age and malaria prevalence in Anopheles gambiae s.l. from Guinea, Sci. Rep. 9 (2019) 1–12. doi:10.1038/s41598-019-45261-5.

[72]    M. Viana, A. Hughes, J. Matthiopoulos, H. Ranson, H.M. Ferguson, Delayed mortality effects cut the malaria transmission potential of insecticide-resistant mosquitoes, Proc. Natl. Acad. Sci. U. S. A. 113 (2016) 8975–8980. doi:10.1073/pnas.1603431113.

[73]    M. Kristan, J. Lines, A. Nuwa, C. Ntege, S.R. Meek, T.A. Abeku, Exposure to deltamethrin affects development of Plasmodium falciparum inside wild pyrethroid resistant Anopheles gambiae s.s. mosquitoes in Uganda, Parasites and Vectors. 9 (2016) 1–9. doi:10.1186/s13071-016-1384-x.

[74]    H. Alout, N.T. Ndam, M.M. Sandeu, I. Djégbe, F. Chandre, R.K. Dabiré, L.S. Djogbénou, V. Corbel, A. Cohuet, Insecticide Resistance Alleles Affect Vector Competence of Anopheles gambiae s.s. for Plasmodium falciparum Field Isolates, PLoS One. 8 (2013) 1–9. doi:10.1371/journal.pone.0063849.

[75]    M. Coetzee, L.L. Koekemoer, Molecular Systematics and Insecticide Resistance in the Major African Malaria Vector Anopheles funestus, Annu. Rev. Entomol. 58 (2013) 393–412. doi:10.1146/annurev-ento-120811-153628.

[76]    C. Strode, S. Donegan, P. Garner, A.A. Enayati, J. Hemingway, The Impact of Pyrethroid Resistance on the Efficacy of Insecticide-Treated Bed Nets against African Anopheline Mosquitoes: Systematic Review and Meta-Analysis, PLoS Med. 11 (2014) 1–32. doi:10.1371/journal.pmed.1001619.