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

Effects of Amino Acid Mutation in Cytochrome P450 (CYP96A146) of Descurainia sophia on the Metabolism and Resistance to Tribenuron-Methyl

Cytochrome P450 monooxygenases (P450s) play important roles in herbicide resistance. In this study, there are four amino acid mutations (F39Y, H163Y, S203A, and V361E) between CYP96A146-S and CYP96A146-R, which were cloned, respectively, from susceptible (S) and tribenuron-methyl-resistant (TR) Descurainia sophia. The Arabidopsis expressing CYP96A146-S or CYP96A146-R showed resistance to tribenuron-methyl, carfentrazone-ethyl, and oxyfluorfen, while Arabidopsis transformed with CYP96A146-R or CYP96A146 with any two or three mutations of H163Y, S203A, or V361E exhibited significantly higher resistance to tribenuron-methyl than Arabidopsis expressing CYP96A146-S. The metabolic rates of tribenuron-methyl were significantly faster in Arabidopsis expressing CYP96A146-R than that with CYP96A146-S. The molecular dynamics simulation demonstrated that amino acid mutations did not affect the domain of the HEM ring, which could significantly enhance the volume of the catalytic pocket in P450 (CYP96A146), thereby increasing the collision rate between the catalytic pocket and tribenuron-methyl. Hence, the amino acid mutations may be one of the mechanisms underlying P450-mediated herbicide resistance.

Two highly selected mutations in the tandemly duplicated CYP6P4a and CYP6P4b genes drive pyrethroid resistance in Anopheles funestus in West Africa

BACKGROUND

Gaining a comprehensive understanding of the genetic mechanisms underlying insecticide resistance in malaria vectors is crucial for optimising the effectiveness of insecticide-based vector control methods and developing diagnostic tools for resistance management. Considering the heterogeneity of metabolic resistance in major malaria vectors, the implementation of tailored resistance management strategies is essential for successful vector control. Here, we provide evidence demonstrating that two highly selected mutations in CYP6P4a and CYP6P4b are driving pyrethroid insecticide resistance in the major malaria vector Anopheles funestus, in West Africa.

RESULTS

Continent-wide polymorphism survey revealed escalated signatures of directional selection of both genes between 2014 and 2021. In vitro insecticide metabolism assays with recombinant enzymes from both genes showed that mutant alleles under selection exhibit higher metabolic efficiency than their wild-type counterparts. Using the GAL4-UAS expression system, transgenic Drosophila flies overexpressing mutant alleles exhibited increased resistance to pyrethroids. These findings were consistent with in silico predictions which highlighted changes in enzyme active site architecture that enhance the affinity of mutant alleles for type I and II pyrethroids. Furthermore, we designed two DNA-based assays for the detection of CYP6P4a-M220I and CYP6P4b-D284E mutations, showing their current confinement to West Africa. Genotype/phenotype correlation analyses revealed that these markers are strongly associated with resistance to types I and II pyrethroids and combine to drastically reduce killing effects of pyrethroid bed nets.

CONCLUSIONS

Overall, this study demonstrated that CYP6P4a and CYP6P4b contribute to pyrethroid resistance in An. funestus and provided two additional insecticide resistance molecular diagnostic tools that would contribute to monitoring and better management of resistance.

Transcriptome Analysis Unveils Molecular Mechanisms of Acaricide Resistance in Two-Spotted Spider Mite Populations on Hops

Broad-spectrum crop protection technologies, such as abamectin and bifenthrin, are globally relied upon to curb the existential threats from economic crop pests such as the generalist herbivore Tetranychus urticae Koch (TSSM). However, the rising cost of discovering and registering new acaricides, particularly for specialty crops, along with the increasing risk of pesticide resistance development, underscores the urgent need to preserve the efficacy of currently registered acaricides. This study examined the overall genetic mechanism underlying adaptation to abamectin and bifenthrin in T. urticae populations from commercial hop fields in the Pacific Northwestern region of the USA. A transcriptomic study was conducted using four populations (susceptible, abamectin-resistant, and two bifenthrin-resistant populations). Differential gene expression analysis revealed a notable disparity, with significantly more downregulated genes than upregulated genes in both resistant populations. Gene ontology enrichment analysis revealed a striking consistency among all three resistant populations, with downregulated genes predominately associated with chitin metabolism. In contrast, upregulated genes in the resistant populations were linked to biological processes, such as peptidase activity and oxidoreductase activity. Proteolytic activity by peptidase enzymes in abamectin- and bifenthrin-resistant TSSM populations may suggest their involvement in acaricide metabolism. These findings provide valuable insights into the molecular mechanisms underlying acaricide resistance in the TSSM. This knowledge can be utilized to develop innovative pesticides and molecular diagnostic tools for effectively monitoring and managing resistant TSSM populations.

Sesamin Exerts an Antioxidative Effect by Activating the Nrf2 Transcription Factor in the Glial Cells of the Central Nervous System in Drosophila Larvae

Sesame seeds are abundant in sesamin, which exerts health-promoting effects such as extending the lifespan of adult Drosophila and suppressing oxidative stress by activating the Nrf2 transcription factor. Here, we investigated whether sesamin activated Nrf2 in larval tissues and induced the expression of Nrf2 target genes. In the sesamin-fed larvae, Nrf2 was activated in the central nervous system (CNS), gut, and salivary glands. The ectopic expression of Keap1 in glial cells inhibited sesamin-induced Nrf2 activation in the whole CNS more than in the neurons, indicating that sesamin activates Nrf2 in glia efficiently. We labeled the astrocytes as well as cortex and surface glia with fluorescence to identify the glial cell types in which Nrf2 was activated; we observed their activation in both cell types. These data suggest that sesamin may stimulate the expression of antioxidative genes in glial cells. Among the 17 candidate Nrf2 targets, the mRNA levels of Cyp6a2 and Cyp6g1 in cytochrome P450 were elevated in the CNS, gut, and salivary glands of the sesamin-fed larvae. However, this elevation did not lead to resistance against imidacloprid, which is detoxified by these enzymes. Our results suggest that sesamin may exert similar health-promoting effects on the human CNS and digestive tissues.

Transcription factor B-H2 regulates CYP9J34 expression conveying deltamethrin resistance in Culex pipiens pallens

BACKGROUND

Mosquitoes are vectors of various diseases, posing significant health threats worldwide. Chemical pesticides, particularly pyrethroids like deltamethrin, are commonly used for mosquito control, but the emergence of resistant mosquito populations has become a concern. In the deltamethrin-resistant (DR) strain of Culex pipiens pallens, the highly expressed cytochrome P450 9 J34 (CYP9J34) gene is believed to play a role in resistance, yet the underlying mechanism remains unclear.

RESULTS

Quantitative polymerase chain reaction with reverse transcription (qRT-PCR) analysis revealed that the expression of CYP9J34 was 14.6-fold higher in DR strains than in deltamethrin-susceptible (DS) strains. The recombinant production of CYP9J34 protein of Cx. pipiens pallens showed that the protein could directly metabolize deltamethrin, yielding the major metabolite 4'-OH deltamethrin. Through dual luciferase reporter assays and RNA interference, the transcription factor homeobox protein B-H2-like (B-H2) was identified to modulate the expression of the CYP9J34 gene, contributing to mosquito resistance to deltamethrin.

CONCLUSIONS

