Diversity Begets Diversity When Diet Drives Snake Venom Evolution, but Evenness Rather Than Richness Is What Counts

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 deep-rooted origin of disulfide-rich spider venom toxins

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 Fast and the Furriest: Investigating the Rate of Selection on Mammalian Toxins

The evolution of venom and the selection pressures that act on toxins have been increasingly researched within toxinology in the last two decades, in part due to the exceptionally high rates of diversifying selection observed in animal toxins. In 2015, Sungar and Moran proposed the 'two-speed' model of toxin evolution linking evolutionary age of a group to the rates of selection acting on toxins but due to a lack of data, mammals were not included as less than 30 species of venomous mammal have been recorded, represented by elusive species which produce small amounts of venom. Due to advances in genomics and transcriptomics, the availability of toxin sequences from venomous mammals has been increasing. Using branch- and site-specific selection models, we present the rates of both episodic and pervasive selection acting upon venomous mammal toxins as a group for the first time. We identified seven toxin groups present within venomous mammals, representing Chiroptera, Eulipotyphla and Monotremata: KLK1, Plasminogen Activator, Desmallipins, PACAP, CRiSP, Kunitz Domain One and Kunitz Domain Two. All but one group (KLK1) was identified by our results to be evolving under both episodic and pervasive diversifying selection with four toxin groups having sites that were implicated in the fitness of the animal by TreeSAAP (Selection on Amino Acid Properties). Our results suggest that venomous mammal ecology, behaviour or genomic evolution are the main drivers of selection, although evolutionary age may still be a factor. Our conclusion from these results indicates that mammalian toxins are following the two-speed model of selection, evolving predominately under diversifying selection, fitting in with other younger venomous taxa like snakes and cone snails-with high amounts of accumulating mutations, leading to more novel adaptions in their toxins.

Vexitoxins: conotoxin-like venom peptides from predatory gastropods of the genus Vexillum

Venoms of predatory marine cone snails are intensely studied because of the biomedical applications of the neuropeptides that they contain, termed conotoxins. Meanwhile some gastropod lineages have independently acquired secretory glands strikingly similar to the venom gland of cone snails, suggesting that they possess similar venoms. Here we focus on the most diversified of these clades, the genus Vexillum. Based on the analysis of a multi-species proteo-transcriptomic dataset, we show that Vexillum species indeed produce complex venoms dominated by highly diversified short cysteine-rich peptides, vexitoxins. Vexitoxins possess the same precursor organization, display overlapping cysteine frameworks and share several common post-translational modifications with conotoxins. Some vexitoxins show sequence similarity to conotoxins and adopt similar domain conformations, including a pharmacologically relevant inhibitory cysteine knot motif. The Vexillum envenomation gland (gL) is a notably more recent evolutionary novelty than the conoidean venom gland. Thus, we hypothesize lower divergence between vexitoxin genes, and their ancestral 'somatic' counterparts compared to that in conotoxins, and we find support for this hypothesis in the evolution of the vexitoxin cluster V027. We use this example to discuss how future studies on vexitoxins can inform the origin of conotoxins, and how they may help to address outstanding questions in venom evolution.

Contextual Constraints: Dynamic Evolution of Snake Venom Phospholipase A2

Venom is a dynamic trait that has contributed to the success of numerous organismal lineages. Predominantly composed of proteins, these complex cocktails are deployed for predation and/or self-defence. Many non-toxic physiological proteins have been convergently and recurrently recruited by venomous animals into their toxin arsenal. Phospholipase A₂ (PLA₂) is one such protein and features in the venoms of many organisms across the animal kingdom, including snakes of the families Elapidae and Viperidae. Understanding the evolutionary history of this superfamily would therefore provide insight into the origin and diversification of venom toxins and the evolution of novelty more broadly. The literature is replete with studies that have identified diversifying selection as the sole influence on PLA₂ evolution. However, these studies have largely neglected the structural/functional constraints on PLA₂s, and the ecology and evolutionary histories of the diverse snake lineages that produce them. By considering these crucial factors and employing evolutionary analyses integrated with a schema for the classification of PLA₂s, we uncovered lineage-specific differences in selection regimes. Thus, our work provides novel insights into the evolution of this major snake venom toxin superfamily and underscores the importance of considering the influence of evolutionary and ecological contexts on molecular evolution.

