Svoboda | Graniru | BBC Russia | Golosameriki | Facebook
Skip to main content
Advertisement
  • Loading metrics

Nicotinic acetylcholine receptor modulator insecticides act on diverse receptor subtypes with distinct subunit compositions

  • Wanjun Lu ,

    Contributed equally to this work with: Wanjun Lu, Zhihan Liu, Xinyu Fan, Xinzhong Zhang

    Roles Data curation

    Affiliation Ministry of Agriculture Key Laboratory of Molecular Biology of Crop Pathogens and Insects, Institute of Insect Sciences, Zhejiang University, Hangzhou, China

  • Zhihan Liu ,

    Contributed equally to this work with: Wanjun Lu, Zhihan Liu, Xinyu Fan, Xinzhong Zhang

    Roles Data curation

    Affiliation Ministry of Agriculture Key Laboratory of Molecular Biology of Crop Pathogens and Insects, Institute of Insect Sciences, Zhejiang University, Hangzhou, China

  • Xinyu Fan ,

    Contributed equally to this work with: Wanjun Lu, Zhihan Liu, Xinyu Fan, Xinzhong Zhang

    Roles Data curation

    Affiliation Ministry of Agriculture Key Laboratory of Molecular Biology of Crop Pathogens and Insects, Institute of Insect Sciences, Zhejiang University, Hangzhou, China

  • Xinzhong Zhang ,

    Contributed equally to this work with: Wanjun Lu, Zhihan Liu, Xinyu Fan, Xinzhong Zhang

    Roles Data curation, Formal analysis, Methodology

    Affiliation Tea Research Institute, Chinese Academy of Agricultural Sciences, Hangzhou, China

  • Xiaomu Qiao,

    Roles Formal analysis, Funding acquisition, Investigation, Methodology, Writing – original draft

    Affiliation Ministry of Agriculture Key Laboratory of Molecular Biology of Crop Pathogens and Insects, Institute of Insect Sciences, Zhejiang University, Hangzhou, China

  • Jia Huang

    Roles Conceptualization, Formal analysis, Funding acquisition, Investigation, Project administration, Supervision, Writing – original draft, Writing – review & editing

    [email protected]

    Affiliation Ministry of Agriculture Key Laboratory of Molecular Biology of Crop Pathogens and Insects, Institute of Insect Sciences, Zhejiang University, Hangzhou, China

Abstract

Insect nicotinic acetylcholine receptors (nAChRs) are pentameric ligand-gated ion channels mainly expressed in the central nervous system of insects. They are the directed targets of many insecticides, including neonicotinoids, which are the most widely used insecticides in the world. However, the development of resistance in pests and the negative impacts on bee pollinators affect the application of insecticides and have created a demand for alternatives. Thus, it is very important to understand the mode of action of these insecticides, which is not fully understood at the molecular level. In this study, we systematically examined the susceptibility of ten Drosophila melanogaster nAChR subunit mutants to eleven insecticides acting on nAChRs. Our results showed that there are several subtypes of nAChRs with distinct subunit compositions that are responsible for the toxicity of different insecticides. At least three of them are the major molecular targets of seven structurally similar neonicotinoids in vivo. Moreover, spinosyns may act exclusively on the α6 homomeric pentamers but not any other nAChRs. Behavioral assays using thermogenetic tools further confirmed the bioassay results and supported the idea that receptor activation rather than inhibition leads to the insecticidal effects of neonicotinoids. The present findings reveal native nAChR subunit interactions with various insecticides and have important implications for the management of resistance and the development of novel insecticides targeting these important ion channels.

Author summary

Neonicotinoids and spinosyns account for approximately 24% and 3% of the world market value of insecticides, respectively. However, the negative effects of neonicotinoids on pollinators have led to the development of novel insecticides, such as sulfoxaflor, flupyradifurone and triflumezopyrim. Although all act via insect nicotinic acetylcholine receptors, their modes of action are not fully understood. Our work shows that these insecticides act on diverse receptor subtypes with distinct subunit compositions. This finding could lead to the development of more selective insecticides to control pests with minimal effects on beneficial insects.

Introduction

Chemical insecticides have been widely used to control pests in the agriculture, horticulture, and forestry industries as well as homes and cities. They have also played a vital role in preventing the spread of human and animal vector-borne diseases. However, insecticide resistance is a serious worldwide problem for invertebrate pest control, and more than 600 different insect and mite species have become resistant to at least one insecticide. In addition, at least one case of resistance to more than 335 insecticides/acaricides has been documented [1]. Therefore, there is great demand for effective insecticide resistance management (IRM) and the development of new pest control compounds. To address both issues, we need to determine the mode of action of insecticides, which is the molecular-level processes underlying the effects of insecticides [2].

A complete understanding of the mode of action of an insecticide requires knowledge of how it affects a specific target site within an organism. Although most insecticides have multiple biological effects, toxicity is usually attributed to a single major effect. For some insecticides, however, the exact molecular targets remain elusive. To ascribe whether a candidate protein is indeed the target for an insecticidal effect in vivo, it is not sufficient to demonstrate an in vitro biochemical interaction between an insecticide and a protein. Genetic evidence demonstrating an effect due to mutation of the candidate target must be obtained before a given protein can be identified as an insecticide target.

Neonicotinoids (acetamiprid, clothianidin, dinotefuran, imidacloprid, nitenpyram, thiacloprid, and thiamethoxam) are remarkably effective at controlling agricultural pests, ectoparasites and arthropod vectors [3]. They are taken up by the roots or leaves and translocated to all parts of the plant due to high systemic activity, making them effectively toxic to a wide range of sap-feeding and foliar-feeding insects. Thus, neonicotinoids account for 24% of the global insecticide market, which is the largest market share of all chemical classes [1]. They act selectively on insect nicotinic acetylcholine receptors (nAChRs) as agonists compared with the mammalian-selective nicotine. Spinosyns are a naturally derived, unique family of macrocyclic lactones that act on insect nAChR in an allosteric fashion. In addition, sulfoximine sulfoxaflor, butenolide flupyradifurone and mesoionic triflumezopyrim are three newly developed insecticides that are also nAChR competitive modulators [4]. It is expected that the market of all of the above nAChRs targeting insecticides that show excellent insect-to-mammalian selectivity will continue to grow. However, the molecular targets of neonicotinoids and other nAChR modulators remain unclear, mainly because the structure and assembly of native nAChRs in insects have not been clarified [5].

