Mining genes involved in insecticide resistance of Liposcelis bostrychophila Badonnel by transcriptome and expression profile analysis

Antioxidant Enzymes and Heat Shock Protein Genes from Liposcelis bostrychophila Are Involved in Stress Defense upon Heat Shock

Simple Summary

Liposcelis bostrychophila is one of the most serious pests of stored commodities among the psocids. Controlling psocids mainly relies on chemical insecticides and heat stress. In fact, L. bostrychophila has developed high levels of resistance or tolerance to heat treatment in grain storage systems. In this study, we evaluated the changes in malondialdehyde (MDA) concentration after different high temperatures. The result showed that MDA is increased slightly overall, but a drastic increase is detected at 42.5 °C for exposure of different times. To further explore the principles of L. bostrychophila in response to heat stress, we tested the changes of superoxide dismutase (SOD), catalase (CAT), peroxidases (POD) and glutathione-S-transferases (GST) activities under different heat treatments and identified four inducible LbHsp70 genes and one LbHsp110 gene. Enzyme activities and transcript levels changed drastically after different heat treatments. These findings contribute to our understanding of the mechanism of L. bostrychophila responding to heat stress and provide baseline information for further understanding the excellent targets of L. bostrychophila.

Abstract

Psocids are a new risk for global food security and safety because they are significant worldwide pests of stored products. Among these psocids, Liposcelis bostrychophila has developed high levels of resistance or tolerance to heat treatment in grain storage systems, and thus has led to investigation of molecular mechanisms underlying heat tolerance in this pest. In this study, the time-related effects of thermal stress treatments at relatively high temperatures on the activities of antioxidant enzymes, including superoxide dismutase (SOD), catalase (CAT), peroxidases (POD), glutathione-S-transferases (GST) and malondialdehyde (MDA), of L. bostrychophila were determined. Thermal stress resulted that L. bostrychophila had a significantly higher MDA concentration at 42.5 °C, which indicated that the heat stress increased lipid peroxidation (LPO) contents and oxidative stress in this psocid pest. Heat stress also resulted in significant elevation of SOD, CAT and GST activities but decreased POD activity. Our data indicates that different antioxidant enzymes contribute to defense mechanisms, counteracting oxidative damage in varying levels. POD play minor roles in scavenging deleterious LPO, while enhanced SOD, CAT and GST activities in response to thermal stress likely play a more important role against oxidative damage. Here, we firstly identified five LbHsps (four LbHsp70s and one LbHsp110) from psocids, and most of these LbHsps (except LbHsp70-1) are highly expressed at fourth instar nymph and adults, and LbHsp70-1 likely presents as a cognate form of HSP due to its non-significant changes of expression. Most LbHsp70s (except LbHsp70-4) are significantly induced at moderate high temperatures (<40 °C) and decreased at extreme high temperatures (40–45 °C), but LbHsp110-1 can be significantly induced at all high temperatures. Results of this study suggest that the LbHsp70s and LbHsp110 genes are involved in tolerance to thermal stress in L. bostrychophila, and antioxidant enzymes and heat shock proteins may be coordinately involved in the tolerance to thermal stress in psocids.

Identification of transcriptional responsive genes to acetic acid, ethanol, and 2-phenylethanol exposure in Drosophila melanogaster

The fruit fly, Drosophila melanogaster, is predominantly found in overripe, rotten, fermenting, or decaying fruits and is constantly exposed to chemical stressors such as acetic acid, ethanol, and 2-phenylethanol. D. melanogaster has been employed as a model system for studying the molecular bases of various types of chemical-induced tolerance. Expression profiling using Illumina sequencing has been performed for identifying changes in gene expression that may be associated with evolutionary adaptation to exposure of acetic acid, ethanol, and 2-phenylethanol. We identified a total of 457 differentially expressed genes that may affect sensitivity or tolerance to three chemicals in the chemical treatment group as opposed to the control group. Gene-set enrichment analysis revealed that the genes involved in metabolism, multicellular organism reproduction, olfaction, regulation of signal transduction, and stress tolerance were over-represented in response to chemical exposure. Furthermore, we also detected a coordinated upregulation of genes in the Toll- and Imd-signaling pathways after the chemical exposure. Quantitative reverse transcription PCR analysis revealed that the expression levels of nine genes within the set of genes identified by RNA sequencing were up- or downregulated owing to chemical exposure. Taken together, our data suggest that such differentially expressed genes are coordinately affected by chemical exposure. Transcriptional analyses after exposure of D. melanogaster with three chemicals provide unique insights into subsequent functional studies on the mechanisms underlying the evolutionary adaptation of insect species to environmental chemical stressors.