Our findings demonstrate that the CYP9J34 protein could directly degrade deltamethrin, and the transcription factor B-H2 could regulate CYP9J34 expression, influencing the resistance of mosquitoes to deltamethrin. © 2023 Society of Chemical Industry.

Genetics, genomics and mechanisms responsible for high levels of pyrethroid resistance in Musca domestica

Insecticide resistance is both economically important and evolutionarily interesting phenomenon. Identification of the mutations responsible for resistance allows for highly sensitive resistance monitoring and allows tools to study the forces (population genetics, fitness costs, etc.) that shape the evolution of resistance. Genes coding for insecticide targets have many well-characterized mutations, but the mutations responsible for enhanced detoxification have proven difficult to identify. We employed multiple strategies to identify the mutations responsible for the extraordinarily high permethrin resistance in the KS17-R strain of house fly (Musca domestica): insecticide synergist assays, linkage analysis, bulk segregant analyses (BSA), transcriptomics and long read DNA (Nanopore) sequencing. The >85,100-fold resistance in KS17-R was partially suppressed by the insecticide synergists piperonyl butoxide and S,S,S-tributylphosphorothionate, but not by diethyl maleate nor by injection. This suggests the involvement of target site insensitivity, CYP-mediated resistance, possibly hydrolase mediated resistance and potentially other unknown factors. Linkage analysis identified chromosomes 1, 2, 3 and 5 as having a role in resistance. BSA mapped resistance loci on chromosomes 3 and 5. The locus on chromosome 3 was centered on the voltage sensitive sodium channel. The locus on chromosome 5 was associated with a duplication of multiple detoxification genes. Transcriptomic analyses and long read DNA sequencing revealed overexpressed CYPs and esterases and identified a complex set of structural variants at the chromosome 5 locus.

The effects of methylphenidate and atomoxetine on Drosophila brain at single-cell resolution and potential drug repurposing for ADHD treatment

