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

Enhanced metabolism and selection of pyrethroid-resistant western corn rootworms (Diabrotica virgifera virgifera LeConte)

Western corn rootworm (WCR) pyrethroid resistance has been previously reported in the United States (US) western Corn Belt, and cross-resistance and synergism studies suggested that both target site insensitivity and enhanced metabolism may be conferring WCR resistance to pyrethroids. The present study aimed to investigate the potential mechanisms of WCR pyrethroid resistance and to estimate the heritability of the resistance trait. Biochemical assays using model substrates and spectrophotometry revealed 2-4-fold higher activity of P450s and esterases in pyrethroid-resistant WCR populations, whereas the biological activity of glutathione S-transferase was similar between populations tested. No mutation in the voltage-gated sodium channel was detected in pyrethroid-resistant WCR individuals by sequencing PCR products containing the para-homologous L1014, T929, and M918 amino acid positions that are commonly associated with target site mutations in other pyrethroid-resistant insects. A pilot estimation of pyrethroid resistance heritability obtained during laboratory selection of a WCR population suggested a major genetic component of the resistance trait and predicted a 10-fold increase in WCR bifenthrin resistance within ~7 generations of insecticide lethal exposure. Results support earlier indirect evidence that enhanced metabolism may be contributing to WCR resistance to pyrethroids and illustrates the potential of WCR pyrethroid resistance evolution.

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

Background

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

Methods

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

Results

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

Conclusions

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

Lipase is associated with deltamethrin resistance in Culex pipiens pallens

The wide application of pyrethroids has led to the rapid development of insecticide resistance in mosquitoes, leading to a rise in mosquito-borne diseases. We previously identified five differentially expressed lipase family genes upon evaluating the transcriptomes of deltamethrin-resistant and deltamethrin-susceptible strains of Culex pipiens pallens. Herein, the gene expression levels were verified by quantitative real-time PCR, and two lipase family genes, lipase A and pancreatic triacylglycerol lipase A, were chosen for further investigations. Using cell viability assays and Centers for Disease Control and Prevention bottle bioassays, lipase A was found to increase the resistance of mosquitoes against deltamethrin both in vitro and in vivo. Our findings indicate that lipase A is involved in conferring deltamethrin resistance in Cx. pipiens pallens.

A cytochrome P450 from the mustard leaf beetles hydroxylates geraniol, a key step in iridoid biosynthesis

Larvae of the leaf beetle Phaedon cochleariae synthesize the iridoid chysomelidial via the mevalonate pathway to repel predators. The normal terpenoid biosynthesis is integrated into the dedicated defensive pathway by the ω-hydroxylation of geraniol to (2E,6E)-2,6-dimethylocta-2,6-diene-1,8-diol (ω-OH-geraniol). Here we identify and characterize the P450 monooxygenase CYP6BH5 as the geraniol hydroxylase using integrated transcriptomics, proteomics and RNA interference (RNAi). In the fat body, 73 cytochrome P450s were identified, and CYP6BH5 was among those that were expressed specifically in fat body. Double stranded RNA mediated knockdown of CYP6BH5 led to a significant reduction of ω-hydroxygeraniol glucoside in the hemolymph and, later, of the chrysomelidial in the defensive secretion. Heterologously expressed CYP6BH5 converted geraniol to ω-OH-geraniol. In addition to geraniol, CYP6BH5 also catalyzes hydroxylation of other monoterpenols, such as nerol and citronellol to the corresponding α,ω-dihydroxy compounds.

Using targeted next-generation sequencing to characterize genetic differences associated with insecticide resistance in Culex quinquefasciatus populations from the southern U.S

Resistance to insecticides can hamper the control of mosquitoes such as Culex quinquefasciatus, known to vector arboviruses such as West Nile virus and others. The strong selective pressure exerted on a mosquito population by the use of insecticides can result in heritable genetic changes associated with resistance. We sought to characterize genetic differences between insecticide resistant and susceptible Culex quinquefasciatus mosquitoes using targeted DNA sequencing. To that end, we developed a panel of 122 genes known or hypothesized to be involved in insecticide resistance, and used an Ion Torrent PGM sequencer to sequence 125 unrelated individuals from seven populations in the southern U.S. whose resistance phenotypes to permethrin and malathion were known from previous CDC bottle bioassay testing. Data analysis consisted of discovering SNPs (Single Nucleotide Polymorphism) and genes with evidence of copy number variants (CNVs) statistically associated with resistance. Ten of the seventeen genes found to be present in higher copy numbers were experimentally validated with real-time PCR. Of those, six, including the gene with the knock-down resistance (kdr) mutation, showed evidence of a ≥ 1.5 fold increase compared to control DNA. The SNP analysis revealed 228 unique SNPs that had significant p-values for both a Fisher’s Exact Test and the Cochran-Armitage Test for Trend. We calculated the population frequency for each of the 64 nonsynonymous SNPs in this group. Several genes not previously well characterized represent potential candidates for diagnostic assays when further validation is conducted.

Survey of the Bradysia odoriphaga Transcriptome Using PacBio Single-Molecule Long-Read Sequencing

The damage caused by Bradysia odoriphaga is the main factor threatening the production of vegetables in the Liliaceae family. However, few genetic studies of B. odoriphaga have been conducted because of a lack of genomic resources. Many long-read sequencing technologies have been developed in the last decade; therefore, in this study, the transcriptome including all development stages of B. odoriphaga was sequenced for the first time by Pacific single-molecule long-read sequencing. Here, 39,129 isoforms were generated, and 35,645 were found to have annotation results when checked against sequences available in different databases. Overall, 18,473 isoforms were distributed in 25 various Clusters of Orthologous Groups, and 11,880 isoforms were categorized into 60 functional groups that belonged to the three main Gene Ontology classifications. Moreover, 30,610 isoforms were assigned into 44 functional categories belonging to six main Kyoto Encyclopedia of Genes and Genomes functional categories. Coding DNA sequence (CDS) prediction showed that 36,419 out of 39,129 isoforms were predicted to have CDS, and 4319 simple sequence repeats were detected in total. Finally, 266 insecticide resistance and metabolism-related isoforms were identified as candidate genes for further investigation of insecticide resistance and metabolism in B. odoriphaga.

Expanding the canon: Non-classical mosquito genes at the interface of arboviral infection

Mosquito transmitted viruses cause significant morbidity and mortality in human populations. Despite the use of insecticides and other measures of vector control, arboviral diseases are on the rise. One potential solution for limiting disease transmission to humans is to render mosquitoes refractory to viral infection through genetic modification. Substantial research effort in Drosophila, Aedes and Anopheles has helped to define the major innate immune pathways, including Toll, IMD, Jak/Stat and RNAi, however we still have an incomplete picture of the mosquito antiviral response. Transcriptional profiles of virus-infected insects reveal a much wider range of pathways activated by the process of infection. Within these lists of genes are unexplored mosquito candidates of viral defense. Wolbachia species are endosymbiotic bacteria that naturally limit arboviral infection in mosquitoes. Our understanding of the Wolbachia-mediated viral blocking mechanism is poor, but it does not appear to operate via the classical immune pathways. Herein, we reviewed the transcriptomic response of mosquitoes to multiple viral species and put forth consensus gene types/families outside the immune canon whose expression responds to infection, including cytoskeleton and cellular trafficking, the heat shock response, cytochromes P450, cell proliferation, chitin and small RNAs. We then examine emerging evidence for their functional role in viral resistance in diverse insect and mammalian hosts and their potential role in Wolbachia-mediated viral blocking. These candidate gene families offer novel avenues for research into the nature of insect viral defense.

Identification and classification of differentially expressed genes in pyrethroid-resistant Culex pipiens pallens

