Genetic mapping of genes conferring permethrin resistance in the malaria vector, Anopheles gambiae

Whole transcriptomic analysis reveals overexpression of salivary gland and cuticular proteins genes in insecticide-resistant Anopheles arabiensis from Western Kenya

BACKGROUND

Effective vector control is key to malaria prevention. However, this is now compromised by increased insecticide resistance due to continued reliance on insecticide-based control interventions. In Kenya, we have observed heterogenous resistance to pyrethroids and organophosphates in Anopheles arabiensis which is one of the most widespread malaria vectors in the country. We investigated the gene expression profiles of insecticide resistant An. arabiensis populations from Migori and Siaya counties in Western Kenya using RNA-Sequencing. Centers for Disease Control and Prevention (CDC) bottle assays were conducted using deltamethrin (DELTA), alphacypermethrin (ACYP) and pirimiphos-methyl (PMM) to determine the resistance status in both sites.

RESULTS

Mosquitoes from Migori had average mortalities of 91%, 92% and 58% while those from Siaya had 85%, 86%, and 30% when exposed to DELTA, ACYP and PMM, respectively. RNA-Seq analysis was done on pools of mosquitoes which survived exposure ('resistant'), mosquitoes that were not exposed, and the insecticide-susceptible An. arabiensis Dongola strain. Gene expression profiles of resistant mosquitoes from both Migori and Siaya showed an overexpression mainly of salivary gland proteins belonging to both the short and long form D7 genes, and cuticular proteins (including CPR9, CPR10, CPR15, CPR16). Additionally, the overexpression of detoxification genes including cytochrome P450s (CYP9M1, CYP325H1, CYP4C27, CYP9L1 and CYP307A1), 2 carboxylesterases and a glutathione-S-transferase (GSTE4) were also shared between DELTA, ACYP, and PMM survivors, pointing to potential contribution to cross resistance to both pyrethroid and organophosphate insecticides.

CONCLUSION

This study provides novel insights into the molecular basis of insecticide resistance in An. arabiensis in Western Kenya and suggests that salivary gland proteins and cuticular proteins are associated with resistance to multiple classes of insecticides.

Predicting the genomic resolution of bulk segregant analysis

Bulk segregant analysis is a technique for identifying the genetic loci that underlie phenotypic trait differences. The basic approach is to compare two pools of individuals from the opposing tails of the phenotypic distribution, sampled from an interbred population. Each pool is sequenced and scanned for alleles that show divergent frequencies between the pools, indicating potential association with the observed trait differences. Bulk segregant analysis has already been successfully applied to the mapping of various quantitative trait loci in organisms ranging from yeast to maize. However, these studies have typically suffered from rather low mapping resolution, and we still lack a detailed understanding of how this resolution is affected by experimental parameters. Here, we use coalescence theory to calculate the expected genomic resolution of bulk segregant analysis for a simple monogenic trait. We first show that in an idealized interbreeding population of infinite size, the expected length of the mapped region is inversely proportional to the recombination rate, the number of generations of interbreeding, and the number of genomes sampled, as intuitively expected. In a finite population, coalescence events in the genealogy of the sample reduce the number of potentially informative recombination events during interbreeding, thereby increasing the length of the mapped region. This is incorporated into our model by an effective population size parameter that specifies the pairwise coalescence rate of the interbreeding population. The mapping resolution predicted by our calculations closely matches numerical simulations and is surprisingly robust to moderate levels of contamination of the segregant pools with alternative alleles. Furthermore, we show that the approach can easily be extended to modifications of the crossing scheme. Our framework will allow researchers to predict the expected power of their mapping experiments, and to evaluate how their experimental design could be tuned to optimize mapping resolution.

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.

Characterization of permethrin resistance in a Musca domestica strain: resistance development, cross-resistance potential and realized heritability

BACKGROUND

Permethrin (a pyrethroid) has been recommended for the management of a number of insect pests, including Musca domestica. Recently, resistance to permethrin has been reported from the Punjab province, Pakistan, which necessitates the need to manage the resistance problem. For this reason, a field strain of M. domestica was selected with permethrin for ten consecutive generations (G1-G10) to determine the risk assessment, the rapidity with which permethrin resistance can develop, and the presence or absence of the cross-resistance phenomenon.

RESULTS

The results revealed that permethrin resistance increased very quickly in response to selection experiments. The resistance ratio (RR) value increased by 16.04- to 605.92-fold and the LD₅₀ value increased by 100.22 to 3787.08 ng fly^-1 after selection experiments. Moreover, the increase in LD₅₀ value remained stable when the selected strain (Perm-SEL) was cultured in the absence of permethrin for the next five generations (G11-G15). The results revealed the possibility of minor cross-resistance in the Perm-SEL strain with β-cyfluthrin and deltamethrin, and lack of cross-resistance with imidacloprid and spinosad. The realized heritability (h² ) value for the Perm-SEL strain was 0.23. Assuming a mean slope value of 1.27 and an h² value of 0.23, then ten to four generations would be estimated for a 10-fold increase in the LD₅₀ values at 50-90% selection intensity, respectively.

CONCLUSION

The results of present study revealed high risk for permethrin resistance development under continuous selection pressure. Permethrin resistance developed very quickly under laboratory conditions, which might reflect the likelihood of resistance evolution in field conditions if permethrin selection pressure holds for a long time. © 2019 Society of Chemical Industry.

Long-Term Population Studies Uncover the Genome Structure and Genetic Basis of Xenobiotic and Host Plant Adaptation in the Herbivore Tetranychus urticae

Pesticide resistance arises rapidly in arthropod herbivores, as can host plant adaptation, and both are significant problems in agriculture. These traits have been challenging to study as both are often polygenic and many arthropods are genetically intractable. Here, we examined the genetic architecture of pesticide resistance and host plant adaptation in the two-spotted spider mite, Tetranychus urticae, a global agricultural pest. We show that the short generation time and high fecundity of T. urticae can be readily exploited in experimental evolution designs for high-resolution mapping of quantitative traits. As revealed by selection with spirodiclofen, an acetyl-CoA carboxylase inhibitor, in populations from a cross between a spirodiclofen-resistant and a spirodiclofen-susceptible strain, and which also differed in performance on tomato, we found that a limited number of loci could explain quantitative resistance to this compound. These were resolved to narrow genomic intervals, suggesting specific candidate genes, including acetyl-CoA carboxylase itself, clustered and copy variable cytochrome P450 genes, and NADPH cytochrome P450 reductase, which encodes a redox partner for cytochrome P450s. For performance on tomato, candidate genomic regions for response to selection were distinct from those responding to the synthetic compound and were consistent with a more polygenic architecture. In accomplishing this work, we exploited the continuous nature of allele frequency changes across experimental populations to resolve the existing fragmented T. urticae draft genome to pseudochromosomes. This improved assembly was indispensable for our analyses, as it will be for future research with this model herbivore that is exceptionally amenable to genetic studies.

