A comprehensive analysis of the Manduca sexta immunotranscriptome

Immune responses of the Asian onion moth, Acrolepiopsis sapporensis, and their genetic factors from RNA-Seq analysis

A nonmodel insect, Acrolepiopsis sapporensis, has been analyzed in immune responses. The total hemocytes in the fifth instar larvae were 2.33 × 106  cells/mL. These hemocytes comprised at least five different types and different relative ratios: 47% granulocytes, 26% plasmatocytes, 11% oenocytoid, 8% prohemocytes, and 5% spherulocytes. Upon bacterial challenge, some of the hemocytes exhibited typical hemocyte-spreading behaviors, such as focal adhesion, and filopodial and lamellipodial cytoplasmic extensions. The hemocyte behaviors induced cellular immune responses demonstrated by nodule formation. In addition, the plasma collected from the immune-challenged larvae exhibited humoral immune responses by bacterial growth inhibition along with enhanced phenoloxidase enzyme activity. These cellular and humoral immune responses were further analyzed by determining the immune-associated genes from a transcriptome generated by RNA-Seq. A total of about 12 Gb sequences led to about 218,116 contigs, which were predicted to encode about 46,808 genes. Comparative expression analysis showed 8392 uniquely expressed genes in the immune-challenged larvae. Differentially expressed gene (DEG) analysis among the commonly expressed genes indicated that 782 genes were upregulated and 548 genes were downregulated in the expressions after bacterial challenge. These immune-associated genes included pattern recognition receptors, immune mediation/signaling genes, and various immune effectors. Specifically, the genetic components of the Toll, IMD, and JAK/STAT immune signaling pathways were included in the DEG database. These results demonstrate the immune responses of A. sapporensis larvae and suggest the genes associated with the immune responses in this nonmodel insect.

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.

Analysis of the Humoral Immunal Response Transcriptome of Ectropis obliqua Infected by Beauveria bassiana

Simple Summary

Ectropis obliqua is a destructive leaf-eating pest that is widely distributed in China’s tea gardens. This pest shows remarkable resistance against multiple insecticides. As an environmentally friendly entomopathogen, Beauveria bassiana has been widely used to prevent agricultural pest infestations. However, the molecular mechanism of B. bassiana against E. obliqua remains unclear. We firstly isolated and identified a highly virulent B. bassiana strain. Using a transcriptome, we analyzed the differences of immune gene expression levels in fat bodies and hemocytes of E. obliqua that were infected by the B. bassiana, which provide molecular insights into the insect–pathogen interaction.

Abstract

Ectropis obliqua is a destructive masticatory pest in China’s tea gardens. Beauveria bassiana as microbial insecticides can effectively control E. obliqua larvae; however, the immune response of this insect infected by B. bassiana are largely unknown. Here, after isolating a highly virulent strain of B. bassiana from E. obliqua, the changes in gene expression among different tissues, including hemocytes and fat bodies, of E. obliqua larvae infected by the entomopathogen were investigated using transcriptome sequencing. A total of 5877 co-expressed differentially expressed genes (DEGs) were identified in hemocytes and fat bodies, of which 5826 were up-regulated in hemocytes and 5784 were up-regulated in fat bodies. We identified 249 immunity-related genes, including pattern recognition receptors, immune effectors, signal modulators, and members of immune pathways. A quantitative real-time PCR analysis confirmed that several pattern recognition receptors were upregulated in hemocytes and fat bodies; however, others were downregulated. The investigated immune effectors (ATT and PPO-1) were suppressed. The results showed that there were tissue differences in the expression of immune genes. This study provides a large number of immunity-related gene sequences from E. obliqua after being infected by B. bassiana, furthering the understanding of the molecular mechanisms of E. obliqua defenses against B. bassiana.

Regulators and signalling in insect antimicrobial innate immunity: Functional molecules and cellular pathways

Insects possess an immune system that protects them from attacks by various pathogenic microorganisms that would otherwise threaten their survival. Immune mechanisms may deal directly with the pathogens by eliminating them from the host organism or disarm them by suppressing the synthesis of toxins and virulence factors that promote the invasion and destructive action of the intruder within the host. Insects have been established as outstanding models for studying immune system regulation because innate immunity can be explored as an integrated system at the level of the whole organism. Innate immunity in insects consists of basal immunity that controls the constitutive synthesis of effector molecules such as antimicrobial peptides, and inducible immunity that is activated after detection of a microbe or its product(s). Activation and coordination of innate immune defenses in insects involve evolutionary conserved immune factors. Previous research in insects has led to the identification and characterization of distinct immune signalling pathways that modulate the response to microbial infections. This work has not only advanced the field of insect immunology, but it has also rekindled interest in the innate immune system of mammals. Here we review the current knowledge on key molecular components of insect immunity and discuss the opportunities they present for confronting infectious diseases in humans.

The innate immune gene Relish and Caudal jointly contribute to the gut immune homeostasis by regulating antimicrobial peptides in Galleria mellonella

Gut microbiota modulates various physiologic processes in insects, such as nutrition, metabolic homeostasis, and pathogen exclusion. Maintaining a normal microbiome is an essential element of the gut homeostasis, requiring an extensive network of regulatory immune responses. The molecular mechanisms driving these various effects and the events leading to the establishment of a normal microbiota in insects are still largely unknown. In this study, the NF-kB (IMD and Toll) signaling pathways in the gut of Galleria mellonella and their roles in the regulation of its gut microbes were assessed. For this, the transcript levels of the IMD pathway (Imd and Relish) and the Toll pathway (Spätzle and Dif/Dorsal) genes were analyzed and the results showed that all the genes were expressed in the gut of G. mellonella. Silencing of Relish resulted in reduced expression levels of the IMD pathway genes and antimicrobial peptides (AMPs) followed by overpopulation of gut bacteria. Antibiotics-treated larvae showed lower expression levels of the IMD and Toll pathway genes followed by lower AMPs expression levels. The expression level of caudal decreased in the antibiotics-treated larvae compared with the controls. Together, these data suggest that the IMD and Toll pathways are active in the gut of G. mellonella. The IMD pathway gene, relish functions in the regulation of gut microbes in this insect model.

