A new transcriptome and transcriptome profiling of adult and larval tissue in the box jellyfish Alatina alata: an emerging model for studying venom, vision and sex

Venom biotechnology: casting light on nature's deadliest weapons using synthetic biology

Venoms are complex chemical arsenals that have evolved independently many times in the animal kingdom. Venoms have attracted the interest of researchers because they are an important innovation that has contributed greatly to the evolutionary success of many animals, and their medical relevance offers significant potential for drug discovery. During the last decade, venom research has been revolutionized by the application of systems biology, giving rise to a novel field known as venomics. More recently, biotechnology has also made an increasing impact in this field. Its methods provide the means to disentangle and study venom systems across all levels of biological organization and, given their tremendous impact on the life sciences, these pivotal tools greatly facilitate the coherent understanding of venom system organization, development, biochemistry, and therapeutic activity. Even so, we lack a comprehensive overview of major advances achieved by applying biotechnology to venom systems. This review therefore considers the methods, insights, and potential future developments of biotechnological applications in the field of venom research. We follow the levels of biological organization and structure, starting with the methods used to study the genomic blueprint and genetic machinery of venoms, followed gene products and their functional phenotypes. We argue that biotechnology can answer some of the most urgent questions in venom research, particularly when multiple approaches are combined together, and with other venomics technologies.

The state of Medusozoa genomics: current evidence and future challenges

Medusozoa is a widely distributed ancient lineage that harbors one-third of Cnidaria diversity divided into 4 classes. This clade is characterized by the succession of stages and modes of reproduction during metagenic lifecycles, and includes some of the most plastic body plans and life cycles among animals. The characterization of traditional genomic features, such as chromosome numbers and genome sizes, was rather overlooked in Medusozoa and many evolutionary questions still remain unanswered. Modern genomic DNA sequencing in this group started in 2010 with the publication of the Hydra vulgaris genome and has experienced an exponential increase in the past 3 years. Therefore, an update of the state of Medusozoa genomics is warranted. We reviewed different sources of evidence, including cytogenetic records and high-throughput sequencing projects. We focused on 4 main topics that would be relevant for the broad Cnidaria research community: (i) taxonomic coverage of genomic information; (ii) continuity, quality, and completeness of high-throughput sequencing datasets; (iii) overview of the Medusozoa specific research questions approached with genomics; and (iv) the accessibility of data and metadata. We highlight a lack of standardization in genomic projects and their reports, and reinforce a series of recommendations to enhance future collaborative research.

Evolution of gene expression across species and specialized zooids in Siphonophora

Siphonophores are complex colonial animals, consisting of asexually produced bodies (zooids) that are functionally specialized for specific tasks, including feeding, swimming, and sexual reproduction. Though this extreme functional specialization has captivated biologists for generations, its genomic underpinnings remain unknown. We use RNA-seq to investigate gene expression patterns in five zooids and one specialized tissue across seven siphonophore species. Analyses of gene expression across species present several challenges, including identification of comparable expression changes on gene trees with complex histories of speciation, duplication, and loss. We examine gene expression within species, conduct classical analyses examining expression patterns between species, and introduce species branch filtering, which allows us to examine the evolution of expression across species in a phylogenetic framework. Within and across species, we identified hundreds of zooid-specific and species-specific genes, as well as a number of putative transcription factors showing differential expression in particular zooids and developmental stages. We found that gene expression patterns tended to be largely consistent in zooids with the same function across species, but also some large lineage-specific shifts in gene expression. Our findings show that patterns of gene expression have the potential to define zooids in colonial organisms. Traditional analyses of the evolution of gene expression focus on the tips of gene phylogenies, identifying large-scale expression patterns that are zooid or species variable. The new explicit phylogenetic approach we propose here focuses on branches (not tips) offering a deeper evolutionary perspective into specific changes in gene expression within zooids along all branches of the gene (and species) trees.

The cnidarian parasite Ceratonova shasta utilizes inherited and recruited venom-like compounds during infection