Identification of glutathione S-transferases in Bemisia tabaci (Hemiptera: Aleyrodidae) and evidence that GSTd7 helps explain the difference in insecticide susceptibility between B. tabaci Middle East-Minor Asia 1 and Mediterranean

The Bemisia tabaci (Gennadius) (Hemiptera:Aleyrodidae) species complex includes invasive and destructive pests of field crops, and the sibling species MEAM1 and MED are its two most damaging members. Previous research indicated that the replacement of Middle East-Minor Asia 1 (MEAM1) by Mediterranean (MED) as the dominant B. tabaci species in China can be mainly attributed to MED's greater tolerance to insecticides. Glutathione S-transferases (GSTs) play important roles in the detoxification of hydrophobic toxic compounds. To increase our understanding of differences in insecticide resistance between B. tabaci MEAM1 and MED, we searched the genomic and transcriptomic databases and identified 23 putative GSTs in both B. tabaci MEAM1 and MED. Through measuring mRNA levels of 18 of the GSTs after B. tabaci MEAM1 and MED adults were exposed to the insecticide imidacloprid, we found that the expression levels were increased more in B. tabaci MED than in MEAM1 (in particular, the expression level of GST-d7 was increased by 4.39-fold relative to the control). Knockdown of GST-d7 in B. tabaci MED but not in B. tabaci MEAM1 resulted in a substantial increase in the mortality of imidacloprid-treated adults. These results indicate that differences in GST-d7 may help explain why insecticide tolerance is greater in B. tabaci MED than in B. tabaci MEAM1.

Human Interventions: Driving Forces of Mosquito Evolution

One of the most common strategies for controlling mosquito-borne diseases relies on the use of chemical pesticides to repel or kill the mosquito vector. Pesticide exposure interferes with several key biological functions in the mosquito and triggers a variety of adaptive responses whose underlying mechanisms are only partially elucidated. The availability of reference genome sequences opens up the possibility of tracking signatures of evolutionary changes, including the most recent, across the genomes of many vector species. In this review, we highlight the recent genomic changes, which contribute to the fascinating adaptation of malaria vectors of the sub-Saharan African region to intensive insecticide-based interventions. We emphasize the operational significance of detailed genomic knowledge for current monitoring and decision-making.

Resistance Status to the Insecticides Temephos, Deltamethrin, and Diflubenzuron in Brazilian Aedes aegypti Populations

Insecticides are still largely applied in public health to control disease vectors. In Brazil, organophosphates (OP) and pyrethroids (PY) are used against Aedes aegypti for years. Since 2009 Insect Growth Regulators (IGR) are also employed in the control of larvae. We quantified resistance to temephos (OP), deltamethrin (PY), and diflubenzuron (IGR) of A. aegypti samples from 12 municipalities distributed throughout the country, collected between 2010 and 2012. High levels of resistance to neurotoxic insecticides were detected in almost all populations: RR₉₅ to temephos varied between 4.0 and 27.1; the lowest RR₉₅ to deltamethrin was 13.1, and values higher than 70.0 were found. In contrast, all samples were susceptible to diflubenzuron (RR₉₅ < 2.3). Biochemical tests performed with larvae and adults discarded the participation of acetylcholinesterase, the OP target, and confirmed involvement of the detoxifying enzymes esterases, mixed function oxidases, and glutathione-S-transferases. The results obtained were discussed taking into account the public chemical control component and the increase in the domestic use of insecticides during dengue epidemic seasons in the evaluated municipalities.

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.

Comparative transcriptome analysis of lufenuron-resistant and susceptible strains of Spodoptera frugiperda (Lepidoptera: Noctuidae)

BACKGROUND

The evolution of insecticide resistance in Spodoptera frugiperda (Lepidoptera: Noctuidae) has resulted in large economic losses and disturbances to the environment and agroecosystems. Resistance to lufenuron, a chitin biosynthesis inhibitor insecticide, was recently documented in Brazilian populations of S. frugiperda. Thus, we utilized large-scale cDNA sequencing (RNA-Seq analysis) to compare the pattern of gene expression between lufenuron-resistant (LUF-R) and susceptible (LUF-S) S. larvae in an attempt to identify the molecular basis behind the resistance mechanism(s) of S. frugiperda to this insecticide.

RESULTS

A transcriptome was assembled using approximately 19.6 million 100 bp-long single-end reads, which generated 18,506 transcripts with a N50 of 996 bp. A search against the NCBI non-redundant database generated 51.1% (9,457) functionally annotated transcripts. A large portion of the alignments were homologous to insects, with the majority (45%) being similar to sequences of Bombyx mori (Lepidoptera: Bombycidae). Moreover, 10% of the alignments were similar to sequences of various species of Spodoptera (Lepidoptera: Noctuidae), with 3% of them being similar to sequences of S. frugiperda. A comparative analysis of the gene expression between LUF-R and LUF-S S. frugiperda larvae identified 940 differentially expressed transcripts (p ≤ 0.05, t-test; fold change ≥ 4). Six of them were associated with cuticle metabolism. Of those, four were overexpressed in LUF-R larvae. The machinery involved with the detoxification process was represented by 35 differentially expressed transcripts; 24 of them belonging to P450 monooxygenases, four to glutathione-S-transferases, six to carboxylases and one to sulfotransferases. RNA-Seq analysis was validated for a number of selected candidate transcripts by using quantitative real time PCR (qPCR).

CONCLUSIONS

The gene expression profile of LUF-R larvae of S. frugiperda differs from LUF-S larvae. In general, gene expression is much higher in resistant larvae when compared to the susceptible ones, particularly for those genes involved with pathways for xenobiotic detoxification, mainly represented by P450 monooxygenases transcripts. Our data indicate that enzymes involved with the detoxification process, and mostly the P450, are one of the resistance mechanisms employed by the LUF-R S. frugiperda larvae against lufenuron.

Identification of a major Quantitative Trait Locus determining resistance to the organophosphate temephos in the dengue vector mosquito Aedes aegypti

