A process of convergent amplification and tissue‐specific expression dominates the evolution of toxin and toxin‐like genes in sea anemones

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.

Diversity analysis of sea anemone peptide toxins in different tissues of Heteractis crispa based on transcriptomics

Peptide toxins found in sea anemones venom have diverse properties that make them important research subjects in the fields of pharmacology, neuroscience and biotechnology. This study used high-throughput sequencing technology to systematically analyze the venom components of the tentacles, column, and mesenterial filaments of sea anemone Heteractis crispa, revealing the diversity and complexity of sea anemone toxins in different tissues. A total of 1049 transcripts were identified and categorized into 60 families, of which 91.0% were proteins and 9.0% were peptides. Of those 1049 transcripts, 416, 291, and 307 putative proteins and peptide precursors were identified from tentacles, column, and mesenterial filaments respectively, while 428 were identified when the datasets were combined. Of these putative toxin sequences, 42 were detected in all three tissues, including 33 proteins and 9 peptides, with the majority of peptides being ShKT domain, β-defensin, and Kunitz-type. In addition, this study applied bioinformatics approaches to predict the family classification, 3D structures, and functional annotation of these representative peptides, as well as the evolutionary relationships between peptides, laying the foundation for the next step of peptide pharmacological activity research.

Evolutionary Analysis of Cnidaria Small Cysteine-Rich Proteins (SCRiPs), an Enigmatic Neurotoxin Family from Stony Corals and Sea Anemones (Anthozoa: Hexacorallia)

Cnidarians (corals, sea anemones, and jellyfish) produce toxins that play central roles in key ecological processes, including predation, defense, and competition, being the oldest extant venomous animal lineage. Cnidaria small cysteine-rich proteins (SCRiPs) were the first family of neurotoxins detected in stony corals, one of the ocean's most crucial foundation species. Yet, their molecular evolution remains poorly understood. Moreover, the lack of a clear classification system has hindered the establishment of an accurate and phylogenetically informed nomenclature. In this study, we extensively surveyed 117 genomes and 103 transcriptomes of cnidarians to identify orthologous SCRiP gene sequences. We annotated a total of 168 novel putative SCRiPs from over 36 species of stony corals and 12 species of sea anemones. Phylogenetic reconstruction identified four distinct SCRiP subfamilies, according to strict discrimination criteria based on well-supported monophyly with a high percentage of nucleotide and amino acids' identity. Although there is a high prevalence of purifying selection for most SCRiP subfamilies, with few positively selected sites detected, a subset of Acroporidae sequences is influenced by diversifying positive selection, suggesting potential neofunctionalizations related to the fine-tuning of toxin potency. We propose a new nomenclature classification system relying on the phylogenetic distribution and evolution of SCRiPs across Anthozoa, which will further assist future proteomic and functional research efforts.

Structure-function relationships in ShKT domain peptides: ShKT-Ts1 from the sea anemone Telmatactis stephensoni

Diverse structural scaffolds have been described in peptides from sea anemones, with the ShKT domain being a common scaffold first identified in ShK toxin from Stichodactyla helianthus. ShK is a potent blocker of voltage-gated potassium channels (KV 1.x), and an analog, ShK-186 (dalazatide), has completed Phase 1 clinical trials in plaque psoriasis. The ShKT domain has been found in numerous other species, but only a tiny fraction of ShKT domains has been characterized functionally. Despite adopting the canonical ShK fold, some ShKT peptides from sea anemones inhibit KV 1.x, while others do not. Mutagenesis studies have shown that a Lys-Tyr (KY) dyad plays a key role in KV 1.x blockade, although a cationic residue followed by a hydrophobic residue may also suffice. Nevertheless, ShKT peptides displaying an ShK-like fold and containing a KY dyad do not necessarily block potassium channels, so additional criteria are needed to determine whether new ShKT peptides might show activity against potassium channels. In this study, we used a combination of NMR and molecular dynamics (MD) simulations to assess the potential activity of a new ShKT peptide. We determined the structure of ShKT-Ts1, from the sea anemone Telmatactis stephensoni, examined its tissue localization, and investigated its activity against a range of ion channels. As ShKT-Ts1 showed no activity against KV 1.x channels, we used MD simulations to investigate whether solvent exposure of the dyad residues may be informative in rationalizing and potentially predicting the ability of ShKT peptides to block KV 1.x channels. We show that either a buried dyad that does not become exposed during MD simulations, or a partially exposed dyad that becomes buried during MD simulations, correlates with weak or absent activity against KV 1.x channels. Therefore, structure determination coupled with MD simulations, may be used to predict whether new sequences belonging to the ShKT family may act as potassium channel blockers.

Venomics Reveals the Venom Complexity of Sea Anemone Heteractis magnifica

The venoms of various sea anemones are rich in diverse toxins, which usually play a dual role in capturing prey and deterring predators. However, the complex components of such venoms have not been well known yet. Here, venomics of integrating transcriptomic and proteomic technologies was applied for the first time to identify putative protein and peptide toxins from different tissues of the representative sea anemone, Heteractis magnifica. The transcriptomic analysis of H. magnifica identified 728 putative toxin sequences, including 442 and 381 from the tentacles and the column, respectively, and they were assigned to 68 gene superfamilies. The proteomic analysis confirmed 101 protein and peptide toxins in the venom, including 91 in the tentacles and 39 in the column. The integrated venomics also confirmed that some toxins such as the ShK-like peptides and defensins are co-expressed in both the tentacles and the column. Meanwhile, a homology analysis was conducted to predict the three-dimensional structures and potential activity of seven representative toxins. Altogether, this venomics study revealed the venom complexity of H. magnifica, which will help deepen our understanding of cnidarian toxins, thereby supporting the in-depth development of valuable marine drugs.

Structural and functional characterisation of Tst2, a novel TRPV1 inhibitory peptide from the Australian sea anemone Telmatactis stephensoni

