The molecular basis of venom resistance in a rattlesnake‐squirrel predator‐prey system

Understanding how interspecific interactions mould the molecular basis of adaptations in coevolving species is a long‐sought goal of evolutionary biology. Venom in predators and venom resistance proteins in prey are coevolving molecular phenotypes, and while venoms are highly complex mixtures it is unclear if prey respond with equally complex resistance traits. Here, we use a novel molecular methodology based on protein affinity columns to capture and identify candidate blood serum resistance proteins (“venom interactive proteins” [VIPs]) in California Ground Squirrels (Otospermophilus beecheyi) that interact with venom proteins from their main predator, Northern Pacific Rattlesnakes (Crotalus o. oreganus). This assay showed that serum‐based resistance is both population‐ and species‐specific, with serum proteins from ground squirrels showing higher binding affinities for venom proteins of local snakes compared to allopatric individuals. Venom protein specificity assays identified numerous and diverse candidate prey resistance VIPs but also potential targets of venom in prey tissues. Many specific VIPs bind to multiple snake venom proteins and, conversely, single venom proteins bind multiple VIPs, demonstrating that a portion of the squirrel blood serum “resistome” involves broad‐based inhibition of nonself proteins and suggests that resistance involves a toxin scavenging mechanism. Analyses of rates of evolution of VIP protein homologues in related mammals show that most of these proteins evolve under purifying selection possibly due to molecular constraints that limit the evolutionary responses of prey to rapidly evolving snake venom proteins. Our method represents a general approach to identify specific proteins involved in co‐evolutionary interactions between species at the molecular level.     

Evolutionary History of the Photolyase/Cryptochrome Superfamily in Eukaryotes

Background

Photolyases and cryptochromes are evolutionarily related flavoproteins, which however perform distinct physiological functions. Photolyases (PHR) are evolutionarily ancient enzymes. They are activated by light and repair DNA damage caused by UV radiation. Although cryptochromes share structural similarity with DNA photolyases, they lack DNA repair activity. Cryptochrome (CRY) is one of the key elements of the circadian system in animals. In plants, CRY acts as a blue light receptor to entrain circadian rhythms, and mediates a variety of light responses, such as the regulation of flowering and seedling growth.

Results

We performed a comprehensive evolutionary analysis of the CRY/PHR superfamily. The superfamily consists of 7 major subfamilies: CPD class I and CPD class II photolyases, (6–4) photolyases, CRY-DASH, plant PHR2, plant CRY and animal CRY. Although the whole superfamily evolved primarily under strong purifying selection (average ω = 0.0168), some subfamilies did experience strong episodic positive selection during their evolution. Photolyases were lost in higher animals that suggests natural selection apparently became weaker in the late stage of evolutionary history. The evolutionary time estimates suggested that plant and animal CRYs evolved in the Neoproterozoic Era (~1000–541 Mya), which might be a result of adaptation to the major climate and global light regime changes occurred in that period of the Earth’s geological history.     

Phylogeny-Based Comparative Analysis of Venom Proteome Variation in a Clade of Rattlesnakes (Sistrurus sp.)

A long-standing question in evolutionary studies of snake venoms is the extent to which phylogenetic divergence and diet can account for between-species differences in venom composition. Here we apply phylogeny-based comparative methods to address this question. We use data on venom variation generated using proteomic techniques for all members of a small clade of rattlesnakes (Sistrurus sp.) and two outgroups for which phylogenetic and diet information is available. We first complete the characterization of venom variation for all members of this clade with a “venomic” analysis of pooled venoms from two members of this genus, S. milarius streckeri and S. m. milarius. These venoms exhibit the same general classes of proteins as those found in other Sistrurus species but differ in their relative abundances of specific protein families. We then test whether there is significant phylogenetic signal in the relative abundances of major venom proteins across species and if diet (measured as percent mammals and lizards among all prey consumed) covaries with venom composition after phylogenetic divergence is accounted for. We found no evidence for significant phylogenetic signal in venom variation: K values for seven snake venom proteins and two composite venom variables [PC 1 and 2]) were all nonsignificant and lower (mean = 0.11+0.06 sd) than mean K values (>0.35) previously reported for a wide range of morphological, life history, physiological and behavioral traits from other species. Finally, analyses based on Phylogenetic Generalized Least Squares (PGLS) methods reveal that variation in abundance of some venom proteins, most strongly CRISP is significantly related to snake diet. Our results demonstrate that venom variation in these snakes is evolutionarily a highly labile trait even among very closely-related taxa and that natural selection acting through diet variation may play a role in molding the relative abundance of specific venom proteins.     

