Cis- and trans-acting variants contribute to survivorship in a naïve Drosophila melanogaster population exposed to ryanoid insecticides

A masked gene concealed hand in glove in the forkhead protein crocodile regulates the predominant detoxification CYP6DA1 in Aphis gossypii Glover

Cytochrome P450-mediated metabolism is an important mechanism of insecticide resistance, most studies show upregulated transcript levels of P450s in resistant insect strains. Our previous studies illustrated that some upregulated P450s were associated with cyantraniliprole resistance, and it is more comprehensive to use the tissue specificity of transcriptomes to compare resistant (CyR) and susceptible (SS) strains. In this study, the expression profiles of P450s in a CyR strain compared with a SS strain in remaining carcass or midgut were investigated by RNA sequencing, and candidate genes were selected for functional study. Drosophila melanogaster bioassays suggested that ectopic overexpression of CYP4CK1, CYP6CY5, CYP6CY9, CYP6CY19, CYP6CZ1 and CYP6DA1 in flies was sufficient to confer cyantraniliprole resistance, among which CYP6DA1 was the predominant contributor to resistance (12.24-fold). RNAi suppression of CYP4CK1, CYP6CY5, CYP6CY9 and CYP6DA1 significantly increased CyR aphid sensitivity to cyantraniliprole. The CYP6DA1 promoter had two predicted binding sites for crocodile (CROC), an intron-free ORF with bidirectional transcription yielding CROC (+) and CROC (-). Y1H, RNAi and EMSA found that CROC (-) was a transcription factor directly regulating CYP6DA1 expression. In conclusion, P450 genes contribute to cyantraniliprole resistance, and the transcription factor CROC (-) regulates the expression of CYP6DA1 in A. gossypii.

The genomic basis of copper tolerance in Drosophila is shaped by a complex interplay of regulatory and environmental factors

BACKGROUND

Escalation in industrialization and anthropogenic activity have resulted in an increase of pollutants released into the environment. Of these pollutants, heavy metals such as copper are particularly concerning due to their bio-accumulative nature. Due to its highly heterogeneous distribution and its dual nature as an essential micronutrient and toxic element, the genetic basis of copper tolerance is likely shaped by a complex interplay of genetic and environmental factors.

RESULTS

In this study, we utilized the natural variation present in multiple populations of Drosophila melanogaster collected across Europe to screen for variation in copper tolerance. We found that latitude and the degree of urbanization at the collection sites, rather than any other combination of environmental factors, were linked to copper tolerance. While previously identified copper-related genes were not differentially expressed in tolerant vs. sensitive strains, genes involved in metabolism, reproduction, and protease induction contributed to the differential stress response. Additionally, the greatest transcriptomic and physiological responses to copper toxicity were seen in the midgut, where we found that preservation of gut acidity is strongly linked to greater tolerance. Finally, we identified transposable element insertions likely to play a role in copper stress response.

CONCLUSIONS

Overall, by combining genome-wide approaches with environmental association analysis, and functional analysis of candidate genes, our study provides a unique perspective on the genetic and environmental factors that shape copper tolerance in natural D. melanogaster populations and identifies new genes, transposable elements, and physiological traits involved in this complex phenotype.

Discovery of malathion resistance QTL in Drosophila melanogaster using a bulked phenotyping approach

