Genomics-based approaches to screening carboxylesterase-like genes potentially involved in malathion resistance in oriental migratory locust (Locusta migratoria manilensis)

Crofton weed derived isomers of ageraphorone as potent antifeedant against Plutella xylostella (L.)

Global agricultural production is significantly hampered by insect pests, and the demand for natural pragmatic pesticides with environmental concern remains unfulfilled. Ageratina adenophora (Spreng.) also known as Crofton weed, is an invasive perennial herbaceous plant that is known to possess multiple bioactive compounds. In our study, two isomers of ageraphorone metabolites i.e, 10 Hα-9-oxo-ageraphorone (10HA) and 10 Hβ-9-oxo-ageraphorone (10HB), were identified from Crofton weed, exhibiting potent antifeedant and larvicidal activities against Plutella xylostella. For antifeedant activity, the median effective concentration (EC50) values for 10HA and 10HB in the choice method were 2279 mg/L and 3233 mg/L, respectively, and for the no choice method, EC50 values were 1721 mg/L and 2394 mg/L, respectively. For larvicidal activity, lethal concentration (LC50) values for 10HA and 10HB were 2421 mg/L and 4109 mg/L at 48 h and 2101 mg/L and 3550 mg/L at 72 h. Furthermore, both in- vivo and in-vitro studies revealed that the isomers 10HA and 10HB exhibited potent detoxifying enzymes inhibition activity such as carboxylesterase and glutathione S-transferases. Molecular docking and MD simulation analysis provide insight into the possible interaction between isomers of ageraphorone metabolites and Carboxylic Ester Hydrolase protein (Gene: pxCCE016b) of P. xylostella, which led to a finding that CarEH protein plays a significant role in the detoxification of the two compounds in P. xylostella. Finally, our findings show that the primary enzymes undergoing inhibition by isomers of ageraphorone metabolites, causing toxicity in insects, are Carboxylesterase and glutathione S-transferase.

The effect of chlorogenic acid, a potential botanical insecticide, on gene transcription and protein expression of carboxylesterases in the armyworm (Mythimna separata)

Chlorogenic acid (CGA) is a potential botanical insecticide metabolite that naturally occurs in various plants. Our previous studies revealed CGA is sufficient to control the armyworm Mythimna separata. In this study, we conducted a proteomic analysis of saliva collected from M. separata following exposure to CGA and found that differentially expressed proteins (DEPs) treated with CGA for 6 h and 24 h were primarily enriched in glutathione metabolism and the pentose phosphate pathway. Notably, we observed six carboxylesterase (CarE) proteins that were enriched at both time points. Additionally, these corresponding genes were expressed at levels 5.05 to 130.25 times higher in our laboratory-selected resistance strains. We also noted a significant increase in the enzyme activity of carboxylesterase following treatments with varying CGA concentrations. Finally, we confirmed that knockdown of MsCarE14, MsCarE28, and MsCCE001h decreased the susceptibility to CGA in resistance strain, indicating three CarE genes play crucial roles in CGA detoxification. This study presents the first report on the salivary proteomics of M. separata, offering valuable insights into the role of salivary proteins. Moreover, the determination of CarE mediated susceptibility change to CGA provides new targets for agricultural pest control and highlights the potential insecticide resistance mechanism for pest resistance management.

Identification, expression profiles and involvement in insecticides tolerance and detoxification of carboxylesterase genes in Bactrocera dorsalis

Carboxylesterases (CarEs) are a multifunctional superfamily of enzymes and play an important role in detoxification of various insecticides in insects. The oriental fruit fly, Bactrocera dorsalis, is one of the most destructive agricultural pests and has developed different degrees of resistance to organophosphates in field. However, the involvement of BdCarEs in tolerance or resistance to other alternative insecticides are still unclear. In the present study, 33 BdCarEs genes were identified based on the genome database of B. dorsalis. Phylogenetic analysis demonstrated that they were classified into nine clades, with abundance of α-esterases. Meanwhile, the sequence characterization and the chromosome distribution were also analyzed. The spatiotemporal expression analysis of BdCarEs genes suggested that the diversity of potential function in different physiological processes. With the exception of BdCarE21, all BdCarEs genes responded to at least one insecticide exposure, and BdCarE20 was found to be up-regulated after exposure to all five tested insecticides individually. Eight BdCarEs genes were overexpressed in MR strain when compared to that in SS strain. Subsequently, knockdown the expression of representative BdCarEs genes significantly increased the susceptibility of the oriental fruit fly to corresponding insecticides, which indicated that the tested BdCarEs genes contributed to one or multiple insecticide detoxification. These findings provide valuable insights into the potential role in respond to tolerance or resistance to insecticides with different mode of action, and will facilitate development of efficiency management strategy for B. dorsalis.

