Abstract Animal venoms are theorized to evolve under the significant influence of positive Darwinian selection in a chemical arms race scenario, where the evolution of venom resistance in prey and the invention of potent venom in the secreting animal exert reciprocal selection pressures. Venom research to date has mainly focused on evolutionarily younger lineages, such as snakes and cone snails, while mostly neglecting ancient clades (e.g., cnidarians, coleoids, spiders and centipedes). By examining genome, venom-gland transcriptome and sequences from the public repositories, we report the molecular evolutionary regimes of several centipede and spider toxin families, which surprisingly accumulated low-levels of sequence variations, despite their long evolutionary histories. Molecular evolutionary assessment of over 3500 nucleotide sequences from 85 toxin families spanning the breadth of the animal kingdom has unraveled a contrasting evolutionary strategy employed by ancient and evolutionarily young clades. We show that the venoms of ancient lineages remarkably evolve under the heavy constraints of negative selection, while toxin families in lineages that originated relatively recently rapidly diversify under the influence of positive selection. We propose that animal venoms mostly employ a ‘two-speed’ mode of evolution, where the major influence of diversifying selection accompanies the earlier stages of ecological specialization (e.g., diet and range expansion) in the evolutionary history of the species–the period of expansion, resulting in the rapid diversification of the venom arsenal, followed by longer periods of purifying selection that preserve the potent toxin pharmacopeia–the period of purification and fixation. However, species in the period of purification may re-enter the period of expansion upon experiencing a major shift in ecology or environment. Thus, we highlight for the first time the significant roles of purifying and episodic selections in shaping animal venoms.

Author Summary While the influence of positive selection in diversifying animal venoms is widely recognized, the role of purifying selection that conserves the amino acid sequence of venom components such as peptide toxins has never been considered. In addition to unraveling the unique strategies of evolution of toxin gene families in centipedes and spiders, which are amongst the first terrestrial venomous lineages, we highlight the significant role of purifying selection in shaping the composition of animal venoms. Analysis of numerous toxin families, spanning the breadth of the animal kingdom, has revealed a striking contrast between the evolution of venom in ancient and evolutionarily young animal groups. Our findings enable the postulation of a new theory of venom evolution. The proposed ‘two-speed’ mode of evolution of venom captures the fascinating evolutionary history and the dynamics of this complex biochemical cocktail.

Citation: Sunagar K, Moran Y (2015) The Rise and Fall of an Evolutionary Innovation: Contrasting Strategies of Venom Evolution in Ancient and Young Animals. PLoS Genet 11(10): e1005596. https://doi.org/10.1371/journal.pgen.1005596 Editor: Hopi E. Hoekstra, Harvard University, UNITED STATES Received: May 28, 2015; Accepted: September 18, 2015; Published: October 22, 2015 Copyright: © 2015 Sunagar, Moran. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited Data Availability: All relevant data are within the paper and its Supporting Information files. Funding: KS was supported by a Marie Skłodowska-Curie Individual Fellowship (654294) (Link: http://ec.europa.eu/research/mariecurieactions/about-mca/actions/index_en.htm). This research was supported by the Israel Science Foundation grant No. 691/14 to YM (Link: http://www.isf.org.il/english/). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing interests: The authors have declared that no competing interests exist.

Introduction Venom is an intriguing evolutionary innovation that is utilized by various animals for predation and/or defense. This complex biochemical cocktail is characterized by a myriad of organic and inorganic molecules, such as proteins, peptides, polyamines and salts that disrupt the normal physiology of the envenomed animal. Evolution of venom has been intensively investigated in more recently diverged lineages (for simplicity, we refer to them as ‘evolutionarily younger’ lineages), such as advanced snakes and cone snails, which originated ~54 [1] and ~33–50 [2, 3] million years ago (MA), respectively. Several venom-encoding genes in these animals have undergone extensive duplications [4, 5] and evolve rapidly under the influence of positive selection [6–10]. In contrast, the evolution of venom in most of the ancient lineages, such as cnidarians (corals, sea anemones, hydroids and jellyfish), coleoids (octopus, squids and cuttlefish), spiders and centipedes, remains understudied, if not completely overlooked. Perhaps the only exhaustively investigated ancient venomous clade are the scorpions, which originated in the Silurian about 430 MA [11, 12]. Moreover, certain potent toxins in species separated by considerable geographic and genetic distance can exhibit remarkable sequence conservation (Fig 1). Yet, research to date has solely focused on how positive selection has expanded the venom arsenal, while completely ignoring the role of negative (purifying) selection. PPT PowerPoint slide

