Summary 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 down-regulation 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.

Parasitoid wasps are diverse and ecologically important insects that use venom to modify their host’s metabolism for the benefit of the parasitoid’s offspring. Thus, the effects of venom can be considered an ‘extended phenotype’ of the wasp. The model parasitoid wasp Nasonia vitripennis has approximately 100 venom proteins, 23 of which do not have sequence similarity to known proteins. Envenomation by N. vitripennis has previously been shown to induce developmental arrest, selective apoptosis and alterations in lipid metabolism in flesh fly hosts. However, the full effects of Nasonia venom are still largely unknown. In this study, we used high throughput RNA sequencing (RNA-Seq) to characterize global changes in Sarcophaga bullata (Diptera) gene expression in response to envenomation by N. vitripennis. Surprisingly, we show that Nasonia venom targets a small subset of S. bullata loci, with ~2% genes being differentially expressed in response to envenomation. Strong upregulation of enhancer of split complex genes provides a potential molecular mechanism that could explain the observed neural cell death and developmental arrest in envenomated hosts. Significant increases in antimicrobial peptides and their corresponding regulatory genes provide evidence that venom could be selectively activating certain immune responses of the hosts. Further, we found differential expression of genes in several metabolic pathways, including glycolysis and gluconeogenesis that may be responsible for the decrease in pyruvate levels found in envenomated hosts. The targeting of Nasonia venom effects to a specific and limited set of genes provides insight into the interaction between the ectoparasitoid wasp and its host.

Venom proteins evolve rapidly, and as a trophic adaptation are excellent models for predator-prey evolutionary studies. The key to a deeper understanding of venom evolution is an integrated approach, combining prey assays with analysis of venom gene expression and venom phenotype. Here, we use such an approach to study venom evolution in the Amazon puffing snake, , a generalist feeder. We identify two novel three-finger toxins: sulditoxin and sulmotoxin 1. These new toxins are not only two of the most abundant venom proteins, but are also functionally intriguing, displaying distinct prey-specific toxicities. Sulditoxin is highly toxic towards lizard prey, but is non-toxic towards mammalian prey, even at greater than 22-fold higher dosage. By contrast, sulmotoxin 1 exhibits the reverse trend. Furthermore, evolutionary analysis and structural modelling show highest sequence variability in the central loop of these proteins, probably driving taxon-specific toxicity. This is, to our knowledge, the first case in which a bimodal and contrasting pattern of toxicity has been shown for proteins in the venom of a single snake in relation to diet. Our study is an example of how toxin gene neofunctionalization can result in a venom system dominated by one protein superfamily and still exhibit flexibility in prey capture efficacy.

Parasitoid wasps inject insect hosts with a cocktail of venoms to manipulate the physiology, development, and immunity of the hosts and to promote development of the parasitoid offspring. The jewel wasp Nasonia vitripennis is a model parasitoid with at least 79 venom proteins. We conducted a high-throughput analysis of Nasonia venom effects on temporal changes of 249 metabolites in pupae of the flesh fly host (Sarcophaga bullata), over a five-day time course. Our results show that venom does not simply arrest the metabolism of the fly host. Rather, it targets specific metabolic processes while keeping hosts alive for at least five days post venom injection by the wasp. We found that venom: (a) Activates the sorbitol biosynthetic pathway while maintaining stable glucose levels, (b) Causes a shift in intermediary metabolism by switching to anaerobic metabolism and blocking the tricarboxylic acid cycle, (c) Arrests chitin biosynthesis that likely reflects developmental arrest of adult fly structures, (d) Elevates the majority of free amino acids, and (e) May be increasing phospholipid degradation. Despite sharing some metabolic effects with cold treatment, diapause, and hypoxia, the venom response is distinct from these conditions. Because Nasonia venom dramatically increases sorbitol levels without changing glucose levels, it could be a useful model for studying the regulation of the sorbitol pathway, which is relevant to diabetes research. Our findings generally support the view that parasitoid venoms are a rich source of bioactive molecules with potential biomedical applications.

No abstract

No abstract

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.

Using an integrated transcriptomic and proteomic approach, we characterized the venom peptidome of the European red ant, Manica rubida. We identified 13 “myrmicitoxins” that share sequence similarities with previously identified ant venom peptides, one of them being identified as an EGF-like toxin likely resulting from a threonine residue modified by O-fucosylation. Furthermore, we conducted insecticidal assays of reversed-phase HPLC venom fractions on the blowfly Lucilia caesar, permitting us to identify six myrmicitoxins (i.e., U3-, U10-, U13-, U20-MYRTX-Mri1a, U10-MYRTX-Mri1b, and U10-MYRTX-Mri1c) with an insecticidal activity. Chemically synthesized U10-MYRTX-Mri1a, -Mri1b, -Mri1c, and U20-MYRTX-Mri1a irreversibly paralyzed blowflies at the highest doses tested (30–125 nmol·g–1). U13-MYRTX-Mri1a, the most potent neurotoxic peptide at 1 h, had reversible effects after 24 h (150 nmol·g–1). Finally, U3-MYRTX-Mri1a has no insecticidal activity, even at up to 55 nmol·g–1. Thus, M. rubida employs a paralytic venom rich in linear insecticidal peptides, which likely act by disrupting cell membranes.