Next Article in Journal
NMR Analyses of the Enzymatic Degradation End-Products of Diabolican: The Secreted EPS of Vibrio diabolicus CNCM I-1629
Previous Article in Journal
Effect of Echinochrome A on Submandibular Gland Dysfunction in Ovariectomized Rats
Previous Article in Special Issue
Transcriptomic Analysis of Four Cerianthid (Cnidaria, Ceriantharia) Venoms
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Marine Drugs
Volume 20
Issue 12
10.3390/md20120730
marinedrugs-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Never, Ever Make an Enemy… Out of an Anemone: Transcriptomic Comparison of Clownfish Hosting Sea Anemone Venoms

by
Alonso Delgado
1,*,
Charlotte Benedict
1,
Jason Macrander
2 and
Marymegan Daly
1
1
Department of Evolution, Ecology and Organismal Biology, The Ohio State University, Columbus, OH 43210, USA
2
Department of Biology, Florida Southern College, Lakeland, FL 33815, USA
*
Author to whom correspondence should be addressed.
Mar. Drugs 2022, 20(12), 730; https://doi.org/10.3390/md20120730
Submission received: 4 October 2022 / Revised: 15 November 2022 / Accepted: 21 November 2022 / Published: 23 November 2022
(This article belongs to the Special Issue Anthozoan Toxins: Using New Approaches to Understand Their Composition, Distribution, and Function)
Browse Figure
Versions Notes

Abstract

:
Sea anemones are predatory marine invertebrates and have diverse venom arsenals. Venom is integral to their biology, and is used in competition, defense, and feeding. Three lineages of sea anemones are known to have independently evolved symbiotic relationships with clownfish, however the evolutionary impact of this relationship on the venom composition of the host is still unknown. Here, we investigate the potential of this symbiotic relationship to shape the venom profiles of the sea anemones that host clownfish. We use transcriptomic data to identify differences and similarities in venom profiles of six sea anemone species, representing the three known clades of clownfish-hosting sea anemones. We recovered 1121 transcripts matching verified toxins across all species, and show that hemolytic and hemorrhagic toxins are consistently the most dominant and diverse toxins across all species examined. These results are consistent with the known biology of sea anemones, provide foundational data on venom diversity of these species, and allow for a review of existing hierarchical structures in venomic studies.

1. Introduction

Cnidaria is the oldest venomous metazoan lineage, exhibiting the greatest functional diversity in venom across animals [1]. Cnidarian venom is stored within glandular cells, as well as within the cnidarian-specific cells called nematocytes that produce microscopic venom delivery structures called nematocysts [2,3]. Nematocysts are a shared feature that unite and diagnose members of the phylum, with structural variation in nematocysts delimiting groups within Cnidaria [4,5]. Venom is an integral component to cnidarian biology as it is used for prey capture, defense, intraspecific aggression, and digestion [6,7,8]. Sea anemones are found in all marine habitats and have widespread ecological success, likely because of the diversity of their ecological and physiological strategies [9,10]. The diverse symbiotic, competitive, and predator-prey interactions of sea anemones certainly rely on venoms [11], although precise correspondence between venoms, ecology, and physiology has yet to be demonstrated (but see [12]). Among the most compelling of the interactions in which venoms play a key role are mutualisms: sea anemones partner with photosynthetic and chemosynthetic microbes (e.g., [13,14]), other marine invertebrates (crustaceans, sponges, gastropods: [15,16,17,18,19]) as well as vertebrates, including clownfish [20].
The clownfish-sea anemone symbiosis has been a model system for research into complex mutualistic relationships. In this mutualism, the clownfish live in, mate, and lay eggs within the area covered by the body and tentacles of the host sea anemone. In at least some of the partnerships, the sea anemone host receives protection from predators by the clownfish [21,22,23] and indirectly benefits through the ammonia provided by the clownfish to the photosymbiotic zooxanthellae inside its tissues [24,25]. Much of the research into this symbiosis has focused on the clownfish, which are all members of the subfamily Amphiprioninae, a clade of ~30 species of clownfishes that has rapidly diversified within the damselfish family Pomacentridae [26,27].
In contrast to clownfish, which represent a single radiation [26,27,28,29], at least three lineages of anemones have independently evolved symbiotic interactions with clownfish [30]. These lineages, designated Stichodactylina (Cryptodendrum, Heteractis magnfica, Stichodactyla, Thalassianthus), Heteractina (Heteractis aurora, Hs. crispa, Hs. malu, and Macrodactyla doreensis), and Entacmaea quadricolor, are all part of superfamily Actinioidea but are not one another’s closest relatives [30].
Although the precise role of venom is unknown, it likely plays a role throughout this mutualistic symbiosis [31]. The venom of host anemones is inferred to be the deterrent to clownfish predators, and thus responsible for at least some of the benefits the fish receive [20,32]. Host anemones do not typically envenomate clownfish while associating with(see for [33] review), however, clownfish are susceptible to the venom of their host species [34]. This interaction is mediated through adaptive behaviors [20] and may change across the lifecycle of the fish [35]. Skin coating on the fish has been implicated as a factor explaining the nematocysts firing supression [31,36,37,38], with at least some of the protective attributes acquired through acclimatization by the fish to its host [39] and others inferred to be innate to the fish [40]. Genome-wide scans of clownfish find positive selection in genes associated with polysaccharides that might provide mechanisms for chemical manipulation of the mucus coat and thus offer some protection against host envenomation [41].
The venom variation across lineages of host sea anemones are poorly known. There are clear differences in the toxicity of venoms from different clownfish-hosting anemones, at least with reference to their mode of action on a crustacean model [42]. Variation in toxin effects does not follow a clear pattern [42], which is unsurprising given the lack of evolutionary relatedness across host anemones [30]. Beyond the ways in which venom might shape this complex and charismatic symbiosis, this system is ripe for investigation into venom constituents that have pharmaceutical potential. This is one of the few instances of clearly coevolutionary interactions between a sea anemone and a vertebrate, which creates interesting potential for specializations for the musculo-nervous system of vertebrates. Evidence of this relevance of the historical connection between host sea anemones and their fish associates for molecules with specificity to other vertebrates comes from Stichodactyla helianthus, the source of the medically promising polypeptide Sea Anemone Toxin ShK [43], and a member of this group.
We investigate the venom in the sea anemones that have co-evolved with clownfish through a comparison of the transcriptomes from clownfish-hosting anemone representatives of each of the three major clades of anemone hosts: Entacmaea; Stichodactylina, and Heteractina. We use transcriptomes because this data source is most widely available for sea anemones, allowing more nuanced comparison with previous studies and other species. We use these accounts of expressed diversity of putative venoms and their function to identify differences and similarities of the venom repertoire of clownfish hosts.

2. Results

2.1. Transcriptome Assemblies

We assembled 20 transcriptomes from six species spanning the three symbiotic lineages using the using the de novo assembly program Trinity v2.2 [44]; Stichodactylina: Cryptodendrum adhaesivum (2), Stichodactyla haddoni (1); Heterodactyla hemprichii (4); Heteractina: Heteractis crispa (3), Macrodactyla doreensis (4); Entacmaea: Entacmaea quadricolor (6). Specific assembly details for concatenated transcriptomes are reported in Table 1, with specifics outlined for each transcriptome in Supplementary Table S1. The number of transcripts across assemblies ranged from 196,433 (S. haddoni) to 1,198,081 (E. quadricolor), with a mean value of ~700,000. The inferred number of genes ranged from 150,054 (S. haddoni) to 581,957 (E. quadricolor). All N50 scores for each species were >1000 bp, except for E. quadricolor (N50 = 545 bp). BUSCO (Benchmarking Universal Single Copy Orthologs) scores for species transcriptomes ranged from 83.6–98.3%.
Transcripts with the presence of a signaling region and with high sequence similarity/bit scores and low E-value to known toxins from ToxProt were identified as candidates and extracted from each transcriptome assembly. Each of these were identified based on PFAM annotation, Tox Prot ID, and NCBI’s non redundant (NR) database BLAST results. These putative toxins were categorized into toxin families and then categorized into functional catagories: auxiliary toxins, allergens and innate immunity toxins, hemostatic and hemorrhagic activity toxins, mixed function enzymes, neurotoxins, pore forming toxins (cytolysins), protease inhibitors, and actiniarian toxins of unknown function. Proportion of functional venom catagories in the total recovered venom arsenal of each of our focal taxa was reported (Figure 1).

2.2. Putative Toxins Inferred from Transcriptomes across Clades

2.2.1. Stichodactylina

We identified 483 putative toxins belonging to 47 toxin gene families (18 actiniarian-specific) from the members of Stichodactylina (C. adhaesivum, Ha. hemprichii, S. haddoni). The greatest diversity and number of putative toxins were recovered from C. adhaesivum, which had 256 putative venom transcripts from 34 toxin gene families in its transcriptome. The transcriptome of C. adhaesivum shows high representation of hemostatic and hemorrhagic toxins (128), auxiliary toxins (40), and neurotoxins (43) which together represent 81% of all putative toxin transcripts. In C. adhaesivum, we found 88 transcripts belonging to 13 actiniarian-specific toxin families, representing 37% of all transcripts. Among the anemone-specific toxin families, C. adhaesivum had a high representation of auxiliary toxin transcripts (39) and neurotoxins (33), which together comprised of 81% of the actiniarian-specific transcripts. Cryptodendron adhaesivum was the sole member of Stichodactylina in which we recover transcripts assigned to the actiniarian-specific toxins sea anemone 8 toxin and sea anemone structural class 9a toxins. We did not recover any mixed enzyme proteins (e.g., PLA2) from the transcriptome of C. adhaesivum.
The 115 putative toxin transcripts we recovered from Ha. hemprichii represent 32 toxin gene families. The transcriptome of Ha. hemprichii shows high representation of hemostatic and hemorrhagic toxins (74) and neurotoxins (18) which together represent 80% of all putative toxin transcripts. Although it contains close to the same number of actiniarian-specific gene families as the other stichodactylines, the transcriptome of Ha. hemprichii had the fewest number of actiniarian-specific toxin transcripts (27/115, 23%). The most abundant of these was the neurotoxin sea anemone type 3 potassium channel toxin (BDS). Despite the poor representation of putative toxin genes in the transcriptome assembly of Ha. hemprichii, this is the sole member of Stichodactylina for which we recovered a transcript for a cnidarian Small Cysteine-Rich Protein (SCRiP) and ShK-like-1 neurotoxin.
A total of 112 putative toxins belonging to 35 toxin gene families were identified from the transcriptome of S. haddoni. The transcriptome of S. haddoni shows high representation of hemostatic and hemorrhagic toxins (74) and toxins with neurotoxin activity (30), which together represent 82% of all putative toxin transcripts. From its transcriptome we identified 46 transcripts belonging to 12 actiniarian-specific toxin families, representing 38% of all toxin transcripts. As was the case for C. adhaesivum, S. haddoni showed strong representation in both the neurotoxin activity transcripts (25) and auxiliary activity transcripts (13), representing 82% of all putative actiniarian toxin transcripts. We did not recover a membrane active protein, actinoporin, or MAC/perforin domain protein from the transcriptome of S. haddoni.

