Deep sequencing-based transcriptome analysis of Plutella xylostella larvae parasitized by Diadegma semiclausum

Determination of suitable reference genes for RT-qPCR normalisation in Bombyx mori (Lepidoptera: Bombycidae) infected by the parasitoid Exorista sorbillans (Diptera, Tachinidae)

The silkworm Bombyx mori (Lepidoptera: Bombycidae) is a lepidopteran model insect of great economic importance. The parasitoid Exorista sorbillans (Diptera, Tachinidae) is the major pest of B. mori and also a promising candidate for biological control. However, the molecular interactions between hosts and dipteran parasitoids have only partially been studied. Gene expression analysis by reverse-transcription quantitative real-time polymerase chain reaction (RT-qPCR) is indispensable to characterise their interactions. Accurate normalisation of RT-qPCR-based gene expression requires the use of reference genes that are constantly expressed irrespective of experimental conditions. In this study, the expression stability of 13 traditionally used reference genes was estimated by five statistical algorithms (ΔCt, geNorm, Normfinder, BestKeeper, and RefFinder) to determine the best reference genes for gene expression studies in different tissues of B. mori under E. sorbillans parasitism. Specifically, TATA-box-binding protein was the best reference gene in epidermis and testis, while elongation factor 1α was the most stable gene in prothoracic gland and midgut. Elongation factor 1γ, ribosomal protein L3, actin A1, ribosomal protein L40, glyceraldehyde-3-phosphate dehydrogenase and eukaryotic translation initiation factor 4A were the most suitable genes in head, silk gland, fat body, haemolymph, Malpighian tubule and ovary, respectively. Our study offers a set of suitable reference genes for gene expression normalisation in B. mori under the parasitic stress of E. sorbillans, which will benefit the in-depth exploration of host-dipteran parasitoid interactions, and also provide insights for further improvements of B. mori resistance against parasitoids and biocontrol efficacy of dipteran parasitoids.

Transcriptional regulation of host insulin signaling pathway genes controlling larval development by Microplitis manilae parasitization

Idiobiont parasitoids using other insects as hosts sabotage the host growth and development to ensure their offspring survival. Numerous studies have discovered that insect development is subtly regulated by the conserved insulin signaling pathway. However, little is known about how wasp parasitization disrupts host development controlled by the insulin signaling pathway. Here we address this study to determine the effect of wasp parasitism on host Spodoptera frugiperda development using the idiobiont parasitoid Microplitis manilae as a model. Upon M. manilae parasitization, the body weight, body length, and food consumption of host insect were dramatically reduced compared to the unparasitized S. frugiperda. We next identified the core genes involved in host insulin signaling pathway and further analyzed the domain organizations of these genes. Phylogenetic reconstruction based on the insulin receptors clustered S. frugiperda together with other noctuidae insects. In the latter study, we profiled the expression patterns of host insulin signaling pathway genes in response to M. manilae parasitization at 2, 24, and 48 h, significant decreases in mRNA levels were recorded in S. frugiperda larvae upon 24 and 48 h parasitization. These current findings substantially add to our understanding of the physiological interaction between parasitoid and host insects, thus contributing to revealing the molecular mechanism of parasitic wasps regulating host development.

Sex dependent transcriptome responses of the diamondback moth, Plutella xylostella L. to cold stress

Temperature has fundamental influences on the performance and distribution of insects. While considerable attention has been devoted to extreme conditions, particularly extreme cold conditions, few studies have investigated effects of mild cold conditions on insects. We examined the transcriptomic changes in mid-fourth instar larvae of both sexes reared at 10 °C and 25 °C to investigate sex-dependent responses of Plutella xylostella to mild cold stress. There were 624 differentially expressed genes (DEGs) in females, the majority of which (n = 386) were down-regulated. In males 3239 genes were differentially expressed and the majority (n = 2341) were up-regulated. Only 280 DEGs were common to both sexes. In females, there were no DEGs encoding heat shock or cold shock proteins, but six of these DEGs were found in males. These differences suggest that females and males might adopt some different strategies to cope with cold stress and/or that they were affected by rearing under cold conditions to different degrees and in different ways. In addition, DEGs encoding antimicrobial peptides, cytochrome P450 monooxygenases, fatty acid-related enzymes, cuticle proteins, myofilament, and hormone-related proteins were found in both sexes under cold stress. The transcriptome study reveals unexpected sex-dependent thermal responses and provides new information of how an insect that does not diapause copes with low temperatures.

Genome, host genome integration, and gene expression in Diadegma fenestrale ichnovirus from the perspective of coevolutionary hosts

Polydnaviruses (PDVs) exhibit species-specific mutualistic relationships with endoparasitoid wasps. PDVs can be categorized into bracoviruses and ichnoviruses, which have independent evolutionary origins. In our previous study, we identified an ichnovirus of the endoparasitoid Diadegma fenestrale and named it DfIV. Here, DfIV virions from the ovarian calyx of gravid female wasps were characterized. DfIV virion particles were ellipsoidal (246.5 nm × 109.0 nm) with a double-layered envelope. Next-generation sequencing of the DfIV genome revealed 62 non-overlapping circular DNA segments (A1-A5, B1-B9, C1-C15, D1-D23, E1-E7, and F1-F3); the aggregate genome size was approximately 240 kb, and the GC content (43%) was similar to that of other IVs (41%-43%). A total of 123 open reading frames were predicted and included typical IV gene families such as repeat element protein (41 members), cysteine motif (10 members), vankyrin (9 members), polar residue-rich protein (7 members), vinnexin (6 members), and N gene (3 members). Neuromodulin N (2 members) was found to be unique to DfIV, along with 45 hypothetical genes. Among the 62 segments, 54 showed high (76%-98%) sequence similarities to the genome of Diadegma semiclausum ichnovirus (DsIV). Three segments, namely, D22, E3, and F2, contained lepidopteran host genome integration motifs with homologous regions of about 36-46 bp between them (Diadegma fenestrale ichnovirus, DfIV and lepidopteran host, Plutella xylostella). Most of the DfIV genes were expressed in the hymenopteran host and some in the lepidopteran host (P. xylostella), parasitized by D. fenestrale. Five segments (A4, C3, C15, D5, and E4) were differentially expressed at different developmental stages of the parasitized P. xylostella, and two segments (C15 and D14) were highly expressed in the ovaries of D. fenestrale. Comparative analysis between DfIV and DsIV revealed that the genomes differed in the number of segments, composition of sequences, and internal sequence homologies.

Identification of Key Genes Involved in Resistance to Early Stage of BmNPV Infection in Silkworms

Bombyx mori nucleopolyhedrovirus (BmNPV) is one of the most serious pathogens restricting the sustainable development of the sericulture industry. Currently, there is no effective treatment for BmNPV infection in silkworms, and the mechanism underlying BmNPV resistance in silkworms is also not clear. In this study, comparative transcriptome analyses were carried out in midguts of two silkworm varieties, namely BaiyuN, which is a resistance variety, and Baiyu, which is a susceptible variety, at five different time points (i.e., 0, 1, 3, 6, and 9 h) post-BmNPV infection to detect the early-stage transcriptional changes in these silkworms. In total, 1911 and 1577 differentially expressed genes (DEGs) were identified in the Baiyu and BaiyuN varieties, respectively, involving a total of 48 metabolic pathways. Of these pathways, eight were shared by the Baiyu and BaiyuN varieties in response to BmNPV infection. Notably, four genes (i.e., BGIBMGA08815, BGIBMGA003935, BGIBMGA003571, BGIBMGA010059) were upregulated in the Baiyu variety while downregulated in the BaiyuN variety. The inhibited expression of these four genes in the resistant variety highlighted their potential roles in the resistance of early-stage viral replication. Thus, our study provided a new avenue for the further study of the mechanism underlying BmNPV infection in silkworms and the potential treatment of BmNPV infection.

Transcriptome analysis of megalurothrips usitatus (Bagnall) identifies olfactory genes with ligands binding characteristics of MusiOBP1 and MusiCSP1

