Molecular characterization of DDT resistance in Anopheles gambiae from Benin

Physico-chemical characterization of Anopheles gambiae s.l. breeding sites and kdr mutations in urban areas of Cotonou and Natitingou, Benin

BACKGROUND

This study aimed to investigate the relationship between the physicochemical characteristics of An. gambiae s.s. and An. coluzzii breeding sites, the susceptibility profiles to commonly used insecticides in public health, and the underlying insecticide resistance mechanisms.

METHODS

Anopheles breeding sites surveys were conducted in Cotonou and Natitingou in September 2020, January and August 2021. Physicochemical properties and bacterial loads were determined in individual breeding sites. The WHO susceptibility assays were carried out using the female of the emerging adult mosquitoes. Anopheles species were identified through PCR techniques. Kdr L1014F/S, N1575Y and G119S mutations were investigated using TaqMan genotyping assays.

RESULTS

Molecular analysis showed that all mosquitoes analyzed in Cotonou were Anopheles coluzzii, while those of Natitingou were Anopheles gambiae s.s. Fecal coliforms were identified as playing a role in this distribution through their significant influence on the presence of An. coluzzii larvae. WHO susceptibility assay indicated a high level of resistance to deltamethrin in the two cities. The resistance levels to deltamethrin were higher in Cotonou (X2 = 31.689; DF = 1; P < 0.0001). There was a suspected resistance to bendiocarb in Cotonou, whereas the mosquito population in Natitingou was resistant. The kdr L1014F mutation was highly observed in both mosquito populations (frequence: 86-91%), while the Ace-1 mutation was found in a small proportion of mosquitoes. In Cotonou, salinity was the only recorded physicochemical parameter that significantly correlated with the resistance of Anopheles mosquitoes to deltamethrin (P < 0.05). In Natitingou, significant correlations were observed between the allelic frequencies of the kdr L1014F mutation and pH, conductivity, and TDS.

CONCLUSION

These results indicate a high level of pyrethroid resistance in the anopheles populations of both Cotonou and Natitingou. Moreover, this study report the involvement of abiotic factors influencing Anopheles susceptibility profile.

Pyrethroid resistance and gene expression profile of a new resistant An. gambiae colony from Uganda reveals multiple resistance mechanisms and overexpression of Glutathione-S-Transferases linked to survival of PBO-pyrethroid combination

Background

The effectiveness of long-lasting insecticidal nets (LLINs) are being threatened by growing resistance to pyrethroids. To restore their efficacy, a synergist, piperonyl butoxide (PBO) which inhibits cytochrome P450s has been incorporated into pyrethroid treated nets. A trial of PBO-LLINs was conducted in Uganda from 2017 and we attempted to characterize mechanisms of resistance that could impact intervention efficacy.

Methods

We established an Anopheles gambiae s.s colony in 2018 using female mosquitoes collected from Busia district in eastern Uganda. We first assessed the phenotypic resistance profile of this colony using WHO tube and net assays using a deltamethrin dose-response approach. The Busia colony was screened for known resistance markers and RT-qPCR targeting 15 genes previously associated with insecticide resistance was performed.

Results

The Busia colony had very high resistance to deltamethrin, permethrin and DDT. In addition, the colony had moderate resistance to alpha-cypermethrin and lambda-cyhalothrin but were fully susceptible to bendiocarb and fenitrothion. Exposure to PBO in combination with permethrin and deltamethrin resulted in higher mortality rates in both net and tube assays, with a higher mortality observed in net assays than tube assays. The kdr marker, Vgsc-995S was at very high frequency (91.7-98.9%) whilst the metabolic markers Coeae1d and Cyp4j5-L43F were at very low (1.3% - 11.5%) and moderate (39.5% - 44.7%) frequencies respectively. Our analysis showed that gene expression pattern in mosquitoes exposed to deltamethrin, permethrin or DDT only were similar in comparison to the susceptible strain and there was significant overexpression of cytochrome P450s, glutathione-s-transferases (GSTs) and carboxyl esterases (COEs). However, mosquitoes exposed to both PBO and pyrethroid strikingly and significantly only overexpressed closely related GSTs compared to unexposed mosquitoes while major cytochrome P450s were underexpressed.

Conclusions

The high levels of pyrethroid resistance observed in Busia appears associated with a wide range of metabolic gene families.

Evolution of the Ace-1 and Gste2 Mutations and Their Potential Impact on the Use of Carbamate and Organophosphates in IRS for Controlling Anopheles gambiae s.l., the Major Malaria Mosquito in Senegal

Widespread of insecticide resistance amongst the species of the Anopheles gambiae complex continues to threaten vector control in Senegal. In this study, we investigated the presence and evolution of the Ace-1 and Gste2 resistance genes in natural populations of Anopheles gambiae s.l., the main malaria vector in Senegal. Using historical samples collected from ten sentinel health districts, this study focused on three different years (2013, 2017, and 2018) marking the periods of shift between the main public health insecticides families (pyrethroids, carbamates, organophosphates) used in IRS to track back the evolutionary history of the resistance mutations on the Ace-1 and Gste2 loci. The results revealed the presence of four members of the Anopheles gambiae complex, with the predominance of An. arabiensis followed by An. gambiae, An. coluzzii, and An. gambiae–coluzzii hybrids. The Ace-1 mutation was only detected in An. gambiae and An. gambiae–coluzzii hybrids at low frequencies varying between 0.006 and 0.02, while the Gste2 mutation was found in all the species with a frequency ranging between 0.02 and 0.25. The Ace-1 and Gste2 genes were highly diversified with twenty-two and thirty-one different haplotypes, respectively. The neutrality tests on each gene indicated a negative Tajima’s D, suggesting the abundance of rare alleles. The presence and spread of the Ace-1 and Gste2 resistance mutations represent a serious threat to of the effectiveness and the sustainability of IRS-based interventions using carbamates or organophosphates to manage the widespread pyrethroids resistance in Senegal. These data are of the highest importance to support the NMCP for evidence-based vector control interventions selection and targeting.

