Mapping a quantitative trait locus (QTL) conferring pyrethroid resistance in the African malaria vector Anopheles funestus

Overexpression and nonsynonymous mutations of UDP-glycosyltransferases are potentially associated with pyrethroid resistance in Anopheles funestus

UDP-glycosyltransferases (UGTs) enzymes are pivotal in insecticide resistance by transforming hydrophobic substrates into more hydrophilic forms for efficient cell elimination. This study provides the first comprehensive investigation of Anopheles funestus UGT genes, their evolution, and their association with pyrethroid resistance. We employed a genome-wide association study using pooled sequencing (GWAS-PoolSeq) and transcriptomics on pyrethroid-resistant An. funestus, along with deep-targeted sequencing of UGTs in 80 mosquitoes Africa-wide. UGT310B2 was consistently overexpressed Africa-wide and significant gene-wise Fst differentiation was observed between resistant and susceptible populations: UGT301C2 and UGT302A3 in Malawi, and UGT306C2 in Uganda. Additionally, nonsynonymous mutations in UGT genes were identified. Gene-wise Tajima's D density curves provide insights into population structures within populations across these countries, supporting previous observations. These findings have important implications for current An. funestus control strategies facilitating the prediction of cross-resistance to other UGT-metabolised polar insecticides, thereby guiding more effective and targeted insecticide resistance management efforts.

Non-Coding RNAs Potentially Involved in Pyrethroid Resistance of Anopheles funestus Population in Western Kenya

Backgrounds The resurgence of Anopheles funestus , a dominant vector of human malaria in western Kenya was partly attributed to insecticide resistance. However, evidence on the molecular basis of pyrethroid resistance in western Kenya is limited. Noncoding RNAs (ncRNAs) form a vast class of RNAs that do not code for proteins and are ubiquitous in the insect genome. Here, we demonstrated that multiple ncRNAs could play a potential role in An. funestus resistance to pyrethroid in western Kenya. Materials and Methods Anopheles funestus mosquitoes were sampled by aspiration methods in Bungoma, Teso, Siaya, Port Victoria and Kombewa in western Kenya. The F1 progenies were exposed to deltamethrin (0.05%), permethrin (0.75%), DDT (4%) and pirimiphos-methyl (0.25%) following WHO test guidelines. A synergist assay using piperonyl butoxide (PBO) (4%) was conducted to determine cytochrome P450s' role in pyrethroid resistance. RNA-seq was conducted on a combined pool of specimens that were resistant and unexposed, and the results were compared with those of the FANG susceptible strain. This approach aimed to uncover the molecular mechanisms underlying pyrethroid resistance. Results Pyrethroid resistance was observed in all the sites with an average mortality rate of 57.6%. Port Victoria had the highest level of resistance to permethrin (MR=53%) and deltamethrin (MR=11%) pyrethroids. Teso had the lowest level of resistance to permethrin (MR=70%) and deltamethrin (MR=87%). Resistance to DDT was observed only in Kombewa (MR=89%) and Port Victoria (MR=85%). A full susceptibility to P-methyl (0.25%) was observed in all the sites. PBO synergist assay revealed high susceptibility (>98%) to the pyrethroids in all the sites except for Port Victoria (MR=96%, n=100). Whole transcriptomic analysis showed that most of the gene families associated with pyrethroid resistance comprised non-coding RNAs (67%), followed by imipenemase (10%),cytochrome P450s (6%), cuticular proteins (5%), olfactory proteins (4%), glutathione S-transferases (3%), UDP-glycosyltransferases (2%), ATP-binding cassettes (2%) and carboxylesterases(1%). Conclusions This study unveils the molecular basis of insecticide resistance in An. funestus in western Kenya, highlighting for the first time the potential role of non-coding RNAs in pyrethroid resistance. Targeting non-coding RNAs for intervention development could help in insecticide resistance management.

Genomic Profiling of Insecticide Resistance in Malaria Vectors: Insights into Molecular Mechanisms

Malaria control faces challenges from widespread insecticide resistance in major Anopheles species. This study, employing a cross-species approach, integrates RNA-Sequencing, whole-genome sequencing, and microarray data to elucidate drivers of insecticide resistance in Anopheles gambiae complex and An. funestus. Findings show an inverse relationship between genetic diversity and gene expression, with highly expressed genes experiencing stronger purifying selection. These genes cluster physically in the genome, revealing potential coordinated regulation. We identified known and novel candidate insecticide resistance genes, enriched in metabolic, cuticular, and behavioural functions. We also present AnoExpress, a Python package, and an online interface for user-friendly exploration of resistance candidate expression. Despite millions of years of speciation, convergent gene expression responses to insecticidal selection pressures are observed across Anopheles species, providing crucial insights for malaria vector control. This study culminates in a rich dataset that allows us to understand molecular mechanisms, better enabling us to combat insecticide resistance effectively.

Genetic Diversity of Cytochrome P450s CYP6M2 and CYP6P4 Associated with Pyrethroid Resistance in the Major Malaria Vectors Anopheles coluzzii and Anopheles gambiae from Yaoundé, Cameroon

Assessing the genetic diversity of metabolic resistance genes, such as cytochrome P450s, helps to understand the dynamics and evolution of resistance in the field. Here, we analyzed the polymorphisms of CYP6M2 and CYP6P4, associated with pyrethroid resistance in Anopheles coluzzii and Anopheles gambiae, to detect potential resistance markers. Field-caught resistant mosquitos and susceptible lab strains were crossed, and F4 was exposed to permethrin for 15 min and 90 min to discriminate highly susceptible (HS) and highly resistant (HR) mosquitos, respectively. Significant permethrin mortality reduction was observed after pre-exposure to PBO, suggesting the gene involvement of P450s. qPCR analysis revealed significant overexpression of CYP6M2 (FC = 19.57 [95% CI 13.96-25.18] for An. coluzzii; 10.16 [7.86-12.46] for An. gambiae) and CYP6P4 (FC = 6.73 [6.15-7.30] An. coluzzii; 23.62 [26.48-20.76] An. gambiae). Full-gene and ≈1 kb upstream were sequenced. For CYP6M2, the upstream region shows low diversity in HR and HS (overall Hd = 0.49, π = 0.018), whereas the full-gene shows allelic-variation but without evidence of ongoing selection. CYP6P4 upstream region showed a lower diversity in HR (Hd = 0.48) than HS (Hd = 0.86) of An. gambiae. These results highlighted that CYP6P4-associated resistance is potentially driven by modification in upstream region. However, further work is needed to determine the real causative variants that will help design rapid detection tools.

Reversion of CYP450 monooxygenase-mediated acetamiprid larval resistance in dengue fever mosquito, Aedes aegypti L

Aedes-borne diseases are on the rampant rise despite continued application of chemical insecticide-based interventions. The appearance of high degree of insecticide resistance in Aedes species and noxious effects on environment and non-targets have raised further concerns. Among new chemical interventions, neonicotinoids are considered a safe and effective approach. The present study investigated the control potency of acetamiprid and development of resistance in Aedes aegypti larvae; and the involvement of CYP450 monooxygenases in inducing resistance. The early fourth instars of Ae. aegypti parent susceptible strain (PS) were selected with acetamiprid for 15 generations (ACSF strain) increasing the resistance to 19.74-fold in ACSF-10 and 36.71-fold in ACSF-15. The ACSF-10 larvae were assayed with acetamiprid combined with piperonyl butoxide (PBO) in three different ratios (1:1, 1:5 and 1:10) and selected for next five generations with 1:10 combination. Selection with synergized acetamiprid (APSF strains) reversed as well as reduced the rate of resistance development resulting in only 1.35-fold resistance in APSF-15. The APSF strains showed %monooxygenase dependency ranging from 86.71 to 96.72%. The estimation of the monooxygenases levels in parent and selected larvae showed increased monooxygenase level in the ACSF strains by 2.42-2.87-fold. The APSF-15 strains exhibited 57.95% lower enzyme production than ACSF-15 strain. The reduction and reversion of resistance by using PBO and the elevated levels of monooxygenases in ACSF and reduction in APSF strains recommend the involvement of CYP450-mediated mechanism in the development of acetamiprid resistance in Ae. aegypti. These studies could help in devising resistance management strategies in order to preserve the efficiency of pre-existing insecticides.

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.

The cytochrome P450 CYP325A is a major driver of pyrethroid resistance in the major malaria vector Anopheles funestus in Central Africa

The overexpression and overactivity of key cytochrome P450s (CYP450) genes are major drivers of metabolic resistance to insecticides in African malaria vectors such as Anopheles funestus s.s. Previous RNAseq-based transcription analyses revealed elevated expression of CYP325A specific to Central African populations but its role in conferring resistance has not previously been demonstrated. In this study, RT-qPCR consistently confirmed that CYP325A is highly over-expressed in pyrethroid-resistant An. funestus from Cameroon, compared with a control strain and insecticide-unexposed mosquitoes. A synergist bioassay with PBO significantly recovered susceptibility for permethrin and deltamethrin indicating P450-based metabolic resistance. Analyses of the coding sequence of CYP325A Africa-wide detected high-levels of polymorphism, but with no predominant alleles selected by pyrethroid resistance. Geographical amino acid changes were detected notably in Cameroon. In silico homology modelling and molecular docking simulations predicted that CYP325A binds and metabolises type I and type II pyrethroids. Heterologous expression of recombinant CYP325A and metabolic assays confirmed that the most-common Cameroonian haplotype metabolises both type I and type II pyrethroids with depletion rate twice that the of the DR Congo haplotype. Analysis of the 1 kb putative promoter of CYP325A revealed reduced diversity in resistant mosquitoes compared to susceptible ones, suggesting a potential selective sweep in this region. The establishment of CYP325A as a pyrethroid resistance metabolising gene further explains pyrethroid resistance in Central African populations of An. funestus. Our work will facilitate future efforts to detect the causative resistance markers in the promoter region of CYP325A to design field applicable DNA-based diagnostic tools.

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.