Identification and characterization of candidate detoxification genes in Pharsalia antennata Gahan (Coleoptera: Cerambycidae)

The wood-boring beetles, including the majority of Cerambycidae, have developed the ability to metabolize a variety of toxic compounds derived from host plants and the surrounding environment. However, detoxification mechanisms underlying the evolutionary adaptation of a cerambycid beetle Pharsalia antennata to hosts and habitats are largely unexplored. Here, we characterized three key gene families in relation to detoxification (cytochrome P450 monooxygenases: P450s, carboxylesterases: COEs and glutathione-S-transferases: GSTs), by combinations of transcriptomics, gene identification, phylogenetics and expression profiles. Illumina sequencing generated 668,701,566 filtered reads in 12 tissues of P. antennata, summing to 100.28 gigabases data. From the transcriptome, 215 genes encoding 106 P450s, 77 COEs and 32 GSTs were identified, of which 107 relatives were differentially expressed genes. Of the identified 215 genes, a number of relatives showed the orthology to those in Anoplophora glabripennis, revealing 1:1 relationships in 94 phylogenetic clades. In the trees, P. antennata detoxification genes mainly clustered into one or two subfamilies, including 64 P450s in the CYP3 clan, 33 COEs in clade A, and 20 GSTs in Delta and Epsilon subclasses. Combining transcriptomic data and PCR approaches, the numbers of detoxification genes expressed in abdomens, antennae and legs were 188, 148 and 141, respectively. Notably, some genes exhibited significantly sex-biased levels in antennae or legs of both sexes. The findings provide valuable reference resources for further exploring xenobiotics metabolism and odorant detection in P. antennata.

The metabolic resistance of Nilaparvata lugens to chlorpyrifos is mainly driven by the carboxylesterase CarE17

The involvement of carboxylesterases (CarEs) in resistance to chlorpyrifos has been confirmed by the synergism analysis in Nilaparvata lugens. However, the function of specific CarE gene in chlorpyrifos resistance and the transcriptional regulatory mechanism are obscure. Herein, the expression patterns of 29 CarE genes in the susceptible and chlorpyrifos-resistant strains were analyzed. Among them, CarE3, CarE17 and CarE19 were overexpressed in the resistant strain, and knockdown of either CarE gene by RNA interference significantly increased the susceptibility to chlorpyrifos. Remarkably, knockdown of CarE17 reduced the enzymatic activity of CarE by 88.63 % and showed a much greater effect on increasing chlorpyrifos toxicity than silencing other two CarE genes. Overexpression of CarE17 in Drosophila melanogaster decreased the toxicity of chlorpyrifos to transgenic fruit flies. Furthermore, the region between - 205 to + 256 of CarE17 promoter sequence showed the highest promoter activity, and 16 transcription factors (TFs) were predicted from this region. Among these TFs, Lim1β and C15 were overexpressed in the resistant strain. Knockdown of either TF resulted in reduced CarE17 expression and a decrease in resistance of N. lugens to chlorpyrifos. These results indicate that the constitutive overexpression of Lim1β and C15 induces CarE17 expression thus conferring chlorpyrifos resistance in N. lugens.

Characterization of the insecticide detoxification carboxylesterase Boest1 from Bradysia odoriphaga Yang et Zhang (Diptera: Sciaridae)

BACKGROUND

In insects, carboxylesterases (CarEs) are enzymes involved in the detoxification of insecticides. However, the molecular mechanism of CarE-mediated insecticide metabolism in Bradysia odoriphaga, a serious agricultural pest, remains unclear. The aim of this study is to investigate the detoxification process of malathion, bifenthrin, and imidacloprid by B. odoriphaga carboxylesterase (Boest1).

RESULTS

An alpha class CarE gene Boest1 was cloned from B. odoriphaga. The results of real-time quantitative polymerase chain reaction showed that Boest1 is up-regulated with age during the larval stage, and the level of transcription of Boest1 is higher in the midgut and Malpighian tubule than in other tissues. The expression level of Boest1 was significantly increased after exposure to malathion and bifenthrin. Recombinant BoEST1 expressed in vitro showed high catalytic activity toward α-naphthyl acetate, which was substantially inhibited by malathion and triphenyl phosphate. The in vitro metabolism assays showed that BoEST1 demonstrates hydrolytic capacity toward malathion and bifenthrin but not imidacloprid. The binding free energy analysis indicates that BoEST1 has a higher affinity for malathion and bifenthrin than imidacloprid.