Cation-selective nAChRs are members of the Cys-loop ligand-gated ion channel superfamily responsible for rapid excitatory neurotransmission. The functional nAChRs are homo- or heteromeric pentamers of structurally related subunits arranged around a central ion-conducting pore [6]. Each subunit has an extracellular N-terminal domain that contains six distinct regions (loops A–F) involved in ligand binding, four C-terminal transmembrane segments (TM1–TM4) and an intracellular loop between TM3 and TM4. nAChRs are divided into α-subunits possessing two adjacent cystine residues in loop C, while those subunits without this motif are termed nonα subunits. In vertebrates, 17 nAChR subunits have been identified, and they can coassemble to generate a diverse family of nAChR subtypes with different pharmacological properties and physiological functions. Insects have fewer nAChR subunits (10–12 subunits) according to the available genome data. Although coimmunoprecipitation studies have indicated potential associations of several subunits, the exact subunit composition of native insect nAChRs remains unknown [5]. Unlike the vertebrate counterparts, heterologous expression of genuine arthropod α and β subunits was not successful until two groups recently found that three ancillary proteins are essential for robust expression of arthropod nAChR heteromers [7, 8]. Thus, for a long time, researchers have used hybrid receptors with insect α subunits and mammalian/avian β subunits to study the interaction of insecticides and receptors. However, such alternatives may not faithfully reflect all features of native nAChRs [9].

In this study, we systematically examined the effects of eleven nAChR-targeting insecticides against ten (seven α and three β) Drosophila melanogaster subunit mutants. We found that there are multiple subtypes of receptors with distinct subunit compositions that are responsible for the toxicity of different insecticides. Artificial activation/inhibition of subunit-expressing neurons also mimicked insecticide poisoning symptoms in pests. Elucidating the molecular targets of these economically important agrochemicals and the assembly of native nAChRs will be very helpful for resistance management and ecotoxicological evaluations of beneficial insects, such as predators and pollinators.

Results

Generation of the nAChRβ1R81T mutant

We obtained all 10 nAChR knockout mutants from Yi Rao’s laboratory and found that the KO of α4 and β1 was homozygous lethal. Previous studies [10, 11] and our results (Tables A-K in S1 Text) both showed that target-site resistance in nAChRs is mostly recessive to semirecessive, with the wild-type allele being dominant. Therefore, we used point mutation alleles of α4 and β1 in all experiments instead of heterozygous KO flies. The α4T227M mutant (redeye, rye) is a dominant-negative allele that causes a reduced sleep phenotype in flies [12]. An R81T mutation of nAChR β1 was found in neonicotinoid-resistant peach aphids and later in cotton aphids [13, 14]; therefore, we introduced a homologous mutation into the β1 locus of Drosophila melanogaster with CRISPR–Cas9–mediated homology-directed repair (HDR). The design of the gRNA target site and HDR template was reported, and the screening of successful R81T knock-in was performed under imidacloprid selection pressure and confirmed by direct DNA sequencing (Fig 1 and Fig A in S1 Text).

thumbnail
Fig 1. Generation of the nAChRβ1R81T mutant by CRISPR/Cas9 genome editing.

(A) Schematic of the nAChRβ1 locus and the sequence of the donor construct. The boxes represent exons, and the coding regions are shown in blue. The gRNA sequence is indicated in red, and the codon for amino acid substitution (CGT to ACT) is highlighted in green. One synonymous mutation (G to A) is also introduced in the PAM region (in yellow) to prevent recleavage from Cas9 after successful integration. (B) Sequence comparison between wild-type flies and flies with point mutations. The nucleotides replaced are highlighted in green and yellow boxes.

https://doi.org/10.1371/journal.pgen.1009920.g001

nAChR mutants showed distinct resistance to multiple insecticides

We tested the effects of 10 nAChR mutants and some heterozygous mutants against 11 insecticides (Fig 2 and Tables A-K in S1 Text). The α1 mutant showed moderate levels of resistance to imidacloprid, thiacloprid, acetamiprid and triflumezopyrim, and the LC50 resistance ratio (RR) was approximately 13.5–88.0. Its heterozygous mutant also showed low levels of resistance to these insecticides. In addition, it showed low but statistically significant increases in RR (2.7–3.7) to thiamethoxam, clothianidine, dinotefuran and nitenpyram. The α2 mutant also showed similar levels of resistance (17.2–48.5 in terms of RR) to imidacloprid, thiacloprid and triflumezopyrim. The α3 mutant showed small RR increases (2.7–5.5) to thiamethoxam, clothianidine, dinotefuran, nitenpyram, sulfoxaflor and flupyradifurone. The α4, α5, α6, α7 and β3 mutants are sensitive to almost all the tested insecticides. The obvious exception is the α6 homozygous mutant, which is resistant to spinetoram and has a RR of 42.8, although the heterozygous mutant is close to the wild type (RR 1.2). The β1 mutant exhibited medium to high resistance to all insecticides (23.9–398.3 in terms of RR) except spinetoram, and its heterozygous mutant showed small RR increases for most insecticides. The resistance profile of the β2 mutant was similar to that of the α1 mutant, with a 13.0- to 84.3-fold increased RR to imidacloprid, thiacloprid, acetamiprid and triflumezopyrim.

thumbnail
Fig 2. Nonlinear log-dose mortality data for tested insecticides against ten Drosophila nAChR homozygous mutants, including eight null alleles and two point mutation alleles (α4T227M and β1R81T).

Mortality (0–1 means 0–100% in terms of percentage) of control and mutant female adults after 48 hours of exposure to increasing concentrations of insecticides. Error bars represent standard deviations.