Listening to your gut: immune challenge to the gut sensitizes body wall nociception in the caterpillar Manduca sexta

Immune-nociceptor connections are found in animals across phyla. Local inflammation and/or damage results in increased nociceptive sensitivity of the affected area. However, in mammals, immune responses far from peripheral nociceptors, such as immune responses in the gut, produce a general increase in peripheral nociceptive sensitivity. This phenomenon has not, to our knowledge, been found in other animal groups. We found that consuming heat-killed pathogens reduced the tactile force needed to induce a defensive strike in the caterpillar Manduca sexta. This increase in the nociceptive sensitivity of the body wall is probably part of the reconfiguration of behaviour and physiology that occurs during an immune response (e.g. sickness behaviour). This increase may help enhance anti-predator behaviour as molecular resources are shifted towards the immune system. This article is part of the Theo Murphy meeting issue 'Evolution of mechanisms and behaviour important for pain'.

Transcriptome and differentially expressed genes of Busseola fusca (Lepidoptera: Noctuidae) larvae challenged with Cry1Ab toxin

Busseola fusca (Fuller) (Lepidoptera: Noctuidae), a major insect pest of maize in sub-Saharan Africa, has developed high levels of non-recessive resistance to Cry1Ab toxin expressed in genetically modified Bt maize. Multiple resistance mechanisms to various Cry toxins have been identified in Lepidoptera, but no study has yet been done to determine the mechanism of Cry1Ab resistance in B. fusca. Therefore, the larval transcriptome of B. fusca was sequenced, de novo assembled and characterized. Differential expression analysis was performed to compare gene expression profiles of Cry toxin challenged and unchallenged neonate larvae to assess the molecular basis of the defence mechanism employed by this insect. Several genes associated with Cry toxin resistance in other lepidopteran pests were detected in B. fusca. Results suggest that differential expression of metabolic and immune-related genes might explain Cry1Ab toxin defence in this pest (supplemental file). Transcript expression profiles of neonates demonstrated that 33.59% and 60.31% of the 131 differentially expressed genes were upregulated and downregulated in the toxin-challenged neonate larvae, respectively. Transcripts were grouped into two subclusters according to the similarity of their expression patterns. Transcripts in subcluster 1 were moderately upregulated in the toxin-challenged neonate larvae, and, conversely, downregulated in the unchallenged neonate larvae. The solute carrier organic anion transporter, which is involved in insecticide detoxification, was upregulated in the toxin-challenged neonate larvae. Conversely, most of the transcripts in subcluster 2 were moderately downregulated in the toxin-challenged neonate larvae, and upregulated for neonates feeding on non-challenged maize. Four unidentified transcripts were extremely down-regulated in the toxin-challenged neonate larvae, and upregulated in the unchallenged neonate larvae. Further studies are recommended to establish if there is a direct correlation between these differentially expressed genes and the observed resistance. Elucidation of such defence mechanisms is crucial for developing insect resistance management strategies to ensure sustainable use of genetically modified maize in Africa. Nevertheless, this is the first study on gene expression profiles of B. fusca strains challenged with Cry toxin. The transcriptome characterized in this study provides a significant resource base for future studies on B. fusca and contributes to understanding some of the gene regulation and signalling networks involved in the defence of B. fusca against Cry1Ab toxin.

Transcriptomic analysis of interactions between Hyphantria cunea larvae and nucleopolyhedrovirus

BACKGROUND

Hyphantria cunea is a destructive forest pest. To utilise H. cunea nucleopolyhedrovirus (HcNPV) for biological control, understanding insect-virus interactions is essential.

RESULTS

Four cDNA libraries were constructed from H. cunea larvae (two HcNPV-infected and two uninfected groups) and 76 004 expressed genes were obtained by next-generation sequencing. Compared with controls, 272 differentially expressed genes (DEGs) were identified in infected groups, including 162 up-regulated and 110 down-regulated genes. Transcription levels of 20 random DEGs were consistent with the Solexa expression profiles obtained by quantitative real-time PCR. DEGs associated with innate immunity were grouped into several categories, including pattern recognition proteins, heat-shock proteins, UDP-glycosyltransferases, cytochrome P450s, antimicrobial peptides and hormonal signalling proteins. Interestingly, up-regulated host genes included farnesoic acid O-methyltransferase, two juvenile hormone (JH) binding proteins, and a circadian clock-controlled protein related to JH regulation. Pathway enrichment analysis indicates that mitogen-activated protein kinase (MAPK) signalling pathways, key candidate genes and important biological pathways may be associated with molecular modification in H. cunea larvae in response to virus stress.

CONCLUSION

These findings provide insight for future research on the molecular mechanisms of HcNPV invasion and anti-HcNPV mechanisms in H. cunea. A better understanding of gene regulation following HcNPV invasion could help to develop the virus as a bio-insecticide. © 2018 Society of Chemical Industry.

Analysis of the Gynaephora qinghaiensis pupae immune transcriptome in response to parasitization by Thektogaster sp

As a pest on the Qinghai-Tibet Plateau, Gynaephora qinghaiensis causes severe damage to grassland vegetation and its pupae are also natural hosts of Thektogaster sp. To successfully parasitize, endoparasitoids generally introduce or secrete multiple parasitic factors into the host body during the spawning stage to suppress the host immune response. To study the parasitic effects of Thektogaster sp. on G. qinghaiensis, a transcriptome analysis of immune-related genes in parasitized and nonparasitized G. qinghaiensis pupae was performed. A total of 371,260,704 clean reads were assembled into 118,144 unigenes with an average length of 884.33 base pairs. Of these, 23,660 unigenes were annotated in at least one database and 94,484 unigenes were not annotated in any databases. These findings indicated that the majority of the genetic resources (79.97% of all unigenes) in Gynaephora should be further explored. Parasitization significantly affected the transcriptional profile of G. qinghaiensis pupae. The present study identified 12,322 differentially expressed genes and 57 immune-related genes were identified in parasitized G. qinghaiensis pupae. Most immune-related genes were downregulated, potentially resulting from the inhibitory effect of Thektogaster sp. on G. qinghaiensis pupae after parasitization. Overall, the transcriptome analysis sheds valuable light on the molecular mechanisms of G. qinghaiensis parasitization by Thektogaster sp. and promotes the development of novel biocontrol strategies for Gynaephora based on immune defense.