The olfactory system is an important component of insect behavior and is vital for survival and reproduction. However, the genomic characterization and molecular basis of the olfactory response of Megalurothrips usitatus remain relatively unknown. RNA sequencing-built developmental transcriptomes of nymphs, pupae, and adult M. usitatus were examined in order to establish the sequence-based background of M. usitatus olfactory responses. A total of 56,669 unigenes were annotated using GO, NR, Pfam, eggNOG, SwissProt, and KEGG. The number of differentially expressed genes between pupae and nymphs, males and nymphs, and females and nymphs were 10,498, 9,235, and 10,964, respectively. One odorant-binding protein (MusiOBP1) and one chemosensory protein (MusiCSP1) were selected from the transcriptome, and their full-length sequences were obtained using RACE PCR. The relative expression of MusiOBP1 was the highest in primordial females, whereas the relative expression of MusiCSP1 was the highest in primordial pupae. The strongest binding ability to the odor-binding protein MusiOBP1 was observed for β-citronellol. 3-Hydroxy-2-methyl-4-pyrone showed the strongest binding affinity to MusiCSP1. Our analysis suggests that MusiOBP1 and MusiCSP1 may play significant roles in mediating M. usitatus host recognition. This research will improve our knowledge of odorant-binding proteins and chemosensory proteins, which will in turn improve our understanding of insect olfactory systems.

Differential temperature responses between Plutella xylostella and its specialist endo-larval parasitoid Diadegma semiclausum-Implications for biological control

Understanding the thermal dynamics of host-parasitoid interactions is crucial to predicting how biological control of pest insects by parasitoids might be affected by geographic location and climate change. We compared performance traits of Plutella xylostella (Lepidoptera: Plutellidae) and its solitary endo-larval parasitoid Diadegma semiclausum (Hymenoptera: Ichneumonidae), over a wide range of constant rearing temperatures (10-30°C). Parasitoids reared at 30°C experienced reductions in pupation rate, pupal mass, egg load, and adult life span when compared with those reared at lower temperatures. Our analyses of the fate of parasitoids and their hosts and intergenerational population growth at different rearing temperatures show that D. semiclausum and P. xylostella respond differently to temperature, leading to divergent outcomes under different temperature conditions. Some parasitoid larvae could not complete development at 30°C, the temperature at which the host biomass was least and the metabolic demands of the parasitoid could be high, suggesting that parasitoid development might be constrained by lack of host resources at higher temperatures. We discuss the potential mechanisms of parasitoid susceptibility to elevated temperatures, which likely explain the pronounced seasonal dynamics of D. semiclausum in subtropical regions and its failure to establish in lowland tropical regions, where P. xylostella is a serious pest. Similar interactions in other host-parasitoid associations would constrain the efficacy of parasitoids as biological control agents as global temperatures increase.

Evaluation of antibacterial activity induced by Staphylococcus aureus and Ent A in the hemolymph of Spodoptera littoralis

The problem of antibiotic resistance considers one of the most dangerous challenges facing the medical field. So, it is necessary to find substitutions to conventional antibiotics. Antimicrobial peptides (AMPs) are a bio-functional derivative that have been observed as one of the important solutions to such upcoming crisis. Owing to their role as the first line of defense against bacteria, fungi, and viruses. This study was conducted to induce the immune response of Spodoptera littoralis larvae by inoculation of sub lethal doses of Staphylococcus aureus and its enterotoxin. Since Staphylococcal enterotoxin A (SEA) considers the major causative agents of Staphylococcal food poisoning, our study oriented to purify and characterize this toxin to provoke its role in yielding AMPs with broad spectrum antimicrobial activity. A great fluctuation was recorded in the biochemical properties of immunized hemolymph not only in the total protein content but also protein banding pattern. Protein bands of ∼22 kDa (attacin-like) and ∼15 kDa (lysozyme-like) were found to be common between the AMPs induced as a result of both treatments. While protein bands of molecular weight ∼70 kDa (phenoloxidase-like) and ∼14 kDa (gloverin-like) were found specific for SEA treatment. Chromatographic analysis using HPLC for the induced AMPs showed different types of amino acids appeared with differences in their quantities and velocities. These peptides exhibited noticeable antimicrobial activity against certain Gram-positive and Gram-negative bacteria. In conclusion, the antimicrobial potential of the antimicrobial peptides (AMP) induced in the larval hemolymph of S. littoralis will be a promising molecule for the development of new therapeutic alternatives.

Diverse Host Immune Responses of Different Geographical Populations of the Coconut Rhinoceros Beetle to Oryctes Rhinoceros Nudivirus (OrNV) Infection

Incursions of the coconut rhinoceros beetle (CRB), Oryctes rhinoceros, into different islands in the South Pacific have been detected in recent years. It has been suggested that this range expansion is related to an O. rhinoceros haplotype reported to show reduced susceptibility to the well-established classical biocontrol agent, Oryctes rhinoceros nudivirus (OrNV). Our understanding of the genetic characteristics which distinguish the population of O. rhinoceros that has recently established in Solomon Islands from other well-established populations across the region is very limited. Here, we hypothesized that the recently established O. rhinoceros population should have greater innate immune responses when challenged by OrNV than those of well-established and native O. rhinoceros populations. We used the RNA sequencing (RNA-Seq) approach to generate gene expression profiles of midgut tissue from OrNV-infected and noninfected individuals collected in the Solomon Islands (recent incursion), Papua New Guinea and Fiji (previously established), and the Philippines (within the native range). The collections included individuals from each of the three major mitochondrial lineages (CRB-G, CRB-PNG, and CRB-S) known to the region, allowing us to explore the specific responses of each haplotype to infection. Although insects from the Philippines and Solomon Islands that were tested belong to the same mitochondrial lineage (CRB-G), their overall responses to infection were different. The number of differentially expressed genes between OrNV-infected and noninfected wild-caught individuals from the four different locations varied from 148 to 252. Persistent OrNV infection caused a high level of induced antimicrobial activity and immune responses in O. rhinoceros, but the direction and magnitude of the responses were population specific. The insects tested from the Solomon Islands displayed extremely high expression of genes which are known to be involved in immune responses (e.g. coleoptericin, cecropin, and serpin). These variations in the host immune system among insects from different geographical regions might be driven by variations in the virulence of OrNV isolates, and this requires further investigation. Overall, our current findings support the importance of immunity in insect pest incursion and an expansion of the pest's geographic range. IMPORTANCE Oryctes rhinoceros nudivirus (OrNV) is a double-stranded DNA (dsDNA) virus which has been used as a biocontrol agent to suppress coconut rhinoceros beetle (CRB) in the Pacific Islands. Recently a new wave of CRB incursions in Oceania is thought to be related to the presence of low-virulence isolates of OrNV or virus-tolerant haplotypes of beetles (CRB-G). Our comparative analysis of OrNV-infected and noninfected CRBs revealed that specific sets of genes were induced by viral infection in the beetles. This induction was much stronger in beetles collected from the Solomon Islands, a newly invaded country, than in individuals collected from within the beetle's native range (the Philippines) or from longer-established populations in its exotic range (Fiji and Papua New Guinea [PNG]). Beetles from the Philippines and the Solomon Islands that were tested in this study all belonged to the CRB-G haplotype, but the country-specific responses of the beetles to OrNV infection were different.

Antibiotics accelerate growth at the expense of immunity

Antibiotics have long been used in the raising of animals for agricultural, industrial or laboratory use. The use of subtherapeutic doses in diets of terrestrial and aquatic animals to promote growth is common and highly debated. Despite their vast application in animal husbandry, knowledge about the mechanisms behind growth promotion is minimal, particularly at the molecular level. Evidence from evolutionary research shows that immunocompetence is resource-limited, and hence expected to trade off with other resource-demanding processes, such as growth. Here, we ask if accelerated growth caused by antibiotics can be explained by genome-wide trade-offs between growth and costly immunocompetence. We explored this idea by injecting broad-spectrum antibiotics into wood tiger moth (Arctia plantaginis) larvae during development. We follow several life-history traits and analyse gene expression (RNA-seq) and bacterial (r16S) profiles. Moths treated with antibiotics show a substantial depletion of bacterial taxa, faster growth rate, a significant downregulation of genes involved in immunity and significant upregulation of growth-related genes. These results suggest that the presence of antibiotics may aid in up-keeping the immune system. Hence, by reducing the resource load of this costly process, bodily resources may be reallocated to other key processes such as growth.