Drosophila melanogaster has proved an effective system with which to understand the evolutionary genetics and molecular mechanisms of insecticide resistance. Insecticide use has left signatures of selection in the fly genome, and both functional and quantitative genetic studies in the system have identified genes and variants associated with resistance. Here, we use D. melanogaster and leverage a bulk phenotyping and pooled sequencing "extreme quantitative trait loci" approach to genetically dissect variation in resistance to malathion, an organophosphate insecticide. We resolve 2 quantitative trait loci, one of which implicates allelic variation at the cytochrome P450 gene Cyp6g1, a strong candidate based on previous work. The second shows no overlap with hits from a previous genome-wide association study for malathion resistance, recapitulating other studies showing that different strategies for complex trait dissection in flies can yield apparently different architectures. Notably, we see no genetic signal at the Ace gene. Ace encodes the target of organophosphate insecticide inhibition, and genome-wide association studies have identified strong Ace-linked associations with resistance in flies. The absence of quantitative trait locus implicating Ace here is most likely because our mapping population does not segregate for several of the known functional polymorphisms impacting resistance at Ace, perhaps because our population is derived from flies collected prior to the widespread use of organophosphate insecticides. Our fundamental approach can be an efficient, powerful strategy to dissect genetic variation in resistance traits. Nonetheless, studies seeking to interrogate contemporary insecticide resistance variation may benefit from deriving mapping populations from more recently collected strains.

Functional Inquiry into ATP-Binding Cassette Transporter Genes Contributing to Spirotetramat Resistance in Aphis gossypii Glover

ATP-binding cassette (ABC) transporters regulate the efflux of a broad spectrum of substrates to extracellular transporting, which play an important role in the detoxification process in arthropods. Here, we described a comprehensive approach to explore the involvement of ABC transporters in spirotetramat resistance in cotton aphids. In this study, synergism bioassays showed 17.05% and 35.42% increases in the toxicity to spirotetramat with the ABC inhibitor verapamil in adult and 3rd instar nymph aphids of the SR strain, respectively. In a competitive assay based on the microinjection of a fluorescent ABC transporter substrate, verapamil (a general ABC inhibitor) and spirotetramat significantly inhibited the elimination of Texas Red. Based on transcriptome data of midguts of spirotetramat-susceptible (SS) and -resistant (SR) strains, the expression levels of ABCB4, ABCB5, ABCF2, MRP11, and MRP12 were significantly upregulated in the SR strain midgut compared to that of the SS strain. Gene functional analysis based on ectopic expression and RNA interference (RNAi) proved that ABCB4, ABCB5, ABCF2, MRP11, and MRP12 were involved in the tolerance of cotton aphids to spirotetramat. Moreover, the upregulated ABCF2, ABCB4, and ABCB5 in the midgut of the SR strain contributed more to the resistance of spirotetramat in in vitro functional analysis. In summary, these results demonstrate that candidate ABC transporter genes in the midgut tissue were involved in spirotetramat resistance, which will help reveal the relationship between ABC transporters and the development of spirotetramat resistance in field populations.

Chemosensory proteins confer adaptation to the ryanoid anthranilic diamide insecticide cyantraniliprole in Aphis gossypii glover

Chemosensory proteins (CSPs) are a class of small transporter proteins expressed only in arthropods with various functions beyond chemoreception. Previous studies have been reported that CSPs are involved in the insecticide resistance. In this study, we found that AgoCSP1, AgoCSP4, and AgoCSP5 were constitutively overexpressed in an insecticide-resistant strain of Aphis gossypii and showed higher expression in broad body tissue (including fat bodies) than in the midgut but without tissue specificity. However, the function of these three upregulated AgoCSPs remains unknown. Here, we investigated the function of AgoCSPs in resistance to the diamide insecticide cyantraniliprole. Suppression of AgoCSP1, AgoCSP4 and AgoCSP5 transcription by RNAi significantly increased the sensitivity of resistant aphids to cyantraniliprole. Molecular docking and competitive binding assays indicated that these AgoCSPs bind moderate with cyantraniliprole. Transgenic Drosophila melanogaster expressing these AgoCSPs in the broad body or midgut showed higher tolerance to cyantraniliprole than control flies with the same genetic background; AgoCSP4 was more effective in broad body tissue, and AgoCSP1 and AgoCSP5 were more effective in the midgut, indicating that broad body and midgut tissues may be involved in the insecticide resistance mediated by the AgoCSPs examined. The present results strongly indicate that AgoCSPs participate in xenobiotic detoxification by sequestering and masking toxic insecticide molecules, providing insights into new factors involved in resistance development in A. gossypii.

