Genomic analysis of detoxification genes in the mosquito Aedes aegypti

Insecticide resistance mechanisms associated with different environments in the malaria vector Anopheles gambiae: a case study in Tanzania

BACKGROUND

Resistance of mosquitoes to insecticides is a growing concern in Africa. Since only a few insecticides are used for public health and limited development of new molecules is expected in the next decade, maintaining the efficacy of control programmes mostly relies on resistance management strategies. Developing such strategies requires a deep understanding of factors influencing resistance together with characterizing the mechanisms involved. Among factors likely to influence insecticide resistance in mosquitoes, agriculture and urbanization have been implicated but rarely studied in detail. The present study aimed at comparing insecticide resistance levels and associated mechanisms across multiple Anopheles gambiae sensu lato populations from different environments.

METHODS

Nine populations were sampled in three areas of Tanzania showing contrasting agriculture activity, urbanization and usage of insecticides for vector control. Insecticide resistance levels were measured in larvae and adults through bioassays with deltamethrin, DDT and bendiocarb. The distribution of An. gambiae sub-species and pyrethroid target-site mutations (kdr) were investigated using molecular assays. A microarray approach was used for identifying transcription level variations associated to different environments and insecticide resistance.

RESULTS

Elevated resistance levels to deltamethrin and DDT were identified in agriculture and urban areas as compared to the susceptible strain Kisumu. A significant correlation was found between adult deltamethrin resistance and agriculture activity. The subspecies Anopheles arabiensis was predominant with only few An. gambiae sensu stricto identified in the urban area of Dar es Salaam. The L1014S kdr mutation was detected at elevated frequency in An gambiae s.s. in the urban area but remains sporadic in An. arabiensis specimens. Microarrays identified 416 transcripts differentially expressed in any area versus the susceptible reference strain and supported the impact of agriculture on resistance mechanisms with multiple genes encoding pesticide targets, detoxification enzymes and proteins linked to neurotransmitter activity affected. In contrast, resistance mechanisms found in the urban area appeared more specific and more related to the use of insecticides for vector control.

CONCLUSIONS

Overall, this study confirmed the role of the environment in shaping insecticide resistance in mosquitoes with a major impact of agriculture activities. Results are discussed in relation to resistance mechanisms and the optimization of resistance management strategies.

Molecular cloning and xenobiotic induction of seven novel cytochrome P450 monooxygenases in Aedes albopictus

Cytochrome P450 monooxygenase (P450) is a superfamily of enzymes that is important in metabolism of endogenous and exogenous compounds. In insects, these enzymes confer resistance to insecticides through its metabolic activities. Members of P450 from family 6 in insects are known to play a role in such function. In this study, we have isolated seven novel family 6 P450 from Aedes albopictus (Skuse) (Diptera: Culicidae), a vector of dengue and chikungunya fever. Induction profile of these seven genes was studied using several insecticides and xenobiotics. It was found that deltamethrin and permethrin did not induce expression of any genes. Another insecticide, temephos, inhibited expression of CYP6P15 for fivefold and twofold for CYP6N29, CYP6Y7, and CYP6Z18. In addition, copper II sulfate induced expression of CYP6M17 and CYP6N28 for up to sixfold. Benzothiazole (BZT), a tire leachate induced the expression of CYP6M17 by fourfold, CYP6N28 by sevenfold, but inhibited the expression of CYP6P15 for threefold and CYP6Y7 for twofold. Meanwhile, piperonyl butoxide (PBO) induced the expression CYP6N28 (twofold), while it inhibited the expression of CYP6P15 (fivefold) and CYP6Y7 (twofold). Remarkably, all seven genes were induced two- to eightfold by acetone in larval stage, but not adult stage. Expression of CYP6N28 was twofold higher, while expression of CYP6P15 was 15-fold lower in adult than larva. The other five P450s were not differentially expressed between the larvae and adult. This finding showed that acetone can be a good inducer of P450 in Ae. albopictus. On the other hand, temephos can act as good suppressor of P450, which may affect its own bioefficacy because it needs to be bioactivated by P450. To the best of our knowledge, this is the first report on acetone-inducible P450 in insects. Further study is needed to characterize the mechanisms involved in acetone induction in P450.

Molecular cloning and expression of CYP9A61: a chlorpyrifos-ethyl and lambda-cyhalothrin-inducible cytochrome P450 cDNA from Cydia pomonella

Cytochrome P450 monooxygenases (CYPs or P450s) play paramount roles in detoxification of insecticides in a number of insect pests. However, little is known about the roles of P450s and their responses to insecticide exposure in the codling moth Cydia pomonella (L.), an economically important fruit pest. Here we report the characterization and expression analysis of the first P450 gene, designated as CYP9A61, from this pest. The full-length cDNA sequence of CYP9A61 is 2071 bp long and its open reading frame (ORF) encodes 538 amino acids. Sequence analysis shows that CYP9A61 shares 51%-60% identity with other known CYP9s and contains the highly conserved substrate recognition site SRS1, SRS4 and SRS5. Quantitative real-time PCR showed that CYP9A61 were 67-fold higher in the fifth instar larvae than in the first instar, and more abundant in the silk gland and fat body than other tissues. Exposure of the 3rd instar larvae to 12.5 mg L(-1) of chlorpyrifos-ethyl for 60 h and 0.19 mg L(-1) of lambda-cyhalothrin for 36 h resulted in 2.20- and 3.47-fold induction of CYP9A61, respectively. Exposure of the 3rd instar larvae to these two insecticides also significantly enhanced the total P450 activity. The results suggested that CYP9A61 is an insecticide-detoxifying P450.

A whole transcriptomal linkage analysis of gene co-regulation in insecticide resistant house flies, Musca domestica

Background

Studies suggest that not only is insecticide resistance conferred via multiple gene up-regulation, but it is mediated through the interaction of regulatory factors. However, no regulatory factors in insecticide resistance have yet been identified, and there has been no examination of the regulatory interaction of resistance genes. Our current study generated the first reference transcriptome from the adult house fly and conducted a whole transcriptome analysis for the multiple insecticide resistant strain ALHF (wild-type) and two insecticide susceptible strains: aabys (with morphological recessive markers) and CS (wild type) to gain valuable insights into the gene interaction and complex regulation in insecticide resistance of house flies, Musca domestica.

Results

Over 56 million reads were used to assemble the adult female M. domestica transcriptome reference and 14488 contigs were generated from the de novo transcriptome assembly. A total of 6159 (43%) of the contigs contained coding regions, among which 1316 genes were identified as being co-up-regulated in ALHF in comparison to both aabys and CS. The majority of these up-regulated genes fell within the SCOP categories of metabolism, general, intra-cellular processes, and regulation, and covered three key detailed function categories: redox detailed function category in metabolism, signal transduction and kinases/phosphatases in regulation, and proteases in intra-cellular processes. The redox group contained detoxification gene superfamilies, including cytochrome P450s, glutathione S-transferases, and esterases. The signal transduction and kinases/phosphatases groups contained gene families of rhodopsin-like GPCRs, adenylate and guanylate cyclases, protein kinases and phosphatases. The proteases group contained genes with digestive, catalytic, and proteinase activities. Genetic linkage analysis with house fly lines comparing different autosomal combinations from ALHF revealed that the up-regulation of gene expression in the three key SCOP detailed function categories occurred mainly through the co-regulation of factors among multiple autosomes, especially between autosomes 2 and 5, suggesting that signaling transduction cascades controlled by GPCRs, protein kinase/phosphates and proteases may be involved in the regulation of resistance P450 gene regulation.

Conclusion

Taken together, our findings suggested that not only is insecticide resistance conferred via multi-resistance mechanisms or up-regulated genes, but it is mediated through the trans and/or cis co-regulations of resistance genes.

Challenges and prospects for dengue and malaria control in Thailand, Southeast Asia

Despite significant advances in the search for potential dengue vaccines and new therapeutic schemes for malaria, the control of these diseases remains difficult. In Thailand, malaria incidence is falling whereas that of dengue is rising, with an increase in the proportion of reported severe cases. In the absence of antiviral therapeutic options for acute dengue, appropriate case management reduces mortality. However, the interruption of transmission still relies on vector control measures that are currently insufficient to curtail the cycle of epidemics. Drug resistance in malaria parasites is increasing, compromising malaria control and elimination. Deficiencies in our knowledge of vector biology and vectorial capacity also hinder public health efforts for vector control. Challenges to dengue and malaria control are discussed, and research priorities identified.

