Contribution of Ldace1 gene to acetylcholinesterase activity in Colorado potato beetle

Leveraging RNA Interference to Impact Insecticide Resistance in the Colorado Potato Beetle, Leptinotarsa decemlineata

The Colorado potato beetle, Leptinotarsa decemlineata Say, is a potato pest that can cause important economic losses to the potato industry worldwide. Diverse strategies have been deployed to target this insect such as biological control, crop rotation, and a variety of insecticides. Regarding the latter, this pest has demonstrated impressive abilities to develop resistance against the compounds used to regulate its spread. Substantial work has been conducted to better characterize the molecular signatures underlying this resistance, with the overarching objective of leveraging this information for the development of novel approaches, including RNAi-based techniques, to limit the damage associated with this insect. This review first describes the various strategies utilized to control L. decemlineata and highlights different examples of reported cases of resistances against insecticides for this insect. The molecular leads identified as potential players modulating insecticide resistance as well as the growing interest towards the use of RNAi aimed at these leads as part of novel means to control the impact of L. decemlineata are described subsequently. Finally, select advantages and limitations of RNAi are addressed to better assess the potential of this technology in the broader context of insecticide resistance for pest management.

The effects of short-term glyphosate-based herbicide exposure on insect gene expression profiles

Glyphosate-based herbicides (GBHs) are the most frequently used herbicides worldwide. The use of GBHs is intended to tackle weeds, but GBHs have been shown to affect the life-history traits and antioxidant defense system of invertebrates found in agroecosystems. Thus far, the effects of GBHs on detoxification pathways among invertebrates have not been sufficiently investigated. We performed two different experiments-1) the direct pure glyphosate and GBH treatment, and 2) the indirect GBH experiment via food-to examine the possible effects of environmentally relevant GBH levels on the survival of the Colorado potato beetle (Leptinotarsa decemlineata) and the expression profiles of their detoxification genes. As candidate genes, we selected four cytochrome P450 (CYP), three glutathione-S-transferase (GST), and two acetylcholinesterase (AChE) genes that are known to be related to metabolic or target-site resistances in insects. We showed that environmentally relevant levels of pure glyphosate and GBH increased the probability for higher mortality in the Colorado potato beetle larvae in the direct experiment, but not in the indirect experiment. The GBHs or glyphosate did not affect the expression profiles of the studied CYP, GST, or AChE genes; however, we found a large family-level variation in expression profiles in both the direct and indirect treatment experiments. These results suggest that the genes selected for this study may not be the ones expressed in response to glyphosate or GBHs. It is also possible that the relatively short exposure time did not affect gene expression profiles, or the response may have already occurred at a shorter exposure time. Our results show that glyphosate products may affect the survival of the herbivorous insect already at lower levels, depending on their sensitivity to pesticides.

The redlegged earth mite draft genome provides new insights into pesticide resistance evolution and demography in its invasive Australian range

Genomic data provide valuable insights into pest management issues such as resistance evolution, historical patterns of pest invasions and ongoing population dynamics. We assembled the first reference genome for the redlegged earth mite, Halotydeus destructor (Tucker, 1925), to investigate adaptation to pesticide pressures and demography in its invasive Australian range using whole-genome pool-seq data from regionally distributed populations. Our reference genome comprises 132 autosomal contigs, with a total length of 48.90 Mb. We observed a large complex of ace genes, which has presumably evolved from a long history of organophosphate selection in H. destructor and may contribute towards organophosphate resistance through copy number variation, target-site mutations and structural variants. In the putative ancestral H. destructor ace gene, we identified three target-site mutations (G119S, A201S and F331Y) segregating in organophosphate-resistant populations. Additionally, we identified two new para sodium channel gene mutations (L925I and F1020Y) that may contribute to pyrethroid resistance. Regional structuring observed in population genomic analyses indicates that gene flow in H. destructor does not homogenize populations across large geographic distances. However, our demographic analyses were equivocal on the magnitude of gene flow; the short invasion history of H. destructor makes it difficult to distinguish scenarios of complete isolation vs. ongoing migration. Nonetheless, we identified clear signatures of reduced genetic diversity and smaller inferred effective population sizes in eastern vs. western populations, which is consistent with the stepping-stone invasion pathway of this pest in Australia. These new insights will inform development of diagnostic genetic markers of resistance, further investigation into the multifaceted organophosphate resistance mechanism and predictive modelling of resistance evolution and spread.