First report of the L1014F kdr mutation in wild populations of Anopheles arabiensis in Cabo Verde, West Africa

BACKGROUND

Due to the lack of vaccines, malaria control mainly involves the control of anopheline vectors (Anopheles spp.) using chemical insecticides. However, the prolonged and indiscriminate use of these compounds has led to the emergence of resistance in Anopheles populations in Africa. Insecticide resistance surveillance programs are less frequent in Cabo Verde than in other African countries. This study aimed to investigate the circulation of the L1014F and L1014S alleles in natural populations of Anopheles arabiensis collected from two sampling sites in the city of Praia, Cabo Verde.

METHODS

Anopheles larvae were collected from the two sampling sites and reared in the laboratory until the adult stage. Mosquitoes were first morphologically identified by classical taxonomy and then by molecular species identification using molecular markers. All Anopheles arabiensis were subjected to PCR analysis to screen for mutations associated to resistance in the Na_v gene.

RESULTS

A total of 105 mosquitoes, all belonging to the Anopheles gambiae complex, were identified by classical taxonomy as well as by molecular taxonomy. Molecular identification showed that 100% of the An. gambiae senso lato specimens analyzed corresponded to An. arabiensis. Analysis of the Na_v gene revealed the presence of L1014S and L1014F alleles with frequencies of 0.10 and 0.19, respectively.

CONCLUSIONS

Our data demonstrated, for the first time, the presence of the L1014F allele in the An. arabiensis population from Cabo Verde, as well as an increase in the frequency of the kdr L1014S allele reported in a previous study. The results of this study demonstrate the need to establish new approaches in vector control programs in Cabo Verde.

Transcriptome analysis and response of three important detoxifying enzymes to Serratia marcescens Bizio (SM1) in Hyphantria cunea (Drury) (Lepidoptera: Noctuidae)

Hyphantria cunea (Drury) (Lepidoptera: Noctuidae) is a main pest of forest trees. In this study, the effects of Serratia marcescens Bizio (SM1) infection on the transcriptome of H. cunea were studied. The expression of 1068 unigenes in the transcriptome of H. cunea infected by S. marcescens was markedly different from that in the control of H. cunea; 474 genes were upregulated, and 594 genes were downregulated in the former. Among them, 8 cytochrome P450s (CYPs), 5 uridine diphosphate-glycosyltransferases (UGTs) and 3 glutathione S-transferases (GSTs) were significantly differentially expressed. Pathway enrichment analysis indicated that these differentially expressed detoxification enzyme genes were mainly involved in the drug metabolism pathway, glutathione metabolism pathway and ABC transporter pathway. Interestingly, we found that five UGTs were related to oestradiol metabolism in the steroid hormone biosynthesis pathway. Furthermore, the real-time fluorescent quantitative PCR results showed that SM1 could induce the expression of CYPs and UGTs, but inhibit the expression of GSTs. This research will identify the response of important detoxification enzymes to S. marcescens, which will provide a theoretical foundation for the development of new immunosuppressants for H. cunea control. Furthermore, H. cunea was performed transcriptome sequencing to explore the key metabolic pathways, signalling pathways and genes affected by S. marcescens, which will clarify the mechanisms of S. marcescens infection of H. cunea. In addition, this study also explored the relationship between H. cunea and S. marcescens, which will provide a theoretical basis for the biological control of H. cunea by using S. marcescens.

Experimental evolution supports the potential of neonicotinoid-pyrethroid combination for managing insecticide resistance in malaria vectors

The introduction of neonicotinoids for managing insecticide resistance in mosquitoes is of high interest as they interact with a biochemical target not previously used in public health. In this concern, Bayer developed a combination of the neonicotinoid clothianidin and the pyrethroid deltamethrin (brand name Fludora Fusion) as a new vector control tool. Although this combination proved to be efficient against pyrethroid-resistant mosquitoes, its ability to prevent the selection of pyrethroid and neonicotinoid resistance alleles was not investigated. In this context, the objective of this work was to study the dynamics and the molecular mechanisms of resistance of An. gambiae to the separated or combined components of this combination. A field-derived An. gambiae line carrying resistance alleles to multiple insecticides at low frequencies was used as a starting for 33 successive generations of controlled selection. Resistance levels to each insecticide and target site mutation frequencies were monitored throughout the selection process. Cross resistance to other public health insecticides were also investigated. RNA-seq was used to compare gene transcription variations and polymorphisms across all lines. This study confirmed the potential of this insecticide combination to impair the selection of resistance as compared to its two separated components. Deltamethrin selection led to the rapid enrichment of the kdr L1014F target-site mutation. Clothianidin selection led to the over-transcription of multiple cytochrome P450s including some showing high homology with those conferring neonicotinoid resistance in other insects. A strong selection signature associated with clothianidin selection was also observed on a P450 gene cluster previously associated with resistance. Within this cluster, the gene CYP6M1 showed the highest selection signature together with a transcription profile supporting a role in clothianidin resistance. Modelling the impact of point mutations selected by clothianidin on CYP6M1 protein structure showed that selection retained a protein variant with a modified active site potentially enhancing clothianidin metabolism. In the context of the recent deployment of neonicotinoids for mosquito control and their frequent usage in agriculture, the present study highlights the benefit of combining them with other insecticides for preventing the selection of resistance and sustaining vector control activities.

