Pyrosequence analysis of expressed sequence tags for Manduca sexta hemolymph proteins involved in immune responses

Preferential binding of DAP-PGs by major peptidoglycan recognition proteins found in cell-free hemolymph of Manduca sexta

Peptidoglycan recognition proteins (PGRPs) detect invading bacteria to trigger or modulate immune responses in insects. While these roles are established in Drosophila, functional studies are not yet achieved at the PGRP family level in other insects. To attain this goal, we selected Manduca sexta PGRP12 and five of the nine secreted PGRPs for recombinant expression and biochemical characterization. We cloned PGRP2-5, 12 and 13 cDNAs, produced the proteins in full (PGRP2-5, 13) or in part (PGRP3s, 12e, 13N, 13C) in Sf9 cells, and tested their bindings of two muramyl pentapeptides by surface plasmon resonance, two soluble peptidoglycans by competitive ELISA, and four insoluble peptidoglycans and eight whole bacteria by a pull-down assay. Preferential binding of meso-diaminopimelic acid-peptidoglycans (DAP-PGs) was observed in all the proteins containing a peptidoglycan binding domain and, since PGRP6, 7 and 9 proteins were hardly detected in cell-free hemolymph, the reportoire of PGRPs (including PGRP1 published previously) in M. sexta hemolymph is likely adapted to mainly detect Gram-negative bacteria and certain Gram-positive bacteria with DAP-PGs located on their surface. After incubation with plasma from naïve larvae, PGRP2, 3f, 4, 5, 13f and 13N considerably stimulated prophenoloxidase activation in the absence of a bacterial elicitor. PGRP3s and 12e had much smaller effects. Inclusion of the full-length PGRPs and their regions in the plasma also led to proHP8 activation, supporting their connections to the Toll pathway, since HP8 is a Spӓtzle-1 processing enzyme in M. sexta. Together, these findings raised concerns on the common belief that the Toll-pathway is specific for Gram-positive bacteria in insects.

iTRAQ-Based Comparative Proteomic Analysis of Larval Midgut From the Beet Armyworm, Spodoptera exigua (Hübner) (Lepidoptera: Noctuidae) Challenged With the Entomopathogenic Bacteria Serratia marcescens

Entomopathogenic bacteria Serratia marcescens is widely used as an environmentally friendly biocontrol agent against various pests, including Spodoptera exigua. Understanding the immune defense mechanism of S. exigua through comparative proteomic analysis can identify the key proteins expressed in response to the microbial infection. Here, we employed the as isobaric tags for relative and absolute quantification (iTRAQ) technique to investigate the effects of S. marcescens on the proteomic expression of S. exigua. Based on the molecular functional analysis, the differentially expressed proteins (DEPs) were mainly involved in the binding process and catalytic activities. Further bioinformatics analysis revealed important DEPs that played a crucial role in innate immunity of S. exigua with recognition (C-type lectin), melanization (propanol oxidase 3, serine protease, Serine-type carboxypeptidase activity, clip domain serine protease 4), antimicrobial activity (lysozyme, lysozyme-like, gloverin, cecropin B), detoxification (acetyl-CoA C-acetyltransferase, 3-dehydroecdysone 3-alpha-reductase, glucuronosyltransferase, glutathione S-transferase) and others. The Quantitative real-time PCR (qRT-PCR) results further indicated the significant upregulation of the immune-related genes in Spodoptera exigua following S. marcescens infection. To the best of our knowledge, this is the first iTRAQ based study to characterize S. marcescens mediated proteomic changes in S. exigua and identified important immune-related DEPs. The results of this study will provide an essential resource for understanding the host-pathogen interactions and the development of novel microbial biopesticides against various pests.

Manduca sexta serpin-12 controls the prophenoloxidase activation system in larval hemolymph

Insect prophenoloxidase activation is coordinated by a serine protease network, which is regulated by serine protease inhibitors of the serpin superfamily. The enzyme system also leads to proteolytic processing of a Spätzle precursor. Binding of Spätzle to a Toll receptor turns on a signaling pathway to induce the synthesis of defense proteins. Previous studies of the tobacco hornworm Manduca sexta have revealed key members of the protease cascade, which generates phenoloxidase for melanogenesis and Spätzle to induce immunity-related genes. Here we provide evidence that M. sexta serpin-12 regulates hemolymph protease-14 (HP14), an initiating protease of the cascade. This inhibitor, unlike the other serpins characterized in M. sexta, has an amino-terminal extension rich in hydrophilic residues and an unusual P1 residue (Leu⁴²⁹) right before the scissile bond cleaved by a target protease. Serpins with similarities to serpin-12, including Drosophila Necrotic, were identified in a wide range of insects including flies, moths, wasps, beetles, and two hemimetabolous species. The serpin-12 mRNA is present at low, constitutive levels in larval fat body and hemocytes and becomes more abundant after an immune challenge. We produced the serpin-12 core domain (serpin-12ΔN) in insect cells and in Escherichia coli and demonstrated its inhibition of human cathepsin G, bovine α-chymotrypsin, and porcine pancreatic elastase. MALDI-TOF analysis of the reaction mixtures confirmed the predicted P1 residue of Leu⁴²⁹. Supplementation of larval plasma samples with the serpin-12ΔN decreased prophenoloxidase activation elicited by microbial cells and reduced the proteolytic activation of the protease precursors of HP6, HP8, PAPs, and other serine protease-related proteins. After incubation of plasma stimulated with peptidoglycan, a 72 kDa protein appeared, which was recognized by polyclonal antibodies against both serpin-12 and HP14, suggesting that a covalent serpin-protease complex formed when serpin-12 inhibited HP14. Together, these data suggest that M. sexta serpin-12 inhibits HP14 to regulate melanization and antimicrobial peptide induction.

Serpin-9 and -13 regulate hemolymph proteases during immune responses of Manduca sexta

Serpins are a superfamily of proteins, most of which inhibit cognate serine proteases by forming inactive acyl-enzyme complexes. In the tobacco hornworm Manduca sexta, serpin-1, -3 through -7 negatively regulate a hemolymph serine protease system that activates precursors of the serine protease homologs (SPHs), phenoloxidases (POs), Spätzles, and other cytokines. Here we report the cloning and characterization of M. sexta serpin-9 and -13. Serpin-9, a 402-residue protein most similar to Drosophila Spn77Ba, has R³⁶⁶ at the P1 position right before the cleavage site; Serpin-13, a 444-residue ortholog of Drosophila Spn28Dc, is longer than the other seven serpins and has R⁴¹⁰ as the P1 residue. Both serpins are mainly produced in fat body and secreted into plasma to function. While their mRNA and protein levels were not up-regulated upon immune challenge, they blocked protease activities and affected proPO activation in hemolymph. Serpin-9 inhibited human neutrophil elastase, cathepsin G, trypsin, and chymotrypsin to different extents; serpin-13 reduced trypsin activity to approximately 10% at a molar ratio of 4:1 (serpin: enzyme). Serpin-9 was cleaved at Arg³⁶⁶ by the enzymes with different specificity, but serpin-13 had four P1 sites (Arg⁴¹⁰ for trypsin-like proteases, Gly⁴⁰⁶ and Ala⁴⁰⁹ for the elastase and Thr⁴⁰⁴ for cathepsin G). Supplementation of induced cell-free hemolymph (IP, P for plasma) with recombinant serpin-9 did not noticeably affect proPO activation, but slightly reduced the PO activity increase after 0-50% ammonium sulfate fraction of the IP had been elicited by bacteria. In comparison, addition of recombinant serpin-13 significantly inhibited proPO activation in IP and the suppression was stronger in the fraction of IP. Serpin-9- and -13-containing protein complexes were isolated from IP using their antibodies. Hemolymph protease-1 precursor (proHP1), HP6 and HP8 were found to be associated with serpin-9, whereas proHP1, HP2 and HP6 were pulled downed with serpin-13. These results indicate that both serpins regulate immune proteases in hemolymph of M. sexta larvae.

An analysis of 67 RNA-seq datasets from various tissues at different stages of a model insect, Manduca sexta

Background

Manduca sexta is a large lepidopteran insect widely used as a model to study biochemistry of insect physiological processes. As a part of its genome project, over 50 cDNA libraries have been analyzed to profile gene expression in different tissues and life stages. While the RNA-seq data were used to study genes related to cuticle structure, chitin metabolism and immunity, a vast amount of the information has not yet been mined for understanding the basic molecular biology of this model insect. In fact, the basic features of these data, such as composition of the RNA-seq reads and lists of library-correlated genes, are unclear. From an extended view of all insects, clear-cut tempospatial expression data are rarely seen in the largest group of animals including Drosophila and mosquitoes, mainly due to their small sizes.

Results

We obtained the transcriptome data, analyzed the raw reads in relation to the assembled genome, and generated heatmaps for clustered genes. Library characteristics (tissues, stages), number of mapped bases, and sequencing methods affected the observed percentages of genome transcription. While up to 40% of the reads were not mapped to the genome in the initial Cufflinks gene modeling, we identified the causes for the mapping failure and reduced the number of non-mappable reads to <8%. Similarities between libraries, measured based on library-correlated genes, clearly identified differences among tissues or life stages. We calculated gene expression levels, analyzed the most abundantly expressed genes in the libraries. Furthermore, we analyzed tissue-specific gene expression and identified 18 groups of genes with distinct expression patterns.

Conclusion

We performed a thorough analysis of the 67 RNA-seq datasets to characterize new genomic features of M. sexta. Integrated knowledge of gene functions and expression features will facilitate future functional studies in this biochemical model insect.

Electronic supplementary material

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

Pathogen-induced maternal effects result in enhanced immune responsiveness across generations