A differential expression of pyrethroid resistance genes in the malaria vector Anopheles funestus across Uganda is associated with patterns of gene flow

BACKGROUND

Insecticide resistance is challenging the effectiveness of insecticide-based control interventions to reduce malaria burden in Africa. Understanding the molecular basis of insecticides resistance and patterns of gene flow in major malaria vectors such as Anopheles funestus are important steps for designing effective resistance management strategies. Here, we investigated the association between patterns of genetic structure and expression profiles of genes involved in the pyrethroid resistance in An. funestus across Uganda and neighboring Kenya.

METHODS

Blood-fed mosquitoes An. funestus were collected across the four localities in Uganda and neighboring Kenya. A Microarray-based genome-wide transcription analysis was performed to identify the set of genes associated with permethrin resistance. 17 microsatellites markers were genotyped and used to establish patterns of genetic differentiation.

RESULTS

Microarray-based genome-wide transcription profiling of pyrethroid resistance in four locations across Uganda (Arua, Bulambuli, Lira, and Tororo) and Kenya (Kisumu) revealed that resistance was mainly driven by metabolic resistance. The most commonly up-regulated genes in pyrethroid resistance mosquitoes include cytochrome P450s (CYP9K1, CYP6M7, CYP4H18, CYP4H17, CYP4C36). However, expression levels of key genes vary geographically such as the P450 CYP6M7 [Fold-change (FC) = 115.8 (Arua) vs 24.05 (Tororo) and 16.9 (Kisumu)]. In addition, several genes from other families were also over-expressed including Glutathione S-transferases (GSTs), carboxylesterases, trypsin, glycogenin, and nucleotide binding protein which probably contribute to insecticide resistance across Uganda and Kenya. Genotyping of 17 microsatellite loci in the five locations provided evidence that a geographical shift in the resistance mechanisms could be associated with patterns of population structure throughout East Africa. Genetic and population structure analyses indicated significant genetic differentiation between Arua and other localities (FST>0.03) and revealed a barrier to gene flow between Arua and other areas, possibly associated with Rift Valley.

CONCLUSION

The correlation between patterns of genetic structure and variation in gene expression could be used to inform future interventions especially as new insecticides are gradually introduced.

Genome-wide survey of cytochrome P450 genes in the salmon louse Lepeophtheirus salmonis (Krøyer, 1837)

Background

The salmon louse (Lepeophtheirus salmonis) infests farmed and wild salmonid fishes, causing considerable economic damage to the salmon farming industry. Infestations of farmed salmon are controlled using a combination of non-medicinal approaches and veterinary drug treatments. While L. salmonis has developed resistance to most available salmon delousing agents, relatively little is known about the molecular mechanisms involved. Members of the cytochrome P450 (CYP) superfamily are typically monooxygenases, some of which are involved in the biosynthesis and metabolism of endogenous compounds, while others have central roles in the detoxification of xenobiotics. In terrestrial arthropods, insecticide resistance can be based on the enhanced expression of CYPs. The reported research aimed to characterise the CYP superfamily in L. salmonis and assess its potential roles in drug resistance.

Methods

Lepeophtheirus salmonis CYPs were identified by homology searches of the genome and transcriptome of the parasite. CYP transcript abundance in drug susceptible and multi-resistant L. salmonis was assessed by quantitative reverse transcription PCR, taking into account both constitutive expression and expression in parasites exposed to sublethal levels of salmon delousing agents, ecdysteroids and environmental chemicals.

Results

The above strategy led to the identification of 25 CYP genes/pseudogenes in L. salmonis, making its CYP superfamily the most compact characterised for any arthropod to date. Lepeophtheirus salmonis possesses homologues of a number of arthropod CYP genes with roles in ecdysteroid metabolism, such as the fruit fly genes disembodied, shadow, shade, spook and Cyp18a1. CYP transcript expression did not differ between one drug susceptible and one multi-resistant strain of L. salmonis. Exposure of L. salmonis to emamectin benzoate or deltamethrin caused the transcriptional upregulation of certain CYPs. In contrast, neither ecdysteroid nor benzo[a]pyrene exposure affected CYP transcription significantly.

Conclusions

The parasite L. salmonis is demonstrated to possess the most compact CYP superfamily characterised for any arthropod to date. The complement of CYP genes in L. salmonis includes conserved CYP genes involved in ecdysteroid biosynthesis and metabolism, as well as drug-inducible CYP genes. The present study does not provide evidence for a role of CYP genes in the decreased susceptibility of the multiresistant parasite strain studied.

Characterisation of Anopheles strains used for laboratory screening of new vector control products

BACKGROUND

Insecticides formulated into products that target Anopheles mosquitos have had an immense impact on reducing malaria cases in Africa. However, resistance to currently used insecticides is spreading rapidly and there is an urgent need for alternative public health insecticides. Potential new insecticides must be screened against a range of characterized mosquito strains to identify potential resistance liabilities. The Liverpool School of Tropical Medicine maintains three susceptible and four resistant Anopheles strains that are widely used for screening for new insecticides. The properties of these strains are described in this paper.

METHODS

WHO tube susceptibility bioassays were used for colony selection and to screen for resistance to the major classes of public health insecticides. Topical and tarsal contact bioassays were used to produce dose response curves to assess resistance intensity. Bioassays with the synergist piperonyl butoxide were also performed. Taqman™ assays were used to screen for known target site resistance alleles (kdr and ace-1). RT-qPCR was used to quantify expression of genes associated with pyrethroid resistance.

RESULTS

Pyrethroid selection pressure has maintained resistance to this class in all four resistant strains. Some carbamate and organophosphate resistance has been lost through lack of exposure to these insecticide classes. The Anopheles gambiae (sensu lato) strains, VK7 2014, Banfora M and Tiassalé 13 have higher levels of pyrethroid resistance than the An. funestus FUMOZ-R strain. Elevated expression of P450s is found in all four strains and the 1014F kdr mutation is present in all three An. gambiae strains at varying frequencies. Tarsal contact data and overexpression of CYP4G16 and SAP2 suggest penetration barriers and/or sequestration also confer resistance in Banfora M.

CONCLUSIONS

Continual selection with deltamethrin has maintained a stable pyrethroid-resistant phenotype over many generations. In conjunction with a standardized rearing regime, this ensures quality control of strains over time allowing for robust product comparison and selection of optimal products for further development. The identification of multiple mechanisms underpinning insecticide resistance highlights the importance of screening new compounds against a range of mosquito strains.

Cis-regulatory CYP6P9b P450 variants associated with loss of insecticide-treated bed net efficacy against Anopheles funestus

Elucidating the genetic basis of metabolic resistance to insecticides in malaria vectors is crucial to prolonging the effectiveness of insecticide-based control tools including long lasting insecticidal nets (LLINs). Here, we show that cis-regulatory variants of the cytochrome P450 gene, CYP6P9b, are associated with pyrethroid resistance in the African malaria vector Anopheles funestus. A DNA-based assay is designed to track this resistance that occurs near fixation in southern Africa but not in West/Central Africa. Applying this assay we demonstrate, using semi-field experimental huts, that CYP6P9b-mediated resistance associates with reduced effectiveness of LLINs. Furthermore, we establish that CYP6P9b combines with another P450, CYP6P9a, to additively exacerbate the reduced efficacy of insecticide-treated nets. Double homozygote resistant mosquitoes (RR/RR) significantly survive exposure to insecticide-treated nets and successfully blood feed more than other genotypes. This study provides tools to track and assess the impact of multi-gene driven metabolic resistance to pyrethroids, helping improve resistance management.

Bed nets treated with insecticides have been instrumental in reducing malaria mortality, but insecticide resistance is on the rise. Here, Mugenzi et al. identify genetic variants in the P450 gene CYP6P9b of Anopheles funestus that associate with insecticide resistance and develop a PCR-based diagnostic assay to help identify pyrethroid-resistant strains.

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.

Candidate-gene based GWAS identifies reproducible DNA markers for metabolic pyrethroid resistance from standing genetic variation in East African Anopheles gambiae

