The insulin signaling pathway in Drosophila melanogaster: A nexus revealing an "Achilles' 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' 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' heel" against pesticide resistance in this pathway.

Dietary antioxidant vitamin C influences the evolutionary path of insecticide resistance in Drosophila melanogaster

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.

Dietary antioxidants impact DDT resistance in Drosophila melanogaster

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.

A review of DDT resistance as it pertains to the 91-C and 91-R strains in Drosophila melanogaster

While insecticide resistance presents a challenge for those intent on controlling insect populations, these challenges have also generated a set of tools that can be used to ask fundamental biological questions about that resistance. Numerous species of insects have evolved resistance to multiple classes of insecticides. Each one of these species and their respective resistant populations represent a potential tool for understanding the molecular basis of the evolution of resistance. However, in-laboratory maintenance of resistant insect populations (and their comparative susceptible populations) suitable for asking the needed set of questions around the molecular consequences of long-term pesticide exposure requires a significant, in places prohibitive, level of resources. Drosophila melanogaster (hereafter referred to as Drosophila) is a model insect system with populations easily selected with pesticides and readily maintainable over decades. Even within Drosophila, however, few populations exist where long-term pesticide selection has occurred along with contrasting non-selected population. As such, the Drosophila 91-C and 91-R populations, which exhibit insecticide resistance to DDT (91-R), compared to a non-selection population (91-C), represent a unique resource for the study of high level DDT resistance. Moreover, with the availability of "omics" technologies over the past several decades, this paired population has emerged as a useful tool for understanding both the molecular basis of pesticide resistance and the molecular consequences of long-term pesticide exposure. In this review, we summarize the studies with these aforementioned populations over the past several decades, addressing what has been learned from these efforts.

Variation in Mitochondria-Derived Transcript Levels Associated With DDT Resistance in the 91-R Strain of Drosophila melanogaster (Diptera: Drosophilidae)

The organochloride insecticide dichlorodiphenyltrichloroethane (DDT) and its metabolites can increase cellular levels of reactive oxygen species (ROS), cause mitochondrial dysfunction, and induce apoptosis. The highly DDT-resistant Drosophila melanogaster Meigen 1830 (Drosophila) strain, 91-R, and its susceptible control, 91-C, were used to investigate functional and structural changes among mitochondrial-derived pathways. Resequencing of mitochondrial genomes (mitogenomes) detected no structural differences between 91-R and 91-C, whereas RNA-seq suggested the differential expression of 221 mitochondrial-associated genes. Reverse transcriptase-quantitative PCR validation of 33 candidates confirmed that transcripts for six genes (Cyp12d1-p, Cyp12a4, cyt-c-d, COX5BL, COX7AL, CG17140) were significantly upregulated and two genes (Dif, Rel) were significantly downregulated in 91-R. Among the upregulated genes, four genes are duplicated within the reference genome (cyt-c-d, CG17140, COX5BL, and COX7AL). The predicted functions of the differentially expressed genes, or known functions of closely related genes, suggest that 91-R utilizes existing ROS regulation pathways of the mitochondria to combat increased ROS levels from exposure to DDT. This study represents, to our knowledge, the initial investigation of mitochondrial genome sequence variants and functional adaptations in responses to intense DDT selection and provides insights into potential adaptations of ROS management associated with DDT selection in Drosophila.

Transcriptomic Analyses of Secreted Proteins From the Salivary Glands of Wheat Midge Larvae

Both the wheat midge (Sitodiplosis mosellana) (Géhin) (Diptera: Cecidomyiidae) and the Hessian fly (Mayetiola destructor) (Say) (Diptera: Cecidomyiidae) belong to a group of insects called gall midges (Diptera: Cecidomyiidae), and both are destructive pests of wheat. From Hessian fly larvae, a large number of genes have been identified to encode secreted salivary gland proteins (SSGPs), which are presumably critical for the insect to feed on and manipulate host plants. For comparison, we conducted an analysis on transcripts encoding SSGPs from the first instar larvae of the wheat midge. In total, 3,500 cDNA clones were sequenced, from which 1,301 high-quality sequences were obtained. Approximately 25% of the cDNAs with high-quality sequences encoded SSGPs. The SSGPs were grouped into 97 groups based on sequence homology. Among the SSGP-encoding transcripts, 206 encoded unique proteins with no sequence similarity to any known protein and 29 encoded proteins similar to known proteins including proteases, serpines, thioesterases, ankyrins, and ferritins. Most (~80%) SSGP-encoding genes appear under strong selection for mutations that generate amino acid changes within the coding region. Identification and characterization of SSGPs in wheat midge larvae provide a foundation for future work to reveal molecular mechanisms behind wheat midge–wheat interactions and the role of these putative effector proteins in insect virulence. Availability of the SSGP transcripts will also facilitate comparative analyses of insect effectors from related species.

The complete mitogenome of the Emerald Ash Borer (EAB), Agrilus planipennis (Insecta: Coleoptera: Buprestidae)

The complete mitogenome of the Emerald Ash Borer (EAB, Agrilus planipennis) was obtained by gleaning mitochondrial sequences from whole-genome Illumina sequencing data. The circular genome has 15,942 base pairs and contains 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), 2 ribosomal RNAs (rRNAs) and an A-T-rich region. All PCGs begin with ATN codons. The nucleotide composition is highly asymmetric (31.65% A, 40.25% T, 17.39% G, 10.71% C), with an overall A-T content of 71.9%. Phylogenetic analysis based on insect mitogenomes indicated that EAB is closely related to other Buprestoidea species, clustering most closely with Chrysochroa fulgidissima.

Transcriptome Analysis of the Emerald Ash Borer (EAB), Agrilus planipennis: De Novo Assembly, Functional Annotation and Comparative Analysis

Background

The Emerald ash borer (EAB), Agrilus planipennis, is an invasive phloem-feeding insect pest of ash trees. Since its initial discovery near the Detroit, US- Windsor, Canada area in 2002, the spread of EAB has had strong negative economic, social and environmental impacts in both countries. Several transcriptomes from specific tissues including midgut, fat body and antenna have recently been generated. However, the relatively low sequence depth, gene coverage and completeness limited the usefulness of these EAB databases.

Methodology and Principal Findings

