RNA-Seq reveals a xenobiotic stress response in the soybean aphid, Aphis glycines, when fed aphid-resistant soybean

Partial funding for Open Access provided by The Ohio State University Open Access Fund.Background: While much recent research has expanded our understanding of the molecular interactions between aphids and their host plants, it is lacking for the soybean aphid, Aphis glycines. Since its North American invasion, A. glycines has become one of the most damaging insect pests on this important crop. Five soybean genes for host plant resistance to A. glycines have been identified, but populations of A. glycines have already adapted to overcome these resistance genes. Understanding the molecular interactions between resistant soybean and A. glycines can provide clues to its adaptation mechanisms. Here, we used RNA-Sequencing to compare and contrast A. glycines gene expression when fed resistant (Rag1) and susceptible soybean. Results: Combining results from a previous A. glycines transcriptome, we generated 64,860 high quality transcripts, totaling 41,151,086 bases. Statistical analysis revealed 914 genes with significant differential expression. Most genes with higher expression in A. glycines on resistant plants (N = 352) were related to stress and detoxification such as cytochrome P450s, glutathione-S-transferases, carboxyesterases, and ABC transporters. A total of 562 genes showed lower transcript abundance in A. glycines on resistant plants. From our extensive transcriptome data, we also identified genes encoding for putative salivary effector proteins (N = 73). Among these, 6 effector genes have lower transcript abundance in A. glycines feeding on resistant soybean. Conclusions: Overall, A. glycines exhibited a pattern typical of xenobiotic challenge, thereby validating antibiosis in Rag1, presumably mediated through toxic secondary metabolites. Additionally, this study identified many A. glycines genes and gene families at the forefront of its molecular interaction with soybean. Further investigation of these genes in other biotypes may reveal adaptation mechanisms to resistant plants

Tissue and Life Stage Specificity of Glutathione S-Transferase Expression in the Hessian Fly, Mayetiola destructor: Implications for Resistance to Host Allelochemicals

Two new Delta and Sigma glutathione S-transferases (GSTs) in the Hessian fly, Mayetiola destructor (Diptera: Cecidomyiidae), were characterized and transcription profiles described. The deduced amino acid sequences for the two M. destructor Delta GSTs (MdesGST-1 and MdesGST-3) showed high similarity with other insect Delta GSTs including the conserved catalytic serine residue. The deduced amino acid sequence for the M. destructor Sigma GST (MdesGST-2) showed high similarity with other insect Sigma GSTs including the conserved glutathione and substrate binding sites. Quantitative tissue expression analysis showed that mRNA levels for MdesGST-1 were predominant in fat body, whereas for MdesGST-2 and MdesGST-3 expression was predominant in the midgut. Temporal expression during development showed peak mRNA levels for MdesGST-1 during larval development, but in the pupal stage for MdesGST-2. MdesGST-3 showed a constitutive expression pattern throughout development. M. destructor feeds on wheat, and expression analysis after feeding indicated that mRNA levels for MdesGST-1 were significantly higher in incompatible interactions in which larvae fed on resistant wheat, while MdesGST-3 was significantly higher in compatible interactions when larvae fed on susceptible wheat. MdesGST-2 showed an equivalent expression pattern during both interactions. These results suggest that the M. destructor Delta GSTs are important in detoxifying wheat allelochemicals during feeding, while Sigma GST participates in metabolism of endogenous substrates

Unbiased Transcriptional Comparisons of Generalist and Specialist Herbivores Feeding on Progressively Defenseless Nicotiana attenuata Plants