The glutathione S-transferase (PxGST2L) may contribute to the detoxification metabolism of chlorantraniliprole in Plutella xylostella(L.)

The diamondback moth (Plutella xylostella L.), is an economic pest of cruciferous plants worldwide, which causes great economic loss to cruciferous plants production. However, the pest has developed resistance to insecticides. One of such insecticides is chlorantraniliprole. The study of the mechanisms underlying resistance is key for the effective management of resistance. In this study, a comparative proteomics approach was used to isolate and identify various proteins that differed between chlorantraniliprole-susceptible and -resistant strains of P. xylostella. Eleven proteins were significantly different and were successfully identified by MALDI-TOF-MS. Metabolism-related proteins accounted for the highest proportion among the eleven different proteins. The function of the PxGST2L protein was validated by RNAi. Knockdown of PxGST2L reduced the GST activity and increased the toxicity of chlorantraniliprole to the diamondback moth. The resistance ratio of diamondback moth to chlorantraniliprole was reduced from 1029 to 505. The results indicated that PxGST2L is partly responsible for chlorantraniliprole insecticide resistance in DBM. Our finding contributes to the understanding of the mechanism underlying resistance to chlorantraniliprole in the DBM, to develop effective resistance management tactics.

Detoxification of plant defensive glucosinolates by an herbivorous caterpillar is beneficial to its endoparasitic wasp

Plant chemical defences impact not only herbivores, but also organisms in higher trophic levels that prey on or parasitize herbivores. While herbivorous insects can often detoxify plant chemicals ingested from suitable host plants, how such detoxification affects endoparasitoids that use these herbivores as hosts is largely unknown. Here, we used transformed plants to experimentally manipulate the major detoxification reaction used by Plutella xylostella (diamondback moth) to deactivate the glucosinolate defences of its Brassicaceae host plants. We then assessed the developmental, metabolic, immune, and reproductive consequences of this genetic manipulation on the herbivore as well as its hymenopteran endoparasitoid Diadegma semiclausum. Inhibition of P. xylostella glucosinolate metabolism by plant-mediated RNA interference increased the accumulation of the principal glucosinolate activation products, the toxic isothiocyanates, in the herbivore, with negative effects on its growth. Although the endoparasitoid manipulated the excretion of toxins by its insect host to its own advantage, the inhibition of herbivore glucosinolate detoxification slowed endoparasitoid development, impaired its reproduction, and suppressed the expression of genes of a parasitoid-symbiotic polydnavirus that aids parasitism. Therefore, the detoxification of plant glucosinolates by an herbivore lowers its toxicity as a host and benefits the parasitoid D. semiclausum at multiple levels.

Transcriptome analysis of blood for the discovery of sex-related genes in ricefield eel Monopterus albus

The blood acts as a transfer channel for a variety of factors in the whole body. The ricefield eel (Monopterus albus) is a protogynous hermaphrodite vertebrate. Until now, no research has reported an analysis of the blood transcriptome during the process of sexual development in the ricefield eel. In this study, the transcriptome sequencing of blood samples from male and female ricefield eels was completed with a total of 34.70 Gb clean data. The clean data of each sample all reached 5.23 GB, and the percent of the Q30 basic group was 88.62% and above. A total of 106,369 unigenes were obtained after assembly, including 13,296 unigenes with a length of more than 1 kb. Further functional annotation analysis showed that there are 28,522 unigenes that can be annotated. The annotations of genes with differential expression revealed that there were 563 genes with significant differential expression in the blood of male and female ricefield eels, including 91 upregulated genes and 472 downregulated genes. Among which, 14 genes may be closely related to sex differentiation, the qPCR was used to confirmed the expression pattern of those genes and result shown that 11 genes were downregulated and 3 genes were upregulated, consistent with the results of our RNA-Seq analysis. This blood transcript dataset will open future research avenues on ricefield eel sex development and differentiation.

Molecular Evidence that Lysiphlebia japonica Regulates the Development and Physiological Metabolism of Aphis gossypii

Lysiphlebia japonica Ashmead (Hymenoptera, Braconidae) is an endophagous parasitoid and Aphis gossypii Glover (Hemiptera, Aphididae) is a major pest in cotton. The relationship between insect host-parasitoids and their hosts involves complex physiological, biochemical and genetic interactions. This study examines changes in the development and physiological metabolism of A. gossypii regulated by L. japonica. Our results demonstrated that both the body length and width increased compared to non-parasitized aphids. We detected significantly increases in the developmental period as well as severe reproductive castration following parasitization by L. japonica. We then used proteomics to characterize these biological changes, and when combined with transcriptomes, this analysis demonstrated that the differential expression of mRNA (up or downregulation) captured a maximum of 48.7% of the variations of protein expression. We assigned these proteins to functional categories that included immunity, energy metabolism and transport, lipid metabolism, and reproduction. We then verified the contents of glycogen and 6-phosphate glucose, which demonstrated that these important energy sources were significantly altered following parasitization. These results uncover the effects on A. gossypii following parasitization by L. japonica, additional insight into the mechanisms behind insect-insect parasitism, and a better understanding of host-parasite interactions.

Variation in Performance and Resistance to Parasitism of Plutella xylostella Populations

Two major ecological factors determine the fitness of an insect herbivore: the ability to overcome plant resistance strategies (bottom-up effects) and the ability to avoid or resist attack by natural enemies such as predators and parasitoids (top-down effects). In response to differences in selection pressure, variation may exist in host-plant adaptation and immunity against parasitism among populations of an insect herbivore. We investigated the variation in larval performance of six different Plutella xylostella populations originating from four continents when feeding on a native Dutch plant species, Brassica rapa. One of the used populations has successfully switched its host plant, and is now adapted to pea. In addition, we determined the resistance to attack by the endoparasitoid Diadegma semiclausum originating from the Netherlands (where it is also native) and measured parasitoid performance as a proxy for host resistance against parasitism. Pupal mortality, immature development times, and adult biomass of P. xylostella differed significantly across populations when feeding on the same host plant species. In addition, parasitism success differed in terms of parasitoid adult emergence and their biomass, but not their development times. Variation among natural populations of insects should be considered more when studying interactions between plants and insects up the food chain.

Regulation of the Larval Transcriptome of Diatraea saccharalis (Lepidoptera: Crambidae) by Maternal and Other Factors of the Parasitoid Cotesia flavipes (Hymenoptera: Braconidae)

Koinobiont endoparasitoid wasps regulate the host's physiology to their own benefit during their growth and development, using maternal, immature and/or derived-tissue weaponry. The tools used to subdue the wasps' hosts interfere directly with host transcription activity. The broad range of host tissues and pathways affected impedes our overall understanding of the host-regulation process during parasitoid development. Next-generation sequencing and de novo transcriptomes are helpful approaches to broad questions, including in non-model organisms. In the present study, we used Illumina sequencing to assemble a de novo reference transcriptome of the sugarcane borer Diatraea saccharalis, to investigate the regulation of host gene expression by the larval endoparasitoid Cotesia flavipes. We obtained 174,809,358 reads and assembled 144,116 transcripts, of which 44,325 were putatively identified as lepidopteran genes and represented a substantial number of pathways that are well described in other lepidopteran species. Comparative transcriptome analyses of unparasitized versus parasitized larvae identified 1,432 transcripts of D. saccharalis that were up-regulated under parasitization by C. flavipes, while 1,027 transcripts were down-regulated. Comparison of the transcriptomes of unparasitized and pseudoparasitized D. saccharalis larvae led to the identification of 1,253 up-regulated transcripts and 972 down-regulated transcripts in the pseudoparasitized larvae. Analysis of the differentially expressed transcripts showed that C. flavipes regulated several pathways, including the Ca+2 transduction signaling pathway, glycolysis/gluconeogenesis, chitin metabolism, and hormone biosynthesis and degradation, as well as the immune system, allowing us to identify key target genes involved in the metabolism and development of D. saccharalis.

Immune response and susceptibility to Cotesia flavipes parasitizing Diatraea saccharalis larvae exposed to and surviving an LC25 dosage of Bacillus thuringiensis