Potentiation Effects of Ficus sycomorus Active Fraction Against Permethrin-Resistant Field-Population of Anopheles coluzzii (Diptera: Culicidae)

Insecticide resistance in mosquitoes is increasing amidst growing cases of global malaria, leading to high fatality in mostly Africa. To overcome the resistance as well as environmental effects of the synthetic insecticides, preliminary insecticidal and botanical potentiating effects of sub-lethal concentration (LC₂₅) Ficus sycomorus active fraction (AFFS) and its synergistic potential with standard insecticide permethrin were evaluated against malarial vector Anopheles coluzzii (Coetzee & Wilkerson) populations. The glutathione-S-transferase (GST) inhibitory activity of the AFFS was also investigated compared to standard GST inhibitor, diethyl meleate (DEM). The WHO standard protocol for adult bioassay was used to expose the adult mosquitoes with sub-lethal concentration (LD₂₅=0.49 mg/ml) of the plants' active fraction and permethrin (0.75%). The permethrin susceptibility screening result showed high level of resistance to permethrin in the field populations of A. coluzzii from Kano with 50.29 ± 2.14% average mortality after exposure to WHO diagnostic dose 0.75% permethrin. Post hoc Fisher's exact test showed that combination of sub-lethal concentration of AFFS with permethrin (mortality=73.02±12.10%; p=0.00352; RR=0.6923 and 95% CI = 0.5358-0.8946) was statistically significant, while the combination of sub-lethal concentration of AFFS with DEM showed no statistical difference (mortality=63.22±5.03; p=1; RR=0.6667 and 95% CI=0.4470-0.8438). This potentiation effect was signified to be additive effects with co-toxicity factor (CTF) of - 12.66. There was significant reduction of GST activities in the AFFS- and permethrin -exposed groups compared to unexposed populations of A. coluzzii (p < 0.05). The AFFS additively potentiate the permethrin activities by inhibiting GSTs, bio-transformational enzymes implicated in pyrethroids resistance. This study finding generally signifies the potential for bio-rational insecticide approach for malarial vector control.

Genetic variation at the Cyp6m2 putative insecticide resistance locus in Anopheles gambiae and Anopheles coluzzii

Background

The emergence of insecticide resistance is a major threat to malaria control programmes in Africa, with many different factors contributing to insecticide resistance in its vectors, Anopheles mosquitoes. CYP6M2 has previously been recognized as an important candidate in cytochrome P450-mediated detoxification in Anopheles. As it has been implicated in resistance against pyrethroids, organochlorines and carbamates, its broad metabolic activity makes it a potential agent in insecticide cross-resistance. Currently, allelic variation within the Cyp6m2 gene remains unknown.

Methods

Here, Illumina whole-genome sequence data from Phase 2 of the Anopheles gambiae 1000 Genomes Project (Ag1000G) was used to examine genetic variation in the Cyp6m2 gene across 16 populations in 13 countries comprising Anopheles gambiae and Anopheles coluzzii mosquitoes. To identify whether these alleles show evidence of selection either through potentially modified enzymatic function or by being linked to variants that change the transcriptional profile of the gene, hierarchical clustering of haplotypes, linkage disequilibrium, median joining networks and extended haplotype homozygosity analyses were performed.

Results

Fifteen missense biallelic substitutions at high frequency (defined as > 5% frequency in one or more populations) are found, which fall into five distinct haplotype groups that carry the main high frequency variants: A13T, D65A, E328Q, Y347F, I359V and A468S. Despite consistent reports of Cyp6m2 upregulation and metabolic activity in insecticide resistant Anophelines, no evidence of directional selection is found occurring on these variants or on the haplotype clusters in which they are found.

Conclusion

These results imply that emerging resistance associated with Cyp6m2 is potentially driven by distant regulatory loci such as transcriptional factors rather than by its missense variants, or that other genes are playing a more significant role in conferring metabolic resistance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-021-03757-4.

First report of N1575Y mutation in Anopheles gambiae in Sierra Leone

The resistance of mosquito vectors to insecticides is one of the biological obstacles in the fight against malaria. Understanding of the status and mechanisms underlying the insecticide resistance in Anopheles gambiae species is necessary for success of vector control efforts. The study aimed to determine the molecular forms of An. gambiae from four districts in Sierra Leone during May and June 2018, and the level of N1575Y mutation. The molecular form identification of adult female An. gambiae mosquitoes reared from larvae were carried out using polymerase chain reaction and sequencing. And the N1575Y mutations were detected using SNaPshot and sequencing. As a result, significant differences were found in the distribution of An. gambiae molecular forms among regions (P < 0.001). And a total of 638 An. gambiae sensu stricto, 106 An. coluzzi, and 4 hybrid individuals were identified. Moreover, the overall N1575Y mutation frequency was 10.2% with no statistical difference among regions (χ² = 3.009, P = 0.390). In addition, no significant differences in N1575Y mutation frequency were found among different An. gambiae molecular forms (P = 0.383). In conclusion, the N1575Y mutation in An. gambiae populations in Sierra Leone was reported for the first time in the present study. It provides key evidence for the necessity of monitoring vector susceptibility levels to insecticides used in this country.

Anopheles metabolic proteins in malaria transmission, prevention and control: a review