Culex pipiens pallens is an important vector that transmits Bancroftian filariasis, Japanese encephalitis and other diseases that pose a serious threat to human health. Extensive and improper use of insecticides has caused insecticide resistance in mosquitoes, which has become an important obstacle to the control of mosquito-borne diseases. It is crucial to investigate the underlying mechanism of insecticide resistance. The aims of this study were to identify genes involved in insecticide resistance based on the resistance phenotype, gene expression profile and single-nucleotide polymorphisms (SNPs) and to screen for major genes controlling insecticide resistance. Using a combination of SNP and transcriptome data, gene expression quantitative trait loci (eQTLs) were studied in deltamethrin-resistant mosquitoes. The most differentially expressed pathway in the resistant group was identified, and a regulatory network was built using these SNPs and the differentially expressed genes (DEGs) in this pathway. The major candidate genes involved in the control of insecticide resistance were analyzed by qPCR, siRNA microinjection and CDC bottle bioassays. A total of 85 DEGs that encoded putative detoxification enzymes (including 61 P450s) were identified in this pathway. The resistance regulatory network was built using SNPs, and these metabolic genes, and a major gene, CYP9AL1, were identified. The functional role of CYP9AL1 in insecticide resistance was confirmed by siRNA microinjection and CDC bottle bioassays. Using the eQTL approach, we identified important genes in pyrethroid resistance that may aid in understanding the mechanism underlying insecticide resistance and in targeting new measures for resistance monitoring and management.

Comparative transcriptome analysis and RNA interference reveal CYP6A8 and SNPs related to pyrethroid resistance in Aedes albopictus

Wide and improper application of pyrethroid insecticides for mosquito control has resulted in widespread resistance in Aedes albopictus mosquitoes, an important dengue vector. Therefore, understanding the molecular regulation of insecticide resistance is urgently needed to provide a basis for developing novel resistance diagnostic methods and vector control approaches. We investigated the transcriptional profiles of deltamethrin-resistant and -susceptible Ae. albopictus by performing paired-end sequencing for RNA expression analysis. The analysis used 24 independent libraries constructed from 12 wild-caught resistant and 12 susceptible Ae. albopictus female adults. A total of 674,503,592 and 612,512,034 reads were obtained, mapped to the Ae. albopictus genome and assembled into 20,091 Ae. albopictus transcripts. A total of 1,130 significantly differentially expressed genes included 874 up-regulated genes and 256 down-regulated genes in the deltamethrin-resistant individuals. These differentially expressed genes code for cytochrome P450s, cuticle proteins, glutathione S-transferase, serine proteases, heat shock proteins, esterase, and others. We selected three highly differentially expressed candidate genes, CYP6A8 and two genes of unknown function (CCG013931 and CCG000656), to test the association between these 3 genes and deltamethrin resistance using RNAi through microinjection in adult mosquitoes and oral feeding in larval mosquitoes. We found that expression knockdown of these three genes caused significant changes in resistance. Further, we detected 1,162 single nucleotide polymorphisms (SNPs) with a frequency difference of more than 50%. Among them, 5 SNPs in 4 cytochrome P450 gene families were found to be significantly associated with resistance in a genotype-phenotype association study using independent field-collected mosquitoes of known resistance phenotypes. Altogether, a combination of novel individually based transcriptome profiling, RNAi, and genetic association study identified both differentially expressed genes and SNPs associated with pyrethroid resistance in Ae. albopictus mosquitoes, and laid a useful foundation for further studies on insecticide resistance mechanisms.

Changes in Neuronal Signaling and Cell Stress Response Pathways are Associated with a Multigenic Response of Drosophila melanogaster to DDT Selection

The adaptation of insect populations to insecticidal control is a continual threat to human health and sustainable agricultural practices, but many complex genomic mechanisms involved in this adaption remain poorly understood. This study applied a systems approach to investigate the interconnections between structural and functional variance in response to dichlorodiphenyltrichloroethane (DDT) within the Drosophila melanogaster strain 91-R. Directional selection in 6 selective sweeps coincided with constitutive gene expression differences in DDT resistant flies, including the most highly upregulated transcript, Unc-115 b, which plays a central role in axon guidance, and the most highly downregulated transcript, the angiopoietin-like CG31832, which is involved in directing vascular branching and dendrite outgrowth but likely may be under trans-regulatory control. Direct functions and protein–protein interactions mediated by differentially expressed transcripts control changes in cell migration, signal transduction, and gene regulatory cascades that impact the nervous system. Although changes to cellular stress response pathways involve 8 different cytochrome P450s, stress response, and apoptosis is controlled by a multifacetted regulatory mechanism. These data demonstrate that DDT selection in 91-R may have resulted in genome-wide adaptations that impacts genetic and signal transduction pathways that converge to modify stress response, cell survival, and neurological functions. This study implicates the involvement of a multigenic mechanism in the adaptation to a chemical insecticide, which impact interconnected regulatory cascades. We propose that DDT selection within 91-R might act systemically, wherein pathway interactions function to reinforce the epistatic effects of individual adaptive changes on an additive or nonadditive basis.

The Effect of Permethrin Resistance on Aedes aegypti Transcriptome Following Ingestion of Zika Virus Infected Blood

Aedes aegypti (L.) is the primary vector of many emerging arboviruses. Insecticide resistance among mosquito populations is a consequence of the application of insecticides for mosquito control. We used RNA-sequencing to compare transcriptomes between permethrin resistant and susceptible strains of Florida Ae. aegypti in response to Zika virus infection. A total of 2459 transcripts were expressed at significantly different levels between resistant and susceptible Ae. aegypti. Gene ontology analysis placed these genes into seven categories of biological processes. The 863 transcripts were expressed at significantly different levels between the two mosquito strains (up/down regulated) more than 2-fold. Quantitative real-time PCR analysis was used to validate the Zika-infection response. Our results suggested a highly overexpressed P450, with AAEL014617 and AAEL006798 as potential candidates for the molecular mechanism of permethrin resistance in Ae. aegypti. Our findings indicated that most detoxification enzymes and immune system enzymes altered their gene expression between the two strains of Ae. aegypti in response to Zika virus infection. Understanding the interactions of arboviruses with resistant mosquito vectors at the molecular level allows for the possible development of new approaches in mitigating arbovirus transmission. This information sheds light on Zika-induced changes in insecticide resistant Ae. aegypti with implications for mosquito control strategies.

Genome-wide and expression-profiling analyses suggest the main cytochrome P450 genes related to pyrethroid resistance in the malaria vector, Anopheles sinensis (Diptera Culicidae)

BACKGROUND

Anopheles sinensis is one of the major malaria vectors. However, pyrethroid resistance in An. sinensis is threatening malaria control. Cytochrome P450-mediated detoxification is an important pyrethroid resistance mechanism that has been unexplored in An. sinensis. In this study, we performed a comprehensive analysis of the An. sinensis P450 gene superfamily with special attention to their role in pyrethroid resistance using bioinformatics and molecular approaches.

RESULTS

Our data revealed the presence of 112 individual P450 genes in An. sinensis, which were classified into four major clans (mitochondrial, CYP2, CYP3 and CYP4), 18 families and 50 subfamilies. Sixty-seven genes formed nine gene clusters, and genes within the same cluster and the same gene family had a similar gene structure. Phylogenetic analysis showed that most of An. sinensis P450s (82/112) had very close 1: 1 orthology with Anopheles gambiae P450s. Five genes (AsCYP6Z2, AsCYP6P3v1, AsCYP6P3v2, AsCYP9J5 and AsCYP306A1) were significantly upregulated in three pyrethroid-resistant populations in both RNA-seq and RT-qPCR analyses, suggesting that they could be the most important P450 genes involved in pyrethroid resistance in An. sinensis.

CONCLUSION

Our study provides insight on the diversity of An. sinensis P450 superfamily and basis for further elucidating pyrethroid resistance mechanism in this mosquito species. © 2018 Society of Chemical Industry.

Gene expression profiles of the Southern house mosquito Culex quinquefasciatus during exposure to permethrin

