skip to main content

Title: Phylogenetically diverse diets favor more complex venoms in North American pitvipers

The role of natural selection in the evolution of trait complexity can be characterized by testing hypothesized links between complex forms and their functions across species. Predatory venoms are composed of multiple proteins that collectively function to incapacitate prey. Venom complexity fluctuates over evolutionary timescales, with apparent increases and decreases in complexity, and yet the causes of this variation are unclear. We tested alternative hypotheses linking venom complexity and ecological sources of selection from diet in the largest clade of front-fanged venomous snakes in North America: the rattlesnakes, copperheads, cantils, and cottonmouths. We generated independent transcriptomic and proteomic measures of venom complexity and collated several natural history studies to quantify dietary variation. We then constructed genome-scale phylogenies for these snakes for comparative analyses. Strikingly, prey phylogenetic diversity was more strongly correlated to venom complexity than was overall prey species diversity, specifically implicating prey species’ divergence, rather than the number of lineages alone, in the evolution of complexity. Prey phylogenetic diversity further predicted transcriptomic complexity of three of the four largest gene families in viper venom, showing that complexity evolution is a concerted response among many independent gene families. We suggest that the phylogenetic diversity of prey measures functionally relevant divergence in the targets of venom, a claim supported by sequence diversity in the coagulation cascade targets of venom. Our results support the general concept that the diversity of species in an ecological community is more important than their overall number in determining evolutionary patterns in predator trait complexity.

