VENOM composition within snake species can show considerable geographical variation¹, an important consideration because bites by conspecific populations may differ in symptomatology and require different treatments^2–5. The underlying causes of this phenomenon have never been explained. Here we present evidence that the variation in the venom of the pitviper Calloselasma rhodostoma (Serpentes: Viperidae) is closely associated with its diet. We also evaluated other possible causes of geographic variation in venom using partial Mantel tests^6–10 and independent contrasts¹¹, but rejected both contemporary gene flow (estimated from geographical proximity) and the phylogenetic relationships (assessed by analysis of mitochondrial DNA) among populations as important influences upon venom evolution. As the primary function of viperid venom is to immobilize and digest prey^12–14 and prey animals vary in their susceptibility to venom^15,16, we suggest that geographical variation in venom composition reflects natural selection for feeding on local prey.

Snakes are limbless predators, and many species use venom to help overpower relatively large, agile prey. Snake venoms are complex protein mixtures encoded by several multilocus gene families that function synergistically to cause incapacitation. To examine venom evolution, we sequenced and interrogated the genome of a venomous snake, the king cobra (Ophiophagus hannah), and compared it, together with our unique transcriptome, microRNA, and proteome datasets from this species, with data from other vertebrates. In contrast to the platypus, the only other venomous vertebrate with a sequenced genome, we find that snake toxin genes evolve through several distinct co-option mechanisms and exhibit surprisingly variable levels of gene duplication and directional selection that correlate with their functional importance in prey capture. The enigmatic accessory venom gland shows a very different pattern of toxin gene expression from the main venom gland and seems to have recruited toxin-like lectin genes repeatedly for new nontoxic functions. In addition, tissue-specific microRNA analyses suggested the co-option of core genetic regulatory components of the venom secretory system from a pancreatic origin. Although the king cobra is limbless, we recovered coding sequences for all Hox genes involved in amniote limb development, with the exception of Hoxd12. Our results provide a unique view of the origin and evolution of snake venom and reveal multiple genome-level adaptive responses to natural selection in this complex biological weapon system. More generally, they provide insight into mechanisms of protein evolution under strong selection. Significance Snake venoms are toxic protein cocktails used for prey capture. To investigate the evolution of these complex biological weapon systems, we sequenced the genome of a venomous snake, the king cobra, and assessed the composition of venom gland expressed genes, small RNAs, and secreted venom proteins. We show that regulatory components of the venom secretory system may have evolved from a pancreatic origin and that venom toxin genes were co-opted by distinct genomic mechanisms. After co-option, toxin genes important for prey capture have massively expanded by gene duplication and evolved under positive selection, resulting in protein neofunctionalization. This diverse and dramatic venom-related genomic response seemingly occurs in response to a coevolutionary arms race between venomous snakes and their prey.

Significance: The toxic composition of snake venom varies between species. Such variation can have major medical implications for the treatment of human snakebite victims. Venom variation is largely attributed to differences in toxin-encoding genes present in the genome or venom gland of snakes. Here, we demonstrate that mechanisms affecting the transcription, translation, and posttranslational modification of toxins also significantly contribute to the diversity of venom protein composition. Venom variation observed between related snake species is therefore the result of a complex interaction between a variety of genetic and postgenomic factors acting on toxin genes. Ultimately, this variation results in significant differences in venom-induced pathology and lethality and can undermine the efficacy of antivenom therapies used to treat human snakebite victims.

The processes that drive the evolution of snake venom variability, particularly the role of diet, have been a topic of intense recent research interest. Here, we test whether extensive variation in venom composition in the medically important viper genus Echis is associated with shifts in diet. Examination of stomach and hindgut contents revealed extreme variation between the major clades of Echis in the proportion of arthropod prey consumed. The toxicity (median lethal dose, LD ₅₀ ) of representative Echis venoms to a natural scorpion prey species was found to be strongly associated with the degree of arthropod feeding. Mapping the results onto a novel Echis phylogeny generated from nuclear and mitochondrial sequence data revealed two independent instances of coevolution of venom toxicity and diet. Unlike venom LD ₅₀ , the speed with which venoms incapacitated and killed scorpions was not associated with the degree of arthropod feeding. The prey-specific venom toxicity of arthropod-feeding Echis may thus be adaptive primarily by reducing venom expenditure. Overall, our results provide strong evidence that variation in snake venom composition results from adaptive evolution driven by natural selection for different diets, and underscores the need for a multi-faceted, integrative approach to the study of the causes of venom evolution.