The Manduca sexta serpinome: Analysis of serpin genes and proteins in the tobacco hornworm

Members of the serpin superfamily of proteins occur in animals, plants, bacteria, archaea and some viruses. They adopt a variety of physiological functions, including regulation of immune system, modulation of apoptosis, hormone transport and acting as storage proteins. Most members of the serpin family are inhibitors of serine proteinases. In this study, we searched the genome of Manduca sexta and identified 32 serpin genes. We analyzed the structure of these genes and the sequences of their encoded proteins. Three M. sexta genes (serpin-1, serpin-15, and serpin-28) have mutually exclusive alternatively spliced exons encoding the carboxyl-terminal reactive center loop of the protein, which is the site of interaction with target proteases. We discovered that MsSerpin-1 has 14 splicing isoforms, including two undiscovered in previous studies. Twenty-eight of the 32 M. sexta serpins include a putative secretion signal peptide and are predicted to be extracellular proteins. Phylogenetic analysis of serpins in M. sexta and Bombyx mori indicates that 17 are orthologous pairs, perhaps carrying out essential physiological functions. Analysis of the reactive center loop and hinge regions of the protein sequences indicates that 16 of the serpin genes encode proteins that may lack proteinase inhibitor activity. Our annotation and analysis of these serpin genes and their transcript profiles should lead to future advances in experimental study of their functions in insect biochemistry.

Manduca sexta hemolymph protease-2 (HP2) activated by HP14 generates prophenoloxidase-activating protease-2 (PAP2) in wandering larvae and pupae

Melanization is a universal defense mechanism of insects against microbial infection. During this response, phenoloxidase (PO) is activated from its precursor by prophenoloxidase activating protease (PAP), the terminal enzyme of a serine protease (SP) cascade. In the tobacco hornworm Manduca sexta, hemolymph protease-14 (HP14) is autoactivated from proHP14 to initiate the protease cascade after host proteins recognize invading pathogens. HP14, HP21, proHP1*, HP6, HP8, PAP1-3, and non-catalytic serine protease homologs (SPH1 and SPH2) constitute a portion of the extracellular SP-SPH system to mediate melanization and other immune responses. Here we report the expression, purification, and functional characterization of M. sexta HP2. The HP2 precursor is synthesized in hemocytes, fat body, integument, nerve and trachea. Its mRNA level is low in fat body of 5th instar larvae before wandering stage; abundance of the protein in hemolymph displays a similar pattern. HP2 exists as an active enzyme in plasma of the wandering larvae and pupae in the absence of an infection. HP14 cleaves proHP2 to yield active HP2. After incubating active HP2 with larval hemolymph, we detected higher levels of PO activity, i.e. an enhancement of proPO activation. HP2 cleaved proPAP2 (but not proPAP3 or proPAP1) to yield active PAP2, responsible for a major increase in IEARpNA hydrolysis. PAP2 activates proPOs in the presence of a cofactor of SPH1 and SPH2. In summary, we have identified a new member of the proPO activation system and reconstituted a pathway of HP14-HP2-PAP2-PO. Since high levels of HP2 mRNA were present in integument and active HP2 in plasma of wandering larvae, HP2 likely plays a role in cuticle melanization during pupation and protects host from microbial infection in a soil environment.

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.

Thioester-containing proteins in the tsetse fly (Glossina) and their response to trypanosome infection

Thioester-containing proteins (TEPs) are conserved proteins with a role in innate immune immunity. In the current study, we characterized the TEP family in the genome of six tsetse fly species (Glossina spp.). Tsetse flies are the biological vectors of several African trypanosomes, which cause sleeping sickness in humans or nagana in livestock. The analysis of the tsetse TEP sequences revealed information about their structure, evolutionary relationships and expression profiles under both normal and trypanosome infection conditions. Phylogenetic analysis of the family showed that tsetse flies harbour a genomic expansion of specific TEPs that are not found in other dipterans. We found a general expression of all TEP genes in the alimentary tract, mouthparts and salivary glands. Glossina morsitans and Glossina palpalis TEP genes display a tissue-specific expression pattern with some that are markedly up-regulated when the fly is infected with the trypanosome parasite. A different TEP response was observed to infection with Trypanosoma brucei compared to Trypanosoma congolense, indicating that the tsetse TEP response is trypanosome-specific. These findings are suggestive for the involvement of the TEP family in tsetse innate immunity, with a possible role in the control of the trypanosome parasite.

Immune functions of insect βGRPs and their potential application

Insects rely completely on the innate immune system to sense the foreign bodies and to mount the immune responses. Germ-line encoded pattern recognition receptors play crucial roles in recognizing pathogen-associated molecular patterns. Among them, β-1,3-glucan recognition proteins (βGRPs) and gram-negative bacteria-binding proteins (GNBPs) belong to the same pattern recognition receptor family, which can recognize β-1,3-glucans. Typical insect βGRPs are comprised of a tandem carbohydrate-binding module in the N-terminal and a glucanase-like domain in the C-terminal. The former can recognize triple-helical β-1,3-glucans, whereas the latter, which normally lacks the enzymatic activity, can recruit adapter proteins to initiate the protease cascade. According to studies, insect βGRPs possess at least three types of functions. Firstly, some βGRPs cooperate with peptidoglycan recognition proteins to recognize the lysine-type peptidoglycans upstream of the Toll pathway. Secondly, some directly recognize fungal β-1,3-glucans to activate the Toll pathway and melanization. Thirdly, some form the 'attack complexes' with other immune effectors to promote the antifungal defenses. The current review will focus on the discovery of insect βGRPs, functions of some well-characterized members, structure-function studies and their potential application.

Response analysis of host Spodoptera exigua larvae to infection by Heliothis virescens ascovirus 3h (HvAV-3h) via transcriptome