High-throughput deep RNA-Sequencing (RNA-Seq) was used to obtain 473.9 million pairs of 100 bp length paired-end reads from various life stages and tissues. These reads were assembled into 88,907 contigs using the Trinity strategy and integrated into 38,160 unigenes after redundant sequences were removed. We annotated 11,229 unigenes by searching against the public nr, Swiss-Prot and COG. The EAB transcriptome assembly was compared with 13 other sequenced insect species, resulting in the prediction of 536 unigenes that are Coleoptera-specific. Differential gene expression revealed that 290 unigenes are expressed during larval molting and 3,911 unigenes during metamorphosis from larvae to pupae, respectively (FDR< 0.01 and log2 FC>2). In addition, 1,167 differentially expressed unigenes were identified from larval and adult midguts, 435 unigenes were up-regulated in larval midgut and 732 unigenes were up-regulated in adult midgut. Most of the genes involved in RNA interference (RNAi) pathways were identified, which implies the existence of a system RNAi in EAB.

Conclusions and Significance

This study provides one of the most fundamental and comprehensive transcriptome resources available for EAB to date. Identification of the tissue- stage- or species- specific unigenes will benefit the further study of gene functions during growth and metamorphosis processes in EAB and other pest insects.

Evidence for trade-offs in detoxification and chemosensation gene signatures in Plutella xylostella

BACKGROUND

Detoxification genes have been associated with insecticide adaptation in the diamondback moth, Plutella xylostella. The link between chemosensation genes and adaptation, however, remains unexplored. To gain a better understanding of the involvement of these genes in insecticide adaptation, the authors exposed lines of P. xylostella to either high uniform (HU) or low heterogeneous (LH) concentrations of permethrin, expecting primarily physiological or behavioral selection respectively. Initially, 454 pyrosequencing was applied, followed by an examination of expression profiles of candidate genes that responded to selection [cytochrome P450 (CYP), glutathione S-transferase (GST), carboxylesterase (CarE), chemosensory protein (CSP) and odorant-binding protein (OBP)] by quantitative PCR in the larvae. Toxicity and behavioral assays were also conducted to document the effects of the two forms of exposure.

RESULTS

Pyrosequencing of the P. xylostella transcriptome from adult heads and third instars produced 198,753 reads with 52,752,486 bases. Quantitative PCR revealed overexpression of CYP4M14, CYP305B1 and CSP8 in HU larvae. OBP13, however, was highest in LH. Larvae from LH and HU lines had up to five- and 752-fold resistance levels respectively, which could be due to overexpression of P450s. However, the behavioral responses of all lines to a series of permethrin concentrations did not vary significantly in any of the generations examined, in spite of the observed upregulation of CSP8 and OBP13.

CONCLUSION

Expression patterns from the target genes provide insights into behavioral and physiological responses to permethrin and suggest a new avenue of research on the role of chemosensation genes in insect adaptation to toxins.

Core RNAi machinery and gene knockdown in the emerald ash borer (Agrilus planipennis)

The RNA interference (RNAi) technology has been widely used in insect functional genomics research and provides an alternative approach for insect pest management. To understand whether the emerald ash borer (Agrilus planipennis), an invasive and destructive coleopteran insect pest of ash tree (Fraxinus spp.), possesses a strong RNAi machinery that is capable of degrading target mRNA as a response to exogenous double-stranded RNA (dsRNA) induction, we identified three RNAi pathway core component genes, Dicer-2, Argonaute-2 and R2D2, from the A. planipennis genome sequence. Characterization of these core components revealed that they contain conserved domains essential for the proteins to function in the RNAi pathway. Phylogenetic analyses showed that they are closely related to homologs derived from other coleopteran species. We also delivered the dsRNA fragment of AplaScrB-2, a β-fructofuranosidase-encoding gene horizontally acquired by A. planipennis as we reported previously, into A. planipennis adults through microinjection. Quantitative real-time PCR analysis on the dsRNA-treated beetles demonstrated a significantly decreased gene expression level of AplaScrB-2 appearing on day 2 and lasting until at least day 6. This study is the first record of RNAi applied in A. planipennis.

Characterization of horizontally transferred β-fructofuranosidase (ScrB) genes in Agrilus planipennis

The emerald ash borer (Agrilus planipennis) is an important invasive insect pest of Fraxinus spp. that feeds on host tissues containing high levels of sucrose. However, little is known about how it digests sucrose. Here, using larval midgut transcriptome data and preliminary genome sequence efforts, two β-fructofuranosidase-encoding ScrB genes, AplaScrB-1 and AplaScrB-2, were identified, and proved to reside within the A. planipennis genome. Homology and phylogenetic analysis revealed that they were acquired by A. planipennis via horizontal gene transfer (HGT) from bacteria, possibly an event independent from that reported in bark beetles (eg ScrB genes). Microsynteny between A. planipennis DNA scaffold #2042940, which hosts AplaScrB-1, and a region in the Tribolium castaneum chromosome LG4 suggested that A. planipennis gained this gene after the separation of Buprestidae and Tenebrionidae. Although both of the putative AplaScrB proteins have conserved β-fructofuranosidase motifs, only AplaScrB-2 was predicted to be a secretory protein. Expression of AplaScrB-1 seemed constitutive during development and in all tissues examined, whereas AplaScrB-2 showed a peak expression in adults and in the midgut. We propose that acquisition of these genes by A. planipennis from bacteria is adaptive, and specifically AplaScrB-2 is involved in breaking down dietary sucrose to obtain energy for development.

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

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.

Molecular characterization and expression analysis of soluble trehalase gene in Aphis glycines, a migratory pest of soybean

In insects, the enzyme trehalase plays a crucial role in energy metabolism, chitin synthesis and possibly during plant-insect interactions. We have characterized a soluble trehalase gene (Tre-1) from cDNA of Aphis glycines, a serious migratory pest of soybean. The full-length cDNA of Tre-1 in A. glycines (AyTre-1) was 2550 bp long with an open reading frame of 1770 bp that encoded for a 589 amino acid residues protein. Sequence assessment and phylogenetic analysis of the putative protein suggested that the selected cDNA belongs to soluble trehalase group. Quantitative PCR (qPCR) analysis in different tissues and developmental stages revealed peak mRNA levels of AyTre-1 in the gut (compared with other tissues assayed) and highest expression in the second instar compared with the other developmental stages assayed. Interestingly, a significantly increased expression of AyTre-1 (1.9-fold, P < 0.05) was observed in the alate morphs compared with that in apterate morphs. However, there was no significant difference in AyTre-1 expression in A. glycines-nymphs fed with resistant and susceptible plants. Expression patterns identified in this study provide a platform to investigate the role of AyTre-1 in physiological activities such as flight and feeding in A. glycines. The characterization of soluble trehalase gene may help to develop novel strategies to manage A. glycines using trehalase inhibitors and using RNA interference for knock-down of AyTre-1 expression.