Biological control using entomopathogens and natural enemies is an ecofriendly method for pest management in agriculture. Biological control agents often can be simultaneously employed and compatibility between agents may improve pest suppression. We investigated the influence of the entomopathogen Bacillus thuringiensis (Bt) on the immune system of the sugarcane borer Diatraea saccharalis (Fabricius, 1794) (Lepidoptera: Crambidae) to determine if such changes impact parasitization by Cotesia flavipes Cameron, 1891 (Hymenoptera: Braconidae). The immune response of surviving D. saccharalis larvae fed with an LC₂₅ dosage of a Bt-based biopesticide (Dipel®) was analyzed (total hemocyte count, hemocyte adhesion, and activities of phenoloxidase and lysozyme). Furthermore, the suitability of surviving Bt-fed larvae as hosts for C. flavipes was assessed by measuring parasitoid attributes such as number and size of teratocytes, weight of pupae, length of adult female tibia and number of emerged adults. Total hemocyte count, but not hemocyte adhesion, total protein content and phenoloxidase activity increased in the hemolymph of non-parasitized Bt-fed larvae (Bt-NP) compared to control larvae (NBt-NP). Lysozyme activity increased only after parasitization without Bt exposure (NBt-P). After parasitization, the immunological parameters activated in Bt-NP larvae decreased to levels at or below those observed in control larvae, showing that C. flavipes can regulate the activated immune response of Bt-fed larvae. The development of C. flavipes was not impaired in Bt-fed larval hosts (Bt-P); no changes were observed for teratocyte size, weight of pupal mass, length of hind tibia and number of adults emerged.

RNA-seq profiles of putative genes involved in specific immune priming in Bombyx mori haemocytes

BACKGROUND

The immune system of many invertebrates, including insects, has been shown to comprise memory, or specific immune priming. However, knowledge of the molecular mechanisms especially the candidate immune-related genes mediated the specificity of the immune priming are still very scarce and fragmentary. We therefore used two closely related Gram-negative pathogenic bacteria (Photorhabdus luminescens TT01 and P. luminescens H06) as the priming agents and employed Illumina/Solexa platform to investigate the transcriptional changes of the haemocytes of Bombyx mori larvae after priming.

RESULTS

In total, 23.0 Gbp of sequence data and 153,331,564 reads were generated, representing 10,496 genes. Approximately 89% of the genes or sequenced reads could be aligned to the silkworm reference genome. The differentially expressed genes (DEGs) of PBS-vs-TT01 (up-regulated expression of TT01 relative to PBS), PBS-vs-H06 (up-regulated expression of H06 relative to PBS) and TT01-vs-H06 (up-regulated expression of H06 relative to TT01) were 707, 159 and 461 respectively. In addition, expression patterns of 25 selected DEGs derived from quantitative real-time polymerase chain reaction (qRT-PCR) were consistent with their transcript abundance changes obtained by transcriptomic analyses. The DEGs are mainly related to pattern recognition receptors (PRRs), antimicrobial peptides (AMPs), signaling molecular, effector molecules, phagosome and spliceosome, indicating that they have participated in the regulation of the specific immune priming in the B. mori larvae.

CONCLUSIONS

The transcriptome profiling data sets from this study will provide valuable resources to better understand the molecular and biological mechanisms regulating the specificity of invertebrates' immune priming. All these will shed light on controlling insect pests or preventing epidemic of infectious diseases in economic invertebrates.

The dual life of ichnoviruses

Ichnoviruses (IVs) are mutualistic, double-stranded DNA viruses playing a key role in the successful parasitism of thousands of endoparasitoid wasp species. IV particles are produced exclusively in the female wasp reproductive tract. They are co-injected along with the parasitoid egg into caterpillar hosts upon parasitization. The expression of viral genes by infected host cells leads to an immunosuppressive state and delayed development of the host, two pathologies that are critical to the successful development of the wasp egg and larva. Ichnovirus is one of the two recognized genera within the family Polydnaviridae (polydnaviruses or PDVs), the other genus being Bracovirus (BV), associated with braconid wasps. IVs are associated with ichneumonid wasps belonging to the subfamilies Campopleginae and Banchinae; attempts to identify IV particles in other ichneumonid subfamilies have so far been unsuccessful. Functional studies targeting IV genes expressed in parasitized hosts, along with investigations of the molecular mechanisms responsible for viral morphogenesis in the female wasp, have resulted in a better understanding of the biology of these atypical viruses.

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.

Comparative transcriptomic analysis provides insights into antibacterial mechanisms of Branchiostoma belcheri under Vibrio parahaemolyticus infection

Amphioxus, a basal chordate, is widely considered to be an existing proxy of the invertebrate ancestor of vertebrates, and it exhibits susceptibility to various pathogen infections and pathogenic mimic challenges. Here, in order to understand more clearly its antibacterial mechanisms, we analyzed the ribosomal RNA (rRNA)-depleted transcriptome of Chinese amphioxus (Branchiostoma belcheri) infected with Vibrio parahaemolyticus (V. p.) via next-generation deep sequencing technology (RNA-seq). We identified a total of 3214 differentially expressed genes (DEGs) by comparing V. p.-infected and control transcriptome libraries, including 2219 significantly up-regulated and 995 significantly down-regulated DEGs in V. p.-infected amphioxus. The DEGs with the top 10 most dramatic expression fold changes after V. p. infection, as well as 53 immune-related DEGs (IRDs) belonging to four primary categories of innate immunity were analyzed further. Through gene ontology (GO) and pathway enrichment analysis, DEGs were found to be primarily related to immune processes, apoptosis, catabolic and metabolic processes, binding and enzyme activity, while pathways involving bacterial infection, immune signaling, immune response, cancer, and apoptosis were overrepresented. We validated the RNA-seq results by detecting the expression levels of 10 IRDs using qRT-PCR, and we surveyed the dynamic variation in gene expression for these IRDs at 0, 6, 12, 24, and 48 h after V. p.

TREATMENT

Subsequently, according to the RNA-seq results, the presence of a primitive Toll-like receptor (TLR)-mediated antibacterial immune signaling pathway was predicted in B. belcheri. This study provides valuable information regarding antibacterial immunity for further research into the evolution of immunity in vertebrates and broadens our understanding of the innate immune response against bacterial invasion in amphioxus.

Transcriptome analysis in different developmental stages of Batocera horsfieldi (Coleoptera: Cerambycidae) and comparison of candidate olfactory genes

The white-striped longhorn beetle Batocera horsfieldi (Coleoptera: Cerambycidae) is a polyphagous wood-boring pest that causes substantial damage to the lumber industry. Moreover olfactory proteins are crucial components to function in related processes, but the B. horsfieldi genome is not readily available for olfactory proteins analysis. In the present study, developmental transcriptomes of larvae from the first instar to the prepupal stage, pupae, and adults (females and males) from emergence to mating were built by RNA sequencing to establish a genetic background that may help understand olfactory genes. Approximately 199 million clean reads were obtained and assembled into 171,664 transcripts, which were classified into 23,380, 26,511, 22,393, 30,270, and 87, 732 unigenes for larvae, pupae, females, males, and combined datasets, respectively. The unigenes were annotated against NCBI’s non-redundant nucleotide and protein sequences, Swiss-Prot, Gene Ontology (GO), Pfam, Clusters of Eukaryotic Orthologous Groups (KOG), and KEGG Orthology (KO) databases. A total of 43,197 unigenes were annotated into 55 sub-categories under the three main GO categories; 25,237 unigenes were classified into 26 functional KOG categories, and 25,814 unigenes were classified into five functional KEGG Pathway categories. RSEM software identified 2,983, 3,097, 870, 2,437, 5,161, and 2,882 genes that were differentially expressed between larvae and males, larvae and pupae, larvae and females, males and females, males and pupae, and females and pupae, respectively. Among them, genes encoding seven candidate odorant binding proteins (OBPs) and three chemosensory proteins (CSPs) were identified. RT-PCR and RT-qPCR analyses showed that BhorOBP3, BhorCSP2, and BhorOBPC1/C3/C4 were highly expressed in the antenna of males, indicating these genes may may play key roles in foraging and host-orientation in B. horsfieldi. Our results provide valuable molecular information about the olfactory system in B. horsfieldi and will help guide future functional studies on olfactory genes.

Transcriptome-wide analysis of immune-responsive microRNAs against poly (I:C) challenge in Branchiostoma belcheri by deep sequencing and bioinformatics