Heliothis virescens ascovirus 3 h (HvAV-3h), a dsDNA insect virus, belonging to the family Ascoviridae, can infect caterpillars of several Noctuidae species by ovipositing parasitoid wasps. In order to provide a comprehensive overview of the interactive responses of host larvae after infection by the ascovirus, a transcriptome analysis of Spodoptera exigua to HvAV-3h was conducted from 6 to 168 hours post infection (hpi). Approximately 101.64 Gb of RNA sequencing (RNA-seq) data obtained from infected and uninfected S. exigua larvae were used to perform a de novo transcriptome assembly, which generated approximately 62,258 S. exigua unigenes. Using differential gene expression analysis, it was determined that the majority of host transcripts were down-regulated beginning at 6 hpi and continuing throughout the infection period, although there was an increase in up-regulated unigene number during the 12 to 72 hpi stage. It is noteworthy that the most abundantly enriched pathways in KEGG annotation were Metabolism terms, indicating that the host larval metabolic mechanisms were highly influenced post HvAV-3h infection. In addition, the host cuticle protein encoding unigenes were highly down-regulated in most of the situations, suggesting that the host larval cuticle synthesis were inhibited by the viral infection.

Comparative genomics analysis of triatomines reveals common first line and inducible immunity-related genes and the absence of Imd canonical components among hemimetabolous arthropods

Background

Insects operate complex humoral and cellular immune strategies to fend against invading microorganisms. The majority of these have been characterized in Drosophila and other dipterans. Information on hemipterans, including Triatominae vectors of Chagas disease remains incomplete and fractionated.

Results

We identified putative immune-related homologs of three Triatominae vectors of Chagas disease, Triatoma pallidipennis, T. dimidiata and T. infestans (TTTs), using comparative transcriptomics based on established immune response gene references, in conjunction with the predicted proteomes of Rhodnius prolixus, Cimex lecticularis and Acyrthosiphon pisum hemimetabolous. We present a compressive description of the humoral and cellular innate immune components of these TTTs and extend the immune information of other related hemipterans. Key homologs of the constitutive and induced immunity genes were identified in all the studied hemipterans.

Conclusions

Our results in the TTTs extend previous observations in other hemipterans lacking several components of the Imd signaling pathway. Comparison with other hexapods, using published data, revealed that the absence of various Imd canonical components is common in several hemimetabolous species.

Electronic supplementary material

The online version of this article (10.1186/s13071-017-2561-2) contains supplementary material, which is available to authorized users.

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.

Transcriptome profiling of the Plutella xylostella (Lepidoptera: Plutellidae) ovary reveals genes involved in oogenesis

BACKGROUND

As a specialized organ, the insect ovary performs valuable functions by ensuring fecundity and population survival. Oogenesis is the complex physiological process resulting in the production of mature eggs, which are involved in epigenetic programming, germ cell behavior, cell cycle regulation, etc. Identification of the genes involved in ovary development and oogenesis is critical to better understand the reproductive biology and screening for the potential molecular targets in Plutella xylostella, a worldwide destructive pest of economically major crops.

RESULTS

Based on transcriptome sequencing, a total of 7.88Gb clean nucleotides was obtained, with 19,934 genes and 1861 new transcripts being identified. Expression profiling indicated that 61.7% of the genes were expressed (FPKM≥1) in the P. xylostella ovary. GO annotation showed that the pathways of multicellular organism reproduction and multicellular organism reproduction process, as well as gamete generation and chorion were significantly enriched. Processes that were most likely relevant to reproduction included the spliceosome, ubiquitin mediated proteolysis, endocytosis, PI3K-Akt signaling pathway, insulin signaling pathway, cAMP signaling pathway, and focal adhesion were identified in the top 20 'highly represented' KEGG pathways. Functional genes involved in oogenesis were further analyzed and validated by qRT-PCR to show their potential predominant roles in P. xylostella reproduction.

CONCLUSIONS

Our newly developed P. xylostella ovary transcriptome provides an overview of the gene expression profiling in this specialized tissue and the functional gene network closely related to the ovary development and oogenesis. This is the first genome-wide transcriptome dataset of P. xylostella ovary that includes a subset of functionally activated genes. This global approach will be the basis for further studies on molecular mechanisms of P. xylostella reproduction aimed at screening potential molecular targets for integrated pest management.

Dynamics of the Interaction between Cotton Bollworm Helicoverpa armigera and Nucleopolyhedrovirus as Revealed by Integrated Transcriptomic and Proteomic Analyses

Over the past decades, Helicoverpa armigera nucleopolyhedrovirus (HearNPV) has been widely used for biocontrol of cotton bollworm, which is one of the most destructive pest insects in agriculture worldwide. However, the molecular mechanism underlying the interaction between HearNPV and host insects remains poorly understood. In this study, high-throughput RNA-sequencing was integrated with label-free quantitative proteomics analysis to examine the dynamics of gene expression in the fat body of H. armigera larvae in response to challenge with HearNPV. RNA sequencing-based transcriptomic analysis indicated that host gene expression was substantially altered, yielding 3,850 differentially expressed genes (DEGs), whereas no global transcriptional shut-off effects were observed in the fat body. Among the DEGs, 60 immunity-related genes were down-regulated after baculovirus infection, a finding that was consistent with the results of quantitative real-time RT-PCR. Gene ontology and functional classification demonstrated that the majority of down-regulated genes were enriched in gene cohorts involved in energy, carbohydrate, and amino acid metabolic pathways. Proteomics analysis identified differentially expressed proteins in the fat body, among which 76 were up-regulated, whereas 373 were significantly down-regulated upon infection. The down-regulated proteins are involved in metabolic pathways such as energy metabolism, carbohydrate metabolism (CM), and amino acid metabolism, in agreement with the RNA-sequence data. Furthermore, correlation analysis suggested a strong association between the mRNA level and protein abundance in the H. armigera fat body. More importantly, the predicted gene interaction network indicated that a large subset of metabolic networks was significantly negatively regulated by viral infection, including CM-related enzymes such as aldolase, enolase, malate dehydrogenase, and triose-phosphate isomerase. Taken together, transcriptomic data combined with proteomic data elucidated that baculovirus established systemic infection of host larvae and manipulated the host mainly by suppressing the host immune response and down-regulating metabolism to allow viral self-replication and proliferation. Therefore, this study provided important insights into the mechanism of host-baculovirus interaction.

The immune properties of Manduca sexta transferrin