Functional investigation of lncRNAs and target cytochrome P450 genes related to spirotetramat resistance in Aphis gossypii Glover

BACKGROUND

Spirotetramat is a tetramic acid derivative insecticide with novel modes of action for controlling Aphis gossypii Glover in the field. Previous studies have shown that long noncoding RNAs (lncRNAs) and cytochrome P450 monooxygenases (P450s) are involved in the detoxification process. However, the functions of lncRNAs in regulating P450 gene expression in spirotetramat resistance in A. gossypii are unknown.

RESULTS

In this study, we found CYP4CJ1, CYP6CY7 and CYP6CY21 expression levels to be significantly upregulated in a spirotetramat-resistant (SR) strain compared with a susceptible (SS) strain. Furthermore, knockdown of CYP4CJ1, CYP6CY7 and CYP6CY21 increased nymph and adult mortality in the SR strain following exposure to spirotetramat. Drosophila ectopically expressing CYP380C6, CYP4CJ1, CYP6DA2, CYP6CY7 and CYP6CY21 showed significantly decreased mortality after spirotetramat exposure, and CYP380C6, CYP4CJ1 and CYP6CY21 are putative targets of six lncRNAs. Silencing of lncRNAs MSTRG.36649.2/5 and MSTRG.71880.1 changed CYP6CY21 and CYP380C6 expression, altering the sensitivity of the SR strain to spirotetramat. Moreover, MSTRG.36649.2/5 did not compete for microRNA (miRNA) binding to regulate CYP6CY21 expression.

CONCLUSION

Our results confirm that CYP380C6, CYP4CJ1, CYP6DA2, CYP6CY7 and CYP6CY21 are potentially involved in the development of spirotetramat resistance in A. gossypii, and MSTRG.36649.2/5 and MSTRG.71880.1 probably regulate CYP6CY21 and CYP380C6 expression other than through the "sponge effect" of competing for miRNA binding. Our results provide a favorable molecular basis for studying cotton aphid P450 genes and lncRNA functions in spirotetramat resistance development.

Chemosensory Proteins Are Associated with Thiamethoxam and Spirotetramat Tolerance in Aphis gossypii Glover

Chemosensory proteins (CSPs) are a class of transporters in arthropods. Deeper research on CSPs showed that CSPs may be involved in some physiological processes beyond chemoreception, such as insect resistance to pesticides. We identified two upregulated CSPs in two resistant strains of Aphis gossypii Glover. To understand their role in the resistance of aphids to pesticides, we performed the functional verification of CSP1 and CSP4 in vivo and in vitro. Results showed that the sensitivity of the thiamethoxam-resistant strain to thiamethoxam increased significantly with the silencing of CSP1 and CSP4 by RNAi (RNA interference), and the sensitivity of the spirotetramat-resistant strain to spirotetramat increased significantly with the silencing of CSP4. Transgenic Drosophila melanogaster expressing CSPs exhibited stronger resistance to thiamethoxam, spirotetramat, and alpha-cypermethrin than the control did. In the bioassay of transgenic Drosophila, CSPs showed different tolerance mechanisms for different pesticides, and the overexpressed CSPs may play a role in processes other than resistance to pesticides. In brief, the present results prove that CSPs are related to the resistance of cotton aphids to insecticides.

MAPK-mediated transcription factor GATAd contributes to Cry1Ac resistance in diamondback moth by reducing PxmALP expression

