Convergent evolution of venom gland transcriptomes across Metazoa

Life in the fastlane? A comparative analysis of gene expression profiles across annual, semi-annual, and non-annual killifishes (Cyprinodontiformes: Nothobranchiidae)

The Turquoise Killifish is an important vertebrate for the study of aging and age-related diseases due to its short lifespan. Within Nothobranchiidae, species possess annual, semi-annual, or non-annual life-histories. We took a comparative approach and examined gene expression profiles (QuantSeq) from 62 individuals from eleven nothobranchid species that span three life-histories. Our results show significant differences in differentially expressed genes (DEGs) across life-histories with non-annuals and semi-annuals being most similar, and annuals being the most distinct. At finer scales, we recovered significant differences in DEGs for DNA repair genes and show that non-annual and semi-annuals share similar gene expression profiles, while annuals are distinct. Most of the GO terms enriched in annuals are related to metabolic processes. However, GO terms, including translation, protein transport, and DNA replication initiation also are enriched in annuals. Non-annuals are enriched in Notch signaling pathway genes and downregulated in the canonical Wnt signaling pathway compared to annual species, which suggests that non-annuals have stronger regulation in cellular processes. This study provides support for congruency in DEGs involved in these life-histories and provides strong evidence that a particular set of candidate genes may be worthy of study to investigate their role in the aging process.

A single-cell atlas of the miracidium larva of Schistosoma mansoni reveals cell types, developmental pathways, and tissue architecture

Schistosoma mansoni is a parasitic flatworm that causes the major neglected tropical disease schistosomiasis. The miracidium is the first larval stage of the life cycle. It swims and infects a freshwater snail, transforms into a mother sporocyst, where its stem cells generate daughter sporocysts that give rise to human-infective cercariae larvae. To understand the miracidium at cellular and molecular levels, we created a whole-body atlas of its ~365 cells. Single-cell RNA sequencing identified 19 transcriptionally distinct cell clusters. In situ hybridisation of tissue-specific genes revealed that 93% of the cells in the larva are somatic (57% neural, 19% muscle, 13% epidermal or tegument, 2% parenchyma, and 2% protonephridia) and 7% are stem. Whereas neurons represent the most diverse somatic cell types, trajectory analysis of the two main stem cell populations indicates that one of them is the origin of the tegument lineage and the other likely contains pluripotent cells. Furthermore, unlike the somatic cells, each of these stem populations shows sex-biased transcriptional signatures suggesting a cell-type-specific gene dosage compensation for sex chromosome-linked loci. The miracidium represents a simple developmental stage with which to gain a fundamental understanding of the molecular biology and spatial architecture of schistosome cells.

Host and venom evolution in parasitoid wasps: does independently adapting to the same host shape the evolution of the venom gland transcriptome?

BACKGROUND

Venoms have repeatedly evolved over 100 occasions throughout the animal tree of life, making them excellent systems for exploring convergent evolutionary novelty. Growing evidence supports that venom evolution is predominantly driven by prey or host-related selection pressures, and the expression patterns of venom glands reflect adaptive evolution. However, it remains elusive whether the evolution of expression patterns in venom glands is likewise a convergent evolution driven by their prey/host species.

RESULTS

We utilized parasitoid wasps that had independently adapted to Drosophila hosts as models to investigate the convergent evolution of venom gland transcriptomes in 19 hymenopteran species spanning ~ 200 million years of evolution. Comparative transcriptome analysis reveals that the global expression patterns among the venom glands of Drosophila parasitoid wasps do not achieve higher similarity compared to non-Drosophila parasitoid wasps. Further evolutionary analyses of expression patterns at the single gene, orthogroup, and Gene Ontology (GO) term levels indicate that some orthogroups/GO terms show correlation with the Drosophila parasitoid wasps. However, these groups rarely include genes highly expressed in venom glands or putative venom genes in the Drosophila parasitoid wasps.

CONCLUSIONS

Our study suggests that convergent evolution may not play a predominant force shaping gene expression levels in the venom gland of the Drosophila parasitoid wasps, offering novel insights into the co-evolution between venom and prey/host.