CONCLUSION

These results suggest that BoEST1 plays a role in the breakdown of insecticides and may be involved in the development of resistance in the Chinese chive pest B. odoriphaga; our findings also provide data for better pest management and perspectives for new pesticides development. © 2021 Society of Chemical Industry.

Functional analysis of a carboxylesterase gene involved in beta-cypermethrin and phoxim resistance in Plutella xylostella (L.)

BACKGROUND

Carboxylesterases (CarEs) are associated with detoxification of xenobiotics, including insecticides, in organism bodies. Overexpression of CarE genes is considered to have an important role in insecticide resistance in insects, however its involvement in multi-insecticide resistance has rarely been reported. This study aimed to assess the function of a CarE gene (PxαE8) in resistance to five insecticides in Plutella xylostella.

RESULTS

Relative expression of PxαE8 in three multi-insecticide-resistant Plutella xylostella populations, GD-2017, GD-2019 and HN, was14.8-, 19.5- and 28.0-fold higher than that in the susceptible population. Exposure to lethal concentrations associated with 25% mortality (LC₂₅ ) of beta-cypermethrin, chlorantraniliprole, metaflumizone, phoxim and tebufenozide could induce the specific activity of CarEs and increase the relative expression of PxαE8. By contrast, knockdown of PxαE8 expression dramatically reduced the activity of CarEs and increased the resistance of P. xylostella (GD-2019) larvae to beta-cypermethrin and phoxim by 47.4% and 45.5%, respectively. Further, a transgenic line of Drosophila melanogaster overexpressing PxαE8 was constructed and the bioassay results showed that the tolerance of transgenic Drosophila to beta-cypermethrin and phoxim was 3.93- and 3.98-fold higher than that of the untransgenic line.

CONCLUSION

These results provide evidence that overexpression of PxαE8 is involved in resistance, at least to beta-cypermethrin and phoxim, in multi-insecticide-resistant P. xylostella populations, which could help in further understanding the molecular mechanisms of multi-insecticide resistance in this pest. © 2020 Society of Chemical Industry.

Expression and kinetic analysis of carboxylesterase LmCesA1 from Locusta migratoria

OBJECTIVE

To investigate the biochemical characterization of the carboxylesterase LmCesA1 from Locusta migratoria.

RESULTS

We expressed recombinant LmCesA1 in Sf9 cells by using the Bac-to-bac baculovirus expression system. Enzyme kinetic assays showed that the K_m values of LmCesA1 for α-naphthyl acetate (α-NA) and β-naphthyl acetate (β-NA) were 0.08 ± 0.01 mM and 0.22 ± 0.03 mM, respectively, suggesting that LmCesA1 has a higher affinity for α-NA. LmCesA1 retained its enzymatic activity during incubations at pH 7-10 and at 10-30 °C. In an inhibition experiment, two organophosphate pesticides (malaoxon and malathion) and one pyrethroid pesticide (deltamethrin) showed different inhibition profiles against purified LmCesA1. Recombinant LmCesA1 activity was significantly inhibited by malaoxon in vitro. UPLC analysis showed that no metabolites were detected.

CONCLUSIONS

These results suggest that overexpression of LmCesA1 enhances malathion sequestration to confer malathion tolerance in L. migratoria.

Integrated analysis of physiological, transcriptomic and metabolomic responses and tolerance mechanism of nitrite exposure in Litopenaeus vannamei

