Differential expression and induction of two Drosophila cytochrome P450 genes near the Rst(2)DDT locus

Comparative Transcriptome Analysis to Reveal Differentially Expressed Cytochrome P450 in Response to Imidacloprid in the Aphid Lion, Chrysoperla zastrowi sillemi (Esben-Petersen)

The aphid lion, Chrysoperla zastrowi sillemi (Neuroptera: Chrysopidae) is a highly effective beneficial predator of many agricultural pests and has developed resistance to several insecticides. Understanding the molecular mechanism of insecticide resistance in the predators is crucial for its effective application in IPM programs. Therefore, transcriptomes of imidacloprid-resistant and susceptible strains have been assessed using RNA-seq. Cytochrome P450 is one of the important gene families involved in xenobiotic metabolism. Hence, our study focused on the CYP gene family where mining, nomenclature, and phylogenetic analysis revealed a total of 95 unique CYP genes with considerable expansion in CYP3 and CYP4 clans. Further, differential gene expression (DGE) analysis revealed ten CYP genes from CYP3 and CYP4 clans to be differentially expressed, out of which nine genes (CYP4419A1, CYP4XK1, CYP4416A10, CYP4416A-fragment8, CYP6YL1, CYP6YH6, CYP9GK-fragment16, CYP9GN2, CYP9GK6) were downregulated and one (CYP9GK3) was upregulated in the resistant strain as compared to the susceptible strain. Expression validation by quantitative real-time PCR (qRT-PCR) is consistent with the DGE results. The expansion and differential expression of CYP genes may be an indicator of the capacity of the predator to detoxify a particular group of insecticides.

CYP4CJ6-mediated resistance to two neonicotinoid insecticides in Sitobion miscanthi (Takahashi)

The wheat aphid Sitobion miscanthi (CWA) is an important harmful pest in wheat fields. Insecticide application is the main method to effectively control wheat aphids. However, CWA has developed resistance to some insecticides due to its extensive application, and understanding resistance mechanisms is crucial for the management of CWA. In our study, a new P450 gene, CYP4CJ6, was identified from CWA and showed a positive response to imidacloprid and thiamethoxam. Transcription of CYP4CJ6 was significantly induced by both imidacloprid and thiamethoxam, and overexpression of CYP4CJ6 in the imidacloprid-resistant strain was also observed. The sensitivity of CWA to these two insecticides was increased after the knockdown of CYP4CJ6. These results indicated that CYP4CJ6 could be associated with CWA resistance to imidacloprid and thiamethoxam. Subsequently, the posttranscriptional regulatory mechanism was assessed, and miR-316 was confirmed to participate in the posttranscriptional regulation of CYP4CJ6. These results are crucial for clarifying the roles of P450 in the resistance of CWA to insecticides.

Xenobiotic responses in insects

Insects have evolved a powerful detoxification system to protect themselves against environmental and anthropogenic xenobiotics including pesticides and nanoparticles. The resulting tolerance to insecticides is an immense problem in agriculture. In this study, we summarize advances in our understanding of insect xenobiotic responses: the detoxification strategies and the regulation mechanisms against xenobiotics including nanoparticles, the problem of response specificity and the potential usefulness of this study field for an elaborate pest management. In particular, we highlight that versatility of the detoxification system relies on the relatively unspecific recognition of a broad range of potential toxic substances that trigger either of various canonical xenobiotic responses signaling pathways, including CncC/Keap1, HR96, AHR/ARNT, GPCR, and MAPK/CREB. However, it has emerged that the actual response to an inducer may nevertheless be specific. There are two nonexclusive possibilities that may explain response specificity: (1) differential cross-talk between the known pathways and (2) additional, yet unidentified regulators and pathways of detoxification. Hence, a deeper and broader understanding of the regulation mechanisms of xenobiotic response in insects in the future might facilitate the development and application of highly efficient and environmentally friendly pest control methods, allowing us to face the challenge of the world population growth.

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.

Post-transcriptional modulation of cytochrome P450s, Cyp6g1 and Cyp6g2, by miR-310s cluster is associated with DDT-resistant Drosophila melanogaster strain 91-R

The role of miRNAs in mediating insecticide resistance remains largely unknown, even for the model species Drosophila melanogaster. Building on prior research, this study used microinjection of synthetic miR-310s mimics into DDT-resistant 91-R flies and observed both a significant transcriptional repression of computationally-predicted endogenous target P450 detoxification genes, Cyp6g1 and Cyp6g2, and also a concomitant increase in DDT susceptibility. Additionally, co-transfection of D. melanogaster S2 cells with dual luciferase reporter constructs validated predictions that miR-310s bind to target binding sites in the 3ʹ untranslated regions (3ʹ-UTR) of both Cyp6g1 and Cyp6g2 in vitro. Findings in the current study provide empirical evidence for a link between reduced miRNA expression and an insecticidal resistance phenotype through reduced targeted post-transcriptional suppression of transcripts encoding proteins involved in xenobiotic detoxification. These insights are important for understanding the breadth of adaptive molecular changes that have contributed to the evolution of DDT resistance in D. melanogaster.

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.

Diverse families of transposable elements affect the transcriptional regulation of stress-response genes in Drosophila melanogaster

Although transposable elements are an important source of regulatory variation, their genome-wide contribution to the transcriptional regulation of stress-response genes has not been studied yet. Stress is a major aspect of natural selection in the wild, leading to changes in the transcriptional regulation of a variety of genes that are often triggered by one or a few transcription factors. In this work, we take advantage of the wealth of information available for Drosophila melanogaster and humans to analyze the role of transposable elements in six stress regulatory networks: immune, hypoxia, oxidative, xenobiotic, heat shock, and heavy metal. We found that transposable elements were enriched for caudal, dorsal, HSF, and tango binding sites in D. melanogaster and for NFE2L2 binding sites in humans. Taking into account the D. melanogaster population frequencies of transposable elements with predicted binding motifs and/or binding sites, we showed that those containing three or more binding motifs/sites are more likely to be functional. For a representative subset of these TEs, we performed in vivo transgenic reporter assays in different stress conditions. Overall, our results showed that TEs are relevant contributors to the transcriptional regulation of stress-response genes.

MicroRNAs as Potential Biomarkers of Insecticide Exposure: A Review