Parental investment theory postulates that adults can accurately perceive cues from their surroundings, anticipate the needs of future offspring based on those cues, and selectively allocate nongenetic resources to their progeny. Such context‐dependent parental contributions can result in phenotypically variable offspring. Consistent with these predictions, we show that bacterially exposed Manduca sexta mothers oviposited significantly more variable embryos (as measured by mass, volume, hatching time, and hatching success) relative to naïve and control mothers. By using an in vivo “clearance of infection” assay, we also show that challenged larvae born to heat‐killed‐ or live‐Serratia‐injected mothers, supported lower microbial loads and cleared the infection faster than progeny of control mothers. Our data support the notion that mothers can anticipate the future pathogenic risks and immunological needs of their unborn offspring, providing progeny with enhanced immune protection likely through transgenerational immune priming. Although the inclusion of live Serratia into oocytes does not appear to be the mechanism by which mothers confer protection to their young, other mechanisms, including epigenetic modifications in the progeny due to maternal pathogenic stress, may be at play. The adaptive nature of maternal effects in the face of pathogenic stress provides insights into parental investment, resource allocation, and life‐history theories and highlights the significant role that pathogen‐induced maternal effects play as generators and modulators of evolutionary change.

Transcriptome profiling of Diachasmimorpha longicaudata towards useful molecular tools for population management

Background

Diachasmimorpha longicaudata (Hymenoptera: Braconidae) is a solitary parasitoid of Tephritidae (Diptera) fruit flies of economic importance currently being mass-reared in bio-factories and successfully used worldwide. A peculiar biological aspect of Hymenoptera is its haplo-diploid life cycle, where females (diploid) develop from fertilized eggs and males (haploid) from unfertilized eggs. Diploid males were described in many species and recently evidenced in D. longicaudata by mean of inbreeding studies. Sex determination in this parasitoid is based on the Complementary Sex Determination (CSD) system, with alleles from at least one locus involved in early steps of this pathway. Since limited information is available about genetics of this parasitoid species, a deeper analysis on D. longicaudata’s genomics is required to provide molecular tools for achieving a more cost effective production under artificial rearing conditions.

Results

We report here the first transcriptome analysis of male-larvae, adult females and adult males of D. longicaudata using 454-pyrosequencing. A total of 469766 reads were analyzed and 8483 high-quality isotigs were assembled. After functional annotation, a total of 51686 unigenes were produced, from which, 7021 isotigs and 20227 singletons had at least one BLAST hit against the NCBI non-redundant protein database. A preliminary comparison of adult female and male evidenced that 98 transcripts showed differential expression profiles, with at least a 10-fold difference. Among the functionally annotated transcripts we detected four sequences potentially involved in sex determination and three homologues to two known genes involved in the sex determination cascade. Finally, a total of 4674SimpleSequence Repeats (SSRs) were in silico identified and characterized.

Conclusion

The information obtained here will significantly contribute to the development of D. longicaudata functional genomics, genetics and population-based genome studies. Thousands of new microsatellite markers were identified as toolkits for population genetics analysis. The transcriptome characterized here is the starting point to elucidate the molecular bases of the sex determination mechanism in this species.

Electronic supplementary material

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

Pyrosequencing the Midgut Transcriptome of the Banana Weevil Cosmopolites sordidus (Germar) (Coleoptera: Curculionidae) Reveals Multiple Protease-Like Transcripts

The banana weevil Cosmopolites sordidus is an important and serious insect pest in most banana and plantain-growing areas of the world. In spite of the economic importance of this insect pest very little genomic and transcriptomic information exists for this species. In the present study, we characterized the midgut transcriptome of C. sordidus using massive 454-pyrosequencing. We generated over 590,000 sequencing reads that assembled into 30,840 contigs with more than 400 bp, representing a significant expansion of existing sequences available for this insect pest. Among them, 16,427 contigs contained one or more GO terms. In addition, 15,263 contigs were assigned an EC number. In-depth transcriptome analysis identified genes potentially involved in insecticide resistance, peritrophic membrane biosynthesis, immunity-related function and defense against pathogens, and Bacillus thuringiensis toxins binding proteins as well as multiple enzymes involved with protein digestion. This transcriptome will provide a valuable resource for understanding larval physiology and for identifying novel target sites and management approaches for this important insect pest.

Changes in the Plasma Proteome of Manduca sexta Larvae in Relation to the Transcriptome Variations after an Immune Challenge: Evidence for High Molecular Weight Immune Complex Formation

Manduca sextais a lepidopteran model widely used to study insect physiological processes, including innate immunity. In this study, we explored the proteomes of cell-free hemolymph from larvae injected with a sterile buffer (C for control) or a mixture of bacteria (I for induced). Of the 654 proteins identified, 70 showed 1.67 to >200-fold abundance increases after the immune challenge; 51 decreased to 0-60% of the control levels. While there was no strong parallel between plasma protein levels and their transcript levels in hemocytes or fat body, the mRNA level changes (i.e.I/C ratios of normalized read numbers) in the tissues concurred with their protein level changes (i.e.I/C ratios of normalized spectral counts) with correlation coefficients of 0.44 and 0.57, respectively. Better correlations support that fat body contributes a more significant portion of the plasma proteins involved in various aspects of innate immunity. Consistently, ratios of mRNA and protein levels were better correlated for immunity-related proteins than unrelated ones. There is a set of proteins whose apparent molecular masses differ considerably from the calculatedMr's, suggestive of posttranslational modifications. In addition, some lowMrproteins were detected in the range of 80 to >300 kDa on a reducing SDS-polyacrylamide gel, indicating the existence of highMrcovalent complexes. We identified 30 serine proteases and their homologs, 11 of which are known members of an extracellular immune signaling network. Along with our quantitative transcriptome data, the protein identification, inducibility, and association provide leads toward a focused exploration of humoral immunity inM. sexta.

Integrated modeling of protein-coding genes in the Manduca sexta genome using RNA-Seq data from the biochemical model insect

The genome sequence of Manduca sexta was recently determined using 454 technology. Cufflinks and MAKER2 were used to establish gene models in the genome assembly based on the RNA-Seq data and other species' sequences. Aided by the extensive RNA-Seq data from 50 tissue samples at various life stages, annotators over the world (including the present authors) have manually confirmed and improved a small percentage of the models after spending months of effort. While such collaborative efforts are highly commendable, many of the predicted genes still have problems which may hamper future research on this insect species. As a biochemical model representing lepidopteran pests, M. sexta has been used extensively to study insect physiological processes for over five decades. In this work, we assembled Manduca datasets Cufflinks 3.0, Trinity 4.0, and Oases 4.0 to assist the manual annotation efforts and development of Official Gene Set (OGS) 2.0. To further improve annotation quality, we developed methods to evaluate gene models in the MAKER2, Cufflinks, Oases and Trinity assemblies and selected the best ones to constitute MCOT 1.0 after thorough crosschecking. MCOT 1.0 has 18,089 genes encoding 31,666 proteins: 32.8% match OGS 2.0 models perfectly or near perfectly, 11,747 differ considerably, and 29.5% are absent in OGS 2.0. Future automation of this process is anticipated to greatly reduce human efforts in generating comprehensive, reliable models of structural genes in other genome projects where extensive RNA-Seq data are available.

Identification and profiling of Manduca sexta microRNAs and their possible roles in regulating specific transcripts in fat body, hemocytes, and midgut

Significance of microRNA-mediated posttranscriptional regulation has been appreciated ever since its discovery. In the tobacco hornworm Manduca sexta, 164 conserved and 16 novel microRNAs have been identified experimentally (Zhang et al., 2012, 2014). To extend the list of microRNAs in this lepidopteran model species and further explore their possible regulatory roles, we constructed and sequenced small RNA libraries of M. sexta fat body, hemocytes and midgut, since transcriptomes of these tissues from the 5th instar larvae had been studied quite extensively. Each library represented a mixture of the same tissues from larvae that were naïve or induced by three different pathogens. From a total of 167 million reads obtained, we identified two new variants of conserved miR-281 and miR-305 and six novel microRNAs. Abundances of all microRNAs were normalized and compared to reveal their differential expression in these three tissues. Star strands of ten microRNAs were present at higher levels than the corresponding mature strands. From a list of tissue-specific transcripts, we predicted target sites in 3'-UTRs using preferentially expressed microRNA groups in each tissue and suggested possible regulatory roles of these microRNAs in energy metabolism, insecticide resistance, and some mitochondrial and immune gene expression. Examining manifold targets, microRNA regulations were suggested of multiple physiological processes. This study has enriched our knowledge of M. sexta microRNAs and how microRNAs potentially coordinate different physiological processes.

Sequence conservation, phylogenetic relationships, and expression profiles of nondigestive serine proteases and serine protease homologs in Manduca sexta

Serine protease (SP) and serine protease homolog (SPH) genes in insects encode a large family of proteins involved in digestion, development, immunity, and other processes. While 68 digestive SPs and their close homologs are reported in a companion paper (Kuwar et al., in preparation), we have identified 125 other SPs/SPHs in Manduca sexta and studied their structure, evolution, and expression. Fifty-two of them contain cystine-stabilized structures for molecular recognition, including clip, LDLa, Sushi, Wonton, TSP, CUB, Frizzle, and SR domains. There are nineteen groups of genes evolved from relatively recent gene duplication and sequence divergence. Thirty-five SPs and seven SPHs contain 1, 2 or 5 clip domains. Multiple sequence alignment and molecular modeling of the 54 clip domains have revealed structural diversity of these regulatory modules. Sequence comparison with their homologs in Drosophila melanogaster, Anopheles gambiae and Tribolium castaneum allows us to classify them into five subfamilies: A are SPHs with 1 or 5 group-3 clip domains, B are SPs with 1 or 2 group-2 clip domains, C, D1 and D2 are SPs with a single clip domain in group-1a, 1b and 1c, respectively. We have classified into six categories the 125 expression profiles of SP-related proteins in fat body, brain, midgut, Malpighian tubule, testis, and ovary at different stages, suggesting that they participate in various physiological processes. Through RNA-Seq-based gene annotation and expression profiling, as well as intragenomic sequence comparisons, we have established a framework of information for future biochemical research of nondigestive SPs and SPHs in this model species.