Background

Cnidarians are the most ancient venomous organisms. They store a cocktail of venom proteins inside unique stinging organelles called nematocysts. When a cnidarian encounters chemical and physical cues from a potential threat or prey animal, the nematocyst is triggered and fires a harpoon-like tubule to penetrate and inject venom into the prey. Nematocysts are present in all Cnidaria, including the morphologically simple Myxozoa, which are a speciose group of microscopic, spore-forming, obligate parasites of fish and invertebrates. Rather than predation or defense, myxozoans use nematocysts for adhesion to hosts, but the involvement of venom in this process is poorly understood. Recent work shows some myxozoans have a reduced repertoire of venom-like compounds (VLCs) relative to free-living cnidarians, however the function of these proteins is not known.

Methods

We searched for VLCs in the nematocyst proteome and a time-series infection transcriptome of Ceratonova shasta, a myxozoan parasite of salmonid fish. We used four parallel approaches to detect VLCs: BLAST and HMMER searches to preexisting cnidarian venom datasets, the machine learning tool ToxClassifier, and structural modeling of nematocyst proteomes. Sequences that scored positive by at least three methods were considered VLCs. We then mapped their time-series expressions in the fish host and analyzed their phylogenetic relatedness to sequences from other venomous animals.

Results

We identified eight VLCs, all of which have closely related sequences in other myxozoan datasets, suggesting a conserved venom profile across Myxozoa, and an overall reduction in venom diversity relative to free-living cnidarians. Expression of the VLCs over the 3-week fish infection varied considerably: three sequences were most expressed at one day post-exposure in the fish’s gills; whereas expression of the other five VLCs peaked at 21 days post-exposure in the intestines, coinciding with the formation of mature parasite spores with nematocysts. Expression of VLC genes early in infection, prior to the development of nematocysts, suggests venoms in C. shasta have been repurposed to facilitate parasite invasion and proliferation within the host. Molecular phylogenetics suggested some VLCs were inherited from a cnidarian ancestor, whereas others were more closely related to sequences from venomous non-Cnidarian organisms and thus may have gained qualities of venom components via convergent evolution. The presence of VLCs and their differential expression during parasite infection enrich the concept of what functions a “venom” can have and represent targets for designing therapeutics against myxozoan infections.

Conservative and Atypical Ferritins of Sponges

Ferritins comprise a conservative family of proteins found in all species and play an essential role in resistance to redox stress, immune response, and cell differentiation. Sponges (Porifera) are the oldest Metazoa that show unique plasticity and regenerative potential. Here, we characterize the ferritins of two cold-water sponges using proteomics, spectral microscopy, and bioinformatic analysis. The recently duplicated conservative HdF1a/b and atypical HdF2 genes were found in the Halisarca dujardini genome. Multiple related transcripts of HpF1 were identified in the Halichondria panicea transcriptome. Expression of HdF1a/b was much higher than that of HdF2 in all annual seasons and regulated differently during the sponge dissociation/reaggregation. The presence of the MRE and HRE motifs in the HdF1 and HdF2 promotor regions and the IRE motif in mRNAs of HdF1 and HpF indicates that sponge ferritins expression depends on the cellular iron and oxygen levels. The gel electrophoresis combined with specific staining and mass spectrometry confirmed the presence of ferric ions and ferritins in multi-subunit complexes. The 3D modeling predicts the iron-binding capacity of HdF1 and HpF1 at the ferroxidase center and the absence of iron-binding in atypical HdF2. Interestingly, atypical ferritins lacking iron-binding capacity were found in genomes of many invertebrate species. Their function deserves further research.

SuMO-Fil: Supervised multi-omic filtering prior to performing network analysis