Insecticides are key weapons for the control of pests. Large scale use of insecticides is harmful to the ecosystem, which is made up of a wide range of species and environments. MicroRNAs (miRNAs) are a class of endogenous single-stranded noncoding small RNAs in length of 20-24 nucleotides (nt), which extensively regulate expression of genes at transcriptional and post-transcriptional levels. The current research on miRNA-induced insecticide resistance reveals that dysregulated miRNAs cause significant changes in detoxification genes, particularly cytochrome P450s. Meanwhile, insecticide-induced changes in miRNAs are related to the decline of honeybees and threatened the development of zebrafish and other animals. Additionally, miRNAs are involved in insecticide-induced cytotoxicity, and dysregulated miRNAs are associated with human occupational and environmental exposure to insecticides. Therefore, miRNAs are valuable novel biomarkers of insecticide exposure, and they are potential factors to explain the toxicological effects of insecticides.

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.

Cytochrome P450s Cyp4p1 and Cyp4p2 associated with the DDT tolerance in the Drosophila melanogaster strain 91-R

Cytochrome P450s are part of a super-gene family that has undergone gene duplication, divergence, over-expression and, in some cases, loss of function. One such case is the 91-R and 91-C strains of common origin, in Drosophila melanogaster, whereby 91-R (DDT resistant strain) overexpresses Cyp4p1 and Cyp4p2 and both genes are lost in 91-C (DDT susceptible strain). In this study, we used a comparative approach to demonstrate that transcription of Cyp4p1 and Cyp4p2 were constitutively up-regulated in the Drosophila melanogaster strain 91-R as compared to another DDT susceptible strain Canton-S which does not have a loss of function of these genes. Furthermore, significantly increased expression of Cyp4p1 and Cyp4p2 was induced in 91-R in response to sublethal DDT exposure, however, such induction did not occur in the DDT treated Canton-S. Additionally, fixed nucleotide variation within putative transcription factor binding sites of Cyp4p1 and Cyp4p2 promoters were observed between 91-R and Canton-S, however, their impact on transcription remains to be determined. Two GAL4/UAS transgenic strains with integrated heat shock-inducible Cyp4p1- or Cyp4p2-RNAi constructs within wild-type genetic backgrounds were developed. Following heat shock induction of Cyp4p1 and Cyp4p2 knockdown, these transgenic lines showed increased DDT mortality as compared to their corresponding non-heat shock controls. These results provide a functional link of Cyp4p1 and Cyp4p2 in conferring tolerance to DDT exposure.

CYP4CJ1-mediated gossypol and tannic acid tolerance in Aphis gossypii Glover

Cytochrome P450 monooxygenases play a key role in herbivorous pest adaptation to host plants by the detoxification against plant allelochemicals. A new P450 gene (CYP4CJ1) was identified from Aphis gossypii, which displayed a positive response to plant allelochemicals. The transcript levels of CYP4CJ1 could be significantly induced by both gossypol and tannic acid. Knockdown of CYP4CJ1 increased the sensitivity of A. gossypii to these two plant allelochemicals. These results suggest that CYP4CJ1 could be involved in the tolerance of A. gossypii to some plant allelochemicals. Subsequently, we examined the regulatory mechanism of CYP4CJ1 based on the transcriptional and post-transcriptional level. A promoter region from -1422 to -1166 of CYP4CJ1 was identified, which was an essential plant allelochemical responsive region. In addition, miR-4133-3p was found to participate in the regulation of CYP4CJ1 post-transcriptionally. Our results suggest that the transcript abundance of CYP4CJ1, following the exposure of A. gossypii to gossypol and tannic acid can be attributed to both the transcriptional and post-transcriptional regulation mechanisms. These results are important for understanding the roles of P450s in the plant allelochemical tolerance of A. gossypii.

Molecular basis of permethrin and DDT resistance in an Anopheles funestus population from Benin

BACKGROUND

Insecticide resistance in Anopheles mosquitoes is threatening the success of malaria control programmes. In order to implement suitable insecticide resistance management strategies, it is necessary to understand the underlying mechanisms involved. To achieve this, the molecular basis of permethrin and DDT resistance in the principal malaria vector, Anopheles funestus from inland Benin (Kpome), was investigated.

RESULTS

Here, using a microarray-based genome-wide transcription and qRT-PCR analysis, we showed that metabolic resistance mechanisms through over-expression of cytochrome P450 and glutathione S-transferase genes (GSTs) are a major contributor to DDT and permethrin resistance in Anopheles funestus from Kpome. The GSTe2 gene was the most upregulated detoxification gene in both DDT- [fold-change (FC: 16.0)] and permethrin-resistant (FC: 18.1) mosquitoes suggesting that upregulation of this gene could contribute to DDT resistance and cross-resistance to permethrin. CYP6P9a and CYP6P9b genes that have been previously associated with pyrethroid resistance were also significantly overexpressed with FC 5.4 and 4.8, respectively, in a permethrin resistant population. Noticeably, the GSTs, GSTd1-5 and GSTd3, were more upregulated in DDT-resistant than in permethrin-resistant Anopheles funestus suggesting these genes are more implicated in DDT resistance. The absence of the L1014F or L1014S kdr mutations in the voltage-gated sodium channel gene coupled with the lack of directional selection at the gene further supported that knockdown resistance plays little role in this resistance.

CONCLUSIONS

The major role played by metabolic resistance to pyrethroids in this An. funestus population in Benin suggests that using novel control tools combining the P450 synergist piperonyl butoxide (PBO), such as PBO-based bednets, could help manage the growing pyrethroid resistance in this malaria vector in Benin.

Comparative CYP-omic analysis between the DDT-susceptible and -resistant Drosophila melanogaster strains 91-C and 91-R

BACKGROUND

Cytochrome P450 monooxygenases (P450s) are involved in the biosynthesis of endogenous intracellular compounds and the metabolism of xenobiotics, including chemical insecticides. We investigated the structural and expression level variance across all P450 genes with respect to the evolution of insecticide resistance under multigenerational dichlorodiphenyltrichloroethane (DDT) selection.

RESULTS