20-hydroxyecdysone transcriptionally regulates humoral immunity in the fat body of Helicoverpa armigera

20-hydroxyecdysone (20E) increases its titre level during the wandering stage and influences innate immunity in many holometabolous insects. However, the function of 20E as an immune-activator or -suppressor needs to be determined. Here, the transcriptome of the peptidoglycan-challenged fat body of the cotton bollworm, Helicoverpa armigera, was analysed using Illumina sequencing technology. Overall, 32 073 unigenes were assembled with a mean length of 643 nucleotides. Gene expression dynamics in the fat body during the wandering stage and of peptidoglycan-challenged individuals were investigated by the digital gene expression system. Pattern recognition receptors [such as peptidoglycan recognition protein B (PGRP B), PGRP S2 precursor, C-type lectin 5, hemolin and β-1,3-glucan recognition protein 2a] and antimicrobial peptides (namely attacin, gloverin, gloverin precursor, gloverin-like, cecropin 2, cecropin D, cecropin D-like and i-type lysozyme) significantly increased their mRNA levels during the wandering stage. 20E treatment significantly induced the expression of these genes. Antibacterial activities were also enhanced during the wandering stage and after 20E injections. Bacillus subtilis peptidoglycan induced the expression of PGRP D, PGRP B, PGRP S2 precursor, gloverin, gloverin precursor, gloverin-like, cecropin 2, cecropin D and lebocin-like genes. These results demonstrate that 20E acts by enhancing humoral immunity in H. armigera.

Transcriptome profiling of pyrethroid resistant and susceptible mosquitoes in the malaria vector, Anopheles sinensis

BACKGROUND

Anopheles sinensis is a major malaria vector in China and other Southeast Asian countries, and it is becoming increasingly resistant to the insecticides used for agriculture, net impregnation, and indoor residual spray. Very limited genomic information on this species is available, which has hindered the development of new tools for resistance surveillance and vector control. We used the 454 GS FLX system and generated expressed sequence tag (EST) databases of various life stages of An. sinensis, and we determined the transcriptional differences between deltamethrin resistant and susceptible mosquitoes.

RESULTS

The 454 GS FLX transcriptome sequencing yielded a total of 624,559 reads (average length of 290 bp) with the pooled An. sinensis mosquitoes across various development stages. The de novo assembly generated 33,411 contigs with average length of 493 bp. A total of 8,057 ESTs were generated with Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation. A total of 2,131 ESTs were differentially expressed between deltamethrin resistant and susceptible mosquitoes collected from the same field site in Jiangsu, China. Among these differentially expressed ESTs, a total of 294 pathways were mapped to the KEGG database, with the predominant ESTs belonging to metabolic pathways. Furthermore, a total of 2,408 microsatellites and 15,496 single nucleotide polymorphisms (SNPs) were identified.

CONCLUSIONS

The annotated EST and transcriptome databases provide a valuable genomic resource for further genetic studies of this important malaria vector species. The differentially expressed ESTs associated with insecticide resistance identified in this study lay an important foundation for further functional analysis. The identified microsatellite and SNP markers will provide useful tools for future population genetic and comparative genomic analyses of malaria vectors.

Semi-quantitative analysis of changes in the plasma peptidome of Manduca sexta larvae and their correlation with the transcriptome variations upon immune challenge

The tobacco hornworm, Manduca sexta, has been used as a biochemical model for studying insect physiological processes. While the transcriptomes of its fat body, hemocytes, midgut, and antennae have been examined in several studies, limited information is available for proteins in tissues, cells, or body fluids of this insect. In keeping pace with the M. sexta genome project, we launched a pilot study to identify differences in the peptidome of cell-free hemolymph samples from larvae injected with buffer or a mixture of bacteria. At 24 h after injection, plasma was collected and treated with 50% acetonitrile to precipitate large proteins. The supernatants, containing peptides (<25 kDa) and other stable proteins (>25 kDa), were digested with trypsin and analyzed by nano-liquid chromatography and nano-electrospray tandem mass spectrometry (nanoLC-MS/MS) on an LTQ Orbitrap XL mass spectrometer. Known M. sexta cDNA sequences and gene transcripts from the draft genome were translated in silico to generate a database of polypeptides (i.e. peptides and proteins) in this species. By searching the database, we identified 268 hemolymph polypeptides, 50 of which showed 1.67-200 fold abundance increases after the immune challenge, as judged by significant changes in normalized spectral counts between the control and induced plasma. These included a total of 33 antimicrobial peptides (attacins, cecropins, defensins, diapausins, gallerimycin, gloverin, lebocins, lysozymes), pattern recognition receptors, and proteinase inhibitors. Although there was no strong parallel (correlation coefficients: -0.13, 0.11, 0.39 and 0.62) between plasma peptide levels and their transcript levels in control or induced hemocytes or fat body, we observed the mRNA level changes in hemocytes and fat body concurred with their peptide level changes with correlation coefficients of 0.67 and 0.76, respectively. These data suggest that fat body contributed a significant portion of the plasma polypeptides involved in various aspects of innate immunity after the bacterial injection.

Characterization of intestinal bacteria in wild and domesticated adult black tiger shrimp (Penaeus monodon)

The black tiger shrimp (Penaeus monodon) is a marine crustacean of economic importance in the world market. To ensure sustainability of the shrimp industry, production capacity and disease outbreak prevention must be improved. Understanding healthy microbial balance inside the shrimp intestine can provide an initial step toward better farming practice and probiotic applications. In this study, we employed a barcode pyrosequencing analysis of V3-4 regions of 16S rRNA genes to examine intestinal bacteria communities in wild-caught and domesticated P. monodon broodstock. Shrimp faeces were removed from intestines prior to further analysis in attempt to identify mucosal bacterial population. Five phyla, Actinobacteria, Fusobacteria, Proteobacteria, Firmicutes and Bacteroidetes, were found in all shrimp from both wild and domesticated environments. The operational taxonomic unit (OTU) was assigned at 97% sequence identity, and our pyrosequencing results identified 18 OTUs commonly found in both groups. Sequences of the shared OTUs were similar to bacteria in three phyla, namely i) Proteobacteria (Vibrio, Photobacterium, Novosphingobium, Pseudomonas, Sphingomonas and Undibacterium), ii) Firmicutes (Fusibacter), and iii) Bacteroidetes (Cloacibacterium). The shared bacterial members in P. monodon from two different habitats provide evidence that the internal environments within the host shrimp also exerts selective pressure on bacterial members. Intestinal bacterial profiles were compared using denaturing gradient gel electrophoresis (DGGE). The sequences from DGGE bands were similar to those of Vibrio and Photobacterium in all shrimp, consistent with pyrosequencing results. This work provides the first comprehensive report on bacterial populations in the intestine of adult black tiger shrimp and reveals some similar bacterial members between the intestine of wild-caught and domesticated shrimp.

Transcriptome analysis in cotton boll weevil (Anthonomus grandis) and RNA interference in insect pests

Cotton plants are subjected to the attack of several insect pests. In Brazil, the cotton boll weevil, Anthonomus grandis, is the most important cotton pest. The use of insecticidal proteins and gene silencing by interference RNA (RNAi) as techniques for insect control are promising strategies, which has been applied in the last few years. For this insect, there are not much available molecular information on databases. Using 454-pyrosequencing methodology, the transcriptome of all developmental stages of the insect pest, A. grandis, was analyzed. The A. grandis transcriptome analysis resulted in more than 500.000 reads and a data set of high quality 20,841 contigs. After sequence assembly and annotation, around 10,600 contigs had at least one BLAST hit against NCBI non-redundant protein database and 65.7% was similar to Tribolium castaneum sequences. A comparison of A. grandis, Drosophila melanogaster and Bombyx mori protein families’ data showed higher similarity to dipteran than to lepidopteran sequences. Several contigs of genes encoding proteins involved in RNAi mechanism were found. PAZ Domains sequences extracted from the transcriptome showed high similarity and conservation for the most important functional and structural motifs when compared to PAZ Domains from 5 species. Two SID-like contigs were phylogenetically analyzed and grouped with T. castaneum SID-like proteins. No RdRP gene was found. A contig matching chitin synthase 1 was mined from the transcriptome. dsRNA microinjection of a chitin synthase gene to A. grandis female adults resulted in normal oviposition of unviable eggs and malformed alive larvae that were unable to develop in artificial diet. This is the first study that characterizes the transcriptome of the coleopteran, A. grandis. A new and representative transcriptome database for this insect pest is now available. All data support the state of the art of RNAi mechanism in insects.

Immune Signaling and Antimicrobial Peptide Expression in Lepidoptera