Characterization of Esterase Genes Involving Malathion Detoxification and Establishment of an RNA Interference Method in Liposcelis bostrychophila

Esterases (ESTs) play important roles in metabolizing various physiologically endogenous and exogenous compounds, and various environmental xenobiotics in insects. The psocid, Liposcelis bostrychophila is a major pest of stored products worldwide and rapidly develops resistance to commonly insecticides. However, the involvement of ESTs in insecticide metabolization and the application of RNAi approach in psocids have not been well elucidated. In this study, we characterized four LbEST genes and investigated the transcriptional levels of these genes at different developmental stages and under different insecticides exposures to assess their potential roles in response to insecticides. The four LbESTs contain a catalytic triad (Ser-His-Glu) linked to an oxyanion hole and acyl pocket involved in substrate stabilization during its hydrolysis. Synergism observed with the esterase-inhibitor DEF suggests the involvement of esterases in malathion detoxification. LbESTs were expressed during the whole of developmental stages, but predominant abundance in the first nymphal instar and adult stage. The mRNA level of three LbEST genes (except for LbEST4) was induced (1.29- to 5.60 fold) in response to malathion or deltamethrin exposures, indicating that these esterases are involved in the detoxification process. Silencing of LbEST1, LbEST2 or LbEST3 through dsRNA feeding led to a higher mortality of psocids upon the malathion treatment compared to controls (1.83 to 2.69-fold), demonstrating that these esterase genes play roles in malathion detoxification in L. bostrychophila. Our study provides new evidence for understanding of the function and regulation mechanism of esterases in L. bostrychophila in insecticide detoxification. The current study also suggests that the present RNAi method could be applied for gene functional studies in psocids.

Lip b 1 is a novel allergenic protein isolated from the booklouse, Liposcelis bostrychophila

BACKGROUND

Booklice, belonging to the order Psocoptera, are small household insect pests that are distributed worldwide. Liposcelis bostrychophila, a common home-inhabiting species of booklouse, infests old books, sheets of paper, and stored food. Recent entomological and serological studies demonstrated that L. bostrychophila accounted for the majority of detectable insects in house dust and could be a potent inducer of respiratory allergy. Our recent proteomic analysis identified a potent allergenic protein from L. bostrychophila, designated Lip b 1, and determined its partial amino acid sequences.

METHODS

Cloning of cDNAs for Lip b 1 was performed by large-scale transcriptome analysis (RNA-seq) and subsequent reverse transcription polymerase chain reaction. The full-length amino acid sequences deduced from Lip b 1 cDNAs were bioinformatically analyzed. The recombinant proteins of glutathione S-transferase (GST)-fused Lip b 1 were analyzed by Western blot and enzyme-linked immunosorbent assay.

RESULTS

Lip b 1 cDNAs encoding two types of 254-amino acid proteins were cloned. The clones shared 87% identity, and the deduced molecular weights and isoelectric points were consistent with those determined in our previous study. The two types of Lip b 1 proteins in the GST-fused form were similarly reactive with sera from allergic patients sensitized with L. bostrychophila.

CONCLUSIONS

Lip b 1 is a novel protein possibly causing booklouse allergy.

Patterns of Gene Expression in Western Corn Rootworm (Diabrotica virgifera virgifera) Neonates, Challenged with Cry34Ab1, Cry35Ab1 and Cry34/35Ab1, Based on Next-Generation Sequencing

With Next Generation Sequencing technologies, high-throughput RNA sequencing (RNAseq) was conducted to examine gene expression in neonates of Diabrotica virgifera virgifera (LeConte) (Western Corn Rootworm, WCR) challenged with individual proteins of the binary Bacillus thuringiensis insecticidal proteins, Cry34Ab1 and Cry35Ab1, and the combination of Cry34/Cry35Ab1, which together are active against rootworm larvae. Integrated results of three different statistical comparisons identified 114 and 1300 differentially expressed transcripts (DETs) in the Cry34Ab1 and Cry34/35Ab1 treatment, respectively, as compared to the control. No DETs were identified in the Cry35Ab1 treatment. Putative Bt binding receptors previously identified in other insect species were not identified in DETs in this study. The majority of DETs (75% with Cry34Ab1 and 68.3% with Cry34/35Ab1 treatments) had no significant hits in the NCBI nr database. In addition, 92 DETs were shared between Cry34Ab1 and Cry34/35Ab1 treatments. Further analysis revealed that the most abundant DETs in both Cry34Ab1 and Cry34/35Ab1 treatments were associated with binding and catalytic activity. Results from this study confirmed the nature of these binary toxins against WCR larvae and provide a fundamental profile of expression pattern of genes in response to challenge of the Cry34/35Ab1 toxin, which may provide insight into potential resistance mechanisms.

DIRProt: a computational approach for discriminating insecticide resistant proteins from non-resistant proteins

BACKGROUND

Insecticide resistance is a major challenge for the control program of insect pests in the fields of crop protection, human and animal health etc. Resistance to different insecticides is conferred by the proteins encoded from certain class of genes of the insects. To distinguish the insecticide resistant proteins from non-resistant proteins, no computational tool is available till date. Thus, development of such a computational tool will be helpful in predicting the insecticide resistant proteins, which can be targeted for developing appropriate insecticides.

RESULTS

Five different sets of feature viz., amino acid composition (AAC), di-peptide composition (DPC), pseudo amino acid composition (PAAC), composition-transition-distribution (CTD) and auto-correlation function (ACF) were used to map the protein sequences into numeric feature vectors. The encoded numeric vectors were then used as input in support vector machine (SVM) for classification of insecticide resistant and non-resistant proteins. Higher accuracies were obtained under RBF kernel than that of other kernels. Further, accuracies were observed to be higher for DPC feature set as compared to others. The proposed approach achieved an overall accuracy of >90% in discriminating resistant from non-resistant proteins. Further, the two classes of resistant proteins i.e., detoxification-based and target-based were discriminated from non-resistant proteins with >95% accuracy. Besides, >95% accuracy was also observed for discrimination of proteins involved in detoxification- and target-based resistance mechanisms. The proposed approach not only outperformed Blastp, PSI-Blast and Delta-Blast algorithms, but also achieved >92% accuracy while assessed using an independent dataset of 75 insecticide resistant proteins.

CONCLUSIONS

This paper presents the first computational approach for discriminating the insecticide resistant proteins from non-resistant proteins. Based on the proposed approach, an online prediction server DIRProt has also been developed for computational prediction of insecticide resistant proteins, which is accessible at http://cabgrid.res.in:8080/dirprot/ . The proposed approach is believed to supplement the efforts needed to develop dynamic insecticides in wet-lab by targeting the insecticide resistant proteins.

The choreography of the chemical defensome response to insecticide stress: insights into the Anopheles stephensi transcriptome using RNA-Seq

Animals respond to chemical stress with an array of gene families and pathways termed “chemical defensome”. In arthropods, despite many defensome genes have been detected, how their activation is arranged during toxic exposure remains poorly understood. Here, we sequenced the transcriptome of Anopheles stephensi larvae exposed for six, 24 and 48 hours to the LD₅₀ dose of the insecticide permethrin to monitor transcriptional changes of defensome genes across time. A total of 177 genes involved in insecticide defense were differentially expressed (DE) in at least one time-point, including genes encoding for Phase 0, I, II, III and antioxidant enzymes and for Heat Shock and Cuticular Proteins. Three major patterns emerged throughout time. First, most of DE genes were down-regulated at all time-points, suggesting a reallocation of energetic resources during insecticide stress. Second, single genes and clusters of genes turn off and on from six to 48 hours of treatment, showing a modulated response across time. Third, the number of up-regulated genes peaked at six hours and then decreased during exposure. Our results give a first picture of how defensome gene families respond against toxicants and provide a valuable resource for understanding how defensome genes work together during insecticide stress.

De novo transcriptome and expression profile analyses of the Asian corn borer (Ostrinia furnacalis) reveals relevant flubendiamide response genes

Background

The Asian corn borer (ACB), Ostrinia furnacalis (Guenée), has become the most damaging insect pest of corn in Asia. However, the lack of genome or transcriptome information heavily hinders our further understanding of ACB in every aspect at a molecular level and on a genome-wide scale. Here, we used the Ion Torrent Personal Genome Machine (PGM) Sequencer to explore the ACB transcriptome and to identify relevant genes in response to flubendiamide, showing high selective activity against ACB.

Results