https://doi.org/10.1371/journal.pgen.1009920.g002

Both the α1 and β1 mutants showed variable resistance to multiple insecticides; thus, we generated a α1/β1 double mutant with recombination. However, the eggs laid by this combined mutant could not hatch and thus were not used in further experiments. A recent paper also generated a β1 R81T Drosophila and found that it has serious defects in reproduction and locomotion [15]; however, the β1 mutant we made here did not show any significant fitness cost (Fig B in S1 Text). The sequences of α5, α6 and α7 are very close and show high similarity to the vertebrate nAChR α7 subunit; however, only the α6 mutant showed resistance to spinetoram. We further generated a α5/α7 double mutant that was still sensitive to spinetoram (Table K in S1 Text), indicating that the α6 homomeric channel could be the sole target for spinosyns.

Hyperactivating/Silencing nAChR-expressing neurons mimics insecticide poisoning symptoms

Insects present similar reactions upon exposure to neonicotinoids, sulfoxaflor, flupyradifurone and spinosyns. Early-onset behaviors include hyperactivity, convulsion, uncoordinated movements, leg extension and tremors. At higher doses, these excitatory symptoms can induce severe tremors and complete paralysis, which ultimately leads to death [1618]. We then wondered whether artificial activation of nAChR-expressing neurons would induce insecticide-like poisoning symptoms. Thus, we used the thermosensitive cation channel Drosophila TRPA1 to acutely hyperstimulate these neurons with all available nAChR KI-Gal4 strains [19]. We found that expressing trpA1 in nAChRα12A-GAL4, nAChRα22A-GAL4, nAChRα32A-GAL4, nAChRα62A-GAL4 and nAChRβ22A-GAL4 neurons strongly induced hyperactivity behavior at 32°C and eventually led to paralysis (Fig 3A and S1 Video), which is similar to the abovementioned symptoms. However, activation of nAChRβ32A-GAL4 neurons did not show any behavioral defects. These results parallel the above bioassay data showing that the deletion of α1, α2, α3, α6 and β2 caused medium to high resistance to these insecticides. Therefore, thermogenetic activation of some nAChR-expressing neurons in a short time window phenocopies the action of insecticides in target pests, which demonstrates that in vivo pharmacological activation of these subunit-containing nAChRs leads to toxicity and ultimately death.

thumbnail
Fig 3. Effects of artificial neuronal activation and inhibition in various nAChR-expressing neurons.

(A) Thermogenetic activation of five nAChR-expressing neurons using UAS-trpA1-induced paralysis behavior. (B) Thermogenetic silencing of five nAChR-expressing neurons using UAS-Shibirets decreased activity. n = 30–50.

https://doi.org/10.1371/journal.pgen.1009920.g003

The poisoning symptoms associated with triflumezopyrim are distinct from other insecticides that act on nAChRs since it inhibits rather than activates insect nAChRs. There are no neuroexcitatory symptoms after treatment with triflumezopyrim; in contrast, triflumezopyrim induces lethargic poisoning characterized by slow but coordinated leg movements, and insects become less responsive to stimuli over time [20]. Thus, we chose to use UAS-Shibirets to inhibit nAChR-expressing neurons [21]. As expected, nAChRα12A-GAL4, nAChRα22A-GAL4 and nAChRβ22A-GAL4 neurons produced “sluggish” behavior rather than hyperactivity (Fig 3B). The flies exhibited almost no translational or rotational body movement (S1 Video). Silencing of nAChRα32A-GAL4 and nAChRα62A-GAL4 neurons also produced similar behaviors, further confirming that the α3- and α6-containing nAChRs cannot be blocked by triflumezopyrim; otherwise, both mutants would show resistance in bioassays.

Expression patterns of nAChRs in KO mutants

We confirmed that the KO coding regions were not detected or barely detectable with real-time PCR quantification (Fig C in S1 Text). There was no large difference in the expression levels of each subunit in these mutant flies except for β3, which was relatively less transcribed than the other genes. For the α1 heterozygous mutant, the mRNA levels of all subunits were almost the same as those of the wild-type control.

Discussion

The Insecticide Resistance Action Committee (IRAC) classifies neonicotinoids, sulfoximines, butenolides and mesoionics according to their chemical similarity relations into subgroups 4A, 4C, 4D and 4E, respectively. However, our results showed that sulfoxaflor and flupyradifurone may mainly act on the same nAChR subtype, which consists of α3 and β1 subunits (Fig 4), although other subunits may also be involved considering genetic redundancy. More importantly, we found that neonicotinoids act on distinct nAChR subtypes and that such selectivity is not dependent on the aromatic heterocyclic (A) or the electron-withdrawing nitro or cyano moiety (X-Y), which is considered the key toxophore. Interestingly, the ring systems and the R2 substituents in the open-chain structures are the determining factors (Fig 4). For example, the α1, α2, β1 and β2 mutants showed similar levels of resistance to imidacloprid and thiacloprid (both have a five-membered ring), indicating that they mainly act on the same α1/α2/β1/β2 pentamer (Fig 4B). This finding is consistent with previous ex vivo recording results [22] and two recent reconstituted studies, which showed that both drugs act as partial agonists on the α1/α2/β1/β2 nAChR [7, 8]. Acetamiprid is structurally similar to thiacloprid with the cyanoimine pharmacophore, although the acyclic configuration changes its molecular target in vivo. It may act on the α1/β1/β2 nAChR, and again, electrophysiological studies have already indicated that acetamiprid is nearly a full agonist [22]. Moreover, its potency on the recombinant louse α1/α2/β1/β2 nAChR is approximately 860-fold lower than that of thiacloprid [8]. Although thiamethoxam has a six-membered ring, it is a prodrug without intrinsic nAChR activity until metabolized to the active form clothianidine in plants and insects [23]. Therefore, thiamethoxam, clothianidine, dinotefuran and nitenpyram can be considered the same type, which has N-methyl substitutions in the R2 position and mainly acts on the α1/α3/β1 nAChR (Fig 4B). Neonicotinoids are traditionally divided into nitroimines (NNO2), nitromethylenes (CHNO2) or cyanoimines (NCN), although our findings indicate that they should be classified according to their major nAChR subtype targets.

