Transcriptomic analysis of resistance and short-term induction response to pyrethroids, in Anopheles coluzzii legs.

BACKGROUND: Insecticide-treated bed nets and indoor residual spraying comprise the major control measures against Anopheles gambiae sl, the dominant vector in sub-Saharan Africa. The primary site of contact with insecticide is through the mosquitoes' legs, which represents the first barrier insecticides have to bypass to reach their neuronal targets. Proteomic changes and leg cuticle modifications have been associated with insecticide resistance that may reduce the rate of penetration of insecticides. Here, we performed a multiple transcriptomic analyses focusing on An. coluzzii legs. RESULTS: Firstly, leg-specific enrichment analysis identified 359 genes including the pyrethroid-binder SAP2 and 2 other chemosensory proteins, along with 4 ABCG transporters previously shown to be leg enriched. Enrichment of gene families included those involved in detecting chemical stimuli, including gustatory and ionotropic receptors and genes implicated in hydrocarbon-synthesis. Subsequently, we compared transcript expression in the legs of a highly resistant strain (VK7-HR) to both a strain with very similar genetic background which has reverted to susceptibility after several generations without insecticide pressure (VK7-LR) and a lab susceptible population (NG). Two hundred thirty-two differentially expressed genes (73 up-regulated and 159 down-regulated) were identified in the resistant strain when compared to the two susceptible counterparts, indicating an over-expression of phase I detoxification enzymes and cuticular proteins, with decrease in hormone-related metabolic processes in legs from the insecticide resistant population. Finally, we analysed the short-term effect of pyrethroid exposure on An. coluzzii legs, comparing legs of 1 h-deltamethrin-exposed An. coluzzii (VK7-IN) to those of unexposed mosquitoes (VK7-HR) and identified 348 up-regulated genes including those encoding for GPCRs, ABC transporters, odorant-binding proteins and members of the divergent salivary gland protein family. CONCLUSIONS: The data on An. coluzzii leg-specific transcriptome provides valuable insights into the first line of defense in pyrethroid resistant and short-term deltamethrin-exposed mosquitoes. Our results suggest that xenobiotic detoxification is likely occurring in legs, while the enrichment of sensory proteins, ABCG transporters and cuticular genes is also evident. Constitutive resistance is primarily associated with elevated levels of detoxification and cuticular genes, while short-term insecticide-induced tolerance is linked with overexpression of transporters, GPCRs and GPCR-related genes, sensory/binding and salivary gland proteins.

Transcriptomic analysis reveals pronounced changes in gene expression due to sub-lethal pyrethroid exposure and ageing in insecticide resistance Anopheles coluzzii.

BACKGROUND: Malaria control is heavily reliant on the use of insecticides that target and kill the adult female Anopheline vector. The intensive use of insecticides of the pyrethroid class has led to widespread resistance in mosquito populations. The intensity of pyrethroid resistance in some settings in Africa means mosquitoes can contact bednets treated with this insecticide class multiple times with minimal mortality effects. Furthermore, both ageing and diel cycle have been shown to have large impacts on the resistance phenotype. Together, these traits may affect other aspects of vector biology controlling the vectorial capacity or fitness of the mosquito. RESULTS: Here we show that sublethal exposure of a highly resistant Anopheles coluzzii population originally from Burkina Faso to the pyrethroid deltamethrin results in large and sustained changes to transcript expression. We identify five clear patterns in the data showing changes to transcripts relating to: DNA repair, respiration, translation, behaviour and oxioreductase processes. Further, we highlight differential regulation of transcripts from detoxification families previously linked with insecticide resistance, in addition to clear down-regulation of the oxidative phosphorylation pathway both indicative of changes in metabolism post-exposure. Finally, we show that both ageing and diel cycle have major effects on known insecticide resistance related transcripts. CONCLUSION: Sub-lethal pyrethroid exposure, ageing and the diel cycle results in large-scale changes in the transcriptome of the major malaria vector Anopheles coluzzii. Our data strongly supports further phenotypic studies on how transcriptional changes such as reduced expression of the oxidative phosphorylation pathway or pyrethroid induced changes to redox state might impact key mosquito traits, such as vectorial capacity and life history traits.

Transcriptomic meta-signatures identified in Anopheles gambiae populations reveal previously undetected insecticide resistance mechanisms.