RNA-sequencing (RNA-seq) and reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR) indicated that the transcript levels of seven P450 genes were significantly up-regulated and three P450 genes were down-regulated in the DDT-resistant strain 91-R, as compared to the control strain 91-C. The overexpression of Cyp6g1 was associated with the presence of an Accord and an HMS-Beagle element insertion in the 5' upstream region in conjunction with copy number variation in the 91-R strain, but not in the 91-C strain. A total of 122 (50.2%) fixed nonsynonymous (amino acid-changing) mutations were found between 91-C and 91-R, and 20 (8.2%) resulted in amino acid changes within functional domains. Three P450 proteins were truncated as a result of premature stop codons and fixed between strains.

CONCLUSION

Our results demonstrate that a combination of changes in P450 protein-coding regions and transcript levels are possibly associated with DDT resistance, and thereby suggest that selection for variant function may occur within this gene family in response to chronic DDT exposure. © 2018 Society of Chemical Industry.

Variation in Mitochondria-Derived Transcript Levels Associated With DDT Resistance in the 91-R Strain of Drosophila melanogaster (Diptera: Drosophilidae)

The organochloride insecticide dichlorodiphenyltrichloroethane (DDT) and its metabolites can increase cellular levels of reactive oxygen species (ROS), cause mitochondrial dysfunction, and induce apoptosis. The highly DDT-resistant Drosophila melanogaster Meigen 1830 (Drosophila) strain, 91-R, and its susceptible control, 91-C, were used to investigate functional and structural changes among mitochondrial-derived pathways. Resequencing of mitochondrial genomes (mitogenomes) detected no structural differences between 91-R and 91-C, whereas RNA-seq suggested the differential expression of 221 mitochondrial-associated genes. Reverse transcriptase-quantitative PCR validation of 33 candidates confirmed that transcripts for six genes (Cyp12d1-p, Cyp12a4, cyt-c-d, COX5BL, COX7AL, CG17140) were significantly upregulated and two genes (Dif, Rel) were significantly downregulated in 91-R. Among the upregulated genes, four genes are duplicated within the reference genome (cyt-c-d, CG17140, COX5BL, and COX7AL). The predicted functions of the differentially expressed genes, or known functions of closely related genes, suggest that 91-R utilizes existing ROS regulation pathways of the mitochondria to combat increased ROS levels from exposure to DDT. This study represents, to our knowledge, the initial investigation of mitochondrial genome sequence variants and functional adaptations in responses to intense DDT selection and provides insights into potential adaptations of ROS management associated with DDT selection in Drosophila.

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.

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.

Insights into DDT Resistance from the Drosophila melanogaster Genetic Reference Panel

Insecticide resistance is considered a classic model of microevolution, where a strong selective agent is applied to a large natural population, resulting in a change in frequency of alleles that confer resistance. While many insecticide resistance variants have been characterized at the gene level, they are typically single genes of large effect identified in highly resistant pest species. In contrast, multiple variants have been implicated in DDT resistance in Drosophila melanogaster; however, only the Cyp6g1 locus has previously been shown to be relevant to field populations. Here we use genome-wide association studies (GWAS) to identify DDT-associated polygenes and use selective sweep analyses to assess their adaptive significance. We identify and verify two candidate DDT resistance loci. A largely uncharacterized gene, CG10737, has a function in muscles that ameliorates the effects of DDT, while a putative detoxifying P450, Cyp6w1, shows compelling evidence of positive selection.

Identification of a novel cytochrome P450 CYP321B1 gene from tobacco cutworm (Spodoptera litura) and RNA interference to evaluate its role in commonly used insecticides

Insect cytochrome P450 monooxygenases (CYPs or P450s) play an important role in detoxifying insecticides leading to resistance in insect populations. A polyphagous pest, Spodoptera litura, has developed resistance to a wide range of insecticides. In the present study, a novel P450 gene, CYP321B1, was cloned from S. litura. The function of CYP321B1 was assessed using RNA interference (RNAi) and monitoring resistance levels for three commonly used insecticides, including chlorpyrifos, β-cypermethrin and methomyl. The full-length complementary DNA sequence of CYP321B1 is 1814 bp long with an open reading frame of 1 488 bp encoding 495 amino acid residues. Quantitative reverse-transcriptase polymerase chain reaction analyses during larval and pupal development indicated that CYP321B1 expression was highest in the midgut of fifth-instar larvae, followed by fat body and cuticle. The expression of CYP321B1 in the midgut was up-regulated by chlorpyrifos, β-cypermethrin and methomyl with both lethal concentration at 15% (LC₁₅ ) (50, 100 and 150 μg/mL, respectively) and 50%(LC₅₀ ) dosages (100, 200 and 300 μg/mL, respectively). Addition of piperonyl butoxide (PBO) significantly increased the toxicity of chlorpyrifos, β-cypermethrin and methomyl to S. litura, suggesting a marked synergism of the three insecticides with PBO and P450-mediated detoxification. RNAi-mediated silencing of CYP321B1 further increased mortality by 25.6% and 38.9% when the fifth-instar larvae were exposed to chlorpyrifos and β-cypermethrin, respectively, at the LC₅₀ dose levels. The results demonstrate that CYP321B1 might play an important role in chlorpyrifos and β-cypermethrin detoxification in S. litura.

Differential Electrophysiological Responses to Odorant Isotopologues in Drosophilid Antennae

Olfaction presents the ultimate challenge to molecular recognition as thousands of molecules have to be recognized by far fewer olfactory receptors. We have presented evidence that Drosophila readily distinguish odorants based on their molecular vibrations using a battery of behavioral assays suggesting engagement of a molecular vibration-sensing component. Here we interrogate electrophysiologically the antennae of four Drosophilids and demonstrate conserved differential response amplitudes to aldehydes, alcohols, ketones, nitriles, and their deuterated isotopologues. Certain deuterated odorants evoked larger electroantennogram (EAG) amplitudes, while the response to the normal odorant was elevated in others. Significantly, benzonitrile isotopologues were not distinguishable as predicted. This suggests that isotopologue-specific EAG amplitudes result from differential activation of specific olfactory receptors. In support of this, odorants with as few as two deuteria evoke distinct EAG amplitudes from their normal isotopologues, and this is independent of the size of the deuterated molecule. Importantly, we find no evidence that these isotopologue-specific amplitudes depend on perireceptor mechanisms or other pertinent physical property of the deuterated odorants. Rather, our results strongly suggest that Drosophilid olfactory receptors are activated by molecular vibrations differentiating similarly sized and shaped odorants in vivo, yielding sufficient differential information to drive behavioral choices.