The benefits of biopesticides and transgenic crops based on the insecticidal Cry-toxins from Bacillus thuringiensis (Bt) are considerably threatened by insect resistance evolution, thus, deciphering the molecular mechanisms underlying insect resistance to Bt products is of great significance to their sustainable utilization. Previously, we have demonstrated that the down-regulation of PxmALP in a strain of Plutella xylostella (L.) highly resistant to the Bt Cry1Ac toxin was due to a hormone-activated MAPK signaling pathway and contributed to the resistance phenotype. However, the underlying transcriptional regulatory mechanism remains enigmatic. Here, we report that the PxGATAd transcription factor (TF) is responsible for the differential expression of PxmALP observed between the Cry1Ac susceptible and resistant strains. We identified that PxGATAd directly activates PxmALP expression via interacting with a non-canonical but specific GATA-like cis-response element (CRE) located in the PxmALP promoter region. A six-nucleotide insertion mutation in this cis-acting element of the PxmALP promoter from the resistant strain resulted in repression of transcriptional activity, affecting the regulatory performance of PxGATAd. Furthermore, silencing of PxGATAd in susceptible larvae reduced the expression of PxmALP and susceptibility to Cry1Ac toxin. Suppressing PxMAP4K4 expression in the resistant larvae transiently recovered both the expression of PxGATAd and PxmALP, indicating that the PxGATAd is a positive responsive factor involved in the activation of PxmALP promoter and negatively regulated by the MAPK signaling pathway. Overall, this study deciphers an intricate regulatory mechanism of PxmALP gene expression and highlights the concurrent involvement of both trans-regulatory factors and cis-acting elements in Cry1Ac resistance development in lepidopteran insects.

Gene expression and regulation are associated with adaptive evolution in living organisms. The rapid evolution of insect resistance to Bt insecticidal Cry toxins is frequently associated with reduced expression of diverse midgut genes that code for Cry-toxin receptors. Nonetheless, our current knowledge about the regulation of gene expression of these pivotal receptor genes in insects is limited. Membrane-bound alkaline phosphatase (mALP) is a known receptor for Cry1Ac toxin in diverse insects and here, we report the transcriptional regulatory mechanism of the PxmALP gene related to Cry1Ac resistance in P. xylostella. We identified a MAPK signaling pathway that negatively regulates the PxGATAd transcriptional factor which is involved in the differential expression of PxmALP via interacting with the PxmALP promoter. Furthermore, a cis-acting element mutation repressing the regulatory activity of PxGATAd for PxmALP expression in the Cry1Ac resistant strain was identified. Our study provides an insight into the precise transcriptional regulatory mechanism that regulates PxmALP expression and is involved in the evolution of Bt Cry1Ac resistance in P. xylostella, which provides a paradigm for decoding the regulation landscape of midgut Cry-toxin receptor genes in insects.

The genomic basis of evolutionary differentiation among honey bees

In contrast to the western honey bee, Apis mellifera, other honey bee species have been largely neglected despite their importance and diversity. The genetic basis of the evolutionary diversification of honey bees remains largely unknown. Here, we provide a genome-wide comparison of three honey bee species, each representing one of the three subgenera of honey bees, namely the dwarf (Apis florea), giant (A. dorsata), and cavity-nesting (A. mellifera) honey bees with bumblebees as an outgroup. Our analyses resolve the phylogeny of honey bees with the dwarf honey bees diverging first. We find that evolution of increased eusocial complexity in Apis proceeds via increases in the complexity of gene regulation, which is in agreement with previous studies. However, this process seems to be related to pathways other than transcriptional control. Positive selection patterns across Apis reveal a trade-off between maintaining genome stability and generating genetic diversity, with a rapidly evolving piRNA pathway leading to genomes depleted of transposable elements, and a rapidly evolving DNA repair pathway associated with high recombination rates in all Apis species. Diversification within Apis is accompanied by positive selection in several genes whose putative functions present candidate mechanisms for lineage-specific adaptations, such as migration, immunity, and nesting behavior.

Resistance Risk Assessment of the Ryanoid Anthranilic Diamide Insecticide Cyantraniliprole in Aphis gossypii Glover