The increasing resistance to currently available insecticides in the malaria vector, Anopheles mosquitoes, hampers their use as an effective vector control strategy for the prevention of malaria transmission. Therefore, there is need for new insecticides and/or alternative vector control strategies, the development of which relies on the identification of possible targets in Anopheles. Some known and promising targets for the prevention or control of malaria transmission exist among Anopheles metabolic proteins. This review aims to elucidate the current and potential contribution of Anopheles metabolic proteins to malaria transmission and control. Highlighted are the roles of metabolic proteins as insecticide targets, in blood digestion and immune response as well as their contribution to insecticide resistance and Plasmodium parasite development. Furthermore, strategies by which these metabolic proteins can be utilized for vector control are described. Inhibitors of Anopheles metabolic proteins that are designed based on target specificity can yield insecticides with no significant toxicity to non-target species. These metabolic modulators combined with each other or with synergists, sterilants, and transmission-blocking agents in a single product, can yield potent malaria intervention strategies. These combinations can provide multiple means of controlling the vector. Also, they can help to slow down the development of insecticide resistance. Moreover, some metabolic proteins can be modulated for mosquito population replacement or suppression strategies, which will significantly help to curb malaria transmission.

Insecticide resistance in indoor and outdoor-resting Anopheles gambiae in Northern Ghana

BACKGROUND

Selection pressure from continued exposure to insecticides drives development of insecticide resistance and changes in resting behaviour of malaria vectors. There is need to understand how resistance drives changes in resting behaviour within vector species. The association between insecticide resistance and resting behaviour of Anopheles gambiae sensu lato (s.l.) in Northern Ghana was examined.

METHODS

F1 progenies from adult mosquitoes collected indoors and outdoors were exposed to DDT, deltamethrin, malathion and bendiocarb using WHO insecticide susceptibility tests. Insecticide resistance markers including voltage-gated sodium channel (Vgsc)-1014F, Vgsc-1014S, Vgsc-1575Y, glutathione-S-transferase epsilon 2 (GSTe2)-114T and acetylcholinesterase (Ace1)-119S, as well as blood meal sources were investigated using PCR methods. Activities of metabolic enzymes, acetylcholine esterase (AChE), non-specific β-esterases, glutathione-S-transferase (GST) and monooxygenases were measured from unexposed F1 progenies using microplate assays.

RESULTS

Susceptibility of Anopheles coluzzii to deltamethrin 24 h post-exposure was significantly higher in indoor (mortality = 5%) than outdoor (mortality = 2.5%) populations (P = 0.02). Mosquitoes were fully susceptible to malathion (mortality: indoor = 98%, outdoor = 100%). Susceptibility to DDT was significantly higher in outdoor (mortality = 9%) than indoor (mortality = 0%) mosquitoes (P = 0.006). Mosquitoes were also found with suspected resistance to bendiocarb but mortality was not statistically different (mortality: indoor = 90%, outdoor = 95%. P = 0.30). Frequencies of all resistance alleles were higher in F1 outdoor (0.11-0.85) than indoor (0.04-0.65) mosquito populations, while Vgsc-1014F in F0 An. gambiae sensu stricto (s.s) was significantly associated with outdoor-resting behaviour (P = 0.01). Activities of non-specific β-esterase enzymes were significantly higher in outdoor than indoor mosquitoes (Mean enzyme activity: Outdoor = : 1.70/mg protein; Indoor = 1.35/mg protein. P < 0.0001). AChE activity was also more elevated in outdoor (0.62/mg protein) than indoor (0.57/mg protein) mosquitoes but this was not significant (P = 0.08). Human blood index (HBI) was predominantly detected in indoor (18%) than outdoor mosquito populations (3%).

CONCLUSIONS

The overall results did not establish that there was a significant preference of resistant malaria vectors to solely rest indoors or outdoors, but varied depending on the resistant alleles present. Phenotypic resistance was higher in indoor than outdoor-resting mosquitoes, but genotypic and metabolic resistance levels were higher in outdoor than the indoor populations. Continued monitoring of changes in resting behaviour within An. gambiae s.l. populations is recommended.

Status of insecticide resistance and its biochemical and molecular mechanisms in Anopheles stephensi (Diptera: Culicidae) from Afghanistan

Background

Insecticide resistance of Anopheles stephensi, the main malaria vector in eastern Afghanistan, has been reported previously. This study describes the biochemical and molecular mechanisms of resistance to facilitate effective vector control and insecticide resistance management.

Methods

Mosquito larvae were collected from the provinces of Kunar, Laghman and Nangarhar from 2014 to 2017. The susceptibility of the reared 3–4 days old adults was tested with deltamethrin 0.05%, bendiocarb 0.1%, malathion 5%, permethrin 0.75% and DDT 4%. Cytochrome P450 content and general esterase, glutathione S-transferase (GST) and acetylcholinesterase (AChE) activities were measured in the three field populations and the results were compared with those of the laboratory susceptible An. stephensi Beech strain. Two separate allele-specific PCR assays were used to identify L1014, L1014F and L1014S mutations in the voltage gated sodium channel gene of An. stephensi. Probit analysis, ANOVA and Hardy–Weinberg equilibrium were used to analyse bioassay, biochemical assay and gene frequency data respectively.

Results

The population of An. stephensi from Kunar was susceptible to bendiocarb, apart from this, all populations were resistant to all the other insecticides tested. The differences between all values for cytochrome P450s, general esterases, GSTs and AChE inhibition rates in the Kunar, Laghman and Nangarhar populations were statistically significant when compared to the Beech strain, excluding GST activities between Kunar and Beech due to the high standard deviation in Kunar. The three different sodium channel alleles [L1014 (wild type), L1014F (kdr west) and L1014S (kdr east)] were all segregated in the Afghan populations. The frequencies of kdr east mutation were 22.9%, 32.7% and 35% in Kunar, Laghman and Nangarhar populations respectively. Kdr west was at the lowest frequency of 4.44%.

Conclusions

Resistance to different groups of insecticides in the field populations of An. stephensi from Kunar, Laghman and Nangarhar Provinces of Afghanistan is caused by a range of metabolic and site insensitivity mechanisms, including esterases, cytochrome P450s and GSTs combined with AChE and sodium channel target site insensitivity. The intensity and frequency of these mechanisms are increasing in these populations, calling for urgent reorientation of vector control programmes and implementation of insecticide resistance management strategies.