Resistance to lambda-cyhalothrin in Spanish field populations of Ceratitis capitata and metabolic resistance mediated by P450 in a resistant strain

BACKGROUND

The withdrawal of malathion in the European Union in 2009 resulted in a large increase in lambda-cyhalothrin applications for the control of the Mediterranean fruit fly, Ceratitis capitata, in Spanish citrus crops.

RESULTS

Spanish field populations of C. capitata have developed resistance to lambda-cyhalothrin (6-14-fold), achieving LC50 values (129-287 ppm) higher than the recommended concentration for field treatments (125 ppm). These results contrast with the high susceptibility to lambda-cyhalothrin found in three Tunisian field populations. We have studied the mechanism of resistance in the laboratory-selected resistant strain W-1Kλ (205-fold resistance). Bioassays with synergists showed that resistance was almost completely suppressed by the P450 inhibitor PBO. The study of the expression of 53 P450 genes belonging to the CYP4, CYP6, CYP9 and CYP12 families in C. capitata revealed that CYP6A51 was overexpressed (13-18-fold) in the resistant strain. The W-1Kλ strain also showed high levels of cross-resistance to etofenprox (240-fold) and deltamethrin (150-fold).

CONCLUSION

Field-evolved resistance to lambda-cyhalothrin has been found in C. capitata. Metabolic resistance mediated by P450 appears to be the main resistance mechanism in the resistant strain W-1Kλ. The levels of cross-resistance found may compromise the effectiveness of other pyrethroids for the control of this species. © 2014 Society of Chemical Industry.

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.

Insecticide-mediated up-regulation of cytochrome P450 genes in the red flour beetle (Tribolium castaneum)

Some cytochrome P450 (CYP) genes are known for their rapid up-regulation in response to insecticide exposures in insects. To date, however, limited information is available with respect to the relationships among the insecticide type, insecticide concentration, exposure duration and the up-regulated CYP genes. In this study, we examined the transcriptional response of eight selected CYP genes, including CYP4G7, CYP4Q4, CYP4BR3, CYP12H1, CYP6BK11, CYP9D4, CYP9Z5 and CYP345A1, to each of four insecticides in the red flour beetle, Tribolium castaneum. Reverse transcription quantitative PCR (RT-qPCR) revealed that CYP4G7 and CYP345A1 can be significantly up-regulated by cypermethrin (1.97- and 2.06-fold, respectively), permethrin (2.00- and 2.03-fold) and lambda-cyhalothrin (1.73- and 1.81-fold), whereas CYP4BR3 and CYP345A1 can be significantly up-regulated by imidacloprid (1.99- and 1.83-fold) when 20-day larvae were exposed to each of these insecticides at the concentration of LC₂₀ for 24 h. Our studies also showed that similar levels of up-regulation can be achieved for CYP4G7, CYP4BR3 and CYP345A1 by cypermethrin, permethrin, lambda-cyhalothrin or imidacloprid with approximately one fourth of LC₂₀ in 6 h. Our study demonstrated that up-regulation of these CYP genes was rapid and only required low concentrations of insecticides, and the up-regulation not only depended on the CYP genes but also the type of insecticides. Our results along with those from previous studies also indicated that there were no specific patterns for predicting the up-regulation of specific CYP gene families based on the insecticide classification.

Identification of a novel cytochrome P450 gene, CYP321E1 from the diamondback moth, Plutella xylostella (L.) and RNA interference to evaluate its role in chlorantraniliprole resistance

Insect cytochrome P450 monooxygenases (P450s) play an important role in catalysis of many reactions leading to insecticides resistance. Our previous studies on transcriptome analysis of chlorantraniliprole-resistant development in the diamondback moth, Plutella xylostella revealed that up-regulation of cytochrome P450s are one of the main factors leading to the development of chlorantraniliprole resistance. Here, we report for the first time a novel cytochrome P450 gene CYP321E1, which belongs to the cytochrome P450 gene family CYP321. Real-time quantitative PCR (RT-qPCR) analyses indicated that CYP321E1 was expressed at all developmental stages of P. xylostella but was highest in the fourth-instar larvae; furthermore, the relatively high expression was observed in the midgut of the fourth-instar larvae, followed by fat bodies and epidermis. The expression of CYP321E1 in P. xylostella was differentially affected by three representative insecticides, including alphamethrin, abamectin and chlorantraniliprole. Among them, the exposure to chlorantraniliprole resulted in the largest transcript level of this cytochrome P450 gene. The findings suggested potential involvement of CYP321E1 in chlorantraniliprole resistance of P. xylostella. To assess the functional link of CYP321E1 to chlorantraniliprole resistance, RNA interference (RNAi)-mediated gene silencing by double stranded RNA (dsRNA) injecting was used. Results revealed that injection delivery of dsRNA can greatly reduce gene expression after 24 h. As a consequence of RNAi, a significant increment in mortality of larvae injected CYP321E1 dsRNA was observed after 24 h of exposure to chlorantraniliprole. These results strongly support our notion that this novel cytochrome P450 gene plays an important role in chlorantraniliprole detoxification in the diamondback moth and is partly responsible for its resistance.

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.

Nrf2/Maf-binding-site-containing functional Cyp6a2 allele is associated with DDT resistance in Drosophila melanogaster

BACKGROUND

Increased insecticide detoxification mediated by cytochrome P450s is a common mechanism of insecticide resistance. Although Cyp6a2 has been observed to be overexpressed in many 4,4'-dichlorodiphenyltrichloroethane (DDT)-resistant strains of Drosophila melanogaster, how Cyp6a2 is regulated and whether its overproduction confers DDT resistance remain elusive.

RESULTS

Molecular analysis identified five Cyp6a2 alleles (Cyp6a2(Canton) (-S-1) , Cyp6a2(Canton) (-S-2) , Cyp6a2(91-C) , Cyp6a2(91-R) and Cyp6a2(Wisconsin) (-) (WD) ) from four D. melanogaster strains, notably differing in the presence or absence of an intact Nrf2/Maf (a transcription factor) binding site in the 5'-promoter core region, a 'G1410' frameshift deletion mutation in the heme-binding region and a long terminal repeat (LTR) of transposable element 17.6 in the 3'-untranslated region (UTR). Linkage analysis confirmed that DDT resistance was genetically linked to a Nrf2/Maf-binding-site-containing, LTR-lacking functional allele of Cyp6a2 (Cyp6a2(91-R) ). The qRT-PCR results showed that overexpression of functional Cyp6a2 was consistently associated with DDT resistance. Luciferase reporter gene assays revealed that an intact Nrf2/Maf binding site in the 5'-promoter core region enhanced the constitutive transcription of Cyp6a2.