A new, widespread venomous mammal species: hemolytic activity of Sorex araneus venom is similar to that of Neomys fodiens venom

BACKGROUND

Venom production has evolved independently many times in the animal kingdom, although it is rare among mammals. Venomous shrews produce toxins in their salivary glands and use their venoms to hunt and store prey. Thus far, the toxicity and composition of shrew venoms have been studied only in two shrew species: the northern short-tailed shrew, Blarina brevicauda, and the Eurasian water shrew, Neomys fodiens. Venom of N. fodiens has potent paralytic activity which enables hunting and storing prey in a comatose state. Here, we assayed the hemolytic effects of extracts from salivary glands of N. fodiens and the common shrew, Sorex araneus, in erythrocytes of Pelophylax sp. frogs. We identified toxins in shrew venom by high-performance liquid chromatography coupled to tandem mass spectrometry.

RESULTS

Our results prove, confirming a suggestion made four centuries ago, that S. araneus is venomous. We also provide the first experimental evidence that shrew venoms produce potent hemolysis in frog erythrocytes. We found significant concentration-dependent effects of venoms of N. fodiens and S. araneus on hemolysis of red blood cells evaluated as hemoglobin release. Treatment of erythrocytes with N. fodiens venom at concentrations of 1.0 and 0.5 mg/ml and with S. araneus venom at concentration of 1.0 mg/ml caused an increased release of hemoglobin. Our findings confirm that hemolytic effects of N. fodiens venom are stronger than those produced by S. araneus venom. We identified four toxins in the venom of N. fodiens: proenkephalin, phospholipase A₂ (PLA₂), a disintegrin and metalloproteinase domain-containing protein (ADAM) and lysozyme C, as well as a non-toxic hyaluronidase. In the venom of S. araneus we found five toxins: proenkephalin, kallikrein 1-related peptidase, beta-defensin, ADAM and lysozyme C. PLA₂ and ADAMs are likely to produce hemolysis in frog erythrocytes.

CONCLUSIONS

Our results clearly show that shrew venoms possess hemolytic action that may allow them to hunt larger prey. Since a member of the numerous genus Sorex is venomous, it is likely that venom production among shrews and other eulipotyphlans may be more widespread than it has previously been assumed.

The rise of genomics in snake venom research: recent advances and future perspectives

Snake venoms represent a danger to human health, but also a gold mine of bioactive proteins that can be harnessed for drug discovery purposes. The evolution of snakes and their venom has been studied for decades, particularly via traditional morphological and basic genetic methods alongside venom proteomics. However, while the field of genomics has matured rapidly over the past 2 decades, owing to the development of next-generation sequencing technologies, snake genomics remains in its infancy. Here, we provide an overview of the state of the art in snake genomics and discuss its potential implications for studying venom evolution and toxinology. On the basis of current knowledge, gene duplication and positive selection are key mechanisms in the neofunctionalization of snake venom proteins. This makes snake venoms important evolutionary drivers that explain the remarkable venom diversification and adaptive variation observed in these reptiles. Gene duplication and neofunctionalization have also generated a large number of repeat sequences in snake genomes that pose a significant challenge to DNA sequencing, resulting in the need for substantial computational resources and longer sequencing read length for high-quality genome assembly. Fortunately, owing to constantly improving sequencing technologies and computational tools, we are now able to explore the molecular mechanisms of snake venom evolution in unprecedented detail. Such novel insights have the potential to affect the design and development of antivenoms and possibly other drugs, as well as provide new fundamental knowledge on snake biology and evolution.

A phylogeny-aware approach reveals unexpected venom components in divergent lineages of cone snails