Genome-wide gene expression profiling reveals that cuticle alterations and P450 detoxification are associated with deltamethrin and DDT resistance in Anopheles arabiensis populations from Ethiopia

BACKGROUND

Vector control is the main intervention in malaria control and elimination strategies. However, the development of insecticide resistance is one of the major challenges for controlling malaria vectors. Anopheles arabiensis populations in Ethiopia showed resistance against both DDT and the pyrethroid deltamethrin. Although an L1014F target-site resistance mutation was present in the voltage gated sodium channel of investigated populations, the levels of resistance indicated the presence of additional resistance mechanisms. In this study, we used genome-wide transcriptome profiling by RNAseq to assess differentially expressed genes between three deltamethrin and DDT resistant An. arabiensis field populations - Asendabo, Chewaka and Tolay - and two susceptible strains - Sekoru and Mozambique.

RESULTS

Both RNAseq analysis and RT-qPCR showed that a glutathione-S-transferase, gstd3, and a cytochrome P450 monooxygenase, cyp6p4, were significantly overexpressed in the group of resistant populations compared to the susceptible strains, suggesting that the enzymes they encode play a key role in metabolic resistance against deltamethrin or DDT. Furthermore, a gene ontology enrichment analysis showed that expression changes of cuticle related genes were strongly associated with insecticide resistance. Although this did not translate in increased thickness of the procuticle, a higher cuticular hydrocarbon content was observed in a resistant population.

CONCLUSION

Our transcriptome sequencing of deltamethrin and DDT resistant An. arabiensis populations from Ethiopia suggests non-target site resistance mechanisms and paves the way for further investigation of the role of cuticle composition in insecticide resistance of malaria vectors. © 2019 Society of Chemical Industry.

Detection of alleles associated with resistance to chemical insecticide in the malaria vector Anopheles arabiensis in Santiago, Cabo Verde

BACKGROUND

Mosquitoes of the Anopheles gambiae complex are the main malaria vectors worldwide. Due to the lack of a vaccine to prevent malaria, the principal way to reduce the impact of this disease relies on the use of chemical insecticides to control its vectors. However, the intensive use of such compounds has led to the emergence of insecticide resistance in several Anopheles populations in Africa. This study aimed to investigate the presence of resistance alleles in an Anopheles arabiensis population from the City of Praia, capital of the Archipelago Cabo Verde, one of the countries on the World Health Organization list of countries that are on a path to eliminate local transmission of malaria.

METHODS

Larvae from the Anopheles genus were collected using a one-pint dipper in three areas of City of Praia. Larvae were fed and maintained until the emergence of adult mosquitoes, and these were morphologically identified. In addition, molecular identification was performed using IGS markers and all An. arabiensis samples were subjected to PCR to screen for mutations associated to resistance in the Ace-1, Na_v and GSTE2 genes.

RESULTS

From a total of 440 mosquitoes collected, 52.3% were morphologically identified as An. gambiae sensu lato (s.l.) and 46.7% as Anopheles pretoriensis. The molecular identification showed that 100% of the An. gambiae s.l. were An. arabiensis. The mutations G119S in the Ace-1 gene and L119F in the GSTE2 gene were screened but not found in any sample. However, sequencing analysis for GSTE2 revealed the presence of 37 haplotypes, 16 polymorphic sites and a high genetic diversity (π = 2.67). The L1014S mutation in the Na_v (voltage-gated sodium channel gene) was detected at a frequency of 7.3%.

CONCLUSION

This is the first study to investigate the circulation of insecticide resistance alleles in An. arabiensis from Cabo Verde. The circulation of the L1014S allele in the population of An. arabiensis in the city of Praia suggests that pyrethroid resistance may arise, be quickly selected, and may affect the process of malaria elimination in Cabo Verde. Molecular monitoring of resistance should continue in order to guide the development of strategies to be used in vector control in the study region.

Insecticide resistance in Anopheles gambiae from the northern Democratic Republic of Congo, with extreme knockdown resistance (kdr) mutation frequencies revealed by a new diagnostic assay

BACKGROUND

Mutations in the voltage-gated sodium channel at codon 1014 confer knock-down resistance (kdr) to pyrethroids in a wide range of insects. Anopheles gambiae exhibits two mutant alleles at codon 1014, serine and phenylalanine; and both are now widespread across Africa. Existing screening methods only allow for one resistant allele to be detected per assay. A new locked nucleic acid (LNA) qPCR assay was developed for the simultaneous detection of both mutant alleles and the wild type allele in a single assay. This tri-allelic detection assay was assessed as part of a study of the insecticide resistance in An. gambiae sensu stricto (s.s.) in the previously un-sampled area of Nord Ubangi, Democratic Republic of the Congo.

METHODS

Samples from three sites were tested for insecticide susceptibility using WHO bioassays, with and without the synergist PBO preceding pyrethroid exposures, and were subsequently analysed for frequency and resistance-association of the Vgsc-1014 and Vgsc-N1575Y mutations. Results from the LNA-kdr 1014 assay were compared to results from standard TaqMan-kdr assays.

RESULTS

Anopheles gambiae sensu lato (s.l.) was by far the predominant vector captured (84%), with only low frequencies of Anopheles funestus s.l. (9%) detected in Nord Ubangi. Molecular identification found An. gambiae s.s. to be the principal vector (99%) although Anopheles coluzzii was detected at very low frequency. Anopheles gambiae were susceptible to the carbamate insecticide bendiocarb, but resistant to DDT and to the pyrethroids permethrin and deltamethrin. Susceptibility to both pyrethroids was partially restored with prior exposure to PBO suggesting likely involvement of metabolic resistance. Anopheles gambiae s.s. was homozygous for kdr resistant alleles with both the L1014F and L1014S mutations present, and the N1575Y polymorphism was present at low frequency. The LNA-kdr assay simultaneously detected both resistant alleles and gave results entirely consistent with those from the two TaqMan-kdr assays.

