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

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.

Comparative analyses of transcriptional responses of Dectes texanus LeConte (Coleoptera: Cerambycidae) larvae fed on three different host plants and artificial diet

Dectes texanus is an important coleopteran pest of soybeans and cultivated sunflowers in the Midwestern United States that causes yield losses by girdling stems of their host plants. Although sunflower and giant ragweed are primary hosts of D. texanus, they began colonizing soybeans approximately 50 years ago and no reliable management method has been established to prevent or reduce losses by this pest. To identify genes putatively involved when feeding soybean, we compared gene expression of D. texanus third-instar larvae fed soybean to those fed sunflower, giant ragweed, or artificial diet. Dectes texanus larvae differentially expressed 514 unigenes when fed on soybean compared to those fed the other diet treatments. Enrichment analyses of gene ontology terms from up-regulated unigenes in soybean-fed larvae compared to those fed both primary hosts highlighted unigenes involved in oxidoreductase and polygalacturonase activities. Cytochrome P450s, carboxylesterases, major facilitator superfamily transporters, lipocalins, apolipoproteins, glycoside hydrolases 1 and 28, and lytic monooxygenases were among the most commonly up-regulated unigenes in soybean-fed larvae compared to those fed their primary hosts. These results suggest that D. texanus larvae differentially expressed unigenes involved in biotransformation of allelochemicals, digestion of plant cell walls and transport of small solutes and lipids when feeding in soybean.

RNA Sequencing Analysis of Metopolophium dirhodum (Walker) (Hemiptera: Aphididae) Reveals the Mechanism Underlying Insecticide Resistance

Xinjiang (XJ) and Ningxia (NX) provinces are important agricultural regions in western China. Aphids are one kind of the most devastating pests in both the provinces. Aphids are typical phloem-feeding insects distributed worldwide and can severely damage crops. In this study, two representative Metopolophium dirhodum (Walker) (Hemiptera: Aphididae) populations were collected from the typical agricultural regions of XJ and NX, respectively for a high-throughput transcriptome sequencing analysis. A total of 5,265 differentially expressed genes (DEGs) were identified. The functional annotation of DEGs and the identification of enriched pathways indicated many of the DEGs are involved in processes related to energy metabolism, development, and insecticide resistance. Furthermore, an investigation of insecticide toxicity revealed the NX population is more resistant to insecticide treatments than the XJ population. Thus, the transcriptome data generated in present study can be used for functional gene characterization relevant to aphid development, metabolism, environmental adaptation, and insecticide resistance.

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.

Transcriptome Sequencing of the Striped Cucumber Beetle, Acalymma vittatum (F.), Reveals Numerous Sex-Specific Transcripts and Xenobiotic Detoxification Genes

Acalymma vittatum (F.), the striped cucumber beetle, is an important pest of cucurbit crops in the contintental United States, damaging plants through both direct feeding and vectoring of a bacterial wilt pathogen. Besides providing basic biological knowledge, biosequence data for A. vittatum would be useful towards the development of molecular biopesticides to complement existing population control methods. However, no such datasets currently exist. In this study, three biological replicates apiece of male and female adult insects were sequenced and assembled into a set of 630,139 transcripts (of which 232,899 exhibited hits to one or more sequences in NCBI NR). Quantitative analyses identified 2898 genes differentially expressed across the male–female divide, and qualitative analyses characterized the insect’s resistome, comprising the glutathione S-transferase, carboxylesterase, and cytochrome P450 monooxygenase families of xenobiotic detoxification genes. In summary, these data provide useful insights into genes associated with sex differentiation and this beetle’s innate genetic capacity to develop resistance to synthetic pesticides; furthermore, these genes may serve as useful targets for potential use in molecular-based biocontrol technologies.

Functional Characterization of an α-Esterase Gene Associated with Malathion Detoxification in Bradysia odoriphaga