Marine gastropods of the genus Conus are renowned for their remarkable diversity and deadly venoms. While Conus venoms are increasingly well studied for their biomedical applications, we know surprisingly little about venom composition in other lineages of Conidae. We performed comprehensive venom transcriptomic profiling for Conasprella coriolisi and Pygmaeconus traillii, first time for both respective genera. We complemented reference-based transcriptome annotation by a de novo toxin prediction guided by phylogeny, which involved transcriptomic data on two additional 'divergent' cone snail lineages, Profundiconus, and Californiconus. We identified toxin clusters (SSCs) shared among all or some of the four analysed genera based on the identity of the signal region-a molecular tag present in toxins. In total, 116 and 98 putative toxins represent 29 and 28 toxin gene superfamilies in Conasprella and Pygmaeconus, respectively; about quarter of these only found by semi-manual annotation of the SSCs. Two rare gene superfamilies, originally identified from fish-hunting cone snails, were detected outside Conus rather unexpectedly, so we further investigated their distribution across Conidae radiation. We demonstrate that both these, in fact, are ubiquitous in Conidae, sometimes with extremely high expression. Our findings demonstrate how a phylogeny-aware approach circumvents methodological caveats of similarity-based transcriptome annotation.

Identification of Novel Toxin Genes from the Stinging Nettle Caterpillar Parasa lepida (Cramer, 1799): Insights into the Evolution of Lepidoptera Toxins

Simple Summary

Many caterpillar species can produce toxins that cause harmful reactions to humans, varying from mild irritation to death. Currently, there is very limited knowledge about caterpillar toxin diversity, because only a few species have been investigated. We used the transcriptome technique to identify candidate toxin genes from the nettle caterpillar Parasa lepida (Cramer, 1799). It is a common pest of oil palm, coconut, and mango in South and South-East Asia, which can cause severe pain and allergic responses to those in contact with them. We reported 168 candidate toxin genes. Most of them are members of the toxin genes families commonly recruited in animal venoms such as serine protease and serine protease inhibitors. However, we identified 21 novel genes encoding knottin-like peptides expressed at a high level in the transcriptome. Their predicted 3D structures are similar to neurotoxins in scorpion and tarantula. Our study suggests that P. lepida venom contains diverse toxin proteins that potentially cause allergic reactions and pain. This study sheds light on the hidden diversity of toxin proteins in caterpillar lineage, which could be future fruitful new drug sources.

Abstract

Many animal species can produce venom for defense, predation, and competition. The venom usually contains diverse peptide and protein toxins, including neurotoxins, proteolytic enzymes, protease inhibitors, and allergens. Some drugs for cancer, neurological disorders, and analgesics were developed based on animal toxin structures and functions. Several caterpillar species possess venoms that cause varying effects on humans both locally and systemically. However, toxins from only a few species have been investigated, limiting the full understanding of the Lepidoptera toxin diversity and evolution. We used the RNA-seq technique to identify toxin genes from the stinging nettle caterpillar, Parasa lepida (Cramer, 1799). We constructed a transcriptome from caterpillar urticating hairs and reported 34,968 unique transcripts. Using our toxin gene annotation pipeline, we identified 168 candidate toxin genes, including protease inhibitors, proteolytic enzymes, and allergens. The 21 P. lepida novel Knottin-like peptides, which do not show sequence similarity to any known peptide, have predicted 3D structures similar to tarantula, scorpion, and cone snail neurotoxins. We highlighted the importance of convergent evolution in the Lepidoptera toxin evolution and the possible mechanisms. This study opens a new path to understanding the hidden diversity of Lepidoptera toxins, which could be a fruitful source for developing new drugs.

The Curious Case of the "Neurotoxic Skink": Scientific Literature Points to the Absence of Venom in Scincidae

In contrast to the clearly documented evolution of venom in many animal lineages, the origin of reptilian venom is highly debated. Historically, venom has been theorised to have evolved independently in snakes and lizards. However, some of the recent works have argued for the common origin of venom in “Toxicofera” reptiles, which include the order Serpentes (all snakes), and Anguimorpha and Iguania lizards. Nevertheless, in both these contrasting hypotheses, the lizards of the family Scincidae are considered to be harmless and devoid of toxic venoms. Interestingly, an unusual clinical case claiming neurotoxic envenoming by a scincid lizard was recently reported in Southern India. Considering its potentially significant medicolegal, conservation and evolutionary implications, we have summarised the scientific evidence that questions the validity of this clinical report. We argue that the symptoms documented in the patient are likely to have resulted from krait envenomation, which is far too frequent in these regions.