Sea anemone venoms are complex mixtures of biologically active compounds, including disulfide-rich peptides, some of which have found applications as research tools, and others as therapeutic leads. Our recent transcriptomic and proteomic studies of the Australian sea anemone Telmatactis stephensoni identified a transcript for a peptide designated Tst2. Tst2 is a 38-residue peptide showing sequence similarity to peptide toxins known to interact with a range of ion channels (NaV, TRPV1, KV and CaV). Recombinant Tst2 (rTst2, which contains an additional Gly at the N-terminus) was produced by periplasmic expression in Escherichia coli, enabling the production of both unlabelled and uniformly 13C,15N-labelled peptide for functional assays and structural studies. The LC-MS profile of the recombinant Tst2 showed a pure peak with molecular mass 6 Da less than that of the reduced form of the peptide, indicating the successful formation of three disulfide bonds from its six cysteine residues. The solution structure of rTst2 was determined using multidimensional NMR spectroscopy and revealed that rTst2 adopts an inhibitor cystine knot (ICK) structure. rTst2 was screened using various functional assays, including patch-clamp electrophysiological and cytotoxicity assays. rTst2 was inactive against voltage-gated sodium channels (NaV) and the human voltage-gated proton (hHv1) channel. rTst2 also did not possess cytotoxic activity when assessed against Drosophila melanogaster flies. However, the recombinant peptide at 100 nM showed >50% inhibition of the transient receptor potential subfamily V member 1 (TRPV1) and slight (∼10%) inhibition of transient receptor potential subfamily A member 1 (TRPA1). Tst2 is thus a novel ICK inhibitor of the TRPV1 channel.

Advances in genome sequencing reveal changes in gene content that contribute to arthropod macroevolution

Current sequencing technology allows for the relatively affordable generation of highly contiguous genomes. Technological advances have made it possible for researchers to investigate the consequences of diverse sorts of genomic variants, such as gene gain and loss. With the extraordinary number of high-quality genomes now available, we take stock of how these genomic variants impact phenotypic evolution. We take care to point out that the identification of genomic variants of interest is only the first step in understanding their impact. Painstaking lab or fieldwork is still required to establish causal relationships between genomic variants and phenotypic evolution. We focus mostly on arthropod research, as this phylum has an impressive degree of phenotypic diversity and is also the subject of much evolutionary genetics research. This article is intended to both highlight recent advances in the field and also to be a primer for learning about evolutionary genetics and genomics.

Spatial Distribution and Biochemical Characterization of Serine Peptidase Inhibitors in the Venom of the Brazilian Sea Anemone Anthopleura cascaia Using Mass Spectrometry Imaging

Sea anemones are known to produce a diverse array of toxins with different cysteine-rich peptide scaffolds in their venoms. The serine peptidase inhibitors, specifically Kunitz inhibitors, are an important toxin family that is believed to function as defensive peptides, as well as prevent proteolysis of other secreted anemone toxins. In this study, we isolated three serine peptidase inhibitors named Anthopleura cascaia peptide inhibitors I, II, and III (ACPI-I, ACPI-II, and ACPI-III) from the venom of the endemic Brazilian sea anemone A. cascaia. The venom was fractionated using RP-HPLC, and the inhibitory activity of these fractions against trypsin was determined and found to range from 59% to 93%. The spatial distribution of the anemone peptides throughout A. cascaia was observed using mass spectrometry imaging. The inhibitory peptides were found to be present in the tentacles, pedal disc, and mesenterial filaments. We suggest that the three inhibitors observed during this study belong to the venom Kunitz toxin family on the basis of their similarity to PI-actitoxin-aeq3a-like and the identification of amino acid residues that correspond to a serine peptidase binding site. Our findings expand our understanding of the diversity of toxins present in sea anemone venom and shed light on their potential role in protecting other venom components from proteolysis.

Spatial metabolomics for symbiotic marine invertebrates

Microbial symbionts frequently localize within specific body structures or cell types of their multicellular hosts. This spatiotemporal niche is critical to host health, nutrient exchange, and fitness. Measuring host-microbe metabolite exchange has conventionally relied on tissue homogenates, eliminating dimensionality and dampening analytical sensitivity. We have developed a mass spectrometry imaging workflow for a soft- and hard-bodied cnidarian animal capable of revealing the host and symbiont metabolome in situ, without the need for a priori isotopic labelling or skeleton decalcification. The mass spectrometry imaging method provides critical functional insights that cannot be gleaned from bulk tissue analyses or other presently available spatial methods. We show that cnidarian hosts may regulate microalgal symbiont acquisition and rejection through specific ceramides distributed throughout the tissue lining the gastrovascular cavity. The distribution pattern of betaine lipids showed that once resident, symbionts primarily reside in light-exposed tentacles to generate photosynthate. Spatial patterns of these metabolites also revealed that symbiont identity can drive host metabolism.

Genomic, functional and structural analyses elucidate evolutionary innovation within the sea anemone 8 toxin family

BACKGROUND

The ShK toxin from Stichodactyla helianthus has established the therapeutic potential of sea anemone venom peptides, but many lineage-specific toxin families in Actiniarians remain uncharacterised. One such peptide family, sea anemone 8 (SA8), is present in all five sea anemone superfamilies. We explored the genomic arrangement and evolution of the SA8 gene family in Actinia tenebrosa and Telmatactis stephensoni, characterised the expression patterns of SA8 sequences, and examined the structure and function of SA8 from the venom of T. stephensoni.

RESULTS

We identified ten SA8-family genes in two clusters and six SA8-family genes in five clusters for T. stephensoni and A. tenebrosa, respectively. Nine SA8 T. stephensoni genes were found in a single cluster, and an SA8 peptide encoded by an inverted SA8 gene from this cluster was recruited to venom. We show that SA8 genes in both species are expressed in a tissue-specific manner and the inverted SA8 gene has a unique tissue distribution. While the functional activity of the SA8 putative toxin encoded by the inverted gene was inconclusive, its tissue localisation is similar to toxins used for predator deterrence. We demonstrate that, although mature SA8 putative toxins have similar cysteine spacing to ShK, SA8 peptides are distinct from ShK peptides based on structure and disulfide connectivity.

CONCLUSIONS

Our results provide the first demonstration that SA8 is a unique gene family in Actiniarians, evolving through a variety of structural changes including tandem and proximal gene duplication and an inversion event that together allowed SA8 to be recruited into the venom of T. stephensoni.

Acontia, a Specialised Defensive Structure, Has Low Venom Complexity in Calliactis polypus