Metabolic resistance to pyrethroid insecticides is widespread in Anopheles mosquitoes and is a major threat to malaria control. DNA markers would aid predictive monitoring of resistance, but few mutations have been discovered outside of insecticide-targeted genes. Isofemale family pools from a wild Ugandan Anopheles gambiae population, from an area where operational pyrethroid failure is suspected, were genotyped using a candidate-gene enriched SNP array. Resistance-associated SNPs were detected in three genes from detoxification superfamilies, in addition to the insecticide target site (the Voltage Gated Sodium Channel gene, Vgsc). The putative associations were confirmed for two of the marker SNPs, in the P450 Cyp4j5 and the esterase Coeae1d by reproducible association with pyrethroid resistance in multiple field collections from Uganda and Kenya, and together with the Vgsc-1014S (kdr) mutation these SNPs explained around 20% of variation in resistance. Moreover, the >20 Mb 2La inversion also showed evidence of association with resistance as did environmental humidity. Sequencing of Cyp4j5 and Coeae1d detected no resistance-linked loss of diversity, suggesting selection from standing variation. Our study provides novel, regionally-validated DNA assays for resistance to the most important insecticide class, and establishes both 2La karyotype variation and humidity as common factors impacting the resistance phenotype.

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.

Genomic Footprints of Selective Sweeps from Metabolic Resistance to Pyrethroids in African Malaria Vectors Are Driven by Scale up of Insecticide-Based Vector Control

Insecticide resistance in mosquito populations threatens recent successes in malaria prevention. Elucidating patterns of genetic structure in malaria vectors to predict the speed and direction of the spread of resistance is essential to get ahead of the ‘resistance curve’ and to avert a public health catastrophe. Here, applying a combination of microsatellite analysis, whole genome sequencing and targeted sequencing of a resistance locus, we elucidated the continent-wide population structure of a major African malaria vector, Anopheles funestus. We identified a major selective sweep in a genomic region controlling cytochrome P450-based metabolic resistance conferring high resistance to pyrethroids. This selective sweep occurred since 2002, likely as a direct consequence of scaled up vector control as revealed by whole genome and fine-scale sequencing of pre- and post-intervention populations. Fine-scaled analysis of the pyrethroid resistance locus revealed that a resistance-associated allele of the cytochrome P450 monooxygenase CYP6P9a has swept through southern Africa to near fixation, in contrast to high polymorphism levels before interventions, conferring high levels of pyrethroid resistance linked to control failure. Population structure analysis revealed a barrier to gene flow between southern Africa and other areas, which may prevent or slow the spread of the southern mechanism of pyrethroid resistance to other regions. By identifying a genetic signature of pyrethroid-based interventions, we have demonstrated the intense selective pressure that control interventions exert on mosquito populations. If this level of selection and spread of resistance continues unabated, our ability to control malaria with current interventions will be compromised.

Malaria control currently relies heavily on insecticide-based vector control interventions. Unfortunately, resistance to insecticides threatens the continued effectiveness of these measures. Metabolic resistance, caused by increased detoxification of insecticides, presents the greatest threat to vector control, yet it remains unclear how these mechanisms are linked to underlying genetic changes driven by the massive selection pressure from these interventions, such as the widespread use of Long Lasting Insecticide Nets (LLINs) across Africa. Therefore, understanding the direction and speed at which this operationally important form of resistance spreads through mosquito populations is essential if we are to get ahead of the ‘resistance curve’ and avert a public health catastrophe. Here, using microsatellite markers, whole genome sequencing and fine-scale sequencing at a major resistance locus, we elucidated the Africa-wide population structure of Anopheles funestus, a major African malaria vector, and detected a strong selective sweep occurring in a genomic region controlling cytochrome P450-based metabolic pyrethroid resistance in this species. Furthermore, we demonstrated that this selective sweep is driven by the scale-up of insecticide-based malaria control in Africa, highlighting the risk that if this level of selection and spread of resistance continues unabated, our ability to control malaria with current interventions will be compromised.

The P450 CYP6Z1 confers carbamate/pyrethroid cross-resistance in a major African malaria vector beside a novel carbamate-insensitive N485I acetylcholinesterase-1 mutation

Carbamates are increasingly used for vector control notably in areas with pyrethroid resistance. However, a cross‐resistance between these insecticides in major malaria vectors such as Anopheles funestus could severely limit available resistance management options. Unfortunately, the molecular basis of such cross‐resistance remains uncharacterized in An. funestus, preventing effective resistance management. Here, using a genomewide transcription profiling, we revealed that metabolic resistance through upregulation of cytochrome P450 genes is driving carbamate resistance. The P450s CYP6P9a,CYP6P9b and CYP6Z1 were the most upregulated detoxification genes in the multiple resistant mosquitoes. However, in silico docking simulations predicted CYP6Z1 to metabolize both pyrethroids and carbamates, whereas CYP6P9a and CYP6P9b were predicted to metabolize only the pyrethroids. Using recombinant enzyme metabolism and inhibition assays, we demonstrated that CYP6Z1 metabolizes bendiocarb and pyrethroids, whereas CYP6P9a and CYP6P9b metabolize only the pyrethroids. Other upregulated gene families in resistant mosquitoes included several cuticular protein genes suggesting a possible reduced penetration resistance mechanism. Investigation of the target‐site resistance in acetylcholinesterase 1 (ace‐1) gene detected and established the association between the new N485I mutation and bendiocarb resistance (odds ratio 7.3; P < 0.0001). The detection of multiple haplotypes in single mosquitoes after cloning suggested the duplication of ace‐1. A TaqMan genotyping of the N485I in nine countries revealed that the mutation is located only in southern Africa with frequency of 10–15% suggesting its recent occurrence. These findings will help in monitoring the spread and evolution of carbamate resistance and improve the design of effective resistance management strategies to control this malaria vector.

Identification, Validation, and Application of Molecular Diagnostics for Insecticide Resistance in Malaria Vectors

Insecticide resistance is a major obstacle to control of Anopheles malaria mosquitoes in sub-Saharan Africa and requires an improved understanding of the underlying mechanisms. Efforts to discover resistance genes and DNA markers have been dominated by candidate gene and quantitative trait locus studies of laboratory strains, but with greater availability of genome sequences a shift toward field-based agnostic discovery is anticipated. Mechanisms evolve continually to produce elevated resistance yielding multiplicative diagnostic markers, co-screening of which can give high predictive value. With a shift toward prospective analyses, identification and screening of resistance marker panels will boost monitoring and programmatic decision making.

Identification of QTLs Conferring Resistance to Deltamethrin in Culex pipiens pallens

Culex pipiens pallens is the most abundant Culex mosquito species in northern China and is an important vector of bancroftian filariasis and West Nile virus. Deltamethrin is an insecticide that is widely used for mosquito control, however resistance to this and other insecticides has become a major challenge in the control of vector-borne diseases that appear to be inherited quantitatively. Furthermore, the genetic basis of insecticide resistance remains poorly understood. In this study, quantitative trait loci (QTL) mapping of resistance to deltamethrin was conducted in F₂ intercross segregation populations using bulked segregation analysis (BSA) and amplified fragment length polymorphism markers (AFLP) in Culex pipiens pallens. A genetic linkage map covering 381 cM was constructed and a total of seven QTL responsible for resistance to deltamethrin were detected by composite interval mapping (CIM), which explained 95% of the phenotypic variance. The major QTL in linkage group 2 accounted for 62% of the variance and is worthy of further study. 12 AFLP markers in the map were cloned and the genomic locations of these marker sequences were determined by applying the Basic Local Alignment Search Tool (BLAST) tool to the genome sequence of the closely related Culex quinquefasciatus. Our results suggest that resistance to deltamethrin is a quantitative trait under the control of a major QTL in Culex pipiens pallens. Cloning of related AFLP markers confirm the potential utility for anchoring the genetic map to the physical map. The results provide insight into the genetic architecture of the trait.

Target-site resistance mutations (kdr and RDL), but not metabolic resistance, negatively impact male mating competiveness in the malaria vector Anopheles gambiae

The implementation of successful insecticide resistance management strategies for malaria control is currently hampered by poor understanding of the fitness cost of resistance on mosquito populations, including their mating competiveness. To fill this knowledge gap, coupled and uncoupled Anopheles gambiae s.l. males (all M form (Anopheles coluzzii)) were collected from mating swarms in Burkina Faso. This multiple insecticide resistant population exhibited high 1014F kdr^R allele frequencies (>60%) and RDL^R (>80%) in contrast to the Ace-1^R allele (<6%). Kdr heterozygote males were more likely to mate than homozygote resistant (OR=2.36; P<0.001), suggesting a negative impact of kdr on An. coluzzii mating ability. Interestingly, heterozygote males were also more competitive than homozygote susceptible (OR=3.26; P=0.006), suggesting a heterozygote advantage effect. Similarly, heterozygote RDL^R/RDL^S were also more likely to mate than homozygote-resistant males (OR=2.58; P=0.007). Furthermore, an additive mating disadvantage was detected in male homozygotes for both kdr/RDL-resistant alleles. In contrast, no fitness difference was observed for the Ace-1 mutation. Comparative microarray-based genome-wide transcription analysis revealed that metabolic resistance did not significantly alter the mating competitiveness of male An. coluzzii mosquitoes. Indeed, no significant difference of expression levels was observed for the main metabolic resistance genes, suggesting that metabolic resistance has a limited impact on male mating competiveness. In addition, specific gene classes/GO terms associated with mating process were detected including sensory perception and peroxidase activity. The detrimental impact of insecticide resistance on mating competiveness observed here suggests that resistance management strategies such as insecticide rotation could help reverse the resistance, if implemented early.

The colonization of pyrethroid resistant strain from wild Anopheles sinensis, the major Asian malaria vector