CONCLUSION

This study provides rare data on insecticide resistance and mechanisms in Anopheles from the centre of Africa, with the first detection of N1575Y. Nord Ubangi populations of An. gambiae s.s. show insecticide resistance mediated by both metabolic mechanisms and Vgsc mutations. The LNA-kdr assay is particularly suitable for use in populations in which both 1014S and 1014F kdr alleles co-occur and provides robust results, with higher throughput and at a quarter of the cost of TaqMan assays.

First report of the presence of L1014S Knockdown-resistance mutation in Anopheles gambiae s.s and Anopheles coluzzii from Togo, West Africa

Background: To optimize the success of insecticide-based malaria control intervention, knowledge of the distribution of Anopheles gambiae species and insecticide resistance mechanisms is necessary. This paper reported an updated data on pyrethroids/DDT resistance in the An. gambiae s.l population from Togo.  Methods: From December 2013 to April 2015, females of indoor-resting An. gambiae s.l were captured in three locations belonging to three different ecological zones. Resistance to DDT, permethrin and deltamethrin was screened in F1 progeny of collected mosquitoes using WHO susceptibility tests. The identification of species of An. gambiae complex and the detection of kdr and ace.1 R allele were carried out using DNA-based molecular techniques. Results:An. gambiae from Kovié and Nangbéto were highly resistant to DDT and permethrin with mortalities rate ranging from 0.83% to 1.58% for DDT and zero to 8.54% for permethrin. Mosquitoes collected in Nangbéto displayed 81.53% mortality with deltamethrin. An. coluzzii and An. gambiae s.s were found in sympatry in Nangbéto and Mango . The allelic frequency of L1014F was high, ranging from 66 to 100% in both An. coluzzii and An. gambiae s.s. For the first time we detected the L1014S allele in both An. coluzzii and An. gambiaes.s. from Togo at the frequency ranging from 5% to 13% in all the sites. The kdr N1575Y was present at various frequencies in both species ranging from 10% to 45%. Both An. gambiae s.s. and An. coluzzii shared the ace1 R mutation in all investigated sites with allelic frequency ranging from 4% to 16%. Conclusion: These results showed that multiple mutations are involved in insecticides resistance in An. gambiae populations from Togo including the kdr L1014F, L1014S, and N1575Y and ace.1 R G119S mutations.

Contributions of cuticle permeability and enzyme detoxification to pyrethroid resistance in the major malaria vector Anopheles gambiae

To tackle the problem of insecticide resistance, all resistance mechanisms need to be studied. This study investigated the involvement of the cuticle in pyrethroid resistance in a strain of Anopheles gambiae, MRS, free of kdr mutations. Bioassays revealed MRS to be resistant to pyrethroids and DDT, indicated by increasing knockdown times and resistance ratios. Moreover, biochemical analysis indicated that metabolic resistance based on enhanced CYP450 activity may also play a role. Insecticide penetration assays showed that there were significantly lower amounts of insecticide in the MRS strain than in the susceptible control. Analysis of the levels of the selected transcripts by qPCR showed that CYP6M2, a major pyrethroid metaboliser, CYP4G16, a gene implicated in resistance via its contribution to the biosynthesis of elevated epicuticular hydrocarbons that delay insecticide uptake, and the cuticle genes CPAP3-E and CPLCX1 were upregulated after insecticide exposure. Other metabolic (CYP6P3, GSTe2) and cuticle (CPLCG3, CPRs) genes were also constitutively upregulated. Microscopic analysis showed that the cuticle layers of the MRS strain were significantly thicker than those of the susceptible strain. This study allowed us to assess the contribution made by the cuticle and metabolic mechanisms to pyrethroid resistance in Anopheles gambiae without target-site mutations.

Chromosomal inversions and ecotypic differentiation in Anopheles gambiae: the perspective from whole-genome sequencing

The molecular mechanisms and genetic architecture that facilitate adaptive radiation of lineages remain elusive. Polymorphic chromosomal inversions, due to their recombination-reducing effect, are proposed instruments of ecotypic differentiation. Here, we study an ecologically diversifying lineage of Anopheles gambiae, known as the Bamako chromosomal form based on its unique complement of three chromosomal inversions, to explore the impact of these inversions on ecotypic differentiation. We used pooled and individual genome sequencing of Bamako, typical (non-Bamako) An. gambiae and the sister species Anopheles coluzzii to investigate evolutionary relationships and genomewide patterns of nucleotide diversity and differentiation among lineages. Despite extensive shared polymorphism and limited differentiation from the other taxa, Bamako clusters apart from the other taxa, and forms a maximally supported clade in neighbour-joining trees based on whole-genome data (including inversions) or solely on collinear regions. Nevertheless, F_ST outlier analysis reveals that the majority of differentiated regions between Bamako and typical An. gambiae are located inside chromosomal inversions, consistent with their role in the ecological isolation of Bamako. Exceptionally differentiated genomic regions were enriched for genes implicated in nervous system development and signalling. Candidate genes associated with a selective sweep unique to Bamako contain substitutions not observed in sympatric samples of the other taxa, and several insecticide resistance gene alleles shared between Bamako and other taxa segregate at sharply different frequencies in these samples. Bamako represents a useful window into the initial stages of ecological and genomic differentiation from sympatric populations in this important group of malaria vectors.

