Investigating molecular basis of lambda-cyhalothrin resistance in an Anopheles funestus population from Senegal

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

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

Recombinant expression and characterization of GSTd3 from a resistant population of Anopheles arabiensis and comparison of DDTase activity with GSTe2

The development of insecticide resistance in malaria vectors is a challenge for the global effort to control and eradicate malaria. Glutathione S-transferases (GSTs) are multifunctional enzymes involved in the detoxification of many classes of insecticides. For mosquitoes, it is known that overexpression of an epsilon GST, GSTe2, confers resistance towards DDT and pyrethroids. In addition to GSTe2, consistent overexpression of a delta class GST, GSTd3, has been observed in insecticide resistant populations of different malaria vector species. However, the functional role of GSTd3 towards DDT resistance has not yet been investigated. Here, we recombinantly expressed both GSTe2 and GSTd3 from Anopheles arabiensis and compared their metabolic activities against DDT. Both AaGSTd3 and AaGSTe2 exhibited CDNB-conjugating and glutathione peroxidase activity and DDT metabolism was observed for both GSTs. However, the DDT dehydrochlorinase activity exhibited by AaGSTe2 was much higher than for AaGSTd3, and AaGSTe2 was also able to eliminate DDE although the metabolite could not be identified. Molecular modeling revealed subtle differences in the binding pocket of both enzymes and a better fit of DDT within the H-site of AaGSTe2. The overexpression but much lower DDT metabolic activity of AaGSTd3, might suggest that AaGSTd3 sequesters DDT. These findings highlight the complexity of insecticide resistance in the major malaria vectors and the difficulties associated with control of the vectors using DDT, which is still used for indoor residual spraying.

Genome-wide Association Study Reveals New Loci Associated With Pyrethroid Resistance in Aedes aegypti

Genome-wide association studies (GWAS) use genetic polymorphism across the genomes of individuals with distinct characteristics to identify genotype-phenotype associations. In mosquitoes, complex traits such as vector competence and insecticide resistance could benefit from GWAS. We used the Aedes aegypti 50k SNP chip to genotype populations with different levels of pyrethroid resistance from Northern Brazil. Pyrethroids are widely used worldwide to control mosquitoes and agricultural pests, and their intensive use led to the selection of resistance phenotypes in many insects including mosquitoes. For Ae. aegypti, resistance phenotypes are mainly associated with several mutations in the voltage-gated sodium channel, known as knockdown resistance (kdr). We phenotyped those populations with the WHO insecticide bioassay using deltamethrin impregnated papers, genotyped the kdr alleles using qPCR, and determined allele frequencies across the genome using the SNP chip. We identified single-nucleotide polymorphisms (SNPs) directly associated with resistance and one epistatic SNP pair. We also observed that the novel SNPs correlated with the known kdr genotypes, although on different chromosomes or not in close physical proximity to the voltage gated sodium channel gene. In addition, pairwise comparison of resistance and susceptible mosquitoes from each population revealed differentiated genomic regions not associated with pyrethroid resistance. These new bi-allelic markers can be used to genotype other populations along with kdr alleles to understand their worldwide distribution. The functional roles of the genes near the newly discovered SNPs require new studies to determine if they act synergistically with kdr alleles or reduce the fitness cost of maintaining resistant alleles.

RNAseq-based gene expression profiling of the Anopheles funestus pyrethroid-resistant strain FUMOZ highlights the predominant role of the duplicated CYP6P9a/b cytochrome P450s

Insecticide-based interventions, notably long-lasting insecticidal nets, against mosquito vectors of malaria are currently threatened by pyrethroid resistance. Here, we contrasted RNAseq-based gene expression profiling of laboratory-resistant (FUMOZ) and susceptible (FANG) strains of the major malaria vector Anopheles funestus. Cytochrome P450 genes were the predominant over-expressed detoxification genes in FUMOZ, with high expression of the duplicated CYP6P9a (fold-change of 82.23 vs FANG) and CYP6P9b (FC 11.15). Other over-expressed P450s belonged to the same cluster of P450s corresponding to the resistance to pyrethroid 1 (rp1) quantitative trait loci (QTL) on chromosome 2R. Several Epsilon class glutathione S-transferases were also over-expressed in FUMOZ, as was the ATP-binding cassette transporter AFUN019220 (ABCA) which also exhibited between-strain alternative splicing events at exon 7. Significant differences in single-nucleotide polymorphism frequencies between strains occurred in resistance QTLs rp1 (CYP6P9a/b and CYP6AA1), rp2 on chromosome 2L (CYP6Z1, CYP6M7, and CYP6Z3), and rp3 on chromosome 3R (CYP9J5, CYP9J4, and CYP9J3). Differences were also detected in CYP4G17 and CYP4G16 genes on the X chromosome, both of which are associated with cuticular resistance in Anopheles gambiae. A close analysis of nonsynonymous diversity at the CYP6P9a/b loci revealed a drastic loss of diversity in FUMOZ with only a single polymorphism and 2 haplotypes vs 18 substitutions and 8 haplotypes in FANG. By contrast, a lowly expressed cytochrome P450 (CYP4C36) did not exhibit diversity differences between strains. We also detected the known pyrethroid resistance conferring amino acid change N384S in CYP6P9b. This study further elucidates the molecular bases of resistance in An. funestus, informing strategies to better manage widespread resistance across Africa.

