Three second chromosome-linked clustered Cyp6 genes show differential constitutive and barbital-induced expression in DDT-resistant and susceptible strains of Drosophila melanogaster

Multiple-P450 Gene Co-Up-Regulation in the Development of Permethrin Resistance in the House Fly, Musca domestica

This paper reports a study conducted at the whole transcriptome level to characterize the P450 genes involved in the development of pyrethroid resistance, utilizing expression profile analyses of 86 cytochrome P450 genes in house fly strains with different levels of resistance to pyrethroids/permethrin. Interactions among the up-regulated P450 genes and possible regulatory factors in different autosomes were examined in house fly lines with different combinations of autosomes from a resistant house fly strain, ALHF. Eleven P450 genes that were significantly up-regulated, with levels > 2-fold those in the resistant ALHF house flies, were in CYP families 4 and 6 and located on autosomes 1, 3 and 5. The expression of these P450 genes was regulated by trans- and/or cis-acting factors, especially on autosomes 1 and 2. An in vivo functional study indicated that the up-regulated P450 genes also conferred permethrin resistance in Drosophila melanogaster transgenic lines. An in vitro functional study confirmed that the up-regulated P450 genes are able to metabolize not only cis- and trans-permethrin, but also two metabolites of permethrin, PBalc and PBald. In silico homology modeling and the molecular docking methodology further support the metabolic capacity of these P450s for permethrin and substrates. Taken together, the findings of this study highlight the important function of multi-up-regulated P450 genes in the development of insecticide resistance in house flies.

Functional Characterization of Drosophila melanogaster CYP6A8 Fatty Acid Hydroxylase

Genomic analysis indicated that the genome of Drosophila melanogaster contains more than 80 cytochrome P450 genes. To date, the enzymatic activity of these P450s has not been extensively studied. Here, the biochemical properties of CYP6A8 were characterized. CYP6A8 was cloned into the pCW vector, and its recombinant enzyme was expressed in Escherichia coli and purified using Ni²⁺-nitrilotriacetate affinity chromatography. Its expression level was approximately 130 nmol per liter of culture. Purified CYP6A8 exhibited a low-spin state in the absolute spectra of the ferric forms. Binding titration analysis indicated that lauric acid and capric acid produced type І spectral changes, with K_d values 28 ± 4 and 144 ± 20 μM, respectively. Ultra-performance liquid chromatography-mass spectrometry analysis showed that the oxidation reaction of lauric acid produced (ω-1)-hydroxylated lauric acid as a major product and ω-hydroxy-lauric acid as a minor product. Steady-state kinetic analysis of lauric acid hydroxylation yielded a k_cat value of 0.038 ± 0.002 min^-1 and a K_m value of 10 ± 2 μM. In addition, capric acid hydroxylation of CYP6A8 yielded kinetic parameters with a k_cat value of 0.135 ± 0.007 min^-1 and a K_m value of 21 ± 4 μM. Because of the importance of various lipids as carbon sources, the metabolic analysis of fatty acids using CYP6A8 in this study can provide an understanding of the biochemical roles of P450 enzymes in many insects, including Drosophila melanogaster.

Transcriptome and population structure of glassy-winged sharpshooters (Homalodisca vitripennis) with varying insecticide resistance in southern California

Background

Homalodisca vitripennis Germar, the glassy-winged sharpshooter, is an invasive insect in California and a critical threat to agriculture through its transmission of the plant pathogen, Xylella fastidiosa. Quarantine, broad-spectrum insecticides, and biological control have been used for population management of H. vitripennis since its invasion and subsequent proliferation throughout California. Recently wide-spread neonicotinoid resistance has been detected in populations of H. vitripennis in the southern portions of California’s Central Valley. In order to better understand potential mechanisms of H. vitripennis neonicotinoid resistance, we performed RNA sequencing on wild-caught insecticide-resistant and relatively susceptible sharpshooters to profile their transcriptome and population structure.

Results

We identified 81 differentially expressed genes with higher expression in resistant individuals. The significant largest differentially expressed candidate gene linked to resistance status was a cytochrome P450 gene with similarity to CYP6A9. Furthermore, we observed an over-enrichment of GO terms representing functions supportive of roles in resistance mechanisms (cytochrome P450s, M13 peptidases, and cuticle structural proteins). Finally, we saw no evidence of broad-scale population structure, perhaps due to H. vitripennis' relatively recent introduction to California or due to the relatively small geographic scale investigated here.

Conclusions

In this work, we characterized the transcriptome of insecticide-resistant and susceptible H. vitripennis and identified candidate genes that may be involved in resistance mechanisms for this species. Future work should seek to build on the transcriptome profiling performed here to confirm the role of the identified genes, particularly the cytochrome P450, in resistance in H. vitripennis. We hope this work helps aid future population management strategies for this and other species with growing insecticide resistance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08939-1.

Poly(ADP-ribose) polymerase 1 in genome-wide expression control in Drosophila

Poly(ADP-ribose) polymerase 1 (PARP-1) is a nuclear enzyme involved in DNA repair and transcription regulation, among other processes. Malignant transformations, tumor progression, the onset of some neuropathies and other disorders have been linked to misregulation of PARP-1 activity. Despite intensive studies during the last few decades, the role of PARP-1 in transcription regulation is still not well understood. In this study, a transcriptomic analysis in Drosophila melanogaster third instar larvae was carried out. A total of 602 genes were identified, showing large-scale changes in their expression levels in the absence of PARP-1 in vivo. Among these genes, several functional gene groups were present, including transcription factors and cytochrome family members. The transcription levels of genes from the same functional group were affected by the absence of PARP-1 in a similar manner. In the absence of PARP-1, all misregulated genes coding for transcription factors were downregulated, whereas all genes coding for members of the cytochrome P450 family were upregulated. The cytochrome P450 proteins contain heme as a cofactor and are involved in oxidoreduction. Significant changes were also observed in the expression of several mobile elements in the absence of PARP-1, suggesting that PARP-1 may be involved in regulating the expression of mobile elements.

A First Glimpse of the Mexican Fruit Fly Anastrepha ludens (Diptera: Tephritidae) Antenna Morphology and Proteome in Response to a Proteinaceous Attractant