thumbnail
Fig 4. Resistance patterns of the tested insecticides on different nAChR mutants (A) and proposed target receptor subtypes for neonicotinoids and others (B).

Various resistance ratios are grouped and represented as different colors in the heatmap. Thiamethoxam is considered a prodrug of clothianidin and thus is not listed in the structural formula.

https://doi.org/10.1371/journal.pgen.1009920.g004

Despite the widespread use of neonicotinoids for almost three decades, the first and only field-evolved target-site resistance mutation (R81T in nAChRβ1) was reported in 2011, and it has only been found in two species to date [13, 14]. This unusual phenomenon is consistent with our findings that the seven neonicotinoids actually act on multiple receptor types in vivo and that only the β1 mutant caused high resistance to all neonicotinoids. New nicotine-mimicking insecticides, such as sulfoxaflor and flupyradifurone, mainly act on another nAChR subtype that is distinct from neonicotinoids (Fig 4), indicating their potential use in insecticide resistance management.

Electrophysiological studies with native tissues or recombinant receptors showed that low concentrations of neonicotinoids can block nAChR, while higher concentrations can activate the receptor [7,24]. Therefore, it is still unclear whether insecticidal activity is the consequence of nAChR inhibition or activation in vivo. We found that transient artificial activation rather than inhibition of nAChR-expressing neurons was sufficient to induce neonicotinoid-like poisoning symptoms in flies (Fig 3). Thus, the overall effect of neonicotinoids is neuronal depolarization by activation of nAChR, which is more physiologically relevant.

Triflumezopyrim is the first member of a new class of mesoionic insecticides that act via inhibition of the orthosteric binding site of the nAChR [20]. We found that the α1/α2/β1/β2 nAChR could be its major target, similar to imidacloprid and thiacloprid, and all these mutants showed high resistance to triflumezopyrim (Fig 4A). This finding is consistent with radioligand binding results showing that triflumezopyrim potently displaced [3H]imidacloprid with a Ki value of 43 nM based on membrane preparations from the aphid [20]. Thermogenetic inhibition of neurons expressing α1, α2 and β2 also mimicked lethargic intoxication symptoms (Fig 3B). Thus, to maintain the durability and effectiveness of this new powerful tool for the control of hopper species in rice, it is critical to avoid repeated use of triflumezopyrim with imidacloprid and thiacloprid.

Spinosyns, including spinosad and spinetoram, have been shown to act on a population of nAChRs that are not targeted by neonicotinoids, and the binding site is also distinct from the orthosteric site [4]. The α6 subunit has been proposed as the main target of spinosyns since the field-evolved resistance to spinosad is associated with loss-of-function mutations of α6 loci in many pest insects [2532]. However, the involvement of other subunits, such as α5 and α7, which are phylogenetically close to α6, has not been clarified (Fig D in S1 Text). Previous reports showed that α5 and α7 can form functional homomeric and heteromeric channels in vitro while α6 can only form heteromeric channels with α5 or α5/α7 together [33, 34]. We then wondered whether there was genetic redundancy among these evolutionarily conserved subunits. We found that the α5, α7 and α5/α7 double mutants were all sensitive to spinetoram (Table K in S1 Text), indicating that spinosyns may exclusively act on the α6 homomeric nAChR, which is consistent with a recent report using spinosad [35]. Thermogenetic activation of α6-expressing neurons also induced spinosyn-like poisoning symptoms in flies.

Our current knowledge about the subunit composition of insect nAChRs is very limited. Immunoprecipitation data with subunit-specific antibodies showed that Drosophila α3 and β1 coassemble within the same receptor complex [36]. Further studies from the same group indicated that α1/α2/β2 and β1/β2 may coassemble into the same receptor complex [37]. Similar studies using the brown planthopper suggested that there are two populations of nAChRs that contain Drosophila-equivalent subunit combinations α1/α2/β1 and α3/β1/β2 [38]. These previous findings are partially confirmed by the present results because α3/β1, α1/α3/β1, α1/β1/β2 and α1/α2/β1/β2 could be the major receptor subtypes for the tested insecticides, indicating that the β1 subunit could be an indispensable component for all heteromecic pentamers (Fig 4). In addition, we noticed that for some insecticides, different subunit mutations contribute in an asymmetrical manner to resistance (Fig 4A). Therefore, there could be functional redundancy between some α-type subunits, and we cannot exclude the existence of other potential receptor subtypes, such as α1/β1 and α3/β1/β2. The diversity of insect nAChRs and their druggability make them an extremely important target for insecticide development.

Growing evidence indicates that sublethal doses of neonicotinoids, such as imidacloprid, thiamethoxam and clothianidin negatively affect wild and managed bees, which are important pollinators in ecosystems and agriculture [3941]. They reduce reproduction and colony development, perhaps by impairing the foraging, homing and nursing behaviors of bees [42]. These severe sublethal effects have led to heavy restrictions on the use of the above three neonicotinoids in Europe to protect bee pollinators [43]. Sulfoxaflor and flupyradifurone are potential alternatives for neonicotinoids; however, their risk to bees is controversial [4446]. Therefore, it is critical to understand the mode of action of these insecticides inside bees. The core groups of nAChR subunits are highly conserved among different insects spanning ~300 million years of evolution [47], which is likely due to their essential roles in the nervous system. Most Drosophila nAChR subunit genes (except α5 and β3) have one-to-one orthologs in other insects, including honeybees and bumblebees (Fig D in S1 Text), and the sequence identities between orthologs are also high (Table L in S1 Text). Thus, the expression, assembly and function of these receptors could be conserved between flies and bees, suggesting that our results will enable further studies about the ecotoxicology and risk assessment of these nAChR modulators.

Materials and methods

Insecticides