Molecular surveillance of resistance to pyrethroids insecticides in Colombian Aedes aegypti populations

Introduction

In Colombia, organochloride, organophosphate, carbamate, and pyrethroid insecticides are broadly used to control Aedes aegypti populations. However, Colombian mosquito populations have shown variability in their susceptibility profiles to these insecticides, with some expressing high resistance levels.

Materials and methods

In this study, we analyzed the susceptibility status of ten Colombian field populations of Ae. aegypti to two pyrethroids; permethrin (type-I pyrethroid) and lambda-cyhalothrin (type-II pyrethroid). In addition, we evaluated if mosquitoes pressured with increasing lambda-cyhalothrin concentrations during some filial generations exhibited altered allelic frequency of these kdr mutations and the activity levels of some metabolic enzymes.

Results

Mosquitoes from all field populations showed resistance to lambda-cyhalothrin and permethrin. We found that resistance profiles could only be partially explained by kdr mutations and altered enzymatic activities such as esterases and mixed-function oxidases, indicating that other yet unknown mechanisms could be involved. The molecular and biochemical analyses of the most pyrethroid-resistant mosquito population (Acacías) indicated that kdr mutations and altered metabolic enzyme activity are involved in the resistance phenotype expression.

Conclusions

In this context, we propose genetic surveillance of the mosquito populations to monitor the emergence of resistance as an excellent initiative to improve mosquito-borne disease control measures.

The main method of preventing Aedes-borne diseases such as dengue, Zika, and chikungunya is by targeting the primary mosquito vector, Aedes aegypti, with insecticides. However, the success of these vector control strategies is jeopardized by the widespread development of insecticide resistance in mosquito populations. Furthermore, the molecular mechanisms of insecticide resistance in Ae. aegypti are still not well understood, resulting in limited resistance mitigation and management strategies. In this paper, we found that resistance to some pyrethroid insecticides in different Colombian cities is associated with three allelic substitutions V419L, V1016I, and F1534C, on the voltage-gated sodium channel gene, known as kdr (‘knock-down resistance’) mutations, with all three mutations present in mosquitoes resistant to pyrethroids. The data also showed that kdr mutations are important in conferring low resistance levels, but after around 10-fold intensity, the allele frequencies don’t change, indicating that other mechanisms contribute to the resistance. Thus, we found that mosquitoes under selective pressure with insecticides present also altered enzymatic activities such as esterases and mixed-function oxidases, indicating that kdr mutations and metabolic enzymes are involved in the resistance expression. The findings on the extent of insecticide resistance and the molecular mechanisms underpinning the problem will impact the surveillance, selection, and rational use of insecticides by local health authorities.

From Global to Local-New Insights into Features of Pyrethroid Detoxification in Vector Mosquitoes

The threat of mosquito-borne diseases continues to be a problem for public health in subtropical and tropical regions of the world; in response, there has been increased use of adulticidal insecticides, such as pyrethroids, in human habitation areas over the last thirty years. As a result, the prevalence of pyrethroid-resistant genetic markers in natural mosquito populations has increased at an alarming rate. This review details recent advances in the understanding of specific mechanisms associated with pyrethroid resistance, with emphasis on features of insecticide detoxification and the interdependence of multiple cellular pathways. Together, these advances add important context to the understanding of the processes that are selected in resistant mosquitoes. Specifically, before pyrethroids bind to their targets on motoneurons, they must first permeate the outer cuticle and diffuse to inner tissues. Resistant mosquitoes have evolved detoxification mechanisms that rely on cytochrome P450s (CYP), esterases, carboxyesterases, and other oxidation/reduction (redox) components to effectively detoxify pyrethroids to nontoxic breakdown products that are then excreted. Enhanced resistance mechanisms have evolved to include alteration of gene copy number, transcriptional and post-transcriptional regulation of gene expression, as well as changes to cellular signaling mechanisms. Here, we outline the variety of ways in which detoxification has been selected in various mosquito populations, as well as key gene categories involved. Pathways associated with potential new genes of interest are proposed. Consideration of multiple cellular pathways could provide opportunities for development of new insecticides.

Investigation of DDT resistance mechanisms in Anopheles funestus populations from northern and southern Benin reveals a key role of the GSTe2 gene

Background

Understanding the molecular basis of insecticide resistance in mosquito, such as Anopheles funestus, is an important step in developing strategies to mitigate the resistance problem. This study aims to assess the role of the GSTe2 gene in DDT resistance and determine the genetic diversity of this gene in An. funestus.

Methods

Gene expression analysis was performed using microarrays and PCR while the potential mutation associated with resistance was determined using sequencing.