Anastrepha ludens is a key pest of mangoes and citrus from Texas to Costa Rica but the mechanisms of odorant perception in this species are poorly understood. Detection of volatiles in insects occurs mainly in the antenna, where molecules penetrate sensillum pores and link to soluble proteins in the hemolymph until reaching specific odor receptors that trigger signal transduction and lead to behavioral responses. Scrutinizing the molecular foundation of odorant perception in A. ludens is necessary to improve biorational management strategies against this pest. After exposing adults of three maturity stages to a proteinaceous attractant, we studied antennal morphology and comparative proteomic profiles using nano-LC-MS/MS with tandem mass tags combined with synchronous precursor selection (SPS)-MS3. Antennas from newly emerged flies exhibited dense agglomerations of olfactory sensory neurons. We discovered 4618 unique proteins in the antennas of A. ludens and identified some associated with odor signaling, including odorant-binding and calcium signaling related proteins, the odorant receptor co-receptor (Orco), and putative odorant-degrading enzymes. Antennas of sexually immature flies exhibited the most upregulation of odor perception proteins compared to mature flies exposed to the attractant. This is the first report where critical molecular players are linked to the odor perception mechanism of A. ludens.

Detoxification Genes Differ Between Cactus-, Fruit-, and Flower-Feeding Drosophila

We use annotated genomes of 14 Drosophila species covering diverse host use phenotypes to test whether 4 gene families that often have detoxification functions are associated with host shifts among species. Bark, slime flux, flower, and generalist necrotic fruit-feeding species all have similar numbers of carboxyl/cholinesterase, glutathione S-transferase, cytochrome P450, and UDP-glucuronosyltransferase genes. However, species feeding on toxic Morinda citrifolia fruit and the fresh fruit-feeding Drosophila suzukii have about 30 and 60 more, respectively. ABC transporters show a different pattern, with the flower-feeding D. elegans and the generalist necrotic fruit and cactus feeder D. hydei having about 20 and >100 more than the other species, respectively. Surprisingly, despite the complex secondary chemistry we find that 3 cactophilic specialists in the mojavensis species cluster have variably fewer genes than any of the other species across all 4 families. We also find 82 positive selection events across the 4 families, with the terminal D. suzukii and M. citrifolia-feeding D. sechellia branches again having the highest number of such events in proportion to their respective branch lengths. Many of the genes involved in these host-use-specific gene number differences or positive selection events lie in specific clades of the gene families that have been recurrently associated with detoxification. Several genes are also found to be involved in multiple duplication and/or positive selection events across the species studied regardless of their host use phenotypes; the most frequently involved are the ABC transporter CG1718, which is not in a specific clade associated with detoxification, and the α-esterase gene cluster, which is.

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.

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.

A novel role of Drosophila cytochrome P450-4e3 in permethrin insecticide tolerance

The exposure of insects to xenobiotics, such as insecticides, triggers a complex defence response necessary for survival. This response includes the induction of genes that encode key Cytochrome P450 monooxygenase detoxification enzymes. Drosophila melanogaster Malpighian (renal) tubules are critical organs in the detoxification and elimination of these foreign compounds, so the tubule response induced by dietary exposure to the insecticide permethrin was examined. We found that expression of the gene encoding Cytochrome P450-4e3 (Cyp4e3) is significantly up-regulated by Drosophila fed on permethrin and that manipulation of Cyp4e3 levels, specifically in the principal cells of the Malpighian tubules, impacts significantly on the survival of permethrin-fed flies. Both dietary exposure to permethrin and Cyp4e3 knockdown cause a significant elevation of oxidative stress-associated markers in the tubules, including H2O2 and lipid peroxidation byproduct, HNE (4-hydroxynonenal). Thus, Cyp4e3 may play an important role in regulating H2O2 levels in the endoplasmic reticulum (ER) where it resides, and its absence triggers a JAK/STAT and NF-κB-mediated stress response, similar to that observed in cells under ER stress. This work increases our understanding of the molecular mechanisms of insecticide detoxification and provides further evidence of the oxidative stress responses induced by permethrin metabolism.

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.

The cytochrome P450 family in the parasitic nematode Haemonchus contortus

•

The Haemonchus contortus genome encodes a large family of cytochrome P450 (CYP) genes.

•

Haemonchus contortus lacks the dramatic CYP family expansions seen in Caenorhabditis elegans.

•

Haemonchus contortus orthologues of C. elegans CYPs share similar expression profiles.

•

The majority of H. contortus CYPs are most highly expressed in larval stages.

•

The parasite intestine is a major site of CYP expression.

Haemonchus contortus, a highly pathogenic and economically important parasitic nematode of sheep, is particularly adept at developing resistance to the anthelmintic drugs used in its treatment and control. The basis of anthelmintic resistance is poorly understood for many commonly used drugs with most research being focused on mechanisms involving drug targets or drug efflux. Altered or increased drug metabolism is a possible mechanism that has yet to receive much attention despite the clear role of xenobiotic metabolism in pesticide resistance in insects. The cytochrome P450s (CYPs) are a large family of drug-metabolising enzymes present in almost all living organisms, but for many years thought to be absent from parasitic nematodes. In this paper, we describe the CYP sequences encoded in the H. contortus genome and compare their expression in different parasite life-stages, sexes and tissues. We developed a novel real-time PCR approach based on partially assembled CYP sequences “tags” and confirmed findings in the subsequent draft genome with RNA-seq. Constitutive expression was highest in larval stages for the majority of CYPs, although higher expression was detected in the adult male or female for a small subset of genes. Many CYPs were expressed in the worm intestine. A number of H. contortus genes share high identity with Caenorhabditis elegans CYPs and the similarity in their expression profiles supports their classification as putative orthologues. Notably, H. contortus appears to lack the dramatic CYP subfamily expansions seen in C. elegans and other species, which are typical of CYPs with exogenous roles. However, a small group of H. contortus genes cluster with the C. elegans CYP34 and CYP35 subfamilies and may represent candidate xenobiotic metabolising genes in the parasite.

Marine invertebrate xenobiotic-activated nuclear receptors: their application as sensor elements in high-throughput bioassays for marine bioactive compounds

Developing high-throughput assays to screen marine extracts for bioactive compounds presents both conceptual and technical challenges. One major challenge is to develop assays that have well-grounded ecological and evolutionary rationales. In this review we propose that a specific group of ligand-activated transcription factors are particularly well-suited to act as sensors in such bioassays. More specifically, xenobiotic-activated nuclear receptors (XANRs) regulate transcription of genes involved in xenobiotic detoxification. XANR ligand-binding domains (LBDs) may adaptively evolve to bind those bioactive, and potentially toxic, compounds to which organisms are normally exposed to through their specific diets. A brief overview of the function and taxonomic distribution of both vertebrate and invertebrate XANRs is first provided. Proof-of-concept experiments are then described which confirm that a filter-feeding marine invertebrate XANR LBD is activated by marine bioactive compounds. We speculate that increasing access to marine invertebrate genome sequence data, in combination with the expression of functional recombinant marine invertebrate XANR LBDs, will facilitate the generation of high-throughput bioassays/biosensors of widely differing specificities, but all based on activation of XANR LBDs. Such assays may find application in screening marine extracts for bioactive compounds that could act as drug lead compounds.

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.