Increasing insecticide resistance in malaria-transmitting vectors represents a public health threat, but underlying mechanisms are poorly understood. Here, a data integration approach is used to analyse transcriptomic data from comparisons of insecticide resistant and susceptible Anopheles populations from disparate geographical regions across the African continent. An unbiased, integrated analysis of this data confirms previously described resistance candidates but also identifies multiple novel genes involving alternative resistance mechanisms, including sequestration, and transcription factors regulating multiple downstream effector genes, which are validated by gene silencing. The integrated datasets can be interrogated with a bespoke Shiny R script, deployed as an interactive web-based application, that maps the expression of resistance candidates and identifies co-regulated transcripts that may give clues to the function of novel resistance-associated genes.

Do bednets including piperonyl butoxide offer additional protection against populations of Anopheles gambiae s.l. that are highly resistant to pyrethroids? An experimental hut evaluation in Burkina Fasov.

Malaria control is dependent on the use of longlasting insecticidal nets (LLINs) containing pyrethroids. A new generation of LLINs containing both pyrethroids and the synergist piperonyl butoxide (PBO) has been developed in response to increasing pyrethroid resistance in African malaria vectors, but questions remain about the performance of these nets in areas where levels of pyrethroid resistance are very high. This study was conducted in two settings in southwest Burkina Faso, Vallée du Kou 5 and Tengrela, where Anopheles gambiae s.l. (Diptera: Culicidae) mortality rates in World Health Organization (WHO) discriminating dose assays were < 14% for permethrin and < 33% for deltamethrin. When mosquitoes were pre-exposed to PBO in WHO tube assays, mortality rates increased substantially but full susceptibility was not restored. Molecular characterization revealed high levels of kdr alleles and elevated levels of P450s previously implicated in pyrethroid resistance. In cone bioassays and experimental huts, PBO LLINs outperformed the pyrethroid-only equivalents from the same manufacturers. Blood feeding rates were 1.6-2.2-fold lower and mortality rates were 1.69-1.78-fold greater in huts with PBO LLINs vs. non-PBO LLINs. This study indicates that PBO LLINs provide greater personal and community-level protection than standard LLINs against highly pyrethroid-resistant mosquito populations.

The Anopheles gambiae ATP-binding cassette transporter family: phylogenetic analysis and tissue localization provide clues on function and role in insecticide resistance.

The role of ATP-binding cassette (ABC) transporters in conferring insecticide resistance has received much attention recently. Here we identify ABC transporters differentially expressed in insecticide-resistant populations of the malaria vector, Anopheles gambiae. Although we found little evidence that the orthologues of the multidrug resistance proteins described in other species are associated with resistance in An. gambiae we did identify a subset of ABC proteins consistently differentially expressed in pyrethroid-resistant populations from across Africa. We present information on the phylogenetic relationship, primary sites of expression and potential role of ABC transporters in mediating the mosquito's response to insecticides. Furthermore we demonstrate that a paralogous group of eight ABCG transporters, clustered on chromosome 3R, are highly enriched in the legs of An. gambiae mosquitoes, consistent with a proposed role for this ABC subfamily in transport of lipids to the outer surface of the cuticle. Finally, antibodies raised against one of the most highly expressed ABC transporters in adult females, ABCG7 (AGAP009850), localized this transporter to the pericardial cells. These data will help prioritize members of this gene family for further localization and functional validation studies to identify the in vivo function of these transporters in the mosquito and determine whether elevated expression of members of this family contribute to insecticide resistance.

Evaluation of a temperate climate mosquito, Ochlerotatus detritus (=Aedes detritus), as a potential vector of Japanese encephalitis virus.

The U.K. has not yet experienced a confirmed outbreak of mosquito-borne virus transmission to people or livestock despite numerous autochthonous epizootic and human outbreaks of mosquito-borne diseases on the European mainland. Indeed, whether or not British mosquitoes are competent to transmit arboviruses has not been established. Therefore, the competence of a local (temperate) British mosquito species, Ochlerotatus detritus (=Aedes detritus) (Diptera: Culicidae) for transmission of a member of the genus Flavivirus, Japanese encephalitis virus (JEV) as a model for mosquito-borne virus transmission was assessed. The JEV competence in a laboratory strain of Culex quinquefasciatus (Diptera: Culicidae), a previously incriminated JEV vector, was also evaluated as a positive control. Ochlerotatus detritus adults were reared from field-collected juvenile stages. In oral infection bioassays, adult females developed disseminated infections and were able to transmit virus as determined by the isolation of virus in saliva secretions. When pooled at 7-21 days post-infection, 13% and 25% of O. detritus were able to transmit JEV when held at 23 °C and 28 °C, respectively. Similar results were obtained for C. quinquefasciatus. To our knowledge, this study is the first to demonstrate that a British mosquito species, O. detritus, is a potential vector of an exotic flavivirus.