Phylum Cnidaria represents a unique group among venomous taxa, with its delivery system organised as individual organelles, known as nematocysts, heterogeneously distributed across morphological structures rather than packaged as a specialised organ. Acontia are packed with large nematocysts that are expelled from sea anemones during aggressive encounters with predatory species and are found in a limited number of species in the superfamily Metridioidea. Little is known about this specialised structure other than the commonly accepted hypothesis of its role in defence and a rudimentary understanding of its toxin content and activity. This study utilised previously published transcriptomic data and new proteomic analyses to expand this knowledge by identifying the venom profile of acontia in Calliactis polypus. Using mass spectrometry, we found limited toxin diversity in the proteome of acontia, with an abundance of a sodium channel toxin type I, and a novel toxin with two ShK-like domains. Additionally, genomic evidence suggests that the proposed novel toxin is ubiquitous across sea anemone lineages. Overall, the venom profile of acontia in Calliactis polypus and the novel toxin identified here provide the basis for future research to define the function of acontial toxins in sea anemones.

Micro and macroevolution of sea anemone venom phenotype

Venom is a complex trait with substantial inter- and intraspecific variability resulting from strong selective pressures acting on the expression of many toxic proteins. However, understanding the processes underlying toxin expression dynamics that determine the venom phenotype remains unresolved. By interspecific comparisons we reveal that toxin expression in sea anemones evolves rapidly and that in each species different toxin family dictates the venom phenotype by massive gene duplication events. In-depth analysis of the sea anemone, Nematostella vectensis, revealed striking variation of the dominant toxin (Nv1) diploid copy number across populations (1-24 copies) resulting from independent expansion/contraction events, which generate distinct haplotypes. Nv1 copy number correlates with expression at both the transcript and protein levels with one population having a near-complete loss of Nv1 production. Finally, we establish the dominant toxin hypothesis which incorporates observations in other venomous lineages that animals have convergently evolved a similar strategy in shaping their venom.

Venom system variation and the division of labor in the colonial hydrozoan Hydractinia symbiolongicarpus

Cnidarians (jellyfish, hydroids, sea anemones, and corals) possess a unique method for venom production, maintenance, and deployment through a decentralized system composed of different types of venom-filled stinging structures called nematocysts. In many species, nematocyst types are distributed heterogeneously across functionally distinct tissues. This has led to a prediction that different nematocyst types contain specific venom components. The colonial hydrozoan, Hydractinia symbiolongicarpus, is an ideal system to study the functional distribution of nematocyst types and their venoms, given that they display a division of labor through functionally distinct polyps within the colony. Here, we characterized the composition and distribution of nematocysts (cnidome) in the different polyp types and show that the feeding polyp (gastrozooid) has a distinct cnidome compared to the reproductive (gonozooid) and predatory polyp (dactylozooid). We generated a nematocyst-specific reporter line to track nematocyst development (nematogenesis) in H. symbiolongicarpus, and were able to confirm that nematogenesis primarily occurs in the mid-region of the gastrozooid and throughout stolons (tubes of epithelia that connect the polyps in the colony). This reporter line enabled us to isolate a nematocyst-specific lineage of cells for de novo transcriptome assembly, annotate venom-like genes (VLGs) and determine differential expression (DE) across polyp types. We show that a majority of VLGs are upregulated in gastrozooids, consistent with it being the primary site of active nematogenesis. However, despite gastrozooids producing more nematocysts, we found a number of VLGs significantly upregulated in dactylozooids, suggesting that these VLGs may be important for prey-capture. Our transgenic Hydractinia reporter line provides an opportunity to explore the complex interplay between venom composition, nematocyst diversity, and ecological partitioning in a colonial hydrozoan that displays a division of labor.

•

Functionally specific polyp types in Hydractinia symbiolongicarpus have distinct cnidomes.

•

We present a nematocyst-targeted transgenic line for H. symbiolongicarpus, showcasing active areas of nematogenesis.

•

105 venom-like genes (VLGs) were annotated from an assembled nematocyst-enriched transcriptome.

•

Several VLGs were significantly upregulated in feeding polyps, consistent with being a site of active nematogenesis.

•

Differential expression analysis suggests that different polyp types express distinct combinations of VLGs.

Identification, diversity and domain structure analysis of mucin and mucin-like genes in sea anemone Actinia tenebrosa

Background

Mucins are part of the glycoprotein family and the main proteinaceous component of mucus. The sea anemone species, Actinia tenebrosa (Phylum Cnidaria) produce large amounts of mucus, which have not been studied in detail. Furthermore, there has only been limited investigation of mucin genes in phylum Cnidaria. Therefore, the aim of current study was to identify and analyse the repertoire mucin genes present in A. tenebrosa and range of other sea anemone species to document their diversity in this group.

Methods

To achieve this aim, we undertook transcriptome sequencing, assembly, and annotation to identify mucin genes in A. tenebrosa.

Results

The results from this study demonstrated a diverse repertoire of mucin proteins, including mucin1-like, mucin4-like, and a range of mucin-like genes in the range of sea anemone species examined. The domain structure of the identified mucin genes was found to be consistent with the conserved domains found in the homologous proteins of vertebrate species. The discovery of a diverse range of mucin genes in sea anemone species provided a basic reference for future mucin studies in cnidarians and could lead to research into their application in the pharmacological, clinical, and cosmetic industries.

Transcriptomic Insights into the Diversity and Evolution of Myxozoa (Cnidaria, Endocnidozoa) Toxin-like Proteins

Myxozoa is a speciose group of endoparasitic cnidarians that can cause severe ecological and economic effects. Their cnidarian affinity is affirmed by genetic relatedness and the presence of nematocysts, historically called “polar capsules”. Previous studies have revealed the presence of toxin-like proteins in myxozoans; however, the diversity and evolution of venom in Myxozoa are not fully understood. Here, we performed a comparative analysis using the newly sequenced transcriptomes of five Myxobolidae species as well as some public datasets. Toxin mining revealed that myxozoans have lost most of their toxin families, while most species retained Kunitz, M12B, and CRISP, which may play a role in endoparasitism. The venom composition of Endocnidozoa (Myxozoa + Polypodium) differs from that of free-living cnidarians and may be influenced by ecological and environmental factors. Phylogenetic analyses showed that toxin families of myxozoans and free-living cnidarians were clustered into different clades. Selection analyses showed that purifying selection was the dominant evolutionary pressure in toxins, while they were still influenced by episodic adaptive selection. This suggests that the potency or specificity of a particular toxin or species might increase. Overall, our findings provide a more comprehensive framework for understanding the diversity and evolution of Myxozoa venoms.