Comparative genomic analysis of carboxylesterase genes in Tenebrio molitor and other four tenebrionids

Carboxylesterases (COEs) have various functions in wide taxons of organisms. In insects, COEs are important enzymes involved in the hydrolysis of a variety of ester-containing xenobiotics, neural signal transmission, pheromone degradation, and reproductive development. Understanding the diversity of COEs is basic to illustrate their functions. In this study, we identified 53, 105, 37, and 39 COEs from the genomes of Tenebrio molitor, Asbolus verucosus, Hycleus cichorii, and H. phaleratus in the superfamily of Tenebrionidea, respectively. Phylogenetic analysis showed that 234 COEs from these four species and those reported in Tribolium castaneum (63) could be divided into 12 clades and three major classes. The α-esterases significantly expanded in T. molitor, A. verucosus, and T. castaneum compared to dipteran and hymenopteran insects. In T. molitor, most COEs showed tissue and stage-specific but not a sex-biased expression. Our results provide insights into the diversity and evolutionary characteristics of COEs in tenebrionids, and lay a foundation for the functional characterization of COEs in the yellow mealworm.

Repeated exposure of fluazinam fungicides affects gene expression profiles yet carries no costs on a nontarget pest

Fungicides are used to control pathogenic fungi of crop species, but they have also been shown to alter behavioral, life history and fitness related traits of nontarget insects. Here, we tested the fungicide effects on feeding behavior, survival and physiology of the nontarget pest insect, the Colorado potato beetle (CPB) (Leptinotarsa decemlineata). Feeding behavior was studied by a choice test of adult beetles, which were allowed to choose between a control and a fungicide (fluazinam) treated potato leaf. Larval survival was recorded after 24 and 72 h exposure to control and fungicide-treated leaves with 2 different concentrations. The adults did not show fungicide avoidance behavior. Similarly, survival of the larvae was not affected by the exposure to fungicides. Finally, to understand the effects of fungicides at the physiological level (gene expression), we tested whether the larval exposure to fungicide alter the expression of 5 metabolic pathway and stress associated genes. Highest concentration and 72-h exposure caused upregulation of 1 cytochrome P450 (CYP9Z14v2) and 1 insecticide resistance gene (Ldace1), whereas metabolic detoxification gene (Ugt1) was downregulated. At 24-h exposure, highest concentration caused downregulation of another common detoxification gene (Gs), while both exposure times to lowest concentration caused upregulation of the Hsp70 stress tolerance gene. Despite these overall effects, there was a considerable amount of variation among different families in the gene expression levels. Even though the behavioral effects of the fungicide treatments were minor, the expression level differences of the studied genes indicate changes on the metabolic detoxifications and stress-related pathways.

Characterization of molecular and kinetic properties of two acetylcholinesterases from the Colorado potato beetle, Leptinotarsa decemlineata

