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

Xenobiotic responses of Drosophila melanogaster to insecticides with different modes of action and entry

Depending on their chemical structure, insecticides enter the insect body either through the cuticle or by ingestion (mode of entry [MoE]), and, naturally, harm or even kill insects through different mechanisms (modes of action). In parallel, they trigger a systemic detoxification response, especially by activation of detoxification gene expression. We monitored the acute genetic alterations of known xenobiotic response target genes against five different insecticides with two most common MoEs (contact toxicity and stomach toxicity), found that: 1. only a few genes were detected responding to acute exposure to insecticides (LD₉₀ ); 2. The expression of cyp12d1 was upregulated in all experiments, except for dichlorodiphenyltrichloroethane exposure, suggesting that cyp12d1 is a general first response gene of the xenobiotic response; 3. The contact and stomach entries did not show any notable difference, both MoEs induced the response of JNK signaling pathway, possibly serving as the driver of the response of cyp12d1 and a few other genes. In conclusion, the changes in gene expression levels were relatively modest and no significant differences were found between the two MoEs, so the insecticide entry route does not seem to have an impact on the detoxification response. However, the two MoEs of the same insecticide showed different efficiencies in our test. Thus, the study of these two MoEs will help to develop more efficient release and management methods for the use of such insecticides.

The Power of Drosophila melanogaster for Modeling Neonicotinoid Effects on Pollinators and Identifying Novel Mechanisms

Neonicotinoids are the most widely used insecticides in the world and are implicated in the widespread population declines of insects including pollinators. Neonicotinoids target nicotinic acetylcholine receptors which are expressed throughout the insect central nervous system, causing a wide range of sub-lethal effects on non-target insects. Here, we review the potential of the fruit fly Drosophila melanogaster to model the sub-lethal effects of neonicotinoids on pollinators, by utilizing its well-established assays that allow rapid identification and mechanistic characterization of these effects. We compare studies on the effects of neonicotinoids on lethality, reproduction, locomotion, immunity, learning, circadian rhythms and sleep in D. melanogaster and a range of pollinators. We also highlight how the genetic tools available in D. melanogaster, such as GAL4/UAS targeted transgene expression system combined with RNAi lines to any gene in the genome including the different nicotinic acetylcholine receptor subunit genes, are set to elucidate the mechanisms that underlie the sub-lethal effects of these common pesticides. We argue that studying pollinators and D. melanogaster in tandem allows rapid elucidation of mechanisms of action, which translate well from D. melanogaster to pollinators. We focus on the recent identification of novel and important sublethal effects of neonicotinoids on circadian rhythms and sleep. The comparison of effects between D. melanogaster and pollinators and the use of genetic tools to identify mechanisms make a powerful partnership for the future discovery and testing of more specific insecticides.

Linkage mapping an indoxacarb resistance locus in Helicoverpa armigera (Lepidoptera: Noctuidae) by genotype-by-sequencing

BACKGROUND

A major challenge to sustainable agricultural pest control is the rapid evolution of insecticide resistance. This is caused by mechanisms that reduce insecticide efficacy. Understanding the genetic mechanisms of resistance is essential for DNA-based monitoring of resistance in field populations. One such insecticide is indoxacarb, an important selective control option for Helicoverpa armigera in a range of crops including grain, horticulture and cotton. Recently, a strain of H. armigera (GY7-39) resistant to indoxacarb (198-fold) was isolated from field-collected moth.

RESULTS

To identify the indoxacarb resistance locus, GY7-39 was backcrossed for six generations to susceptible strain New GR. In each generation, only resistant males were used to cross back to New GR. Genotype-by-sequencing was carried out on 95 H. armigera samples. In total, 13 203 tags with 8697 unique locations on the H. armigera genome were obtained. The indoxacarb resistance locus in strain GY7-39 was mapped to a 2.6 Mbp region on chromosome 16. In this region, two closely linked loci (IndoR1 and IndoR2) were found to be associated with indoxacarb resistant GY7-39.

CONCLUSIONS

We mapped indoxacarb resistance in GY7-39 to two closely linked loci IndoR1 and IndoR2 in a narrowed 2.6 Mbp region of H. armigera chromosome 16. The results provide essential background data for future genetic investigations including fine mapping of the indoxacarb resistance gene and the eventual development of an effective DNA-based diagnostic to support resistance management. © 2019 Society of Chemical Industry.

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.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Changes in Neuronal Signaling and Cell Stress Response Pathways are Associated with a Multigenic Response of Drosophila melanogaster to DDT Selection