Overexpression of multiple detoxification genes in deltamethrin resistant Laodelphax striatellus (Hemiptera: Delphacidae) in China

BACKGROUND

The small brown planthopper (SBPH), Laodelphax striatellus (Fallén), is one of the major rice pests in Asia and has developed resistance to multiple classes of insecticides. Understanding resistance mechanisms is essential to the management of this pest. Biochemical and molecular assays were performed in this study to systematically characterize deltamethrin resistance mechanisms with laboratory-selected resistant and susceptible strains of SBPH.

METHODOLOGY/PRINCIPAL FINDINGS

Deltamethrin resistant strains of SBPH (JH-del) were derived from a field population by continuously selections (up to 30 generations) in the laboratory, while a susceptible strain (JHS) was obtained from the same population by removing insecticide pressure for 30 generations. The role of detoxification enzymes in the resistance was investigated using synergism and enzyme activity assays with strains of different resistant levels. Furthermore, 71 cytochrome P450, 93 esterases and 12 glutathione-S-transferases cDNAs were cloned based on transcriptome data of a field collected population. Semi-quantitative RT-PCR screening analysis of 176 identified detoxification genes demonstrated that multiple P450 and esterase genes were overexpressed (>2-fold) in JH-del strains (G4 and G30) when compared to that in JHS, and the results of quantitative PCR coincided with the semi-quantitative RT-PCR results. Target mutation at IIS3-IIS6 regions encoded by the voltage-gated sodium channel gene was ruled out for conferring the observed resistance.

CONCLUSION/SIGNIFICANCE

As the first attempt to discover genes potentially involved in SBPH pyrethroid resistance, this study putatively identified several candidate genes of detoxification enzymes that were significantly overexpressed in the resistant strain, which matched the synergism and enzyme activity testing. The biochemical and molecular evidences suggest that the high level pyrethroid resistance in L. striatellus could be due to enhanced detoxification rather than target insensitivity. The findings lay a solid ground for further resistance mechanism elucidation studies.

Characterization and expression of cytochrome p450 cDNA (CYP9AT2) in Chironomus riparius fourth instar larvae exposed to multiple xenobiotics

We identified and characterized a CYP9 family gene, CrCYP9AT2, from Chironomus riparius, an eco-toxicologically important model organism. The 1978 base pair (bp) length CrCYP9AT2 cDNA has an open reading frame of 1587 bp encoding a putative 528 amino acid protein. There was 267 bp 5' and 123 bp 3' untranslated region with a polyadenylation signal site (AATAAA). The putative heme-binding cysteine at position 471 and the typical p450 signature sequence of 463-FGIGPRNCIG-473 were also present. The CrCYP9AT2 transcript was present in all life stages with the highest expression in larvae. The modulation of CrCYP9AT2 was studied using real-time polymerase chain reaction after 24h exposure to cadmium chloride, benzo(a)pyrene; bisphenol A; nonylphenol; chlorpyrifos and ethinylestradiol. Significant up-regulation of CrCYP9AT2 gene was observed after exposure to Cd, B(a)P and CP. However, CrCYP9AT2 was significantly down-regulated after exposure to BPA, NP and EE.

The central role of mosquito cytochrome P450 CYP6Zs in insecticide detoxification revealed by functional expression and structural modelling

The resistance of mosquitoes to chemical insecticides is threatening vector control programmes worldwide. Cytochrome P450 monooxygenases (CYPs) are known to play a major role in insecticide resistance, allowing resistant insects to metabolize insecticides at a higher rate. Among them, members of the mosquito CYP6Z subfamily, like Aedes aegypti CYP6Z8 and its Anopheles gambiae orthologue CYP6Z2, have been frequently associated with pyrethroid resistance. However, their role in the pyrethroid degradation pathway remains unclear. In the present study, we created a genetically modified yeast strain overexpressing Ae. aegypti cytochrome P450 reductase and CYP6Z8, thereby producing the first mosquito P450-CPR (NADPH-cytochrome P450-reductase) complex in a yeast recombinant system. The results of the present study show that: (i) CYP6Z8 metabolizes PBAlc (3-phenoxybenzoic alcohol) and PBAld (3-phenoxybenzaldehyde), common pyrethroid metabolites produced by carboxylesterases, producing PBA (3-phenoxybenzoic acid); (ii) CYP6Z8 transcription is induced by PBAlc, PBAld and PBA; (iii) An. gambiae CYP6Z2 metabolizes PBAlc and PBAld in the same way; (iv) PBA is the major metabolite produced in vivo and is excreted without further modification; and (v) in silico modelling of substrate-enzyme interactions supports a similar role of other mosquito CYP6Zs in pyrethroid degradation. By playing a pivotal role in the degradation of pyrethroid insecticides, mosquito CYP6Zs thus represent good targets for mosquito-resistance management strategies.

Cytochrome P450-encoding genes from the Heliconius genome as candidates for cyanogenesis

Cytochrome P450s are important both in the metabolism of xenobiotics and the production of compounds such as cyanogenic glucosides, which insects use in their defence. In the present study, we use transcriptomic and genomic information to isolate and name P450-encoding genes from the butterfly Heliconius melpomene. We classify each of the putative genes into its appropriate superfamily and compare the distribution of P450s across sequenced insects. We also identify homologues of two P450s known to be involved in cyanogenesis in the six-spot Burnet moth, Zygaena filipendulae. Classification of Heliconius P450s should be an important step in the dissection of their role in the exploitation of their host plant, the passion vine Passiflora.

Zero effect of Bt rice on expression of genes coding for digestion, detoxification and immune responses and developmental performances of brown planthopper Nilaparvata lugens (Stål)

Transgenic Cry1Ac, Cry2Aa and Cry1Ca (Bt toxins) rice lines are well developed to manage lepidopteron pests in China. The impact of transgenic Bt rice on the non-target Brown Planthopper (BPH) has become an essential part of environmental risk assessment, however, scanty evidence is found addressing on developmental and molecular responses of BPH to the ingestion of Bt protein from transgenic rice. The focus of the current study is to examine the developmental characteristics and the expression profiles of gene in relation to digestion, detoxification and immune responses were examined. Our study strongly revealed that the tested Bt rice strains have no unfavorable effect on fecundity, survival and growth of BPH. Furthermore, each of the tested genes did not exhibit distinct expression pattern responding to non Bt parental cultivar, thus, it could be concluded that Bt rice have no detrimental effects on the physiological processes of digestion, detoxification and immune responses of BPH.

RNA interference revealed the roles of two carboxylesterase genes in insecticide detoxification in Locusta migratoria

Carboxylesterases (CarEs) play key roles in metabolism of specific hormones and detoxification of dietary and environmental xenobiotics in insects. We sequenced and characterized CarE cDNAs putatively derived from two different genes named LmCesA1 and LmCesA2 from the migratory locust, Locusta migratoria, one of the most important agricultural pests in the world. The full-length cDNAs of LmCesA1 (1892 bp) and LmCesA2 (1643 bp) encode 543 and 501 amino acid residues, respectively. The two deduced CarEs share a characteristic α/β-hydrolase structure, including a catalytic triad composed of Ser-Glu (Asp)-His and a consensus sequence GQSAG, which suggests that both CarEs are biologically active. Phylogenetic analysis grouped both LmCesA1 and LmCesA2 into clade A which has been suggested to be involved in dietary detoxification. Both transcripts were highly expressed in all the nymphal and adult stages, but only slightly expressed in eggs. Analyses of tissue-dependent expression and in situ hybridization revealed that both transcripts were primarily expressed in gastric caeca. RNA interference (RNAi) of LmCesA1 and LmCesA2 followed by a topical application of carbaryl or deltamethrin did not lead to a significantly increased mortality with either insecticide. However, RNAi of LmCesA1 and LmCesA2 increased insect mortalities by 20.9% and 14.5%, respectively, when chlorpyrifos was applied. These results suggest that these genes might not play a significant role in detoxification of carbaryl and deltamethrin but are most likely to be involved in detoxification of chlorpyrifos in L. migratoria.

Permethrin induction of multiple cytochrome P450 genes in insecticide resistant mosquitoes, Culex quinquefasciatus