Results

Low expression level of GSTe2 gene was recorded in Burkina-Faso samples with a fold change of 3.3 while high expression (FC 35.6) was recorded in southern Benin in Pahou (FC 35.6) and Kpome (FC 13.3). The sequencing of GSTe2 gene in six localities showed that L119F-GSTe2 mutation is almost getting fixed in highly DDT-resistant Benin (Pahou, Kpome, Doukonta) and Nigeria (Akaka Remo) mosquitoes with a low mutation rate observed in Tanongou (Benin) and Burkina-Faso mosquitoes.

Conclusion

This study shows the key role of the GSTe2 gene in DDT resistant An. funestus in Benin. Polymorphism analysis of this gene across Benin revealed possible barriers to gene flow, which could impact the design and implementation of resistance management strategies in the country.

Mugenzi LM, Tripet F, Konaté L, Diabate A, Dia I, Gaye O, Faye O, Niang EHA, S. Wondji C. Insecticide Resistance Profiling of Anopheles coluzzii and Anopheles gambiae Populations in the Southern Senegal: Role of Target Sites and Metabolic Resistance Mechanisms

The emergence and spread of insecticide resistance among the main malaria vectors is threatening the effectiveness of vector control interventions in Senegal. The main drivers of this resistance in the Anopheles gambiae complex (e.g., An. gambiae and Anopheles coluzzii) remains poorly characterized in Senegal. Here we characterized the main target site and metabolic resistances mechanisms among the An. gambiae and An. coluzzii populations from their sympatric and allopatric or predominance area in Senegal. Larvae and pupae of An. gambiae s.l. were collected, reared to adulthood, and then used for insecticides susceptibility and synergist assays using the WHO (World Health Organisation) test kits for adult mosquitoes. The TaqMan method was used for the molecular characterization of the main target site insecticide resistance mechanisms (Vgsc-1014F, Vgsc-1014S, N1575Y and G119S). A RT-qPCR (Reverse Transcriptase-quantitative Polymerase Chaine Reaction) was performed to estimate the level of genes expression belonging to the CYP450 (Cytochrome P450) family. Plasmodium infection rate was investigated using TaqMan method. High levels of resistance to pyrethroids and DDT and full susceptibility to organophosphates and carbamates where observed in all three sites, excepted a probable resistance to bendiocarb in Kedougou. The L1014F, L1014S, and N1575Y mutations were found in both species. Pre-exposure to the PBO (Piperonyl butoxide) synergist induced a partial recovery of susceptibility to permethrin and full recovery to deltamethrin. Subsequent analysis of the level of genes expression, revealed that the CYP6Z1 and CYP6Z2 genes were over-expressed in wild-resistant mosquitoes compared to the reference susceptible strain (Kisumu), suggesting that both the metabolic resistance and target site mutation involving kdr mutations are likely implicated in this pyrethroid resistance. The presence of both target-site and metabolic resistance mechanisms in highly pyrethroid-resistant populations of An. gambiae s.l. from Senegal threatens the effectiveness and the sustainability of the pyrethroid-based tools and interventions currently deployed in the country. The Kdr-west mutation is widely widespread in An. coluzzii sympatric population. PBO or Duo nets and IRS (Indoor Residual Spraying) with organophosphates could be used as an alternative measure to sustain malaria control in the study area.

Influence of GST- and P450-based metabolic resistance to pyrethroids on blood feeding in the major African malaria vector Anopheles funestus

Insecticide resistance genes are often associated with pleiotropic effects on various mosquito life-history traits. However, very little information is available on the impact of insecticide resistance on blood feeding process in mosquitoes. Here, using two recently detected DNA-based metabolic markers in the major malaria vector, An. funestus, we investigated how metabolic resistance genes could affect the blood meal intake. After allowing both the field F1 and lab F8 Anopheles funestus strains to feed on the human arm for 30 minutes, we assessed the association between key parameters of blood meal process including, probing time, feeding duration, blood feeding success, blood meal size, and markers of glutathione S-transferase (L119F-GSTe2) and cytochrome P450 (CYP6P9a_R)-mediated metabolic resistance. None of the parameters of blood meal process was associated with L119F-GSTe2 genotypes. By contrast, for CYP6P9a_R, homozygous resistant mosquitoes were significantly more able to blood-feed than homozygous susceptible (OR = 3.3; CI 95%: 1.4-7.7; P = 0.01) mosquitoes. Moreover, the volume of blood meal ingested by CYP6P9a-SS mosquitoes was lower than that of CYP6P9a-RS (P<0.004) and of CYP6P9a-RR (P<0.006). This suggests that CYP6P9a gene is inked with the feeding success and blood meal size of An. funestus. However, no correlation was found in the expression of CYP6P9a and that of genes encoding for salivary proteins involved in blood meal process. This study suggests that P450-based metabolic resistance may influence the blood feeding process of Anopheles funestus mosquito and consequently its ability to transmit malaria parasites.

