Parasitization by Cotesia chilonis influences gene expression in fatbody and hemocytes of Chilo suppressalis

Immune system modulation & virus transmission during parasitism identified by multi-species transcriptomics of a declining insect biocontrol system

BACKGROUND

The Argentine stem weevil (ASW, Listronotus bonariensis) is a significant pasture pest in Aotearoa New Zealand, primarily controlled by the parasitoid biocontrol agent Microctonus hyperodae. Despite providing effective control of ASW soon after release, M. hyperodae parasitism rates have since declined significantly, with ASW hypothesised to have evolved resistance to its biocontrol agent. While the parasitism arsenal of M. hyperodae has previously been investigated, revealing many venom components and an exogenous novel DNA virus Microctonus hyperodae filamentous virus (MhFV), the effects of said arsenal on gene expression in ASW during parasitism have not been examined. In this study, we performed a multi-species transcriptomic analysis to investigate the biology of ASW parasitism by M. hyperodae, as well as the decline in efficacy of this biocontrol system.

RESULTS

The transcriptomic response of ASW to parasitism by M. hyperodae involves modulation of the weevil's innate immune system, flight muscle components, and lipid and glucose metabolism. The multispecies approach also revealed continued expression of venom components in parasitised ASW, as well as the transmission of MhFV to weevils during parasitism and some interrupted parasitism attempts. Transcriptomics did not detect a clear indication of parasitoid avoidance or other mechanisms to explain biocontrol decline.

CONCLUSIONS

This study has expanded our understanding of interactions between M. hyperodae and ASW in a biocontrol system of critical importance to Aotearoa-New Zealand's agricultural economy. Transmission of MhFV to ASW during successful and interrupted parasitism attempts may link to a premature mortality phenomenon in ASW, hypothesised to be a result of a toxin-antitoxin system. Further research into MhFV and its potential role in ASW premature mortality is required to explore whether manipulation of this viral infection has the potential to increase biocontrol efficacy in future.

Multi-Omic Identification of Venom Proteins Collected from Artificial Hosts of a Parasitoid Wasp

Habrobracon hebetor is a parasitoid wasp capable of infesting many lepidopteran larvae. It uses venom proteins to immobilize host larvae and prevent host larval development, thus playing an important role in the biocontrol of lepidopteran pests. To identify and characterize its venom proteins, we developed a novel venom collection method using an artificial host (ACV), i.e., encapsulated amino acid solution in paraffin membrane, allowing parasitoid wasps to inject venom. We performed protein full mass spectrometry analysis of putative venom proteins collected from ACV and venom reservoirs (VRs) (control). To verify the accuracy of proteomic data, we also collected venom glands (VGs), Dufour's glands (DGs) and ovaries (OVs), and performed transcriptome analysis. In this paper, we identified 204 proteins in ACV via proteomic analysis; compared ACV putative venom proteins with those identified in VG, VR, and DG via proteome and transcriptome approaches; and verified a set of them using quantitative real-time polymerase chain reaction. Finally, 201 ACV proteins were identified as potential venom proteins. In addition, we screened 152 and 148 putative venom proteins identified in the VG transcriptome and the VR proteome against those in ACV, and found only 26 and 25 putative venom proteins, respectively, were overlapped with those in ACV. Altogether, our data suggest proteome analysis of ACV in combination with proteome-transcriptome analysis of other organs/tissues will provide the most comprehensive identification of true venom proteins in parasitoid wasps.

An endoparasitoid uses its egg surface proteins to regulate its host immune response

With proteomic analysis, we identified 379 egg surface proteins from an endoparasitoid, Cotesia chilonis. Proteins containing conserved enzymatic domains constitute a large proportion of egg surface components. Some proteins, such as superoxidase dismutase, homolog of C. rubecula 32-kDa protein, and immunoevasive protein-2A, are classical parasitism factors that have known functions in host immunity regulation. Melanization assays revealed that a novel egg surface protein, C. chilonis egg surface serpin domain-containing protein had the same function as a C. chilonis venom serpin, as both suppressed host melanization in a dose-dependent manner. C. chilonis egg surface serpin domain-containing protein is mainly transcribed in C. chilonis oocytes with follicular cells, and it is located on both the anterior and posterior sides of the mature egg surface. Additionally, we used LC-MS/MS to identify 586 binding proteins sourced from C. suppressalis plasma located on the eggshell surface of C. chilonis, which included some immunity-related proteins. These results not only indicate that C. chilonis uses its egg surface proteins to reduce the immune response of its host but also imply that endoparasitoid egg surface proteins might be a new parasitism factor involved in host immune regulation.

