Nasonia vitripennis venom causes targeted gene expression changes in its fly host

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.

Comprehensive isoform-level analysis reveals the contribution of alternative isoforms to venom evolution and repertoire diversity

Animal venom systems have emerged as valuable models for investigating how novel polygenic phenotypes may arise from gene evolution by varying molecular mechanisms. However, a significant portion of venom genes produce alternative mRNA isoforms that have not been extensively characterized, hindering a comprehensive understanding of venom biology. In this study, we present a full-length isoform-level profiling workflow integrating multiple RNA sequencing technologies, allowing us to reconstruct a high-resolution transcriptome landscape of venom genes in the parasitoid wasp Pteromalus puparum Our findings demonstrate that more than half of the venom genes generate multiple isoforms within the venom gland. Through mass spectrometry analysis, we confirm that alternative splicing contributes to the diversity of venom proteins, acting as a mechanism for expanding the venom repertoire. Notably, we identified seven venom genes that exhibit distinct isoform usages between the venom gland and other tissues. Furthermore, evolutionary analyses of venom serpin3 and orcokinin further reveal that the co-option of an ancient isoform and a newly evolved isoform, respectively, contributes to venom recruitment, providing valuable insights into the genetic mechanisms driving venom evolution in parasitoid wasps. Overall, our study presents a comprehensive investigation of venom genes at the isoform level, significantly advancing our understanding of alternative isoforms in venom diversity and evolution and setting the stage for further in-depth research on venoms.

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.

Host Transcriptome Analysis of Spodoptera frugiperda Larvae Parasitized by Microplitis manilae

It has been extensively found that parasitoids manipulate host physiology to benefit the survival and development of their offspring. However, the underlying regulatory mechanisms have not received much attention. To reveal the effects of parasitization of the larval solitary endoparasitoid Microplitis manilae (Hymenoptera: Braconidae) on host Spodoptera frugiperda (Lepidoptera: Noctuidae), one of the most destructive agricultural pests in China, deep-sequencing-based transcriptome analysis was conducted to compare the host gene expression levels after 2 h, 24 h, and 48 h parasitization. A total of 1861, 962, and 108 differentially expressed genes (DEGs) were obtained from the S. frugiperda larvae at 2 h, 24 h, and 48 h post-parasitization, respectively, compared with unparasitized controls. The changes in host gene expressions were most likely caused by the injection of wasp parasitic factors, including PDVs, that were injected along with the eggs during oviposition. Based on the functional annotations in GO and KEGG databases, we revealed that most DEGs were implicated in host metabolism and immunity. Further analysis of the common DEGs in three comparisons between the unparasitized and parasitized groups identified four genes, including one unknown and three prophenoloxidase (PPO) genes. Moreover, 46 and 7 common DEGs involved in host metabolism and immunity were identified at two or three time points after parasitization, respectively. Among these, most DEGs showed increased expressions at 2 h post-wasp parasitization while exhibiting significantly decreased expression levels at 24 h post-parasitization, demonstrating the expression regulations of M. manilae parasitization on host metabolism and immune-related genes. Further qPCR verification in 20 randomly selected DEGs confirmed the accuracy and reproducibility of the gene expression profiles generated from RNA-seq. This study reveals the molecular regulatory network about how host insects respond to wasp parasitism, laying a solid foundation for revealing the physiological manipulation of wasp parasitization on host insects, which facilitates the development of biological control practices for parasitoids.

Genome of the parasitoid wasp Cotesia chilonis sheds light on amino acid resource exploitation

Background

A fundamental feature of parasitism is the nutritional exploitation of host organisms by their parasites. Parasitoid wasps lay eggs on arthropod hosts, exploiting them for nutrition to support larval development by using diverse effectors aimed at regulating host metabolism. However, the genetic components and molecular mechanisms at the basis of such exploitation, especially the utilization of host amino acid resources, remain largely unknown. To address this question, here, we present a chromosome-level genome assembly of the parasitoid wasp Cotesia chilonis and reconstruct its amino acid biosynthetic pathway.

Results

Analyses of the amino acid synthetic pathway indicate that C. chilonis lost the ability to synthesize ten amino acids, which was confirmed by feeding experiments with amino acid-depleted media. Of the ten pathways, nine are known to have been lost in the common ancestor of animals. We find that the ability to synthesize arginine was also lost in C. chilonis because of the absence of two key genes in the arginine synthesis pathway. Further analyses of the genomes of 72 arthropods species show that the loss of arginine synthesis is common in arthropods. Metabolomic analyses by UPLC-MS/MS reveal that the temporal concentrations of arginine, serine, tyrosine, and alanine are significantly higher in host (Chilo suppressalis) hemolymph at 3 days after parasitism, whereas the temporal levels of 5-hydroxylysine, glutamic acid, methionine, and lysine are significantly lower. We sequence the transcriptomes of a parasitized host and non-parasitized control. Differential gene expression analyses using these transcriptomes indicate that parasitoid wasps inhibit amino acid utilization and activate protein degradation in the host, likely resulting in the increase of amino acid content in host hemolymph.

Conclusions

We sequenced the genome of a parasitoid wasp, C. chilonis, and revealed the features of trait loss in amino acid biosynthesis. Our work provides new insights into amino acid exploitation by parasitoid wasps, and this knowledge can specifically be used to design parasitoid artificial diets that potentially benefit mass rearing of parasitoids for pest control.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01313-3.

A venom protein of ectoparasitoid Pachycrepoideus vindemiae, PvG6PDH, contributes to parasitism by inhibiting host glucose-6-phosphate metabolism

To achieve successful development, female parasitoids, while laying eggs, introduce various virulence factors, mainly venoms, into host insects to manipulate their physiology. Although numerous studies have been conducted to characterize the components of venoms that regulate host immune responses, few systematic investigations have been conducted on the roles of venom proteins in host metabolic regulation. In this investigation, we characterized a novel venom protein in Pachycrepoideus vindemiae called glucose-6-phosphate dehydrogenase (PvG6PDH) and showed it has a vital role in regulating host carbohydrate metabolism. PvG6PDH encodes 510 amino acids and features a signal peptide and two conserved "G6PDH" domains. Multiple sequence alignment showed it has high amino acid identity with G6PDH from other pteromalids, and quantitative polymerase chain reaction analysis and immunofluorescent staining demonstrated a significantly higher expression of PvG6PDH in the venom apparatus compared with the carcass. We report that PvG6PDH contributes to parasitism by inhibiting the glucose-6-phosphate (G6P) metabolism of host Drosophila melanogaster, as demonstrated by PvG6PDH injection and RNA interference analysis. Further tests revealed that the accumulation of host G6P was caused by the transcriptional inhibition of G6P-metabolism-related genes. These findings greatly contribute to our understanding of venom-mediated host metabolic regulation, further laying the foundation for the development of venom proteins as biological agents for pest control.