Multi-omic analyses that integrate many high-dimensional datasets often present significant deficiencies in statistical power and require time consuming computations to execute the analytical methods. We present SuMO-Fil to remedy against these issues which is a pre-processing method for Supervised Multi-Omic Filtering that removes variables or features considered to be irrelevant noise. SuMO-Fil is intended to be performed prior to downstream analyses that detect supervised gene networks in sparse settings. We accomplish this by implementing variable filters based on low similarity across the datasets in conjunction with low similarity with the outcome. This approach can improve accuracy, as well as reduce run times for a variety of computationally expensive downstream analyses. This method has applications in a setting where the downstream analysis may include sparse canonical correlation analysis. Filtering methods specifically for cluster and network analysis are introduced and compared by simulating modular networks with known statistical properties. The SuMO-Fil method performs favorably by eliminating non-network features while maintaining important biological signal under a variety of different signal settings as compared to popular filtering techniques based on low means or low variances. We show that the speed and accuracy of methods such as supervised sparse canonical correlation are increased after using SuMO-Fil, thus greatly improving the scalability of these approaches.

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.

A Review of Toxins from Cnidaria

Cnidarians have been known since ancient times for the painful stings they induce to humans. The effects of the stings range from skin irritation to cardiotoxicity and can result in death of human beings. The noxious effects of cnidarian venoms have stimulated the definition of their composition and their activity. Despite this interest, only a limited number of compounds extracted from cnidarian venoms have been identified and defined in detail. Venoms extracted from Anthozoa are likely the most studied, while venoms from Cubozoa attract research interests due to their lethal effects on humans. The investigation of cnidarian venoms has benefited in very recent times by the application of omics approaches. In this review, we propose an updated synopsis of the toxins identified in the venoms of the main classes of Cnidaria (Hydrozoa, Scyphozoa, Cubozoa, Staurozoa and Anthozoa). We have attempted to consider most of the available information, including a summary of the most recent results from omics and biotechnological studies, with the aim to define the state of the art in the field and provide a background for future research.

A comparative genomics study of neuropeptide genes in the cnidarian subclasses Hexacorallia and Ceriantharia

Background

Nervous systems originated before the split of Proto- and Deuterostomia, more than 600 million years ago. Four animal phyla (Cnidaria, Placozoa, Ctenophora, Porifera) diverged before this split and studying these phyla could give us important information on the evolution of the nervous system. Here, we have annotated the neuropeptide preprohormone genes of twenty species belonging to the subclass Hexacorallia or Ceriantharia (Anthozoa: Cnidaria), using thirty-seven publicly accessible genome or transcriptome databases. Studying hexacorals is important, because they are versatile laboratory models for development (e.g., Nematostella vectensis) and symbiosis (e.g., Exaiptasia diaphana) and also are prominent reef-builders.

Results

We found that each hexacoral or ceriantharian species contains five to ten neuropeptide preprohormone genes. Many of these preprohormones contain multiple copies of immature neuropeptides, which can be up to 50 copies of identical or similar neuropeptide sequences. We also discovered preprohormones that only contained one neuropeptide sequence positioned directly after the signal sequence. Examples of them are neuropeptides that terminate with the sequence RWamide (the Antho-RWamides). Most neuropeptide sequences are N-terminally protected by pyroglutamyl (pQ) or one or more prolyl residues, while they are C-terminally protected by an amide group. Previously, we isolated and sequenced small neuropeptides from hexacorals that were N-terminally protected by an unusual L-3-phenyllactyl group. In our current analysis, we found that these N-phenyllactyl-peptides are derived from N-phenylalanyl-peptides located directly after the signal sequence of the preprohormone. The N-phenyllactyl- peptides appear to be confined to the hexacorallian order Actiniaria and do not occur in other cnidarians. On the other hand, (1) the neuropeptide Antho-RFamide (pQGRFamide); (2) peptides with the C-terminal sequence GLWamide; and (3) tetrapeptides with the X₁PRX₂amide consensus sequence (most frequently GPRGamide) are ubiquitous in Hexacorallia.

Conclusions