Amphioxus is a key experimental animal for studying the evolution of vertebrate immune system. However, we still do not know about the roles of microRNAs (miRNAs) under viral stress in amphioxus. In this study, we sequenced six small RNA libraries (three biological replicates were included in the treatments challenged by the viral mimic, poly (I:C) (pIC) and control groups, respectively) from Branchiostoma belcheri. A total of 151 known miRNAs, 197 new miRNAs (named novel_mir, including nine conserved miRNAs) were identified by deep sequencing from the six libraries. We primarily focused on differentially expressed miRNAs (DEMs) after pIC challenge. Next, we screened a total of 77 DEMs, including 27 down- and 50 up-regulated DEMs in response to pIC challenge. Furthermore, we used real-time quantitative PCR (qRT-PCR) to verify the expression levels of 10 randomly selected DEMs. Target genes likely regulated by DEMs were predicted, and functional enrichment analyses of these targets were performed using bioinformatics approach. MiRNA targets of DEMs are primarily involved in immune response, diseases, cancer and regulation process, and could be largely linked to 14 immune-related signaling pathways, including NF-kappa B, NOD-like receptor, RIG-I-like receptor and endocytosis. The present study for the first time explores key regulatory roles of miRNAs in the innate antiviral immune response in amphioxus, and will provide insight into the molecular basis of antiviral immunity and evolution of immune-related miRNAs.

The Entomopathogenic Fungi Isaria fumosorosea Plays a Vital Role in Suppressing the Immune System of Plutella xylostella: RNA-Seq and DGE Analysis of Immunity-Related Genes

Most, if not all, entomopathogenic fungi have been used as alternative control agents to decrease the insect resistance and harmful effects of the insecticides on the environment. Among them, Isaria fumosorosea has also shown great potential to control different insect pests. In the present study, we explored the immune response of P. xylostella to the infection of I. fumosorosea at different time points by using RNA-Sequencing and differential gene expression technology at the genomic level. To gain insight into the host-pathogen interaction at the genomic level, five libraries of P. xylostella larvae at 12, 18, 24, and 36 h post-infection and a control were constructed. In total, 161 immunity-related genes were identified and grouped into four categories; immune recognition families, toll and Imd pathway, melanization, and antimicrobial peptides (AMPs). The results of differentially expressed immunity-related genes depicted that 15, 13, 53, and 14 up-regulated and 38, 51, 56, and 49 were down-regulated in P. xylostella at 12, 18, 24, and 36 h post-treatment, respectively. RNA-Seq results of immunity-related genes revealed that the expression of AMPs was reduced after treatment with I. fumosorosea. To validate RNA-Seq results by RT-qPCR, 22 immunity-related genes were randomly selected. In conclusion, our results demonstrate that I. fumosorosea has the potential to suppress the immune response of P. xylostella and can become a potential biopesticide for controlling P. xylostella.

The developmental transcriptome of the bamboo snout beetle Cyrtotrachelus buqueti and insights into candidate pheromone-binding proteins

Cyrtotrachelus buqueti is an extremely harmful bamboo borer, and the larvae of this pest attack clumping bamboo shoots. Pheromone-binding proteins (PBPs) play an important role in identifying insect sex pheromones, but the C. buqueti genome is not readily available for PBP analysis. Developmental transcriptomes of eggs, larvae from the first instar to the prepupal stage, pupae, and adults (females and males) from emergence to mating were built by RNA sequencing (RNA-Seq) in the present study to establish a sequence background of C. buqueti to help understand PBPs. Approximately 164.8 million clean reads were obtained and annotated into 108,854 transcripts. These were assembled into 24,338, 21,597, 24,798, 21,886, 24,642, and 83,115 unigenes for eggs, larvae, pupae, females, males, and the combined datasets, respectively. Unigenes were annotated against NCBI non-redundant protein sequences, NCBI non-redundant nucleotide sequences, Gene Ontology (GO), Protein family, Clusters of Orthologous Groups of Proteins/ Clusters of Eukaryotic Orthologous Groups (KOG), Swiss-Prot, and KEGG Orthology databases. A total of 17,213 unigenes were annotated into 55 sub-categories belonging to three main GO categories; 10,672 unigenes were classified into 26 functional categories by KOG classification, and 8,063 unigenes were classified into five functional KEGG categories. RSEM software for RNA sequencing showed that 4,816, 3,176, 3,661, 2,898, 4,316, 8,019, 7,273, 5,922, 5,844, and 4,570 genes were differentially expressed between larvae and males, larvae and eggs, larvae and pupae, larvae and females, males and females, males and eggs, males and pupae, females and eggs, females and pupae, and eggs and pupae, respectively. Of these, three were confirmed to be significantly differentially expressed between larvae, females, and males. Furthermore, PBP Cbuq7577_g1 was highly expressed in the antenna of males. A comprehensive sequence resource of a desirable quality was constructed from developmental transcriptomes of C. buqueti eggs, larvae, pupae, and adults. This work enriches the genomic data of C. buqueti, and facilitates our understanding of its metamorphosis, development, and response to environmental change. The identified candidate PBP Cbuq7577_g1 might play a crucial role in identifying sex pheromones, and could be used as a targeted gene to control C. buqueti numbers by disrupting sex pheromone communication.

Characterization of joining sites of a viral histone H4 on host insect chromosomes

A viral histone H4 (CpBV-H4) is encoded in a polydnavirus, Cotesia plutellae bracovirus (CpBV). It plays a crucial role in parasitism of an endoparasitoid wasp, C. plutellae, against diamondback moth, Plutella xylostella, by altering host gene expression in an epigenetic mode by its N-terminal tail after joining host nucleosomes. Comparative transcriptomic analysis between parasitized and nonparasitized P. xylostella by RNA-Seq indicated that 1,858 genes were altered at more than two folds in expression levels at late parasitic stage, including 877 up-regulated genes and 981 down-regulated genes. Among parasitic factors altering host gene expression, CpBV-H4 alone explained 16.3% of these expressional changes. To characterize the joining sites of CpBV-H4 on host chromosomes, ChIP-Seq (chromatin immunoprecipitation followed by deep sequencing) was applied to chromatins extracted from parasitized larvae. It identified specific 538 ChIP targets. Joining sites were rich (60.2%) in AT sequence. Almost 40% of ChIP targets included short nucleotide repeat sequences presumably recognizable by transcriptional factors and chromatin remodeling factors. To further validate these CpBV-H4 targets, CpBV-H4 was transiently expressed in nonparasitized host at late larval stage and subjected to ChIP-Seq. Two kinds of ChIP-Seqs shared 51 core joining sites. Common targets were close (within 1 kb) to genes regulated at expression levels by CpBV-H4. However, other host genes not close to CpBV-H4 joining sites were also regulated by CpBV-H4. These results indicate that CpBV-H4 joins specific chromatin regions of P. xylostella and controls about one sixth of the total host genes that were regulated by C. plutellae parasitism in an epigenetic mode.

An endoparasitoid wasp influences host DNA methylation

Parasitism by endoparasitoid wasps changes the expression of various host genes, and alters host immune and developmental processes. However, it is not clearly understood how parasitism changes host gene expression in a whole genome scale. This study focused on an epigenetic control of Cotesia plutellae, an endoparasitoid wasp, against its host, Plutella xylostella. Two DNA methyltransferases (DNMT-1 and DNMT-2) are encoded in the genome of P. xylostella. In addition, methyl-binding domain proteins (MBDs) and DNA demethylation factor, ten-eleven translation protein (TET) are encoded. DNA methylation of P. xylostella genomic DNA was confirmed by restriction digestion with Gla I specific to 5-methylcytosine. DNA methylation intensity in parasitized (P) larvae was decreased compared to that in nonparasitized (NP) larvae, especially at late parasitic stage, at which expression levels of both DNMT-1 and DNMT-2 were also decreased. DNA demethylation of P. xylostella was confirmed in both NP and P larvae by restriction digestion with PvuRts1I recognizing 5-hydroxymethyl cytosine. Parasitism also suppressed expression levels of TET and MBDs. Treatment of 5-aza-2′-deoxycytidine (AZA) reduced DNA methylation intensity of NP larvae, causing suppression of hemocyte-spreading behavior and delay of immature development. RNA interference of DNMT-1 or DNMT-2 mimicked the adverse effects of AZA.