The adaptation of insect populations to insecticidal control is a continual threat to human health and sustainable agricultural practices, but many complex genomic mechanisms involved in this adaption remain poorly understood. This study applied a systems approach to investigate the interconnections between structural and functional variance in response to dichlorodiphenyltrichloroethane (DDT) within the Drosophila melanogaster strain 91-R. Directional selection in 6 selective sweeps coincided with constitutive gene expression differences in DDT resistant flies, including the most highly upregulated transcript, Unc-115 b, which plays a central role in axon guidance, and the most highly downregulated transcript, the angiopoietin-like CG31832, which is involved in directing vascular branching and dendrite outgrowth but likely may be under trans-regulatory control. Direct functions and protein–protein interactions mediated by differentially expressed transcripts control changes in cell migration, signal transduction, and gene regulatory cascades that impact the nervous system. Although changes to cellular stress response pathways involve 8 different cytochrome P450s, stress response, and apoptosis is controlled by a multifacetted regulatory mechanism. These data demonstrate that DDT selection in 91-R may have resulted in genome-wide adaptations that impacts genetic and signal transduction pathways that converge to modify stress response, cell survival, and neurological functions. This study implicates the involvement of a multigenic mechanism in the adaptation to a chemical insecticide, which impact interconnected regulatory cascades. We propose that DDT selection within 91-R might act systemically, wherein pathway interactions function to reinforce the epistatic effects of individual adaptive changes on an additive or nonadditive basis.

Identification of Differentially Expressed miRNAs in Colorado Potato Beetles (Leptinotarsa decemlineata (Say)) Exposed to Imidacloprid

The Colorado potato beetle (Leptinotarsa decemlineata (Say)) is a significant pest of potato plants that has been controlled for more than two decades by neonicotinoid imidacloprid. L. decemlineata can develop resistance to this agent even though the molecular mechanisms underlying this resistance are not well characterized. MicroRNAs (miRNAs) are short ribonucleic acids that have been linked to response to various insecticides in several insect models. Unfortunately, the information is lacking regarding differentially expressed miRNAs following imidacloprid treatment in L. decemlineata. In this study, next-generation sequencing and quantitative real-time polymerase chain reaction (qRT-PCR) were used to identify modulated miRNAs in imidacloprid-treated versus untreated L. decemlineata. This approach identified 33 differentially expressed miRNAs between the two experimental conditions. Of interest, miR-282 and miR-989, miRNAs previously shown to be modulated by imidacloprid in other insects, and miR-100, a miRNA associated with regulation of cytochrome P450 expression, were significantly modulated in imidacloprid-treated beetles. Overall, this work presents the first report of a miRNA signature associated with imidacloprid exposure in L. decemlineata using a high-throughput approach. It also reveals interesting miRNA candidates that potentially underly imidacloprid response in this insect pest.

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

Discovery of Organophosphate Resistance-Related Genes Associated With Well-known Resistance Mechanisms of Plutella xylostella (L.) (Lepidoptera: Plutellidae) by RNA-Seq

Pesticide resistance poses many challenges for pest control, particularly for destructive pests such as diamondback moths ( Plutella xylostella ). Organophosphates have been used in the field since the 1950s, leading to selection for resistance-related gene variants and the development of resistance to new insecticides in the diamondback moth. Identifying actual and potential genes involved in resistance could offer solutions for control. This study established resistant diamondback moth strains from two different collections using mevinphos. Two sets of transcriptome sequencing (RNA-Seq) data were generated for pairs of mevinphos-resistant versus susceptible (wild-type) strains. One susceptible strain containing 14 giga base pairs was assembled into a reference-based assembly using published scaffold sequences as reference. Differential expression data between resistant and susceptible strains revealed 944 transcripts (803 with annotations) showing upregulation and 427 transcripts (150 with annotations) showing downregulation. Around 6.8% of the differential expression transcripts (65) could be categorized as associated with well-known resistance mechanisms such as penetration, detoxification, and behavior response; of these 65 transcripts, 38 showed upregulation, and 12 relating to penetration were upregulated when the transcripts of 19 cytochrome P450s, 2 zeta-class glutathione S-transferases, and 4 ATP-binding cassette transporters showed upregulation. In addition, 11 groups of transcripts related to olfactory perception appeared to be downregulated in trade-off situations. Quantitative polymerase chain reaction expression results were consistent with RNA-Seq data. Possible roles of these differentially expressed genes in resistance mechanisms are discussed in this study.

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

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

Spirotetramat resistance adaption analysis of Aphis gossypii Glover by transcriptomic survey