PowerPoint slide PNG larger image

larger image TIFF original image Download: Fig 1. Remarkable sequence conservation in distantly related toxins. Sequence alignments of widely separated sea anemone and scorpion neurotoxins are depicted. Sampled locations of these toxins are indicated on the map. Identical positions in sequence alignments are shown in blue, while differing amino acids are shown in brown. Uniprot IDs of sequences are: 1) B1NWR0; 2) P01532; 3) P0C5F4; 4) P29187; 5) E2S062; 6) Q7YXD3; 7) D5HR48; 8) P01484 and 9) D5HR56. https://doi.org/10.1371/journal.pgen.1005596.g001 Phylum Cnidaria consists of animals such as sea anemones, jellyfish, corals and hyrdroids that originated in the Ediacaran Period, approximately 600 MA [13–15]. They are characterized by unique stinging organelles called nematocysts with which they inject venom. Cnidaria represents the oldest venomous lineage known and includes some of the most notorious animals, such as the sea wasp (Chironex fleckeri), a species of box jellyfish. Coleoids, which first appeared in the Early Devonian 380–390 MA [16], represent yet another neglected lineage of ancient venomous animals. Although the venomous nature of coleoids was established as early as 1888 [17], their venoms have received scant attention from toxinological research [17–20]. Centipedes are amongst the oldest living terrestrial venomous animals, with the fossil record extending back to ~420 MA [21]. All ~3,300 species of centipedes (class Chilopoda) described to date belong to five extant orders: Craterostigmomorpha, Geophilomorpha, Lithobiomorpha, Scolopendromorpha and Scutigeromorpha. They inject venom into victims via modified first pair of trunk limbs (forcipules) and use venom for predation and defense. Venoms of certain centipedes can cause excruciating pain, paresthesia, edema, necrosis [22, 23] and can be fatal to mammals as large as dogs [24]. Yet, only a handful of centipede toxins have been pharmacologically characterized to date. Similarly, despite their remarkable ability to target a diversity of ion channels, only toxins from certain medically significant species of spiders have been investigated to date [25]. Thus, the evolutionary history and phyletic distribution of venom from these aforementioned ancient lineages, which represent the first venomous animal groups, remain understudied [18, 26–28]. It should be noted that the divergence times of these lineages can be safely assumed to be equivalent to the time of origin of venom in those respective lineages, as all of the examined lineages are (i) venomous, (ii) do not share between them a common venomous ancestor, and/or (iii) for most of them the fossil data clearly indicates the presence of a venom delivery apparatus [29–36]. By examining a large number of nucleotide sequences from a diversity of species, we report for the first time the molecular evolutionary histories of a number of venom protein families in centipedes, spiders and Toxicofera (clade of venomous snakes and lizards) lizards. In contrast to the rapid evolution of venom in evolutionarily younger lineages, we report an unusually high conservation of venom in centipedes and spiders, despite their long evolutionary histories. Moreover, molecular evolutionary assessments of toxin-encoding genes distributed across the tree of life, has unraveled a surprisingly strong influence of negative selection on the venoms of ancient animals. Our findings reveal contrasting trajectories of venom evolution in ancient and evolutionarily young clades, and emphasize the significant roles of purifying and episodic selections in shaping animal venoms. Further, these results enabled the postulation of a new model of venom evolution that captures their evolutionary dynamics, and the rise and fall in evolutionary rates of animal venoms.