Transcriptomic analysis reveals potential candidate pathways and genes involved in toxin biosynthesis in true toads

Synthesized chemical defenses have broadly evolved across countless taxa and are important in shaping evolutionary and ecological interactions within ecosystems. However, the underlying genomic mechanisms by which these organisms synthesize and utilize their toxins are relatively unknown. Herein, we use comparative transcriptomics to uncover potential toxin synthesizing genes and pathways, as well as interspecific patterns of toxin synthesizing genes across ten species of North American true toads (Bufonidae). Upon assembly and annotation of the ten transcriptomes, we explored patterns of relative gene expression and possible protein-protein interactions across the species to determine what genes and/or pathways may be responsible for toxin synthesis. We also tested our transcriptome dataset for signatures of positive selection to reveal how selection may be acting upon potential toxin producing genes. We assembled high quality transcriptomes of the bufonid parotoid gland, a tissue not often investigated in other bufonid related RNAseq studies. We found several genes involved in metabolic and biosynthetic pathways (e.g. steroid biosynthesis, terpenoid backbone biosynthesis, isoquinoline biosynthesis, glucosinolate biosynthesis) that were functionally enriched and/or relatively expressed across the ten focal species that may be involved in the synthesis of alkaloid and steroid toxins, as well as other small metabolic compounds that cause distastefulness in bufonids. We hope that our study lays a foundation for future studies to explore the genomic underpinnings and specific pathways of toxin synthesis in toads, as well as at the macroevolutionary scale across numerous taxa that produce their own defensive toxins.

Venoms for all occasions: The functional toxin profiles of different anatomical regions in sea anemones are related to their ecological function

The phylum Cnidaria is the oldest extant venomous group and is defined by the presence of nematocysts, specialized organelles responsible for venom production and delivery. Although toxin peptides and the cells housing nematocysts are distributed across the entire animal, nematocyte and venom profiles have been shown to differ across morphological structures in actiniarians. In this study, we explore the relationship between patterns of toxin expression and the ecological roles of discrete anatomical structures in Telmatactis stephensoni. Specifically, using a combination of proteomic and transcriptomic approaches, we examined whether there is a direct correlation between the functional similarity of regions and the similarity of their associated toxin expression profiles. We report that the regionalization of toxin production is consistent with the partitioning of the ecological roles of venom across envenomating structures, and that three major functional regions are present in T. stephensoni: tentacles, epidermis and gastrodermis. Additionally, we find that most structures that serve similar functions not only have comparable putative toxin profiles but also similar nematocyst types. There was no overlap in the putative toxins identified using proteomics and transcriptomics, but the expression patterns of specific milked venom peptides were conserved across RNA-sequencing and mass spectrometry imaging data sets. Furthermore, based on our data, it appears that acontia of T. stephensoni may be transcriptionally inactive and only mature nematocysts are present in the distal portions of the threads. Overall, we find that the venom profile of different anatomical regions in sea anemones varies according to its ecological functions.

The Tentacular Spectacular: Evolution of Regeneration in Sea Anemones

Sea anemones vary immensely in life history strategies, environmental niches and their ability to regenerate. While the sea anemone Nematostella vectensis is the starlet of many key regeneration studies, recent work is emerging on the diverse regeneration strategies employed by other sea anemones. This manuscript will explore current molecular mechanisms of regeneration employed by non-model sea anemones Exaiptasia diaphana (an emerging model species for coral symbiosis studies) and Calliactis polypus (a less well-studied species) and examine how these species compare to the model sea anemone N. vectensis. We summarize the field of regeneration within sea anemones, within the greater context of phylum Cnidaria and in other invertebrate models of regeneration. We also address the current knowledge on two key systems that may be implemented in regeneration: the innate immune system and developmental pathways, including future aspects of work and current limitations.

Tentacle Morphological Variation Coincides with Differential Expression of Toxins in Sea Anemones

Phylum Cnidaria is an ancient venomous group defined by the presence of cnidae, specialised organelles that serve as venom delivery systems. The distribution of cnidae across the body plan is linked to regionalisation of venom production, with tissue-specific venom composition observed in multiple actiniarian species. In this study, we assess whether morphological variants of tentacles are associated with distinct toxin expression profiles and investigate the functional significance of specialised tentacular structures. Using five sea anemone species, we analysed differential expression of toxin-like transcripts and found that expression levels differ significantly across tentacular structures when substantial morphological variation is present. Therefore, the differential expression of toxin genes is associated with morphological variation of tentacular structures in a tissue-specific manner. Furthermore, the unique toxin profile of spherical tentacular structures in families Aliciidae and Thalassianthidae indicate that vesicles and nematospheres may function to protect branched structures that host a large number of photosynthetic symbionts. Thus, hosting zooxanthellae may account for the tentacle-specific toxin expression profiles observed in the current study. Overall, specialised tentacular structures serve unique ecological roles and, in order to fulfil their functions, they possess distinct venom cocktails.

Phylogenetic and Selection Analysis of an Expanded Family of Putatively Pore-Forming Jellyfish Toxins (Cnidaria: Medusozoa)

Many jellyfish species are known to cause a painful sting, but box jellyfish (class Cubozoa) are a well-known danger to humans due to exceptionally potent venoms. Cubozoan toxicity has been attributed to the presence and abundance of cnidarian-specific pore-forming toxins called jellyfish toxins (JFTs), which are highly hemolytic and cardiotoxic. However, JFTs have also been found in other cnidarians outside of Cubozoa, and no comprehensive analysis of their phylogenetic distribution has been conducted to date. Here, we present a thorough annotation of JFTs from 147 cnidarian transcriptomes and document 111 novel putative JFTs from over 20 species within Medusozoa. Phylogenetic analyses show that JFTs form two distinct clades, which we call JFT-1 and JFT-2. JFT-1 includes all known potent cubozoan toxins, as well as hydrozoan and scyphozoan representatives, some of which were derived from medically relevant species. JFT-2 contains primarily uncharacterized JFTs. Although our analyses detected broad purifying selection across JFTs, we found that a subset of cubozoan JFT-1 sequences are influenced by gene-wide episodic positive selection compared with homologous toxins from other taxonomic groups. This suggests that duplication followed by neofunctionalization or subfunctionalization as a potential mechanism for the highly potent venom in cubozoans. Additionally, published RNA-seq data from several medusozoan species indicate that JFTs are differentially expressed, spatially and temporally, between functionally distinct tissues. Overall, our findings suggest a complex evolutionary history of JFTs involving duplication and selection that may have led to functional diversification, including variability in toxin potency and specificity.