Nitrite accumulation in aquatic environments is a potential risk factor that disrupts multiple physiological functions in aquatic animals. In this study, the physiology, transcriptome and metabolome of the control group (LV-C), nitrite-tolerance group (LV-NT) and nitrite-sensitive group (LV-NS) were investigated to identify the stress responses and mechanisms underlying the nitrite tolerance of Litopenaeus vannamei. After LV-NT and LV-NS were subjected to nitrite stress, the hemocyanin contents were significantly decreased, and hepatopancreas showed severe histological damage compared with LV-C. Likewise, the antioxidant enzymes were also significantly changed after nitrite exposure. The transcriptome data revealed differentially expressed genes associated with immune system, cytoskeleton remodeling and apoptosis in LV-NT and LV-NS. The combination of transcriptomic and metabolomic analysis revealed nitrite exposure disturbed metabolism processes in L. vannamei, including amino acid metabolism, nucleotide metabolism and lipid metabolism. The multiple comparative analysis implicated that higher nitrite tolerance of LV-NT than LV-NS may be attributed to enhanced hypoxia inducible factor-1α expression to regulate energy supply and gaseous exchange. Moreover, LV-NT showed higher antioxidative ability, detoxification gene expression and enhanced fatty acids contents after nitrite exposure in relative to LV-NS. Collectively, all these results will greatly provide new insights into the molecular mechanisms underlying the stress responses and tolerance of nitrite exposure in L. vannamei.

Identification and biochemical characterization of carboxylesterase 001G associated with insecticide detoxification in Helicoverpa armigera

Carboxylesterases (CarEs) are a major class of detoxification enzymes involved in insecticide resistance in various insect species. In this study, a novel CarE 001G was isolated from the cotton bollworm Helicoverpa armigera, one of the most destructive agricultural insect pests. The open reading frame of 001G has 2244 nucleotides and putatively encodes 747 amino acid residues. The deduced CarE possessed the highly conserved catalytic triads(Ser-Glu-His) and pentapeptide motifs (Gly-X-Ser-X-Gly), suggesting 001G is biologically active. The truncated 001G was successfully expressed in Escherichia coli, and the recombinant proteins were purified and tested. The enzyme kinetic assay showed the purified proteins could catalyze two model substrates, α-naphthyl acetate and β-naphthyl acetate, with a kcat of 8.8 and 2.3 s^-1, a Km of 9.6 and 16.2 μM, respectively. The inhibition study with pyrethroid, organophosphate and neonicotinoid insecticides showed different inhibition profile against the purified CarE. The HPLC assay demonstrated that the purified proteins were able to metabolize β-cypermethrin, λ-cyhalothrin and fenvalerate insecticides, exhibiting respective specific activities of 1.7, 1.4 and 0.5 nM/min/mg protein. However, the purified proteins were not able to metabolize the chlorpyrifos, parathion-methyl, paraoxon-ethyl and imidacloprid. The modeling and docking analyses consistently demonstrated that the pyrethroid molecule fits snugly into the catalytic pocket of the CarE 001G. Collectively, our results suggest that 001G may play a role in pyrethroids detoxification in H. armigera.

RNA interference technology in crop protection against arthropod pests, pathogens and nematodes

Scientists have made significant progress in understanding and unraveling several aspects of double-stranded RNA (dsRNA)-mediated gene silencing during the last two decades. Now that the RNA interference (RNAi) mechanism is well understood, it is time to consider how to apply the acquired knowledge to agriculture and crop protection. Some RNAi-based products are already available for farmers and more are expected to reach the market soon. Tailor-made dsRNA as an active ingredient for biopesticide formulations is considered a raw material that can be used for diverse purposes, from pest control and bee protection against viruses to pesticide resistance management. The RNAi mechanism works at the messenger RNA (mRNA) level, exploiting a sequence-dependent mode of action, which makes it unique in potency and selectivity compared with conventional agrochemicals. Furthermore, the use of RNAi in crop protection can be achieved by employing plant-incorporated protectants through plant transformation, but also by non-transformative strategies such as the use of formulations of sprayable RNAs as direct control agents, resistance factor repressors or developmental disruptors. In this review, RNAi is presented in an agricultural context (discussing products that have been launched on the market or will soon be available), and we go beyond the classical presentation of successful examples of RNAi in pest-insect control and comprehensively explore its potential for the control of plant pathogens, nematodes and mites, and to fight against diseases and parasites in beneficial insects. Moreover, we also discuss its use as a repressor for the management of pesticide-resistant weeds and insects. Finally, this review reports on the advances in non-transformative dsRNA delivery and the production costs of dsRNA, and discusses environmental considerations. © 2017 Society of Chemical Industry.

Overexpression of Mn-superoxide dismutase in Oxya chinensis mediates increased malathion tolerance