Discussion The significance of purifying selection in the evolution venom Animal venoms are assumed to rapidly diversify under the unabated influence of positive Darwinian selection. They have been theorized to undergo a chemical arms race with prey animals, where the evolution of venom resistance in prey and the invention of efficient toxins in the predatory venomous animal exert reciprocal selection pressure [64], as postulated in the Red Queen hypothesis of Van Valen [65]. While the influence of positive selection is widely recognized, the role of purifying selection in shaping animal venoms has rarely been considered. Investigation of a large number of toxin-encoding gene families in this study has revealed a significant influence of negative selection on venom. Whilst positive selection increases the diversity of venom proteins, purifying selection probably aids in preserving the potency of the venom by filtering out mutations that negatively affect toxin efficiency. However, rare mutations that increase the potency of the venom arsenal (e.g., evolution of novel biochemical activity or increased binding efficiency) are likely to be propagated and preserved in the population. In the absence of a conservatory evolutionary force, neutral or positive selection could modify key residues and result in the reduction of potency or, for worse, the complete loss of bioactivity, which could severely decrease the fitness of the animal. Thus, purifying selection pressure appears to be vital for sustaining the potency and, consequently, shaping the animal venom arsenal. Certain toxins are more constrained than others It has been recently demonstrated that PFTs in Cnidaria, which bind to cell membranes and punch holes, evolve under the heavy constraints of negative selection [26, 66]. The lack of variation in this group of toxins, which includes several unrelated toxin types (e.g., aerolysin-related toxins in sea anemones, independently recruited hydralysins in hydroids, actinoporins and jellyfish toxins), was theorized to be a result of their complex multi-subunit packaging [67] and their ability to attack highly conserved molecular targets, such as cell membranes [26]. Toxins that undergo oligomerization in other classes of animals have also been noted to evolve relatively slowly as a result of structural constraints like the need to conserve sites involved in subunit interaction. While most 3FTxs in snake venoms diversify rapidly, κ-3FTxs, which undergo dimerization, were found to accumulate relatively fewer variations [59]. Similarly, toxins that may function in a ‘non-specific’ manner may also experience negative selection. Here, non-specificity of action is defined as the ability to target regions in a structural/biochemical property dependent (e.g., surface electrostatic charge) and target motif independent manner. For example, cytotoxic 3FTxs and β-defensin toxins—two very potent snake venom proteins, induce cytotoxicity by non-specifically binding to negatively charged cell membranes using hydrophobicity [68] and positively charged molecular surface [69], respectively. As a result, unlike most snake venom components, these proteins remain evolutionarily constrained [59, 62]. Similarly, scorpion lipolytic toxins were also theorized to be evolutionarily constrained because of their non-specific mechanism of action [58]. We found that β-PFTs in centipede venoms, which are similar to the aerolysin-like toxins, evolve under the significant influence of negative selection (Table 1). The lack of variation in this group of toxins may suggest that they either undergo oligomerization like their aerolysin homologues in other lineages or the possibility that they may employ a non-specific mechanism of action. A plot of site-specific ω against their respective amino acid positions reveals the extreme conservation of such toxin types that employ unique strategies for causing toxicity in envenomed animals (S2 Fig). As it allows the targeting of a wide variety of animals, the strategy of exerting toxic action non-specifically or by targeting highly conserved molecular sites, appears to be advantageous and follows a contrastingly different evolutionary regime in comparison to toxins that specialize in attacking highly plastic molecular receptors. The rise and fall of venom evolution A comparison of evolutionary regimes of ancient and evolutionarily younger lineages suggests a fascinating strategy of venom evolution. When venomous animals venture into novel ecological niches, they encounter new types of prey and predatory animals. Consequently, in order to adapt and conquer niches, they would need to fine-tune venom proteins to efficiently target these new animals. Several sites detected as episodically adaptive—i.e., sites that experience short bursts of adaptive selection, in these ancient clades may be reflective of such shifts in ecology. We propose that these earlier periods in the evolutionary history of a venomous species are accompanied by the significant influence of diversifying selection on the venom arsenal, which would expand the range of target sites and/or result in the origination of novel biochemical activities. This is particularly advantageous, since novel toxins generated may facilitate the efficient and rapid incapacitation of newly encountered prey and predatory animals. The period of expansion is followed by longer periods of purification, where the significant influence of negative selection preserves the potency of the toxin. Whenever there is a major shift in ecology or environment, the aforementioned stages of evolution repeat. Thus, we propose that venom-encoding genes mostly employ a ‘two-speed’ mode of evolution, where episodic diversifying selection accompanies the earlier stages of ecological specialization (e.g., diet and range expansion), resulting in the rapid diversification of the venom arsenal, followed by a longer period of purification and fixation that ensure the sustainability of venom potency. The low sequence variation in venom-encoding genes of ancient clades could be reflective of such long periods of purification and fine-tuning. In contrast, advanced snakes and cone snails, being evolutionarily very young, could still be undergoing the period of expansion and, consequently, exhibit a pronounced signature of positive Darwinian selection. However, it should be noted that the ‘two-speed’ model of evolution is likely applicable to venoms that serve predominantly predatory roles. Due to limited toxin sequence information from venoms that are employed for non-predatory functions (e.g., intraspecific competition in platypus, exclusively defensive roles in fishes, etc.), it remains to be seen whether they too follow our proposed evolutionary model. To conclude, in addition to unraveling the evolutionary regimes of toxin families in centipedes and spiders, which are amongst the first terrestrial venomous lineages, our findings highlight the pivotal roles of purifying and episodic selections in shaping animal venoms. Our findings enabled the postulation of a new theory of venom evolution in the animal kingdom that emphasizes the dynamic nature of these complex biochemical cocktails.