We found GRFamide, GLWamide, and X₁PRX₂amide peptides in all tested Hexacorallia. Previously, we discovered these three neuropeptide classes also in Cubozoa, Scyphozoa, and Staurozoa, indicating that these neuropeptides originated in the common cnidarian ancestor and are evolutionarily ancient. In addition to these ubiquitous neuropeptides, other neuropeptides appear to be confined to specific cnidarian orders or subclasses.

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.

Cassiosomes are stinging-cell structures in the mucus of the upside-down jellyfish Cassiopea xamachana

Snorkelers in mangrove forest waters inhabited by the upside-down jellyfish Cassiopea xamachana report discomfort due to a sensation known as stinging water, the cause of which is unknown. Using a combination of histology, microscopy, microfluidics, videography, molecular biology, and mass spectrometry-based proteomics, we describe C. xamachana stinging-cell structures that we term cassiosomes. These structures are released within C. xamachana mucus and are capable of killing prey. Cassiosomes consist of an outer epithelial layer mainly composed of nematocytes surrounding a core filled by endosymbiotic dinoflagellates hosted within amoebocytes and presumptive mesoglea. Furthermore, we report cassiosome structures in four additional jellyfish species in the same taxonomic group as C. xamachana (Class Scyphozoa; Order Rhizostomeae), categorized as either motile (ciliated) or nonmotile types. This inaugural study provides a qualitative assessment of the stinging contents of C. xamachana mucus and implicates mucus containing cassiosomes and free intact nematocytes as the cause of stinging water.

Cheryl L Ames, Anna Klompen et al. describe cassiosomes, stinging cell structures in the mucus of the upside-down jellyfish Cassiopea xamachana. They show that these motile cell masses consist of an outer epithelial layer largely composed of nematocytes surrounding centralized clusters of endosymbiotic dinoflagellates.

Transcriptome Characterization of Reverse Development in Turritopsis dohrnii (Hydrozoa, Cnidaria)

Medusae of Turritopsis dohrnii undergo reverse development in response to physical damage, adverse environmental conditions, or aging. Senescent, weakened or damaged medusae transform into a cluster of poorly differentiated cells (known as the cyst stage), which metamorphose back into a preceding life cycle stage, the polyp. During the metamorphosis, cell transdifferentiation occurs. The cyst represents the intermediate stage between a reverting medusa and a healthy polyp, during which cell transdifferentiation and tissue reorganization take place. Here we characterize and compare the transcriptomes of the polyp and newborn medusa stages of T. dohrnii with that of the cyst, to identify biological networks potentially involved in the reverse development and transdifferentiation processes. The polyp, medusa and cyst of T. dohrnii were sequenced through Illumina RNA-sequencing and assembled using a de novo approach, resulting in 92,569, 74,639 and 86,373 contigs, respectively. The transcriptomes were annotated and comparative analyses among the stages identified biological networks that were significantly over-and under-expressed in the cyst as compared to the polyp and medusa stages. Biological processes that occur at the cyst stage such as telomerase activity, regulation of transposable elements and DNA repair systems, and suppression of cell signaling pathways, mitotic cell division and cellular differentiation and development may be involved in T. dohrnii's reverse development and transdifferentiation. Our results are the first attempt to understand T. dohrnii's life-cycle reversal at the genetic level, and indicate possible avenues of future research on developmental strategies, cell transdifferentiation, and aging using T. dohrnii as a non-traditional in vivo system.

A non-bilaterian perspective on the development and evolution of animal digestive systems

Digestive systems and extracellular digestion are key animal features, but their emergence during early animal evolution is currently poorly understood. As the last common ancestor of non-bilaterian animal groups (sponges, ctenophores, placozoans and cnidarians) dates back to the beginning of animal life, their study and comparison provides important insights into the early evolution of digestive systems and functions. Here, I have compiled an overview of the development and cell biology of digestive tissues in non-bilaterian animals. I will highlight the fundamental differences between extracellular and intracellular digestive processes, and how these are distributed among animals. Cnidarians (e.g. sea anemones, corals, jellyfish), the phylogenetic outgroup of bilaterians (e.g. vertebrates, flies, annelids), occupy a key position to reconstruct the evolution of bilaterian gut evolution. A major focus will therefore lie on the development and cell biology of digestive tissues in cnidarians, especially sea anemones, and how they compare to bilaterian gut tissues. In that context, I will also review how a recent study on the gastrula fate map of the sea anemone Nematostella vectensis challenges our long-standing conceptions on the evolution of cnidarian and bilaterian germ layers and guts.

