Novel phenotypes are commonly associated with gene duplications and neofunctionalization, less documented are the cases of phenotypic maintenance through the recruitment of novel genes. Proteolysis is the primary toxic character of many snake venoms, and ADAM metalloproteinases, named Snake Venom Metalloproteinases (SVMPs), are largely recognized as the major effectors of this phenotype. However, by investigating original transcriptomes from 58 species of advanced snakes (Caenophidia) across their phylogeny, we discovered that a different enzyme, matrix metalloproteinase (MMP), is actually the dominant venom component in three tribes (Tachymenini, Xenodontini, and Conophiini) of rear-fanged snakes (Dipsadidae). Proteomic and functional analyses of these venoms further indicate that MMPs are likely playing an ‘SVMP-like’ function in the proteolytic phenotype. A detailed look into the venom-specific sequences revealed a new highly expressed MMP subtype, named snake venom MMP (svMMP), which originated independently on at least three occasions from an endogenous MMP-9. We further show that by losing ancillary non-catalytic domains present in its ancestors, svMMPs followed an evolutionary path toward a simplified structure during their expansion in the genomes, thus paralleling what has been proposed for the evolution of their Viperidae counterparts, the SVMPs. Moreover, we inferred an inverse relationship between the expression of svMMPs and SVMPs along the evolutionary history of Xenodontinae, pointing out that one type of enzyme may be substituting for the other, while the general (metallo)proteolytic phenotype is maintained. These results provide rare evidence on how relevant phenotypic traits can be optimized via natural selection on non-homologous genes, yielding alternate biochemical components.

Snakes of the Philodryadini tribe are included in the Dipsadidae family, which is a diverse group of rear-fanged snakes widespread in different ecological conditions, including habitats and diet. However, little is known about the composition and effects of their venoms despite their relevance for understanding the evolution of these snakes or even their impact on the occasional cases of human envenoming. In this study, we integrated venom gland transcriptomics, venom proteomics and functional assays to characterize the venoms from eight species of the Philodryadini tribe, which includes the genus Philodryas, Chlorosoma and Xenoxybelis. The most abundant components identified in the venoms were snake venom metalloproteinases (SVMPs), cysteine-rich secretory proteins (CRISPs), C-type lectins (CTLs), snake endogenous matrix metalloproteinases type 9 (seMMP-9) and snake venom serinoproteinases (SVSPs). These protein families showed a variable expression profile in each genus. SVMPs were the most abundant components in Philodryas, while seMMP-9 and CRISPs were the most expressed in Chlorosoma and Xenoxybelis, respectively. Lineage-specific differences in venom composition were also observed among Philodryas species, whereas P. olfersii presented the highest amount of SVSPs and P. agassizii was the only species to express significant amounts of 3FTx. The variability observed in venom composition was confirmed by the venom functional assays. Philodryas species presented the highest SVMP activity, whereas Chlorosoma species showed higher levels of gelatin activity, which may correlate to the seMMP-9 enzymes. The variability observed in the composition of these venoms may be related to the tribe phylogeny and influenced by their diets. In the presented study, we expanded the set of venomics studies of the Philodryadini tribe, which paves new roads for further studies on the evolution and ecology of Dipsadidae snakes.

No abstract

Snake venom harbors a wide and diverse array of enzymatic and nonenzymatic toxic components, allowing them to exert myriad effects on their prey. However, they appear to trend toward a few optimal compositional scaffolds, dominated by four major toxin classes: SVMPs, SVSPs, 3FTxs and PLA2s. Nevertheless, the latter appears to be restricted to vipers and elapids, as it has never been reported as a major venom component in rear-fanged species. Here, by investigating the original transcriptomes from 19 species distributed in eight genera from the Pseudoboini tribe (Dipsadidae: Xenodontinae) and screening among seven additional tribes of Dipsadidae and three additional families of advanced snakes, we discovered that a novel type of venom, PLA2, resembling a PLA2-IIE, has been recruited to the venom of some species of the Pseudoboini tribe, where it is a major component. Proteomic and functional analyses of these venoms further indicate that these PLA2s play a relevant role in the venoms from this tribe. Moreover, we reconstructed the phylogeny of PLA2s across different snake groups and show that different types of these toxins have been recruited in at least five independent events in caenophidian snakes. Additionally, we present the first compositional profiling of Pseudoboini venoms. Our results demonstrate how relevant phenotypic traits are convergently recruited by different means and from homologous and nonhomologous genes in phylogenetically and ecologically divergent snake groups, possibly optimizing venom composition to overcome diverse adaptative landscapes.