Background: Understanding the processes that drive the evolution of snake venom is a topic of great research interest in molecular and evolutionary toxinology. Recent studies suggest that ontogenetic changes in venom composition are genetically controlled rather than environmentally induced. However, the molecular mechanisms underlying these changes remain elusive. Here we have explored the basis and level of regulation of the ontogenetic shift in the venom composition of the Central American rattlesnake,Crotalus s. simususing a combined proteomics and transcriptomics approach.Results: Proteomic analysis showed that the ontogenetic shift in the venom composition ofC. s. simusis essentially characterized by a gradual reduction in the expression of serine proteinases and PLA₂molecules, particularly crotoxin, a β-neurotoxic heterodimeric PLA₂, concominantly with an increment of PI and PIII metalloproteinases at age 9–18 months. Comparison of the transcriptional activity of the venom glands of neonate and adultC. s. simusspecimens indicated that their transcriptomes exhibit indistinguisable toxin family profiles, suggesting that the elusive mechanism by which shared transcriptomes generate divergent venom phenotypes may operate post-transcriptionally. Specifically, miRNAs with frequency count of 1000 or greater exhibited an uneven distribution between the newborn and adult datasets. Of note, 590 copies of a miRNA targeting crotoxin B-subunit was exclusively found in the transcriptome of the adult snake, whereas 1185 copies of a miRNA complementary to a PIII-SVMP mRNA was uniquely present in the newborn dataset. These results support the view that age-dependent changes in the concentration of miRNA modulating the transition from a crotoxin-rich to a SVMP-rich venom from birth through adulhood can potentially explain what is observed in the proteomic analysis of the ontogenetic changes in the venom composition ofC. s. simus.Conclusions: Existing snake venom toxins are the result of early recruitment events in the Toxicofera clade of reptiles by which ordinary genes were duplicated, and the new genes selectively expressed in the venom gland and amplified to multigene families with extensive neofunctionalization throughout the approximately 112–125 million years of ophidian evolution. Our findings support the view that understanding the phenotypic diversity of snake venoms requires a deep knowledge of the mechanisms regulating the transcriptional and translational activity of the venom gland. Our results suggest a functional role for miRNAs. The impact of specific miRNAs in the modulation of venom composition, and the integration of the mechanisms responsible for the generation of these miRNAs in the evolutionary landscape of the snake's venom gland, are further challenges for future research.

Snakebite envenoming is a neglected tropical disease that may claim over 100,000 human lives annually worldwide. Snakebite occurs as the result of an interaction between a human and a snake that elicits either a defensive response from the snake or, more rarely, a feeding response as the result of mistaken identity. Snakebite envenoming is therefore a biological and, more specifically, an ecological problem. Snake venom itself is often described as a “cocktail”, as it is a heterogenous mixture of molecules including the toxins (which are typically proteinaceous) responsible for the pathophysiological consequences of envenoming. The primary function of venom in snake ecology is pre-subjugation, with defensive deployment of the secretion typically considered a secondary function. The particular composition of any given venom cocktail is shaped by evolutionary forces that include phylogenetic constraints associated with the snake's lineage and adaptive responses to the snake's ecological context, including the taxa it preys upon and by which it is predated upon. In the present article, we describe how conceptual frameworks from ecology and evolutionary biology can enter into a mutually enlightening relationship with clinical toxinology by enabling the consideration of snakebite envenoming from an “ecological stance”. We detail the insights that may emerge from such a perspective and highlight the ways in which the high-fidelity descriptive knowledge emerging from applications of -omics era technologies – “venomics” and “antivenomics” – can combine with evolutionary explanations to deliver a detailed understanding of this multifactorial health crisis.

The metacommunity concept provides a theoretical framework that aims at explaining organism distributions by a combination of environmental filtering, dispersal, and drift. However, few works have attempted a multitaxon approach and even fewer have compared two distant biogeographical regions using the same methodology. We tested the expectation that temperate (mediterranean‐climate) pond metacommunities would be more influenced by environmental and spatial processes than tropical ones, because of stronger environmental gradients and a greater isolation of waterbodies. However, the pattern should be different among groups of organisms depending on their dispersal abilities. We surveyed 30 tropical and 32 mediterranean temporary ponds from Costa Rica and Spain, respectively, and obtained data on 49 environmental variables. We characterized the biological communities of bacteria and archaea (from the water column and the sediments), phytoplankton, zooplankton, benthic invertebrates, amphibians and birds, and estimated the relative role of space and environment on metacommunity organization for each group and region, by means of variation partitioning using generalized additive models. Purely environmental effects were important in both tropical and mediterranean ponds, but stronger in the latter, probably due to their larger limnological heterogeneity. Spatially correlated environment and pure spatial effects were greater in the tropics, related to higher climatic heterogeneity and dispersal processes (e.g., restriction, surplus) acting at different scales. The variability between taxonomic groups in the contribution of spatial and environmental factors to metacommunity variation was very wide, but higher in active, compared with passive, dispersers. Higher environmental effects were observed in mediterranean passive dispersers, and higher spatial effects in tropical passive dispersers. The unexplained variation was larger in the tropical setting, suggesting a higher role for stochastic processes, unmeasured environmental factors, or biotic interactions in the tropics, although this difference affected some actively dispersing groups (insects and birds) more than passive dispersers. These results, despite our limitations in comparing only two regions, provide support, for a wide variety of aquatic organisms, for the classic view of stronger abiotic niche constraints in temperate areas compared with the tropics. The heterogeneous response of taxonomic groups between regions also points to a stronger influence of regional context than organism adaptations on metacommunity organization.

The latitudinal diversity gradient predicts that tropical regions should have higher alpha, beta, and gamma diversity than temperate areas. However, only a few studies have assessed the temporal variability of the different components of diversity across climatic regions. In this study, we compare, using a spatial and temporal approach, the diversity of multiple taxa inhabiting tropical and Mediterranean temporary ponds. We sampled the biological communities of each set of ponds on three occasions during the same hydrological year. Under a spatial framework, we analyzed, alpha, beta, and gamma diversities. With a temporal approach, we compared the coefficients of variation in alpha diversity for each local community, and temporal beta diversity. Differences between regions and sampling periods were tested using generalized linear mixed models. We found higher gamma and alpha diversity in the tropical ponds, as expected given the latitudinal differences between them. However, phytoplankton and microinvertebrates from the Mediterranean region, matched or even exceeded tropical alpha diversity on some occasions. Spatial beta diversity did not differ between regions, and it showed lower values at the middle or the end of the hydroperiod in bacteria, micro‐ and macroinvertebrates and amphibians. Thus, processes homogenizing and heterogenising pond metacommunities must be balanced in both studied regions. Temporal variation in alpha and beta diversity was similar for ponds in both regions, except for macroinvertebrates and amphibians, suggesting differential effects on community variation observable only in animals with longer life‐spans, at our temporal scale of analysis.