Phylogeny-based comparative analysis of venom proteome variation in a clade of rattlesnakes (Sistrurus sp.)

• Published in: PloS one 2013-06, Vol.8 (6), p.e67220-e67220
• Main Authors: Gibbs, H Lisle, Sanz, Libia, Sovic, Michael G, Calvete, Juan J
• Format: Article
• Language: English
• Subjects: Abundance; Agkistrodon contortrix; Animals; Biology; Bothrops jararaca; Chromatography, High Pressure Liquid; Chromatography, Reverse-Phase; Comparative analysis; Crotalid Venoms - metabolism; Crotalus - metabolism; Crotalus atrox; Diet; Divergence; Evolution; Evolution, Molecular; Genes; Genetic Loci - genetics; Hypotheses; Least-Squares Analysis; Life history; Lizards; Natural selection; Organismal biology; Peptides; Phylogenetics; Phylogeny; Prey; Protein families; Proteins; Proteome - metabolism; Proteomics; Relative abundance; Sistrurus catenatus; Sistrurus miliarius; Snakes; Taxa; Variation; Venom; Venom toxins; Viperidae
• ISSN: 1932-6203; 1932-6203
• DOI: 10.1371/journal.pone.0067220

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.