A resistant strain of the cotton aphid (SR) developed 441.26-fold and 11.97-fold resistance to spirotetramat for adult aphids and nymphs, respectively, compared with the susceptible (SS) strain. Solexa sequencing technology was employed to identify differentially expressed genes (DEGs) in the spirotetramat-resistant cotton aphid. Respective totals of 22,430,522 and 21,317,732 clean reads were obtained from SR and SS cDNA libraries and assembled into 35,222 non-redundant (Nr) consensus sequences. A total of 14,913, 9,220, 7,922, 4,314 and 4,686 sequences were annotated using Nr, Swiss-Prot, Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) and Clusters of Orthologous Groups (COG), respectively. Compared with the SS strain, the SR strain had 1287 significantly changed unigenes, of which 130 genes were up-regulated and 1157 genes were down-regulated (P ≤ 0.001). Among these genes, 440 unigenes were annotated, consisting of 114 up-regulated and 326 down-regulated genes. The expression levels of heat shock protein 70 (Hsp70) and UDP-glucuronosyltransferase were significantly up-regulated in the SR strain compared to the SS strain. The genes encoding cuticle proteins, salivary glue protein, fibroin heavy chain, energy ATP synthase, and cytochrome c oxidase were dramatically decreased. Among the DEGs, cytochrome P450 6A2 (c20965.graph_c0) was the only P450 gene up-regulated in the SR strain. The expression levels of 10 DEGs were confirmed by real-time qPCR, and the trends in gene expression observed by qPCR matched those of the Solexa expression profiles. The acetyl-CoA carboxylase (ACC) genes in the SR and SS libraries both contain four single nucleotide polymorphisms (SNPs), with three common SNPs: 1227 (C/T), 1811 (A/T: F/Y) and 3759 (C/T); however, 7540 (A/T) and 108 (G/A) occurred solely in the SS and SR strains, respectively.

Effect of Low Doses (5-40 cGy) of Gamma-irradiation on Lifespan and Stress-related Genes Expression Profile in Drosophila melanogaster

Studying of the effects of low doses of γ-irradiation is a crucial issue in different areas of interest, from environmental safety and industrial monitoring to aerospace and medicine. The goal of this work is to identify changes of lifespan and expression stress-sensitive genes in Drosophila melanogaster, exposed to low doses of γ-irradiation (5 – 40 cGy) on the imaginal stage of development. Although some changes in life extensity in males were identified (the effect of hormesis after the exposure to 5, 10 and 40 cGy) as well as in females (the effect of hormesis after the exposure to 5 and 40 cGy), they were not caused by the organism “physiological” changes. This means that the observed changes in life expectancy are not related to the changes of organism physiological functions after the exposure to low doses of ionizing radiation. The identified changes in gene expression are not dose-dependent, there is not any proportionality between dose and its impact on expression. These results reflect nonlinear effects of low dose radiation and sex-specific radio-resistance of the postmitotic cell state of Drosophila melanogaster imago.

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.

Effects of neonicotinoids and fipronil on non-target invertebrates

We assessed the state of knowledge regarding the effects of large-scale pollution with neonicotinoid insecticides and fipronil on non-target invertebrate species of terrestrial, freshwater and marine environments. A large section of the assessment is dedicated to the state of knowledge on sublethal effects on honeybees (Apis mellifera) because this important pollinator is the most studied non-target invertebrate species. Lepidoptera (butterflies and moths), Lumbricidae (earthworms), Apoidae sensu lato (bumblebees, solitary bees) and the section "other invertebrates" review available studies on the other terrestrial species. The sections on freshwater and marine species are rather short as little is known so far about the impact of neonicotinoid insecticides and fipronil on the diverse invertebrate fauna of these widely exposed habitats. For terrestrial and aquatic invertebrate species, the known effects of neonicotinoid pesticides and fipronil are described ranging from organismal toxicology and behavioural effects to population-level effects. For earthworms, freshwater and marine species, the relation of findings to regulatory risk assessment is described. Neonicotinoid insecticides exhibit very high toxicity to a wide range of invertebrates, particularly insects, and field-realistic exposure is likely to result in both lethal and a broad range of important sublethal impacts. There is a major knowledge gap regarding impacts on the grand majority of invertebrates, many of which perform essential roles enabling healthy ecosystem functioning. The data on the few non-target species on which field tests have been performed are limited by major flaws in the outdated test protocols. Despite large knowledge gaps and uncertainties, enough knowledge exists to conclude that existing levels of pollution with neonicotinoids and fipronil resulting from presently authorized uses frequently exceed the lowest observed adverse effect concentrations and are thus likely to have large-scale and wide ranging negative biological and ecological impacts on a wide range of non-target invertebrates in terrestrial, aquatic, marine and benthic habitats.

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 whole transcriptomal linkage analysis of gene co-regulation in insecticide resistant house flies, Musca domestica

Background

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

Results

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

Conclusion

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

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

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