We obtained 35,430 unigenes, with an average length of 716 bp, representing a dramatic expansion of existing cDNA sequences available for ACB. These sequences were annotated with Non-redundant Protein (Nr), Gene Ontology (GO), Clusters of Orthologous Groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) to better understand their functions. A total of 31 cytochrome P450 monooxygenases (P450s), 27 carboxyl/cholinesterases (CCEs) and 19 glutathione S-transferases (GSTs) were manually curated to construct phylogenetic trees, and 25 unigenes encoding target proteins (acetylcholinesterase, nicotinic acetylcholine receptor, gamma-aminobutyric acid receptor, glutamate-gated chloride channel, voltage-gated sodium channel and ryanodine receptor) were identified. In addition, we compared and validated the differentially expressed unigenes upon flubendiamide treatment, revealing that the genes for detoxification enzymes (P450s and esterase), calcium signaling pathways and muscle control pathways (twitchin and tropomyosin), immunoglobulin (hemolin), chemosensory protein and heat shock protein 70 were significantly overexpressed in response to flubendiamide, while the genes for cuticular protein, protease and oxidoreductase showed much lower expression levels.

Conclusion

The obtained transcriptome information provides large genomic resources available for further studies of ACB. The differentially expressed gene data will elucidate the molecular mechanisms of ACB in response to the novel diamide insecticide, flubendiamide. In particular, these findings will facilitate the identification of the genes involved in insecticide resistance and the development of new compounds to control the ACB.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3431-6) contains supplementary material, which is available to authorized users.

How heterogeneous is the involvement of ABC transporters against insecticides?

Understanding the molecular mechanisms underlying cellular defense against xenobiotic compounds is a main research issue in medical and veterinary entomology, as insecticide/acaricide resistance is a major threat in the control of arthropods. ABC transporters are recognized as a component of the detoxifying mechanism in arthropods. We investigated the possible involvement of ABC transporters in defense to the organophosphate insecticide temephos in the malarial vector Anopheles stephensi. We performed bioassays on larvae of An. stephensi, using insecticide alone and in combination with ABC-transporter inhibitors, to assess synergism between these compounds. Next, we investigated the expression profiles of six ABC transporter genes in larvae exposed to temephos. Surprisingly, neither bioassays nor gene expression analyses provided any evidence for a major role of ABC transporters in defense against temephos in An. stephensi. We thus decided to review existing literature to generate a record of other studies that failed to reveal a role for ABC transporters against particular insecticides/acaricides. A review of the scientific literature led to the recovery of 569 papers about ABC transporters; among these, 50 involved arthropods, and 10 reported negative results. Our study on An. stephensi and accompanying literature review highlight the heterogeneity that exists in ABC transporter involvement in defense/resistance mechanisms in arthropods.

Transcriptome profiling of the testis reveals genes involved in spermatogenesis and marker discovery in the oriental fruit fly, Bactrocera dorsalis

The testis is a highly specialized tissue that plays a vital role in ensuring fertility by producing spermatozoa, which are transferred to the female during mating. Spermatogenesis is a complex process, resulting in the production of mature sperm, and involves significant structural and biochemical changes in the seminiferous epithelium of the adult testis. The identification of genes involved in spermatogenesis of Bactrocera dorsalis (Hendel) is critical for a better understanding of its reproductive development. In this study, we constructed a cDNA library of testes from male B. dorsalis adults at different ages, and performed de novo transcriptome sequencing to produce a comprehensive transcript data set, using Illumina sequencing technology. The analysis yielded 52 016 732 clean reads, including a total of 4.65 Gb of nucleotides. These reads were assembled into 47 677 contigs (average 443 bp) and then clustered into 30 516 unigenes (average 756 bp). Based on BLAST hits with known proteins in different databases, 20 921 unigenes were annotated with a cut-off E-value of 10(-5). The transcriptome sequences were further annotated using the Clusters of Orthologous Groups, Gene Orthology and the Kyoto Encyclopedia of Genes and Genomes databases. Functional genes involved in spermatogenesis were analysed, including cell cycle proteins, metalloproteins, actin, and ubiquitin and antihyperthermia proteins. Several testis-specific genes were also identified. The transcripts database will help us to understand the molecular mechanisms underlying spermatogenesis in B. dorsalis. Furthermore, 2913 simple sequence repeats and 151 431 single nucleotide polymorphisms were identified, which will be useful for investigating the genetic diversity of B. dorsalis in the future.