Superoxide dismutase (SOD) is the first line of defense against oxidative damage. Malathion is an organophosphate insecticide and can induce the production of reactive oxygen species (ROS) and cause the intracellular oxidative stress. The present study was undertaken to examine the effects of malathion on SODs activity and their transcriptional levels in Oxya chinensis (Thunberg) (Orthoptera: Acrididae). The results showed that total SOD and MnSOD activities increased as a dose-dependent manner while CuZnSOD activity has no significant changes after malathion treatments. Total SOD and MnSOD activities were the highest at the concentration of 0.8 μg μL^-1 malathion treatment and increased significantly about 1.81- and 2.48-fold compared with the control, respectively. Increased mRNA expression of MnSOD, ecCuZnSOD1, and ecCuZnSOD2 were observed after malathion treatments. Moreover, the alteration of MnSOD transcript was similar to the profiles of MnSOD activity. These results suggested that the up-regulation expression of MnSOD transcript led to the increase of MnSOD activity in order to eliminate the excessive ROS caused by malathion. In addition, we evaluated the role of individual SOD gene in malathion stress by using RNAi and recombinant SOD proteins. The results showed that ROS contents increased significantly after the silencing of MnSOD and ecCuZnSOD1 genes. The OD values of the E. coli cells transformed with pET-28a-OcMnSOD plasmid were 1.13-1.31-fold and 1.08-1.33-fold higher than those of cells with pET-28a plasmids under 0.4 and 0.8 μg μL^-1 malathion treatments, respectively. These findings indicated that MnSOD exerted an important role in defense oxidative stress caused by malathion.

Ligand-binding characterization of simulated β-adrenergic-like octopamine receptor in Schistocerca gregaria via progressive structure simulation

It is important to design insecticides having both low drug resistance and less undesirable toxicity for desert locust control. Specific GPCRs of Schistocerca gregaria, especially β-adrenergic-like octopamine receptor (SgOctβR), can be considered as its potential effective insecticide targets. However, either the unavailability of SgOctβR's structure or the inadequate capability of its sequence lead the development of insecticide for Schistocerca gregaria meets its plateau. To relax this difficulty, this paper develops a promising progressive structure simulation from SgOctβR's sequence, to its predicted structure of SgOctβR in vacuum, to its conformation as well as its complex with endogenous ligand octopamine in a solvent-membrane system. The combined approach of multiple sequence alignment, static structural characterization, and dynamic process of conformational change during binding octopamine reveal three important aspects. The first one is the characterization of SgOctβR's active pocket, including the attending secondary structure elements, its hydrophobic residues and nonpolar surface. The second one is the interaction with octopamine, especially the involved hydrogen bonds and an aromatic stacking of pi-pi interactions. The third one is the potential binding sites, including six highly conserved residues and one highly variable residue for locust insecticide design. This work is definitely helpful for the further structure-based drug design for efficient and eco-friendly insecticides, as well as site-directed mutagenesis biochemical research of SgOctβR.

Transcriptome Analysis and Identification of Major Detoxification Gene Families and Insecticide Targets in Grapholita Molesta (Busck) (Lepidoptera: Tortricidae)

The oriental fruit moth, Grapholita molesta (Busck) (Lepidoptera: Tortricidae), is an important pest of most stone and pome fruits and causes serious damage to the fruit industry worldwide. This insect pest has been primarily controlled through the application of insecticides; as a result, G. molesta has developed resistance to many different types of insecticides. To identify detoxification genes, we have, de novo, sequenced the transcriptome of G. molesta (Lepidoptera: Tortricidae) and yielded 58,970 unigenes of which 26,985 unigenes matched to known proteins. In total, 2,040 simple sequence repeats have been identified. The comprehensive transcriptome data set has permitted us to identify members of important gene families related to detoxification in G. molesta, including 77 unigenes of putative cytochrome P450s, 28 of glutathione S-transferases, 46 of Carboxylesterases, and 31 of insecticide targets. Orthologs of some of these unigenes have shown to play a pivotal role in insecticide resistance in other insect species and those unigenes likely have similar functions in G. molesta.

Microarray-based annotation of the gut transcriptome of the migratory locust, Locusta migratoria

The migratory locust, Locusta migratoria, is a serious agricultural pest and important insect model in the study of insect digestion and feeding behaviour. The gut is one of the primary interfaces between the insect and its environment. Nevertheless, knowledge on the gut transcriptome of L. migratoria is still very limited. Here, 48 802 expressed sequence tags were extracted from publicly available databases and their expression in larval gut and/or brain tissue was determined using microarray hybridization. Our data show 2765 transcripts predominantly or exclusively expressed in the gut. Many transcripts had putative functions closely related to the physiological functions of the gut as a muscular digestive organ and as the first barrier against microorganisms and a wide range of toxins. By means of a ranking procedure based on the relative signal intensity, we estimated 15% of the transcripts to show high expression levels, the highest belonging to diverse digestive enzymes and muscle-related proteins. We also found evidence for very high expression of an allergen protein, which could have important implications, as locusts form a traditional food source in various parts of the world, and were also recently added to the list of insects fit for human consumption in Europe. Interestingly, many highly expressed sequences have as yet unknown functions. Taken together, the present data provide significant insight into locust larval gut physiology, and will be valuable for future studies on the insect gut.