Differential transcription profiles in Aedes aegypti detoxification genes after temephos selection.

The mosquito Aedes aegypti is the main vector of Dengue and Yellow Fever flaviviruses. The organophosphate insecticide temephos is a larvicide that is used globally to control Ae. aegypti populations; many of which have in turn evolved resistance. Target site alteration in the acetylcholine esterase of this species has not being identified. Instead, we tracked changes in transcription of metabolic detoxification genes using the Ae. aegypti 'Detox Chip' microarray during five generations of temephos selection. We selected for temephos resistance in three replicates in each of six collections, five from Mexico, and one from Peru. The response to selection was tracked in terms of lethal concentrations. Uniform upregulation was seen in the epsilon class glutathione-S-transferase (eGST) genes in strains from Mexico prior to laboratory selection, while eGSTs in the Iquitos Peru strain became upregulated after five generations of temephos selection. While expression of many carboxyl/cholinesterase esterase (CCE) genes increased with selection, no single esterase was consistently upregulated and this same pattern was noted in the cytochrome P450 monooxygenase (CYP) genes and in other genes involved in reduction or oxidation of xenobiotics. Bioassays using glutathione-S-transferase (GST), CCE and CYP inhibitors suggest that various CCEs instead of GSTs are the main metabolic mechanism conferring resistance to temephos. We show that temephos-selected strains show no cross resistance to permethrin and that genes associated with temephos selection are largely independent of those selected with permethrin in a previous study.

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.

Transcription of detoxification genes after permethrin selection in the mosquito Aedes aegypti.

Changes in gene expression before, during and after five generations of permethrin laboratory selection were monitored in six strains of Aedes aegypti: five F(2)-F(3) collections from the Yucatán Peninsula of Mexico and one F(2) from Iquitos, Peru. Three biological replicate lines were generated for each strain. The response to selection was measured as changes in the lethal and knockdown permethrin concentrations (LC(50), KC(50)) and in the frequency of the Ile1,016 substitution in the voltage-gated sodium channel (para) gene. Changes in expression of 290 metabolic detoxification genes were measured using the 'Aedes Detox' microarray. Selection simultaneously increased the LC(50), KC(50) and Ile1,016 frequency. There was an inverse relationship between Ile1,016 frequency and the numbers of differentially transcribed genes. The Iquitos strain lacked the Ile1,016 allele and 51 genes were differentially transcribed after selection as compared with 10-18 genes in the Mexican strains. Very few of the same genes were differentially transcribed among field strains but 10 cytochrome P(450) genes were upregulated in more than one strain. Laboratory adaptation to permethrin in Ae. aegypti is genetically complex and largely conditioned by geographic origin and pre-existing target site insensitivity in the para gene. The lack of uniformity in the genes that responded to artificial selection as well as differences in the direction of their responses challenges the assumption that one or a few genes control permethrin metabolic resistance. Attempts to identify one or a few metabolic genes that are predictably associated with permethrin adaptation may be futile.

Temephos resistance and esterase activity in the mosquito Aedes aegypti in Havana, Cuba increased dramatically between 2006 and 2008.

Aedes aegypti (L.) (Diptera: Culicidae) control programmes in Cuba rely on the application of the organophosphate temephos for larval control. Hence, the monitoring of resistance to this insecticide is an essential component of such programmes. Here, 15 field populations from different municipalities of Havana City were assayed for resistance to temephos. High levels of resistance were detected in all strains and resistance ratios were highly correlated with esterase activity (P = 0.00001). Populations from three municipalities were tested in both 2006 and 2008; resistance and esterase activities both significantly increased during this 2-year period. Synergist studies demonstrated that neither glutathione transferases nor monooxygenases were associated with the increase in resistance to temephos in this period. The duration of the efficacy of commercial formulations of temephos in controlling Ae. aegypti populations in Havana City was reduced by the high level of temephos resistance observed; hence these data are of clear operational significance for the dengue control programme in Cuba. New integrated strategies to avoid further increases in temephos resistance in Cuba are necessary.

