Juan Calvete

Understanding the origin and maintenance of phenotypic variation, particularly across a continuous spatial distribution, represents a key challenge in evolutionary biology. For this, animal venoms represent ideal study systems: they are complex, variable, yet easily quantifiable molecular phenotypes with a clear function. Rattlesnakes display tremendous variation in their venom composition, mostly through strongly dichotomous venom strategies, which may even coexist within a single species. Here, through dense, widespread population-level sampling of the Mojave rattlesnake, Crotalus scutulatus, we show that genomic structural variation at multiple loci underlies extreme  geographical variation in venom composition, which is maintained despite extensive gene flow. Unexpectedly, neither diet composition nor neutral population structure explain venom variation. Instead, venom divergence is strongly correlated with environmental conditions. Individual toxin genes correlate with distinct environmental factors, suggesting that different selective pressures can act on individual loci independently of their co-expression patterns or genomic proximity. Our results challenge common assumptions about diet composition as the key selective driver of snake venom evolution and emphasize how the interplay between genomic architecture and local scale spatial heterogeneity in selective pressures may facilitate the retention of adaptive functional polymorphisms across a continuous space.

Venomous snakes often display extensive variation in venom composition both between and within species. However, the mechanisms underlying the distribution of different toxins and venom types among populations and taxa remain insufficiently known. Rattlesnakes (Crotalus, Sistrurus) display extreme inter-and intraspecific variation in venom composition, centered particularly on the presence or absence of presynaptically neurotoxic phospholipases A 2 such as Mojave toxin (MTX). Interspecific hybridization has been invoked as a mechanism to explain the distribution of these toxins across rattlesnakes, with the implicit assumption that they are adaptively advantageous. Here, we test the potential of adaptive hybridization as a mechanism for venom evolution by assessing the distribution of genes encoding the acidic and basic subunits of Mojave toxin across a hybrid zone between MTX-positive Crotalus scutulatus and MTX-negative C. viridis in southwestern New Mexico, USA. Analyses of morphology, mitochondrial and single copy-nuclear genes document extensive admixture within a narrow hybrid zone. The genes encoding the two MTX subunits are strictly linked, and found in most hybrids and backcrossed individuals, but not in C. viridis away from the hybrid zone. Presence of the genes is invariably associated with presence of the corresponding toxin in the venom. We conclude that introgression of highly lethal neurotoxins through hybridization is not necessarily favored by natural selection in rattlesnakes, and that even extensive hybridization may not lead to introgression of these genes into another species.