Box, stalked, and upside-down? Draft genomes from diverse jellyfish (Cnidaria, Acraspeda) lineages: Alatina alata (Cubozoa), Calvadosia cruxmelitensis (Staurozoa), and Cassiopea xamachana (Scyphozoa)

BACKGROUND

Anthozoa, Endocnidozoa, and Medusozoa are the 3 major clades of Cnidaria. Medusozoa is further divided into 4 clades, Hydrozoa, Staurozoa, Cubozoa, and Scyphozoa-the latter 3 lineages make up the clade Acraspeda. Acraspeda encompasses extraordinary diversity in terms of life history, numerous nuisance species, taxa with complex eyes rivaling other animals, and some of the most venomous organisms on the planet. Genomes have recently become available within Scyphozoa and Cubozoa, but there are currently no published genomes within Staurozoa and Cubozoa.

FINDINGS

Here we present 3 new draft genomes of Calvadosia cruxmelitensis (Staurozoa), Alatina alata (Cubozoa), and Cassiopea xamachana (Scyphozoa) for which we provide a preliminary orthology analysis that includes an inventory of their respective venom-related genes. Additionally, we identify synteny between POU and Hox genes that had previously been reported in a hydrozoan, suggesting this linkage is highly conserved, possibly dating back to at least the last common ancestor of Medusozoa, yet likely independent of vertebrate POU-Hox linkages.

CONCLUSIONS

These draft genomes provide a valuable resource for studying the evolutionary history and biology of these extraordinary animals, and for identifying genomic features underlying venom, vision, and life history traits in Acraspeda.

Global Neuropeptide Annotations From the Genomes and Transcriptomes of Cubozoa, Scyphozoa, Staurozoa (Cnidaria: Medusozoa), and Octocorallia (Cnidaria: Anthozoa)

During animal evolution, ancestral Cnidaria and Bilateria diverged more than 600 million years ago. The nervous systems of extant cnidarians are strongly peptidergic. Neuropeptides have been isolated and sequenced from a few model cnidarians, but a global investigation of the presence of neuropeptides in all cnidarian classes has been lacking. Here, we have used a recently developed software program to annotate neuropeptides in the publicly available genomes and transcriptomes from members of the classes Cubozoa, Scyphozoa, and Staurozoa (which all belong to the subphylum Medusozoa) and contrasted these results with neuropeptides present in the subclass Octocorallia (belonging to the class Anthozoa). We found three to six neuropeptide preprohormone genes in members of the above-mentioned cnidarian classes or subclasses, each coding for several (up to thirty-two) similar or identical neuropeptide copies. Two of these neuropeptide preprohormone genes are present in all cnidarian classes/subclasses investigated, so they are good candidates for being among the first neuropeptide genes evolved in cnidarians. One of these primordial neuropeptide genes codes for neuropeptides having the C-terminal sequence GRFamide (pQGRFamide in Octocorallia; pQWLRGRFamide in Cubozoa and Scyphozoa; pQFLRGRFamide in Staurozoa). The other primordial neuropeptide gene codes for peptides having RPRSamide or closely resembling amino acid sequences. In addition to these two primordial neuropeptide sequences, cnidarians have their own class- or subclass-specific neuropeptides, which probably evolved to serve class/subclass-specific needs. When we carried out phylogenetic tree analyses of the GRFamide or RPRSamide preprohormones from cubozoans, scyphozoans, staurozoans, and octocorallia, we found that their phylogenetic relationships perfectly agreed with current models of the phylogeny of the studied cnidarian classes and subclasses. These results support the early origins of the GRFamide and RPRSamide preprohormone genes.