Metabolic enzymes associated with xenobiotic and chemosensory responses in Nasonia vitripennis.

The numbers of glutathione S-transferase, cytochrome P450 and esterase genes in the genome of the hymenopteran parasitoid Nasonia vitripennis are about twice those found in the genome of another hymenopteran, the honeybee Apis mellifera. Some of the difference is associated with clades of these families implicated in xenobiotic resistance in other insects and some is in clades implicated in hormone and pheromone metabolism. The data support the hypothesis that the eusocial behaviour of the honeybee and the concomitant homeostasis of the nest environment may obviate the need for as many gene/enzyme systems associated with xenobiotic metabolism as are found in other species, including N. vitripennis, that are thought to encounter a wider range of potentially toxic xenobiotics in their diet and habitat.

Evidence of multiple pyrethroid resistance mechanisms in the malaria vector Anopheles gambiae sensu stricto from Nigeria.

Pyrethroid insecticide resistance in Anopheles gambiae sensu stricto is a major concern to malaria vector control programmes. Resistance is mainly due to target-site insensitivity arising from a single point mutation, often referred to as knockdown resistance (kdr). Metabolic-based resistance mechanisms have also been implicated in pyrethroid resistance in East Africa and are currently being investigated in West Africa. Here we report the co-occurrence of both resistance mechanisms in a population of An. gambiae s.s. from Nigeria. Bioassay, synergist and biochemical analysis carried out on resistant and susceptible strains of An. gambiae s.s. from the same geographical area revealed >50% of the West African kdr mutation in the resistant mosquitoes but <3% in the susceptible mosquitoes. Resistant mosquitoes synergized using pyperonyl butoxide before permethrin exposure showed a significant increase in mortality compared with the non-synergized. Biochemical assays showed an increased level of monooxygenase but not glutathione-S-transferase or esterase activities in the resistant mosquitoes. Microarray analysis using the An. gambiae detox-chip for expression of detoxifying genes showed five over-expressed genes in the resistant strain when compared with the susceptible one. Two of these, CPLC8 and CPLC#, are cuticular genes not implicated in pyrethroid metabolism in An. gambiae s.s, and could constitute a novel set of candidate genes that warrant further investigation.

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.

Over expression of a cytochrome P450 (CYP6P9) in a major African malaria vector, Anopheles Funestus, resistant to pyrethroids.

Anopheles funestus Giles is one of the major African malaria vectors. It has previously been implicated in a major outbreak of malaria in KwaZulu/Natal, South Africa, during the period 1996 to 2000. The re-emergence of this vector was associated with monooxygenase-based resistance to pyrethroid insecticides. We have identified a gene from the monooxygenase CYP6 family, CYP6P9, which is over expressed in a pyrethroid resistant strain originating from Mozambique. Quantitative Real-Time PCR shows that this gene is highly over expressed in the egg and adult stages of the resistant strain relative to the susceptible strain but the larval stages showed almost no difference in expression between strains. This gene is genetically linked to a major locus associated with pyrethroid resistance in this A. funestus population.

A mutation in the voltage-gated sodium channel gene associated with pyrethroid resistance in Latin American Aedes aegypti.

Pyrethroids are commonly used as mosquito adulticides and evolution of resistance to these compounds is a major threat to public health. 'Knockdown resistance' to pyrethroids (kdr) is frequently caused by nonsynonymous mutations in the voltage-gated sodium channel transmembrane protein (para) that reduce pyrethroid binding. Early detection of kdr is critical to the development of resistance management strategies in mosquitoes including Aedes aegypti, the most prevalent vector of dengue and yellow fever viruses. Brengues et al. described seven novel mutations in hydrophobic segment 6 of domain II of para in Ae. aegypti. Assays on larvae from strains bearing these mutations indicated reduced nerve sensitivity to permethrin inhibition. Two of these occurred in codons Iso1011 and Val1016 in exons 20 and 21 respectively. A transition in the third position of Iso1011 encoded a Met1011 replacement and a transversion in the second position of Val1016 encoded a Gly1016 replacement. We have screened this same region in 1318 mosquitoes in 32 additional strains; 30 from throughout Latin America. While the Gly1016 allele was never detected in Latin America, we found two new mutations in these same codons. A transition in the first position of codon 1011 encodes a Val replacement while a transition in the first position of codon 1016 encodes an Iso replacement. We developed PCR assays for these four mutations that can be read either on an agarose gel or as a melting curve. Selection experiments, one with deltamethrin on a field strain from Santiago de Cuba and another with permethrin on a strain from Isla Mujeres, Mexico rapidly increased the frequency of the Iso1016 allele. Bioassays of F(3) offspring arising from permethrin susceptible Val1016 homozygous parents and permethrin resistant Iso1016 homozygous parents show that Iso1016 segregates as a recessive allele in conferring kdr. Analysis of segregation between alleles at the 1011 and 1016 codons in the F(3) showed a high rate of recombination even though the two codons are only separated by a ~250 bp intron. The tools and information presented provide a means for early detection and characterization of kdr that is critical to the development of strategies for resistance management.