Organophosphate insecticides (OP) have extensively been used to control mosquitoes, such as the vector Aedes aegypti. Unfortunately, OP resistance has hampered control programs worldwide. We used Quantitative Trait Locus (QTL) mapping to evaluate temephos resistance in two F1 intercross populations derived from crosses between a resistant Ae. aegypti strain (RecR) and two susceptible strains (MoyoD and Red). A single major effect QTL was identified on chromosome 2 of both segregating populations, named rtt1 (resistance to temephos 1). Bioinformatics analyses identified a cluster of carboxylesterase genes (CCE) within the rtt1 interval. qRT-PCR demonstrated that different CCEs were up-regulated in F2 resistant individuals from both crosses. However, none exceeded the 2-fold expression. Primary mechanisms for temephos resistance may vary between Ae. aegypti populations, yet also appear to support previous findings suggesting that multiple linked esterase genes may contribute to temephos resistance in the RecR strain as well as other populations.

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.

RNA-seq analyses of changes in the Anopheles gambiae transcriptome associated with resistance to pyrethroids in Kenya: identification of candidate-resistance genes and candidate-resistance SNPs

Background

The extensive use of pyrethroids for control of malaria vectors, driven by their cost, efficacy and safety, has led to widespread resistance. To favor their sustainable use, the World Health Organization (WHO) formulated an insecticide resistance management plan, which includes the identification of the mechanisms of resistance and resistance surveillance. Recognized physiological mechanisms of resistance include target site mutations in the para voltage-gated sodium channel, metabolic detoxification and penetration resistance. Such understanding of resistance mechanisms has allowed the development of resistance monitoring tools, including genotyping of the kdr mutation L1014F/S in the para gene.

Methods

The sequence-based technique RNA-seq was applied to study changes in the transcriptome of deltamethrin-resistant and -susceptible Anopheles gambiae mosquitoes from the Western Province of Kenya. The resulting gene expression profiles were compared to data in the most recent literature to derive a list of candidate resistance genes. RNA-seq data were analyzed also to identify sequence polymorphisms linked to resistance.

Results

A total of five candidate-resistance genes (AGAP04177, AGAP004572, AGAP008840, AGAP007530 and AGAP013036) were identified with altered expression between resistant and susceptible mosquitoes from West and East Africa. A change from G to C at position 36043997 of chromosome 3R resulting in A101G of the sulfotransferase gene AGAP009551 was significantly associated with the resistance phenotype (odds ratio: 5.10). The kdr L1014S mutation was detected at similar frequencies in both phenotypically resistant and susceptible mosquitoes, suggesting it is no longer fully predictive of the resistant phenotype.

Conclusions

Overall, these results support the conclusion that resistance to pyrethroids is a complex and evolving phenotype, dependent on multiple gene functions including, but not limited to, metabolic detoxification. Functional convergence among metabolic detoxification genes may exist, with the role of each gene being modulated by the life history and selection pressure on mosquito populations. As a consequence, biochemical assays that quantify overall enzyme activity may be a more suitable method for predicting metabolic resistance than gene-based assays.

Electronic supplementary material

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

Identification and characterisation of multiple glutathione S-transferase genes from the diamondback moth, Plutella xylostella

BACKGROUND

The diamondback moth (DBM), Plutella xylostella, is one of the most harmful insect pests on crucifer crops worldwide. In this study, 19 cDNAs encoding glutathione S-transferases (GSTs) were identified from the genomic and transcriptomic database for DBM (KONAGAbase) and further characterized.

RESULTS

Phylogenetic analysis showed that the 19 GSTs were classified into six different cytosolic classes, including four in delta, six in epsilon, three in omega, two in sigma, one in theta and one in zeta. Two GSTs were unclassified. RT-PCR analysis revealed that most GST genes were expressed in all developmental stages, with higher expression in the larval stages. Six DBM GSTs were expressed at the highest levels in the midgut tissue. Twelve purified recombinant GSTs showed varied enzymatic properties towards 1-chloro-2,4-dinitrobenzene and glutathione, whereas rPxGSTo2, rPxGSTz1 and rPxGSTu2 had no activity. Real-time quantitative PCR revealed that expression levels of the 19 DBM GST genes were varied and changed after exposure to acephate, indoxacarb, beta-cypermethrin and spinosad. PxGSTd3 was significantly overexpressed, while PxGSTe3 and PxGSTs2 were significantly downregulated by all four insecticide exposures.

CONCLUSION

The changes in DBM GST gene expression levels exposed to different insecticides indicate that they may play individual roles in tolerance to insecticides and xenobiotics.

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.

Phylogenetic analysis of the GST family in Anopheles (Nyssorhynchus) darlingi

Anopheles darlingi Root, 1926 and Anopheles gambiae (Diptera: Culicidae) are the most important human malaria vectors in South America and Africa, respectively. The two species are estimated to have diverged 100 million years ago. Studies on the phylogenetics and evolution of gene sequences, such as glutathione S-transferase (GST) in disease-transmitting mosquitoes are scarce. The sigma class GST (KC890767) from the transcriptome of An. darlingi captured in the Brazilian Amazon was studied by in silico hybridization, and mapped to chromosome 3 of An. gambiae. The sigma class GST of An. darlingi was used for phylogenetic analyses to understand the GST base composition of the most recent common ancestor between An. darlingi, Anopheles gambiae, Aedes aegypti and Culex quinquefasciatus. The GST (KC890767) of An. darlingi was studied to generate the main divergence branches using a Neighbor-Joining and bootstrapping approaches to confirm confidence levels on the tree nodes that separate the An. darlingi and other mosquito species. The results showed divergence between An. gambiae, Ae. Aegypti, Cx. quinquefasciatus, and Phlebotomus papatasi as outgroup, and the homology relationship between sigma class GST of An. darlingi and GSTS1_1 gene of An. gambiae was valuable for phylogenetic and evolutionary studies.

The role of vector control in stopping the transmission of malaria: threats and opportunities

Malaria control, and that of other insect borne diseases such as dengue, is heavily dependent on our ability to control the mosquito populations that transmit these diseases. The major push over the last decade to reduce the global burden of malaria has been driven by the distribution of pyrethroid insecticide-treated bednets and an increase in coverage of indoor residual spraying (IRS). This has reduced malaria deaths by a third. Progress towards the goal of reducing this further is threatened by lack of funding and the selection of drug and insecticide resistance. When malaria control was initially scaled up, there was little pyrethroid resistance in the major vectors, today there is no country in Africa where the vectors remain fully susceptible to pyrethroids. The first pyrethroid resistance mechanisms to be selected produced low-level resistance which had little or no operational significance. More recently, metabolically based resistance has been selected, primarily in West Africa, which in some mosquito populations produces more than 1000-fold resistance. As this spreads the effectiveness of pyrethroid-based bednets and IRS will be compromised. New public health insecticides are not readily available. The pipeline of agrochemical insecticides that can be re-purposed for public health dried up 30 years ago when the target product profile for agricultural insecticides shifted from broad spectrum, stable, contact-acting insecticides to narrow spectrum stomach poisons that could be delivered through the plant. A public–private partnership, the Innovative Vector Control Consortium, was established in 2005 to stimulate the development of new public health pesticides. Nine potential new classes of chemistry are in the pipeline, with the intention of developing three into new insecticides. While this has been successfully achieved, it will still take 6–9 years for new insecticides to reach the market. Careful management of the resistance situation in the interim will be needed if current gains in malaria control are not to be reversed.