The expression of some insect P450 genes can be induced by both exogenous and endogenous compounds and there is evidence to suggest that multiple constitutively overexpressed P450 genes are co-responsible for the development of resistance to permethrin in resistant mosquitoes. This study characterized the permethrin induction profiles of P450 genes known to be constitutively overexpressed in resistant mosquitoes, Culex quinquefasciatus. The gene expression in 7 of the 19 P450 genes CYP325K3v1, CYP4D42v2, CYP9J45, (CYP) CPIJ000926, CYP325G4, CYP4C38, CYP4H40 in the HAmCqG8 strain, increased more than 2-fold after exposure to permethrin at an LC50 concentration (10 ppm) compared to their acetone treated counterpart; no significant differences in the expression of these P450 genes in susceptible S-Lab mosquitoes were observed after permethrin treatment. Eleven of the fourteen P450 genes overexpressed in the MAmCqG6 strain, CYP9M10, CYP6Z12, CYP9J33, CYP9J43, CYP9J34, CYP306A1, CYP6Z15, CYP9J45, CYPPAL1, CYP4C52v1, CYP9J39, were also induced more than doubled after exposure to an LC50 (0.7 ppm) dose of permethrin. No significant induction in P450 gene expression was observed in the susceptible S-Lab mosquitoes after permethrin treatment except for CYP6Z15 and CYP9J39, suggesting that permethrin induction of these two P450 genes are common to both susceptible and resistant mosquitoes while the induction of the others are specific to insecticide resistant mosquitoes. These results demonstrate that multiple P450 genes are co-up-regulated in insecticide resistant mosquitoes through both constitutive overexpression and induction mechanisms, providing additional support for their involvement in the detoxification of insecticides and the development of insecticide resistance.

Insecticide resistance and, efficacy of space spraying and larviciding in the control of dengue vectors Aedes aegypti and Aedes albopictus in Sri Lanka

Unprecedented incidence of dengue has been recorded in Sri Lanka in recent times. Source reduction and use of insecticides in space spraying/fogging and larviciding, are the primary means of controlling the vector mosquitoes Aedes aegypti and Ae. albopictus in the island nation. A study was carried out to understand insecticide cross-resistance spectra and mechanisms of insecticide resistance of both these vectors from six administrative districts, i.e. Kandy, Kurunegala, Puttalam, Gampaha, Ratnapura and Jaffna, of Sri Lanka. Efficacy of the recommended dosages of frequently used insecticides in space spraying and larviciding in dengue vector control programmes was also tested. Insecticide bioassay results revealed that, in general, both mosquito species were highly resistant to DDT but susceptible to propoxur and malathion except Jaffna Ae. aegypti population. Moderate resistance to malathion shown by Jaffna Ae. aegypti population correlated with esterase and malathion carboxylesterase activities of the population. High levels of acetylcholinesterase (AChE) insensitivity in the absence of malathion and propoxur resistance may be due to non-synaptic forms of AChE proteins. Moderate pyrethroid resistance in the absence of high monooxygenase levels indicated the possible involvement of 'kdr' type resistance mechanism in Sri Lankan dengue vectors. Results of the space spraying experiments revealed that 100% mortality at a 10 m distance and >50% mortality at a 50 m distance can be achieved with malathion, pesguard and deltacide even in a ground with dense vegetation. Pesguard and deltacide spraying gave 100% mortality up to 50 m distance in open area and areas with little vegetation. Both species gave >50% mortalities for deltacide at a distance of 75 m in a dense vegetation area. Larval bioassays conducted in the laboratory showed that a 1 ppm temephos solution can maintain a larval mortality rate of 100% for ten months, and the mortality rate declined to 0% in the eleventh month. In the field, where 1 ppm concentration is gradually decreased with water usage, 100% mortality was observed only for the first four months, <50% mortality for the next two months, and 0% mortality was observed eight months after the application of temephos. Deltacide can be effectively used for space spraying programmes in Sri Lanka. Larval control can be successfully achieved through temephos with public participation.

Molecular cloning and expression in vitro of a carboxylesterase gene from the Glanville fritillary butterfly (Melitaea cinxia)

Carboxylesterase (EC 3.1.1.1) is a member of the carboxyl/cholinesterase (CCE) superfamily, which is widely distributed in animals, plants and microorganisms. This enzyme has been known to be associated with insecticide resistance and detoxification. Although CCEs have been extensively studied in insects, including lepidopterans, the research on butterflies, a major subgroup in Lepidoptera, is still poor. In the present study, we cloned a CCE gene (McCCE1) from the Glanville fritillary butterfly (Melitaea cinxia, Lepidoptera: Nymphalidae). The full-length cDNA encoding McCCE1 was 1786 bp, containing a 1641 bp open reading frame encoding 546 amino acids, a 38 bp 5'-untranslated region (5'-UTR), and a 107 bp 3'-UTR with a poly(A) tail. The functionally conserved amino acids in McCCE1 shared the 55% identity with the cytoplasmic esterase CCE017a in Helicoverpa armigera (Lepidoptera: Noctuidae), which has been associated with detoxification. Assays in vitro showed that the recombinant McCCE1 could hydrolyze α- and β-naphthyl acetate. Thus, the present study adds to the body of knowledge concerning the detoxification of pesticides by lepidopterans.

Evidence for gene duplication in the voltage-gated sodium channel gene of Aedes aegypti

BACKGROUND AND OBJECTIVES

Mutations in the voltage-gated sodium channel gene (NaV), known as kdr mutations, are associated with pyrethroid and DDT insecticide resistance in a number of species. In the mosquito dengue vector Aedes aegypti, besides kdr, other polymorphisms allowed grouping AaNaV sequences as type 'A' or 'B'. Here, we point a series of evidences that these polymorphisms are actually involved in a gene duplication event.

METHODOLOGY

Four series of methods were employed: (i) genotypying, with allele-specific PCR (AS-PCR), of two AaNaV sites that can harbor kdr mutations (Ile1011Met and Val1016Ile), (ii) cloning and sequencing of part of the AaNaV gene, (iii) crosses with specific lineages and analysis of the offspring genotypes and (iv) copy number variation assays, with TaqMan quantitative real-time PCR.

RESULTS

kdr mutations in 1011 and 1016 sites were present only in type 'A' sequences, but never in the same haplotype. In addition, although the 1011Met-mutant allele is widely disseminated, no homozygous (1011Met/Met) was detected. Sequencing revealed three distinct haplotypes in some individuals, raising the hypothesis of gene duplication, which was supported by the genotype frequencies in the offspring of specific crosses. Furthermore, it was estimated that a laboratory strain selected for insecticide resistance had 5-fold more copies of the sodium channel gene compared with a susceptible reference strain.

CONCLUSIONS AND IMPLICATIONS

The AaNaV duplication here found might be a recent adaptive response to the intense use of insecticides, maintaining together wild-type and mutant alleles in the same organism, conferring resistance and reducing some of its deleterious effects.

Analysis of the Olive Fruit Fly Bactrocera oleae Transcriptome and Phylogenetic Classification of the Major Detoxification Gene Families

The olive fruit fly Bactrocera oleae has a unique ability to cope with olive flesh, and is the most destructive pest of olives worldwide. Its control has been largely based on the use of chemical insecticides, however, the selection of insecticide resistance against several insecticides has evolved. The study of detoxification mechanisms, which allow the olive fruit fly to defend against insecticides, and/or phytotoxins possibly present in the mesocarp, has been hampered by the lack of genomic information in this species. In the NCBI database less than 1,000 nucleotide sequences have been deposited, with less than 10 detoxification gene homologues in total. We used 454 pyrosequencing to produce, for the first time, a large transcriptome dataset for B. oleae. A total of 482,790 reads were assembled into 14,204 contigs. More than 60% of those contigs (8,630) were larger than 500 base pairs, and almost half of them matched with genes of the order of the Diptera. Analysis of the Gene Ontology (GO) distribution of unique contigs, suggests that, compared to other insects, the assembly is broadly representative for the B. oleae transcriptome. Furthermore, the transcriptome was found to contain 55 P450, 43 GST-, 15 CCE- and 18 ABC transporter-genes. Several of those detoxification genes, may putatively be involved in the ability of the olive fruit fly to deal with xenobiotics, such as plant phytotoxins and insecticides. In summary, our study has generated new data and genomic resources, which will substantially facilitate molecular studies in B. oleae, including elucidation of detoxification mechanisms of xenobiotic, as well as other important aspects of olive fruit fly biology.

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.

Genomic insights into the glutathione S-transferase gene family of two rice planthoppers, Nilaparvata lugens (Stål) and Sogatella furcifera (Horváth) (Hemiptera: Delphacidae)

BACKGROUND

Glutathione S-transferase (GST) genes control crucial traits for the metabolism of various toxins encountered by insects in host plants and the wider environment, including insecticides. The planthoppers Nilaparvata lugens and Sogatella furcifera are serious specialist pests of rice throughout eastern Asia. Their capacity to rapidly adapt to resistant rice varieties and to develop resistance to various insecticides has led to severe outbreaks over the last decade.