Snake venomics of Crotalus tigris: the minimalist toxin arsenal of the deadliest Neartic rattlesnake venom. Evolutionary Clues for generating a pan-specific antivenom against crotalid type II venoms

We report the proteomic and antivenomic characterization of Crotalus tigris venom. This venom exhibits the highest lethality for mice among rattlesnakes and the simplest toxin proteome reported to date. The venom proteome of C. tigris comprises 7-8 gene products from 6 toxin families; the presynaptic β-neurotoxic heterodimeric PLA(2), Mojave toxin, and two serine proteinases comprise, respectively, 66 and 27% of the C. tigris toxin arsenal, whereas a VEGF-like protein, a CRISP molecule, a medium-sized disintegrin, and 1-2 PIII-SVMPs each represent 0.1-5% of the total venom proteome. This toxin profile really explains the systemic neuro- and myotoxic effects observed in envenomated animals. In addition, we found that venom lethality of C. tigris and other North American rattlesnake type II venoms correlates with the concentration of Mojave toxin A-subunit, supporting the view that the neurotoxic venom phenotype of crotalid type II venoms may be described as a single-allele adaptation. Our data suggest that the evolutionary trend toward neurotoxicity, which has been also reported for the South American rattlesnakes, may have resulted by pedomorphism. The ability of an experimental antivenom to effectively immunodeplete proteins from the type II venoms of C. tigris, Crotalus horridus , Crotalus oreganus helleri, Crotalus scutulatus scutulatus, and Sistrurus catenatus catenatus indicated the feasibility of generating a pan-American anti-Crotalus type II antivenom, suggested by the identification of shared evolutionary trends among South and North American Crotalus species.     

Phenotypic integration in the feeding system of the eastern diamondback rattlesnake (Crotalus adamanteus)

Selection can vary geographically across environments and temporally over the lifetime of an individual. Unlike geographic contexts, where different selective regimes can act on different alleles, age‐specific selection is constrained to act on the same genome by altering age‐specific expression. Snake venoms are exceptional traits for studying ontogeny because toxin expression variation directly changes the phenotype; relative amounts of venom components determine, in part, venom efficacy. Phenotypic integration is the dependent relationship between different traits that collectively produce a complex phenotype and, in venomous snakes, may include traits as diverse as venom, head shape and fang length. We examined the feeding system of the eastern diamondback rattlesnake (Crotalus adamanteus) across environments and over the lifetime of individuals and used a genotype–phenotype map approach, protein expression data and morphological data to demonstrate that: (i) ontogenetic effects explained more of the variation in toxin expression variation than geographic effects, (ii) both juveniles and adults varied geographically, (iii) toxin expression variation was a result of directional selection and (iv) different venom phenotypes covaried with morphological traits also associated with feeding in temporal (ontogenetic) and geographic (functional) contexts. These data are the first to demonstrate, to our knowledge, phenotypic integration between multiple morphological characters and a biochemical phenotype across populations and age classes. We identified copy number variation as the mechanism driving the difference in the venom phenotype associated with these morphological differences, and the parallel mitochondrial, venom and morphological divergence between northern and southern clades suggests that each clade may warrant classification as a separate evolutionarily significant unit.     