Investigating the molecular basis of multiple insecticide resistance in a major malaria vector Anopheles funestus (sensu stricto) from Akaka-Remo, Ogun State, Nigeria

Background

Understanding the mechanisms used by Anopheles mosquitoes to survive insecticide exposure is key to manage existing insecticide resistance and develop more suitable insecticide-based malaria vector control interventions as well as other alternative integrated tools. To this regard, the molecular basis of permethrin, DDT and dieldrin resistance in Anopheles funestus (sensu stricto) at Akaka-Remo was investigated.

Methods

Bioassays were conducted on 3–5-day-old adult An. funestus (s.s.) mosquitoes for permethrin, DDT and dieldrin susceptibility test. The molecular mechanisms of mosquito resistance to these insecticides were investigated using microarray and reverse transcriptase PCR techniques. The voltage-gated sodium channel region of mosquitoes was also screened for the presence of knockdown resistance mutations (kdr west and east) by sequencing method.

Results

Anopheles funestus (s.s.) population was resistant to permethrin (mortality rate of 68%), DDT (mortality rate of 10%) and dieldrin (mortality rate of 8%) insecticides. Microarray and RT-PCR analyses revealed the overexpression of glutathione S-transferase genes, cytochrome P450s, esterase, trypsin and cuticle proteins in resistant mosquitoes compared to control. The GSTe2 was the most upregulated detoxification gene in permethrin-resistant (FC = 44.89), DDT-resistant (FC = 57.39) and dieldrin-resistant (FC = 41.10) mosquitoes compared to control population (FC = 22.34). The cytochrome P450 gene, CYP6P9b was also upregulated in both permethrin- and DDT-resistant mosquitoes. The digestive enzyme, trypsin (hydrolytic processes) and the cuticle proteins (inducing cuticle thickening leading to reduced insecticides penetration) also showed high involvement in insecticide resistance, through their overexpression in resistant mosquitoes compared to control. The kdr east and west were absent in all mosquitoes analysed, suggesting their non-involvement in the observed mosquito resistance.

Conclusions

The upregulation of metabolic genes, especially the GSTe2 and trypsin, as well as the cuticle proteins is driving insecticide resistance of An. funestus (s.s.) population. However, additional molecular analyses, including functional metabolic assays of these genes as well as screening for a possible higher cuticular hydrocarbon and lipid contents, and increased procuticle thickness in resistant mosquitoes are needed to further describe their distinct roles in mosquito resistance.

Anopheles arabiensis and Anopheles funestus biting patterns in Dielmo, an area of low level exposure to malaria vectors

Background

In Dielmo, Senegal, the widespread use of long-lasting insecticidal nets has decreased both the incidence of malaria and the density of the Anopheles population. However, persistent low-level malaria transmission may hamper efforts to eliminate the disease. Therefore, continuous monitoring of the vector population is needed in order to improve knowledge of Anopheles biting behaviour and to readjust control interventions.

Methods

In 2015, Anopheles were collected every month for a whole year and each specimen was identified using morphological and molecular techniques. The biting pattern of each species was analysed according to night (7 pm–7am) and morning (7am–11am) periods, the place of biting and the season. The ELISA CSP technique was used to assess the Plasmodium falciparum sporozoite rate to evaluate the entomological inoculation rate (EIR).

Results

Anopheles arabiensis and Anopheles funestus sensu stricto were found to be the main vectors biting humans. Overall, the biting rate was low, at 3.84bites per night (bpn) and 1.27 bites per morning (bpm), respectively (IRR = 3.04, CI [1.84–5.00], p < 0.001). The EIR was 2.51 and 5.03 infectious bites per year during the night and morning, respectively. During the night, the An. arabiensis and An. funestus biting rate was 1.81 bpn and 1.71 bpn, respectively (IRR = 0.95, CI [0.46–1.92], p = 0.88). During the morning, their density decreased to 0.51 bpm and 0.73 bpm for An. arabiensis and An. funestus, respectively (IRR = 1.47, CI [0.58–3.71], p = 0.41). During the night and the morning, no specific trend of indoor or outdoor biting was observed in the dry and rainy season for both vectors.

Conclusion

This study highlighted low level Anopheles nocturnal and diurnal biting and the associated risk of malaria transmission. It showed also the influence of the season on the indoor and outdoor biting pattern, indicating that the human population could be exposed all year round to a low level of Anopheles bites. Control programmes should increase awareness of the use of bed nets throughout the year and promote the development and implementation of complimentary tools to target Anopheles biting shortly after dawn when people are still indoors and outside the bed nets.

Exploring the Mechanisms of Multiple Insecticide Resistance in a Highly Plasmodium-Infected Malaria Vector Anopheles funestus Sensu Stricto from Sahel of Northern Nigeria

