Ontogenetic changes in venom composition have been described in Bothrops snakes, but only a few studies have attempted to identify the targeted paralogues or the molecular mechanisms involved in modifications of gene expression during ontogeny. In this study, we decoded B. jararacussu venom gland transcripts from six specimens of varying sizes and analyzed the variability in the composition of independent venom proteomes from 19 individuals. We identified 125 distinct putative toxin transcripts, and of these, 73 were detected in venom proteomes and only 10 were involved in the ontogenetic changes. Ontogenetic variability was linearly related to snake size and did not correspond to the maturation of the reproductive stage. Changes in the transcriptome were highly predictive of changes in the venom proteome. The basic myotoxic phospholipases A2 (PLA2s) were the most abundant components in larger snakes, while in venoms from smaller snakes, PIII-class SVMPs were the major components. The snake venom metalloproteinases (SVMPs) identified corresponded to novel sequences and conferred higher pro-coagulant and hemorrhagic functions to the venom of small snakes. The mechanisms modulating venom variability are predominantly related to transcriptional events and may consist of an advantage of higher hematotoxicity and more efficient predatory function in the venom from small snakes.

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.

We present here a comprehensive study on Bothrops jararacussu snake venom composition based on data obtained from 19 individuals of different sex, size, reproductive status, or geographical location. We carried out the first next-generation sequencing transcriptome (NGS) of venom glands from six individuals, which allowed a label-free characterization and quantification of venom toxins at the isoform level and evidencing the great diversity of PIII-class SVMPs in this venom. Venom variability was modulated mostly at transcriptional levels of a limited number of paralogues. As result, coagulant and hemorrhagic venoms from small snakes might be related to predatory function, while the venom from large snakes seems more related to immobilizing by high tissue-damaging myotoxic activity.