Revised annotation of Plutella xylostella microRNAs and their genome-wide target identification

The diamondback moth, Plutella xylostella, is the most devastating pest of brassica crops worldwide. Although 128 mature microRNAs (miRNAs) have been annotated from this species in miRBase, there is a need to extend and correct the current P. xylostella miRNA repertoire as a result of its recently improved genome assembly and more available small RNA sequence data. We used our new ultra-deep sequence data and bioinformatics to re-annotate the P. xylostella genome for high confidence miRNAs with the correct 5p and 3p arm features. Furthermore, all the P. xylostella annotated genes were also screened to identify potential miRNA binding sites using three target-predicting algorithms. In total, 203 mature miRNAs were annotated, including 33 novel miRNAs. We identified 7691 highly confident binding sites for 160 pxy-miRNAs. The data provided here will facilitate future studies involving functional analyses of P. xylostella miRNAs as a platform to introduce novel approaches for sustainable management of this destructive pest.

Translational Control of Host Gene Expression by a Cys-Motif Protein Encoded in a Bracovirus

Translational control is a strategy that various viruses use to manipulate their hosts to suppress acute antiviral response. Polydnaviruses, a group of insect double-stranded DNA viruses symbiotic to some endoparasitoid wasps, are divided into two genera: ichnovirus (IV) and bracovirus (BV). In IV, some Cys-motif genes are known as host translation-inhibitory factors (HTIF). The genome of endoparasitoid wasp Cotesia plutellae contains a Cys-motif gene (Cp-TSP13) homologous to an HTIF known as teratocyte-secretory protein 14 (TSP14) of Microplitis croceipes. Cp-TSP13 consists of 129 amino acid residues with a predicted molecular weight of 13.987 kDa and pI value of 7.928. Genomic DNA region encoding its open reading frame has three introns. Cp-TSP13 possesses six conserved cysteine residues as other Cys-motif genes functioning as HTIF. Cp-TSP13 was expressed in Plutella xylostella larvae parasitized by C. plutellae. C. plutellae bracovirus (CpBV) was purified and injected into non-parasitized P. xylostella that expressed Cp-TSP13. Cp-TSP13 was cloned into a eukaryotic expression vector and used to infect Sf9 cells to transiently express Cp-TSP13. The synthesized Cp-TSP13 protein was detected in culture broth. An overlaying experiment showed that the purified Cp-TSP13 entered hemocytes. It was localized in the cytosol. Recombinant Cp-TSP13 significantly inhibited protein synthesis of secretory proteins when it was added to in vitro cultured fat body. In addition, the recombinant Cp-TSP13 directly inhibited the translation of fat body mRNAs in in vitro translation assay using rabbit reticulocyte lysate. Moreover, the recombinant Cp-TSP13 significantly suppressed cellular immune responses by inhibiting hemocyte-spreading behavior. It also exhibited significant insecticidal activities by both injection and feeding routes. These results indicate that Cp-TSP13 is a viral HTIF.

De novo sequencing and transcriptome analysis of venom glands of endoparasitoid Aenasius arizonensis (Girault) (=Aenasius bambawalei Hayat) (Hymenoptera, Encyrtidae)

Aenasius bambawalei Hayat (Encyrtidae: Hymenoptera) has been synonymized with Aenasius arizonensis (Girault) is a small, newly discovered endoparasitoid of the cotton mealybug Phenacoccuss solenopsis Tinsley (Pseudococcidae: Hemiptera), which completes its life cycle inside the body of its host and it is a potential insect control tool. Despite the acquired knowledge regarding host-parasitoid interaction, little information is available on the factors of parasitoid origin able to modulate mealybug physiology. The components of A. arizonensis venom have not been well studied but venom from other parasitoids and wasps contain biologically active proteins that have potential applications in pest management or may be of medicinal importance. To provide an insight into the transcripts expressed in the venom gland of A. arizonensis, a transcriptomic database was developed utilizing high throughput RNA sequencing approaches to analyze the genes expressed in venom glands of this endoparasitic wasp. The resulting A. arizonensis RNA sequences were assembled de-novo with contigs then blasted against the NCBI non-redundant sequence database. Contigs which matched database sequences were mostly homologous to genes from hymenopteran parasitoids such as Nasonia vitripennis, Copidosoma floridanum, Fopius arsenus and Pteromalas puparium. Further analysis of the A. arizonensis database was then performed which focused on selected genes encoding proteins potentially involved in host developmental arrest, disrupting the host immune system, host paralysis, and transcripts that support these functions. Sequenced mRNAS predicted to encode full length ORFs of Calreticulin, Serine Protease Precursor and Arginine kinase proteins were identified and the tissue specific expression of these putative venom genes was analyzed by RT-PCR. In addition, results also demonstrate that de novo transcriptome assembly allows useful venom gene expression analysis in a species lacking a genome sequence database and may provide useful information for devising control tools for insect pests and other applications.

De novo sequencing and transcriptome analysis of female venom glands of ectoparasitoid Bracon hebetor (Say.) (Hymenoptera: Braconidae)

Venom is a key-factor in the regulation of host physiology by parasitic Hymenoptera and a potentially rich source of novel bioactive substances for biotechnological applications. The limited study of venom from the ectoparasitoid Bracon hebetor, a tiny wasp that attacks larval pest insects of field and stored products and is thus a potential insect control agent, has not described the full complement and composition of these biomolecules. To have a comprehensive picture of genes expressed in the venom glands of B. hebetor, a venom gland transcriptome was assembled by using next generation sequencing technologies followed by de novo assemblies of the 10.81 M sequence reads yielded 22,425 contigs, of which 10,581 had significant BLASTx hits to know genes. The majority of hits were to Diachasma alloeum, an ectoparasitoid from same taxonomic family, as well as other wasps. Gene ontology grouped the sequences into molecular functions in which catalytic activity with 42.2% was maximum, cellular components in which cells with 33.8% and biological processes among which metabolic process with 30% had the most representatives. In this study, we highlight the most abundant sequences, and those that are likely to be functional components of the venom for parasitization. Full length ORFs of Calreticulin, Venom Acid Phosphatase Acph-1 like protein and arginine kinase proteins were isolated and their tissue specific expression was studied by RT-PCR. Our report is the first to characterize components of the B. hebetor venom glands that may be useful for developing control tools for insect pests and other applications.

De Novo Assembly and Developmental Transcriptome Analysis of the Small White Butterfly Pieris rapae

The small white butterfly Pieris rapae is one of the most destructive pests of Brassicaceae. Yet little is understood about its genes involved in development. To facilitate research on P. rapae, we sequenced the transcriptome of P. rapae during six developmental stages, including the egg, three larval stages, the pupa, and the adult. In total, 240 million high-quality reads were obtained. De novo assembly generated 96,069 unigenes with an average length of 1353 nt. Of these, 31,629 unigenes had homologs as determined by a blastx search against the NR database with a cut-off e-value of 10⁻⁵. Clusters of Orthologous Groups of proteins (COG), Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were conducted to functionally annotate those genes. Then, 849 genes involved in seven canonical development signaling pathway were identified, including dozens of key genes such as Hippo, Notch, and JAK2. A total of 21,883 differentially expressed (cut-off of 2-fold) unigenes were detected across the developmental stages, most of which were found between the egg and first larval stages. Interestingly, only 34 differentially expressed unigenes, most of which are cuticle protein related genes, were detected with a cut-off of 2¹⁰-fold. Furthermore, we identified 32 heat shock protein (Hsp) genes that were expressed with complete open reading frames. Based on phylogenetic trees of the Hsp genes, we found that Hsp genes with close evolutionary relationships had similar expression pattern. Additionally, partial pattern recognition receptors genes were found to be developmental regulated. This study provides comprehensive sequence resources for P. rapae and numerous differential expressed genes, and these findings will lay the foundation for future functional genomics studies on this species.

Deep Sequencing-Based Transcriptome Analysis Reveals the Regulatory Mechanism of Bemisia tabaci (Hemiptera: Aleyrodidae) Nymph Parasitized by Encarsia sophia (Hymenoptera: Aphelinidae)