An ancient, conserved gene regulatory network led to the rise of oral venom systems

Although oral venom systems are ecologically important characters, how they originated is still unclear. In this study, we show that oral venom systems likely originated from a gene regulatory network conserved across amniotes. This network, which we term the “metavenom network,” comprises over 3,000 housekeeping genes coexpressed with venom and play a role in protein folding and modification. Comparative transcriptomics revealed that the network is conserved between venom glands of snakes and salivary glands of mammals. This suggests that while these tissues have evolved different functions, they share a common regulatory core, that persisted since their common ancestor. We propose several evolutionary mechanisms that can utilize this common regulatory core to give rise to venomous animals from their nonvenomous ancestors.

Oral venom systems evolved multiple times in numerous vertebrates enabling the exploitation of unique predatory niches. Yet how and when they evolved remains poorly understood. Up to now, most research on venom evolution has focused strictly on the toxins. However, using toxins present in modern day animals to trace the origin of the venom system is difficult, since they tend to evolve rapidly, show complex patterns of expression, and were incorporated into the venom arsenal relatively recently. Here we focus on gene regulatory networks associated with the production of toxins in snakes, rather than the toxins themselves. We found that overall venom gland gene expression was surprisingly well conserved when compared to salivary glands of other amniotes. We characterized the “metavenom network,” a network of ∼3,000 nonsecreted housekeeping genes that are strongly coexpressed with the toxins, and are primarily involved in protein folding and modification. Conserved across amniotes, this network was coopted for venom evolution by exaptation of existing members and the recruitment of new toxin genes. For instance, starting from this common molecular foundation, Heloderma lizards, shrews, and solenodon, evolved venoms in parallel by overexpression of kallikreins, which were common in ancestral saliva and induce vasodilation when injected, causing circulatory shock. Derived venoms, such as those of snakes, incorporated novel toxins, though still rely on hypotension for prey immobilization. These similarities suggest repeated cooption of shared molecular machinery for the evolution of oral venom in mammals and reptiles, blurring the line between truly venomous animals and their ancestors.

Insights into how development and life-history dynamics shape the evolution of venom

Venomous animals are a striking example of the convergent evolution of a complex trait. These animals have independently evolved an apparatus that synthesizes, stores, and secretes a mixture of toxic compounds to the target animal through the infliction of a wound. Among these distantly related animals, some can modulate and compartmentalize functionally distinct venoms related to predation and defense. A process to separate distinct venoms can occur within and across complex life cycles as well as more streamlined ontogenies, depending on their life-history requirements. Moreover, the morphological and cellular complexity of the venom apparatus likely facilitates the functional diversity of venom deployed within a given life stage. Intersexual variation of venoms has also evolved further contributing to the massive diversity of toxic compounds characterized in these animals. These changes in the biochemical phenotype of venom can directly affect the fitness of these animals, having important implications in their diet, behavior, and mating biology. In this review, we explore the current literature that is unraveling the temporal dynamics of the venom system that are required by these animals to meet their ecological functions. These recent findings have important consequences in understanding the evolution and development of a convergent complex trait and its organismal and ecological implications.

A disulfide-stabilised helical hairpin fold in acrorhagin I: An emerging structural motif in peptide toxins

Acrorhagin I (U-AITX-Aeq5a) is a disulfide-rich peptide identified in the aggressive organs (acrorhagi) of the sea anemone Actinia equina. Previous studies (Toxicon 2005, 46:768-74) found that the peptide is toxic in crabs, although the structural and functional properties of acrorhagin I have not been reported. In this work, an Escherichia coli (BL21 strain) expression system was established for the preparation of ¹³C,¹⁵N-labelled acrorhagin I, and the solution structure was determined using NMR spectroscopy. Structurally, acrorhagin I is similar to B-IV toxin from the marine worm Cerebratulus lacteus (PDB id 1VIB), with a well-defined helical hairpin structure stabilised by four intramolecular disulfide bonds. The recombinant peptide was tested in patch-clamp electrophysiology assays against voltage-gated potassium and sodium channels, and in bacterial and fungal growth inhibitory assays and haemolytic assays. Acrorhagin I was not active against any of the ion channels tested and showed no activity in functional assays, indicating that this peptide may possess a different biological function. Metal ion interaction studies using NMR spectroscopy showed that acrorhagin I bound zinc and nickel, suggesting that its function might be modulated by metal ions or that it may be involved in regulating metal ion levels and their transport. The similarity between the structure of acrorhagin I and that of B-IV toxin from a marine worm suggests that this fold may prove to be a recurring motif in disulfide-rich peptides from marine organisms.

How the Toxin got its Toxicity

Venom systems are functional and ecological traits, typically used by one organism to subdue or deter another. A predominant subset of their constituent molecules—“toxins”—share this ecological function and are therefore molecules that mediate interactions between organisms. Such molecules have been referred to as “exochemicals.” There has been debate within the field of toxinology concerning the evolutionary pathways leading to the “recruitment” of a gene product for a toxic role within venom. We review these discussions and the evidence interpreted in support of alternate pathways, along with many of the most popular models describing the origin of novel molecular functions in general. We note that such functions may arise with or without gene duplication occurring and are often the consequence of a gene product encountering a novel “environment,” i.e., a range of novel partners for molecular interaction. After stressing the distinction between “activity” and “function,” we describe in detail the results of a recent study which reconstructed the evolutionary history of a multigene family that has been recruited as a toxin and argue that these results indicate that a pluralistic approach to understanding the origin of novel functions is advantageous. This leads us to recommend that an expansive approach be taken to the definition of “neofunctionalization”—simply the origins of a novel molecular function by any process—and “recruitment”—the “weaponization” of a molecule via the acquisition of a toxic function in venom, by any process. Recruitment does not occur at the molecular level or even at the level of gene expression, but only when a confluence of factors results in the ecological deployment of a physiologically active molecule as a toxin. Subsequent to recruitment, the evolutionary regime of a gene family may shift into a more dynamic form of “birth-and-death.” Thus, recruitment leads to a form of “downwards causation,” in which a change at the ecological level at which whole organisms interact leads to a change in patterns of evolution at the genomic level.

Some like it hot: population-specific adaptations in venom production to abiotic stressors in a widely distributed cnidarian