A Wolf in Another Wolf's Clothing: Post-Genomic Regulation Dictates Venom Profiles of Medically-Important Cryptic Kraits in India

The Common Krait (Bungarus caeruleus) shares a distribution range with many other ‘phenotypically-similar’ kraits across the Indian subcontinent. Despite several reports of fatal envenomings by other Bungarus species, commercial Indian antivenoms are only manufactured against B. caeruleus. It is, therefore, imperative to understand the distribution of genetically distinct lineages of kraits, the compositional differences in their venoms, and the consequent impact of venom variation on the (pre)clinical effectiveness of antivenom therapy. To address this knowledge gap, we conducted phylogenetic and comparative venomics investigations of kraits in Southern and Western India. Phylogenetic reconstructions using mitochondrial markers revealed a new species of krait, Romulus’ krait (Bungarus romulusi sp. nov.), in Southern India. Additionally, we found that kraits with 17 mid-body dorsal scale rows in Western India do not represent a subspecies of the Sind Krait (B. sindanus walli) as previously believed, but are genetically very similar to B. sindanus in Pakistan. Furthermore, venom proteomics and comparative transcriptomics revealed completely contrasting venom profiles. While the venom gland transcriptomes of all three species were highly similar, venom proteomes and toxicity profiles differed significantly, suggesting the prominent role of post-genomic regulatory mechanisms in shaping the venoms of these cryptic kraits. In vitro venom recognition and in vivo neutralisation experiments revealed a strong negative impact of venom variability on the preclinical performance of commercial antivenoms. While the venom of B. caeruleus was neutralised as per the manufacturer’s claim, performance against the venoms of B. sindanus and B. romulusi was poor, highlighting the need for regionally-effective antivenoms in India.

Adaptive venom evolution and toxicity in octopods is driven by extensive novel gene formation, expansion, and loss

BACKGROUND

Cephalopods represent a rich system for investigating the genetic basis underlying organismal novelties. This diverse group of specialized predators has evolved many adaptations including proteinaceous venom. Of particular interest is the blue-ringed octopus genus (Hapalochlaena), which are the only octopods known to store large quantities of the potent neurotoxin, tetrodotoxin, within their tissues and venom gland.

FINDINGS

To reveal genomic correlates of organismal novelties, we conducted a comparative study of 3 octopod genomes, including the Southern blue-ringed octopus (Hapalochlaena maculosa). We present the genome of this species and reveal highly dynamic evolutionary patterns at both non-coding and coding organizational levels. Gene family expansions previously reported in Octopus bimaculoides (e.g., zinc finger and cadherins, both associated with neural functions), as well as formation of novel gene families, dominate the genomic landscape in all octopods. Examination of tissue-specific genes in the posterior salivary gland revealed that expression was dominated by serine proteases in non-tetrodotoxin-bearing octopods, while this family was a minor component in H. maculosa. Moreover, voltage-gated sodium channels in H. maculosa contain a resistance mutation found in pufferfish and garter snakes, which is exclusive to the genus. Analysis of the posterior salivary gland microbiome revealed a diverse array of bacterial species, including genera that can produce tetrodotoxin, suggestive of a possible production source.

CONCLUSIONS

We present the first tetrodotoxin-bearing octopod genome H. maculosa, which displays lineage-specific adaptations to tetrodotoxin acquisition. This genome, along with other recently published cephalopod genomes, represents a valuable resource from which future work could advance our understanding of the evolution of genomic novelty in this family.

Toxin-like neuropeptides in the sea anemone Nematostella unravel recruitment from the nervous system to venom

