Alterations of the acetylcholinesterase enzyme in the oriental fruit fly Bactrocera dorsalis are correlated with resistance to the organophosphate insecticide fenitrothion

Reduced insecticide susceptibility of the wheat aphid Sitobion miscanthi after infection by the secondary bacterial symbiont Hamiltonella defensa

BACKGROUND

Bacterial symbionts in insects, especially aphids, have a major influence on host adaptation. We previously showed that infection with the secondary symbiont Hamiltonella defensa increases the fitness of the wheat aphid Sitobion miscanthi, yielding increases in fitness parameters such as adult weight and offspring number. However, whether H. defensa affects the sensitivity of host aphids to insecticides remains unknown.

RESULTS

We tested the effects of H. defensa on host aphid susceptibility to the insecticides chlorpyrifos methyl, imidacloprid, cyantraniliprole and acetamiprid. Our results showed that compared with Hamiltonella-free aphid clones, Hamiltonella-infected aphid clones exhibited lower sensitivity to most of the tested insecticides at low concentrations. Quantitative polymerase chain reaction showed that the density of H. defensa in the infected clones was slightly decreased at 24 h but then sharply increased until the late stage after treatment with the different insecticides. H. defensa in the host aphids was detected by fluorescence in situ hybridization and was localized to the aphid hindgut. Levels of the detoxification enzymes acetylcholinesterase, glutathione transferase and carboxylesterase were significantly higher in Hamiltonella-infected clones than in Hamiltonella-free clones.

CONCLUSIONS

The findings indicated that infection with H. defensa reduced aphid susceptibility to the investigated insecticides at low concentrations, potentially by increasing detoxification enzyme activity in the host. Therefore, symbiont-mediated insecticide resistance should be taken into account when performing resistance-monitoring studies. Studies of symbiont-mediated insecticide resistance may enhance our understanding of the emergence of insecticide resistance in agricultural systems. © 2020 Society of Chemical Industry.

Jaburetox, a urease-derived peptide: Effects on enzymatic pathways of the cockroach Nauphoeta cinerea

Jaburetox is a recombinant peptide derived from one of the Canavalia ensiformis urease isoforms. This peptide induces several toxic effects on insects of different orders, including interference on muscle contractility in cockroaches, modulation of UDP-N-acetylglucosamine pyrophosphorylase (UAP) and nitric oxide synthase (NOS) activities in the central nervous system of triatomines, as well as activation of the immune system in Rhodnius prolixus. When injected, the peptide is lethal for R. prolixus and Triatoma infestans. Here, we evaluated Jaburetox toxicity to Nauphoeta cinerea cockroaches, exploring the effects on the central nervous system through the activities of UAP, NOS, acid phosphatases (ACP), and acetylcholinesterase (AChE). The results indicated that N. cinerea is not susceptible to the lethal effect of the peptide. Moreover, both in vivo and in vitro treatments with Jaburetox inhibited NOS activity, without modifying the protein levels. No alterations on ACP activity were observed. In addition, the enzyme activity of UAP only had its activity affected at 18 hr after injection. The peptide increased the AChE activity, suggesting a mechanism involved in overcoming the toxic effects. In conclusion, our findings indicate that Jaburetox affects the nitrinergic signaling as well as the AChE and UAP activities and establishes N. cinerea as a Jaburetox-resistant model for future comparative studies.

A glutathione S-transferase (BdGSTd9) participates in malathion resistance via directly depleting malathion and its toxic oxide malaoxon in Bactrocera dorsalis (Hendel)

BACKGROUND

The oriental fruit fly, Bactrocera dorsalis (Hendel), is a widespread agricultural pest that has evolved resistance to many commonly used insecticides including malathion. Glutathione S-transferases (GSTs) are multifunctional enzymes that metabolize insecticides directly or indirectly. The specific mechanism used by GSTs to confer malathion resistance in B. dorsalis is unclear.

RESULTS