BACKGROUND

In cnidarians, antagonistic interactions with predators and prey are mediated by their venom, whose synthesis may be metabolically expensive. The potentially high cost of venom production has been hypothesized to drive population-specific variation in venom expression due to differences in abiotic conditions. However, the effects of environmental factors on venom production have been rarely demonstrated in animals. Here, we explore the impact of specific abiotic stresses on venom production of distinct populations of the sea anemone Nematostella vectensis (Actiniaria, Cnidaria) inhabiting estuaries over a broad geographic range where environmental conditions such as temperatures and salinity vary widely.

RESULTS

We challenged Nematostella polyps with heat, salinity, UV light stressors, and a combination of all three factors to determine how abiotic stressors impact toxin expression for individuals collected across this species' range. Transcriptomics and proteomics revealed that the highly abundant toxin Nv1 was the most downregulated gene under heat stress conditions in multiple populations. Physiological measurements demonstrated that venom is metabolically costly to produce. Strikingly, under a range of abiotic stressors, individuals from different geographic locations along this latitudinal cline modulate differently their venom production levels.

CONCLUSIONS

We demonstrate that abiotic stress results in venom regulation in Nematostella. Together with anecdotal observations from other cnidarian species, our results suggest this might be a universal phenomenon in Cnidaria. The decrease in venom production under stress conditions across species coupled with the evidence for its high metabolic cost in Nematostella suggests downregulation of venom production under certain conditions may be highly advantageous and adaptive. Furthermore, our results point towards local adaptation of this mechanism in Nematostella populations along a latitudinal cline, possibly resulting from distinct genetics and significant environmental differences between their habitats.

Transcriptomic Analysis of Four Cerianthid (Cnidaria, Ceriantharia) Venoms

Tube anemones, or cerianthids, are a phylogenetically informative group of cnidarians with complex life histories, including a pelagic larval stage and tube-dwelling adult stage, both known to utilize venom in stinging-cell rich tentacles. Cnidarians are an entirely venomous group that utilize their proteinaceous-dominated toxins to capture prey and defend against predators, in addition to several other ecological functions, including intraspecific interactions. At present there are no studies describing the venom for any species within cerianthids. Given their unique development, ecology, and distinct phylogenetic-placement within Cnidaria, our objective is to evaluate the venom-like gene diversity of four species of cerianthids from newly collected transcriptomic data. We identified 525 venom-like genes between all four species. The venom-gene profile for each species was dominated by enzymatic protein and peptide families, which is consistent with previous findings in other cnidarian venoms. However, we found few toxins that are typical of sea anemones and corals, and furthermore, three of the four species express toxin-like genes closely related to potent pore-forming toxins in box jellyfish. Our study is the first to provide a survey of the putative venom composition of cerianthids and contributes to our general understanding of the diversity of cnidarian toxins.

Venom Systems as Models for Studying the Origin and Regulation of Evolutionary Novelties

A central goal in biology is to determine the ways in which evolution repeats itself. One of the most remarkable examples in nature of convergent evolutionary novelty is animal venom. Across diverse animal phyla, various specialized organs and anatomical structures have evolved from disparate developmental tissues to perform the same function, that is, produce and deliver a cocktail of potent molecules to subdue prey or predators. Venomous organisms therefore offer unique opportunities to investigate the evolutionary processes of convergence of key adaptive traits, and the molecular mechanisms underlying the emergence of novel genes, cells, and tissues. Indeed, some venomous species have already proven to be highly amenable as models for developmental studies, and recent work with venom gland organoids provides manipulatable systems for directly testing important evolutionary questions. Here, we provide a synthesis of the current knowledge that could serve as a starting point for the establishment of venom systems as new models for evolutionary and molecular biology. In particular, we highlight the potential of various venomous species for the study of cell differentiation and cell identity, and the regulatory dynamics of rapidly evolving, highly expressed, tissue-specific, gene paralogs. We hope that this review will encourage researchers to look beyond traditional study organisms and consider venom systems as useful tools to explore evolutionary novelties.

Characterising Functional Venom Profiles of Anthozoans and Medusozoans within Their Ecological Context

This review examines the current state of knowledge regarding toxins from anthozoans (sea anemones, coral, zoanthids, corallimorphs, sea pens and tube anemones). We provide an overview of venom from phylum Cnidaria and review the diversity of venom composition between the two major clades (Medusozoa and Anthozoa). We highlight that the functional and ecological context of venom has implications for the temporal and spatial expression of protein and peptide toxins within class Anthozoa. Understanding the nuances in the regulation of venom arsenals has been made possible by recent advances in analytical technologies that allow characterisation of the spatial distributions of toxins. Furthermore, anthozoans are unique in that ecological roles can be assigned using tissue expression data, thereby circumventing some of the challenges related to pharmacological screening.

Tentacle Transcriptomes of the Speckled Anemone (Actiniaria: Actiniidae: Oulactis sp.): Venom-Related Components and Their Domain Structure

Cnidarians are one of the oldest known animal lineages (ca. 700 million years), with a unique envenomation apparatus to deliver a potent mixture of peptides and proteins. Some peptide toxins from cnidarian venom have proven therapeutic potential. Here, we use a transcriptomic/proteomic strategy to identify sequences with similarity to known venom protein families in the tentacles of the endemic Australian 'speckled anemone' (Oulactis sp.). Illumina RNASeq data were assembled de novo. Annotated sequences in the library were verified by cross-referencing individuals' transcriptomes or protein expression evidence from LC-MS/MS data. Sequences include pore-forming toxins, phospholipases, peptidases, neurotoxins (sodium and potassium channel modulators), cysteine-rich secretory proteins and defensins (antimicrobial peptides). Fewer than 4% of the sequences in the library occurred across the three individuals examined, demonstrating high sequence variability of an individual's arsenal. We searched for actinoporins in Oulactis sp. to assess sequence similarity to the only described toxins (OR-A and -G) for this genus and examined the domain architecture of venom-related peptides and proteins. The novel putative actinoporin of Oulactis sp. has a greater similarity to other species in the Actiniidae family than to O. orientalis. Venom-related sequences have an architecture that occurs in single, repeat or multi-domain combinations of venom-related (e.g. ShK-like) and non-venom (e.g. whey acid protein) domains. This study has produced the first transcriptomes for an endemic Australian sea anemone species and the genus Oulactis, while identifying nearly 400 novel venom-related peptides and proteins for future structural and functional analyses and venom evolution studies.