The following were purchased commercially: imidacloprid (600 g/LSC, Bayer CropScience, Germany), thiamethoxam (70%GZ, Syngenta, China), clothianidin (48%SC, HeNan Hansi crop protection, China), dinotefuran (20%SG, Mitsui Chemicals, Japan), nitenpyram (30%WG, ZinGrow, China), acetamiprid (20%SP, Noposion, China), thiacloprid (40%SC, Limin Chemical, China), sulfoxaflor (22%SC, Dow AgroSciences, USA), flupyradifurone (17%SC, Bayer CropScience, Germany), triflumezopyrim (10%SC, DuPont, USA), spinetoram (60 g/LSC, Dow AgroSciences, USA) and Triton X-100 (Sangon Biotech, China).

Fly strains

Flies were maintained and reared on conventional cornmeal agar molasses medium at 25 ± 1°C and 60% ± 10% humidity with a photoperiod of 12 hours light:12 hours night. For experiments using UAS-trpA1 and UAS-Shibirets transgenes, flies were reared at 21°C. The following strains were sourced from the Bloomington Stock Center (Indiana University): vas-cas (#51323), UAS-trpA1 (#26263), and UAS-Shibirets (44222). All nAChR KO mutants and KI-Gal4 strains were gifts from Dr. Yi Rao (Deng et al., 2019) (Peking University). The w1118 strain used for outcrossing was used as a wild type for the insecticide bioassays.

We generated the nAChRβ1R81T mutant by CRISPR/Cas9 genome editing. The gRNA sequence (3 L:4433329~4433352, ATCAAACGTTTGGTTAACTTTAG) was designed with flyCRISPR Target Finder (https://flycrispr.org/target-finder/) and cloned into the pDCC6 plasmid (Addgene #59985). A 110 bp ssODN (single-strand oligodeoxynucleotide) was custom-synthesized as the donor template to replace the targeted genomic region. This ssODN contained three nucleotide changes, with two (CG to AC) conferring the R81T mutation and one synonymous mutation (G to A) to prevent recleavage from Cas9 after incorporation. Both the gRNA plasmid and ssODN were microinjected into the embryos of vas-cas flies (BL #51323). The crossing and selection scheme is shown in Fig A in S1 Text.

Insecticide bioassays

Three- to five-day-old and uniformly sized adult females were used in insecticide bioassays to assess the susceptibility of different fly strains. The testing method was modified from the IRAC susceptibility test method 026 (https://irac-online.org/methods/). Briefly, the required serial dilutions of insecticide solution were prepared in 200 g/L sucrose using formulated insecticides. Approximately 5 ml of insecticide solution is required for each concentration. A piece of dental wick (1 cm) was placed in a standard Drosophila vial and treated with 800 μL of 20% aqueous sucrose with or without insecticide. Ten flies of each genotype were transferred into vials, with 3–6 vials for each concentration, and each genotype was repeated at least 3 times for every tested insecticide. The vials were kept upside down until all flies became active to avoid flies becoming trapped in the dental wick. The bioassay was assessed after 48 h, and dead flies and seriously affected flies displaying no coordinated movement that were unable to walk up the vial or get to their feet were cumulatively scored as ‘affected’. The LC50 values were calculated by probit analysis using Polo Plus software (LeOra Software, Berkeley, CA, USA). Nonlinear log dose–response curves were generated in GraphPad Prism 8.21 (GraphPad Software Inc., La Jolla, CA, USA).

Thermogenetic activation and silencing assays

Flies for TRPA1-mediated thermogenetic activation and Shibire-mediated silencing experiments were collected upon eclosion and reared in vials containing standard food medium at 21°C for 5–8 days. For thermogenetic activation with the UAS-trpA1 transgene, 10 flies were transferred to new empty vials by gentle inspiration, and then the assays were performed at 23°C and 32°C for 10 minutes. Each genotype was repeated for at least 5 times. The percentage of paralysis behavior in which the animal lies on its back with little effective movement of the legs and wings was measured. For the silencing assays, the UAS-Shibirets transgene was used and flies were also transferred to fly vials at 23°C and 32°C for 10 minutes.

Real-time quantitative PCR

The relative transcription levels of nAChRs in different KO mutants were examined using real-time quantitative PCR performed with a CFX96TM Real-Time PCR System (Bio–Rad, Hercules, USA). Total RNA was isolated with TRIzol reagent according to the manufacturer’s instructions. Residual genomic DNA was removed by RQ1 RNase-Free DNase (Promega, Madison, USA). Total RNA was reverse transcribed to cDNA with EasyScript First-Strand cDNA Synthesis SuperMix (Transgene, Beijing, China). qPCR with gene-specific primers was performed with ChamQ Universal SYBR qPCR Master Mix (Vazyme, Nanjing, China) to investigate the relative expression levels of different nAChRs. RpL32 (ribosomal protein L32) was used as an internal control. The relative expression of nAChRs was normalized to the reference (RpL32) using the 2-ΔΔCT method. The primers used are listed in Table M in S1 Text.

Fecundity and development assays

Ten pairs of freshly emerged couples of wild-type control and β1R81T mutant were transferred into vials containing normal food for 72 hours. These files were then transferred into a new dish that was used for the egg-laying assay. The numbers of eggs laid in each dish were recorded after 24 hours. To calculate the larval to pupal rate, 60 second-instar larvae were collected and transferred into a new vial as one group. The numbers of pupae in each vial were recorded after 7 days in an incubator. Each genotype was repeated at least three times in duplicate.

Climbing assay

Approximately three-day-old male flies were collected with CO2 anesthesia into groups of 10 and then allowed to recover for 2 days. A climbing tube consisted of two vials with 90 mm height and 20 mm diameter. The flies were filmed for 30 seconds with a SONY HDR-CX900E camera. The climbing index (percentage of flies in the upper half of the vial) was determined at 5 second intervals after the flies had been tapped down to the bottom of the vials.

Phylogenetic analysis

The following representative species from a variety of orders were selected: Apis mellifera (honey bee), Tribolium castaneum (red flour beetle), Myzus persicae (green peach aphid), Bombyx mori (silk worm), Bombus terrestris (bumble bee) and Drosophila melanogaster (fruit fly). To identify the orthologs of the D. melanogaster nAChR subunits, we searched against NCBI non-redundant protein database using BLASTP. We renamed M. persicae and B. terrestris nAChR subunit proteins according to their closest orthologs. All the amino acid sequences were aligned by Clustal X. A neighbor-joining tree was performed by MEGA 11 with default parameters, 1000 bootstrap replications, and substitution with JTT model and visualized by Evolview (https://www.evolgenius.info//evolview/).

Supporting information

S1 Video. Effects of thermogenetic activation and inhibition in nAChRα1- expressing neurons.

The following transgenes were used: nAChRα12A-GAL4 > UAS-trpA1; nAChRα12A-GAL4 > UAS-Shibirets. Other nAChR KI-Gal4 strains like nAChRα22A-GAL4, nAChRα32A-GAL4, nAChRα62A-GAL4 and nAChRβ22A-GAL4 also produced similar behaviors when stimulated under 32°C.

https://doi.org/10.1371/journal.pgen.1009920.s001

(MP4)

S1 Text.

Fig A in S1 Text. The crossing schemes to establish the nAChRβ1R81T knock-in line. The HDR event was isolated by imidacloprid selection and confirmed by PCR. The vas-Cas9 (3XP3 RFP) was removed by the absence of red fluorescence in eyes. Fig B in S1 Text. Effects of nAChRβ1R81T point mutation on number of eggs laid (A), pupation rate of larvae (B) and negative geotaxis behavior (C). Fig C in S1 Text. Expression patterns of the nAChR genes in different KO mutants. Fig D in S1 Text. Phylogenetic relationships of core groups of nAChR subunits from 6 representative insect species including Apis mellifera (honey bee), Tribolium castaneum (red flour beetle), Myzus persicae (green peach aphid), Bombyx mori (silk worm), Bombus terrestris (bumble bee) and Drosophila melanogaster (fruit fly). The colorful dots at the nodes of the branches represent the values of bootstrap support for each branch. The D. melanogaster FMRFamide receptor (DmFR) was used as an outgroup. The sequence accession numbers are shown in Table N in S1 Text. Table A in S1 Text. Log dose probit mortality data and resistance ratios for imidacloprid. Table B in S1 Text. Log dose probit mortality data and resistance ratios for thiacloprid. Table C in S1 Text. Log dose probit mortality data and resistance ratios for acetamiprid. Table D in S1 Text. Log dose probit mortality data and resistance ratios for thiamethoxam. Table E in S1 Text. Log dose probit mortality data and resistance ratios for clothianidin. Table F in S1 Text. dose probit mortality data and resistance ratios for dinotefuran. Table G in S1 Text. Log dose probit mortality data and resistance ratios for nitenpyram. Table H in S1 Text. Log dose probit mortality data and resistance ratios for flupyradifurone. Table I in S1 Text. Log dose probit mortality data and resistance ratios for sulfoxaflor. Table J in S1 Text. Log dose probit mortality data and resistance ratios for triflumezopyrim. Table K in S1 Text. Log dose probit mortality data and resistance ratios for spinetoram. Table L in S1 Text. Sequence identities between Drosophila nAChR subunits and corresponding orthologs in other insects. Table M in S1 Text. Primers used in qPCR analysis. Table N in S1 Text. The accession numbers of sequences used in Fig D in S1 Text.

https://doi.org/10.1371/journal.pgen.1009920.s002

(DOCX)

Acknowledgments

We thank Yi Rao (Peking University) and the Bloomington Drosophila Stock Center for fly stocks.

References

  1. 1. Sparks TC, Crossthwaite AJ, Nauen R, Banba S, Cordova D, Earley F, et al. Insecticides, biologics and nematicides: Updates to IRAC’s mode of action classification—a tool for resistance management. Pestic Biochem Physiol. 2020;167:104587. Epub 2020/06/13. pmid:32527435.
  2. 2. Casida JE, Bryant RJ. The ABCs of pesticide toxicology: amounts, biology, and chemistry. Toxicol Res (Camb). 2017;6(6):755–63. pmid:30090540.
  3. 3. Jeschke P, Nauen R, Beck ME. Nicotinic Acetylcholine Receptor Agonists: A Milestone for Modern Crop Protection. Angew Chem Int Edit. 2013;52(36):9464–85. pmid:23934864.
  4. 4. Crossthwaite AJ, Bigot A, Camblin P, Goodchild J, Lind RJ, Slater R, et al. The invertebrate pharmacology of insecticides acting at nicotinic acetylcholine receptors. J Pestic Sci. 2017;42(3):67–83. Epub 2017/08/20. pmid:30363948.
  5. 5. Matsuda K, Ihara M, Sattelle DB. Neonicotinoid Insecticides: Molecular Targets, Resistance, and Toxicity. Annu Rev Pharmacol Toxicol. 2020;60:241–55. Epub 2020/01/10. pmid:31914891.
  6. 6. Casida JE. Neonicotinoids and Other Insect Nicotinic Receptor Competitive Modulators: Progress and Prospects. Annu Rev Entomol. 2018;63:125–44. Epub 2018/01/13. pmid:29324040.
  7. 7. Ihara M, Furutani S, Shigetou S, Shimada S, Niki K, Komori Y, et al. Cofactor-enabled functional expression of fruit fly, honeybee, and bumblebee nicotinic receptors reveals picomolar neonicotinoid actions. Proc Natl Acad Sci U S A. 2020;117(28):16283–91. Epub 2020/07/03. pmid:32611810.
  8. 8. Rufener L, Kaur K, Sarr A, Aaen SM, Horsberg TE. Nicotinic acetylcholine receptors: Ex-vivo expression of functional, non-hybrid, heteropentameric receptors from a marine arthropod, Lepeophtheirus salmonis. PLoS Pathog. 2020;16(7):e1008715. Epub 2020/07/28. pmid:32716968.
  9. 9. Ihara M, Buckingham SD, Matsuda K, Sattelle DB. Modes of Action, Resistance and Toxicity of Insecticides Targeting Nicotinic Acetylcholine Receptors. Curr Med Chem. 2017;24(27):2925–34. Epub 2017/02/09. pmid:28176635.
  10. 10. Zimmer CT, Garrood WT, Puinean AM, Eckel-Zimmer M, Williamson MS, Davies TG, et al. A CRISPR/Cas9 mediated point mutation in the alpha 6 subunit of the nicotinic acetylcholine receptor confers resistance to spinosad in Drosophila melanogaster. Insect Biochem Mol Biol. 2016;73:62–9. Epub 2016/04/28. pmid:27117524.
  11. 11. Mottet C, Fontaine S, Caddoux L, Brazier C, Maheo F, Simon JC, et al. Assessment of the Dominance Level of the R81T Target Resistance to Two Neonicotinoid Insecticides in Myzus persicae (Hemiptera: Aphididae). J Econ Entomol. 2016;109(5):2182–9. Epub 2016/08/09. pmid:27498842.
  12. 12. Shi M, Yue Z, Kuryatov A, Lindstrom JM, Sehgal A. Identification of Redeye, a new sleep-regulating protein whose expression is modulated by sleep amount. Elife. 2014;3:e01473. Epub 2014/02/06. pmid:24497543.
  13. 13. Bass C, Puinean AM, Andrews M, Cutler P, Daniels M, Elias J, et al. Mutation of a nicotinic acetylcholine receptor beta subunit is associated with resistance to neonicotinoid insecticides in the aphid Myzus persicae. Bmc Neuroscience. 2011;12. Artn 51 pmid:21627790
  14. 14. Shi XG, Zhu YK, Xia XM, Qiao K, Wang HY, Wang KY. The mutation in nicotinic acetylcholine receptor beta 1 subunit may confer resistance to imidacloprid in Aphis gossypii (Glover). J Food Agric Environ. 2012;10(2):1227–30.
  15. 15. Homem RA, Buttery B, Richardson E, Tan Y, Field LM, Williamson MS, et al. Evolutionary trade-offs of insecticide resistance—The fitness costs associated with target-site mutations in the nAChR of Drosophila melanogaster. Molecular Ecology. 2020;29(14):2661–75. pmid:32510730
  16. 16. Salgado VL. Studies on the mode of action of spinosad: Insect symptoms and physiological correlates. Pesticide Biochemistry and Physiology. 1998;60(2):91–102.
  17. 17. Watson GB, Loso MR, Babcock JM, Hasler JM, Letherer TJ, Young CD, et al. Novel nicotinic action of the sulfoximine insecticide sulfoxaflor. Insect Biochem Molec. 2011;41(7):432–9. pmid:21296156
  18. 18. Nauen R, Jeschke P, Velten R, Beck ME, Ebbinghaus-Kintscher U, Thielert W, et al. Flupyradifurone: a brief profile of a new butenolide insecticide. Pest Management Science. 2015;71(6):850–62. pmid:25351824
  19. 19. Hamada FN, Rosenzweig M, Kang K, Pulver SR, Ghezzi A, Jegla TJ, et al. An internal thermal sensor controlling temperature preference in Drosophila. Nature. 2008;454(7201):217–20. pmid:18548007.
  20. 20. Cordova D, Benner EA, Schroeder ME, Holyoke CW Jr., Zhang W, Pahutski TF, et al. Mode of action of triflumezopyrim: A novel mesoionic insecticide which inhibits the nicotinic acetylcholine receptor. Insect Biochem Mol Biol. 2016;74:32–41. Epub 2016/05/01. pmid:27130855.
  21. 21. Kitamoto T. Conditional disruption of synaptic transmission induces male-male courtship behavior in Drosophila. Proc Natl Acad Sci U S A. 2002;99(20):13232–7. pmid:12239352.
  22. 22. Tan J, Galligan JJ, Hollingworth RM. Agonist actions of neonicotinoids on nicotinic acetylcholine receptors expressed by cockroach neurons. Neurotoxicology. 2007;28(4):829–42. pmid:17561262
  23. 23. Nauen R, Ebbinghaus-Kintscher U, Salgado VL, Kaussmann M. Thiamethoxam is a neonicotinoid precursor converted to clothianidin in insects and plants. Pesticide Biochemistry and Physiology. 2003;76(2):55–69.
  24. 24. Salgado VL. Selective actions of insecticides on desensitizing and non-desensitizing nicotinic acetylcholine receptors in cockroach (Periplaneta americana) neurons. Pest Management Science. 2021;77(8):3663–72. pmid:33821538
  25. 25. Bao WX, Narai Y, Nakano A, Kaneda T, Murai T, Sonoda S. Spinosad resistance of melon thrips, Thrips palmi, is conferred by G275E mutation in alpha6 subunit of nicotinic acetylcholine receptor and cytochrome P450 detoxification. Pestic Biochem Physiol. 2014;112:51–5. Epub 2014/06/30. pmid:24974117.
  26. 26. Baxter SW, Chen M, Dawson A, Zhao JZ, Vogel H, Shelton AM, et al. Mis-spliced transcripts of nicotinic acetylcholine receptor alpha6 are associated with field evolved spinosad resistance in Plutella xylostella (L.). PLoS Genet. 2010;6(1):e1000802. Epub 2010/01/12. pmid:20062520.
  27. 27. Puinean AM, Lansdell SJ, Collins T, Bielza P, Millar NS. A nicotinic acetylcholine receptor transmembrane point mutation (G275E) associated with resistance to spinosad in Frankliniella occidentalis. J Neurochem. 2013;124(5):590–601. Epub 2012/09/29. pmid:23016960.
  28. 28. Rinkevich FD, Chen M, Shelton AM, Scott JG. Transcripts of the nicotinic acetylcholine receptor subunit gene Pxylalpha6 with premature stop codons are associated with spinosad resistance in diamondback moth, Plutella xylostella. Invert Neurosci. 2010;10(1):25–33. Epub 2010/05/26. pmid:20499126.
  29. 29. Silva WM, Berger M, Bass C, Williamson M, Moura DM, Ribeiro LM, et al. Mutation (G275E) of the nicotinic acetylcholine receptor alpha6 subunit is associated with high levels of resistance to spinosyns in Tuta absoluta (Meyrick) (Lepidoptera: Gelechiidae). Pestic Biochem Physiol. 2016;131:1–8. Epub 2016/06/07. pmid:27265820.
  30. 30. Urena E, Guillem-Amat A, Couso-Ferrer F, Beroiz B, Perera N, Lopez-Errasquin E, et al. Multiple mutations in the nicotinic acetylcholine receptor Ccalpha6 gene associated with resistance to spinosad in medfly. Sci Rep. 2019;9(1):2961. Epub 2019/03/01. pmid:30814521.
  31. 31. Wang HT, Tsai CL, Chen ME. Nicotinic acetylcholine receptor subunit alpha6 associated with spinosad resistance in Rhyzopertha dominica (Coleoptera: Bostrichidae). Pestic Biochem Physiol. 2018;148:68–73. Epub 2018/06/13. pmid:29891379.
  32. 32. Wang J, Wang X, Lansdell SJ, Zhang J, Millar NS, Wu Y. A three amino acid deletion in the transmembrane domain of the nicotinic acetylcholine receptor alpha6 subunit confers high-level resistance to spinosad in Plutella xylostella. Insect Biochem Mol Biol. 2016;71:29–36. Epub 2016/02/09. pmid:26855198.
  33. 33. Lansdell SJ, Collins T, Goodchild J, Millar NS. The Drosophila nicotinic acetylcholine receptor subunits Dalpha5 and Dalpha7 form functional homomeric and heteromeric ion channels. BMC Neurosci. 2012;13:73. Epub 2012/06/26. pmid:22727315.
  34. 34. Watson GB, Chouinard SW, Cook KR, Geng C, Gifford JM, Gustafson GD, et al. A spinosyn-sensitive Drosophila melanogaster nicotinic acetylcholine receptor identified through chemically induced target site resistance, resistance gene identification, and heterologous expression. Insect Biochem Mol Biol. 2010;40(5):376–84. Epub 2009/12/01. pmid:19944756.
  35. 35. Perry T, Chen W, Ghazali R, Yang YT, Christesen D, Martelli F, et al. Role of nicotinic acetylcholine receptor subunits in the mode of action of neonicotinoid, sulfoximine and spinosyn insecticides in Drosophila melanogaster. Insect Biochem Mol Biol. 2021;131:103547. Epub 2021/02/07. pmid:33548485.
  36. 36. Chamaon K, Schulz R, Smalla KH, Seidel B, Gundelfinger ED. Neuronal nicotinic acetylcholine receptors of Drosophila melanogaster: the alpha-subunit dalpha3 and the beta-type subunit ARD co-assemble within the same receptor complex. FEBS Lett. 2000;482(3):189–92. Epub 2000/10/12. pmid:11024458.
  37. 37. Chamaon K, Smalla KH, Thomas U, Gundelfinger ED. Nicotinic acetylcholine receptors of Drosophila: three subunits encoded by genomically linked genes can co-assemble into the same receptor complex. Journal of Neurochemistry. 2002;80(1):149–57. pmid:11796753
  38. 38. Li J, Shao Y, Ding Z, Bao H, Liu Z, Han Z, et al. Native subunit composition of two insect nicotinic receptor subtypes with differing affinities for the insecticide imidacloprid. Insect Biochem Mol Biol. 2010;40(1):17–22. Epub 2009/12/17. pmid:20005950.
  39. 39. Tsvetkov N, Samson-Robert O, Sood K, Patel HS, Malena DA, Gajiwala PH, et al. Chronic exposure to neonicotinoids reduces honey bee health near corn crops. Science. 2017;356(6345):1395–7. Epub 2017/07/01. pmid:28663503.
  40. 40. Woodcock BA, Bullock JM, Shore RF, Heard MS, Pereira MG, Redhead J, et al. Country-specific effects of neonicotinoid pesticides on honey bees and wild bees. Science. 2017;356(6345):1393–5. Epub 2017/07/01. pmid:28663502.
  41. 41. Rundlof M, Andersson GK, Bommarco R, Fries I, Hederstrom V, Herbertsson L, et al. Seed coating with a neonicotinoid insecticide negatively affects wild bees. Nature. 2015;521(7550):77–80. pmid:25901681.
  42. 42. Crall JD, Switzer CM, Oppenheimer RL, Ford Versypt AN, Dey B, Brown A, et al. Neonicotinoid exposure disrupts bumblebee nest behavior, social networks, and thermoregulation. Science. 2018;362(6415):683–6. Epub 2018/11/10. pmid:30409882.
  43. 43. Cressey D. Neonics Vs Bees. Nature. 2017;551(7679):156–8. pmid:29120435
  44. 44. Siviter H, Brown MJF, Leadbeater E. Sulfoxaflor exposure reduces bumblebee reproductive success. Nature. 2018;561(7721):109–12. pmid:30111837.
  45. 45. Haas J, Zaworra M, Glaubitz J, Hertlein G, Kohler M, Lagojda A, et al. A toxicogenomics approach reveals characteristics supporting the honey bee (Apis mellifera L.) safety profile of the butenolide insecticide flupyradifurone. Ecotoxicol Environ Saf. 2021;217:112247. Epub 2021/04/27. pmid:33901780.
  46. 46. Tosi S, Nieh JC, Brandt A, Colli M, Fourrier J, Giffard H, et al. Long-term field-realistic exposure to a next-generation pesticide, flupyradifurone, impairs honey bee behaviour and survival. Commun Biol. 2021;4(1):805. Epub 2021/06/30. pmid:34183763.
  47. 47. Jones AK, Brown LA, Sattelle DB. Insect nicotinic acetylcholine receptor gene families: from genetic model organism to vector, pest and beneficial species. Invert Neurosci. 2007;7(1):67–73. Epub 2007/01/12. pmid:17216517.