Cyantraniliprole targets the ryanodine receptor and shows cross-spectrum activity against a broad range of chewing and sucking pests. In this study, a cyantraniliprole-resistant cotton aphid strain (CyR) developed resistance 17.30-fold higher than that of a susceptible (SS) strain. Bioassay results indicated that CyR developed increased cross-resistance to cyfluthrin, α-cypermethrin, imidacloprid, and acephate. In CyR, piperonyl butoxide synergistically increased the toxicity of cyantraniliprole, α-cypermethrin, and cyfluthrin. The cytochrome P450 activities in the CyR strain were significantly higher than those in the SS strain. The mRNA expression of CYP6CY7, CYP6CY12, CYP6CY21, CYP6CZ1, CYP6DA1, and CYP6DC1 in the CYP3 clade, and CYP380C6, CYP380C12, CYP380C44, CYP4CJ1, and CYP4CJ5 in the CYP4 clade, was significantly higher in CyR than in SS. The depletion of the most abundant CYP380C6 transcript by RNAi also significantly increased the sensitivity of CyR to cyantraniliprole. Transgenic expression of CYP380C6, CYP6CY7, CYP6CY21, and CYP4CJ1 in Drosophila melanogaster suggested that the expression of CYP380C6 and CYP4CJ1 was sufficient to confer cyantraniliprole resistance, with CYP380C6 being the most effective, and that CYP380C6, CYP6CY7, and CYP6CY21 were related to α-cypermethrin cross-resistance. These results indicate the involvement of P450 genes in cyantraniliprole resistance and pyrethroid cross-resistance and provide an overall view of the metabolic factors involved in resistance development.

Specific Identification Method based on PCR for Drosophila melanogaster

D. melanogaster is one of the most harmful citrus fruit flies having a large number of host plants. The molecular diagnostic method has been created for identification the D. melanogaster from another non-quarantine species Drosophila spp. The proposed method for differentiation is to use the mitochondrial DNA cytochrome oxidase I gene region 709-bp. We amplified samples of DNA with primers Droso-S391 and Droso-A381 by D. melanogaster, D. suzukii, and D. Simulans collections in the laboratory samples from many countries and contrasted with sequences of other GenBank Drosophila taxa. The findings of a polymerase chain reaction (PCR) based on DNA sequence polymorphisms showed that these primers accurately identify the area of the gene as well as the unique primers of Drosophila melanogaster.

Functional analysis of cytochrome P450 genes linked with acetamiprid resistance in melon aphid, Aphis gossypii

Cytochrome P450 monooxygenases (P450s) are highly conserved multifunctional enzymes that play crucial roles in insecticide resistance development. In this study, the molecular mechanisms of P450s in acetamiprid resistance development to melon aphid, Aphis gossypii was investigated. Acetamiprid resistant (32.64-fold resistance) population (Ace-R) of A. gossypii was established by continuous selection with acetamiprid for 24 generations. Quantitative Real Time PCR was carried out to analyze the expression of P450 genes in both acetamiprid resistant (Ace-R) and susceptible (Ace-S) strains. Result showed that nine genes (CYP6CY14, CYP6DC1, CYP6CZ1, CYP6DD1, CYP6CY5, CYP6CY9, CYP6DA1, CYP6CY18, and CYP6CY16) of CYP3 clade, four genes (CYP302A1, CYP315A1, CYP301A1, and CYP314A1) of CYP2 clade, two genes (CYP4CK1, CYP4G51) of CYP4 clade and three genes (CYP306A1, CYP305E1, CYP307A1) of mitochondrial clade (Mito clad) were significantly up-regulated, in Ace-R compared to Ace-S strain. Whilst CYP4CJ2 gene from (CYP4 clade) was significantly down-regulated in Ace-R strain. Furthermore, RNA interference-mediated knockdown of CYP6CY14, CYP6DC1, and CYP6CZ1 genes significantly increased the sensitivity of Ace-R strain to acetamiprid. Taken together, this study showed that P450 genes especially CYP6CY14, CYP6DC1 and CYP6CZ1 are potentially involved in acetamiprid resistance development in A. gossypii. This study could be useful to understand the molecular basis of acetamiprid resistance mechanism in A. gossypii.

Genomic and transcriptomic analyses in Drosophila suggest that the ecdysteroid kinase-like (EcKL) gene family encodes the 'detoxification-by-phosphorylation' enzymes of insects