Computational evidence for the detoxifying mechanism of epsilon class glutathione transferase toward the insecticide DDT

A combined quantum mechanics/molecular mechanics (QM/MM) computation of the detoxifying mechanism of an epsilon class glutathione transferases (GSTs) toward organochlorine insecticide DDT, 1,1,1-trichloro-2,2-bis(p-chlorophenyl)ethane, has been carried out. The exponential average barrier of the proton transfer mechanism is 15.2 kcal/mol, which is 27.6 kcal/mol lower than that of the GS-DDT conjugant mechanism. It suggests that the detoxifying reaction proceeds via a proton transfer mechanism where GSH acts as a cofactor rather than a conjugate. The study reveals that the protein environment has a strong effect on the reaction barrier. The experimentally proposed residues Arg112, Glu116 and Phe120 were found to have a strong influence on the detoxifying reaction. The influence of residues Pro13, Cys15, His53, Ile55, Glu67, Ser68, Phe115, and Leu119 was detected as well. It is worth noticing that Ile55 facilitates the detoxifying reaction most. On the basis of the structure of DDT, structure 2, (BrC6H4)2CHCCl3, is the best candidate among all the tested structures in resisting the detoxification of enzyme agGSTe2.

Co-occurrence and distribution of East (L1014S) and West (L1014F) African knock-down resistance in Anopheles gambiae sensu lato population of Tanzania

Objective

Insecticide resistance molecular markers can provide sensitive indicators of resistance development in Anopheles vector populations. Assaying these makers is of paramount importance in the resistance monitoring programme. We investigated the presence and distribution of knock-down resistance (kdr) mutations in Anopheles gambiae s.l. in Tanzania.

Methods

Indoor-resting Anopheles mosquitoes were collected from 10 sites and tested for insecticide resistance using the standard WHO protocol. Polymerase chain reaction-based molecular diagnostics were used to genotype mosquitoes and detect kdr mutations.

Results

The An. gambiae tested were resistance to lambdacyhalothrin in Muheza, Arumeru and Muleba. Out of 350 An. gambiae s.l. genotyped, 35% were An. gambiae s.s. and 65% An. arabiensis. L1014S and L1014F mutations were detected in both An. gambiae s.s. and An. arabiensis. L1014S point mutation was found at the allelic frequency of 4–33%, while L1014F was at the allelic frequency 6–41%. The L1014S mutation was much associated with An. gambiae s.s. (χ² = 23.41; P < 0.0001) and L1014F associated with An. arabiensis (χ² = 11.21; P = 0.0008). The occurrence of the L1014S allele was significantly associated with lambdacyhalothrin resistance mosquitoes (Fisher exact P < 0.001).

Conclusion

The observed co-occurrence of L1014S and L1014F mutations coupled with reports of insecticide resistance in the country suggest that pyrethroid resistance is becoming a widespread phenomenon among our malaria vector populations. The presence of L1014F mutation in this East African mosquito population indicates the spreading of this gene across Africa. The potential operational implications of these findings on malaria control need further exploration.

Objectif

Les marqueurs moléculaires de la résistance aux insecticides peuvent fournir des indicateurs sensibles du développement de la résistance dans les populations de vecteurs Anopheles. Le test de ces indicateurs est d'une importance énorme dans le programme de surveillance de la résistance. Nous avons étudié la présence et la répartition des mutations de résistance knockdown (kdr) chez Anopheles gambiae s.l. en Tanzanie.

Méthodes

Des anophèles d'intérieur, au repos ont été collectées dans 10 sites et testées pour la résistance aux insecticides en utilisant le protocole standard de l'OMS. Les diagnostics moléculaires basés sur la PCR ont été utilisés pour le génotypage des moustiques et la détection des génotypes kdr.

Résultats

Les An. gambiae testées étaient résistantes à la lambdacyhalothrine à Muheza, Arumeru et Muleba. Sur 350 An. gambiae s.l. génotypées, 35% étaient An. gambiae s.s. et 65% étaient An. arabiensis. Les mutations L1014S et L1014F ont été détectées à la fois chez An. gambiae s.s. et An. arabiensis. La mutation ponctuelle L1014S a été trouvée à la fréquence allélique de 4 à 33%, tandis que L1014F était à la fréquence allélique de 6 à 14%. La mutation L1014S a été fortement associée à An. gambiae s.s. (Chi carré = 23,41; P<0,0001) et L1014F était associée à An. arabiensis (chi carré = 11,21; P = 0,0008). L'allèle L1014S était significativement associé aux moustiques résistants à la lambdacyhalothrine (Fisher P exact <0,001).

Conclusion

La cooccurrence des mutations L1014S et L1014F couplées à des rapports sur la résistance aux insecticides suggèrent que la résistance aux pyréthrinoïdes est en train de devenir un phénomène répandu dans les populations de vecteurs du paludisme en Tanzanie. La présence de la mutation L1014F dans cette population de moustiques en Afrique de l'Est indique la propagation de ce gène à travers l'Afrique. L'investigation des implications opérationnelles potentielles de ces résultats sur le contrôle du paludisme devraient être approfondie.

Objetivo

Los marcadores moleculares de resistencia a insecticidas pueden ser indicadores sensibles del desarrollo de resistencias en las poblaciones de los vectores Anopheles. Evaluar dichos marcadores es crucial para los programas de monitorización de resistencias. Hemos investigado la presencia y la distribución de las mutaciones de resistencia knockdown (kdr) en Anopheles gambiae s.l. en Tanzania.

Métodos

Se recolectaron mosquitos Anopheles intradomiciliarios de 10 lugares diferentes y se evaluaron en busca de resistencia a insecticidas utilizando el protocolo estándar de la OMS. Mediante un diagnóstico molecular basado en la PCR se genotiparon los mosquitos y se detectaron los genotipos kdr.

Resultados