Many lepidopteran insects are agricultural pests that affect stored grains, food and fiber crops. These insects have negative ecological and economic impacts since they lower crop yield, and pesticides are expensive and can have off-target effects on beneficial arthropods. A better understanding of lepidopteran immunity will aid in identifying new targets for the development of specific insect pest management compounds. A fundamental aspect of immunity, and therefore a logical target for control, is the induction of antimicrobial peptide (AMP) expression. These peptides insert into and disrupt microbial membranes, thereby promoting pathogen clearance and insect survival. Pathways leading to AMP expression have been extensively studied in the dipteran Drosophila melanogaster. However, Diptera are an important group of pollinators and pest management strategies that target their immune systems is not recommended. Recent advances have facilitated investigation of lepidopteran immunity, revealing both conserved and derived characteristics. Although the general pathways leading to AMP expression are conserved, specific components of these pathways, such as recognition proteins have diverged. In this review we highlight how such comparative immunology could aid in developing pest management strategies that are specific to agricultural insect pests.

Nutrigenomics in Arma chinensis: transcriptome analysis of Arma chinensis fed on artificial diet and Chinese oak silk moth Antheraea pernyi pupae

Background

The insect predator, Arma chinensis, is capable of effectively controlling many pests, such as Colorado potato beetle, cotton bollworm, and mirid bugs. Our previous study demonstrated several life history parameters were diminished for A. chinensis reared on an artificial diet compared to a natural food source like the Chinese oak silk moth pupae. The molecular mechanisms underlying the nutritive impact of the artificial diet on A. chinensis health are unclear. So we utilized transcriptome information to better understand the impact of the artificial diet on A. chinensis at the molecular level.

Methodology/Principal Findings

Illumina HiSeq2000 was used to sequence 4.79 and 4.70 Gb of the transcriptome from pupae-fed and artificial diet-fed A. chinensis libraries, respectively, and a de novo transcriptome assembly was performed (Trinity short read assembler). This resulted in 112,029 and 98,724 contigs, clustered into 54,083 and 54,169 unigenes for pupae-fed and diet-fed A. chinensis, respectively. Unigenes from each sample’s assembly underwent sequence splicing and redundancy removal to acquire non-redundant unigenes. We obtained 55,189 unigenes of A. chinensis, including 12,046 distinct clusters and 43,143 distinct singletons. Unigene sequences were aligned by BLASTx to nr, Swiss-Prot, KEGG and COG (E-value <10⁻⁵), and further aligned by BLASTn to nt (E-value <10⁻⁵), retrieving proteins of highest sequence similarity with the given unigenes along with their protein functional annotations. Totally, 22,964, 7,898, 18,069, 15,416, 8,066 and 5,341 unigenes were annotated in nr, nt, Swiss-Prot, KEGG, COG and GO, respectively. We compared gene expression variations and found thousands of genes were differentially expressed between pupae-fed and diet-fed A. chinensis.

Conclusions/Significance

Our study provides abundant genomic data and offers comprehensive sequence information for studying A. chinensis. Additionally, the physiological roles of the differentially expressed genes enable us to predict effects of some dietary ingredients and subsequently propose formulation improvements to artificial diets.

A comprehensive analysis of the Manduca sexta immunotranscriptome

As a biochemical model, Manduca sexta has substantially contributed to our knowledge on insect innate immunity. The RNA-Seq approach was implemented in three studies to examine tissue immunotranscriptomes of this species. With the latest and largest focusing on highly regulated process- and tissue-specific genes, we further analyzed the same set of data using BLAST2GO to explore functional aspects of the larval fat body (F) and hemocyte (H) transcriptomes with (I) or without (C) immune challenge. Using immunity-related sequences from other insects, we found 383 homologous contigs and compared them with those discovered based on relative abundance changes. The major overlap of the two lists validated our previous research designed for gene discovery and transcript profiling in organisms lacking sequenced genomes. By concatenating the contigs, we established a repertoire of 232 immunity-related genes encoding proteins for pathogen recognition (16%), signal transduction (53%), microbe killing (13%) and others (18%). We examined their transcript levels along with attribute classifications and detected prominent differences in nine of the 30 level 2 gene ontology (GO) categories. The increase in extracellular proteins (155%) was consistent with the highly induced synthesis of defense molecules (e.g., antimicrobial peptides) in fat body after the immune challenge. We identified most members of the putative Toll, IMD, MAPK-JNK-p38 and JAK-STAT pathways and small changes in their mRNA levels. Together, these findings set the stage for on-going analysis of the M. sexta immunogenome.

Identification of plasma proteinase complexes with serpin-3 in Manduca sexta

Extracellular serine proteinase cascades stimulate prophenoloxidase (proPO) activation and antimicrobial peptide production in insect innate immune responses. Serpins in plasma regulate such cascades by selective inhibition of proteinases, in reactions which result in the formation of covalent serpin-proteinase complexes. We carried out experiments to identify plasma proteinases that are inhibited by Manduca sexta serpin-3, an immune-inducible serpin known to regulate proPO activation. Immunoaffinity chromatography, using antiserum to serpin-3, yielded serpin-3 complexes with proteinases identified by immunoblot analysis as prophenoloxidase-activating proteinase (PAP)-1, PAP-2, PAP-3, and hemolymph proteinase 8 (HP8). HP8 can cleave and activate the Toll ligand, Spätzle, leading to synthesis of antimicrobial peptides. Analysis by mass spectrometry of tryptic peptides derived from the serpin-3 complexes confirmed the presence of PAP-1, PAP-3, and HP8. Purified recombinant serpin-3 and active HP8 formed an SDS-stable complex in vitro. Identification of serpin-3-proteinase complexes in plasma provides insight into proteinase targets of serpin-3 and extends the understanding of serpin/proteinase function in the immune response of M. sexta.

Transcriptome of the plant virus vector Graminella nigrifrons, and the molecular interactions of maize fine streak rhabdovirus transmission

BACKGROUND

Leafhoppers (HEmiptera: Cicadellidae) are plant-phloem feeders that are known for their ability to vector plant pathogens. The black-faced leafhopper (Graminella nigrifrons) has been identified as the only known vector for the Maize fine streak virus (MFSV), an emerging plant pathogen in the Rhabdoviridae. Within G. nigrifrons populations, individuals can be experimentally separated into three classes based on their capacity for viral transmission: transmitters, acquirers and non-acquirers. Understanding the molecular interactions between vector and virus can reveal important insights in virus immune defense and vector transmission.

RESULTS

RNA sequencing (RNA-Seq) was performed to characterize the transcriptome of G. nigrifrons. A total of 38,240 ESTs of a minimum 100 bp were generated from two separate cDNA libraries consisting of virus transmitters and acquirers. More than 60% of known D. melanogaster, A. gambiae, T. castaneum immune response genes mapped to our G. nigrifrons EST database. Real time quantitative PCR (RT-qPCR) showed significant down-regulation of three genes for peptidoglycan recognition proteins (PGRP - SB1, SD, and LC) in G. nigrifrons transmitters versus control leafhoppers.

CONCLUSIONS

Our study is the first to characterize the transcriptome of a leafhopper vector species. Significant sequence similarity in immune defense genes existed between G. nigrifrons and other well characterized insects. The down-regulation of PGRPs in MFSV transmitters suggested a possible role in rhabdovirus transmission. The results provide a framework for future studies aimed at elucidating the molecular mechanisms of plant virus vector competence.

Global Transcriptional Analysis of Olfactory Genes in the Head of Pine Shoot Beetle, Tomicus yunnanensis

The most important proteins involved in olfaction include odorant binding protein (OBP), chemosensory protein (CSP), olfactory receptor (OR), and gustatory receptor (GR). Despite that the exhaustive genomic analysis has revealed a large number of olfactory genes in a number of model insects, it is still poorly understood for most nonmodel species. This is mostly due to the reason that the small antenna is challenging for collection. We can generally isolate one or few genes at a time by means of the traditional method. Here, we present the large-scale identifying members of the main olfactory genes from the head of Tomicus yunnanensis using Illumina sequencing. In a single run, we obtained over 51.8 million raw reads. These reads were assembled into 57,142 unigenes. Nearly 29,384 of them were functionally annotated in the NCBI nonredundant database. By depth analysis of the data, 11 OBPs, 8 CSPs, 18 ORs, and 8 GRs were retrieved. Sequences encoding full length proteins were further characterised for one OBP and two CSPs. The obtained olfactory genes provide a major resource in further unraveling the molecular mechanisms of T. yunnanensis chemoperception. This study indicates that the next generation sequencing is an attractive approach for efficient identification of olfactory genes from insects, for which the genome sequence is unavailable.

Identification and developmental profiling of conserved and novel microRNAs in Manduca sexta

MicroRNAs (miRNAs) are a group of small RNAs involved in translation inhibition or mRNA degradation. Due to its large size, Manduca sexta has long been used as a model to study insect physiology and biochemistry. While transcriptome studies have greatly enriched our knowledge on M. sexta structural genes, little is known about posttranscriptional regulation by miRNAs in this lepidopteran species. We constructed four small RNA libraries from embryos, 4th instar feeding larvae, pupae, and adults, obtained 21 million reads of 18-31 nucleotides by Illumina sequencing, and found 163 conserved and 13 novel miRNAs. By searching the M. sexta genome assembly, we identified precursors of 82 conserved miRNAs, 76 of which had mapped reads in one or more of these libraries. After normalization, we compared numbers of miRNA and miRNA-star reads in these libraries and observed abundance changes during development. Interestingly, mse-miR-281-star, mse-miR-31-star, mse-miR-965-star, mse-miR-9a-star, mse-miR-9b-star, mse-miR-2a-star, mse-miR-92b-star and mse-miR-279c-star are either more abundant or maintained at similar levels compared to respective mature miRNA strand. Expression profiling of the first set of miRNAs provided insights to their possible involvement in developmental regulation. This study will aid in the annotation of miRNA genes in the genome.

Pyrosequencing the Bemisia tabaci transcriptome reveals a highly diverse bacterial community and a robust system for insecticide resistance

Background