Multiple insecticide resistance in an infected population of the malaria vector Anopheles funestus in Benin

Background

Knowledge on the spread and distribution of insecticide resistance in major malaria vectors such as Anopheles funestus is key to implement successful resistance management strategies across Africa. Here, by assessing the susceptibility status of an inland population of An. funestus Giles (Kpome) and investigating molecular basis of resistance, we show that multiple resistance and consistent plasmodium infection rate are present in Anopheles funestus populations from Kpome.

Methods

The insecticide susceptibility level of collected Anopheles funestus was assessed. Synergist (PBO) was used to screen resistance mechanisms. The TaqMan technique was used for genotyping of insecticide resistant alleles and detecting plasmodium infection levels. The nested PCR was used to further assess the plasmodium infection rate.

Results

The TaqMan analysis of plasmodial infections revealed an infection rate (18.2 %) of An. funestus in this locality. The WHO bioassays revealed a multiple phenotypic resistance profile for An. funestus in Kpome. This population is highly resistant to pyrethroids (permethrin and deltamethrin), organochlorines (DDT), and carbamates (bendiocarb). A reduced susceptibility was observed with dieldrin. Mortalities did not vary after pre-exposure to PBO for DDT indicating that cytochrome P450s play little role in DDT resistance in Kpome. In contrast, we noticed, a significant increase in mortalities when PBO was combined to permethrin suggesting the direct involvement of P450s in pyrethroid resistance. A high frequency of the L119F-GSTe2 DDT resistance marker was observed in the wild DDT resistant population (9 %RS and 91 %RR) whereas the A296S mutation was detected at a low frequency (1 %RS and 99 %SS).

Conclusion

The presence of multiple resistance in An. funestus populations in the inland locality of Kpome is established in this study as recently documented in the costal locality of Pahou. Data from both localities suggest that resistance could be widespread in Benin and this highlights the need for further studies to assess the geographical distribution of insecticide resistance across Benin and neighboring countries as well as a more comprehensive analysis of the resistance mechanisms involved.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1723-y) contains supplementary material, which is available to authorized users.

Dynamics of pyrethroid resistance in malaria vectors in southern Benin following a large scale implementation of vector control interventions

Background

Large-scale implementation of Indoor Residual Spraying and Insecticide Treated Nets has been implemented in Plateau Department, Benin between 2011 and 2014. The purpose of this study was to monitor the frequency and mechanisms of pyrethroid resistance in malaria vectors following the implementation of vector control tools for malaria prevention.

Methods

Anopheles larvae were collected in 13 villages twice a year from 2012 to 2014. WHO tube tests were used to assess the phenotypic resistance of each population to 0.05 % deltamethrin. Sibling species within Anopheles gambiae complex were identified by PCR techniques. Taqman and biochemical assays were performed to identify the presence of kdr mutations in individual mosquitoes and to detect any increase in the activity of enzymes putatively involved in insecticide metabolism (oxidases, esterase and glutathione-S-transferases). Quantitative real time PCR was used to measure the expression of three metabolic genes involved in pyrethroid resistance (CYP6P3, CYP6M2 and GSTD3).

Results

Anopheles populations showed < 90 % mortality to deltamethrin in all villages and at all time points. The 1014 F kdr allele frequency was close to fixation (> 0.9) over the sampling periods in both An. gambiae and An. coluzzii. Biochemical assays showed higher activities of alpha esterase and GST in field malaria vector populations compared to susceptible mosquitoes. qPCR assays showed a significant increase of CYP6P3, CYP6M2 GSTD3 expression in An. gambiae after a three-year implementation of LLINs.

Conclusion

The study confirmed that deltamethrin resistance is widespread in malaria vectors in Southern Benin. We suspect that the increase in deltamethrin resistance between 2012 and 2014 resulted from an increased expression of metabolic detoxification genes (CYP6M2 and CYP6P3) rather than from kdr mutations. It is urgent to evaluate further the impact of metabolic resistance on the efficacy of vector control interventions using pyrethroid insecticides.

Independent origins of loss-of-function mutations conferring oxamniquine resistance in a Brazilian schistosome population

Molecular surveillance provides a powerful approach to monitoring the resistance status of parasite populations in the field and for understanding resistance evolution. Oxamniquine was used to treat Brazilian schistosomiasis patients (mid-1970s to mid-2000s) and several cases of parasite infections resistant to treatment were recorded. The gene underlying resistance (SmSULT-OR) encodes a sulfotransferase required for intracellular drug activation. Resistance has a recessive basis and occurs when both SmSULT-OR alleles encode for defective proteins. Here we examine SmSULT-OR sequence variation in a natural schistosome population in Brazil ∼40years after the first use of this drug. We sequenced SmSULT-OR from 189 individual miracidia (1-11 per patient) recovered from 49 patients, and tested proteins expressed from putative resistance alleles for their ability to activate oxamniquine. We found nine mutations (four non-synonymous single nucleotide polymorphisms, three non-coding single nucleotide polymorphisms and two indels). Both mutations (p.E142del and p.C35R) identified previously were recovered in this field population. We also found two additional mutations (a splice site variant and 1bp coding insertion) predicted to encode non-functional truncated proteins. Two additional substitutions (p.G206V, p.N215Y) tested had no impact on oxamniquine activation. Three results are of particular interest: (i) we recovered the p.E142del mutation from the field: this same deletion is responsible for resistance in an oxamniquine selected laboratory parasite population; (ii) frequencies of resistance alleles are extremely low (0.27-0.8%), perhaps due to fitness costs associated with carriage of these alleles; (iii) that four independent resistant alleles were found is consistent with the idea that multiple mutations can generate loss-of-function alleles.