Molecular characterization of genes encoding inward rectifier potassium (Kir) channels in the bed bug (Cimex lectularius)

The molecular genetics of inward-rectifier potassium (Kir) channels in insects is poorly understood. To date, Kir channel genes have been characterized only from a few representative dipterans (i.e., fruit flies and mosquitoes). The goal of the present study was to characterize Kir channel cDNAs in a hemipteran, the bed bug (Cimex lectularius). Using our previously reported bed bug transcriptome (RNA-seq), we identified two cDNAs that encode putative Kir channels. One was a full-length cDNA that encodes a protein belonging to the insect 'Kir3' clade, which we designate as 'ClKir3'. The other was a partial cDNA that encodes a protein with similarity to both the insect 'Kir1' and 'Kir2' clades, which we designate as 'ClKir1/2'. Quantitative real-time PCR analysis revealed that ClKir1/2 and ClKir3 exhibited peak expression levels in late-instar nymphs and early-instar nymphs, respectively. Furthermore, ClKir3, but not ClKir1/2, showed tissue-specific expression in Malpighian tubules of adult bed bugs. Lastly, using an improved procedure for delivering double-stranded RNA (dsRNA) to male and female bed bugs (via the cervical membrane) we demonstrate rapid and systemic knockdown of ClKir3 transcripts. In conclusion, we demonstrate that the bed bug possesses at least two genes encoding Kir channels, and that RNAi is possible for at least Kir3, thereby offering a potential approach for elucidating the roles of Kir channel genes in bed bug physiology.

Identification of odor-processing genes in the emerald ash borer, Agrilus planipennis

Background

Insects rely on olfaction to locate food, mates, and suitable oviposition sites for successful completion of their life cycle. Agrilus planipennis Fairmaire (emerald ash borer) is a serious invasive insect pest that has killed tens of millions of North American ash (Fraxinus spp) trees and threatens the very existence of the genus Fraxinus. Adult A. planipennis are attracted to host volatiles and conspecifics; however, to date no molecular knowledge exists on olfaction in A. planipennis. Hence, we undertook an antennae-specific transcriptomic study to identify the repertoire of odor processing genes involved in A. planipennis olfaction.

Methodology and Principal Findings

We acquired 139,085 Roche/454 GS FLX transcriptomic reads that were assembled into 30,615 high quality expressed sequence tags (ESTs), including 3,249 isotigs and 27,366 non-isotigs (contigs and singletons). Intriguingly, the majority of the A. planipennis antennal transcripts (59.72%) did not show similarity with sequences deposited in the non-redundant database of GenBank, potentially representing novel genes. Functional annotation and KEGG analysis revealed pathways associated with signaling and detoxification. Several odor processing genes (9 odorant binding proteins, 2 odorant receptors, 1 sensory neuron membrane protein and 134 odorant/xenobiotic degradation enzymes, including cytochrome P450s, glutathione-S-transferases; esterases, etc.) putatively involved in olfaction processes were identified. Quantitative PCR of candidate genes in male and female A. planipennis in different developmental stages revealed developmental- and sex-biased expression patterns.

Conclusions and Significance

The antennal ESTs derived from A. planipennis constitute a rich molecular resource for the identification of genes potentially involved in the olfaction process of A. planipennis. These findings should help in understanding the processing of antennally-active compounds (e.g. 7-epi-sesquithujene) previously identified in this serious invasive pest.

Transcriptome analysis of the salivary glands of potato leafhopper, Empoasca fabae

The potato leafhopper, Empoasca fabae, is a pest of economic crops in the United States and Canada, where it causes damage known as hopperburn. Saliva, along with mechanical injury, leads to decreases in gas exchange rates, stunting and chlorosis. Although E. fabae saliva is known to induce plant responses, little knowledge exists of saliva composition at the molecular level. We subjected the salivary glands of E. fabae to Roche 454-pyrosequencing which resulted significant number (30,893) of expressed sequence tags including 2805 contigs and 28,088 singletons. A high number of sequences (78%) showed similarity to other insect species in GenBank, including Triboliumcastaneum, Drosophilamelanogaster and Acrythosiphonpisum. KEGG analysis predicted the presence of pathways for purine and thiamine metabolic, biosynthesis of secondary metabolites, drug metabolism, and lysine degradation. Pfam analysis showed a high number of cellulase and carboxylesterase protein domains. Expression analysis of candidate genes (alpha amylase, lipase, pectin lyase, etc.) among different tissues revealed tissue-specific expression of digestive enzymes in E. fabae. This is the first study to characterize the sialotranscriptome of E. fabae and the first for any species in the family of Cicadellidae. Due to the status of these insects as economic pests, knowledge of which genes are active in the salivary glands is important for understanding their impact on host plants.

Permanent genetic resources added to Molecular Ecology Resources Database 1 August 2012 - 30 September 2012

This article documents the addition of 83 microsatellite marker loci and 96 pairs of single-nucleotide polymorphism (SNP) sequencing primers to the Molecular Ecology Resources Database. Loci were developed for the following species: Bembidion lampros, Inimicus japonicus, Lymnaea stagnalis, Panopea abbreviata, Pentadesma butyracea, Sycoscapter hirticola and Thanatephorus cucumeris (anamorph: Rhizoctonia solani). These loci were cross-tested on the following species: Pentadesma grandifolia and Pentadesma reyndersii. This article also documents the addition of 96 sequencing primer pairs and 88 allele-specific primers or probes for Plutella xylostella.

Characterization of a chitin synthase encoding gene and effect of diflubenzuron in soybean aphid, Aphis glycines

Chitin synthases are critical enzymes for synthesis of chitin and thus for subsequent growth and development in insects. We identified the cDNA of chitin synthase gene (CHS) in Aphis glycines, the soybean aphid, which is a serious pest of soybean. The full-length cDNA of CHS in A. glycines (AyCHS) was 5802 bp long with an open reading frame of 4704 bp that encoded for a 1567 amino acid residues protein. The predicted AyCHS protein had a molecular mass of 180.05 kDa and its amino acid sequence contained all the signature motifs (EDR, QRRRW and TWGTR) of chitin synthases. The quantitative real-time PCR (qPCR) analysis revealed that AyCHS was expressed in all major tissues (gut, fat body and integument); however, it had the highest expression in integument (~3.5 fold compared to gut). Interestingly, the expression of AyCHS in developing embryos was nearly 7 fold higher compared to adult integument, which probably is a reflection of embryonic molts in hemimetabolus insects. Expression analysis in different developmental stages of A. glycines revealed a consistent AyCHS expression in all stages. Further, through leaf dip bioassay, we tested the effect of diflubenzuron (DFB, Dimilin ®), a chitin-synthesis inhibitor, on A. glycines' survival, fecundity and body weight. When fed with soybean leaves previously dipped in 50 ppm DFB solution, A. glycines nymphs suffered significantly higher mortality compared to control. A. glycines nymphs feeding on diflubenzuron treated leaves showed a slightly enhanced expression (1.67 fold) of AyCHS compared to nymphs on untreated leaves. We discussed the potential applications of the current study to develop novel management strategies using chitin-synthesis inhibitors and using RNAi by knocking down AyCHS expression.