CONCLUSION

The results suggest that the Nrf2/Maf binding-site-containing functional Cyp6a2 allele is associated with DDT resistance in the D. melanogaster strains under study.

The mechanisms underlying α-amanitin resistance in Drosophila melanogaster: a microarray analysis

The rapid evolution of toxin resistance in animals has important consequences for the ecology of species and our economy. Pesticide resistance in insects has been a subject of intensive study; however, very little is known about how Drosophila species became resistant to natural toxins with ecological relevance, such as α-amanitin that is produced in deadly poisonous mushrooms. Here we performed a microarray study to elucidate the genes, chromosomal loci, molecular functions, biological processes, and cellular components that contribute to the α-amanitin resistance phenotype in Drosophila melanogaster. We suggest that toxin entry blockage through the cuticle, phase I and II detoxification, sequestration in lipid particles, and proteolytic cleavage of α-amanitin contribute in concert to this quantitative trait. We speculate that the resistance to mushroom toxins in D. melanogaster and perhaps in mycophagous Drosophila species has evolved as cross-resistance to pesticides, other xenobiotic substances, or environmental stress factors.

The molecular genetics of insecticide resistance

The past 60 years have seen a revolution in our understanding of the molecular genetics of insecticide resistance. While at first the field was split by arguments about the relative importance of mono- vs. polygenic resistance and field- vs. laboratory-based selection, the application of molecular cloning to insecticide targets and to the metabolic enzymes that degrade insecticides before they reach those targets has brought out an exponential growth in our understanding of the mutations involved. Molecular analysis has confirmed the relative importance of single major genes in target-site resistance and has also revealed some interesting surprises about the multi-gene families, such as cytochrome P450s, involved in metabolic resistance. Identification of the mutations involved in resistance has also led to parallel advances in our understanding of the enzymes and receptors involved, often with implications for the role of these receptors in humans. This Review seeks to provide an historical perspective on the impact of molecular biology on our understanding of resistance and to begin to look forward to the likely impact of rapid advances in both sequencing and genome-wide association analysis.

Genomic regions harboring insecticide resistance-associated Cyp genes are enriched by transposable element fragments carrying putative transcription factor binding sites in two sibling Drosophila species

In the present study, an in silico analysis was performed to identify transposable element (TE) fragments inserted in Cyps with functions associated with resistance to insecticides and developmental regulation as well as in neighboring genes in two sibling species, Drosophila melanogaster and Drosophila simulans. The Cyps associated with insecticide resistance and their neighboring non-Cyp genes have accumulated a greater number of TE fragments than the other Cyps or a random sample of genes, predominantly in the 5'-flanking regions. Most of the insertions were due to DNA transposons, with DNAREP1 fragments being the most common. These fragments carry putative binding sites for transcription factors, which reinforces the hypothesis that DNAREP1 may influence gene regulation and play a role in the adaptation of the Drosophila species.

Identification of cytochrome P450 monooxygenase genes and their expression profiles in cyhalothrin-treated Colorado potato beetle, Leptinotarsa decemlineata

Based on a Leptinotarsa decemlineata transcriptome dataset and the GenBank sequences, a total of 74 cytochrome P450 monooxygenase genes (Cyps) were identified. These genes fell into CYP2 clan, mitochondrial clan, CYP3 clan and CYP4 clan, and were classified into 19 families and 35 subfamilies according to standard nomenclature. Two new families were discovered in CYP4 clan, and were named CYP412 and CYP413 respectively. Four new families that were recently discovered in Tribolium castaneum, including mitochondrial family CYP353, CYP3 clan families CYP345 and CYP347, and CYP4 clan family CYP350, were also found in L. decemlineata. The phylogenetic trees of CYPs from L. decemlineata and other representative insect species were constructed, and these trees provided evolutionary insight for the genetic distance. Our results facilitate further researches to understand the functions and evolution of L. decemlineata Cyp genes. In order to find cyhalothrin-inducible Cyp genes, the expression levels of Cyps belonging to CYP12, CYP6, CYP9 and CYP4 families were determined by quantitative reverse transcriptase-PCR in cyhalothrin-treated and control fourth-instar larvae. Nine Cyp genes, i.e., Cyp12H2, Cyp6BH2, Cyp6BJ1, Cyp6BQ17, Cyp6EG1, Cyp6EH1, Cyp6EJ1 Cyp4BN13v1 and Cyp4BN15, were highly expressed in cyhalothrin-treated larvae. These CYPs are the candidates that are involved in cyhalothrin detoxification.

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

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

Imidacloprid does not induce Cyp genes involved in insecticide resistance of a mutant Drosophila melanogaster line

Certain xenobiotics have the capacity to induce the expression of genes involved in various biological phenomena, including insecticide resistance. The induction potential of different chemicals, among them different insecticides, has been documented for a number of insect species. In this study, we have analyzed the induction potential of Imidacloprid, a widely used member of the neonicotinoid insecticide family. Genes Cyp6g1 and Cyp6a2, known to be involved in the resistance of mutant Drosophila melanogaster line MiT[W⁻]3R2 to Imidacloprid and DDT were included in the analyzed sample. We find that Imidacloprid does not induce expression of the analyzed genes.

Genome-wide analysis of genes associated with moderate and high DDT resistance in Drosophila melanogaster

BACKGROUND

Moderate to high DDT resistance in generally associated with overexpression of multiple genes and therefore has been considered to be polygenic. However, very little information is available about the molecular mechanisms that insect populations employ when evolving increased levels of resistance. The presence of common regulatory motifs among resistance-associated genes may help to explain how and why certain suites of genes are preferentially represented in genomic-scale analyses.

RESULTS

A set of commonly differentially expressed genes associated with DDT resistance in the fruit fly was identified on the basis of genome-wide microarray analysis followed by qRT-PCR verification. More genes were observed to be overtranscribed in the highly resistant strain (91-R) than in the moderately resistant strain (Wisconsin) and susceptible strain (Canton-S). Furthermore, possible transcription factor binding sites that occurred in coexpressed resistance-associated genes were discovered by computational motif discovery methods.

