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

Novel cytochrome P450s, CYP6BB1 and CYP6P10, from the salt marsh mosquito Aedes sollicitans (Walker) (Diptera: Culicidae)

Based on the conserved heme-binding region and the charge pair consensus of insect cytochrome P450s, two novel full-length P450 cDNAs, CYP6BB1 and CYP6P10, were cloned from the salt marsh mosquito Aedes sollicitans (Walker). CYP6BB1 and CYP6P10 had open reading frames of 1,518 and 1,521 nucleotides encoding 506 and 507 amino acid residue proteins, respectively. Several alleles with amino acid substitutions were found both in CYP6BB1 and CYP6P10. The deduced proteins are typical microsomal P450s sharing signature sequences with other insect CYP6 P450s. Sequence analysis showed that both CYP6BB1 and CYP6P10 shared highest sequence identities with P450 CYP6P4, 56% and 65%, respectively. Phylogenetic analysis showed both CYP6BB1 and CYP6P10 were grouped into the clade containing several P450s from subfamily CYP6P. Real-time RT-PCR analysis showed CYP6BB1 but not CYP6P10 transcription in females was significantly increased 24 h after a blood meal. Neither CYP6BB1 nor CYP6P10 were life stage or gender specific. Protein expression experiments are needed to determine the functions of these proteins.

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.

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.

Molecular marker systems in insects: current trends and future avenues

Insects comprise the largest species composition in the entire animal kingdom and possess a vast undiscovered genetic diversity and gene pool that can be better explored using molecular marker techniques. Current trends of application of DNA marker techniques in diverse domains of insect ecological studies show that mitochondrial DNA (mtDNA), microsatellites, random amplified polymorphic DNA (RAPD), expressed sequence tags (EST) and amplified fragment length polymorphism (AFLP) markers have contributed significantly for progresses towards understanding genetic basis of insect diversity and for mapping medically and agriculturally important genes and quantitative trait loci in insect pests. Apart from these popular marker systems, other novel approaches including transposon display, sequence-specific amplification polymorphism (S-SAP), repeat-associated polymerase chain reaction (PCR) markers have been identified as alternate marker systems in insect studies. Besides, whole genome microarray and single nucleotide polymorphism (SNP) assays are becoming more popular to screen genome-wide polymorphisms in fast and cost effective manner. However, use of such methodologies has not gained widespread popularity in entomological studies. The current study highlights the recent trends of applications of molecular markers in insect studies and explores the technological advancements in molecular marker tools and modern high throughput genotyping methodologies that may be applied in entomological researches for better understanding of insect ecology at molecular level.

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.

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.

Evidence for a functional interaction between the Bari1 transposable element and the cytochrome P450 cyp12a4 gene in Drosophila melanogaster

Previous studies of the genomic distribution of the transposon Bari1 in Drosophila melanogaster have revealed an element which is fixed at division 91F in over 90 lab and natural populations. Here we report about the structural and transcriptional features of the insertion site which was studied in sublines isolated from an exceptional Drosophila line polymorphic for the presence/absence of Bari1 at 91F. The insert is located at the 3' end of the cyp12a4 gene that belongs to the cytochrome P450 family. In flies with the insert the transcript of this gene encompasses 18 nucleotides of the transposon, it is shorter and is about tenfold more abundant compared to flies devoid of it. Although the hypothetical selective agent remains unknown, these data are suggestive of a selective advantage brought about by the Bari1 insert and are reminiscent of recent evidence for functional mutagenesis of cyp6g1, another P450 gene, brought about by Accord and Doc transposable elements in D. melanogaster and Drosophila simulans.

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

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.

Point mutations associated with insecticide resistance in the Drosophila cytochrome P450 Cyp6a2 enable DDT metabolism

Three point mutations R335S, L336V and V476L, distinguish the sequence of a cytochrome P450 CYP6A2 variant assumed to be responsible for 1,1,1-trichloro-2,2-bis-(4'-chlorophenyl)ethane (DDT) resistance in the RDDT(R) strain of Drosophila melanogaster. To determine the impact of each mutation on the function of CYP6A2, the wild-type enzyme (CYP6A2wt) of Cyp6a2 was expressed in Escherichia coli as well as three variants carrying a single mutation, the double mutant CYP6A2vSV and the triple mutant CYP6A2vSVL. All CYP6A2 variants were less stable than the CYP6A2wt protein. Two activities enhanced in the RDDT(R) strain were measured with all recombinant proteins, namely testosterone hydroxylation and DDT metabolism. Testosterone was hydroxylated at the 2beta position with little quantitative variation among the variants. In contrast, metabolism of DDT was strongly affected by the mutations. The CYP6A2vSVL enzyme had an enhanced metabolism of DDT, producing dicofol, dichlorodiphenyldichloroethane and dichlorodiphenyl acetic acid. The apparent affinity of the enzymes CYP6A2wt and CYP6A2vSVL for DDT and testosterone was not significantly different as revealed by the type I difference spectra. Sequence alignments with CYP102A1 provided clues to the positions of the amino acids mutated in CYP6A2. These mutations were found spatially clustered in the vicinity of the distal end of helix I relative to the substrate recognition valley. Thus this area, including helix J, is important for the structure and activity of CYP6A2. Furthermore, we show here that point mutations in a cytochrome P450 can have a prominent role in insecticide resistance.

Transcriptome analysis of the cowpea weevil bruchid: identification of putative proteinases and alpha-amylases associated with food breakdown

We describe here the first systematic work to discover insect genes involved in food breakdown using a cDNA library enriched for gut-expressed transcripts from Callosobruchus maculatus. A total of 1056 clones were screened for cDNA insert-containing plasmids, and 503 nonredundant open reading frames were discovered. Twenty-three inferred genes potentially involved in digestive processes in cowpea weevil were identified, including proteinases and amylases. The predicted catalytic sites were identified in the inferred cysteine and aspartic acid proteinases, and in alpha-amylases. Transcriptome analysis of the cowpea bruchid will potentially permit gene discovery in other beetles, an insect order of major economic and ecological importance that is poorly represented in genomic databases.

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.