Insecticide resistance is a major obstacle to the management of disease-vectoring mosquitoes worldwide. The genetic changes and detoxification genes involved in insecticide resistance have been extensively studied in populations of insecticide-resistant mosquitoes, however few studies have focused on the resistance genes upregulated upon insecticide exposure and the possible regulation pathways involved in insecticide resistance. To characterize the changes in gene expression during insecticide exposure, and to investigate the possible connection of known regulation pathways with insecticide resistance, we conducted RNA-Seq analysis of a highly permethrin-resistant strain of Culex quinquefasciatus following permethrin exposure. Gene expression profiles revealed a total of 224 upregulated and 146 downregulated genes when compared to a blank acetone carrier treated control, respectively, suggesting that there were multiple, but specific genes involved in permethrin resistance. Functional enrichment analysis showed that the upregulated genes contained multiple detoxification genes including a glutathione S-transferase and multiple cytochrome P450 genes, as well as several immune-related genes, while the downregulated genes consisted primarily of proteases and carbohydrate metabolism and transport. Further analysis showed that permethrin exposure resulted in a decrease in the expression of serum storage proteins and likely represented a delay in the development of the fourth instar possibly due to a decrease in feeding. This effect was more pronounced in an insecticide-resistant strain than in an insecticide-susceptible strain and may represent a behavioral mechanism of insecticide resistance in Culex mosquitoes.

Expansion of cytochrome P450 and cathepsin genes in the generalist herbivore brown marmorated stink bug

BACKGROUND

The brown marmorated stink bug (Halyomorpha halys) is an invasive pest in North America which causes severe economic losses on tree fruits, ornamentals, vegetables, and field crops. The H. halys is an extreme generalist and this feeding behaviour may have been a major contributor behind its establishment and successful adaptation in invasive habitats of North America. To develop an understanding into the mechanism of H. halys' generalist herbivory, here we specifically focused on genes putatively facilitating its adaptation on diverse host plants.

RESULTS

We generated over 142 million reads via sequencing eight RNA-Seq libraries, each representing an individual H. halys adult. The de novo assembly contained 79,855 high quality transcripts, totalling 39,600,178 bases. Following a comprehensive transcriptome analysis, H. halys had an expanded suite of cytochrome P450 and cathepsin-L genes compared to other insects. Detailed characterization of P450 genes from the CYP6 family, known for herbivore adaptation on host plants, strongly hinted towards H. halys-specific expansions involving gene duplications. In subsequent RT-PCR experiments, both P450 and cathepsin genes exhibited tissue-specific or distinct expression patterns which supported their principal roles of detoxification and/or digestion in a particular tissue.

CONCLUSIONS

Our analysis into P450 and cathepsin genes in H. halys offers new insights into potential mechanisms for understanding generalist herbivory and adaptation success in invasive habitats. Additionally, the large-scale transcriptomic resource developed here provides highly useful data for gene discovery; functional, population and comparative genomics as well as efforts to assemble and annotate the H. halys genome.

Imidacloprid is hydroxylated by Laodelphax striatellus CYP6AY3v2

Laodelphax striatellus (Fallén) is one of the most destructive pests of rice, and has developed high resistance to imidacloprid. Our previous work indicated a strong association between imidacloprid resistance and the overexpression of a cytochrome P450 gene CYP6AY3v2 in a L. striatellus imidacloprid resistant strain (Imid-R). In this study, a transgenic Drosophila melanogaster line that overexpressed the L. striatellus CYP6AY3v2 gene was established and was found to confer increased levels of imidacloprid resistance. Furthermore, CYP6AY3v2 was co-expressed with D. melanogaster cytochrome P450 reductase (CPR) in Spodoptera frugiperda 9 (SF9) cells. A carbon monoxide difference spectra analysis indicated that CYP6AY3v2 was expressed predominately in its cytochrome P450 (P450) form, which is indicative of a good-quality functional enzyme. The recombinant CYP6AY3v2 protein efficiently catalysed the model substrate P-nitroanisole to p-nitrophenol with a maximum velocity (V_max ) of 60.78 ± 3.93 optical density (mOD)/min/mg protein. In addition, imidacloprid itself was metabolized by the recombinant CYP6AY3v2/nicotinamide adenine dinucleotide 2'-phosphate reduced tetrasodium salt (NADPH) CPR microsomes in in vitro assays (catalytic constant (K_cat ) = 0.34 pmol/min/pmol P450, michaelis constant (K_m ) = 41.98 μM), and imidacloprid depletion and metabolite peak formation were with a time dependence. The data provided direct evidence that CYP6AY3v2 is capable of hydroxylation of imidacloprid and conferring metabolic resistance in L. striatellus.

Regulation of P450-mediated permethrin resistance in Culex quinquefasciatus by the GPCR/Gαs/AC/cAMP/PKA signaling cascade

This study explores the role of G-protein-coupled receptor-intracellular signaling in the development of P450-mediated insecticide resistance in mosquitoes, Culex quinquefasciatus, focusing on the essential function of the GPCRs and their downstream effectors of Gs alpha subunit protein (Gαs) and adenylyl cyclase (ACs) in P450-mediated insecticide resistance of Culex mosquitoes. Our RNAi-mediated functional study showed that knockdown of Gαs caused the decreased expression of the downstream effectors of ACs and PKAs in the GPCR signaling pathway and resistance P450 genes, whereas knockdown of ACs decreased the expression of PKAs and resistance P450 genes. Knockdown of either Gαs or ACs resulted in an increased susceptibility of mosquitoes to permethrin. These results add significantly to our understanding of the molecular basis of resistance P450 gene regulation through GPCR/Gαs/AC/cAMP-PKA signaling pathways in the insecticide resistance of mosquitoes. The temporal and spatial dynamic analyses of GPCRs, Gαs, ACs, PKAs, and P450s in two insecticide resistant mosquito strains revealed that all the GPCR signaling pathway components tested, namely GPCRs, Gαs, ACs and PKAs, were most highly expressed in the brain for both resistant strains, suggesting the role played by these genes in signaling transduction and regulation. The resistance P450 genes were mainly expressed in the brain, midgut and malpighian tubules (MTs), suggesting their critical function in the central nervous system and importance for detoxification. The temporal dynamics analysis for the gene expression showed a diverse expression profile during mosquito development, indicating their initially functional importance in response to exposure to insecticides during their life stages.

•

GPCR, Gs alpha subunit protein (Gαs), adenylyl cyclase (ACs), and protein kinase A (PKAs) regulates resistance P450 gene expression and the development of insecticide resistance in mosquitoes, Culex quinquefasciatus.

•

GPCR, Gαs, ACs, and PKAs, are highly expressed in the brain of mosquitoes, corresponding their role in signaling transduction, and regulation.

•

GPCR, Gαs, ACs, PKAs, and P450s are expressed in the different life stages of mosquitoes, revealing their functional importance in response to exposure to insecticides during mosquito life stages.

The role of miR-2∼13∼71 cluster in resistance to deltamethrin in Culex pipiens pallens

Excessive and continuous application of deltamethrin has resulted in the development of deltamethrin resistance among mosquitoes, which becomes a major obstacle for mosquito control. In a previous study, differentially expressed miRNAs between deltamethrin-susceptible (DS) strain and deltamethrin-resistant (DR) strain using illumina sequencing in Culex pipiens pallens were identified. In this study, we applied RNAi and the Centers for Disease Control and Prevention (CDC) bottle bioassay to investigate the relationship between miR-2∼13∼71 cluster (miR-2, miR-13 and miR-71) and deltamethrin resistance. We used quantitative real-time PCR (qRT-PCR) to measure expression levels of miR-2∼13∼71 clusters. MiR-2∼13∼71 cluster was down regulated in adult female mosquitoes from the DR strain and played important roles in deltamethrin resistance through regulating target genes, CYP9J35 and CYP325BG3. Knocking down CYP9J35 and CYP325BG3 resulted in decreased mortality of DR mosquitoes. This study provides the first evidence that miRNA clusters are associated with deltamethrin resistance in mosquitoes. Moreover, we investigated the regulatory networks formed between miR-2∼13∼71 cluster and its target genes, which provide a better understanding of the mechanism involved in deltamethrin resistance.

The function of two P450s, CYP9M10 and CYP6AA7, in the permethrin resistance of Culex quinquefasciatus