CLP gene family, a new gene family of Cotesia vestalis bracovirus inhibits melanization of Plutella xylostella hemolymph

Polydnaviruses (PDVs) are obligatory symbionts of parasitoid wasps and play an important role in suppressing host immune defenses. Although PDV genes that inhibit host melanization are known in Microplitis bracovirus, the functional homologs in Cotesia bracoviruses remain unknown. Here, we find that Cotesia vestalis bracovirus (CvBV) can inhibit hemolymph melanization of its host, Plutella xylostella larvae, during the early stages of parasitization, and that overexpression of highly expressed CvBV genes reduced host phenoloxidase activity. Furthermore, CvBV-7-1 in particular reduced host phenoloxidase activity within 12 h, and the injection of anti-CvBV-7-1 antibody increased the melanization of parasitized host larvae. Further analyses showed that CvBV-7-1 and three homologs from other Cotesia bracoviruses possessed a C-terminal leucine/isoleucine-rich region and had a similar function in inhibiting melanization. Therefore, a new family of bracovirus genes was proposed and named as C-terminal Leucine/isoleucine-rich Protein (CLP). Ectopic expression of CvBV-7-1 in Drosophila hemocytes increased susceptibility to bacterial repression of melanization and reduced the melanotic encapsulation of parasitized D. melanogaster by the parasitoid Leptopilina boulardi. The formation rate of wasp pupae and the eclosion rate of C. vestalis were affected when the function of CvBV-7-1 was blocked. Our findings suggest that CLP genes from Cotesia bracoviruses encoded proteins that contain a C-terminal leucine/isoleucine-rich region and function as melanization inhibitors during the early stage of parasitization, which is important for successful parasitization.

Comparison of morphology, development and expression patterns of hsf and hsp11.0 of Cotesia chilonis under normal and high temperature

Cotesia chilonis (Munakata) is the dominant parasitic wasp of the rice pest, Chilo suppressalis (Walker), and is a valuable parasitic wasp for the prevention and control of C. suppressalis. In this study, developmental indicators and expression of Cchsp11.0 (heat shock protein 11.0) and Cchsf (heat shock factor) were compared for C. chilonis at 27 °C and 36 °C. Developmental duration, morphology, emergence rate, and number of C. chilonis offspring were shortened at 36 °C while the ratio of females to males increased. Cchsp11.0 and Cchsf were highly expressed in the 1^st instar stage at 36 °C, and Cchsp11.0 expression gradually decreased as C. chilonis matured; Cchsf expression was not correlated with Cchsp11.0 expression. Compared with 27 °C, the expression pattern of Cchsp11.0 and Cchsf was also not consistent, and Cchsp11.0 expression increased significantly at the adult stage. In conclusion, mildly high temperatures impact growth, development and reproduction of C. chilonis and stimulate the expression of Cchsp11.0 and Cchsf, and Cchsp11.0 and Cchsf play different roles in different developmental stages of C. chilonis at normal and high temperature.

Towards unravelling Wolbachia global exchange: a contribution from the Bicyclus and Mylothris butterflies in the Afrotropics

Background

Phylogenetically closely related strains of maternally inherited endosymbiotic bacteria are often found in phylogenetically divergent, and geographically distant insect host species. The interspecies transfer of the symbiont Wolbachia has been thought to have occurred repeatedly, facilitating its observed global pandemic. Few ecological interactions have been proposed as potential routes for the horizontal transfer of Wolbachia within natural insect communities. These routes are however likely to act only at the local scale, but how they may support the global distribution of some Wolbachia strains remains unclear.

Results

Here, we characterize the Wolbachia diversity in butterflies from the tropical forest regions of central Africa to discuss transfer at both local and global scales. We show that numerous species from both the Mylothris (family Pieridae) and Bicyclus (family Nymphalidae) butterfly genera are infected with similar Wolbachia strains, despite only minor interclade contacts across the life cycles of the species within their partially overlapping ecological niches. The phylogenetic distance and differences in resource use between these genera rule out the role of ancestry, hybridization, and shared host-plants in the interspecies transfer of the symbiont. Furthermore, we could not identify any shared ecological factors to explain the presence of the strains in other arthropod species from other habitats, or even ecoregions.