The Nigerian Government is scaling up the distribution of insecticide-treated bed nets for malaria control, but the lack of surveillance data, especially in the Sudan/Sahel region of the country, may hinder targeting priority populations. Here, the vectorial role and insecticide resistance profile of a population of a major malaria vector Anopheles funestus sensu stricto from Sahel of Nigeria was characterised. An. funestus s.s. was the only vector found, with a high human blood index (100%) and a biting rate of 5.3/person/night. High Plasmodium falciparum infection was discovered (sporozoite rate = 54.55%). The population is resistant to permethrin (mortality = 48.30%, LT₅₀ = 65.76 min), deltamethrin, DDT (dichlorodiphenyltrichloroethane) and bendiocarb, with mortalities of 29.44%, 56.34% and 54.05%, respectively. Cone-bioassays established loss of efficacy of the pyrethroid-only long-lasting insecticidal nets (LLINs); but 100% recovery of susceptibility was obtained for piperonylbutoxide (PBO)-containing PermaNet^®3.0. Synergist bioassays with PBO and diethyl maleate recovered susceptibility, implicating CYP450s (permethrin mortality = 78.73%, χ² = 22.33, P < 0.0001) and GSTs (DDT mortality = 81.44%, χ² = 19.12, P < 0.0001). A high frequency of 119F GSTe2 mutation (0.84) was observed (OR = 16, χ² = 3.40, P = 0.05), suggesting the preeminent role of metabolic resistance. These findings highlight challenges associated with deployment of LLINs and indoor residual spraying (IRS) in Nigeria.

Insecticide Resistance Status and Mechanisms of Anopheles sinensis (Diptera: Culicidae) in Wenzhou, an Important Coastal Port City in China

Although scaled-up interventions and effective control efforts have drastically reduced malaria morbidity and mortality, malaria remains a serious threat to public health worldwide. Anopheles sinensis Wiedemann 1828 is a historically important vector of Plasmodium vivax (Haemosporida: Plasmodiidae) malaria in China. Insecticide resistance has become a major obstacle to vector-borne disease control. However, little is known about the insecticide resistance of An. sinensis in Wenzhou, an important coastal port city in Zhejiang province, China. The aim of this study was to examine insecticide resistance and mechanisms in An. sinensis field mosquito populations. Evidence of multiple insecticide resistance was found in An. sinensis adult female populations. Medium to high frequencies of target site kdr together with fixed ace-1 mutations was detected in both the Ruian and Yongjia populations. Both populations showed an association between kdr L1014 mutation and resistance phenotype when tested against deltamethrin and DDT. Significantly different metabolic enzyme activities were found between the susceptible laboratory strain and field-collected mosquitoes from both Ruian and Yongjia. Both field collected An. sinensis populations exhibited significantly higher P450 enzyme activity compared with the laboratory strain, while the field-collected resistant mosquitoes exhibited various GST and COE enzyme activities. These results indicate multiple resistance mechanisms in An. sinensis field populations. Effective implementation of insecticide resistance management strategies is urgently needed. The data collected in this study will be valuable for modeling insecticide resistance spread and vector-control interventions.

A marker of glutathione S-transferase-mediated resistance to insecticides is associated with higher Plasmodium infection in the African malaria vector Anopheles funestus

Metabolic resistance to insecticides is threatening malaria control in Africa. However, the extent to which it impacts malaria transmission remains unclear. Here, we investigated the association between a marker of glutathione S-transferase mediated metabolic resistance and Plasmodium infection in field population of Anopheles funestus s.s. in comparison to the A296S-RDL target site mutation. The 119F-GSTe2 resistant allele was present in southern (Obout) (56%) and central (Mibellon) (25%) regions of Cameroon whereas the 296S-RDL resistant allele was detected at 98.5% and 15% respectively. The whole mosquito Plasmodium and sporozoite infection rates were 57% and 14.8% respectively in Obout (n = 508) and 19.7% and 5% in Mibellon (n = 360). No association was found between L119F-GSTe2 genotypes and whole mosquito infection status. However, when analyzing oocyst and sporozoite infection rates separately, the resistant homozygote 119F/F genotype was significantly more associated with Plasmodium infection in Obout than both heterozygote (OR = 2.5; P = 0.012) and homozygote susceptible (L/L119) genotypes (OR = 2.10; P = 0.013). In contrast, homozygote RDL susceptible mosquitoes (A/A296) were associated more frequently with Plasmodium infection than other genotypes (OR = 4; P = 0.03). No additive interaction was found between L119F and A296S. Sequencing of the GSTe2 gene showed no association between the polymorphism of this gene and Plasmodium infection. Glutathione S-transferase metabolic resistance is potentially increasing the vectorial capacity of resistant An. funestus mosquitoes. This could result in a possible exacerbation of malaria transmission in areas of high GSTe2-based metabolic resistance to insecticides.

Overexpression of Two Members of D7 Salivary Genes Family is Associated with Pyrethroid Resistance in the Malaria Vector Anopheles Funestus s.s. but Not in Anopheles Gambiae in Cameroon