Cytochrome P450 monooxygenases play a critical role in insecticide resistance by allowing resistant insects to metabolize insecticides. Previous studies revealed that two P450 genes, CYP9M10 and CYP6AA7, are not only up-regulated but also induced in resistant Culex mosquitoes. In this study, CYP9M10 and CYP6AA7 were separately co-expressed with cytochrome P450 reductase (CPR) in insect Spodoptera frugiperda (Sf9) cells using a baculovirus-mediated expression system and the enzymatic activity and metabolic ability of CYP9M10/CPR and CYP6AA7/CPR to permethrin and its metabolites, including 3-phenoxybenzoic alcohol (PBOH) and 3-phenoxybenzaldehyde (PBCHO), characterized. PBOH and PBCHO, both of which are toxic to Culex mosquito larvae, can be further metabolized by CYP9M10/CPR and CYP6AA7/CPR, with the ultimate metabolite identified here as PBCOOH, which is considerably less toxic to mosquito larvae. A cell-based MTT (3-[4,5-dimethylthiazol-2-yl]-2,5-diphenyltetrazolium bromide) cytotoxicity assay revealed that Sf9 cells expressing CYP9M10/CPR or CYP6AA7/CPR increased the cell line's tolerance to permethrin, PBOH, and PBCHO. This study confirms the important role played by CYP9M10 and CYP6AA7 in the detoxification of permethrin and its metabolites PBOH and PBCHO.

Metabolic resistance in Nilaparvata lugens to etofenprox, a non-ester pyrethroid insecticide

Etofenprox, a non-ester pyrethroid insecticide, will be registered to control rice pests such as the brown planthopper (BPH, Nilaparvata lugens Stål) in mainland China. Insecticide resistance is always a problem to the effective control of N. lugens by chemical insecticides. An etofenprox resistance selection of N. lugens was performed in order to understand the related mechanisms. Through successive selection by etofenprox for 16 generations in the laboratory, an etofenprox-resistant strain (G16) with the resistance ratio (RR) of 422.3-fold was obtained. The resistance was partly synergised (2.68-fold) with the metabolic inhibitor PBO, suggesting a role for P450 monooxygenases. In this study, 11 P450 genes were significantly up-regulated in G16, among which eight genes was above 2.0-fold higher than that in US16, a population with the same origin of G16 but without contacting any insecticide in the laboratory. The expression level of four genes (CYP6AY1, CYP6FU1 and CYP408A1 from Clade 3, and CYP425A1 from Clade 4) were above 4.0-fold when compared to US16. RNA interference (RNAi) was performed to evaluate the importance of the selected P450s in etofenprox resistance. When CYP6FU1, CYP425A1 or CYP6AY1 was interfered, the susceptibility was significantly recovered in both G16 and US16, while the knockdown of CYP408A1 or CYP353D1 did not cause significant changes in etofenprox susceptibility. We supposed that CYP6FU1 was the most important P450 member for etofenprox resistance because of the highest expression level and the most noticeable effects on resistance ratios following RNAi.

Metabolic imidacloprid resistance in the brown planthopper, Nilaparvata lugens, relies on multiple P450 enzymes

Target insensitivity contributing to imidacloprid resistance in Nilaparvata lugens has been reported to occur either through point mutations or quantitative change in nicotinic acetylcholine receptors (nAChRs). However, the metabolic resistance, especially the enhanced detoxification by P450 enzymes, is the major mechanism in fields. From one field-originated N. lugens population, an imidacloprid resistant strain G25 and a susceptible counterpart S25 were obtained to analyze putative roles of P450s in imidacloprid resistance. Compared to S25, over-expression of twelve P450 genes was observed in G25, with ratios above 5.0-fold for CYP6AY1, CYP6ER1, CYP6CS1, CYP6CW1, CYP4CE1 and CYP425B1. RNAi against these genes in vivo and recombinant tests on the corresponding proteins in vitro revealed that four P450s, CYP6AY1, CYP6ER1, CYP4CE1 and CYP6CW1, played important roles in imidacloprid resistance. The importance of the four P450s was not equal at different stages of resistance development based on their over-expression levels, among which CYP6ER1 was important at all stages, and that the others might only contribute at certain stages. The results indicated that, to completely reflect roles of P450s in insecticide resistances, their over-expression in resistant individuals, expression changes at the stages of resistance development, and catalytic activities against insecticides should be considered. In this study, multiple P450s, CYP6AY1, CYP6ER1, CYP4CE1 and CYP6CW1, have proven to be important in imidacloprid resistance.

Transcriptome and Difference Analysis of Fenpropathrin Resistant Predatory Mite, Neoseiulus barkeri (Hughes)

Several fenpropathrin-resistant predatory mites have been reported. However, the molecular mechanism of the resistance remains unknown. In the present study, the Neoseiulus barkeri (N. barkeri) transcriptome was generated using the Illumina sequencing platform, 34,211 unigenes were obtained, and 15,987 were manually annotated. After manual annotation, attentions were attracted to resistance-related genes, such as voltage-gated sodium channel (VGSC), cytochrome P450s (P450s), and glutathione S-transferases (GSTs). A polymorphism analysis detected two point mutations (E1233G and S1282G) in the linker region between VGSC domain II and III. In addition, 43 putative P450 genes and 10 putative GST genes were identified from the transcriptome. Among them, two P450 genes, NbCYP4EV2 and NbCYP4EZ1, and four GST genes, NbGSTd01, NbGSTd02, NbGSTd03 and NbGSTm03, were remarkably overexpressed 3.64-46.69-fold in the fenpropathrin resistant strain compared to that in the susceptible strain. These results suggest that fenpropathrin resistance in N. barkeri is a complex biological process involving many genetic changes and provide new insight into the N. barkeri resistance mechanism.

Molecular Cloning, Expression Pattern and Polymorphisms of NADPH-Cytochrome P450 Reductase in the Bird Cherry-Oat Aphid Rhopalosiphum padi (L.)

NADPH–cytochrome P450 reductase (CPR) plays an important role in the cytochrome P450 (CYP)-mediated metabolism of endogenous and exogenous substrates. CPR has been found to be associated with insecticide metabolism and resistance in many insects. However, information regarding CPR in the bird cherry-oat aphid, Rhopalosiphum padi, is unavailable. In the current study, a full-length cDNA (2,476 bp) of CPR (RpCPR) encoding 681 amino acids was cloned from R. padi. Nucleotide sequence and deduced amino acid sequence analysis showed that RpCPR exhibits characteristics of classical CPRs and shares high identities with those of other insects, especially with the pea aphid, Acyrthosiphon pisum. The mRNA of RpCPR was expressed at all developmental stages, with the highest expression level found in the second instar and the lowest in adult. Expression levels of RpCPR in isoprocarb-resistant and imidacloprid-resistant strains were 3.74- and 3.53-fold higher, respectively, than that of a susceptible strain. RpCPR expression could also be induced by low concentrations (LC₃₀) of isoprocarb and imidacloprid. Moreover, we sequenced the open reading frame (ORF) of RpCPR from 167 field samples collected in 11 geographical populations. Three hundred and thirty-four SNPs were detected, of which, 65 were found in more than two individuals. One hundred and ninety-four missense mutations were present in the amino acid sequence, of which, the P484S mutant had an allele frequency of 35.1%. The present results suggest that RpCPR may play an important role in the P450-mediated insecticide resistance of R. padi to isoprocarb and imidacloprid and possibly other insecticides. Meanwhile, RpCPRmaintains high genetic diversity in natural individuals, which provides the possibility of studying potential correlations between variants and certain special physiological characters.

Transcriptome Analysis of an Insecticide Resistant Housefly Strain: Insights about SNPs and Regulatory Elements in Cytochrome P450 Genes

Background

Insecticide resistance in the housefly, Musca domestica, has been investigated for more than 60 years. It will enter a new era after the recent publication of the housefly genome and the development of multiple next generation sequencing technologies. The genetic background of the xenobiotic response can now be investigated in greater detail. Here, we investigate the 454-pyrosequencing transcriptome of the spinosad-resistant 791spin strain in relation to the housefly genome with focus on P450 genes.

Results