Transcriptional Analysis of The Adaptive Digestive System of The Migratory Locust in Response to Plant Defensive Protease Inhibitors

Herbivorous insects evolved adaptive mechanisms to compensate for the presence of plant defensive protease inhibitors (PI) in their food. The underlying regulatory mechanisms of these compensatory responses remain largely elusive. In the current study, we investigated the initiation of this adaptive response in the migratory locust, Locusta migratoria, via microarray analysis of gut tissues. Four hours after dietary uptake of PIs, 114 and 150 transcripts were respectively found up- or downregulated. The results suggest a quick trade-off between compensating for potential loss of digestive activity on the one hand, and stress tolerance, defense, and structural integrity of the gut on the other hand. We additionally addressed the role of a group of related upregulated hexamerin-like proteins in the PI-induced response. Simultaneous knockdown of corresponding transcripts by means of RNA interference resulted in a reduced capacity of the locust nymphs to cope with the effects of PI. Moreover, since insect hexamerins have been shown to bind Juvenile Hormone (JH), we also investigated the effect of JH on the proteolytic digestion in L. migratoria. Our results indicate that JH has a stimulatory effect on the expression of three homologous chymotrypsin genes, while knocking down the JH receptor (methoprene tolerant) led to opposite effects.

Insecticide resistance monitoring and metabolic mechanism study of the green peach aphid, Myzus persicae (Sulzer) (Hemiptera: Aphididae), in Chongqing, China

Myzus persicae (Sulzer) is one of the most important agricultural pests in China, which caused serious losses every year. For resistance monitoring, twelve populations of this pest were collected from tobacco field in Chongqing, China, and their sensitivities to 4 insecticides were tested. Results showed that only WL (RR=6.51) and FJ (RR=6.03) populations have developed minor resistance to imidacloprid, and the others still remained susceptible. One population (NC) has reached a high resistance level to cyhalothrin (RR=41.28), five populations showed medium level (10.36≤RR≤20.45), and the other six remained susceptible (0.39≤RR≤3.53). As regards carbosulfan, three populations have developed medium resistance, four populations showed only minor resistance, and the other five (0.81≤RR≤3.97) were still susceptible. Population SZ developed a medium level (RR=14.83) to phoxim, the other 11 were susceptible (0.29≤RR≤2.41). To analysis the potential resistance mechanism, inhibition effects of synergists and detoxifying enzyme activities were detected. The results indicated that the MFO was the most important detoxifying enzyme conferring imidacloprid resistance, and CarE was most important to cyhalothrin, carbosulfan and phoxim. Our study provided a comprehensive survey of insecticide resistance of M. persicae in Chongqing, and suggested that different counties should take corresponding management to delay the insecticide resistance development and prolong the usefulness of insecticides.

De Novo Assembly and Characterization of the Transcriptome of Grasshopper Shirakiacris shirakii

Background: The grasshopper Shirakiacris shirakii is an important agricultural pest and feeds mainly on gramineous plants, thereby causing economic damage to a wide range of crops. However, genomic information on this species is extremely limited thus far, and transcriptome data relevant to insecticide resistance and pest control are also not available. Methods: The transcriptome of S. shirakii was sequenced using the Illumina HiSeq platform, and we de novo assembled the transcriptome. Results: Its sequencing produced a total of 105,408,878 clean reads, and the de novo assembly revealed 74,657 unigenes with an average length of 680 bp and N50 of 1057 bp. A total of 28,173 unigenes were annotated for the NCBI non-redundant protein sequences (Nr), NCBI non-redundant nucleotide sequences (Nt), a manually-annotated and reviewed protein sequence database (Swiss-Prot), Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. Based on the Nr annotation results, we manually identified 79 unigenes encoding cytochrome P450 monooxygenases (P450s), 36 unigenes encoding carboxylesterases (CarEs) and 36 unigenes encoding glutathione S-transferases (GSTs) in S. shirakii. Core RNAi components relevant to miroRNA, siRNA and piRNA pathways, including Pasha, Loquacious, Argonaute-1, Argonaute-2, Argonaute-3, Zucchini, Aubergine, enhanced RNAi-1 and Piwi, were expressed in S. shirakii. We also identified five unigenes that were homologous to the Sid-1 gene. In addition, the analysis of differential gene expressions revealed that a total of 19,764 unigenes were up-regulated and 4185 unigenes were down-regulated in larvae. In total, we predicted 7504 simple sequence repeats (SSRs) from 74,657 unigenes. Conclusions: The comprehensive de novo transcriptomic data of S. shirakii will offer a series of valuable molecular resources for better studying insecticide resistance, RNAi and molecular marker discovery in the transcriptome.

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.

