Praveen Mamidala

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.

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

BackgroundBed 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.ResultsWe 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.ConclusionsWe 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.

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.

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.

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.

Identification and Validation of Reference Genes for Quantitative Real-Time Polymerase Chain Reaction in Cimex lectularius

ABSTRACT 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; &agr;-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.

Validation of reference genes for gene expression studies in the emerald ash borer (Agrilus planipennis)

Abstract  The Emerald ash borer (EAB, Agrilus planipennis Fairmaire) an exotic invasive insect pest has killed millions of ash trees (Fraxinus spp.) across North America and threatens billions more. We validated six A. planipennis reference genes (actin, ACT; beta tubulin, β‐TUB; glyceraldehyde‐3‐phosphate dehydrogenase, GAPDH; ribosomal protein, RPL7; translation elongation factor 1α, TEF‐1α; and ubiquitin, UBQ) using geNorm, Normfinder and BestKeeper for accurate determination of target messenger RNA levels in gene expression studies. The stability of the six reference genes was evaluated in different larval tissues, developmental stages and two treatments of A. planipennis using quantitative real‐time polymerase chain reaction. Although there was no consistent ranking observed among the reference genes across the samples, the overall analysis revealed TEF‐1α as the most stable reference gene. GAPDH and ACT showed least stability for all the samples studied. We conclude that TEF‐1α is the most appropriate reference gene for gene expression studies in A. planipennis. Results obtained can be applicable for transcript profiling in other invasive insect pests. Further, these validated reference genes could also serve as the basis for selection of candidate reference genes in any given insect system post‐validation.

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.

Evaluation of Reference Genes for Expression Studies in Ash (Fraxinus spp.)

Ash (Fraxinus spp.) is a dominant tree species in North America, in both managed and natural landscapes. However, due to the rapid invasion by the emerald ash borer (Agrilus planipennis), an exotic invasive insect pest, millions of North American ash trees have been killed. Real-time quantitative polymerase chain reaction (RTq-PCR) is widely used for validating transcript levels in gene expression studies for which a good reference gene is mandatory. In the current study, we evaluated the stability of ten reference genes in at least five different tissues (phloem, roots, shoots, immature leaves, and mature leaves), and two developmental stages (young and old) among three ash species including the resistant Asian Manchurian ash (F. mandshurica) and two susceptible North American ash species (green—F. pennsylvanica and white—F. americana). Of the examined genes, the translation elongation factor alpha (eEF1α) was observed to be most stable and thus is recommended for RTq-PCR based gene expression studies in Fraxinus species. To our knowledge, this is the first report on the stability of reference genes across ash species (in different tissues and during development).

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.