The de novo assembly of clean reads gave 35,834 contigs consisting of 21,780 sequences of the spinosad resistant strain. The 3,648 sequences were annotated with an enzyme code EC number and were mapped to 124 KEGG pathways with metabolic processes as most highly represented pathway. One hundred and twenty contigs were annotated as P450s covering 44 different P450 genes of housefly. Eight differentially expressed P450s genes were identified and investigated for SNPs, CpG islands and common regulatory motifs in promoter and coding regions. Functional annotation clustering of metabolic related genes and motif analysis of P450s revealed their association with epigenetic, transcription and gene expression related functions. The sequence variation analysis resulted in 12 SNPs and eight of them found in cyp6d1. There is variation in location, size and frequency of CpG islands and specific motifs were also identified in these P450s. Moreover, identified motifs were associated to GO terms and transcription factors using bioinformatic tools.

Conclusion

Transcriptome data of a spinosad resistant strain provide together with genome data fundamental support for future research to understand evolution of resistance in houseflies. Here, we report for the first time the SNPs, CpG islands and common regulatory motifs in differentially expressed P450s. Taken together our findings will serve as a stepping stone to advance understanding of the mechanism and role of P450s in xenobiotic detoxification.

Identification of protease m1 zinc metalloprotease conferring resistance to deltamethrin by characterization of an AFLP marker in Culex pipiens pallens

BACKGROUND

Continuous and excessive application of deltamethrin (DM) has resulted in the rapid development of insecticide resistance in Culex pipiens pallens. The quantitative trait loci (QTL) responsible for resistance to DM had previously been detected in Cx. pipiens pallens. But locating the QTLs on the chromosomes remained difficult. An available approach is to first characterize DNA molecular markers linked with the phenotype, and then identify candidate genes.

METHODS

In this study, the amplified fragment length polymorphism (AFLP) marker L3A8.177 associated with the QTL, was characterized. We searched for potential candidate genes in the flank region of L3A8.177 in the genome sequence of the closely related Cx. pipiens quinquefasciatus and conducted mRNA expression analysis of the candidate gene via quantitative real-time PCR. Then the relationship between DM resistance and the candidate gene was identified using RNAi and American CDC Bottle Bioassay in vivo. We also cloned the ORF sequences of the candidate gene from both susceptible and resistant mosquitoes.

RESULTS

The genes CYP6CP1 and protease m1 zinc metalloprotease were in the flank region of L3A8.177 and had significantly different expression levels between susceptible and resistant strains. Protease m1 zinc metalloprotease was significantly up-regulated in the susceptible strains compared with the resistant and remained over-expressed in the susceptible field-collected strains. For deduced amino acid sequences of protease m1 zinc metalloprotease, there was no difference between susceptible and resistant mosquitoes. Knockdown of protease m1 zinc metalloprotease not only decreased the sensitivity of mosquitoes to DM in the susceptible strain but also increased the expression of CYP6CP1, suggesting the role of protease m1 zinc metalloprotease in resistance may be involved in the regulation of the P450 gene expression.

CONCLUSION

Our study represents an example of candidate genes derived from the AFLP marker associated with the QTL and provides the first evidence that protease m1 zinc metalloprotease may play a role in the regulation of DM resistance in Cx. pipiens pallens.

Insecticide Resistance and Management Strategies in Urban Ecosystems

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

A G-protein-coupled receptor regulation pathway in cytochrome P450-mediated permethrin-resistance in mosquitoes, Culex quinquefasciatus

Rhodopsin-like G protein-coupled receptors (GPCRs) are known to be involved in the GPCR signal transduction system and regulate many essential physiological processes in organisms. This study, for the first time, revealed that knockdown of the rhodopsin-like GPCR gene in resistant mosquitoes resulted in a reduction of mosquitoes’ resistance to permethrin, simultaneously reducing the expression of two cAMP-dependent protein kinase A genes (PKAs) and four resistance related cytochrome P450 genes. The function of rhodopsin-like GPCR was further confirmed using transgenic lines of Drosophila melanogaster, in which the tolerance to permethrin and the expression of Drosophila resistance P450 genes were both increased. The roles of GPCR signaling pathway second messenger cyclic adenosine monophosphate (cAMP) and downstream effectors PKAs in resistance were investigated using cAMP production inhibitor Bupivacaine HCl and the RNAi technique. Inhibition of cAMP production led to significant decreases in both the expression of four resistance P450 genes and two PKA genes and mosquito resistance to permethrin. Knockdown of the PKA genes had shown the similar effects on permethrin resistance and P450 gene expression. Taken together, our studies revealed, for the first time, the role of the GPCR/cAMP/PKA-mediated regulatory pathway governing P450 gene expression and P450-mediated resistance in Culex mosquitoes.

Genomic Analysis of Detoxification Supergene Families in the Mosquito Anopheles sinensis

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

Comparative transcriptome analyses of deltamethrin-susceptible and -resistant Culex pipiens pallens by RNA-seq

The widespread and improper use of pyrethroid insecticides, such as deltamethrin, has resulted in the evolution of resistance in many mosquito species, including Culex pipiens pallens. With the development of high-throughput sequencing, it is possible to massively screen pyrethroid resistance-associated gene. In this study, we used Illumina-Solexa transcriptome sequencing to identify genes that are expressed differently in deltamethrin-susceptible and -resistant strains of Culex pipiens pallens as a critical knowledge base for further studies. A total of 4,961,197,620 base pairs and 55,124,418 reads were sequenced, mapped to the Culex quinquefasciatus genome and assembled into 17,679 known genes. We recorded 1826 significantly differentially expressed genes (DEGs). Among them, 1078 genes were up-regulated and 748 genes were down-regulated in the deltamethrin-resistant strain compared to -susceptible strain. These DEGs contained cytochrome P450 s, cuticle proteins, UDP-glucuronosyltransferases, lipases, serine proteases, heat shock proteins, esterases and others. Among the 1826 DEGs, we found that the transcriptional levels of CYP6AA9 in the laboratory populations was elevated as the levels of deltamethrin resistance increased. Moreover, the expression levels of the CYP6AA9 were significantly higher in the resistant strains than the susceptible strains in three different field populations. We further confirmed the association between the CYP6AA9 gene and deltamethrin resistance in mosquitoes by RNA interfering (RNAi). Altogether, we explored massive potential pyrethroid resistance-associated genes and demonstrated that CYP6AA9 participated in the pyrethroid resistance in mosquitoes.

De novo sequencing and characterization of the Bradysia odoriphaga (Diptera: Sciaridae) larval transcriptome

The most serious pestilent threat to the Chinese chive, Allium tuberosum Rottle ex Spreng (Liliaceae) is the Bradysia odoriphaga Yang and Zhang. There is limited genetic research focused on B. odoriphaga, partially due to the lack of genomic resources. The advent of high-throughput sequencing technologies has enabled generation of genomic resources in a short time frame and at minimal costs. In this study, we performed, for the first time, de novo transcriptome sequencing of the B. odoriphaga. Here, 16,829 unigenes were assembled from the total reads, 12,024 of these unigenes were annotated in the NCBI NR protein database, and 9784 were annotated in the Swiss-Prot database. Of these annotated unigenes, 7903 and 5060 unigenes have been assigned to gene ontology categories and clusters of orthologous groups, respectively. Furthermore, 8647 unigenes were mapped to 257 pathways using the Kyoto Encyclopedia of Genes and Genomes Pathway database. We found that 408 unigenes were related to insecticide resistance and metabolism. In addition, 23,122 simple sequence repeats (SSRs) were identified in 11,009 unigenes, and 100 PCR primers of SSR loci were used to validate the assembly quality and polymorphisms. These results provide a good platform for further investigations into the insecticide resistance of B. odoriphaga. Finally, the SSRs identified in B. odoriphaga may be a useful genomic resource.

Venom allergen 5 is Associated With Deltamethrin Resistance in Culex pipiens pallens (Diptera: Culicidae)