Los An. gambiae evaluados eran resistentes a lambdacialotrina en Muheza, Arumeru y Muleba. De 350 An. gambiae s.l. genotipados, 35% eran An. gambiae s.s. y 65% eran An. arabiensis. Se detectaron mutaciones L1014S y L1014F tanto en An. gambiae s.s. como en An. arabiensis. La mutación puntual L1014S se encontró con una frecuencia alélica de 4-33%, mientras que L1014F tenía una frecuencia alélica de 6-14%. La mutación L1014S estaba ampliamente asociada a An. gambiae s.s. (Chi-Cuadrado = 23.41; P < 0.0001) y la L1014F estaba asociada con An. arabiensis (Chi-Square = 11.21; P = 0.0008). El alelo L1014S estaba significativamente asociado con mosquitos resistentes a la lambdacialotrina (P < 0.001).

Conclusión

La simultaneidad de mutaciones de L1014S y L1014F junto con informes de resistencia a los insecticidas sugiere que la resistencia a piretroides se está convirtiendo en un fenómeno común entre las poblaciones del vector de la malaria en Tanzania. La presencia de la mutación L1014F en estas poblaciones del Este de África indican la diseminación del gen a lo largo del continente africano. Determinar las implicaciones potenciales a nivel operativo de estos hallazgos sobre el control de la malaria requiere de más estudios.

DDT and pyrethroid resistance in Anopheles arabiensis from South Africa

Background

Pyrethroid resistance has been well documented in Anopheles arabiensis, one of the major African malaria vectors, and the predominant malaria vector in South Africa.

Methods

In this study, the genetic basis of pyrethroid resistance in a selected laboratory strain of An. arabiensis from South Africa was investigated using a custom-made microarray, known as the An. gambiae detoxification chip.

Results

A large number of P450 genes were over-transcribed, as well as a suite of redox genes and glutathione S-transferases. The five genes that showed the highest level of gene transcription when compared with an insecticide susceptible strain were: CYP6AG2, CYPZ1, TPX2, CYPZ2 and CYP6P1.

Conclusions

Permethrin resistance in South African An. arabiensis is associated with increased transcription of multiple genes, and a large proportion of these genes were also previously recorded as over-transcribed in another An. arabiensis strain selected for resistance to DDT with cross-resistance to deltamethrin. The deltamethrin resistance developed de novo in the DDT-selected strain and is most likely due to increased transcription of those genes associated with DDT resistance. However, of particular interest was the fact that the strain selected for resistance to pyrethroids did not develop de novo resistance to DDT. These differences are compared and discussed.

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.

Marked biological differences between insecticide resistant and susceptible strains of Anopheles funestus infected with the murine parasite Plasmodium berghei

Background

Anopheles funestus is one of the major malaria vectors in Africa but research on this species has been restricted due to the lack of viable laboratory colonies. The vectorial capacity of natural populations of An. funestus is well known but its ability to host Plasmodium in the laboratory and the development cycle of the parasite within this mosquito species was, until very recently, unknown. In this study we compared laboratory strains of An. funestus that were resistant and susceptible to pyrethroid insecticides, for their receptiveness to infection with Plasmodium berghei and compared development times with other vector species available in our laboratory.

Methods

The murine parasite P. berghei was used to infect two base An. funestus colonies (FANG and FUMOZ) and two selected sub-colonies with different degrees of pyrethroid resistance (FUMOZ-BS susceptible and FUMOZ-R resistant). Results were compared with the G3 strain of An. gambiae.

Results

While all colonies were able to support the parasite, the development time in An. funestus was generally longer than that recorded in the laboratory strain of An. gambiae. Infected females were able to initiate new rounds of infection when feeding on healthy mice. The pyrethroid resistant strain FUMOZ-R supported the lowest numbers of oocysts and sporozoites while the insecticide susceptible strain FUMOZ-BS produced one of the highest sporozoite indices ever documented in P. berghei research. The FUMOZ base colony, exhibiting partial insecticide resistance was the median in terms of infection intensity. The oocyst number in all colonies did not fully correlate with the sporozoite index, indicating possible factors influencing the sporozoites’ transit from the midgut to the salivary glands.

Conclusions

The presence of both insecticide resistance and limited parasite infection phenotypes in the same individuals suggests there may be association between the two mechanisms, but further elucidation is required.

Country-level operational implementation of the Global Plan for Insecticide Resistance Management

Malaria control is reliant on the use of long-lasting pyrethroid-impregnated nets and/or indoor residual spraying (IRS) of insecticide. The rapid selection and spread of operationally significant pyrethroid resistance in African malaria vectors threatens our ability to sustain malaria control. Establishing whether resistance is operationally significant is technically challenging. Routine monitoring by bioassay is inadequate, and there are limited data linking resistance selection with changes in disease transmission. The default is to switch insecticides when resistance is detected, but limited insecticide options and resistance to multiple insecticides in numerous locations make this approach unsustainable. Detailed analysis of the resistance situation in Anopheles gambiae on Bioko Island after pyrethroid resistance was detected in this species in 2004, and the IRS program switched to carbamate bendiocarb, has now been undertaken. The pyrethroid resistance selected is a target-site knock-down resistance kdr-form, on a background of generally elevated metabolic activity, compared with insecticide-susceptible A. gambiae, but the major cytochrome P450-based metabolic pyrethroid resistance mechanisms are not present. The available evidence from bioassays and infection data suggests that the pyrethroid resistance mechanisms in Bioko malaria vectors are not operationally significant, and on this basis, a different, long-lasting pyrethroid formulation is now being reintroduced for IRS in a rotational insecticide resistance management program. This will allow control efforts to be sustained in a cost-effective manner while reducing the selection pressure for resistance to nonpyrethroid insecticides. The methods used provide a template for evidence-based insecticide resistance management by malaria control programs.

Impact of environment on mosquito response to pyrethroid insecticides: facts, evidences and prospects

By transmitting major human diseases such as malaria, dengue fever and filariasis, mosquito species represent a serious threat worldwide in terms of public health, and pose a significant economic burden for the African continent and developing tropical regions. Most vector control programmes aiming at controlling life-threatening mosquitoes rely on the use of chemical insecticides, mainly belonging to the pyrethroid class. However, resistance of mosquito populations to pyrethroids is increasing at a dramatic rate, threatening the efficacy of control programmes throughout insecticide-treated areas, where mosquito-borne diseases are still prevalent. In the absence of new insecticides and efficient alternative vector control methods, resistance management strategies are therefore critical, but these require a deep understanding of adaptive mechanisms underlying resistance. Although insecticide resistance mechanisms are intensively studied in mosquitoes, such adaptation is often considered as the unique result of the selection pressure caused by insecticides used for vector control. Indeed, additional environmental parameters, such as insecticides/pesticides usage in agriculture, the presence of anthropogenic or natural xenobiotics, and biotic interactions between vectors and other organisms, may affect both the overall mosquito responses to pyrethroids and the selection of resistance mechanisms. In this context, the present work aims at updating current knowledge on pyrethroid resistance mechanisms in mosquitoes and compiling available data, often from different research fields, on the impact of the environment on mosquito response to pyrethroids. Key environmental factors, such as the presence of urban or agricultural pollutants and biotic interactions between mosquitoes and their microbiome are discussed, and research perspectives to fill in knowledge gaps are suggested.