Two homologous carboxylesterase genes from Locusta migratoria with different tissue expression patterns and roles in insecticide detoxification

Carboxylesterases (CarEs) play a crucial role in detoxification of xenobiotics and resistance to insecticides in insects. In this study, two cDNAs of CarE genes (LmCesA4 and LmCesA5) were sequenced from the migratory locust, Locusta migratoria. The cDNAs of LmCesA4 and LmCesA5 putatively encoded 538 and 470 amino acid residues, respectively. The deduced amino acid sequences of the two CarE genes showed 45.0% identities, possessed highly conserved catalytic triads (Ser-Glu-His), and clustered in phylogenetic analysis. These results suggest that they are homologous genes. Both CarE genes were expressed throughout the developmental stages. However, LmCesA4 was predominately expressed in the midgut (including the gastric caeca) and fat bodies, whereas LmCesA5 was mainly expressed in the gastric caeca. The in situ hybridization results showed that the transcripts of the two genes were localized in apical and basal regions of the columnar cells in the gastric caeca. Gene silencing followed by insecticide bioassay increased the mortalities of deltamethrin-, malathion-, and carbaryl-treated locusts by 29.5%, 31.0% and 20.4%, respectively, after the locusts were injected with LmCesA4 double-stranded RNA (dsRNA). In contrast, the injection of LmCesA5 dsRNA did not significantly increase the susceptibility of the locusts to any of these insecticides. These results suggest that these genes not only show different tissue expression patterns but also play different roles in insecticide detoxification.

RNA interference: Applications and advances in insect toxicology and insect pest management

Since its discovery, RNA interference (RNAi) has revolutionized functional genomic studies due to its sequence-specific nature of post-transcriptional gene silencing. In this paper, we provide a comprehensive review of the recent literature and summarize the current knowledge and advances in the applications of RNAi technologies in the field of insect toxicology and insect pest management. Many recent studies have focused on identification and validation of the genes encoding insecticide target proteins, such as acetylcholinesterases, ion channels, Bacillus thuringiensis receptors, and other receptors in the nervous system. RNAi technologies have also been widely applied to reveal the role of genes encoding cytochrome P450 monooxygenases, carboxylesterases, and glutathione S-transferases in insecticide detoxification and resistance. More recently, studies have focused on understanding the mechanism of insecticide-mediated up-regulation of detoxification genes in insects. As RNAi has already shown great potentials for insect pest management, many recent studies have also focused on host-induced gene silencing, in which several RNAi-based transgenic plants have been developed and tested as proof of concept for insect pest management. These studies indicate that RNAi is a valuable tool to address various fundamental questions in insect toxicology and may soon become an effective strategy for insect pest management.

Oral delivery mediated RNA interference of a carboxylesterase gene results in reduced resistance to organophosphorus insecticides in the cotton Aphid, Aphis gossypii Glover

Background

RNA interference (RNAi) is an effective tool to examine the function of individual genes. Carboxylesterases (CarE, EC 3.1.1.1) are known to play significant roles in the metabolism of xenobiotic compounds in many insect species. Previous studies in our laboratory found that CarE expression was up-regulated in Aphis gossypii (Glover) (Hemiptera: Aphididae) adults of both omethoate and malathion resistant strains, indicating the potential involvement of CarE in organophosphorus (OP) insecticide resistance. Functional analysis (RNAi) is therefore warranted to investigate the role of CarE in A. gossypii to OPs resistance.

Result