Background

Anopheles sinensis is one of the most important malaria vectors in Asian countries. The rapid spread of insecticide resistance has become a major obstacle for insecticide-based strategies for vector control. Therefore, it is necessary to prepare an insecticide-resistant strain of An. sinensis to further understand the insecticide resistance mechanisms in this species to facilitate genetic approaches to targeting the insecticide-resistant population of this important malaria vector.

Methods

An. sinensis mosquitoes were collected from regions where pyrethroid resistance had been reported. The mosquitoes were subjected to continuous pyrethroid selection after species confirmation, and the forced copulation method was used to increase the mating rate. In addition, the knockdown-resistance (kdr) mutation frequencies of each generation of An. sinensis were measured; and the metabolic enzyme activities of cytochrome P450 monoxygenases (P450s) and glutathione S-transferases (GSTs) were detected.

Results

The identification of field-captured An. sinensis was confirmed by both morphological and molecular methods. The population of An. sinensis exhibited stable resistance to pyrethroid after continuous generations of pyrethroid selection in the laboratory with high kdr mutation frequencies; and elevated levels of both P450s and GSTs were significantly found in field selected populations comparing with the laboratory susceptible strain. So far, the colonised strain has reached its eleventh generation and culturing well in the laboratory.

Conclusions

We colonised a pyrethroid-resistant population of An. sinensis in the laboratory, which provides a fundamental model for genetic studies of this important malaria vector.

Mutations in the voltage-gated sodium channel gene of anophelines and their association with resistance to pyrethroids - a review

Constant and extensive use of chemical insecticides has created a selection pressure and favored resistance development in many insect species worldwide. One of the most important pyrethroid resistance mechanisms is classified as target site insensitivity, due to conformational changes in the target site that impair a proper binding of the insecticide molecule. The voltage-gated sodium channel (NaV) is the target of pyrethroids and DDT insecticides, used to control insects of medical, agricultural and veterinary importance, such as anophelines. It has been reported that the presence of a few non-silent point mutations in the NaV gene are associated with pyrethroid resistance, termed as 'kdr' (knockdown resistance) for preventing the knockdown effect of these insecticides. The presence of these mutations, as well as their effects, has been thoroughly studied in Anopheles mosquitoes. So far, kdr mutations have already been detected in at least 13 species (Anopheles gambiae, Anopheles arabiensis, Anopheles sinensis, Anopheles stephensi, Anopheles subpictus, Anopheles sacharovi, Anopheles culicifacies, Anopheles sundaicus, Anopheles aconitus, Anopheles vagus, Anopheles paraliae, Anopheles peditaeniatus and Anopheles albimanus) from populations of African, Asian and, more recently, American continents. Seven mutational variants (L1014F, L1014S, L1014C, L1014W, N1013S, N1575Y and V1010L) were described, with the highest prevalence of L1014F, which occurs at the 1014 site in NaV IIS6 domain. The increase of frequency and distribution of kdr mutations clearly shows the importance of this mechanism in the process of pyrethroid resistance. In this sense, several species-specific and highly sensitive methods have been designed in order to genotype individual mosquitoes for kdr in large scale, which may serve as important tolls for monitoring the dynamics of pyrethroid resistance in natural populations. We also briefly discuss investigations concerning the course of Plasmodium infection in kdr individuals. Considering the limitation of insecticides available for employment in public health campaigns and the absence of a vaccine able to brake the life cycle of the malaria parasites, the use of pyrethroids is likely to remain as the main strategy against mosquitoes by either indoor residual spraying (IR) and insecticide treated nets (ITN). Therefore, monitoring insecticide resistance programs is a crucial need in malaria endemic countries.

High frequency of kdr L1014F is associated with pyrethroid resistance in Anopheles coluzzii in Sudan savannah of northern Nigeria

BACKGROUND

Malaria burden is high in Nigeria, yet information on the major mosquito vectors is lacking especially in the Sudan savannah region of the country. In order to facilitate the design of future insecticide-based control interventions in the region, this study has established the resistance profile of An. gambiae s.l. populations in two northern Nigeria locations and assessed the contribution of target site resistance mutations.

METHODS

Larval collection was conducted in two localities in Sudan savannah (Bunkure and Auyo) of northern Nigeria between 2009 and 2011, from which resulting adult, female mosquitoes were used for insecticides bioassays with deltamethrin, lambda-cyhalothrin, DDT and malathion. The mosquitoes were identified to species level and molecular forms and then genotyped for the presence of L1014F-kdr, L1014S-kdr and ace-1R mutations.

RESULTS

WHO bioassays revealed that An. gambiae s.l. from both localities were highly resistant to lambda-cyhalothrin and DDT, but only moderately resistant to deltamethrin. Full susceptibility was observed to malathion. An. gambiae, M form (now An. coluzzii), was predominant over An. arabiensis in Auyo and was more resistant to lambda-cyhalothrin than An. arabiensis. No 'S' form (An. gambiae s.s.) was detected. A high frequency of 1014 F mutation (80.1%) was found in An. coluzzii in contrast to An. arabiensis (13.5%). The presence of the 1014 F kdr allele was significantly associated with resistance to lambda-cyhalothrin in An. coluzzii (OR = 9.85; P < 0.001) but not in An. arabiensis. The L1014S-kdr mutation was detected in a single An. arabiensis mosquito while no ace-1R mutation was found in any of the mosquitoes analysed.

CONCLUSIONS

The predominance of An. coluzzii and its resistance profile to main insecticides described in this study can guide the implementation of appropriate vector control interventions in this region of Nigeria where such information was previously lacking.

The dynamics of pyrethroid resistance in Anopheles arabiensis from Zanzibar and an assessment of the underlying genetic basis

Background

The emergence of pyrethroid resistance in the malaria vector, Anopheles arabiensis, threatens to undermine the considerable gains made towards eliminating malaria on Zanzibar. Previously, resistance was restricted to the island of Pemba while mosquitoes from Unguja, the larger of the two islands of Zanzibar, were susceptible. Here, we characterised the mechanism(s) responsible for resistance on Zanzibar using a combination of gene expression and target-site mutation assays.

Methods

WHO resistance bioassays were conducted using 1-5d old adult Anopheles gambiae s.l. collected between 2011 and 2013 across the archipelago. Synergist assays with the P450 inhibitor piperonyl-butoxide were performed in 2013. Members of the An. gambiae complex were PCR-identified and screened for target-site mutations (kdr and Ace-1). Gene expression in pyrethroid resistant An. arabiensis from Pemba was analysed using whole-genome microarrays.

Results

Pyrethroid resistance is now present across the entire Zanzibar archipelago. Survival to the pyrethroid lambda-cyhalothrin in bioassays conducted in 2013 was 23.5-54.3% on Unguja and 32.9-81.7% on Pemba. We present evidence that resistance is mediated, in part at least, by elevated P450 monoxygenases. Whole-genome microarray scans showed that the most enriched gene terms in resistant An. arabiensis from Pemba were associated with P450 activity and synergist assays with PBO completely restored susceptibility to pyrethroids in both islands. CYP4G16 was the most consistently over-expressed gene in resistant mosquitoes compared with two susceptible strains from Unguja and Dar es Salaam. Expression of this P450 is enriched in the abdomen and it is thought to play a role in hydrocarbon synthesis. Microarray and qPCR detected several additional genes putatively involved in this pathway enriched in the Pemba pyrethroid resistant population and we hypothesise that resistance may be, in part, related to alterations in the structure of the mosquito cuticle. None of the kdr target-site mutations, associated with pyrethroid/DDT resistance in An. gambiae elsewhere in Africa, were found on the islands.

Conclusion

The consequences of this resistance phenotype are discussed in relation to future vector control strategies on Zanzibar to support the ongoing malaria elimination efforts on the islands.

Development of reference transcriptomes for the major field insect pests of cowpea: a toolbox for insect pest management approaches in west Africa

Cowpea is a widely cultivated and major nutritional source of protein for many people that live in West Africa. Annual yields and longevity of grain storage is greatly reduced by feeding damage caused by a complex of insect pests that include the pod sucking bugs, Anoplocnemis curvipes Fabricius (Hemiptera: Coreidae) and Clavigralla tomentosicollis Stål (Hemiptera: Coreidae); as well as phloem-feeding cowpea aphids, Aphis craccivora Koch (Hemiptera: Aphididae) and flower thrips, Megalurothrips sjostedti Trybom (Thysanoptera: Thripidae). Efforts to control these pests remain a challenge and there is a need to understand the structure and movement of these pest populations in order to facilitate the development of integrated pest management strategies (IPM). Molecular tools have the potential to help facilitate a better understanding of pest populations. Towards this goal, we used 454 pyrosequencing technology to generate 319,126, 176,262, 320,722 and 227,882 raw reads from A. curvipes, A. craccivora, C. tomentosicollis and M. sjostedti, respectively. The reads were de novo assembled into 11,687, 7,647, 10,652 and 7,348 transcripts for A. curvipes, A. craccivora, C. tomentosicollis and M. sjostedti, respectively. Functional annotation of the resulting transcripts identified genes putatively involved in insecticide resistance, pathogen defense and immunity. Additionally, sequences that matched the primary aphid endosymbiont, Buchnera aphidicola, were identified among A. craccivora transcripts. Furthermore, 742, 97, 607 and 180 single nucleotide polymorphisms (SNPs) were respectively predicted among A. curvipes, A. craccivora, C. tomentosicollis and M. sjostedti transcripts, and will likely be valuable tools for future molecular genetic marker development. These results demonstrate that Roche 454-based transcriptome sequencing could be useful for the development of genomic resources for cowpea pest insects in West Africa.

