Search for: All records

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. 

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. 

Abstract Despite the medical significance to humans and important ecological roles filled by vipers, few high-quality genomic resources exist for these snakes outside of a few genera of pitvipers. Here we sequence, assemble, and annotate the genome of Fea’s Viper (Azemiops feae). This taxon is distributed in East Asia and belongs to a monotypic subfamily, sister to the pitvipers. The newly sequenced genome resulted in a 1.56 Gb assembly, a contig N50 of 1.59 Mb, with 97.6% of the genome assembly in contigs >50 Kb, and a BUSCO completeness of 92.4%. We found that A. feae venom is primarily composed of phospholipase A2 (PLA2) proteins expressed by genes that likely arose from lineage-specific PLA2 gene duplications. Additionally, we show that renin, an enzyme associated with blood pressure regulation in mammals and known from the venoms of two viper species including A. feae, is expressed in the venom gland at comparative levels to known toxins and is present in the venom proteome. The cooption of this gene as a toxin may be more widespread in viperids than currently known. To investigate the historical population demographics of A. feae, we performed coalescent-based analyses and determined that the effective population size has remained stable over the last 100 kyr. This suggests Quaternary glacial cycles likely had minimal influence on the demographic history of A. feae. This newly assembled genome will be an important resource for studying the genomic basis of phenotypic evolution and understanding the diversification of venom toxin gene families. 

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. 

Ecologically divergent selection can lead to the evolution of reproductive isolation through the process of ecological speciation, but the balance of responsible evolutionary forces is often obscured by an inadequate assessment of demographic history and the genetics of traits under selection. Snake venoms have emerged as a system for studying the genetic basis of adaptation because of their genetic tractability and contributions to fitness, and speciation in venomous snakes can be associated with ecological diversification such as dietary shifts and corresponding venom changes. Here, we explored the neurotoxic (type A)–hemotoxic (type B) venom dichotomy and the potential for ecological speciation among Timber Rattlesnake (Crotalus horridus) populations. Previous work identified the genetic basis of this phenotypic difference, enabling us to characterize the roles geography, history, ecology, selection, and chance play in determining when and why new species emerge or are absorbed. We identified significant genetic, proteomic, morphological, and ecological/environmental differences at smaller spatial scales, suggestive of incipient ecological speciation between type A and type B C. horridus. Range-wide analyses, however, rejected the reciprocal monophyly of venom type, indicative of varying intensities of introgression and a lack of reproductive isolation across the range. Given that we have now established the phenotypic distributions and ecological niche models of type A and B populations, genome-wide data are needed and capable of determining whether type A and type B C. horridus represent distinct, reproductively isolated lineages due to incipient ecological speciation or differentiated populations within a single species.