Phosphorylation is a phase II detoxification reaction that, among animals, occurs near exclusively in insects, but the enzymes responsible have never been cloned or otherwise identified. We propose the hypothesis that members of the arthropod-specific ecdysteroid kinase-like (EcKL) gene family encode detoxicative kinases. To test this hypothesis, we annotated the EcKL gene family in 12 species of Drosophila and explored their evolution within the genus. Many ancestral EcKL clades are evolutionarily unstable and have experienced repeated gene gain and loss events, while others are conserved as single-copy orthologs. Leveraging multiple published gene expression datasets from D. melanogaster, and using the cytochrome P450s-a classical detoxification family-as a test case, we demonstrate relationships between xenobiotic induction, detoxification tissue-enriched expression and evolutionary instability in the EcKLs and the P450s. We devised a systematic method for identifying candidate detoxification genes in large gene families that is concordant with experimentally determined functions of P450 genes in D. melanogaster. Applying this method to the EcKLs suggested a significant proportion of these genes play roles in detoxification, and that the EcKLs may constitute a detoxification gene family in insects. Additionally, we estimate that between 11 and 16 uncharacterised D. melanogaster P450s are strong detoxification candidates. Lastly, we also found previously unreported genomic and transcriptomic variation in a number of EcKLs and P450s associated with toxic stress phenotypes using a targeted phenome-wide association study (PheWAS) approach in D. melanogaster, presenting multiple future avenues of research for detoxification genetics in this species.

Genetic Diversity in Invasive Populations of Argentine Stem Weevil Associated with Adaptation to Biocontrol

Modified, agricultural landscapes are susceptible to damage by insect pests. Biological control of pests is typically successful once a control agent has established, but this depends on the agent’s capacity to co-evolve with the host. Theoretical studies have shown that different levels of genetic variation between the host and the control agent will lead to rapid evolution of resistance in the host. Although this has been reported in one instance, the underlying genetics have not been studied. To address this, we measured the genetic variation in New Zealand populations of the pasture pest, Argentine stem weevil (Listronotus bonariensis), which is controlled with declining effectiveness by a parasitoid wasp, Microctonus hyperodae. We constructed a draft reference genome of the weevil, collected samples from a geographical survey of 10 sites around New Zealand, and genotyped them using a modified genotyping-by-sequencing approach. New Zealand populations of Argentine stem weevil have high levels of heterozygosity and low population structure, consistent with a large effective population size and frequent gene flow. This implies that Argentine stem weevils were able to evolve more rapidly than their biocontrol agent, which reproduces asexually. These findings show that monitoring genetic diversity in biocontrol agents and their targets is critical for long-term success of biological control.

Transcriptome analysis of Drosophila melanogaster laboratory strains of different geographical origin after long-term laboratory maintenance

Positive selection may be the main factor of the between-population divergence in gene expression. Expression profiles of two Drosophila melanogaster laboratory strains of different geographical origin and long-term laboratory maintenance were analyzed using microchip arrays encompassing probes for 18,500 transcripts. The Russian strain D18 and the North American strain Canton-S were compared. A set of 223 known or putative genes demonstrated significant changes in expression levels between these strains. Differentially expressed genes (DEG) were enriched in response to DDT (p = .0014), proteolysis (p = 2.285E-5), transmembrane transport (p = 1.03E-4), carbohydrate metabolic process (p = .0317), protein homotetramerization (p = .0444), and antibacterial humoral response (p = 425E-4). The expression in subset of genes from different categories was verified by qRT-PCR. Analysis of transcript abundance between Canton-S and D18 strains allowed to select several genes to estimate their participation in latitude adaptation. Expression of selected genes was analyzed in five D. melanogaster lines of different geographic origins by qRT-PCR, and we found two candidate genes that may be associated with latitude adaptation in adult flies-smp-30 and Cda9. Quite possible that several alleles of these genes may be important for insect survival in the environments of global warming. It is interesting that the number of genes involved in local adaptation demonstrates expression level appropriate to their geographical origin even after decades of laboratory maintenance.