Carboxylesterases (CarEs) are a multigene superfamily of metabolic enzymes involved in metabolic detoxification of xenobiotics. In this study, an α-esterase gene (BoαE1) was identified from Bradysia odoriphaga. Phylogenetic analysis classified BoαE1 into the α-esterase clade. Developmental expression analysis indicated that BoαE1 was significantly expressed in the second to fourth larval stages. Tissue-specific expression analysis indicated that BoαE1 was highly expressed in the larval midgut. After exposure to LC₃₀ of malathion, the CarE activity of B. odoriphaga was induced and the transcriptional level of BoαE1 was significantly up-regulated. Silencing of BoαE1 significantly increased the susceptibility of B. odoriphaga larvae to malathion. Inhibition assays in vitro indicated that malathion significantly inhibited BoαE1 activity. GC-MS assay showed that BoαE1 possesses hydrolase activity toward malathion and participates in the detoxification of malathion. These results strongly suggest that BoαE1 plays a crucial role in detoxification of malathion in B. odoriphaga.

De novo transcriptome analysis and identification of genes associated with immunity, detoxification and energy metabolism from the fat body of the tephritid gall fly, Procecidochares utilis

The fat body, a multifunctional organ analogous to the liver and fat tissue of vertebrates, plays an important role in insect life cycles. The fat body is involved in protein storage, energy metabolism, elimination of xenobiotics, and production of immunity regulator-like proteins. However, the molecular mechanism of the fat body's physiological functions in the tephritid stem gall-forming fly, Procecidochares utilis, are still unknown. In this study, we performed transcriptome analysis of the fat body of P. utilis using Illumina sequencing technology. In total, 3.71 G of clean reads were obtained and assembled into 30,559 unigenes, with an average length of 539 bp. Among those unigenes, 21,439 (70.16%) were annotated based on sequence similarity to proteins in NCBI's non-redundant protein sequence database (Nr). Sequences were also compared to NCBI's non-redundant nucleotide sequence database (Nt), a manually curated and reviewed protein sequence database (SwissProt), and KEGG and gene ontology annotations were applied to better understand the functions of these unigenes. A comparative analysis was performed to identify unigenes related to detoxification, immunity and energy metabolism. Many unigenes involved in detoxification were identified, including 50 unigenes of putative cytochrome P450s (P450s), 18 of glutathione S-transferases (GSTs), 35 of carboxylesterases (CarEs) and 26 of ATP-binding cassette (ABC) transporters. Many unigenes related to immunity were identified, including 17 putative serpin genes, five peptidoglycan recognition proteins (PGRPs) and four lysozyme genes. In addition, unigenes potentially involved in energy metabolism, including 18 lipase genes, five fatty acid synthase (FAS) genes and six elongases of very long chain fatty acid (ELOVL) genes, were identified. This transcriptome improves our genetic understanding of P. utilis and the identification of a numerous transcripts in the fat body of P. utilis offer a series of valuable molecular resources for future studies on the functions of these genes.

Latrophilin mediates insecticides susceptibility and fecundity through two carboxylesterases, esterase4 and esterase6, in Tribolium castaneum

Latrophilin (LPH) is known as an adhesion G-protein-coupled receptor which involved in multiple physiological processes in organisms. Previous studies showed that lph not only involved the susceptibility to anticholinesterase insecticides but also affected fecundity in Tribolium castaneum. However, its regulatory mechanisms in these biological processes are still not clear. Here, we identified two potential downstream carboxylesterase (cce) genes of Tclph, esterase4 and esterase6, and further characterized their interactions with Tclph. After treatment of T. castaneum larvae with carbofuran or dichlorvos insecticides, the transcript levels of Tcest4 and Tcest6 were significantly induced from 12 to 72 h. RNAi against Tcest4 or Tcest6 led to the higher mortality compared with the controls after the insecticides treatment, suggesting that these two genes play a vital role in detoxification of insecticides in T. castaneum. Furthermore, with insecticides exposure to Tclph knockdown beetles, the expression of Tcest4 was upregulated but Tcest6 was downregulated, indicating that beetles existed a compensatory response against the insecticides. Additionally, RNAi of Tcest6 resulted in 43% reductions in female egg laying and completely inhibited egg hatching, which showed the similar phenotype as that of Tclph knockdown. These results indicated that Tclph affected fecundity by positively regulating Tcest6 expression. Our findings will provide a new insight into the molecular mechanisms of Tclph involved in physiological functions in T. castaneum.

Functional characterization of the Hyles euphorbiae hawkmoth transcriptome reveals strong expression of phorbol ester detoxification and seasonal cold hardiness genes

Background

The European spurge hawkmoth, Hyles euphorbiae (Lepidoptera, Sphingidae), has been intensively studied as a model organism for insect chemical ecology, cold hardiness and evolution of species delineation. To understand species isolation mechanisms at a molecular level, this study aims at determining genetic factors underlying two adaptive ecological trait candidates, phorbol ester (TPA) detoxification and seasonal cold acclimation.