Insecticide Resistance and Management Strategies in Urban Ecosystems

The increased urbanization of a growing global population makes imperative the development of sustainable integrated pest management (IPM) strategies for urban pest control. This emphasizes pests that are closely associated with the health and wellbeing of humans and domesticated animals. Concurrently there are regulatory requirements enforced to minimize inadvertent exposures to insecticides in the urban environment. Development of insecticide resistance management (IRM) strategies in urban ecosystems involves understanding the status and mechanisms of insecticide resistance and reducing insecticide selection pressure by combining multiple chemical and non-chemical approaches. In this review, we will focus on the commonly used insecticides and molecular and physiological mechanisms underlying insecticide resistance in six major urban insect pests: house fly, German cockroach, mosquitoes, red flour beetle, bed bugs and head louse. We will also discuss several strategies that may prove promising for future urban IPM programs.

Impact of Insecticide Resistance on the Effectiveness of Pyrethroid-Based Malaria Vectors Control Tools in Benin: Decreased Toxicity and Repellent Effect

Since the first evidence of pyrethroids resistance in 1999 in Benin, mutations have rapidly increased in mosquitoes and it is now difficult to design a study including a control area where malaria vectors are fully susceptible. Few studies have assessed the after effect of resistance on the success of pyrethroid based prevention methods in mosquito populations. We therefore assessed the impact of resistance on the effectiveness of pyrethroids based indoor residual spraying (IRS) in semi-field conditions and long lasting insecticidal nets (LLINs) in laboratory conditions. The results observed showed low repulsion and low toxicity of pyrethroids compounds in the test populations. The toxicity of pyrethroids used in IRS was significantly low with An. gambiae s.l (< 46%) but high for other predominant species such as Mansonia africana (93% to 97%). There were significant differences in terms of the repellent effect expressed as exophily and deterrence compared to the untreated huts (P<0.001). Furthermore, mortality was 23.71% for OlyseNet® and 39.06% for PermaNet®. However, with laboratory susceptible "Kisumu", mortality was 100% for both nets suggesting a resistance within the wild mosquito populations. Thus treatment with pyrethroids at World Health Organization recommended dose will not be effective at reducing malaria in the coming years. Therefore it is necessary to study how insecticide resistance decreases the efficacy of particular pyrethroids used in pyrethroid-based vector control so that a targeted approach can be adopted.

Genomic Analysis of Detoxification Supergene Families in the Mosquito Anopheles sinensis

Anopheles sinensis is an important malaria vector in China and other Southeast Asian countries, and the emergence of insecticide resistance in this mosquito poses a serious threat to the efficacy of malaria control programs. The recently published An. sinensis genome and transcriptome provide an opportunity to understand the molecular mechanisms of insecticide resistance. Analysis of the An. sinensis genome revealed 174 detoxification genes, including 93 cytochrome P450s (P450s), 31 glutathione-S-transferases (GSTs), and 50 choline/carboxylesterases (CCEs). The gene number was similar to that in An. gambiae, but represented a decrease of 29% and 42% compared with Aedes aegypti and Culex quinquefasciatus, respectively. The considerable contraction in gene number in Anopheles mosquitoes mainly occurred in two detoxification supergene families, P450s and CCEs. The available An. sinensis transcriptome was also re-analyzed to further identify key resistance-associated detoxification genes. Among 174 detoxification genes, 124 (71%) were detected. Several candidate genes overexpressed in a deltamethrin-resistant strain (DR-strain) were identified as belonging to the CYP4 or CYP6 family of P450s and the Delta GST class. These generated data provide a basis for identifying the resistance-associated genes of An. sinensis at the molecular level.

The kdr-bearing haplotype and susceptibility to Plasmodium falciparum in Anopheles gambiae: genetic correlation and functional testing

Background

Members of the Anophelesgambiae species complex are primary vectors of human malaria in Africa. It is known that a large haplotype shared between An. gambiae and Anophelescoluzzii by introgression carries point mutations of the voltage-gated sodium channel gene para, including the L1014F kdr mutation associated with insensitivity to pyrethroid insecticides. Carriage of L1014F kdr is also correlated with higher susceptibility to infection with Plasmodium falciparum. However, the genetic mechanism and causative gene(s) underlying the parasite susceptibility phenotype are not known.

Methods

Mosquitoes from the wild Burkina Faso population were challenged by feeding on natural P. falciparum gametocytes. Oocyst infection phenotypes were determined and were tested for association with SNP genotypes. Candidate genes in the detected locus were prioritized and RNAi-mediated gene silencing was used to functionally test for gene effects on P. falciparum susceptibility.

Results

A genetic locus, Pfin6, was identified that influences infection levels of P. falciparum in mosquitoes. The locus segregates as a ~3 Mb haplotype carrying 65 predicted genes including the para gene. The haplotype carrying the kdr allele of para is linked to increased parasite infection prevalence, but many single nucleotide polymorphisms on the haplotype are also equally linked to the infection phenotype. Candidate genes in the haplotype were prioritized and functionally tested. Silencing of para did not influence P. falciparum infection, while silencing of a predicted immune gene, serine protease ClipC9, allowed development of significantly increased parasite numbers.

Conclusions

Genetic variation influencing Plasmodium infection in wild Anopheles is linked to a natural ~3 megabase haplotype on chromosome 2L that carries the kdr allele of the para gene. Evidence suggests that para gene function does not directly influence parasite susceptibility, and the association of kdr with infection may be due to tight linkage of kdr with other gene(s) on the haplotype. Further work will be required to determine if ClipC9 influences the outcome of P. falciparum infection in nature, as well as to confirm the absence of a direct influence by para.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-015-0924-8) contains supplementary material, which is available to authorized users.