The whitefly Bemisia tabaci is a genetically diverse complex with multiple cryptic species, and some are the most destructive invasive pests of many ornamentals and crops worldwide. Encarsia sophia is an autoparasitoid wasp that demonstrated high efficiency as bio-control agent of whiteflies. However, the immune mechanism of B. tabaci parasitization by E. sophia is unknown. In order to investigate immune response of B. tabaci to E. Sophia parasitization, the transcriptome of E. sophia parasitized B. tabaci nymph was sequenced by Illumina sequencing. De novo assembly generated 393,063 unigenes with average length of 616 bp, in which 46,406 unigenes (15.8% of all unigenes) were successfully mapped. Parasitization by E. sophia had significant effects on the transcriptome profile of B. tabaci nymph. A total of 1482 genes were significantly differentially expressed, of which 852 genes were up-regulated and 630 genes were down-regulated. These genes were mainly involved in immune response, development, metabolism and host signaling pathways. At least 52 genes were found to be involved in the host immune response, 33 genes were involved in the development process, and 29 genes were involved in host metabolism. Taken together, the assembled and annotated transcriptome sequences provided a valuable genomic resource for further understanding the molecular mechanism of immune response of B. tabaci parasitization by E. sophia.

Glyceraldehyde-3-phosphate dehydrogenase is a mediator of hemocyte-spreading behavior and molecular target of immunosuppressive factor CrV1

Cellular immunity is accompanied by hemocyte-spreading behavior, which undergoes cytoskeletal rearrangement. Polydnaviral factor CpBV-CrV1 can inhibit the hemocyte-spreading behavior and suppress host immune response of Plutella xylostella. However, host target molecule of CpBV-CrV1 that inhibits the hemocyte behavior has not been identified yet. This study used a pull-down approach to identify the target molecule of CpBV-CrV1. A protein bound to CpBV-CrV1 was co-precipitated and identified to be glyceraldehyde-3-phosphate dehydrogenase (GAPDH) by LC-MS/MS analysis. RNA interference (RNAi) specific to GAPDH of P. xylostella was found to be able to inhibit the hemocyte-spreading behavior, while RNAi treatments with other glycolytic genes had no effect on the spreading behavior. An addition of recombinant CpBV-CrV1 on hemocyte monolayer interrupted the association between GAPDH and α-tubulin in the cytoplasm. Overlay of mutant proteins (Y492A or Y501A with tyrosine to alanine at putative GAPDH-binding site) of CpBV-CrV1 on hemocyte monolayer revealed that they could enter hemocytes unlike a mutant in the N-terminal coiled-coil domain. However, they failed to inhibit the hemocyte-spreading behavior without any binding affinity to GAPDH. These results suggest that GAPDH plays a critical role in hemocyte-spreading behavior during immune challenge as a molecular target of viral factor CpBV-CrV1.

Venom gland components of the ectoparasitoid wasp, Anisopteromalus calandrae

The wasp Anisopteromalus calandrae is a small ectoparasitoid that attacks stored product pest beetle larvae that develop inside grain kernels, and is thus a potential insect control tool. The components of A. calandrae venom have not been studied, but venom from other organisms contains proteins with potential applications, such as pest management tools and treatments for human diseases. We dissected female A. calandrae and collected venom and associated glands. Using high throughput sequencing, a venom gland transcriptome was assembled that contained 45,432 contigs, 25,726 of which had BLASTx hits. The majority of hits were to Nasonia vitripennis, an ectoparasitoid from the same taxonomic family, as well as other bees, wasps, and ants. Gene ontology grouped sequences into eleven molecular functions, among which binding and catalytic activity had the most representatives. In this study, we highlighted the most abundant sequences, including those that are likely the functional components of the venom. Specifically, we focused on genes encoding proteins potentially involved in host developmental arrest, disrupting the host immune system, host paralysis, and transcripts that support these functions. Our report is the first to characterize components of the A. calandrae venom gland that may be useful as control tools for insect pests and other applications.

De novo transcriptome sequencing and analysis of Coccinella septempunctata L. in non-diapause, diapause and diapause-terminated states to identify diapause-associated genes

BACKGROUND

The most common ladybird beetle, Coccinella septempunctata L., is an excellent predator of crop pests such as aphids and white flies, and it shows a wide range of adaptability, a large appetite and a high reproductive ability. Diapause research plays an important role in the artificial propagation and shelf-life extension of insect products. Although this lady beetle's regulatory, physiological and biochemical characteristics in the diapause period are well understood, the molecular mechanism of diapause remains unknown. Therefore, we collected female adults in three different states, i.e., non-diapause, diapause and diapause termination, for transcriptome sequencing.

RESULTS

After transcriptome sequencing using the Illumina HiSeq 2500 platform with pretreatment, a total of 417.6 million clean reads from nine samples were filtered using the program FASTX (version 0.0). Additionally, 106,262 contigs were assembled into 82,820 unigenes with an average length of 921 bp and an N50 of 1,241 bp. All of the unigenes were annotated through BLASTX alignment against the Nr or UniProt database, and 37,872 unigenes were matched. We performed further analysis of these unigenes using the Clusters of Orthologous Groups of proteins (COG), Gene Ontology (GO), and the Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. Through pairwise comparisons of the non-diapause (ND), diapause (D), and diapause-terminated (DT) groups, 3,501 and 1,427 differentially expressed genes (DEGs) were identified between D and ND and between DT and D, respectively. Moreover, 443 of the DEGs were specifically expressed during the diapause period (i.e., DEGs that were expressed at the highest or lowest levels during diapause compared with the other stages). GO function and KEGG pathway enrichment were performed on all DEGs and showed that RNA-directed DNA polymerase activity and fatty acid metabolism were significantly affected. Furthermore, eight specific expressed genes were selected for validation using qRT-PCR. Among these eight genes, seven genes were up-regulated, and one gene was down-regulated; the change trends of the eight genes were the same between the qRT-PCR and RNA-seq analysis results.

CONCLUSIONS

In this study, a new method for collecting and identifying diapause insects was described. We generated a vast quantity of transcriptome data from C. septempunctata L., providing a resource for gene function research. The diapause-associated genes that we identified establish a foundation for future studies on the molecular mechanisms of diapause.

THE ENDOPARASITOID Pteromalus puparum INFLUENCES HOST GENE EXPRESSION WITHIN FIRST HOUR OF PARASITIZATION

The small cabbage butterfly, Pieris rapae, is an important pest of cruciferous corps, and Pteromalus puparum is a predominant pupal endoparasitoid wasp of this butterfly. For successful development of parasitoid offspring, female parasitoids usually introduce one or several kinds of maternal factors into the hemocoels during oviposition to suppress host immunity. To investigate the early changes in host immune-related genes following parasitization, we analyzed transcriptomes of parasitized and unparasitized, control, host pupae. Approximately 17.7 and 19.3 million paired-end reads were generated from nonparasitized and parasitized host pupae, and assembled de novo into 45,639 transcripts and 27,659 nonredundant unigenes. The average unigene length was 790 bp. A total 18,377 of 27,659 unigenes were annotated and we identified 557 differentially expressed unigenes in host pupae at 1 h after parasitization, of which 21 were immune-related. Parasitization led to downregulation of most pattern recognition receptors and upregulation of all serine protease inhibitors. The transcirptomic profile of P. rapae is considerably affected by parasitization. This study provides valuable sources for future investigations of the molecular interaction between P. puparum and its host P. rapae.

Genome wide discovery of long intergenic non-coding RNAs in Diamondback moth (Plutella xylostella) and their expression in insecticide resistant strains

Long non-coding RNAs (lncRNAs) play important roles in genomic imprinting, cancer, differentiation and regulation of gene expression. Here, we identified 3844 long intergenic ncRNAs (lincRNA) in Plutella xylostella, which is a notorious pest of cruciferous plants that has developed field resistance to all classes of insecticides, including Bacillus thuringiensis (Bt) endotoxins. Further, we found that some of those lincRNAs may potentially serve as precursors for the production of small ncRNAs. We found 280 and 350 lincRNAs that are differentially expressed in Chlorpyrifos and Fipronil resistant larvae. A survey on P. xylostella midgut transcriptome data from Bt-resistant populations revealed 59 altered lincRNA in two resistant strains compared with the susceptible population. We validated the transcript levels of a number of putative lincRNAs in deltamethrin-resistant larvae that were exposed to deltamethrin, which indicated that this group of lincRNAs might be involved in the response to xenobiotics in this insect. To functionally characterize DBM lincRNAs, gene ontology (GO) enrichment of their associated protein-coding genes was extracted and showed over representation of protein, DNA and RNA binding GO terms. The data presented here will facilitate future studies to unravel the function of lincRNAs in insecticide resistance or the response to xenobiotics of eukaryotic cells.

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.