2.2.2. Heteractina

In our representative members of Heteractina, Hs. crispa and M. doreensis, we recovered 310 transcripts assigned to 37 toxin gene families, with 13 of these belonging to actiniarian-specific gene families. In the transcriptome of Hs. crispa, we recovered 175 putative toxin transcripts belonging to 35 gene families. The transcriptome of Hs. crispa was heavily represented by hemostatic and hemorrhagic toxins (90/175, 51%) which made up half of the total number of transcripts recovered. Transcripts identified as venom prothrombin activator (32), ryncolin (13), and veficolin-1 (10) were the most abundant hemostatic and hemorrhagic toxins. We identified 59 transcripts (34%) that belong to 12 actiniarian-specific gene families. Most common among these are auxiliary activity, astacin-like metalloprotease toxin M12A (13 transcripts) and nematocyst expressed protein, NEP-6 (13 transcripts), which together represent 54% of all putative actiniarian toxin transcripts. Among Heteractina, phospholipase A2 (PLA2) was found only in Hs. crispa, which showed more putative neurotoxin activity proteins, cysteine-rich venom proteins, and potassium channel toxins (type II & III) than M. doreensis.
We identified 135 putative toxin transcripts in the transcriptome assembly for M. doreensis; these belong to 25 toxin gene families. The transcriptome of M. doreensis is rich in hemostatic and hemorrhagic toxins (63/135, 46%) and auxiliary proteins (25/135, 18%), representing 64% of all transcripts recovered. Transcripts identified with venom protease (11) and Snaclec (10) gene families were the most represented transcripts among hemostatic and hemorrhagic toxins. We recovered 48 transcripts (38%) that were assigned to 12 actiniarian-specific gene families. The transcriptome of M. doreensis includes notable representation of transcripts identified as belonging to membrane active activity protein MAC/perforin domain (7 transcripts) and venom Kunitz-type protease inhibitor (9 transcripts) gene families. Among Heteractina, M. doreensis. was the only species to recover a transcript from the sea anemone 8 toxin family (1). We did not recover SCRiP, ShK-like-1 neurotoxin, sea anemone short toxin (type III), type II sodium channel toxins, or sea anemone structural class 9a toxins from the transcriptomes of Heteractina studied here.

2.2.3. Entacmaea

Entacmaea quadricolor is the only member of Entacmaea to host clownfish. We found 328 putative toxin transcripts in its transcriptome, representing 37 toxin gene families. The transcriptome of E. quadricolor includes hemostatic and hemorrhagic toxins (222/328, 67%), with the true venom lectin (26), ryncolin (35), and veficolin-1 (45) families heavily represented. Of the 328 transcripts, 79 (24%) belonged to 13 actiniarian-specific gene families, with strong representation of the auxiliary activity proteins astacin-like metalloprotease toxin M12A (22) and nematocyst expressed protein NEP-6 (21). Despite having the second-highest number of actiniarian-specific putative toxin transcripts, we recovered no neurotoxins (including no transcripts similar to any voltage gated potassium or sodium toxins) in the transcriptome of E. quadricolor.

2.3. Diversity of Toxins within Functional Groups

A primary objective of this study was to identify functional groupings of toxins and characterize their distribution and diversity across actinioidean lineages that have independently associated with clownfish. We found multiple clade-specific putative toxins in Stichodactylina, and one unique putative toxin each in Heteractina and Entacmaea. Collectively, these unique putative toxins were low in abundance and in every functional class except mixed enzymes and pore-forming toxins. Additionally, knowing that the clownfish-hosting anemones represent independent associations in terms of the host lineage [30], we have looked at shared transcripts within a gene family (Table 2) and also at groups of functionally similar toxins that may include genes from many gene families to understand whether there are similarities in the composition of the venom arsenal at a functional level, if not at a phylogenetic one.

2.3.1. Hemostatic and Hemorrhagic Toxins

Hemostatic and hemorrhagic toxins were the most diverse and prevalent types of putative venom genes recovered in this study. These toxins are associated with blood coagulation, inflammation, myotoxicity, platelet aggregation, and homeostasis interference and are ubiquitous in venom [48] and well documented among cnidarians [49,50,51,52].
Several hemostatic or hemorrhagic toxins were found across all studied taxa. We recovered abundant transcripts of hemotoxin venom prothrombin activator all species: C. adhaesivum (25), Ha. hemprichii (15), S. haddoni (5), M. doreensis (8), Hs. crispa (32), and E. quadricolor (48). Transcripts that closely match veficolin toxins that interfere in platelet aggregation [53] were recovered in all species, with the highest abundance in E. quadricolor (45). Transcripts matching the anti-coagulant ryncolin were found in all species and highly expressed in E. quadricolor (35). Hemorrhagic toxins similar to zinc metalloproteinase/disintegrins were recovered in all species. Peptide isomerase heavy chain transcripts were recovered across all species and heavily expressed in C. adhaesivum (17) and E. quadricolor (15).
We recovered transcripts matching American bumble bee venom protease activators [54,55] and Snaclec family, and c-type snake lectins [56] across all species. Transcripts similar to Snake venom lectins were recovered in all species but Ha. hemprichii. In all species but M. doreensis, we recovered transcripts closely matching blarina toxins (BLTX), which have dilatory effects on blood vessel walls in short tailed shrew [57]. Transcripts highly similar to snake venom coagulation factor V & X transcripts were also found in all taxa except M. doreensis. The procoagulant hemolin (EC 3.4.21.) toxin was not recovered in S. haddoni or Hs. crispa, but was found in C. adhaesivum, E. quadricolor, and Ha. hemprichii. Cobra venom factors (VF) from the venom complement C3 homolog family [58] were not recovered in S. haddoni or any members of Heteractina.
Putative toxin transcripts for snake venom coagulation factor V & X were recovered across all species except M. doreensis. We recovered transcripts with high similarity to two types of toad fish toxins: galactose-specific lectin nattectin [59,60], a pro-inflammatory activity toxin which also induces neutrophil mobilization, was recovered in all species but Ha. hemprichi. The proteolytic and myxotoxic toxin, natterin-4 [61] was recovered only in Hs. crispa (4) and E. quadricolor (1).

2.3.2. Neurotoxins/Protease Inhibitors

Cnidarian neurotoxins are a set of diverse and well-characterized sodium and potassium channels toxins used to both immobilize the prey and to defend against predators [62]. They are widely used in the development of pharmaceuticals and bioinsecticides [63,64,65,66]. In this study, neurotoxins represented the second-most diverse type of toxins recovered in focal taxa. We found sea anemone type I and III potassium channel toxin subfamily (KTx1, KTx3) transcripts across all clades, with transcript abundances ranging from 1–9 (Table 3). The only species in which we found no KTx1 transcripts is Ha. hemprichii.
The cnidarian-specific neuropeptide SCRiP [8], originally described in reef building corals, was recovered in only Ha. hemprichii (1) and E. quadricolor (1). Cysteine-rich venom proteins were recovered in C. adhaesivum and in both members of Heteractina. Sea anemone structural class 9a transcripts were recovered in the stichodactylines C. adhaesivum (1) and S. haddoni (1) and in E. quadricolor (3).
In all species, we recovered transcripts matching to the conopeptide p-like superfamily. In cone snails, these proteins act both as neurotoxin and serine protease inhibitor [67]; they were found in high relative abundance in both C. adhaesivum (10) and E. quadricolor (8). We found transcripts closely matching neurotoxic u-scoloptoxin [68] in S. haddoni (1) and Hs. crispa (2).
Kunitz-domain peptides block ion channels and inhibit protease, leading to blood coagulation and inflammation [69]. Transcripts matching the venom Kunitz-type family were found across all species, with Heteractina member M. doreensis recovering the greatest number of transcripts (9). Sea anemone type II potassium channel toxin (KTx2) is a neurotoxin that has Kunitz domain; transcripts belonging to this group were abundant in all species, being most abundant in C. adhaesivum (11) and E. quadricolor (9).
Only in Stichodactylina did we recover transcripts identified as sea anemone sodium channel toxin type II (NaTx2). Both S. haddoni and Ha. hemprichii had transcripts matching the sea anemone short toxin (potassium type III) family. In Ha. hemprichii, we found the only ShK-like-1 neurotoxin transcript, similar to what has previously been described as type II from the model actiniarian Nematostella vectensis [8,70]. A transcript closely matching the protein omega-theraphotoxin-Pm1a was only recovered in S. haddoni. Additionally, a single transcript closely matching psalmotoxin-1 (PcTx1) was only recovered in S. haddoni. Psalmotoxins have pharmacological neuroprotection applications [71].

2.3.3. Auxiliary Venom

Auxiliary toxin proteins act as venom stabilizers and are known to induce hemorrhaging and cause necrosis by degrading the extracellular matrix [72]. We found cnidarian-specific nematocyst expressed protein (NEP-6) transcripts, first characterized in the nematocysts of the starlet sea anemone, Nematostella vectensis [73], across all species. The greatest number of copies of NEP-6 were in C. adhaesivum (24) and E. quadricolor (21).
Across all species, we recovered astacin-like metalloprotease toxin (M12A). We recover transcripts closely matching the cysteine-type endopeptidase inhibitor cystatin-2 in all taxa except the stichodactylines C. adhaesivum and Ha. hemprichii.
Crypotdendrum adhaesivum was the only species in which we recovered auxiliary venom protein 302. Venom protein 302 has been recovered in many taxa, including the zoantharian cnidarians Protopalythoa variabilis and model actiniarian Exaiptasia pallida (see [49,74]). In C. adhaesivum we recovered a singular copy of venom protein 59.1, a protein associated with insulin-like activity [75]. We also found a reticulocalbin-2 in Ha. hemprichii, similar to a taipoxin-associated calcium-binding protein [76].

2.3.4. Mixed Enzyme

PLA2s have been recruited to venom function by almost all venomous animals [77], acting as a multifunctional toxin with myotoxic, neurotoxic, and hemotoxic activity. We found PLA2s in the transcriptomes of S. haddoni (2), Hs. crispa (5), and E. quadricolor (2).