D7 family proteins are among the most expressed salivary proteins in mosquitoes. They facilitate blood meal intake of the mosquito by scavenging host amines that induce vasoconstriction, platelet aggregation and pain. Despite this important role, little information is available on the impact of insecticide resistance on the regulation of D7 proteins and consequently on the blood feeding success. In this study, real-time quantitative polymerase chain reaction (qPCR) analyses were performed to investigate how pyrethroid resistance could influence the expression of genes encoding D7 family proteins in Anopheles gambiae and Anopheles funestus s.s. mosquitoes from Elon in the Central Cameroon. Out of 328 collected mosquitoes, 256 were identified as An. funestus sl and 64 as An. gambiae sl. Within the An. funestus group, An. funestus s.s. was the most abundant species (95.95%) with An. rivulorum, An. parensis and An. rivulorum-like also detected. All An. gambiae s.l mosquitoes were identified as An. gambiae. High levels of pyrethroid resistance were observed in both An. gambiae and An. funestus mosquitoes. RT-qPCR analyses revealed a significant overexpression of two genes encoding D7 proteins, D7r3 and D7r4, in pyrethroids resistant An. funestus. However, no association was observed between the polymorphism of these genes and their overexpression. In contrast, overall D7 salivary genes were under-expressed in pyrethroid resistant An. gambiae. This study provides preliminary evidences that pyrethroid resistance could influence blood meal intake through over-expression of D7 proteins although future studies will help establishing potential impact on vectorial capacity.

Fitness Costs of the Glutathione S-Transferase Epsilon 2 (L119F-GSTe2) Mediated Metabolic Resistance to Insecticides in the Major African Malaria Vector Anopheles Funestus

Metabolic resistance to insecticides threatens malaria control. However, little is known about its fitness cost in field populations of malaria vectors, thus limiting the design of suitable resistance management strategies. Here, we assessed the association between the glutathione S-transferase GSTe2-mediated metabolic resistance and life-traits of natural populations of Anopheles funestus. A total of 1200 indoor resting blood-fed female An. funestus (F₀) were collected in Mibellon, Cameroon (2016/2017), and allowed to lay eggs individually. Genotyping of F1 mosquitoes for the L119F-GSTE2 mutation revealed that L/L119-homozygote susceptible (SS) mosquitoes significantly laid more eggs than heterozygotes L119F-RS (odds ratio (OR) = 2.06; p < 0.0001) and homozygote resistant 119F/F-RR (OR = 2.93; p < 0.0001). L/L119-SS susceptible mosquitoes also showed the higher ability for oviposition than 119F/F-RR resistant (OR = 2.68; p = 0.0002) indicating a reduced fecundity in resistant mosquitoes. Furthermore, L119F-RS larvae developed faster (nine days) than L119F-RR and L119F-SS (11 days) (X² = 11.052; degree of freedom (df) = 4; p = 0.02) suggesting a heterozygote advantage effect for larval development. Interestingly, L/L119-SS developed faster than 119F/F-RR (OR = 5.3; p < 0.0001) revealing an increased developmental time in resistant mosquitoes. However, genotyping and sequencing revealed that L119F-RR mosquitoes exhibited a higher adult longevity compared to RS (OR > 2.2; p < 0.05) and SS (OR > 2.1; p < 0.05) with an increased frequency of GSTe2-resistant haplotypes in mosquitoes of D30 after adult emergence. Additionally, comparison of the expression of GSTe2 revealed a significantly increased expression from D1-D30 after emergence of adults (Anova test (F) = 8; df= 3; p = 0.008). The negative association between GSTe2 and some life traits of An. funestus could facilitate new resistance management strategies. However, the increased longevity of GSTe2-resistant mosquitoes suggests that an increase in resistance could exacerbate malaria transmission.

Molecular basis of permethrin and DDT resistance in an Anopheles funestus population from Benin

BACKGROUND

Insecticide resistance in Anopheles mosquitoes is threatening the success of malaria control programmes. In order to implement suitable insecticide resistance management strategies, it is necessary to understand the underlying mechanisms involved. To achieve this, the molecular basis of permethrin and DDT resistance in the principal malaria vector, Anopheles funestus from inland Benin (Kpome), was investigated.

RESULTS

Here, using a microarray-based genome-wide transcription and qRT-PCR analysis, we showed that metabolic resistance mechanisms through over-expression of cytochrome P450 and glutathione S-transferase genes (GSTs) are a major contributor to DDT and permethrin resistance in Anopheles funestus from Kpome. The GSTe2 gene was the most upregulated detoxification gene in both DDT- [fold-change (FC: 16.0)] and permethrin-resistant (FC: 18.1) mosquitoes suggesting that upregulation of this gene could contribute to DDT resistance and cross-resistance to permethrin. CYP6P9a and CYP6P9b genes that have been previously associated with pyrethroid resistance were also significantly overexpressed with FC 5.4 and 4.8, respectively, in a permethrin resistant population. Noticeably, the GSTs, GSTd1-5 and GSTd3, were more upregulated in DDT-resistant than in permethrin-resistant Anopheles funestus suggesting these genes are more implicated in DDT resistance. The absence of the L1014F or L1014S kdr mutations in the voltage-gated sodium channel gene coupled with the lack of directional selection at the gene further supported that knockdown resistance plays little role in this resistance.