Validation of reference genes for gene expression studies in Aphis glycines (Hemiptera: Aphididae)

Quantitative real-time polymerase chain reaction (qRT-PCR) is a common and robust tool for accurate quantification of mRNA transcripts. To normalize results, a housekeeping gene ([HKG], reference gene or endogenous control gene) is mandatory. Soybean aphid, Aphis glycines Matsumura (Hemiptera: Aphididae), is a significant soybean, Glycine max (L.) Merr., pest, yet gene expression and functional genomics studies are hindered by a lack of stable HKGs. We evaluated seven potential HKGs (SDFS, succinate dehydrogenase flavoprotein subunit; EF1a, elongation factor-la; HEL, helicase; GAPDH, glyceraldehyde-3 phosphate dehydrogenase; RPS9, ribosomal protein S9; TBP, TATA-box binding protein; and UBQ, ubiquitin-conjugating protein) to determine the most efficient HKGs that have stable expression among tissues, developmental stages, and aphids fed on susceptible and host plant-resistant soybean. HKG stability was determined using GeNorm and NormFinder. Results from three different experimental conditions revealed high stability of TBP compared with the other HKGs profiled across the samples assayed. RPS9 showed stable expression among aphids on susceptible and resistant plants, whereas EF1a showed stable expression in tissues and developmental stages. Therefore, we recommend the TBP as a suitable HKG for efficient normalization among treatments, tissues, and developmental stages of A. glycines. In addition, RPS9 may be used for host-plant resistance experiments and EF1a could be considered for testing differential expression across tissues or developmental stages. These results will enable a more accurate and reliable normalization of qRT-PCR data in A. glycines.

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

BACKGROUND

Bed bugs (Cimex lectularius) are hematophagous nocturnal parasites of humans that have attained high impact status due to their worldwide resurgence. The sudden and rampant resurgence of C. lectularius has been attributed to numerous factors including frequent international travel, narrower pest management practices, and insecticide resistance.

RESULTS

We performed a next-generation RNA sequencing (RNA-Seq) experiment to find differentially expressed genes between pesticide-resistant (PR) and pesticide-susceptible (PS) strains of C. lectularius. A reference transcriptome database of 51,492 expressed sequence tags (ESTs) was created by combining the databases derived from de novo assembled mRNA-Seq 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 S-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.

CONCLUSIONS

We developed significant molecular resources for C. lectularius putatively involved in metabolic resistance as well as those participating in other modes of insecticide resistance. RNA-Seq profiles of PR strains combined with tissue-specific profiles and molecular docking revealed multi-level insecticide resistance in C. lectularius. 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 C. lectularius.

Cloning and characterization of mariner-like elements in the soybean aphid, Aphis glycines Matsumura

Soybean aphid, Aphis glycines Matsumura (Hemiptera: Aphididae), is currently the most important insect pest of soybean (Glycine max (L.) Merr.) in the United States and causes significant economic damage worldwide, but little is known about the aphid at the molecular level. Mariner-like transposable elements (MLEs) are ubiquitous within the genomes of arthropods and various other invertebrates. In this study, we report the cloning of MLEs from the soybean aphid genome using degenerate PCR primers designed to amplify conserved regions of mariner transposases. Two of the ten sequenced clones (designated as Agmar1 and Agmar2) contained partial but continuous open reading frames, which shared high levels of homology at the protein level with other mariner transposases from insects and other taxa. Phylogenetic analysis revealed Agmar1 to group within the irritans subfamily of MLEs and Agmar2 within the mellifera subfamily. Southern blot analysis and quantitative PCR analysis indicated a low copy number for Agmar1-like elements within the soybean aphid genome. These results suggest the presence of at least two different putative mariner-like transposases encoded by the soybean aphid genome. Both Agmar1 and Agmar2 could play influential roles in the architecture of the soybean aphid genome. Transposable elements are also thought to potentially mediate resistance in insects through changes in gene amplification and mutations in coding sequences. Finally, Agmar1 and Agmar2 may represent useful genetic tools and provide insights on A. glycines adaptation.

Identification and validation of reference genes for quantitative real-time polymerase chain reaction in Cimex lectularius

Quantitative real-time polymerase chain reaction (qRT-PCR) has emerged as robust methodology for gene expression studies, but reference genes are crucial for accurate normalization. Commonly used reference genes are housekeeping genes that are thought to be nonregulated; however, their expression can be unstable across different experimental conditions. We report the identification and validation of suitable reference genes in the bed bug, Cimex lectularius, by using qRT-PCR. The expression stability of eight reference genes in different tissues (abdominal cuticle, midgut, Malpighian tubules, and ovary) and developmental stages (early instar nymphs, late instar nymphs, and adults) of pesticide-susceptible and pesticide-exposed C. lectularius were analyzed using geNorm, NormFinder, and BestKeeper. Overall expression analysis of the eight reference genes revealed significant variation among samples, indicating the necessity of validating suitable reference genes for accurate quantification of mRNA transcripts. Ribosomal protein (RPL18) exhibited the most stable gene expression across all the tissue and developmental-stage samples; a-tubulin revealed the least stability across all of the samples examined. Thus, we recommend RPL18 as a suitable reference gene for normalization in gene expression studies of C. lectularius.

Antioxidant genes of the emerald ash borer (Agrilus planipennis): gene characterization and expression profiles