Molecular population genetics of the NADPH cytochrome P450 reductase (CPR) gene in Anopheles minimus

Development of insecticide resistance (IR) in mosquito vectors is a primary huddle to malaria control program. Since IR has genetic basis, and genes constantly evolve with response to environment for adaptation to organisms, it is important to know evolutionary pattern of genes conferring IR in malaria vectors. The mosquito Anopheles minimus is a major malaria vector of the Southeast (SE) Asia and India and is susceptible to all insecticides, and thus of interest to know if natural selection has shaped variations in the gene conferring IR. If not, the DNA fragment of such a gene could be used to infer population structure and demography of this species of malaria vector. We have therefore sequenced a ~569 bp DNA segment of the NADPH cytochrome P450 reductase (CPR) gene (widely known to confer IR) in 123 individuals of An. minimus collected in 10 different locations (eight Indian, one Thai and one Vietnamese). Two Indian population samples were completely mono-morphic in the CPR gene. In general, low genetic diversity was found with no evidence of natural selection in this gene. The data were therefore analyzed to infer population structure and demography of this species. The 10 populations could be genetically differentiated into four different groups; the samples from Thailand and Vietnam contained high nucleotide diversity. All the 10 populations conform to demographic equilibrium model with signature of past population expansion in four populations. The results in general indicate that the An. minimus mosquitoes sampled in the two SE Asian localities contain several genetic characteristics of being parts of the ancestral population.

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 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.

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

Background

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

Results

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

Conclusion

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

Use of mutagenesis, genetic mapping and next generation transcriptomics to investigate insecticide resistance mechanisms

Insecticide resistance is a worldwide problem with major impact on agriculture and human health. Understanding the underlying molecular mechanisms is crucial for the management of the phenomenon; however, this information often comes late with respect to the implementation of efficient counter-measures, particularly in the case of metabolism-based resistance mechanisms. We employed a genome-wide insertional mutagenesis screen to Drosophila melanogaster, using a Minos-based construct, and retrieved a line (MiT[w⁻]3R2) resistant to the neonicotinoid insecticide Imidacloprid. Biochemical and bioassay data indicated that resistance was due to increased P450 detoxification. Deep sequencing transcriptomic analysis revealed substantial over- and under-representation of 357 transcripts in the resistant line, including statistically significant changes in mixed function oxidases, peptidases and cuticular proteins. Three P450 genes (Cyp4p2, Cyp6a2 and Cyp6g1) located on the 2R chromosome, are highly up-regulated in mutant flies compared to susceptible Drosophila. One of them (Cyp6g1) has been already described as a major factor for Imidacloprid resistance, which validated the approach. Elevated expression of the Cyp4p2 was not previously documented in Drosophila lines resistant to neonicotinoids. In silico analysis using the Drosophila reference genome failed to detect transcription binding factors or microRNAs associated with the over-expressed Cyp genes. The resistant line did not contain a Minos insertion in its chromosomes, suggesting a hit-and-run event, i.e. an insertion of the transposable element, followed by an excision which caused the mutation. Genetic mapping placed the resistance locus to the right arm of the second chromosome, within a ∼1 Mb region, where the highly up-regulated Cyp6g1 gene is located. The nature of the unknown mutation that causes resistance is discussed on the basis of these results.

Transcriptional regulation of xenobiotic detoxification in Drosophila

Living organisms, from bacteria to humans, display a coordinated transcriptional response to xenobiotic exposure, inducing enzymes and transporters that facilitate detoxification. Several transcription factors have been identified in vertebrates that contribute to this regulatory response. In contrast, little is known about this pathway in insects. Here we show that the Drosophila Nrf2 (NF-E2-related factor 2) ortholog CncC (cap 'n' collar isoform-C) is a central regulator of xenobiotic detoxification responses. A binding site for CncC and its heterodimer partner Maf (muscle aponeurosis fibromatosis) is sufficient and necessary for robust transcriptional responses to three xenobiotic compounds: phenobarbital (PB), chlorpromazine, and caffeine. Genetic manipulations that alter the levels of CncC or its negative regulator, Keap1 (Kelch-like ECH-associated protein 1), lead to predictable changes in xenobiotic-inducible gene expression. Transcriptional profiling studies reveal that more than half of the genes regulated by PB are also controlled by CncC. Consistent with these effects on detoxification gene expression, activation of the CncC/Keap1 pathway in Drosophila is sufficient to confer resistance to the lethal effects of the pesticide malathion. These studies establish a molecular mechanism for the regulation of xenobiotic detoxification in Drosophila and have implications for controlling insect populations and the spread of insect-borne human diseases.

Regulation of cytochrome P450 expression in Drosophila: Genomic insights

Genomic tools such as the availability of the Drosophila genome sequence, the relative ease of stable transformation, and DNA microarrays have made the fruit fly a powerful model in insecticide toxicology research. We have used transgenic promoter-GFP constructs to document the detailed pattern of induced Cyp6a2 gene expression in larval and adult Drosophila tissues. We also compared various insecticides and xenobiotics for their ability to induce this cytochrome P450 gene, and show that the pattern of Cyp6a2 inducibility is comparable to that of vertebrate CYP2B genes, and different from that of vertebrate CYP1A genes, suggesting a degree of evolutionary conservation for the "phenobarbital-type" induction mechanism. Our results are compared to the increasingly diverse reports on P450 induction that can be gleaned from whole genome or from "detox" microarray experiments in Drosophila. These suggest that only a third of the genomic repertoire of CYP genes is inducible by xenobiotics, and that there are distinct subsets of inducers / induced genes, suggesting multiple xenobiotic transduction mechanisms. A relationship between induction and resistance is not supported by expression data from the literature. The relative abundance of expression data now available is in contrast to the paucity of studies on functional expression of P450 enzymes, and this remains a challenge for our understanding of the toxicokinetic aspects of insecticide action.

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.

Characterization of two novel cytochrome P450 genes in insecticide-resistant house-flies