Males of the parasitoid wasp, Nasonia vitripennis, can identify which fly hosts contain females

The reproductive success of a male is limited by the number of females it can mate with. Thus, males deploy elaborate strategies to maximize access to females. In Nasonia, which are parasitoids of cyclorrhaphous flies, such reproductive strategies are thought to be restricted to competition among males for access to females in the natal patch. This study investigates whether additional strategies are present, especially the capability to identify which fly hosts contain adult females inside. Behavioural assays revealed that only one out of the four species, N. vitripennis, can distinguish which hosts specifically have adult female wasps, indicating a species-specific reproductive strategy. Results of gas chromatography-mass spectrometry analyses and behavioural data suggest that female-signature cuticular hydrocarbons (CHCs) are used as chemical cues, possibly emanating from within the host puparium. Further assays indicated that N. vitripennis males can also detect differences in the intensities of female-signature CHCs, giving them the capability to seek out hosts with maximum number of females. This study uncovers a previously unknown reproductive strategy in one of the most widely studied parasitoid wasps.

Tissue-specific gene expression shows a cynipid wasp repurposes oak host gene networks to create a complex and novel parasite-specific organ

Every organism on Earth depends on interactions with other organisms to survive. In each of these interactions, an organism must utilize the limited toolbox of genes and proteins it possesses to successfully manipulate or cooperate with another species, but it can also co-opt the genome machinery of its partner to expand its available tools. Insect-induced plant galls are an extreme example of this, wherein an insect hijacks the plant's genome to direct the initiation and development of galls consisting of plant tissue. However, previous transcriptomic studies have not evaluated individual tissues within a gall to determine the full extent to which a galling insect manipulates its host plant. Here we demonstrate that the cynipid wasp Dryocosmus quercuspalustris creates a complex parasite-specific organ from red oak tissue via massive changes in host gene expression. Our results show that the gall wasp is not merely modifying oak leaf tissue but creating extensive changes in gene expression between galled and ungalled tissue (differential expression in 28% of genes) and distinct gall tissue types (20% of genes). The outer gall tissue shows increases in various plant defence systems, which is consistent with its predicted functional role of protecting the wasp larva. The inner larval capsule shows suppression of large parts of the plant innate immune system and evidence for the wasp utilizing the plant's RNA interference mechanisms, which may be a potential mechanism for the wasp's control on gall growth.

Genome of the pincer wasp Gonatopus flavifemur reveals unique venom evolution and a dual adaptation to parasitism and predation

Background

Hymenoptera comprise extremely diverse insect species with extensive variation in their life histories. The Dryinidae, a family of solitary wasps of Hymenoptera, have evolved innovations that allow them to hunt using venom and a pair of chelae developed from the fore legs that can grasp prey. Dryinidae larvae are also parasitoids of Auchenorrhyncha, a group including common pests such as planthoppers and leafhoppers. Both of these traits make them effective and valuable for pest control, but little is yet known about the genetic basis of its dual adaptation to parasitism and predation.

Results

We sequenced and assembled a high-quality genome of the dryinid wasp Gonatopus flavifemur, which at 636.5 Mb is larger than most hymenopterans. The expansion of transposable elements, especially DNA transposons, is a major contributor to the genome size enlargement. Our genome-wide screens reveal a number of positively selected genes and rapidly evolving proteins involved in energy production and motor activity, which may contribute to the predatory adaptation of dryinid wasp. We further show that three female-biased, reproductive-associated yellow genes, in response to the prey feeding behavior, are significantly elevated in adult females, which may facilitate the egg production. Venom is a powerful weapon for dryinid wasp during parasitism and predation. We therefore analyze the transcriptomes of venom glands and describe specific expansions in venom Idgf-like genes and neprilysin-like genes. Furthermore, we find the LWS2-opsin gene is exclusively expressed in male G. flavifemur, which may contribute to partner searching and mating.

Conclusions

Our results provide new insights into the genome evolution, predatory adaptation, venom evolution, and sex-biased genes in G. flavifemur, and present genomic resources for future in-depth comparative analyses of hymenopterans that may benefit pest control.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01081-6.

Comparative transcriptome analysis reveals a potential mechanism for host nutritional manipulation after parasitization by Leptopilina boulardi

Parasitoids have been extensively found to manipulate nutrient amounts of their hosts to benefit their own development and survival, but the underlying mechanisms are largely unknown. Leptopilina boulardi (Hymenoptera: Figitidae) is a larval-pupal endoparasitoid wasp of Drosophila melanogaster whose survival relies on the nutrients provided by its Drosophila host. Here, we used RNA-seq to compare the gene expression levels of the host midgut at 24 h and 48 h post L. boulardi parasitization. We obtained 95 and 191 differentially expressed genes (DEGs) in the parasitized host midgut at 24 h and 48 h post L. boulardi parasitization, respectively. A KEGG analysis revealed that several metabolic pathways were significantly enriched in the upregulated DEGs, and these pathways included "starch and sucrose metabolism" and "galactose metabolism". A functional annotation analysis showed that four classes of genes involved in carbohydrate digestion process had increased expression levels in the midgut post L.boulardi parasitization than nonparasitized groups: glucosidase, mannosidase, chitinase and amylase. Genes involved in protein digestion process were also found among the DEGs, and most of these genes, which belonged to the metallopeptidase and serine-type endopeptidase families, were found at higher expression levels in the parasitized host midgut comparing with nonparasitized hosts. Moreover, some immune genes, particularly those involved in the Toll and Imd pathways, also exhibited high expression levels after L.boulardi parasitization. Our study provides large-scale transcriptome data and identifies sets of DEGs between parasitized and nonparasitized host midgut tissues at 24 h and 48 h post L. boulardi parasitization. These resources help improve our understanding of how parasitoid infection affects the nutrient components in the hosts.

The evolutionary dynamics of venom toxins made by insects and other animals