2.3.5. Pore Forming Toxins

Pore forming toxins are common in cnidarian venoms [78] and disrupt or penetrate cell membranes via lysis [79]. Pore forming toxins that form α-helical barrel structures are referred to as α-PFTs or actinoporins [80]. We recovered actinoporin transcripts in all species except for S. haddoni. We found actinoporins in greatest abundance in the heteractine Hs. crispa (5). Pore forming toxins that form β-barrel pores are referred to as (β-PFTs) or Membrane Attack Complex/Perforins (MAC/PF) and are found across eukaryotes [81]. As was the case for actinoporins, we recovered transcripts for MAC/PFs in all species except S. haddoni, with highest abundance (7) in the heteractine M. doreensis.

2.3.6. Allergens and Innate Immunity

Transcripts classified as allergens and innate immunity genes are common in the venom arsenals of arthropod and reptile toxins [82,83] and have been recovered across Cnidaria [50,84]. The stichodactyline C. adhaesivum was the only species in this study in which we recover a transcript closely matching snake venom serine protease Dav-PA (3). This protein has fibrinogenolytic, esterolysis, and amidolytic activities [85]. In E. quadricolor, we recovered a single transcript closely matching a techylectin-like protein known to lead to platelet aggregation and blood coagulation [53].
Across all transcriptomes we recovered transcripts for the cysteine-rich venom protein, venom allergen 5. Venom serine protease Bi-VSP, a multifunctional enzyme which was previously found in both wasps and in bees [54,55] was also recovered across all species. Entacmaea quadricolor was the only species in which we did not recover the silk moth venom serine protease and Ha. hemprichii was the only species in which we did not recover transcripts that align to the gene for venom serine carboxypeptidase previously identified in bees and inferred to be involved in the degradation of [86]. We recovered transcripts matching the complement-activating protein Austrelaps superbus venom factor 1 [87] in all taxa except S. haddoni and M. doreensis.

2.3.7. Venom Transcrits of Unknown Function

We find transcripts similar to those associated with venom or toxins in other lineages, but for which precise function or mode of action is unknown. Several of these are part of toxin gene families that are part of physiological pathways and systems other than venom. For example, Endothelial Growth Factors or EGF domain peptides like OMEGA-stichotoxin-Shd4a may have both toxic and EGF activity. We found OMEGA-stichotoxins in all species studied here. We found transcripts matching the sea anemone 8 toxin family [88,89,90] across three species: C. adhaesivum, M. doreensis, and E. quadricolor.

3. Discussion

Actiniarians are predatory invertebrates that use venom to catch prey. Previous research on actinarian venoms has focused on describing and characterizing these toxins broadly, and those studies have shown that Actinarian venom composition is largely neurotoxin-rich [3,91] in contrast to the venom of other cnidarians, such as the enzymatically rich toxin arsenal of medusozoans [3,92,93]. Given the presumed defensive function of neurotoxins for their symbiotic clownfish, we expected greater prevalence of neurotoxins in the venom transcriptome of clownfish-hosting sea anemones [20]. Instead, from our findings we conclude that, when looking at the whole venom arsenal reconstructed from transcriptomes, the venom composition of clownfish-hosting sea anemones largely emphasizes predatory rather than defensive behavior.
The prevalence of hemostatic and hemorrhagic transcripts in the expressed venom of clownfish hosting actiniarians is the rationale for interpreting their venom as largely predatory. Within this functional group, we found key gene families which were highly abundant and consistent across all species: for example, transcripts that match platelet aggregation disruptors veficolin-1 and ryncolin were found in high numbers across all species, as were multiple members of the subfamily of venom prothrombin activators. The ecology of clownfish-hosting anemones is likely relevant to understanding their venom [94] suggests that the hemostatic and hemorrhagic “predatory” venom is dual function and important in defense of the fish. Dual function for predatory venom is common [94] and is likely relevant here because the sea anemones themselves are expected to rely minimally on predation, using zooxanthellae to meet most of their nutrition [20]. Future studies of the target and effect of these hemostatic and hemorrhagic venoms should evaluate whether it is more attuned to the diet of the sea anemone host or the predators of the clownfish.
The prevalence of hemostatic and hemorrhagic transcripts in the expressed venom of clownfish hosting actiniarians is the rationale for interpreting their venom as largely predatory. Within this functional group, we found key gene families which were highly abundant and consistent across all species: for example, transcripts that match platelet aggregation disruptors veficolin-1 and ryncolin were found in high numbers across all species, as were multiple members of the subfamily of venom prothrombin activators. The ecology of clownfish-hosting anemones is likely relevant to understanding their venom [94] and may explain differences from related but ecologically dissimilar species. Dual function for predatory venom is common [94], and hemolytic and hemorrhagic toxins may be important in the defense of the fish. Venom from a species which is invoved in complex evolutionary interactions is likely to have its venom composition shaped by those interactions [1,48]. Sea anemones are generalist predators and thus the breadth of predators of clownfish, rather than the breadth of the host anemone’s diet, would be expected to shape the venom in these species.
The dynamic interactions which occur between clownfish and their venomous hosts places clownfish-hosting sea anemones in a unique category of venomous animals with pharmaceutical potential [62,95]. Among sea anemone-derived toxins, one (a ShK derivative called dalazatide) has passed initial pharmaceutical testing and is currently being used in human trials to combat multiple sclerosis and autoimmune disorders [96,97]. Given the long evolutionary history of Cnidaria and symbiotic association with potential vertebrate targets, clownfish-hosting anemones present a unique system in which toxins with pharmaceutical potential have undergone multualistic co-evolution with vertebrates, rather than co-evolving as predators or prey. Although it is unlikely every putative toxin gene identified here is used in the sea anemone venom arsenal, these accounts of the putative toxin repertore of clownshigh-hosting anemones serve as a starting point to identify additional toxin candidates with pharmaceutical potential.
Beyond what has been previously characterized in these species, we found eight genes affiliated with toxins that have not been reported in actiniarians. These span functional groups and include members of gene families previously reported in cnidarians and those known only from other venomous lineages. The insulin-like venom proteins 302 and 59.1 that we identify in the transcriptome of C. adhaesivum and Ha. hemprichii, respectively, have not been reported in any other actiniarians, but have been found in other cnidarians [51,98,99,100]. Insulin-like peptides reported from the actiniarians Oulactis sp., N. vectensis, and Ex. pallida [101,102] are cono-insulins either matching to human insulins or cone snail toxins; the relationships between these cono-insulin-like peptides and insulin-like venom proteins 302 or 59.1 are unknown. We did not recover any cono-insulin-like peptides using our pipeline. Like the insulin-like peptides, reticulocalbin is known in cnidarians [98] but not in actiniarians. In contrast, the neurotoxins omega-theraphotoxin-Pm1a and psalmotoxin-1, originally described in tarantulas, have not been reported in any cnidariana. Although hemotoxins were diverse and abundant in the transcriptomes of all taxa (Table 3), snake venom phosphodiesterase was recovered only in Hs. crispa and snake venom serine protease (SVSP) Dav-PA was recovered only in C. adhaesivum.
No toxin transcript previously reported in anemones was uniquely identified in any of the study species, although some of the more diverse toxins were present in low copy number, occured in only a couple of taxa, or were restricted to a clade (Table 3). Based on the diversity of toxins reported to date for actiniarians, we expect to find Type 1, 2, and 3 sodium channel toxins; types 1, 3, 4, and 5 potassium channel toxins; cytolysin 1 and 3; NEPs; SCRiPS; acrorhagins; and actinoporins (see [4]). The fact that we did not recover all of these in all taxa (Table 3) indicates that transcripts of those genes are either not expressed or expressed at very low levels when the data were collected, or that they are absent from the venom of these species altogether. Of the types expected to be present, acrorhagins are an interesting omission: although these toxins can kill crustaceans [103], they are associated with intraspecific aggression in actinioideans and metridiodeans by virtue of their pattern of expression in tissues involved in self/not self-interactions [104,105,106] and so are expected to be part of the shared repertoire of both clades [4]. However, no clownfish-hosting anemone has structures associated with the use of acrorhagins in other taxa, and none are known to engage in intraspecific competitive interactions. Thus, the lack of expression may be expected based on function, if not on phylogeny. More information about the function of these toxins and a genomic perspective on acrorhagins in species that do and do not engage in intraspecific competition is needed to understand the evolution of these toxin types.
Species-specificity of a toxin implies either a related lineage loss or pseudogenization, or an independent recruitment event involving a protein gaining a toxic function. Because each transcriptome is a snapshot in the expression profile, absence may simply reflect the developmental or physiological state of the animal sampled and should be contextualized with an understanding of how these factors impact venome [12,107,108,109,110]. Although our use of a standard method for assembly, search, and identification of putative venom transcripts mitigates some problems with comparing venom across transcriptomes, given the lack of standardization and replication of transcriptome generation, absent toxins should be interpreted cautiously [111].However, unique occurrences and absences are provocative, indicating novelties that can (and should) be searched for deliberately in subsequent transcriptomic or proteomic studies.
Although venom is dynamic and expression can change across time, space, and clade, we are encouraged and intrigued to see stability across functional groups in lineages that share similar ecologies. All studied transcriptomes were dominated by transcripts matching those known to have hemostatic and hemorrhagic function. Consequently, we infer that this functional group is of primary importance in these species. The fraction of toxins identified as hemostatic or hemorrhagic per transcriptome ranged from 46–66% of total toxins recovered across our study species. These were the most diverse type of toxin in terms of the number of toxin families (Table 3). The second-most prevalent toxin type varies by clade, auxiliary toxins second in Heteractina and neurotoxins being second in both stychodactylines and in Entacmaea. Macrodactyla doreensis showed a relatively high abundance of pore forming toxins, consistent with Ashwood et al. [45], which may indicate the particular importance of pore forming toxins in the venom arsenal of M. doreensis.
A few published accounts can help contextualize what we report for the clownfish-hosting anemones. In cerithiarians [98], hemostatic and hemorrhagic toxins are the dominant constituents. In Nematostella [110] and in Anthopleura elegantissima [105], neurotoxins are the most abundant reported toxin, and potentially the most abundant among ancestral actinarians [110]. The proteome of mucus of envenomating Anthopleura dowii is also especially rich in neurotoxins, although its tentacle transcriptome is relatively enriched in hemorrhagic toxins (constituting >50% of the reported toxins) [112], which may be due to differences in molecular stability, tissue specific expression (with mucus representing a pooled sample), or individual-level differences in physiological state and expression at any point in time. In N. vectensis, the voltage gated neurotoxins are made in the ectodermal gland cells whereas the hemorrhagic toxin metalloprotease is made in the nematocyte [73]. This offers a possible explanation for and the relative abundance and control of abundance for various kinds of toxins, since ectodermal gland cells and nematocytes occur in different densities and arrangements in different parts of the animal.
Transcriptome assembly and completeness is a major determinant of the venom repertoire that can be inferred [111,113]. Because of differences in assembly parameters and growth of databases against which the transcriptomes are queried, we expect to recover transcripts that could not be interpreted or compared in earlier assemblies (Table 1). Furthermore, methodological differences in transcriptome queries make it hard to compare results across studies. Nonetheless, despite differences in assembly and match criteria, our results are broadly consistent with previous studies of cnidarian venoms in our recovery of an array of transcripts associated with the types of venoms found in actiniarians: auxiliary, pore forming, mixed enzyme, neurotoxins, protease inhibitor activities and anemone toxins with unknown functions. In addition to the categories typically used for describing actiniarian venom, we identified allergens and innate immunity toxins and hemostatic and hemorrhagic activity toxins to better reflect functionality of venom in these species.
In general, across the focal taxa, we find that most putative toxin groups are represented by a small number of transcripts, so a greater number of transcripts in the assembly broadly translates into a greater number of gene families. Increasing the number of taxa studied within a lineage did not necessarily increase the perceived venom diversity or the number of toxin families reported for that lineage (Table 2). For example, in the single representative of Entacmaea (E. quadricolor), we recovered more putative toxins and gene families than in the two species of Heteractina combined. Although the number of actiniarian-specific toxins varies considerably, the number of actiniarian-specific gene families varies relatively little across samples (range is 12–13, see Table 2). The overall similarity in the percentage of the recovered toxins that belong to actiniarian-specific gene families obscures diversity across taxa in terms of which gene families are present. Furthermore, the proportion of actiniarian-specific transcripts is wide-ranging, with more than double the proportion of actiniarian-specific transcripts in M. doreensis than in E. quadricolor or in Ha. hemprichii.
To properly evaluate the pharmaceutical potential for sea anemones that host vertebrate symbionts, it is important to consider all possible toxin candidates found within a transcriptome or proteome. Further exploration is needed with to properly evaluate the medical potential with with these understudied proteins, which have been shaped by evolution over millions of years.