Background Herbivore feeding elicits dramatic increases in defenses, most of which require jasmonate (JA) signaling, and against which specialist herbivores are thought to be better adapted than generalist herbivores. Unbiased transcriptional analyses of how neonate larvae cope with these induced plant defenses are lacking. Methodology/Principal Findings We created cDNA microarrays for Manduca sexta and Heliothis virescens separately, by spotting normalized midgut-specific cDNA libraries created from larvae that fed for 24 hours on MeJA-elicited wild-type (WT) Nicotiana attenuata plants. These microarrays were hybridized with labeled probes from neonates that fed for 24 hours on WT and isogenic plants progressively silenced in JA-mediated defenses (N: nicotine; N/PI: N and trypsin protease inhibitors; JA: all JA-mediated defenses). H. virescens neonates regulated 16 times more genes than did M. sexta neonates when they fed on plants silenced in JA-mediated defenses, and for both species, the greater the number of defenses silenced in the host plant (JA > N/PI > N), the greater were the number of transcripts regulated in the larvae. M. sexta larvae tended to down-regulate while H. virescens larvae up- and down-regulated transcripts from the same functional categories of genes. M. sexta larvae regulated transcripts in a diet-specific manner, while H. virescens larvae regulated a similar suite of transcripts across all diet types. Conclusions/Significance The observations are consistent with the expectation that specialists are better adapted than generalist herbivores to the defense responses elicited in their host plants by their feeding. While M. sexta larvae appear to be better adapted to N. attenuata's defenses, some of the elicited responses remain effective defenses against both herbivore species. The regulated genes provide novel insights into larval adaptations to N. attenuata's induced defenses, and represent potential targets for plant-mediated RNAi to falsify hypotheses about the process of adaptation

The insulin signaling pathway in \u3cem\u3eDrosophila melanogaster\u3c/em\u3e: A nexus revealing an Achilles\u27 heel in DDT resistance

Insecticide resistance is an ongoing challenge in agriculture and disease vector control. Here, we demonstrate a novel strategy to attenuate resistance. We used genomics tools to target fundamental energy-associated pathways and identified a potential “Achilles\u27 heel” for resistance, a resistance-associated protein that, upon inhibition, results in a substantial loss in the resistance phenotype. Specifically, we compared the gene expression profiles and structural variations of the insulin/insulin-like growth factor signaling (IIS) pathway genes in DDT-susceptible (91-C) and -resistant (91-R) Drosophila melanogaster (Drosophila) strains. A total of eight and seven IIS transcripts were up- and down-regulated, respectively, in 91-R compared to 91-C. A total of 114 nonsynonymous mutations were observed between 91-C and 91-R, of which 51.8% were fixed. Among the differentially expressed transcripts, phosphoenolpyruvate carboxykinase (PEPCK), down-regulated in 91-R, encoded the greatest number of amino acid changes, prompting us to perform PEPCK inhibitor–pesticide exposure bioassays. The inhibitor of PEPCK, hydrazine sulfate, resulted in a 161- to 218-fold decrease in the DDT resistance phenotype (91-R) and more than a 4- to 5-fold increase in susceptibility in 91-C. A second target protein, Glycogen synthase kinase 3β (GSK3β-PO), had one amino acid difference between 91-C and 91-R, and the corresponding transcript was also down-regulated in 91-R. A GSK3β-PO inhibitor, lithium chloride, likewise reduced the resistance but to a lesser extent than did hydrazine sulfate for PEPCK. We demonstrate the potential role of IIS genes in DDT resistance and the potential discovery of an “Achilles\u27 heel” against pesticide resistance in this pathway

RNA-Seq and molecular docking reveal multi-level pesticide resistance in the bed bug