The mosquito, Culex pipiens pallens (L.), is an important vector of encephalitis and filariasis in northern China. The control of these mosquitoes occurs primarily via the use of pyrethroid insecticides, such as deltamethrin. The widespread and improper application of pyrethroid has resulted in the evolution of pyrethroid resistance amongst many mosquito populations, including Cx. pipiens pallens. Previous studies using high-throughput transcriptome sequencing have identified that the venom allergen 5 gene is differentially expressed between deltamethrin-susceptible and deltamethrin-resistant Cx. pipiens pallens. In this study, quantitative real-time polymerase chain reaction analyses revealed that venom allergen 5 was significantly overexpressed in adult females of both deltamethrin-resistant laboratory populations and two field populations. The transcriptional level of venom allergen 5 in the laboratory populations was elevated as the levels of deltamethrin resistance increased. Full-length cDNAs of the venom allergen 5 gene were cloned from Cx. pipiens pallens, and contained an open reading frame of 765 bp, encoding a protein with 254 amino acids. The deduced amino acid sequence shared 100% identity with the ortholog in Culex quinquefasciatus Say. The overexpression of venom allergen 5 decreased the susceptibility of mosquito cells to deltamethrin, while knockdown of this gene by RNAi increased the susceptibility of mosquitoes to deltamethrin. This study provides the first evidence of the association between the venom allergen 5 gene and deltamethrin resistance in mosquitoes.

Social evolution. Genomic signatures of evolutionary transitions from solitary to group living

The evolution of eusociality is one of the major transitions in evolution, but the underlying genomic changes are unknown. We compared the genomes of 10 bee species that vary in social complexity, representing multiple independent transitions in social evolution, and report three major findings. First, many important genes show evidence of neutral evolution as a consequence of relaxed selection with increasing social complexity. Second, there is no single road map to eusociality; independent evolutionary transitions in sociality have independent genetic underpinnings. Third, though clearly independent in detail, these transitions do have similar general features, including an increase in constrained protein evolution accompanied by increases in the potential for gene regulation and decreases in diversity and abundance of transposable elements. Eusociality may arise through different mechanisms each time, but would likely always involve an increase in the complexity of gene networks.

Cytochrome P450s--Their expression, regulation, and role in insecticide resistance

P450s are known to be critical for the detoxification and/or activation of xenobiotics such as drugs and pesticides and overexpression of P450 genes can significantly affect the disposition of xenobiotics in the tissues of organisms, altering their pharmacological/toxicological effects. In insects, P450s play an important role in detoxifying exogenous compounds such as insecticides and plant toxins and their overexpression can result in increased levels of P450 proteins and P450 activities. This has been associated with enhanced metabolic detoxification of insecticides and has been implicated in the development of insecticide resistance in insects. Multiple P450 genes have been found to be co-overexpressed in individual insect species via several constitutive overexpression and induction mechanisms, which in turn are co-responsible for high levels of insecticide resistance. Many studies have also demonstrated that the transcriptional overexpression of P450 genes in resistant insects is regulated by trans and/or cis regulatory genes/factors. Taken together, these earlier findings suggest not only that insecticide resistance is conferred via multi-resistance P450 genes, but also that it is mediated through the interaction of regulatory genes/factors and resistance genes. This chapter reviews our current understanding of how the molecular mechanisms of P450 interaction/gene regulation govern the development of insecticide resistance in insects and our progress along the road to a comprehensive characterization of P450 detoxification-mediated insecticide resistance.

Pyrethroid resistance in Culex pipiens mosquitoes

Mosquitoes within the Culex pipiens complex are widely distributed and important in the transmission of many human diseases. Insecticides, pyrethroids in particular, remain a mainstay for control of these important vectors. In this paper we review what is known about the levels, mechanisms and fitness costs of pyrethroid resistance in Cx. pipiens. Pyrethroid resistance in Cx. pipiens is a global problem, and resistance ratios of up to 7000-fold have been found in larvae of field collected mosquitoes. However, there is considerable variation between populations, indicating significant geographic heterogeneity of the resistance. The two major mechanisms of resistance to pyrethroids in Culex are mutations in Vssc (target site insensitivity) and overexpression of cytochrome P450(s) (increased detoxification). The most frequently reported Vssc mutation is L1014F (i.e. kdr), which has been found throughout the world. The L1014S mutation has been found in Cx. p. pallens from Japan and China, and in Cx. p. pipiens from China. The L1014C mutation has only been reported for Cx. p. pipens molestus from China and the V1016G mutation has only been reported from Saudi Arabia. Studies on the P450s of Cx. pipiens have identified several that are overexpressed (measured as transcript levels) in pyrethroid resistant strains. CYP9M10 is consistently overexpressed in pyrethroid resistant Cx. pipiens from at least seven countries, suggesting this P450 might be of global importance in resistance. Both CYP9M10-mediated pyrethroid resistance and kdr have fitness costs in the absence of insecticides under certain environmental conditions. Research needs and future directions are discussed.

Increased expression of CSP and CYP genes in adult silkworm females exposed to avermectins

We analyzed 20 chemosensory protein (CSP) genes of the silkworm Bombyx mori. We found a high number of retrotransposons inserted in introns. We then analyzed expression of the 20 BmorCSP genes across tissues using quantitative real-time polymerase chain reaction (PCR). Relatively low expression levels of BmorCSPs were found in the gut and fat body tissues. We thus tested the effects of endectocyte insecticide abamectin (B1a and B1b avermectins) on BmorCSP gene expression. Quantitative real-time PCR experiments showed that a single brief exposure to insecticide abamectin increased dramatically CSP expression not only in the antennae but in most tissues, including gut and fat body. Furthermore, our study showed coordinate expression of CSPs and metabolic cytochrome P450 enzymes in a tissue-dependent manner in response to the insecticide. The function of CSPs remains unknown. Based on our results, we suggest a role in detecting xenobiotics that are then detoxified by cytochrome P450 anti-xenobiotic enzymes.

Characterization and expression of the cytochrome P450 gene family in diamondback moth, Plutella xylostella (L.)

Cytochrome P450 monooxygenases are present in almost all organisms and can play vital roles in hormone regulation, metabolism of xenobiotics and in biosynthesis or inactivation of endogenous compounds. In the present study, a genome-wide approach was used to identify and analyze the P450 gene family of diamondback moth, Plutella xylostella, a destructive worldwide pest of cruciferous crops. We identified 85 putative cytochrome P450 genes from the P. xylostella genome, including 84 functional genes and 1 pseudogene. These genes were classified into 26 families and 52 subfamilies. A phylogenetic tree constructed with three additional insect species shows extensive gene expansions of P. xylostella P450 genes from clans 3 and 4. Gene expression of cytochrome P450s was quantified across multiple developmental stages (egg, larva, pupa and adult) and tissues (head and midgut) using P. xylostella strains susceptible or resistant to insecticides chlorpyrifos and fiprinol. Expression of the lepidopteran specific CYP367s predominantly occurred in head tissue suggesting a role in either olfaction or detoxification. CYP340s with abundant transposable elements and relatively high expression in the midgut probably contribute to the detoxification of insecticides or plant toxins in P. xylostella. This study will facilitate future functional studies of the P. xylostella P450s in detoxification.

Identification of proteins associated with pyrethroid resistance by iTRAQ-based quantitative proteomic analysis in Culex pipiens pallens

Background

Mosquito control based on chemical insecticides is considered as an important element in the current global strategies for the control of mosquito-borne diseases. Unfortunately, the development of pyrethroid resistance in important vector mosquito species jeopardizes the effectiveness of insecticide-based mosquito control. To date, the mechanisms of pyrethroid resistance are still unclear. Recent advances in proteomic techniques can facilitate to identify pyrethroid resistance-associated proteins at a large-scale for improving our understanding of resistance mechanisms, and more importantly, for seeking some genetic markers used for monitoring and predicting the development of resistance.

Methods

We performed a quantitative proteomic analysis between a deltamethrin-susceptible strain and a deltamethrin-resistant strain of laboratory population of Culex pipiens pallens using isobaric tags for relative and absolute quantitation (iTRAQ) analysis. Gene Ontology (GO) analysis was used to find the relative processes that these differentially expressed proteins were involved in. One differentially expressed protein was chosen to confirm by Western blot in the laboratory and field populations of Cx. pipiens pallens.

Results