Conclusion

Only the systematic surveys of the Wolbachia strains from entire species communities may offer the material currently lacking for understanding how Wolbachia may transfer between highly different and unrelated hosts, as well as across environmental scales.

Supplementary information

Supplementary information accompanies this paper at 10.1186/s12866-020-02011-2.

Amino acid synthesis loss in parasitoid wasps and other hymenopterans

Insects utilize diverse food resources which can affect the evolution of their genomic repertoire, including leading to gene losses in different nutrient pathways. Here, we investigate gene loss in amino acid synthesis pathways, with special attention to hymenopterans and parasitoid wasps. Using comparative genomics, we find that synthesis capability for tryptophan, phenylalanine, tyrosine, and histidine was lost in holometabolous insects prior to hymenopteran divergence, while valine, leucine, and isoleucine were lost in the common ancestor of Hymenoptera. Subsequently, multiple loss events of lysine synthesis occurred independently in the Parasitoida and Aculeata. Experiments in the parasitoid Cotesia chilonis confirm that it has lost the ability to synthesize eight amino acids. Our findings provide insights into amino acid synthesis evolution, and specifically can be used to inform the design of parasitoid artificial diets for pest control.

Characterization of an inducible HSP70 gene in Chilo suppressalis and expression in response to environmental and biological stress

The highly conserved heat shock protein 70 (HSP70) contributes to survival at a cellular level and greatly enhances stress tolerance in many organisms. In this study, we isolate and characterize Cshsp702, which encodes an inducible form of HSP70 in the rice stem borer, Chilo suppressalis. Cshsp702 does not contain introns; the translational product is comprised of 629 amino acids with an isoelectric point of 5.69. Real-time quantitative PCR revealed that Cshsp702 was expressed at maximal levels in hemocytes and was minimally expressed in the midgut. Expression of Cshsp702 in response to a range of temperatures (-11 to 43 °C) indicated significant induction by extreme cold and hot temperatures, with maximum expression after 2 h at 42 °C. The induction of Cshsp702 in response to the endoparasite Cotesia chilonis was also studied; interestingly, Cshsp702 expression in C. suppressalis was significantly induced at 24 h and 5 days, which correspond to predicted times of C. chilonis feeding and growth, respectively. The potential induction of Cshsp702 as an inflammatory response due to parasitic stress is discussed. In conclusion, Cshsp702 is induced in response to both environmental and biotic stress and plays an important role in the physiological adaptation of C. suppressalis.

Selection of valid reference genes for quantitative real-time PCR in Cotesia chilonis (Hymenoptera: Braconidae) exposed to different temperatures

In quantitative real-time PCR (qRT-PCR), data are normalized using reference genes, which helps to control for internal differences and reduce error among samples. In this study, the expression profiles of eight candidate housekeeping genes, 18S ribosomal (18S rRNA), elongation factor (EF1), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), ribosomal protein L10 (RPL10), ribosomal protein L17 (RPL17), histone 3 (H3), arginine kinase (AK), amd β-Actin (ACTB), were evaluated in the parasitic wasp Cotesia chilonis in response to different temperatures. Specifically, the performance and stabilities of these genes were compared in adult wasps maintained in a growth condition at 27°C (normal storage conditions) and in adults obtained from pupae refrigerated at 4°C for five days (cold storage conditions). Data were analyzed using the ΔCt method, BestKeeper, NormFinder, and geNorm. The optimal numbers and stabilities of reference genes varied between the two temperature treatments (27°C and 4°C). In samples stored at normal developmental temperature (27°C), the requirement for normalization in response to low temperature exposures was three genes (18S, H3, AK), whereas normalization in response to high temperature exposures required only two reference genes (GAPDH, ACTB). In samples stored at cold temperature (4°C), for low temperature exposures two reference genes (RPL17, RPL10) were required for standardization, while following high temperature exposures three reference genes (18S, H3, ACTB) were needed. This study strengthens understanding of the selection of reference genes before qRT-PCR analysis in C. chilonis. The reference genes identified here will facilitate further investigations of the biological characteristics of this important parasitoid.

An Ovarian Protein Involved in Passive Avoidance of an Endoparasitoid To Evade Its Host Immune Response

