Among extant reptiles only two lineages are known to have evolved venom delivery systems, the advanced snakes and helodermatid lizards (Gila Monster and Beaded Lizard)¹. Evolution of the venom system is thought to underlie the impressive radiation of the advanced snakes (2,500 of 3,000 snake species)^2,3,4,5. In contrast, the lizard venom system is thought to be restricted to just two species and to have evolved independently from the snake venom system¹. Here we report the presence of venom toxins in two additional lizard lineages (Monitor Lizards and Iguania) and show that all lineages possessing toxin-secreting oral glands form a clade, demonstrating a single early origin of the venom system in lizards and snakes. Construction of gland complementary-DNA libraries and phylogenetic analysis of transcripts revealed that nine toxin types are shared between lizards and snakes. Toxinological analyses of venom components from the Lace Monitor Varanus varius showed potent effects on blood pressure and clotting ability, bioactivities associated with a rapid loss of consciousness and extensive bleeding in prey. The iguanian lizard Pogona barbata retains characteristics of the ancestral venom system, namely serial, lobular non-compound venom-secreting glands on both the upper and lower jaws, whereas the advanced snakes and anguimorph lizards (including Monitor Lizards, Gila Monster and Beaded Lizard) have more derived venom systems characterized by the loss of the mandibular (lower) or maxillary (upper) glands. Demonstration that the snakes, iguanians and anguimorphs form a single clade provides overwhelming support for a single, early origin of the venom system in lizards and snakes. These results provide new insights into the evolution of the venom system in squamate reptiles and open new avenues for biomedical research and drug design using hitherto unexplored venom proteins.

Venom has only been recently discovered to be a basal trait of the Anguimorpha lizards. Consequently, very little is known about the timings of toxin recruitment events, venom protein molecular evolution, or even the relative physical diversifications of the venom system itself. A multidisciplinary approach was used to examine the evolution across the full taxonomical range of this ∼130 million-year-old clade. Analysis of cDNA libraries revealed complex venom transcriptomes. Most notably, three new cardioactive peptide toxin types were discovered (celestoxin, cholecystokinin, and YY peptides). The latter two represent additional examples of convergent use of genes in toxic arsenals, both having previously been documented as components of frog skin defensive chemical secretions. Two other novel venom gland-overexpressed modified versions of other protein frameworks were also recovered from the libraries (epididymal secretory protein and ribonuclease). Lectin, hyaluronidase, and veficolin toxin types were sequenced for the first time from lizard venoms and shown to be homologous to the snake venom forms. In contrast, phylogenetic analyses demonstrated that the lizard natriuretic peptide toxins were recruited independently of the form in snake venoms. The de novo evolution of helokinestatin peptide toxin encoding domains within the lizard venom natriuretic gene was revealed to be exclusive to the helodermatid/anguid subclade. New isoforms were sequenced for cysteine-rich secretory protein, kallikrein, and phospholipase A₂ toxins. Venom gland morphological analysis revealed extensive evolutionary tinkering. Anguid glands are characterized by thin capsules and mixed glands, serous at the bottom of the lobule and mucous toward the apex. Twice, independently this arrangement was segregated into specialized serous protein-secreting glands with thick capsules with the mucous lobules now distinct (Heloderma and the Lanthanotus/Varanus clade). The results obtained highlight the importance of utilizing evolution-based search strategies for biodiscovery and emphasize the largely untapped drug design and development potential of lizard venoms.

The binding of compounds to nicotinic acetylcholine receptors is of great interest in biomedical research. However, progress in this area is hampered by the lack of a high-throughput, cost-effective, and taxonomically flexible platform. Current methods are low-throughput, consume large quantities of sample, or are taxonomically limited in which targets can be tested. We describe a novel assay which utilizes a label-free bio-layer interferometry technology, in combination with adapted mimotope peptides, in order to measure ligand binding to the orthosteric site of nicotinic acetylcholine receptor alpha-subunits of diverse organisms. We validated the method by testing the evolutionary patterns of a generalist feeding species (Acanthophis antarcticus), a fish specialist species (Aipysurus laevis), and a snake specialist species (Ophiophagus hannah) for comparative binding to the orthosteric site of fish, amphibian, lizard, snake, bird, marsupial, and rodent alpha-1 nicotinic acetylcholine receptors. Binding patterns corresponded with diet, with the Acanthophis antarcticus not showing bias towards any particular lineage, while Aipysurus laevis showed selectivity for fish, and Ophiophagus hannah a selectivity for snake. To validate the biodiscovery potential of this method, we screened Acanthophis antarcticus and Tropidolaemus wagleri venom for binding to human alpha-1, alpha-2, alpha-3, alpha-4, alpha-5, alpha-6, alpha-7, alpha-9, and alpha-10. While A. antarcticus was broadly potent, T. wagleri showed very strong but selective binding, specifically to the alpha-1 target which would be evolutionarily selected for, as well as the alpha-5 target which is of major interest for drug design and development. Thus, we have shown that our novel method is broadly applicable for studies including evolutionary patterns of venom diversification, predicting potential neurotoxic effects in human envenomed patients, and searches for novel ligands of interest for laboratory tools and in drug design and development.