The molecular and biochemical properties of two acetylcholinesterases (LdAChE1 and LdAChE2) from the Colorado potato beetle, Leptinotarsa decemlineata, were investigated in this study. Polyacrylamide gel electrophoresis in conjunction with western blotting with LdAChE1- or LdAChE2-specific antibodies suggested that LdAChE1 exists in a soluble form, whereas LdAChE2 exists in both soluble and amphiphilic forms with a glycophosphatidylinositol anchor. Both LdAChEs exist as homodimers with each monomer connected with a disulfide bond. LdAChE1 was the most highly expressed in the thorax followed by the head, leg, and abdomen, whereas LdAChE2 was the most highly expressed in the head, followed by the thorax, leg, and abdomen. The overall expression levels of LdAChE1, however, were higher than those of LdAChE2 in all examined tissues. Kinetic analysis using recombinant LdAChE1 and LdAChE2 showed that LdAChE2 has a 4.8-fold higher catalytic efficiency toward acetylthiocholine iodide compared to LdAChE1. LdAChE2 was more sensitive to organophosphorus and carbamate insecticides than LdAChE1. The addition of irreversibly phosphorylated LdAChE1 via paraoxon titration significantly reduced LdAChE2 inhibition by insecticides and glycoalkaloids, suggesting a sequestration role of soluble LdAChE1 in the chemical defense against xenobiotics. Taken together, LdAChE2 may be the main enzyme for synaptic transmission, thus serving as a toxicologically more relevant target, whereas the soluble LdAChE1 may function as a bioscavenger.

Sequence variation and regulatory variation in acetylcholinesterase genes contribute to insecticide resistance in different populations of Leptinotarsa decemlineata

Although insect herbivores are known to evolve resistance to insecticides through multiple genetic mechanisms, resistance in individual species has been assumed to follow the same mechanism. While both mutations in the target site insensitivity and increased amplification are known to contribute to insecticide resistance, little is known about the degree to which geographic populations of the same species differ at the target site in a response to insecticides. We tested structural (e.g., mutation profiles) and regulatory (e.g., the gene expression of Ldace1 and Ldace2, AChE activity) differences between two populations (Vermont, USA and Belchow, Poland) of the Colorado potato beetle, Leptinotarsa decemlineata in their resistance to two commonly used groups of insecticides, organophosphates, and carbamates. We established that Vermont beetles were more resistant to azinphos-methyl and carbaryl insecticides than Belchow beetles, despite a similar frequency of resistance-associated alleles (i.e., S291G) in the Ldace2 gene. However, the Vermont population had two additional amino acid replacements (G192S and F402Y) in the Ldace1 gene, which were absent in the Belchow population. Moreover, the Vermont population showed higher expression of Ldace1 and was less sensitive to AChE inhibition by azinphos-methyl oxon than the Belchow population. Therefore, the two populations have evolved different genetic mechanisms to adapt to organophosphate and carbamate insecticides.

Can Indirect Herbicide Exposure Modify the Response of the Colorado Potato Beetle to an Organophosphate Insecticide?

Organisms live in complex multivariate environments. In agroecosystems, this complexity is often human-induced as pest individuals can be exposed to many xenobiotics simultaneously. Predicting the effects of multiple stressors can be problematic, as two or more stressors can have interactive effects. Our objective was to investigate whether indirect glyphosate-based herbicide (GBH) exposure of the host plant has interactive effects in combination with an insecticide (azinphos-methyl) on an invasive pest Colorado potato beetle (Leptinotarsa decemlineata Say). We tested the effects of GBH and insecticide on the survival, insecticide target genes expression (acetylcholinesterase genes) and oxidative status biomarkers (glutathione S-transferase [GST], glucose-6-phosphate dehydrogenase [G6PDH], glutathione reductase homolog [GR], glutathione peroxidase homolog [GPx], total glutathione [totGSH], glutathione reduced-oxidized [GSH: GSSG], catalase [CAT], superoxide dismutase [SOD], lipid hydroperoxides). We found that exposure to indirect GBH has no single or interactive effects in combination with the insecticide on larval survival. However, prior exposure to GBH inhibits Ldace1 gene expression by 0.55-fold, which is the target site for the organophosphate and carbamate insecticides. This difference disappears when individuals are exposed to both GBH and insecticide, suggesting an antagonistic effect. On the other hand, oxidative status biomarker scores (PCAs of GPx, GR, and CAT) were decreased when exposed to both stressors, indicating a synergistic effect. Overall, we found that indirect GBH exposure can have both antagonistic and synergistic effects in combination with an insecticide, which should be considered when aiming for an ecologically relevant risk assessment of multiple human-induced stressors.