Through a combination of transcriptomic and proteomic analyses, we identified 817 secreted ovarian proteins from an endoparasitoid wasp, Cotesia chilonis, of which five proteins are probably involved in passive evasion. The results of an encapsulation assay revealed that one of these passive evasion-associated proteins (Crp32B), a homologue of a 32-kDa protein (Crp32) from C. rubecula, could protect resin beads from being encapsulated by host hemocytes in a dose-dependent manner. Crp32B is transcribed in ovarian cells, nurse cells, follicular cells, and oocytes, and the protein is located throughout the ovary and on the egg surface. Moreover, Crp32B has antigenic similarity to several host components. These results indicate that C. chilonis may use molecular mimicry as a mechanism to avoid host cellular immune response.

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.

Transferrin Family Genes in the Brown Planthopper, Nilaparvata lugens (Hemiptera: Delphacidae) in Response to Three Insecticides

Transferrins are involved in iron metabolism, immunity, xenobiotics tolerance, and development in eukaryotic organisms including insects. However, little is known about the relationship between transferrins and insecticide toxicology and resistance. Three transferrin family genes, NlTsf1, NlTsf2, and NlTsf3, of the brown planthopper, Nilaparvata lugens (Stål) (Hemiptera: Delphacidae)a major insect pest of rice field in Asia, had been identified and characterized in this study. Quantitative polymerase chain reaction results demonstrated that NlTsf1 was significantly higher than the other two genes in different tissues. All of them were expressed at higher levels in abdomen and head than in antenna, leg, stylet, and thorax. Compared with the control, the expression of three N. lugens transferrin family genes decreased dramatically 24 h after treatment with buprofezin, pymetrozine and imidacloprid.

A Horizontally Transferred Autonomous Helitron Became a Full Polydnavirus Segment in Cotesia vestalis

Bracoviruses associate symbiotically with thousands of parasitoid wasp species in the family Braconidae, working as virulence gene vectors, and allowing the development of wasp larvae within hosts. These viruses are composed of multiple DNA circles that are packaged into infective particles, and injected together with wasp's eggs during parasitization. One of the viral segments of Cotesia vestalis bracovirus contains a gene that has been previously described as a helicase of unknown origin. Here, we demonstrate that this gene is a Rep/Helicase from an intact Helitron transposable element that covers the viral segment almost entirely. We also provide evidence that this element underwent at least two horizontal transfers, which appear to have occurred consecutively: first from a Drosophila host ancestor to the genome of the parasitoid wasp C. vestalis and its bracovirus, and then from C. vestalis to a lepidopteran host (Bombyx mori). Our results reinforce the idea of parasitoid wasps as frequent agents of horizontal transfers in eukaryotes. Additionally, this Helitron-bracovirus segment is the first example of a transposable element that effectively became a whole viral circle.

ACE: an efficient and sensitive tool to detect insecticide resistance-associated mutations in insect acetylcholinesterase from RNA-Seq data

Background

Insecticide resistance is a substantial problem in controlling agricultural and medical pests. Detecting target site mutations is crucial to manage insecticide resistance. Though PCR-based methods have been widely used in this field, they are time-consuming and inefficient, and typically have a high false positive rate. Acetylcholinesterases (Ace) is the neural target of the widely used organophosphate (OP) and carbamate insecticides. However, there is not any software available to detect insecticide resistance associated mutations in RNA-Seq data at present.

Results

A computational pipeline ACE was developed to detect resistance mutations of ace in insect RNA-Seq data. Known ace resistance mutations were collected and used as a reference. We constructed a Web server for ACE, and the standalone software in both Linux and Windows versions is available for download. ACE was used to analyse 971 RNA-Seq data from 136 studies in 7 insect pests. The mutation frequency of each RNA-Seq dataset was calculated. The results indicated that the resistance frequency was 30%–44% in an eastern Ugandan Anopheles population, thus suggesting this resistance-conferring mutation has reached high frequency in these mosquitoes in Uganda. Analyses of RNA-Seq data from the diamondback moth Plutella xylostella indicated that the G227A mutation was positively related with resistance levels to organophosphate or carbamate insecticides. The wasp Nasonia vitripennis had a low frequency of resistant reads (<5%), but the agricultural pests Chilo suppressalis and Bemisia tabaci had a high resistance frequency. All ace reads in the 30 B. tabaci RNA-Seq data were resistant reads, suggesting that insecticide resistance has spread to very high frequency in B. tabaci.

Conclusions