The sea anemone Nematostella vectensis (Anthozoa, Cnidaria) is a powerful model for characterizing the evolution of genes functioning in venom and nervous systems. Although venom has evolved independently numerous times in animals, the evolutionary origin of many toxins remains unknown. In this work, we pinpoint an ancestral gene giving rise to a new toxin and functionally characterize both genes in the same species. Thus, we report a case of protein recruitment from the cnidarian nervous to venom system. The ShK-like1 peptide has a ShKT cysteine motif, is lethal for fish larvae and packaged into nematocysts, the cnidarian venom-producing stinging capsules. Thus, ShK-like1 is a toxic venom component. Its paralog, ShK-like2, is a neuropeptide localized to neurons and is involved in development. Both peptides exhibit similarities in their functional activities: They provoke contraction in Nematostella polyps and are toxic to fish. Because ShK-like2 but not ShK-like1 is conserved throughout sea anemone phylogeny, we conclude that the two paralogs originated due to a Nematostella-specific duplication of a ShK-like2 ancestor, a neuropeptide-encoding gene, followed by diversification and partial functional specialization. ShK-like2 is represented by two gene isoforms controlled by alternative promoters conferring regulatory flexibility throughout development. Additionally, we characterized the expression patterns of four other peptides with structural similarities to studied venom components and revealed their unexpected neuronal localization. Thus, we employed genomics, transcriptomics, and functional approaches to reveal one venom component, five neuropeptides with two different cysteine motifs, and an evolutionary pathway from nervous to venom system in Cnidaria.

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.

Causes and Consequences of Snake Venom Variation

Snake venoms are mixtures of toxins that vary extensively between and within snake species. This variability has serious consequences for the management of the world’s 1.8 million annual snakebite victims. Advances in ‘omic’ technologies have empowered toxinologists to comprehensively characterize snake venom compositions, unravel the molecular mechanisms that underpin venom variation, and elucidate the ensuing functional consequences. In this review, we describe how such mechanistic processes have resulted in suites of toxin isoforms that cause diverse pathologies in human snakebite victims and we detail how variation in venom composition can result in treatment failure. Finally, we outline current therapeutic approaches designed to circumvent venom variation and deliver next-generation treatments for the world’s most lethal neglected tropical disease.

Structural venomics reveals evolution of a complex venom by duplication and diversification of an ancient peptide-encoding gene

Spiders are one of the most successful venomous animals, with more than 48,000 described species. Most spider venoms are dominated by cysteine-rich peptides with a diverse range of pharmacological activities. Some spider venoms contain thousands of unique peptides, but little is known about the mechanisms used to generate such complex chemical arsenals. We used an integrated transcriptomic, proteomic, and structural biology approach to demonstrate that the lethal Australian funnel-web spider produces 33 superfamilies of venom peptides and proteins. Twenty-six of the 33 superfamilies are disulfide-rich peptides, and we show that 15 of these are knottins that contribute >90% of the venom proteome. NMR analyses revealed that most of these disulfide-rich peptides are structurally related and range in complexity from simple to highly elaborated knottin domains, as well as double-knot toxins, that likely evolved from a single ancestral toxin gene.

Lack of Signal for the Impact of Conotoxin Gene Diversity on Speciation Rates in Cone Snails

Understanding why some groups of organisms are more diverse than others is a central goal in macroevolution. Evolvability, or the intrinsic capacity of lineages for evolutionary change, is thought to influence disparities in species diversity across taxa. Over macroevolutionary time scales, clades that exhibit high evolvability are expected to have higher speciation rates. Cone snails (family: Conidae, $>$900 spp.) provide a unique opportunity to test this prediction because their toxin genes can be used to characterize differences in evolvability between clades. Cone snails are carnivorous, use prey-specific venom (conotoxins) to capture prey, and the genes that encode venom are known and diversify through gene duplication. Theory predicts that higher gene diversity confers a greater potential to generate novel phenotypes for specialization and adaptation. Therefore, if conotoxin gene diversity gives rise to varying levels of evolvability, conotoxin gene diversity should be coupled with macroevolutionary speciation rates. We applied exon capture techniques to recover phylogenetic markers and conotoxin loci across 314 species, the largest venom discovery effort in a single study. We paired a reconstructed timetree using 12 fossil calibrations with species-specific estimates of conotoxin gene diversity and used trait-dependent diversification methods to test the impact of evolvability on diversification patterns. Surprisingly, we did not detect any signal for the relationship between conotoxin gene diversity and speciation rates, suggesting that venom evolution may not be the rate-limiting factor controlling diversification dynamics in Conidae. Comparative analyses showed some signal for the impact of diet and larval dispersal strategy on diversification patterns, though detection of a signal depended on the dataset and the method. If our results remain true with increased taxonomic sampling in future studies, they suggest that the rapid evolution of conid venom may cause other factors to become more critical to diversification, such as ecological opportunity or traits that promote isolation among lineages.