Isolation of CarE genes from the Chinese white pine beetle Dendroctonus armandi (Curculionidae: Scolytinae) and their response to host chemical defense

BACKGROUND

Bark beetles rely on detoxifying enzymes to resist the defensive terpenoids of their host trees. Research on carboxylesterases (CarEs) has focused on their multiple functions in the metabolic detoxification of pesticides and plant allelochemicals, drug resistance, and juvenile hormone and pheromone degradation.

RESULT

We identified eight new CarE genes in the Chinese white pine beetle (Dendroctonus armandi) and carried out bioinformatics analysis on the deduced full-length amino acid sequences. Differential transcript levels of CarE genes were observed between sexes; within these levels, significant differences were found among the different development stages, and between insects fed on the phloem of Pinus armandi and exposed to five stimuli [(-)-α-pinene, (-)-β-pinene, (+)-3-carene, limonene and turpentine] at 8 and 24 h.

CONCLUSION

Transcription levels of CarE genes suggest some relationship with the detoxification of terpenoids released by host trees. The functions of bark beetle esterase are mainly in hydrolyzing the host chemical defense and degrading odorant molecules during host selection and colonization. © 2018 Society of Chemical Industry.

Effects of a glyphosate-based herbicide on survival and oxidative status of a non-target herbivore, the Colorado potato beetle (Leptinotarsa decemlineata)

Glyphosate is the globally most used herbicide against a wide range of weeds. Glyphosate has been considered safe to animals as it mainly targets physiological pathways in plants. However, recent toxicological studies have revealed that glyphosate can cause various toxic effects also on animals. In this study, we investigated the direct toxic effects of a glyphosate-based herbicide (GBH, Roundup® Bio) on 1) survival and 2) oxidative status of a non-target herbivore by using Colorado potato beetles (Leptinotarsa decemlineata), originating from Poland and USA, as model species. Larvae were randomly divided into three groups: 1) high concentration (100% Roundup Bio, 360 g/l), 2) low concentration (1.5% Roundup Bio) and 3) control group (water). Larvae were exposed to Roundup for different time periods: 2 h, 24 h, 48 h, 72 h and 96 h. Larval survival decreased in the group treated with high concentration of GBH compared to controls, whereas the low concentration group did not differ from the control group. GBH treatment had no association with oxidative status biomarkers (i.e. catalase, superoxide dismutase, glutathione-S-transferase, glutathione and glutathione related enzymes), but increased lipid hydroperoxide levels after 2 h exposure, suggesting increased oxidative damage soon after the exposure. Larvae of different origin also differed in their oxidative status, indicating population-dependent differences in antioxidant defence system. Environmentally relevant concentrations of GBH are not likely to affect larval survival, but high concentrations can reduce survival and increase oxidative damage of non-target herbivores. Also, populations of different origin and pesticide usage history can differ in their tolerance to GBH.

Molecular Characterization and Expression Analysis of Two Acetylcholinesterase Genes From the Small White Butterfly Pieris rapae (Lepidoptera: Pieridae)

Acetylcholinesterases (AChEs) are essential for the hydrolysis of the neurotransmitter acetylcholine and play crucial roles in the termination of neurotransmission. AChEs are encoded by the ace genes. However, the ace genes from the small white butterfly, Pieris rapae (L.) (Lepidoptera: Pieridae), remained uncharacterized. In this study, two aces (Prace1 and Prace2) were identified from P. rapae. Prace1 encoded a PrAChE1 protein consisting of 694 amino acid residues, and Prace2 encoded the 638-amino-acid PrAChE2. The two identified PrAChEs both had features typical of AChEs, including the catalytic triad, choline-binding sites, an oxyanion hole, an acyl pocket, a peripheral anionic subsite, an FGESAG motif and 14 conserved aromatic amino acids. Phylogenetic analysis showed that Prace1 and Prace2 were clustered into two distinct groups: ace1 and ace2, respectively. The two Praces were distributed on different genomic scaffolds: Prace1 on scaffold 156 and Prace2 on scaffold 430. Additionally, Prace1 consisted of three exons and two introns, whereas Prace2 consisted of six exons and five introns. One amino acid mutation (Gly324Ala) in PrAChE1 and two (Ser291Gly and Ser431Phe) in PrAChE2 were consistent with mutations in other insect AChEs that are associated with insecticide insensitivity. Both Prace1 and Prace2 were highly expressed at the fifth-instar larval stage and in the larval head, and the transcriptional levels of Prace1 were significantly higher than those of Prace2 in all of the tested life stages and tissues. This is the first report characterizing two ace genes in P. rapae. The results pave the way for functional study of these genes.