To the best of our knowledge, the ACE pipeline is the first tool to detect resistance mutations from RNA-Seq data, and it facilitates the full utilization of large-scale genetic data obtained by using next-generation sequencing.

Electronic supplementary material

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

Molecular characterization and expression profiles of nicotinic acetylcholine receptors in the rice striped stem borer, Chilo suppressalis (Lepidoptera: Crambidae)

Nicotinic acetylcholine receptors (nAChRs) are members of the cys-loop ligand-gated ion channel (cysLGIC) superfamily, mediating fast synaptic cholinergic transmission in the central nervous system in insects. Insect nAChRs are the molecular targets of economically important insecticides, such as neonicotinoids and spinosad. Identification and characterization of the nAChR gene family in the rice striped stem borer, Chilo suppressalis, could provide beneficial information about this important receptor gene family and contribute to the investigation of the molecular modes of insecticide action and resistance for current and future chemical control strategies. We searched our C. suppressalis transcriptome database using Bombyx mori nAChR sequences in local BLAST searches and obtained the putative nAChR subunit complementary DNAs (cDNAs) via reverse transcription polymerase chain reaction (RT-PCR) and rapid amplification of cDNA ends methods. Similar to B. mori, C. suppressalis possesses 12 nAChR subunits, including nine α-type and three β-type subunits. Quantitative RT-PCR analysis revealed the expression profiles of the nAChR subunits in various tissues, including the brain, subesophageal ganglion, thoracic ganglion, abdominal ganglion, hemocytes, fat body, foregut, midgut, hindgut and Malpighian tubules. Developmental expression analyses showed clear differential expression of nAChR subunits throughout the C. suppressalis life cycle. The identification of nAChR subunits in this study will provide a foundation for investigating the diverse roles played by nAChRs in C. suppressalis and for exploring specific target sites for chemicals that control agricultural pests while sparing beneficial species.

Protein Discovery: Combined Transcriptomic and Proteomic Analyses of Venom from the Endoparasitoid Cotesia chilonis (Hymenoptera: Braconidae)

Many species of endoparasitoid wasps provide biological control services in agroecosystems. Although there is a great deal of information on the ecology and physiology of host/parasitoid interactions, relatively little is known about the protein composition of venom and how specific venom proteins influence physiological systems within host insects. This is a crucial gap in our knowledge because venom proteins act in modulating host physiology in ways that favor parasitoid development. Here, we identified 37 possible venom proteins from the polydnavirus-carrying endoparasitoid Cotesia chilonis by combining transcriptomic and proteomic analyses. The most abundant proteins were hydrolases, such as proteases, peptidases, esterases, glycosyl hydrolase, and endonucleases. Some components are classical parasitoid venom proteins with known functions, including extracellular superoxide dismutase 3, serine protease inhibitor and calreticulin. The venom contains novel proteins, not recorded from any other parasitoid species, including tolloid-like proteins, chitooligosaccharidolytic β-N-acetylglucosaminidase, FK506-binding protein 14, corticotropin-releasing factor-binding protein and vascular endothelial growth factor receptor 2. These new data generate hypotheses and provide a platform for functional analysis of venom components.

Pharmacological characterization of dopamine receptors in the rice striped stem borer, Chilo suppressalis

Dopamine is an important neurotransmitter and neuromodulator in both vertebrates and invertebrates and is the most abundant monoamine present in the central nervous system of insects. A complement of functionally distinct dopamine receptors mediate the signal transduction of dopamine by modifying intracellular Ca²⁺ and cAMP levels. In the present study, we pharmacologically characterized three types of dopamine receptors, CsDOP1, CsDOP2 and CsDOP3, from the rice striped stem borer, Chilo suppressalis. All three receptors show considerable sequence identity with orthologous dopamine receptors. The phylogenetic analysis also clusters the receptors within their respective groups. Transcript levels of CsDOP1, CsDOP2 and CsDOP3 were all expressed at high levels in the central nervous system, indicating their important roles in neural processes. After heterologous expression in HEK 293 cells, CsDOP1, CsDOP2 and CsDOP3 were dose-dependently activated by dopamine and synthetic dopamine receptor agonists. They can also be blocked by different series of antagonists. This study offers important information on three dopamine receptors from C. suppressalis that will provide the basis for forthcoming studies investigating their roles in behaviors and physiology, and facilitate the development of new insecticides for pest control.

Differential responses of Helicoverpa armigera C-type immunlectin genes to the endoparasitoid Campoletis chlorideae