BdGSTd9 was identified from B. dorsalis and was expressed at twice the level in a malathion-resistant strain (MR) than in a susceptible strain (MS). By using RNAi of BdGSTd9, the toxicity of malathion against MR was increased. Protein modelling and docking of BdGSTd9 with malathion and malaoxon indicated key amino acid residues for direct binding in the active site. In vitro assays with engineered Sf9 cells overexpressing BdGSTd9 demonstrated lower cytotoxicity of malathion. High performance liquid chromatography (HPLC) analysis indicated that malathion could be broken down significantly by BdGSTd9, and it also could deplete the malathion metabolite malaoxon, which possesses a higher toxicity to B. dorsalis. Taken together, the BdGSTd9 of B. dorsalis could not only deplete malathion, but also react with malaoxon and therefore enhance malathion resistance.

CONCLUSION

BdGSTd9 is a component of malathion resistance in B. dorsalis. It acts by depleting both malathion and malaoxon. © 2020 Society of Chemical Industry.

Purification and characterization of acetylcholinesterase in Rhynchophorus ferrugineus (Olivier) (Coleoptera: Curculionidae)

Red palm weevil (RPW), Rhynchophorus ferrugineus, is one of the most destructive pests of cultivated palm trees. The application of synthetic insecticides is currently a main strategy for RPW control. In this study we estimated the distribution of acetylcholinesterase (AChE), as a detoxifying enzyme and the target site of inhibition by insecticides, using ASChI as substrate in different organs of the pest including whole gut, cuticle, fat body, head and haemolymph. The activity ranged from 314.9 to 3868 U in individual organs while the specific activity ranged from 99 to 340.8 U/mg proteins; the cuticle had the highest enzyme level. During larval development, the 11th instar larvae had the highest enzyme content with 5630 U in the cuticle, with a specific activity of 140 U/mg protein. The two major AChE isoenzymes were purified by chromatography on gel filtration and ion exchange columns. They had specific activities of 3504.3 and 2979 U/mg protein, molecular weights of 33 and 54 kDa and activation energies of 8.3 and 4.4 kcal/mol, respectively. Both isoenzymes had monomeric forms, optimum activity at pH 8.0 and 40 °C, were completely inhibited by Hg²⁺ and Cu² and showed similar trends towards the inhibitors eserine, BW284C51 and iso-OMPA. The catalytic properties were compared with those previously recorded for different insect species. This work will pave the way for more studies for improving the understanding of insecticide resistance and developing the field application of synthetic insecticides for controlling R. ferrugineus to ensure successful application.

Acetylcholinesterases from entomopathogenic nematode Heterorhabditid bacteriophora: Susceptibility to insecticides and immunological characteristics

Acetylcholinesterases (AChEs) from the infective juveniles (IJs) of entomopathogenic nematode (EPN) have been investigated with respect to their susceptibility to insecticides and immunological characteristics, aiming at nominating the most compatible insecticide(s) to be used in conjunction with the most insecticide-tolerant EPN strain before incorporation in integrated pest management (IPM) programs. The inhibition kinetics of two purified AChE isoenzymes, AChEAII and AChEBI isolated from Heterorhabditid bacteriophora EM2 strain, by different insecticides revealed that the insensitivity to inhibition by such insecticides could be arranged in a descending order as; methomyl>carbofuran>acetamiprid>oxamyl>malathion. Except for malathion, the insecticides competitively inhibited AChEs with Ki values ranging from 0.1 to 15mM and IC₅₀ values from 1.25 to 23mM. The two AChE isoforms are several folds less sensitive to inhibition by methomyl and carbofuran compared to those previously reported for other insect species. AChEBI was used as an immunogen to raise anti-AChEBI antisera in rabbits. The prepared antisera cross-reacted with AChEs of five different heterorhabditid nematode strains implying the presence of common epitopes shared along all the examined strains. Such studies could aid in the rational selection of the compatible insecticide(s) and the prepared polyclonal anti-AChE antisera would be a valuable immunodiagnostic tool for evaluating the most insecticide-tolerant EPN strain(s) in IPM programs.

Functional characterization of an α-esterase gene involving malathion detoxification in Bactrocera dorsalis (Hendel)

Extensive use of insecticides in many orchards has prompted resistance development in the oriental fruit fly, Bactrocera dorsalis (Hendel). In this study, a laboratory selected strain of B. dorsalis (MR) with a 21-fold higher resistance to malathion was used to examine the resistance mechanisms to this organophosphate insecticide. Carboxylesterase (CarE) was found to be involved in malathion resistance in B. dorsalis from the synergism bioassay by CarE-specific inhibitor triphenylphosphate (TPP). Molecular studies further identified a previously uncharacterized α-esterase gene, BdCarE2, that may function in the development of malathion resistance in B. dorsalis via gene upregulation. This gene is predominantly expressed in the Malpighian tubules, a key insect tissue for detoxification. The transcript levels of BdCarE2 were also compared between the MR and a malathion-susceptible (MS) strain of B. dorsalis, and it was significantly more abundant in the MR strain. No sequence mutation or gene copy changes were detected between the two strains. Functional studies using RNA interference (RNAi)-mediated knockdown of BdCarE2 significantly increased the malathion susceptibility in the adult files. Furthermore, heterologous expression of BdCarE2 combined with cytotoxicity assay in Sf9 cells demonstrated that BdCarE2 could probably detoxify malathion. Taken together, the current study bring new molecular evidence supporting the involvement of CarE-mediated metabolism in resistance development against malathion in B. dorsalis and also provide bases on functional analysis of insect α-esterase associated with insecticide resistance.

The activity of detoxifying enzymes in the infective juveniles of Heterorhabditis bacteriophora strains: Purification and characterization of two acetylcholinesterases

The infectivity and detoxifying enzyme activities including glutathione-S-transferase (GST), acetylcholinesterase (AChE) and carboxylesterase (CaE) are investigated in the infective juveniles (IJs) of six different strains of Heterorhabditis bacteriophora as a biocontrol agent against insect pests. The specific activities ranged from 10.8-29.8 and 50-220units/mg protein for GST and AChE, respectively; and from 24.7-129 and 22.6-77.3units/mg protein for CaE as estimated by P-nitrophenyl and α-naphthyl acetates, respectively. H. bacteriophora EM2 strain has the highest infectivity and the highest enzymatic activities as well. AChE is the predominant detoxifying enzyme that might imply its major role in the detoxification of insecticide(s). The isoenzyme pattern demonstrated two major slow-moving isoforms in all EPN strains examined. Purification of two AChE isoforms, AChEAII and AChEBI, from H. bacteriophora EM2 strain is performed by ammonium sulfate precipitation, gel filtration on Sephacryl S-200 and chromatography on DEAE-Sepharose. AChEAII and AChEBII have specific activities of 1207 and 1560unit/mg protein, native molecular weights of 180 and 68kDa, and are found in dimeric and monomeric forms, respectively. Both isoforms showed optimum activity at pH8.5 and 35°C. AChEBI exhibited higher thermal stability and higher activation energy than AChEAII. The enzymatic activities of purified AChEs are completely inhibited by Hg(+2) and Ni(+2) and greatly enhanced by Mn(+2). The substrate specificity, the relative efficiency of substrates hydrolysis, substrate inhibition and inhibition by BW284C51, but not by iso-OMPA, clearly indicated that they are true AChEs; their properties are compared with those recorded for insects as target hosts for H. bacteriophora EM2.

Discovery of genes related to formothion resistance in oriental fruit fly (Bactrocera dorsalis) by a constrained functional genomics analysis

Artificial selection can provide insights into how insecticide resistance mechanisms evolve in populations. The underlying basis of such phenomena can involve complex interactions of multiple genes, and the resolution of this complexity first necessitates confirmation that specific genes are involved in resistance mechanisms. Here, we used a novel approach invoking a constrained RNA sequencing analysis to refine the discovery of specific genes involved in insecticide resistance. Specifically, for gene discovery, an additional constraint was added to the traditional comparisons of susceptible vs. resistant flies by the incorporation of a line in which insecticide susceptibility was 'recovered' within a resistant line by the removal of insecticide stress. In our analysis, the criterion for the classification of any gene as related to insecticide resistance was based on evidence for differential expression in the resistant line as compared with both the susceptible and recovered lines. The incorporation of this additional constraint reduced the number of differentially expressed genes putatively involved in resistance to 464, compared with more than 1000 that had been identified previously using this same species. In addition, our analysis identified several key genes involved in metabolic detoxification processes that showed up-regulated expression. Furthermore, the involvement of acetylcholinesterase, a known target for modification in insecticide resistance, was associated with three key nonsynonymous amino acid substitutions within our data. In conclusion, the incorporation of an additional constraint using a 'recovered' line for gene discovery provides a higher degree of confidence in genes identified to be involved in insecticide resistance phenomena.

Characterization of Bactrocera dorsalis serine proteases and evidence for their indirect role in insecticide tolerance

The oriental fruit fly Bactrocera dorsalis (Hendel) causes devastating losses to agricultural crops world-wide and is considered to be an economically important pest. Little is known about the digestive enzymes such as serine proteases (SPs) in B. dorsalis, which are important both for energy supply and mitigation of fitness cost associated with insecticide tolerance. In this study, we identified five SP genes in the midgut of B. dorsalis, and the alignments of their deduced amino acid sequences revealed the presence of motifs conserved in the SP superfamily. Phylogenetic analyses with known SPs from other insect species suggested that three of them were trypsin-like proteases. Analyses of the expression profiles among the different developmental stages showed that all five genes were most abundant in larvae than in other stages. When larvae were continuously fed on diet containing 0.33 μg/g β-Cypermethrin, expression of all five genes were upregulated in the midgut but the larval development was delayed. Biochemical assays were consistent with the increased protease activity exhibited by SPs in the midgut after treatment with β-Cypermethrin. Taken together, these findings provide evidence for the hypothesis that enhanced SP activity may play an indirect role in relieving the toxicity stress of insecticide in B. dorsalis.

Toxicological and biochemical characterizations of AChE in phosalone-susceptible and resistant populations of the common pistachio psyllid, Agonoscena pistaciae

The toxicological and biochemical characteristics of acetylcholinesterases (AChE) in nine populations of the common pistachio psyllid, Agonoscena pistaciae Burckhardt and Lauterer (Hemiptera: Psyllidae), were investigated in Kerman Province, Iran. Nine A. pistaciae populations were collected from pistachio orchards, Pistacia vera L. (Sapindales: Anacardiaceae), located in Rafsanjan, Anar, Bam, Kerman, Shahrbabak, Herat, Sirjan, Pariz, and Paghaleh regions of Kerman province. The previous bioassay results showed these populations were susceptible or resistant to phosalone, and the Rafsanjan population was most resistant, with a resistance ratio of 11.3. The specific activity of AChE in the Rafsanjan population was significantly higher than in the susceptible population (Bam). The affinity (K(M)) and hydrolyzing efficiency (Vmax) of AChE on acetylthiocholine iodide, butyrylthiocholine iodide, and propionylthiocholine odide as artificial substrates were clearly lower in the Bam population than that in the Rafsanjan population. These results indicated that the AChE of the Rafsanjan population had lower affinity to these substrates than that of the susceptible population. The higher Vmax value in the Rafsanjan population compared to the susceptible population suggests a possible over expression of AChE in the Rafsanjan population. The in vitro inhibitory effect of several organophosphates and carbamates on AChE of the Rafsanjan and Bam populations was determined. Based on I50, the results showed that the ratios of AChE insensitivity of the resistant to susceptible populations were 23 and 21.7-fold to monocrotophos and phosphamidon, respectively. Whereas, the insensitivity ratios for Rafsanjan population were 0.86, 0.8, 0.78, 0.46, and 0.43 for carbaryl, eserine, propoxur, m-tolyl methyl carbamate, and carbofuran, respectively, suggesting negatively correlated sensitivity to organophosphate-insensitive AChE. Therefore, AChE from the Rafsanjan population showed negatively correlated sensitivity, being insensitive to phosphamidon and monocrotophos and sensitive to N-methyl carbamates.

In silico cloning and annotation of genes involved in the digestion, detoxification and RNA interference mechanism in the midgut of Bactrocera dorsalis [Hendel (Diptera: Tephritidae)]

As the second largest organ in insects, the insect midgut is the major tissue involved in the digestion of food and detoxification of xenobiotics, such as insecticides, and the first barrier and target for oral RNA interference (RNAi). In this study, we performed a midgut-specific transcriptome analysis in the oriental fruit fly, Bactrocera dorsalis, an economically important worldwide pest, with many populations showing high levels of insecticide resistance. Using high-throughput sequencing, 52 838 060 short reads were generated and assembled to 25 236 unigenes with a mean length of 758 bp. Interestingly, 34 unique sequences encoding digestion enzymes were newly described and these included aminopeptidase and trypsin, genes associated with Bacillus thuringiensis resistance and fitness cost. Second, 41 transcripts were annotated to particular detoxification genes such as glutathione S-transferases, carboxylesterases and cytochrome P450s, and the subsequent phylogenetic analysis indicated homology with tissue-specific and insecticide resistance-related genes of Drosophila melanogaster. Third, we identified the genes involved in the mechanism of RNAi and the uptake of double-stranded RNA. The sequences encoding Dicer-2, R2D2, AGO2, and Eater were confirmed, but SID and SR-CI were absent in the midgut transcriptome. In conclusion, the results provide basic molecular information to better understand the mechanisms of food digestion, insecticide resistance and oral RNAi in this important pest insect in agriculture. Specific genes in these systems can be used in the future as potential targets for pest control, for instance, with RNAi technology.

Testing insecticidal activity of novel chemically synthesized siRNA against Plutella xylostella under laboratory and field conditions

BACKGROUND

Over the last 60 years, synthetic chemical pesticides have served as a main tactic in the field of crop protection, but their availability is now declining as a result of the development of insect resistance. Therefore, alternative pest management agents are needed. However, the demonstration of RNAi gene silencing in insects and its successful usage in disrupting the expression of vital genes opened a door to the development of a variety of novel, environmentally sound approaches for insect pest management.

METHODOLOGY/PRINCIPAL FINDINGS

Six small interfering RNAs (siRNAs) were chemically synthesized and modified according to the cDNA sequence of P. xylostella acetylcholine esterase genes AChE1 and AChE2. All of them were formulated and used in insecticide activity screening against P. xylostella. Bioassay data suggested that Si-ace1_003 and Si-ace2_001 at a concentration of 3 µg cm(-2) displayed the best insecticidal activity with 73.7% and 89.0%, mortality, respectively. Additional bioassays were used to obtain the acute lethal concentrations of LC50 and LC90 for Si-ace2_001, which were 53.66 µg/ml and 759.71 µg/ml, respectively. Quantitative Real-time PCR was used to confirm silencing and detected that the transcript levels of P. xylostella AChE2 (PxAChE2) were reduced by 5.7-fold compared to the control group. Consequently, AChE activity was also reduced by 1.7-fold. Finally, effects of the siRNAs on treated plants of Brassica oleracea and Brassica alboglabra were investigated with different siRNA doses. Our results showed that Si-ace2_001 had no negative effects on plant morphology, color and growth of vein under our experimental conditions.

CONCLUSIONS

The most important finding of this study is the discovery that chemically synthesized and modified siRNA corresponding to P. xylostella AChE genes cause significant mortality of the insect both under laboratory and field conditions, which provides a novel strategy to control P. xylostella and to develop bio-pesticides based on the RNA interference technology.

Biochemical and molecular characterisation of acetylcholinesterase in four field populations of Bactrocera dorsalis (Hendel) (Diptera: Tephritidae)

BACKGROUND

The oriental fruit fly, Bactrocera dorsalis, is a major pest that infects fruits and agricultural products worldwide. The latest resistance monitoring of B. dorsalis from mainland China has identified high levels of resistance to insecticides. In this study, the biochemical and molecular characteristics of acetylcholinesterase (AChE) in four field populations of B. dorsalis are investigated.

RESULTS

Among the four populations, the DG population and its purified AChE were found to be the least susceptible to malathion and five inhibitors, whereas the KM population and its purified AChE were the most susceptible. The highest catalytic activity of purified AChE was found for the KM population, and the catalytic activity of the DG population was the lowest. Among developmental stages, the AChE purified from larvae was found to be the most insusceptible to inhibitors, but its catalytic activity was the highest. Sequence analysis of the cDNA encoding AChE showed that some residue differences existed. However, no significant differences in expression levels of the AChE gene among populations and developmental stages were detected.

CONCLUSION

The results suggest that the decrease in susceptibility of B. dorsalis was mainly caused by decrease in AChE activity, and they provide a broad view on the relation between AChE and resistance.

Proposed model for in vitro interaction between fenitrothion and DNA, by using competitive fluorescence, (31)P NMR, (1)H NMR, FT-IR, CD and molecular modeling

In this work we proposed a model for in vitro interaction of fenitrothion (FEN) with calf thymus-DNA by combination of multispectroscopic and two dimensional molecular modeling (ONIOM) methods. The circular dichroism results showed that FEN changes the conformation of B-DNA and caused some changes to C-DNA form. The FT-IR results confirmed a partial intercalation between FEN and edges of all base pairs. The competitive fluorescence, using methylene blue as fluorescence probe, in the presence of increasing amounts of FEN, revealed that FEN is able to release the non-intercalated methylene blue from the DNA. The weak chemical shift and peak broadening of (1)H NMR spectrum of FEN in the presence of DNA confirmed a non-intercalation mode. The (31)P NMR showed that FEN interacts more with DNA via its -NO2 moiety. The ONIOM, based on the hybridization of QM/MM (DFT, 6.31++G (d,p)/UFF) methodology, was also performed by Gaussian 2003 package. The results revealed that the interaction is base sequence dependent, and FEN interacts more with AT base sequences.

Discovery of genes related to insecticide resistance in Bactrocera dorsalis by functional genomic analysis of a de novo assembled transcriptome

Insecticide resistance has recently become a critical concern for control of many insect pest species. Genome sequencing and global quantization of gene expression through analysis of the transcriptome can provide useful information relevant to this challenging problem. The oriental fruit fly, Bactrocera dorsalis, is one of the world's most destructive agricultural pests, and recently it has been used as a target for studies of genetic mechanisms related to insecticide resistance. However, prior to this study, the molecular data available for this species was largely limited to genes identified through homology. To provide a broader pool of gene sequences of potential interest with regard to insecticide resistance, this study uses whole transcriptome analysis developed through de novo assembly of short reads generated by next-generation sequencing (NGS). The transcriptome of B. dorsalis was initially constructed using Illumina's Solexa sequencing technology. Qualified reads were assembled into contigs and potential splicing variants (isotigs). A total of 29,067 isotigs have putative homologues in the non-redundant (nr) protein database from NCBI, and 11,073 of these correspond to distinct D. melanogaster proteins in the RefSeq database. Approximately 5,546 isotigs contain coding sequences that are at least 80% complete and appear to represent B. dorsalis genes. We observed a strong correlation between the completeness of the assembled sequences and the expression intensity of the transcripts. The assembled sequences were also used to identify large numbers of genes potentially belonging to families related to insecticide resistance. A total of 90 P450-, 42 GST-and 37 COE-related genes, representing three major enzyme families involved in insecticide metabolism and resistance, were identified. In addition, 36 isotigs were discovered to contain target site sequences related to four classes of resistance genes. Identified sequence motifs were also analyzed to characterize putative polypeptide translational products and associate them with specific genes and protein functions.

Molecular characterizations of two cytochrome P450 genes encoding CYP6A41 and CYP6EK1 from the oriental fruit fly, Bactrocera dorsalis (Diptera: Tephritidae)

Two P450 genes encoding CYP6A41 and CYP6EK1 were cloned from the oriental fruit fly using polymerase chain reaction (PCR) and rapid amplification of cDNA ends (RACE) techniques. CYP6A41 and CYP6EK1 contained open reading frames of 1,530 and 1,524 nucleotides that encode 510 and 508 amino acid residues, respectively. The putative proteins shared 44% identity with each other. Phylogenetic analysis showed that CYP6A41 and CYP6EK1 were most closely related to Ceratitis capitata CYP6A10 and CYP6A subfamily. Expression patterns of the two genes in different geographical populations (Yunnan, Hainan, Dongguang, and Guangzhou), developmental stages (eggs, larvae, pupae, and adults), and tissues (midguts, fat bodies, and Malpighian tubules) were analyzed by real-time quantitative PCR (RT-qPCR) methods. The results showed that the expression levels of CYP6EK1 were significantly different among the four populations, but were not different for CYP6A41. Both the expressions of CYP6A41 and CYP6EK1 were development specific and had significantly higher levels in the larval stage. The expression of CYP6A41 did not vary among the midgut, fat body, or Malpighian tubules; however, CYP6EK1 expression was higher in the Malpighian tubules. The results suggest that CYP6A41 and CYP6EK1 might be involved in detoxification of xenobiotic compounds that were harmful to larval flies or development. Moreover, high expression of CYP6EK1 in the Malpighian tubules also implied participation in detoxification.