CONCLUSION

A glucocorticoid receptor (GR)-like putative transcription factor binding motif (TFBM) was observed to be associated with genes commonly differentially transcribed in both the 91-R and Wisconsin lines of DDT-resistant Drosophila.

Soluble and membrane-bound Drosophila melanogaster CYP6G1 expressed in Escherichia coli: purification, activity, and binding properties toward multiple pesticides

Cytochrome P450 CYP6G1 has been implicated in the resistance of Drosophila melanogaster to numerous pesticides. While in vivo and in vitro studies have provided insight to the diverse functions of this enzyme, direct studies on the isolated CYP6G1 enzyme have not been possible due to the need for a source of recombinant enzyme. In the current study, the Cyp6g1 gene was isolated from D. melanogaster and re-engineered for heterologous expression in Escherichia coli. Approximately 460 nmol L⁻¹ of P450 holoenzyme were obtained in 500 mL cultures. The recombinant enzyme was located predominantly within the bacterial cytosol. A two-step purification protocol using Ni-chelate affinity chromatography followed by removal of detergent on a hydroxyapatite column produced essentially homogenous enzyme from both soluble and membrane fractions. Recombinant CYP6G1 exhibited p-nitroanisole O-dealkylation activity but was not active against eleven other typical P450 marker substrates. Substrate-induced binding spectra and IC₅₀ values for inhibition of p-nitroanisole O-dealkylation were obtained for a wide selection of pesticides, namely DDT, imidacloprid, chlorfenvinphos, malathion, endosulfan, dieldrin, dicyclanil, lufenuron and carbaryl, supporting previous in vivo and in vitro studies on Drosophila that have suggested that the enzyme is involved in multi-pesticide resistance in insects.

Integrated analysis of cytochrome P450 gene superfamily in the red flour beetle, Tribolium castaneum

Background

The functional and evolutionary diversification of insect cytochrome P450s (CYPs) shaped the success of insects. CYPs constitute one of the largest and oldest gene superfamilies that are found in virtually all aerobic organisms. Because of the availability of whole genome sequence and well functioning RNA interference (RNAi), the red flour beetle, Tribolium castaneum serves as an ideal insect model for conducting functional genomics studies. Although several T. castaneum CYPs had been functionally investigated in our previous studies, the roles of the majority of CYPs remain largely unknown. Here, we comprehensively analyzed the phylogenetic relationship of all T. castaneum CYPs with genes in other insect species, investigated the CYP6BQ gene cluster organization, function and evolution, as well as examined the mitochondrial CYPs gene expression patterns and intron-exon organization.

Results

A total 143 CYPs were identified and classified into 26 families and 59 subfamilies. The phylogenetic trees of CYPs among insects across taxa provided evolutionary insight for the genetic distance and function. The percentage of singleton (33.3%) in T. castaneum CYPs is much less than those in Drosophila melanogaster (52.5%) and Bombyx mori (51.2%). Most members in the largest CYP6BQ gene cluster may make contribution to deltamethrin resistance in QTC279 strain. T. castaneum genome encodes nine mitochondrial CYPs, among them CYP12H1 is only expressed in the final instar larval stage. The intron-exon organizations of these mitochondrial CYPs are highly diverse.

Conclusion

Our studies provide a platform to understand the evolution and functions of T. castaneum CYP gene superfamily which will help reveal the strategies employed by insects to cope with their environment.

Resistance of Australian Helicoverpa armigera to fenvalerate is due to the chimeric P450 enzyme CYP337B3

Worldwide, increasing numbers of insects have evolved resistance to a wide range of pesticides, which hampers their control in the field and, therefore, threatens agriculture. Members of the carboxylesterase and cytochrome P450 monooxygenase superfamilies are prominent candidates to confer metabolic resistance to pyrethroid insecticides. Both carboxylesterases and P450 enzymes have been shown to be involved in pyrethroid resistance in Australian Helicoverpa armigera, the noctuid species possessing by far the most reported resistance cases worldwide. However, specific enzymes responsible for pyrethroid resistance in field populations of this species have not yet been identified. Here, we show that the resistance toward fenvalerate in an Australian strain of H. armigera is due to a unique P450 enzyme, CYP337B3, which arose from unequal crossing-over between two parental P450 genes, resulting in a chimeric enzyme. CYP337B3 is capable of metabolizing fenvalerate into 4'-hydroxyfenvalerate, which exhibits no toxic effect on susceptible larvae; enzymes from the parental P450 genes showed no detectable fenvalerate metabolism. Furthermore, a polymorphic H. armigera strain could be bred into a susceptible line possessing the parental genes CYP337B1 and CYP337B2 and a resistant line possessing only CYP337B3. The exclusive presence of CYP337B3 in resistant insects of this strain confers a 42-fold resistance to fenvalerate. Thus, in addition to previously documented genetic mechanisms of resistance, recombination can also generate selectively advantageous variants, such as this chimeric P450 enzyme with an altered substrate specificity leading to a potent resistance mechanism.

Evolutionary toxicogenomics: diversification of the Cyp12d1 and Cyp12d3 genes in Drosophila species

Gene duplication and divergence are overwhelmingly considered to be the primary mechanisms by which cytochrome P450 monooxygenases (P450s) have radiated into a large and diverse gene superfamily. To address how environmental stress drives the fixation and diversification of gene duplications, we have analyzed Cyp12d1 and Cyp12d3, a pair of duplicated genes found in the sequenced Drosophila genomes of the melanogaster group. The paralog Cyp12d3, which is not found in Drosophila melanogaster, is basal to the melanogaster group, after it split from the obscura group (ca. 50 mya), and has a significant signature of positive selection in two species (D. sechellia and D. ananassae). Examination of the Cyp12d1 region in D. melanogaster wildtype and isoline populations revealed variation both in copy number and sequence, including splice-site variations, which certainly alter gene function. Further investigations of several strains have identified three cases in which differences in the Cyp12d1 gene region are associated with the differences in transcript abundance and transcriptional responses to the environmental stresses that have not been seen for other detoxificative loci. Together, these data highlight the value of using both macro- and microevolutionary approaches in studying the duplication and divergence events associated with detoxification genes and lay important groundwork for future studies in the field of evolutionary toxicogenomics, which uses the principles of phylogenetic analysis to predict possible enzymatic functions.

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

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

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Characterization of five CYP4 genes from Asian citrus psyllid and their expression levels in Candidatus Liberibacter asiaticus-infected and uninfected psyllids

