# Sublethal Exposure Effects of the Neonicotinoid Clothianidin Strongly Modify the Brain Transcriptome and Proteome in the Male Moth Agrotis ipsilon

Camille Meslin, Françoise Bozzolan, Virginie Braman, Solenne Chardonnet, Cédric Pionneau, Marie-Christine François, Dany Severac, Christophe Gadenne, Sylvia Anton, Martine Maibèche, et al.

<figure><img src="image_1.png" /></figure>

To cite this version:

Camille Meslin, Françoise Bozzolan, Virginie Braman, Solenne Chardonnet, Cédric Pionneau, et al. <a href="https://hal.inrae.fr/hal-03138710">https://hal.inrae.fr/hal-03138710</a><br/>Submitted on 11 Feb 2021

HAL Id: hal-03138710

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Sublethal Exposure Effects of the Neonicotinoid Clothianidin Strongly Modify the Brain Transcriptome and Proteome in the Male Moth *Agrotis ipsilon*

Camille Meslin<fnref n="1" />, Françoise Bozzolan<fnref n="1" />, Virginie Braman<fnref n="1" />, Solenne Chardonnet<fnref n="2" />, Cédric Pionneau<fnref n="2" />, Marie-Christine François<fnref n="1" />, Dany Severac<fnref n="3" />, Christophe Gadenne<fnref n="4" />, Sylvia Anton<fnref n="4" />, Martine Maibèche<fnref n="1" />, Emmanuelle Jacquin-Joly<fnref n="1" /> and David Siaussat <fnref n="1" />,*,†

**Citation:** Meslin, C.; Bozzolan, F.; Braman, V.; Chardonnet, S.; Pionneau, C.; François, M.-C.; Severac, D.; Gadenne, C.; Anton, S.; Maibèche, M.; et al. Sublethal Exposure Effects of the Neonicotinoid Clothianidin Strongly Modify the Brain Transcriptome and Proteome in the Male Moth *Agrotis ipsilon*. *Insects* **2021**, *12*, 152. https://doi.org/10.3390/insects12020152

Academic Editors: Ian M. Scott and David Siaussat
Received: 6 January 2021
Accepted: 4 February 2021
Published: 11 February 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

**Simple Summary:** Insect pest management relies mainly on neurotoxic insecticides, including neonicotinoids such as clothianidin. Low doses of insecticides can stimulate various life traits in target pest insects, whereas negative effects are expected. We recently showed that treatments with different low doses of clothianidin could modify behavioral and neuronal sex pheromone responses in the male moth, *Agrotis ipsilon*. In this study, we showed that clothianidin disrupted 1229 genes and 49 proteins at the molecular level, including numerous enzymes of detoxification and neuronal actors, which could explain the acclimatization in pest insects to the insecticide-contaminated environment.

**Abstract:** Insect pest management relies mainly on neurotoxic insecticides, including neonicotinoids such as clothianidin. The residual accumulation of low concentrations of these insecticides can have positive effects on target pest insects by enhancing various life traits. Because pest insects often rely on sex pheromones for reproduction and olfactory synaptic transmission is cholinergic, neonicotinoid residues could indeed modify chemical communication. We recently showed that treatments with low doses of clothianidin could induce hormetic effects on behavioral and neuronal sex pheromone responses in the male moth, *Agrotis ipsilon*. In this study, we used high-throughput RNAseq and proteomic analyses from brains of *A. ipsilon* males that were intoxicated with a low dose of clothianidin to investigate the molecular mechanisms leading to the observed hormetic effect. Our results showed that clothianidin induced significant changes in transcript levels and protein quantity in the brain of treated moths: 1229 genes and 49 proteins were differentially expressed upon clothianidin exposure. In particular, our analyses highlighted a regulation in numerous enzymes as a possible detoxification response to the insecticide and also numerous changes in neuronal processes, which could act as a form of acclimatization to the insecticide-contaminated environment, both leading to enhanced neuronal and behavioral responses to sex pheromone.

<footnote n="1">¹ Département Ecologie Sensorielle, Institut d'Ecologie et des Sciences de l'Environnement de Paris (iEES-Paris), Sorbonne Université, INRAE, CNRS, IRD, UPEC, Université de Paris, 75005 Paris, France or 78026 Versailles, France; Camille.Meslin@inrae.fr (C.M.); Francoise.Bozzolan@sorbonne-universite.fr (F.B.); virginie.braman@upmc.fr (V.B.); marie-christine.francois@inrae.fr (M.-C.F.); martine.maibeche@sorbonne-universite.fr (M.M.); emmanuelle.joly@inrae.fr (E.J.-J.)</footnote>

<footnote n="2">² Sorbonne Université, INSERM, UMS 37 PASS, Plateforme Post-Génomique de la Pitié-Salpêtrière (P3S), 75013 Paris, France; solenne.chardonnet@sorbonne-universite.fr (S.C.); cedric.pionneau@sorbonne-universite.fr (C.P.)</footnote>

<footnote n="3">³ MGX, BioCampus Montpellier, CNRS, INSERM, Université de Montpellier, 34000 Montpellier, France; Dany.Severac@mgx.cnrs.fr</footnote>

<footnote n="4">⁴ Institut de Génétique Environnement et Protection des Plantes IGEPP, INRAE, Institut Agro, Univ Rennes, 49045 Angers, France; antongadenne@orange.fr (C.G.); sylvia.anton@inrae.fr (S.A.)</footnote>

<footnote>* Correspondence: david.siaussat@sorbonne-universite.fr</footnote>

<footnote>† Present address: Faculté des sciences et ingéniérie, Campus Pierre et Marie Curie, Sorbonne Université, 4 Place Jussieu—Tours 44-45, 3ème Étage, Boîte Courrier 237, 75005 Paris, France.</footnote>
**Keywords:** pest insect; clothianidin; proteomics; transcriptomics; hormesis; *Agrotis ipsilon*

# 1. Introduction

Although integrated pest management strategies are increasingly being developed [1], the majority of treatments for pest insects still rely exclusively on the use of neurotoxic chemicals, such as neonicotinoid insecticides [2]. These molecules, including the widely used last-generation insecticide clothianidin, are known to disrupt synaptic transmission through their action on nicotinic acetylcholine receptors [3,4]. The widespread use of these neurotoxic insecticides raises numerous issues, such as residual accumulation in the environment [5], which is present for many years even after stopping treatments, and negative effects on physiology and behavior of non-target insects such as honeybees [6]. Indeed, there is growing evidence that sublethal or low doses of neonicotinoids impact insect physiology and thus vital behaviors such as reproduction or searching for food [7,8].

Contrary to these effects disturbing insect physiology and behavior, low doses of insecticides can also elicit hormetic effects—i.e., enhance certain physiological and behavioral traits. Hormesis is defined as a biphasic response following exposure to a given toxicant with beneficial effects at low-dose exposure and adverse effects at high-dose exposure [9]. Numerous examples of this toxicological phenomenon have been reported for many types of biological and pathological processes in microorganisms, plants, and mammals including humans [10,11]. In insects, insecticide-induced hormesis in developmental and reproductive life traits (such as growth stimulation, enhanced pupation, decrease in pupal mortality, increased fecundity and longevity, and increase in oviposition) has likewise been observed following treatments with different insecticides, including neonicotinoids, carbamates, and organophosphates [12–14]. In addition to the mentioned life traits, insecticides also interfere with chemical communication in insects: they can, for example, disrupt the behavioral response of pest insects to sex pheromones or food odors [15–18]. Recent results from the black cutworm, *Agrotis ipsilon* (Hufnagel) (Lepidoptera: Noctuidae), show that low doses of the neonicotinoid clothianidin induce a biphasic effect on pheromone-guided behavior with a hormetic-like effect [17]. In another moth species, the cotton leafworm *Spodoptera littoralis* (Boisduval) (Lepidoptera: Noctuidae), the same effect has been observed for deltamethrin, and detailed investigations revealed that these effects were restricted to the response of males to the sex pheromone produced by con-specific females, whereas no effect of the same doses was observed on behavioral responses to food odor [16].

Although evidence is accumulating that the concept of hormesis is valid for many insecticides used to control pest insects, very little is known about the molecular and cellular mechanisms leading to this phenomenon [14]. Most likely, the effects are due to changes in the expression of genes involved in basic cellular and physiological processes. The recent progresses in “omic” approaches have opened up possibilities to decipher such changes, especially in non-model species such as crop pest insects. It is now possible to acquire a large amount of data on gene expression (transcriptomics, e.g., RNAseq) or protein levels (proteomics) on a given species and in different experimental conditions, pinpointing any mechanisms or possible regulatory pathways involved, without any a priori study. However, most of such studies have mainly focused on gene expression alone (see, for example, [19]), and rare are those that also take into account protein levels [20].

The well described olfactory plasticity and available molecular data on neuromodulatory mechanisms in *A. ipsilon* males [21,22] make it an excellent model with which to study the molecular mechanisms responsible for the hormetic action of pesticides on behavior and on the central nervous system.

In our previous study, the lethal dose 50 (LD50: dose resulting in 50% mortality in tested insects) was found to be 69 ng/moth in male *A. ipsilon* [17]. At higher doses, intoxicated insects exhibited trembling and incapacity to move before dying. In the low lethal
dose range, 5 ng of clothianidin decreased the proportion of males able to fly at the exception of the LD20 dose (10 ng), which induces a hormetic-like effect. Indeed, at this dose, we observed an improved orientation behavior of males in response to female pheromone, whereas no clothianidin effect was observed on behavioral responses to plant odor [17]. We also discovered that this dose of clothianidin modifies the pheromone response thresholds of central neurons of *A. ipsilon*, but not those of peripheral olfactory receptor neurons [23]. This correlates with the changes in behavioral responses after clothianidin treatment and suggests the antennal lobe—i.e., the part of the brain that processes responses to odorants—as the neural substrate involved in clothianidin-induced behavioral changes. We therefore attempted to explore the molecular mechanisms underlying the hormetic effects within the brain, using combined transcriptomic and proteomic analyses of male *A. ipsilon* brains originating from individuals orally treated with a LD20 dose of clothianidin dissolved in dimethyl sulfoxide (DMSO), compared with untreated individuals and with individuals treated with the solvent DMSO alone in the same way as in Rabhi et al. [17,23]. Concerning oral treatments with insecticides, even though the gut epithelium is supposed to have a barrier role with the presence of a high number and quantity of detoxification enzymes, several studies have shown that the ingestion of pesticides can induce strong effects in the brain, such as the degeneration of neural tissue or the molecular alteration of neuronal actors and degradation/detoxification enzymes [24–26]. Because obvious physiological and behavioral effects in *A. ipsilon* with oral insecticide treatments were observed, we kept the same approach for our molecular analysis.