Animal venoms are recognised as unique biological systems in which to study molecular evolution. Venom use has evolved numerous times among the insects, and insects today use venom to capture prey, defend themselves from predators, or to subdue and modulate host responses during parasitism. However, little is known about most insect venom toxins or the mode and tempo by which they evolve. Here, I review the evolutionary dynamics of insect venom toxins, and argue that insects offer many opportunities to examine novel aspects of toxin evolution. The key questions addressed are: How do venomous animals evolve from non-venomous animals, and how does this path effect the composition and pharmacology of the venom? What genetic processes (gene duplication, co-option, neofunctionalisation) are most important in toxin evolution? What kinds of selection pressures are acting on toxin-encoding genes and their cognate targets in envenomated animals? The emerging evidence highlights that venom composition and pharmacology adapts quickly in response to changing selection pressures resulting from new ecological interactions, and that such evolution occurs through a stunning variety of genetic mechanisms. Insects offer many opportunities to investigate the evolutionary dynamics of venom toxins due to their evolutionary history rich in venom-related adaptations, and their quick generation time and suitability for culture in the laboratory.

Identification of the Ricin-B-Lectin LdRBLk in the Colorado Potato Beetle and an Analysis of Its Expression in Response to Fungal Infections

Ricin-B-lectins (RBLs) have been identified in many groups of organisms, including coleopterans insects, particularly the Colorado potato beetle Leptinotarsa decemlineata (LdRBLs). We hypothesized that one of these LdRBLs (LdRBLk) may be involved in the immune response to fungal infections. We performed a theoretical analysis of the structure of this protein. Additionally, the expression levels of the LdRBlk gene were measured in L. decemlineata in response to infections with the fungi Metarhizium robertsii and Beauveria bassiana. The expression levels of LdRBlk in the L. decemlineata cuticle and fat body were increased in response to both infections. The induction of LdRBlk expression was dependent on the susceptibility of larvae to the fungi. Upregulation of the LdRBlk gene was also observed in response to other stresses, particularly thermal burns. Elevation of LdRBlk expression was frequently observed to be correlated with the expression of the antimicrobial peptide attacin but was not correlated with hsp90 regulation. Commercially available β-lectin of ricin from Ricinuscommunis was observed to inhibit the germination of conidia of the fungi. We suggest that LdRBLk is involved in antifungal immune responses in the Colorado potato beetle, either exerting fungicidal properties directly or acting as a modulator of the immune response.

Characterization of a cell death-inducing endonuclease-like venom protein from the parasitoid wasp Pteromalus puparum (Hymenoptera: Pteromalidae)

BACKGROUND

Parasitoid wasps are valuable natural enemies for controlling pests. To ensure successful parasitism, these wasps inject venoms along with their eggs that are deposited either into or on their hosts. Parasitoid venoms regulate host behaviors, development, metabolism and immune responses. Pteromalus puparum is a pupal endoparasitoid that parasitizes a number of butterflies, including the worldwide pest cabbage butterfly, Pieris rapae. Venom from P. puparum has a variety of effects on host hemocytes, including alteration of absolute and relative hemocyte counts, and inhibition of hemocyte spreading and encapsulation. In particular, P. puparum venom causes hemocyte cell death in vivo and in vitro.

RESULTS

Using assay-guided chromatography, a cell death-inducing venom fraction was identified and defined as P. puparum endonuclease-like venom protein (PpENVP). It belongs to the DNA/RNA nonspecific endonuclease family, which contains two conserved endonuclease activation sites. We analyzed its expression profiles and demonstrated that PpENVP inhibits gene expression in transfected cells relying on two activation sites. However, RNA interference of PpENVP did not significantly reduce P. puparum venom cytotoxicity, suggesting that PpENVP may not be the sole cytotoxic factor present.

CONCLUSION

Our results provide novel insight into the function of the P. puparum venom cocktail and identify a promising insecticide candidate endonuclease that targets insect hemocytes.

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.

Metamorphosis-related changes in the free fatty acid profiles of Sarcophaga (Liopygia) argyrostoma (Robineau-Desvoidy, 1830)

The flies of the Sarcophagidae, widespread throughout the temperate zone, are of great significance in Medicine, Veterinary science, Forensics and Entomotoxicology. Lipids are important elements of cell and organelle membranes and a source of energy for embryogenesis, metamorphosis and flight. Cuticular lipids protect from desiccation and act as recognition cues for species, nest mates and castes, and are a source of various pheromones. The free fatty acid (FFA) profile of cuticular and internal extracts of Sarcophaga (Liopygia) argyrostoma (Robineau-Desvoidy, 1830) larvae, pupae and adults was determined by gas chromatography-mass spectrometry (GC-MS). The larvae, pupae and adults contained FFAs from C5:0 to C28:0. The extracts differed quantitatively and qualitatively from each other: C18:1 > C16:1 > C16:0 > C18:0 predominated in the cuticular and internal extracts from the larvae and adults, while 18:1 > C16:0 > C16:1 > C18:0 predominated in the pupae. The FFA profile of the cuticle varies considerably between each development stage: C23:0 and C25:0 are only present in larvae, C28:0 in the pupal cuticle, and C12:1 and C18:3 in internal extracts from adults. The mechanisms underlying this diversity are discussed herein.

Genome Report: Whole Genome Sequence and Annotation of the Parasitoid Jewel Wasp Nasonia giraulti Laboratory Strain RV2X[u]

Jewel wasps in the genus of Nasonia are parasitoids with haplodiploidy sex determination, rapid development and are easy to culture in the laboratory. They are excellent models for insect genetics, genomics, epigenetics, development, and evolution. Nasonia vitripennis (Nv) and N. giraulti (Ng) are closely-related species that can be intercrossed, particularly after removal of the intracellular bacterium Wolbachia, which serve as a powerful tool to map and positionally clone morphological, behavioral, expression and methylation phenotypes. The Nv reference genome was assembled using Sanger, PacBio and Nanopore approaches and annotated with extensive RNA-seq data. In contrast, Ng genome is only available through low coverage resequencing. Therefore, de novo Ng assembly is in urgent need to advance this system. In this study, we report a high-quality Ng assembly using 10X Genomics linked-reads with 670X sequencing depth. The current assembly has a genome size of 259,040,977 bp in 3,160 scaffolds with 38.05% G-C and a 98.6% BUSCO completeness score. 97% of the RNA reads are perfectly aligned to the genome, indicating high quality in contiguity and completeness. A total of 14,777 genes are annotated in the Ng genome, and 72% of the annotated genes have a one-to-one ortholog in the Nv genome. We reported 5 million Ng-Nv SNPs which will facility mapping and population genomic studies in Nasonia In addition, 42 Ng-specific genes were identified by comparing with Nv genome and annotation. This is the first de novo assembly for this important species in the Nasonia model system, providing a useful new genomic toolkit.