CONCLUSIONS

The major role played by metabolic resistance to pyrethroids in this An. funestus population in Benin suggests that using novel control tools combining the P450 synergist piperonyl butoxide (PBO), such as PBO-based bednets, could help manage the growing pyrethroid resistance in this malaria vector in Benin.

16S Metagenomic Comparison of Plasmodium falciparum-Infected and Noninfected Anopheles gambiae and Anopheles funestus Microbiota from Senegal

In the context of the pre-elimination of malaria, biological control may provide an alternative or additional tool to current malaria control strategies. During their various stages of development, mosquitoes undergo subsequent changes in their associated microbiota, depending on their environment and nutritional status. Although Anopheles gambiae s.l. and Anopheles funestus are the two major malaria vectors in Senegal, the composition of their microbiota is not yet well known. In this study, we explored the microbiota of mosquitoes naturally infected or not by Plasmodium falciparum (Pf) using the 16S ribosomal RNA gene-based bacterial metagenomic approach. In both vector species, the microbiota was more diverse in Pf-infected samples than in the noninfected ones, although the total number of reads appeared to be higher in noninfected mosquitoes. Overall, the microbiota was different between the two vector species. Noteworthy, the bacterial microbiota was significantly different between Pf-positive and Pf-negative groups whatever the species, but was similar between individuals of the same infection status within a species. Overall, the phylum of Proteobacteria was the most predominant in both species, with bacteria of the genus Burkholderia outweighing the others in noninfected vectors. The presence of some specific bacterial species such as Asaia bogorensis, Enterobacter cloacae, Burkholderia fungorum, and Burkholderia cepacia was also observed in Pf-free samples only. These preliminary observations pave the way for further characterization of the mosquito microbiota to select promising bacterial candidates for potential use in an innovative approach to controlling malaria and overcoming the challenges to achieving a malaria-free world.

Genome-wide transcriptional analyses in Anopheles mosquitoes reveal an unexpected association between salivary gland gene expression and insecticide resistance

BACKGROUND

To combat malaria transmission, the Ugandan government has embarked upon an ambitious programme of indoor residual spraying (IRS) with a carbamate class insecticide, bendiocarb. In preparation for this campaign, we characterized bendiocarb resistance and associated transcriptional variation among Anopheles gambiae s.s. mosquitoes from two sites in Uganda.

RESULTS

Gene expression in two mosquito populations displaying some resistance to bendiocarb (95% and 79% An. gambiae s.l. WHO tube bioassay mortality in Nagongera and Kihihi, respectively) was investigated using whole-genome microarrays. Significant overexpression of several genes encoding salivary gland proteins, including D7r2 and D7r4, was detected in mosquitoes from Nagongera. In Kihihi, D7r4, two detoxification-associated genes (Cyp6m2 and Gstd3) and an epithelial serine protease were among the genes most highly overexpressed in resistant mosquitoes. Following the first round of IRS in Nagongera, bendiocarb-resistant mosquitoes were collected, and real-time quantitative PCR analyses detected significant overexpression of D7r2 and D7r4 in resistant mosquitoes. A single nucleotide polymorphism located in a non-coding transcript downstream of the D7 genes was found at a significantly higher frequency in resistant individuals. In silico modelling of the interaction between D7r4 and bendiocarb demonstrated similarity between the insecticide and serotonin, a known ligand of D7 proteins. A meta-analysis of published microarray studies revealed a recurring association between D7 expression and insecticide resistance across Anopheles species and locations.

CONCLUSIONS

A whole-genome microarray approach identified an association between novel insecticide resistance candidates and bendiocarb resistance in Uganda. In addition, a single nucleotide polymorphism associated with this resistance mechanism was discovered. The use of such impartial screening methods allows for discovery of resistance candidates that have no previously-ascribed function in insecticide binding or detoxification. Characterizing these novel candidates will broaden our understanding of resistance mechanisms and yield new strategies for combatting widespread insecticide resistance among malaria vectors.

Insecticide resistance and target site mutations (G119S ace-1 and L1014F kdr) of Culex pipiens in Morocco

Background

Control of the mosquito vector Culex pipiens with insecticides is the main way to control arboviruses that the species can transmit such as West Nile virus (WNV) and Rift Valley fever virus (RVFV). However, its efficiency has been hampered by the emergence of insecticide resistance. Little is known about the insecticide-resistance status and underlying resistance mechanisms of field-collected populations of Cx. pipiens in Morocco.

Methods

Mosquito adults from Mohammadia city in Morocco were reared from immature stages. The level of their susceptibility to insecticides was assessed using standard WHO bioassay. The two forms of the Cx. pipiens complex and their hybrids were identified by a multiplex PCR. Identified mosquitoes were then tested for the presence of the G119S ace-1 and L1014F kdr mutations using PCR-RFLP and PCR assays, respectively.

Results