A flavin-dependent monooxgenase confers resistance to chlorantraniliprole in the diamondback moth, Plutella xylostella

The diamondback moth, Plutella xylostella, is a damaging pest of cruciferous crops, and has evolved resistance to many of the insecticides used for control, including members of the diamide class. Previous work on the molecular basis of resistance to diamides has documented mutations in the target-site, the ryanodine receptor, in resistant populations of P. xylostella worldwide. In contrast the role of metabolic resistance to this insecticide class is significantly less clear. Here we show that overexpression of a flavin-dependent monooxgenase (FMO) confers resistance to the diamide chlorantraniliprole in P. xylostella. Transcriptome profiling of diamide resistant strains, with and without target-site resistance, revealed constitutive over-expression of several transcripts encoding detoxification enzymes compared to susceptible strains. Two of these, CYP6BG1, and PxFMO2 were particularly highly overexpressed (33,000 and 14,700-fold, respectively) in a resistant strain (HAW) lacking target-site resistance. After 17 generations without diamide selection the resistance of the HAW strain fell by 52-fold and the expression of PxFMO2 by > 1300-fold, however, the expression of CYP6BG1 declined by only 3-fold. Generation of transgenic Drosophila melanogaster expressing these genes demonstrated that PxFMO2, but not CYP6BG1, confers resistance in vivo. Overexpression of PxFMO2 in the HAW strain is associated with mutations, including a putative transposable element insertion, in the promoter of this gene. These enhance the expression of a reporter gene when expressed in a lepidopteran cell line suggesting they are, at least in part, responsible for the overexpression of PxFMO2 in the resistant strain. Our results provide new evidence that insect FMOs can be recruited to provide resistance to synthetic insecticides.

Sublethal larval exposure to imidacloprid impacts adult behaviour in Drosophila melanogaster

Pesticides are now chronically found in numerous ecosystems incurring widespread toxic effects on multiple organisms. For insects, the larvae are very exposed to pesticide pollution and the acute effect of insecticides on larvae has been characterized in a range of species. However, the carry-on effects in adults of sublethal exposure occurring in larvae are not well characterized. Here, we use a collection of strains of Drosophila melanogaster differing in their larval resistance to a commonly used insecticide, imidacloprid, and we test the effect of larval exposure on behavioural traits at the adult stage. Focusing on locomotor activity and on courtship and mating behaviour, we observed a significant carry-on effect of imidacloprid exposure. The heritability of activity traits measured in flies exposed to imidacloprid was higher than measured in controls and in these, courtship traits were genetically less correlated from mating success. Altogether, we did not observe a significant effect of the larval insecticide resistance status on adult behavioural traits, suggesting that selection for resistance in larvae does not involve repeatable behavioural changes in adults. This lack of correlation between larval resistance and adult behaviour also suggests that resistance at the larval stage does not necessarily result in increased behavioural resilience at a later life stage. These findings imply that selection for resistance in larvae as well as for behavioural resilience to sublethal exposure in adult will combine and impose a greater evolutionary constraint. Our conclusions further substantiate the need to encompass multiple trait measures and life stages in toxicological assays to properly assess the environmental impact of pesticides.

Drosophila melanogaster as a powerful tool for studying insect toxicology

Insecticides are valuable and widely used tools for the control of pest insects. Despite the use of synthetic insecticides for >50 years, we continue to have a limited understanding of the genes that influence the key steps of the poisoning process. Major barriers for improving our understanding of insecticide toxicity have included a narrow range of tools and/or a large number of candidate genes that could be involved in the poisoning process. Herein, we discuss the numerous tools and resources available in Drosophila melanogaster that could be brought to bear to improve our understanding of the processes determining insecticide toxicity. These include unbiased approaches such as forward genetic screens, population genetic methods and candidate gene approaches. Examples are provided to showcase how D. melanogaster has been successfully used for insecticide toxicology studies in the past, and ideas for future studies using this valuable insect are discussed.