A Systematic Study of RNAi Effects and dsRNA Stability in Tribolium castaneum and Acyrthosiphon pisum, Following Injection and Ingestion of Analogous dsRNAs

RNA interference (RNAi) effects in insects are highly variable and may be largely dependent upon the stability of introduced double-stranded RNAs to digestion by nucleases. Here, we report a systematic comparison of RNAi effects in susceptible red flour beetle (Tribolium castaneum) and recalcitrant pea aphid (Acyrthosiphon pisum) following delivery of dsRNAs of identical length targeting expression of V-type ATPase subunit E (VTE) and inhibitor of apoptosis (IAP) genes. Injection and ingestion of VTE and IAP dsRNAs resulted in up to 100% mortality of T. castaneum larvae and sustained suppression (>80%) of transcript levels. In A. pisum, injection of VTE but not IAP dsRNA resulted in up to 65% mortality and transient suppression (ca. 40%) of VTE transcript levels. Feeding aphids on VTE dsRNA reduced growth and fecundity although no evidence for gene suppression was obtained. Rapid degradation of dsRNAs by aphid salivary, haemolymph and gut nucleases contrasted with stability in T. castaneum larvae where it appears that exo-nuclease activity is responsible for relatively slow digestion of dsRNAs. This is the first study to directly compare RNAi effects and dsRNA stability in receptive and refractory insect species and provides further evidence that dsRNA susceptibility to nucleases is a key factor in determining RNAi efficiency.

Identification and characterization of ace1-type acetylcholinesterase in insecticide-resistant and -susceptible Propylaea japonica (Thunberg)

Characterization and gene cloning of acetylecholinesterase (AChE) in the insecticide-resistant (R) and -susceptible (S) insects have been reported in the past. However, the studies focused mostly on herbivorous pests, rather than predacious species, such as ladybird beetles. Using R and S Propylaea japonica (thunberg), a full-length cDNA sequence (2928 bp) of the ace1-type AChE gene was determined for the first time. The ace1 encoding a protein of 645 amino acids contained typical conserved motifs, such as FGESAG domains, catalytic triad, acyl pocket, oxyanino hole, choline binding site, peripheral anionic site, omega loop and conserved aromatic residues. R P. japonica displayed 50-times greater resistance to chlorpyrifos or mathamidophos with a significantly lower AChE sensitivity to paraoxon, malaoxon, chlorpyrifos or methamidophos than its S counterpart. Five amino acids in the ace1 of R P. japonica differed from those found in S P. japonica. One of them, F358S, located in the acyl-binding pocket, might play a crucial role in the resistance of the insect to organophosphates (OPs). Whereas, K493E and I538V, which were close to some of the conserved aromatic amino acids (i.e., H509, Y511, and W499) in the gorge, and G571R and T576A near C593 that formed the disulfide bonds with C471, might also involve in the change of insecticide resistance in P. japonica. AChE insensitivity and amino acid replacements, particularly F358S, might be the determining factors in the alteration of OPs-resistance in P. japonica.

A model species for agricultural pest genomics: the genome of the Colorado potato beetle, Leptinotarsa decemlineata (Coleoptera: Chrysomelidae)