The endoparasitoid Psyllaephagus arenarius benefits from ectoparasitic venom via multiparasitism with the ectoparasitoid Tamarixia lyciumi

As solitary nymphal parasitoids of Paratrioza sinica, the ectoparasitoid Tamarixia lyciumi and the endoparasitoid Psyllaephagus arenarius act as effective biocontrol agents. Thus, it is necessary to facilitate mass productions of both species. Despite showing an excellent parasitic ability, Ps. arenarius is often trapped fatally inside 5th-instar nymphs of Pa. sinica due to strong host immunity. To improve the emergence rate of Ps. arenarius, we evaluated whether Ps. arenarius could utilize T. lyciumi venom via multiparasitism, so the parasitism characteristics of both species were examined between separate-existence (monoparasitism only) and co-existence (mono- and multiparasitism) systems. Further, the parasitism characteristics of Ps. arenarius on venom-injected hosts with/without T. lyciumi eggs were tested to further identify the facilitator. The results showed the parasitism rate of T. lyciumi was increased while that of Ps. arenarius did not change from separate-existence to co-existence systems. The intrinsic performances of two species in monoparasitism did not differ between separate- and co-existence systems. From monoparasitism (separate-existence) to multiparasitism (co-existence), no differences were detected in the intrinsic performances of T. lyciumi, but those of Ps. arenarius were greatly improved. After T. lyciumi venom injection, the parasitism characteristics of Ps. arenarius did not differ between venom-injected hosts with T. lyciumi eggs and those without, further indicating Ps. arenarius benefited from the venom of T. lyciumi females rather than T. lyciumi egg/larval secretions. Instead of negative effects, multiparasitism with ectoparasitoids improves endoparasitoids due to ectoparasitic venom. The study increases host resource utilization and provides creative ways for mass production of endoparasitoids.

A chromosome-level genome assembly of the parasitoid wasp Pteromalus puparum

Parasitoid wasps represent a large proportion of hymenopteran species. They have complex evolutionary histories and are important biocontrol agents. To advance parasitoid research, a combination of Illumina short-read, PacBio long-read and Hi-C scaffolding technologies was used to develop a high-quality chromosome-level genome assembly for Pteromalus puparum, which is an important pupal endoparasitoid of caterpillar pests. The chromosome-level assembly has aided in studies of venom and detoxification genes. The assembled genome size is 338 Mb with a contig N50 of 38.7 kb and a scaffold N50 of 1.16 Mb. Hi-C analysis assembled scaffolds onto five chromosomes and raised the scaffold N50 to 65.8 Mb, with more than 96% of assembled bases located on chromosomes. Gene annotation was assisted by RNA sequencing for the two sexes and four different life stages. Analysis detected 98% of the BUSCO (Benchmarking Universal Single-Copy Orthologs) gene set, supporting a high-quality assembly and annotation. In total, 40.1% (135.6 Mb) of the assembly is composed of repetitive sequences, and 14,946 protein-coding genes were identified. Although venom genes play important roles in parasitoid biology, their spatial distribution on chromosomes was poorly understood. Mapping has revealed venom gene tandem arrays for serine proteases, pancreatic lipase-related proteins and kynurenine-oxoglutarate transaminases, which have amplified in the P. puparum lineage after divergence from its common ancestor with Nasonia vitripennis. In addition, there is a large expansion of P450 genes in P. puparum. These examples illustrate how chromosome-level genome assembly can provide a valuable resource for molecular, evolutionary and biocontrol studies of parasitoid wasps.

Lipidomics reveals how the endoparasitoid wasp Pteromalus puparum manipulates host energy stores for its young

Endoparasitoid wasps inject venom along with their eggs to adjust the physiological and nutritional environment inside their hosts to benefit the development of their offspring. In particular, wasp venoms are known to modify host lipid metabolism, lipid storage in the fat body, and release of lipids into the hemolymph, but how venoms accomplish these functions remains unclear. Here, we use an UPLC-MS-based lipidomics approach to analyze the identities and concentrations of lipids in both fat body and hemolymph of host cabbage butterfly (Pieris rapae) infected by the pupal endoparasitoid Pteromalus puparum. During infection, host fat body levels of highly unsaturated, soluble triacylglycerides (TAGs) increased while less unsaturated, less soluble forms decreased. Furthermore, in infected host hemolymph, overall levels of TAG and phospholipids (the major component of cell membranes) increased, suggesting that fat body cells are destroyed and their contents are dispersed. Altogether, these data suggest that wasp venom induces host fat body TAGs to be transformed into lower melting point (more liquid) forms and released into the host hemolymph following infection, allowing simple absorption and nutritional acquisition by wasp larvae. Finally, cholesteryl esters (CEs, a dietary lipid derived from cholesterol) increased in host hemolymph following infection with no concomitant decrease in host cholesterol, implying that the wasp may provide this necessary food resource to its offspring via its venom. This study provides novel insight into how parasitoid infection alters lipid metabolism in insect hosts, and begins to uncover the wasp venom proteins responsible for host physiological changes and offspring development.

Ant behaviour and brain gene expression of defending hosts depend on the ecological success of the intruding social parasite

The geographical mosaic theory of coevolution predicts that species interactions vary between locales. Depending on who leads the coevolutionary arms race, the effectivity of parasite attack or host defence strategies will explain parasite prevalence. Here, we compare behaviour and brain transcriptomes of Temnothorax longispinosus ant workers when defending their nest against an invading social parasite, the slavemaking ant Temnothorax americanus. A full-factorial design allowed us to test whether behaviour and gene expression are linked to parasite pressure on host populations or to the ecological success of parasite populations. Albeit host defences had been shown before to covary with local parasite pressure, we found parasite success to be much more important. Our chemical and behavioural analyses revealed that parasites from high prevalence sites carry lower concentrations of recognition cues and are less often attacked by hosts. This link was further supported by gene expression analysis. Our study reveals that host-parasite interactions are strongly influenced by social parasite strategies, so that variation in parasite prevalence is determined by parasite traits rather than the efficacy of host defence. Gene functions associated with parasite success indicated strong neuronal responses in hosts, including long-term changes in gene regulation, indicating an enduring impact of parasites on host behaviour. This article is part of the theme issue 'The coevolutionary biology of brood parasitism: from mechanism to pattern'.

The Pupal Ectoparasitoid Pachycrepoideus vindemmiae Regulates Cellular and Humoral Immunity of Host Drosophila melanogaster