The yellow-bellied sea snake, Pelamis platura, is the most broadly distributed snake species. Despite being endowed with a highly lethal venom, a proteomic analysis of its toxin composition was unavailable. The venoms of specimens collected in Golfo de Papagayo and Golfo Dulce (Costa Rica), where two distinctive color morphs occur, were chromatographically compared. The latter inhabits a fjord-like gulf where the transit of oceanic sea snakes into and from the basin is restricted, thus possibly affecting gene flow. RP-HPLC evidenced a conserved venom protein profile in both populations, despite their divergent color phenotypes. Following a trend observed in other sea snakes, P. platura venom is relatively simple, being composed of proteins of the three-finger toxin (3FTx), phospholipase A2 (PLA2), cysteine-rich secretory protein (CRISP), 5'-nucleotidase, and metalloproteinase families. The first three groups represent 49.9%, 32.9%, and 9.1% of total venom protein, respectively. The most abundant component (~26%) is pelamitoxin (P62388), a short-chain 3FTx, followed by a major basic PLA2 (~20%) and a group of three isoforms of CRISPs (~9%). Whereas isolated pelamitoxin was highly lethal to mice, neither the PLA2 nor the CRISP fraction caused death. However, the PLA2 rapidly increased plasma creatine kinase activity after intramuscular injection, indicating its myotoxic action. Differing from myotoxic PLA2s of viperids, this PLA2 was not cytolytic to murine myogenic cells in vitro, suggesting possible differences in its mechanism of action. The median lethal dose (LD50) estimates for P. platura crude venom in mice and in three species of fishes did not differ significantly. The sea snake antivenom manufactured by CSL Ltd. (Australia), which uses Enhydrina schistosa as immunogen, cross-recognized the three major components of P. platura venom and, accordingly, neutralized the lethal activity of crude venom and pelamitoxin, therefore being of potential usefulness in the treatment of envenomations by this species. Integrative analyses of animal venoms that combine the power of proteomics (venomics) with the characterization of their functional and immunological properties are significantly expanding knowledge on these remarkable bioweapons, both from a basic and a medical perspective. Costa Rica harbors a unique population of the yellow-bellied sea snake, Pelamis platura, that is restricted to a fjord-like gulf (Golfo Dulce). This population differs markedly from oceanic populations found elsewhere along the Pacific coast of this country, by presenting a patternless bright yellow coloration, instead of the typical bicolored or tricolored pattern of this species. It has been suggested that the dominance of this yellow-morph in Golfo Dulce might reflect gene flow restrictions, caused by the oceanographic conditions at this location. The present study demonstrates that the remarkable phenotypic variation between the two color morphs inhabiting Golfo Dulce and Golfo de Papagayo, respectively, is not associated with differences in the expression of venom components, as shown by their conserved RP-HPLC profiles. Proteomic analysis revealed the relatively simple toxin composition of P. platura venom, which contains three predominant types of proteins: three-finger toxins (protein abundance: 49.9%), phospholipases A2 (32.9%), and cysteine-rich secretory proteins (9.1%), together with few minor components. Further, the involvement of these most abundant proteins in the toxic effects of the venom, and their cross-recognition and neutralization by a sea snake antivenom produced against the venom of Enhydrina schistosa, were analyzed.

Bothriechis nigroviridis is an arboreal Neotropical pitviper found in Costa Rica and Panamá. A previous proteomic profiling of its venom revealed the presence of proteins with homology to the A and B subunits of crotoxin/Mojave toxin, a heterodimeric phospholipase A2 (PLA2) complex only described in rattlesnake venoms (genera Crotalus and Sistrurus). The native crotoxin-like heterodimer, named nigroviriditoxin, and its A and B subunits were isolated in the present work, and the complete amino acid sequence of the B subunit was determined. The purified A and B components were demonstrated to form a complex when reconstituted under native conditions. Nigroviriditoxin presents features similar to crotoxin, albeit displaying lower toxicity: the A component decreases the PLA2 activity of the B component, and increases its lethal potency in mice. Also in similarity to crotoxin B, nigroviriditoxin B induces myonecrosis. Its 122 amino acid sequence presents 81% identity with crotoxin B. Accordingly, nigroviriditoxin B was cross-recognized by equine antibodies from a Crotalus durissus terrificus antivenom. Phylogenetic analysis shows that the novel PLA2 from B. nigroviridis venom is basal to the branch including all the homologous PLA2 enzymes described in rattlesnakes, and more distant from PLA2s from Bothriechis species. Nigroviriditoxin is the first heterodimeric PLA2 complex found in a non-rattlesnake, Neotropical viperid venom, which displays structural, functional, and immunochemical similarities to crotoxin. The present findings are compatible with the existence of the particular structural trait of crotoxin-like molecules in New World pitvipers before the split of the Meso-South American and the Nearctic clades.

Seeds of Bauhinia forficata species were submited to biochemical characterization concerning fatty acids analysis, protein fractionization, and hemaglutinanting activity. The seed elementary analysis showed a high protein and lipids contents with 21.24% and 19.45% respectively. The more abundant fatty acid was linoleic acid with 46.47% of the lipidic fraction. With the exception of prolamins, the different proteic fractions (albumin, globulins, acid and basic glutelins) showed hemaglutinanting activity against rabbit red cells no treated and treated with proteolitic enzymes. The fraction acid glutelin showed the higher specific hemaglutinanting activity (1072.25 UH/mg P) against rabbit blood pre-treated with trypsin. Glutamin (16.20%) and Valin (11.07%) were the more abundant amino acids in the seeds. Therefore, B. forficata represent a possible optional source of food because exhibit a high energetic values.