Two novel P450 cDNAs, CYP6A36 and CYP6A37, were isolated from house-flies. Putative protein sequences deduced from the cDNA sequences shared 58% identity. Predicted protein sequences of CYP6A36 and CYP6A37 from pyrethroid-resistant ALHF house-flies were identical to their corresponding orthologues in susceptible aabys flies. Expression of CYP6A36 was developmentally regulated with significant overexpression in ALHF compared with susceptible CS flies. Overexpression of CYP6A36 was detected in ALHF abdomen, where the primary detoxification organs of the midgut and fat body are located. CYP6A37, however, was expressed in all tested life stages, with no significant difference in expression between ALHF and CS. Genetic linkage analysis located CYP6A36 on autosome 5; overexpression of CYP6A36 was linked to the factors on autosomes 1 and 2, corresponding to the linkage of P450-mediated resistance in ALHF. This evidence suggests the importance of CYP6A36 in detoxification of insecticides and evolution of insecticide resistance in ALHF.

Genomic analysis of detoxification genes in the mosquito Aedes aegypti

Annotation of the recently determined genome sequence of the major dengue vector, Aedes aegypti, reveals an abundance of detoxification genes. Here, we report the presence of 235 members of the cytochrome P450, glutathione transferase and carboxy/cholinesterase families in Ae. aegypti. This gene count represents an increase of 58% and 36% compared with the fruitfly, Drosophila melanogaster, and the malaria mosquito, Anopheles gambiae, respectively. The expansion is not uniform within the gene families. Secure orthologs can be found across the insect species for enzymes that have presumed or proven biosynthetic or housekeeping roles. In contrast, subsets of these gene families that are associated with general xenobiotic detoxification, in particular the CYP6, CYP9 and alpha esterase families, have expanded in Ae. aegypti. In order to identify detoxification genes associated with resistance to insecticides we constructed an array containing unique oligonucleotide probes for these genes and compared their expression level in insecticide resistant and susceptible strains. Several candidate genes were identified with the majority belonging to two gene families, the CYP9 P450s and the Epsilon GSTs. This 'Ae. aegypti Detox Chip' will facilitate the implementation of insecticide resistance management strategies for arboviral control programmes.

Transposable elements are enriched within or in close proximity to xenobiotic-metabolizing cytochrome P450 genes

Background

Transposons, i.e. transposable elements (TEs), are the major internal spontaneous mutation agents for the variability of eukaryotic genomes. To address the general issue of whether transposons mediate genomic changes in environment-adaptation genes, we scanned two alleles per each of the six xenobiotic-metabolizing Helicoverpa zea cytochrome P450 loci, including CYP6B8, CYP6B27, CYP321A1, CYP321A2, CYP9A12v3 and CYP9A14, for the presence of transposon insertions by genome walking and sequence analysis. We also scanned thirteen Drosophila melanogaster P450s genes for TE insertions by in silico mapping and literature search.

Results

Twelve novel transposons, including LINEs (long interspersed nuclear elements), SINEs (short interspersed nuclear elements), MITEs (miniature inverted-repeat transposable elements), one full-length transib-like transposon, and one full-length Tcl-like DNA transpson, are identified from the alleles of the six H. zea P450 genes. The twelve transposons are inserted into the 5'flanking region, 3'flanking region, exon, or intron of the six environment-adaptation P450 genes. In D. melanogaster, seven out of the eight Drosophila P450s (CYP4E2, CYP6A2, CYP6A8, CYP6A9, CYP6G1, CYP6W1, CYP12A4, CYP12D1) implicated in insecticide resistance are associated with a variety of transposons. By contrast, all the five Drosophila P450s (CYP302A1, CYP306A1, CYP307A1, CYP314A1 and CYP315A1) involved in ecdysone biosynthesis and developmental regulation are free of TE insertions.

Conclusion

These results indicate that TEs are selectively retained within or in close proximity to xenobiotic-metabolizing P450 genes.

Mapping a quantitative trait locus (QTL) conferring pyrethroid resistance in the African malaria vector Anopheles funestus

BACKGROUND

Pyrethroid resistance in Anopheles funestus populations has led to an increase in malaria transmission in southern Africa. Resistance has been attributed to elevated activities of cytochrome P450s but the molecular basis underlying this metabolic resistance is unknown. Microsatellite and SNP markers were used to construct a linkage map and to detect a quantitative trait locus (QTL) associated with pyrethroid resistance in the FUMOZ-R strain of An. funestus from Mozambique.

RESULTS

By genotyping 349 F2 individuals from 11 independent families, a single major QTL, rp1, at the telomeric end of chromosome 2R was identified. The rp1 QTL appears to present a major effect since it accounts for more than 60% of the variance in susceptibility to permethrin. This QTL has a strong additive genetic effect with respect to susceptibility. Candidate genes associated with pyrethroid resistance in other species were physically mapped to An. funestus polytene chromosomes. This showed that rp1 is genetically linked to a cluster of CYP6 cytochrome P450 genes located on division 9 of chromosome 2R and confirmed earlier reports that pyrethroid resistance in this strain is not associated with target site mutations (knockdown resistance).

CONCLUSION

We hypothesize that one or more of these CYP6 P450s clustered on chromosome 2R confers pyrethroid resistance in the FUMOZ-R strain of An. funestus.

A comparison of Drosophila melanogaster detoxification gene induction responses for six insecticides, caffeine and phenobarbital

Modifications of metabolic pathways are important in insecticide resistance evolution. Mutations leading to changes in expression levels or substrate specificities of cytochrome P450 (P450), glutathione-S-transferase (GST) and esterase genes have been linked to many cases of resistance with the responsible enzyme shown to utilize the insecticide as a substrate. Many studies show that the substrates of enzymes are capable of inducing the expression of those enzymes. We investigated if this was the case for insecticides and the enzymes responsible for their metabolism. The induction responses for P450s, GSTs and esterases to six different insecticides were investigated using a custom designed microarray in Drosophila melanogaster. Even though these gene families can all contribute to insecticide resistance, their induction responses when exposed to insecticides are minimal. The insecticides spinosad, diazinon, nitenpyram, lufenuron and dicyclanil did not induce any P450, GST or esterase gene expression after a short exposure to high lethal concentrations of insecticide. DDT elicited the low-level induction of one GST and one P450. These results are in contrast to induction responses we observed for the natural plant compound caffeine and the barbituate drug phenobarbital, both of which highly induced a number of P450 and GST genes under the same short exposure regime. Our results indicate that, under the insecticide exposure conditions we used, constitutive over-expression of metabolic genes play more of a role in insect survival than induction of members of these gene families.