Toxins from scratch? Diverse, multimodal gene origins in the predatory robber fly Dasypogon diadema indicate a dynamic venom evolution in dipteran insects

BACKGROUND

Venoms and the toxins they contain represent molecular adaptations that have evolved on numerous occasions throughout the animal kingdom. However, the processes that shape venom protein evolution are poorly understood because of the scarcity of whole-genome data available for comparative analyses of venomous species.

RESULTS

We performed a broad comparative toxicogenomic analysis to gain insight into the genomic mechanisms of venom evolution in robber flies (Asilidae). We first sequenced a high-quality draft genome of the hymenopteran hunting robber fly Dasypogon diadema, analysed its venom by a combined proteotranscriptomic approach, and compared our results with recently described robber fly venoms to assess the general composition and major components of asilid venom. We then applied a comparative genomics approach, based on 1 additional asilid genome, 10 high-quality dipteran genomes, and 2 lepidopteran outgroup genomes, to reveal the evolutionary mechanisms and origins of identified venom proteins in robber flies.

CONCLUSIONS

While homologues were identified for 15 of 30 predominant venom protein in the non-asilid genomes, the remaining 15 highly expressed venom proteins appear to be unique to robber flies. Our results reveal that the venom of D. diadema likely evolves in a multimodal fashion comprising (i) neofunctionalization after gene duplication, (ii) expression-dependent co-option of proteins, and (iii) asilid lineage-specific orphan genes with enigmatic origin. The role of such orphan genes is currently being disputed in evolutionary genomics but has not been discussed in the context of toxin evolution. Our results display an unexpected dynamic venom evolution in asilid insects, which contrasts the findings of the only other insect toxicogenomic evolutionary analysis, in parasitoid wasps (Hymenoptera), where toxin evolution is dominated by single gene co-option. These findings underpin the significance of further genomic studies to cover more neglected lineages of venomous taxa and to understand the importance of orphan genes as possible drivers for venom evolution.

Dynamics of venom composition across a complex life cycle

Little is known about venom in young developmental stages of animals. The appearance of toxins and stinging cells during early embryonic stages in the sea anemone Nematostella vectensis suggests that venom is already expressed in eggs and larvae of this species. Here, we harness transcriptomic, biochemical and transgenic tools to study venom production dynamics in Nematostella. We find that venom composition and arsenal of toxin-producing cells change dramatically between developmental stages of this species. These findings can be explained by the vastly different interspecific interactions of each life stage, as individuals develop from a miniature non-feeding mobile planula to a larger sessile polyp that predates on other animals and interact differently with predators. Indeed, behavioral assays involving prey, predators and Nematostella are consistent with this hypothesis. Further, the results of this work suggest a much wider and dynamic venom landscape than initially appreciated in animals with a complex life cycle.

Some animals produce a mixture of toxins, commonly known as venom, to protect themselves from predators and catch prey. Cnidarians – a group of animals that includes sea anemones, jellyfish and corals – have stinging cells on their tentacles that inject venom into the animals they touch.

The sea anemone Nematostella goes through a complex life cycle. Nematostella start out life in eggs. They then become swimming larvae, barely visible to the naked eye, that do not feed. Adult Nematostella are cylindrical, stationary ‘polyps’ that are several inches long. They use tentacles at the end of their tube-like bodies to capture small aquatic animals. Sea anemones therefore change how they interact with predators and prey at different stages of their life. Most research on venomous animals focuses on adults, so until now it was not clear whether the venom changes along their maturation.