Proteins from shrews' venom glands play a role in gland functioning and venom production

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 gland organogenesis in the common house spider

Venom is a remarkable innovation found across the animal kingdom, yet the evolutionary origins of venom systems in various groups, including spiders, remain enigmatic. Here, we investigated the organogenesis of the venom apparatus in the common house spider, Parasteatoda tepidariorum. The venom apparatus consists of a pair of secretory glands, each connected to an opening at the fang tip by a duct that runs through the chelicerae. We performed bulk RNA-seq to identify venom gland-specific markers and assayed their expression using RNA in situ hybridisation experiments on whole-mount time-series. These revealed that the gland primordium emerges during embryonic stage 13 at the chelicera tip, progresses proximally by the end of embryonic development and extends into the prosoma post-eclosion. The initiation of expression of an important toxin component in late postembryos marks the activation of venom-secreting cells. Our selected markers also exhibited distinct expression patterns in adult venom glands: sage and the toxin marker were expressed in the secretory epithelium, forkhead and sum-1 in the surrounding muscle layer, while Distal-less was predominantly expressed at the gland extremities. Our study provides the first comprehensive analysis of venom gland morphogenesis in spiders, offering key insights into their evolution and development.

The venom and telopodal defence systems of the centipede Lithobius forficatus are functionally convergent serial homologues

BACKGROUND

Evolution of novelty is a central theme in evolutionary biology, yet studying the origins of traits with an apparently discontinuous origin remains a major challenge. Venom systems are a well-suited model for the study of this phenomenon because they capture several aspects of novelty across multiple levels of biological complexity. However, while there is some knowledge on the evolution of individual toxins, not much is known about the evolution of venom systems as a whole. One way of shedding light on the evolution of new traits is to investigate less specialised serial homologues, i.e. repeated traits in an organism that share a developmental origin. This approach can be particularly informative in animals with repetitive body segments, such as centipedes.

RESULTS

Here, we investigate morphological and biochemical aspects of the defensive telopodal glandular organs borne on the posterior legs of venomous stone centipedes (Lithobiomorpha), using a multimethod approach, including behavioural observations, comparative morphology, proteomics, comparative transcriptomics and molecular phylogenetics. We show that the anterior venom system and posterior telopodal defence system are functionally convergent serial homologues, where one (telopodal defence) represents a model for the putative early evolutionary state of the other (venom). Venom glands and telopodal glandular organs appear to have evolved from the same type of epidermal gland (four-cell recto-canal type) and while the telopodal defensive secretion shares a great degree of compositional overlap with centipede venoms in general, these similarities arose predominantly through convergent recruitment of distantly related toxin-like components. Both systems are composed of elements predisposed to functional innovation across levels of biological complexity that range from proteins to glands, demonstrating clear parallels between molecular and morphological traits in the properties that facilitate the evolution of novelty.

CONCLUSIONS

The evolution of the lithobiomorph telopodal defence system provides indirect empirical support for the plausibility of the hypothesised evolutionary origin of the centipede venom system, which occurred through functional innovation and gradual specialisation of existing epidermal glands. Our results thus exemplify how continuous transformation and functional innovation can drive the apparent discontinuous emergence of novelties on higher levels of biological complexity.

Distinct regulatory networks control toxin gene expression in elapid and viperid snakes

BACKGROUND

Venom systems are ideal models to study genetic regulatory mechanisms that underpin evolutionary novelty. Snake venom glands are thought to share a common origin, but there are major distinctions between venom toxins from the medically significant snake families Elapidae and Viperidae, and toxin gene regulatory investigations in elapid snakes have been limited. Here, we used high-throughput RNA-sequencing to profile gene expression and microRNAs between active (milked) and resting (unmilked) venom glands in an elapid (Eastern Brown Snake, Pseudonaja textilis), in addition to comparative genomics, to identify cis- and trans-acting regulation of venom production in an elapid in comparison to viperids (Crotalus viridis and C. tigris).

RESULTS