CarE expression in omethoate resistant individuals of Aphis gossypii was dramatically suppressed following ingestion of dsRNA-CarE. The highest knockdown efficiency (33%) was observed at 72 h after feeding when dsRNA-CarE concentration was 100 ng/µL. The CarE activities from the CarE knockdown aphids were consistent with the correspondingly significant reduction in CarE expression. The CarE activity in the individuals of control aphids was concentrated in the range of 650–900 mOD/per/min, while in the individuals of dsRNA-CarE-fed aphids, the CarE activity was concentrated in the range of 500–800 mOD/per/min. In vitro inhibition experiments also demonstrated that total CarE activity in the CarE knockdown aphids decreased significantly as compared to control aphids. Bioassay results of aphids fed dsRNA-CarE indicated that suppression of CarE expression increased susceptibility to omethoate in individuals of the resistant aphid strains.

Conclusion

The results of this study not only suggest that ingestion of dsRNA through artificial diet could be exploited for functional genomic studies in cotton aphids, but also indicate that CarE can be considered as a major target of organophosphorus insecticide (OPs) resistance in A. gossypii. Further, our results suggest that the CarE would be a propitious target for OPs resistant aphid control, and insect-resistant transgenic plants may be obtained through plant RNAi-mediated silencing of insect CarE expression.

Molecular and functional characterization of cDNAs putatively encoding carboxylesterases from the migratory locust, Locusta migratoria

Carboxylesterases (CarEs) belong to a superfamily of metabolic enzymes encoded by a number of genes and are widely distributed in microbes, plants and animals including insects. These enzymes play important roles in detoxification of insecticides and other xenobiotics, degradation of pheromones, regulation of neurodevelopment, and control of animal development. In this study, we characterized a total of 39 full-length cDNAs putatively encoding different CarEs from the migratory locust, Locusta migratoria, one of the most severe insect pests in many regions of the world, and evaluated the role of four CarE genes in insecticide detoxification. Our phylogenetic analysis grouped the 39 CarEs into five different clades including 20 CarEs in clade A, 3 in D, 13 in E, 1 in F and 2 in I. Four CarE genes (LmCesA3, LmCesA20, LmCesD1, LmCesE1), representing three different clades (A, D and E), were selected for further analyses. The transcripts of the four genes were detectable in all the developmental stages and tissues examined. LmCesA3 and LmCesE1 were mainly expressed in the fat bodies and Malpighian tubules, whereas LmCesA20 and LmCesD1 were predominately expressed in the muscles and hemolymph, respectively. The injection of double-stranded RNA (dsRNA) synthesized from each of the four CarE genes followed by the bioassay with each of four insecticides (chlorpyrifos, malathion, carbaryl and deltamethrin) increased the nymphal mortalities by 37.2 and 28.4% in response to malathion after LmCesA20 and LmCesE1 were silenced, respectively. Thus, we proposed that both LmCesA20 and LmCesE1 played an important role in detoxification of malathion in the locust. These results are expected to help researchers reveal the characteristics of diverse CarEs and assess the risk of insecticide resistance conferred by CarEs in the locust and other insect species.

RNA interference revealed the roles of two carboxylesterase genes in insecticide detoxification in Locusta migratoria

Carboxylesterases (CarEs) play key roles in metabolism of specific hormones and detoxification of dietary and environmental xenobiotics in insects. We sequenced and characterized CarE cDNAs putatively derived from two different genes named LmCesA1 and LmCesA2 from the migratory locust, Locusta migratoria, one of the most important agricultural pests in the world. The full-length cDNAs of LmCesA1 (1892 bp) and LmCesA2 (1643 bp) encode 543 and 501 amino acid residues, respectively. The two deduced CarEs share a characteristic α/β-hydrolase structure, including a catalytic triad composed of Ser-Glu (Asp)-His and a consensus sequence GQSAG, which suggests that both CarEs are biologically active. Phylogenetic analysis grouped both LmCesA1 and LmCesA2 into clade A which has been suggested to be involved in dietary detoxification. Both transcripts were highly expressed in all the nymphal and adult stages, but only slightly expressed in eggs. Analyses of tissue-dependent expression and in situ hybridization revealed that both transcripts were primarily expressed in gastric caeca. RNA interference (RNAi) of LmCesA1 and LmCesA2 followed by a topical application of carbaryl or deltamethrin did not lead to a significantly increased mortality with either insecticide. However, RNAi of LmCesA1 and LmCesA2 increased insect mortalities by 20.9% and 14.5%, respectively, when chlorpyrifos was applied. These results suggest that these genes might not play a significant role in detoxification of carbaryl and deltamethrin but are most likely to be involved in detoxification of chlorpyrifos in L. migratoria.

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.