Cnidarian peptide neurotoxins: a new source of various ion channel modulators or blockers against central nervous systems disease

Cnidaria provide the largest source of bioactive peptides for new drug development. The venoms contain enzymes, potent pore-forming toxins and neurotoxins. The neurotoxins can immobilize predators rapidly when discharged via modifying sodium-channel-gating or blocking the potassium channel during the repolarization stage. Most cnidarian neurotoxins remain conserved under the strong influence of negative selection. Neuroactive peptides targeting the central nervous system through affinity with ion channels could provide insight leading to drug treatment of neurological diseases, which arise from ion channel dysfunctions. Although marine resources offer thousands of possible peptides, only one peptide derived from Cnidaria: ShK-186, also named dalazatide, has reached the pharmaceutical market. This review focuses on neuroprotective agents derived from cnidarian neurotoxic peptides.

Venomix: a simple bioinformatic pipeline for identifying and characterizing toxin gene candidates from transcriptomic data

The advent of next-generation sequencing has resulted in transcriptome-based approaches to investigate functionally significant biological components in a variety of non-model organism. This has resulted in the area of “venomics”: a rapidly growing field using combined transcriptomic and proteomic datasets to characterize toxin diversity in a variety of venomous taxa. Ultimately, the transcriptomic portion of these analyses follows very similar pathways after transcriptome assembly often including candidate toxin identification using BLAST, expression level screening, protein sequence alignment, gene tree reconstruction, and characterization of potential toxin function. Here we describe the Python package Venomix, which streamlines these processes using common bioinformatic tools along with ToxProt, a publicly available annotated database comprised of characterized venom proteins. In this study, we use the Venomix pipeline to characterize candidate venom diversity in four phylogenetically distinct organisms, a cone snail (Conidae; Conus sponsalis), a snake (Viperidae; Echis coloratus), an ant (Formicidae; Tetramorium bicarinatum), and a scorpion (Scorpionidae; Urodacus yaschenkoi). Data on these organisms were sampled from public databases, with each original analysis using different approaches for transcriptome assembly, toxin identification, or gene expression quantification. Venomix recovered numerically more candidate toxin transcripts for three of the four transcriptomes than the original analyses and identified new toxin candidates. In summary, we show that the Venomix package is a useful tool to identify and characterize the diversity of toxin-like transcripts derived from transcriptomic datasets. Venomix is available at: https://bitbucket.org/JasonMacrander/Venomix/.

Evidence for an Alternative Mechanism of Toxin Production in the Box Jellyfish Alatina alata

Cubozoans (box jellyfish) have a reputation as the most venomous animals on the planet. Herein, we provide a review of cubozoan prey capture and digestion informed by the scientific literature. Like all cnidarians, box jellyfish envenomation originates from structures secreted within nematocyte post-Golgi vesicles called nematocysts. When tentacles come in contact with prey or would-be predators, a cocktail of toxins is rapidly deployed from nematocysts via a long spiny tubule that serves to immobilize the target organism. The implication has long been that toxin peptides and proteins making up the venom within the nematocyst capsule are secreted directly by nematocytes during nematogenesis. However, our combined molecular and morphological analysis of the venomous box jellyfish Alatina alata suggests that gland cells with possible dual roles in secreting toxins and toxic-like enzymes are found in the gastric cirri. These putative gland cell assemblages might be functionally important internally (digestion of prey) as well as externally (envenomation) in cubozoans. Despite the absence of nematocysts in the gastric cirri of mature A. alata medusae, this area of the digestive system appears to be the region of the body where venom-implicated gene products are found in highest abundance, challenging the idea that in cnidarians venom is synthesized exclusively in, or nearby, nematocysts. In an effort to uncover evidence for a central area enriched in gland cells associated with the gastric cirri we provide a comparative description of the morphology of the digestive structures of A. alata and Carybdea box jellyfish species. Finally, we conduct a multi-faceted analysis of the gene ontology terms associated with venom-implicated genes expressed in the tentacle/pedalium and gastric cirri, with a particular emphasis on zinc metalloprotease homologs and genes encoding other bioactive proteins that are abundant in the A. alata transcriptome.