High expression of Cyp6g1, a cytochrome P450 gene, does not necessarily confer DDT resistance in Drosophila melanogaster

Cytochrome P450 monooxygenases, a family of detoxifying enzymes, are thought to confer resistance to various insecticides including DDT. Daborn et al. [Daborn, P., Yen, J.L., Bogwitz, M., Le Goff, G., Feil, et al. 2002. A single p450 allele associated with insecticide resistance in Drosophila. Science 297, 2253-2256.] suggested that the Accord transposable element causes overexpression of a Cyp6g1 allele, which has spread globally and is the basis of DDT resistance in Drosophila melanogaster populations. To determine whether the same phenomenon also operates in other Drosophila strains, we investigated 91-R, 91-C, ry(506), Wisconsin, Canton-SH and Hikone-RH strains. While the LC(50) values for the 91-R and Wisconsin strains are 8348 microg and 447 microg of DDT, respectively, values for the other four strains range between 0.74 to 20.9 microg. As expected, the susceptible ry(506) and 91-C strains have about 16-33-fold lower levels of CYP6G1 mRNA than the resistant 91-R and Wisconsin strains. Surprisingly, CYP6G1 mRNA and protein levels in the Canton-SH and Hikone-RH strains are as high as in the two resistant strains, yet they are as susceptible as the 91-C strain. The susceptible phenotype of the Canton-SH and Hikone-RH strains is not due to mutation in the Cyp6g1 gene; sequence analysis showed that Cyp6g1 alleles of resistant and susceptible strains are very similar and cannot be classified into resistant and susceptible alleles. As observed by others, we also found that only the 5'-upstream DNA of overexpressing alleles of Cyp6g1 has an insertional DNA, which is similar to Accord and Ninja elements. To examine the role of Cyp6g1 in DDT resistance, we substituted the Cyp6g1 allele of the 91-R strain with the allele from the susceptible 91-C strain via recombination and synthesized three recombinant lines. All three lines lacked Accord insertion and showed low expression of Cyp6g1 like the 91-C strain, yet they were as highly resistant as the 91-R strain. We conclude a strain may not have to have Accord insertion in the Cyp6g1 gene and the Cyp6g1 itself may not have to be overexpressed for DDT resistance to occur.

Genome-wide analysis of phenobarbital-inducible genes in Drosophila melanogaster

An oligoarray analysis was conducted to determine the differential expression of genes due to phenobarbital exposure in Drosophila melanogaster (w(1118) strain) third instar larvae. Seventeen genes were observed to be induced with increased expression by a statistical analysis of microarrays approach with a q < or = 0.05. At q < or = 0.12, four more genes (Cyp12d1, DmGstd4, and two genes with unknown function) were found to be up-regulated, and 11 genes with unknown function were found to be down-regulated. Fifteen of these genes, Cyp4d14, Cyp6a2, Cyp6a8, Cyp12d1, Cyp6d5, Cyp6w1, CG2065, DmGstd6, DmGstd7, Amy-p/Amy-d, Ugt86Dd, GC5724, Jheh1, Jheh2 and CG11893, were verified using quantitative real time polymerase chain reaction. Some of these genes have been shown to be over-transcribed in metabolically DDT-resistant Drosophila strains.

The DHR96 nuclear receptor regulates xenobiotic responses in Drosophila

Exposure to xenobiotics such as plant toxins, pollutants, or prescription drugs triggers a defense response, inducing genes that encode key detoxification enzymes. Although xenobiotic responses have been studied in vertebrates, little effort has been made to exploit a simple genetic system for characterizing the molecular basis of this coordinated transcriptional response. We show here that approximately 1000 transcripts are significantly affected by phenobarbital treatment in Drosophila. We also demonstrate that the Drosophila ortholog of the human SXR and CAR xenobiotic receptors, DHR96, plays a role in this response. A DHR96 null mutant displays increased sensitivity to the sedative effects of phenobarbital and the pesticide DDT as well as defects in the expression of many phenobarbital-regulated genes. Metabolic and stress-response genes are also controlled by DHR96, implicating its role in coordinating multiple response pathways. This work establishes a new model system for defining the genetic control of xenobiotic stress responses.

Toxicity of aflatoxin B1 to Helicoverpa zea and bioactivation by cytochrome P450 monooxygenases

Infestation of corn (Zea mays) by corn earworm (Helicoverpa zea) predisposes the plant to infection by Aspergillus fungi and concomitant contamination with the carcinogenic mycotoxin aflatoxin B1 (AFB1). Although effects of ingesting AFB1 are well documented in livestock and humans, the effects on insects that naturally encounter this mycotoxin are not as well defined. Toxicity of AFB1 to different stages of H. zea (first, third, and fifth instars) was evaluated with artificial diets containing varying concentrations. Although not acutely toxic at low concentrations (1-20 ng/g), AFB1 had significant chronic effects, including protracted development, increased mortality, decreased pupation rate, and reduced pupal weight. Sensitivity varied with developmental stage; whereas intermediate concentrations (200 ng/g) caused complete mortality in first instars, this same concentration had no detectable adverse effects on larvae encountering AFB1 in fifth instar. Fifth instars consuming AFB1 at higher concentrations (1 microg/g), however, displayed morphological deformities at pupation. That cytochrome P450 monooxygenases (P450s) are involved in the bioactivation of aflatoxin in this species is evidenced by the effects of piperonyl butoxide (PBO), a known P450 inhibitor, on toxicity; whereas no fourth instars pupated in the presence of 1 mug/g AFB1 in the diet, the presence of 0.1% PBO increased the pupation rate to 71.7%. Pupation rates of both fourth and fifth instars on diets containing 1 mug/g AFB1 also increased significantly in the presence of PBO. Effects of phenobarbital, a P450 inducer, on AFB1 toxicity were less dramatic than those of PBO. Collectively, these findings indicate that, as in many other vertebrates and invertebrates, toxicity of AFB1 to H. zea results from P450-mediated metabolic bioactivation.

Induction of two cytochrome P450 genes, Cyp6a2 and Cyp6a8, of Drosophila melanogaster by caffeine in adult flies and in cell culture