Transcriptome analysis of the oriental fruit fly (Bactrocera dorsalis)

Background

The oriental fruit fly, Bactrocera dorsalis (Hendel), is one of the most economically important pests in the world, causing serious damage to fruit production. However, lack of genetic information on this organism is an obstacle to understanding the mechanisms behind its development and its ability to resist insecticides. Analysis of the B. dorsalis transcriptome and its expression profile data is essential to extending the genetic information resources on this species, providing a shortcut that will support studies on B. dorsalis.

Methodology/Principal Findings

We performed de novo assembly of a transcriptome using short read sequencing technology (Illumina). The results generated 484,628 contigs, 70,640 scaffolds, and 49,804 unigenes. Of those unigenes, 27,455 (55.13%) matched known proteins in the NCBI database, as determined by BLAST search. Clusters of orthologous groups (COG), gene orthology (GO), and the Kyoto Encyclopedia of Genes and Genomes (KEGG) annotations were performed to better understand the functions of these unigenes. Genes related to insecticide resistance were analyzed in additional detail. Digital gene expression (DGE) libraries showed differences in gene expression profiles at different developmental stages (eggs, third-instar larvae, pupae, and adults). To confirm the DGE results, the expression profiles of six randomly selected genes were analyzed.

Conclusion/Significance

This transcriptome greatly improves our genetic understanding of B. dorsalis and makes a huge number of gene sequences available for further study, including both genes of known importance and genes of unknown function. The DGE data provide comprehensive insight into gene expression profiles at different developmental stages. This facilitates the study of the role of each gene in the developmental process and in insecticide resistance.

Purification and biochemical characterization of glutathione S-transferases from four field populations of Bactrocera dorsalis (Hendel) (Diptera: Tephritidae)

Glutathione S-transferases (GSTs) are a group of detoxification enzymes that catalyze the nucleophilic addition of glutathione to a wide variety of endogenous and xenobiotic compounds. In this study, GSTs were purified from four field populations of Bactrocera dorsalis with different insecticide susceptibilities by glutathione-agarose affinity chromatography. The populations were collected from Dongguan (DG) and Guangzhou (GZ) of the Guangdong Province, Haikou of the Hainan province (HN), and Kunming of the Yunnan province (YN), China. Differences in GST characteristics among the four populations were studied using purified enzyme samples through comparative SDS-PAGE, kinetic, and inhibition experiments. The specific activities of the purified enzymes were similar, but the purification yield of the GZ population (31.54%) was the lowest. SDS-PAGE analysis showed only one band at approximately 23 kDa for these four populations. Kinetic analyses showed that the affinities of the purified GSTs from the GZ and YN populations for 1-chloro-2.4-dinitrobenzene (CDNB) were much higher than those of GSTs from the other two populations, whereas the HN population had the highest catalytic capability in terms of V(max) value. The optimum temperature for CDNB conjugation was 37 °C and the optimum pH was 7.5 in all four populations. Inhibition kinetics showed that ethacrynic acid, diethyl maleate, tetraethylthiuram disulfide, curcumin, bromosulfalein, and β-cypermethrin had excellent inhibitory effects on GSTs in the four populations of B. dorsalis, but the low inhibitory effects of malathion and avermectin did not differ between populations. These results suggest that GSTs may have a role in detoxification of β-cypermethrin in B. dorsalis.

Comparative studies of acetylcholinesterase purified from three field populations of Liposcelis entomophila (enderlein) (psocoptera: liposcelididae)

Acetylcholinesterace (AChE) is known to be the major target for organophophate and carbamate insecticides and biomolecular changes to AChE have been demonstrated to be an important mechanism for insecticide resistance in many insect species. In this study, AChE from three field populations of Liposcelis entomophila (Enderlein) (Psocoptera: Liposcelididae) was purified by affinity chromatography and subsequently characterized by its Michaelis-Menten kinetics to determine if detectable changes to AChE have occurred. Bioassays revealed that the potential resistance threat of psocids in Sichuan Province (GH) was greater than either Hubei Province (WH) or Chongqing Municipality (BB). Compared to the other two populations, the WH population possessed the highest specific activity of purified AChE. Kinetic analyses indicated that the purified AChE from GH population expressed a significantly lower affinity to the substrate and a higher catalytic activity toward acetylthiocholine iodide (ATChI) (i.e., higher K(m) and V(max) values) than BB and WH populations. In vitro studies of AChE suggest that five inhibitors (aldicarb, eserine, BW284C51, omethoate, and propoxur) all possess strong inhibitory effects with eserine having the strongest inhibitory effect against purified AChE. According to bimolecular rate constants (k(i)), the purified AChE from GH population was least sensitive to all inhibitors except for omethoate. The differences in AChE among the three populations may be partially attributed to the differences in pesticide application and control practices for psocids among the three locations.