A Dipteran's Novel Sucker Punch: Evolution of Arthropod Atypical Venom with a Neurotoxic Component in Robber Flies (Asilidae, Diptera)

Predatory robber flies (Diptera, Asilidae) have been suspected to be venomous due to their ability to overpower well-defended prey. However, details of their venom composition and toxin arsenal remained unknown. Here, we provide a detailed characterization of the venom system of robber flies through the application of comparative transcriptomics, proteomics and functional morphology. Our results reveal asilid venoms to be dominated by peptides and non-enzymatic proteins, and that the majority of components in the crude venom is represented by just ten toxin families, which we have named Asilidin1-10. Contrary to what might be expected for a liquid-feeding predator, the venoms of robber flies appear to be rich in novel peptides, rather than enzymes with a putative pre-digestive role. The novelty of these peptides suggests that the robber fly venom system evolved independently from hematophagous dipterans and other pancrustaceans. Indeed, six Asilidins match no other venom proteins, while three represent known examples of peptide scaffolds convergently recruited to a toxic function. Of these, members of Asilidin1 closely resemble cysteine inhibitor knot peptides (ICK), of which neurotoxic variants occur in cone snails, assassin bugs, scorpions and spiders. Synthesis of one of these putative ICKs, U-Asilidin₁-Mar1a, followed by toxicity assays against an ecologically relevant prey model revealed that one of these likely plays a role as a neurotoxin involved in the immobilization of prey. Our results are fundamental to address these insights further and to understand processes that drive venom evolution in dipterans as well as other arthropods.

Medusa: A Review of an Ancient Cnidarian Body Form

Medusae (aka jellyfish) have multiphasic life cycles and a propensity to adapt to, and proliferate in, a plethora of aquatic habitats, connecting them to a number of ecological and societal issues. Now, in the midst of the genomics era, affordable next-generation sequencing (NGS) platforms coupled with publically available bioinformatics tools present the much-anticipated opportunity to explore medusa taxa as potential model systems. Genome-wide studies of medusae would provide a remarkable opportunity to address long-standing questions related to the biology, physiology, and nervous system of some of the earliest pelagic animals. Furthermore, medusae have become key targets in the exploration of marine natural products, in the development of marine biomarkers, and for their application to the biomedical and robotics fields. Presented here is a synopsis of the current state of medusa research, highlighting insights provided by multi-omics studies, as well as existing knowledge gaps, calling upon the scientific community to adopt a number of medusa taxa as model systems in forthcoming research endeavors.

Comparative proteomics reveals recruitment patterns of some protein families in the venoms of Cnidaria

Cnidarians are probably the oldest group of animals to be venomous, yet our current picture of cnidarian venom evolution is highly imbalanced due to limited taxon sampling. High-throughput tandem mass spectrometry was used to determine venom composition of the scyphozoan Chrysaora lactea and two cubozoans Tamoya haplonema and Chiropsalmus quadrumanus. Protein recruitment patterns were then compared against 5 other cnidarian venom proteomes taken from the literature. A total of 28 putative toxin protein families were identified, many for the first time in Cnidaria. Character mapping analysis revealed that 17 toxin protein families with predominantly cytolytic biological activities were likely recruited into the cnidarian venom proteome before the lineage split between Anthozoa and Medusozoa. Thereafter, venoms of Medusozoa and Anthozoa differed during subsequent divergence of cnidarian classes. Recruitment and loss of toxin protein families did not correlate with accepted phylogenetic patterns of Cnidaria. Selective pressures that drive toxin diversification independent of taxonomic positioning have yet to be identified in Cnidaria and now warrant experimental consideration.