Although there is conservation in high-level mechanistic pathways regulating venom production (unfolded protein response, Notch signaling and cholesterol homeostasis), there are differences in the regulation of histone methylation enzymes, transcription factors, and microRNAs in venom glands from these two snake families. Histone methyltransferases and transcription factor (TF) specificity protein 1 (Sp1) were highly upregulated in the milked elapid venom gland in comparison to the viperids, whereas nuclear factor I (NFI) TFs were upregulated after viperid venom milking. Sp1 and NFI cis-regulatory elements were common to toxin gene promoter regions, but many unique elements were also present between elapid and viperid toxins. The presence of Sp1 binding sites across multiple elapid toxin gene promoter regions that have been experimentally determined to regulate expression, in addition to upregulation of Sp1 after venom milking, suggests this transcription factor is involved in elapid toxin expression. microRNA profiles were distinctive between milked and unmilked venom glands for both snake families, and microRNAs were predicted to target a diversity of toxin transcripts in the elapid P. textilis venom gland, but only snake venom metalloproteinase transcripts in the viperid C. viridis venom gland. These results suggest differences in toxin gene posttranscriptional regulation between the elapid P. textilis and viperid C. viridis.

CONCLUSIONS

Our comparative transcriptomic and genomic analyses between toxin genes and isoforms in elapid and viperid snakes suggests independent toxin regulation between these two snake families, demonstrating multiple different regulatory mechanisms underpin a venomous phenotype.

Proteotransciptomics of the Most Popular Host Sea Anemone Entacmaea quadricolor Reveals Not All Toxin Genes Expressed by Tentacles Are Recruited into Its Venom Arsenal

While the unique symbiotic relationship between anemonefishes and sea anemones is iconic, it is still not fully understood how anemonefishes can withstand and thrive within the venomous environment of their host sea anemone. In this study, we used a proteotranscriptomics approach to elucidate the proteinaceous toxin repertoire from the most common host sea anemone, Entacmaea quadricolor. Although 1251 different toxin or toxin-like RNA transcripts were expressed in E. quadricolor tentacles (0.05% of gene clusters, 1.8% of expression) and 5375 proteins were detected in milked venom, only 4% of proteins detected in venom were putative toxins (230), and they only represent on average 14% of the normalised protein expression in the milked venom samples. Thus, most proteins in milked venom do not appear to have a toxin function. This work raises the perils of defining a dominant venom phenotype based on transcriptomics data alone in sea anemones, as we found that the dominant venom phenotype differs between the transcriptome and proteome abundance data. E. quadricolor venom contains a mixture of toxin-like proteins of unknown and known function. A newly identified toxin protein family, Z3, rich in conserved cysteines of unknown function, was the most abundant at the RNA transcript and protein levels. The venom was also rich in toxins from the Protease S1, Kunitz-type and PLA2 toxin protein families and contains toxins from eight venom categories. Exploring the intricate venom toxin components in other host sea anemones will be crucial for improving our understanding of how anemonefish adapt to the venomous environment.

ToxCodAn-Genome: an automated pipeline for toxin-gene annotation in genome assembly of venomous lineages

BACKGROUND

The rapid development of sequencing technologies resulted in a wide expansion of genomics studies using venomous lineages. This facilitated research focusing on understanding the evolution of adaptive traits and the search for novel compounds that can be applied in agriculture and medicine. However, the toxin annotation of genomes is a laborious and time-consuming task, and no consensus pipeline is currently available. No computational tool currently exists to address the challenges specific to toxin annotation and to ensure the reproducibility of the process.

RESULTS

Here, we present ToxCodAn-Genome, the first software designed to perform automated toxin annotation in genomes of venomous lineages. This pipeline was designed to retrieve the full-length coding sequences of toxins and to allow the detection of novel truncated paralogs and pseudogenes. We tested ToxCodAn-Genome using 12 genomes of venomous lineages and achieved high performance on recovering their current toxin annotations. This tool can be easily customized to allow improvements in the final toxin annotation set and can be expanded to virtually any venomous lineage. ToxCodAn-Genome is fast, allowing it to run on any personal computer, but it can also be executed in multicore mode, taking advantage of large high-performance servers. In addition, we provide a guide to direct future research in the venomics field to ensure a confident toxin annotation in the genome being studied. As a case study, we sequenced and annotated the toxin repertoire of Bothrops alternatus, which may facilitate future evolutionary and biomedical studies using vipers as models.