Columbus-Shenkar, Sachkova et al. genetically modified Nematostella so that the cells that produce distinct venom components were labeled with different fluorescent markers. The composition of the venom could then be linked to how the anemones interacted with their fish and shrimp predators at each life stage.

The results of the experiments showed that Nematostella mothers pass on a toxin to their eggs that makes them unpalatable to predators. Larvae then produce high levels of other toxins that allow them to incapacitate or kill potential predators. Adults have a different mix of toxins that likely help them capture prey.

Venom is often studied because the compounds it contains have the potential to be developed into new drugs. The jellyfish and coral relatives of Nematostella may also produce different venoms at different life stages. This means that there are likely to be many toxins that we have not yet identified in these animals. As some jellyfish venoms are very active on humans and reef corals have a pivotal role in ocean ecology, further research into the venoms produced at different life stages could help us to understand and preserve marine ecosystems, as well as having medical benefits.

Divergence of the Venom Exogene Repertoire in Two Sister Species of Turriconus

The genus Conus comprises approximately 700 species of venomous marine cone snails that are highly efficient predators of worms, snails, and fish. In evolutionary terms, cone snails are relatively young with the earliest fossil records occurring in the Lower Eocene, 55 Ma. The rapid radiation of cone snail species has been accompanied by remarkably high rates of toxin diversification. To shed light on the molecular mechanisms that accompany speciation, we investigated the toxin repertoire of two sister species, Conus andremenezi and Conus praecellens, that were until recently considered a single variable species. A total of 196 and 250 toxin sequences were identified in the venom gland transcriptomes of C. andremenezi and C. praecellens belonging to 25 and 29 putative toxin gene superfamilies, respectively. Comparative analysis with closely (Conus tribblei and Conus lenavati) and more distantly related species (Conus geographus) suggests that speciation is associated with significant diversification of individual toxin genes (exogenes) whereas the expression pattern of toxin gene superfamilies within lineages remains largely conserved. Thus, changes within individual toxin sequences can serve as a sensitive indicator for recent speciation whereas changes in the expression pattern of gene superfamilies are likely to reflect more dramatic differences in a species' interaction with its prey, predators, and competitors.

Molecular Diversity and Gene Evolution of the Venom Arsenal of Terebridae Predatory Marine Snails

Venom peptides from predatory organisms are a resource for investigating evolutionary processes such as adaptive radiation or diversification, and exemplify promising targets for biomedical drug development. Terebridae are an understudied lineage of conoidean snails, which also includes cone snails and turrids. Characterization of cone snail venom peptides, conotoxins, has revealed a cocktail of bioactive compounds used to investigate physiological cellular function, predator-prey interactions, and to develop novel therapeutics. However, venom diversity of other conoidean snails remains poorly understood. The present research applies a venomics approach to characterize novel terebrid venom peptides, teretoxins, from the venom gland transcriptomes of Triplostephanus anilis and Terebra subulata. Next-generation sequencing and de novo assembly identified 139 putative teretoxins that were analyzed for the presence of canonical peptide features as identified in conotoxins. To meet the challenges of de novo assembly, multiple approaches for cross validation of findings were performed to achieve reliable assemblies of venom duct transcriptomes and to obtain a robust portrait of Terebridae venom. Phylogenetic methodology was used to identify 14 teretoxin gene superfamilies for the first time, 13 of which are unique to the Terebridae. Additionally, basic local algorithm search tool homology-based searches to venom-related genes and posttranslational modification enzymes identified a convergence of certain venom proteins, such as actinoporin, commonly found in venoms. This research provides novel insights into venom evolution and recruitment in Conoidean predatory marine snails and identifies a plethora of terebrid venom peptides that can be used to investigate fundamental questions pertaining to gene evolution.

Production and packaging of a biological arsenal: evolution of centipede venoms under morphological constraint