The Colorado potato beetle is one of the most challenging agricultural pests to manage. It has shown a spectacular ability to adapt to a variety of solanaceaeous plants and variable climates during its global invasion, and, notably, to rapidly evolve insecticide resistance. To examine evidence of rapid evolutionary change, and to understand the genetic basis of herbivory and insecticide resistance, we tested for structural and functional genomic changes relative to other arthropod species using genome sequencing, transcriptomics, and community annotation. Two factors that might facilitate rapid evolutionary change include transposable elements, which comprise at least 17% of the genome and are rapidly evolving compared to other Coleoptera, and high levels of nucleotide diversity in rapidly growing pest populations. Adaptations to plant feeding are evident in gene expansions and differential expression of digestive enzymes in gut tissues, as well as expansions of gustatory receptors for bitter tasting. Surprisingly, the suite of genes involved in insecticide resistance is similar to other beetles. Finally, duplications in the RNAi pathway might explain why Leptinotarsa decemlineata has high sensitivity to dsRNA. The L. decemlineata genome provides opportunities to investigate a broad range of phenotypes and to develop sustainable methods to control this widely successful pest.

RNA interference of three up-regulated transcripts associated with insecticide resistance in an imidacloprid resistant population of Leptinotarsa decemlineata

The Colorado potato beetle, Leptinotarsa decemlineata (Say), is a major agricultural pest of potatoes in the Central Sands production region of Wisconsin. Previous studies have shown that populations of L. decemlineata have become resistant to many classes of insecticides, including the neonicotinoid insecticide, imidacloprid. Furthermore, L. decemlineata has multiple mechanisms of resistance to deal with a pesticide insult, including enhanced metabolic detoxification by cytochrome p450s and glutathione S-transferases. With recent advances in the transcriptomic analysis of imidacloprid susceptible and resistant L. decemlineata populations, it is possible to investigate the role of candidate genes involved in imidacloprid resistance. A recently annotated transcriptome analysis of L. decemlineata was obtained from select populations of L. decemlineata collected in the Central Sands potato production region, which revealed a subset of mRNA transcripts constitutively up-regulated in resistant populations. We hypothesize that a portion of the up-regulated transcripts encoding for genes within the resistant populations also encode for pesticide resistance and can be suppressed to re-establish a susceptible phenotype. In this study, a discrete set of three up-regulated targets were selected for RNA interference experiments using a resistant L. decemlineata population. Following the successful suppression of transcripts encoding for a cytochrome p450, a cuticular protein, and a glutathione synthetase protein in a select L. decemlineata population, we observed reductions in measured resistance to imidacloprid that strongly suggest these genes control essential steps in imidacloprid metabolism in these field populations.

EXPRESSION AND EFFECTS OF MUTANT Bombyx mori ACETYLCHOLINESTRASE1 IN BmN CELLS

The main mechanism of toxicity of organophosphate (OP) and carbamate (CB) insecticides is their irreversible binding and inhibition of acetylcholinestrase (AChE), encoded by ace1 (acetylcholinestrase gene 1), leading to eventual death of insects. Mutations in AChE may significantly reduce insects susceptibility to these pesticides. Bombyx mori is an important beneficial insect, and no OP- or CB-resistant strains have been generated. In this study, wild-type ace1 (wace1) and mutant ace1 (mace1) were introduced into BmN cells, confirmed by screening and identification. The expression of wace1 and mace1 in the cells was confirmed by Western blot and their expression levels were about 21-fold higher than the endogenous ace1 level. The activities of AChE in wace1 and mace1 transgenic cells were 10.6 and 20.2% higher compared to control cells, respectively. mace1 transgenic cells had higher remaining activity than wace1 transgenic cells under the treatment of physostigmine (a reversible cholinesterase inhibitor) and phoxim (an OP acaricide). The results showed that ace1 transgene can significantly improve ace1 expression, and ace1 mutation at a specific site can reduce the sensitivity to AChE inhibitors. Our study provides a new direction for the exploration of the relationship between AChE mutations and drug resistance.