METHODOLOGY/PRINCIPAL FINDINGS

Using the genome sequence of N. lugens, we identified for the first time the complete GST gene family of a delphacid insect whilst nine GST gene orthologs were identified from the closely related species S. furcifera. Nilaparvata lugens has 11 GST genes belonging to six cytosolic subclasses and a microsomal class, many fewer than seen in other insects with known genomes. Sigma is the largest GST subclass, and the intron-exon pattern deviates significantly from that of other species. Higher GST gene expression in the N. lugens adult migratory form reflects the higher risk of this life stage in encountering the toxins of non-host plants. After exposure to a sub-lethal dose of four insecticides, chlorpyrifos, imidacloprid, buprofezin or beta-cypermethrin, more GST genes were upregulated in S. furcifera than in N. lugens. RNA interference targeting two N. lugens GST genes, NlGSTe1 and NlGSTm2, significantly increased the sensitivity of fourth instar nymphs to chlorpyrifos but not to beta-cypermethrin.

CONCLUSIONS/SIGNIFICANCE

This study provides the first elucidation of the nature of the GST gene family in a delphacid species, offering new insights into the evolution of metabolic enzyme genes in insects. Further, the use of RNA interference to identify the GST genes induced by insecticides illustrates likely mechanisms for the tolerance of these insects.

Molecular mechanisms associated with increased tolerance to the neonicotinoid insecticide imidacloprid in the dengue vector Aedes aegypti

Mosquitoes are vectors of several major human diseases and their control is mainly based on the use of chemical insecticides. Resistance of mosquitoes to organochlorines, organophosphates, carbamates and pyrethroids led to a regain of interest for the use of neonicotinoid insecticides in vector control. The present study investigated the molecular basis of neonicotinoid resistance in the mosquito Aedes aegypti. A strain susceptible to insecticides was selected at the larval stage with imidacloprid. After eight generations of selection, larvae of the selected strain (Imida-R) showed a 5.4-fold increased tolerance to imidacloprid while adult tolerance level remained low. Imida-R larvae showed significant cross-tolerance to other neonicotinoids but not to pyrethroids, organophosphates and carbamates. Transcriptome profiling identified 344 and 108 genes differentially transcribed in larvae and adults of the Imida-R strain compared to the parental strain. Most of these genes encode detoxification enzymes, cuticle proteins, hexamerins as well as other proteins involved in cell metabolism. Among detoxification enzymes, cytochrome P450 monooxygenases (CYPs) and glucosyl/glucuronosyl transferases (UDPGTs) were over-represented. Bioassays with enzyme inhibitors and biochemical assays confirmed the contribution of P450s with an increased capacity of the Imida-R microsomes to metabolize imidacloprid in presence of NADPH. Comparison of substrate recognition sites and imidacloprid docking models of six CYP6s over-transcribed in the Imida-R strain together with Bemisia tabaci CYP6CM1vQ and Drosophila melanogaster CYP6G1, both able to metabolize imidacloprid, suggested that CYP6BB2 and CYP6N12 are good candidates for imidacloprid metabolism in Ae. aegypti. The present study revealed that imidacloprid tolerance in mosquitoes can arise after few generations of selection at the larval stage but does not lead to a significant tolerance of adults. As in other insects, P450-mediated insecticide metabolism appears to play a major role in imidacloprid tolerance in mosquitoes.

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

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

Pyrethroid resistance in Aedes aegypti and Aedes albopictus from Port-au-Prince, Haiti

In Port-au-Prince, Haiti, the status of insecticide resistance has not recently been evaluated for Aedes aegypti (L) and Aedes albopictus (Skuse) populations. No prophylactics exist for dengue, so prevention is only through vector control methods. An earthquake occurred in Haiti on January 12, 2010, with a magnitude of 7.0 Mw that devastated the area. Dengue became a major concern for the humanitarian relief workers that entered the country. Bottle bioassays were conducted in the field on adult mosquitoes reared from larvae collected from the grounds of the U.S. Embassy and from an adjacent neighborhood in eastern Port-au-Prince, Haiti. At the CDC, Fort Collins, CO, bioassays, molecular, and biochemical assays were performed on mosquitoes reared from field-collected eggs. A small percentage of the population was able to survive the diagnostic dose in bioassays run in Haiti. Mosquitoes tested at the CDC demonstrated no phenotypic resistance. A variety of factors could be responsible for the discrepancies between the field and lab data, but temperature and larval nutrition are probably most important. Knowledge of localized resistance and underlying mechanisms helps in making rational decisions in selection of appropriate and effective insecticides in the event of a dengue outbreak.

Using Drosophila melanogaster to validate metabolism-based insecticide resistance from insect pests

Identifying molecular mechanisms of insecticide resistance is important for preserving insecticide efficacy, developing new insecticides and implementing insect control. The metabolic detoxification of insecticides is a widespread resistance mechanism. Enzymes with the potential to detoxify insecticides are commonly encoded by members of the large cytochrome P450, glutathione S-transferase and carboxylesterase gene families, all rapidly evolving in insects. Here, we demonstrate that the model insect Drosophila melanogaster is useful for functionally validating the role of metabolic enzymes in conferring metabolism-based insecticide resistance. Alleles of three well-characterized genes from different pest insects were expressed in transgenic D. melanogaster : a carboxylesterase gene (αE7) from the Australian sheep blowfly Lucilia cuprina, a glutathione S-transferase gene (GstE2) from the mosquito Anopheles gambiae and a cytochrome P450 gene (Cyp6cm1) from the whitefly Bemisia tabaci. For all genes, expression in D. melanogaster resulted in insecticide resistance phenotypes mirroring those observed in resistant populations of the pest species. Using D. melanogaster to assess the potential for novel metabolic resistance mechanisms to evolve in pest species is discussed.

Transcriptomic and phylogenetic analysis of Culex pipiens quinquefasciatus for three detoxification gene families

Background

The genomes of three major mosquito vectors of human diseases, Anopheles gambiae, Aedes aegypti, and Culex pipiens quinquefasciatus, have been previously sequenced. C. p. quinquefasciatus has the largest number of predicted protein-coding genes, which partially results from the expansion of three detoxification gene families: cytochrome P450 monooxygenases (P450), glutathione S-transferases (GST), and carboxyl/cholinesterases (CCE). However, unlike An. gambiae and Ae. aegypti, which have large amounts of gene expression data, C. p. quinquefasciatus has limited transcriptomic resources. Knowledge of complete gene expression information is very important for the exploration of the functions of genes involved in specific biological processes. In the present study, the three detoxification gene families of C. p. quinquefasciatus were analyzed for phylogenetic classification and compared with those of three other dipteran insects. Gene expression during various developmental stages and the differential expression responsible for parathion resistance were profiled using the digital gene expression (DGE) technique.

Results

A total of 302 detoxification genes were found in C. p. quinquefasciatus, including 71 CCE, 196 P450, and 35 cytosolic GST genes. Compared with three other dipteran species, gene expansion in Culex mainly occurred in the CCE and P450 families, where the genes of α-esterases, juvenile hormone esterases, and CYP325 of the CYP4 subfamily showed the most pronounced expansion on the genome. For the five DGE libraries, 3.5-3.8 million raw tags were generated and mapped to 13314 reference genes. Among 302 detoxification genes, 225 (75%) were detected for expression in at least one DGE library. One fourth of the CCE and P450 genes were detected uniquely in one stage, indicating potential developmentally regulated expression. A total of 1511 genes showed different expression levels between a parathion-resistant and a susceptible strain. Fifteen detoxification genes, including 2 CCEs, 6 GSTs, and 7 P450s, were expressed at higher levels in the resistant strain.

Conclusions

The results of the present study provide new insights into the functions and evolution of three detoxification gene families in mosquitoes and comprehensive transcriptomic resources for C. p. quinquefasciatus, which will facilitate the elucidation of molecular mechanisms underlying the different biological characteristics of the three major mosquito vectors.

Multiple sodium channel variants in the mosquito Culex quinquefasciatus