4. Materials and Methods

We analyzed transcriptomes of six species of sea anemones from NCBI (Table 1) using all available data, including transcriptomes from specific tissues rather than the whole body (Table S1). We compiled transcriptomes for Cryptodendrum adhaesivum from run IDs SRR14115232, SRR14115233; Stichodactyla haddoni from the run Id SRR5397293, Heterodactyla hemprichii from run IDs SRR14115226, SRR14115227, SRR14115228, for Heteractis crispa from run IDI SRR1950632, SRR1950633, SRR1950656; for SRR14115229; for Macrodactyla doreensis for run IDs SRR14115222, SRR14115223, SRR14115224, SRR14115225; and for Entacmaea quadricolor from run IDs ERR2045166, ERR2045167, ERR2045168, RR2045169, ERR2045170, EERR204517. Transcriptome construction followed a de novo assembly using Trinity v2.2 [114] using the trimmomatic option, on the Ohio State Supercomputer Center (OSC) running on 2 nodes with 30 CPU each. Transcriptome completeness was determined via BUSCO v5.3.1 [115] against the metazoa_odb10 lineage dataset for each transcriptome. Each transcriptome was bioinformatically annotated for putative venom transcripts. Results per species were then pooled for each species where there were multiple transcriptomes.
Protein-coding regions were predicted from assembled transcriptomes using TransDecoder v5.5.0 (https://transdecoder.github.io, accessed on November 2021), with a minimum sequence length of 50 amino acids. All transcriptomes were queried against the Tox-prot animal venom database using the key word (“Cnidaria,” downloaded November 2021) and the NCBI Protein Database using the key words (“Cnidaria AND ((Toxin) OR (Venom)),” downloaded 23 November 2021) by using blastp from NCBI BLAST + v.2.8.1 [116] using an e-value cutoff of 0.001. Predicted protein-coding region were searched using hmmscan with an e-value cutoff of 0.001 from HMMER 3.1b2 [117] against hidden markov model (HMM) profiles alignments of all the cnidarian venom protein classes from VenomZone (venomzone.expasy.org, accessed on November 2021). The results from all three queries (ToxProt, cnidaria specific-NCBI, and hmmsearch) were combined for each species and only complete coding sequences used for downstream analysis.
Because venoms are secreted proteins and peptides, signal peptides were predicted from the sequences of the combined queries using the SignalP v5.0 server (https://services.healthtech.dtu.dk/service.php?SignalP-5.0, accessed on 23 November 2021) [118]. Sequences with a 70% >probability of having a signaling region were clustered using CD-HIT v.4.6.8 [119] with a cutoff of 0.95. Only the top hit from each cluster was used in downstream analysis. To ensure completeness of venom recovery, sequences which passed these thresholds were used for a reciprocal search against a concatenated transcriptome form each species. The resulting datasets (signal peptide present, with redundant sequences removed) were screened via Blastp with an e-value cutoff of 1 × 10−5 against Tox-Prot animal venom database and the NCBI non-redundant protein sequences database (NR DB) (downloaded March 2022), and via hmmsearch (PFAM) search with an evalue cutoff of 1 × 10−5 against Pfam (downloaded March 2022). Manual curation of query hits was performed to confirm that annotations from ToxProt and NR DB matched the venom domain from Pfam. The final list of Toxin candidates for each species were classified into protein families and grouped by toxin function. Pipeline can be visualized in Figure S1 (Supplementary Material).
Select venoms were further evaluated by creating alignments using MUSCLE [120] on the EMBL-EBI web server [121] (https://www.ebi.ac.uk/Tools/msa/muscle/, accessed on 23 April 2022) aligning to proteins found on UniProt [122]. Alignment visualization was created on Geneious Prime v2022.2.2 (https://www.geneious.com, accessed on 23 April 2022) (Figures S2–S5).

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/md20120730/s1, Figure S1: Bioinformatic pipeline for the annotation of venom-like genes, Figures S2–S5: alignments of: NaTx1&2, SHk-like-1, venom 302, sea anemone short type 3 toxins to their gene families, Table S1: Individual transcriptome trinity statistics and BUSCO scores, Table S2: List of Identified venom transcripts (Excel).

Author Contributions

Conceptualization, J.M. and M.D.; Data curation, A.D.; Formal analysis, A.D. and C.B.; Investigation, A.D.; Methodology, A.D. and M.D.; Writing—original draft, A.D., C.B., J.M. and M.D.; Writing—review & editing, A.D., C.B., J.M. and M.D.. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We are grateful to Anna M. L. Klompen, Micheal Broe, and Katie Statile for discussions about the study. We would also like to acknowledge the reviewers who improved the initial submition of this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Casewell, N.R.; Wüster, W.; Vonk, F.J.; Harrison, R.A.; Fry, B.G. Complex Cocktails: The Evolutionary Novelty of Venoms. Trends Ecol. Evol. 2013, 28, 219–229. [Google Scholar] [CrossRef]
  2. Moran, Y.; Genikhovich, G.; Gordon, D.; Wienkoop, S.; Zenkert, C.; Özbek, S.; Technau, U.; Gurevitz, M. Neurotoxin Localization to Ectodermal Gland Cells Uncovers an Alternative Mechanism of Venom Delivery in Sea Anemones. Proc. R. Soc. B Biol. Sci. 2012, 279, 1351–1358. [Google Scholar] [CrossRef] [Green Version]
  3. Rachamim, T.; Morgenstern, D.; Aharonovich, D.; Brekhman, V.; Lotan, T.; Sher, D. The Dynamically Evolving Nematocyst Content of an Anthozoan, a Scyphozoan, and a Hydrozoan. Mol. Biol. Evol. 2015, 32, 740–753. [Google Scholar] [CrossRef] [Green Version]
  4. Daly, M. Functional and genetic diversity of toxins in sea anemones. In Evolution of Venomous Animals and Their Toxins; Gopalakrishnakone, P., Malhotra, A., Eds.; Toxinology; Springer: Dordrecht, The Netherlands, 2016; pp. 1–18. ISBN 978-94-007-6727-0. [Google Scholar]
  5. Fautin, D.G. Structural Diversity, Systematics, and Evolution of Cnidae. Toxicon 2009, 54, 1054–1064. [Google Scholar] [CrossRef] [PubMed]
  6. Schlesinger, A.; Zlotkin, E.; Kramarsky-Winter, E.; Loya, Y. Cnidarian Internal Stinging Mechanism. Proc. R. Soc. B Biol. Sci. 2009, 276, 1063–1067. [Google Scholar] [CrossRef] [Green Version]
  7. Nevalainen, T.J.; Peuravuori, H.J.; Quinn, R.J.; Llewellyn, L.E.; Benzie, J.A.H.; Fenner, P.J.; Winkel, K.D. Phospholipase A2 in Cnidaria. Comp. Biochem. Physiol. B Biochem. Mol. Biol. 2004, 139, 731–735. [Google Scholar] [CrossRef]
  8. Jouiaei, M.; Sunagar, K.; Gross, A.F.; Scheib, H.; Alewood, P.F.; Moran, Y.; Fry, B.G. Evolution of an Ancient Venom: Recognition of a Novel Family of Cnidarian Toxins and the Common Evolutionary Origin of Sodium and Potassium Neurotoxins in Sea Anemone. Mol. Biol. Evol. 2015, 32, 1598–1610. [Google Scholar] [CrossRef] [Green Version]
  9. Suggett, D.J.; Hall-Spencer, J.M.; Rodolfo-Metalpa, R.; Boatman, T.G.; Payton, R.; Tye Pettay, D.; Johnson, V.R.; Warner, M.E.; Lawson, T. Sea Anemones May Thrive in a High CO2 World. Glob. Change Biol. 2012, 18, 3015–3025. [Google Scholar] [CrossRef] [PubMed]
  10. Daly, M.; Rack, F.; Zook, R. Edwardsiella andrillae, a New Species of Sea Anemone from Antarctic Ice. PLoS ONE 2013, 8, e83476. [Google Scholar] [CrossRef] [Green Version]
  11. Menezes, C.; Thakur, N.L. Sea Anemone Venom: Ecological Interactions and Bioactive Potential. Toxicon 2022, 208, 31–46. [Google Scholar] [CrossRef]
  12. Sachkova, M.Y.; Macrander, J.; Surm, J.M.; Aharoni, R.; Menard-Harvey, S.S.; Klock, A.; Leach, W.B.; Reitzel, A.M.; Moran, Y. Some like It Hot: Population-Specific Adaptations in Venom Production to Abiotic Stressors in a Widely Distributed Cnidarian. BMC Biol. 2020, 18, 121. [Google Scholar] [CrossRef]
  13. LaJeunesse, T.C.; Parkinson, J.E.; Gabrielson, P.W.; Jeong, H.J.; Reimer, J.D.; Voolstra, C.R.; Santos, S.R. Systematic Revision of Symbiodiniaceae Highlights the Antiquity and Diversity of Coral Endosymbionts. Curr. Biol. 2018, 28, 2570–2580. [Google Scholar] [CrossRef] [Green Version]
  14. Goffredi, S.K.; Motooka, C.; Fike, D.A.; Gusmão, L.C.; Tilic, E.; Rouse, G.W.; Rodríguez, E. Mixotrophic Chemosynthesis in a Deep-Sea Anemone from Hydrothermal Vents in the Pescadero Basin, Gulf of California. BMC Biol. 2021, 19, 8. [Google Scholar] [CrossRef]
  15. Williams, J.D.; McDermott, J.J. Hermit Crab Biocoenoses: A Worldwide Review of the Diversity and Natural History of Hermit Crab Associates. J. Exp. Mar. Biol. Ecol. 2004, 305, 1–128. [Google Scholar] [CrossRef]
  16. Mercier, A.; Hamel, J.-F. Nature and Role of Newly Described Symbiotic Associations between a Sea Anemone and Gastropods at Bathyal Depths in the NW Atlantic. J. Exp. Mar. Biol. Ecol. 2008, 358, 57–69. [Google Scholar] [CrossRef]
  17. Sanamyan, N.P.; Sanamyan, K.E.; Tabachnik, K.R. The First Species of Actiniaria, Spongiactis japonica gen.n., sp.n. (Cnidaria: Anthozoa), an Obligate Symbiont of a Glass Sponge. Invertebr. Zool. 2012, 9, 127–141. [Google Scholar] [CrossRef]
  18. Izumi, T.; Ise, Y.; Yanagi, K.; Shibata, D.; Ueshima, R. First Detailed Record of Symbiosis Between a Sea Anemone and Homoscleromorph Sponge, With a Description of Tempuractis rinkai Gen. et Sp. Nov. (Cnidaria: Anthozoa: Actiniaria: Edwardsiidae). Zoolog. Sci. 2018, 35, 188–198. [Google Scholar] [CrossRef] [Green Version]
  19. Vader, W.; Tandberg, A.H.S. Amphipods and Sea Anemones, an Update. J. Crustac. Biol. 2020, 40, 872–878. [Google Scholar] [CrossRef]
  20. Fautin, D.G. The Anemonefish Symbiosis: What Is Known and What Is Not. Symbiosis 1991, 10, 23–26. [Google Scholar]
  21. Dunn, D.F. The Clownfish Sea Anemones: Stichodactylidae (Coelenterata: Actiniaria) and Other Sea Anemones Symbiotic with Pomacentrid Fishes. Trans. Am. Philos. Soc. 1981, 71, 3–115. [Google Scholar] [CrossRef]
  22. Godwin, J.; Fautin, D.G. Defense of Host Actinians by Anemonefishes. Copeia 1992, 1992, 902–908. [Google Scholar] [CrossRef]
  23. Holbrook, S.J.; Schmitt, R.J. Growth, Reproduction and Survival of a Tropical Sea Anemone (Actiniaria): Benefits of Hosting Anemonefish. Coral Reefs 2005, 24, 67–73. [Google Scholar] [CrossRef]
  24. Porat, D.; Chadwick-Furman, N.E. Effects of Anemonefish on Giant Sea Anemones: Ammonium Uptake, Zooxanthella Content and Tissue Regeneration. Mar. Freshw. Behav. Physiol. 2005, 38, 43–51. [Google Scholar] [CrossRef]
  25. Alan Verde, E.; Cleveland, A.; Lee, R.W. Nutritional Exchange in a Tropical Tripartite Symbiosis II: Direct Evidence for the Transfer of Nutrients from Host Anemone and Zooxanthellae to Anemonefish. Mar. Biol. 2015, 162, 2409–2429. [Google Scholar] [CrossRef]
  26. Litsios, G.; Sims, C.A.; Wüest, R.O.; Pearman, P.B.; Zimmermann, N.E.; Salamin, N. Mutualism with Sea Anemones Triggered the Adaptive Radiation of Clownfishes. BMC Evol. Biol. 2012, 12, 212. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. McCord, C.L.; Nash, C.M.; Cooper, W.J.; Westneat, M.W. Phylogeny of the Damselfishes (Pomacentridae) and Patterns of Asymmetrical Diversification in Body Size and Feeding Ecology. PLoS ONE 2021, 16, e0258889. [Google Scholar] [CrossRef] [PubMed]
  28. Frédérich, B.; Sorenson, L.; Santini, F.; Slater, G.J.; Alfaro, M.E. Iterative Ecological Radiation and Convergence during the Evolutionary History of Damselfishes (Pomacentridae). Am. Nat. 2013, 181, 94–113. [Google Scholar] [CrossRef]
  29. Elliott, J.K.; Lougheed, S.C.; Bateman, B.; McPhee, L.K.; Boag, P.T. Molecular Phylogenetic Evidence for the Evolution of Specialization in Anemonefishes. Proc. R. Soc. Lond. B Biol. Sci. 1999, 266, 677–685. [Google Scholar] [CrossRef] [Green Version]
  30. Titus, B.M.; Benedict, C.; Laroche, R.; Gusmão, L.C.; Van Deusen, V.; Chiodo, T.; Meyer, C.P.; Berumen, M.L.; Bartholomew, A.; Yanagi, K.; et al. Phylogenetic Relationships among the Clownfish-Hosting Sea Anemones. Mol. Phylogenet. Evol. 2019, 139, 106526. [Google Scholar] [CrossRef]
  31. Lubbock, R. The Clownfish/Anemone Symbiosis: A Problem of Cellular Recognition. Parasitology 1981, 82, 159–173. [Google Scholar] [CrossRef]
  32. Allen, G.R. Anemonefishes and Their Amazing Partnership. Aust. Nat. Hist. 1975, 18, 274–277. [Google Scholar]
  33. Mebs, D. Chemical Biology of the Mutualistic Relationships of Sea Anemones with Fish and Crustaceans. Toxicon 2009, 54, 1071–1074. [Google Scholar] [CrossRef]
  34. Mebs, D. Anemonefish Symbiosis: Vulnerability and Resistance of Fish to the Toxin of the Sea Anemone. Toxicon 1994, 32, 1059–1068. [Google Scholar] [CrossRef] [PubMed]
  35. Elliott, J.K.; Mariscal, R.N. Coexistence of Nine Anemonefish Species: Differential Host and Habitat Utilization, Size and Recruitment. Mar. Biol. 2001, 138, 23–36. [Google Scholar] [CrossRef]
  36. Lubbock, R.; Smith, D.C. Why Are Clownfishes Not Stung by Sea Anemones? Proc. R. Soc. Lond. B Biol. Sci. 1980, 207, 35–61. [Google Scholar] [CrossRef]
  37. Mariscal, R.N. The Nature of the Symbiosis between Indo-Pacific Anemone Fishes and Sea Anemones. Mar. Biol. 1970, 6, 58–65. [Google Scholar] [CrossRef]
  38. Davenport, D.; Norris, K.S. Observations on the Symbiosis of the Sea Anemone Stoichactis and the Pomacentrid Fish, Amphiprion percula. Biol. Bull. 1958, 115, 397–410. [Google Scholar] [CrossRef]
  39. Elliott, J.K.; Mariscal, R.N.; Roux, K.H. Do Anemonefishes Use Molecular Mimicry to Avoid Being Stung by Host Anemones? J. Exp. Mar. Biol. Ecol. 1994, 179, 99–113. [Google Scholar] [CrossRef]
  40. Elliott, J.K.; Mariscal, R.N. Acclimation or Innate Protection of Anemonefishes from Sea Anemones? Copeia 1997, 1997, 284–289. [Google Scholar] [CrossRef]
  41. Marcionetti, A.; Rossier, V.; Roux, N.; Salis, P.; Laudet, V.; Salamin, N. Insights into the Genomics of Clownfish Adaptive Radiation: Genetic Basis of the Mutualism with Sea Anemones. Genome Biol. Evol. 2019, 11, 869–882. [Google Scholar] [CrossRef] [Green Version]
  42. Nedosyko, A.M.; Young, J.E.; Edwards, J.W.; Silva, K.B. da Searching for a Toxic Key to Unlock the Mystery of Anemonefish and Anemone Symbiosis. PLoS ONE 2014, 9, e98449. [Google Scholar] [CrossRef]
  43. Norton, R.S.; Pennington, M.W.; Wulff, H. Potassium Channel Blockade by the Sea Anemone Toxin ShK for the Treatment of Multiple Sclerosis and Other Autoimmune Diseases. Curr. Med. Chem. 2004, 11, 3041–3052. [Google Scholar] [CrossRef]
  44. Haas, B.J.; Papanicolaou, A.; Yassour, M.; Grabherr, M.; Blood, P.D.; Bowden, J.; Couger, M.B.; Eccles, D.; Li, B.; Lieber, M.; et al. De Novo Transcript Sequence Reconstruction from RNA-Seq: Reference Generation and Analysis with Trinity. Nat. Protoc. 2013, 8, 1494–1512. [Google Scholar] [CrossRef]
  45. Ashwood, L.M.; Mitchell, M.L.; Madio, B.; Hurwood, D.A.; King, G.F.; Undheim, E.A.B.; Norton, R.S.; Prentis, P.J. Tentacle Morphological Variation Coincides with Differential Expression of Toxins in Sea Anemones. Toxins 2021, 13, 452. [Google Scholar] [CrossRef] [PubMed]
  46. Madio, B.; Undheim, E.A.B.; King, G.F. Revisiting Venom of the Sea Anemone Stichodactyla haddoni: Omics Techniques Reveal the Complete Toxin Arsenal of a Well-Studied Sea Anemone Genus. J. Proteom. 2017, 166, 83–92. [Google Scholar] [CrossRef] [PubMed]
  47. Macrander, J.; Broe, M.; Daly, M. Tissue-Specific Venom Composition and Differential Gene Expression in Sea Anemones. Genome Biol. Evol. 2016, 8, 2358–2375. [Google Scholar] [CrossRef] [Green Version]
  48. Fry, B.G.; Roelants, K.; Champagne, D.E.; Scheib, H.; Tyndall, J.D.A.; King, G.F.; Nevalainen, T.J.; Norman, J.A.; Lewis, R.J.; Norton, R.S.; et al. The Toxicogenomic Multiverse: Convergent Recruitment of Proteins Into Animal Venoms. Annu. Rev. Genom. Hum. Genet. 2009, 10, 483–511. [Google Scholar] [CrossRef] [Green Version]
  49. Huang, C.; Morlighem, J.-É.R.; Zhou, H.; Lima, É.P.; Gomes, P.B.; Cai, J.; Lou, I.; Pérez, C.D.; Lee, S.M.; Rádis-Baptista, G. The Transcriptome of the Zoanthid Protopalythoa variabilis (Cnidaria, Anthozoa) Predicts a Basal Repertoire of Toxin-like and Venom-Auxiliary Polypeptides. Genome Biol. Evol. 2016, 8, 3045–3064. [Google Scholar] [CrossRef] [Green Version]
  50. Ponce, D.; Brinkman, D.L.; Potriquet, J.; Mulvenna, J. Tentacle Transcriptome and Venom Proteome of the Pacific Sea Nettle, Chrysaora fuscescens (Cnidaria: Scyphozoa). Toxins 2016, 8, 102. [Google Scholar] [CrossRef] [Green Version]
  51. Liao, Q.; Feng, Y.; Yang, B.; Lee, S.M.-Y. Cnidarian Peptide Neurotoxins: A New Source of Various Ion Channel Modulators or Blockers against Central Nervous Systems Disease. Drug Discov. Today 2019, 24, 189–197. [Google Scholar] [CrossRef]
  52. Li, R.; Yu, H.; Xue, W.; Yue, Y.; Liu, S.; Xing, R.; Li, P. Jellyfish Venomics and Venom Gland Transcriptomics Analysis of Stomolophus meleagris to Reveal the Toxins Associated with Sting. J. Proteom. 2014, 106, 17–29. [Google Scholar] [CrossRef]
  53. OmPraba, G.; Chapeaurouge, A.; Doley, R.; Devi, K.R.; Padmanaban, P.; Venkatraman, C.; Velmurugan, D.; Lin, Q.; Kini, R.M. Identification of a Novel Family of Snake Venom Proteins Veficolins from Cerberus rynchops Using a Venom Gland Transcriptomics and Proteomics Approach. J. Proteome Res. 2010, 9, 1882–1893. [Google Scholar] [CrossRef]
  54. Choo, Y.M.; Lee, K.S.; Yoon, H.J.; Kim, B.Y.; Sohn, M.R.; Roh, J.Y.; Je, Y.H.; Kim, N.J.; Kim, I.; Woo, S.D.; et al. Dual Function of a Bee Venom Serine Protease: Prophenoloxidase-Activating Factor in Arthropods and Fibrin(Ogen)Olytic Enzyme in Mammals. PLoS ONE 2010, 5, e10393. [Google Scholar] [CrossRef]
  55. Choo, Y.M.; Lee, K.S.; Yoon, H.J.; Qiu, Y.; Wan, H.; Sohn, M.R.; Sohn, H.D.; Jin, B.R. Antifibrinolytic Role of a Bee Venom Serine Protease Inhibitor That Acts as a Plasmin Inhibitor. PLoS ONE 2012, 7, e32269. [Google Scholar] [CrossRef] [Green Version]
  56. Zhong, S.-R.; Jin, Y.; Wu, J.-B.; Chen, R.-Q.; Jia, Y.-H.; Wang, W.-Y.; Xiong, Y.-L.; Zhang, Y. Characterization and Molecular Cloning of Dabocetin, a Potent Antiplatelet C-Type Lectin-like Protein from Daboia russellii siamensis Venom. Toxicon 2006, 47, 104–112. [Google Scholar] [CrossRef]
  57. Kita, M.; Nakamura, Y.; Okumura, Y.; Ohdachi, S.D.; Oba, Y.; Yoshikuni, M.; Kido, H.; Uemura, D. Blarina Toxin, a Mammalian Lethal Venom from the Short-Tailed Shrew Blarina brevicauda: Isolation and Characterization. Proc. Natl. Acad. Sci. USA 2004, 101, 7542–7547. [Google Scholar] [CrossRef] [Green Version]
  58. Rokyta, D.R.; Lemmon, A.R.; Margres, M.J.; Aronow, K. The Venom-Gland Transcriptome of the Eastern Diamondback Rattlesnake (Crotalus adamanteus). BMC Genom. 2012, 13, 312. [Google Scholar] [CrossRef] [Green Version]
  59. Lopes-Ferreira, M.; Magalhães, G.S.; Fernandez, J.H.; Junqueira-de-Azevedo, I.d.L.M.; Le Ho, P.; Lima, C.; Valente, R.H.; Moura-da-Silva, A.M. Structural and Biological Characterization of Nattectin, a New C-Type Lectin from the Venomous Fish Thalassophryne nattereri. Biochimie 2011, 93, 971–980. [Google Scholar] [CrossRef] [Green Version]
  60. Junqueira-de-Azevedo, I.d.L.M.; Ho, P.L. A Survey of Gene Expression and Diversity in the Venom Glands of the Pitviper Snake Bothrops insularis through the Generation of Expressed Sequence Tags (ESTs). Gene 2002, 299, 279–291. [Google Scholar] [CrossRef]
  61. Magalhães, G.S.; Lopes-Ferreira, M.; Junqueira-de-Azevedo, I.d.L.M.; Spencer, P.J.; Araújo, M.S.; Portaro, F.C.V.; Ma, L.; Valente, R.H.; Juliano, L.; Fox, J.W.; et al. Natterins, a New Class of Proteins with Kininogenase Activity Characterized from Thalassophryne nattereri Fish Venom. Biochimie 2005, 87, 687–699. [Google Scholar] [CrossRef]
  62. Prentis, P.J.; Pavasovic, A.; Norton, R.S. Sea Anemones: Quiet Achievers in the Field of Peptide Toxins. Toxins 2018, 10, 36. [Google Scholar] [CrossRef] [Green Version]
  63. Bosmans, F.; Tytgat, J. Sea Anemone Venom as a Source of Insecticidal Peptides Acting on Voltage-Gated Na+ Channels. Toxicon 2007, 49, 550–560. [Google Scholar] [CrossRef]
  64. Upadhyay, S.K.; Eckel-Mahan, K.L.; Mirbolooki, M.R.; Tjong, I.; Griffey, S.M.; Schmunk, G.; Koehne, A.; Halbout, B.; Iadonato, S.; Pedersen, B.; et al. Selective Kv1.3 Channel Blocker as Therapeutic for Obesity and Insulin Resistance. Proc. Natl. Acad. Sci. USA 2013, 110, E2239–E2248. [Google Scholar] [CrossRef] [Green Version]
  65. Yan, F.; Cheng, X.; Ding, X.; Yao, T.; Chen, H.; Li, W.; Hu, S.; Yu, Z.; Sun, Y.; Zhang, Y.; et al. Improved Insecticidal Toxicity by Fusing Cry1Ac of Bacillus Thuringiensis with Av3 of Anemonia viridis. Curr. Microbiol. 2014, 68, 604–609. [Google Scholar] [CrossRef]
  66. Chi, V.; Pennington, M.W.; Norton, R.S.; Tarcha, E.J.; Londono, L.M.; Sims-Fahey, B.; Upadhyay, S.K.; Lakey, J.T.; Iadonato, S.; Wulff, H.; et al. Development of a Sea Anemone Toxin as an Immunomodulator for Therapy of Autoimmune Diseases. Toxicon 2012, 59, 529–546. [Google Scholar] [CrossRef] [Green Version]
  67. Watkins, M.; Hillyard, D.R.; Olivera, B.M. Genes Expressed in a Turrid Venom Duct: Divergence and Similarity to Conotoxins. J. Mol. Evol. 2006, 62, 247–256. [Google Scholar] [CrossRef]
  68. Ward, M.J.; Rokyta, D.R. Venom-Gland Transcriptomics and Venom Proteomics of the Giant Florida Blue Centipede, Scolopendra viridis. Toxicon 2018, 152, 121–136. [Google Scholar] [CrossRef]
  69. Undheim, E.A.B.; Sunagar, K.; Herzig, V.; Kely, L.; Low, D.H.W.; Jackson, T.N.W.; Jones, A.; Kurniawan, N.; King, G.F.; Ali, S.A.; et al. A Proteomics and Transcriptomics Investigation of the Venom from the Barychelid Spider Trittame loki (Brush-Foot Trapdoor). Toxins 2013, 5, 2488–2503. [Google Scholar] [CrossRef] [Green Version]
  70. Moran, Y.; Gurevitz, M. When Positive Selection of Neurotoxin Genes Is Missing. FEBS J. 2006, 273, 3886–3892. [Google Scholar] [CrossRef]
  71. Baconguis, I.; Gouaux, E. Structural Plasticity and Dynamic Selectivity of Acid-Sensing Ion Channel–Spider Toxin Complexes. Nature 2012, 489, 400–405. [Google Scholar] [CrossRef] [Green Version]
  72. Laustsen, A.H. Toxin Synergism in Snake Venoms. Toxin Rev. 2016, 35, 165–170. [Google Scholar] [CrossRef] [Green Version]
  73. Moran, Y.; Praher, D.; Schlesinger, A.; Ayalon, A.; Tal, Y.; Technau, U. Analysis of Soluble Protein Contents from the Nematocysts of a Model Sea Anemone Sheds Light on Venom Evolution. Mar. Biotechnol. 2013, 15, 329–339. [Google Scholar] [CrossRef]
  74. Baumgarten, S.; Simakov, O.; Esherick, L.Y.; Liew, Y.J.; Lehnert, E.M.; Michell, C.T.; Li, Y.; Hambleton, E.A.; Guse, A.; Oates, M.E.; et al. The Genome of Aiptasia, a Sea Anemone Model for Coral Symbiosis. Proc. Natl. Acad. Sci. USA 2015, 112, 11893–11898. [Google Scholar] [CrossRef] [Green Version]
  75. Ruiming, Z.; Yibao, M.; Yawen, H.; Zhiyong, D.; Yingliang, W.; Zhijian, C.; Wenxin, L. Comparative Venom Gland Transcriptome Analysis of the Scorpion Lychas mucronatus Reveals Intraspecific Toxic Gene Diversity and New Venomous Components. BMC Genom. 2010, 11, 452. [Google Scholar] [CrossRef] [Green Version]
  76. Dodds, D.; Schlimgen, A.K.; Lu, S.-Y.; Perin, M.S. Novel Reticular Calcium Binding Protein Is Purified on Taipoxin Columns. J. Neurochem. 1995, 64, 2339–2344. [Google Scholar] [CrossRef]
  77. Six, D.A.; Dennis, E.A. The Expanding Superfamily of Phospholipase A2 Enzymes: Classification and Characterization. Biochim. Biophys. Acta BBA—Mol. Cell Biol. Lipids 2000, 1488, 1–19. [Google Scholar] [CrossRef]
  78. Podobnik, M.; Anderluh, G. Pore-Forming Toxins in Cnidaria. Semin. Cell Dev. Biol. 2017, 72, 133–141. [Google Scholar] [CrossRef]
  79. Parker, M.W.; Feil, S.C. Pore-Forming Protein Toxins: From Structure to Function. Prog. Biophys. Mol. Biol. 2005, 88, 91–142. [Google Scholar] [CrossRef]
  80. Anderluh, G.; Maček, P. Dissecting the Actinoporin Pore-Forming Mechanism. Structure 2003, 11, 1312–1313. [Google Scholar] [CrossRef] [Green Version]
  81. Anderluh, G.; Kisovec, M.; Kraševec, N.; Gilbert, R.J.C. Distribution of MACPF/CDC proteins. In MACPF/CDC Proteins—Agents of Defence, Attack and Invasion; Anderluh, G., Gilbert, R., Eds.; Subcellular Biochemistry; Springer: Dordrecht, The Netherlands, 2014; pp. 7–30. ISBN 978-94-017-8881-6. [Google Scholar]
  82. Tibballs, J.; Yanagihara, A.A.; Turner, H.C.; Winkel, K. Immunological and Toxinological Responses to Jellyfish Stings. Inflamm. Allergy—Drug Targets Former. Curr. Drug Targets—Inflamm. Allergy 2011, 10, 438–446. [Google Scholar] [CrossRef] [Green Version]
  83. Horiike, T.; Nagai, H.; Kitani, S. Identification of Allergens in the Box Jellyfish Chironex yamaguchii That Cause Sting Dermatitis. Int. Arch. Allergy Immunol. 2015, 167, 73–82. [Google Scholar] [CrossRef] [PubMed]
  84. Tibballs, J. Australian chirodropid cubozoan jellyfish envenomation. In Clinical Toxinology; Gopalakrishnakone, P., Faiz, S.M.A., Gnanathasan, C.A., Habib, A.G., Fernando, R., Yang, C.-C., Vogel, C.-W., Tambourgi, D.V., Seifert, S.A., Eds.; Toxinology; Springer: Dordrecht, The Netherlands, 2017; pp. 1–19. ISBN 978-94-007-6288-6. [Google Scholar]
  85. Wang, Y.-M.; Wang, S.-R.; Tsai, I.-H. Serine Protease Isoforms of Deinagkistrodon acutus Venom: Cloning, Sequencing and Phylogenetic Analysis. Biochem. J. 2001, 354, 161–168. [Google Scholar] [CrossRef]
  86. Weinstock, G.M.; Robinson, G.E.; Gibbs, R.A.; Weinstock, G.M.; Weinstock, G.M.; Robinson, G.E.; Worley, K.C.; Evans, J.D.; Maleszka, R.; Robertson, H.M.; et al. Insights into Social Insects from the Genome of the Honeybee Apis mellifera. Nature 2006, 443, 931–949. [Google Scholar] [CrossRef] [Green Version]
  87. Bernheimer, A.W.; Weinstein, S.A.; Linder, R. Isoelectric Analysis of Some Australian Elapid Snake Venoms with Special Reference to Phospholipase B and Hemolysis. Toxicon 1986, 24, 841–849. [Google Scholar] [CrossRef]
  88. Sabourault, C.; Ganot, P.; Deleury, E.; Allemand, D.; Furla, P. Comprehensive EST Analysis of the Symbiotic Sea Anemone, Anemonia Viridis. BMC Genom. 2009, 10, 333. [Google Scholar] [CrossRef] [Green Version]
  89. Kozlov, S.; Grishin, E. The Mining of Toxin-like Polypeptides from EST Database by Single Residue Distribution Analysis. BMC Genom. 2011, 12, 88. [Google Scholar] [CrossRef] [Green Version]
  90. Oliveira, J.S.; Fuentes-Silva, D.; King, G.F. Development of a Rational Nomenclature for Naming Peptide and Protein Toxins from Sea Anemones. Toxicon 2012, 60, 539–550. [Google Scholar] [CrossRef]
  91. Hines, D.E.; Pawlik, J.R. Assessing the Antipredatory Defensive Strategies of Caribbean Non-Scleractinian Zoantharians (Cnidaria): Is the Sting the Only Thing? Mar. Biol. 2012, 159, 389–398. [Google Scholar] [CrossRef]
  92. Ashwood, L.M.; Norton, R.S.; Undheim, E.A.B.; Hurwood, D.A.; Prentis, P.J. Characterising Functional Venom Profiles of Anthozoans and Medusozoans within Their Ecological Context. Mar. Drugs 2020, 18, 202. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  93. Jaimes-Becerra, A.; Gacesa, R.; Doonan, L.B.; Hartigan, A.; Marques, A.C.; Okamura, B.; Long, P.F. “Beyond Primary Sequence”—Proteomic Data Reveal Complex Toxins in Cnidarian Venoms. Integr. Comp. Biol. 2019, 59, 777–785. [Google Scholar] [CrossRef]
  94. Schendel, V.; Rash, L.D.; Jenner, R.A.; Undheim, E.A.B. The Diversity of Venom: The Importance of Behavior and Venom System Morphology in Understanding Its Ecology and Evolution. Toxins 2019, 11, 666. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  95. Pennington, M.W.; Czerwinski, A.; Norton, R.S. Peptide Therapeutics from Venom: Current Status and Potential. Bioorg. Med. Chem. 2018, 26, 2738–2758. [Google Scholar] [CrossRef]
  96. Tarcha, E.J.; Chi, V.; Muñoz-Elías, E.J.; Bailey, D.; Londono, L.M.; Upadhyay, S.K.; Norton, K.; Banks, A.; Tjong, I.; Nguyen, H.; et al. Durable Pharmacological Responses from the Peptide ShK-186, a Specific Kv1.3 Channel Inhibitor That Suppresses T Cell Mediators of Autoimmune Disease. J. Pharmacol. Exp. Ther. 2012, 342, 642–653. [Google Scholar] [CrossRef] [Green Version]
  97. Tarcha, E.J.; Olsen, C.M.; Probst, P.; Peckham, D.; Muñoz-Elías, E.J.; Kruger, J.G.; Iadonato, S.P. Safety and Pharmacodynamics of Dalazatide, a Kv1.3 Channel Inhibitor, in the Treatment of Plaque Psoriasis: A Randomized Phase 1b Trial. PLoS ONE 2017, 12, e0180762. [Google Scholar] [CrossRef] [PubMed]
  98. Klompen, A.M.L.; Macrander, J.; Reitzel, A.M.; Stampar, S.N. Transcriptomic Analysis of Four Cerianthid (Cnidaria, Ceriantharia) Venoms. Mar. Drugs 2020, 18, 413. [Google Scholar] [CrossRef]
  99. Li, R.; Yu, H.; Yue, Y.; Liu, S.; Xing, R.; Chen, X.; Li, P. Combined Proteomics and Transcriptomics Identifies Sting-Related Toxins of Jellyfish Cyanea nozakii. J. Proteom. 2016, 148, 57–64. [Google Scholar] [CrossRef]
  100. Brinkman, D.L.; Jia, X.; Potriquet, J.; Kumar, D.; Dash, D.; Kvaskoff, D.; Mulvenna, J. Transcriptome and Venom Proteome of the Box Jellyfish Chironex fleckeri. BMC Genom. 2015, 16, 407. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  101. Fu, J.; He, Y.; Peng, C.; Tang, T.; Jin, A.; Liao, Y.; Shi, Q.; Gao, B. Transcriptome Sequencing of the Pale Anemones (Exaiptasia diaphana) Revealed Functional Peptide Gene Resources of Sea Anemone. Front. Mar. Sci. 2022, 9, 856501. [Google Scholar] [CrossRef]
  102. Mitchell, M.L.; Hossain, M.A.; Lin, F.; Pinheiro-Junior, E.L.; Peigneur, S.; Wai, D.C.C.; Delaine, C.; Blyth, A.J.; Forbes, B.E.; Tytgat, J.; et al. Identification, Synthesis, Conformation and Activity of an Insulin-like Peptide from a Sea Anemone. Biomolecules 2021, 11, 1785. [Google Scholar] [CrossRef]
  103. Krishnarjuna, B.; Sunanda, P.; Villegas–Moreno, J.; Csoti, A.; Morales, R.A.V.; Wai, D.C.C.; Panyi, G.; Prentis, P.; Norton, R.S. A Disulfide-Stabilised Helical Hairpin Fold in Acrorhagin I: An Emerging Structural Motif in Peptide Toxins. J. Struct. Biol. 2021, 213, 107692. [Google Scholar] [CrossRef]
  104. Honma, T.; Minagawa, S.; Nagai, H.; Ishida, M.; Nagashima, Y.; Shiomi, K. Novel Peptide Toxins from Acrorhagi, Aggressive Organs of the Sea Anemone Actinia equina. Toxicon 2005, 46, 768–774. [Google Scholar] [CrossRef]
  105. Macrander, J.; Brugler, M.R.; Daly, M. A RNA-Seq Approach to Identify Putative Toxins from Acrorhagi in Aggressive and Non-Aggressive Anthopleura elegantissima Polyps. BMC Genom. 2015, 16, 221. [Google Scholar] [CrossRef]
  106. Minagawa, S.; Sugiyama, M.; Ishida, M.; Nagashima, Y.; Shiomi, K. Kunitz-Type Protease Inhibitors from Acrorhagi of Three Species of Sea Anemones. Comp. Biochem. Physiol. B Biochem. Mol. Biol. 2008, 150, 240–245. [Google Scholar] [CrossRef] [PubMed]
  107. Columbus-Shenkar, Y.Y.; Sachkova, M.Y.; Macrander, J.; Fridrich, A.; Modepalli, V.; Reitzel, A.M.; Sunagar, K.; Moran, Y. Dynamics of Venom Composition across a Complex Life Cycle. eLife 2018, 7, e35014. [Google Scholar] [CrossRef] [PubMed]
  108. Leach, W.B.; Reitzel, A.M. Transcriptional Remodelling upon Light Removal in a Model Cnidarian: Losses and Gains in Gene Expression. Mol. Ecol. 2019, 28, 3413–3426. [Google Scholar] [CrossRef]
  109. Sachkova, M.Y.; Singer, S.A.; Macrander, J.; Reitzel, A.M.; Peigneur, S.; Tytgat, J.; Moran, Y. The Birth and Death of Toxins with Distinct Functions: A Case Study in the Sea Anemone Nematostella. Mol. Biol. Evol. 2019, 36, 2001–2012. [Google Scholar] [CrossRef] [Green Version]
  110. Smith, E.G.; Surm, J.M.; Macrander, J.; Simhi, A.; Amir, G.; Sachkova, M.Y.; Lewandowska, M.; Reitzel, A.M.; Moran, Y. Dominant Toxin Hypothesis: Unravelling the Venom Phenotype across Micro and Macroevolution. bioRxiv 2022. [Google Scholar] [CrossRef]
  111. Verdes, A.; Taboada, S.; Hamilton, B.R.; Undheim, E.A.B.; Sonoda, G.G.; Andrade, S.C.S.; Morato, E.; Marina, A.I.; Cárdenas, C.A.; Riesgo, A. Evolution, Expression Patterns, and Distribution of Novel Ribbon Worm Predatory and Defensive Toxins. Mol. Biol. Evol. 2022, 39, msac096. [Google Scholar] [CrossRef]
  112. Ramírez-Carreto, S.; Vera-Estrella, R.; Portillo-Bobadilla, T.; Licea-Navarro, A.; Bernaldez-Sarabia, J.; Rudiño-Piñera, E.; Verleyen, J.J.; Rodríguez, E.; Rodríguez-Almazán, C. Transcriptomic and Proteomic Analysis of the Tentacles and Mucus of Anthopleura dowii Verrill, 1869. Mar. Drugs 2019, 17, 436. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  113. Holding, M.L.; Margres, M.J.; Mason, A.J.; Parkinson, C.L.; Rokyta, D.R. Evaluating the Performance of De Novo Assembly Methods for Venom-Gland Transcriptomics. Toxins 2018, 10, 249. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  114. Grabherr, M.G.; Haas, B.J.; Yassour, M.; Levin, J.Z.; Thompson, D.A.; Amit, I.; Adiconis, X.; Fan, L.; Raychowdhury, R.; Zeng, Q.; et al. Full-Length Transcriptome Assembly from RNA-Seq Data without a Reference Genome. Nat. Biotechnol. 2011, 29, 644–652. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  115. Simão, F.A.; Waterhouse, R.M.; Ioannidis, P.; Kriventseva, E.V.; Zdobnov, E.M. BUSCO: Assessing Genome Assembly and Annotation Completeness with Single-Copy Orthologs. Bioinformatics 2015, 31, 3210–3212. [Google Scholar] [CrossRef]
  116. Boratyn, G.M.; Camacho, C.; Cooper, P.S.; Coulouris, G.; Fong, A.; Ma, N.; Madden, T.L.; Matten, W.T.; McGinnis, S.D.; Merezhuk, Y.; et al. BLAST: A More Efficient Report with Usability Improvements. Nucleic Acids Res. 2013, 41, W29–W33. [Google Scholar] [CrossRef] [Green Version]
  117. Finn, R.D.; Clements, J.; Arndt, W.; Miller, B.L.; Wheeler, T.J.; Schreiber, F.; Bateman, A.; Eddy, S.R. HMMER Web Server: 2015 Update. Nucleic Acids Res. 2015, 43, W30–W38. [Google Scholar] [CrossRef]
  118. Armenteros, J.J.A.; Salvatore, M.; Emanuelsson, O.; Winther, O.; von Heijne, G.; Elofsson, A.; Nielsen, H. Detecting Sequence Signals in Targeting Peptides Using Deep Learning. Life Sci. Alliance 2019, 2, e201900429. [Google Scholar] [CrossRef] [Green Version]
  119. Fu, L.; Niu, B.; Zhu, Z.; Wu, S.; Li, W. CD-HIT: Accelerated for Clustering the next-Generation Sequencing Data. Bioinformatics 2012, 28, 3150–3152. [Google Scholar] [CrossRef] [PubMed]
  120. Edgar, R.C. MUSCLE: Multiple Sequence Alignment with High Accuracy and High Throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef] [Green Version]
  121. Madeira, F.; Pearce, M.; Tivey, A.R.N.; Basutkar, P.; Lee, J.; Edbali, O.; Madhusoodanan, N.; Kolesnikov, A.; Lopez, R. Search and Sequence Analysis Tools Services from EMBL-EBI in 2022. Nucleic Acids Res. 2022, 50, W276–W279. [Google Scholar] [CrossRef] [PubMed]
  122. The UniProt Consortium UniProt: A Worldwide Hub of Protein Knowledge. Nucleic Acids Res. 2019, 47, D506–D515. [CrossRef]
Marinedrugs 20 00730 g001 550
Figure 1. (a) Stacked bar graph showing the proportion of toxin functional catagories across each of the focal taxa. (b) table showing the percentage of toxin categories per species.
Figure 1. (a) Stacked bar graph showing the proportion of toxin functional catagories across each of the focal taxa. (b) table showing the percentage of toxin categories per species.
Marinedrugs 20 00730 g001
Table
Table 1. Assembly metrics for transcriptome assemblies. The values given in the original assemblies are italicized and in parentheses, dash (-) indicates previously missing metric.
Table 1. Assembly metrics for transcriptome assemblies. The values given in the original assemblies are italicized and in parentheses, dash (-) indicates previously missing metric.
CladeTaxon [Source]TranscriptsGenesN50BUSCO
StichodactylinaC. adhaesivum [45]905,882 (628,469)484,240 (451,132)1308
(609)		96.2
(93.9)
S. haddoni [46]196,433 (1,800,000)150,054 (269,628)1038
(-)		83.6
(-)
Ha. hemprichii [45]600,056 (101,150)165,714 (74,496)1055
(1370)		94.5
(92)
HeteractinaHs. crispa [47]655,116
(-)		581,957
(-)		545
(-)		91.8
(-)
M. doreensis [45]624,291
(-)		218,128
(-)		1328
(-)		91.4
(-)
EntacmaeaE. quadricolor
[PRJEB21970]		1,198,081
(-)		296,968
(-)		1088
(-)		98.3
(-)
Table
Table 2. Number of toxin hits across focal taxa and across families (in bold). (T = number of putative toxins; TF = Number of Toxin Families; AT = Number of Actiniarian Toxins; ATF = Number of Actiniarian Toxin Families).
Table 2. Number of toxin hits across focal taxa and across families (in bold). (T = number of putative toxins; TF = Number of Toxin Families; AT = Number of Actiniarian Toxins; ATF = Number of Actiniarian Toxin Families).
TTFATATF
C. adhaesivum256348813
S. haddoni112334612
Ha. hemprichii115322713
Stychodactylina4834716118
M. doreensis135255612
Hs. crispa175355912
Heteractina3103711513
Entacmaea quadricolor
(Entacmaea)		328377913
Total11215135518
Table
Table 3. Putative toxin diversity among focal taxa sorted by various venom categories.
Table 3. Putative toxin diversity among focal taxa sorted by various venom categories.
Venom CategorySpeciesVenom IDCASHHHMDHCEQ
Auxiliarysea anemoneAstacin-like metalloprotease toxin M12A1552151322
Nematocyst expressed protein 6248491321
snakeCystatin-2-2-111
scorpionVenom protein 3021-----
Venom protein 59.1-1----
Pore Formingsea anemoneActinoporin family1-1152
MAC/PF1-1711
Neurotoxinsea anemoneCnidaria small cysteine-rich protein (SCRiP) family--1--1
Cysteine-rich venom protein8--25-
Neurotoxin ShK-like1--1---
Sea anemone short toxin (type III) family (Delta-actitoxin-Avd2b1)-11---
Sea anemone structural class 9a11---3
Sea anemone type 1 potassium channel toxin subfamily43-111
Sea anemone type 3 potassium channel toxin subfamily (BDS)786154
sodium channel toxin Type I-3-22-
sodium channel toxin Type II326---
sea snailConopeptide P-like superfamily1023258
spiderOmega-theraphotoxin-Pm1a, Neurotoxin 10-1----
Psalmotoxin-1 (PcTx1)-1----
centipedeScoloptoxin-1--2-
Mixed Enzymesea anemonePhospholipase A2-2--52
Protease Inhibitorsea anemoneVenom Kunitz-type family733927
Allergen and Innate ImmunitysnakeA. superbus venom factor 12-2-13
snake venom serine protease Dav-PA3-----
Snake venom serine protease salmobin---21-
Venom serine carboxypeptidase82-675
spiderTechylectin-like protein-----1
antVenom allergen 5231365
mothVenom serine protease71321-
beeVenom serine protease Bi-VSP713213
Hemostatic and Hemorrhagic ToxinshrewBlarina toxin1314-39
snakeCoagulation factor V--1-31
Coagulation factor X-activating enzyme41--34
Cobra venom factor-12--1
Ryncolin148133335
Snaclec family77110136
True venom lectin family43-5326
Veficolin-11791491045
Venom factor (VF)3-4--1
Venom phosphodiesterase-- -1-
Venom prothrombin activator2551583248
Zinc metalloproteinase/disintegrin-like (M12B)726639
fishGalactose-specific lectin nattectin86-614
Natterin-4----41
lizardVenom protease64511412
spiderVenom peptide isomerase heavy chain17563715
unknownsea anemoneEGF domain peptide family631215
Sea anemone 8 toxin family1--1-1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Delgado, A.; Benedict, C.; Macrander, J.; Daly, M. Never, Ever Make an Enemy… Out of an Anemone: Transcriptomic Comparison of Clownfish Hosting Sea Anemone Venoms. Mar. Drugs 2022, 20, 730. https://doi.org/10.3390/md20120730

AMA Style

Delgado A, Benedict C, Macrander J, Daly M. Never, Ever Make an Enemy… Out of an Anemone: Transcriptomic Comparison of Clownfish Hosting Sea Anemone Venoms. Marine Drugs. 2022; 20(12):730. https://doi.org/10.3390/md20120730

Chicago/Turabian Style

Delgado, Alonso, Charlotte Benedict, Jason Macrander, and Marymegan Daly. 2022. "Never, Ever Make an Enemy… Out of an Anemone: Transcriptomic Comparison of Clownfish Hosting Sea Anemone Venoms" Marine Drugs 20, no. 12: 730. https://doi.org/10.3390/md20120730

APA Style

Delgado, A., Benedict, C., Macrander, J., & Daly, M. (2022). Never, Ever Make an Enemy… Out of an Anemone: Transcriptomic Comparison of Clownfish Hosting Sea Anemone Venoms. Marine Drugs, 20(12), 730. https://doi.org/10.3390/md20120730

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop