Microarray analysis of cytochrome P450 mediated insecticide resistance in Drosophila

Up-regulation of CYP6G4 mediated by a CncC/maf binding-site-containing insertion confers resistance to multiple classes of insecticides in the house fly Musca domestica

Populations of many insect species have evolved a variety of resistance mechanisms in response to insecticide selection. Current knowledge about mutations responsible for insecticide resistance is largely achieved from studies on target-site resistance, while much less is known about metabolic resistance. Although it is well known that P450 monooxygenases are one of the major players involved in insecticide metabolism and resistance, understanding mutation(s) responsible for CYP-mediated resistance has been a big challenge. In this study, we used the house fly to pursue a better understanding of P450 mediated insecticide resistance at the molecular level. Metabolism studies illustrated that CYP6G4 had a broad-spectrum metabolic activity in metabolizing insecticides. Population genotyping revealed that the CYP6G4v1 allele harboring a DNA insertion (MdIS1) had been selected in many house fly populations on different continents. Dual luciferase reporter assays identified that the MdIS1 contained a CncC/Maf binding site, and electrophoretic mobility shift assay confirmed that transcription factor CncC was involved in the MdIS1-mediated regulation. This study highlights the common involvement of the CncC pathway in adaptive evolution, and provides an interesting case supportive of parallel evolution in P450-mediated insecticide resistance in insects.

Resistance to dicyclanil and imidacloprid in the sheep blowfly, Lucilia cuprina, in Australia

BACKGROUND

The sheep blowfly, Lucila cuprina, is a myiasis-causing parasite responsible for significant production losses and welfare issues for the Australian sheep industry. Control relies largely on the use of insecticides. The pyrimidine compound, dicyclanil, is the predominant control chemical, although other insecticides also are used, including imidacloprid, ivermectin, cyromazine and spinosad. We investigated in vitro resistance patterns and mechanisms in field-collected blowfly strains.

RESULTS

The Walgett 2019 strain showed significant levels of resistance to both dicyclanil and imidacloprid, with resistance factors at the IC₅₀ of 26- and 17-fold, respectively, in in vitro bioassays. Co-treatment with the cytochrome P450 inhibitor, aminobenzotriazole, resulted in significant levels of synergism for dicyclanil and imidacloprid (synergism ratios of 7.2- and 6.1-fold, respectively), implicating cytochrome P450 in resistance to both insecticides. Cyp12d1 transcription levels were increased up to 40-fold throughout the larval life stages in the resistant strain compared to a reference susceptible strain, whereas transcription levels of some other cyp genes (6g1, 4d1, 28d1) did not differ between the strains. Similar resistance levels also were observed in flies collected from the same property in two subsequent years.

CONCLUSION

This study indicates that in vitro resistance to both dicyclanil and imidacloprid in this field-collected blowfly strain is likely mediated by cytochrome P450, with Cyp12d1 implicated as the enzyme responsible; however, it remains possible that another P450 also may be involved. A common resistance mechanism for the two drugs has important implications for drug rotation strategies designed to prolong the useful life of flystrike control chemicals. © 2022 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Transcriptional and functional motifs defining renal function revealed by single-nucleus RNA sequencing

Recent advances in single-cell sequencing provide a unique opportunity to gain novel insights into the diversity, lineage, and functions of cell types constituting a tissue/organ. Here, we performed a single-nucleus study of the adult Drosophila renal system, consisting of Malpighian tubules and nephrocytes, which shares similarities with the mammalian kidney. We identified 11 distinct clusters representing renal stem cells, stellate cells, regionally specific principal cells, garland nephrocyte cells, and pericardial nephrocytes. Characterization of the transcription factors specific to each cluster identified fruitless (fru) as playing a role in stem cell regeneration and Hepatocyte nuclear factor 4 (Hnf4) in regulating glycogen and triglyceride metabolism. In addition, we identified a number of genes, including Rho guanine nucleotide exchange factor at 64C (RhoGEF64c), Frequenin 2 (Frq2), Prip, and CG1093 that are involved in regulating the unusual star shape of stellate cells. Importantly, the single-nucleus dataset allows visualization of the expression at the organ level of genes involved in ion transport and junctional permeability, providing a systems-level view of the organization and physiological roles of the tubules. Finally, a cross-species analysis allowed us to match the fly kidney cell types to mouse kidney cell types and planarian protonephridia, knowledge that will help the generation of kidney disease models. Altogether, our study provides a comprehensive resource for studying the fly kidney.

Comparative transcriptome and RNA interference reveal CYP6DC1 and CYP380C47 related to lambda-cyhalothrin resistance in Rhopalosiphum padi

The bird-cherry-oat aphid, Rhopalosiphum padi, is a serious agricultural pest of Triticeae crops, and pyrethroids are the most widely used chemical pesticides for the control of the aphid. Our previous studies found that some R. padi field populations have developed resistance against pyrethroids; an M918L target-site mutation of the voltage gated sodium channel was present in the pyrethroid resistant individuals, while the high-level resistance to lambda-cyhalothrin revealed the presence of other mechanisms in the pest. Here, we conducted genome-wide transcriptional analysis for the lambda-cyhalothrin susceptible (SS) and resistant (LC-RR) strains of R. padi. Results indicated that 2457 genes were differently expressed between the SS and LC-RR strains. In the LC-RR, a total of 1265 and 1192 genes were up- and down-regulated, respectively. KEGG analysis implicated enrichment of P450 involved in insecticide metabolic pathways in the resistant transcriptome. qRT-PCR results confirmed that two P450 genes (CYP6DC1 and CYP380C47) were significantly overexpressed in the LC-RR individuals. Furthermore, RNA interference (RNAi) of CYP6DC1 or CYP380C47 significantly increased mortality of R. padi exposure to lambda-cyhalothrin. These results suggest that the overexpression of CYP6DC1 and CYP380C47 contributed to the lambda-cyhalothrin resistance in the pest. This study provides knowledge for further analyzing the molecular mechanism of resistance to pyrethroids in R. padi.

Experimental evolution supports the potential of neonicotinoid-pyrethroid combination for managing insecticide resistance in malaria vectors

The introduction of neonicotinoids for managing insecticide resistance in mosquitoes is of high interest as they interact with a biochemical target not previously used in public health. In this concern, Bayer developed a combination of the neonicotinoid clothianidin and the pyrethroid deltamethrin (brand name Fludora Fusion) as a new vector control tool. Although this combination proved to be efficient against pyrethroid-resistant mosquitoes, its ability to prevent the selection of pyrethroid and neonicotinoid resistance alleles was not investigated. In this context, the objective of this work was to study the dynamics and the molecular mechanisms of resistance of An. gambiae to the separated or combined components of this combination. A field-derived An. gambiae line carrying resistance alleles to multiple insecticides at low frequencies was used as a starting for 33 successive generations of controlled selection. Resistance levels to each insecticide and target site mutation frequencies were monitored throughout the selection process. Cross resistance to other public health insecticides were also investigated. RNA-seq was used to compare gene transcription variations and polymorphisms across all lines. This study confirmed the potential of this insecticide combination to impair the selection of resistance as compared to its two separated components. Deltamethrin selection led to the rapid enrichment of the kdr L1014F target-site mutation. Clothianidin selection led to the over-transcription of multiple cytochrome P450s including some showing high homology with those conferring neonicotinoid resistance in other insects. A strong selection signature associated with clothianidin selection was also observed on a P450 gene cluster previously associated with resistance. Within this cluster, the gene CYP6M1 showed the highest selection signature together with a transcription profile supporting a role in clothianidin resistance. Modelling the impact of point mutations selected by clothianidin on CYP6M1 protein structure showed that selection retained a protein variant with a modified active site potentially enhancing clothianidin metabolism. In the context of the recent deployment of neonicotinoids for mosquito control and their frequent usage in agriculture, the present study highlights the benefit of combining them with other insecticides for preventing the selection of resistance and sustaining vector control activities.

Molecular bases of P450-mediated resistance to the neonicotinoid insecticide imidacloprid in the mosquito Ae. aegypti

Resistance to chemical insecticides including pyrethroids, the main insecticide class used against mosquitoes, has re-kindled interest in the use of neonicotinoids. In this context, the present study aimed to characterize the molecular basis of neonicotinoid resistance in the mosquito Aedes aegypti. Resistance mechanisms were studied by combining transcriptomic and genomic data obtained from a laboratory strain selected at the larval stage after 30 generations of exposure to imidacloprid (Imida-R line). After thirty generations of selection, larvae of the Imida-R line showed an 8-fold increased resistance to imidacloprid and a significant cross-tolerance to the pyrethroids permethrin and deltamethrin. Cross-resistance to pyrethroids was only observed in adults when larvae were previously exposed to imidacloprid suggesting a low but inducible expression of resistance alleles at the adult stage. Resistance of the Imida-R line was associated with a slower larval development time in females. Multiple detoxification enzymes were over-transcribed in larvae in association with resistance including the P450s CYP6BB2, CYP9M9 and CYP6M11 previously associated with pyrethroid resistance. Some of them together with their redox partner NADPH P450 reductase were also affected by non-synonymous mutations associated with resistance. Combining genomic and transcriptomic data allowed identifying promoter variations associated with the up-regulation of CYP6BB2 in the resistant line. Overall, these data confirm the key role of P450s in neonicotinoid resistance in Ae. aegypti and their potential to confer cross-resistance to pyrethroids, raising concerns about the use of neonicotinoids for resistance management in this mosquito species.

Resistance in the Genus Spodoptera: Key Insect Detoxification Genes

The genus Spodoptera (Lepidoptera: Noctuidae) includes species that are among the most important crop pests in the world. These polyphagous species are able to feed on many plants, including corn, rice and cotton. In addition to their ability to adapt to toxic compounds produced by plants, they have developed resistance to the chemical insecticides used for their control. One of the main mechanisms developed by insects to become resistant involves detoxification enzymes. In this review, we illustrate some examples of the role of major families of detoxification enzymes such as cytochromes P450, carboxyl/cholinesterases, glutathione S-transferases (GST) and transporters such as ATP-binding cassette (ABC) transporters in insecticide resistance. We compare available data for four species, Spodoptera exigua, S. frugiperda, S. littoralis and S. litura. Molecular mechanisms underlying the involvement of these genes in resistance will be described, including the duplication of the CYP9A cluster, over-expression of GST epsilon or point mutations in acetylcholinesterase and ABCC2. This review is not intended to be exhaustive but to highlight the key roles of certain genes.

Transcriptome analysis unravels RNAi pathways genes and putative expansion of CYP450 gene family in cotton leafhopper Amrasca biguttula (Ishida)

Cotton Leafhopper, Amrasca biguttula is an important pest of cotton and okra in the Indian subcontinent. Presently limited genomic/transcriptomic information is available for this insect in any of open source databases. The present study reports the first assembled and annotated de novo transcriptome of cotton leafhopper. Out of 75,551 transcripts, 39,613 CDS (Coding Sequence) were predicted with 35,282 showing positive blast hits with NCBI nr database. The Gene ontology (GO) analysis annotated 7431 CDS  with KEGG pathway categorizing these CDS into 22 different functional groups. The majority of CDS were annotated in signal transduction and transport catabolism pathways. The sequence data was screened for RNAi pathway genes and presence of 37 transcripts associated with this process confirmed the existence of robust RNAi machinery. The role of core RNAi machinery genes (Dicer-2, Ago-2, Piwi and Staufen) has been validated through dsRNA feeding studies. The data resource has also been used to identify potential RNAi targets and genes associated with insecticide detoxification specifically CYP 450 family. The current study provides a useful sequence resource which can be used to initiate molecular studies in this insect with emphasis on insecticide resistance, RNAi and functional genomics.

Molecular characterization of western corn rootworm pyrethroid resistance

BACKGROUND

Western corn rootworm (WCR) pyrethroid resistance has been confirmed in the western US Corn Belt. Toxicological and biochemical studies indicated that multiple mechanisms of resistance might be involved in the resistance trait, such as enhanced metabolism and/or kdr target-site mutation(s) in the voltage-gated sodium channels. To characterize the mechanisms of WCR pyrethroid resistance at the molecular level, pairwise comparisons were made between RNA-Seq data collected from pyrethroid-resistant and -susceptible WCR populations. Gene expression levels and sodium channel sequences were evaluated.

RESULTS

Seven transcripts exhibited significantly different expression (q ≤ 0.05) when comparing field-collected pyrethroid-resistant (R-Field) and -susceptible (S-Field) WCR populations. Three of the differentially expressed transcripts were P450s overexpressed in R-Field (9.2-26.2-fold). A higher number (99) of differentially expressed transcripts was found when comparing laboratory-derived pyrethroid-resistant (R-Lab) and -susceptible (S-Lab) WCR populations. Eight of the significant transcripts were P450s overexpressed in R-Lab (2.7-39.8-fold). This study did not detect kdr mutations in pyrethroid-resistant WCR populations. Other differentially expressed transcripts that may play a role in WCR pyrethroid resistance are discussed.

CONCLUSION

This study revealed that P450-mediated metabolism is likely to be a major mechanism of WCR pyrethroid resistance, which could affect the efficacy of other insecticides sharing similar metabolic pathways. Additionally, results suggested that although laboratory selection of a pyrethroid-resistant WCR population may help to characterize resistance mechanisms, a field-selected population provided rare and perhaps major variants corresponding to the resistance trait.

Molecular basis of resistance to organophosphate insecticides in the New World screw-worm fly

Background

The emergence of insecticide resistance is a fast-paced example of the evolutionary process of natural selection. In this study, we investigated the molecular basis of resistance in the myiasis-causing fly Cochliomyia hominivorax (Diptera: Calliphoridae) to dimethyl-organophosphate (OP) insecticides.

Methods

By sequencing the RNA from surviving larvae treated with dimethyl-OP (resistant condition) and non-treated larvae (control condition), we identified genes displaying condition-specific polymorphisms, as well as those differentially expressed.

Results

Both analyses revealed that resistant individuals have altered expression and allele-specific expression of genes involved in proteolysis (specifically serine-endopeptidase), olfactory perception and cuticle metabolism, among others. We also confirmed that resistant individuals carry almost invariably the Trp251Ser mutation in the esterase E3, known to confer OP and Pyrethroid resistance. Interestingly, genes involved in metabolic and detoxifying processes (notably cytochrome P450s) were found under-expressed in resistant individuals. An exception to this were esterases, which were found up-regulated.

Conclusions

These observations suggest that reduced penetration and aversion to dimethyl-OP contaminated food may be important complementary strategies of resistant individuals. The specific genes and processes found are an important starting point for future functional studies. Their role in insecticide resistance merits consideration to better the current pest management strategies.

Dietary antioxidant vitamin C influences the evolutionary path of insecticide resistance in Drosophila melanogaster

Herbivorous insects encounter a variety of toxic environmental substances ranging from ingested plant defensive compounds to human-introduced insecticidal agents. Dietary antioxidants are known to reduce the negative physiological impacts of toxins in mammalian systems through amelioration of reactive oxygen-related cellular damage. The analogous impacts to insects caused by multigenerational exposure to pesticides and the effects on adaptive responses within insect populations, however, are currently unknown. To address these research gaps, we used Drosophila as a model system to explore adaptive phenotypic responses to acute dichlorodiphenyltrichloroethane (DDT) exposure in the presence of the dietary antioxidant vitamin C and to examine the structural genomic consequences of this exposure. DDT resistance increased significantly among four replicates exposed to a low concentration of DDT for 10 generations. In contrast, dietary intake of vitamin C significantly reduced DDT resistance after mutigenerational exposure to the same concentration of DDT. As to the genomic consequences, no significant differences were predicted in overall nucleotide substitution rates across the genome between any of the treatments. Despite this, replicates exposed to a low concentration of DDT without vitamin C showed the highest number of synonymous and non-synonymous variants (3196 in total), followed by the DDT plus vitamin C (1174 in total), and vitamin C alone (728 in total) treatments. This study demonstrates the potential role of diet (specifically, antioxidant intake) on adaptive genome responses, and thus on the evolution of pesticide resistance within insect populations.

The overexpression of three cytochrome P450 genes CYP6CY14, CYP6CY22 and CYP6UN1 contributed to metabolic resistance to dinotefuran in melon/cotton aphid, Aphis gossypii Glover

Dinotefuran, the third-generation neonicotinoid, has been applied against melon/cotton aphid Aphis gossypii Glover in China. The risk of resistance development, cross-resistance pattern and potential resistance mechanism of dinotefuran in A. gossypii were investigated. A dinotefuran-resistant strain of A. gossypii (DinR) with 74.7-fold resistance was established by continuous selection using dinotefuran. The DinR strain showed a medium level of cross resistance to thiamethoxam (15.2-fold), but no cross resistance to imidacloprid. The synergism assay indicated that piperonyl butoxide and triphenyl phosphate showed synergistic effects on dinotefuran toxicity to the DinR strain with a synergistic ratio of 8.3 and 2.5, respectively, while diethyl maleate showed no synergistic effect. The activities of cytochrome P450 monooxygenase and carboxylesterase were significantly higher in DinR strain than in susceptible strain (SS). Moreover, the gene expression results showed that CYP6CY14, CYP6CY22 and CYP6UN1 were significantly overexpressed in DinR strain compared with SS strain. The expression of CYP6CY14 was 5.8-fold higher in DinR strain than in SS strain. Additionally, the transcription of CYP6CY14, CYP6CY22 and CYP6UN1 in A. gossypii showed dose- and time-dependent responses to dinotefuran exposure. Furthermore, knockdown of CYP6CY14, CYP6CY22 and CYP6UN1 via RNA interference (RNAi) significantly increased mortality of A. gossypii, when A. gossypii was treated with dinotefuran. These results demonstrated the overexpression of CYP6CY14, CYP6CY22 and CYP6UN1 contributed to dinotefuran resistance in A. gossypii.

Toxicogenomics provides insights to toxicity pathways of neonicotinoids to aquatic insect, Chironomus dilutus

Neonicotinoid insecticides have posed a great threat to non-target organisms, yet the mechanisms underlying their toxicity are not well characterized. Major modes of action (MoAs) of imidacloprid were analyzed in an aquatic insect Chironomus dilutus. Lethal and sublethal outcomes were assessed in the midges after 96-h exposure to imidacloprid. Global transcriptomic profiles were determined using de novo RNA-sequencing to more holistically identify toxicity pathways. Transcriptional 10% biological potency values derived from ranked KEGG pathways and GO terms were 0.02 (0.01-0.08) (mean (95% confidence interval) and 0.05 (0.04-0.06) μg L^-1, respectively, which were more sensitive than those from phenotypic traits (10% lethal concentration: 0.44 (0.23-0.79) μg L^-1; 10% burrowing behavior concentration: 0.30 (0.22-0.43) μg L^-1). Major MoAs of imidacloprid in aquatic species were identified as follows: the activation of nicotinic acetylcholine receptors (nAChRs) induced by imidacloprid impaired organisms' nerve system through calcium ion homeostasis imbalance and mitochondrial dysfunction, which posed oxidative stress and DNA damage and eventually caused death of organisms. The current investigation highlighted that imidacloprid affected C. dilutus at environmentally relevant concentrations, and elucidated toxicity pathways derived from gene alteration to individual outcomes, calling for more attention to toxicity of neonicotinoids to aquatic organisms.

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.

The role of G protein-coupled receptor-related genes in cytochrome P450-mediated resistance of the house fly, Musca domestica (Diptera: Muscidae), to imidacloprid

Ninety-four putative G protein-coupled receptors (GPCRs) were identified in the Musca domestica genome. They were annotated and compared with their homologues in Drosophila melanogaster. Phylogenetic analyses of the GPCRs from both species revealed that several family members shared a closer relationship based on the domain architecture. The expression profiles of these genes were examined by quantitative real-time PCR amongst three strains of the house fly, a near-isogenic line strain with imidacloprid resistance (N-IRS), the corresponding susceptible strain (CSS) and another strain derived from field populations with imidacloprid resistance (IRS). We found that five GPCR genes were upregulated in the N-IRS and eight GPCR genes were upregulated in the IRS strains compared to the CSS strain. The transgenic lines of D. melanogaster with the GPCR genes (LOC101899380 in the N-IRS strain and LOC101895664 in the IRS strain) exhibited significantly increased tolerance to imidacloprid, and higher expression of cytochrome P450 genes. Bioinformatic analysis of LOC101899380 was carried out based on its full-length nucleic acid sequence and putative amino acid sequence, and it was named Methuselah-like10 (Mthl10) owing to its homology with D. melanogaster Mthl10. A cell-base cell counting kit-8 toxicity assay demonstrated that the expression of the GPCR gene LOC101899380 in Spodoptera frugiperda (Sf9) cells using a baculovirus-mediated expression system can elevate the cell tolerance to imidacloprid, indirectly supporting the hypothesis that the GPCR gene LOC101899380 plays some role in imidacloprid resistance. These results should be useful for furthering understanding of the regulatory pathway by which house flies develop resistance.

A review of DDT resistance as it pertains to the 91-C and 91-R strains in Drosophila melanogaster

While insecticide resistance presents a challenge for those intent on controlling insect populations, these challenges have also generated a set of tools that can be used to ask fundamental biological questions about that resistance. Numerous species of insects have evolved resistance to multiple classes of insecticides. Each one of these species and their respective resistant populations represent a potential tool for understanding the molecular basis of the evolution of resistance. However, in-laboratory maintenance of resistant insect populations (and their comparative susceptible populations) suitable for asking the needed set of questions around the molecular consequences of long-term pesticide exposure requires a significant, in places prohibitive, level of resources. Drosophila melanogaster (hereafter referred to as Drosophila) is a model insect system with populations easily selected with pesticides and readily maintainable over decades. Even within Drosophila, however, few populations exist where long-term pesticide selection has occurred along with contrasting non-selected population. As such, the Drosophila 91-C and 91-R populations, which exhibit insecticide resistance to DDT (91-R), compared to a non-selection population (91-C), represent a unique resource for the study of high level DDT resistance. Moreover, with the availability of "omics" technologies over the past several decades, this paired population has emerged as a useful tool for understanding both the molecular basis of pesticide resistance and the molecular consequences of long-term pesticide exposure. In this review, we summarize the studies with these aforementioned populations over the past several decades, addressing what has been learned from these efforts.

Development of a goosegrass (Eleusine indica) draft genome and application to weed science research

BACKGROUND

Genomes are vital to the study of genomics, population genetics, and evolution of species. To date, only one genome (Echinochloa crus-galli) for C4 annual weedy grass species has been sequenced. Research was conducted to develop a draft genome of goosegrass (Eleusine indica; 2n = 2x = 18), one of the most common and troublesome weeds in the world.

RESULTS

A draft assembly of an approximately 492 Mb whole-genome sequence of goosegrass was obtained by de novo assembly of paired-end and mate-paired reads generated by Illumina sequencing of total genomic DNA. The genome was assembled into 24,072 scaffolds with N50 = 233,459 bp. More than 99% of transcriptome sequences were mapped to the goosegrass draft genome, and 95% of the commonly conserved plant genes were present. The assembled genome contains 25,467 unique protein-coding genes. Genes associated with herbicide resistance were obtained and variant calling allowed the detection of 754,409 single nucleotide polymorphisms. In addition, we also report 115,417 simple sequence repeats which can be deployed in population genetics and phylogenetic analysis.

CONCLUSION

This is the first report of genome sequence of goosegrass. Our assembly was able to identify all major herbicide-resistance related genes and develop a useful tool for other genomic and evolutionary analysis. © 2019 Society of Chemical Industry.

Evaluating the evolution and function of the dynamic Venom Y protein in ectoparasitoid wasps

Venom of the parasitoid wasp Nasonia vitripennis changes the metabolism and gene expression in its fly host Sarcophaga bullata to induce developmental arrest, suppression of the immune response and various other venom effects. Yet, the venom of ectoparasitoid wasps has not been fully characterized. A major component of N. vitripennis venom is an uncharacterized, high-expressing protein referred to as Venom Y. Here we describe the evolutionary history and possible functions of this venom protein. We found that Venom Y is a relatively young gene that has duplicated to form two distinct paralogue groups. A copy of Venom Y has been recruited as a venom protein in at least five wasp species. Functional analysis found that Venom Y affects detoxification and immunity genes in envenomated fly hosts. Many of these genes are fat-body specific, suggesting that Venom Y may have a targeted effect on fat body tissue. We also show that Venom Y may mitigate negative effects of other venom proteins. Finally, protein sequencing indicates that Venom Y is post-translationally modified. This study contributes to elucidating parasitoid venom by using RNA interference knockdown to investigate venom protein function in the context of the whole venom cocktail.

Expression Analysis Reveals the Association of Several Genes with Pupal Diapause in Bactrocera minax (Diptera: Tephritidae)

The Chinese citrus fly, Bactrocera minax, is a devastating pest of citrus, which enters the obligatory diapause in overwintering pupae to resist harsh environmental conditions. However, little is known about the molecular mechanisms underlying pupal diapause. The previous transcriptomic analysis revealed that a large number of genes were regulated throughout the pupal stage. Of these genes, 12 and six ones that are remarkably up- and downregulated, respectively, specifically in intense diapause were manually screened out in present study. To validate the expression of these genes throughout the pupal stage, the quantitative real-time PCR (qRT-PCR) was conducted, and the genes displaying different expression patterns with those of previous study were excluded. Then, the expressions of remaining genes were compared between diapause-destined and non-diapause-destined pupae to reveal their association with diapause using qRT-PCR and semiquantitative PCR. Finally, five genes, TTLL3B, Cyp6a9, MSTA, Fru, and UC2, and two genes, KSPI and LYZ1, were demonstrated to be positively and negatively associated with diapause, respectively. These findings provide a solid foundation for the further investigation of molecular mechanisms underlying B. minax pupal diapause.

Global transcriptome profiling and functional analysis reveal that tissue-specific constitutive overexpression of cytochrome P450s confers tolerance to imidacloprid in palm weevils in date palm fields

Background

Cytochrome P450-dependent monooxygenases (P450s), constituting one of the largest and oldest gene superfamilies found in many organisms from bacteria to humans, play a vital role in the detoxification and inactivation of endogenous toxic compounds. The use of various insecticides has increased over the last two decades, and insects have developed resistance to most of these compounds through the detoxifying function of P450s. In this study, we focused on the red palm weevil (RPW), Rhynchophorus ferrugineus, the most devastating pest of palm trees worldwide, and demonstrated through functional analysis that upregulation of P450 gene expression has evolved as an adaptation to insecticide stress arising from exposure to the neonicotinoid-class systematic insecticide imidacloprid.

Results

Based on the RPW global transcriptome analysis, we identified 101 putative P450 genes, including 77 likely encoding protein coding genes with ubiquitous expression. A phylogenetic analysis revealed extensive functional and species-specific diversification of RPW P450s, indicating that multiple CYPs actively participated in the detoxification process. We identified highly conserved paralogs of insect P450s that likely play a role in the development of resistance to imidacloprid: Drosophila Cyp6g1 (CYP6345J1) and Bemisia tabaci CYP4C64 (CYP4LE1). We performed a toxicity bioassay and evaluated the induction of P450s, followed by the identification of overexpressed P450s, including CYP9Z82, CYP6fra5, CYP6NR1, CYP6345J1 and CYP4BD4, which confer cross-resistance to imidacloprid. In addition, under imidacloprid insecticide stress in a date palm field, we observed increased expression of various P450 genes, with CYP9Z82, CYP4BD4, CYP6NR1 and CYP6345J1 being the most upregulated detoxification genes in RPWs. Expression profiling and cluster analysis revealed P450 genes with multiple patterns of induction and differential expression. Furthermore, we used RNA interference to knock down the overexpressed P450s, after which a toxicity bioassay and quantitative expression analysis revealed likely candidates involved in metabolic resistance against imidacloprid in RPW. Ingestion of double-stranded RNA (dsRNA) successfully knocked down the expression of CYP9Z82, CYP6NR1 and CYP345J1 and demonstrated that silencing of CYP345J1 and CYP6NR1 significantly decreased the survival rate of adult RPWs treated with imidacloprid, indicating that overexpression of these two P450s may play an important role in developing tolerance to imidacloprid in a date palm field.

Conclusion

Our study provides useful background information on imidacloprid-specific induction and overexpression of P450s, which may enable the development of diagnostic tools/markers for monitoring the spread of insecticide resistant RPWs. The observed trend of increasing tolerance to imidacloprid in the date palm field therefore indicated that strategies for resistance management are urgently needed.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5837-4) contains supplementary material, which is available to authorized users.

Structural Variants and Selective Sweep Foci Contribute to Insecticide Resistance in the Drosophila Genetic Reference Panel

Patterns of nucleotide polymorphism within populations of Drosophila melanogaster suggest that insecticides have been the selective agents driving the strongest recent bouts of positive selection. However, there is a need to explicitly link selective sweeps to the particular insecticide phenotypes that could plausibly account for the drastic selective responses that are observed in these non-target insects. Here, we screen the Drosophila Genetic Reference Panel with two common insecticides; malathion (an organophosphate) and permethrin (a pyrethroid). Genome-wide association studies map survival on malathion to two of the largest sweeps in the D. melanogaster genome; Ace and Cyp6g1. Malathion survivorship also correlates with lines which have high levels of Cyp12d1, Jheh1 and Jheh2 transcript abundance. Permethrin phenotypes map to the largest cluster of P450 genes in the Drosophila genome, however in contrast to a selective sweep driven by insecticide use, the derived allele seems to be associated with susceptibility. These results underscore previous findings that highlight the importance of structural variation to insecticide phenotypes: Cyp6g1 exhibits copy number variation and transposable element insertions, Cyp12d1 is tandemly duplicated, the Jheh loci are associated with a Bari1 transposable element insertion, and a Cyp6a17 deletion is associated with susceptibility.

Identification and interaction of multiple genes resulting in DDT resistance in the 91-R strain of Drosophila melanogaster by RNAi approaches

4,4'-dichlorodiphenyltrichloroethane (DDT) has been re-recommended by the World Health Organization for malaria mosquito control. Previous DDT use has resulted in resistance, and with continued use resistance will likely increase in terms of level and extent. Drosophila melanogaster is a model dipteran with a well annotated genome allowing both forward and reverse genetic manipulation, numerous studies done on insecticide resistance mechanisms, and is related to malaria mosquitoes allowing for extrapolation. The 91-R strain of D. melanogaster is highly resistant to DDT (>1500-fold) and recently, reduced penetration, increased detoxification, and direct excretion have been identified as resistance mechanisms. Their interactions, however, remain unclear. Use of Gal4/UAS-RNAi transgenic lines of D. melanogaster allowed for the targeted knockdown of genes putatively involved in DDT resistance and has identified the role of several cuticular proteins (Cyp4g1 and Lcp1), cytochrome P450 monooxygenases (Cyp6g1 and Cyp12d1), and ATP binding cassette transporters (Mdr50, Mdr65, and Mrp1) involved in decreased sensitivity to DDT. These above findings have been further validated in 91-R flies using a nanoparticle-enhanced RNAi strategy, directly implication these genes in DDT resistance in 91-R flies.

Evidence for activation of nitenpyram by a mitochondrial cytochrome P450 in Drosophila melanogaster

BACKGROUND

Nitenpyram is a member of the economically important neonicotinoid class of insecticides. The in vivo metabolism of nitenpyram is not well characterised, but cytochrome P450 activity is the major mechanism of resistance to neonicotinoids identified in insect pests, and P450s metabolise other neonicotinoids including imidacloprid.

RESULTS

Here, we used the GAL4-UAS targeted expression system to direct RNA interference (RNAi) against the cytochrome P450 redox partners to interrupt P450 functions in specific tissues in Drosophila melanogaster. RNAi of the mitochondrial redox partner defective in the avoidance of repellents (dare) in the digestive tissues reduced nitenpyram mortality, suggesting an activation step in the metabolism of nitenpyram carried out by a mitochondrial P450. RNAi of the mitochondrial cytochrome P450 Cyp12a5, which is expressed in the digestive tissues, resulted in the same phenotype, and transgenic overexpression of Cyp12a5 increased nitenpyram sensitivity.

CONCLUSION

These results suggest that in vivo metabolism of nitenpyram by the mitochondrial P450 CYP12A5 results in the formation of a product with higher toxicity than the parent compound. © 2018 Society of Chemical Industry.

An overview of functional genomic tools in deciphering insecticide resistance

In this short review, we highlight three functional genomic technologies that have recently been contributing to the understanding of the molecular mechanisms underpinning insecticide resistance: the GAL4/UAS system, a molecular tool used to express genes of interest in a spatiotemporal controlled manner; the RNAi system, which is used to knock-down gene expression; and the most recently developed gene editing tool, CRISPR/Cas9, which can be used to knock-out and knock-in sequences of interest.

Baseline susceptibility of spotted wing Drosophila (Drosophila suzukii) to four key insecticide classes

BACKGROUND

The invasive drosophilid pest, Drosophila suzukii Matsumura, is affecting berry production in most fruit-producing regions of the world. Chemical control is the dominant management approach, creating concern for insecticide resistance in this pest. We compared the insecticide susceptibility of D. suzukii populations collected from conventional, organic or insecticide-free blueberry sites.

RESULTS

The sensitivity of D. suzukii to malathion and spinetoram declined slightly across the 3 years of monitoring, whereas it was more consistent for methomyl and zeta-cypermethrin. The sensitivity of D. suzukii to all four insecticides (LC₅₀ and LC₉₀ values) did not differ significantly among the blueberry fields using different management practices.

CONCLUSIONS

The baseline sensitivity of D. suzukii has been characterized, allowing future comparisons if field failures of chemical control are reported. The concentration achieving high control indicates that effective levels of control can still be achieved with field rates of these four insecticides. However, declining susceptibility of some populations of D. suzukii to some key insecticides highlights the need for resistance monitoring. © 2017 Society of Chemical Industry.

Genetic Diversity in Insect Metal Tolerance

Insects encounter a variety of metals in their environment, many of which are required at some concentration for normal organismal homeostasis, but essentially all of which are toxic at higher concentrations. Insects have evolved a variety of genetic, and likely epigenetic, mechanisms to deal with metal stress. A recurring theme in all these systems is complexity and diversity; even simple, single gene, cases are complex. Of the known gene families, the metallothioneins are perhaps the best understood and provide good examples of how diverse metal response is. Interestingly, there is considerable diversity across taxa in these metal-responsive systems, including duplications to form small gene families and complex expression of single loci. Strikingly, different species have evolved different mechanisms to cope with the same, or similar, stress suggesting both independent derivation of, and plasticity in, the pathways involved. It is likely that some metal-response systems evolved early in evolutionary time and have been conserved, while others have diverged, and still others evolved more recently and convergently. In addition to conventional genetics, insects likely respond to environmental metal through a variety of epigenetic systems, but direct tests are lacking. Ultimately, it is likely that classical genetic and epigenetic factors interact in regulating insect metal responses. In light of this diversity across species, future studies including a broad-based examination of gene expression in non-model species in complex environments will likely uncover additional genes and genetic and epigenetic mechanisms.

Cap n collar transcription factor regulates multiple genes coding for proteins involved in insecticide detoxification in the red flour beetle, Tribolium castaneum

In invertebrates, a heterodimer of xenobiotic transcription factors, cap n collar C isoform (CncC) and muscle aponeurosis fibromatosis (Maf) mediate cellular defense. In insects, these proteins regulate expression of genes involved in insecticide detoxification. In the current study, we performed sequencing of cDNA copied from RNA isolated from Tribolium castaneum pyrethroid resistant strain (QTC279) beetles injected with CncC or green fluorescence protein (GFP, control) dsRNA. Differential expression analysis of sequences identified 662 genes that showed a decrease and 91 genes that showed an increase in expression (p value ≤ 0.01 and log2 fold change of ≥ 1.5) in CncC knockdown insects when compared to their expression in control insects. We selected a subset of 27 downregulated genes and verified their differential expression using qRT-PCR. This subset of 27 genes included 21 genes with a predicted function in xenobiotic detoxification. RNAi and insecticide bioassays were employed to study the function of six of these genes coding for CYP4G7, CYP4G14, GST-1 and four ABC transporters, ABCA-UB, ABCA-A1 and ABCA-A1L and ABCA-9B involved in all three phases of insecticide detoxification. These data suggest that CncC regulates genes coding for proteins involved in detoxification of insecticides.

Resistance evolution in Drosophila: the case of CYP6G1

The massive use of DDT as an insecticide between 1940 and 1970 has resulted in the emergence of a resistant population of insects. One of the main metabolic mechanisms developed by resistant insects involves detoxification enzymes such as cytochrome P450s. These enzymes can metabolise the insecticide to render it less toxic and facilitate its elimination from the organism. The P450 Cyp6g1 was identified as the major factor responsible for DDT resistance in Drosophila melanogaster field populations. In this article, we review the data available for this gene since it was associated with resistance in 2002. The knowledge gained on Cyp6g1 allows a better understanding of the evolution of insecticide resistance mechanisms and highlights the major role of transposable elements in evolutionary processes. © 2016 Society of Chemical Industry.

Intralocus sexual conflict and insecticide resistance

The BA allele of the Drosophila cytochrome P450 gene Cyp6g1 confers resistance to a range of insecticides. It is also subject to intralocus sexual conflict when introgressed into the Canton-S background, whose collection predates the widespread use of insecticides. In this genetic background, the allele confers a pleiotropic fitness benefit to females but a cost to males, and exhibits little sexual dimorphism in conferred insecticide resistance. It is unclear whether these sexually antagonistic effects also exist in current populations that have naturally evolved with insecticides, where genetic modifiers that offset male costs might be expected to evolve. Here, we explore these issues using Drosophila melanogaster caught recently from an Australian population in which the BA allele naturally segregates. While we find increased fecundity in insecticide-resistant BA females and no consistent evidence of fitness costs in males, experimental evolution indicates balancing selection at the locus. We suggest that this apparent discrepancy may be due to reduced investment in reproduction in resistant males. Our results at the population level are consistent with previous work, and suggest that individual-level fitness assays do not always capture sexually antagonistic fitness effects that emerge in a population context.

Identification of transcripts involved in digestion, detoxification and immune response from transcriptome of Empoasca vitis (Hemiptera: Cicadellidae) nymphs

Tea production has been significantly impacted by the false-eye leafhopper, Empoasca vitis (Göthe), around Asia. To identify the key genes which are responsible for nutrition absorption, xenobiotic metabolism and immune response, the transcriptome of either alimentary tracts or bodies minus alimentary tract of E. vitis was sequenced and analyzed. Over 31 million reads were obtained from Illumina sequencing. De novo sequence assembly resulted in 52,182 unigenes with a mean size of 848nt. The assembled unigenes were then annotated using various databases. Transcripts of at least 566 digestion-, 224 detoxification-, and 288 immune-related putative genes in E. vitis were identified. In addition, relative expression of highly abundant transcripts was verified through quantitative real-time PCR. Results from this investigation provide genomic information about E. vitis, which will be helpful in further study of E. vitis biology and in the development of novel strategies to control this devastating pest.

Trends and Challenges in Pesticide Resistance Detection

Pesticide resistance is a crucial factor to be considered when developing strategies for the minimal use of pesticides while maintaining pesticide efficacy. This goal requires monitoring the emergence and development of resistance to pesticides in crop pests. To this end, various methods for resistance diagnosis have been developed for different groups of pests. This review provides an overview of biological, biochemical, and molecular methods that are currently used to detect and quantify pesticide resistance. The agronomic, technical, and economic advantages and drawbacks of each method are considered. Emerging technologies are also described, with their associated challenges and their potential for the detection of resistance mechanisms likely to be selected by current and future plant protection methods.

Collaborative contribution of six cytochrome P450 monooxygenase genes to fenpropathrin resistance in Tetranychus cinnabarinus (Boisduval)

Cytochrome P450 monooxygenases (P450s), as an important family of detoxification enzymes, participate in the metabolism of agrochemicals in almost all agricultural pests and play important roles in the development of insecticide resistance. Two P450 genes (CYP389B1 and CYP392A26) were identified and their expression patterns were investigated in our previous study. In this study, four more P450 gene sequences (CYP391A1, CYP384A1, CYP392D11 and CYP392A28) from the Clan 2, Clan 3 and Clan 4 families were identified and characterized. Quantitative PCR analysis showed that these four P450 genes were highly expressed in a fenpropathrin-resistant (FeR) strain of Tetranychus cinnabarinus. In addition, their expressions were much more sensitive to fenpropathrin induction in the FeR strain than the susceptible strain. Gene-silencing experiments via double-stranded RNA feeding were carried out. The results showed that mRNA levels of these six P450 genes were reduced in the FeR strain and the activities of P450s were decreased. Consequently mite susceptibilities to fenpropathrin were increased. Interestingly, silencing all six P450 genes simultaneously had an even greater effect on resistance than silencing them individually. This study increases our understanding of the molecular mechanisms of insecticide detoxification, suggesting that the overexpression of these six P450 genes might play important roles in fenpropathrin resistance in T. cinnabarinus collaboratively.

Over-expression of multiple cytochrome P450 genes in fenvalerate-resistant field strains of Helicoverpa armigera from north of China

Pyrethroid resistance was one of the main reasons for control failure of cotton bollworm Helicoverpa armigera (Hübner) in China. The promotion of Bt crops decreased the application of chemical insecticides in controlling H.armigera. However, the cotton bollworm still kept high levels of resistance to fenvalerate. In this study, the resistance levels of 8 field-collected strains of H. armigera from north of China to 4 insecticides, as well as the expression levels of related P450 genes were investigated. The results of bioassay indicated that the resistance levels to fenvalerate in the field strains varied from 5.4- to 114.7-fold, while the resistance levels to lambda-cyhalothrin, phoxim and methomyl were low, which were ranged from 1.5- to 5.2-, 0.2- to 1.6-, and 2.9- to 8.3- fold, respectively, compared to a susceptible strain. Synergistic experiment showed that PBO was the most effective synergist in increasing the sensitivity of H. armigera to fenvalerate, suggesting that P450 enzymes were involved in the pyrethroid resistance in the field strains. The results of quantitative RT-PCR indicated that eight P450 genes (CYP332A1, CYP4L11, CYP4L5, CYP4M6, CYP4M7, CYP6B7, CYP9A12, CYP9A14) were all significantly overexpressed in Hejian1 and Xiajin1 strains of H. armigera collected in 2013, and CYP4L5 was significantly overexpressed in all the 6 field strains collected in 2014. CYP332A1, CYP6B7 and CYP9A12 had very high overexpression levels in all the field strains, indicating their important roles in fenvalerate resistance. The results suggested that multiple P450 genes were involved in the high-level fenvalerate-resistance in different field strains of H. armigera collected from north of China.

Transcriptional Differences between Diapausing and Non-Diapausing D. montana Females Reared under the Same Photoperiod and Temperature

BACKGROUND

A wide range of insects living at higher latitudes enter diapause at the end of the warm season, which increases their chances of survival through harsh winter conditions. In this study we used RNA sequencing to identify genes involved in adult reproductive diapause in a northern fly species, Drosophila montana. Both diapausing and non-diapausing flies were reared under a critical day length and temperature, where about half of the emerging females enter diapause enabling us to eliminate the effects of varying environmental conditions on gene expression patterns of the two types of female flies.

RESULTS

RNA sequencing revealed large differences between gene expression patterns of diapausing and non-diapausing females, especially in genes involved with metabolism, fatty acid biosynthesis, and metal and nucleotide binding. Differently expressed genes included several gene groups, including myosin, actin and cytochromeP450 genes, which have been previously associated with diapause. This study also identified new candidate genes, including some involved in cuticular hydrocarbon synthesis or regulation (desat1 and desat2), and acyl-CoA Δ11-desaturase activity (CG9747), and few odorant-binding protein genes (e.g. Obp44A). Also, several transposable elements (TEs) showed differential expression between the two female groups motivating future research on their roles in diapause.

CONCLUSIONS

Our results demonstrate that the adult reproductive diapause in D. montana involves changes in the expression level of a variety of genes involved in key processes (e.g. metabolism and fatty acid biosynthesis) which help diapausing females to cope with overwintering. This is consistent with the view that diapause is a complex adaptive phenotype where not only sexual maturation is arrested, but also changes in adult physiology are required in order to survive over the winter.

Mechanisms of TiO2 NPs-induced phoxim metabolism in silkworm (Bombyx mori) fat body

Silkworm is an important economic insect. Abuse of organophosphorus pesticides in recent years often leads to poisoning of silkworms, which significantly affects sericulture development by reducing silk production. Previous studies have shown that TiO2 NPs can effectively mitigate the damages caused by organophosphorus pesticides in silk glands and nerve tissues. The fat body is an important metabolic detoxification organ of silkworms, but it is unknown whether TiO2 NPs affect pesticide metabolism in fat body. In this study, we characterized the transcription of antioxidant genes and enzyme activity in fat body after TiO2 NPs and phoxim treatments using transcriptome sequencing, real-time PCR, and enzyme activity assay. Transcriptome sequencing detected 10 720, 10 641, 10 403, and 10 489 genes for control group, TiO2 NPs group, phoxim group, and TiO2 NPs+phoxim group, respectively. The TiO2 NPs+phoxim group had 705 genes with significantly differential expression (FDR<0.001), among which the antioxidant genes thioredoxin reductase 1 and glutathione S-transferase omega 3 were significantly upregulated. In phoxim group, the expression levels of superoxide dismutase (SOD), catalase (CAT), glutathione S-transferase delta (GSTd), and thioredoxin peroxidase (TPx) were increased by 1.365 -fold, 1.335 -fold, 1.642 -fold, and 1.765 -fold, respectively. The level changes of SOD, CAT, GSTd, and TPx were validated by real time PCR. The contents of reactive oxygen species (ROS), malondialdehyde (MDA), and hydrogen peroxide (H2O2) were increased by 1.598 -fold, 1.946 -fold, and 1.506 -fold, respectively, indicating that TiO2 NPs treatment can relieve phoxim-induced oxidative stress. To clarify the mechanism of TiO2 NPs's effect, the transcription levels of P450 gene family were measured for the TiO2 NPs+phoxim group; the expression levels of CYP4M5, CYP6AB4, CYP6A8, and CYP9G3 were elevated by 2.784 -fold, 3.047 -fold, 2.254 -fold, and 4.253 -fold, respectively, suggesting that high expression of P450 family genes can enhance the metabolism of phoxim in the fat body. The results of this study indicated that TiO2 NPs treatment promoted the transcriptional expression of the P450 family genes to improve the fat body's ability to metabolize phoxim and reduce phoxim-induced oxidative stress. This may be the main mechanism of TiO2 NPs' mitigation of phoxim-induced damages in the fat body.

Splice form variant and amino acid changes in MDR49 confers DDT resistance in transgenic Drosophila

The ATP-binding cassette (ABC) transporters represent a superfamily of proteins that have important physiological roles in both prokaryotes and eukaryotes. In insects, ABC transporters have previously been implicated in insecticide resistance. The 91-R strain of Drosophila melanogaster has been intensely selected with DDT over six decades. A recent selective sweeps analysis of 91-R implicated the potential role of MDR49, an ABC transporter, in DDT resistance, however, to date the details of how MDR49 may play a role in resistance have not been elucidated. In this study, we investigated the impact of structural changes and an alternative splicing event in MDR49 on DDT-resistance in 91-R, as compared to the DDT susceptible strain 91-C. We observed three amino acid differences in MDR49 when 91-R was compared with 91-C, and only one isoform (MDR49B) was implicated in DDT resistance. A transgenic Drosophila strain containing the 91-R-MDR49B isoform had a significantly higher LD₅₀ value as compared to the 91-C-MDR49B isoform at the early time points (6 h to 12 h) during DDT exposure. Our data support the hypothesis that the MDR49B isoform, with three amino acid mutations, plays a role in the early aspects of DDT resistance in 91-R.

Quantitative trait locus mapping and functional genomics of an organophosphate resistance trait in the western corn rootworm, Diabrotica virgifera virgifera

The western corn rootworm, Diabrotica virgifera virgifera, is an insect pest of corn and population suppression with chemical insecticides is an important management tool. Traits conferring organophosphate insecticide resistance have increased in frequency amongst D. v. virgifera populations, resulting in the reduced efficacy in many corn-growing regions of the USA. We used comparative functional genomic and quantitative trait locus (QTL) mapping approaches to investigate the genetic basis of D. v. virgifera resistance to the organophosphate methyl-parathion. RNA from adult methyl-parathion resistant and susceptible adults was hybridized to 8331 microarray probes. The results predicted that 11 transcripts were significantly up-regulated in resistant phenotypes, with the most significant (fold increases ≥ 2.43) being an α-esterase-like transcript. Differential expression was validated only for the α-esterase (ST020027A20C03), with 11- to 13-fold greater expression in methyl-parathion resistant adults (P < 0.05). Progeny with a segregating methyl-parathion resistance trait were obtained from a reciprocal backcross design. QTL analyses of high-throughput single nucleotide polymorphism genotype data predicted involvement of a single genome interval. These data suggest that a specific carboyxesterase may function in field-evolved corn rootworm resistance to organophosphates, even though direct linkage between the QTL and this locus could not be established.

Functional characterization of NADPH-cytochrome P450 reductase from Bactrocera dorsalis: Possible involvement in susceptibility to malathion

NADPH cytochrome P450 reductase (CPR) is essential for cytochrome P450 catalysis, which is important in the detoxification and activation of xenobiotics. In this study, two transcripts of Bactrocera dorsalis CPR (BdCPR) were cloned, and the deduced amino-acid sequence had an N-terminus membrane anchor for BdCPR-X1 and three conserved binding domains (FMN, FAD, and NADP), as well as an FAD binding motif and catalytic residues for both BdCPR-X1 and BdCPR-X2. BdCPR-X1 was detected to have the high expression levels in adults and in Malpighian tubules, fat bodies, and midguts of adults, but BdCPR-X2 expressed lowly in B. dorsalis. The levels of BdCPRs were similar in malathion-resistant strain compared to susceptible strain. However, injecting adults with double-stranded RNA against BdCPR significantly reduced the transcript levels of the mRNA, and knockdown of BdCPR increased adult susceptibility to malathion. Expressing complete BdCPR-X1 cDNA in Sf9 cells resulted in high activity determined by cytochrome c reduction and these cells had higher viability after exposure to malathion than control. The results suggest that BdCPR could affect the susceptibility of B. dorsalis to malathion and eukaryotic expression of BdCPR would lay a solid foundation for further investigation of P450 in B. dorsalis.

Transcription analysis of neonicotinoid resistance in Mediterranean (MED) populations of B. tabaci reveal novel cytochrome P450s, but no nAChR mutations associated with the phenotype

Background

Bemisia tabaci is one of the most damaging agricultural pests world-wide. Although its control is based on insecticides, B. tabaci has developed resistance against almost all classes of insecticides, including neonicotinoids.

Results

We employed an RNA-seq approach to generate genome wide expression data and identify genes associated with neonicotinoid resistance in Mediterranean (MED) B. tabaci (Q1 biotype). Twelve libraries from insecticide resistant and susceptible whitefly populations were sequenced on an Illumina Next-generation sequencing platform, and genomic sequence information of approximately 73 Gbp was generated.

A reference transcriptome was built by de novo assembly and functionally annotated. A total of 146 P450s, 18 GSTs and 23 CCEs enzymes (unigenes) potentially involved in the detoxification of xenobiotics were identified, along with 78 contigs encoding putative target proteins of six different insecticide classes. Ten unigenes encoding nicotinic Acetylcholine Receptors (nAChR), the target of neoinicotinoids, were identified and phylogenetically classified. No nAChR polymorphism potentially related with the resistant phenotypes, was observed among the studied strains.

DE analysis revealed that among the 550 differentially (logFC > 1) over-transcribed unigenes, 52 detoxification enzymes were over expressed including unigenes with orthologues in P450s, GSTs, CCE and UDP-glucuronosyltransferases.

Eight P450 unigenes belonging to clades CYP2, CYP3 and CYP4 were highly up-regulated (logFC > 2) including CYP6CM1, a gene already known to confer imidacloprid resistance in B. tabaci. Using quantitative qPCRs, a larger screening of field MED B. tabaci from Crete with known neonicotinoid phenotype was performed to associate expression levels of P450s with resistance levels. Expression levels of five P450s, including CYP6CM1, were found associated with neonicotinoid resistance. However, a significant correlation was found only in CYP303 and CYP6CX3, with imidacloprid and acetamiprid respectively.

Conclusion

Our work has generated new toxicological data and genomic resources which will significantly enrich the available dataset and substantially facilitate the molecular studies in MED B. tabaci. No evidence of target site neonicotinoid resistance has been found. Eight P450 unigenes, including CYP6CM1, were found significantly over-expressed in resistant B. tabaci. This study suggests at least two novel P450s (CYP303 and CYP6CX3) as candidates for their functional characterization as detoxification mechanisms of neonicotinoid resistance in B. tabaci.

Electronic supplementary material

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

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

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

Expression profile analysis of silkworm P450 family genes after phoxim induction

Silkworm (Bombyx mori) is an important economic insect and a model species for Lepidopteran. Each year, O,O-diethyl O-(alpha-cyanobenzylideneamino) phosphorothioate (phoxim) pesticide poisoning in China results in huge economic losses in sericulture. Silkworm fat body is the main organ for nutrient storage, energy supply, intermediary metabolism, and detoxification. Microarray analysis of silkworm Cytochrome P450 detoxification enzyme genes revealed that all tested P450 4 (CYP4) family genes are expressed in the fat body. Quantitative Real-time PCR (QRT-PCR) was used to detect the expression of CYP4 family genes in silkworm fat body 0, 24, 48, and 72 h after phoxim exposure. The expression levels of silkworm molting hormone synthesis-related genes started to change 24 h after phoxim exposure, with those of CYP302A1, CYP306A1, and CYP314A1 being elevated by 1.38-, 1.33-, and 2.10-fold, respectively. The CYP18A1 gene that participates in steroid hormone inactivation and the CYP15C1 gene that participates in the epoxidation during the synthesis of juvenile hormone (JH) from methyl farnesoate (MF) were increased by 3.85- and 7.82-fold, respectively. Phylogenetic analysis indicated that these endogenous hormone metabolism-related genes belong to CYP mito clan and clan 2, and that phoxim exposure may affect silkworm development and metamorphosis. The CYP4, CYP6, and CYP9 families all showed some degrees of increases in gene expression; among them, CYP49A1, CYP4L6, CYP6AB4, CYP9G3, CYP9A19, and CYP9A22's transcription levels were significantly upregulated to 12.77-, 2.64-, 2.42-, 4.06-, 3.32-, and 2.98-fold, respectively, of the control levels. In the fat body, CYP49A1, CYP6AB4, CYP9A19, and CYP9A22 were constantly expressed at high levels after 24, 48, and 72 h of phoxim treatments; according to phylogenetic analysis, these genes belong to detoxification-related clan 3 and clan 4 CYP families. These genes may participate in the metabolism of phoxim in silkworm fat body. The results obtained in this study provide a basis for future in-depth investigations of insect P450 family genes in metabolic detoxification.

Cytochrome P450 CYP4DE1 and CYP6CW3v2 contribute to ethiprole resistance in Laodelphax striatellus (Fallén)

Laodelphax striatellus Fallén (Hemiptera: Delphacidae), a destructive pest of rice, has developed high resistance to multiple insecticides, threatening the success of pest management programmes. The present study investigated ethiprole resistance mechanisms in a field population that is highly resistant to ethiprole. That population was used to establish a laboratory population that was subjected to further selection to produce a resistant strain. Target genes were cloned and compared between the resistant and the susceptible strains, the role of detoxification enzymes was examined, and the relative expression levels of 71 detoxification enzyme genes were tested using quantitative real time (RT)-PCR. The laboratory selection enhanced the resistance from 107-fold to 180-fold. The Rdl-type target site mutation seldom occurred in the resistant strain and is unlikely to represent the major mechanism underlying the observed resistance. Of the three important detoxification enzymes, only P450 monooxygenase was found to be associated with ethiprole resistance. Moreover, two genes, CYP4DE1 and CYP6CW3v2, were found to be overexpressed in the resistant strain. Furthermore, gene-silencing via a double-stranded RNA feeding test was carried out, and the results showed that the mRNA levels of CYP4DE1 and CYP6CW3v2 were reduced in the resistant strain, whereas ethiprole susceptibility was increased. These results suggest that CYP4DE1 and CYP6CW3v2 play an important role in ethiprole resistance in L. striatellus.

RNAi validation of resistance genes and their interactions in the highly DDT-resistant 91-R strain of Drosophila melanogaster

4,4'-dichlorodiphenyltrichloroethane (DDT) has been re-recommended by the World Health Organization for malaria mosquito control. Previous DDT use has resulted in resistance, and with continued use resistance will increase in terms of level and extent. Drosophila melanogaster is a model dipteran that has many available genetic tools, numerous studies done on insecticide resistance mechanisms, and is related to malaria mosquitoes allowing for extrapolation. The 91-R strain of D. melanogaster is highly resistant to DDT (>1500-fold), however, there is no mechanistic scheme that accounts for this level of resistance. Recently, reduced penetration, increased detoxification, and direct excretion have been identified as resistance mechanisms in the 91-R strain. Their interactions, however, remain unclear. Use of UAS-RNAi transgenic lines of D. melanogaster allowed for the targeted knockdown of genes putatively involved in DDT resistance and has validated the role of several cuticular proteins (Cyp4g1 and Lcp1), cytochrome P450 monooxygenases (Cyp6g1 and Cyp12d1), and ATP binding cassette transporters (Mdr50, Mdr65, and Mrp1) involved in DDT resistance. Further, increased sensitivity to DDT in the 91-R strain after intra-abdominal dsRNA injection for Mdr50, Mdr65, and Mrp1 was determined by a DDT contact bioassay, directly implicating these genes in DDT efflux and resistance.