Transferrins are secreted proteins that bind iron. The well-studied transferrins are mammalian serum transferrin, which is involved in iron transport, and mammalian lactoferrin, which functions as an immune protein. Lactoferrin and lactoferrin-derived peptides have bactericidal activity, and the iron-free form of lactoferrin has bacteriostatic activity due to its ability to sequester iron. Insect transferrin is similar in sequence to both serum transferrin and lactoferrin, and its functions are not well-characterized; however, many studies of insect transferrin indicate that it has some type of immune function. The goal of this study was to determine the specific immune functions of transferrin from Manduca sexta (tobacco hornworm). We verified that transferrin expression is upregulated in response to infection in M. sexta larvae and determined that the concentration of transferrin in hemolymph increases from 2 μM to 10 μM following an immune challenge. It is also present in molting fluid and prepupal midgut fluid, two extracellular fluids with immune capabilities. No immune-induced proteolytic cleavage of transferrin in hemolymph was observed; therefore, M. sexta transferrin does not appear to be a source of antimicrobial peptides. Unlike iron-saturated lactoferrin, iron-saturated transferrin had no detectable antibacterial activity. In contrast, 1 μM iron-free transferrin inhibited bacterial growth, and this inhibition was blocked by supplementing the culture medium with 1 μM iron. Our results suggest that M. sexta transferrin does not have bactericidal activity, but that it does have a bacteriostatic function that depends on its iron sequestering ability. This study supports the hypothesis that insect transferrin participates in an iron withholding strategy to protect insects from infectious bacteria.

The parasitic wasp Cotesia congregata uses multiple mechanisms to control host (Manduca sexta) behaviour

Some parasites alter the behaviour of their hosts. The larvae of the parasitic wasp Cotesia congregata develop within the body of the caterpillar Manduca sexta During the initial phase of wasp development, the host's behaviour remains unchanged. However, once the wasps begin to scrape their way out of the caterpillar, the caterpillar host stops feeding and moving spontaneously. We found that the caterpillar also temporarily lost sensation around the exit hole created by each emerging wasp. However, the caterpillars regained responsiveness to nociception in those areas within 1 day. The temporary reduction in skin sensitivity is probably important for wasp survival because it prevents the caterpillar from attacking the emerging wasp larvae with a defensive strike. We also found that expression of plasmatocyte spreading peptide (PSP) and spätzle genes increased in the fat body of the host during wasp emergence. This result supports the hypothesis that the exiting wasps induce a cytokine storm in their host. Injections of PSP suppressed feeding, suggesting that an augmented immune response may play a role in the suppression of host feeding. Injection of wasp larvae culture media into non-parasitized caterpillars reduced feeding, suggesting that substances secreted by the wasp larvae may help alter host behaviour.

Multifaceted biological insights from a draft genome sequence of the tobacco hornworm moth, Manduca sexta

Manduca sexta, known as the tobacco hornworm or Carolina sphinx moth, is a lepidopteran insect that is used extensively as a model system for research in insect biochemistry, physiology, neurobiology, development, and immunity. One important benefit of this species as an experimental model is its extremely large size, reaching more than 10 g in the larval stage. M. sexta larvae feed on solanaceous plants and thus must tolerate a substantial challenge from plant allelochemicals, including nicotine. We report the sequence and annotation of the M. sexta genome, and a survey of gene expression in various tissues and developmental stages. The Msex_1.0 genome assembly resulted in a total genome size of 419.4 Mbp. Repetitive sequences accounted for 25.8% of the assembled genome. The official gene set is comprised of 15,451 protein-coding genes, of which 2498 were manually curated. Extensive RNA-seq data from many tissues and developmental stages were used to improve gene models and for insights into gene expression patterns. Genome wide synteny analysis indicated a high level of macrosynteny in the Lepidoptera. Annotation and analyses were carried out for gene families involved in a wide spectrum of biological processes, including apoptosis, vacuole sorting, growth and development, structures of exoskeleton, egg shells, and muscle, vision, chemosensation, ion channels, signal transduction, neuropeptide signaling, neurotransmitter synthesis and transport, nicotine tolerance, lipid metabolism, and immunity. This genome sequence, annotation, and analysis provide an important new resource from a well-studied model insect species and will facilitate further biochemical and mechanistic experimental studies of many biological systems in insects.

Characterization and regulation of expression of an antifungal peptide from hemolymph of an insect, Manduca sexta

Insects secrete antimicrobial peptides as part of the innate immune response. Most antimicrobial peptides from insects have antibacterial but not antifungal activity. We have characterized an antifungal peptide, diapausin-1 from hemolymph of a lepidopteran insect, Manduca sexta (tobacco hornworm). Diapausin-1 was isolated by size exclusion chromatography from hemolymph plasma of larvae that were previously injected with a yeast, Saccharomyces cerevisiae. Fractions containing activity against S. cerevisiae were analyzed by SDS-PAGE and MALDI-TOF MS/MS and found to contain a 45-residue peptide that was encoded by sequences identified in M. sexta transcriptome and genome databases. A cDNA for diapausin-1 was cloned from cDNA prepared from fat body RNA. Diapausin-1 is a member of the diapausin family of peptides, which includes members known to have antifungal activity. The M. sexta genome contains 14 genes with high similarity to diapausin-1, each with 6 conserved Cys residues. Diapausin-1 was produced as a recombinant protein in Escherichia coli. Purified recombinant diapausin-1 was active against S. cerevisiae, with IC50 of 12 μM, but had no detectable activity against bacteria. Spores of some plant fungal pathogens treated with diapausin-1 had curled germination tubes or reduced and branched hyphal growth. Diapausin-1 mRNA level in fat body strongly increased after larvae were injected with yeast or with Micrococcus luteus. In addition, diapausin-1 mRNA levels increased in midgut and fat body at the wandering larval stage prior to pupation, suggesting developmental regulation of the gene. Our results indicate that synthesis of diapausin-1 is part of an antifungal innate immune response to infection in M. sexta.

Immune memory in invertebrates