To examine whether caffeine, the most widely used xenobiotic compound, would induce insect cytochrome P450 or CYP gene expression, upstream DNA fragments of Cyp6a2 (0.12, 0.26, 0.52 and 0.98-kb) and Cyp6a8 (0.06, 0.1, 0.2, 0.5 and 0.8-kb) genes of Drosophila melanogaster were individually fused to the firefly luciferase (luc) reporter gene. Promoter activities of these constructs were examined in Drosophila SL-2 cells using luciferase assays. Activity of 0.2- and 0.8-kb upstream DNA of Cyp6a8 was also measured in transgenic female flies. When these flies were treated with 2 mM pure caffeine or Vivarin caffeine, both DNA fragments showed a 4-5-fold induction of promoter activity. Endogenous Cyp6a8 and Cyp6a2 genes in these flies also showed caffeine-induced expression. In addition, both 0.2- and 0.8-kb DNAs showed differential basal and caffeine-induced activity in head, ovaries, gut, cuticle plus fat body and malpighian tubules. However, in all tissues 0.8-kb DNA always showed higher basal and caffeine-induced activities compared to the 0.2-kb DNA, suggesting that the additional DNA present in the 0.8-kb fragment has sequences that enhance both activities. In SL-2 cells, all reporter constructs of each Cyp6 gene showed significantly higher basal activity than the empty vector. Sequences that boost basal activity are located in -265/-129 and -983/-522 DNA of Cyp6a2, and -199/-109 and -491/-199 DNA of Cyp6a8 genes. While the 0.12- and 0.1-kb upstream DNAs of Cyp6a2 and Cyp6a8 genes respectively did not show caffeine-inducibility in SL-2 cells, the longest upstream DNA of each gene gave the highest level of induction. Caffeine-responsive sequences are not clustered at one place; they appear to be dispersed in -983/-126 and -761/-109 regions of Cyp6a2 and Cyp6a8 genes which also contain many binding sites for activator protein 1 (AP1) and cyclic AMP response element binding protein (CRE-BP). Significance of these binding sites in caffeine-inducibility has been discussed.

Gene expression profiles of Drosophila melanogaster exposed to an insecticidal extract of Piper nigrum

Black pepper, Piper nigrum L. (Piperaceae), has insecticidal properties and could potentially be utilized as an alternative to synthetic insecticides. Piperine extracted from P. nigrum has a biphasic effect upon cytochrome P450 monooxygenase activity with an initial suppression followed by induction. In this study, an ethyl acetate extract of P. nigrum seeds was tested for insecticidal activity toward adult Musca domestica and Drosophila melanogaster. The effect of this same P. nigrum extract upon differential gene expression in D. melanogaster was investigated using cDNA microarray analysis of 7380 genes. Treatment of D. melanogaster with P. nigrum extract led to a greater than 2-fold upregulation of transcription of the cytochrome P450 phase I metabolism genes Cyp 6a8, Cyp 9b2, and Cyp 12d1 as well as the glutathione-S-transferase phase II metabolism gene Gst-S1. These data suggests a complex effect of P. nigrum upon toxin metabolism.

Cloning, Functional Expression, and Characterization of CYP709C1, the First Sub-terminal Hydroxylase of Long Chain Fatty Acid in Plants

We cloned and characterized CYP709C1, a new plant cytochrome P450 belonging to the P450 family, that so far has no identified function except for clustering with a fatty acid metabolizing clade of P450 enzymes. We showed here that CYP709C1 is capable of hydroxylating fatty acids at the omega-1 and omega-2 positions. This work was performed after recoding and heterologous expression of a full-length cDNA isolated from a wheat cDNA library in an engineered yeast strain. Investigation on substrate specificity indicates that CYP709C1 metabolizes different fatty acids varying in their chain length (C12 to C18) and unsaturation. CYP709C1 is the first identified plant cytochrome P450 that can catalyze sub-terminal hydroxylation of C18 fatty acids. cis-9,10-Epoxystearic acid is metabolized with the highest efficiency, i.e. K((m)(app)) of 8 microM and V(max(app)) of 328 nmol/min/nmol P450. This, together with the fact that wheat possesses a microsomal peroxygenase able to synthesize this compound from oleic acid, strongly suggests that it is a physiological substrate. Hydroxylated fatty acids are implicated in plant defense events. We postulated that CYP709C1 could be involved in plant defense by producing such compounds. This receives support from the observation that (i) sub-terminal hydroxylation of 9,10-epoxystearic acid is induced (15-fold after 3 h) in microsomes of wheat seedlings treated with the stress hormone methyl jasmonate and (ii) CYP709C1 is enhanced at the transcriptional level by this treatment. CYP709C1 transcript also accumulated after treatment with a combination of the safener naphthalic acid anhydride and phenobarbital. This indicates a possible detoxifying function for CYP709C1 that we discussed.

Expression of Cyp6g1 and Cyp12d1 in DDT resistant and susceptible strains of Drosophila melanogaster

The Rst(2)DDT locus (loci) in Drosophila is associated with the over-expression of two cytochrome P450 genes, Cyp6g1 and Cyp12d1. Using northern and western blot analysis we observed the expression pattern of these two genes in two DDT susceptible (Canton-S and 91-C) and three DDT resistant strains (Wisconsin, 91-R and Hikone-R). In Canton-S and 91-R, the CYP6G1 protein was constitutively expressed throughout development. In the Wisconsin strain, CYP6G1 was not expressed in third instar larvae unless the larvae are exposed to DDT. CYP12D1 protein was only expressed in adults. Cyp12d1 mRNA is induced in DDT resistant strains post-exposure to DDT and the expression patterns of Cyp12d1 mRNA varied across DDT resistant strains. Our data support the hypothesis that there is evolutionary plasticity in the expression patterns of P450s associated with metabolic pesticide resistance.

Profiling of abundant proteins associated with dichlorodiphenyltrichloroethane resistance inDrosophila melanogaster

Dichlorodiphenyltrichloroethane (DDT) metabolism-based resistance in Drosophila melanogaster is a complex metabolic system associated with the transcription of detoxification related genes, ion transport, lipid and sugar metabolism pathways. However, little is known about the differences regarding the proteome of field- and laboratory-selected resistant Drosophila genotypes. We investigated the impact of DDT resistance in the abundant proteome of field- and laboratory- selected resistant Drosophila using a two-dimensional gel electrophoresis DDT reference map. Proteomic profiling was performed in two DDT susceptible genotypes (Canton-S and 91-C) and three DDT resistant lines (Rst(2)DDT(91-R), Rst(2)DDT(Wisconsin) and Rst(2)DDT(Hikone-R)). Protein spots were stained with Coomassie blue and compared using PDQuest software. Selected protein spots were cut out and analyzed using matrix assisted laser desorption-time of flight mass spectrometry. Querying the NCBInr. 10.21.2003 database with mass spectrometric data yielded the identity of 21 differentially translated proteins in Rst(2)DDT(91-R), Rst(2)DDT(Wisconsin) and Canton-S representing proteins putatively involved in biochemical pathways such as glycolysis and gluconeogenesis, the pentose phosphate pathway, the Krebs cycle and fatty acid oxidation. We hypothesize that both strategies are aimed to use of the pentose phosphate pathway to increase glucose utilization while Rst(2)DDT(91-R) relies primarily on glycolysis to produce reduced NADP and increase DDT detoxification. DDT exposure in Canton-S induced six proteins, while four proteins were repressed in Rst(2)DDT(Hikone-R). Our data suggest that insecticide resistance appears to impact different metabolic pathways in Drosophila genotypes selected with the same pesticide (DDT).