We identified 30 differentially expressed proteins assigned into 10 different categories, including oxidoreductase activity, transporter activity, catalytic activity, structural constituent of cuticle and hypothetical proteins. GO analysis revealed that 25 proteins were sub-categorized into 35 hierarchically-structured GO classifications. Western blot results showed that CYP6AA9 as one of the up-regulated proteins was confirmed to be overexpressed in the deltamethrin-resistant strains compared with the deltamethrin-susceptible strains both in the laboratory and field populations.

Conclusions

This is the first study to use modern proteomic tools for identifying pyrethroid resistance-related proteins in Cx. pipiens. The present study brought to light many proteins that were not previously thought to be associated with pyrethroid resistance, which further expands our understanding of pyrethroid resistance mechanisms. CYP6AA9 was overexpressed in the deltamethrin-resistant strains, indicating that CYP6AA9 may be involved in pyrethroid resistance and may be used as a potential genetic marker to monitor and predict the pyrethroid resistance level of field populations.

Electronic supplementary material

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

Insecticide resistance in mosquitoes: impact, mechanisms, and research directions

Mosquito-borne diseases, the most well known of which is malaria, are among the leading causes of human deaths worldwide. Vector control is a very important part of the global strategy for management of mosquito-associated diseases, and insecticide application is the most important component in this effort. However, mosquito-borne diseases are now resurgent, largely because of the insecticide resistance that has developed in mosquito vectors and the drug resistance of pathogens. A large number of studies have shown that multiple, complex resistance mechanisms-in particular, increased metabolic detoxification of insecticides and decreased sensitivity of the target proteins-or genes are likely responsible for insecticide resistance. Gene overexpression and amplification, and mutations in protein-coding-gene regions, have frequently been implicated as well. However, no comprehensive understanding of the resistance mechanisms or regulation involved has yet been developed. This article reviews current knowledge of the molecular mechanisms, genes, gene interactions, and gene regulation governing the development of insecticide resistance in mosquitoes and discusses the potential impact of the latest research findings on the basic and practical aspects of mosquito resistance research.

Inducible Expression of Mu-Class Glutathione S-Transferases Is Associated with Fenpropathrin Resistance in Tetranychus cinnabarinus

The carmine spider mite, Tetranychus cinnabarinus (Boisduval), is a serious pest on a variety of economically important crops widely distributed in China, and its resistance to acaricides has quickly developed. In this study, we fully sequenced 13 GST genes of T. cinnabarinus (TcGSTs). The phylogenetic tree showed that five of them belonged to the delta class and the other eight belonged to the mu class. The alignment of gene sequences and comparison of gene expressions between a fenpropathrin-resistant strain (FR) and a susceptible strain (SS) showed that neither point mutation nor overexpression was detected in TcGSTs. However, when challenged by a sublethal dose of fenpropathrin, the mRNA levels of three GSTs from the mu class (TCGSTM2, TCGSTM3, and TCGSTM8) highly increased in FR, while in SS, the expression of these genes was still at the same level under the treatment. In conclusion, specific TcGSTs were identified that were inducible to stimulation by fenpropathrin, and proved that TcGSTs in FR were not constantly expressed at a high level, but could react much more quickly under the stress of fenpropathrin than SS.

In silico models for predicting vector control chemicals targeting Aedes aegypti

Human arboviral diseases have emerged or re-emerged in numerous countries worldwide due to a number of factors including the lack of progress in vaccine development, lack of drugs, insecticide resistance in mosquitoes, climate changes, societal behaviours, and economical constraints. Thus, Aedes aegypti is the main vector of the yellow fever and dengue fever flaviviruses and is also responsible for several recent outbreaks of the chikungunya alphavirus. As for the other mosquito species, the A. aegypti control relies heavily on the use of insecticides. However, because of increasing resistance to the different families of insecticides, reduction of Aedes populations is becoming increasingly difficult. Despite the unquestionable utility of insecticides in fighting mosquito populations, there are very few new insecticides developed and commercialized for vector control. This is because the high cost of the discovery of an insecticide is not counterbalanced by the 'low profitability' of the vector control market. Fortunately, the use of quantitative structure-activity relationship (QSAR) modelling allows the reduction of time and cost in the discovery of new chemical structures potentially active against mosquitoes. In this context, the goal of the present study was to review all the existing QSAR models on A. aegypti. The homology and pharmacophore models were also reviewed. Specific attention was paid to show the variety of targets investigated in Aedes in relation to the physiology and ecology of the mosquito as well as the diversity of the chemical structures which have been proposed, encompassing man-made and natural substances.

Permethrin resistance variation and susceptible reference line isolation in a field population of the mosquito, Culex quinquefasciatus (Diptera: Culicidae)

This study examines the genetic variations and mechanisms involved in the development of permethrin resistance in individual mosquitoes from a field population of Culex quinquefasciatus, HAmCq(G0) , and characterizes susceptible reference lines of mosquitoes with a similar genetic background to the field HAmCq(G0) strain. Six upregulated cytochrome P450 genes, CYP9M10, CYP9J34, CYP6P14, CYP9J40, CYP6AA7, and CYP4C52v1, previously identified as being upregulated in the larvae of resistant HAmCq(G8) mosquitoes were examined in the larvae of 3 strains (susceptible S-Lab, parental HAmCq(G0) and permethrin-selected highly resistant HAmCq(G8) ) and 8 HAmCq(G0) single-egg raft colonies, covering a range of levels of susceptibility/resistance to permethrin and exhibiting different variations in the expression of A and/or T alleles at the L-to-F kdr locus of the sodium channel. The 2 lines with the lowest tolerance to permethrin and bearing solely the susceptible A allele at the L-to-F kdr locus of the sodium channels, from colonies Cx_SERC5 and Cx_SERC8, showed lower or similar levels of all 6 of the P450 genes tested compared with the S-Lab strain, suggesting that these 2 lines could be used as the reference mosquitoes in future studies characterizing insecticide resistance in HAmCq mosquitoes. This study also provides a detailed investigation of the mechanisms involved in insecticide resistance in individuals within a population: individuals with elevated levels of resistance to permethrin all displayed one or more potential resistance mechanisms-either elevated levels of P450 gene expression, or L-to-F mutations in the sodium channel, or both.

Role of G-protein-coupled receptor-related genes in insecticide resistance of the mosquito, Culex quinquefasciatus

G-protein-coupled receptors regulate signal transduction pathways and play diverse and pivotal roles in the physiology of insects, however, the precise function of GPCRs in insecticide resistance remains unclear. Using quantitative RT-PCR and functional genomic methods, we, for the first time, explored the function of GPCRs and GPCR-related genes in insecticide resistance of mosquitoes, Culex quinquefasciatus. A comparison of the expression of 115 GPCR-related genes at a whole genome level between resistant and susceptible Culex mosquitoes identified one and three GPCR-related genes that were up-regulated in highly resistant Culex mosquito strains, HAmCq^G8 and MAmCq^G6, respectively. To characterize the function of these up-regulated GPCR-related genes in resistance, the up-regulated GPCR-related genes were knockdown in HAmCq^G8 and MAmCq^G6 using RNAi technique. Knockdown of these four GPCR-related genes not only decreased resistance of the mosquitoes to permethrin but also repressed the expression of four insecticide resistance-related P450 genes, suggesting the role of GPCR-related genes in resistance is involved in the regulation of resistance P450 gene expression. This results help in understanding of molecular regulation of resistance development in Cx. quinquefasciatus.

Biochemical and molecular characterisation and cross-resistance in field and laboratory chlorpyrifos-resistant strains of Laodelphax striatellus (Hemiptera: Delphacidae) from eastern China

BACKGROUND

Laboratory selection is often employed in resistance mechanism studies because field-derived populations commonly do not have high enough resistance for such studies. In the present study, a field-collected Laodelphax striatellus population from eastern China was laboratory selected for chlorpyrifos resistance and susceptibility, and the developed strains, along with a field population, were studied for cross-resistance and resistance mechanisms at biochemical and molecular levels.

RESULTS