This approach allowed us to perform an unbiased search for genes and corresponding proteins, which could be modulated by the insecticide treatments and eventually by the solvent treatment. We expected to identify the genes and proteins involved in different physiological processes, among which would be those possibly responsible for the modified neuronal odor responses [23].

## 2. Materials and Methods

### 2.1. Insects

Adult males of *A. ipsilon* Hufnagel originated from a laboratory colony. Wild insects were introduced into the colony each spring to maintain the genetic diversity and overall health of the colony. Insects were reared on an artificial diet [27] in individual cups at 22 °C under a 16 h:8 h light:dark photoperiod until pupation. Pupae were sexed, and males and females were kept separately at 22 °C in an inversed light–dark cycle (16 h:8 h light:dark photoperiod). Newly emerged adults were removed every day from hatching containers and were given access to a 20% sucrose solution ad libitum. The day of emergence was considered as Day 0.

### 2.2. Chemicals

For the stock solution, 25 µg of clothianidin (99% purity, Sigma-Aldrich, Saint-Quentin Fallavier, France) were dissolved in 100 µL of 2% DMSO. The solution was diluted to 1 ng/µL clothianidin with a 20% sucrose solution, with the final concentration of DMSO being 0.008%. Fresh solutions were prepared weekly as needed from frozen aliquots of the stock solution. All the chemicals were purchased from Sigma-Aldrich (Saint-Quentin Fallavier, France), unless stated otherwise.

### 2.3. Clothianidin Intoxication

Oral acute clothianidin intoxications were performed as described previously [17]. Briefly, 4-day-old males were each fed with 10 µL of the freshly prepared clothianidin solution, resulting in a dose of 10 ng of insecticide per insect. We chose 4-day-old males for treatments because the effects of low doses of clothianidin were assessed 24 h after treatment (day 5), which corresponds to the tested age in wind tunnel and electrophysiological experiments and to the optimal age for behavioral responses to sex pheromone in
wind tunnel experiments [27]. Males were individually treated with clothianidin before the onset of the scotophase under a ventilated fume hood, as described previously [17]. In order to control for the effect of DMSO (clothianidin solvent), animals were exposed to DMSO alone in the same way at the same final concentration (0.008%). Once the males had ingested the presented solution, they were placed in individual plastic cups and transferred to a climate chamber different from the rearing room for 24 h until the dissection of the brains for proteomic and transcriptomic analyses. Control and treated males originated from at least 3 different generations. Brains from the 3 groups of *A. ipsilon* males (not treated, DMSO-treated, 10 ng of clothianidin-treated) were independently dissected at day 5 after their emergence in Ringer solution and immediately dipped in Eppendorf vials kept in liquid nitrogen, then stored at -80 °C until further treatment. For RNAseq, 50 brains from a given condition were cumulated to make one sample, and this was replicated 3 times for each condition. In total, 9 samples of 50 brains each were prepared. For proteomics, the same procedure was performed to produce 15 samples (5 replicates for each condition), except that 30 to 50 brains were used for each sample.

## 2.4. *Transcriptomic Analyses*

Total RNA extractions were performed on the 9 brain samples using TRIzol (Invitrogen, Paris, France), according to the manufacturer's protocol. Total RNA samples were stored at -80 °C until library preparation and sequencing. All the samples were processed at the MGX platform (Montpellier, France). All 9 libraries were prepared separately using the TruSeq Stranded mRNA Sample Preparation Kit (Illumina, Paris, France) according to the manufacturer's protocol and sequenced on an Illumina HiSeq2000 to generate paired-end reads of 150 bp. After trimming off the adaptor sequences, raw reads were processed in terms of both their quality and length using Cutadapt [28]. Reads were scanned and trimmed off when a quality score <30 was encountered. Reads with a length <20 bp were discarded. Clean Illumina single-end reads from a previous round of *A. ipsilon* brain sequencing [21] were added for the de novo assembly of the transcriptome, generating 734,263,081 clean paired-end reads and 86,325,883 clean single-end reads that were used for the transcriptome reconstruction using the MIRA assembler v4.0.2 with default parameters [29]. MIRA generated 514,857 contigs, and multiple filtration steps were then applied to reduce the complexity of the de novo transcriptome. First, only contigs with a length >200 bp were kept. Second, CD-HIT [30,31] was applied with default parameters to lower the redundancy. All the Illumina reads were then mapped to the new transcriptome, and only the contigs with an expression >1 fragment per kilobase of exon per million fragments mapped (FPKM) were kept. Finally, only contigs with an open reading frame >30 amino acids were kept, resulting in a final *A. ipsilon* brain transcriptome of 17,986 contigs. The completeness of the transcriptome was assessed using BUSCO v3.0.2 [32] and the Insecta gene reference set. The functional annotation of the contigs was carried out by (1) blastp against the nr database (NR-2016-12-09) and blastx against the Uniprot-sprot database to capture BLAST homologies, (2) running HMMER to identify protein domains [33], (3) running SignalP [34] to predict signal peptides, and (4) running TMHMM v2.0 to predict the transmembrane regions [35]. Gene Ontologies (GO) were mapped to each transcript according to the annotation of their best blast hit by blastp and blastx and assigned to 12,627 contigs. GO Slim annotations were used in order to give a broad overview of the ontology content. Enrichment or depletion for GO categories was determined in comparison to the whole GO-annotated transcriptome using the Fisher exact test and was considered significant when the FDR (False Discovery Rate) was <0.1.
## 2.5. Abundance Estimation and Differential Expression Analysis

All the clean reads from the 9 samples generated in this study were mapped on the assembly using a Bowtie aligner [36]. Transcript abundance was estimated for each sample using RNA-Seq by Expectation Maximization (RSEM) [37] and was measured as the FPKM values. RNAseq counts were normalized between the different samples and replicates using the trimmed mean of M-values normalization method (TMM) [38]. After that step, a quality check was performed to determine if the biological replicates were well correlated for each condition. That quality check revealed that for each condition, one sample did not correlate with the two others. These outliers (DMSO1, clothianidin2 and Control3) were removed from further analyses and the correlations between samples were checked again (Supplementary Data 1). Differentially expressed transcripts were identified using edgeR within the Bioconductor package [39] by taking into account two biological replicates per condition. Genes were considered differentially expressed for an FDR (False Discovery Rate) <0.10.

## 2.6. Proteomics Analysis

Three biological replicates were prepared for each condition and analyzed using a proteomics approach [40]. Proteins were extracted using a polytron in 2D buffer (Urea 7 M; Thiourea 2 M; CHAPS 1%; SB3-10 0.5%; Triton-X100 0.5%; Isobutanol 10%, 25 mM Tris pH 8.8). Protein content was assessed using the Bradford quick start protein assay (BioRad Paris, France). An internal standard was prepared by pooling an equal amount of all samples. Protein labelling was performed using the 3Dye Cy2/3/5 fluor labelling (FluoProbe, Intechim, Paris, France) with 400 pmol of CyDye for 50 µg of protein, incubated 30 min on ice, then quenched with 0.35 mM of Lysine for 10 min. Cy3 and Cy5 were used to label individual samples whereas Cy2 was used to label the internal standard. Following labelling, samples were stored at -80 °C until use. 2DE was run on 24 cm gels in two complementary pH ranges, pH 5-8 and pH 6-9, using commercial strips (GE Healthcare, Paris, France). Each strip was run with 50 µg of internal standard labelled with Cy2 and 50 µg of two different samples labelled with Cy3 and Cy5, respectively. Strips were passively rehydrated overnight directly with the samples diluted in a rehydration buffer (Urea 7 M; Thiourea 2 M; CHAPS 1%; SD3-10 0.5%; Triton-X100 0.5%; Isobutanol 10%, 40 mM DTT, 0.5% ampholites) or for a basic pH range of 6-9 strips, with a rehydration buffer only (Urea 7 M, Thiourea 2 M, CHAPS 4%, Triton X-100 0.05%, glycerol 5% and Destreak 10 mg/mL). Isoelectrofocalisation (IEF) of the pH range 5-8 strips was performed on an Ettan IPGphor (GE Healthcare) as follows: 2 h at 50 V, 2 h at 200 V, 2 h gradient from 200 V to 1000 V, 2 h at 1000 V, 2 h gradient from 1000 V to 10,000 V, 7 h at 10,000 V. For basic pH range 6-9 strips, samples were incorporated by cup-loading during IEF, as described by [41]. The IEF running program was: 7 h at 50 V, 2 h at 500 V, 2 h gradient from 500 V to 1000 V, 7 h gradient from 1000 V to 10,000 V, 2.5 h at 10,000 V. Strips were incubated for 15 min in equilibration buffer (Urea 6 M, Tris pH 8.8 75 mM, Glycerol 26%, SDS 2%) supplemented with 65 mM of DTT, then for 20 min in equilibration buffer supplemented with 135 mM of iodoacetamide. The second dimension was run in 12% acrylamide gels at 30 V for the first hour, then 150 V and 12 mA per gel in a Tris-Glycine buffer. Gel images were acquired on a scanner Ettan DIGE Imager (GE Healthcare). Images were analyzed using Progenesis SameSpots 3.2.3107.24565 (Nonlinear Dynamics). Spots were automatically detected and matched, then manually validated. Relative quantification was performed in each individual gel against its own internal standard, and signal intensities were normalized between all gels based on the internal standard signal of one reference gel. Statistical analyses between each pair of sample groups were performed using ANOVA. Significant spots were selected when the fold change was >1.6 with a *p* value < 0.05 or a fold change >1.3 with a *p* value < 0.01 in at least one of the two comparisons: not treated/clothianidin-treated and DMSO-treated/clothianidin-treated (see Table S1 for details).
Spots of interest were excised, cut into 1 mm³ cubes and destained in potassium ferricyanure 15 mM and sodium thiosulfate 50 mM. The gel pieces were submitted to reduction by incubation in 10 mM of DTT in 50 mM of ammonium bicarbonate (AmBic, Simga aldrci-, Paris, France) for 30 min at 56 °C, then to alkylation by incubation in 50 mM of iodoacetamide in 50 mM of AmBic for 30 min at room temperature. Proteins were digested with 200 ng of trypsin per spot, overnight at 37 °C, and tryptic peptides were collected. The gel pieces were washed twice in 60% acetonitrile (ACN) and 0.1% trifluoroacetic acid (TFA) for 10 min in an ultrasonic bath. Extracted peptides were concentrated in a speed vacuum dryer and resuspended in 5% ACN with 0.1% formic acid and stored at -80 °C until MS analysis. Peptide mixtures were analyzed by LC-MS/MS on a U3000 nanoLC (Thermo, Paris, France) coupled to an HCTultra ion trap (Bruker). Peptides were concentrated and desalted for 5 min on a pre-column RP-C18 (5 mm, 300 µm i.d., 100 Å, Thermo) with mobile phase A (2% ACN/0.1% formic acid) at a flow rate of 20 µL/min, then separated on an analytical column RP-C18 (15 cm, 75 µm i.d., 100 Å, Dionex, Intershim, Paris, France) at a flow rate of 300 nL/min. Elution gradient was run from 2% to 30% of solvent B (95% ACN/0.1% formic acid) in 35 min and 30% to 40% in 5 min. The ion trap was used in the positive mode with the selection of 8 precursors from each MS spectrum for fragmentation by collision-induced dissociation (CID). The capillary voltage was set at 2 kV. Full scan spectra were acquired in the mass range 250 to 1600 m/z and MS/MS spectra were acquired from 100 to 2800 m/z with singly charged ion exclusion, a dynamic exclusion of 30 s, and an isolation width of 4 Da. Raw data were processed using Data Analysis 3.4 (Bruker). Mgf files were generated with a maximum of 5000 compounds with a signal intensity threshold of 100,000 (AU) and spectra deconvolution. Protein identification was performed with X!TandemPipeline 3.4.0 (with X!Tandem search engine software 2015.04.01.1) using a database resulting from the translation of *A. ipsilon* transcriptomic data. Trypsin was selected as the enzyme with a maximum of 1 missed cleavage. Carbamidomethylation of Cys was set as a fixed modification, oxidation of Met as variable modification; MS and MS/MS tolerance at 0.5 Da. At least 2 unique peptides with a *p* value < 0.05 were required for protein validation.