The immunological interaction between Drosophila melanogaster and its larval parasitoids has been thoroughly investigated, however, little is known about the interaction between the host and its pupal parasitoids. Pachycrepoideus vindemmiae, a pupal ectoparasitoid of D. melanogaster, injects venom into its host while laying eggs on the puparium, which regulates host immunity and interrupts host development. To resist the invasion of parasitic wasps, various immune defense strategies have been developed in their hosts as a consequence of co-evolution. In this study, we mainly focused on the host immunomodulation by P. vindemmiae and thoroughly investigated cellular and humoral immune response, including cell adherence, cell viability, hemolymph melanization and the Toll, Imd, and JAK/STAT immune pathways. Our results indicated that venom had a significant inhibitory effect on lamellocyte adherence and induced plasmatocyte cell death. Venom injection and in vitro incubation strongly inhibited hemolymph melanization. More in-depth investigation revealed that the Toll and Imd immune pathways were immediately activated upon parasitization, followed by the JAK/STAT pathway, which was activated within the first 24 h post-parasitism. These regulatory effects were further validated by qPCR. Our present study manifested that P. vindemmiae regulated the cellular and humoral immune system of host D. melanogaster in many aspects. These findings lay the groundwork for studying the immunological interaction between D. melanogaster and its pupal parasitoid.

Evaluating the evolution and function of the dynamic Venom Y protein in ectoparasitoid wasps

Venom of the parasitoid wasp Nasonia vitripennis changes the metabolism and gene expression in its fly host Sarcophaga bullata to induce developmental arrest, suppression of the immune response and various other venom effects. Yet, the venom of ectoparasitoid wasps has not been fully characterized. A major component of N. vitripennis venom is an uncharacterized, high-expressing protein referred to as Venom Y. Here we describe the evolutionary history and possible functions of this venom protein. We found that Venom Y is a relatively young gene that has duplicated to form two distinct paralogue groups. A copy of Venom Y has been recruited as a venom protein in at least five wasp species. Functional analysis found that Venom Y affects detoxification and immunity genes in envenomated fly hosts. Many of these genes are fat-body specific, suggesting that Venom Y may have a targeted effect on fat body tissue. We also show that Venom Y may mitigate negative effects of other venom proteins. Finally, protein sequencing indicates that Venom Y is post-translationally modified. This study contributes to elucidating parasitoid venom by using RNA interference knockdown to investigate venom protein function in the context of the whole venom cocktail.

Bioinformatic analysis suggests potential mechanisms underlying parasitoid venom evolution and function

Hymenopteran parasitoid wasps are a diverse collection of species that infect arthropod hosts and use factors found in their venoms to manipulate host immune responses, physiology, and behaviour. Whole parasitoid venoms have been profiled using proteomic approaches, and here we present a bioinformatic characterization of the venom protein content from Ganaspis sp. 1, a parasitoid that infects flies of the genus Drosophila. We find evidence that diverse evolutionary processes including multifunctionalization, co-option, gene duplication, and horizontal gene transfer may be acting in concert to drive venom gene evolution in Ganaspis sp.1. One major role of parasitoid wasp venom is host immune evasion. We previously demonstrated that Ganaspis sp. 1 venom inhibits immune cell activation in infected Drosophila melanogaster hosts, and our current analysis has uncovered additional predicted virulence functions. Overall, this analysis represents an important step towards understanding the composition and activity of parasitoid wasp venoms.

Genome and Ontogenetic-Based Transcriptomic Analyses of the Flesh Fly, Sarcophaga bullata

The flesh fly, Sarcophaga bullata, is a widely-used model for examining the physiology of insect diapause, development, stress tolerance, neurobiology, and host-parasitoid interactions. Flies in this taxon are implicated in myiasis (larval infection of vertebrates) and feed on carrion, aspects that are important in forensic studies. Here we present the genome of S. bullata, along with developmental- and reproduction-based RNA-Seq analyses. We predict 15,768 protein coding genes, identify orthology in relation to closely related flies, and establish sex and developmental-specific gene sets based on our RNA-Seq analyses. Genomic sequences, predicted genes, and sequencing data sets have been deposited at the National Center for Biotechnology Information. Our results provide groundwork for genomic studies that will expand the flesh fly’s utility as a model system.

Insights into the venom protein components of Microplitis mediator, an endoparasitoid wasp

Endoparasitoid wasps deliver a variety of maternal factors, such as venom proteins, viruses, and virus-like particles, from their venom and calyx fluid into hosts and thereby regulate the hosts' immune response, metabolism and development. The endoparasitoid, Microplitis mediator, is used as an important biological agent for controlling the devastating pest Helicoverpa armigera. In this study, using an integrated transcriptomic and proteomic analysis approach, we identified 75 putative venom proteins in M. mediator. The identified venom components were consistent with other known parasitoid wasps' venom proteins, including metalloproteases, serine protease inhibitors, and glycoside hydrolase family 18 enzymes. The metalloprotease and serpin family showed extensive gene duplications in venom apparatus. Isobaric tags for relative and absolute quantitation (iTRAQ) based quantitative proteomics revealed 521 proteins that were differentially expressed at 6 h and 24 h post-parasitism, including 10 wasp venom proteins that were released into the host hemolymph. Further analysis indicated that 511 differentially expressed proteins (DEP) from the host are primarily involved in the immune response, material metabolism, and extracellular matrix receptor interaction. Taken together, our results on parasitoid wasp venoms have the potential to enhance the application of endoparasitoid wasps for controlling insect pest.

Venom is beneficial but not essential for development and survival of Nasonia

Parasitoid wasps sting and inject venom in arthropod hosts, which alters host metabolism and development while keeping the host alive for several days, presumably to induce benefits for the parasitoid young.Here we investigate the consequences of host envenomation on development and fitness of wasp larvae in the ectoparasitoid Nasonia vitripennis, by comparing wasps reared on live unstung, previously stung, and cold-killed hosts. Developmental arrest and suppression of host response to larvae are major venom effects that occur in both stung and cold-killed hosts, but not unstung hosts; while cold-killed hosts lack venom effects that require a living host. Thus, cold-killed hosts mimic some of the effects of venom, but not others.Eggs placed on live unstung hosts have significantly higher mortality during development, however successfully developing wasps from these hosts have similar lifetime fecundity to wasps from cold-killed or stung hosts. Therefore, although venom is beneficial, it is not required for wasp survival.While wasps developing on cold-killed versus stung hosts have similar fitness, multiple generations of rearing on cold-killed hosts results in significant fitness reductions of wasps.We conclude that the largest benefits of venom are induction of host developmental arrest and suppression of host response to larva (e.g. immune responses), although more subtle benefits may accrue across generations, or under stressful conditions.