How many genetic options for evolving insecticide resistance in heliothine and spodopteran pests?

The widely accepted paradigm for the development of insecticide resistance in field populations of insects is of selection for one or a very few genes of major effect. Limited genetic mapping data for organophosphate and pyrethroid resistance in heliothine and spodopteran pests generally agrees with this paradigm. However, other biochemical and transcriptomic data suggest a more complex set of changes in multiple P450 and esterase gene/enzyme systems in resistant strains of these species. We discuss possible explanations for this paradox, including the likely embedding of these genes in regulatory cascades and emerging evidence for their arrangement in large clusters of closely related genes. We conclude that there could indeed be an unusually large number of genetic options for evolving resistance in these species.

Dissecting the mechanisms responsible for the multiple insecticide resistance phenotype in Anopheles gambiae s.s., M form, from Vallée du Kou, Burkina Faso

With the exception of target site mutations, insecticide resistance mechanisms in the principle malaria vector Anopheles gambiae, remains largely uncharacterized in Burkina Faso.

Here we detected high prevalence of resistance in Vallée du Kou (VK) to pyrethroids, DDT and dieldrin, moderate level for carbamates and full susceptibility to organophosphates. High frequencies of L1014F kdr (75%) and Rdl (87%) mutations were observed showing strong correlation with pyrethroids/DDT and dieldrin resistance. The frequency of ace1^R mutation was low even in carbamate resistant mosquitoes. Microarray analysis identified genes significantly over-transcribed in VK. These include the cytochrome P450 genes, CYP6P3 and CYP6Z2, previously associated with pyrethroid resistance. Gene Ontology (GO) enrichment analysis suggested that elevated neurotransmitter activity is associated with resistance, with the over-transcription of target site resistance genes such as acetylcholinesterase and the GABA receptor. A rhodopsin receptor gene previously associated with pyrethroid resistance in Culex pipiens pallens was also over-transcribed in VK.

This study highlights the complex network of mechanisms conferring multiple resistance in malaria vectors and such information should be taken into account when designing and implementing resistance control strategies.

Genetic mapping identifies a major locus spanning P450 clusters associated with pyrethroid resistance in kdr-free Anopheles arabiensis from Chad

Prevention of malaria transmission throughout much of Africa is dependent on bednets that are impregnated with pyrethroid insecticides. Anopheles arabiensis is the major malaria vector in Chad and efforts to control this vector are threatened by the emergence of pyrethroid resistance. WHO bioassays revealed that An. arabiensis from Ndjamena is resistant to pyrethroids and dichlorodiphenyltrichloroethane (DDT) but fully susceptible to carbamates and organophosphates. No 1014F or 1014S kdr alleles were detected in this population. To determine the mechanisms that are responsible for resistance, genetic crosses were established between the Ndja strain and an insecticide susceptible population from Mozambique. Resistance was inherited as an autosomal trait and quantitative trait locus (QTL) mapping identified a single major locus on chromosome 2R, which explained 24.4% of the variance in resistance. This QTL is enriched in P450 genes including 25 cytochrome P450s in total. One of these, Cyp6p4 is 22-fold upregulated in the Ndja strain compared with the susceptible. Piperonyl butoxide (PBO) synergist and biochemical assays further support a role for P450s in conferring pyrethroid resistance in this population.

Role of cytochrome P450s in insecticide resistance: impact on the control of mosquito-borne diseases and use of insecticides on Earth

The fight against diseases spread by mosquitoes and other insects has enormous environmental, economic and social consequences. Chemical insecticides remain the first line of defence but the control of diseases, especially malaria and dengue fever, is being increasingly undermined by insecticide resistance. Mosquitoes have a large repertoire of P450s (over 100 genes). By pinpointing the key enzymes associated with insecticide resistance we can begin to develop new tools to aid the implementation of control interventions and reduce their environmental impact on Earth. Recent technological advances are helping us to build a functional profile of the P450 determinants of insecticide metabolic resistance in mosquitoes. Alongside, the cross-responses of mosquito P450s to insecticides and pollutants are also being investigated. Such research will provide the means to produce diagnostic tools for early detection of P450s linked to resistance. It will also enable the design of new insecticides with optimized efficacy in different environments.

Directionally selected cytochrome P450 alleles are driving the spread of pyrethroid resistance in the major malaria vector Anopheles funestus

Pyrethroid insecticides are critical for malaria control in Africa. However, resistance to this insecticide class in the malaria vector Anopheles funestus is spreading rapidly across Africa, threatening the success of ongoing and future malaria control programs. The underlying resistance mechanisms driving the spread of this resistance in wild populations remain largely unknown. Here, we show that increased expression of two tandemly duplicated P450 genes, CYP6P9a and CYP6P9b, is the main mechanism driving pyrethroid resistance in Malawi and Mozambique, two southern African countries where this insecticide class forms the mainstay of malaria control. Genome-wide transcription analysis using microarray and quantitative RT-PCR consistently revealed that CYP6P9a and CYP6P9b are the two genes most highly overexpressed (>50-fold; q < 0.01) in permethrin-resistant mosquitoes. Transgenic expression of CYP6P9a and CYP6P9b in Drosophila melanogaster demonstrated that elevated expression of either of these genes confers resistance to both type I (permethrin) and type II (deltamethrin) pyrethroids. Functional characterization of recombinant CYP6P9b confirmed that this protein metabolized both type I (permethrin and bifenthrin) and type II (deltamethrin and Lambda-cyhalothrin) pyrethroids but not DDT. Variability analysis identified that a single allele of each of these genes is predominantly associated with pyrethroid resistance in field populations from both countries, which is suggestive of a single origin of this resistance that has since spread across the region. Urgent resistance management strategies should be implemented in this region to limit a further spread of this resistance and minimize its impact on the success of ongoing malaria control programs.

Comparative transcriptome analyses of deltamethrin-resistant and -susceptible Anopheles gambiae mosquitoes from Kenya by RNA-Seq

Malaria causes more than 300 million clinical cases and 665,000 deaths each year, and the majority of the mortality and morbidity occurs in sub-Saharan Africa. Due to the lack of effective vaccines and wide-spread resistance to antimalarial drugs, mosquito control is the primary method of malaria prevention and control. Currently, malaria vector control relies on the use of insecticides, primarily pyrethroids. The extensive use of insecticides has imposed strong selection pressures for resistance in the mosquito populations. Consequently, resistance to pyrethroids in Anopheles gambiae, the main malaria vector in sub-Saharan Africa, has become a major obstacle for malaria control. A key element of resistance management is the identification of resistance mechanisms and subsequent development of reliable resistance monitoring tools. Field-derived An. gambiae from Western Kenya were phenotyped as deltamethrin-resistant or -susceptible by the standard WHO tube test, and their expression profile compared by RNA-seq. Based on the current annotation of the An. gambiae genome, a total of 1,093 transcripts were detected as significantly differentially accumulated between deltamethrin-resistant and -susceptible mosquitoes. These transcripts are distributed over the entire genome, with a large number mapping in QTLs previously linked to pyrethorid resistance, and correspond to heat-shock proteins, metabolic and transport functions, signal transduction activities, cytoskeleton and others. The detected differences in transcript accumulation levels between resistant and susceptible mosquitoes reflect transcripts directly or indirectly correlated with pyrethroid resistance. RNA-seq data also were used to perform a de-novo Cufflinks assembly of the An. gambiae genome.

Positional cloning of rp2 QTL associates the P450 genes CYP6Z1, CYP6Z3 and CYP6M7 with pyrethroid resistance in the malaria vector Anopheles funestus

Pyrethroid resistance in Anopheles funestus is threatening malaria control in Africa. Elucidation of underlying resistance mechanisms is crucial to improve the success of future control programs. A positional cloning approach was used to identify genes conferring resistance in the uncharacterised rp2 quantitative trait locus (QTL) previously detected in this vector using F6 advanced intercross lines (AIL). A 113 kb BAC clone spanning rp2 was identified and sequenced revealing a cluster of 15 P450 genes and one salivary protein gene (SG7-2). Contrary to A. gambiae, AfCYP6M1 is triplicated in A. funestus, while AgCYP6Z2 orthologue is absent. Five hundred and sixty-five new single nucleotide polymorphisms (SNPs) were identified for genetic mapping from rp2 P450s and other genes revealing high genetic polymorphisms with one SNP every 36 bp. A significant genotype/phenotype association was detected for rp2 P450s but not for a cluster of cuticular protein genes previously associated with resistance in A. gambiae. QTL mapping using F6 AIL confirms the rp2 QTL with an increase logarithm of odds score of 5. Multiplex gene expression profiling of 15 P450s and other genes around rp2 followed by individual validation using qRT–PCR indicated a significant overexpression in the resistant FUMOZ-R strain of the P450s AfCYP6Z1, AfCYP6Z3, AfCYP6M7 and the glutathione-s-transferase GSTe2 with respective fold change of 11.2, 6.3, 5.5 and 2.8. Polymorphisms analysis of AfCYP6Z1 and AfCYP6Z3 identified amino acid changes potentially associated with resistance further indicating that these genes are controlling the pyrethroid resistance explained by the rp2 QTL. The characterisation of this rp2 QTL significantly improves our understanding of resistance mechanisms in A. funestus.