WHO bioassays indicated that Cx. pipiens was resistant to all tested insecticides: lambda-cyhalothrin (49% mortality), permethrin (63% mortality), DDT (16% mortality), malation (52% mortality) and bendiocarb (39% mortality). The frequency of the 119S allele was almost identical in the pipiens form and hybrids (0.11 and 0.15, respectively) whereas it remained low in the molestus form (0.03). No significant correlation was observed between the G119S allele and the resistance phenotype to two tested insecticides (malathion and bendiocarb). The frequency of the L1014F allele was identical in the pipiens form and hybrids (0.44) whereas it was low in the molestus form (0.36) but no significant difference was detected (χ² = 1.46, df = 1, P = 0.225). The presence of the L1014F kdr mutation was significantly associated with resistance to three tested insecticides in pipiens form (P = 0.0019, P = 0.0023 and P = 0.023, respectively, to lambda-cyhalothrin, permethrin and DDT) whereas no significant correlation was observed between the L1014F kdr mutation and resistance phenotype in molestus form and hybrids to the three tested insecticides.

Conclusion

These findings showed that wild populations of Cx. pipiens have developed resistance against the main insecticide families with different modes of action: organochlorines (DDT), organophosphates (malathion), carbamates (bendiocarb), pyrethroids (lambda-cyhalothrin, permethrin). Therefore, urgent action should be taken to manage the resistance in this species to maintain the effectiveness of arbovirus control.

Investigating knockdown resistance (kdr) mechanism against pyrethroids/DDT in the malaria vector Anopheles funestus across Africa

Background

Understanding the molecular basis of insecticide resistance is key to improve the surveillance and monitoring of malaria vector populations under control. In the major malaria vector Anopheles funestus, little is currently known about the role of the knockdown resistance (kdr) mechanism. Here, we investigated the presence and contribution of knockdown resistance (kdr) to pyrethroids/DDT resistance observed in Anopheles funestus across Africa.

Results

Pyrosequencing genotyping and sequencing of the voltage gated sodium channel (VGSC) gene did not detect the common L1014F mutation in field collected An. funestus across Africa. Amplification and cloning of the full-length of the sodium channel gene in pyrethroid resistant mosquitoes revealed evidences of alternative splicing events with three transcripts of 2092, 2061 and 2117 amino acids (93% average similarity to An. gambiae). Several amino acid changes were detected close to the domain II of the protein such as L928R, F938 W, I939S, L802S and T1008 M. However, all these mutations are found at low frequency and their role in pyrethroid resistance could not be established. The presence of the exclusive alternative splicing at exon 19 was not associated with resistance phenotype. Analysis of patterns of genetic diversity of the VGSC gene revealed a high polymorphism level of this gene across Africa with no evidence of directional selection suggesting a limited role for knockdown resistance in pyrethroid resistance in An. funestus. Patterns of genetic differentiation correlate with previous observations of the existence of barriers to gene flow Africa-wide with southern population significantly differentiated from other regions.

Conclusion

Despite an apparent limited role of knockdown resistance in An. funestus, it is necessary to continue to monitor the contribution of the mutations detected here as increasing selection from insecticide-based interventions may change the dynamic in field populations as previously observed in other vectors.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0539-x) contains supplementary material, which is available to authorized users.

Genome-Wide Transcription and Functional Analyses Reveal Heterogeneous Molecular Mechanisms Driving Pyrethroids Resistance in the Major Malaria Vector Anopheles funestus Across Africa

Pyrethroid resistance in malaria vector, An. funestus is increasingly reported across Africa, threatening the sustainability of pyrethroid-based control interventions, including long lasting insecticidal nets (LLINs). Managing this problem requires understanding of the molecular basis of the resistance from different regions of the continent, to establish whether it is being driven by a single or independent selective events. Here, using a genome-wide transcription profiling of pyrethroid resistant populations from southern (Malawi), East (Uganda), and West Africa (Benin), we investigated the molecular basis of resistance, revealing strong differences between the different African regions. The duplicated cytochrome P450 genes (CYP6P9a and CYP6P9b) which were highly overexpressed in southern Africa are not the most upregulated in other regions, where other genes are more overexpressed, including GSTe2 in West (Benin) and CYP9K1 in East (Uganda). The lack of directional selection on both CYP6P9a and CYP6P9b in Uganda in contrast to southern Africa further supports the limited role of these genes outside southern Africa. However, other genes such as the P450 CYP9J11 are commonly overexpressed in all countries across Africa. Here, CYP9J11 is functionally characterized and shown to confer resistance to pyrethroids and moderate cross-resistance to carbamates (bendiocarb). The consistent overexpression of GSTe2 in Benin is coupled with a role of allelic variation at this gene as GAL4-UAS transgenic expression in Drosophila flies showed that the resistant 119F allele is highly efficient in conferring both DDT and permethrin resistance than the L119. The heterogeneity in the molecular basis of resistance and cross-resistance to insecticides in An. funestus populations throughout sub-Saharan African should be taken into account in designing resistance management strategies.