Bemisia tabaci (Gennadius) is a phloem-feeding insect poised to become one of the major insect pests in open field and greenhouse production systems throughout the world. The high level of resistance to insecticides is a main factor that hinders continued use of insecticides for suppression of B. tabaci. Despite its prevalence, little is known about B. tabaci at the genome level. To fill this gap, an invasive B. tabaci B biotype was subjected to pyrosequencing-based transcriptome analysis to identify genes and gene networks putatively involved in various physiological and toxicological processes.

Methodology and Principal Findings

Using Roche 454 pyrosequencing, 857,205 reads containing approximately 340 megabases were obtained from the B. tabaci transcriptome. De novo assembly generated 178,669 unigenes including 30,980 from insects, 17,881 from bacteria, and 129,808 from the nohit. A total of 50,835 (28.45%) unigenes showed similarity to the non-redundant database in GenBank with a cut-off E-value of 10–5. Among them, 40,611 unigenes were assigned to one or more GO terms and 6,917 unigenes were assigned to 288 known pathways. De novo metatranscriptome analysis revealed highly diverse bacterial symbionts in B. tabaci, and demonstrated the host-symbiont cooperation in amino acid production. In-depth transcriptome analysis indentified putative molecular markers, and genes potentially involved in insecticide resistance and nutrient digestion. The utility of this transcriptome was validated by a thiamethoxam resistance study, in which annotated cytochrome P450 genes were significantly overexpressed in the resistant B. tabaci in comparison to its susceptible counterparts.

Conclusions

This transcriptome/metatranscriptome analysis sheds light on the molecular understanding of symbiosis and insecticide resistance in an agriculturally important phloem-feeding insect pest, and lays the foundation for future functional genomics research of the B. tabaci complex. Moreover, current pyrosequencing effort greatly enriched the existing whitefly EST database, and makes RNAseq a viable option for future genomic analysis.

Pyrosequencing-based expression profiling and identification of differentially regulated genes from Manduca sexta, a lepidopteran model insect

Although Manduca sexta has significantly contributed to our knowledge on a variety of insect physiological processes, the lack of its genome sequence hampers the large-scale gene discovery, transcript profiling, and proteomic analysis in this biochemical model species. Here we report our implementation of the RNA-Seq cDNA sequencing approach based on massively parallel pyrosequencing, which allows us to categorize transcripts based on their relative abundances and to discover process- or tissue-specifically regulated genes simultaneously. We obtained 1,821,652 reads with an average length of 289 bp per read from fat body and hemocytes of naïve and microbe-injected M. sexta larvae. After almost all (92.1%) of these reads were assembled into 19,020 contigs, we identified 528 contigs whose relative abundances increased at least 5- and 8-fold in fat body and hemocytes, respectively, after the microbial challenge. Polypeptides encoded by these contigs include pathogen recognition receptors, extracellular and intracellular signal mediators and regulators, antimicrobial peptides, and proteins with no known sequence but likely participating in defense in novel ways. We also found 250 and 161 contigs that were preferentially expressed in fat body and hemocytes, respectively. Furthermore, we integrated data from our previous study and generated a sequence database to support future gene annotation and proteomic analysis in M. sexta. In summary, we have successfully established a combined approach for gene discovery and expression profiling in organisms lacking known genome sequences.

Analyses of Expressed Sequence Tags from Chironomus riparius Using Pyrosequencing : Molecular Ecotoxicology Perspective

OBJECTIVES

Chironomus riparius, a non-biting midge (Chironomidae, Diptera), is extensively used as a model organism in aquatic ecotoxicological studies, and considering the potential of C. riparius larvae as a bio-monitoring species, little is known about its genome sequences. This study reports the results of an Expressed Sequence Tags (ESTs) sequencing project conducted on C. riparius larvae using 454 pyrosequencing.

METHODS

To gain a better understanding of C. riparius transcriptome, we generated ESTs database of C. ripairus using pyrosequencing method.

RESULTS

Sequencing runs, using normalized cDNA collections from fourth instar larvae, yielded 20,020 expressed sequence tags, which were assembled into 8,565 contigs and 11,455 singletons. Sequence analysis was performed by BlastX search against the National Center for Biotechnology Information (NCBI) nucleotide (nr) and uniprot protein database. Based on the gene ontology classifications, 24% (E-value ≤1(-5)) of the sequences had known gene functions, 24% had unknown functions and 52% of sequences did not match any known sequences in the existing database. Sequence comparison revealed 81% of the genes have homologous genes among other insects belonging to the order Diptera providing tools for comparative genome analyses. Targeted searches using these annotations identified genes associated with essential metabolic pathways, signaling pathways, detoxification of toxic metabolites and stress response genes of ecotoxicological interest.

CONCLUSIONS

The results obtained from this study would eventually make ecotoxicogenomics possible in a truly environmentally relevant species, such as, C. riparius.

Applications of next generation sequencing in molecular ecology of non-model organisms

As most biologists are probably aware, technological advances in molecular biology during the last few years have opened up possibilities to rapidly generate large-scale sequencing data from non-model organisms at a reasonable cost. In an era when virtually any study organism can 'go genomic', it is worthwhile to review how this may impact molecular ecology. The first studies to put the next generation sequencing (NGS) to the test in ecologically well-characterized species without previous genome information were published in 2007 and the beginning of 2008. Since then several studies have followed in their footsteps, and a large number are undoubtedly under way. This review focuses on how NGS has been, and can be, applied to ecological, population genetic and conservation genetic studies of non-model species, in which there is no (or very limited) genomic resources. Our aim is to draw attention to the various possibilities that are opening up using the new technologies, but we also highlight some of the pitfalls and drawbacks with these methods. We will try to provide a snapshot of the current state of the art for this rapidly advancing and expanding field of research and give some likely directions for future developments.

A comprehensive transcriptome and immune-gene repertoire of the lepidopteran model host Galleria mellonella

Background

The larvae of the greater wax moth Galleria mellonella are increasingly used (i) as mini-hosts to study pathogenesis and virulence factors of prominent bacterial and fungal human pathogens, (ii) as a whole-animal high throughput infection system for testing pathogen mutant libraries, and (iii) as a reliable host model to evaluate the efficacy of antibiotics against human pathogens. In order to compensate for the lack of genomic information in Galleria, we subjected the transcriptome of different developmental stages and immune-challenged larvae to next generation sequencing.

Results

We performed a Galleria transcriptome characterization on the Roche 454-FLX platform combined with traditional Sanger sequencing to obtain a comprehensive transcriptome. To maximize sequence diversity, we pooled RNA extracted from different developmental stages, larval tissues including hemocytes, and from immune-challenged larvae and normalized the cDNA pool. We generated a total of 789,105 pyrosequencing and 12,032 high-quality Sanger EST sequences which clustered into 18,690 contigs with an average length of 1,132 bases. Approximately 40% of the ESTs were significantly similar (E ≤ e^-03) to proteins of other insects, of which 45% have a reported function. We identified a large number of genes encoding proteins with established functions in immunity related sensing of microbial signatures and signaling, as well as effector molecules such as antimicrobial peptides and inhibitors of microbial proteinases. In addition, we found genes known as mediators of melanization or contributing to stress responses. Using the transcriptomic data, we identified hemolymph peptides and proteins induced upon immune challenge by 2D-gelelectrophoresis combined with mass spectrometric analysis.

Conclusion

Here, we have developed extensive transcriptomic resources for Galleria. The data obtained is rich in gene transcripts related to immunity, expanding remarkably our knowledge about immune and stress-inducible genes in Galleria and providing the complete sequences of genes whose primary structure have only partially been characterized using proteomic methods. The generated data provide for the first time access to the genetic architecture of immunity in this model host, allowing us to elucidate the molecular mechanisms underlying pathogen and parasite response and detailed analyses of both its immune responses against human pathogens, and its coevolution with entomopathogens.

Transcriptomic analysis of the housefly (Musca domestica) larva using massively parallel pyrosequencing

To explore the transcriptome of Musca domestica larvae and to identify unique sequences, we used massively parallel pyrosequencing on the Roche 454-FLX platform to generate a substantial EST dataset of this fly. As a result, we obtained a total of 249,555 ESTs with an average read length of 373 bp. These reads were assembled into 13,206 contigs and 20,556 singletons. Using BlastX searches of the Swissprot and Nr databases, we were able to identify 4,814 contigs and 8,166 singletons as unique sequences. Subsequently, the annotated sequences were subjected to GO analysis and the search results showed a majority of the query sequences were assignable to certain gene ontology terms. In addition, functional classification and pathway assignment were performed by KEGG and 2,164 unique sequences were mapped into 184 KEGG pathways in total. As the first attempt on large-scale RNA sequencing of M. domestica, this general picture of the transcriptome can establish a fundamental resource for further research on functional genomics.

Immunity in lepidopteran insects

Lepidopteran insects provide important model systems for innate immunity of insects, particularly for cell biology of hemocytes and biochemical analyses of plasma proteins. Caterpillars are also among the most serious agricultural pests, and understanding of their immune systems has potential practical significance. An early response to infection in lepidopteran larvae is the activation of hemocyte adhesion, leading to phagocytosis, nodule formation, or encapsulation. Plasmatocytes and granular cells are the hemocyte types involved in these responses. Infectious microorganisms are recognized by binding of hemolymph plasma proteins to microbial surface components. This "pattern recognition" triggers phagocytosis and nodule formation, activation of prophenoloxidase and melanization and the synthesis of antimicrobial proteins that are secreted into the hemolymph. Many hemolymph proteins that function in such innate immune responses of insects were first discovered in lepidopterans. Microbial proteinases and nucleic acids released from lysed host cells may also activate lepidopteran immune responses. Hemolymph antimicrobial peptides and proteins can reach high concentrations and may have activity against a broad spectrum of microorganisms, contributing significantly to clearing of infections. Serine proteinase cascade pathways triggered by microbial components interacting with pattern recognition proteins stimulate activation of the cytokine Spätzle, which initiates the Toll pathway for expression of antimicrobial peptides. A proteinase cascade also results inproteolytic activation of phenoloxidase and production of melanin coatings that trap and kill parasites and pathogens. The proteinases in hemolymph are regulated by specific inhibitors, including members of the serpin superfamily. New developments in lepidopteran functional genomics should lead to much more complete understanding of the immune systems of this insect group.