The genomic and transcriptomic analyses of serine proteases and their homologs in an endoparasitoid, Pteromalus puparum

In insects, serine proteases (SPs) and serine protease homologs (SPHs) constitute a large family of proteins involved in multiple physiological processes such as digestion, development, and immunity. Here we identified 145 SPs and 38 SPHs in the genome of an endoparasitoid, Pteromalus puparum. Gene duplication and tandem repeats were observed in this large SPs/SPHs family. We then analyzed the expression profiles of SP/SPH genes in response to different microbial infections (Gram-positive bacterium Micrococcus luteus, Gram-negative bacterium Escherichia coli, and entomopathogenic fungus Beauveria bassiana), as well as in different developmental stages and tissues. Some SPs/SPHs also displayed distinct expression patterns in venom gland, suggesting their specific physiological functions as venom proteins. Our finding lays groundwork for further research of SPs and SPHs expressed in the venom glands.

Induction of stress- and immune-associated genes in the Indian meal moth Plodia interpunctella against envenomation by the ectoparasitoid Bracon hebetor

Envenomation is an important process in parasitism by parasitic wasps; it suppresses the immune and development of host insects. However, the molecular mechanisms of host responses to envenomation are not yet clear. This study aimed to determine the transcription-level responses of the Indian meal moth Plodia interpunctella against envenomation of the ectoparasitoid Bracon hebetor. Quantitative real-time reverse-transcription PCR was used to determine the transcriptional changes of 13 selected genes, which are associated with development, metabolism, stress, or immunity, in the feeding and wandering fifth instar larvae over a 4-day period after envenomation. The effects of envenomation on the feeding-stage larvae were compared with those of starvation in the transcriptional levels of the 13 genes. Most selected genes were altered in their expression by either envenomation or starvation. In particular, a heat shock protein, hsp70, was highly upregulated in envenomated larvae in both the feeding and wandering stages as well as in starved larvae. Further, some genes were upregulated by envenomation in a stage-specific manner. For example, hsp25 was upregulated after envenomation in the feeding larvae, but hsp90 and an immune-associated gene, hemolin, were upregulated in the wandering larvae. However, both envenomation and starvation resulted in the downregulation of genes associated with development and metabolism. Taken together, P. interpunctella upregulated stress- and immune-responsive genes, but downregulated genes associated with development and metabolism after envenomation. This study provides important information for understanding the molecular mechanisms of host responses to parasitism.

Host regulation by the ectophagous parasitoid wasp Bracon nigricans

The host regulation process has been widely investigated in endophagous parasitoid wasps, which in most cases finely interact with living hosts (i.e. koinobiont parasitoids). In contrast, only very limited information is available for ectophagous parasitoids that permanently paralyze and rapidly suppress their victims (i.e. idiobiont parasitoids). Here we try to fill this research gap by investigating the host regulation by Bracon nigricans, an ectophagous idiobiont wasp species. Parasitism, mainly by venom action, is able to redirect host metabolism in order to enhance its nutritional suitability for the developing parasitoid larvae and to provide the required metabolic support to host tissues. The observed alterations of the host titers of haemolymph proteins, carbohydrates and acylglycerols are associated with a parasitoid-induced mobilization of nutrients stored in the fat body. This tissue undergoes a controlled degradation mediated by a close surface interaction with haemocytes, where a cathepsin L activity is localized, as demonstrated by immunolocalization, biochemical and transcriptional data. B. nigricans parasitism does not markedly influence the survival of haemocytes, even though a persistent suppression of the immune competence is observed in parasitized hosts, which show a reduced capacity to encapsulate and melanize non-self objects. These immune alterations likely allow a more efficient food uptake and use by the ectophagous larvae. The obtained results indicate that the host regulation process in basal lineages of parasitic Hymenoptera is more complex than expected and shares functional similarities with adaptive strategies occurring in derived koinobiont species.

Transcriptional Profiling of Midguts Prepared from Trypanosoma/T. congolense-Positive Glossina palpalis palpalis Collected from Two Distinct Cameroonian Foci: Coordinated Signatures of the Midguts' Remodeling As T. congolense-Supportive Niches

Our previous transcriptomic analysis of Glossina palpalis gambiensis experimentally infected or not with Trypanosoma brucei gambiense aimed to detect differentially expressed genes (DEGs) associated with infection. Specifically, we selected candidate genes governing tsetse fly vector competence that could be used in the context of an anti-vector strategy, to control human and/or animal trypanosomiasis. The present study aimed to verify whether gene expression in field tsetse flies (G. p. palpalis) is modified in response to natural infection by trypanosomes (T. congolense), as reported when insectary-raised flies (G. p. gambiensis) are experimentally infected with T. b. gambiense. This was achieved using the RNA-seq approach, which identified 524 DEGs in infected vs. non-infected tsetse flies, including 285 downregulated genes and 239 upregulated genes (identified using DESeq2). Several of these genes were highly differentially expressed, with log2 fold change values in the vicinity of either +40 or −40. Downregulated genes were primarily involved in transcription/translation processes, whereas encoded upregulated genes governed amino acid and nucleotide biosynthesis pathways. The BioCyc metabolic pathways associated with infection also revealed that downregulated genes were mainly involved in fly immunity processes. Importantly, our study demonstrates that data on the molecular cross-talk between the host and the parasite (as well as the always present fly microbiome) recorded from an experimental biological model has a counterpart in field flies, which in turn validates the use of experimental host/parasite couples.

An alternative pathway to eusociality: Exploring the molecular and functional basis of fortress defense

Some animals express a form of eusociality known as "fortress defense," in which defense rather than brood care is the primary social act. Aphids are small plant-feeding insects, but like termites, some species express division of labor and castes of aggressive juvenile "soldiers." What is the functional basis of fortress defense eusociality in aphids? Previous work showed that the acquisition of venoms might be a key innovation in aphid social evolution. We show that the lethality of aphid soldiers derives in part from the induction of exaggerated immune responses in insects they attack. Comparisons between closely related social and nonsocial species identified a number of secreted effector molecules that are candidates for immune modulation, including a convergently recruited protease described in unrelated aphid species with venom-like functions. These results suggest that aphids are capable of antagonizing conserved features of the insect immune response, and provide new insights into the mechanisms underlying the evolution of fortress defense eusociality in aphids.

The Evolution of Venom by Co-option of Single-Copy Genes