Phylogeny and Patterns of Diversity of Goat mtDNA Haplogroup A Revealed by Resequencing Complete Mitogenomes

We sequenced to near completion the entire mtDNA of 28 Sardinian goats, selected to represent the widest possible diversity of the most widespread mitochondrial evolutionary lineage, haplogroup (Hg) A. These specimens were reporters of the diversity in the island but also elsewhere, as inferred from their affiliation to each of 11 clades defined by D-loop variation. Two reference sequences completed the dataset. Overall, 206 variations were found in the full set of 30 sequences, of which 23 were protein-coding non-synonymous single nucleotide substitutions. Many polymorphic sites within Hg A were informative for the reconstruction of its internal phylogeny. Bayesian and network clustering revealed a general similarity over the entire molecule of sequences previously assigned to the same D-loop clade, indicating evolutionarily meaningful lineages. Two major sister groupings emerged within Hg A, which parallel distinct geographical distributions of D-loop clades in extant stocks. The pattern of variation in protein-coding genes revealed an overwhelming role of purifying selection, with the quota of surviving variants approaching neutrality. However, a simple model of relaxation of selection for the bulk of variants here reported should be rejected. Non-synonymous diversity of Hg''s A, B and C denoted that a proportion of variants not greater than that allowed in the wild was given the opportunity to spread into domesticated stocks. Our results also confirmed that a remarkable proportion of pre-existing Hg A diversity became incorporated into domestic stocks. Our results confirm clade A11 as a well differentiated and ancient lineage peculiar of Sardinia.     

Assembling an arsenal: origin and evolution of the snake venom proteome inferred from phylogenetic analysis of toxin sequences

We analyzed the origin and evolution of snake venom toxin families represented in both viperid and elapid snakes by means of phylogenetic analysis of the amino acid sequences of the toxins and related nonvenom proteins. Out of eight toxin families analyzed, five provided clear evidence of recruitment into the snake venom proteome before the diversification of the advanced snakes (Kunitz-type protease inhibitors, CRISP toxins, galactose-binding lectins, M12B peptidases, nerve growth factor toxins), and one was equivocal (cystatin toxins). In two others (phospholipase A(2) and natriuretic toxins), the nonmonophyly of venom toxins demonstrates that presence of these proteins in elapids and viperids results from independent recruitment events. The ANP/BNP natriuretic toxins are likely to be basal, whereas the CNP/BPP toxins are Viperidae only. Similarly, the lectins were recruited twice. In contrast to the basal recruitment of the galactose-binding lectins, the C-type lectins were shown to be Viperidae only, with the alpha-chains and beta-chains resulting from an early duplication event. These results provide strong additional evidence that venom evolved once, at the base of the advanced snake radiation, rather than multiple times in different lineages, with these toxins also present in the venoms of the "colubrid" snake families. Moreover, they provide a first insight into the composition of the earliest ophidian venoms and point the way toward a research program that could elucidate the functional context of the evolution of the snake venom proteome.     

Identification and characterization of a taxon-specific three-finger toxin from the venom of the Green Vinesnake (Oxybelis fulgidus; family Colubridae)

Snake venoms contain a variety of protein and peptide toxins, and the three-finger toxins (3FTxs) are among the best characterized family of venom proteins. The compact nature and highly conserved molecular fold of 3FTxs, together with their abundance in many venoms, has contributed to their utility in structure-function studies. Although many target the nicotinic acetylcholine receptor of vertebrate skeletal muscle, often binding with nanomolar K_ds, several non-conventional 3FTxs show pronounced taxon-specific neurotoxic effects. Here we describe the purification and characterization of fulgimotoxin, a monomeric 3FTx from the venom of Oxybelis fulgidus, a neotropical rear-fanged snake. Fulgimotoxin retains the canonical 5 disulfides of the non-conventional 3FTxs and is highly neurotoxic to lizards; however, mice are unaffected, demonstrating that this toxin is taxon-specific in its effects. Analysis of structural features of fulgimotoxin and other colubrid venom 3FTxs indicate the presence of a “colubrid toxin motif” (CYTLY) and a second conserved segment (WAVK) found in Boiga and Oxybelis taxon-specific 3FTxs, both in loop II. Because specific residues in loop II conventional α-neurotoxic 3FTxs are intimately associated with receptor binding, we hypothesize that this loop, with its highly conserved substitutions, confers taxon-specific neurotoxicity. These findings underscore the importance of rear-fanged snake venoms for understanding the evolution of toxin molecules and demonstrate that even among well-characterized toxin families, novel structural and functional motifs may be found.     