Previously, we reported that Candidatus Liberibacter asiaticus (Las)-infected Diaphorina citri are characterized by lower levels of cytochrome P450 monooxygenases than uninfected counterparts. In the present study, we investigated expression levels of family 4 cytochrome P450 (CYP4) genes in Las-infected and uninfected D.citri adults. Five novel CYP4 genes (CYP4C67, CYP4DA1, CYP4C68, CYP4DB1 and CYP4G70) were identified. Four of the five CYP4 genes were expressed at significantly higher levels in uninfected than Las-infected males, whereas only one was expressed at significantly higher levels in uninfected than Las-infected females. These results suggest that levels of cytochrome P450 monooxygenases in D.citri may be linked to expression levels of these CYP4 genes. Expression of all five CYP4 genes was induced by exposure of D.citri to imidacloprid, suggesting their possible involvement in metabolism of this toxin. Higher expression of the five CYP4 genes was found in nymphs than adults, which is congruent with previous results indicating higher levels of cytochrome P450 monooxygenases in nymphs than adults. These five CYP4 genes may be promising candidates for RNA-interference to silence overexpression of genes associated with insecticide resistance in D.citri. These newly identified genes may also serve as DNA-based screening markers for cytochrome P450-mediated insecticide resistance in field populations of D.citri.

Brief exposures of human body lice to sublethal amounts of ivermectin over-transcribes detoxification genes involved in tolerance

Transcriptional profiling results, using our non-invasive induction assay {short exposure intervals (2-5 h) to sublethal amounts of insecticides [< lethal concentration 3% (LC(3)) at 24 h] administered by stress-reducing means (contact vs. immersion screen) and with induction assessed in a time frame when tolerance is still present [~lethal concentration 90% (LC(90)) in 2-4 h]}, showed that ivermectin-induced detoxification genes from body lice are identified by quantitative real-time PCR analyses. Of the cytochrome P450 monooxygenase and ATP binding cassette transporter genes induced by ivermectin, CYP6CJ1, CYP9AG1, CYP9AG2 and PhABCC4 were respectively most significantly over-expressed, had high basal expression levels and were most closely related to genes from other organisms that metabolized insecticides, including ivermectin. Injection of double-stranded RNAs (dsRNAs) against either CYP9AG2 or PhABCC4 into non-induced female lice reduced their respective transcript level and resulted in increased sensitivity to ivermectin, indicating that these two genes are involved in the xenobiotic metabolism of ivermectin and in the production of tolerance.

Homology modelling of Drosophila cytochrome P450 enzymes associated with insecticide resistance

BACKGROUND

Overexpression of the cytochrome P450 gene Cyp6g1 confers resistance against DDT and a broad range of other insecticides in Drosophila melanogaster Meig. In the absence of crystal structures of CYP6G1 or complexes with its substrates, structural studies rely on homology modelling and ligand docking to understand P450-substrate interactions.

RESULTS

Homology models are presented for CYP6G1, a P450 associated with resistance to DDT and neonicotinoids, and two other enzymes associated with insecticide resistance in D. melanogaster, CYP12D1 and CYP6A2. The models are based on a template of the X-ray structure of the phylogenetically related human CYP3A4, which is known for its broad substrate specificity. The model of CYP6G1 has a much smaller active site cavity than the template. The cavity is also 'V'-shaped and is lined with hydrophobic residues, showing high shape and chemical complementarity with the molecular characteristics of DDT. Comparison of the DDT-CYP6G1 complex and a non-resistant CYP6A2 homology model implies that tight-fit recognition of this insecticide is important in CYP6G1. The active site can accommodate differently shaped substrates ranging from imidacloprid to malathion but not the pyrethroids permethrin and cyfluthrin.

CONCLUSION

The CYP6G1, CYP12D1 and CYP6A2 homology models can provide a structural insight into insecticide resistance in flies overexpressing P450 enzymes with broad substrate specificities.

Decreased detoxification genes and genome size make the human body louse an efficient model to study xenobiotic metabolism

The human body louse, Pediculus humanus humanus, has one of the smallest insect genomes, containing ∼10 775 annotated genes. Annotation of detoxification [cytochrome P450 monooxygenase (P450), glutathione-S-transferase (GST), esterase (Est) and ATP-binding cassette transporter (ABC transporter)] genes revealed that they are dramatically reduced in P. h. humanus compared to other insects except for Apis mellifera. There are 37 P450, 13 GST and 17 Est genes present in P. h. humanus, approximately half the number found in Drosophila melanogaster and Anopheles gambiae. The number of putatively functional ABC transporter genes in P. h. humanus and Ap. mellifera are the same (36) but both have fewer than An. gambiae (44) or Dr. melanogaster (65). The reduction of detoxification genes in P. h. humanus may be a result of this louse's simple life history, in which it does not encounter a wide variety of xenobiotics. Neuronal component genes are highly conserved across different insect species as expected because of their critical function. Although reduced in number, P. h. humanus still retains at least a minimum repertoire of genes known to confer metabolic or toxicokinetic resistance to xenobiotics (eg Cyp3 clade P450s, Delta GSTs, B clade Ests and B/C subfamily ABC transporters), suggestive of its high potential for resistance development.

Copy number variation and transposable elements feature in recent, ongoing adaptation at the Cyp6g1 locus

The increased transcription of the Cyp6g1 gene of Drosophila melanogaster, and consequent resistance to insecticides such as DDT, is a widely cited example of adaptation mediated by cis-regulatory change. A fragment of an Accord transposable element inserted upstream of the Cyp6g1 gene is causally associated with resistance and has spread to high frequencies in populations around the world since the 1940s. Here we report the existence of a natural allelic series at this locus of D. melanogaster, involving copy number variation of Cyp6g1, and two additional transposable element insertions (a P and an HMS-Beagle). We provide evidence that this genetic variation underpins phenotypic variation, as the more derived the allele, the greater the level of DDT resistance. Tracking the spatial and temporal patterns of allele frequency changes indicates that the multiple steps of the allelic series are adaptive. Further, a DDT association study shows that the most resistant allele, Cyp6g1-[BP], is greatly enriched in the top 5% of the phenotypic distribution and accounts for ∼16% of the underlying phenotypic variation in resistance to DDT. In contrast, copy number variation for another candidate resistance gene, Cyp12d1, is not associated with resistance. Thus the Cyp6g1 locus is a major contributor to DDT resistance in field populations, and evolution at this locus features multiple adaptive steps occurring in rapid succession.

