Parasitoid wasps represent a group of parasitic insects with high species diversity that have played a pivotal role in biological control and evolutionary studies. Over the past 20 years, developments in genomics have greatly enhanced our understanding of the biology of these species. Technological leaps in sequencing have facilitated the improvement of genome quality and quantity, leading to the availability of hundreds of parasitoid wasp genomes. Here, we summarize recent progress in parasitoid wasp genomics, focusing on the evolution of genome size (GS) and the genomic basis of several key traits. We also discuss the contributions of genomics in studying venom evolution and endogenization of viruses. Finally, we advocate for increased sequencing and functional research to better understand parasitoid biology and enhance biological control.

Venomous marine gastropods of the superfamily Conoidea possess a rich arsenal of toxins, including neuroactive toxins. Venom adaptations might have played a fundamental role in the radiation of conoideans; nevertheless, there is still no knowledge about the venom of the most diversified family of the group: Raphitomidae Bellardi, 1875. In this study, transcriptomes were produced from the carcase, salivary glands, and proximal and distal venom ducts of the northeastern Atlantic species Raphitoma purpurea (Montagu, 1803). Using a gut barcoding approach, we were also able to report, for the first time, molecular evidence of a vermivorous diet for the genus. Transcriptomic analyses revealed over a hundred putative venom components (PVC), including 69 neurotoxins. Twenty novel toxin families, including some with high levels of expansion, were discovered. No significant difference was observed between the distal and proximal venom duct secretions. Peptides related to cone snail toxins (Cerm06, Pgam02, and turritoxin) and other venom-related proteins (disulfide isomerase and elevenin) were retrieved from the salivary glands. These salivary venom components may constitute ancestral adaptations for venom production in conoideans. Although often neglected, salivary gland secretions are of extreme importance for understanding the evolutionary history of conoidean venom.

Gene duplication followed by nucleotide differentiation is one of the simplest mechanisms to develop new functions for genes. However, the evolutionary processes underlying the divergence of multigene families remain controversial. We used multigene families found within the diversity of toxic proteins in centipede venom to test two hypotheses related to venom evolution: the two-speed mode of venom evolution and the rapid accumulation of variation in exposed residues (RAVER) model. The two-speed mode of venom evolution proposes that different types of selection impact ancient and younger venomous lineages with negative selection being the predominant form in ancient lineages and positive selection being the dominant form in younger lineages. The RAVER hypothesis proposes that, instead of different types of selection acting on different ages of venomous lineages, the different types of selection will selectively contribute to amino acid variation based on whether the residue is exposed to the solvent where it can potentially interact directly with toxin targets. This hypothesis parallels the longstanding understanding of protein evolution that suggests that residues found within the structural or active regions of the protein will be under negative or purifying selection, and residues that do not form part of these areas will be more prone to positive selection. To test these two hypotheses, we compared the venom of 26 centipedes from the order Scolopendromorpha from six currently recognized species from across North America using both transcriptomics and proteomics. We first estimated their phylogenetic relationships and uncovered paraphyly among the genus Scolopendra and evidence for cryptic diversity among currently recognized species. Using our phylogeny, we then characterized the diverse venom components from across the identified clades using a combination of transcriptomics and proteomics. We conducted selection-based analyses in the context of predicted three-dimensional properties of the venom proteins and found support for both hypotheses. Consistent with the two-speed hypothesis, we found a prevalence of negative selection across all proteins. Consistent with the RAVER hypothesis, we found evidence of positive selection on solvent-exposed residues, with structural and less-exposed residues showing stronger signal for negative selection. Through the use of phylogenetics, transcriptomics, proteomics, and selection-based analyses, we were able to describe the evolution of venom from an ancient venomous lineage and support principles of protein evolution that directly relate to multigene family evolution.

Venom production has evolved independently many times in the animal kingdom, although it is rare among mammals. Venomous shrews produce venom in their submandibular salivary glands and use it for food acquisition. Only a few toxins have been identified in shrew venoms thus far, and their modes of action require investigation. The biological and molecular processes relating to venom production and gland functioning also remain unknown. To address this gap, we investigated protein content in extracts from venom glands of two shrew species, Neomys fodiens and Sorex araneus, and interpreted their biological functions. Applying a proteomic approach coupled with Gene Ontology enrichment analysis, we identified 313 and 187 putative proteins in venom glands of N. fodiens and S. araneus, respectively. A search of the UniProt database revealed that most of the proteins found in both shrew species were involved in metabolic processes and stress response, while GO enrichment analysis revealed more stress-related proteins in the glands of S. araneus. Molecules that regulate molecule synthesis, cell cycles, and cell divisions are necessary to enable venom regeneration and ensure its effectiveness in predation and food hoarding. The presence of proteins involved in stress response may be the result of shrews\' high metabolic rate and the costs of venom replenishment. Some proteins are likely to promote toxin spreading during envenomation and, due to their proteolytic action, reinforce venom toxicity. Finally, finding numerous proteins involved in immune response suggests a potential role of shrew venom gland secretions in protection against pathogens. These findings open up new perspectives for studying biological functions of molecules from shrew venom glands and extend our knowledge on the functioning of eulipotyphlan venom systems. Because the majority of existing and putative venomous mammals use oral venom systems to inject venom into target species, the methods presented here provide a promising avenue for confirming or discovering new taxa of venomous mammals.