Voltage-gated sodium channels are the target sites of both DDT and pyrethroid insecticides. The importance of alternative splicing as a key mechanism governing the structural and functional diversity of sodium channels and the resulting development of insecticide and acaricide resistance is widely recognized, as shown by the extensive research on characterizing alternative splicing and variants of sodium channels in medically and agriculturally important insect species. Here we present the first comparative study of multiple variants of the sodium channel transcripts in the mosquito Culex quinquefasciatus. The variants were classified into two categories, CxNa-L and CxNa-S based on their distinguishing sequence sizes of ~6.5 kb and ~4.0 kb, respectively, and generated via major extensive alternative splicing with minor small deletions/ insertions in susceptible S-Lab, low resistant HAmCq(G0), and highly resistant HAmCq(G8)Culex strains. Four alternative Cx-Na-L splice variants were identified, including three full length variants with three optional exons (2, 5, and 21i) and one with in-frame-stop codons. Large, multi-exon-alternative splices were identified in the CxNa-S category. All CxNa-S splicing variants in the S-Lab and HAmCq(G0) strains contained in-frame stop codons, suggesting that any resulting proteins would be truncated. The ~1000 to ~3000-fold lower expression of these splice variants with stop codons compared with the CxNa-L splicing variants may support the lower importance of these variants in S-Lab and HAmCq(G0). Interestingly, two alternative splicing variants of CxNa-S in HAmCq(G8) included entire ORFs but lacked exons 5 to18 and these two variants had much higher expression levels in HAmCq(G8) than in S-Lab and HAmCq(G0). These results provide a functional basis for further characterizing how alternative splicing of a voltage-gated sodium channel contributes to diversity in neuronal signaling in mosquitoes in response to pyrethroids, and possibly indicates the role of these variants in the development of pyrethroid resistance.

The transcriptome profile of the mosquito Culex quinquefasciatus following permethrin selection

To gain valuable insights into the gene interaction and the complex regulation system involved in the development of insecticide resistance in mosquitoes Culex quinquefasciatus, we conducted a whole transcriptome analysis of Culex mosquitoes following permethrin selection. Gene expression profiles for the lower resistant parental mosquito strain HAmCq^G0 and their permethrin-selected high resistant offspring HAmCq^G8 were compared and a total of 367 and 3982 genes were found to be up- and down-regulated, respectively, in HAmCq^G8, indicating that multiple genes are involved in response to permethrin selection. However, a similar overall cumulative gene expression abundance was identified between up- and down-regulated genes in HAmCq^G8 mosquitoes following permethrin selection, suggesting a homeostatic response to insecticides through a balancing of the up- and down-regulation of the genes. While structural and/or cuticular structural functions were the only two enriched GO terms for down-regulated genes, the enriched GO terms obtained for the up-regulated genes occurred primarily among the catalytic and metabolic functions where they represented three functional categories: electron carrier activity, binding, and catalytic activity. Interestingly, the functional GO terms in these three functional categories were overwhelmingly overrepresented in P450s and proteases/serine proteases. The important role played by P450s in the development of insecticide resistance has been extensively studied but the function of proteases/serine proteases in resistance is less well understood. Hence, the characterization of the functions of these proteins, including their digestive, catalytic and proteinase activities; regulation of signaling transduction and protein trafficking, immunity and storage; and their precise function in the development of insecticide resistance in mosquitoes will provide new insights into how genes are interconnected and regulated in resistance.

Identification and expression profiles of nine glutathione S-transferase genes from the important rice phloem sap-sucker and virus vector Laodelphax striatellus (Fallén) (Hemiptera: Delphacidae)

BACKGROUND

Glutathione S-transferases (GSTs) have received considerable attention in insects for their roles in insecticide resistance. Laodelphax striatellus (Fallén) is a serious rice pest. L. striatellus outbreaks occur frequently throughout eastern Asia. A key problem in controlling this pest is its rapid adaptation to numerous insecticides. In this research, nine cDNAs encoding GSTs in L. striatellus were cloned and characterised.

RESULTS

The cloned GSTs of L. striatellus belonged to six cytosolic classes and a microsomal subgroup. Exposure to sublethal concentrations of each of the six insecticides, DDT, chlorpyrifos, fipronil, imidacloprid, buprofezin and beta-cypermethrin, quickly induced (6 h) up-expression of LsGSTe1. The expression of LsGSTs2 was increased by chlorpyrifos, fipronil and beta-cypermethrin. Furthermore, exposure of L. striatellus to fipronil, imidacloprid, buprofezin and beta-cypermethrin increased the expression of the LsGSTm gene after 24 or 48 h.

CONCLUSION

This work is the first identification of GST genes from different GST groups in Auchenorrhyncha species and their induction characteristics with insecticide types and time. The elevated expression of GST genes induced by insecticides might be related to the enhanced tolerance of this insect to insecticides and xenobiotics.

Expression profile of genes during resistance reversal in a temephos selected strain of the dengue vector, Aedes aegypti

Background

The mosquito Aedes aegypti is one of the most important disease vectors because it transmits two major arboviruses, dengue and yellow fever, which cause significant global morbidity and mortality. Chemical insecticides form the cornerstone of vector control. The organophosphate temephos a larvicide recommended by WHO for controlling Ae. aegypti, however, resistance to this compound has been reported in many countries, including Brazil.

Methodology/Principal Findings

The aim of this study was to identify genes implicated in metabolic resistance in an Ae. aegypti temephos resistant strain, named RecR, through microarray analysis. We utilized a custom ‘Ae. aegypti detox chip’ and validated microarray data through RT-PCR comparing susceptible and resistant individuals. In addition, we analyzed gene expression in 4^th instar larvae from a reversed susceptible strain (RecRev), exposed and unexposed to temephos. The results obtained revealed a set of 13 and 6 genes significantly over expressed in resistant adult mosquitoes and larvae, respectively. One of these genes, the cytochrome P450 CYP6N12, was up-regulated in both stages. RT-PCR confirmed the microarray results and, additionally, showed no difference in gene expression between temephos exposed and unexposed RecRev mosquitoes. This suggested that the differences in the transcript profiles among the strains are heritable due to a selection process and are not caused by immediate insecticide exposure. Reversal of temephos resistance was demonstrated and, importantly, there was a positive correlation between a decrease in the resistance ratio and an accompanying decrease in the expression levels of previously over expressed genes. Some of the genes identified here have also been implicated in metabolic resistance in other mosquito species and insecticide resistant populations of Ae. aegypti.

Conclusions/Significance

The identification of gene expression signatures associated to insecticide resistance and their suppression could greatly aid the development of improved strategies of vector control.

Gene amplification, ABC transporters and cytochrome P450s: unraveling the molecular basis of pyrethroid resistance in the dengue vector, Aedes aegypti

Background

Pyrethroid insecticides are widely utilized in dengue control. However, the major vector, Aedes aegypti, is becoming increasingly resistant to these insecticides and this is impacting on the efficacy of control measures. The near complete transcriptome of two pyrethroid resistant populations from the Caribbean was examined to explore the molecular basis of this resistance.

Principal Findings

Two previously described target site mutations, 1016I and 1534C were detected in pyrethroid resistant populations from Grand Cayman and Cuba. In addition between two and five per cent of the Ae. aegypti transcriptome was differentially expressed in the resistant populations compared to a laboratory susceptible population. Approximately 20 per cent of the genes over-expressed in resistant mosquitoes were up-regulated in both Caribbean populations (107 genes). Genes with putative monooxygenase activity were significantly over represented in the up-regulated subset, including five CYP9 P450 genes. Quantitative PCR was used to confirm the higher transcript levels of multiple cytochrome P450 genes from the CYP9J family and an ATP binding cassette transporter. Over expression of two genes, CYP9J26 and ABCB4, is due, at least in part, to gene amplification.

Significance

These results, and those from other studies, strongly suggest that increases in the amount of the CYP9J cytochrome P450s are an important mechanism of pyrethroid resistance in Ae. aegypti. The genetic redundancy resulting from the expansion of this gene family makes it unlikely that a single gene or mutation responsible for pyrethroid resistance will be identified in this mosquito species. However, the results from this study do pave the way for the development of new pyrethroid synergists and improved resistance diagnostics. The role of copy number polymorphisms in detoxification and transporter genes in providing protection against insecticide exposure requires further investigation.

Dengue is the most rapidly spreading arboviral infection of humans and each year there are 50–100 million cases of dengue fever. There is no vaccine or drug to prevent dengue infection so control of the mosquitoes that transmit this virus is the only option to reduce transmission. Removing mosquito habitats close to human homes can be effective but in reality most dengue control programmes rely on a small number of chemical insecticides. Therefore, when the mosquito vectors develop resistance to the available insecticides, dengue control is jeopardized. In this study we examined the causes of resistance to the insecticide class most commonly used in mosquito control, the pyrethroids. We found that a group of genes, which have been implicated in detoxifying these insecticides in other populations of dengue vectors, were highly over expressed in both these Caribbean populations and we investigated the molecular basis of this increased expression. The next steps, which will be a considerable challenge, are to utilize this information to develop effective means of restoring insecticide susceptibility in dengue vectors.