Identification and validation of a gene causing cross-resistance between insecticide classes in Anopheles gambiae from Ghana

In the last decade there have been marked reductions in malaria incidence in sub-Saharan Africa. Sustaining these reductions will rely upon insecticides to control the mosquito malaria vectors. We report that in the primary African malaria vector, Anopheles gambiae sensu stricto, a single enzyme, CYP6M2, confers resistance to two classes of insecticide. This is unique evidence in a disease vector of cross-resistance associated with a single metabolic gene that simultaneously reduces the efficacy of two of the four classes of insecticide routinely used for malaria control. The gene-expression profile of a highly DDT-resistant population of A. gambiae s.s. from Ghana was characterized using a unique whole-genome microarray. A number of genes were significantly overexpressed compared with two susceptible West African colonies, including genes from metabolic families previously linked to insecticide resistance. One of the most significantly overexpressed probe groups (false-discovery rate-adjusted P < 0.0001) belonged to the cytochrome P450 gene CYP6M2. This gene is associated with pyrethroid resistance in wild A. gambiae s.s. populations) and can metabolize both type I and type II pyrethroids in recombinant protein assays. Using in vitro assays we show that recombinant CYP6M2 is also capable of metabolizing the organochlorine insecticide DDT in the presence of solubilizing factor sodium cholate.

Exploring mechanisms of multiple insecticide resistance in a population of the malaria vector Anopheles funestus in Benin

Background

The insecticide resistance status of the malaria vector Anopheles funestus and the underlying resistance mechanisms remain uncharacterised in many parts of Africa, notably in Benin, West Africa. To fill this gap in our knowledge, we assessed the susceptibility status of a population of this species in Pahou, Southern Benin and investigated the potential resistance mechanisms.

Methodology/Principal Findings

WHO bioassays revealed a multiple resistance profile for An. funestus in Pahou. This population is highly resistant to DDT with no mortality in females after 1h exposure to 4%DDT. Resistance was observed against the Type I pyrethroid permethrin and the carbamate bendiocarb. A moderate resistance was detected against deltamethrin (type II pyrethroids). A total susceptibility was observed against malathion, an organophosphate. Pre-exposure to PBO did not change the mortality rates for DDT indicating that cytochrome P450s play no role in DDT resistance in Pahou. No L1014F kdr mutation was detected but a correlation between haplotypes of two fragments of the Voltage-Gated Sodium Channel gene and resistance was observed suggesting that mutations in other exons may confer the knockdown resistance in this species. Biochemical assays revealed elevated levels of GSTs and cytochrome mono-oxygenases in Pahou. No G119S mutation and no altered acetylcholinesterase gene were detected in the Pahou population. qPCR analysis of five detoxification genes revealed that the GSTe2 is associated to the DDT resistance in this population with a significantly higher expression in DDT resistant samples. A significant over-expression of CYP6P9a and CYP6P9b previously associated with pyrethroid resistance was also seen but at a lower fold change than in southern Africa.

Conclusion

The multiple insecticide resistance profile of this An. funestus population in Benin shows that more attention should be paid to this important malaria vector for the implementation and management of current and future malaria vector control programs in this country.

Identification and distribution of a GABA receptor mutation conferring dieldrin resistance in the malaria vector Anopheles funestus in Africa

Growing problems of pyrethroid resistance in Anopheles funestus have intensified efforts to identify alternative insecticides. Many agrochemicals target the GABA receptors, but cross-resistance from dieldrin resistance may preclude their introduction. Dieldrin resistance was detected in An. funestus populations from West (Burkina Faso) and central (Cameroon) Africa, but populations from East (Uganda) and Southern Africa (Mozambique and Malawi) were fully susceptible to this insecticide. Partial sequencing of the dieldrin target site, the γ-aminobutyric acid (GABA) receptor, identified two amino acid substitutions, A296S and V327I. The A296S mutation has been associated with dieldrin resistance in other species. The V327I mutations was detected in the resistant sample from Burkina Faso and Cameroon and consistently associated with the A296S substitution. The full-length of the An. funestus GABA-receptor gene, amplified by RT-PCR, generated a sequence of 1674 bp encoding 557 amino acid of the protein in An. funestus with 98% similarity to that of Anopheles gambiae. Two diagnostic assays were developed to genotype the A296S mutation (pyrosequencing and PCR-RFLP), and use of these assays revealed high frequency of the resistant allele in Burkina Faso (60%) and Cameroon (82%), moderate level in Benin (16%) while low frequency or absence of the mutation was observed respectively in Uganda (7.5%) or 0% in Malawi and Mozambique. The distribution of the Rdl(R) mutation in An. funestus populations in Africa suggests extensive barriers to gene flow between populations from different regions.

Cytogenetic map for Anopheles nili: application for population genetics and comparative physical mapping

Anopheles nili is one of the major malaria vectors in Africa with a wide geographic distribution. However, the taxonomic and population genetic studies on this species are scarce. New research tools are urgently needed to genetically characterize this important malaria vector. In this study, a high-resolution cytogenetic map was developed for An. nili polytene chromosomes. Chromosomes were straightened and subdivided into 46 numbered divisions according to the banding pattern. Population analysis of An. nili females collected in Burkina Faso revealed the presence of two highly polymorphic inversions on the 2R chromosomal arm. A statistically significant departure from Hardy-Weinberg equilibrium due to a deficit in heterozygotes was detected for inversion 2Rb. To determine chromosome homologies and gene order conservation between An. nili and other major malaria vectors, PCR probes based on the An. gambiae coding sequences were mapped to An. nili chromosomes. Comparative mapping demonstrated that An. nili chromosomes have an An. stephensi-like arm association and that whole-arm translocations and paracentric inversions were the major types of rearrangement in evolution of these mosquitoes. The minimum number of fixed inversions among An. nili, An. gambiae, and An. stephensi was calculated using the Multiple Genome Rearrangements (MGR), Genome Rearrangements In Man and Mouse (GRIMM), and Sorting Permutation by Reversals and block-INterchanGes (SPRING) programs. The data suggest that the An. nili is, at least, as diverged from An. gambiae as An. stephensi. We provide evidence that 2La/a arrangement of An. gambiae is present in outgroup species An. nili and An. stephensi confirming the ancestral status of the 2La inversion in the An. gambiae complex. Availability of the new polytene chromosome map, polymorphic inversions, and physically mapped DNA markers for An. nili will further stimulate population genetic, taxonomic, and genomic studies of this neglected malaria vector.

Cuticle thickening associated with pyrethroid resistance in the major malaria vector Anopheles funestus

Background

Malaria in South Africa is primarily transmitted by Anopheles funestus Giles. Resistance to pyrethroid insecticides in An. funestus in northern Kwazulu/Natal, South Africa, and in neighbouring areas of southern Mozambique enabled populations of this species to increase their ranges into areas where pyrethroids were being exclusively used for malaria control. Pyrethroid resistance in southern African An. funestus is primarily conferred by monooxygenase enzyme metabolism. However, selection for this resistance mechanism is likely to have occurred in conjunction with other factors that improve production of the resistance phenotype. A strong candidate is cuticle thickening. This is because thicker cuticles lead to slower rates of insecticide absorption, which is likely to increase the efficiency of metabolic detoxification.

Results

Measures of mean cuticle thickness in laboratory samples of female An. funestus were obtained using scanning electron microscopy (SEM). These females were drawn from a laboratory colony carrying the pyrethroid resistance phenotype at a stable rate, but not fixed. Prior to cuticle thickness measurements, these samples were characterised as either more or less tolerant to permethrin exposure in one experiment, and either permethrin resistant or susceptible in another experiment. There was a significant and positive correlation between mean cuticle thickness and time to knock down during exposure to permethrin. Mean cuticle thickness was significantly greater in those samples characterised either as more tolerant or resistant to permethrin exposure compared to those characterised as either less tolerant or permethrin susceptible. Further, insecticide susceptible female An. funestus have thicker cuticles than their male counterparts.

Conclusion

Pyrethroid tolerant or resistant An. funestus females are likely to have thicker cuticles than less tolerant or susceptible females, and females generally have thicker cuticles than males. In pyrethroid resistant An. funestus, this increase in cuticle thickness is likely to have developed as an auxiliary to the primary mode of pyrethroid resistance which is based on enzyme-mediated detoxification.

Pyrethroid resistance in an Anopheles funestus population from Uganda

BACKGROUND

The susceptibility status of Anopheles funestus to insecticides remains largely unknown in most parts of Africa because of the difficulty in rearing field-caught mosquitoes of this malaria vector. Here we report the susceptibility status of the An. funestus population from Tororo district in Uganda and a preliminary characterisation of the putative resistance mechanisms involved.

METHODOLOGY/PRINCIPAL FINDINGS