Target site insensitivity mutations in the AChE enzyme confer resistance to organophosphorous insecticides in Leptinotarsa decemlineata (Say)

In the present study, we demonstrated the use and optimization of the tetra-primer ARMS-PCR procedure to detect and analyze the frequency of the R30K and I392T mutations in resistant field populations of CPB. The R30K mutation was detected in 72%, 84%, 52% and 64% of Bahar, Dehpiaz, Aliabad and Yengijeh populations, respectively. Overall frequencies of the I392T mutation were 12%, 8% and 16% of Bahar, Aliabad and Yengijeh populations, respectively. No I392T point mutation was found among samples from Dehpiaz field population. Moreover, only 31% and 2% of samples from the resistant field populations were homozygous for R30K and I392T mutations, respectively. No individual simultaneously had both I392T and S291G/R30K point mutations. The incidence of individuals with both S291G and R30K point mutations in the samples from Bahar, Dehpiaz, Aliabad, and Yengijeh populations were 31.5%, 44.7%, 41.6%, and 27.3% respectively. Genotypes determined by the tetra-primer ARMS-PCR method were consistent with those determined by PCR sequencing. There was no significant correlation between the mutation frequencies and resistance levels in the resistant populations, indicating that other mutations may contribute to this variation. Polymorphism in the partial L. decemlineata cDNA AChE gene Ldace2 of four field populations was identified by direct sequencing of PCR-amplified fragments. Among 45 novel mutations detected in this study, T29P mutation was found across all four field populations that likely contribute to the AChE insensitivity. Site-directed mutagenesis and protein expression experiments are needed for a more complete evaluation.

Identification of carboxylesterase genes and their expression profiles in the Colorado potato beetle Leptinotarsa decemlineata treated with fipronil and cyhalothrin

Based on the Leptinotarsa decemlineata transcriptome dataset and the GenBank sequences, 70 novel carboxylesterases and 2 acetylcholinesterases were found. The 72 members belong to a multifunctional carboxylesterase/cholinesterase superfamily (CCE). A phylogenetic tree including the 72 LdCCEs and the CCEs from Tribolium castaneum, Drosophila melanogaster and Apis mellifera revealed that all CCEs fell into three main phylogenetic groups: dietary/detoxification, hormone/semiochemical processing, and neurodevelopmental classes. Numbers of L. decemlineata CCEs in the three classes were 52, 12 and 8, respectively. The dietary/detoxification class includes two clades: coleopteran xenobiotic metabolizing and α-esterase type CCEs. CCEs in the two clades have independently expanded in L. decemlineata. The hormone/semiochemical processing class has three clades: integument CCEs, β- and pheromone CCEs and juvenile hormone CCEs. Integument CCEs in L. decemlineata have also expanded. The neurodevelopmental CCEs are implicated the most ancient class, containing acetylcholinesterase, neuroligin, neurotactin, glutactin, gliotactin and others. Among the 70 novel CCE genes, KM220566, KM220530, KM220576, KM220527 and KM220541 were fipronil-inducible, and KM220578, KM220566, KM220542, KM220564, KM220561, KM220554, KM220527, KM220538 and KM220541 were cyhalothrin-inducible. They were the candidates involving in insecticide detoxification. Moreover, our results also provided a platform to understand the functions and evolution of L. decemlineata CCE genes.

Colorado potato beetle (Coleoptera) gut transcriptome analysis: expression of RNA interference-related genes