CONCLUSIONS

ToxCodAn-Genome is suitable to perform toxin annotation in the genome of venomous species and may help to improve the reproducibility of further studies. ToxCodAn-Genome and the guide are freely available at https://github.com/pedronachtigall/ToxCodAn-Genome.

Web of venom: exploration of big data resources in animal toxin research

Research on animal venoms and their components spans multiple disciplines, including biology, biochemistry, bioinformatics, pharmacology, medicine, and more. Manipulating and analyzing the diverse array of data required for venom research can be challenging, and relevant tools and resources are often dispersed across different online platforms, making them less accessible to nonexperts. In this article, we address the multifaceted needs of the scientific community involved in venom and toxin-related research by identifying and discussing web resources, databases, and tools commonly used in this field. We have compiled these resources into a comprehensive table available on the VenomZone website (https://venomzone.expasy.org/10897). Furthermore, we highlight the challenges currently faced by researchers in accessing and using these resources and emphasize the importance of community-driven interdisciplinary approaches. We conclude by underscoring the significance of enhancing standards, promoting interoperability, and encouraging data and method sharing within the venom research community.

Morphological, functional, compositional and transcriptional constraints shape the distinct venom profiles of the assassin bug Sycanus croceovittatus

Predatory bugs employ a salivary venom apparatus to generate complex venoms for capturing and digesting prey. The venom apparatus consists of different glands for the production of distinct venom sets, but the underlying mechanisms behind this process remain poorly understood. Here we present a comprehensive analysis of the morphological, functional, compositional and transcriptional characteristics of venoms derived from posterior main gland (PMG), anterior main gland (AMG), and accessory gland (AG) of the assassin bug Sycanus croceovittatus. Structural observations revealed the intricate constructions of the venom apparatus, enabling the production and storage of three distinct venom sets in anatomically varied glands and allowing them to be modulated in a context-dependent manner upon utilization. There were remarkable differences in the biological activities exhibited by PMG, AMG, and AG venoms. Proteotranscriptomic analysis demonstrated that these venoms displayed compositional heterogeneity at both the quantity and variety levels of proteins. Transcriptional profiles of the identified venom proteins revealed gland-specific or biased expression patterns. These findings indicate that the divergence in venom profiles among different glands arises from morphological, functional, compositional and transcriptional constraints on the venom apparatus, reflecting remarkable morphogenesis and regulatory gene networks responsible for the compartmentalized production of venom proteins in different glands.

Slowly Making Sense: A Review of the Two-Step Venom System within Slow (Nycticebus spp.) and Pygmy Lorises (Xanthonycticebus spp.)

Since the early 2000s, studies of the evolution of venom within animals have rapidly expanded, offering new revelations on the origins and development of venom within various species. The venomous mammals represent excellent opportunities to study venom evolution due to the varying functional usages, the unusual distribution of venom across unrelated mammals and the diverse variety of delivery systems. A group of mammals that excellently represents a combination of these traits are the slow (Nycticebus spp.) and pygmy lorises (Xanthonycticebus spp.) of south-east Asia, which possess the only confirmed two-step venom system. These taxa also present one of the most intriguing mixes of toxic symptoms (cytotoxicity and immunotoxicity) and functional usages (intraspecific competition and ectoparasitic defence) seen in extant animals. We still lack many pieces of the puzzle in understanding how this venom system works, why it evolved what is involved in the venom system and what triggers the toxic components to work. Here, we review available data building upon a decade of research on this topic, focusing especially on why and how this venom system may have evolved. We discuss that research now suggests that venom in slow lorises has a sophisticated set of multiple uses in both intraspecific competition and the potential to disrupt the immune system of targets; we suggest that an exudate diet reveals several toxic plants consumed by slow and pygmy lorises that could be sequestered into their venom and which may help heal venomous bite wounds; we provide the most up-to-date visual model of the brachial gland exudate secretion protein (BGEsp); and we discuss research on a complement component 1r (C1R) protein in saliva that may solve the mystery of what activates the toxicity of slow and pygmy loris venom. We conclude that the slow and pygmy lorises possess amongst the most complex venom system in extant animals, and while we have still a lot more to understand about their venom system, we are close to a breakthrough, particularly with current technological advances.

Discovering Venom-Derived Drug Candidates Using Differential Gene Expression

Venoms are a diverse and complex group of natural toxins that have been adapted to treat many types of human disease, but rigorous computational approaches for discovering new therapeutic activities are scarce. We have designed and validated a new platform-named VenomSeq-to systematically identify putative associations between venoms and drugs/diseases via high-throughput transcriptomics and perturbational differential gene expression analysis. In this study, we describe the architecture of VenomSeq and its evaluation using the crude venoms from 25 diverse animal species and 9 purified teretoxin peptides. By integrating comparisons to public repositories of differential expression, associations between regulatory networks and disease, and existing knowledge of venom activity, we provide a number of new therapeutic hypotheses linking venoms to human diseases supported by multiple layers of preliminary evidence.

Intersexual Differences in the Gene Expression of Phoneutria depilata (Araneae, Ctenidae) Toxins Revealed by Venom Gland Transcriptome Analyses

The wandering spider, Phoneutria depilata, is one of Colombia's most active nocturnal arthropod predators of vertebrates and invertebrates. Its venom has been a relevant subject of study in the last two decades. However, the scarcity of transcriptomic data for the species limits our knowledge of the distinct components present in its venom for linking the mainly neurotoxic effects of the spider venom to a particular molecular target. The transcriptome of the P. depilata venom gland was analyzed to understand the effect of different diets or sex and the impact of these variables on the composition of the venom. We sequenced venom glands obtained from ten males and ten females from three diet treatments: (i) invertebrate: Tenebrio molitor, (ii) vertebrate: Hemidactylus frenatus, and (iii) mixed (T. molitor + H. frenatus). Of 17,354 assembled transcripts from all samples, 65 transcripts relating to venom production differed between males and females. Among them, 36 were classified as neurotoxins, 14 as serine endopeptidases, 11 as other proteins related to venom production, three as metalloprotease toxins, and one as a venom potentiator. There were no differences in transcripts across the analyzed diets, but when considering the effect of diets on differences between the sexes, 59 transcripts were differentially expressed. Our findings provide essential information on toxins differentially expressed that can be related to sex and the plasticity of the diet of P. depilata and thus can be used as a reference for venomics of other wandering spider species.

Single-Cell Heterogeneity in Snake Venom Expression Is Hardwired by Co-Option of Regulators from Progressively Activated Pathways

The ubiquitous cellular heterogeneity underlying many organism-level phenotypes raises questions about what factors drive this heterogeneity and how these complex heterogeneous systems evolve. Here, we use single-cell expression data from a Prairie rattlesnake (Crotalus viridis) venom gland to evaluate hypotheses for signaling networks underlying snake venom regulation and the degree to which different venom gene families have evolutionarily recruited distinct regulatory architectures. Our findings suggest that snake venom regulatory systems have evolutionarily co-opted trans-regulatory factors from extracellular signal-regulated kinase and unfolded protein response pathways that specifically coordinate expression of distinct venom toxins in a phased sequence across a single population of secretory cells. This pattern of co-option results in extensive cell-to-cell variation in venom gene expression, even between tandemly duplicated paralogs, suggesting this regulatory architecture has evolved to circumvent cellular constraints. While the exact nature of such constraints remains an open question, we propose that such regulatory heterogeneity may circumvent steric constraints on chromatin, cellular physiological constraints (e.g., endoplasmic reticulum stress or negative protein-protein interactions), or a combination of these. Regardless of the precise nature of these constraints, this example suggests that, in some cases, dynamic cellular constraints may impose previously unappreciated secondary constraints on the evolution of gene regulatory networks that favors heterogeneous expression.

Genomic and transcriptomic analyses support a silk gland origin of spider venom glands

BACKGROUND

Spiders comprise a hyperdiverse lineage of predators with venom systems, yet the origin of functionally novel spider venom glands remains unclear. Previous studies have hypothesized that spider venom glands originated from salivary glands or evolved from silk-producing glands present in early chelicerates. However, there is insufficient molecular evidence to indicate similarity among them. Here, we provide comparative analyses of genome and transcriptome data from various lineages of spiders and other arthropods to advance our understanding of spider venom gland evolution.

RESULTS

We generated a chromosome-level genome assembly of a model spider species, the common house spider (Parasteatoda tepidariorum). Module preservation, GO semantic similarity, and differentially upregulated gene similarity analyses demonstrated a lower similarity in gene expressions between the venom glands and salivary glands compared to the silk glands, which questions the validity of the salivary gland origin hypothesis but unexpectedly prefers to support the ancestral silk gland origin hypothesis. The conserved core network in the venom and silk glands was mainly correlated with transcription regulation, protein modification, transport, and signal transduction pathways. At the genetic level, we found that many genes in the venom gland-specific transcription modules show positive selection and upregulated expressions, suggesting that genetic variation plays an important role in the evolution of venom glands.

CONCLUSIONS

This research implies the unique origin and evolutionary path of spider venom glands and provides a basis for understanding the diverse molecular characteristics of venom systems.

Adaptation and evolution of the sea anemone Alvinactis sp. to deep-sea hydrothermal vents: A comparison using transcriptomes

Sea anemones are diverse and ecologically successful members of Anthozoa. They are often found in intertidal and shallow waters, although a few of them inhabit harsher living conditions, such as deep-sea hydrothermal vents. Here, we sequenced the transcriptome of the vent sea anemone Alvinactis sp., which was collected from Edmond vent along the central Indian Ocean ridge at a depth of 3275 m, to explore the molecular mechanisms related to adaptation to vents. Compared with another deep-sea anemone (Paraphelliactis xishaensis) and five shallow water sea anemones, a total of 117 positively selected genes and 46 significantly expanded gene families were found in Alvinactis sp. specifically that may be related to its vent-specific aspect of adaptation. In addition, 127 positively selected genes and 23 significantly expanded gene families that were found in both Alvinactis sp. and P. xishaensis. Among these, vent-specific adaptations of Alvinactis sp. may involve genetic alterations in peroxisome, ubiquitin-mediated protein degradation, oxidative phosphorylation, and endocytosis, and its deep-sea adaptation may involve changes in genetic information processing. Differentially expressed genes between Alvinactis sp. and the deep-sea anemone P. xishaensis were enriched in a variety of pathways related to adaptation, such as energy metabolism, genetic information processing, endocytosis, and peroxisomes. Overall, we provided the first transcriptome of sea anemones that inhabit vents, which enriches our knowledge of deep-sea hydrothermal vent adaptation and the diversity of sea anemones.

Evolution, expression patterns and distribution of novel ribbon worm predatory and defensive toxins

Ribbon worms are active predators that use an eversible proboscis to inject venom into their prey and defend themselves with toxic epidermal secretions. Previous work on nemertean venom has largely focused on just a few species and has not investigated the different predatory and defensive secretions in detail. Consequently, our understanding of the composition and evolution of ribbon worm venoms is still very limited. Here, we present a comparative study of nemertean venom combining RNA-seq differential gene expression analyses of venom-producing tissues, tandem mass spectrometry-based proteomics of toxic secretions, and mass spectrometry imaging of proboscis sections, to shed light onto the composition and evolution of predatory and defensive toxic secretions in Antarctonemertes valida. Our analyses reveal a wide diversity of putative defensive and predatory toxins with tissue-specific gene expression patterns and restricted distributions to the mucus and proboscis proteomes respectively, suggesting that ribbon worms produce distinct toxin cocktails for predation and defense. Our results also highlight the presence of numerous lineage-specific toxins, indicating that venom evolution is highly divergent across nemerteans, producing toxin cocktails that might be finely tuned to subdue different prey. Our data also suggest that the hoplonemertean proboscis is a highly specialized predatory organ that seems to be involved in a variety of biological functions besides predation, including secretion and sensory perception. Overall, our results advance our knowledge into the diversity and evolution of nemertean venoms and highlight the importance of combining different types of data to characterize toxin composition in understudied venomous organisms.