The study of insecticide resistance has greatly enriched our understanding of the genetic basis of adaptation, because it represents some of the most intense selection pressures acting on any natural population of eukaryote. Thus it can inform us about the limits of natural selection, both in terms of the number and type of mutations that can arise and also in terms of the rate at which these spread throughout populations. Fifty years ago, studies in Drosophila melanogaster indicated that many genes contributed to DDT resistance. Subsequent research into the Hikone-R strain indicated much of the resistance in this particular strain could be attributed to a single gene known as Cyp6g1. Here we show that there have been successive DDT resistance mutations occurring at the Cyp6g1 locus. They include an increase in gene copy number and the insertion of transposable elements into the regulatory regions of the Cyp6g1 gene. These mutations have swept to high frequencies in natural populations since World War II, when insecticides were first used. D. melanogaster is not a pest and has not been targeted by insecticides, and yet profound changes are occurring within its genome in response to man-made chemicals in the environment.

Metabolism of methoxychlor by the P450-monooxygenase CYP6G1 involved in insecticide resistance of Drosophila melanogaster after expression in cell cultures of Nicotiana tabacum

Cytochrome P450 monooxygenase CYP6G1 of Drosophila melanogaster was heterologously expressed in a cell suspension culture of Nicotiana tabacum. This in vitro system was used to study the capability of CYP6G1 to metabolize the insecticide methoxychlor (=1,1,1-trichloro-2,2-bis(4-methoxyphenyl)ethane, 1) against the background of endogenous enzymes of the corresponding non-transgenic culture. The Cyp6g1-transgenic cell culture metabolized 96% of applied methoxychlor (45.8 microg per assay) within 24 h by demethylation and hydroxylation mainly to trishydroxy and catechol methoxychlor (16 and 17%, resp.). About 34% of the metabolism and the distinct formation of trishydroxy and catechol methoxychlor were due to foreign enzyme CYP6G1. Furthermore, methoxychlor metabolism was inhibited by 43% after simultaneous addition of piperonyl butoxide (458 microg), whereas inhibition in the non-transgenic culture amounted to 92%. Additionally, the rate of glycosylation was reduced in both cultures. These results were supported by the inhibition of the metabolism of the insecticide imidacloprid (6; 20 microg, 24 h) in the Cyp6g1-transgenic culture by 82% in the presence of piperonyl butoxide (200 microg). Due to CYP6G1 being responsible for imidacloprid resistance of Drosophila or being involved in DDT resistance, it is likely that CYP6G1 conveys resistance to methoxychlor (1). Furthermore, treating Drosophila with piperonyl butoxide could weaken the observed resistance phenomena.

Metabolic enzymes associated with xenobiotic and chemosensory responses in Nasonia vitripennis

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

Functional analysis of the cis-acting elements responsible for the induction of the Cyp6a8 and Cyp6g1 genes of Drosophila melanogaster by DDT, phenobarbital and caffeine

Many Drosophila cytochrome P450 or Cyp genes are induced by caffeine and phenobarbital (PB). To understand the induction mechanism, we created Drosophila S2 cell lines stably transformed with different luciferase reporter plasmids carrying upstream DNAs of Cyp6a8 allele of the resistant 91-R strain, and the 1.1-kb upstream DNAs of Cyp6g1 of the 91-R and the susceptible 91-C strains. Following 24 h treatment with dichlorodiphenyltrichloroethane (DDT), caffeine or PB, luciferase activity of all cell lines was determined. Results showed that the 0.1-kb DNA of Cyp6a8 and the upstream DNAs of Cyp6g1 from both strains are not induced by these chemicals in S2 cells. However, the 0.2-, 0.5- and 0.8-kb DNAs of Cyp6a8 showed 13-24-, 4-5- and 2.2-2.7-fold induction with caffeine, PB and DDT, respectively. These DNAs also showed a 2-3-fold synergistic effect of caffeine and PB but not of caffeine and DDT. The results suggest that the cis-regulatory elements for all three chemicals are located within the -11/-199 DNA of Cyp6a8. Furthermore, caffeine and PB inductions appear to be mediated via different cis-elements, whereas caffeine and DDT induction may involve common regulatory elements. These stably transformed cell lines should help understand the mechanism of resistance-associated Cyp gene overexpression in Drosophila.

Two duplicated P450 genes are associated with pyrethroid resistance in Anopheles funestus, a major malaria vector

Pyrethroid resistance in Anopheles funestus is a potential obstacle to malaria control in Africa. Tools are needed to detect resistance in field populations. We have been using a positional cloning approach to identify the major genes conferring pyrethroid resistance in this vector. A quantitative trait locus (QTL) named rp1 explains 87% of the genetic variance in pyrethroid susceptibility in two families from reciprocal crosses between susceptible and resistant strains. Two additional QTLs of minor effect, rp2 and rp3, were also detected. We sequenced a 120-kb BAC clone spanning the rp1 QTL and identified 14 protein-coding genes and one putative pseudogene. Ten of the 14 genes encoded cytochrome P450s, and expression analysis indicated that four of these P450s were differentially expressed between susceptible and resistant strains. Furthermore, two of these genes, CYP6P9 and CYP6P4, which are 25 and 51 times overexpressed in resistant females, are tandemly duplicated in the BAC clone as well as in laboratory and field samples, suggesting that P450 gene duplication could contribute to pyrethroid resistance in An. funestus. Single nucleotide polymorphisms (SNPs) were identified within CYP6P9 and CYP6P4, and genotyping of the progeny of the genetic crosses revealed a maximum penetrance value f(2) = 1, confirming that these SNPs are valid resistance markers in the laboratory strains. This serves as proof of principle that a DNA-based diagnostic test could be designed to trace metabolic resistance in field populations. This will be a major advance for insecticide resistance management in malaria vectors, which requires the early detection of resistance alleles.