Phytophagous insects frequently encounter reactive oxygen species (ROS) from exogenous and endogenous sources. To overcome the effect of ROS, insects have evolved a suite of antioxidant defense genes such as superoxide dismutase (SOD), catalase (CAT), ascorbate peroxidase (APX) and glutathione peroxidase (GPX). The emerald ash borer (Agrilus planipennis Fairmaire), an exotic invasive insect pest from Asia has killed millions of ash trees and continues to invade North America at a rapid pace. From an on-going expressed sequence tag (EST) project of A. planipennis larval tissues, we identified ESTs coding for a Cu-Zn SOD (ApSOD1), a CAT (ApCAT1) and a GPX (ApGPX1). A multiple sequence alignment of the derived A. planipennis sequences revealed high homology with other insect sequences at the amino acid level. Phylogenetic analysis of ApSOD1 grouped it with Cu-Zn SODs of other insect taxa. Quantitative real time PCR (qRT-PCR) analysis in different larval tissues (midgut, fat body, Malpighian tubule and cuticle) revealed high mRNA levels of ApCAT1 in the midgut. Interestingly, high mRNA levels for both ApSOD1 and ApGPX1 were observed in the Malpighian tubules. Assay of mRNA levels in developmental stages (larva, prepupa and adults) by qRT-PCR indicated high transcript levels of ApCAT1 and ApGPX1 in larval and prepupal stages with a decline in adults. On the other hand, the transcript levels of ApSOD1 were observed to be constitutive in all the developmental stages assayed. Results obtained reflect a plausible role of these A. planipennis antioxidant genes in quenching ROS from both diet (ash allelochemicals) as well as endogenous sources. These studies further help in understanding the adaptation/invasiveness of A. planipennis.

Metabolic Resistance in Bed Bugs

Blood-feeding insects have evolved resistance to various insecticides (organochlorines, pyrethroids, carbamates, etc.) through gene mutations and increased metabolism. Bed bugs (Cimex lectularius) are hematophagous ectoparasites that are poised to become one of the major pests in households throughout the United States. Currently, C. lectularius has attained a high global impact status due to its sudden and rampant resurgence. Resistance to pesticides is one factor implicated in this phenomenon. Although much emphasis has been placed on target sensitivity, little to no knowledge is available on the role of key metabolic players (e.g., cytochrome P450s and glutathione S-transferases) towards pesticide resistance in C. lectularius. In this review, we discuss different modes of resistance (target sensitivity, penetration resistance, behavioral resistance, and metabolic resistance) with more emphasis on metabolic resistance.

Transcriptomic signatures of ash (Fraxinus spp.) phloem

BACKGROUND

Ash (Fraxinus spp.) is a dominant tree species throughout urban and forested landscapes of North America (NA). The rapid invasion of NA by emerald ash borer (Agrilus planipennis), a wood-boring beetle endemic to Eastern Asia, has resulted in the death of millions of ash trees and threatens billions more. Larvae feed primarily on phloem tissue, which girdles and kills the tree. While NA ash species including black (F. nigra), green (F. pennsylvannica) and white (F. americana) are highly susceptible, the Asian species Manchurian ash (F. mandshurica) is resistant to A. planipennis perhaps due to their co-evolutionary history. Little is known about the molecular genetics of ash. Hence, we undertook a functional genomics approach to identify the repertoire of genes expressed in ash phloem.

METHODOLOGY AND PRINCIPAL FINDINGS

Using 454 pyrosequencing we obtained 58,673 high quality ash sequences from pooled phloem samples of green, white, black, blue and Manchurian ash. Intriguingly, 45% of the deduced proteins were not significantly similar to any sequences in the GenBank non-redundant database. KEGG analysis of the ash sequences revealed a high occurrence of defense related genes. Expression analysis of early regulators potentially involved in plant defense (i.e. transcription factors, calcium dependent protein kinases and a lipoxygenase 3) revealed higher mRNA levels in resistant ash compared to susceptible ash species. Lastly, we predicted a total of 1,272 single nucleotide polymorphisms and 980 microsatellite loci, among which seven microsatellite loci showed polymorphism between different ash species.

CONCLUSIONS AND SIGNIFICANCE

The current transcriptomic data provide an invaluable resource for understanding the genetic make-up of ash phloem, the target tissue of A. planipennis. These data along with future functional studies could lead to the identification/characterization of defense genes involved in resistance of ash to A. planipennis, and in future ash breeding programs for marker development.

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.

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.

Decreased detoxification genes and genome size make the human body louse an efficient model to study xenobiotic metabolism

The human body louse, Pediculus humanus humanus, has one of the smallest insect genomes, containing ∼10 775 annotated genes. Annotation of detoxification [cytochrome P450 monooxygenase (P450), glutathione-S-transferase (GST), esterase (Est) and ATP-binding cassette transporter (ABC transporter)] genes revealed that they are dramatically reduced in P. h. humanus compared to other insects except for Apis mellifera. There are 37 P450, 13 GST and 17 Est genes present in P. h. humanus, approximately half the number found in Drosophila melanogaster and Anopheles gambiae. The number of putatively functional ABC transporter genes in P. h. humanus and Ap. mellifera are the same (36) but both have fewer than An. gambiae (44) or Dr. melanogaster (65). The reduction of detoxification genes in P. h. humanus may be a result of this louse's simple life history, in which it does not encounter a wide variety of xenobiotics. Neuronal component genes are highly conserved across different insect species as expected because of their critical function. Although reduced in number, P. h. humanus still retains at least a minimum repertoire of genes known to confer metabolic or toxicokinetic resistance to xenobiotics (eg Cyp3 clade P450s, Delta GSTs, B clade Ests and B/C subfamily ABC transporters), suggestive of its high potential for resistance development.

The gut transcriptome of a gall midge, Mayetiola destructor

The Hessian fly, Mayetiola destructor, is a serious pest of wheat and an experimental organism for the study of gall midge-plant interactions. In addition to food digestion and detoxification, the gut of Hessian fly larvae is also an important interface for insect-host interactions. Analysis of the genes expressed in the Hessian fly larval gut will enhance our understanding of the overall gut physiology and may also lead to the identification of critical molecules for Hessian fly-host plant interactions. Over 10,000 Expressed Sequence Tags (ESTs) were generated and assembled into 2007 clusters. The most striking feature of the Hessian fly larval transcriptome is the existence of a large number of transcripts coding for so-called small secretory proteins (SSP) with amino acids less than 250. Eleven of the 30 largest clusters were SSP transcripts with the largest cluster containing 11.3% of total ESTs. Transcripts coding for diverse digestive enzymes and detoxification proteins were also identified. Putative digestive enzymes included trypsins, chymotrypsins, cysteine proteases, aspartic protease, endo-oligopeptidase, aminopeptidases, carboxypeptidases, and alpha-amylases. Putative detoxification proteins included cytochrome P450s, glutathione S-transferases, peroxidases, ferritins, a catalase, peroxiredoxins, and others. This study represents the first global analysis of gut transcripts from a gall midge. The identification of a large number of transcripts coding for SSPs, digestive enzymes, detoxification proteins in the Hessian fly larval gut provides a foundation for future studies on the functions of these genes.