Method

A draft transcriptome of H. euphorbiae was generated using Illumina sequencing, providing the first genomic resource for the hawkmoth subfamily Macroglossinae. RNA expression levels in tissues of experimental TPA feeding larvae and cooled pupae was compared to levels in control larvae and pupae using 26 bp RNA sequence tag libraries (DeepSuperSAGE). Differential gene expression was assessed by homology searches of the tags in the transcriptome.

Results

In total, 389 and 605 differentially expressed transcripts for detoxification and cold hardiness, respectively, could be identified and annotated with proteins. The majority (22 of 28) of differentially expressed detox transcripts of the four 'drug metabolism' enzyme groups (cytochrome P450 (CYP), carboxylesterases (CES), glutathione S-transferases (GST) and lipases) are up-regulated. Triacylglycerol lipase was significantly over proportionally annotated among up-regulated detox transcripts. We record several up-regulated lipases, GSTe2, two CESs, CYP9A21, CYP6BD6 and CYP9A17 as candidate genes for further H. euphorbiae TPA detoxification analyses. Differential gene expression of the cold acclimation treatment is marked by metabolic depression with enriched Gene Ontology terms among down-regulated transcripts almost exclusively comprising metabolism, aerobic respiration and dissimilative functions. Down-regulated transcripts include energy expensive respiratory proteins like NADH dehydrogenase, cytochrome oxidase and ATP synthase. Gene expression patterns show shifts in carbohydrate metabolism towards cryoprotectant production. The Glycolysis enzymes, G1Pase, A1e, Gpi and an Akr isoform are up-regulated. Glycerol, an osmolyte which lowers the body liquid supercooling point, appears to be the predominant polyol cryoprotectant in H. euphorbiae diapause pupae. Several protein candidates involved in glucose, glycerol, myo-inositol and potentially sorbitol and trehalose synthesis were identified.

Conclusions

A majority of differently expressed transcripts unique for either detoxification or cold hardiness indicates highly specialized functional adaptation which may have evolved from general cell metabolism and stress response.The transcriptome and extracted candidate biomarkers provide a basis for further gene expression studies of physiological processes and adaptive traits in H. euphorbiae.

Functional characterization of carboxylesterase gene mutations involved in Aphis gossypii resistance to organophosphate insecticides

Carboxylesterases (CarEs) play an important role in detoxifying insecticides in insects. Over-expression and structural modification of CarEs have been implicated in the development of organophosphate (OP) insecticide resistance in insects. A previous study identified four nonsynonymous mutations (resulting in four amino acid residue substitutions) in the open reading frame of the carboxylesterase gene of resistant cotton aphids compared to the omethoate susceptible strain, which has possibly influenced the development of resistance to omethoate (a systemic OP insecticide). The current study further characterized the function of these mutations, both alone and in combination, in the hydrolysis of OP insecticides. The metabolism results suggest that the combination of four mutations, mainly existing in the laboratory-selected OP-resistant cotton aphid population, increased the OP hydrolase activity (approximately twofold) at the cost of detectable carboxylesterase activity. The functional studies of single or multiple mutations suggest the positive effect of H104R, A128V and T333P on the acquisition of OP hydrolase activity, especially the combination of H104R with A128V or T333P. K484R substitution decreased both the OP hydrolase activity and the CarE activity, indicating that this mutation primarily drives the negative effect on the acquisition of OP hydrolase activity amongst these four mutations in the resistant strain. The modelling and docking results are basically consistent with the metabolic results, which strongly suggest that the structural gene modification is the molecular basis for the OP resistance in this laboratory-selected cotton aphid strain.

Functional characterization of BdB1, a well-conserved carboxylesterase among tephritid fruit flies associated with malathion resistance in Bactrocera dorsalis (Hendel)