Selective sweep analysis in the genomes of the 91-R and 91-C Drosophila melanogaster strains reveals few of the 'usual suspects' in dichlorodiphenyltrichloroethane (DDT) resistance

Adaptation of insect phenotypes for survival after exposure to xenobiotics can result from selection at multiple loci with additive genetic effects. To the authors' knowledge, no selective sweep analysis has been performed to identify such loci in highly dichlorodiphenyltrichloroethane (DDT) resistant insects. Here we compared a highly DDT resistant phenotype in the Drosophila melanogaster (Drosophila) 91-R strain to the DDT susceptible 91-C strain, both of common origin. Whole genome re-sequencing data from pools of individuals was generated separately for 91-R and 91-C, and mapped to the reference Drosophila genome assembly (v. 5.72). Thirteen major and three minor effect chromosome intervals with reduced nucleotide diversity (π) were identified only in the 91-R population. Estimates of Tajima's D (D) showed corresponding evidence of directional selection in these same genome regions of 91-R, however, no similar reductions in π or D estimates were detected in 91-C. An overabundance of non-synonymous proteins coding to synonymous changes were identified in putative open reading frames associated with 91-R. Except for NinaC and Cyp4g1, none of the identified genes were the 'usual suspects' previously observed to be associated with DDT resistance. Additionally, up-regulated ATP-binding cassette transporters have been previously associated with DDT resistance; however, here we identified a structurally altered MDR49 candidate resistance gene. The remaining fourteen genes have not previously been shown to be associated with DDT resistance. These results suggest hitherto unknown mechanisms of DDT resistance, most of which have been overlooked in previous transcriptional studies, with some genes having orthologs in mammals.

A glycine insertion in the estrogen-related receptor (ERR) is associated with enhanced expression of three cytochrome P450 genes in transgenic Drosophila melanogaster

Insecticide-resistant Drosophila melanogaster strains represent a resource for the discovery of the underlying molecular mechanisms of cytochrome P450 constitutive over-expression, even if some of these P450s are not directly involved in the resistance phenotype. For example, in select 4,4'-dichlorodiphenyltrichloroethane (DDT) resistant strains the glucocorticoid receptor-like (GR-like) potential transcription factor binding motifs (TFBMs) have previously been shown to be associated with constitutively differentially-expressed cytochrome P450s, Cyp12d1, Cyp6g2 and Cyp9c1. However, insects are not known to have glucocorticoids. The only ortholog to the mammalian glucocorticoid receptor (GR) in D. melanogaster is an estrogen-related receptor (ERR) gene, which has two predicted alternative splice isoforms (ERRa and ERRb). Sequencing of ERRa and ERRb in select DDT susceptible and resistant D. melanogaster strains has revealed a glycine (G) codon insertion which was only observed in the ligand binding domain of ERR from the resistant strains tested (ERR-G). Transgenic flies, expressing the ERRa-G allele, constitutively over-expressed Cyp12d1, Cyp6g2 and Cyp9c1. Only Cyp12d1 and Cyp6g2 were over-expressed in the ERRb-G transgenic flies. Phylogenetic studies show that the G-insertion appeared to be located in a less conserved domain in ERR and this insertion is found in multiple species across the Sophophora subgenera.

Characterization of the transcriptome of the Asian gypsy moth Lymantria dispar identifies numerous transcripts associated with insecticide resistance

Although the Asian gypsy moth Lymantria dispar causes extensive forest damage worldwide, little is known regarding the genes involved in its development or response to insecticides. Accordingly, characterization of the transcriptome of L. dispar larvae would promote the development of toxicological methods for its control. RNA-seq analysis of L. dispar larvae messenger RNA (mRNA) generated 62,063 unigenes with N50 of 993 bp, from which 23,975 unique sequences (E-value < 10(-5)) were identified using a BLASTx search of the NCBI non-redundant (nr) database. Using functional classification in the Gene Ontology (GO) and Clusters of Orthologous Groups (COG) databases, 7,309 indentified sequences were categorized into 51 functional groups and 8,079 sequences were categorized into 25 functional groups, respectively. Moreover, we identified a large number of transcripts encoding known insecticide targets, or proteins involved in the metabolism of insecticides. Reads per kilobase of unigene length per million mapped reads (RPKM) analysis identified 39 high abundance transcripts, of which 27 exhibited significantly altered expression patterns across the egg, larvae, pupae, male and female adult stages. Our study provides the most comprehensive transcriptomic sequence resource for L. dispar, which will form the basis for future identification of candidate insecticide resistance genes in L. dispar.

Identification of promoter polymorphisms in the cytochrome P450 CYP6AY1 linked with insecticide resistance in the brown planthopper, Nilaparvata lugens

Imidacloprid resistance in the brown planthopper, Nilaparvata lugens, is primarily the result of the over-expression of cytochrome P450 monooxygenases. Here, a field-collected strain of N. lugens was shown to be highly resistant to both imidacloprid and buprofezin. Insecticide exposure and quantitative real-time PCR revealed that its resistance was mainly associated with a cytochrome P450 gene, CYP6AY1. CYP6AY1 is known to metabolize imidacloprid but its effect on buprofezin is unclear. In the 5'-untranslated region of CYP6AY1, a novel alternative splicing was detected. After a 1990-bp promoter region was cloned, its basal luciferase activity was assessed. Furthermore, genotyping studies identified 12 variations in the promoter region that discriminated between the field-collected and control strain. Finally, survival bioassays revealed a single nucleotide polymorphism and an insertion-deletion polymorphism linked to buprofezin and imidacloprid resistance. Mutagenesis of these sites enhanced the promoter activity of CYP6AY1. These results suggest that promoter polymorphisms may affect P450-mediated multiple insecticide resistance of pests.

Identification of candidate odorant degrading gene/enzyme systems in the antennal transcriptome of Drosophila melanogaster

The metabolism of volatile signal molecules by odorant degrading enzymes (ODEs) is crucial to the ongoing sensitivity and specificity of chemoreception in various insects, and a few specific esterases, cytochrome P450s, glutathione S-transferases (GSTs) and UDP-glycosyltransferases (UGTs) have previously been implicated in this process. Significant progress has been made in characterizing ODEs in Lepidoptera but very little is known about them in Diptera, including in Drosophila melanogaster, a major insect model. We have therefore carried out a transcriptomic analysis of the antennae of D. melanogaster in order to identify candidate ODEs. Virgin male and female and mated female antennal transcriptomes were determined by RNAseq. As with the Lepidoptera, we found that many esterases, cytochrome P450 enzymes, GSTs and UGTs are expressed in D. melanogaster antennae. As olfactory genes generally show selective expression in the antennae, a comparison to previously published transcriptomes for other tissues has been performed, showing preferential expression in the antennae for one esterase, JHEdup, one cytochrome P450, CYP308a1, and one GST, GSTE4. These largely uncharacterized enzymes are now prime candidates for ODE functions. JHEdup was expressed heterologously and found to have high catalytic activity against a chemically diverse group of known ester odorants for this species. This is a finding consistent with an ODE although it might suggest a general role in clearing several odorants rather than a specific role in clearing a particular odorant. Our findings do not preclude the possibility of odorant degrading functions for other antennally expressed esterases, P450s, GSTs and UGTs but, if so, they suggest that these enzymes also have additional functions in other tissues.

Whole-genome expression analysis in the third instar larval midgut of Drosophila melanogaster

Survival of insects on a substrate containing toxic substances such as plant secondary metabolites or insecticides is dependent on the metabolism or excretion of those xenobiotics. The primary sites of xenobiotic metabolism are the midgut, Malpighian tubules, and fat body. In general, gene expression in these organs is reported for the entire tissue by online databases, but several studies have shown that gene expression within the midgut is compartmentalized. Here, RNA sequencing is used to investigate whole-genome expression in subsections of third instar larval midguts of Drosophila melanogaster. The data support functional diversification in subsections of the midgut. Analysis of the expression of gene families that are implicated in the metabolism of xenobiotics suggests that metabolism may not be uniform along the midgut. These data provide a starting point for investigating gene expression and xenobiotic metabolism and other functions of the larval midgut.