The Rapid Regenerative Response of a Model Sea Anemone Species Exaiptasia pallida Is Characterised by Tissue Plasticity and Highly Coordinated Cell Communication

Regeneration of a limb or tissue can be achieved through multiple different pathways and mechanisms. The sea anemone Exaiptasia pallida has been observed to have excellent regenerative proficiency, but this has not yet been described transcriptionally. In this study, we examined the genetic expression changes during a regenerative timecourse and reported key genes involved in regeneration and wound healing. We found that the major response was an early (within the first 8 h) upregulation of genes involved in cellular movement and cell communication, which likely contribute to a high level of tissue plasticity resulting in the rapid regeneration response observed in this species. We find the immune system was only transcriptionally active in the first 8 h post-amputation and conclude, in accordance with previous literature, that the immune system and regeneration have an inverse relationship. Fifty-nine genes (3.8% of total) differentially expressed during regeneration were identified as having no orthologues in other species, indicating that regeneration in E. pallida may rely on the activation of species-specific novel genes. Additionally, taxonomically restricted novel genes, including species-specific novels, and highly conserved genes were identified throughout the regenerative timecourse, showing that both may work in concert to achieve complete regeneration.

The genome of the sea anemone Actinia equina (L.): Meiotic toolkit genes and the question of sexual reproduction

The beadlet anemone Actinia equina (L.) (Cnidaria: Anthozoa: Actiniaria: Actiniidae) is one of the most familiar organisms of the North European intertidal zone. Once considered a single, morphologically variable species across northern Europe, it is now recognised as one member of a variable species complex. Previous studies of distribution, aggression, allozymes and mitochondrial DNA suggest that the diversity in form and colour within A. equina may hide still unrecognised species diversity. To empower further study of A. equina population genetics and systematics, we sequenced (PacBio Sequel) the genome of a single A. equina individual to produce a high-quality genome assembly (contig N₅₀ = 492,607 bp, 1485 contigs, number of protein coding genes = 47,671, 97% BUSCO completeness). There is debate as to whether A. equina reproduces solely asexually, since no reliable, consistent evidence of sexual reproduction has been found. To gain further insight, we examined the genome for evidence of a 'meiotic toolkit' - genes believed to be found consistently in sexually reproducing organisms - and demonstrate that the A. equina genome appears not to have this full complement. Additionally, Smudgeplot analysis, coupled with high haplotype diversity, indicates this genome assembly to be of ambiguous ploidy, suggesting that A. equina may not be diploid. The suggested polyploid nature of this species coupled with the deficiency in meiotic toolkit genes, indicates that further field and laboratory studies of this species is warranted to understand how this species reproduces and what role ploidy may play in speciation within this speciose genus.

A Versatile and Robust Serine Protease Inhibitor Scaffold from Actinia tenebrosa

Serine proteases play pivotal roles in normal physiology and a spectrum of patho-physiological processes. Accordingly, there is considerable interest in the discovery and design of potent serine protease inhibitors for therapeutic applications. This led to concerted efforts to discover versatile and robust molecular scaffolds for inhibitor design. This investigation is a bioprospecting study that aims to isolate and identify protease inhibitors from the cnidarian Actinia tenebrosa. The study isolated two Kunitz-type protease inhibitors with very similar sequences but quite divergent inhibitory potencies when assayed against bovine trypsin, chymostrypsin, and a selection of human sequence-related peptidases. Homology modeling and molecular dynamics simulations of these inhibitors in complex with their targets were carried out and, collectively, these methodologies enabled the definition of a versatile scaffold for inhibitor design. Thermal denaturation studies showed that the inhibitors were remarkably robust. To gain a fine-grained map of the residues responsible for this stability, we conducted in silico alanine scanning and quantified individual residue contributions to the inhibitor's stability. Sequences of these inhibitors were then used to search for Kunitz homologs in an A. tenebrosa transcriptome library, resulting in the discovery of a further 14 related sequences. Consensus analysis of these variants identified a rich molecular diversity of Kunitz domains and expanded the palette of potential residue substitutions for rational inhibitor design using this domain.

The Diversity of Venom: The Importance of Behavior and Venom System Morphology in Understanding Its Ecology and Evolution

Venoms are one of the most convergent of animal traits known, and encompass a much greater taxonomic and functional diversity than is commonly appreciated. This knowledge gap limits the potential of venom as a model trait in evolutionary biology. Here, we summarize the taxonomic and functional diversity of animal venoms and relate this to what is known about venom system morphology, venom modulation, and venom pharmacology, with the aim of drawing attention to the importance of these largely neglected aspects of venom research. We find that animals have evolved venoms at least 101 independent times and that venoms play at least 11 distinct ecological roles in addition to predation, defense, and feeding. Comparisons of different venom systems suggest that morphology strongly influences how venoms achieve these functions, and hence is an important consideration for understanding the molecular evolution of venoms and their toxins. Our findings also highlight the need for more holistic studies of venom systems and the toxins they contain. Greater knowledge of behavior, morphology, and ecologically relevant toxin pharmacology will improve our understanding of the evolution of venoms and their toxins, and likely facilitate exploration of their potential as sources of molecular tools and therapeutic and agrochemical lead compounds.

The draft genome of Actinia tenebrosa reveals insights into toxin evolution

Sea anemones have a wide array of toxic compounds (peptide toxins found in their venom) which have potential uses as therapeutics. To date, the majority of studies characterizing toxins in sea anemones have been restricted to species from the superfamily, Actinioidea. No highly complete draft genomes are currently available for this superfamily, however, highlighting our limited understanding of the genes encoding toxins in this important group. Here we have sequenced, assembled, and annotated a draft genome for Actinia tenebrosa. The genome is estimated to be approximately 255 megabases, with 31,556 protein-coding genes. Quality metrics revealed that this draft genome matches the quality and completeness of other model cnidarian genomes, including Nematostella, Hydra, and Acropora. Phylogenomic analyses revealed strong conservation of the Cnidaria and Hexacorallia core-gene set. However, we found that lineage-specific gene families have undergone significant expansion events compared with shared gene families. Enrichment analysis performed for both gene ontologies, and protein domains revealed that genes encoding toxins contribute to a significant proportion of the lineage-specific genes and gene families. The results make clear that the draft genome of A. tenebrosa will provide insight into the evolution of toxins and lineage-specific genes, and provide an important resource for the discovery of novel biological compounds.