Venom represents one of the most extreme manifestations of a chemical arms race. Venoms are complex biochemical arsenals, often containing hundreds to thousands of unique protein toxins. Despite their utility for prey capture, venoms are energetically expensive commodities, and consequently it is hypothesized that venom complexity is inversely related to the capacity of a venomous animal to physically subdue prey. Centipedes, one of the oldest yet least-studied venomous lineages, appear to defy this rule. Although scutigeromorph centipedes produce less complex venom than those secreted by scolopendrid centipedes, they appear to rely heavily on venom for prey capture. We show that the venom glands are large and well developed in both scutigerid and scolopendrid species, but that scutigerid forcipules lack the adaptations that allow scolopendrids to inflict physical damage on prey and predators. Moreover, we reveal that scolopendrid venom glands have evolved to accommodate a much larger number of secretory cells and, by using imaging mass spectrometry, we demonstrate that toxin production is heterogeneous across these secretory units. We propose that the differences in venom complexity between centipede orders are largely a result of morphological restrictions of the venom gland, and consequently there is a strong correlation between the morphological and biochemical complexity of this unique venom system. The current data add to the growing body of evidence that toxins are not expressed in a spatially homogenous manner within venom glands, and they suggest that the link between ecology and toxin evolution is more complex than previously thought.

A Polychaete's powerful punch: venom gland transcriptomics of Glycera reveals a complex cocktail of toxin homologs

Glycerids are marine annelids commonly known as bloodworms. Bloodworms have an eversible proboscis adorned with jaws connected to venom glands. Bloodworms prey on invertebrates, and it is known that the venom glands produce compounds that can induce toxic effects in animals. Yet, none of these putative toxins has been characterized on a molecular basis. Here we present the transcriptomic profiles of the venom glands of three species of bloodworm, Glycera dibranchiata, Glycera fallax and Glycera tridactyla, as well as the body tissue of G. tridactyla. The venom glands express a complex mixture of transcripts coding for putative toxin precursors. These transcripts represent 20 known toxin classes that have been convergently recruited into animal venoms, as well as transcripts potentially coding for Glycera-specific toxins. The toxins represent five functional categories: Pore-forming and membrane-disrupting toxins, neurotoxins, protease inhibitors, other enzymes, and CAP domain toxins. Many of the transcripts coding for putative Glycera toxins belong to classes that have been widely recruited into venoms, but some are homologs of toxins previously only known from the venoms of scorpaeniform fish and monotremes (stonustoxin-like toxin), turrid gastropods (turripeptide-like peptides), and sea anemones (gigantoxin I-like neurotoxin). This complex mixture of toxin homologs suggests that bloodworms employ venom while predating on macroscopic prey, casting doubt on the previously widespread opinion that G. dibranchiata is a detritivore. Our results further show that researchers should be aware that different assembly methods, as well as different methods of homology prediction, can influence the transcriptomic profiling of venom glands.

The first venomous crustacean revealed by transcriptomics and functional morphology: remipede venom glands express a unique toxin cocktail dominated by enzymes and a neurotoxin

Animal venoms have evolved many times. Venomous species are especially common in three of the four main groups of arthropods (Chelicerata, Myriapoda, and Hexapoda), which together represent tens of thousands of species of venomous spiders, scorpions, centipedes, and hymenopterans. Surprisingly, despite their great diversity of body plans, there is no unambiguous evidence that any crustacean is venomous. We provide the first conclusive evidence that the aquatic, blind, and cave-dwelling remipede crustaceans are venomous and that venoms evolved in all four major arthropod groups. We produced a three-dimensional reconstruction of the venom delivery apparatus of the remipede Speleonectes tulumensis, showing that remipedes can inject venom in a controlled manner. A transcriptomic profile of its venom glands shows that they express a unique cocktail of transcripts coding for known venom toxins, including a diversity of enzymes and a probable paralytic neurotoxin very similar to one described from spider venom. We screened a transcriptomic library obtained from whole animals and identified a nontoxin paralog of the remipede neurotoxin that is not expressed in the venom glands. This allowed us to reconstruct its probable evolutionary origin and underlines the importance of incorporating data derived from nonvenom gland tissue to elucidate the evolution of candidate venom proteins. This first glimpse into the venom of a crustacean and primitively aquatic arthropod reveals conspicuous differences from the venoms of other predatory arthropods such as centipedes, scorpions, and spiders and contributes valuable information for ultimately disentangling the many factors shaping the biology and evolution of venoms and venomous species.