Recognition of microbial molecular patterns and stimulation of prophenoloxidase activation by a β-1,3-glucanase-related protein in Manduca sexta larval plasma

Detection of pathogenic invaders is the essential first step of a successful defense response in multicellular organisms. In this study, we have identified a new member of the β-1,3-glucanase-related protein superfamily from the tobacco hornworm Manduca sexta. This protein, designated microbe binding protein (MBP), is 61% identical in sequence to Bombyx mori Gram-negative bacteria binding protein, but only 34-36% identical to M. sexta β-1,3-glucan recognition protein-1 and 2. Its mRNA levels were strongly up-regulated in hemocytes and fat body of immune challenged larvae, along with an increase in concentration of the plasma protein. We expressed M. sexta MBP in a baculovirus-insect cell system. The purified protein associated with intact bacteria and fungi. It specifically bound to lipoteichoic acid, lipopolysaccharide, diaminopimelic acid-type peptidoglycans (DAP-PGs) from Escherichia coli and Bacillus subtilis, but less so to laminarin or Lys-type PG from Staphylococcus aureus. The complex binding pattern was influenced by other plasma factors and additional microbial surface molecules. After different amounts of MBP had been incubated with larval plasma on ice, a concentration-dependent increase in phenoloxidase (PO) activity occurred in the absence of any microbial elicitor. The activity increase was also observed in the mixture of plasma and a bacterial or fungal cell wall component. The prophenoloxidase (proPO) activation became more prominent when DAP-PGs, Micrococcus luteus Lys-PG, or lipoteichoic acid was included in the mixture of MBP and plasma. Statistic analysis suggested that a synergistic enhancement of proPO activation was caused by an interaction between MBP and these elicitors, but not S. aureus Lys-PG, lipopolysaccharide, curdlan, or laminarin. These data indicate that M. sexta MBP is a component of the surveillance mechanism and, by working together with other pattern recognition molecules and serine proteinases, triggers the proPO activation system.

Analysis of expressed sequence tags from Maize mosaic rhabdovirus-infected gut tissues of Peregrinus maidis reveals the presence of key components of insect innate immunity

The corn planthopper, Peregrinus maidis, causes direct feeding damage to plants and transmits Maize mosaic rhabdovirus (MMV) in a persistent-propagative manner. MMV must cross several insect tissue layers for successful transmission to occur, and the gut serves as an important barrier for rhabdovirus transmission. In order to facilitate the identification of proteins that may interact with MMV either by facilitating acquisition or responding to virus infection, we generated and analysed the gut transcriptome of P. maidis. From two normalized cDNA libraries, we generated a P. maidis gut transcriptome composed of 20,771 expressed sequence tags (ESTs). Assembly of the sequences yielded 1860 contigs and 14,032 singletons, and biological roles were assigned to 5793 (36%). Comparison of P. maidis ESTs with other insect amino acid sequences revealed that P. maidis shares greatest sequence similarity with another hemipteran, the brown planthopper Nilaparvata lugens. We identified 202 P. maidis transcripts with putative homology to proteins associated with insect innate immunity, including those implicated in the Toll, Imd, JAK/STAT, Jnk and the small-interfering RNA-mediated pathways. Sequence comparisons between our P. maidis gut EST collection and the currently available National Center for Biotechnology Information EST database collection for Ni. lugens revealed that a pathogen recognition receptor in the Imd pathway, peptidoglycan recognition protein-long class (PGRP-LC), is present in these two members of the family Delphacidae; however, these recognition receptors are lacking in the model hemipteran Acyrthosiphon pisum. In addition, we identified sequences in the P. maidis gut transcriptome that share significant amino acid sequence similarities with the rhabdovirus receptor molecule, acetylcholine receptor (AChR), found in other hosts. This EST analysis sheds new light on immune response pathways in hemipteran guts that will be useful for further dissecting innate defence response pathways to rhabdovirus infection.

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.

Pyrosequencing and characterization of immune response genes from the American dog tick, Dermacentor variabilis (L.)

Ticks continue to be a threat to animal and human health, and new and novel control strategies are needed for ticks and tick-borne pathogens. The characterization of the tick-pathogen interface and the tick immune response to microbial infections is fundamental toward the formulation of new control strategies for ticks and the pathogens they transmit. Our overall hypothesis for this research is that the tick immune system manages the maintenance of pathogens. Therefore, discovery of tick immune response genes may provide targets for novel control strategies directed toward reducing vector competency and pathogen transmission. In these studies, 454 pyrosequencing, a high-throughput genomic sequencing method was used to discover tick genes expressed in response to bacterial and fungal infections. Expressed sequence tags (ESTs) were analysed from Dermacentor variabilis ticks that had been injected with bacteria (Escherichia coli, Bacillus subtilis, Micrococcus luteus) or fungi (Saccharomyces cerevisiae and Candida albicans) and ticks that were naturally infected with the intracellular bacterium, Anaplasma marginale. By this approach, ESTs were assembled into 5995 contigs. Contigs fell into the five main functional categories of metabolism, genetic information processing, environmental information processing, cellular processes and human diseases. We identified more than 30 genes that are likely to encode for proteins involved in tick immune function. We further analysed by reverse transcriptase PCR (RT-PCR) the expression of 22 of these genes in each of our bacterial or fungal treatment groups and found that seven were up-regulated. Up-regulation of these seven genes was confirmed for bacterial, but not fungal treatment by quantitative PCR (qPCR). One of these products was novel, encoding a new tick defensin. Our results clearly demonstrate the complexities of the tick immune system and mark new directions for further study and characterization of proteins that modulate microbial infections in the American dog tick.

Functional analysis of four processing products from multiple precursors encoded by a lebocin-related gene from Manduca sexta

Antimicrobial peptides (AMPs) are a crucial component of the natural immune system in insects. Five types of AMPs have been identified in the tobacco hornworm Manduca sexta, including attacin, cecropin, moricin, gloverin, and lebocin. Here we report the isolation of lebocin-related cDNA clones and antibacterial activity of their processed protein products. The 17 cDNA sequences are composed of a constant 5' end and a variable 3' region containing 3-16 copies of an 81-nucleotide repeat. The sequence of the corresponding gene isolated from a M. sexta genomic library and Southern blotting results indicated that the gene lacks introns and exists as a single copy in the genome. The genomic sequence contained 13 complete and one partial copy of the 81-nucleotide repeat. Northern blot analysis revealed multiple transcripts with major size differences. The mRNA level of M. sexta lebocin increased substantially in fat body after larvae had been injected with bacteria. The RXXR motifs in the protein sequences led us to postulate that the precursors are processed by an intracellular convertase to form four bioactive peptides. To test this hypothesis, we chemically synthesized the peptides and examined their antibacterial activity. Peptide 1 killed Gram-positive and Gram-negative bacteria. Peptide 2, similar in sequence to a Galleria mellonella AMP, did not affect the bacterial growth. Peptide 3 was inactive but peptide 3 with an extra Arg at the carboxyl terminus was active against Escherichia coli at a high minimum inhibitory concentration. Peptide 4, encoded by the 81-bp repeat, was inactive in the antibacterial tests. The hypothesis that posttranslational processing of the precursor proteins produces multiple bioactive peptides for defense purposes was validated by identification of peptides 1, 2, and 3 from larval hemolymph via liquid chromatography and tandem mass spectrometry. Comparison with the orthologs from other lepidopteran insects indicates that the same mechanism may be used to generate several functional products from a single precursor.

Uncovering the evolutionary origin of plant molecular processes: comparison of Coleochaete (Coleochaetales) and Spirogyra (Zygnematales) transcriptomes

The large and diverse land plant lineage is nested within a clade of fresh water green algae, the charophytes. Collection of genome-scale data for land plants and other organisms over the past decade has invigorated the field of evolutionary biology. One of the core questions in the field asks: how did a colonization event by a green algae over 450 mya lead to one of the most successful lineages on the tree of life? This question can best be answered using the comparative method, the first step of which is to gather genome-scale data across closely related lineages to land plants. Before sequencing an entire genome it is useful to first gather transcriptome data: it is less expensive, it targets the protein coding regions of the genome, and provides support for gene models for future genome sequencing. We built Expressed Sequence Tag (EST) libraries for two charophyte species, Coleochaete orbicularis (Coleochaetales) and Spirogyra pratensis (Zygnematales). We used both Sanger sequencing and next generation 454 sequencing to cover as much of the transcriptome as possible.

RESULTS

Our sequencing effort for Spirogyra pratensis yielded 9,984 5' Sanger reads plus 598,460 GS FLX Standard 454 sequences; Coleochaete orbicularis yielded 4,992 5' Sanger reads plus 673,811 GS FLX Titanium 454 sequences. After clustering S. pratensis yielded 12,000 unique transcripts, or unigenes, and C. orbicularis yielded 19,000. Both transcriptomes were very plant-like, i.e. most of the transcripts were more similar to streptophytes (land plants + charophyte green algae) than to other green algae in the sister group chlorophytes. BLAST results of several land plant genes hypothesized to be important in early land plant evolution resulted in high quality hits in both transcriptomes revealing putative orthologs ripe for follow-up studies.