Venom systems are complex traits that have independently emerged multiple times in diverse plant and animal phyla. Within each venomous lineage there typically exists interspecific variation in venom composition where several factors have been proposed as drivers of variation, including phylogeny and diet. Understanding these factors is of broad biological interest and has implications for the development of antivenom therapies and venom-based drug discovery. Because of their high species richness and the presence of several major evolutionary prey shifts, venomous marine cone snails (genus Conus) provide an ideal system to investigate drivers of interspecific venom variation. Here, by analyzing the venom gland expression profiles of ∼3,000 toxin genes from 42 species of cone snail, we elucidate the role of prey-specific selection pressures in shaping venom variation. By analyzing overall venom composition and individual toxin structures, we demonstrate that the shifts from vermivory to piscivory in Conus are complemented by distinct changes in venom composition independent of phylogeny. In vivo injections of venom from piscivorous cone snails in fish further showed a higher potency compared with venom of nonpiscivores demonstrating a selective advantage. Together, our findings provide compelling evidence for the role of prey shifts in directing the venom composition of cone snails and expand our understanding of the mechanisms of venom variation and diversification.

Venom plays a crucial role in the defense and predation of venomous animals. Spiders (Araneae) are among the most successful predators and have a fascinating venom composition. Their venom mainly contains disulfide-rich peptides and large proteins. Here, we analyzed spider venom protein families, utilizing transcriptomic and genomic data, and highlighted their similarities and differences. We show that spiders have specific combinations of toxins for better predation and defense, typically comprising a core toxin expressed alongside several auxiliary toxins. Among them, the CAP superfamily is widely distributed and highly expressed in web-building Araneoidea spiders. Our analysis of evolutionary relationships revealed four subfamilies (subA-subD) of the CAP superfamily that differ in structure and potential functions. CAP proteins are composed of a conserved CAP domain and diverse C-terminal domains. CAP subC shares similar domains with the snake ion channel regulator svCRISP proteins, while CAP subD possesses a sequence similar to that of insect venom allergen 5 (Ag5). Furthermore, we show that gene duplication and selective expression lead to increased expression of CAP subD, making it a core member of the CAP superfamily. This study sheds light on the functional diversity of CAP subfamilies and their evolutionary history, which has important implications for fully understanding the composition of spider venom proteins and the core toxin components of web-building spiders.

Venomous marine gastropods of the family Conidae are among the most diversified predators in marine realm-in large due to their complex venoms. Besides being a valuable source of bioactive neuropeptides conotoxins, cone-snails venoms are an excellent model for molecular evolution studies, addressing origin of key innovations. However, these studies are handicapped by scarce current knowledge on the tissues involved in venom production, as it is generally assumed the sole prerogative of the venom gland (VG). The role of other secretory glands that are present in all Conus species (salivary gland, SG) or only in some species (accessory salivary gland, ASG) remains poorly understood. Here, for the first time, we carry out a detailed analysis of the VG, SG, and ASG transcriptomes in the vermivorous Conus virgo. We detect multiple transcripts clusters in both the SG and ASG, whose annotations imply venom-related functions. Despite the subsets of transcripts highly-expressed in the VG, SG, and ASG being very distinct, SG expresses an L-, and ASG-Cerm08-, and MEFRR- superfamily conotoxins, all previously considered specific for VG. We corroborate our results with the analysis of published SG and VG transcriptomes from unrelated fish-hunting C. geographus, and C. striatus, possibly fish-hunting C. rolani, and worm-hunting Conus quercinus. In spite of low expression levels of conotoxins, some other specific clusters of putative venom-related peptides are present and may be highly expressed in the SG of these species. Further functional studies are necessary to determine the role that these peptides play in envenomation. In the meantime, our results show importance of routine multi-tissue sampling both for accurate interpretation of tissue-specific venom composition in cone-snails, and for better understanding origin and evolution of venom peptides genes.

Snake venoms are primarily used to subjugate prey, and consequently, their evolution has been shown to be predominantly driven by diet-related selection pressure. Venoms tend to be more lethal to prey than non-prey species (except in cases of toxin resistance), prey-specific toxins have been identified, and preliminary work has demonstrated an association between the diversity of diet classes and that of toxicological activities of whole venom. However, venoms are complex mixtures of many toxins, and it remains unclear how toxin diversity is driven by diet. Prey-specific toxins do not encompass the molecular diversity of venoms, and whole venom effects could be driven by one, few, or all components, so the link between diet and venom diversity remains minimally understood. Here, we collated a database of venom composition and diet records and used a combination of phylogenetic comparative methods and two quantitative diversity indices to investigate whether and how diet diversity relates to the toxin diversity of snake venoms. We reveal that venom diversity is negatively related to diet diversity using Shannon\'s index but positively related using Simpson\'s index. Since Shannon\'s index predominantly considers the number of prey/toxins, whereas Simpson\'s index more strongly reflects evenness, we provide insights into how the diet-venom diversity link is driven. Specifically, species with low diet diversity tend to have venoms dominated by a few abundant (possibly specialised) toxin families, whereas species with diverse diets tend to \'hedge their bets\' by having venoms with a more even composition of different toxin classes.