Detoxification enzymes associated with insecticide resistance in laboratory strains of Anopheles arabiensis of different geographic origin

Background

The use of insecticides to control malaria vectors is essential to reduce the prevalence of malaria and as a result, the development of insecticide resistance in vector populations is of major concern. Anopheles arabiensis is one of the main African malaria vectors and insecticide resistance in this species has been reported in a number of countries. The aim of this study was to investigate the detoxification enzymes that are involved in An. arabiensis resistance to DDT and pyrethroids.

Methods

The detoxification enzyme profiles were compared between two DDT selected, insecticide resistant strains of An. arabiensis, one from South Africa and one from Sudan, using the An. gambiae detoxification chip, a boutique microarray based on the major classes of enzymes associated with metabolism and detoxification of insecticides. Synergist assays were performed in order to clarify the roles of over-transcribed detoxification genes in the observed resistance phenotypes. In addition, the presence of kdr mutations in the colonies under investigation was determined.

Results

The microarray data identifies several genes over-transcribed in the insecticide selected South African strain, while in the Sudanese population, only one gene, CYP9L1, was found to be over-transcribed. The outcome of the synergist experiments indicate that the over-transcription of detoxification enzymes is linked to deltamethrin resistance, while DDT and permethrin resistance are mainly associated with the presence of the L1014F kdr mutation.

Conclusions

These data emphasise the complexity associated with resistance phenotypes and suggest that specific insecticide resistance mechanisms cannot be extrapolated to different vector populations of the same species.

Insecticide resistance mechanisms in the green peach aphid Myzus persicae (Hemiptera: Aphididae) I: A transcriptomic survey

BACKGROUND

Insecticide resistance is one of the best examples of rapid micro-evolution found in nature. Since the development of the first synthetic insecticide in 1939, humans have invested considerable effort to stay ahead of resistance phenotypes that repeatedly develop in insects. Aphids are a group of insects that have become global pests in agriculture and frequently exhibit insecticide resistance. The green peach aphid, Myzus persicae, has developed resistance to at least seventy different synthetic compounds, and different insecticide resistance mechanisms have been reported worldwide.

METHODOLOGY/PRINCIPAL FINDINGS

To further characterize this resistance, we analyzed genome-wide transcriptional responses in three genotypes of M. persicae, each exhibiting different resistance mechanisms, in response to an anti-cholinesterase insecticide. The sensitive genotype (exhibiting no resistance mechanism) responded to the insecticide by up-regulating 183 genes primarily ones related to energy metabolism, detoxifying enzymes, proteins of extracellular transport, peptidases and cuticular proteins. The second genotype (resistant through a kdr sodium channel mutation), up-regulated 17 genes coding for detoxifying enzymes, peptidase and cuticular proteins. Finally, a multiply resistant genotype (carrying kdr and a modified acetylcholinesterase), up-regulated only 7 genes, appears not to require induced insecticide detoxification, and instead down-regulated many genes.

CONCLUSIONS/SIGNIFICANCE

This study suggests strongly that insecticide resistance in M. persicae is more complex that has been described, with the participation of a broad array of resistance mechanisms. The sensitive genotype exhibited the highest transcriptional plasticity, accounting for the wide range of potential adaptations to insecticides that this species can evolve. In contrast, the multiply resistant genotype exhibited a low transcriptional plasticity, even for the expression of genes encoding enzymes involved in insecticide detoxification. Our results emphasize the value of microarray studies to search for regulated genes in insects, but also highlights the many ways those different genotypes can assemble resistant phenotypes depending on the environmental pressure.

Novel cytochrome P450 genes, CYP6EB1 and CYP6EC1, are over-expressed in acrinathrin-resistant Frankliniella occidentalis (Thysanoptera: Thripidae)

Control of Frankliniella occidentalis (Pergande) is a serious problem for agriculture all over the world because of the limited range of insecticides that are available. Insecticide resistance in F. occidentalis has been reported for all major insecticide groups. Our previous studies showed that cytochrome P450-mediated detoxification is a major mechanism responsible for insecticide resistance in this pest. Degenerate polymerase chain reaction was used to identify P450 genes that might be involved in acrinathrin resistance, in a laboratory population of F. occidentalis. Associated sequences were classified as belonging to the CYP4 and CYP6 families. Real-time quantitative polymerase chain reaction analyses revealed that two genes, CYP6EB1 and CYP6EC1, were over-expressed in adults and L2 larvae of the resistant population, when compared with the susceptible population, suggesting their possible involvement in resistance to acrinathrin.

Phylogenomics of the benzoxazinoid biosynthetic pathway of Poaceae: gene duplications and origin of the Bx cluster

Background

The benzoxazinoids 2,4-dihydroxy-1,4-benzoxazin-3-one (DIBOA) and 2,4-dihydroxy-7- methoxy-1,4-benzoxazin-3-one (DIMBOA), are key defense compounds present in major agricultural crops such as maize and wheat. Their biosynthesis involves nine enzymes thought to form a linear pathway leading to the storage of DI(M)BOA as glucoside conjugates. Seven of the genes (Bx1-Bx6 and Bx8) form a cluster at the tip of the short arm of maize chromosome 4 that includes four P450 genes (Bx2-5) belonging to the same CYP71C subfamily. The origin of this cluster is unknown.

Results

We show that the pathway appeared following several duplications of the TSA gene (α-subunit of tryptophan synthase) and of a Bx2-like ancestral CYP71C gene and the recruitment of Bx8 before the radiation of Poaceae. The origins of Bx6 and Bx7 remain unclear. We demonstrate that the Bx2-like CYP71C ancestor was not committed to the benzoxazinoid pathway and that after duplications the Bx2-Bx5 genes were under positive selection on a few sites and underwent functional divergence, leading to the current specific biochemical properties of the enzymes. The absence of synteny between available Poaceae genomes involving the Bx gene regions is in contrast with the conserved synteny in the TSA gene region.

Conclusions

These results demonstrate that rearrangements following duplications of an IGL/TSA gene and of a CYP71C gene probably resulted in the clustering of the new copies (Bx1 and Bx2) at the tip of a chromosome in an ancestor of grasses. Clustering favored cosegregation and tip chromosomal location favored gene rearrangements that allowed the further recruitment of genes to the pathway. These events, a founding event and elongation events, may have been the key to the subsequent evolution of the benzoxazinoid biosynthetic cluster.

A cytochrome p450 conserved in insects is involved in cuticle formation

The sequencing of numerous insect genomes has revealed dynamic changes in the number and identity of cytochrome P450 genes in different insects. In the evolutionary sense, the rapid birth and death of many P450 genes is observed, with only a small number of P450 genes showing orthology between insects with sequenced genomes. It is likely that these conserved P450s function in conserved pathways. In this study, we demonstrate the P450 gene, Cyp301a1, present in all insect genomes sequenced to date, affects the formation of the adult cuticle in Drosophila melanogaster. A Cyp301a1 piggyBac insertion mutant and RNAi of Cyp301a1 both show a similar cuticle malformation phenotype, which can be reduced by 20-hydroxyecdysone, suggesting that Cyp301a1 is an important gene involved in the formation of the adult cuticle and may be involved in ecdysone regulation in this tissue.

Pyrosequencing the Bemisia tabaci transcriptome reveals a highly diverse bacterial community and a robust system for insecticide resistance

Background

Bemisia tabaci (Gennadius) is a phloem-feeding insect poised to become one of the major insect pests in open field and greenhouse production systems throughout the world. The high level of resistance to insecticides is a main factor that hinders continued use of insecticides for suppression of B. tabaci. Despite its prevalence, little is known about B. tabaci at the genome level. To fill this gap, an invasive B. tabaci B biotype was subjected to pyrosequencing-based transcriptome analysis to identify genes and gene networks putatively involved in various physiological and toxicological processes.

Methodology and Principal Findings