CONCLUSIONS

Two main conclusions were drawn from this study. One illustrates the utility of next generation sequencing for transcriptome studies: larger scale data collection at a lower cost enabled us to cover a considerable portion of the transcriptome for both species. And, two, that the charophyte green algal transcriptoms are remarkably plant-like, which gives them the unique capacity to be major players for future evolutionary genomic studies addressing origin of land plant questions.

Pyrosequencing the Manduca sexta larval midgut transcriptome: messages for digestion, detoxification and defence

The tobacco hornworm Manduca sexta is an important model for insect physiology but genomic and transcriptomic data are currently lacking. Following a recent pyrosequencing study generating immune related expressed sequence tags (ESTs), here we use this new technology to define the M. sexta larval midgut transcriptome. We generated over 387,000 midgut ESTs, using a combination of Sanger and 454 sequencing, and classified predicted proteins into those involved in digestion, detoxification and immunity. In many cases the depth of 454 pyrosequencing coverage allowed us to define the entire cDNA sequence of a particular gene. Many new M. sexta genes are described including up to 36 new cytochrome P450s, some of which have been implicated in the metabolism of host plant-derived nicotine. New lepidopteran gene families such as the beta-fructofuranosidases, previously thought to be restricted to Bombyx mori, are also described. An unexpectedly high number of ESTs were involved in immunity, for example 39 contigs encoding serpins, and the increasingly appreciated role of the midgut in insect immunity is discussed. Similar studies of other tissues will allow for a tissue by tissue description of the M. sexta transcriptome and will form an essential complimentary step on the road to genome sequencing and annotation.

Next generation transcriptomes for next generation genomes using est2assembly

Background

The decreasing costs of capillary-based Sanger sequencing and next generation technologies, such as 454 pyrosequencing, have prompted an explosion of transcriptome projects in non-model species, where even shallow sequencing of transcriptomes can now be used to examine a range of research questions. This rapid growth in data has outstripped the ability of researchers working on non-model species to analyze and mine transcriptome data efficiently.

Results

Here we present a semi-automated platform 'est2assembly' that processes raw sequence data from Sanger or 454 sequencing into a hybrid de-novo assembly, annotates it and produces GMOD compatible output, including a SeqFeature database suitable for GBrowse. Users are able to parameterize assembler variables, judge assembly quality and determine the optimal assembly for their specific needs. We used est2assembly to process Drosophila and Bicyclus public Sanger EST data and then compared them to published 454 data as well as eight new insect transcriptome collections.

Conclusions

Analysis of such a wide variety of data allows us to understand how these new technologies can assist EST project design. We determine that assembler parameterization is as essential as standardized methods to judge the output of ESTs projects. Further, even shallow sequencing using 454 produces sufficient data to be of wide use to the community. est2assembly is an important tool to assist manual curation for gene models, an important resource in their own right but especially for species which are due to acquire a genome project using Next Generation Sequencing.

454 pyrosequencing based transcriptome analysis of Zygaena filipendulae with focus on genes involved in biosynthesis of cyanogenic glucosides

Background

An essential driving component in the co-evolution of plants and insects is the ability to produce and handle bioactive compounds. Plants produce bioactive natural products for defense, but some insects detoxify and/or sequester the compounds, opening up for new niches with fewer competitors. To study the molecular mechanism behind the co-adaption in plant-insect interactions, we have investigated the interactions between Lotus corniculatus and Zygaena filipendulae. They both contain cyanogenic glucosides which liberate toxic hydrogen cyanide upon breakdown. Moths belonging to the Zygaena family are the only insects known, able to carry out both de novo biosynthesis and sequestration of the same cyanogenic glucosides as those from their feed plants. The biosynthetic pathway for cyanogenic glucoside biosynthesis in Z. filipendulae proceeds using the same intermediates as in the well known pathway from plants, but none of the enzymes responsible have been identified. A genomics strategy founded on 454 pyrosequencing of the Z. filipendulae transcriptome was undertaken to identify some of these enzymes in Z. filipendulae.

Results

Comparisons of the Z. filipendulae transcriptome with the sequenced genomes of Bombyx mori, Drosophila melanogaster, Tribolium castaneum, Apis mellifera and Anopheles gambiae indicate a high coverage of the Z. filipendulae transcriptome. 11% of the Z. filipendulae transcriptome sequences were assigned to Gene Ontology categories. Candidate genes for enzymes functioning in the biosynthesis of cyanogenic glucosides (cytochrome P450 and family 1 glycosyltransferases) were identified based on sequence length, number of copies and presence/absence of close homologs in D. melanogaster, B. mori and the cyanogenic butterfly Heliconius. Examination of biased codon usage, GC content and selection on gene candidates support the notion of cyanogenesis as an "old" trait within Ditrysia, as well as its origins being convergent between plants and insects.

Conclusion

Pyrosequencing is an attractive approach to gain access to genes in the biosynthesis of bio-active natural products from insects and other organisms, for which the genome sequence is not known. Based on analysis of the Z. filipendulae transcriptome, promising gene candidates for biosynthesis of cyanogenic glucosides was identified, and the suitability of Z. filipendulae as a model system for cyanogenesis in insects is evident.

Proteolytic activation and function of the cytokine Spätzle in the innate immune response of a lepidopteran insect, Manduca sexta

The innate immune response of insects includes induced expression of genes encoding a variety of antimicrobial peptides. The signaling pathways that stimulate this gene expression have been well characterized by genetic analysis in Drosophila melanogaster, but are not well understood in most other insect species. One such pathway involves proteolytic activation of a cytokine called Spätzle, which functions in dorsal-ventral patterning in early embryonic development and in the antimicrobial immune response in larvae and adults. We have investigated the function of Spätzle in a lepidopteran insect, Manduca sexta, in which hemolymph proteinases activated during immune responses have been characterized biochemically. Two cDNA isoforms for M. sexta Spätzle-1 differ because of alternative splicing, resulting in a 10 amino acid residue insertion in the pro-region of proSpätzle-1B that is not present in proSpätzle-1A. The proSpätzle-1A cDNA encodes a 32.7 kDa polypeptide that is 23% and 44% identical to D. melanogaster and Bombyx mori Spätzle-1, respectively. Recombinant proSpätzle-1A was a disulfide-linked homodimer. M. sexta hemolymph proteinase 8 cleaved proSpätzle-1A to release Spätzle-C108, a dimer of the C-terminal 108 residue cystine-knot domain. Injection of Spätzle-C108, but not proSpätzle-1A, into larvae stimulated expression of several antimicrobial peptides and proteins, including attacin-1, cecropin-6, moricin, lysozyme, and the immunoglobulin domain protein hemolin, but did not significantly affect the expression of two bacteria-inducible pattern recognition proteins, immulectin-2 and beta-1,3-glucan recognition protein-2. The results of this and other recent studies support a model for a pathway in which the clip-domain proteinase pro-hemolymph proteinase 6 becomes activated in plasma upon exposure to Gram-negative or Gram-positive bacteria or to beta-1,3-glucan. Hemolymph proteinase 6 then activates pro-hemolymph proteinase 8, which in turn activates Spätzle-1. The resulting Spätzle-C108 dimer is likely to function as a ligand to activate a Toll pathway in M. sexta as a response to a wide variety of microbial challenges, stimulating a broad response to infection. Structured digital abstract * MINT-7295125: Spätzle 1A (uniprotkb:C8BMD1) and Spätzle 1A (uniprotkb:C8BMD1) bind (MI:0407) by comigration in gel electrophoresis (MI:0807).

Extensive conserved synteny of genes between the karyotypes of Manduca sexta and Bombyx mori revealed by BAC-FISH mapping

Background

Genome sequencing projects have been completed for several species representing four highly diverged holometabolous insect orders, Diptera, Hymenoptera, Coleoptera, and Lepidoptera. The striking evolutionary diversity of insects argues a need for efficient methods to apply genome information from such models to genetically uncharacterized species. Constructing conserved synteny maps plays a crucial role in this task. Here, we demonstrate the use of fluorescence in situ hybridization with bacterial artificial chromosome probes as a powerful tool for physical mapping of genes and comparative genome analysis in Lepidoptera, which have numerous and morphologically uniform holokinetic chromosomes.

Methodology/Principal Findings

We isolated 214 clones containing 159 orthologs of well conserved single-copy genes of a sequenced lepidopteran model, the silkworm, Bombyx mori, from a BAC library of a sphingid with an unexplored genome, the tobacco hornworm, Manduca sexta. We then constructed a BAC-FISH karyotype identifying all 28 chromosomes of M. sexta by mapping 124 loci using the corresponding BAC clones. BAC probes from three M. sexta chromosomes also generated clear signals on the corresponding chromosomes of the convolvulus hawk moth, Agrius convolvuli, which belongs to the same subfamily, Sphinginae, as M. sexta.

Conclusions/Significance

Comparison of the M. sexta BAC physical map with the linkage map and genome sequence of B. mori pointed to extensive conserved synteny including conserved gene order in most chromosomes. Only a few rearrangements, including three inversions, three translocations, and two fission/fusion events were estimated to have occurred after the divergence of Bombycidae and Sphingidae. These results add to accumulating evidence for the stability of lepidopteran genomes. Generating signals on A. convolvuli chromosomes using heterologous M. sexta probes demonstrated that BAC-FISH with orthologous sequences can be used for karyotyping a wide range of related and genetically uncharacterized species, significantly extending the ability to develop synteny maps for comparative and functional genomics.