In the search for new methods of pest control, the potential of RNA interference (RNAi) is being explored. Because the gut is the first barrier for the uptake of double-stranded (ds)RNA, pyrosequencing of the gut transcriptome is a powerful tool for obtaining the necessary sequences for specific dsRNA-mediated pest control. In the present study, a dataset representing the gut transcriptome of the Colorado potato beetle (CPB; Leptinotarsa decemlineata) was generated and analysed for the presence of RNAi-related genes. Almost all selected genes that were implicated in silencing efficiency at different levels in the RNAi pathway (core machinery, associated intracellular factors, dsRNA uptake, antiviral RNAi, nucleases), which uses different types of small RNA (small interfering RNA, microRNA and piwi-RNA), were expressed in the CPB gut. Although the database is of lower quality, the majority of the RNAi genes are also found to be present in the gut transcriptome of the tobacco hornworm [TH; Manduca sexta (19 out of 35 genes analysed)]. The high quality of the CPB transcriptome database will lay the foundation for future gene expression and functional studies regarding the gut and RNAi.

Limitations of RNAi of α6 nicotinic acetylcholine receptor subunits for assessing the in vivo sensitivity to spinosad

Spinosad is a widely used insecticide that exerts its toxic effect primarily through interactions with the nicotinic acetylcholine receptor. The α6 nicotinic acetylcholine receptor subunit is involved in spinosad toxicity as demonstrated by the high levels of resistance observed in strains lacking α6. RNAi was performed against the Dα6 nicotinic acetylcholine receptor subunit in Drosophila melanogaster using the Gal4-UAS system to examine if RNAi would yield results similar to those of Dα6 null mutants. These Dα6-deficient flies were subject to spinosad contact bioassays to evaluate the role of the Dα6 nicotinic acetylcholine receptor subunit on spinosad sensitivity. The expression of Dα6 was reduced 60%-75% as verified by quantitative polymerase chain reaction. However, there was no change in spinosad sensitivity in D. melanogaster. We repeated RNAi experiments in Tribolium castaneum using injection of dsRNA for Tcasα6. RNAi of Tcasα6 did not result in changes in spinosad sensitivity, similar to results obtained with D. melanogaster. The lack of change in spinosad sensitivity in both D. melanogaster and T. castaneum using two routes of dsRNA administration shows that RNAi may not provide adequate conditions to study the role of nicotinic acetylcholine receptor subunits on insecticide sensitivity due to the inability to completely eliminate expression of the α6 subunit in both species. Potential causes for the lack of change in spinosad sensitivity are discussed.

Pre-invasion history and demography shape the genetic variation in the insecticide resistance-related acetylcholinesterase 2 gene in the invasive Colorado potato beetle

Background

Invasive pest species offers a unique opportunity to study the effects of genetic architecture, demography and selection on patterns of genetic variability. Invasive Colorado potato beetle (Leptinotarsa decemlineata) populations have experienced a rapid range expansion and intense selection by insecticides. By comparing native and invasive beetle populations, we studied the origins of organophosphate (OP) resistance-associated mutations in the acetylcholinesterase 2 (AChE2) gene, and the role of selection and demography on its genetic variability.

Results

Analysis of three Mexican, two US and five European populations yielded a total of 49 haplotypes. Contrary to the expectations all genetic variability was associated with a point mutation linked to insecticide resistance (S291G), this mutation was found in 100% of Mexican, 95% of US and 71% of European beetle sequences analysed. Only two susceptible haplotypes, genetically very differentiated, were found, one in US and one in Europe. The genetic variability at the AChE2 gene was compared with two other genes not directly affected by insecticide selection, diapause protein 1 and juvenile hormone esterase. All three genes showed reduction in genetic variability indicative of a population bottleneck associated with the invasion.

Conclusions

Stochastic effects during invasion explain most of the observed patterns of genetic variability at the three genes investigated. The high frequency of the S291G mutation in the AChE2 gene among native populations suggests this mutation is the ancestral state and thus, either a pre-adaptation of the beetle for OP resistance or the AChE2 is not the major gene conferring OP resistance. The long historical association with host plant alkaloids together with recombination may have contributed to the high genetic variation at this locus. The genetic diversity in the AChE2 locus of the European beetles, in turn, strongly reflects founder effects followed by rapid invasion. Our results suggest that despite the long history of insecticide use in this species, demographic events together with pre-invasion history have been strongly influential in shaping the genetic diversity of the AChE2 gene in the invasive beetle populations.