Using Roche 454 pyrosequencing, 857,205 reads containing approximately 340 megabases were obtained from the B. tabaci transcriptome. De novo assembly generated 178,669 unigenes including 30,980 from insects, 17,881 from bacteria, and 129,808 from the nohit. A total of 50,835 (28.45%) unigenes showed similarity to the non-redundant database in GenBank with a cut-off E-value of 10–5. Among them, 40,611 unigenes were assigned to one or more GO terms and 6,917 unigenes were assigned to 288 known pathways. De novo metatranscriptome analysis revealed highly diverse bacterial symbionts in B. tabaci, and demonstrated the host-symbiont cooperation in amino acid production. In-depth transcriptome analysis indentified putative molecular markers, and genes potentially involved in insecticide resistance and nutrient digestion. The utility of this transcriptome was validated by a thiamethoxam resistance study, in which annotated cytochrome P450 genes were significantly overexpressed in the resistant B. tabaci in comparison to its susceptible counterparts.

Conclusions

This transcriptome/metatranscriptome analysis sheds light on the molecular understanding of symbiosis and insecticide resistance in an agriculturally important phloem-feeding insect pest, and lays the foundation for future functional genomics research of the B. tabaci complex. Moreover, current pyrosequencing effort greatly enriched the existing whitefly EST database, and makes RNAseq a viable option for future genomic analysis.

Phenotypic plasticity of the Drosophila transcriptome

Phenotypic plasticity is the ability of a single genotype to produce different phenotypes in response to changing environments. We assessed variation in genome-wide gene expression and four fitness-related phenotypes of an outbred Drosophila melanogaster population under 20 different physiological, social, nutritional, chemical, and physical environments; and we compared the phenotypically plastic transcripts to genetically variable transcripts in a single environment. The environmentally sensitive transcriptome consists of two transcript categories, which comprise ∼15% of expressed transcripts. Class I transcripts are genetically variable and associated with detoxification, metabolism, proteolysis, heat shock proteins, and transcriptional regulation. Class II transcripts have low genetic variance and show sexually dimorphic expression enriched for reproductive functions. Clustering analysis of Class I transcripts reveals a fragmented modular organization and distinct environmentally responsive transcriptional signatures for the four fitness-related traits. Our analysis suggests that a restricted environmentally responsive segment of the transcriptome preserves the balance between phenotypic plasticity and environmental canalization.

Unlike Mendelian traits, where the genotype allows a direct prediction of the phenotype, predicting phenotypic values is not straightforward for complex traits, which arise from multiple segregating genes and their interactions with the environment. Here, a single genotype can often express different phenotypes in different environments. Such phenotypic plasticity is the counterpoint to “environmental canalization,” whereby genotypes produce the same phenotype in different environments. Whereas phenotypic plasticity allows organisms to respond rapidly to changing environments, environmental canalization buffers phenotypes against environmental perturbations. The balance between plasticity and robustness is crucial for optimal fitness, but the genetic basis for phenotypic plasticity is poorly defined. Here, we present the most comprehensive analysis to date of variation in genome-wide gene expression of an outbred Drosophila melanogaster population under 20 different environments. We find that a restricted environmentally responsive segment of the transcriptome (∼15%) preserves the balance between phenotypic plasticity and environmental canalization. Environmentally plastic transcripts can be grouped into two categories. Class I transcripts are genetically variable and associated with detoxification, metabolism, proteolysis, heat shock proteins, and transcriptional regulation. Class II transcripts have low genetic variance and show sexually dimorphic expression enriched for reproductive functions. Despite low genetic variance these transcripts evolve rapidly.

Pinpointing P450s associated with pyrethroid metabolism in the dengue vector, Aedes aegypti: developing new tools to combat insecticide resistance

BACKGROUND

Pyrethroids are increasingly used to block the transmission of diseases spread by Aedes aegypti such as dengue and yellow fever. However, insecticide resistance poses a serious threat, thus there is an urgent need to identify the genes and proteins associated with pyrethroid resistance in order to produce effective counter measures. In Ae. aegypti, overexpression of P450s such as the CYP9J32 gene have been linked with pyrethroid resistance. Our aim was to confirm the role of CYP9J32 and other P450s in insecticide metabolism in order to identify potential diagnostic resistance markers.

METHODOLOGY/PRINCIPAL FINDINGS