## 3. Results

### 3.1. Brain Transcriptome Assembly and Annotation

We obtained 1462 complete and single-copy (83.3%), 81 complete and duplicated (4.9%), 99 fragmented (6%), and 97 missing BUSCO genes (5.8%). The high number of complete genes that were reconstructed shows the depth and the completeness of the final transcriptome. The combination of both previous and actual data (deeper sequencing, longer reads, paired-end libraries) improves the overall statistics of the de novo assembly, as described in Table 1. Compared with the previously published transcriptome [30], the median contig length is doubled, leading to more complete genes, as reflected in the BUSCO statistics. Contigs are also less fragmented, and fewer genes are missing in the new transcriptome. GO annotations were assigned to 12,627 contigs (70.2% of the transcriptome) by using the annotation of their best blastp and blastx hits.

<table><caption>**Table 1.** *A. ipsilon* brain transcriptome assembly statistics.</caption><thead><tr><td></td><td>Diesner et al.</td><td>This Study</td></tr></thead><tbody><tr><td>Number of contigs</td><td>50,750</td><td>17,986</td></tr><tr><td>Median contig length (nt)</td><td>618 bp</td><td>1233 bp</td></tr><tr><td>Complete BUSCOs</td><td>1062 (64.0%)</td><td>1462 (88.2%)</td></tr><tr><td>Complete and single-copy BUSCOs</td><td>891 (53.7%)</td><td>1381 (83.3%)</td></tr><tr><td>Complete and duplicated BUSCOs</td><td>171 (10.3%)</td><td>81 (4.9%)</td></tr><tr><td>Fragmented BUSCOs</td><td>309 (18.6%)</td><td>99 (6.0%)</td></tr><tr><td>Missing BUSCOs</td><td>287 (17.4%)</td><td>97 (5.8%)</td></tr></tbody></table>
3.2. *DMSO Exposure Affects Both A. ipsilon Brain Transcriptome and Proteome Profiles*

The transcriptomic analysis revealed that 3938 transcripts were deregulated after DMSO exposure compared to the control group (Table S2). Even though no significant effect of DMSO had been shown on flight activity and pheromone responses in our earlier study [17], the GO analysis of the genes showed that many biological processes were al- tered following DMSO exposure (Figure 1). The “small molecule metabolic process” cat- egory was enriched for the genes deregulated by DMSO while all the other categories were depleted, suggesting that DMSO was slowing down cell metabolism. DMSO also affected the brain proteome, as 51 protein spots presented a significant difference in in- tensity in the comparison between DMSO and control samples out of the 1760 that were detected after gel migration (Figure 2 and grey circle in Figure 3).

<figure><img src="image_4.png" /></figure>