Reconstructing the Origins of the Somatostatin and Allatostatin-C Signaling Systems Using the Accelerated Evolution of Biodiverse Cone Snail Toxins

Somatostatin and its related peptides (SSRPs) form an important family of hormones with diverse physiological roles. The ubiquitous presence of SSRPs in vertebrates and several invertebrate deuterostomes suggests an ancient origin of the SSRP signaling system. However, the existence of SSRP genes outside of deuterostomes has not been established, and the evolutionary history of this signaling system remains poorly understood. Our recent discovery of SSRP-like toxins (consomatins) in venomous marine cone snails (Conus) suggested the presence of a related signaling system in mollusks and potentially other protostomes. Here, we identify the molluscan SSRP-like signaling gene that gave rise to the consomatin family. Following recruitment into venom, consomatin genes experienced strong positive selection and repeated gene duplications resulting in the formation of a hyperdiverse family of venom peptides. Intriguingly, the largest number of consomatins was found in worm-hunting species (>400 sequences), indicating a homologous system in annelids, another large protostome phylum. Consistent with this, comprehensive sequence mining enabled the identification of SSRP-like sequences (and their corresponding orphan receptor) in annelids and several other protostome phyla. These results established the existence of SSRP-like peptides in many major branches of bilaterians and challenge the prevailing hypothesis that deuterostome SSRPs and protostome allatostatin-C are orthologous peptide families. Finally, having a large set of predator–prey SSRP sequences available, we show that although the cone snail’s signaling SSRP-like genes are under purifying selection, the venom consomatin genes experience rapid directional selection to target receptors in a changing mix of prey.     

Recruitment and diversification of an ecdysozoan family of neuropeptide hormones for black widow spider venom expression

Venoms have attracted enormous attention because of their potent physiological effects and dynamic evolution, including the convergent recruitment of homologous genes for venom expression. Here we provide novel evidence for the recruitment of genes from the Crustacean Hyperglycemic Hormone (CHH) and arthropod Ion Transport Peptide (ITP) superfamily for venom expression in black widow spiders. We characterized latrodectin peptides from venom gland cDNAs from the Western black widow spider (Latrodectus hesperus), the brown widow (Latrodectus geometricus) and cupboard spider (Steatoda grossa). Phylogenetic analyses of these sequences with homologs from other spider, scorpion and wasp venom cDNAs, as well as CHH/ITP neuropeptides, show latrodectins as derived members of the CHH/ITP superfamily. These analyses suggest that CHH/ITP homologs are more widespread in spider venoms, and were recruited for venom expression in two additional arthropod lineages. We also found that the latrodectin 2 gene and nearly all CHH/ITP genes include a phase 2 intron in the same position, supporting latrodectin's placement within the CHH/ITP superfamily. Evolutionary analyses of latrodectins suggest episodes of positive selection along some sequence lineages, and positive and purifying selection on specific codons, supporting its functional importance in widow venom. We consider how this improved understanding of latrodectin evolution informs functional hypotheses regarding its role in black widow venom as well as its potential convergent recruitment for venom expression across arthropods.     

Molecular Phylogeny and Evolution of the Proteins Encoded by Coleoid (Cuttlefish, Octopus, and Squid) Posterior Venom Glands