The cDNAs encoding the Zn+2-metalloproteases (SVMPs) of Trimeresurus gracilis (abbreviated as Tgc), a pitviper endemic to Taiwan, were cloned from venom glands and sequenced. The amino-acid sequences of five novel SVMPs, including one P-III, three P-II and one P-I class enzymes, were thus deduced and subjected to BLAST-analyses. The P-III enzyme (designated as Tgc-PIII) is structurally most similar to the PIII-SVMPs of New World pitvipers but not similar to the PIII-SVMP of Ovophis okinavensis. Sequence-similarity analysis of 22 homologous PIII-SVMPs reveal three major structural subtypes of the pitviper PIII-SVMPs, which possibly have different substrate specificities. In addition, Tgc-PIII and the PI-class SVMP (named Tgc-MP) were isolated from the venom and verified by mass spectrometry. All the three deduced sequences of PII-SVMPs (Tgc-PIIs) contain an abnormal Zn+2-binding-site in their catalytic-domain, and an identical \"long-disintegrin\" domain. The predicted 85-residues disintegrin, gracilisin, bears high similarities to some long-disintegrins of the New-World pitvipers and salmosin3. By BLAST search and comparison, Tgc-MP is 96% similar to okinalysin, the hemorrhagic PI-SVMP of O. okinavensis, rather than any other PI-SVMPs in the databanks. Our results confirm the fast evolution of Tgc-SVMPs as well as their structural similarities to different SVMP-classes of the New-World pitvipers and of O. okinavensis, respectively. The implications of our findings are discussed along with our previous sequence comparisons of venom phospholipases A2 and ten venom serine proteases of Tgc.

Spider venoms are a complex concoction of enzymes, polyamines, inorganic salts, and disulfide-rich peptides (DRPs). Although DRPs are widely distributed and abundant, their bevolutionary origin has remained elusive. This knowledge gap stems from the extensive molecular divergence of DRPs and a lack of sequence and structural data from diverse lineages. By evaluating DRPs under a comprehensive phylogenetic, structural and evolutionary framework, we have not only identified 78 novel spider toxin superfamilies but also provided the first evidence for their common origin. We trace the origin of these toxin superfamilies to a primordial knot - which we name \'Adi Shakti\', after the creator of the Universe according to Hindu mythology - 375 MYA in the common ancestor of Araneomorphae and Mygalomorphae. As the lineages under evaluation constitute nearly 60% of extant spiders, our findings provide fascinating insights into the early evolution and diversification of the spider venom arsenal. Reliance on a single molecular toxin scaffold by nearly all spiders is in complete contrast to most other venomous animals that have recruited into their venoms diverse toxins with independent origins. By comparatively evaluating the molecular evolutionary histories of araneomorph and mygalomorph spider venom toxins, we highlight their contrasting evolutionary diversification rates. Our results also suggest that venom deployment (e.g. prey capture or self-defense) influences evolutionary diversification of DRP toxin superfamilies.
The majority of spiders rely on their venom to defend themselves, to hunt, or both. Armed with this formidable weapon, they have managed to conquer every continent besides Antarctica since they first emerged about 495 million years ago. A closer look at spider venoms hints at an intriguing evolutionary history which has been rarely examined so far. The venom of other animals, such as snakes or scorpions, is usually formed of a wide range of unrelated toxins; in contrast, spiders rely on a single class of proteins, known as disulfide-rich peptides, to create their deadly venom cocktail. This family of molecules is impressively diverse, with each peptide having a distinct structure and mode of action. Its origins, however, have remained elusive. To fill this knowledge gap, Shaikh and Sunagar scanned the sequences of all disulfide-rich peptides generated to date, bringing together a dataset that includes 60% of all modern-day spiders. The analyses allowed the identification of 78 new superfamilies of spider toxins. They also revealed that all existing peptides originate from a single molecule, which Shaikh and Sunagar named after the powerful Hindu goddess Adi Shakti. This ancestral toxin was present 375 million years ago in the last common ancestor of modern-day spiders. The work also highlighted that disulfide-rich peptides evolved under different pressures in various groups of spiders; this may be because some species primarily use their venom for hunting, and others for defence. While the ‘hunters’ may need to constantly acquire toxins with new roles and structures to keep their edge over their prey, those that rely on venom to protect themselves may instead benefit from relying on tried-and-tested toxins useful against a range of infrequent predators. Finally, the analyses revealed that the disulphide-rich peptides of Mygalomorphae tarantulas, which form one of the three major groups of spiders, are much more diverse than the related toxins in other spiders. The underlying reason for this difference is still unclear. Several life-saving drugs currently on the market are based on toxins first identified in the venoms of snakes, cone sails or lizards. Similar discoveries could be unlocked by better understanding the range of deadly molecules used by spiders, and how these came to be.