**Figure 1.** DMSO has an impact on many cellular processes. Gene Ontology Enrichment analysis showing the GO terms of transcripts deregulated after DMSO exposure compared to the control group (test set, black) and that exhibit statistically significant differences with the whole transcrip- tome (reference set, grey) (Fisher's exact test, FDR < 0.1). The X-axis shows the percentage of se- quences in each dataset. The Y-axis shows the GO terms. The expression of multiple genes with GO terms reflecting basic cellular processes is affected by DMSO exposure. The GO category “small molecule metabolic process” is enriched in genes deregulated by DMSO (black), while all the other categories are depleted compared to the whole transcript set (light grey).

<figure><img src="image_5.png" /></figure>

**Figure 2.** Proteomic profiles showing protein abundance between control and DMSO-treated brain groups. (A). Bidimen- sional gel electrophoresis separation of the brain of *Agrotis ipsilon* showing about 1292 protein spots after a pre-migration in the pH range 5–8 (Colloidal Coomassie Blue G250). (B). Bidimensional gel electrophoresis separation of the brain of *Agrotis ipsilon* showing about 468 protein spots after a pre-migration in the pH range 6–9 (Colloidal Coomassie Blue G250).
Numbers and colored spots were defined by the software in order to visually identify and number the spots presenting a
difference in protein quantity between two conditions.

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Differential expression analysis between the 3 experimental conditions for both transcriptomic and proteomic analysis. Venn diagram of transcriptomic and proteomic analyses of *A. ipsilon* brains in the 3 experimental conditions (controls, DMSO-treated, and clothianidin-treated). For transcriptomic results (values in bold), differentially expressed genes are indicated (FDR < 0.01). For proteomic results (values in italics), spots showing a significant difference in protein ratio values ≤-1.25 or ≥1.25 (Student T-Test) are indicated. The grey area indicates proteins/transcripts for which a DMSO effect was observed, while the blue area indicates a clothianidin effect. Blue numbers outside of the Venn diagram are the total numbers of transcripts and proteins for which the expression is modified following clothianidin exposure without the DMSO effect.</figcaption></figure>

3.3. *Clothianidin Exposure Affects Protein Biosynthesis and Metabolic Processes*

In order to identify cellular processes related to clothianidin exposure only, we con-
sidered the transcripts and proteins that had their expression modified following the clo-
thianidin exposure (in blue, Figure 3) and excluded transcripts and proteins showing a
DMSO effect (in grey, Figure 3). By doing so, 1229 transcripts and 49 protein spots were
identified and considered differentially expressed in relation with clothianidin exposure
(in blue, Figure 3).

Among these 1229 deregulated transcripts upon clothianidin exposure, 707 were up-
regulated while 522 were downregulated compared to both controls (control and DMSO,
Table S3). Gene ontology enrichment analyses were performed on both categories of tran-
scripts and revealed that clothianidin exposure led to (1) the downregulation of transla-
tion (GO categories “translation”, “ribosome biogenesis”, “ribonucleoprotein complex as-
sembly” as well as “small molecule metabolic processes” and “generation of precursor
metabolites and energy”) (Figure 4) and to (2) the upregulation of genes involved in par-
ticular in two categories that are relevant to our study: “signal transduction” and “neuro-
logical system process” (Figure 4).

Among the 49 spots showing a difference in intensity due to clothianidin, 35 spots
presented a significant difference of intensity in the comparison between clothianidin-
treated and control samples and 4 in the comparison between clothianidin-treated and
DMSO samples (Figure 3; *p* < 0.05, Student Test). Ten spots presented a significant differ-
ence of intensity in both comparisons. Out of these 49 spots, 27 proteins were identified
by LC-MS/MS procedure and placed in groups based on their known function (Figure 5
and Table 2): 6 proteins belonged to energy metabolism; 4 proteins to protein biosynthesis;
4 proteins to proteolysis processes; 3 proteins to transport processes; 3 to valine, thymine,
and pyrimidine metabolism; 2 to neuronal processes; 2 to isoprenoid biosynthetic pro-
cesses; 2 to xenobiotic enzyme families; and 1 to glutamate metabolism. The peptide frag-
ment of an unknown protein was recognized during the LC-MS/MS procedure, but no
correspondence was found in online databases and none could be deduced from the brain
transcriptome. Clothianidin treatment induced a significant decrease in the quantity of 21
of the identified proteins (from 17 to 63% in comparison to controls, Figure 5), mainly affecting energy metabolism, neuronal and isoprenoid biosynthetic processes, xenobiotic metabolizing enzymes, and glutamate metabolism. A significant increase in 6 proteins (from 17 to 32% in comparison to controls) was reported mainly in proteolysis and, to a lesser extent, in energy metabolism; transport; and valine, thymine, and pyrimidine metabolism.

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Exposure to clothianidin deregulates gene expression in the brain of *A. ipsilon*. Gene Ontology Enrichment analysis showing the GO terms of genes deregulated by clothianidin exposure compared to DMSO and controls (test set, blue) and that exhibit statistical significant differences with GO terms of to the whole transcriptome (reference set, grey) (Fisher exact test, FDR < 0.1). The X-axis shows the percentage of sequences in each dataset. The Y-axis shows the GO terms. Upregulated genes are shown on the right of the figure, downregulated genes on the left. The two categories that are relevant to our study, “signal transduction” and “neurological system process”, are highlighted with a red box.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Functional categories of proteins displaying a change in their clothianidin/control level ratio. Proteins of control individuals with or without DMSO and clothianidin-treated individuals showing a significant variation in quantity were identified by LC-MS/MS procedures and classified into 10 functional categories. The name of the protein and its reference in the online database are indicated for each identified protein.</figcaption></figure>
<table><caption>**Table 2.** Characterization of differentially abundant proteins in the brain of *A. ipsilon* following clothianidin exposure.</caption><thead><tr><th>Functional Categories</th><th>Protein Name</th><th>Uniprot Reference</th><th>Species</th><th>Peptide Number</th><th>p Value*</th><th>Protein Quantity in Comparison to Control (%)</th></tr></thead><tbody><tr><td>Glutamate metabolism</td><td>Glutamate dehydrogenase</td><td>A0A0G3VFZ9_CHISP</td><td>Chilo suppressalis</td><td>4</td><td>0.009</td><td>-17.1509</td></tr><tr><td rowspan="2">Xenobiotic enzymes</td><td>Dimeric dihydrodiol dehydrogenase</td><td>S4P4A7_9NEOP</td><td>Pararge aegaria</td><td>8</td><td>0.026</td><td>-62.4242</td></tr><tr><td>Aldehyde dehydrogenase X. mitochondrial-like</td><td>AL1B1_HUMAN</td><td>Amyelois transitella</td><td>6</td><td>0.019</td><td>-18.8077</td></tr><tr><td rowspan="2">Isoprenoid biosynthetic process</td><td>Geranylgeranyl diphosphate synthase</td><td>U3RD44_CHOFU</td><td>Choristoneura fumiferana</td><td>15</td><td>0.010</td><td>-21.7245</td></tr><tr><td>Farnesyl diphosphate synthase 2</td><td>A5A7A5_BOMMO</td><td>Bombyx mori</td><td>7</td><td>0.026</td><td>-62.4242</td></tr><tr><td rowspan="3">Neuronal processes</td><td>Palmitoyl-protein Thioesterase 2</td><td>KGM_213554</td><td>Danaus plexipus</td><td>4</td><td>0.019</td><td>-18.8077</td></tr><tr><td>Fatty alcohol acetyltransferase</td><td>A0A088M9K3_AGRSE</td><td>Agrotis segetum</td><td>13</td><td>0.006</td><td>-43.1842</td></tr><tr><td>N-ethylmaleimide sensitive fusion protein</td><td>Q5G1P2_HELM</td><td>Helicoverpa armigera</td><td>11</td><td>0.035</td><td>-15.9046</td></tr><tr><td rowspan="3">Valine, thymine and pyrimidine metabolism</td><td>Methylmalonate-semialdehyde dehydrogenase [acylating]. mitochondrial isoform X2</td><td>MMSA_DROME</td><td>Papilio polytes</td><td>2</td><td>0.019</td><td>-18.8077</td></tr><tr><td>3-hydroxyisobutyryl-CoA hydrolase. mitochondrial</td><td>A0A068FL83_AGRSE</td><td>Agrotis segetum</td><td>7</td><td>0.026</td><td>-62.4242</td></tr><tr><td>Dihydropyrimidinase isoform X2</td><td>A0A170WSE7_TRIIIF</td><td>Bombyx mori</td><td>3</td><td>0.028</td><td>16.72897</td></tr><tr><td rowspan="3">Transport</td><td>Transferrin</td><td>A7IT76_SPOLT</td><td>Spodoptera litura</td><td>2</td><td>0.004</td><td>-18.3537</td></tr><tr><td>Apolipophorins</td><td>APLP_DROME</td><td>Papilio polytes</td><td>3</td><td>0.028</td><td>16.72897</td></tr><tr><td>Rabenosyn-5</td><td>A0A194QJD8_PA PXU</td><td>Papilio xuthus</td><td>6</td><td>0.019</td><td>-18.8077</td></tr><tr><td rowspan="3">Proteolysis</td><td>Protease m1 zinc metalloprotease</td><td>G6D6V6_DANPL</td><td>Danaus plexipus</td><td>23</td><td>0.017</td><td>28.09224</td></tr><tr><td>Protease m1 zinc metalloprotease</td><td>G6D6V6_DANPL</td><td>Danaus plexipus</td><td>20</td><td>0.030</td><td>19.18782</td></tr><tr><td>Mitochondrial processing peptidase beta subunit</td><td>I4DK27_PAPXU</td><td>Papilio xuthus</td><td>9</td><td>0.006</td><td>-43.1842</td></tr><tr><td rowspan="4">Protein biosynthesis</td><td>Elongation factor 1-a</td><td>L7WID6_SPOLT</td><td>Spodoptera litura</td><td>4</td><td>0.029</td><td>-23.5498</td></tr><tr><td>Elongation factor 1 alpha</td><td>E0D8P8 LOCM</td><td>Locusta migratoria</td><td>4</td><td>0.029</td><td>-23.5498</td></tr><tr><td>UDP-N-acetylglucosamine-peptide N-acetylglucosaminyltransferase 110 kDa subunit</td><td>OGT1_HUMAN</td><td>Bombyx mori</td><td>2</td><td>0.013</td><td>-16.5289</td></tr><tr><td>Translation initiation factor 2 gamma subunit</td><td>Q684K4_9NEOP</td><td>Scoliopteryx li-batrix</td><td>3</td><td>0.029</td><td>-23.5498</td></tr><tr><td rowspan="6">Energy metabolism</td><td>ATP synthase subunit b. mitochondrial</td><td>A0A194PTU2_PA PXU</td><td>Papilio xuthus</td><td>2</td><td>0.034</td><td>-33.3046</td></tr><tr><td>ATP synthase</td><td>Q2F5T4_BOMMO</td><td>Bombyx mori</td><td>5</td><td>0.020</td><td>-18.8077</td></tr><tr><td>Aconitate hydratase. mitochondrial</td><td>Q86GF8_ANTYA</td><td>Antheraea ya-mamai</td><td>5</td><td>0.004</td><td>-18.3537</td></tr><tr><td>Pyruvate kinase</td><td>H9IZ23_BOMMO</td><td>Bombyx mori</td><td>5</td><td>0.009</td><td>-17.1509</td></tr><tr><td>NADH-ubiquinone oxidoreductase 75 kDa subunit. mitochondrial</td><td>NDUS1_DROME</td><td>Bombyx mori</td><td>28</td><td>0.028</td><td>16.72897</td></tr><tr><td>UDP-glucose pyrophosphorylase</td><td>A0A0S1MMM3_A NTPE</td><td>Antheraea pernyi</td><td>5</td><td>0.009</td><td>-17.1509</td></tr><tr><td>Unknown</td><td>unknown</td><td></td><td></td><td>17</td><td>0.003</td><td>-28.2986</td></tr></tbody></table>

Legends: ↗ more abundant in clothianidin-treated group compared to controls; ↘ less abundant in clothianidin-treated group compared to controls; *p* value of the Student's *t*-test analysis.
### 3.4. Clothianidin Exposure Affects Detoxification Enzymes

Transcriptomic and proteomic analyses revealed that 10 genes and 2 proteins related to detoxification enzymes were regulated by the clothianidin treatment within *A. ipsilon* brains. Among them, one cytochrome P450 (CYPs) was upregulated, whereas five other CYPs, one Glutathione-S-transferase (GST), one aldo-keto reductase (AKR), and one carboxylesterase (CCE) were downregulated (Tables 3 and S3). The quantity of two xenobiotic metabolizing enzymes also decreased following clothianidin exposure: one aldehyde dehydrogenase (ALDH) and a dimeric dihydrodiol dehydrogenase (Figure 5).

<table><caption>**Table 3.** Clothianidin exposure affects the transcript abundance of detoxification enzymes.</caption><thead><tr><th>*</th><th>Transcript ID</th><th>logFC<br>Clothianidin vs. DMSO</th><th>logFC<br>Clothianidin vs. Control</th><th>Annotation by Blast Research</th></tr></thead><tbody><tr><td rowspan="2">Up</td><td>DN9632_c0_g1_i1</td><td>1.35</td><td rowspan="2">1.13</td><td>Cytochrome P450 4C1</td></tr><tr><td>DN9385_c0_g1_i1</td><td>0.78</td><td>UDP-glucuronosyltransferase 2B7</td></tr><tr><td rowspan="10">Down</td><td>DN27963_c0_g1_i1</td><td></td><td>-1.49</td><td>Probable cytochrome P450 6a13</td></tr><tr><td>DN20937_c0_g1_i1</td><td></td><td>-1.42</td><td>Cytochrome P450 4c3</td></tr><tr><td>DN37703_c0_g1_i1</td><td></td><td>-1.16</td><td>Cytochrome P450 6B7</td></tr><tr><td>DN34383_c0_g1_i1</td><td></td><td>-0.96</td><td>Cytochrome P450 4C1</td></tr><tr><td>DN20657_c0_g1_i1</td><td></td><td>-0.69</td><td>Aldo-keto reductase AKR2E4</td></tr><tr><td>DN5740_c0_g1_i1</td><td></td><td>-0.63</td><td>Glutathione S-transferase 1</td></tr><tr><td>DN4625_c0_g1_i1</td><td></td><td>-0.35</td><td>Cytochrome P450 CYP12A2</td></tr><tr><td>DN4767_c0_g1_i1</td><td></td><td>-0.50</td><td>Venom carboxylesterase-6</td></tr></tbody></table>

### 3.5. Clothianidin Exposure Affects Neuronal Processes

The transcriptomic analysis revealed significant differences in the expression of genes that are involved in synaptic function after clothianidin exposure (Tables 4 and S3). Indeed, we observed an upregulation of the gene expression of a neurobeachin and synaptotagmin and a downregulation of the synaptic vesicle glycoprotein 2B. Other neuronal processes are also disrupted by clothianidin exposure according to the transcriptomic analysis (Tables 4 and S3). Genes involved in the regulation of exchanges through the neuronal membrane have also their expression modified, such as the tyrosine-protein phosphatase 69D, or the flow of Ca²⁺ and/or Na⁺ through the activity of a transient receptor potential cation channel (trpm), a sodium channel protein (para type), a chloride channel, a voltage-dependent calcium channel (type A), or a sodium/calcium exchanger (type 3). For all these genes, we noted a significant upregulation of gene expression in insects exposed to clothianidin. Interestingly, we also observed the upregulation of other regulatory systems of neuronal activities (Tables 4 and S3), such as, for example, the Gamma-amino-butyric acid (GABA) system. Proteomic analyses revealed that the quantities of three proteins involved in neuronal processes were decreased after clothianidin exposure: a N-ethylmaleimide sensitive fusion protein, also known as N-ethylmaleimide-sensitive factor (NSF); a Palmitoyl-protein Thioesterase 2 (PPT2); and a fatty alcohol acetyltransferase (FAA) (Figure 5).

<table><caption>**Table 4.** Clothianidin exposure affects the transcript abundance of neuronal processes.</caption><thead><tr><th rowspan="2">GO Function</th><th rowspan="2">* Transcript ID</th><th>logFC</th><th>logFC</th><th rowspan="2">Annotation by Blast Research</th></tr><tr><th>Clothianidin vs. DMSO</th><th>Clothianidin vs. Control</th></tr></thead><tbody><tr><td rowspan="4">Synaptic function</td><td>Up DN6306_c0_g1_i1</td><td>0.89</td><td></td><td>Neurobeachin</td></tr><tr><td>Up DN2580_c0_g1_i1</td><td>0.42</td><td></td><td>Synaptotagmin 1</td></tr><tr><td>Up DN4891_c0_g1_i1</td><td>0.33</td><td></td><td>Neurocalcin homolog</td></tr><tr><td>Down DN10953_c0_g1_i1</td><td>-0.53</td><td></td><td>Synaptic vesicle glycoprotein 2B</td></tr><tr><td rowspan="3">Exchanges across the neuronal membrane</td><td>Up DN21572_c0_g1_i1</td><td>1.14</td><td></td><td>Tyrosine-protein phosphatase 69D</td></tr><tr><td>Up DN19704_c0_g1_i1</td><td>0.46</td><td></td><td>Transient receptor potential cation channel trpm</td></tr><tr><td>Up DN6537_c0_g1_i1</td><td>0.56</td><td></td><td>Sodium channel protein para</td></tr></tbody></table>

<footnote>* Up- or downregulated following clothianidin exposure.</footnote>
<table><tbody><tr><td></td><td></td><td>DN12028_c0_g1_i1</td><td>0.53</td><td>Chloride channel protein 2</td></tr><tr><td></td><td></td><td>DN27405_c0_g1_i1</td><td>1.21</td><td>Voltage-dependent calcium channel type D subunit alpha-1</td></tr><tr><td></td><td></td><td>DN1850_c0_g1_i1</td><td>1.59</td><td>Sodium/calcium exchanger 3</td></tr><tr><td rowspan="2">Neuronal activity</td><td>Up</td><td>DN6443_c0_g1_i1</td><td>0.49</td><td>Sodium- and chloride-dependent GABA transporter 1</td></tr><tr><td>Down</td><td>DN392_c0_g1_i1</td><td>-0.31</td><td>Sodium- and chloride-dependent GABA transporter 1</td></tr><tr><td rowspan="5">Neuronal processes</td><td></td><td>DN20357_c0_g1_i1</td><td>0.64</td><td>Mushroom body large-type Kenyon cell-specific protein 1</td></tr><tr><td></td><td>DN7325_c0_g1_i1</td><td>0.96</td><td>Krueppel homolog 1</td></tr><tr><td>Up</td><td>DN25883_c0_g1_i1</td><td>1.36</td><td>Octopamine receptor beta-3R</td></tr><tr><td></td><td>DN49002_c0_g1_i1</td><td>0.84</td><td>Glutamate-gated chloride channel</td></tr><tr><td></td><td>DN18143_c0_g1_i1</td><td>0.48</td><td>Neuropeptide CCHamide-2 receptor</td></tr><tr><td rowspan="2">Neurogenesis/axonal growth</td><td>Up</td><td>DN33688_c0_g1_i1</td><td>0.88</td><td>Lachesin</td></tr><tr><td></td><td>DN10822_c0_g1_i1</td><td>0.49</td><td>Protein sickie</td></tr></tbody></table>

* Up- or downregulated expression following clothianidin exposure.

## 4. Discussion

In this study, we showed that exposure to a low dose of clothianidin strongly modified gene and protein expression in the brain of *A. ipsilon* males. However, this is also the case for DMSO. DMSO is often used as a solvent or vehicle solution since this product is a good solvent for hydrophobic molecules such as clothianidin. The selection of appropriate vehicles or solvents to administer the compounds of interest is determinant for the quality of the results. There is now more and more evidence that solvents affect organisms. This also includes DMSO. Whether in vertebrates or invertebrates, DMSO effects can be deleterious, with the observation of some cases of cell toxicity [42], sterility [43], or neurotoxicity [44], but also beneficial, with some report of the neutralization or attenuation of pathologies [45] or with no observed effect [46]. Our results confirm that DMSO, a solvent used in many studies, even at very low doses, can influence gene and protein expression and that these effects need to be taken into account in differential studies.

After the exclusion of transcripts and proteins showing a DMSO effect, we thus identified 1229 transcripts and 49 protein spots considered differentially expressed in relation to clothianidin exposure. Interestingly, although we were able to identify common regulated categories such as neuronal processes using transcriptomics and proteomics, we did not identify the same genes/proteins. Additionally, genes of these categories were upregulated while protein levels were decreased following clothianidin exposure. While this may be surprising at a first glance, multiple studies showed that levels of mRNA and proteins are often not well correlated. This is due to numerous reasons, such as translational efficiency, mRNA structure, and stability [47] but also due to the difference in sensitivity between proteomics and transcriptomics methods. This work is another evidence to support the fact that these two high-throughput omics approaches are complementary for the identification of precise and complete pathways [48].

We have been able to find a wide variety of families of molecular actors and more specifically detoxification enzymes which show that clothianidin, after having crossed the digestive barrier, or its bio-transformed metabolites induce a detoxifying system response in the brain. All these enzyme super families are crucial in insects for their diverse roles on endogenous and both natural or synthetic exogenous compounds. In particular, insect CYPs, CCEs and GSTs play a prominent role in xenobiotic metabolism and many of them are involved in insecticide detoxification/resistance mechanisms [49]. Metabolic detoxification can be divided into phase I processes, consisting mainly of hydrolysis and oxidation reactions (by CYPs, AKRs or ALDHs), and phase II processes, involving the conjugation of phase I products with endogenous compounds (mainly by UGTs and GSTs), leading to the production of hydrophilic metabolites excreted outside the cells.

Neonicotinoids are complex molecules that could be metabolized by the combined action of various phase I and II enzymes [50,51]. The involvement of CYPs, CCEs, and GSTs in the neonicotinoid metabolism and/or neonicotinoid resistance in insects has been well documented [50,52]. The induction of these enzyme genes after neonicotinoid exposure and/or enhanced neonicotinoid detoxification associated with these enzymes were
receptor potential cation channel (trpm) was realized in rat primary cortical neurons ex-
posed to rotenone, an insecticide disrupting the energetic metabolism and inducing oxi-
dative stress [75]. The authors observed an increase in the expression of the trpm2 isoform
after exposure to rotenone, showing that there is a link between dysfunction of TRP chan-
nels, calcium dyshomeostasis and oxidative stress induced by insecticides.

Interestingly, we also observed the upregulation of the beta subunit of the Gamma-
aminobutyric acid (GABA) receptor. The Gamma-aminobutyric acid (GABA) system is
often associated to insecticide response since it appeared to be modulated from the reab-
sorption of the neurotransmitter (i.e., GABA) with a sodium- and chloride-dependent
GABA transporter, to its reception by the GABA receptor, a receptor that may be involved
in resistance to other insecticides, such as cyclodiene in *D. melanogaster* [76]. In fact, as we
observed mainly an upregulation of genes involved in neuronal processes, we can hy-
pothesize that this is the consequence of a general acclimatization of the neuronal system
to a low lethal dose (LD<sub>20</sub>) of clothianidin. The upregulation of mushroom body large-type
Kenyon cell-specific protein 1, Krueppel homolog 1, Octopamine receptor beta-3R, Gluta-
mate-gated chloride channel or Neuropeptide CCHamide-2 receptor could also support
this hypothesis. These actors are known to allow hormonal or neuropeptide modulation
of neuronal activity, and they were previously described to be regulated by insecticide
exposure [77,78]. Finally, some of our results also suggest that the low lethal dose expo-
sure induced some neurogenesis or axonal growth, since we observed a significant in-
crease in lachesin, an immunoglobulin superfamily protein whose expression correlates
with neurogenesis [79] and of the SICKIE protein, which regulates F-actin mediated ax-
onal growth in Drosophila mushroom body neurons [80].

Finally, we were able to observe effects on candidates still little studied but which potentially seem to play an important role in neuronal functioning. It is the case of PPT2 which is involved in the Palmitoylation cycle, a post-translation modification that can occur on secreted, cytoplasmic, and integral membrane proteins. PPT1 and PPT2 are the two primary thioesterases involved in removing palmitoyl groups during the lysosomal degradation process. While PPT1 is associated with neuronal disorders [81], little is known about the function of PPT2, despite its crucial role in the search for therapeutic solutions against neurodegenerative diseases in humans [82]. The difference of identity percentages between the two proteins and the inability of PPT2 to rescue the neural disruption associated with loss of PPT1 support distinct functions and substrates for these two thioesterases [83,84]. In insects, even less information is available on these proteins. In *D. melanogaster*, PPT1 plays a role in neuronal development and function [85,86]. However, the function of PPT2 remains elusive. Another interesting candidates is FAA which is often associated with sex pheromone biosynthesis, which takes place in pheromonal glands of insects. The FAA we found with our proteomic approach is very similar to the one found in a study on the role of enzymes in the pheromone biosynthesis in *Agrotis segetum* [87]. However, this FAA, as with the 33 other acetyltransferases found in the transcriptome of the pheromonal glands, appeared not to be involved in pheromone biosynthesis. Unfortunately, there is currently no information about the role of FAA in the nervous system in insects or vertebrates. Interestingly, this protein, as PPT proteins, was classified within a family of proteins having a palmitoyl-(protein) hydrolase activity (GO molecular function—UniProtKB—A0A088M9W1).

Whereas many genes/proteins involved in neuronal modulation were activated by
the low dose of clothianidin, we did not find changes in the expression of the acetylcholine
receptors themselves, the prime target for neonicotinoids such as clothianidin. Indeed,
other ways of neuronal modulation were activated. However, we cannot conclude
whether this large range of effects is due to a direct effect of clothianidin, or of its bio-
transformed metabolites, or of a general response due to the stress induced by the insec-
ticide exposure (i.e., oxidative stress, increase in metabolic cost...). Nevertheless, these
neuronal changes could provide an explanatory element for the behavioral and neuronal
sensitivity modifications that we observed previously in clothianidin-treated males [17,23,88].

The results of our study are consistent with the few studies showing anatomical or molecular modifications in the brain after the ingestion of pesticides in insects [24]. Studies on the capacity and limits of the digestive system to manage exposure to pesticides still need to progress since many processes remain unclear (passage of molecules, production of toxic metabolites, etc.). These effects highlight the fact that the digestive system does not represent a strict barrier to toxic compounds and that indeed insecticide effects occur in the brain after oral application. High-throughput RNAseq and proteomic analyses in our study showed a regulation of the expression of numerous enzymes as a possible detoxification response to the insecticide and also numerous changes in neuronal processes, which could act as a form of acclimatization to the insecticide-polluted environment. Functional studies will now be needed to investigate how up- or down-regulation of the differentially expressed genes and proteins might be involved in the behavioral and neuronal hormetic effects of increased sex pheromone responses in male *A. ipsilon*. An increased number of ecotoxicological studies, including the present one, now combine data from multiple omics techniques (Table 5). The complementarity of all these techniques permits more generally the identification of numerous targets that could be used for functional studies to disentangle the molecular mechanisms of pollutant effects on insects and other invertebrates.

<table><caption>**Table 5.** Ecotoxicological studies using transcriptomics and proteomics methods.</caption><thead><tr><th colspan="3">Proteins</th><th colspan="3">Transcripts</th><th rowspan="2">Molecules or Condition</th><th rowspan="2">Species</th><th rowspan="2">References</th></tr><tr><th>Number of Detected Spots or Proteins</th><th>Proteins Showing Significant Variation in Quantity</th><th>Up</th><th>Down</th><th>Number of Analyzed Unigenes</th><th>Up</th><th>Down</th></tr></thead><tbody><tr><td>1760</td><td>49</td><td>6</td><td>23</td><td>17,986</td><td>2292</td><td>1646</td><td>Clothianidin</td><td>Agrotis ipsilon</td><td>Present study</td></tr><tr><td rowspan="3">700</td><td rowspan="3">12</td><td rowspan="3">nc</td><td rowspan="3">nc</td><td>nc</td><td>41</td><td>56</td><td>Imidacloprid</td><td rowspan="3">Mytilus galloprovincialis</td><td rowspan="3">[20]</td></tr><tr><td>nc</td><td>43</td><td>37</td><td>Thiacloprid</td></tr><tr><td>nc</td><td>26</td><td>23</td><td>Mix</td></tr><tr><td rowspan="3">1005</td><td rowspan="3">52</td><td rowspan="3">38</td><td rowspan="3">14</td><td>29,146 to 31,467</td><td>646 to 658</td><td>284 to 533</td><td>Clothianidin</td><td>Bradysia odoriphaga</td><td>[66]</td></tr><tr><td>35,222</td><td>349</td><td>271</td><td>Thiamethoxam</td><td>Aphis gossypii</td><td>[89]</td></tr><tr><td>11,150 to 11,426</td><td>225</td><td>384</td><td>Thiamethoxam</td><td>Apis mellifera</td><td>[64]</td></tr><tr><td>821</td><td>143</td><td>35</td><td>108</td><td>nc</td><td>664</td><td>674</td><td>Thiamethoxam</td><td>Bemisia tabaci</td><td>[65]</td></tr><tr><td rowspan="2">&gt;1300</td><td rowspan="2">28</td><td rowspan="2">14</td><td rowspan="2">14</td><td rowspan="2"></td><td rowspan="2"></td><td rowspan="2"></td><td>Mesoionic pyr-ido[1,2]-pyrimidinone</td><td>Aphis craccivora</td><td>[90]</td></tr><tr><td>resistant strain to imidacloprid</td><td>Myzus persicae</td><td>[91]</td></tr><tr><td>1470</td><td>155</td><td>96</td><td>59</td><td></td><td></td><td></td><td>Nicotine</td><td>Apis mellifera</td><td>[92]</td></tr></tbody></table>

## 5. Conclusions

Our study offers potential molecular explanations for the hormetic response that we previously observed in adult *A. ipsilon* males. Testing and understanding the effects of low doses is particularly important because these doses can have unexpected stimulating effects on pests, for which alternative management solutions are sought. In the context of integrated control measures, this information is therefore particularly crucial in order to avoid potentially counterproductive strategies (e.g., reduction in the quantities used) or to choose new alternative control methods.

**Supplementary Materials:** The following are available online at www.mdpi.com/2075-4450/12/2/152/s1: Supplementary Data 1. Fragment counts were used to build a correlation matrix between all conditions of the experiments to perform a quality check. The biological replicates DMSO1, clothianidin2, and Control3 did not correlate with the other biological replicates of their exposure groups and were removed from subsequent analysis (A). After the removal of these three samples, the two remaining biological replicates for each exposure group are correlated to each
other (B). Table S1: Proteomics results for the two pH migrations. The spots for which the ANOVA result is <0.05 and the ratio is ≥1.3 were coloured in green for the CTRL/Clothia and Clothia/DMSO comparisons and in red for the CTRL/DMSO comparison. Table S2: List of deregulated genes following DMSO exposure. Table S3: List of deregulated genes following clothianidin exposure.

**Author Contributions:** C.M.: Bioinformatic analysis and data curation, formal analysis, prepared graphs and statistical analyses, writing—review and editing; F.B.: Methodology; V.B.: Methodology; S.C.: Methodology, prepared graphs and statistical analyses; C.P.: Methodology, prepared graphs and statistical analyses; M.-C.F.: Methodology; D.S. (Dany Severac): Methodology, rearing; C.G.: Rearing and sample collection, Conceptualization, Formal analysis, review and editing; S.A.: Conceptualization, Formal analysis, review and editing; M.M.: Conceptualization, review and editing; E.J.-J.: Conceptualization, Formal analysis, review and editing; D.S. (David Siaussat): Conceptualization, data curation, formal analysis, prepared graphs and statistical analyses, writing—review and editing. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by ANR program, PHEROTOX.

**Data Availability Statement:** Data available on request due to restrictions eg privacy or ethical.

**Acknowledgments:** The authors thank Cyril Le Corre for help with insect rearing, Lydie Garnier for insect rearing and moth treatments. This work was supported by French National Founding Agency grant ANR-12 ADAP-0012-01.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Veres, A.; Wyckhuys, K.A.G.; Kiss, J.; Tóth, F.; Burgio, G.; Pons, X.; Avilla, C.; Vidal, S.; Razinger, J.; Bazok, R.; et al. An update of the Worldwide Integrated Assessment (WIA) on systemic pesticides. Part 4: Alternatives in major cropping systems. *Environ. Sci. Pollut. Res. Int.* **2020**, *27*, 29867–29899, doi:10.1007/s11356-020-09279-x.

2. Jeschke, P.; Nauen, R.; Schindler, M.; Elbert, A. Overview of the status and global strategy for neonicotinoids. *J. Agric. Food Chem.* **2011**, *59*, 2897–2908, doi:10.1021/jf101303g.

3. Casida, J.E.; Durkin, K.A. Neuroactive Insecticides: Targets, Selectivity, Resistance, and Secondary Effects. *Annu. Rev. Entomol.* **2013**, *58*, 99–117, doi:10.1146/annurev-ento-120811-153645.

4. Matsuda, K.; Buckingham, S.D.; Kleier, D.; Rauh, J.J.; Grauso, M.; Sattelle, D.B. Neonicotinoids: Insecticides acting on insect nicotinic acetylcholine receptors. *Trends Pharmacol. Sci.* **2001**, *22*, 573–580, doi:10.1016/s0165-6147(00)01820-4.

5. Goulson, D. REVIEW: An overview of the environmental risks posed by neonicotinoid insecticides. *J. Appl. Ecol.* **2013**, *50*, 977–987, doi:10.1111/1365-2664.12111@10.1111/.

6. Blacquière, T.; Smagghe, G.; Van Gestel, C.A.M.; Mommaerts, V. Neonicotinoids in bees: A review on concentrations, side-effects and risk assessment. *Ecotoxicology* **2012**, *21*, 973–992, doi:10.1007/s10646-012-0863-x.

7. Desneux, N.; Decourtye, A.; Delpuech, J.-M. The Sublethal Effects of Pesticides on Beneficial Arthropods. *Annu. Rev. Entomol.* **2007**, *52*, 81–106, doi:10.1146/annurev.ento.52.110405.091440.

8. Lürling, M.; Scheffer, M. Info-disruption: Pollution and the transfer of chemical information between organisms. *Trends Ecol. Evol.* **2007**, *22*, 374–379, doi:10.1016/j.tree.2007.04.002.

9. Calabrese, E.J.; Baldwin, L.A. Defining hormesis. *Hum. Exp. Toxicol.* **2002**, *21*, 91–97., doi:10.1191/0960327102ht217oa.

10. Calabrese, E.J. The Frequency of U-Shaped Dose Responses in the Toxicological Literature. *Toxicol. Sci.* **2001**, *62*, 330–338, doi:10.1093/toxsci/62.2.330.

11. Calabrese, E.J. The road to linearity: Why linearity at low doses became the basis for carcinogen risk assessment. *Arch. Toxicol.* **2009**, doi:10.1007/s00204-009-0412-4.

12. Cutler, G.C. Insects, Insecticides and Hormesis: Evidence and Considerations for Study. *Dose-Response* **2012**, *11*, 154–177, doi:10.2203/dose-response.12-008.Cutler.

13. Ding, J.; Zhao, Y.; Zhang, Z.; Xu, C.; Mu, W. Sublethal and Hormesis Effects of Clothianidin on the Black Cutworm (Lepidoptera: Noctuidae). *J. Econ. Entomol.* **2018**, *111*, 2809–2816, doi:10.1093/jee/toy254.

14. Guedes, R.N.C.; Cutler, G.C. Insecticide-induced hormesis and arthropod pest management. *Pest Manag. Sci.* **2014**, *70*, 690–697, doi:10.1002/ps.3669.

15. Dewer, Y.; Pottier, M.A.; Lalouette, L.; Maria, A.; Dacher, M.; Belzunces, L.P.; Kairo, G.; Renault, D.; Maibeche, M.; Siaussat, D. Behavioral and metabolic effects of sublethal doses of two insecticides, chlorpyrifos and methomyl, in the Egyptian cotton leafworm, Spodoptera littoralis (Boisduval) (Lepidoptera: Noctuidae). *Environ. Sci. Pollut. Res.* **2016**, *23*, 3086–3096, doi:10.1007/s11356-015-5710-1.

16. Lalouette, L.; Pottier, M.A.; Wycke, M.A.; Boitard, C.; Bozzolan, F.; Maria, A.; Demondion, E.; Chertemps, T.; Lucas, P.; Renault, D.; et al. Unexpected effects of sublethal doses of insecticide on the peripheral olfactory response and sexual behavior in a pest insect. *Environ. Sci. Pollut. Res.* **2016**, *23*, 3073–3085, doi:10.1007/s11356-015-5923-3.
17. Rabhi, K.K.; Esancy, K.; Voisin, A.; Crespin, L.; Le Corre, J.; Tricoire-Leignel, H.; Anton, S.; Gadenne, C. Unexpected effects of low doses of a neonicotinoid insecticide on behavioral responses to sex pheromone in a pest insect. *PLoS ONE*. **2014**, *9*, doi:10.1371/journal.pone.0114411.

18. Tricoire-Leignel, H.; Thany, S.H.; Gadenne, C.; Anton, S. Pest insect olfaction in an insecticide-contaminated environment: Info-disruption or hormesis effect. *Front Physiol*. **2012**, *3*, 58.

19. Enders, L.S.; Rault, L.C.; Heng-moss, T.M.; Siegfried, B.D.; Miller, J. Transcriptional responses of soybean aphids to sublethal insecticide exposure. *Insect Biochem. Mol. Biol*. **2019**, *103285*, doi:10.1016/j.ibmb.2019.103285.

20. Dondero, F.; Negri, A.; Boatti, L.; Marsano, F.; Mignone, F.; Viarengo, A. Transcriptomic and proteomic effects of a neonicotinoid insecticide mixture in the marine mussel (Mytilus galloprovincialis, Lam.). *Sci. Total Environ*. **2010**, *408*, 3775-3786, doi:10.1016/j.scitotenv.2010.03.040.

21. Diesner, M.; Gallot, A.; Binz, H.; Gaertner, C.; Vitecek, S.; Kahnt, J.; Schachtner, J.; Jacquin-Joly, E.; Gadenne, C. Mating-Induced Differential Peptidomics of Neuropeptides and Protein Hormones in Agrotis ipsilon Moths. *J. Proteome Res*. **2018**, *17*, 1397-1414, doi:10.1021/acs.jproteome.7b00779.

22. Gadenne, C.; Barrozo, R.B.; Anton, S. Plasticity in Insect Olfaction: To Smell or Not to Smell? *Annu. Rev. Entomol*. **2016**, *61*, 317-333, doi:10.1146/annurev-ento-010715-023523.

23. Rabhi, K.K.; Deisig, N.; Demondion, E.; Le Corre, J.; Robert, G.; Tricoire-Leignel, H.; Lucas, P.; Gadenne, C.; Anton, S. Low doses of a neonicotinoid insecticide modify pheromone response thresholds of central but not peripheral olfactory neurons in a pest insect. *Proc. R. Soc. B Biol. Sci*. **2016**, *283*, doi:10.1098/rspb.2015.2987.

24. Arthidoro de Castro, M.B.; Martinez, L.C.; Cossolin, J.F.S.; Serra, R.S.; Serrão, J.E. Cytotoxic effects on the midgut, hypopharyngeal, glands and brain of Apis mellifera honey bee workers exposed to chronic concentrations of lambda-cyhalothrin. *Chemosphere* **2020**, *248*, 126075, doi:10.1016/j.chemosphere.2020.126075. Epub 2020 Jan 30.

25. Moreira, D.R.; Sinópolis Gigliolli, A.A.; Falco, J.R.P.; Julio, A.H.F.; Volnistem, E.A.; Chagas, F.D.; Toledo, V.A.A.; Ruvolo-Takasusuki, M.C.C. Toxicity and effects of the neonicotinoid thiamethoxam on Scaptotrigona bipunctata Lepeletier, 1836 (Hymenoptera: Apidae). *Environ. Toxicol*. **2018**, *33*, 463-475, doi:10.1002/tox.22533.

26. Surendra Nath, B.; Surendra Kumar, R.P. Toxic impact of organophosphorus insecticides on acetylcholinesterase activity in the silkworm, Bombyx mori L. *Ecotoxicol. Environ. Saf.* **1999**, *42*, 157-162, doi:10.1006/eesa.1998.1735.

27. Gadenne, C.; Renou, M.; Sreng, L. Hormonal control of pheromone responsiveness in the male black cutworm Agrotis ipsilon. *Experientia* **1993**, *49*, 721-724, doi:10.1007/BF01923960.

28. Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. *EMBnet J.* **2011**, *17*, 10, doi:10.14806/ej.17.1.200.

29. Chevreux, B.; Suhai, S. Genome Sequence Assembly Using Trace Signals and Additional Sequence Information. *Comput. Sci. Biol. Proc. Ger. Conf. Bioinf.* **1999**, *99*, 45-56.

30. Fu, L.; Niu, B.; Zhu, Z.; Wu, S.; Li, W. CD-HIT: Accelerated for clustering the next-generation sequencing data. *Bioinformatics* **2012**, *28*, 3150-3152, doi:10.1093/bioinformatics/bts565.

31. Li, W.; Godzik, A. Cd-hit: A fast program for clustering and comparing large sets of protein or nucleotide sequences. *Bioinformatics* **2006**, *22*, 1658-1659.

32. Simão, F.A.; Waterhouse, R.M.; Ioannidis, P.; Kriventseva, E.V.; Zdobnov, E.M. BUSCO: Assessing genome assembly and annotation completeness with single-copy orthologs. *Bioinformatics* **2015**, *31*, 3210-3212, doi:10.1093/bioinformatics/btv351.

33. Eddy, S.R. A Probabilistic Model of Local Sequence Alignment That Simplifies Statistical Significance Estimation. *PLoS Comput. Biol.* **2008**, *4*, e1000069, doi:10.1371/journal.pcbi.1000069.

34. Petersen, T.N.; Brunak, S.; von Heijne, G.; Nielsen, H. SignalP 4.0: Discriminating signal peptides from transmembrane regions. *Nat. Methods* **2011**, *8*, 785-786.

35. Krogh, A.; Larsson, B.; Fau-von Heijne, G.; von Heijne, G.; Fau-Sonnhammer, E.L.; Sonnhammer, E.L.; Larsson, B.; von Heijne, G.; Sonnhammer, E.L. Predicting transmembrane protein topology with a hidden Markov model: Application to complete genomes. *J. Mol. Biol.* **2001**, *305*, 567-580, doi:10.1006/jmbi.2000.4315.

36. Langmead, B.; Trapnell, C.; Pop, M.; Salzberg, S.L. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. *Genome Biol.* **2009**, *10*, R25, doi:10.1186/gb-2009-10-3-r25.

37. Li, B.; Dewey, C.N. RSEM: Accurate transcript quantification from RNA-Seq data with or without a reference genome. *BMC Bioinform.* **2011**, *12*, 323, doi:10.1186/1471-2105-12-323.

38. Robinson, M.D.; Oshlack, A. A scaling normalization method for differential expression analysis of RNA-seq data. *Genome Biol.* **2009**, *10*, R25, doi:10.1186/gb-2010-11-3-r25.

39. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. *Bioinformatics* **2010**, *26*, 139-140, doi:10.1093/bioinformatics/btp616.

40. Bigot, L.; Shaik, H.A.; Bozzolan, F.; Party, V.; Lucas, P.; Debernard, S.; Siaussat, D. Peripheral regulation by ecdysteroids of olfactory responsiveness in male Egyptian cotton leaf worms, Spodoptera littoralis. *Insect Biochem. Mol. Biol.* **2012**, *42*, 22-31, doi:10.1016/j.ibmb.2011.10.003. Epub 2011 Oct 25.

41. Dépagne, J.; Chevalier, F. Technical updates to basic proteins focalization using IPG strips. *Proteome Sci.* **2012**, *10*, 54, doi:10.1186/1477-5956-10-54.
42. Ruiz-Delgado, G.J.; Mancías-Guerra, C.; Tamez-Gómez, E.L.; Rodríguez-Romo, L.N.; López-Otero, A.; Hernández-Arizpe, A.; Gómez-Almaguer, D.; Ruiz-Argüelles, G.J. Dimethyl Sulfoxide-Induced Toxicity in Cord Blood Stem Cell Transplantation: Report of Three Cases and Review of the Literature. *Acta Haematol.* **2009**, *122*, 1–5, doi:10.1159/00022767.

43. Harbo, J.R. Sterility in honey bees caused by dimethyl sulfoxide. *J. Hered.* **1986**, *77*, 129–130, doi:10.1093/oxfordjournals.jhered.a110187.

44. Milchreit, K.; Ruhnke, H.; Wegener, J.; Bienefeld, K. Effects of an insect growth regulator and a solvent on honeybee (*Apis mellifera* L.) brood development and queen viability. *Ecotoxicology* **2016**, *25*, 530–537, doi:10.1007/s10646-016-1611-4.

45. Jacob, S.W.; de la Torre, J.C. Pharmacology of dimethyl sulfoxide in cardiac and CNS damage. *Pharmacol. Rep.* **2009**, doi:10.1016/S1734-1140(09)70026-X.

46. Drummond, J.; Williamson, S.M.; Fitchett, A.E.; Wright, G.A.; Judge, S.J. Spontaneous honeybee behaviour is altered by persistent organic pollutants. *Ecotoxicology* **2017**, *26*, 141–150, doi:10.1007/s10646-016-1749-0.

47. Haider, S.; Pal, R. Integrated Analysis of Transcriptomic and Proteomic Data. *Curr. Genom.* **2013**, *14*, 91–110, doi:10.2174/1389202911314020003.

48. Kumar, D.; Bansal, G.; Narang, A.; Basak, T.; Abbas, T.; Dash, D. Integrating transcriptome and proteome profiling: Strategies and applications. *Proteomics* **2016**, *16*, 2533–2544, doi:10.1002/pmic.201600140.

49. Li, Z.; Yang, D.; Ding, J.; Zhang, J.; Pan, L.; Wang, X.; Yin, Q.; Gao, Y.; Chen, Z.; Chen, D. Resistance of mosquitoes to insecticide. Chinese. *J. Hyg. Insectic. Equip.* **2007**, *13*, 256–258.

50. Casida, J.E. Neonicotinoid metabolism: Compounds, substituents, pathways, enzymes, organisms, and relevance. *J. Agric. Food Chem.* **2011**, *59*, 2923–2931.

51. Iwasa, T.; Motoyama, N.; Ambrose, J.T.; Roe, R.M. Mechanism for the differential toxicity of neonicotinoid insecticides in the honey bee, *Apis mellifera*. *Crop Prot.* **2004**, *23*, 371–378.

52. Karunker, I.; Benting, J.; Lueke, B.; Ponge, T.; Nauen, R.; Roditakis, E.; Vontas, J.; Gorman, K.; Denholm, I.; Morin, S. Over-expression of cytochrome P450 CYP6CM1 is associated with high resistance to imidacloprid in the B and Q biotypes of Bemisia tabaci (Hemiptera: Aleyrodidae). *Insect Biochem. Mol. Biol.* **2008**, *38*, 634–644.

53. Daborn, P.; Boundy, S.; Yen, J.; Pittendrigh, B. DDT resistance in Drosophila correlates with Cyp6g1 over-expression and confers cross-resistance to the neonicotinoid imidacloprid. *Mol. Genet. Genom.* **2001**, *266*, 556–563.

54. Magesh, V.; Zhu, Z.; Tang, T.; Chen, S.; Li, L.; Wang, L.; Varma, K.K.; Wu, Y. Toxicity of Neonicotinoids to honey bees and detoxification mechanism in honey bees. *IOSR J. Environ. Sci. Toxicol. Food Technol.* **2017**, *11*, 102–110.

55. Yang, N.; Xie, W.; Yang, X.; Wang, S.; Wu, Q.; Li, R.; Pan, H.; Liu, B.; Shi, X.; Fang, Y. Transcriptomic and proteomic responses of sweetpotato whitefly, Bemisia tabaci, to thiamethoxam. *PLoS ONE* **2013**, *8*, e61820.

56. Chen, C.; Shi, X.; Desneux, N.; Han, P.; Gao, X. Detection of insecticide resistance in Bradysia odoriphaga Yang et Zhang (Diptera: Sciaridae) in China. *Ecotoxicology* **2017**, *26*, 868–875.

57. Alyokhin, A.; Baker, M.; Mota-Sanchez, D.; Dively, G.; Grafius, E. Colorado potato beetle resistance to insecticides. *Am. J. Potato Res.* **2008**, *85*, 395–413.

58. Yang, Y.X.; Yu, N.; Zhang, J.H.; Zhang, Y.X.; Liu, Z.W. Induction of P450 genes in Nilaparvata lugens and Sogatella furcifera by two neonicotinoid insecticides. *Insect Sci.* **2017**, *25*, 401–408, doi:10.1111/1744-7917.12440. Epub 2017 Mar 20.

59. Penning, T.M. The aldo-keto reductases (AKRs): Overview. *Chem. Biol. Interact.* **2015**, *234*, 236–246.

60. Chen, J.; Webb, T.J.; Powls, R.; Rees, H.H. Purification and Characterisation of Haemolymph 3-Dehydroecdysone 3β-Reductase in Relation to Ecdysteroid Biosynthesis in the Cotton Leafworm Spodoptera littoralis. *Eur. J. Biochem.* **1996**, *242*, 394–401.

61. Yamamoto, K.; Wilson, D.K. Identification, characterization, and crystal structure of an aldo-keto reductase (AKR2E4) from the silkworm Bombyx mori. *Arch. Biochem. Biophys.* **2013**, *538*, 156–163.

62. Yamamoto, K.; Yamada, N. Identification of a diazinon-metabolizing glutathione S-transferase in the silkworm, Bombyx mori. *Sci. Rep.* **2016**, *6*, 30073.

63. Marchitti, S.A.; Brocker, C.; Stagos, D.; Vasiliou, V. Non-P450 aldehyde oxidizing enzymes: The aldehyde dehydrogenase superfamily. *Expert Opin. Drug Metab. Toxicol.* **2008**, *4*, 697–720.

64. Fry, J.D.; Saweikis, M. Aldehyde dehydrogenase is essential for both adult and larval ethanol resistance in Drosophila melanogaster. *Genet. Res. (Camb)* **2006**, *87*, 87–92.

65. Lumjuan, N.; Wicheer, J.; Leelapat, P.; Choochote, W.; Somboon, P. Identification and characterisation of Aedes aegypti aldehyde dehydrogenases involved in pyrethroid metabolism. *PLoS ONE* **2014**, *9*, e102746.

66. Doorn, J.A.; Florang, V.R.; Schamp, J.H.; Vanle, B.C. Aldehyde dehydrogenase inhibition generates a reactive dopamine metabolite autotoxic to dopamine neurons. *Parkinsonism Relat. Disord.* **2014**, *20*, S73–S75.

67. Fitzmaurice, A.G.; Rhodes, S.L.; Cockburn, M.; Ritz, B.; Bronstein, J.M. Aldehyde dehydrogenase variation enhances effect of pesticides associated with Parkinson disease. *Neurology* **2014**, *82*, 419–426.

68. Volders, K.; Scholz, S.; Slabbaert, J.R.; Nagel, A.C.; Verstreken, P.; Creemers, J.W.M.; Callaerts, P.; Schwärzel, M. Drosophila rugose is a functional homolog of mammalian Neurobeachin and affects synaptic architecture, brain morphology, and associative learning. *J. Neurosci.* **2012**, *32*, 15193–15204.

69. Akbergenova, Y.; Cunningham, K.L.; Zhang, Y.V.; Weiss, S.; Littleton, J.T. Characterization of developmental and molecular factors underlying release heterogeneity at Drosophila synapses. *Elife* **2018**, *7*, e38268.

70. Faurobert, E.; Chen, C.-K.; Hurley, J.B.; Teng, D.H.-F. Drosophila neurocalcin, a fatty acylated, Ca-binding protein that associates with membranes and inhibits in vitro phosphorylation of bovine rhodopsin. *J. Biol. Chem.* **1996**, *271*, 10256–10262.
71. Gokhale, A.; Mullin, A.P.; Zlatic, S.A.; Easley, C.A.; Merritt, M.E.; Raj, N.; Larimore, J.; Gordon, D.E.; Peden, A.A.; Sanyal, S.; et al. The N-ethylmaleimide-sensitive factor and dysbindin interact to modulate synaptic plasticity. *J. Neurosci.* **2015**, *35*, 7643–7653, doi:10.1523/JNEUROSCI.4724-14.2015.

72. Li, T.; Tian, Y.; Li, Q.; Chen, H.; Lv, H.; Xie, W.; Han, J. The neurexin/N-ethylmaleimide-sensitive factor (NSF) interaction regulates short term synaptic depression. *J. Biol. Chem.* **2015**, *290*, 17656–17667, doi:10.1074/jbc.M115.644583.

73. Nunes, P.; Haines, N.; Kuppuswamy, V.; Fleet, D.J.; Stewart, B.A. Synaptic vesicle mobility and presynaptic F-actin are disrupted in a N-ethylmaleimide-sensitive factor allele of Drosophila. *Mol. Biol. Cell* **2006**, *17*, 4709–4719, doi:10.1091/mbc.E06-03-0253.

74. Caballero, J.P.; Murillo, L.; List, O.; Bastiat, G.; Flochlay-Sigognault, A.; Guerino, F.; Lefrançois, C.; Lautram, N.; Lapied, B.; Apaire-Marchais, V. Nanoencapsulated deltamethrin as synergistic agent potentiates insecticide effect of indoxacarb through an unusual neuronal calcium-dependent mechanism. *Pestic. Biochem. Physiol.* **2019**, *157*, 1–12.

75. Roedding, A.S.; Tong, S.Y.; Au-Yeung, W.; Li, P.P.; Warsh, J.J. Chronic oxidative stress modulates TRPC3 and TRPM2 channel expression and function in rat primary cortical neurons: Relevance to the pathophysiology of bipolar disorder. *Brain Res.* **2013**, *1517*, 16–27.

76. Ffrench-Constant, R.H.; Rocheleau, T.A. Drosophila γ-aminobutyric acid receptor gene Rdl shows extensive alternative splicing. *J. Neurochem.* **1993**, *60*, 2323–2326.

77. Jonsson, N.N.; Klafke, G.; Corley, S.W.; Tidwell, J.; Berry, C.M.; Koh-Tan, H.H. Molecular biology of amitraz resistance in cattle ticks of the genus Rhipicephalus. *Front. Biosci.* **2018**, *23*, 796–810.

78. Wang, X.; Puinean, A.M.; Williamson, M.S.; Smelt, C.L.C.; Millar, N.S.; Wu, Y. Mutations on M3 helix of Plutella xylostella glutamate-gated chloride channel confer unequal resistance to abamectin by two different mechanisms. *Insect Biochem. Mol. Biol.* **2017**, *86*, 50–57.

79. Karlstrom, R.O.; Wilder, L.P.; Bastiani, M.J. Lachesin: An immunoglobulin superfamily protein whose expression correlates with neurogenesis in grasshopper embryos. *Development* **1993**, *118*, 509–522.

80. Abe, T.; Yamazaki, D.; Murakami, S.; Hiroi, M.; Nitta, Y.; Maeyama, Y.; Tabata, T. The NAV2 homolog Sickie regulates F-actin-mediated axonal growth in Drosophila mushroom body neurons via the non-canonical Rac-Cofilin pathway. *Development* **2014**, *141*, 4716–4728.

81. Jones, M.A.; Amr, S.; Ferebee, A.; Huynh, P.; Rosenfeld, J.A.; Miles, M.F.; Davies, A.G.; Korey, C.A.; Warrick, J.M.; Shiang, R.; et al. Genetic studies in Drosophila and humans support a model for the concerted function of CISD2, PPT1 and CLN3 in disease. *Biol. Open* **2014**, *3*, 342–352, doi:10.1242/bio.20147559.

82. Koster, K.P.; Yoshii, A. Depalmitoylation by palmitoyl-protein thioesterase 1 in neuronal health and degeneration. *Front. Synaptic Neurosci.* **2019**, doi:10.3389/fnsyn.2019.00025.

83. Bannan, B.A.; Van Etten, J.; Kohler, J.A.; Tsoi, Y.; Hansen, N.M.; Sigmon, S.; Fowler, E.; Buff, H.; Williams, T.S.; Ault, J.G.; et al. The drosophila protein palmitoylome characterizing palmitoyl-thioesterases and DHHC palmitoyl-transferases. *Fly (Austin)* **2008**, *2*, 198–214, doi:10.4161/fly.6621.

84. Aby, E.; Gumps, K.; Roth, A.; Sigmon, S.; Jenkins, S.E.; Kim, J.J.; Kramer, N.J.; Parfitt, K.D.; Korey, C.A. Mutations in palmitoyl-protein thioesterase 1 alter exocytosis and endocytosis at synapses in Drosophila larvae. *Fly (Austin)* **2013**, *7*, 267–279.

85. Chu-LaGraff, Q.; Blanchette, C.; O'Hern, P.; Denefrio, C. The Batten disease Palmitoyl Protein Thioesterase 1 gene regulates neural specification and axon connectivity during Drosophila embryonic development. *PLoS ONE* **2010**, *5*, e14402.

86. Hickey, A.J.; Chotkowski, H.L.; Singh, N.; Ault, J.G.; Korey, C.A.; MacDonald, M.E.; Glaser, R.L. Palmitoyl-protein thioesterase 1 deficiency in Drosophila melanogaster causes accumulation of abnormal storage material and reduced life span. *Genetics* **2006**, *172*, 2379–2390.

87. Ding, B.J.; Löfstedt, C. Analysis of the Agrotis segetum pheromone gland transcriptome in the light of Sex pheromone biosynthesis. *BMC Genom.* **2015**, *16*, doi:10.1186/s12864-015-1909-2.

88. Abrieux, A.; Mhamdi, A.; Rabhi, K.K.; Egon, J.; Debernard, S.; Duportets, L.; Tricoire-Leignel, H.; Anton, S.; Gadenne, C. An Insecticide Further Enhances Experience-Dependent Increased Behavioural Responses to Sex Pheromone in a Pest Insect. *PLoS ONE* **2016**, *11*, e0167469, doi:10.1371/journal.pone.0167469.

89. Pan, Y.; Tian, F.; Wei, X.; Wu, Y.; Gao, X.; Xi, J.; Shang, Q. Thiamethoxam resistance in Aphis gossypii glover relies on multiple UDP-glucuronosyltransferases. *Front. Physiol.* **2018**, *9*, doi:10.3389/fphys.2018.00322.

90. Pan, J.; Yu, L.; Liu, D.; Hu, D. Synthesis and Insecticidal Activity of Mesoionic Pyrido [1, 2-α] pyrimidinone Derivatives Containing a Neonicotinoid Moiety. *Molecules* **2018**, *23*, 1217.

91. Meng, J.; Zhang, C.; Chen, X.; Cao, Y.; Shang, S. Differential protein expression in the susceptible and resistant Myzus persicae (Sulzer) to imidacloprid. *Pestic. Biochem. Physiol.* **2014**, *115*, 1–8.

92. Rand, E.E. du, Smit, S.; Beukes, M.; Apostolides, Z.; Pirk, C.W.W.; Nicolson, S.W. Detoxification mechanisms of honey bees (Apis mellifera) resulting in tolerance of dietary nicotine. *Sci. Rep.* **2015**, *5*, 11779.