Transcriptional analysis of insecticide resistance in Anopheles stephensi using cross-species microarray hybridization.

A large scale microarray (20k MMC1) from the African malaria vector Anopheles gambiae was used to monitor gene expression in insecticide resistant and susceptible strains of the Asian mosquito Anopheles stephensi. Heterologous hybridization at slightly reduced stringency yielded approximately 7000 significant signals. Thirty-six putative genes were differentially transcribed between the pyrethroid-resistant (DUB-R) and the susceptible (BEECH) strains. The expression profiles of selected transcripts were verified by real-time PCR. A gene putatively involved in the thickening of the adult cuticle showed the most striking up-regulation in DUB-R. A more specialized microarray containing 231 An. gambiae genes putatively involved in insecticide detoxification was used to further analyse classical insecticide resistance genes. Three glutathione S-transferase (GST) transcripts, one esterase and a cytochrome P450 were up-regulated in the resistant strain, while two peroxidases were down-regulated.

A deficit of detoxification enzymes: pesticide sensitivity and environmental response in the honeybee.

The honeybee genome has substantially fewer protein coding genes ( approximately 11 000 genes) than Drosophila melanogaster ( approximately 13 500) and Anopheles gambiae ( approximately 14 000). Some of the most marked differences occur in three superfamilies encoding xenobiotic detoxifying enzymes. Specifically there are only about half as many glutathione-S-transferases (GSTs), cytochrome P450 monooxygenases (P450s) and carboxyl/cholinesterases (CCEs) in the honeybee. This includes 10-fold or greater shortfalls in the numbers of Delta and Epsilon GSTs and CYP4 P450s, members of which clades have been recurrently associated with insecticide resistance in other species. These shortfalls may contribute to the sensitivity of the honeybee to insecticides. On the other hand there are some recent radiations in CYP6, CYP9 and certain CCE clades in A. mellifera that could be associated with the evolution of the hormonal and chemosensory processes underpinning its highly organized eusociality.

Differential expression of the detoxification genes in the different life stages of the malaria vector Anopheles gambiae.

The diverse habitats and diets encountered during the life cycle of an Anopheles mosquito have necessitated the development of extensive families of detoxification enzymes. Expansion of the three detoxification enzyme families (cytochrome P450s, carboxylesterases and glutathione transfereases), has occurred in mosquitoes compared with Drosophila, however, very little is known regarding the developmental expression of theses genes. Using a custom made microarray we determined the expression profile of the detoxification genes in adults, larvae and pupae of the malaria vector A. gambiae. The expression of approximately one quarter of these genes was developmentally regulated. The expression profile of each of these genes and the information this data provides on putative functions of the mosquito detoxification enzymes is discussed.

Anopheles gambiae P450 reductase is highly expressed in oenocytes and in vivo knockdown increases permethrin susceptibility.

We describe an in vivo model for investigation of detoxification mechanisms of the mosquito Anopheles gambiae, important for the development of malaria control programmes. Cytochrome P450s are involved in metabolic insecticide resistance and require NADPH cytochrome P450 reductase (CPR) to function. Here we demonstrate that the major sites of adult mosquito CPR expression are oenocytes, mid-gut epithelia and head appendages. High CPR expression was also evident in Drosophila oenocytes indicating a general functional role in these insect cells. RNAi mediated knockdown drastically reduced CPR expression in oenocytes, and to a lesser extent in mid-gut epithelia; the head was unaffected. These flies showed enhanced sensitivity to permethrin, demonstrating a key role for abdominal/mid-gut P450s in pyrethroid metabolism, aiding the development of insecticides.