A new forced egg laying technique used in this study significantly increased the numbers of field-caught females laying eggs and generated more than 4000 F1 adults. WHO bioassays indicated that An. funestus in Tororo is resistant to pyrethroids (62% mortality after 1 h exposure to 0.75% permethrin and 28% mortality to 0.05% deltamethrin). Suspected DDT resistance was also observed with 82% mortality. However this population is fully susceptible to bendiocarb (carbamate), malathion (organophosphate) and dieldrin with 100% mortality observed after exposure to each of these insecticides. Sequencing of a fragment of the sodium channel gene containing the 1014 codon conferring pyrethroid/DDT resistance in An. gambiae did not detect the L1014F kdr mutation but a correlation between haplotypes and resistance phenotype was observed indicating that mutations in other exons may be conferring the knockdown resistance in this species. Biochemical assays suggest that resistance in this population is mediated by metabolic resistance with elevated level of GSTs, P450s and pNPA compared to a susceptible strain of Anopheles gambiae. RT-PCR further confirmed the involvement of P450s with a 12-fold over-expression of CYP6P9b in the Tororo population compared to the fully susceptible laboratory colony FANG.

CONCLUSION

This study represents the first report of pyrethroid/DDT resistance in An. funestus from East Africa. With resistance already reported in southern and West Africa, this indicates that resistance in An. funestus may be more widespread than previously assumed and therefore this should be taken into account for the implementation and management of vector control programs in Africa.

High level of pyrethroid resistance in an Anopheles funestus population of the Chokwe District in Mozambique

BACKGROUND

Although Anopheles funestus is difficult to rear, it is crucial to analyse field populations of this malaria vector in order to successfully characterise mechanisms of insecticide resistance observed in this species in Africa. In this study we carried out a large-scale field collection and rearing of An. funestus from Mozambique in order to analyse its susceptibility status to insecticides and to broadly characterise the main resistance mechanisms involved in natural populations.

METHODOLOGY/PRINCIPAL FINDINGS

3,000 F(1) adults were obtained through larval rearing. WHO susceptibility assays indicated a very high resistance to pyrethroids with no mortality recorded after 1 h 30 min exposure and less than 50% mortality at 3 h 30 min. Resistance to the carbamate, bendiocarb was also noted, with 70% mortality after 1h exposure. In contrast, no DDT resistance was observed, indicating that no kdr-type resistance was involved. The sequencing of the acetylcholinesterase gene indicated the absence of the G119S and F455W mutations associated with carbamate and organophosphate resistance. This could explain the absence of malathion resistance in this population. Both biochemical assays and quantitative PCR implicated up-regulated P450 genes in pyrethroid resistance, with GSTs playing a secondary role. The carbamate resistance observed in this population is probably conferred by the observed altered AChE with esterases also involved.

CONCLUSION/SIGNIFICANCE

The high level of pyrethroid resistance in this population despite the cessation of pyrethroid use for IRS in 1999 is a serious concern for resistance management strategies such as rotational use of insecticides. As DDT has now been re-introduced for IRS, susceptibility to DDT needs to be closely monitored to prevent the appearance and spread of resistance to this insecticide.

Genome landscape and evolutionary plasticity of chromosomes in malaria mosquitoes

BACKGROUND

Nonrandom distribution of rearrangements is a common feature of eukaryotic chromosomes that is not well understood in terms of genome organization and evolution. In the major African malaria vector Anopheles gambiae, polymorphic inversions are highly nonuniformly distributed among five chromosomal arms and are associated with epidemiologically important adaptations. However, it is not clear whether the genomic content of the chromosomal arms is associated with inversion polymorphism and fixation rates.

METHODOLOGY/PRINCIPAL FINDINGS

To better understand the evolutionary dynamics of chromosomal inversions, we created a physical map for an Asian malaria mosquito, Anopheles stephensi, and compared it with the genome of An. gambiae. We also developed and deployed novel Bayesian statistical models to analyze genome landscapes in individual chromosomal arms An. gambiae. Here, we demonstrate that, despite the paucity of inversion polymorphisms on the X chromosome, this chromosome has the fastest rate of inversion fixation and the highest density of transposable elements, simple DNA repeats, and GC content. The highly polymorphic and rapidly evolving autosomal 2R arm had overrepresentation of genes involved in cellular response to stress supporting the role of natural selection in maintaining adaptive polymorphic inversions. In addition, the 2R arm had the highest density of regions involved in segmental duplications that clustered in the breakpoint-rich zone of the arm. In contrast, the slower evolving 2L, 3R, and 3L, arms were enriched with matrix-attachment regions that potentially contribute to chromosome stability in the cell nucleus.

CONCLUSIONS/SIGNIFICANCE

These results highlight fundamental differences in evolutionary dynamics of the sex chromosome and autosomes and revealed the strong association between characteristics of the genome landscape and rates of chromosomal evolution. We conclude that a unique combination of various classes of genes and repetitive DNA in each arm, rather than a single type of repetitive element, is likely responsible for arm-specific rates of rearrangements.

Overexpression of cytochrome P450 genes in pyrethroid-resistant Culex quinquefasciatus

JPal-per strain of Culex quinquefasciatus exhibits extremely high resistance against pyrethroids in larvae, though the resistance is greatly lower in adults. Increased microsome monooxygenase metabolism is one of the major factors of the larval resistance in this strain. We cloned 46 novel cytochrome P450 cDNAs from JPal-per strain. An oligonucleotide microarray was designed for the novel 46 genes plus 16 previously reported P450 genes along with other non-P450 gene probes. Of these, five P450 genes were upregulated (>2.5-fold) in JPal-per larvae as compared with a susceptible strain. The expression ratios for the highest three among the five P450 genes screened in the microarray analysis, CYP9M10, CYP4H34 and CYP6Z10, were further validated by qPCR as 264-, 8.3-, and 3.9-fold, respectively. In JPal-per, the transcription levels of CYP9M10 and CYP4H34 showed a similar stage-dependent pattern as a high expression level during the larvfrom Ogasawara Islands in Japanal stage dramatically decreases in the adult stage. This larval specific overexpression manner of the two genes was consistent with the characteristic of stage-dependent resistance of JPal-per strain previously reported, suggesting that the two P450s, CYP9M10 and CYP4H34, are involved in pyrethroid detoxification in JPal-per strain.

An initial linkage map of the West Nile Virus vector Culex tarsalis

We have constructed the first genetic linkage map for the North American arboviral vector Culex tarsalis. 120 F(2) offspring from a cross between two colonies were genotyped using 25 microsatellites and six inter-simple sequence repeat (ISSR) markers. We resolved four linkage groups which likely correspond to two full-length chromosomes and two arms of the final chromosome. The longest linkage group contains the sex locus and corresponds to chromosome 3. Recombination rates around the sex locus were dramatically higher in females compared to males. The majority of microsatellite loci share sequence identity with regions of the Culex quinquefasciatus genome, whose assembly should aid in anchoring linkage groups to physical chromosomes. This map will aid in identification of loci involved with variable phenotypes in C. tarsalis including WNV susceptibility.

Isolation and characterization of microsatellite markers from the Mediterranean fruit fly, Ceratitis capitata: cross-species amplification in other Tephritidae species reveals a varying degree of transferability

The Mediterranean fruit fly, Ceratitis capitata, is a pest of major economic importance and has become a model for the development of SIT control programs for insect pests. Significant information has been accumulated on classical and population genetics of this species during the past 2 decades. However, the availability of molecular markers is limited. Here, we present the isolation and characterization of 159 microsatellite clones and the development of 108 polymorphic microsatellite markers for this insect pest. Mapping by in situ hybridization to polytene chromosomes of 21 microsatellite clones enriched the cytogenetic map that was previously constructed by our group. The enriched map provides a large number of STSs for future genome mapping projects. Cross-species amplification of these microsatellite loci in 12 Tephritidae species and sequence analysis of several amplification products indicated a varying degree of transferability and their possible usefulness as molecular and genetic markers in these species where genetic and molecular tools are limited.

The cytochrome P450 (CYP) gene superfamily in Daphnia pulex

BACKGROUND

Cytochrome P450s (CYPs) in animals fall into two categories: those that synthesize or metabolize endogenous molecules and those that interact with exogenous chemicals from the diet or the environment. The latter form a critical component of detoxification systems.

RESULTS

Data mining and manual curation of the Daphnia pulex genome identified 75 functional CYP genes, and three CYP pseudogenes. These CYPs belong to 4 clans, 13 families, and 19 subfamilies. The CYP 2, 3, 4, and mitochondrial clans are the same four clans found in other sequenced protostome genomes. Comparison of the CYPs from D. pulex to the CYPs from insects, vertebrates and sea anemone (Nematostella vectensis) show that the CYP2 clan, and to a lesser degree, the CYP4 clan has expanded in Daphnia pulex, whereas the CYP3 clan has expanded in insects. However, the expansion of the Daphnia CYP2 clan is not as great as the expansion observed in deuterostomes and the nematode C. elegans. Mapping of CYP tandem repeat regions demonstrated the unusual expansion of the CYP370 family of the CYP2 clan. The CYP370s are similar to the CYP15s and CYP303s that occur as solo genes in insects, but the CYP370s constitute approximately 20% of all the CYP genes in Daphnia pulex. Lastly, our phylogenetic comparisons provide new insights into the potential origins of otherwise mysterious CYPs such as CYP46 and CYP19 (aromatase).