The importance of mosquito behavioural adaptations to malaria control in Africa

Over the past decade the use of long-lasting insecticidal nets (LLINs), in combination with improved drug therapies, indoor residual spraying (IRS), and better health infrastructure, has helped reduce malaria in many African countries for the first time in a generation. However, insecticide resistance in the vector is an evolving threat to these gains. We review emerging and historical data on behavioral resistance in response to LLINs and IRS. Overall the current literature suggests behavioral and species changes may be emerging, but the data are sparse and, at times unconvincing. However, preliminary modeling has demonstrated that behavioral resistance could have significant impacts on the effectiveness of malaria control. We propose seven recommendations to improve understanding of resistance in malaria vectors. Determining the public health impact of physiological and behavioral insecticide resistance is an urgent priority if we are to maintain the significant gains made in reducing malaria morbidity and mortality.

Identifying permethrin resistance loci in malaria vectors by genetic mapping

Identification of the major loci responsible for insecticide resistance in malaria vectors would aid the development and implementation of effective resistance management strategies, which are urgently needed to tackle the growing threat posed by resistance to the limited insecticides available for malaria control. Genome-wide association studies in the major malaria vector, Anopheles gambiae, have been hindered by the high degree of within-population structuring and very low levels of linkage disequilibrium hence we revisited the use of quantitative trait loci (QTL) mapping to study resistance phenotypes in this vector species. Earlier work, identified two major QTL associated with pyrethroid resistance in A. gambiae s.s. from East Africa using genetic crossing of laboratory-colonized resistant and susceptible strains. In this study, we report the results from genetic mapping of pyrethroid resistance in three isofemale pedigrees established from wild-caught female A. gambiae s.s. mosquitoes from Benin. We identified two QTL on chromosomes 2L and 3R in these field populations, in similar genomic locations to the QTL identified in laboratory strains. The relative merits of two alternative study designs are discussed and suggestions made for future genetic mapping studies of insecticide resistance in mosquitoes.

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.

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.

Cytochrome P450 6M2 from the malaria vector Anopheles gambiae metabolizes pyrethroids: Sequential metabolism of deltamethrin revealed

Resistance to pyrethroid insecticides in the malaria vector Anopheles gambiae is a major threat to malaria control programmes. Cytochome P450-mediated detoxification is an important resistance mechanism. CYP6M2 is over-expressed in wild populations of permethrin resistant A. gambiae but its role in detoxification is not clear. CYP6M2 was expressed in Escherichia coli and a structural model was produced to examine its role in pyrethroid metabolism. Both permethrin and deltamethrin were metabolized. Rates were enhanced by A. gambiae cytochrome b(5) with kinetic parameters of K(M)=11±1μM and k(cat)=6.1±0.4 per min for permethrin (1:1 cis-trans) and K(M)=2.0±0.3μM and k(cat)=1.2±0.1 per min for deltamethrin. Mass spectrometry and NMR analysis identified 4'-hydroxy deltamethrin and hydroxymethyl deltamethrin as major and minor deltamethrin metabolites respectively. Secondary breakdown products included cyano(3-hydroxyphenyl)methyl deltamethrate and deltamethric acid. CYP6M2 was most highly transcribed in the midgut and Malpighian tubules of adult A. gambiae, consistent with a role in detoxification. Our data indicates that CYP6M2 plays an important role in metabolic resistance to pyrethroids and thus an important target for the design of new tools to combat malaria.

Identification and characterisation of ten glutathione S-transferase genes from oriental migratory locust, Locusta migratoria manilensis (Meyen)

BACKGROUND

Synthetic pyrethroids are the primary insecticides that are widely used for controlling Locusta migratoria manilensis (Meyen), a major pest in eastern and southern Asia and the Pacific region. In this paper, ten cDNAs encoding glutathione S-transferases (GSTs) were sequenced and characterised in L. migratoria manilensis. The effects of deltamethrin on the ten GST gene expressions were studied.

RESULTS

Phylogenetic analysis revealed nine GSTs in three different classes, including seven in sigma, one in delta and one in theta. The remaining GST (LmGSTu1) was unclassified. RT-PCR analysis showed that most GST genes were expressed in all tissues examined, including the foregut, midgut, gastric caecum, hindgut, Malpighian tubules, fat bodies, muscles, spermaries and ovaries, except that LmGSTs2, LmGSTs4, LmGSTs7 and LmGSTu1 were expressed in several tissues. LmGSTu1 appeared to be the only gene whose expressions could not be detected in eggs. Real-time quantitative PCR showed that deltamethrin at 0.08 and/or 0.12 µg mL⁻¹ increased almost all ten GST gene expressions in third-instar nymph locusts. However, deltamethrin at 0.16 and/or 0.2 µg mL⁻¹ decreased the expressions of LmGSTd1, LmGSTs1, LmGSTs5 and LmGSTs6.

CONCLUSION

The increases in GST gene expressions after deltamethrin exposure in L. migratoria manilensis might result in its elevating tolerance to other insecticides and xenobiotics.

Cross-amplified microsatellites in the European cherry fly, Rhagoletis cerasi: medium polymorphic-highly informative markers

The European cherry fruit fly, Rhagoletis cerasi (Diptera: Tephritidae), is a major pest of cherries in Europe and parts of Asia. Despite its big economic significance, there is a lack of studies on the genetic structure of its natural populations. Knowledge about an insect pest on molecular, genetic and population levels facilitates the development of environmentally friendly control methods. In this study, we present the development of 13 microsatellite markers for R. cerasi, through cross-species amplification. These markers have been used for the genotyping of 130 individuals from five different sampling sites in Greece. Our results indicate that (i) cross-species amplification is a versatile and rapid tool for developing microsatellite markers in Rhagoletis spp., (ii) the microsatellite markers presented here constitute an important tool for population studies on this pest, and (iii) there is clear structuring of natural European cherry fly populations.

How can we identify parasite genes that underlie antimalarial drug resistance?