In this study, we report for the first time a detailed evaluation of the phylogenetic history and molecular evolution of the major coleoid toxins: CAP, carboxypeptidase, chitinase, metalloprotease GON-domain, hyaluronidase, pacifastin, PLA2, SE-cephalotoxin and serine proteases, with the carboxypeptidase and GON-domain documented for the first time in the coleoid venom arsenal. We show that although a majority of sites in these coleoid venom-encoding genes have evolved under the regime of negative selection, a very small proportion of sites are influenced by the transient selection pressures. Moreover, nearly 70 % of these episodically adapted sites are confined to the molecular surface, highlighting the importance of variation of the toxin surface chemistry. Coleoid venoms were revealed to be as complex as other venoms that have traditionally been the recipient of the bulk of research efforts. The presence of multiple peptide/protein types in coleoids similar to those present in other animal venoms identifies a convergent strategy, revealing new information as to what characteristics make a peptide/protein type amenable for recruitment into chemical arsenals. Coleoid venoms have significant potential not only for understanding fundamental aspects of venom evolution but also as an untapped source of novel toxins for use in drug design and discovery.     

Hemorrhagic and mojave toxins in the venoms of the offspring of two mojave rattlesnakes (Crotalus scutulatus scutulatus)

• 1.1. The venoms of two Mojave rattlesnakes and those of their offsprings were analyzed for Mojave toxin and hemorrhagic toxin.
• 2.2. The venom of one female, collected in Pima County, Arizona, and the venoms of her six offspring contained hemorrhagic toxin but not Mojave toxin (venom B).
• 3.3. The venom of the second female, captured in El Paso County, Texas, contained both toxins (A + B venom). Of her 10 offspring, five contained venom with both toxins, two had hemorrhagic toxin only, and three contained neither toxin.
• 4.4. Venoms that caused hemorrhage also inactivated complement. A pool of the venoms of the venom B offspring was less toxic than adult pooled venom A.

     

Snake venomics of Central American pitvipers: clues for rationalizing the distinct envenomation profiles of Atropoides nummifer and Atropoides picadoi

We report the proteomic characterization of the Central American pitvipers Atropoides nummifer and Atropoides picadoi. The crude venoms were fractionated by reverse-phase high-performance liquid chromatography (HPLC), followed by analysis of each chromatographic fraction by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE), N-terminal sequencing, matrix-assisted laser desorption ionization-time of flight (MALDI-TOF) mass fingerprinting, and collision-induced dissociation-tandem mass spectrometry (CID-MS/MS) of tryptic peptides. Each venom contained a number of bradykinin-potentiating peptides and around 25-27 proteins of molecular masses in the range of 7-112 kDa, belonging to only nine different toxin families (disintegrin, DC fragment, snake venom vascular endothelial growth factor, phospholipases A2, serine protease, cysteine-rich secretory proteins, C-type lectins, L-amino acid oxidase, and Zn2+-dependent metalloproteases), albeit distinctly distributed among the two Atropoides species. In addition, A. nummifer expresses low amounts of a three-finger toxin not detected in the venom of A. picadoi. The major toxins of A. nummifer belong to the PLA2 (relative abundance, 36.5%) and the serine proteinase (22%) families, whereas the most abundant A. picadoi toxins are Zn2+-dependent metalloproteinases (66.4%). We estimate that the similarity of venom proteins between the two Atropoides taxa may be around 14-16%. The high degree of differentiation in the venom proteome among congeneric taxa emphasizes unique aspects of venom composition of related species of Atropoides snakes and points to a strong role for adaptive diversification via natural selection as a cause of this distinctiveness. On the other hand, their distinct venom toxin compositions provide clues for rationalizing the low hemorrhagic, coagulant, and defibrinating activities and the high myotoxic and proteolytic effects evoked by A. nummifer snakebite in comparison to other crotaline snake venoms and the high hemorrhagic activity of A. picadoi.