A 158.58-fold chlorpyrifos-resistant strain (JH-chl) and a chlorpyrifos-susceptible strain (JHS) were established after laboratory selection of 25 generations. Cross-resistance to deltamethrin, diazinon, methomyl, carbosulfan, acephate and imidacloprid were detected in JH-chl and a field-collected strain (JHF). Synergism and enzyme activity data suggested potential involvement of P450s and esterases in JH-chl as well as AChE alteration. Furthermore, CYP6AY3v2, CYP306A2v2, CYP353D1v2 and LSCE36 genes were significantly overexpressed in JH-chl (6.87-12.14-fold). Feeding of dsRNAs reduced the expression of the four target genes (35.6-56.8%) and caused significant adult mortality (75.21-88.45%), implying resistance reduction. However, mechanism(s) conferring chlorpyrifos resistance in JHF were unclear.

CONCLUSION

In contrast to previous reports, multiple overexpressed detoxification genes were potentially associated with chlorpyrifos resistance, as confirmed by RNAi feeding tests. Chlorpyrifos resistance exhibits cross-resistance with insecticides in the same and different classes.

Transcriptome profiling of pyrethroid resistant and susceptible mosquitoes in the malaria vector, Anopheles sinensis

BACKGROUND

Anopheles sinensis is a major malaria vector in China and other Southeast Asian countries, and it is becoming increasingly resistant to the insecticides used for agriculture, net impregnation, and indoor residual spray. Very limited genomic information on this species is available, which has hindered the development of new tools for resistance surveillance and vector control. We used the 454 GS FLX system and generated expressed sequence tag (EST) databases of various life stages of An. sinensis, and we determined the transcriptional differences between deltamethrin resistant and susceptible mosquitoes.

RESULTS

The 454 GS FLX transcriptome sequencing yielded a total of 624,559 reads (average length of 290 bp) with the pooled An. sinensis mosquitoes across various development stages. The de novo assembly generated 33,411 contigs with average length of 493 bp. A total of 8,057 ESTs were generated with Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation. A total of 2,131 ESTs were differentially expressed between deltamethrin resistant and susceptible mosquitoes collected from the same field site in Jiangsu, China. Among these differentially expressed ESTs, a total of 294 pathways were mapped to the KEGG database, with the predominant ESTs belonging to metabolic pathways. Furthermore, a total of 2,408 microsatellites and 15,496 single nucleotide polymorphisms (SNPs) were identified.

CONCLUSIONS

The annotated EST and transcriptome databases provide a valuable genomic resource for further genetic studies of this important malaria vector species. The differentially expressed ESTs associated with insecticide resistance identified in this study lay an important foundation for further functional analysis. The identified microsatellite and SNP markers will provide useful tools for future population genetic and comparative genomic analyses of malaria vectors.

Transcriptional analysis of four family 4 P450s in a Puerto Rico strain of Aedes aegypti (Diptera: Culicidae) compared with an Orlando strain and their possible functional roles in permethrin resistance

A field strain of Aedes aegypti (L.) was collected from Puerto Rico in October 2008. Based on LD50 values by topical application, the Puerto Rico strain was 73-fold resistant to permethrin compared with a susceptible Orlando strain. In the presence of piperonyl butoxide, the resistance of Puerto Rico strain of Ae. aegypti was reduced to 15-fold, suggesting that cytochrome P450-mediated detoxification is involved in the resistance of the Puerto Rico strain to permethrin. To determine the cytochrome P450s that might play a role in the resistance to permethrin, the transcriptional levels of 164 cytochrome P450 genes in the Puerto Rico strain were compared with that in the Orlando strain. Of the 164 cytochrome P450s, 33 were significantly (P < 0.05) up-regulated, including cytochrome P450s in families four, six, and nine. Multiple studies have investigated the functionality of family six and nine cytochrome P450s, therefore, we focused on the up-regulated family 4 cytochrome P450s. To determine whether up-regulation of the four cytochrome P450s had any functional role in permethrin resistance, transgenic Drosophila melanogaster Meigen lines overexpressing the four family 4 P450 genes were generated, and their ability to survive exposure to permethrin was evaluated. When exposed to 5 microg per vial permethrin, transgenic D. melanogaster expressing CYP4D24, CYP4H29, CYP4J15v1, and CYP4H33 had a survival rate of 60.0 +/- 6.7, 29.0 +/- 4.4, 64.4 +/- 9.7, and 11.0 +/- 4.4%, respectively. However, none of the control flies survived the permethrin exposure at the same concentration. Similarly, none of the transgenic D. melanogaster expressing CYP4J15v1 or CYP4H33 ?5 survived when they were exposed to permethrin at 10 microg per vial. However, transgenic D. melanogaster expressing CYP4D24 and CYP4H29 had a survival rate of 37.8 +/- 4.4 and 2.2 +/- 2.2%, respectively. Taken together, our results suggest that CYP4D24 might play an important role in cytochrome P450-mediated resistance to permethrin.

Genomic and bioinformatic analysis of NADPH-cytochrome P450 reductase in Anopheles stephensi (Diptera: Culicidae)

The cytochrome P450 monooxygenase (P450) enzyme system is a major mechanism of xenobiotic biotransformation. The nicotinamide adenine dinucleotide phosphate (NADPH)-cytochrome P450 reductase (CPR) is required for transfer of electrons from NADPH to P450. One CPR gene was identified in the genome of the malaria-transmitting mosquito Anopheles stephensi Liston (Diptera: Culicidae). The gene encodes a polypeptide containing highly conserved flavin mononucleotide-, flavin adenine dinucleotide-, and NADPH-binding domains, a unique characteristic of the reductase. Phylogenetic analysis revealed that the A. stephensi and other known mosquito CPRs belong to a monophyletic group distinctly separated from other insects in the same order, Diptera. Amino acid residues of CPRs involved in binding of P450 and cytochrome c are conserved between A. stephensi and the Norway rat Rattus norvegicus Berkenhout (Rodentia: Muridae). However, gene structure particularly within the coding region is evidently different between the two organisms. Such difference might arise during the evolution process as also seen in the difference of P450 families and isoforms found in these organisms. CPR in the mosquito A. stephensi is expected to be active and serve as an essential component of the P450 system.

Detoxification and stress response genes expressed in a western North American bumble bee, Bombus huntii (Hymenoptera: Apidae)

Background

The Hunt bumble bee (Bombus huntii Greene, Hymenoptera: Apidae) is a holometabolous, social insect important as a pollinator in natural and agricultural ecosystems in western North America. Bumble bees spend a significant amount of time foraging on a wide variety of flowering plants, and this activity exposes them to both plant toxins and pesticides, posing a threat to individual and colony survival. Little is known about what detoxification pathways are active in bumble bees, how the expression of detoxification genes changes across life stages, or how the number of detoxification genes expressed in B. huntii compares to other insects.

Results

We found B. huntii expressed at least 584 genes associated with detoxification and stress responses. The expression levels of some of these genes, such as those encoding the cytochrome P450s, glutathione S-transferases (GSTs) and glycosidases, vary among different life stages to a greater extent than do other genes. We also found that the number of P450s, GSTs and esterase genes expressed by B. huntii is similar to the number of these genes found in the genomes of other bees, namely Bombus terrestris, Bombus impatiens, Apis mellifera and Megachile rotundata, but many fewer than are found in the fly Drosophila melanogaster.

Conclusions

Bombus huntii has transcripts for a large number of detoxification and stress related proteins, including oxidation and reduction enzymes, conjugation enzymes, hydrolytic enzymes, ABC transporters, cadherins, and heat shock proteins. The diversity of genes expressed within some detoxification pathways varies among the life stages and castes, and we typically identified more genes in the adult females than in larvae, pupae, or adult males, for most pathways. Meanwhile, we found the numbers of detoxification and stress genes expressed by B. huntii to be more similar to other bees than to the fruit fly. The low number of detoxification genes, first noted in the honey bee, appears to be a common phenomenon among bees, and perhaps results from their symbiotic relationship with plants. Many flowering plants benefit from pollinators, and thus offer these insects rewards (such as nectar) rather than defensive plant toxins.

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.

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.