The use of venom in intraspecific aggression is uncommon and venom-transmitting structures specifically used for intraspecific competition are found in few lineages of venomous taxa. Next-generation transcriptome sequencing allows robust characterization of venom diversity and exploration of functionally unique tissues. Using a tissue-specific RNA-seq approach, we investigate the venom composition and gene ontology diversity of acrorhagi, specialized structures used in intraspecific competition, in aggressive and non-aggressive polyps of the aggregating sea anemone Anthopleura elegantissima (Cnidaria: Anthozoa: Hexacorallia: Actiniaria: Actiniidae). Collectively, we generated approximately 450,000 transcripts from acrorhagi of aggressive and non-aggressive polyps. For both transcriptomes we identified 65 candidate sea anemone toxin genes, representing phospholipase A2s, cytolysins, neurotoxins, and acrorhagins. When compared to previously characterized sea anemone toxin assemblages, each transcriptome revealed greater within-species sequence divergence across all toxin types. The transcriptome of the aggressive polyp had a higher abundance of type II voltage gated potassium channel toxins/Kunitz-type protease inhibitors and type II acrorhagins. Using toxin-like proteins from other venomous taxa, we also identified 612 candidate toxin-like transcripts with signaling regions, potentially unidentified secretory toxin-like proteins. Among these, metallopeptidases and cysteine rich (CRISP) candidate transcripts were in high abundance. Furthermore, our gene ontology analyses identified a high prevalence of genes associated with “blood coagulation” and “positive regulation of apoptosis”, as well as “nucleoside: sodium symporter activity” and “ion channel binding”. The resulting assemblage of expressed genes may represent synergistic proteins associated with toxins or proteins related to the morphology and behavior exhibited by the aggressive polyp. We implement a multifaceted approach to investigate the assemblage of expressed genes specifically within acrorhagi, specialized structures used only for intraspecific competition. By combining differential expression, phylogenetic, and gene ontology analyses, we identify several candidate toxins and other potentially important proteins in acrorhagi of A. elegantissima. Although not all of the toxins identified are used in intraspecific competition, our analysis highlights some candidates that may play a vital role in intraspecific competition. Our findings provide a framework for further investigation into components of venom used exclusively for intraspecific competition in acrorhagi-bearing sea anemones and potentially other venomous animals.

Sea anemones (Cnidaria, Anthozoa, and Actiniaria) use toxic peptides to incapacitate and immobilize prey and to deter potential predators. Their toxin arsenal is complex, targeting a variety of functionally important protein complexes and macromolecules involved in cellular homeostasis. Among these, actinoporins are one of the better characterized toxins; these venom proteins form a pore in cellular membranes containing sphingomyelin. We used a combined bioinformatic and phylogenetic approach to investigate how actinoporins have evolved across three superfamilies of sea anemones (Actinioidea, Metridioidea, and Actinostoloidea). Our analysis identified 90 candidate actinoporins across 20 species. We also found clusters of six actinoporin-like genes in five species of sea anemone (Nematostella vectensis, Stomphia coccinea, Epiactis japonica, Heteractis crispa, and Diadumene leucolena); these actinoporin-like sequences resembled actinoporins but have a higher sequence similarity with toxins from fungi, cone snails, and Hydra. Comparative analysis of the candidate actinoporins highlighted variable and conserved regions within actinoporins that may pertain to functional variation. Although multiple residues are involved in initiating sphingomyelin recognition and membrane binding, there is a high rate of replacement for a specific tryptophan with leucine (W112L) and other hydrophobic residues. Residues thought to be involved with oligomerization were variable, while those forming the phosphocholine (POC) binding site and the N-terminal region involved with cell membrane penetration were highly conserved.

The mitochondrial genome of basal animals is generally more slowly evolving than that of bilaterians. This difference in rate complicates the study of relationships among members of these lineages and the discovery of cryptic species or the testing of morphological species concepts within them. We explore the properties of mitochondrial and nuclear ribosomal genes in the cnidarian order Actiniaria, using both an ordinal- and familial-scale sample of taxa. Although the markers do not show significant incongruence, they differ in their phylogenetic informativeness and the kinds of relationships they resolve. Among the markers studied here, the fragments of 12S rDNA and 18S rDNA most effectively recover well-supported nodes; those of 16S rDNA and 28S rDNA are less effective. The general patterns we observed are similar to those in other hexacorallians, although Actiniaria alone show saturation of transitions for ordinal-scale analyses.