Molecular characterization of mariner-like elements in emerald ash borer, Agrilus planipennis (Coleoptera, Polyphaga)

Emerald ash borer (EAB, Agrilus planipennis), an exotic invasive pest, has killed millions of ash trees (Fraxinus spp.) in North America and continues to threaten the very survival of the entire Fraxinus genus. Despite its high-impact status, to date very little knowledge exists for this devastating insect pest at the molecular level. Mariner-like elements (MLEs) are transposable elements, which are ubiquitous in occurrence in insects and other invertebrates. Because of their low specificity and broad host range, they can be used for epitope-tagging, gene mapping, and in vitro mutagenesis. The majority of the known MLEs are inactive due to in-frame shifts and stop codons within the open reading frame (ORF). We report on the cloning and characterization of two MLEs in A. planipennis genome (Apmar1 and Apmar2). Southern analysis indicated a very high copy number for Apmar1 and a moderate copy number for Apmar2. Phylogenetic analysis revealed that both elements belong to the irritans subfamily. Based on the high copy number for Apmar1, the full-length sequence was obtained using degenerate primers designed to the inverted terminal repeat (ITR) sequences of irritans MLEs. The recovered nucleotide sequence for Apmar1 consisted of 1,292 bases with perfect ITRs, and an ORF of 1,050 bases encoding a putative transposase of 349 amino acids. The deduced amino acid sequence of Apmar1 contained the conserved regions of mariner transposases including WVPHEL and YSPDLAP, and the D,D(34)D motif. Both Apmar1 and Apmar2 could represent useful genetic tools and provide insights on EAB adaptation.

Combining next-generation sequencing strategies for rapid molecular resource development from an invasive aphid species, Aphis glycines

Background

Aphids are one of the most important insect taxa in terms of ecology, evolutionary biology, genetics and genomics, and interactions with endosymbionts. Additionally, many aphids are serious pest species of agricultural and horticultural plants. Recent genetic and genomic research has expanded molecular resources for many aphid species, including the whole genome sequencing of the pea aphid, Acrythosiphon pisum. However, the invasive soybean aphid, Aphis glycines, lacks in any significant molecular resources.

Methodology/Principal Findings

Two next-generation sequencing technologies (Roche-454 and Illumina GA-II) were used in a combined approach to develop both transcriptomic and genomic resources, including expressed genes and molecular markers. Over 278 million bp were sequenced among the two methods, resulting in 19,293 transcripts and 56,688 genomic sequences. From this data set, 635 SNPs and 1,382 microsatellite markers were identified. For each sequencing method, different soybean aphid biotypes were used which revealed potential biotype specific markers. In addition, we uncovered 39,822 bp of sequence that were related to the obligatory endosymbiont, Buchnera aphidicola, as well as sequences that suggest the presence of Hamiltonella defensa, a facultative endosymbiont.

Conclusions and Significance

Molecular resources for an invasive, non-model aphid species were generated. Additionally, the power of next-generation sequencing to uncover endosymbionts was demonstrated. The resources presented here will complement ongoing molecular studies within the Aphididae, including the pea aphid whole genome, lead to better understanding of aphid adaptation and evolution, and help provide novel targets for soybean aphid control.

Bowman-Birk inhibitor affects pathways associated with energy metabolism in Drosophila melanogaster

Bowman-Birk inhibitor (BBI) is toxic when fed to certain insects, including the fruit fly, Drosophila melanogaster. Dietary BBI has been demonstrated to slow growth and increase insect mortality by inhibiting the digestive enzymes trypsin and chymotrypsin, resulting in a reduced supply of amino acids. In mammals, BBI influences cellular energy metabolism. Therefore, we tested the hypothesis that dietary BBI affects energy-associated pathways in the D. melanogaster midgut. Through microarray and metabolomic analyses, we show that dietary BBI affects energy utilization pathways in the midgut cells of D. melanogaster. In addition, ultrastructure studies indicate that microvilli are significantly shortened in BBI-fed larvae. These data provide further insights into the complex cellular response of insects to dietary protease inhibitors.

Analysis of gene expression in emerald ash borer (Agrilus planipennis) using quantitative real time-PCR

Emerald ash borer (EAB, Agrilus planipennis) is an exotic invasive pest, which has killed millions of ash trees (Fraxinus spp) in North America. EAB continues to spread rapidly and attacks ash trees of different ages, from saplings to mature trees. However, to date very little or no molecular knowledge exists for EAB. We are interested in deciphering the molecular-based physiological processes at the tissue level that aid EAB in successful colonization of ash trees. In this report we show the effective use of quantitative real-time PCR (qRT-PCR) to ascertain mRNA levels in different larval tissues (including midgut, fat bodies and cuticle) and different developmental stages (including 1(st)-, 2(nd)-, 3(rd)-, 4(th)-instars, prepupae and adults) of EAB. As an example, a peritrophin gene (herein named, AP-PERI1) is exemplified as the gene of interest and a ribosomal protein (AP-RP1) as the internal control. Peritrophins are important components of the peritrophic membrane/matrix (PM), which is the lining of the insect gut. The PM has diverse functions including digestion and mechanical protection to the midgut epithelium.

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.

Transcriptional signatures in response to wheat germ agglutinin and starvation in Drosophila melanogaster larval midgut

One function of plant lectins such as wheat germ agglutinin is to serve as defences against herbivorous insects. The midgut is one critical site affected by dietary lectins. We observed marked cellular, structural and gene expression changes in the midguts of Drosophila melanogaster third instar larvae that were fed wheat germ agglutinin. Some of these changes were similar to those observed in the midguts of starved D. melanogaster. Dietary wheat germ agglutinin caused shortening, branching, swelling, distortion and in some cases disintegration of the midgut microvilli. Starvation was accompanied primarily by shortening of the microvilli. Microarray analyses revealed that dietary wheat germ agglutinin evoked differential expression of 61 transcripts; seven of these were also differentially expressed in starved D. melanogaster. The differentially transcribed gene clusters in wheat germ agglutinin-fed larvae were associated with (1) cytoskeleton organization; (2) digestive enzymes; (3) detoxification reactions; and (4) energy metabolism. Four possible transcription factor binding motifs were associated with the differentially expressed genes. One of these exhibited substantial similarity to MyoD, a transcription factor binding motif associated with cellular structures in mammals. These results are consistent with the hypothesis that wheat germ agglutinin caused a starvation-like effect and structural changes of midgut cells of D. melanogaster third-instar larvae.