Genomic reaction norms: using integrative biology to understand molecular mechanisms of phenotypic plasticity

Phenotypic plasticity is the development of different phenotypes from a single genotype, depending on the environment. Such plasticity is a pervasive feature of life, is observed for various traits and is often argued to be the result of natural selection. A thorough study of phenotypic plasticity should thus include an ecological and an evolutionary perspective. Recent advances in large-scale gene expression technology make it possible to also study plasticity from a molecular perspective, and the addition of these data will help answer long-standing questions about this widespread phenomenon. In this review, we present examples of integrative studies that illustrate the molecular and cellular mechanisms underlying plastic traits, and show how new techniques will grow in importance in the study of these plastic molecular processes. These techniques include: (i) heterologous hybridization to DNA microarrays; (ii) next generation sequencing technologies applied to transcriptomics; (iii) techniques for studying the function of noncoding small RNAs; and (iv) proteomic tools. We also present recent studies on genetic model systems that uncover how environmental cues triggering different plastic responses are sensed and integrated by the organism. Finally, we describe recent work on changes in gene expression in response to an environmental cue that persist after the cue is removed. Such long-term responses are made possible by epigenetic molecular mechanisms, including DNA methylation. The results of these current studies help us outline future avenues for the study of plasticity.

Applications of new sequencing technologies for transcriptome analysis

Transcriptome analysis has been a key area of biological inquiry for decades. Over the years, research in the field has progressed from candidate gene-based detection of RNAs using Northern blotting to high-throughput expression profiling driven by the advent of microarrays. Next-generation sequencing technologies have revolutionized transcriptomics by providing opportunities for multidimensional examinations of cellular transcriptomes in which high-throughput expression data are obtained at a single-base resolution.

A 454 sequencing approach for large scale phylogenomic analysis of the common emperor scorpion (Pandinus imperator)

In recent years, phylogenetic tree reconstructions that rely on multiple gene alignments that had been deduced from expressed sequence tags (ESTs) have become a popular method in molecular systematics. Here, we present a 454 pyrosequencing approach to infer the transcriptome of the Emperor scorpion Pandinus imperator. We obtained 428,844 high-quality reads (mean length=223+/-50 b) from total cDNA, which were assembled into 8334 contigs (mean length 422+/-313 bp) and 26,147 singletons. About 1200 contigs were successfully annotated by BLAST and orthology search. Specific analyses of eight distinct hemocyanin sequences provided further proof for the quality of the 454 reads and the assembly process. The P. imperator sequences were included in a concatenated alignment of 149 orthologous genes of 67 metazoan taxa that covers 39,842 amino acids. After removal of low-quality regions, 11,168 positions were employed for phylogenetic reconstructions. Using Bayesian and maximum likelihood methods, we obtained strongly supported monophyletic Ecdysozoa, Arthropoda (excluding Tardigrada), Euarthropoda, Pancrustacea and Hexapoda. We also recovered the Myriochelata (Chelicerata+Myriapoda). Within the chelicerates, Pycnogonida form the sister group of Euchelicerata. However, Arachnida were found paraphyletic because the Acari (mites and ticks) were recovered as sister group of a clade comprising Xiphosura, Scorpiones and Araneae. In summary, we have shown that 454 pyrosequencing is a cost-effective method that provides sufficient data and coverage depth for gene detection and multigene-based phylogenetic analyses.

Chromatographic and electrophoretic resolution of proteins and protein complexes from the larval midgut microvilli of Manduca sexta

The microvillar proteome of Manduca sexta larval midguts was analyzed by subjecting brush border membrane vesicles (BBMV) to two different two-dimensional approaches: (i) Anion exchange chromatography followed by SDS-PAGE and (ii) Blue Native-PAGE followed by SDS-PAGE. The first technique was superior to conventional 2-D gel electrophoresis in resolving the most abundant proteins associated with the midgut microvilli. Twenty of them were successfully identified as digestive enzymes, binding targets of the insecticidal Cry1A toxins from Bacillus thuringiensis (Bt), and signal transduction proteins. A homolog of the chlorophyllide A binding protein from the silkworm and several aminopeptidases N represent the most abundant proteins associated with the BBMV. The second technique revealed protein oligomeric complexes associated with midgut microvilli in vivo. Two such complexes contained subunits of the vacuolar ATP synthase complex, and one was an oligomer of the chlorophyllide A binding protein. An additional complex consisted of homo- or hetero-tetramers of three different aminopeptidases N (APNs). As APNs are well-known binding partners of Cry1A toxins, their quaternary structure has implications for Bt toxin mode of action. Both techniques provide a useful complement to conventional 2-D gel electrophoresis in analyzing the complex proteome of the microvillar membrane fraction.

Expressed sequence tags from larval gut of the European corn borer (Ostrinia nubilalis): exploring candidate genes potentially involved in Bacillus thuringiensis toxicity and resistance

Background

Lepidoptera represents more than 160,000 insect species which include some of the most devastating pests of crops, forests, and stored products. However, the genomic information on lepidopteran insects is very limited. Only a few studies have focused on developing expressed sequence tag (EST) libraries from the guts of lepidopteran larvae. Knowledge of the genes that are expressed in the insect gut are crucial for understanding basic physiology of food digestion, their interactions with Bacillus thuringiensis (Bt) toxins, and for discovering new targets for novel toxins for use in pest management. This study analyzed the ESTs generated from the larval gut of the European corn borer (ECB, Ostrinia nubilalis), one of the most destructive pests of corn in North America and the western world. Our goals were to establish an ECB larval gut-specific EST database as a genomic resource for future research and to explore candidate genes potentially involved in insect-Bt interactions and Bt resistance in ECB.

Results

We constructed two cDNA libraries from the guts of the fifth-instar larvae of ECB and sequenced a total of 15,000 ESTs from these libraries. A total of 12,519 ESTs (83.4%) appeared to be high quality with an average length of 656 bp. These ESTs represented 2,895 unique sequences, including 1,738 singletons and 1,157 contigs. Among the unique sequences, 62.7% encoded putative proteins that shared significant sequence similarities (E-value ≤ 10^-3)with the sequences available in GenBank. Our EST analysis revealed 52 candidate genes that potentially have roles in Bt toxicity and resistance. These genes encode 18 trypsin-like proteases, 18 chymotrypsin-like proteases, 13 aminopeptidases, 2 alkaline phosphatases and 1 cadherin-like protein. Comparisons of expression profiles of 41 selected candidate genes between Cry1Ab-susceptible and resistant strains of ECB by RT-PCR showed apparently decreased expressions in 2 trypsin-like and 2 chymotrypsin-like protease genes, and 1 aminopeptidase genes in the resistant strain as compared with the susceptible strain. In contrast, the expression of 3 trypsin- like and 3 chymotrypsin-like protease genes, 2 aminopeptidase genes, and 2 alkaline phosphatase genes were increased in the resistant strain. Such differential expressions of the candidate genes may suggest their involvement in Cry1Ab resistance. Indeed, certain trypsin-like and chymotrypsin-like proteases have previously been found to activate or degrade Bt protoxins and toxins, whereas several aminopeptidases, cadherin-like proteins and alkaline phosphatases have been demonstrated to serve as Bt receptor proteins in other insect species.

Conclusion

We developed a relatively large EST database consisting of 12,519 high-quality sequences from a total of 15,000 cDNAs from the larval gut of ECB. To our knowledge, this database represents the largest gut-specific EST database from a lepidopteran pest. Our work provides a foundation for future research to develop an ECB gut-specific DNA microarray which can be used to analyze the global changes of gene expression in response to Bt protoxins/toxins and the genetic difference(s) between Bt- resistant and susceptible strains. Furthermore, we identified 52 candidate genes that may potentially be involved in Bt toxicity and resistance. Differential expressions of 15 out of the 41 selected candidate genes examined by RT-PCR, including 5 genes with apparently decreased expression and 10 with increased expression in Cry1Ab-resistant strain, may help us conclusively identify the candidate genes involved in Bt resistance and provide us with new insights into the mechanism of Cry1Ab resistance in ECB.

Functions of Manduca sexta hemolymph proteinases HP6 and HP8 in two innate immune pathways

Serine proteinases in insect plasma have been implicated in two types of immune responses; that is, activation of prophenoloxidase (proPO) and activation of cytokine-like proteins. We have identified more than 20 serine proteinases in hemolymph of the tobacco hornworm, Manduca sexta, but functions are known for only a few of them. We report here functions of two additional M. sexta proteinases, hemolymph proteinases 6 and 8 (HP6 and HP8). HP6 and HP8 are each composed of an amino-terminal clip domain and a carboxyl-terminal proteinase domain. HP6 is an apparent ortholog of Drosophila Persephone, whereas HP8 is most similar to Drosophila and Tenebrio spätzle-activating enzymes, all of which activate the Toll pathway. proHP6 and proHP8 are expressed constitutively in fat body and hemocytes and secreted into plasma, where they are activated by proteolytic cleavage in response to infection. To investigate activation and biological activity of HP6 and HP8, we purified recombinant proHP8, proHP6, and mutants of proHP6 in which the catalytic serine was replaced with alanine, and/or the activation site was changed to permit activation by bovine factor Xa. HP6 was found to activate proPO-activating proteinase (proPAP1) in vitro and induce proPO activation in plasma. HP6 was also determined to activate proHP8. Active HP6 or HP8 injected into larvae induced expression of antimicrobial peptides and proteins, including attacin, cecropin, gloverin, moricin, and lysozyme. Our results suggest that proHP6 becomes activated in response to microbial infection and participates in two immune pathways; activation of PAP1, which leads to proPO activation and melanin synthesis, and activation of HP8, which stimulates a Toll-like pathway.