There are many evidences that insect carboxylesterase possess important physiological roles in xenobiotic metabolism and are implicated in the detoxification of organophosphate (OP) insecticides. Despite the ongoing resistance development in the oriental fruit fly, Bactrocera dorsalis (Hendel), the molecular basis of carboxylesterase and its ability to confer OP resistance remain largely obscure. This study was initiated to provide a better understanding of carboxylesterase-mediated resistance mechanism in a tephritid pest fly. Here, we narrow this research gap by demonstrating a well-conserved esterase B1 gene, BdB1, mediates malathion resistance development via gene upregulation with the use of a laboratory selected malathion-resistant strain (MR) of B. dorsalis. No sequence mutation of BdB1 was detected between MR and the susceptible strain (MS) of B. dorsalis. BdB1 is predominantly expressed in the midgut, a key insect tissue for detoxification. As compared with transcripts in MS, BdB1 was significantly more abundant in multiple tissues in the MR. RNA interference (RNAi)-mediated knockdown of BdB1 significantly increased malathion susceptibility. Furthermore, heterologous expression along with cytotoxicity assay revealed BdB1 could probably have the function of malathion detoxification.

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.

Overexpression of two α-esterase genes mediates metabolic resistance to malathion in the oriental fruit fly, Bactrocera dorsalis (Hendel)

Esterase has been reported to be involved in malathion resistance in the oriental fruit fly, Bactrocera dorsalis (Hendel). However, the underlying molecular mechanism of the esterase-mediated resistance remains largely unknown in this species. Here, with the use of a strain selected for malathion resistance in the laboratory (MR), we found that two overexpressed α-esterase genes, namely BdCarE4 and BdCarE6, predominant in the adult midgut and fat body, function in conferring malathion resistance in B. dorsalis. Notably, these two genes were found to be mostly close to the esterase E3, which are usually implicated in detoxifying organophosphate insecticides. The transcript levels of BdCarE4 and BdCarE6 were investigated and compared between the MR and a susceptible (MS) strain of B. dorsalis. Both genes were significantly up-regulated in the MR strain, which was consistent with the enhanced esterase activity in the MR strain. However, no changes in either the coding sequence or gene copy number were observed between the two strains. Subsequently, heterologous expression combined with cytotoxicity assay in Sf9 cells demonstrated that BdCarE4 and BdCarE6 can probably detoxify malathion. Furthermore, RNA interference-mediated knockdown of each of these two genes significantly increased malathion susceptibility in the MR strain adults. In conclusion, these results expand our molecular understanding of the important role of α-esterases during the development of resistance to organophosphorous insecticides in B. dorsalis.

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.

Key enzymes and proteins of crop insects as candidate for RNAi based gene silencing

RNA interference (RNAi) is a mechanism of homology dependent gene silencing present in plants and animals. It operates through 21-24 nucleotides small RNAs which are processed through a set of core enzymatic machinery that involves Dicer and Argonaute proteins. In recent past, the technology has been well appreciated toward the control of plant pathogens and insects through suppression of key genes/proteins of infecting organisms. The genes encoding key enzymes/proteins with the great potential for developing an effective insect control by RNAi approach are actylcholinesterase, cytochrome P450 enzymes, amino peptidase N, allatostatin, allatotropin, tryptophan oxygenase, arginine kinase, vacuolar ATPase, chitin synthase, glutathione-S-transferase, catalase, trehalose phosphate synthase, vitellogenin, hydroxy-3-methylglutaryl coenzyme A reductase, and hormone receptor genes. Through various studies, it is demonstrated that RNAi is a reliable molecular tool which offers great promises in meeting the challenges imposed by crop insects with careful selection of key enzymes/proteins. Utilization of RNAi tool to target some of these key proteins of crop insects through various approaches is described here. The major challenges of RNAi based insect control such as identifying potential targets, delivery methods of silencing trigger, off target effects, and complexity of insect biology are very well illustrated. Further, required efforts to address these challenges are also discussed.

Systemic RNA interference in locusts: reverse genetics and possibilities for locust pest control

RNA interference (RNAi) is a biological process triggered by double stranded (ds)RNA that results in sequence-dependent mRNA degradation. Because of its high specificity, this post-transcriptional gene silencing mechanism is a widely used tool for reverse genetics in several insect species. In particular, locusts possess a very robust and sensitive RNAi response that has already been exploited to investigate a diverse range of important physiological processes. These orthopteran insects constitute important model organisms in several areas of entomology, but they can also become voracious swarming pests that threaten the agricultural production in large parts of the world. In comparison to the widely applied chemical insecticides, the RNAi-technology could contribute to the development of a novel generation of insecticides, with high species-specificity. In this article, we discuss the potential of the RNAi-technology in loss of function studies in locusts, as well as to control locust populations.