Characterization and expression analysis of a gene encoding a secreted lipase-like protein expressed in the salivary glands of the larval Hessian fly, Mayetiola destructor (Say)

In a salivary gland transcriptomics study we identified a cDNA with a full-length open reading frame for a gene (MdesL1) encoding a lipase-like protein expressed in the salivary glands of the larval Hessian fly, Mayetiola destructor (Say). Fluorescent in situ hybridization on salivary polytenes positioned MdesL1 on the long arm of Autosome 1. BLASTp and conserved domain searches revealed the deduced amino acid sequence contained a lipase superfamily domain with similarity to lipases and phospholipases from other insects. A secretion signal peptide was identified at the amino terminus of the deduced amino acid sequence. Analysis of the transcript of MdesL1 in larval Hessian fly tissues by quantitative real-time PCR (qPCR) revealed the greatest abundance was in salivary glands. Analysis of transcript levels during development showed the greatest level was detected in feeding 1st-instar and early 2nd-instar larvae. Transcript levels increased dramatically over time in larvae feeding on susceptible wheat but were detected at low levels in larvae feeding on resistant wheat. These data suggest the protein encoded by MdesL1 is likely secreted into host-plant cells during larval feeding and could be involved in extra-oral digestion and changes in host-cell permeability or in generating a second messenger in a host-cell-signaling cascade.

Transcriptomic profiles of Drosophila melanogaster third instar larval midgut and responses to oxidative stress

Oligoarray analysis was used to determine the number and nature of genes expressed in third instar Drosophila melanogaster larval midguts. The majority of transcripts were associated with protein synthesis and metabolism. Serine proteases were the main proteolytic enzymes detected. Some 40% of the cytochrome P450 genes and 74% of the glutathione S transferases (GSTs) in the genome of D. melanogaster were observed to be expressed in the midgut by oligoarray analysis. We also identified potential transcription factor binding motifs (TFBMs) of P450s, GSTs and carboxylesterases. Many of the midgut-expressed GST genes contained candidate TFBMs homologous to TFBMs in mammals that have been associated with responses to oxidative stress. We also investigated the response of GSTs in the midgut to dietary H2O2, which showed a dosage-based differential response.

Molecular characterization and responsive expression of a defender against apoptotic cell death homologue from the Hessian fly, Mayetiola destructor

Apoptosis or programmed cell death is an active process occurring in multicellular organisms to maintain growth and development. The Hessian fly, Mayetiola destructor, is rapidly emerging as a model insect species to study insect-plant interactions and to decipher some exceptional physiological phenomena. In this study, we report the characterization and expression profiles of a putative Hessian fly defender against apoptotic cell death (DAD1) homologue designated MdesDAD1. The deduced amino acid sequence of MdesDAD1 revealed significant similarity (75% identity, 9e-42) to other insect and non-insect DAD1 sequences. Phylogenetic analysis grouped MdesDAD1 within a sub-clade consisting of other insect DAD1 homologues. Quantitative analysis indicated constitutive levels of MdesDAD1 mRNA in all the tissues examined but an altered expression pattern during development, wherein the highest mRNA levels observed were prior to pupation. Most interestingly, MdesDAD1 transcript was found to be up-regulated during incompatible (larvae reared on resistant wheat) Hessian fly/wheat interactions compared to compatible (larvae reared on susceptible wheat) interactions. These results suggest MdesDAD1 to have a putative role in the inhibition of unwanted apoptosis triggered during development and in incompatible Hessian fly/wheat interactions. The results obtained provide clues to plausible insect and host-plant factors that could be responsible for the induction of MdesDAD1.

Antioxidant defense response in a galling insect

Herbivorous insect species are constantly challenged with reactive oxygen species (ROS) generated from endogenous and exogenous sources. ROS produced within insects because of stress and prooxidant allelochemicals produced by host plants in response to herbivory require a complex mode of antioxidant defense during insect/plant interactions. Some insect herbivores have a midgut-based defense against the suite of ROS encountered. Because the Hessian fly (Mayetiola destructor) is the major insect pest of wheat worldwide, and an emerging model for all gall midges, we investigated its antioxidant responses during interaction with its host plant. Quantitative data for two phospholipid glutathione peroxidases (MdesPHGPX-1 and MdesPHGPX-2), two catalases (MdesCAT-1 and MdesCAT-2), and two superoxide dismutases (MdesSOD-1 and MdesSOD-2) revealed high levels of all of the mRNAs in the midgut of larvae on susceptible wheat (compatible interaction). During development of the Hessian fly on susceptible wheat, a differential expression pattern was observed for all six genes. Analysis of larvae on resistant wheat (incompatible interaction) compared with larvae on susceptible wheat showed increased levels of mRNAs in larvae on resistant wheat for all of the antioxidant genes except MdesSOD-1 and MdesSOD-2. We postulate that the increased mRNA levels of MdesPHGPX-1, MdesPHGPX-2, MdesCAT-1, and MdesCAT-2 reflect responses to ROS encountered by larvae while feeding on resistant wheat seedlings and/or ROS generated endogenously in larvae because of stress/starvation. These results provide an opportunity to understand the cooperative antioxidant defense responses in the Hessian fly/wheat interaction and may be applicable to other insect/plant interactions.

cDNA cloning and transcriptional expression of a peritrophin-like gene in the Hessian fly, Mayetiola destructor [Say]

One of the well-studied components of the insect gut is the peritrophic matrix (PM). This semipermeable structure primarily functions in digestion, and protection against invasive microorganisms and mechanical damage. We report the cDNA cloning and transcription profiles of a peritrophin-A like gene (designated MdesPERI-A1) in the Hessian fly Mayetiola destructor. The predicted amino acid sequence of MdesPERI-A1 revealed a putative secretion signal peptide at its amino terminus, similarity to peritrophins from other insects including dipterans, and the presence of two chitin binding domains each containing six cysteine residues. Quantitative expression analysis of MdesPERI-A1 mRNA in different larval tissues revealed the transcript to be predominantly present in the midgut (597.9-fold) compared to other tissues assayed including salivary glands and fat bodies. Spatial expression patterns during development showed a peak expression of MdesPERI-A1 in the feeding second-instars (146-fold) and a decline in expression in the pupal and adult stages. Transcription profiling of MdesPERI-A1 during compatible (larvae on susceptible plants) and incompatible (larvae on resistant plants) interactions with wheat revealed a greater level (1.7-fold) of MdesPERI-A1 transcript in larvae on resistant plants in the initial time point examined. However, MdesPERI-A1 expression declined in larvae on resistant plants at the later time points.

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.