<p>Abstract</p> <p>Background</p> <p>Bed bugs (<it>Cimex lectularius</it>) are hematophagous nocturnal parasites of humans that have attained high impact status due to their worldwide resurgence. The sudden and rampant resurgence of <it>C. lectularius </it>has been attributed to numerous factors including frequent international travel, narrower pest management practices, and insecticide resistance.</p> <p>Results</p> <p>We performed a next-generation RNA sequencing (RNA-<it>Seq</it>) experiment to find differentially expressed genes between pesticide-resistant (PR) and pesticide-susceptible (PS) strains of <it>C. lectularius</it>. A reference transcriptome database of 51,492 expressed sequence tags (ESTs) was created by combining the databases derived from <it>de novo </it>assembled mRNA-<it>Seq </it>tags (30,404 ESTs) and our previous 454 pyrosequenced database (21,088 ESTs). The two-way GLMseq analysis revealed ~15,000 highly significant differentially expressed ESTs between the PR and PS strains. Among the top 5,000 differentially expressed ESTs, 109 putative defense genes (cuticular proteins, cytochrome P450s, antioxidant genes, ABC transporters, glutathione <it>S</it>-transferases, carboxylesterases and acetyl cholinesterase) involved in penetration resistance and metabolic resistance were identified. Tissue and development-specific expression of P450 CYP3 clan members showed high mRNA levels in the cuticle, Malpighian tubules, and midgut; and in early instar nymphs, respectively. Lastly, molecular modeling and docking of a candidate cytochrome P450 (CYP397A1V2) revealed the flexibility of the deduced protein to metabolize a broad range of insecticide substrates including DDT, deltamethrin, permethrin, and imidacloprid.</p> <p>Conclusions</p> <p>We developed significant molecular resources for <it>C. lectularius </it>putatively involved in metabolic resistance as well as those participating in other modes of insecticide resistance. RNA-<it>Seq </it>profiles of PR strains combined with tissue-specific profiles and molecular docking revealed multi-level insecticide resistance in <it>C. lectularius</it>. Future research that is targeted towards RNA interference (RNAi) on the identified metabolic targets such as cytochrome P450s and cuticular proteins could lay the foundation for a better understanding of the genetic basis of insecticide resistance in <it>C. lectularius</it>.</p

Dietary Antioxidants Impact DDT Resistance in \u3cem\u3eDrosophila melanogaster\u3c/em\u3e

Insects experience a diversity of subtoxic and/or toxic xenobiotics through exposure to pesticides and, in the case of herbivorous insects, through plant defensive compounds in their diets. Many insects are also concurrently exposed to antioxidants in their diets. The impact of dietary antioxidants on the toxicity of xenobiotics in insects is not well understood, in part due to the challenge of developing appropriate systems in which doses and exposure times (of both the antioxidants and the xenobiotics) can be controlled and outcomes can be easily measured. However, in Drosophila melanogaster, a well-established insect model system, both dietary factors and pesticide exposure can be easily controlled. Additionally, the mode of action and xenobiotic metabolism of dichlorodiphenyltrichloroethane (DDT), a highly persistent neurotoxic organochlorine insecticide that is detected widely in the environment, have been well studied in DDT-susceptible and -resistant strains. Using a glass-vial bioassay system with blue diet as the food source, seven compounds with known antioxidant effects (ascorbic acid, β-carotene, glutathione, α-lipoic acid, melatonin, minocycline, and serotonin) were orally tested for their impact on DDT toxicity across three strains of D. melanogaster: one highly susceptible to DDT (Canton-S), one mildly susceptible (91-C), and one highly resistant (91-R). Three of the antioxidants (serotonin, ascorbic acid, and β-carotene) significantly impacted the toxicity of DDT in one or more strains. Fly strain and gender, antioxidant type, and antioxidant dose all affected the relative toxicity of DDT. Our work demonstrates that dietary antioxidants can potentially alter the toxicity of a xenobiotic in an insect population

Dietary antioxidant vitamin C influences the evolutionary path of insecticide resistance in \u3cem\u3eDrosophila melanogaster\u3c/em\u3e

Herbivorous insects encounter a variety of toxic environmental substances ranging from ingested plant defensive compounds to human-introduced insecticidal agents. Dietary antioxidants are known to reduce the negative physiological impacts of toxins in mammalian systems through amelioration of reactive oxygen-related cellular damage. The analogous impacts to insects caused by multigenerational exposure to pesticides and the effects on adaptive responses within insect populations, however, are currently unknown. To address these research gaps, we used Drosophila as a model system to explore adaptive phenotypic responses to acute dichlorodiphenyltrichloroethane (DDT) exposure in the presence of the dietary antioxidant vitamin C and to examine the structural genomic consequences of this exposure. DDT resistance increased significantly among four replicates exposed to a low concentration of DDT for 10 generations. In contrast, dietary intake of vitamin C significantly reduced DDT resistance after mutigenerational exposure to the same concentration of DDT. As to the genomic consequences, no significant differences were predicted in overall nucleotide substitution rates across the genome between any of the treatments. Despite this, replicates exposed to a low concentration of DDT without vitamin C showed the highest number of synonymous and non-synonymous variants (3196 in total), followed by the DDT plus vitamin C (1174 in total), and vitamin C alone (728 in total) treatments. This study demonstrates the potential role of diet (specifically, antioxidant intake) on adaptive genome responses, and thus on the evolution of pesticide resistance within insect populations