Expansion of Hydra actinoporin-like toxin (HALT) gene family: Expression divergence and functional convergence evolved through gene duplication

Hydra actinoporin-like toxin 1 (HALT-1) was previously shown to cause cytolysis and haemolysis in a number of human cells and has similar functional properties to the actinoporins equinatoxin and sticholysin. In addition to HALT-1, five other HALTs (HALTs 2, 3, 4, 6 and 7) were also isolated from Hydra magnipapillata and expressed as recombinant proteins in this study. We demonstrated that recombinant HALTs have cytolytic activity on HeLa cells but each exhibited a different range of toxicity. All six recombinant HALTs bound to sulfatide, while rHALT-1 and rHALT-3 bound to two additional sphingolipids, lysophosphatidic acid and sphingosine-1-phosphate as indicated by the protein-lipid overlay assay. When either tryptophan¹³³ or tyrosine¹²⁹ of HALT-1 was mutated, the mutant protein lost binding to sulfatide, lysophosphatidic acid and sphingosine-1-phosphate. As further verification of HALTs' binding to sulfatide, we performed ELISA for each HALT. To determine the cell-type specific gene expression of seven HALTs in Hydra, we searched for individual HALT expression in the single-cell RNA-seq data set of Single Cell Portal. The results showed that HALT-1, 4 and 7 were expressed in differentiating stenoteles. HALT-1 and HALT-6 were expressed in the female germline during oogenesis. HALT-2 was strongly expressed in the gland and mucous cells in the endoderm. Information on HALT-3 and HALT-5 could not be found in the single-cell data set. Our findings show that subfunctionalisation of gene expression following duplication enabled HALTs to become specialized in various cell types of the interstitial cell lineage.

Transcriptomic and Proteomic Analysis of the Tentacles and Mucus of Anthopleura dowii Verrill, 1869

Sea anemone venom contains a complex and diverse arsenal of peptides and proteins of pharmacological and biotechnological interest, however, only venom from a few species has been explored from a global perspective to date. In the present study, we identified the polypeptides present in the venom of the sea anemone Anthopleura dowii Verrill, 1869 through a transcriptomic and proteomic analysis of the tentacles and the proteomic profile of the secreted mucus. In our transcriptomic results, we identified 261 polypeptides related to or predicted to be secreted in the venom, including proteases, neurotoxins that could act as either potassium (K⁺) or sodium (Na⁺) channels inhibitors, protease inhibitors, phospholipases A2, and other polypeptides. Our proteomic data allowed the identification of 156 polypeptides—48 exclusively identified in the mucus, 20 in the tentacles, and 88 in both protein samples. Only 23 polypeptides identified by tandem mass spectrometry (MS/MS) were related to the venom and 21 exclusively identified in the mucus, most corresponding to neurotoxins and hydrolases. Our data contribute to the knowledge of evolutionary and venomic analyses of cnidarians, particularly of sea anemones.

Transcriptional characterisation of the Exaiptasia pallida pedal disc

BACKGROUND

Biological adhesion (bioadhesion), enables organisms to attach to surfaces as well as to a range of other targets. Bioadhesion evolved numerous times independently and is ubiquitous throughout the kingdoms of life. To date, investigations have focussed on various taxa of animals, plants and bacteria, but the fundamental processes underlying bioadhesion and the degree of conservation in different biological systems remain poorly understood. This study had two aims: 1) To characterise tissue-specific gene regulation in the pedal disc of the model cnidarian Exaiptasia pallida, and 2) to elucidate putative genes involved in pedal disc adhesion.

RESULTS

Five hundred and forty-seven genes were differentially expressed in the pedal disc compared to the rest of the animal. Four hundred and twenty-seven genes were significantly upregulated and 120 genes were significantly downregulated. Forty-one condensed gene ontology terms and 19 protein superfamily classifications were enriched in the pedal disc. Eight condensed gene ontology terms and 11 protein superfamily classifications were depleted. Enriched superfamilies were consistent with classifications identified previously as important for the bioadhesion of unrelated marine invertebrates. A host of genes involved in regulation of extracellular matrix generation and degradation were identified, as well as others related to development and immunity. Ab initio prediction identified 173 upregulated genes that putatively code for extracellularly secreted proteins.

CONCLUSION

The analytical workflow facilitated identification of genes putatively involved in adhesion, immunity, defence and development of the E. pallida pedal disc. When defence, immunity and development-related genes were identified, those remaining corresponded most closely to formation of the extracellular matrix (ECM), implicating ECM in the adhesion of anemones to surfaces. This study therefore provides a valuable high-throughput resource for the bioadhesion community and lays a foundation for further targeted research to elucidate bioadhesion in the Cnidaria.

The Birth and Death of Toxins with Distinct Functions: A Case Study in the Sea Anemone Nematostella

The cnidarian Nematostella vectensis has become an established lab model, providing unique opportunities for venom evolution research. The Nematostella venom system is multimodal: involving both nematocytes and ectodermal gland cells, which produce a toxin mixture whose composition changes throughout the life cycle. Additionally, their modes of interaction with predators and prey vary between eggs, larvae, and adults, which is likely shaped by the dynamics of the venom system.

Nv1 is a major component of adult venom, with activity against arthropods (through specific inhibition of sodium channel inactivation) and fish. Nv1 is encoded by a cluster of at least 12 nearly identical genes that were proposed to be undergoing concerted evolution. Surprisingly, we found that Nematostella venom includes several Nv1 paralogs escaping a pattern of general concerted evolution, despite belonging to the Nv1-like family. Here, we show two of these new toxins, Nv4 and Nv5, are lethal for zebrafish larvae but harmless to arthropods, unlike Nv1. Furthermore, unlike Nv1, the newly identified toxins are expressed in early life stages. Using transgenesis and immunostaining, we demonstrate that Nv4 and Nv5 are localized to ectodermal gland cells in larvae.

The evolution of Nv4 and Nv5 can be described either as neofunctionalization or as subfunctionalization. Additionally, the Nv1-like family includes several pseudogenes being an example of nonfunctionalization and venom evolution through birth-and-death mechanism. Our findings reveal the evolutionary history for a toxin radiation and point toward the ecological function of the novel toxins constituting a complex cnidarian venom.