Evidence for innate immune memory (or 'priming') in invertebrates has been accumulating over the last years. We here provide an in-depth review of the current state of evidence for immune memory in invertebrates, and in particular take a phylogenetic viewpoint. Invertebrates are a very heterogeneous group of animals and accordingly, evidence for the phenomenon of immune memory as well as the hypothesized molecular underpinnings differ largely for the diverse invertebrate taxa. The majority of research currently focuses on Arthropods, while evidence from many other groups of invertebrates is fragmentary or even lacking. We here concentrate on immune memory that is induced by pathogenic challenges, but also extent our view to a non-pathogenic context, i.e. allograft rejection, which can also show forms of memory and can inform us about general principles of specific self-nonself recognition. We discuss definitions of immune memory and a number of relevant aspects such as the type of antigens used, the route of exposure, and the kinetics of reactions following priming.

Characterization of a dual-CRD galectin in the silkworm Bombyx mori

Galectins (S-type lectins) are an ancient family of lectins with the β-galactoside binding activity. In mammals, galectins play essential roles in many biological processes, such as development, immune homeostasis and tumor progression. However, few studies have been devoted to their functions in insects. Here, we characterized the only dual-CRD galectin in the silkworm Bombyx mori (BmGalectin-4). BmGalectin-4 cDNA possesses an open reading frame of 1089 bp, which encodes a putative galectin of 363 amino acids containing tandem carbohydrate recognition domains (CRDs). BmGalectin-4 was expressed in various tissues but the protein was most abundant in fertilized eggs. Its transcript level in fertilized eggs was upregulated upon bacterial challenge. Recombinant BmGalectin-4 purified from Escherichia coli bound to bacterial cell wall components and bacterial cells. In addition, the recombinant protein induced bacterial agglutination, but did not have antibacterial activity against selected microorganisms. Taken together, our results suggest that BmGalectin-4 may function as a pattern recognition receptor primarily in silkworm fertilized eggs.

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.

Transcriptomic insight into the immune defenses in the ghost moth, Hepialus xiaojinensis, during an Ophiocordyceps sinensis fungal infection

Hepialus xiaojinensis is an economically important species of Lepidopteran insect. The fungus Ophiocordyceps sinensis can infect its larvae, which leads to mummification after 5-12 months, providing a valuable system with which to study interactions between the insect hosts and pathogenic fungi. However, little sequence information is available for this insect. A time-course analysis of the fat body transcriptome was performed to explore the host immune response to O. sinensis infection. In total, 50,164 unigenes were obtained by assembling the reads from two high-throughput approaches: 454 pyrosequencing and Illumina Hiseq2000. Hierarchical clustering and functional examination revealed four major gene clusters. Clusters 1-3 included transcripts markedly induced by the fungal infection within 72 h. Cluster 4, with a lower number of transcripts, was suppressed during the early phase of infection but returned to normal expression levels sometime before 1 year. Based on sequence similarity to orthologs known to participate in immune defenses, 258 candidate immunity-related transcripts were identified, and their functions were hypothesized. The genes were more primitive than those in other Lepidopteran insects. In addition, lineage-specific family expansion of the clip-domain serine proteases and C-type lectins were apparent and likely caused by selection pressures. Global expression profiles of immunity-related genes indicated that H. xiaojinensis was capable of a rapid response to an O. sinensis challenge; however, the larvae developed tolerance to the fungus after prolonged infection, probably due to immune suppression. Specifically, antimicrobial peptide mRNAs could not be detected after chronic infection, because key components of the Toll pathway (MyD88, Pelle and Cactus) were downregulated. Taken together, this study provides insights into the defense system of H. xiaojinensis, and a basis for understanding the molecular aspects of the interaction between the host and the entomopathogen.

The plastic response of Manduca sexta to host and non-host plants

Specialist insect herbivores have evolved efficient ways to adapt to the major defenses of their host plants. Although Manduca sexta, specialized on Solanaceous plants, has become a model organism for insect molecular biology, little is known about its adaptive responses to the chemical defenses of its hosts. To study larval performance and transcriptomic responses to host and non-host plants, we conducted developmental assays and replicated RNAseq experiments with Manduca larvae fed on different Solanaceous plants as well as on a Brassicaceous non-host plant, Brassica napus. Manduca larvae developed fastest on Nicotiana attenuata, but no significant differences in performance were found on larvae fed on other Solanaceae or the non-host B. napus. The RNAseq experiments revealed that Manduca larvae display plastic responses at the gene expression level, and transcriptional signatures specific to the challenges of each host- and non-host plant. Our observations are not consistent with expectations that specialist herbivores would perform poorly on non-host plants. Instead, our findings demonstrate the ability of this specialized insect herbivore to efficiently use a larger repertoire of host plants than it utilizes in the field.

Phylogenetic analysis and expression profiling of the pattern recognition receptors: Insights into molecular recognition of invading pathogens in Manduca sexta

Pattern recognition receptors (PRRs) detect microbial pathogens and trigger innate immune responses. Previous biochemical studies have elucidated the physiological functions of eleven PRRs in Manduca sexta but our understanding of the recognition process is still limited, lacking genomic perspectives. While 34 C-type lectin-domain proteins and 16 Toll-like receptors are reported in the companion papers, we present here 120 other putative PRRs identified through the genome annotation. These include 76 leucine-rich repeat (LRR) proteins, 14 peptidoglycan recognition proteins, 6 EGF/Nim-domain proteins, 5 β-1,3-glucanase-related proteins, 4 galectins, 4 fibrinogen-related proteins, 3 thioester proteins, 5 immunoglobulin-domain proteins, 2 hemocytins, and 1 Reeler. Sequence alignment and phylogenetic analysis reveal the evolution history of a diverse repertoire of proteins for pathogen recognition. While functions of insect LRR proteins are mostly unknown, their structure diversification is phenomenal: In addition to the Toll homologs, 22 LRR proteins with a signal peptide are expected to be secreted; 18 LRR proteins lacking signal peptides may be cytoplasmic; 36 LRRs with a signal peptide and a transmembrane segment may be non-Toll receptors on the surface of cells. Expression profiles of the 120 genes in 52 tissue samples reflect complex regulation in various developmental stages and physiological states, including some likely by Rel family transcription factors via κB motifs in the promoter regions. This collection of information is expected to facilitate future biochemical studies detailing their respective roles in this model insect.

The immune signaling pathways of Manduca sexta