Drosophila melanogaster CYP6A8, an insect P450 that catalyzes lauric acid (ω-1)-hydroxylation

Only a handful of P450 genes have been functionally characterized from the approximately 90 recently identified in the genome of Drosophila melanogaster. Cyp6a8 encodes a 506-amino acid protein with 53.6% amino acid identity with CYP6A2. CYP6A2 has been shown to catalyze the metabolism of several insecticides including aldrin and heptachlor. CYP6A8 is expressed at many developmental stages as well as in adult life. CYP6A8 was produced in Saccharomyces cerevisiae and enzymatically characterized after catalytic activity was reconstituted with D. melanogaster P450 reductase and NADPH. Although several saturated or non-saturated fatty acids were not metabolized by CYP6A8, lauric acid (C12:0), a short-chain unsaturated fatty acid, was oxidized by CYP6A8 to produce 11-hydroxylauric acid with an apparent V(max) of 25 nmol/min/nmol P450. This is the first report showing that a member of the CYP6 family catalyzes the hydroxylation of lauric acid. Our data open new prospects for the CYP6 P450 enzymes, which could be involved in important physiological functions through fatty acid metabolism.

Conserved regulatory elements in the promoters of two allelochemical-inducible cytochrome P450 genes differentially regulate transcription

CYP6B4, a cytochrome P450 gene from the tiger swallowtail Papilio glaucus, is transcriptionally induced in the midgut by dietary furanocoumarins, plant allelochemicals that can crosslink DNA in their UV-activated form. The CYP6B4 promoter contains an overlapping EcRE/ARE/XRE-xan element similar to that used for basal and xanthotoxin-inducible expression of the CYP6B1 promoter from the black swallowtail Papilio polyxenes. Transfection of the CYP6B4 promoter:CAT reporter construct into Sf9 cells demonstrates that the basal and xanthotoxin-inducible expression levels observed reflect the relative expression levels of this gene in the midguts of tiger swallowtail larvae. Transfections of mutant CYP6B4 promoter constructs into Sf9 cells indicate that the EcRE/ARE/XRE-xan element is necessary for CYP6B4 induction by xanthotoxin but not for its minimal basal expression. In addition to these elements, the CYP6B4 and CYP6B1 promoters also contain putative XRE-AhR elements identical to the aryl hydrocarbon response elements present in mammalian phase I detoxification genes. Transfections of CYP6B4 and CYP6B1 promoters containing EcRE/ARE/XRE-xan and XRE-AhR elements indicate that both are induced significantly by benzo(alpha)pyrene, an aryl hydrocarbon widespread in the environment, as well as by xanthotoxin, an allelochemical encountered in their hostplants.

Genetic mapping of genes conferring permethrin resistance in the malaria vector, Anopheles gambiae

Resistance to permethrin in an East African population of the major malaria vector, Anopheles gambiae is multifactorial. A mutated sodium channel allele and enhanced insecticide metabolism contribute to the resistance phenotype. We used microsatellite markers to scan the genome for quantitative trait loci (QTL) associated with permethrin resistance. Two major and one minor QTL were identified. The first QTL, rtp1, colocalizes with the sodium channel gene on chromosome 2L thus further supporting the importance of mutations in this gene in conferring permethrin resistance. The second two loci are located on the third chromosome and one of these, rtp2, flanks a large cluster of cytochrome P450 genes. Further detailed mapping of these regions will help elucidate the molecular mechanisms of metabolic resistance to insecticides.

The molecular basis of insecticide resistance in mosquitoes

Insecticide resistance is an inherited characteristic involving changes in one or more insect gene. The molecular basis of these changes are only now being fully determined, aided by the availability of the Drosophila melanogaster and Anopheles gambiae genome sequences. This paper reviews what is currently known about insecticide resistance conferred by metabolic or target site changes in mosquitoes.

Genome-wide transcription profile of field- and laboratory-selected dichlorodiphenyltrichloroethane (DDT)-resistant Drosophila

Genome-wide microarray analysis (Affymetrix array) was used (i) to determine whether only one gene, the cytochrome P450 enzyme Cyp6g1, is differentially transcribed in dichlorodiphenyltrichloroethane (DDT)-resistant vs. -susceptible Drosophila; and (ii) to profile common genes differentially transcribed across a DDT-resistant field isolate [Rst(2)DDT(Wisconsin)] and a laboratory DDT-selected population [Rst(2)DDT(91-R)]. Statistical analysis (ANOVA model) identified 158 probe sets that were differentially transcribed among Rst(2)DDT(91-R), Rst(2)DDT(Wisconsin), and the DDT-susceptible genotype Canton-S (P < 0.01). The cytochrome P450 Cyp6a2 and the diazepam-binding inhibitor gene (Dbi) were over transcribed in the two DDT-resistant genotypes when compared to the wild-type Drosophila, and this difference was significant at the most stringent statistical level, a Bonferroni correction. The list of potential candidates differentially transcribed also includes 63 probe sets for which molecular function ontology annotation of the probe sets did not exist. A total of four genes (Cyp6a2, Dbi, Uhg1, and CG11176) were significantly different (P < 5.6 e(-06)) between Rst(2)DDT(91-R) and Canton-S. Additionally, two probe sets encoding Cyp12d1 and Dbi were significantly different between Rst(2)DDT(Wisconsin) and Canton-S after a Bonferroni correction. Fifty-two probe sets, including those associated with pesticide detoxification, ion transport, signal transduction, RNA transcription, and lipid metabolism, were commonly expressed in both resistant lines but were differentially transcribed in Canton-S. Our results suggest that more than Cyp6g1 is overtranscribed in field and laboratory DDT-resistant genotypes, and the number of commonalities suggests that similar resistance mechanisms may exist between laboratory- and field-selected DDT-resistant fly lines.