Tissue-Specific Transcriptomics of the Exotic Invasive Insect Pest Emerald Ash Borer (Agrilus planipennis)

BACKGROUND: The insect midgut and fat body represent major tissue interfaces that deal with several important physiological functions including digestion, detoxification and immune response. The emerald ash borer (Agrilus planipennis), is an exotic invasive insect pest that has killed millions of ash trees (Fraxinus spp.) primarily in the Midwestern United States and Ontario, Canada. However, despite its high impact status little knowledge exists for A. planipennis at the molecular level. METHODOLOGY AND PRINCIPAL FINDINGS: Newer-generation Roche-454 pyrosequencing was used to obtain 126,185 reads for the midgut and 240,848 reads for the fat body, which were assembled into 25,173 and 37,661 high quality expressed sequence tags (ESTs) for the midgut and the fat body of A. planipennis larvae, respectively. Among these ESTs, 36% of the midgut and 38% of the fat body sequences showed similarity to proteins in the GenBank nr database. A high number of the midgut sequences contained chitin-binding peritrophin (248)and trypsin (98) domains; while the fat body sequences showed high occurrence of cytochrome P450s (85) and protein kinase (123) domains. Further, the midgut transcriptome of A. planipennis revealed putative microbial transcripts encoding for cell-wall degrading enzymes such as polygalacturonases and endoglucanases. A significant number of SNPs (137 in midgut and 347 in fat body) and microsatellite loci (317 in midgut and 571 in fat body) were predicted in the A. planipennis transcripts. An initial assessment of cytochrome P450s belonging to various CYP clades revealed distinct expression patterns at the tissue level. CONCLUSIONS AND SIGNIFICANCE: To our knowledge this study is one of the first to illuminate tissue-specific gene expression in an invasive insect of high ecological and economic consequence. These findings will lay the foundation for future gene expression and functional studies in A. planipennis

Transcriptomics of the Bed Bug (Cimex lectularius)

BACKGROUND: Bed bugs (Cimex lectularius) are blood-feeding insects poised to become one of the major pests in households throughout the United States. Resistance of C. lectularius to insecticides/pesticides is one factor thought to be involved in its sudden resurgence. Despite its high-impact status, scant knowledge exists at the genomic level for C. lectularius. Hence, we subjected the C. lectularius transcriptome to 454 pyrosequencing in order to identify potential genes involved in pesticide resistance. METHODOLOGY AND PRINCIPAL FINDINGS: Using 454 pyrosequencing, we obtained a total of 216,419 reads with 79,596,412 bp, which were assembled into 35,646 expressed sequence tags (3902 contigs and 31744 singletons). Nearly 85.9% of the C. lectularius sequences showed similarity to insect sequences, but 44.8% of the deduced proteins of C. lectularius did not show similarity with sequences in the GenBank non-redundant database. KEGG analysis revealed putative members of several detoxification pathways involved in pesticide resistance. Lamprin domains, Protein Kinase domains, Protein Tyrosine Kinase domains and cytochrome P450 domains were among the top Pfam domains predicted for the C. lectularius sequences. An initial assessment of putative defense genes, including a cytochrome P450 and a glutathione-S-transferase (GST), revealed high transcript levels for the cytochrome P450 (CYP9) in pesticide-exposed versus pesticide-susceptible C. lectularius populations. A significant number of single nucleotide polymorphisms (296) and microsatellite loci (370) were predicted in the C. lectularius sequences. Furthermore, 59 putative sequences of Wolbachia were retrieved from the database. CONCLUSIONS: To our knowledge this is the first study to elucidate the genetic makeup of C. lectularius. This pyrosequencing effort provides clues to the identification of potential detoxification genes involved in pesticide resistance of C. lectularius and lays the foundation for future functional genomics studies