Signal transduction pathways and their coordination are critically important for proper functioning of animal immune systems. Our knowledge of the constituents of the intracellular signaling network in insects mainly comes from genetic analyses in Drosophila melanogaster. To facilitate future studies of similar systems in the tobacco hornworm and other lepidopteran insects, we have identified and examined the homologous genes in the genome of Manduca sexta. Based on 1:1 orthologous relationships in most cases, we hypothesize that the Toll, Imd, MAPK-JNK-p38 and JAK-STAT pathways are intact and operative in this species, as are most of the regulatory mechanisms. Similarly, cellular processes such as autophagy, apoptosis and RNA interference probably function in similar ways, because their mediators and modulators are mostly conserved in this lepidopteran species. We have annotated a total of 186 genes encoding 199 proteins, studied their domain structures and evolution, and examined their mRNA levels in tissues at different life stages. Such information provides a genomic perspective of the intricate signaling system in a non-drosophiline insect.

Transcriptomic response of Manduca sexta immune tissues to parasitization by the bracovirus associated wasp Cotesia congregata

During oviposition, Cotesia congregata parasitoid wasps inject into their host, Manduca sexta, some biological factors such as venom, ovarian fluid and a symbiotic polydnavirus (PDV) named Cotesia congregata bracovirus (CcBV). During parasitism, complex interactions occur between wasp-derived factors and host targets that lead to important modifications in host physiology. In particular, the immune response leading to wasp egg encapsulation is inhibited allowing wasp survival. To date, the regulation of host genes during the interaction had only been studied for a limited number of genes. In this study, we analysed the global impact of parasitism on host gene regulation 24 h post oviposition by high throughput 454 transcriptomic analyses of two tissues known to be involved in the host immune response (hemocytes and fat body). To identify specific effects of parasitism on host transcription at this time point, transcriptomes were obtained from non-treated and parasitized larvae, and also from larvae injected with heat-killed bacteria and double stimulated larvae that were parasitized prior to bacterial challenge. Results showed that, immune challenge by bacteria leads to induction of certain antimicrobial peptide (AMP) genes in M. sexta larvae whether they were parasitized or not prior to bacterial challenge. These results show that at 24 h post oviposition pathways leading to expression of AMP genes are not all inactivated suggesting wasps are in an antiseptic environment. In contrast, at this time point genes involved in phenoloxidase activation and cellular immune responses were globally down-regulated after parasitism in accordance with the observed inhibition of wasp egg encapsulation.

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.

Structural features, evolutionary relationships, and transcriptional regulation of C-type lectin-domain proteins in Manduca sexta

C-type lectins (CTLs) are a large family of Ca(2+)-dependent carbohydrate-binding proteins recognizing various glycoconjugates and functioning primarily in immunity and cell adhesion. We have identified 34 CTLDP (for CTL-domain protein) genes in the Manduca sexta genome, which encode proteins with one to three CTL domains. CTL-S1 through S9 (S for simple) have one or three CTL domains; immulectin-1 through 19 have two CTL domains; CTL-X1 through X6 (X for complex) have one or two CTL domains along with other structural modules. Nine simple CTLs and seventeen immulectins have a signal peptide and are likely extracellular. Five complex CTLs have both an N-terminal signal peptide and a C-terminal transmembrane region, indicating that they are membrane anchored. Immulectins exist broadly in Lepidoptera and lineage-specific gene duplications have generated three clusters of fourteen genes in the M. sexta genome, thirteen of which have similar expression patterns. In contrast to the family expansion, CTL-S1∼S6, S8, and X1∼X6 have 1:1 orthologs in at least four lepidopteran/dipteran/coleopteran species, suggestive of conserved functions in a wide range of holometabolous insects. Structural modeling suggests the key residues for Ca(2+)-dependent or independent binding of certain carbohydrates by CTL domains. Promoter analysis identified putative κB motifs in eighteen of the CTL genes, which did not have a strong correlation with immune inducibility in the mRNA or protein levels. Together, the gene identification, sequence comparisons, structure modeling, phylogenetic analysis, and expression profiling establish a solid foundation for future studies of M. sexta CTL-domain proteins.

A genome-wide analysis of antimicrobial effector genes and their transcription patterns in Manduca sexta

Antimicrobial proteins/peptides (AMPs) are effectors of innate immune systems against pathogen infection in multicellular organisms. Over half of the AMPs reported so far come from insects, and these effectors act in concert to suppress or kill bacteria, fungi, viruses, and parasites. In this work, we have identified 86 AMP genes in the Manduca sexta genome, most of which seem likely to be functional. They encode 15 cecropins, 6 moricins, 6 defensins, 3 gallerimycins, 4 X-tox splicing variants, 14 diapausins, 15 whey acidic protein homologs, 11 attacins, 1 gloverin, 4 lebocins, 6 lysozyme-related proteins, and 4 transferrins. Some of these genes (e.g. attacins, cecropins) constitute large clusters, likely arising after rounds of gene duplication. We compared the amino acid sequences of M. sexta AMPs with their homologs in other insects to reveal conserved structural features and phylogenetic relationships. Expression data showed that many of them are synthesized in fat body and midgut during the larval-pupal molt. Certain genes contain one or more predicted κB binding sites and other regulatory elements in their promoter regions, which may account for the dramatic mRNA level increases in fat body and hemocytes after an immune challenge. Consistent with these strong mRNA increases, many AMPs become highly abundant in the larval plasma at 24 h after the challenge, as demonstrated in our previous peptidomic study. Taken together, these data suggest the existence of a large repertoire of AMPs in M. sexta, whose expression is up-regulated via immune signaling pathways to fight off invading pathogens in a coordinated manner.

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.

Evolutionary genetics of insect innate immunity

Patterns of evolution in immune defense genes help to understand the evolutionary dynamics between hosts and pathogens. Multiple insect genomes have been sequenced, with many of them having annotated immune genes, which paves the way for a comparative genomic analysis of insect immunity. In this review, I summarize the current state of comparative and evolutionary genomics of insect innate immune defense. The focus is on the conserved and divergent components of immunity with an emphasis on gene family evolution and evolution at the sequence level; both population genetics and molecular evolution frameworks are considered.