Microarray analysis of cytochrome P450 mediated insecticide resistance in Drosophila

Insecticide resistance in laboratory selected Drosophila strains has been associated with upregulation of a range of different cytochrome P450s, however in recent field isolates of D. melanogaster resistance to DDT and other compounds is conferred by one P450 gene, Cyp6g1. Using microarray analysis of all Drosophila P450 genes, here we show that different P450 genes such as Cyp12d1 and Cyp6a8 can also be selected using DDT in the laboratory. We also show, however, that a homolog of Cyp6g1 is over-expressed in a field resistant strain of D. simulans. In order to determine why Cyp6g1 is so widely selected in the field we examine the pattern of cross-resistance of both resistant strains and transgenic flies over-expressing Cyp6g1 alone. We show that all three DDT selected P450s can confer resistance to the neonicotinoid imidacloprid but that Cyp6a8 confers no cross-resistance to malathion. Transgenic flies over-expressing Cyp6g1 also show cross-resistance to other neonicotinoids such as acetamiprid and nitenpyram. We suggest that the broad level of cross-resistance shown by Cyp6g1 may have facilitated its selection as a resistance gene in natural Drosophila populations.

Forty years of solitude: life-history divergence and behavioural isolation between laboratory lines of Drosophila melanogaster

The study of the early stages of speciation can benefit from examination of differences between populations of known history that have been separated for a short time, such as a few thousands of generations. We asked whether two lines of Drosophila melanogaster that were isolated more than 40 years ago have evolved differences in life-history characters, or have begun to evolve behavioural or postzygotic isolation. One line, which is resistant to DDT, showed lower egg production and a shorter lifespan than a susceptible line. These differences are not a pleiotropic effect of resistance because they are not attributable to the chromosome that contains the resistance factors. The two lines have begun to become behaviourally isolated. Again, the isolation is not attributable to genes on the chromosome that contains resistance factors. The lines show only prezygotic isolation; there is no evidence of reduced fitness of F1 or F2 hybrids. These lines and others like them, should be excellent subjects for analyses of genetic changes that could lead to speciation.

Molecular analysis of multiple cytochrome P450 genes from the malaria vector, Anopheles gambiae

Cytochrome P450s are a superfamily of haemoproteins, important in the metabolism of endogenous compounds and xenobiotics. As a first step to elucidating the role of this family in insecticide resistance in the malaria mosquito, Anopheles gambiae, we have cloned and mapped multiple P450 genes. Sixteen cDNAs encoding full-length P450s were cloned and physically mapped to the mosquito's polytene chromosomes. Fourteen of these encode putative CYP6 proteins and two encode P450s belonging to the CYP9 class. Eighteen new A. gambiae Cyp4 P450 genes were identified using degenerate PCR primers, cDNAs were detected for ten and in situ locations for thirteen members of this gene family.

Cyp6a8 of Drosophila melanogaster: gene structure, and sequence and functional analysis of the upstream DNA

In Drosophila, the insecticide resistant 91-R strain is an overproducer and susceptible 91-C and ry(506) strains are the underproducers of CYP6A8 mRNA encoded by a cytochrome P450 gene, Cyp6a8. Low expression of Cyp6a8 in the underproducer strains is due to a downregulatory effect of a putative repressor locus, which is thought to be mutant in the overproducer strain. In the present investigation, organization of Cyp6a8 and promoter activity of its upstream DNA were analyzed. Cyp6a8 has two introns of which intron II is similar to the introns of other insect CYP genes with respect to its length and position. Intron I is only 36 bp long and lacks consensus splice sites. It is also in-frame with the CYP6A8 open reading frame. Therefore, inefficient splicing of intron I may produce two isoforms of CYP6A8. Analysis of Cyp6a8 upstream DNA of the overproducer 91-R strain showed that DNA sequences between -199 and -761 bp are required for the highest constitutive and barbital-induced expression of Cyp6a8. This region has six barbie boxes and binding sites for various transcription factors. Promoter activity of the -11/-761 DNA of the overproducer 91-R strain was found to be 4-fold lower in the genome of underproducer ry(506) strain, which is wild type for the putative repressor gene, than in the genome of F1 hybrids of 91-R and ry(506) strains. These results suggest that -11/-761 Cyp6a8 DNA of the 91-R strain can respond to the active repressor present in the hybrid genome and further support our previous findings that overexpression of Cyp6a8 is a result of mutation of a repressor gene rather than mutation of the cis-regulatory sequences.

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

Previous studies have shown that the major metabolic resistance locus in the insecticide-resistant Drosophila line Rst(2)DDT(Wisconsin) maps between the markers cn and vg on chromosome 2. Six cytochrome P450 genes exist in this region. We investigated the expression levels of these P450 genes in DDT-resistant and -susceptible fly lines. We report: (i) DDT resistance is significant (> 30-fold) and dominant, (ii) resistance is reduced by the cytochrome P450 inhibitor PBO, (iii) there is constitutive over-expression relative to susceptible flies of two genes encoding cytochrome P450 enzymes within the cn-vg region (CYP6G1 = 4.3-fold; CYP12D1 = 6.0-fold), and (iv) exposure to DDT results in an increased expression of only one of these two P450 genes (CYP12D1 > or = 6-fold above constitutive resistant fly baselines).

Overproduction of a P450 that metabolizes diazinon is linked to a loss-of-function in the chromosome 2 ali-esterase (MdalphaE7) gene in resistant house flies

Up-regulation of detoxifying enzymes in insecticide-resistant strains of the house fly is a common mechanism for metabolic resistance. However, the molecular basis of this increased insecticide metabolism is not well understood. In the multiresistant Rutgers strain, several cytochromes P450 and glutathione S-transferases are constitutively overexpressed at the transcriptional level. Overexpression is the result of trans-regulation, and a regulatory gene has been located on chromosome 2. A Gly137 to Asp point mutation in alphaE7 esterase gene, leading to the loss of carboxylesterase activity, has been associated with organophosphate resistance in the house fly and the sheep blowfly. We show here that purified recombinant CYP6A1 is able to detoxify diazinon with a high efficiency. We also show that either the Gly137 to Asp point mutation in alphaE7 esterase gene or a deletion at this locus confer resistance and overproduction of the CYP6A1 protein. Based on these findings, we propose it is the absence of the wild-type Gly137 allele of the alphaE7 gene that releases the transcriptional repression of genes coding for detoxification enzymes such as CYP6A1, thereby leading to metabolic resistance to diazinon.