Expression patterns of antibacterial genes in the Hessian fly

We report on the transcriptional patterns of three antibacterial genes, a defensin (MdesDEF-1), a diptericin (MdesDIP-1) and a lysozyme (MdesLYS-1), during development in Hessian fly, Mayetiola destructor. Quantitative analysis by real-time PCR of mRNA levels in different tissues revealed a predominance of the transcripts for all three genes in the midgut, while analysis during development revealed greatest abundance in mRNA during the 3rd-instar. An evaluation of the midgut lumen revealed the presence of a diverse bacterial flora in larvae maintained on susceptible wheat. Further, the titer of bacteria in the midgut increased approximately 250-fold from the 1st-instar through the 2nd-instar. However, no detectable titer of bacteria was observed from the midgut lumen of larvae maintained on resistant plants. PCR amplicons produced using primers designed to conserved regions of the Pseudomonas 16S rRNA gene supported taxonomic identification for some of the bacteria comprising the midgut flora as belonging to the genus Pseudomonas. Analysis of mRNA for the Hessian fly antibacterial genes in larvae feeding on susceptible and resistant plants revealed an increase in the transcript level for MdesDEF-1 in 1st-instar larvae on susceptible plants, while the transcript levels for MdesDIP-1 and MdesLYS-1 were constant. Results suggest the transcriptional patterns of the Hessian fly antibacterial genes observed could be associated with the developing midgut bacterial flora present in larvae feeding on susceptible wheat as well as microbial challenge encountered at other stages in development.

Gene-for-gene defense of wheat against the Hessian fly lacks a classical oxidative burst

Genetic similarities between plant interactions with microbial pathogens and wheat interactions with Hessian fly larvae prompted us to investigate defense and counterdefense mechanisms. Plant oxidative burst, a rapid increase in the levels of active oxygen species (AOS) within the initial 24 h of an interaction with pathogens, commonly is associated with defenses that are triggered by gene-for-gene recognition events similar to those involving wheat and Hessian fly larvae. RNAs encoded by Hessian fly superoxide dismutase (SOD) and catalase (CAT) genes, involved in detoxification of AOS, increased in first-instar larvae during both compatible and incompatible interactions. However, mRNA levels of a wheat NADPH oxidase (NOX) gene that generates superoxide (O2-) did not increase. In addition, inhibiting wheat NOX enzyme with diphenyleneiodonium did not result in increased survival of avirulent larvae. However, nitro blue tetrazolium staining indicated that basal levels of O2- are present in both uninfested and infested wheat tissue. mRNA encoded by wheat genes involved in detoxification of the cellular environment, SOD, CAT, and glutathione-S-transferase did not increase in abundance. Histochemical staining with 3,3-diaminobenzidine revealed no increases in wheat hydrogen peroxide (H2O2) during infestation that were correlated with the changes in larval SOD and CAT mRNA. However, treatment with 2',7'-dichlorofluorescin demonstrated the presence of basal levels of H2O2 in the elongation zone of both infested and uninfested plants. The accumulation of a wheat flavanone 3-hydroxylase mRNA did show some parallels with larval gene mRNA profiles. These results suggested that larvae encounter stresses imposed by mechanisms other than an oxidative burst in wheat seedlings.

Characterization of a serine carboxypeptidase in the salivary glands and fat body of the orange wheat blossom midge, Sitodiplosis mosellana (Diptera: Cecidomyiidae)

A full-length cDNA encoding a serine carboxypeptidase (designated SmSCP-1) was recovered from an ongoing salivary gland EST project of the wheat midge. The deduced 461-amino acid sequence had a putative signal sequence at the amino terminus, indicating it was a secreted protein. The protein shared homology with serine carboxypeptidases from other insects, mammals, plants, and yeasts. SmSCP-1 mRNA was expressed in all stages of development and detected in salivary gland and fat body tissues but not in midgut tissue. Expression analysis and quantitative real-time PCR assays in male and female wheat midges and the fat body tissue of adult midges revealed that SmSCP-1 was up-regulated nearly four-fold in the female midges compared to males and nearly two-fold in female fat body compared to male fat body. The wheat midge serine carboxypeptidase (SmSCP-1) most likely has a dual function. As a secreted digestive enzyme, it could play a role in mobilizing host-plant seed reserves for feeding larvae and as expressed in fat body could function as an exopeptidase in degradation of vitellogenin and/or in post-translational processing of other enzymes.

Differential expression of two cytochrome P450 genes in compatible and incompatible Hessian fly/wheat interactions

We have recovered two Hessian fly cytochrome P450 cDNAs from an ongoing midgut EST project. CYP6AZ1 and CYP6BA1 represent two new subfamilies within the CYP6 family. The deduced amino acid sequences for CYP6AZ1 and CYP6BA1 show conserved structural and functional domains of insect P450s. Expression analysis with reverse transcription-polymerase chain reaction (RT-PCR) indicated that CYP6AZ1 is midgut specific and induced during active larval feeding, whereas CYP6BA1 was expressed in all tissues and developmental stages examined. Further expression analysis of CYP6AZ1 with RT-PCR in compatible and incompatible Hessian fly/wheat interactions suggested that CYP6AZ1 may be required for larval feeding in compatible interactions. These results should lead to a better understanding of the Hessian fly/wheat interaction with emphasis on the larval midgut as a critical interface with its host plant.

Molecular cloning and characterization of two digestive serine proteases from the Hessian fly, Mayetiola destructor

Full-length cDNA and genomic sequences for two genes (designated mdesprot-I and mdesprot-II) encoding digestive serine proteases in Hessian fly, Mayetiola destructor, have been cloned and characterized. The deduced amino acid sequences revealed similarity with trypsin-like digestive serine proteases from other Dipterans. Both mdesprot-I and mdesprot-II encoded proteins with secretion signal peptides at the N-terminals, indicating the proteins are secreted proteases that should function as midgut digestive proteases. A cytological analysis with fluorescent in situ hybridization revealed the cytological localization of mdesprot-I and mdesprot-II on the long arm of Autosome 2. Results are discussed in the context of the efficacy of potential protease inhibitors to develop Hessian fly resistant wheat through genetic engineering approaches.