more » « less
Award ID(s):
1638872 1638879 1822417 1638902
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; « less
Publisher / Repository:
Proceedings of the National Academy of Sciences
Date Published:
Journal Name:
Proceedings of the National Academy of Sciences
Page Range / eLocation ID:
Article No. e2015579118
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Background The explosive radiation and diversification of the advanced snakes (superfamily Colubroidea) was associated with changes in all aspects of the shared venom system. Morphological changes included the partitioning of the mixed ancestral glands into two discrete glands devoted for production of venom or mucous respectively, as well as changes in the location, size and structural elements of the venom-delivering teeth. Evidence also exists for homology among venom gland toxins expressed across the advanced snakes. However, despite the evolutionary novelty of snake venoms, in-depth toxin molecular evolutionary history reconstructions have been mostly limited to those types present in only two front-fanged snake families, Elapidae and Viperidae. To have a broader understanding of toxins shared among extant snakes, here we first sequenced the transcriptomes of eight taxonomically diverse rear-fanged species and four key viperid species and analysed major toxin types shared across the advanced snakes. Results Transcriptomes were constructed for the following families and species: Colubridae - Helicops leopardinus , Heterodon nasicus , Rhabdophis subminiatus ; Homalopsidae – Homalopsis buccata ; Lamprophiidae - Malpolon monspessulanus , Psammophis schokari , Psammophis subtaeniatus , Rhamphiophis oxyrhynchus ; and Viperidae – Bitis atropos , Pseudocerastes urarachnoides , Tropidolaeumus subannulatus , Vipera transcaucasiana . These sequences were combined with those from available databases of other species in order to facilitate a robust reconstruction of the molecular evolutionary history of the key toxin classes present in the venom of the last common ancestor of the advanced snakes, and thus present across the full diversity of colubroid snake venoms. In addition to differential rates of evolution in toxin classes between the snake lineages, these analyses revealed multiple instances of previously unknown instances of structural and functional convergences. Structural convergences included: the evolution of new cysteines to form heteromeric complexes, such as within kunitz peptides (the beta-bungarotoxin trait evolving on at least two occasions) and within SVMP enzymes (the P-IIId trait evolving on at least three occasions); and the C-terminal tail evolving on two separate occasions within the C-type natriuretic peptides, to create structural and functional analogues of the ANP/BNP tailed condition. Also shown was that the de novo evolution of new post-translationally liberated toxin families within the natriuretic peptide gene propeptide region occurred on at least five occasions, with novel functions ranging from induction of hypotension to post-synaptic neurotoxicity. Functional convergences included the following: multiple occasions of SVMP neofunctionalised in procoagulant venoms into activators of the clotting factors prothrombin and Factor X; multiple instances in procoagulant venoms where kunitz peptides were neofunctionalised into inhibitors of the clot destroying enzyme plasmin, thereby prolonging the half-life of the clots formed by the clotting activating enzymatic toxins; and multiple occasions of kunitz peptides neofunctionalised into neurotoxins acting on presynaptic targets, including twice just within Bungarus venoms. Conclusions We found novel convergences in both structural and functional evolution of snake toxins. These results provide a detailed roadmap for future work to elucidate predator–prey evolutionary arms races, ascertain differential clinical pathologies, as well as documenting rich biodiscovery resources for lead compounds in the drug design and discovery pipeline. 
    more » « less
  2. Yoder, Anne (Ed.)
    Abstract Understanding the joint roles of protein sequence variation and differential expression during adaptive evolution is a fundamental, yet largely unrealized goal of evolutionary biology. Here, we use phylogenetic path analysis to analyze a comprehensive venom-gland transcriptome dataset spanning three genera of pitvipers to identify the functional genetic basis of a key adaptation (venom complexity) linked to diet breadth (DB). The analysis of gene-family-specific patterns reveals that, for genes encoding two of the most important venom proteins (snake venom metalloproteases and snake venom serine proteases), there are direct, positive relationships between sequence diversity (SD), expression diversity (ED), and increased DB. Further analysis of gene-family diversification for these proteins showed no constraint on how individual lineages achieved toxin gene SD in terms of the patterns of paralog diversification. In contrast, another major venom protein family (PLA2s) showed no relationship between venom molecular diversity and DB. Additional analyses suggest that other molecular mechanisms—such as higher absolute levels of expression—are responsible for diet adaptation involving these venom proteins. Broadly, our findings argue that functional diversity generated through sequence and expression variations jointly determine adaptation in the key components of pitviper venoms, which mediate complex molecular interactions between the snakes and their prey. 
    more » « less
  3. Traits for prey acquisition form the phenotypic interface of predator–prey interactions. In venomous predators, morphological variation in venom delivery apparatus like fangs and stingers may be optimized for dispatching prey. Here, we determine how a single dimension of venom injection systems evolves in response to variation in the size, climatic conditions and dietary ecology of viperid snakes. We measured fang length in more than 1900 museum specimens representing 199 viper species (55% of recognized species). We find both phylogenetic signal and within-clade variation in relative fang length across vipers suggesting both general taxonomic trends and potential adaptive divergence in fang length. We recover positive evolutionary allometry and little static allometry in fang length. Proportionally longer fangs have evolved in larger species, which may facilitate venom injection in more voluminous prey. Finally, we leverage climatic and diet data to assess the global correlates of fang length. We find that models of fang length evolution are improved through the inclusion of both temperature and diet, particularly the extent to which diets are mammal-heavy diets. These findings demonstrate how adaptive variation can emerge among components of complex prey capture systems. 
    more » « less
  4. Abstract Convergent evolution is pervasive in nature, but it is poorly understood how various constraints and natural selection limit the diversity of evolvable phenotypes. Here, we analyze the transcriptome across fruiting body development to understand the independent evolution of complex multicellularity in the two largest clades of fungi—the Agarico- and Pezizomycotina. Despite >650 My of divergence between these clades, we find that very similar sets of genes have convergently been co-opted for complex multicellularity, followed by expansions of their gene families by duplications. Over 82% of shared multicellularity-related gene families were expanding in both clades, indicating a high prevalence of convergence also at the gene family level. This convergence is coupled with a rich inferred repertoire of multicellularity-related genes in the most recent common ancestor of the Agarico- and Pezizomycotina, consistent with the hypothesis that the coding capacity of ancestral fungal genomes might have promoted the repeated evolution of complex multicellularity. We interpret this repertoire as an indication of evolutionary predisposition of fungal ancestors for evolving complex multicellular fruiting bodies. Our work suggests that evolutionary convergence may happen not only when organisms are closely related or are under similar selection pressures, but also when ancestral genomic repertoires render certain evolutionary trajectories more likely than others, even across large phylogenetic distances. 
    more » « less
  5. Abstract

    The proximate and ultimate mechanisms underlying scaling relationships as well as their evolutionary consequences remain an enigmatic issue in evolutionary biology. Here, I investigate the evolution of wing allometries in the Schizophora, a group of higher Diptera that radiated about 65 million years ago, by studying static allometries in five species using multivariate approaches. Despite the vast ecological diversity observed in contemporary members of the Schizophora and independent evolutionary histories throughout most of the Cenozoic, size‐related changes represent a major contributor to overall variation in wing shape, both within and among species. Static allometries differ between species and sexes, yet multivariate allometries are correlated across species, suggesting a shared developmental programme underlying size‐dependent phenotypic plasticity. Static allometries within species also correlate with evolutionary divergence across 33 different families (belonging to 11 of 13 superfamilies) of the Schizophora. This again points towards a general developmental, genetic or evolutionary mechanism that canalizes or maintains the covariation between shape and size in spite of rapid ecological and morphological diversification during the Cenozoic. I discuss the putative roles of developmental constraints and natural selection in the evolution of wing allometry in the Schizophora.

    more » « less