Changes in the hemolymph protein profiles in Galleria mellonella infected with Bacillus thuringiensis involve apolipophorin III. The effect of heat shock

This report concerns the effect of heat shock on host-pathogen interaction in Galleria mellonella infected with Bacillus thuringiensis. We show enhanced activity against Gram-positive bacteria in the hemolymph of larvae pre-exposed to heat shock before infection with B. thuringiensis. Heat shock influenced the protein pattern in the hemolymph of infected larvae: more peptides with a molecular weight below 10 kDa were detected in comparison with nonshocked animals. Additionally, we noticed that the amount of apolipophorin III (apoLp-III) in the hemolymph decreased transiently following infection, which was considerably higher in larvae pre-exposed to heat shock. On the other hand, its expression in the fat body showed a consequent infection-induced decline, observed equally in shocked and nonshocked animals. This suggests that the amount of apoLp-III in the hemolymph of G. mellonella larvae is regulated at multiple levels. We also report that this protein is more resistant to degradation in the hemolymph of larvae pre-exposed to heat shock in comparison to nonshocked larvae. Two-dimensional analysis revealed the presence of three isoforms of apoLp-III, all susceptible to proteolytic degradation. However, one of them was the most abundant, both in the protease-treated and untreated hemolymph. Taking into consideration that, in general, apoLp-III has a stimulative effect on different immune-related hemolymph proteins and peptides, the reported findings bring us closer to understanding the effect of heat shock on the resistance of G. mellonella to infection.

Dscam and pancrustacean immune memory - a review of the evidence

Evidence is accumulating for a memory-like phenomenon in the immune defence of invertebrates. Down syndrome cell adhesion molecule (Dscam) has been proposed as a key candidate for a somatically diversified receptor system in the crustaceans and insects (Pancrustacea) that could enable challenge-specific protection. However, what is the evidence for an involvement of Dscam in pancrustacean immune memory, and in particular specificity? Here we review the current state of the art, and discuss hypotheses of how Dscam could be involved in immunity. We conclude that while there is increasing evidence for the involvement of Dscam in pancrustacean immunity, crucial experiments to address whether it plays a role in specificity upon secondary encounter with a pathogen still remain to be done.

The structure and function of thioester-containing proteins in arthropods

Thioester-containing proteins (TEPs) form an ancient and diverse family of secreted proteins that play central roles in the innate immune response. Two families of TEPs, complement factors and α₂-macroglobulins, have been known and studied in vertebrates for many years, but only in the last decade have crystal structures become available. In the same period, the presence of two additional classes of TEPs has been revealed in arthropods. In this review, we discuss the common structural features TEPs and how this knowledge can be applied to the many arthropod TEPs of unknown function. TEPs perform a wide variety of functions that are driven by different quaternary structures and protein-protein interactions between a common set of folded domains. A common theme is regulated conformational change triggered by proteolysis. Structure-function analysis of the diverse arthropod TEPs may identify not just new mechanisms in innate immunity but also interfaces between immunity, development and cell death.

Identification of conserved and novel microRNAs in Manduca sexta and their possible roles in the expression regulation of immunity-related genes

The tobacco hornworm Manduca sexta has served as a model for insect biochemical and physiological research for decades. However, knowledge of the posttranscriptional regulation of gene expression by microRNAs is still rudimentary in this species. Our previous study (Zhang et al., 2012) identified 163 conserved and 13 novel microRNAs in M. sexta, most of which were present at low levels in pupae. To identify additional M. sexta microRNAs and more importantly to examine their possible roles in the expression regulation of immunity-related genes, we constructed four small RNA libraries using fat body and hemocytes from naïve or bacteria-injected larvae and obtained 32.9 million reads of 18-31 nucleotides by Illumina sequencing. Mse-miR-929 and mse-miR-1b (antisense microRNA of mse-miR-1) were predicted in the previous study and now found to be conserved microRNAs in the tissue samples. We also found four novel microRNAs, two of which result from a gene cluster. Mse-miR-281-star, mse-miR-965-star, mse-miR-31-star, and mse-miR-9b-star were present at higher levels than their respective mature strands. Abundance changes of microRNAs were observed after the immune challenge. Based on the quantitative data of mRNA levels in control and induced fat body and hemocytes as well as the results of microRNA target site prediction, we suggest that certain microRNAs and microRNA*s regulate gene expression for pattern recognition, prophenoloxidase activation, cellular responses, antimicrobial peptide synthesis, and conserved intracellular signal transduction (Toll, IMD, JAK-STAT, MAPK-JNK-p38, and small interfering RNA pathways). In summary, this study has enriched our knowledge on M. sexta microRNAs and how some of them may participate in the expression regulation of immunity-related genes.

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.

Transcriptome immune analysis of the invasive beetle Octodonta nipae (Maulik) (Coleoptera: Chrysomelidae) parasitized by Tetrastichus brontispae Ferrière (Hymenoptera: Eulophidae)

The beetle Octodonta nipae (Maulik) (Coleoptera: Chrysomelidae) is a serious invasive insect pest of palm plants in southern China, and the endoparasitoid Tetrastichus brontispae Ferrière (Hymenoptera: Eulophidae) is a natural enemy of this pest that exhibits great ability in the biocontrol of O. nipae. For successful parasitism, endoparasitoids often introduce or secrete various virulence factors to suppress host immunity. To investigate the effects of parasitization by T. brontispae on the O. nipae immune system, the transcriptome of O. nipae pupae was analyzed with a focus on immune-related genes through Illumina sequencing. De novo assembly generated 49,919 unigenes with a mean length of 598 bp. Of these genes, 27,490 unigenes (55.1% of all unigenes) exhibited clear homology to known genes in the NCBI nr database. Parasitization had significant effects on the transcriptome profile of O. nipae pupae, and most of these differentially expressed genes were down-regulated. Importantly, the expression profiles of immune-related genes were significantly regulated after parasitization. Taken together, these transcriptome sequencing efforts shed valuable light on the host (O. nipae) manipulation mechanisms induced by T. brontispae, which will pave the way for the development of novel immune defense-based management strategies of O. nipae, and provide a springboard for further molecular analyses, particularly of O. nipae invasion.

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.