CONCLUSION

Overall, the cladoceran, D. pulex has a wide range of CYPs with the same clans as insects and nematodes, but with distinct changes in the size and composition of each clan.

Two duplicated P450 genes are associated with pyrethroid resistance in Anopheles funestus, a major malaria vector

Pyrethroid resistance in Anopheles funestus is a potential obstacle to malaria control in Africa. Tools are needed to detect resistance in field populations. We have been using a positional cloning approach to identify the major genes conferring pyrethroid resistance in this vector. A quantitative trait locus (QTL) named rp1 explains 87% of the genetic variance in pyrethroid susceptibility in two families from reciprocal crosses between susceptible and resistant strains. Two additional QTLs of minor effect, rp2 and rp3, were also detected. We sequenced a 120-kb BAC clone spanning the rp1 QTL and identified 14 protein-coding genes and one putative pseudogene. Ten of the 14 genes encoded cytochrome P450s, and expression analysis indicated that four of these P450s were differentially expressed between susceptible and resistant strains. Furthermore, two of these genes, CYP6P9 and CYP6P4, which are 25 and 51 times overexpressed in resistant females, are tandemly duplicated in the BAC clone as well as in laboratory and field samples, suggesting that P450 gene duplication could contribute to pyrethroid resistance in An. funestus. Single nucleotide polymorphisms (SNPs) were identified within CYP6P9 and CYP6P4, and genotyping of the progeny of the genetic crosses revealed a maximum penetrance value f(2) = 1, confirming that these SNPs are valid resistance markers in the laboratory strains. This serves as proof of principle that a DNA-based diagnostic test could be designed to trace metabolic resistance in field populations. This will be a major advance for insecticide resistance management in malaria vectors, which requires the early detection of resistance alleles.

Field-caught permethrin-resistant Anopheles gambiae overexpress CYP6P3, a P450 that metabolises pyrethroids

Insects exposed to pesticides undergo strong natural selection and have developed various adaptive mechanisms to survive. Resistance to pyrethroid insecticides in the malaria vector Anopheles gambiae is receiving increasing attention because it threatens the sustainability of malaria vector control programs in sub-Saharan Africa. An understanding of the molecular mechanisms conferring pyrethroid resistance gives insight into the processes of evolution of adaptive traits and facilitates the development of simple monitoring tools and novel strategies to restore the efficacy of insecticides. For this purpose, it is essential to understand which mechanisms are important in wild mosquitoes. Here, our aim was to identify enzymes that may be important in metabolic resistance to pyrethroids by measuring gene expression for over 250 genes potentially involved in metabolic resistance in phenotyped individuals from a highly resistant, wild A. gambiae population from Ghana. A cytochrome P450, CYP6P3, was significantly overexpressed in the survivors, and we show that the translated enzyme metabolises both alpha-cyano and non-alpha-cyano pyrethroids. This is the first study to demonstrate the capacity of a P450 identified in wild A. gambiae to metabolise insecticides. The findings add to the understanding of the genetic basis of insecticide resistance in wild mosquito populations.

Expression of the cytochrome P450s, CYP6P3 and CYP6M2 are significantly elevated in multiple pyrethroid resistant populations of Anopheles gambiae s.s. from Southern Benin and Nigeria

Background

Insecticide resistance in Anopheles mosquitoes is threatening the success of malaria control programmes. This is particularly true in Benin where pyrethroid resistance has been linked to the failure of insecticide treated bed nets. The role of mutations in the insecticide target sites in conferring resistance has been clearly established. In this study, the contribution of other potential resistance mechanisms was investigated in Anopheles gambiae s.s. from a number of localities in Southern Benin and Nigeria. The mosquitoes were sampled from a variety of breeding sites in a preliminary attempt to investigate the role of contamination of mosquito breeding sites in selecting for resistance in adult mosquitoes.

Results

All mosquitoes sampled belonged to the M form of An. gambiae s.s. There were high levels of permethrin resistance in an agricultural area (Akron) and an urban area (Gbedjromede), low levels of resistance in mosquito samples from an oil contaminated site (Ojoo) and complete susceptibility in the rural Orogun location. The target site mutation kdrW was detected at high levels in two of the populations (Akron f = 0.86 and Gbedjromede f = 0.84) but was not detected in Ojoo or Orogun. Microarray analysis using the Anopheles gambiae detox chip identified two P450s, CYP6P3 and CYP6M2 up regulated in all three populations, the former was expressed at particularly high levels in the Akron (12.4-fold) and Ojoo (7.4-fold) populations compared to the susceptible population. Additional detoxification and redox genes were also over expressed in one or more populations including two cuticular pre-cursor genes which were elevated in two of the three resistant populations.

Conclusion

Multiple resistance mechanisms incurred in the different breeding sites contribute to resistance to permethrin in Benin. The cytochrome P450 genes, CYP6P3 and CYP6M2 are upregulated in all three resistant populations analysed. Several additional potential resistance mechanisms were also identified that warrant further investigation. Metabolic genes were over expressed irrespective of the presence of kdr, the latter resistance mechanism being absent in one resistant population. The discovery that mosquitoes collected from different types of breeding sites display differing profiles of metabolic genes at the adult stage may reflect the influence of a range of xenobiotics on selecting for resistance in mosquitoes.

The effect of a single blood meal on the phenotypic expression of insecticide resistance in the major malaria vector Anopheles funestus

BACKGROUND

Anopheles funestus is a major malaria vector in southern Africa. Vector control relies on the use of insecticide chemicals to significantly reduce the number of malaria vectors by targeting that portion of the female population that takes blood meals and subsequently rests indoors. It has been suggested that the intake of a blood meal may assist female mosquitoes to tolerate higher doses of insecticide through vigour tolerance. It is hypothesized that during the process of blood digestion, detoxification mechanisms required for the neutralizing of harmful components in the blood meal may also confer an increased ability to tolerate insecticide intoxication through increased enzyme regulation.

METHODS

Bottle bioassays using a range of concentrations of the pyrethroid insecticide permethrin were performed on pyrethroid susceptible and resistant laboratory strains of An. funestus in order to detect differences in insecticide susceptibility following a single blood meal. Based on these results, a discriminating dosage was identified (double the lowest dosage that resulted in 100% mortality of the susceptible strain). Blood-fed and unfed females drawn from the resistant strain of An. funestus were then assayed against this discriminating dose, and the percentage mortality for each sample was scored and compared.

RESULTS

In the insecticide dose response assays neither the fully susceptible nor the resistant strain of An. funestus showed any significant difference in insecticide susceptibility following a blood meal, regardless of the stage of blood meal digestion. A significant increase in the level of resistance was however detected in the resistant An. funestus strain following a single blood meal, based on exposure to a discriminating dose of permethrin.

CONCLUSION

The fully susceptible An. funestus strain did not show any significant alteration in susceptibility to insecticide following a blood meal suggesting that vigour tolerance through increased body mass (and increased dilution of internalized insecticide) does not play a significant role in tolerance to insecticide intoxication. The increase in insecticide tolerance in the pyrethroid resistant strain of An. funestus following a blood meal suggests that insecticide detoxification mechanisms involved in insecticide resistance are stimulated by the presence of a blood meal prior to insecticide exposure, leading to enhanced expression of the resistance phenotype. This finding may be significant in terms of the methods used to control indoor resting populations of An. funestus if the mass killing effect of insecticide application proves increasingly inadequate against blood-feeding females already carrying the insecticide resistance phenotype.

Quantitative trait loci mapping of genome regions controlling permethrin resistance in the mosquito Aedes aegypti

The mosquito Aedes aegypti is the principal vector of dengue and yellow fever flaviviruses. Permethrin is an insecticide used to suppress Ae. aegypti adult populations but metabolic and target site resistance to pyrethroids has evolved in many locations worldwide. Quantitative trait loci (QTL) controlling permethrin survival in Ae. aegypti were mapped in an F(3) advanced intercross line. Parents came from a collection of mosquitoes from Isla Mujeres, México, that had been selected for permethrin resistance for two generations and a reference permethrin-susceptible strain originally from New Orleans. Following a 1-hr permethrin exposure, 439 F(3) adult mosquitoes were phenotyped as knockdown resistant, knocked down/recovered, or dead. For QTL mapping, single nucleotide polymorphisms (SNPs) were identified at 22 loci with potential antixenobiotic activity including genes encoding cytochrome P450s (CYP), esterases (EST), or glutathione transferases (GST) and at 12 previously mapped loci. Seven antixenobiotic genes mapped to chromosome I, six to chromosome II, and nine to chromosome III. Two QTL of major effect were detected on chromosome III. One corresponds with a SNP previously associated with permethrin resistance in the para sodium channel gene and the second with the CCEunk7o esterase marker. Additional QTL but of relatively minor effect were also found. These included two sex-linked QTL on chromosome I affecting knockdown and recovery and a QTL affecting survival and recovery. On chromosome II, one QTL affecting survival and a second affecting recovery were detected. The patterns confirm that mutations in the para gene cause target-site insensitivity and are the major source of permethrin resistance but that other genes dispersed throughout the genome contribute to recovery and survival of mosquitoes following permethrin exposure.