This article outlines genome-scale approaches that can be used to identify mutations in malaria (Plasmodium) parasites that underlie drug resistance and contribute to treatment failure. These approaches include genetic mapping by linkage or genome-wide association studies, drug selection and characterization of resistant mutants, and the identification of genome regions under strong recent selection. While these genomic approaches can identify candidate resistance loci, genetic manipulation is needed to demonstrate causality. We therefore also describe the growing arsenal of available transfection approaches for direct incrimination of mutations suspected to play a role in resistance. Our intention is both to review past progress and highlight promising approaches for future investigations.

Patterns and seasonality of malaria transmission in the forest-savannah transitional zones of Ghana

Background

Knowledge of the local pattern of malaria transmission and the effect of season on transmission is essential for the planning and evaluation of malaria interventions. Therefore, entomological surveys were carried out in the forest-savannah transitional belt of Ghana (Kintampo) from November 2003 to November 2005 in preparation for drug and vaccine trials.

Results

A total of 23,406 mosquitoes were caught from 919 traps over the two-year period (November 2003 to November 2005): 54.3% were Culicines, 36.2% Anopheles funestus, and 9.4% Anopheles gambiae. Infection rates with Plasmodium falciparum were 4.7% and 1.5% for Anopheles gambiae and Anopheles funestus, respectively. Entomological inoculation rates (EIRs) were 269 infective bites per person per year in the first year (November 2003-October 2004) and 231 the following year (November 2004-November 2005). Polymerase Chain Reaction (PCR) analysis detected only Anopheles gambiae s.s. Nineteen mosquitoes were tested by PCR in the wet season; 16 were S-molecular form, 2 M-molecular form and 1 hybrid (S/M). In the dry season, sixteen mosquitoes were tested; 11 S-molecular form, 2 M-molecular form and 3 S/M hybrids. The frequency of knock down resistance (kdr) genotypes F(R) was 0.60.

Conclusion

The dynamics and seasonal abundance of malaria vectors in the Kintampo area was influenced by micro-ecology, rainfall and temperature patterns. Transmission patterns did not differ significantly between the two years (2004 and 2005) and both Anopheles gambiae and Anopheles funestus were identified as effective vectors. EIR estimates in 2004/2005 were between 231 and 269 infective bites per person per year. The information provided by the study will help in planning intensified malaria control activities as well as evaluating the impact of malaria interventions in the middle belt of Ghana.

Major effect genes or loose confederations? The development of insecticide resistance in the malaria vector Anopheles gambiae

Insecticide use in public health and agriculture presents a dramatic adaptive challenge to target and non-target insect populations. The rapid development of genetically modulated resistance to insecticides is postulated to develop in two distinct ways: By selection for single major effect genes or by selection for loose confederations in which several factors, not normally associated with each other, inadvertently combine their effects to produce resistance phenotypes. Insecticide resistance is a common occurrence and has been intensively studied in the major malaria vector Anopheles gambiae, providing a useful model for examining how insecticide resistance develops and what pleiotropic effects are likely to emerge as a consequence of resistance. As malaria vector control becomes increasingly reliant on successfully managing insecticide resistance, the characterisation of resistance mechanisms and their pleiotropic effects becomes increasingly important.

Does kdr genotype predict insecticide-resistance phenotype in mosquitoes?

Several groups are developing and applying DNA-based technologies to monitor insecticide-based disease control programmes. However, several recent papers have concluded that the knockdown resistance (kdr) genotype-phenotype correlation that is observed in a wide variety of taxa might not hold in all mosquitoes. In this article, we review the evidence to support this putative breakdown and argue that the conclusion follows from unreliable data or the unparsimonious interpretation of data. We assert that the link between kdr genotype and DDT- and pyrethroid-susceptibility phenotype is clear. However, we emphasize that kdr genotype might explain only a portion of heritable variation in resistance and that diagnostic assays to test the importance of other resistance mechanisms in field populations are required.

High, clustered, nucleotide diversity in the genome of Anopheles gambiae revealed through pooled-template sequencing: implications for high-throughput genotyping protocols

Background

Association mapping approaches are dependent upon discovery and validation of single nucleotide polymorphisms (SNPs). To further association studies in Anopheles gambiae we conducted a major resequencing programme, primarily targeting regions within or close to candidate genes for insecticide resistance.

Results

Using two pools of mosquito template DNA we sequenced over 300 kbp across 660 distinct amplicons of the An. gambiae genome. Comparison of SNPs identified from pooled templates with those from individual sequences revealed a very low false positive rate. False negative rates were much higher and mostly resulted from SNPs with a low minor allele frequency. Pooled-template sequencing also provided good estimates of SNP allele frequencies. Allele frequency estimation success, along with false positive and negative call rates, improved significantly when using a qualitative measure of SNP call quality. We identified a total of 7062 polymorphic features comprising 6995 SNPs and 67 indels, with, on average, a SNP every 34 bp; a high rate of polymorphism that is comparable to other studies of mosquitoes. SNPs were significantly more frequent in members of the cytochrome p450 mono-oxygenases and carboxy/cholinesterase gene-families than in glutathione-S-transferases, other detoxification genes, and control genomic regions. Polymorphic sites showed a significantly clustered distribution, but the degree of SNP clustering (independent of SNP frequency) did not vary among gene families, suggesting that clustering of polymorphisms is a general property of the An. gambiae genome.

Conclusion

The high frequency and clustering of SNPs has important ramifications for the design of high-throughput genotyping assays based on allele specific primer extension or probe hybridisation. We illustrate these issues in the context of the design of Illumina GoldenGate assays.

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

The Yin and Yang of linkage disequilibrium: mapping of genes and nucleotides conferring insecticide resistance in insect disease vectors

Genetic technologies developed in the last 20 years have lead to novel and exciting methods to identify genes and specific nucleotides within genes that control phenotypes in field collected organisms. In this review we define and explain two of these methods: linkage disequilibrium (LD) mapping and quantitative trait nucleotide (QTN) mapping. The power to detect valid genotype-phenotype associations with LD or QTN mapping depends critically on the extent to which segregating sites in a genome assort independently. LD mapping depends on markers being in disequilibrium with the genes that condition expression of the phenotype. In contrast, QTN mapping depends critically upon most proximal loci being at equilibrium. We show that both patterns actually exist in the genome of Anapheles gambiae, the most important malaria vector in sub-Saharan Africa while segregating sites appear to be largely in equilibrium throughout the genome of Aedes aegypti, the vector of Dengue and Yellow fever flaviviruses. We discuss additional approaches that will be needed to identify genes and nucleotides that control phenotypes in field collected organisms, focusing specifically on ongoing studies of genes conferring resistance to insecticides.