The C-type lectins mediate nonself recognition in insects. The previous studies focused on host immunlectin response to bacterial infection; however, the molecular basis of immunlectin reactions to endoparasitoids has not been elucidated. The present study investigated the effect of parasitization by Campoletis chlorideae on hemagglutination activity (HA; defined as the ability of lectin to agglutinate erythrocytes or other cells), and transcriptional expression of C-type immunlectin genes in the larval host, Helicoverpa armigera. Parasitization induced four- to eightfold higher HA in the parasitized larvae, compared to nonparasitized larvae at days 2 and 6 postparasitization (PP), however inhibited HA at other days PP. Eight C-type lectins were differentially expressed in different host developmental stages, from feeding to wandering stage. The mRNA levels of HaCTL1, HaCTL3, HaCTL4, and HaCTL5 were upregulated and HaCTL2 and HaCTL7 were downregulated. Tissue analysis showed that HaCTLs were mainly expressed in fat body or hemocytes, while HaCTL5 was highly expressed in testes. The effects of parasitization on the lectin expression patterns differed. Lectins except HaCTL6 or HaCTL5 were significantly down- or upregulated in parasitized larvae at day 4 or 6 PP compared with that of nonparasitized larvae. We infer from our results that C-type immunlectins are involved in host-parasitoid interactions, and parasitization alter host immunlectin levels both in inhibiting and promoting host immune defenses to endoparasitoids. These immunlectin genes indicated an altered physiological status of the host insect, depending on developmental stage, tissue, and parasitization.

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.

Effects of the endoparasitoid Cotesia chilonis (Hymenoptera: Braconidae) parasitism, venom, and calyx fluid on cellular and humoral immunity of its host Chilo suppressalis (Lepidoptera: Crambidae) larvae

The larval endoparasitoid Cotesia chilonis injects venom and bracoviruses into its host Chilo suppressalis during oviposition. Here we study the effects of the polydnavirus (PDV)-carrying endoparasitoid C. chilonis (Hymenoptera: Braconidae) parasitism, venom and calyx fluid on host cellular and humoral immunity, specifically hemocyte composition, cellular spreading, encapsulation and melanization. Total hemocyte counts (THCs) were higher in parasitized larvae than in unparasitized larvae in the late stages following parasitization. While both plasmatocyte and granulocyte fractions and hemocyte mortality did not differ between parasitized and unparasitized hosts, in vitro spreading behavior of hemocytes was inhibited significantly by parasitism throughout the course of parasitoid development. C. chilonis parasitism suppressed the encapsulation response and melanization in the early stages. Venom alone did not alter cellular immune responses, including effects on THCs, mortality, hemocyte composition, cell spreading and encapsulation, but venom did inhibit humoral immunity by reducing melanization within 6h after injection. In contrast to venom, calyx fluid had a significant effect on cell spreading, encapsulation and melanization from 6h after injection. Dose-response injection studies indicated the effects of venom and calyx fluid synergized, showing a stronger and more persistent reduction in immune system responses than the effect of either injected alone.

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.

De novo assembly and characterization of central nervous system transcriptome reveals neurotransmitter signaling systems in the rice striped stem borer, Chilo suppressalis

Background

Neurotransmitter signaling systems play crucial roles in multiple physiological and behavioral processes in insects. Genome wide analyses of de novo transcriptome sequencing and gene specific expression profiling provide rich resources for studying neurotransmitter signaling pathways. The rice striped stem borer, Chilo suppressalis is a destructive rice pest in China and other Asian countries. The characterization of genes involved in neurotransmitter biosynthesis and transport could identify potential targets for disruption of the neurochemical communication and for crop protection.

Results

Here we report de novo sequencing of the C. suppressalis central nervous system transcriptome, identification and expression profiles of genes putatively involved in neurotransmitter biosynthesis, packaging, and recycling/degradation. A total of 54,411 unigenes were obtained from the transcriptome analysis. Among these unigenes, we have identified 32 unigenes (31 are full length genes), which encode 21 enzymes and 11 transporters putatively associated with biogenic aminergic signaling, acetylcholinergic signaling, glutamatergic signaling and GABAergic signaling. RT-PCR and qRT-PCR results indicated that 12 enzymes were highly expressed in the central nervous system and all the transporters were expressed at significantly high levels in the central nervous system. In addition, the transcript abundances of enzymes and transporters in the central nervous system were validated by qRT-PCR. The high expression levels of these genes suggest their important roles in the central nervous system.