Genome-wide characterization and expression profiling of immune genes in the diamondback moth, Plutella xylostella (L.)

The diamondback moth, Plutella xylostella (L.), is a destructive pest that attacks cruciferous crops worldwide. Immune responses are important for interactions between insects and pathogens and information on these underpins the development of strategies for biocontrol-based pest management. Little, however, is known about immune genes and their regulation patterns in P. xylostella. A total of 149 immune-related genes in 20 gene families were identified through comparison of P. xylostella genome with the genomes of other insects. Complete and conserved Toll, IMD and JAK-STAT signaling pathways were found in P. xylostella. Genes involved in pathogen recognition were expanded and more diversified than genes associated with intracellular signal transduction. Gene expression profiles showed that the IMD pathway may regulate expression of antimicrobial peptide (AMP) genes in the midgut, and be related to an observed down-regulation of AMPs in experimental lines of insecticide-resistant P. xylostella. A bacterial feeding study demonstrated that P. xylostella could activate different AMPs in response to bacterial infection. This study has established a framework of comprehensive expression profiles that highlight cues for immune regulation in a major pest. Our work provides a foundation for further studies on the functions of P. xylostella immune genes and mechanisms of innate immunity.

Specificities of ichnoviruses associated with campoplegine wasps: genome, genes and role in host-parasitoid interaction

Ichnoviruses (IVs), unique symbiotic viruses carried by ichneumonid campoplegine wasps, derive from integration of a paleo-ichnovirus into an ancestral wasp genome. The modern 'genome' is composed of both regions that are amplified, circularized and encapsidated into viral particles and non-encapsidated viral genomic regions involved in particle morphogenesis. Packaged genomes include multiple circular dsDNAs encoding many genes mostly organized in gene families. Virus particles are assembled in specialized ovarian cells from which they exit into the oviduct lumen; mature virions are injected during oviposition into the insect host. Expression of viral proteins in infected cells correlates with physiological alterations of the host enabling success of parasitism.

A comparison of next generation sequencing technologies for transcriptome assembly and utility for RNA-Seq in a non-model bird

De novo assembled transcriptomes, in combination with RNA-Seq, are powerful tools to explore gene sequence and expression level in organisms without reference genomes. Investigators must first choose which high throughput sequencing platforms will provide data most suitable for their experimental goals. In this study, we explore the utility of 454 and Illumina sequences for de novo transcriptome assembly and downstream RNA-Seq applications in a reproductive gland from a non-model bird species, the Japanese quail (Coturnix japonica). Four transcriptomes composed of either pure 454 or Illumina reads or mixtures of read types were assembled and evaluated for the same cost. Illumina assemblies performed best for de novo transcriptome characterization in terms of contig length, transcriptome coverage, and complete assembly of gene transcripts. Improvements over the Hybrid assembly were marginal, with the exception that the addition of 454 data significantly increased the number of genes annotated. The Illumina assembly provided the best reference to align an independent set of RNA-Seq data as ∼84% of reads mapped to single genes in the transcriptome. Contigs constructed solely from 454 data may impose problems for RNA-Seq as our 454 transcriptome revealed a high number of indels and many ambiguously mapped reads. Correcting the 454 transcriptome with Illumina reads was an effective strategy to deal with indel and frameshift errors inherent to the 454 transcriptome, but at the cost of transcriptome coverage. In the absence of a reference genome, we find that Illumina reads alone produced a high quality transcriptome appropriate for RNA-Seq gene expression analyses.

Antiviral activity of the inducible humoral immunity and its suppression by eleven BEN family members encoded in Cotesia plutellae bracovirus

Upon parasitization by some endoparasitoids, polydnaviruses (PDVs) play a crucial role in inducing host immunosuppression. This study reports a novel immunosuppressive activity against humoral immune responses by BEN family genes encoded in Cotesia plutellae bracovirus (CpBV). A total of 11 BEN family members are encoded in 10 different CpBV DNA segments. When the CpBV segments were individually injected, specific BEN genes were expressed and suppressed the expression of antimicrobial peptide (AMP) and prophenoloxidase genes following bacterial challenge. The suppressive activities of the BEN genes were reversed by injection of the double-stranded RNA (dsRNA) specific to each BEN gene. The suppression of the AMP gene expressions by the BEN genes was also confirmed using an inhibition zone assay against Gram-positive and Gram-negative bacterial growth. The significance of the suppressive activity of BEN genes against humoral immune responses was analyzed in terms of suppression of antiviral activity by the host humoral immunity. When CpBV was incubated with the plasma obtained from the larvae challenged with bacteria, the immunized plasma severely impaired the expression activity of the viral genes. However, an expression of BEN gene significantly rescued the viral gene expression by suppressing humoral immune response. These results suggest that BEN family genes of CpBV play a crucial role in defending the antiviral response of the parasitized Plutella xylostella by inhibiting humoral immune responses.

Extensive transcription analysis of the Hyposoter didymator Ichnovirus genome in permissive and non-permissive lepidopteran host species

Ichnoviruses are large dsDNA viruses that belong to the Polydnaviridae family. They are specifically associated with endoparasitic wasps of the family Ichneumonidae and essential for host parasitization by these wasps. We sequenced the Hyposoter didymator Ichnovirus (HdIV) encapsidated genome for further analysis of the transcription pattern of the entire set of HdIV genes following the parasitization of four different lepidopteran host species. The HdIV genome was found to consist of at least 50 circular dsDNA molecules, carrying 135 genes, 98 of which formed 18 gene families. The HdIV genome had general features typical of Ichnovirus (IV) genomes and closely resembled that of the IV carried by Hyposoter fugitivus. Subsequent transcriptomic analysis with Illumina technology during the course of Spodoptera frugiperda parasitization led to the identification of a small subset of less than 30 genes with high RPKM values in permissive hosts, consisting with these genes encoding crucial virulence proteins. Comparisons of HdIV expression profiles between host species revealed differences in transcript levels for given HdIV genes between two permissive hosts, S. frugiperda and Pseudoplusia includens. However, we found no evident intrafamily gene-specific transcription pattern consistent with the presence of multigenic families within IV genomes reflecting an ability of the wasps concerned to exploit different host species. Interestingly, in two non-permissive hosts, Mamestra brassiccae and Anticarsia gemmatalis (most of the parasitoid eggs were eliminated by the host cellular immune response), HdIV genes were generally less strongly transcribed than in permissive hosts. This suggests that successful parasitism is dependent on the expression of given HdIV genes exceeding a particular threshold value. These results raise questions about the mecanisms involved in regulating IV gene expression according to the nature of the lepidopteran host species encountered.

Parasitization by Scleroderma guani influences protein expression in Tenebrio molitor pupae

Ectoparasitoid wasps deposit their eggs onto the surface and inject venom into their hosts. Venoms are chemically complex and they exert substantial impact on hosts, including permanent or temporary paralysis and developmental arrest. These visible venom effects are due to changes in expression of genes encoding physiologically relevant proteins. While the influence of parasitization on gene expression in several lepidopterans has been reported, the molecular details of parasitoid/beetle relationships remain mostly unknown. This shortcoming led us to pose the hypothesis that envenomation by the ectoparasitic ant-like bethylid wasp Scleroderma guani leads to changes in protein expression in the yellow mealworm beetle Tenebrio molitor. We tested our hypothesis by comparing the proteomes of non-parasitized and parasitized host pupae using iTRAQ-based proteomics. We identified 41 proteins that were differentially expressed (32↑- and 9↓-regulated) in parasitized pupae. We assigned these proteins to functional categories, including immunity, stress and detoxification, energy metabolism, development, cytoskeleton, signaling and others. We recorded parallel changes in mRNA levels and protein abundance in 14 selected proteins following parasitization. Our findings support our hypothesis by documenting changes in protein expression in parasitized hosts.