Characterization of inhibitors and substrates of Anopheles gambiae CYP6Z2

Three CYP6Z genes are linked to a major pyrethroid resistance locus in the mosquito Anopheles gambiae. We have expressed CYP6Z2 in Escherichia coli and produced a structural model in order to examine its role in detoxification. E. coli membranes co-expressing CYP6Z2 and An. gambiae P450 reductase (AgCPR) catalysed the dealkylation of benzyloxyresorufin with kinetic parameters K(m) = 0.13 microM; K(cat) = 1.5 min(-1). The IC(50) values of a wide range of compounds were measured. Pyrethroids cypermethrin and permethrin produced low IC(50) values, but were not metabolized. Plant flavanoids were the most potent inhibitors. Several compounds were shown to be substrates, suggesting that CYP6Z2 has broad substrate specificity and plays an important chemo-protective role during the herbivorous phase of the life-cycle.

Molecular evolution of glutathione S-transferases in the genus Drosophila

As classical phase II detoxification enzymes, glutathione S-transferases (GSTs) have been implicated in insecticide resistance and may have evolved in response to toxins in the niche-defining feeding substrates of Drosophila species. We have annotated the GST genes of the 12 Drosophila species with recently sequenced genomes and analyzed their molecular evolution. Gene copy number variation is attributable mainly to unequal crossing-over events in the large delta and epsilon clusters. Within these gene clusters there are also GST genes with slowly diverging orthologs. This implies that they have their own unique functions or have spatial/temporal expression patterns that impose significant selective constraints. Searches for positively selected sites within the GSTs identified G171K in GSTD1, a protein that has previously been shown to be capable of metabolizing the insecticide DDT. We find that the same radical substitution (G171K) in the substrate-binding domain has occurred at least three times in the Drosophila radiation. Homology-modeling places site 171 distant from the active site but adjacent to an alternative DDT-binding site. We propose that the parallel evolution observed at this site is an adaptive response to an environmental toxin and that sequencing of historical alleles suggests that this toxin was not a synthetic insecticide.

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

BACKGROUND

Pyrethroid resistance in Anopheles funestus populations has led to an increase in malaria transmission in southern Africa. Resistance has been attributed to elevated activities of cytochrome P450s but the molecular basis underlying this metabolic resistance is unknown. Microsatellite and SNP markers were used to construct a linkage map and to detect a quantitative trait locus (QTL) associated with pyrethroid resistance in the FUMOZ-R strain of An. funestus from Mozambique.

RESULTS

By genotyping 349 F2 individuals from 11 independent families, a single major QTL, rp1, at the telomeric end of chromosome 2R was identified. The rp1 QTL appears to present a major effect since it accounts for more than 60% of the variance in susceptibility to permethrin. This QTL has a strong additive genetic effect with respect to susceptibility. Candidate genes associated with pyrethroid resistance in other species were physically mapped to An. funestus polytene chromosomes. This showed that rp1 is genetically linked to a cluster of CYP6 cytochrome P450 genes located on division 9 of chromosome 2R and confirmed earlier reports that pyrethroid resistance in this strain is not associated with target site mutations (knockdown resistance).

CONCLUSION

We hypothesize that one or more of these CYP6 P450s clustered on chromosome 2R confers pyrethroid resistance in the FUMOZ-R strain of An. funestus.

Transcription profiling of a recently colonised pyrethroid resistant Anopheles gambiae strain from Ghana

BACKGROUND

Mosquito resistance to the pyrethroid insecticides used to treat bednets threatens the sustainability of malaria control in sub-Saharan Africa. While the impact of target site insensitivity alleles is being widely discussed the implications of insecticide detoxification--though equally important--remains elusive. The successful development of new tools for malaria intervention and management requires a comprehensive understanding of insecticide resistance, including metabolic resistance mechanisms. Although three enzyme families (cytochrome P450s, glutathione S-transferases and carboxylesterases) have been widely associated with insecticide detoxification the role of individual enzymes is largely unknown.

RESULTS

Here, constitutive expression patterns of genes putatively involved in conferring pyrethroid resistance was investigated in a recently colonised pyrethroid resistant Anopheles gambiae strain from Odumasy, Southern Ghana. RNA from the resistant strain and a standard laboratory susceptible strain, of both sexes was extracted, reverse transcribed and labelled with either Cy3- or Cy5-dye. Labelled cDNA was co-hybridised to the detox chip, a custom-made microarray containing over 230 A. gambiae gene fragments predominantly from enzyme families associated with insecticide resistance. After hybridisation, Cy3- and Cy5-signal intensities were measured and compared gene by gene. In both females and males of the resistant strain the cytochrome P450s CYP6Z2 and CYP6M2 are highly over-expressed along with a member of the superoxide dismutase (SOD) gene family.

CONCLUSION

These genes differ from those found up-regulated in East African strains of pyrethroid resistant A. gambiae and constitute a novel set of candidate genes implicated in insecticide detoxification. These data suggest that metabolic resistance may have multiple origins in A. gambiae, which has strong implications for the management of resistance.

Gene expression in insecticide resistant and susceptible Anopheles gambiae strains constitutively or after insecticide exposure

A microarray containing approximately 20 000 expressed sequence tags (ESTs; 11 760 unique EST clusters) from the malaria vector, Anopheles gambiae, was used to monitor differences in global gene expression in two insecticide resistant and one susceptible strains. Statistical analysis identified 77 ESTs that were differentially transcribed among the three strains. These include the cytochrome P450 CYP314A1, over-transcribed in the DDT resistant ZAN/U strain, and many genes that belong to families not usually associated with insecticide resistance, such as peptidases, sodium/calcium exchangers and genes implicated in lipid and carbohydrate metabolism. Short-term (6 and 10 h) effects of exposure of the pyrethroid resistant RSP strain to permethrin were also detected. Several genes belonging to enzyme families already implicated in insecticide or xenobiotic detoxification were induced, including the carboxylesterase COEAE2F gene and members of the UDP-glucuronosyl transferase and nitrilase families.

The molecular basis of insecticide resistance in mosquitoes

Insecticide resistance is an inherited characteristic involving changes in one or more insect gene. The molecular basis of these changes are only now being fully determined, aided by the availability of the Drosophila melanogaster and Anopheles gambiae genome sequences. This paper reviews what is currently known about insecticide resistance conferred by metabolic or target site changes in mosquitoes.

Malaria: past problems and future prospects. After more than a decade of neglect, malaria is finally black on the agenda for both biomedical research and public health politics

After more than a decade of neglect, malaria is finally back on the agenda for both biomedical research and public health politics