The stimulant methylphenidate (MPH) and the non-stimulant atomoxetine (ATX) are frequently used for the treatment of attention-deficit/hyperactivity disorder (ADHD); however, the function of these drugs in different types of brain cells and their effects on related genes remain largely unknown. To address these questions, we built a pipeline for the simultaneous examination of the activity behavior and transcriptional responses of Drosophila melanogaster at single-cell resolution following drug treatment. We selected the Drosophila with significantly increased locomotor activities (hyperactivity-like behavior) following the administration of each drug in comparison with the control (same food as the drug-treated groups with 5% sucrose, yeast, and blue food dye solution) using EasyFlyTracker. Subsequently, single cell RNA sequencing (scRNASEQ) was used to capture the transcriptome of 82,917 cells, unsupervised clustering analysis of which yielded 28 primary cell clusters representing the major cell types in adult Drosophila brain. Indeed, both neuronal and glial cells responded to MPH and ATX. Further analysis of differentially expressed genes (DEGs) revealed distinct transcriptional changes associated with these two drugs, such as two well-studied dopamine receptor genes (Dop2R and DopEcR) were responsive to MPH but not to ATX at their optimal doses, in addition to genes involved in dopamine metabolism pathways such as Syt1, Sytalpha, Syt7, and Ih in different cell types. More importantly, MPH also suppressed the expression of genes encoding other neurotransmitter receptors and synaptic signaling molecules in many cell types, especially those for Glu and GABA, while the responsive effects of ATX were much weaker. In addition to monoaminergic neuronal transmitters, other neurotransmitters have also shown a similar pattern with respect to a stronger effect associated with MPH than with ATX. Moreover, we identified four distinct glial cell subtypes responsive to the two drugs and detected a greater number of differentially expressed genes associated with ensheathing and astrocyte-like glia. Furthermore, our study provides a rich resource of candidate target genes, supported by drug set enrichment analysis (P = 2.10E-4; hypergeometric test), for the further exploration of drug repurposing. The whole list of candidates can be found at ADHDrug ( http://adhdrug.cibr.ac.cn/ ). In conclusion, we propose a fast and cost-efficient pipeline to explore the underlying molecular mechanisms of ADHD drug treatment in Drosophila brain at single-cell resolution, which may further facilitate drug repurposing applications.

The role of P450 enzymes in malaria and other vector-borne infectious diseases

Vector-borne infectious diseases are still an important global health problem. Malaria is the most important among them, mainly pediatric, life-threatening disease. Malaria and other vector-borne disorders caused by parasites, bacteria, and viruses have a strong impact on public health and significant economic costs. Most vector-borne diseases could be prevented by vector control, with attention to the ecological and biodiversity conservation aspects. Chemical control with pesticides and insecticides is widely used as a measure of prevention although increasing resistance to insecticides is a serious issue in vector control. Metabolic resistance is the most common mechanism and poses a big challenge. Insect enzyme systems, including monooxygenase CYP P450 enzymes, are employed by vectors mainly to metabolize insecticides thus causing resistance. The discovery and application of natural specific inhibitors/blockers of vector P450 enzymes as synergists for commonly used pesticides will contribute to the "greening" of insecticides. Besides vector CYPs, host CYP enzymes could also be exploited to fight against vector-borne diseases: using mostly their detoxifying properties and involvement in the immune response. Here, we review published research data on P450 enzymes from all players in vector-borne infections, that is, pathogens, vectors, and hosts, regarding the potential role of CYPs in disease. We discuss strategies on how to exploit cytochromes P450 in vector-borne disease control.

A glutamate-gated chloride channel as the mite-specific target-site of dicofol and other diphenylcarbinol acaricides

Dicofol has been widely used to control phytophagous mites. Although dicofol is chemically related to DDT, its mode of action has remained elusive. Here, we mapped dicofol resistance in the spider mite Tetranychus urticae to two genomic regions. Each region harbored a glutamate-gated chloride channel (GluCl) gene that contained a mutation-G314D or G326E-known to confer resistance against the unrelated acaricide abamectin. Using electrophysiology assays we showed that dicofol and other diphenylcarbinol acaricides-bromopropylate and chlorobenzilate-induce persistent currents in Xenopus oocytes expressing wild-type T. urticae GluCl3 receptors and potentiate glutamate responses. In contrast, the G326E substitution abolished the agonistic activity of all three compounds. Assays with the wild-type Drosophila GluClα revealed that this receptor was unresponsive to dicofol. Homology modeling combined with ligand-docking confirmed the specificity of electrophysiology assays. Altogether, this work elucidates the mode of action of diphenylcarbinols as mite-specific agonists of GluCl.

A single point mutation in the Bemisia tabaci cytochrome-P450 CYP6CM1 causes enhanced resistance to neonicotinoids

The tobacco whitefly, Bemisia tabaci, is a polyphagous crop pest which causes high levels of economic damage across the globe. Insecticides are often required for the effective control of this species, among which the neonicotinoid class have been particularly widely used. Deciphering the mechanisms responsible for resistance to these chemicals is therefore critical to maintain control of B. tabaci and limit the damage it causes. An important mechanism of resistance to neonicotinoids in B. tabaci is the overexpression of the cytochrome P450 gene CYP6CM1 which leads to the enhanced detoxification of several neonicotinoids. In this study we show that qualitative changes in this P450 dramatically alter its metabolic capacity to detoxify neonicotinoids. CYP6CM1 was significantly over-expressed in two strains of B. tabaci which displayed differing levels of resistance to the neonicotinoids imidacloprid and thiamethoxam. Sequencing of the CYP6CM1 coding sequence from these strains revealed four different alleles encoding isoforms carrying several amino acid changes. Expression of these alleles in vitro and in vivo provided compelling evidence that a mutation (A387G), present in two of the CYP6CM1 alleles, results in enhanced resistance to several neonicotinoids. These data demonstrate the importance of both qualitative and quantitative changes in genes encoding detoxification enzymes in the evolution of insecticide resistance and have applied implications for resistance monitoring programs.

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.

Cytochrome P450 polymorphism: From evolution to clinical use

The cytochromes P450s can be divided in two groups, those of high importance for endogenous functions being evolutionary quite stable and those participating in detoxification of drugs and other xenobiotics having less important endogenous functions. In the latter group extensive genetic diversity has been allowed and in addition this is of high importance for survival in different environments. The genetic polymorphisms in these genes have evolved to some extent based on dietary restrictions and environmental factors and have not been subject of conservation due to less importance for survival. In cases of high dietary selection events, gene multiplication and amplification events have been seen. The different variants in genes encoding drug metabolizing enzymes can be used as genetic biomarkers (pharmacogenomic labels) for adjustment of drug treatment leading to less adverse drug reactions and better response. Indeed, this has improved the use of personalized medicine, although the missing heredity seen based on twin studies indicates that there are indeed many more genetic variants to be discovered before one can achieve a satisfactory relationship between genotype and phenotype with respect to drug metabolism and toxicity.

The Role of Cytochrome P450s in Insect Toxicology and Resistance

Insect cytochrome P450 monooxygenases (P450s) perform a variety of important physiological functions, but it is their role in the detoxification of xenobiotics, such as natural and synthetic insecticides, that is the topic of this review. Recent advances in insect genomics and postgenomic functional approaches have provided an unprecedented opportunity to understand the evolution of insect P450s and their role in insect toxicology. These approaches have also been harnessed to provide new insights into the genomic alterations that lead to insecticide resistance, the mechanisms by which P450s are regulated, and the functional determinants of P450-mediated insecticide resistance. In parallel, an emerging body of work on the role of P450s in defining the sensitivity of beneficial insects to insecticides has been developed. The knowledge gained from these studies has applications for the management of P450-mediated resistance in insect pests and can be leveraged to safeguard the health of important beneficial insects.

Why are Drosophila larvae more sensitive to avermectin than adults?

The insects have different physiological and morphological characteristics in various developmental stages. The difference in the characteristics may be related to the different sensitivity of insects to insecticides. In avermectin resistant strain screening assay, we found that the Drosophila larvae displayed a higher sensitivity to the insecticidal effect of avermectin, compared with adults. In this study, we found that the Drosophila larvae have relatively thicker chitin layer, faster avermectin metabolism and lower P-glycoprotein (P-gp) level, when compared with the adults. Besides, the expression levels of the molecular targets of avermectin, glutamate-gated chloride channel and γ-aminobutyric acid (GABA)-gated chloride channel, are lower in the larval stage than the adult. These results suggested that lower P-gp level in the body especially in brain may be the major reason for the higher sensitivity of Drosophila larvae to the insecticide. In summary, these results shed new light on the concept that different developmental stages of insects display different sensitivity to the same insecticide, which also provided a physiological explanation of the relevant mechanism of the difference of sensitivity of insect at its larval and adult stages to insecticide.

Transcriptomic and proteomic analysis of pyrethroid resistance in the CKR strain of Aedes aegypti

Aedes aegypti is an important vector of human viral diseases. This mosquito is distributed globally and thrives in urban environments, making it a serious risk to human health. Pyrethroid insecticides have been the mainstay for control of adult A. aegypti for decades, but resistance has evolved, making control problematic in some areas. One major mechanism of pyrethroid resistance is detoxification by cytochrome P450 monooxygenases (CYPs), commonly associated with the overexpression of one or more CYPs. Unfortunately, the molecular basis underlying this mechanism remains unknown. We used a combination of RNA-seq and proteomic analysis to evaluate the molecular basis of pyrethroid resistance in the highly resistant CKR strain of A. aegypti. The CKR strain has the resistance mechanisms from the well-studied Singapore (SP) strain introgressed into the susceptible Rockefeller (ROCK) strain genome. The RNA-seq and proteomics data were complimentary; each offering insights that the other technique did not provide. However, transcriptomic results did not quantitatively mirror results of the proteomics. There were 10 CYPs which had increased expression of both transcripts and proteins. These CYPs appeared to be largely trans-regulated, except for some CYPs for which we could not rule out gene duplication. We identified 65 genes and lncRNAs as potentially being responsible for elevating the expression of CYPs in CKR. Resistance was associated with multiple loci on chromosome 1 and at least one locus on chromosome 3. We also identified five CYPs that were overexpressed only as proteins, suggesting that stabilization of CYP proteins could be a mechanism of resistance. Future studies to increase the resolution of the resistance loci, and to examine the candidate genes and lncRNAs identified here will greatly enhance our understanding of CYP-mediated resistance in A. aegypti.

The overexpression and variant of CYP6G4 associated with propoxur resistance in the housefly, Musca domestica L

BACKGROUND

The control of the housefly, Musca domestica, heavily relies on the application of insecticides. Propoxur, a carbamate, has been widely used for vector control. The housefly populations with high propoxur resistance display point mutations and overexpression of acetylcholinesterase. However, the roles of cytochrome P450 monoxygenases (P450s), as important detoxification enzymes, remain poorly understand in the housefly resistant to propoxur.

RESULTS

P450s were implied to contribute to propoxur resistance based on the synergism of piperonyl butoxide (PBO) and the increase of P450 enzyme activity in the near-isogenic line propoxur resistant strain (N-PRS). Five P450 genes (CYP6G4, CYP6A25, CYP304A1, CYP6D3, and CYP6A1) by RNA-sequencing comparison were significantly up-regulated in the N-PRS strain with >1035-fold resistance to propoxur. A total of 13 non-synonymous mutations of three P450 genes (CYP6G4, CYP6D3, and CYP6D8) were found in the N-PRS strain. The amino acid substitutions of CYP6D3 and CYP6D8 were probably not resistance-associated single nucleotide polymorphisms (SNPs) because they were also found in the aabys susceptible strain. However, CYP6G4 variant in the N-PRS strain was not found in the aabys strain. The conjoint analysis of mutations and a series of genetic crosses exhibited that the housefly propoxur resistance was strongly associated with the mutations of CYP6G4 gene.

CONCLUSION

Our results suggested that a combination of up-regulated transcript levels and mutations of CYP6G4 contributed to propoxur resistance in the housefly. © 2021 Society of Chemical Industry.

The contribution of detoxification pathways to pyrethroid resistance in Hyalella azteca

Chronic exposure to pyrethroid insecticides can result in strong selective pressures on non-target species in aquatic systems and drive the evolution of resistance and population-level changes. Characterizing the underlying mechanisms of resistance is essential to better understanding the potential consequences of contaminant-driven microevolution. The current study found that multiple mechanisms enhance the overall tolerance of Hyalella azteca to the pyrethroid permethrin. In H. azteca containing mutations in the voltage-gated sodium channel (VGSC), both adaptation and acclimation played a role in mitigating the adverse effects of pyrethroid exposures. Pyrethroid resistance is primarily attributed to the heritable mutation at a single locus of the VGSC, resulting in reduced target-site sensitivity. However, additional pyrethroid tolerance was conferred through enhanced enzyme-mediated detoxification. Cytochrome P450 monooxygenases (CYP450) and general esterases (GE) significantly contributed to the detoxification of permethrin in H. azteca. Over time, VGSC mutated H. azteca retained most of their pyrethroid resistance, though there was some increased sensitivity from parent to offspring when reared in the absence of pyrethroid exposure. Permethrin median lethal concentrations (LC50s) declined from 1809 ng/L in parent (P₀) individuals to 1123 ng/L in the first filial (F₁) generation, and this reduction in tolerance was likely related to alterations in acclimation mechanisms, rather than changes to target-site sensitivity. Enzyme bioassays indicated decreased CYP450 and GE activity from P₀ to F₁, whereas the VGSC mutation was retained. The permethrin LC50s in resistant H. azteca were still two orders-of-magnitude higher than non-resistant populations indicating that the largest proportion of resistance was maintained through the inherited VGSC mutation. Thus, the noted variation in tolerance in H. azteca is likely associated with inducible traits controlling enzyme pathways. A better understanding of the mechanistic and genomic basis of acclimation is necessary to more accurately predict the ecological and evolutionary consequences of contaminant-driven change in H. azteca.

From Global to Local-New Insights into Features of Pyrethroid Detoxification in Vector Mosquitoes

The threat of mosquito-borne diseases continues to be a problem for public health in subtropical and tropical regions of the world; in response, there has been increased use of adulticidal insecticides, such as pyrethroids, in human habitation areas over the last thirty years. As a result, the prevalence of pyrethroid-resistant genetic markers in natural mosquito populations has increased at an alarming rate. This review details recent advances in the understanding of specific mechanisms associated with pyrethroid resistance, with emphasis on features of insecticide detoxification and the interdependence of multiple cellular pathways. Together, these advances add important context to the understanding of the processes that are selected in resistant mosquitoes. Specifically, before pyrethroids bind to their targets on motoneurons, they must first permeate the outer cuticle and diffuse to inner tissues. Resistant mosquitoes have evolved detoxification mechanisms that rely on cytochrome P450s (CYP), esterases, carboxyesterases, and other oxidation/reduction (redox) components to effectively detoxify pyrethroids to nontoxic breakdown products that are then excreted. Enhanced resistance mechanisms have evolved to include alteration of gene copy number, transcriptional and post-transcriptional regulation of gene expression, as well as changes to cellular signaling mechanisms. Here, we outline the variety of ways in which detoxification has been selected in various mosquito populations, as well as key gene categories involved. Pathways associated with potential new genes of interest are proposed. Consideration of multiple cellular pathways could provide opportunities for development of new insecticides.

A Near-Chromosome Level Genome Assembly of Anopheles stephensi

Malaria remains a major healthcare risk to growing economies like India, and a chromosome-level reference genome of Anopheles stephensi is critical for successful vector management and understanding of vector evolution using comparative genomics. We report chromosome-level assemblies of an Indian strain, STE2, and a Pakistani strain SDA-500 by combining draft genomes of the two strains using a homology-based iterative approach. The resulting assembly IndV3/PakV3 with L50 of 9/12 and N50 6.3/6.9 Mb had scaffolds long enough for building 90% of the euchromatic regions of the three chromosomes, IndV3s/PakV3s, using low-resolution physical markers and enabled the generation of the next version of genome assemblies, IndV4/PakV4, using HiC data. We have validated these assemblies using contact maps against publicly available HiC raw data from two strains including STE2 and another lab strain of An. stephensi from UCI and compare the quality of the assemblies with other assemblies made available as preprints since the submission of the manuscript. We show that the IndV3s and IndV4 assemblies are sensitive in identifying a homozygous 2Rb inversion in the UCI strain and a 2Rb polymorphism in the STE2 strain. Multiple tandem copies of CYP6a14, 4c1, and 4c21 genes, implicated in insecticide resistance, lie within this inversion locus. Comparison of assembled genomes suggests a variation of 1 in 81 positions between the UCI and STE2 lab strains, 1 in 82 between SDA-500 and UCI strain, and 1 in 113 between SDA-500 and STE2 strains of An. stephensi, which are closer than 1 in 68 variations among individuals from two other lab strains sequenced and reported here. Based on the developmental transcriptome and orthology of all the 54 olfactory receptors (ORs) to those of other Anopheles species, we identify an OR with the potential for host recognition in the genus Anopheles. A comparative analysis of An. stephensi genomes with the completed genomes of a few other Anopheles species suggests limited inter-chromosomal gene flow and loss of synteny within chromosomal arms even among the closely related species.

Abiotic and past climatic conditions drive protein abundance variation among natural populations of the caddisfly Crunoecia irrorata

Deducing impacts of environmental change on species and the populations they form in nature is an important goal in contemporary ecology. Achieving this goal is hampered by our limited understanding of the influence of naturally occurring environmental variation on the molecular systems of ecologically relevant species, as the pathways underlying fitness-affecting plastic responses have primarily been studied in model organisms and under controlled laboratory conditions. Here, to test the hypothesis that proteome variation systematically relates to variation in abiotic conditions, we establish such relationships by profiling the proteomes of 24 natural populations of the spring-dwelling caddisfly Crunoecia irrorata. We identified protein networks whose abundances correlated with environmental (abiotic) gradients such as in situ pH, oxygen- and nitrate concentrations but also climatic data such as past thermal minima and temperature seasonality. Our analyses suggest that variations in abiotic conditions induce discrete proteome responses such as the differential abundance of proteins associated with cytoskeletal function, heat-shock proteins and proteins related to post-translational modification. Identifying these drivers of proteome divergence characterizes molecular "noise", and positions it as a background against which molecular signatures of species' adaptive responses to stressful conditions can be identified.

Variation among 532 genomes unveils the origin and evolutionary history of a global insect herbivore

The diamondback moth, Plutella xylostella is a cosmopolitan pest that has evolved resistance to all classes of insecticide, and costs the world economy an estimated US $4-5 billion annually. We analyse patterns of variation among 532 P. xylostella genomes, representing a worldwide sample of 114 populations. We find evidence that suggests South America is the geographical area of origin of this species, challenging earlier hypotheses of an Old-World origin. Our analysis indicates that Plutella xylostella has experienced three major expansions across the world, mainly facilitated by European colonization and global trade. We identify genomic signatures of selection in genes related to metabolic and signaling pathways that could be evidence of environmental adaptation. This evolutionary history of P. xylostella provides insights into transoceanic movements that have enabled it to become a worldwide pest.

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.

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

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

Transcript-Level Analysis of Detoxification Gene Mutation-Mediated Chlorpyrifos Resistance in Laodelphax striatellus (Hemiptera: Delphacidae)

Enhanced detoxification and target mutations that weaken insecticide binding ability are major mechanisms of insecticide resistance. Among these, over-expression or site mutations of carboxylesterase (CarE), cytochrome P450s (CYP450), and glutathione-S-transferase (GST) were the main form responsible for insecticide detoxification; however, transcript-level analysis of the relationship of detoxification gene mutations with chlorpyrifos (an organophosphorus insecticide) resistance is scarce thus far. In this study, multiple sites exhibiting polymorphisms within three detoxification genes were firstly examined via sequencing among different chlorpyrifos-resistant and susceptible individuals of Laodelphax striatellus. For example, the mutation frequencies of A374V in LsCarE16 were 83, 33, and 3%, S277A in LsCarE24 were 88, 28, and 3%, E36K in LsCYP426A1 were 100, 65, and 0% for chlorpyrifos-resistant, resistant decay, and susceptible individuals, respectively. Analysis also found expression levels of GSTd1, GSTt1, GSTs2, CYP4DE1U1, and CYP425B1 are coordinated with chlorpyrifos resistance levels; moreover, we found the deficiencies of 43S and 44A as well as two point mutations of E60D and Q61H at N-terminal region of the OP potential target acetylcholinesterase (AChE) in high resistant but not in low-chlorpyrifos resistant individuals. The results above all demonstrated the dynamic evolutionary process of insecticide resistance and revealed some resistance factors that only played roles at certain resistance level; high insecticide resistance in this example is the result of synergistic impact from multiple resistance factors.

UDP-glycosyltransferase genes and their association and mutations associated with pyrethroid resistance in Anopheles sinensis (Diptera: Culicidae)

BACKGROUND

UDP-glycosyltransferase (UGT) is an important biotransformation superfamily of enzymes. They catalyze the transfer of glycosyl residues from activated nucleotide sugars to acceptor hydrophobic molecules, and function in several physiological processes, including detoxification, olfaction, cuticle formation, pigmentation. The diversity, classification, scaffold location, characteristics, phylogenetics, and evolution of the superfamily of genes at whole genome level, and their association and mutations associated with pyrethroid resistance are still little known.

METHODS

The present study identified UGT genes in Anopheles sinensis genome, classified UGT genes in An. sinensis, Anopheles gambiae, Aedes aegypti and Drosophila melanogaster genomes, and analysed the scaffold location, characteristics, phylogenetics, and evolution of An. sinensis UGT genes using bioinformatics methods. The present study also identified the UGTs associated with pyrethroid resistance using three field pyrethroid-resistant populations with RNA-seq and RT-qPCR, and the mutations associated with pyrethroid resistance with genome re-sequencing in An. sinensis.

RESULTS

There are 30 putative UGTs in An. sinensis genome, which are classified into 12 families (UGT301, UGT302, UGT306, UGT308, UGT309, UGT310, UGT313, UGT314, UGT315, UGT36, UGT49, UGT50) and further into 23 sub-families. The UGT308 is significantly expanded in gene number compared with other families. A total of 119 UGTs from An. sinensis, An. gambiae, Aedes aegypti and Drosophila melanogaster genomes are classified into 19 families, of which seven are specific for three mosquito species and seven are specific for Drosophila melanogaster. The UGT308 and UGT302 are proposed to main families involved in pyrethroid resistance. The AsUGT308D3 is proposed to be the essential UGT gene for the participation in biotransformation in pyrethroid detoxification process, which is possibly regulated by eight SNPs in its 3' flanking region. The UGT302A3 is also associated with pyrethroid resistance, and four amino acid mutations in its coding sequences might enhance its catalytic activity and further result in higher insecticide resistance.

CONCLUSIONS

This study provides the diversity, phylogenetics and evolution of UGT genes, and potential UGT members and mutations involved in pyrethroid resistance in An. sinensis, and lays an important basis for the better understanding and further research on UGT function in defense against insecticide stress.

Protein-protein interaction network analysis of insecticide resistance molecular mechanism in Drosophila melanogaster

The problem of resistance has not been solved fundamentally at present, because the development speed of new insecticides can not keep pace with the development speed of resistance, and the lack of understanding of molecular mechanism of resistance. Here we collected seed genes and their interacting proteins involved in insecticide resistance molecular mechanism in Drosophila melanogaster by literature mining and the String database. We identified a total of 528 proteins and 13514 protein-protein interactions. The protein interaction network was constructed by String and Pajek, and we analyzed the topological properties, such as degree centrality and eigenvector centrality. Proteasome complexes and drug metabolism-cytochrome P450 were an enrichment by Gene Ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. This is the first time to explore the insecticide resistance molecular mechanism of D. melanogaster by the methods and tools of network biology, it can provide the bioinformatic foundation for further understanding the mechanisms of insecticide resistance.

Signatures of Insecticide Selection in the Genome of Drosophila melanogaster

Resistance to insecticides has evolved in multiple insect species, leading to increased application rates and even control failures. Understanding the genetic basis of insecticide resistance is fundamental for mitigating its impact on crop production and disease control. We performed a GWAS approach with the Drosophila Genetic Reference Panel (DGRP) to identify the mutations involved in resistance to two widely used classes of insecticides: organophosphates (OPs, parathion) and pyrethroids (deltamethrin). Most variation in parathion resistance was associated with mutations in the target gene Ace, while most variation in deltamethrin resistance was associated with mutations in Cyp6a23, a gene encoding a detoxification enzyme never previously associated with resistance. A "nested GWAS" further revealed the contribution of other loci: Dscam1 and trpl were implicated in resistance to parathion, but only in lines lacking WolbachiaCyp6a17, the paralogous gene of Cyp6a23, and CG7627, an ATP-binding cassette transporter, were implicated in deltamethrin resistance. We observed signatures of recent selective sweeps at all of these resistance loci and confirmed that the soft sweep at Ace is indeed driven by the identified resistance mutations. Analysis of allele frequencies in additional population samples revealed that most resistance mutations are segregating across the globe, but that frequencies can vary substantially among populations. Altogether, our data reveal that the widely used OP and pyrethroid insecticides imposed a strong selection pressure on natural insect populations. However, it remains unclear why, in Drosophila, resistance evolved due to changes in the target site for OPs, but due to a detoxification enzyme for pyrethroids.

Prokaryotic functional expression and activity comparison of three CYP9A genes from the polyphagous pest Helicoverpa armigera

Cytochrome P450s (CYPs or P450s) have been long recognized as very important enzymes in the metabolism of xenobiotic and endogenous compounds, but only a few CYPs have been functionally characterized in insects. The effort in functional characterization of insect P450s is heavily hindered by technical difficulties in preparing active, individual P450 enzymes directly from the target insect. In this paper, we describe the functional expression of two additional pyrethroid resistance-associated CYP9A genes (CYP9A12 and CYP9A17) from the polyphagous pest Helicoverpa armigera in the facile Escherichia coli. The functionality of E. coli produced CYP9A12, CYP9A14, and CYP9A17 was investigated and activities of these CYP9As were compared against three probe substrates after reconstitution with NADPH-dependent cytochrome P450 reductase. The results showed that active forms of CYP9A12 and CYP9A17 were expressed in E. coli with a content of about 1.0-1.5 nmol mg-1 protein in membrane preparations. In vitro assays showed that CYP9A14 was capable of catalyzing O-dealkylation of methoxyresorufin (MROD), ethoxyresorufin (EROD), and benzyloxyresorufin (BROD), while CYP9A12 and CYP9A17 exhibited only MROD and EROD activities. Kinetic studies demonstrated that CYP9A14 had the greatest k cat/K m value for MROD, and CYP9A17 for EROD, while the lowest k cat/K m values for both MROD and EROD were observed for CYP9A12. The distinct biochemical traits suggest that the three paralogous CYP9As may play different roles in xenobiotic metabolism in this important pest.

The R81T mutation in the nicotinic acetylcholine receptor of Aphis gossypii is associated with neonicotinoid insecticide resistance with differential effects for cyano- and nitro-substituted neonicotinoids

The cotton aphid, Aphis gossypii Glover, is one of the most agriculturally important insect pests. Neonicotinoid insecticides and sulfoxaflor have generally shown excellent control of A. gossypii, but these aphids have recently developed resistance against neonicotinoid insecticides. We previously characterized a field-collected A. gossypii Kushima clone that showed higher resistance to nitro-substituted neonicotinoids, such as imidacloprid, than to cyano-substituted neonicotinoids, such as acetamiprid. This Kushima clone harbors the R81T mutation in the nicotinic acetylcholine receptor (nAChR) β1 subunit; this mutation is the source of neonicotinoid insecticide resistance. In the present study, electrophysiological analyses and molecular modeling were employed to investigate the differential effects of the R81T mutation on cyano- and nitro-substituted neonicotinoids and sulfoxaflor. We isolated full-length coding sequences of A. gossypii nAChR α1, α2, and β1 subunits. When co-expressed in Xenopus laevis oocytes with chicken β2 nAChR, A. gossypii α1 evoked inward currents in a concentration-dependent manner in response to acetylcholine (ACh) and showed sensitivity to neonicotinoid and sulfoxaflor. Additionally, the chicken β2 T77R+E79V (equivalent double mutant of R81T) mutation resulted in a lower effect to cyano-substituted neonicotinoids and sulfoxaflor than to nitro-substituted neonicotinoids. Electrophysiological data and nAChR homology modeling analysis suggested that the Kushima clone exhibited different levels of resistance to cyano- and nitro-substituted neonicotinoid insecticides.

Insights into DDT Resistance from the Drosophila melanogaster Genetic Reference Panel

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

Studies on Aphis gossypii cytochrome P450s CYP6CY22 and CYP6CY13 using an in vitro system

A field-collected Aphis gossypii clone [Kushima resistant (KR) clone] was resistant to neonicotinoid insecticides (23.8- to 394-fold). RNA-seq and next-generation sequence analyses were conducted to identify nine cytochrome P450 (CYP) genes that were significantly upregulated in the KR clone as compared with those in the insecticide-susceptible clone. A. gossypii P450s were transiently and efficiently expressed in S2 cell to show that CYP6CY22 (c21228) and CYP6CY13 (c21368), which were the most upregulated of the nine P450s in the KR clone, did not degrade sulfoxaflor, a new class of insecticides acting on insect nAChRs, but markedly metabolized all of the neonicotinoids tested. Hence, P450s are likely to underpin neonicotinoid resistance in other aphids as well in the future, and the P450 expression protocol established here will prompt studies on P450-medidated insecticide resistance and structural analyses of relevant metabolites.

α-amanitin resistance in Drosophila melanogaster: A genome-wide association approach

We investigated the mechanisms of mushroom toxin resistance in the Drosophila Genetic Reference Panel (DGRP) fly lines, using genome-wide association studies (GWAS). While Drosophila melanogaster avoids mushrooms in nature, some lines are surprisingly resistant to α-amanitin—a toxin found solely in mushrooms. This resistance may represent a pre-adaptation, which might enable this species to invade the mushroom niche in the future. Although our previous microarray study had strongly suggested that pesticide-metabolizing detoxification genes confer α-amanitin resistance in a Taiwanese D. melanogaster line Ama-KTT, none of the traditional detoxification genes were among the top candidate genes resulting from the GWAS in the current study. Instead, we identified Megalin, Tequila, and widerborst as candidate genes underlying the α-amanitin resistance phenotype in the North American DGRP lines, all three of which are connected to the Target of Rapamycin (TOR) pathway. Both widerborst and Tequila are upstream regulators of TOR, and TOR is a key regulator of autophagy and Megalin-mediated endocytosis. We suggest that endocytosis and autophagy of α-amanitin, followed by lysosomal degradation of the toxin, is one of the mechanisms that confer α-amanitin resistance in the DGRP lines.

Multiple cis-acting elements involved in up-regulation of a cytochrome P450 gene conferring resistance to deltamethrin in smal brown planthopper, Laodelphax striatellus (Fallén)

As well as arising from single point mutations in binding sites or detoxifying enzymes, it is likely that insecticide resistance mechanisms are frequently controlled by multiple genetic factors, resulting in resistance being inherited as a quantitative trait. However, empirical evidence for this is still rare. Here we analyse the causes of up-regulation of CYP6FU1, a monoxygenase implicated in resistance to deltamethrin in the rice pest Laodelphax striatellus. The 5'-flanking region of this gene was cloned and sequenced from individuals of a susceptible and a resistant strain. A luminescent reporter assay was used to evaluate different 5'-flanking regions and their fragments for promoter activity. Mutations enhancing promoter activity in various fragments were characterized, singly and in combination, by site mutation recovery. Nucleotide diversity in flanking sequences was greatly reduced in deltamethrin-resistant insects compared to susceptible ones. Phylogenetic sequence analysis found that CYP6FU1 had five different types of 5'-flanking region. All five types were present in a susceptible strain but only a single type showing the highest promoter activity was present in a resistant strain. Four cis-acting elements were identified whose influence on up-regulation was much more pronounced in combination than when present singly. Of these, two were new transcription factor (TF) binding sites produced by mutations, another one was also a new TF binding site alternated from an existing one, and the fourth was a unique transcription start site. These results demonstrate that multiple cis-acting elements are involved in up-regulating CYP6FU1 to generate a resistance phenotype.

Pyrethroid resistance in Aedes aegypti and Aedes albopictus: Important mosquito vectors of human diseases

Aedes aegypti and A. albopictus mosquitoes are vectors of important human disease viruses, including dengue, yellow fever, chikungunya and Zika. Pyrethroid insecticides are widely used to control adult Aedes mosquitoes, especially during disease outbreaks. Herein, we review the status of pyrethroid resistance in A. aegypti and A. albopictus, mechanisms of resistance, fitness costs associated with resistance alleles and provide suggestions for future research. The widespread use of pyrethroids has given rise to many populations with varying levels of resistance worldwide, albeit with substantial geographical variation. In adult A. aegypti and A. albopictus, resistance levels are generally lower in Asia, Africa and the USA, and higher in Latin America, although there are exceptions. Susceptible populations still exist in several areas of the world, particularly in Asia and South America. Resistance to pyrethroids in larvae is also geographically widespread. The two major mechanisms of pyrethroid resistance are increased detoxification due to P450-monooxygenases, and mutations in the voltage sensitive sodium channel (Vssc) gene. Several P450s have been putatively associated with insecticide resistance, but the specific P450s involved are not fully elucidated. Pyrethroid resistance can be due to single mutations or combinations of mutations in Vssc. The presence of multiple Vssc mutations can lead to extremely high levels of resistance. Suggestions for future research needs are presented.

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

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

Allelic Variation of Cytochrome P450s Drives Resistance to Bednet Insecticides in a Major Malaria Vector

Scale up of Long Lasting Insecticide Nets (LLINs) has massively contributed to reduce malaria mortality across Africa. However, resistance to pyrethroid insecticides in malaria vectors threatens its continued effectiveness. Deciphering the detailed molecular basis of such resistance and designing diagnostic tools is critical to implement suitable resistance management strategies. Here, we demonstrated that allelic variation in two cytochrome P450 genes is the most important driver of pyrethroid resistance in the major African malaria vector Anopheles funestus and detected key mutations controlling this resistance. An Africa-wide polymorphism analysis of the duplicated genes CYP6P9a and CYP6P9b revealed that both genes are directionally selected with alleles segregating according to resistance phenotypes. Modelling and docking simulations predicted that resistant alleles were better metabolizers of pyrethroids than susceptible alleles. Metabolism assays performed with recombinant enzymes of various alleles confirmed that alleles from resistant mosquitoes had significantly higher activities toward pyrethroids. Additionally, transgenic expression in Drosophila showed that flies expressing resistant alleles of both genes were significantly more resistant to pyrethroids compared with those expressing the susceptible alleles, indicating that allelic variation is the key resistance mechanism. Furthermore, site-directed mutagenesis and functional analyses demonstrated that three amino acid changes (Val¹⁰⁹Ile, Asp³³⁵Glu and Asn³⁸⁴Ser) from the resistant allele of CYP6P9b were key pyrethroid resistance mutations inducing high metabolic efficiency. The detection of these first DNA markers of metabolic resistance to pyrethroids allows the design of DNA-based diagnostic tools to detect and track resistance associated with bednets scale up, which will improve the design of evidence-based resistance management strategies.

Scale up of Long Lasting Insecticide Nets has massively reduced malaria mortality across Africa. However, resistance to pyrethroid insecticides in malaria vectors threatens its continued effectiveness. Here, we established that allelic variation in two CYP450s is the most important driver of pyrethroid resistance in the major African vector Anopheles funestus and detected key mutations controlling this resistance. The duplicated P450s CYP6P9a and CYP6P9b are directionally selected across Africa with alleles segregating according to resistance phenotypes. Alleles from resistant mosquitoes present significantly higher metabolic activities toward pyrethroids compared with alleles from susceptible mosquitoes. Furthermore, transgenic flies over-expressing resistant alleles of both genes were significantly more resistant to pyrethroids. Three amino acid changes from the resistant CYP6P9b allele are the key pyrethroid resistance mutations which induce high metabolic efficiency. The detection of these first DNA markers of metabolic resistance to pyrethroids allows the design of diagnostic tools to detect and track resistance.

Identifying genomic changes associated with insecticide resistance in the dengue mosquito Aedes aegypti by deep targeted sequencing

The capacity of mosquitoes to resist insecticides threatens the control of diseases such as dengue and malaria. Until alternative control tools are implemented, characterizing resistance mechanisms is crucial for managing resistance in natural populations. Insecticide biodegradation by detoxification enzymes is a common resistance mechanism; however, the genomic changes underlying this mechanism have rarely been identified, precluding individual resistance genotyping. In particular, the role of copy number variations (CNVs) and polymorphisms of detoxification enzymes have never been investigated at the genome level, although they can represent robust markers of metabolic resistance. In this context, we combined target enrichment with high-throughput sequencing for conducting the first comprehensive screening of gene amplifications and polymorphisms associated with insecticide resistance in mosquitoes. More than 760 candidate genes were captured and deep sequenced in several populations of the dengue mosquito Ae. aegypti displaying distinct genetic backgrounds and contrasted resistance levels to the insecticide deltamethrin. CNV analysis identified 41 gene amplifications associated with resistance, most affecting cytochrome P450s overtranscribed in resistant populations. Polymorphism analysis detected more than 30,000 variants and strong selection footprints in specific genomic regions. Combining Bayesian and allele frequency filtering approaches identified 55 nonsynonymous variants strongly associated with resistance. Both CNVs and polymorphisms were conserved within regions but differed across continents, confirming that genomic changes underlying metabolic resistance to insecticides are not universal. By identifying novel DNA markers of insecticide resistance, this study opens the way for tracking down metabolic changes developed by mosquitoes to resist insecticides within and among populations.

Long-Term Resistance of Drosophila melanogaster to the Mushroom Toxin Alpha-Amanitin

Insect resistance to toxins exerts not only a great impact on our economy, but also on the ecology of many species. Resistance to one toxin is often associated with cross-resistance to other, sometimes unrelated, chemicals. In this study, we investigated mushroom toxin resistance in the fruit fly Drosophila melanogaster (Meigen). This fruit fly species does not feed on mushrooms in nature and may thus have evolved cross-resistance to α-amanitin, the principal toxin of deadly poisonous mushrooms, due to previous pesticide exposure. The three Asian D. melanogaster stocks used in this study, Ama-KTT, Ama-MI, and Ama-KLM, acquired α-amanitin resistance at least five decades ago in their natural habitats in Taiwan, India, and Malaysia, respectively. Here we show that all three stocks have not lost the resistance phenotype despite the absence of selective pressure over the past half century. In response to α-amanitin in the larval food, several signs of developmental retardation become apparent in a concentration-dependent manner: higher pre-adult mortality, prolonged larva-to-adult developmental time, decreased adult body size, and reduced adult longevity. In contrast, female fecundity nearly doubles in response to higher α-amanitin concentrations. Our results suggest that α-amanitin resistance has no fitness cost, which could explain why the resistance has persisted in all three stocks over the past five decades. If pesticides caused α-amanitin resistance in D. melanogaster, their use may go far beyond their intended effects and have long-lasting effects on ecosystems.

A novel cytochrome P450 CYP6AB14 gene in Spodoptera litura (Lepidoptera: Noctuidae) and its potential role in plant allelochemical detoxification

Cytochrome P450 monooxygenases (P450s) play a prominent role in the adaptation of insects to host plant chemical defenses. To investigate the potential role of P450s in adaptation of the lepidopteran pest Spodoptera litura to host plant allelochemicals, an expressed sequence data set derived from 6th instar midgut tissues was first mined. One sequence identified from the S. litura 6th instar midgut EST database was determined by phylogenetic analysis to belong to the CYP6AB P450 subfamily, and named CYP6AB14. Dietary supplementation of S. litura larvae with either xanthotoxin (XAN), coumarin (COU) and flavone (FLA) led to elevated CYP6AB14 transcript levels in both midgut and fat body tissues. Injection of CYP6AB14-derived double-stranded RNA (dsRNA) into S. litura individuals significantly reduced CYP6AB14 transcript levels, and resulted in increased developmental abnormalities and higher mortality rates among XAN, COU and FLA-fed larvae. Our results strongly suggest a key role for CYP6AB14 in plant allelochemical detoxification in S. litura.

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

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

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

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

The cytochrome P450 family in the parasitic nematode Haemonchus contortus

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.

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.

Mechanisms of pyrethroid resistance in the dengue mosquito vector, Aedes aegypti: target site insensitivity, penetration, and metabolism

Aedes aegypti is the major vector of yellow and dengue fevers. After 10 generations of adult selection, an A. aegypti strain (SP) developed 1650-fold resistance to permethrin, which is one of the most widely used pyrethroid insecticides for mosquito control. SP larvae also developed 8790-fold resistance following selection of the adults. Prior to the selections, the frequencies of V1016G and F1534C mutations in domains II and III, respectively, of voltage-sensitive sodium channel (Vssc, the target site of pyrethroid insecticide) were 0.44 and 0.56, respectively. In contrast, only G1016 alleles were present after two permethrin selections, indicating that G1016 can more contribute to the insensitivity of Vssc than C1534. In vivo metabolism studies showed that the SP strain excreted permethrin metabolites more rapidly than a susceptible SMK strain. Pretreatment with piperonyl butoxide caused strong inhibition of excretion of permethrin metabolites, suggesting that cytochrome P450 monooxygenases (P450s) play an important role in resistance development. In vitro metabolism studies also indicated an association of P450s with resistance. Microarray analysis showed that multiple P450 genes were over expressed during the larval and adult stages in the SP strain. Following quantitative real time PCR, we focused on two P450 isoforms, CYP9M6 and CYP6BB2. Transcription levels of these P450s were well correlated with the rate of permethrin excretion and they were certainly capable of detoxifying permethrin to 4'-HO-permethrin. Over expression of CYP9M6 was partially due to gene amplification. There was no significant difference in the rate of permethrin reduction from cuticle between SP and SMK strains.

The molecular evolution of cytochrome P450 genes within and between drosophila species

We map 114 gene gains and 74 gene losses in the P450 gene family across the phylogeny of 12 Drosophila species by examining the congruence of gene trees and species trees. Although the number of P450 genes varies from 74 to 94 in the species examined, we infer that there were at least 77 P450 genes in the ancestral Drosophila genome. One of the most striking observations in the data set is the elevated loss of P450 genes in the Drosophila sechellia lineage. The gain and loss events are not evenly distributed among the P450 genes-with 30 genes showing no gene gains or losses whereas others show as many as 20 copy number changes among the species examined. The P450 gene clades showing the fewest number of gene gain and loss events tend to be those evolving with the most purifying selection acting on the protein sequences, although there are exceptions, such as the rapid rate of amino acid replacement observed in the single copy phantom (Cyp306a1) gene. Within D. melanogaster, we observe gene copy number polymorphism in ten P450 genes including multiple cases of interparalog chimeras. Nonallelic homologous recombination (NAHR) has been associated with deleterious mutations in humans, but here we provide a second possible example of an NAHR event in insect P450s being adaptive. Specifically, we find that a polymorphic Cyp12a4/Cyp12a5 chimera correlates with resistance to an insecticide. Although we observe such interparalog exchange in our within-species data sets, we have little evidence of it between species, raising the possibility that such events may occur more frequently than appreciated but are masked by subsequent sequence change.

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

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

A P450 gene associated with robust resistance to DDT in ciliated protozoan, Tetrahymena thermophila by efficient degradation

Analysis of metabolic mechanisms of dichlorodiphenyltrichloroethane (DDT) accumulation and degradation in microorganisms, which could be used to reduce its hazard to higher organisms at the higher in the food chain, have not been investigated. Robust resistance to DDT (grows well in 256 mg/L DDT) and a surprising ability to degrade DDT (more than 70% DDT within 4h) were found in the ciliated protozoan Tetrahymena thermophila. A P450 gene (CYP5013C2) was found to respond specifically to DDT treatment. In the presence of 256 mg/L DDT, cells with overexpressing CYP5013C2 (p450-OE) grew faster and degraded DDT more efficiently than wild-type (WT) cells, while cells with CYP5013C2 partially knocked down (p450-KD) grew slower and exhibited reduced ability to degrade DDT compared to WT cells. Both dichlorodiphenyldichloroethylene (DDE) and dichlorodiphenyldichloroethane (DDD) were detected in cells after exposure to DDT, and the concentration of DDD in the p450-OE strain gradually decreased from 0.5 to 4h. Thus, we argue that this P450 gene (CYP5013C2), by efficiently degrading DDT to DDD, is associated with robust resistance to DDT in Tetrahymena, and that a strain overexpressing this gene has the potential to serve as bioreactor that degrades environmental DDT.

Mining gene expression data for pollutants (dioxin, toluene, formaldehyde) and low dose of gamma-irradiation

General and specific effects of molecular genetic responses to adverse environmental factors are not well understood. This study examines genome-wide gene expression profiles of Drosophila melanogaster in response to ionizing radiation, formaldehyde, toluene, and 2,3,7,8-tetrachlorodibenzo-p-dioxin. We performed RNA-seq analysis on 25,415 transcripts to measure the change in gene expression in males and females separately. An analysis of the genes unique to each treatment yielded a list of genes as a gene expression signature. In the case of radiation exposure, both sexes exhibited a reproducible increase in their expression of the transcription factors sugarbabe and tramtrack. The influence of dioxin up-regulated metabolic genes, such as anachronism, CG16727, and several genes with unknown function. Toluene activated a gene involved in the response to the toxins, Cyp12d1-p; the transcription factor Fer3's gene; the metabolic genes CG2065, CG30427, and CG34447; and the genes Spn28Da and Spn3, which are responsible for reproduction and immunity. All significantly differentially expressed genes, including those shared among the stressors, can be divided into gene groups using Gene Ontology Biological Process identifiers. These gene groups are related to defense response, biological regulation, the cell cycle, metabolic process, and circadian rhythms. KEGG molecular pathway analysis revealed alteration of the Notch signaling pathway, TGF-beta signaling pathway, proteasome, basal transcription factors, nucleotide excision repair, Jak-STAT signaling pathway, circadian rhythm, Hippo signaling pathway, mTOR signaling pathway, ribosome, mismatch repair, RNA polymerase, mRNA surveillance pathway, Hedgehog signaling pathway, and DNA replication genes. Females and, to a lesser extent, males actively metabolize xenobiotics by the action of cytochrome P450 when under the influence of dioxin and toluene. Finally, in this work we obtained gene expression signatures pollutants (dioxin, toluene), low dose of gamma-irradiation and common molecular pathways for different kind of stressors.

Molecular cloning and xenobiotic induction of seven novel cytochrome P450 monooxygenases in Aedes albopictus

Cytochrome P450 monooxygenase (P450) is a superfamily of enzymes that is important in metabolism of endogenous and exogenous compounds. In insects, these enzymes confer resistance to insecticides through its metabolic activities. Members of P450 from family 6 in insects are known to play a role in such function. In this study, we have isolated seven novel family 6 P450 from Aedes albopictus (Skuse) (Diptera: Culicidae), a vector of dengue and chikungunya fever. Induction profile of these seven genes was studied using several insecticides and xenobiotics. It was found that deltamethrin and permethrin did not induce expression of any genes. Another insecticide, temephos, inhibited expression of CYP6P15 for fivefold and twofold for CYP6N29, CYP6Y7, and CYP6Z18. In addition, copper II sulfate induced expression of CYP6M17 and CYP6N28 for up to sixfold. Benzothiazole (BZT), a tire leachate induced the expression of CYP6M17 by fourfold, CYP6N28 by sevenfold, but inhibited the expression of CYP6P15 for threefold and CYP6Y7 for twofold. Meanwhile, piperonyl butoxide (PBO) induced the expression CYP6N28 (twofold), while it inhibited the expression of CYP6P15 (fivefold) and CYP6Y7 (twofold). Remarkably, all seven genes were induced two- to eightfold by acetone in larval stage, but not adult stage. Expression of CYP6N28 was twofold higher, while expression of CYP6P15 was 15-fold lower in adult than larva. The other five P450s were not differentially expressed between the larvae and adult. This finding showed that acetone can be a good inducer of P450 in Ae. albopictus. On the other hand, temephos can act as good suppressor of P450, which may affect its own bioefficacy because it needs to be bioactivated by P450. To the best of our knowledge, this is the first report on acetone-inducible P450 in insects. Further study is needed to characterize the mechanisms involved in acetone induction in P450.

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.