We have expressed CYP9J32 in Escherichia coli and show that the enzyme can metabolize the pyrethroids permethrin and deltamethrin. In addition, three other Ae. aegypti P450s (CYP9J24, CYP9J26, CYP9J28) were found capable of pyrethroid metabolism, albeit with lower activity. Both Ae. aegypti and Anopheles gambiae P450s (CYP's 6M2, 6Z2, 6P3) were screened against fluorogenic and luminescent substrates to identify potential diagnostic probes for P450 activity. Luciferin-PPXE was preferentially metabolised by the three major pyrethroid metabolisers (CYP9J32, CYP6M2 and CYP6P3), identifying a potential diagnostic substrate for these P450s.

CONCLUSIONS/SIGNIFICANCE

P450s have been identified with the potential to confer pyrethroid resistance in Ae.aegypti. It is recommended that over expression of these enzymes should be monitored as indicators of resistance where pyrethroids are used.

Insecticide resistance in the dengue vector Aedes aegypti from Martinique: distribution, mechanisms and relations with environmental factors

Dengue is an important mosquito borne viral disease in Martinique Island (French West Indies). The viruses responsible for dengue are transmitted by Aedes aegypti, an indoor day-biting mosquito. The most effective proven method for disease prevention has been by vector control by various chemical or biological means. Unfortunately insecticide resistance has already been observed on the Island and recently showed to significantly reduce the efficacy of vector control interventions. In this study, we investigated the distribution of resistance and the underlying mechanisms in nine Ae. aegypti populations. Statistical multifactorial approach was used to investigate the correlations between insecticide resistance levels, associated mechanisms and environmental factors characterizing the mosquito populations. Bioassays revealed high levels of resistance to temephos and deltamethrin and susceptibility to Bti in the 9 populations tested. Biochemical assays showed elevated detoxification enzyme activities of monooxygenases, carboxylesterases and glutathione S-tranferases in most of the populations. Molecular screening for common insecticide target-site mutations, revealed the presence of the "knock-down resistance" V1016I Kdr mutation at high frequency (>87%). Real time quantitative RT-PCR showed the potential involvement of several candidate detoxification genes in insecticide resistance. Principal Component Analysis (PCA) performed with variables characterizing Ae. aegypti from Martinique permitted to underline potential links existing between resistance distribution and other variables such as agriculture practices, vector control interventions and urbanization. Insecticide resistance is widespread but not homogeneously distributed across Martinique. The influence of environmental and operational factors on the evolution of the resistance and mechanisms are discussed.

De novo intestine-specific transcriptome of the brown planthopper Nilaparvata lugens revealed potential functions in digestion, detoxification and immune response

The brown planthopper (Nilaparvata lugens, BPH) is the most serious rice plant pests in Asia. In this study, we performed transcriptome-wide analysis on BPH intestine. We obtained more than 26 million sequencing reads that were then assembled into 53,553 unigenes with a mean size of 388 bp. Based on similarity search with the nucleotide sequences available at NCBI, BPH intestine-specific transcriptome analysis identified 21,405 sequences. Assembled sequences were annotated with gene description, gene ontology and clusters of orthologous group terms. The digestion-, defense- and xenobiotic metabolism-related genes were abundantly detected in the transcripts from BPH intestine. Many novel genes including 33 digestion-related genes, 25 immune responsive genes and 27 detoxification-related genes are first reported here. We investigated the gene expression patterns at the transcript levels in different tissues by quantitative real-time PCR analysis, which revealed that some genes had intestine-specific expression, implicating their potential significance for BPH management.

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.

Comparative analysis of the UDP-glycosyltransferase multigene family in insects

UDP-glycosyltransferases (UGT) catalyze the conjugation of a range of diverse small lipophilic compounds with sugars to produce glycosides, playing an important role in the detoxification of xenobiotics and in the regulation of endobiotics in insects. Recent progress in genome sequencing has enabled an assessment of the extent of the UGT multigene family in insects. Here we report over 310 putative UGT genes identified from genomic databases of eight different insect species together with a transcript database from the lepidopteran Helicoverpa armigera. Phylogenetic analysis of the insect UGTs showed Order-specific gene diversification and inter-species conservation of this multigene family. Only one family (UGT50) is found in all insect species surveyed (except the pea aphid) and may be homologous to mammalian UGT8. Three families (UGT31, UGT32, and UGT305) related to Lepidopteran UGTs are unique to baculoviruses. A lepidopteran sub-tree constructed with 40 H. armigera UGTs and 44 Bombyx mori UGTs revealed that lineage-specific expansions of some families in both species appear to be driven by diversification in the N-terminal substrate binding domain, increasing the range of compounds that could be detoxified or regulated by glycosylation. By comparison of the deduced protein sequences, several important domains were predicted, including the N-terminal signal peptide, UGT signature motif, and C-terminal transmembrane domain. Furthermore, several conserved residues putatively involved in sugar donor binding and catalytic mechanism were also identified by comparison with human UGTs. Many UGTs were expressed in fat body, midgut, and Malpighian tubules, consistent with functions in detoxification, and some were expressed in antennae, suggesting a role in pheromone deactivation. Transcript variants derived from alternative splicing, exon skipping, or intron retention produced additional UGT diversity. These findings from this comparative study of two lepidopteran UGTs as well as other insects reveal a diversity comparable to this gene family in vertebrates, plants and fungi and show the magnitude of the task ahead, to determine biochemical function and physiological relevance of each UGT enzyme.

RNA-Seq and molecular docking reveal multi-level pesticide resistance in the bed bug

BACKGROUND

Bed bugs (Cimex lectularius) are hematophagous nocturnal parasites of humans that have attained high impact status due to their worldwide resurgence. The sudden and rampant resurgence of C. lectularius has been attributed to numerous factors including frequent international travel, narrower pest management practices, and insecticide resistance.

RESULTS

We performed a next-generation RNA sequencing (RNA-Seq) experiment to find differentially expressed genes between pesticide-resistant (PR) and pesticide-susceptible (PS) strains of C. lectularius. A reference transcriptome database of 51,492 expressed sequence tags (ESTs) was created by combining the databases derived from de novo assembled mRNA-Seq tags (30,404 ESTs) and our previous 454 pyrosequenced database (21,088 ESTs). The two-way GLMseq analysis revealed ~15,000 highly significant differentially expressed ESTs between the PR and PS strains. Among the top 5,000 differentially expressed ESTs, 109 putative defense genes (cuticular proteins, cytochrome P450s, antioxidant genes, ABC transporters, glutathione S-transferases, carboxylesterases and acetyl cholinesterase) involved in penetration resistance and metabolic resistance were identified. Tissue and development-specific expression of P450 CYP3 clan members showed high mRNA levels in the cuticle, Malpighian tubules, and midgut; and in early instar nymphs, respectively. Lastly, molecular modeling and docking of a candidate cytochrome P450 (CYP397A1V2) revealed the flexibility of the deduced protein to metabolize a broad range of insecticide substrates including DDT, deltamethrin, permethrin, and imidacloprid.

CONCLUSIONS

We developed significant molecular resources for C. lectularius putatively involved in metabolic resistance as well as those participating in other modes of insecticide resistance. RNA-Seq profiles of PR strains combined with tissue-specific profiles and molecular docking revealed multi-level insecticide resistance in C. lectularius. Future research that is targeted towards RNA interference (RNAi) on the identified metabolic targets such as cytochrome P450s and cuticular proteins could lay the foundation for a better understanding of the genetic basis of insecticide resistance in C. lectularius.

Molecular characterizations of two cytochrome P450 genes encoding CYP6A41 and CYP6EK1 from the oriental fruit fly, Bactrocera dorsalis (Diptera: Tephritidae)

Two P450 genes encoding CYP6A41 and CYP6EK1 were cloned from the oriental fruit fly using polymerase chain reaction (PCR) and rapid amplification of cDNA ends (RACE) techniques. CYP6A41 and CYP6EK1 contained open reading frames of 1,530 and 1,524 nucleotides that encode 510 and 508 amino acid residues, respectively. The putative proteins shared 44% identity with each other. Phylogenetic analysis showed that CYP6A41 and CYP6EK1 were most closely related to Ceratitis capitata CYP6A10 and CYP6A subfamily. Expression patterns of the two genes in different geographical populations (Yunnan, Hainan, Dongguang, and Guangzhou), developmental stages (eggs, larvae, pupae, and adults), and tissues (midguts, fat bodies, and Malpighian tubules) were analyzed by real-time quantitative PCR (RT-qPCR) methods. The results showed that the expression levels of CYP6EK1 were significantly different among the four populations, but were not different for CYP6A41. Both the expressions of CYP6A41 and CYP6EK1 were development specific and had significantly higher levels in the larval stage. The expression of CYP6A41 did not vary among the midgut, fat body, or Malpighian tubules; however, CYP6EK1 expression was higher in the Malpighian tubules. The results suggest that CYP6A41 and CYP6EK1 might be involved in detoxification of xenobiotic compounds that were harmful to larval flies or development. Moreover, high expression of CYP6EK1 in the Malpighian tubules also implied participation in detoxification.

Genome analysis of cytochrome P450s and their expression profiles in insecticide resistant mosquitoes, Culex quinquefasciatus

Here we report a study of the 204 P450 genes in the whole genome sequence of larvae and adult Culex quinquefasciatus mosquitoes. The expression profiles of the P450 genes were compared for susceptible (S-Lab) and resistant mosquito populations, two different field populations of mosquitoes (HAmCq and MAmCq), and field parental mosquitoes (HAmCq(G0) and MAmCq(G0)) and their permethrin selected offspring (HAmCq(G8) and MAmCq(G6)). While the majority of the P450 genes were expressed at a similar level between the field parental strains and their permethrin selected offspring, an up- or down-regulation feature in the P450 gene expression was observed following permethrin selection. Compared to their parental strains and the susceptible S-Lab strain, HAmCq(G8) and MAmCq(G6) were found to up-regulate 11 and 6% of total P450 genes in larvae and 7 and 4% in adults, respectively, while 5 and 11% were down-regulated in larvae and 4 and 2% in adults. Although the majority of these up- and down-regulated P450 genes appeared to be developmentally controlled, a few were either up- or down-regulated in both the larvae and adult stages. Interestingly, a different gene set was found to be up- or down-regulated in the HAmCq(G8) and MAmCq(G6) mosquito populations in response to insecticide selection. Several genes were identified as being up- or down-regulated in either the larvae or adults for both HAmCq(G8) and MAmCq(G6); of these, CYP6AA7 and CYP4C52v1 were up-regulated and CYP6BY3 was down-regulated across the life stages and populations of mosquitoes, suggesting a link with the permethrin selection in these mosquitoes. Taken together, the findings from this study indicate that not only are multiple P450 genes involved in insecticide resistance but up- or down-regulation of P450 genes may also be co-responsible for detoxification of insecticides, insecticide selection, and the homeostatic response of mosquitoes to changes in cellular environment.

Pyrosequencing the midgut transcriptome of the brown planthopper, Nilaparvata lugens

The brown planthopper, Nilaparvata lugens, is a serious pest threatening rice production across the world. To identify the main features of the gene expression and the key components of the midgut of N. lugens responsible for nutrition, xenobiotic metabolism and the immune response, we used pyrosequencing to sample the transcriptome. More than 190,000 clean sequences were generated, which led to about 30,000 unique sequences. Sequence analysis indicated that genes with abundant transcripts in the midgut of N. lugens were mainly sugar hydrolyases and transporters, proteases and detoxification-related proteins. Based on the sequence information, we cloned the candidate sucrase gene; this enzyme is likely to interact with the perimicrovillar membrane through its highly hydrophobic C-terminal region. Many proteases were identified, which supported the hypothesis that N. lugens uses the proteolysis system for digestion. Scores of detoxification genes were newly identified, including cytochrome P450s, glutathione S-transferases, caroxylesterases. A wealth of new transcripts possibly participating in the immune response were described as well. The gene encoding a peptidoglycan recognition protein was cloned. Unlike in Acyrthosiphon pisum, the immunodeficiency pathway may be present in N. lugens. This is the first global analysis of midgut transcriptome from N. lugens.

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

Background

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

Methodology/Principal Findings

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

Conclusion

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

Multiple insecticide resistance in Aedes aegypti (Diptera: Culicidae) populations compromises the effectiveness of dengue vector control in French Guiana

In French Guiana, pyrethroids and organophosphates have been used for many years against Aedes aegypti. We aimed to establish both the resistance level of Ae. aegypti and the ultra low volume spray efficacy to provide mosquito control services with practical information to implement vector control and resistance management. Resistance to deltamethrin and fenitrothion was observed. In addition, the profound loss of efficacy of AquaK'othrine® and the moderate loss of efficacy of Paluthion® 500 were recorded. Fenitrothion remained the most effective candidate for spatial application in French Guiana until its removal in December 2010. Further investigation of the mechanism of resistance to deltamethrin demonstrated the involvement of mixed-function oxidases and, to a lesser extent, of carboxylesterases. However, these observations alone cannot explain the level of insecticide resistance we observed during tube and cage tests.