The classic model for the evolution of novel gene function is through gene duplication followed by evolution of a new function by one of the copies (neofunctionalization) [1, 2]. However, other modes have also been found, such as novel genes arising from non-coding DNA, chimeric fusions, and lateral gene transfers from other organisms [3-7]. Here we use the rapid turnover of venom genes in parasitoid wasps to study how new gene functions evolve. In contrast to the classic gene duplication model, we find that a common mode of acquisition of new venom genes in parasitoid wasps is co-option of single-copy genes from non-venom progenitors. Transcriptome and proteome sequencing reveal that recruitment and loss of venom genes occur primarily by rapid cis-regulatory expression evolution in the venom gland. Loss of venom genes is primarily due to downregulation of expression in the gland rather than gene death through coding sequence degradation. While the majority of venom genes have specialized expression in the venom gland, recent losses of venom function occur primarily among genes that show broader expression in development, suggesting that they can more readily switch functional roles. We propose that co-option of single-copy genes may be a common but relatively understudied mechanism of evolution for new gene functions, particularly under conditions of rapid evolutionary change.

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.

The venom gland transcriptome of the parasitoid wasp Nasonia vitripennis highlights the importance of novel genes in venom function

Background

Prior to egg laying the parasitoid wasp Nasonia vitripennis envenomates its pupal host with a complex mixture of venom peptides. This venom induces several dramatic changes in the host, including developmental arrest, immunosuppression, and altered metabolism. The diverse and potent bioactivity of N. vitripennis venom provides opportunities for the development of novel acting pharmaceuticals based on these molecules. However, currently very little is known about the specific functions of individual venom peptides or what mechanisms underlie the hosts response to envenomation. Many of the venom peptides also lack bioinformatically derived annotations because no homologs can be identified in the sequences databases. The RNA interference system of N. vitripennis provides a method for functional characterisation of venom protein encoding genes, however working with the current list of 79 candidates represents a daunting task. For this reason we were interested in determining the expression levels of venom encoding genes in the venom gland, as this information could be used to rank candidates for further study. To do this we carried out deep transcriptome sequencing of the venom gland and ovary tissue and used RNA-seq to rank the venom protein encoding genes by expression level. The generation of a specific venom gland transcriptome dataset also provides further opportunities to investigate novel features of this specialised organ.

Results

RNA-seq revealed that the highest expressed venom encoding gene in the venom gland was ‘Venom protein Y’. The highest expressed annotated gene in this tissue was serine protease Nasvi2EG007167, which has previously been implicated in the apoptotic activity of N. vitripennis venom. As expected the RNA-seq confirmed that venom encoding genes are almost exclusively expressed in the venom gland relative to the neighbouring ovary tissue. Novel genes appear to perform key roles in N. vitripennis venom function, with over half of the 15 highest expressed venom encoding loci lacking bioinformatic annotations. The high throughput sequencing data also provided evidence for the existence of an additional 472 previously undescribed transcribed regions in the N. vitripennis genome. Finally, metatranscriptomic analysis of the venom gland transcriptome finds little evidence for the role of Wolbachia in the venom system.

Conclusions

The expression level information provided here for the N. vitripennis venom protein encoding genes represents a valuable dataset that can be used by the research community to rank candidates for further functional characterisation. These candidates represent bioactive peptides valuable in the development of new pharmaceuticals.

Electronic supplementary material

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

Effects of Scleroderma sichuanensis Xiao (Hymenoptera: Bethylidae) venom and parasitism on nutritional content regulation in host Tenebrio molitor L. (Coleoptera: Tenebrionidae)

To explore the mechanisms by which the wasp Scleroderma sichuanensis Xiao regulates the physiology and biochemistry of its host, effects of S. sichuanensis venom and parasitism on host the Tenebrio molitor L. pupae were examined. Significant differences in nutritional content were noted between parasitized and non-parasitized pupae and between venom- and phosphate buffered saline-injected pupae. When pupae were injected with venom, the fat body could not be disintegrated into granules; however, when pupae were parasitized, fat-body disintegration occurred. Electrophoresis showed no differences in hemolymph protein content between parasitized pupae and those injected with venom, indicating that the wasp did not have narrow-spectrum peptides. These findings confirmed that S. sichuanensis was a typical idiobiont ectoparasitoid wasp, and that nutrient regulation was similar between idiobiont and koinobiont wasps. The strong similarities between the two treatments suggest that venom injection is a major factor responsible for changes in host nutrient content. The wasp fed mainly on reducing sugars, free amino acids, and fat-body tissues; larval fat bodies were derived from hemolymph and from host tissue. Our findings suggest that lipid catabolism might be accelerated, and that lipid biosynthesis might be inhibited, when host pupae are parasitized or injected with venom. In addition to venom, physiological and biochemical changes that occur during the parasitic process might be caused by venom, ovarian proteins, saliva, or secretions.

Laterally Transferred Gene Recruited as a Venom in Parasitoid Wasps

Parasitoid wasps use venom to manipulate the immunity and metabolism of their host insects in a variety of ways to provide resources for their offspring. Yet, how genes are recruited and evolve to perform venom functions remain open questions. A recently recognized source of eukaryotic genome innovation is lateral gene transfer (LGT). Glycoside hydrolase family 19 (GH19) chitinases are widespread in bacteria, microsporidia, and plants where they are used in nutrient acquisition or defense, but have previously not been known in metazoans. In this study, a GH19 chitinase LGT is described from the unicellular microsporidia/Rozella clade into parasitoid wasps of the superfamily Chalcidoidea, where it has become recruited as a venom protein. The GH19 chitinase is present in 15 species of chalcidoid wasps representing four families, and phylogenetic analysis indicates that it was laterally transferred near or before the origin of Chalcidoidea (∼95 Ma). The GH19 chitinase gene is highly expressed in the venom gland of at least seven species, indicating a role in the complex host manipulations performed by parasitoid wasp venom. RNAi knockdown in the model parasitoid Nasonia vitripennis reveals that-following envenomation-the GH19 chitinase induces fly hosts to upregulate genes involved in an immune response to fungi. A second, independent LGT of GH19 chitinase from microsporidia into mosquitoes was also found, also supported by phylogenetic reconstructions. Besides these two LGT events, GH19 chitinase is not found in any other sequenced animal genome, or in any fungi outside the microsporidia/Rozella clade.

Venom of Parasitoid Pteromalus puparum Impairs Host Humoral Antimicrobial Activity by Decreasing Host Cecropin and Lysozyme Gene Expression