Evaluation of endogenous references for gene expression profiling in different tissues of the oriental fruit fly Bactrocera dorsalis (Diptera: Tephritidae)

Background

Quantitative real-time reverse transcriptase PCR (RT-qPCR) has been widely used for quantification of mRNA as a way to determine key genes involved in different biological processes. For accurate gene quantification analysis, normalization of RT-qPCR data is absolutely essential. To date, normalization is most frequently achieved by the use of internal controls, often referred to as reference genes. However, several studies have shown that the reference genes used for the quantification of mRNA expression can be affected by the experimental set-up or cell type resulting in variation of the expression level of these key genes. Therefore, the evaluation of reference genes is critical for gene expression profiling, which is often neglected in gene expression studies of insects. For this purpose, ten candidate reference genes were investigated in three different tissues (midgut, Malpighian tubules, and fat body) of the oriental fruit fly, Bactrocera dorsalis (Hendel).

Results

Two different programs, geNorm and Normfinder, were used to analyze the data. According to geNorm, α-TUB + ACT5 are the most appropriate reference genes for gene expression profiling across the three different tissues in the female flies, while ACT3 + α-TUB are considered as the best for males. Furthermore, we evaluated the stability of the candidate reference genes to determine the sexual differences in the same tissue. In the midgut and Malpighian tubules, ACT2 + α-TUB are the best choice for both males and females. However, α-TUB + ACT1 are the best pair for fat body. Meanwhile, the results calculated by Normfinder are quite the same as the results with geNorm; α-TUB is always one of the most stable genes in each sample validated by the two programs.

Conclusions

In this study, we validated the suitable reference genes for gene expression profiling in different tissues of B. dorsalis. Moreover, appropriate reference genes were selected out for gene expression profiling of the same tissues taking the sexual differences into consideration. This work not only formed a solid basis for future gene expression study in B. dorsalis, but also will serve as a resource to screen reference genes for gene expression studies in any other insects.

Molecular cloning and characterization of an acetylcholinesterase cDNA in the brown planthopper, Nilaparvata lugens

A full cDNA encoding an acetylcholinesterase (AChE, EC 3.1.1.7) was cloned and characterized from the brown planthopper, Nilaparvata lugens Stål (Hemiptera: Delphacidae). The complete cDNA (2467 bp) contains a 1938-bp open reading frame encoding 646 amino acid residues. The amino acid sequence of the AChE deduced from the cDNA consists of 30 residues for a putative signal peptide and 616 residues for the mature protein with a predicted molecular weight of 69,418. The three residues (Ser242, Glu371, and His485) that putatively form the catalytic triad and the six Cys that form intra-subunit disulfide bonds are completely conserved, and 10 out of the 14 aromatic residues lining the active site gorge of the AChE are also conserved. Northern blot analysis of poly(A)+ RNA showed an approximately 2.6-kb transcript, and Southern blot analysis revealed there likely was just a single copy of this gene in N. lugens. The deduced protein sequence is most similar to AChE of Nephotettix cincticeps with 83% amino acid identity. Phylogenetic analysis constructed with 45 AChEs from 30 species showed that the deduced N. lugens AChE formed a cluster with the other 8 insect AChE2s. Additionally, the hypervariable region and amino acids specific to insect AChE2 also existed in the AChE of N. lugens. The results revealed that the AChE cDNA cloned in this work belongs to insect AChE2 subgroup, which is orthologous to Drosophila AChE. Comparison of the AChEs between the susceptible and resistant strains revealed a point mutation, Gly185Ser, is likely responsible for the insensitivity of the AChE to methamidopho in the resistant strain.