Conclusions

Our study identified genes potentially involved in neurotransmitter biosynthesis and transport in C. suppressalis and these genes could serve as targets to interfere with neurotransmitter production. This study presents an opportunity for the development of specific and environmentally safe insecticides for pest control.

Electronic supplementary material

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

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.

Transcriptome analysis of an endoparasitoid wasp Cotesia chilonis (Hymenoptera: Braconidae) reveals genes involved in successful parasitism

For successful parasitization, parasitiods usually depend on the chemosensory cues for the selection of hosts, as well as a variety of virulence factors introduced into their hosts to overcome host immunity and prevent rejection of progeny development. In bracovirus-carrying wasps, the symbiotic polydnaviruses act in manipulating development and immunity of hosts. The endoparasitoid Cotesia chilonis carrying bracovirus as a key host immunosuppressive factor is a superior endoparasitoid of rice stem borer, Chilo suppressalis. So far, genomic information for C. chilonis is not available and transcriptomic data may provide valuable resources for global studying on physiological processes of C. chilonis, including chemosensation and parasitism at molecular level. Here, we performed RNA-seq to characterize the transcriptome of C. chilonis adults. We obtained 27,717,892 reads, assembled into 38,318 unigenes with a mean size of 690 bp. Approximately, 62.1% of the unigenes were annotated using NCBI databases. A large number of chemoreception-related genes encoding proteins including odorant receptors, gustatory receptors, odorant-binding proteins, chemosensory proteins, transient receptor potential ion channels, and sensory neuron membrane proteins were identified in silico. Totally, 72 transcripts possessing high identities with the bracovirus-related genes were identified. We investigated the mRNA expression levels of several transcripts at different developmental stages (including egg, larva, pupae, and adult) by quantitative real-time PCR analysis. The results revealed that some genes had adult-specific expression, indicating their potential significance for mating and parasitism. Overall, these results provide comprehensive insights into transcriptomic data of a polydnavirus-carrying parasitoid of a rice pest.

Functional annotation of Cotesia congregata bracovirus: identification of viral genes expressed in parasitized host immune tissues

Bracoviruses (BVs) from the Polydnaviridae family are symbiotic viruses used as biological weapons by parasitoid wasps to manipulate lepidopteran host physiology and induce parasitism success. BV particles are produced by wasp ovaries and injected along with the eggs into the caterpillar host body, where viral gene expression is necessary for wasp development. Recent sequencing of the proviral genome of Cotesia congregata BV (CcBV) identified 222 predicted virulence genes present on 35 proviral segments integrated into the wasp genome. To date, the expressions of only a few selected candidate virulence genes have been studied in the caterpillar host, and we lacked a global vision of viral gene expression. In this study, a large-scale transcriptomic analysis by 454 sequencing of two immune tissues (fat body and hemocytes) of parasitized Manduca sexta caterpillar hosts allowed the detection of expression of 88 CcBV genes expressed 24 h after the onset of parasitism. We linked the expression profiles of these genes to several factors, showing that different regulatory mechanisms control viral gene expression in the host. These factors include the presence of signal peptides in encoded proteins, diversification of promoter regions, and, more surprisingly, gene position on the proviral genome. Indeed, most genes for which expression could be detected are localized in particular proviral regions globally producing higher numbers of circles. Moreover, this polydnavirus (PDV) transcriptomic analysis also reveals that a majority of CcBV genes possess at least one intron and an arthropod transcription start site, consistent with an insect origin of these virulence genes.

IMPORTANCE

Bracoviruses (BVs) are symbiotic polydnaviruses used by parasitoid wasps to manipulate lepidopteran host physiology, ensuring wasp offspring survival. To date, the expressions of only a few selected candidate BV virulence genes have been studied in caterpillar hosts. We performed a large-scale analysis of BV gene expression in two immune tissues of Manduca sexta caterpillars parasitized by Cotesia congregata wasps. Genes for which expression could be detected corresponded to genes localized in particular regions of the viral genome globally producing higher numbers of circles. Our study thus brings an original global vision of viral gene expression and paves the way to the determination of the regulatory mechanisms enabling the expression of BV genes in targeted organisms, such as major insect pests. In addition, we identify sequence features suggesting that most BV virulence genes were acquired from insect genomes.