Insect host/parasitoid interactions are co-evolved systems in which host defenses are balanced by parasitoid mechanisms to disable or hide from host immune effectors. Here, we report that Pteromalus puparum venom impairs the antimicrobial activity of its host Pieris rapae. Inhibition zone results showed that bead injection induced the antimicrobial activity of the host hemolymph but that venom inhibited it. The cDNAs encoding cecropin and lysozyme were screened. Relative quantitative PCR results indicated that all of the microorganisms and bead injections up-regulated the transcript levels of the two genes but that venom down-regulated them. At 8 h post bead challenge, there was a peak in the transcript level of the cecropin gene, whereas the peak of lysozyme gene occurred at 24 h. The transcripts levels of the two genes were higher in the granulocytes and fat body than in other tissues. RNA interference decreased the transcript levels of the two genes and the antimicrobial activity of the pupal hemolymph. Venom injections similarly silenced the expression of the two genes during the first 8 h post-treatment in time- and dose-dependent manners, after which the silence effects abated. Additionally, recombinant cecropin and lysozyme had no significant effect on the emergence rate of pupae that were parasitized by P. puparum females. These findings suggest one mechanism of impairing host antimicrobial activity by parasitoid venom.

Insights into the venom composition and evolution of an endoparasitoid wasp by combining proteomic and transcriptomic analyses

Parasitoid wasps are abundant and diverse hymenopteran insects that lay their eggs into the internal body (endoparasitoid) or on the external surface (ectoparasitoid) of their hosts. To make a more conducive environment for the wasps’ young, both ecto- and endoparasitoids inject venoms into the host to modulate host immunity, metabolism and development. Endoparasitoids have evolved from ectoparasitoids independently in different hymenopteran lineages. Pteromalus puparum, a pupal endoparasitoid of various butterflies, represents a relatively recent evolution of endoparasitism within pteromalids. Using a combination of transcriptomic and proteomic approaches, we have identified 70 putative venom proteins in P. puparum. Most of them show higher similarity to venom proteins from the related ectoparasitoid Nasonia vitripennis than from other more distantly related endoparasitoids. In addition, 13 venom proteins are similar to venoms of distantly related endoparasitoids but have no detectable venom matches in Nasonia. These venom proteins may have a role in adaptation to endoparasitism. Overall, these results lay the groundwork for more detailed studies of venom function and adaptation to the endoparasitic lifestyle.

The expanding genetic toolbox of the wasp Nasonia vitripennis and its relatives

The parasitoid wasp Nasonia represents a genus of four species that is emerging as a powerful genetic model system that has made and will continue to make important contributions to our understanding of evolutionary biology, development, ecology, and behavior. Particularly powerful are the haplodiploid genetics of the system, which allow some of the advantages of microbial genetics to be applied to a complex multicellular eukaryote. In addition, fertile, viable hybrids can be made among the four species in the genus. This makes Nasonia exceptionally well suited for evolutionary genetics approaches, especially when combined with its haploid genetics and tractability in the laboratory. These features are complemented by an expanding array of genomic, transcriptomic, and functional resources, the application of which has already made Nasonia an important model system in such emerging fields as evolutionary developmental biology and microbiomics. This article describes the genetic and genomic advantages of Nasonia wasps and the resources available for their genetic analysis.

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.

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.

Origin and evolution of the enhancer of split complex

BACKGROUND

The Enhancer of split complex is an unusual gene complex found in Arthropod genomes. Where known this complex of genes is often regulated by Notch cell signalling and is critically important for neurogenesis. The Enhancer of split complex is made up of two different classes of genes, basic helix-loop-helix-orange domain transcription factors and bearded class genes. The association of these genes has been detected in the genomes of insects and crustaceans.

RESULTS

Tracing the evolution of the Enhancer of split complex in recently sequenced Arthropod genomes indicates that enhancer of split basic helix-loop-helix orange domain genes arose before the common ancestor of insects and Crustacea, and before the formation of the complex. Throughout insect and crustacean evolution, a four-gene cluster has been present with lineage specific gene losses and duplications. The complex can be found in the vast majority of genomes, but appears to be missing from the genomes of chalcid wasps, raising questions as to how they carry out neurogenesis in the absence of these crucial genes.

CONCLUSIONS

The enhancer of split complex arose in the common ancestor of Crustacea and insects, probably through the linkage of a basic helix-loop-helix orange domain gene and a bearded class gene. The complex has been maintained, with variations, throughout insect and crustacean evolution indicating some function of the complex, such as coordinate regulation, may maintain its structure through evolutionary time.

Ecological venomics: How genomics, transcriptomics and proteomics can shed new light on the ecology and evolution of venom

Animal venom is a complex cocktail of bioactive chemicals that traditionally drew interest mostly from biochemists and pharmacologists. However, in recent years the evolutionary and ecological importance of venom is realized as this trait has direct and strong influence on interactions between species. Moreover, venom content can be modulated by environmental factors. Like many other fields of biology, venom research has been revolutionized in recent years by the introduction of systems biology approaches, i.e., genomics, transcriptomics and proteomics. The employment of these methods in venom research is known as 'venomics'. In this review we describe the history and recent advancements of venomics and discuss how they are employed in studying venom in general and in particular in the context of evolutionary ecology. We also discuss the pitfalls and challenges of venomics and what the future may hold for this emerging scientific field.

A new approach for investigating venom function applied to venom calreticulin in a parasitoid wasp

A new method is developed to investigate functions of venom components, using venom gene RNA interference knockdown in the venomous animal coupled with RNA sequencing in the envenomated host animal. The vRNAi/eRNA-Seq approach is applied to the venom calreticulin component (v-crc) of the parasitoid wasp Nasonia vitripennis. Parasitoids are common, venomous animals that inject venom proteins into host insects, where they modulate physiology and metabolism to produce a better food resource for the parasitoid larvae. vRNAi/eRNA-Seq indicates that v-crc acts to suppress expression of innate immune cell response, enhance expression of clotting genes in the host, and up-regulate cuticle genes. V-crc KD also results in an increased melanization reaction immediately following envenomation. We propose that v-crc inhibits innate immune response to parasitoid venom and reduces host bleeding during adult and larval parasitoid feeding. Experiments do not support the hypothesis that v-crc is required for the developmental arrest phenotype observed in envenomated hosts. We propose that an important role for some venom components is to reduce (modulate) the exaggerated effects of other venom components on target host gene expression, physiology, and survival, and term this venom mitigation. A model is developed that uses vRNAi/eRNA-Seq to quantify the contribution of individual venom components to total venom phenotypes, and to define different categories of mitigation by individual venoms on host gene expression. Mitigating functions likely contribute to the diversity of venom proteins in parasitoids and other venomous organisms.