Methods Genome searches Toxin homologues were identified in the recently published genome of the coastal European centipede Strigamia maritima [70] by querying amino acid sequences of each toxin type against all six reading frames using the tblastn tool [71]. Evolutionary analyses Translated nucleotide sequences were aligned using MUSCLE 3.8 [72]. The best-fit model of nucleotide substitution for individual toxin datasets was determined according to the Akaike’s information criterion using jModeltest 2.1 [73] and model-averaged parameter estimates were used for the reconstruction of trees. Phylogenetic trees were built using PhyML 3.0 [74], where node support was evaluated with 1,000 bootstrapping replicates. Maximum-likelihood (ML) models [75] implemented in Codeml of the PAML package [76] were utilized to identify the influence of natural selection on toxin families[6]. As no a priori expectation exists, we compared likelihood values for a pair of models with different assumed ω distributions: M7 (β) versus M8 (β and ω) [77]. Only when the alternate model (M8) shows a better fit than the null model (M7) in the likelihood ratio test (LRT), are its results considered significant. LRT is estimated as twice the difference in ML values between the nested models, and is compared with the χ2 distribution with the appropriate degree of freedom—the difference in the number of parameters of the two models. Further, we used the Bayes empirical Bayes (BEB) approach [78] in M8 to detect amino acids under positive selection by calculating the posterior probability (PP) that a particular site belongs to a given selection class (neutral, conserved, or highly variable). Sites with PP ≥ 95% of belonging to the ‘‘ω > 1 class” are inferred to be positively selected. HyPhy’s [79] FUBAR approach [80] was used to detect sites evolving under pervasive diversifying and purifying selection pressures. MEME [81] was also employed to identify episodically diversifying sites. Sequence alignments used for selection assessments have been made available as a zipped file (S1 File; see S3 Table for accession list). Nucleotide substitution saturation was tested using DAMBE 5.5.9 [82] using the recommended protocol [83].

Acknowledgments We are thankful to Juan Calvete and Fernando Gonzalez Candelas (the University of València) for sharing data, and to Ariel Chipman and Joseph Heller (the Hebrew University) for their advice on centipede and molluscan natural history, respectively. We are also thankful to Volker Herzig and Sandy Pineda Gonzalez (the University of Queensland) for their advice on spider taxonomy and toxin nomenclature. We are thankful to Nicholas Casewell (Liverpool School of Tropical Medicine) and Sebastien Dutertre (The University of Montpellier) for discussions.

Author Contributions Conceived and designed the experiments: KS YM. Performed the experiments: KS. Analyzed the data: KS. Contributed reagents/materials/analysis tools: KS YM. Wrote the paper: KS YM.